Need help with Salome Meca and Code Satune

Hi all,

I'm a total beginner at CFD/FEA stuff and I need some help with Salome Meca. I'm trying to export a mesh I've created that has certain faces grouped together (as an inlet, outlet, and symmetry groups) and I want to export the mesh to run in Code Saturne. How do I let Code Saturne know which faces correspond to which boundary conditions? The wizard doesn't seem to detect anything.

Hi claws,

Thank you for your help. I've run into a few more problems now though, with code_saturne telling me I have an error in my boundary conditions. Basically all I'm trying to simulate is laminar flow with a standard object in the middle (in my case a cone). Basically it's a box with one face an inlet, the opposite face an outlet, and the rest as symmetry. I created the mesh using Netgen 1D-2D-3D algorithm. I checked the error file and all it said was there was a problem with my boundary conditions, and to look at the end of my calculation log, which I cannot find. Could you point me in the right direction?

Hi Claws,

Thanks for the advice so far. I've found the listing file and I've tried making the mesh finer so the boundary condition error doesn't occur anymore, but that did not work. I'm not sure what the problem is. I've attached the listing file here, could you help me out by telling me what is actually going on, and how to avoid this error? Thank you!



command:
/home/zixin/tmp_Saturne/CONESTUDY2.CONE2.11232228/cs_solver --param cone2

***************************************************************

(R)
Code_Saturne

Version 2.0.2


Copyright (C) 1998-2011 EDF S.A., France

build Sat 01 Oct 2011 12:32:08 PM EDT
MPI version 2.1 (Open MPI)


The Code_Saturne CFD tool is free software;
you can redistribute it and/or modify it under the terms
of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License,
or (at your option) any later version.

The Code_Saturne CFD tool is distributed in the hope that
it will be useful, but WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A
PARTICULAR PURPOSE. See the GNU General Public License
for more details.

***************************************************************


Local case configuration:

Date: Sat 23 Nov 2013 10:28:35 PM EST
System: Linux 2.6.32-52-generic
Machine: zixin-laptop
Processor: Intel(R) Core(TM) i7 CPU Q 720 @ 1.60GHz
Memory: 5917 MB
User: zixin (Zixin Nie)
Directory: /home/zixin/tmp_Saturne/CONESTUDY2.CONE2.11232228
MPI ranks: 1

===============================================================

CALCULATION PREPARATION
=======================


===========================================================




Reading file: preprocessor_output

No error detected during the data verification
(interface, usini1 and others).


===========================================================

CALCULATION PARAMETERS SUMMARY
==============================

---

** DIMENSIONS

---

--- Geometry
NPRFML = 1 (Nb max. family properties )
NFML = 4 (Nb families )

--- Physics
NPHAS = 1 (Nb phases )
NVAR = 4 (Nb variables )
NSCAL = 0 (Nb scalars )
NSCAUS = 0 (Nb user scalars )
NSCAPP = 0 (Nb specific physics scalars )
NPROCE = 6 (Nb cell properties )
NPROFA = 1 (Nb internal face properties )
NPROFB = 2 (Nb boundary face properties )

---

** PHYSICAL PROPERTIES

---

GX = 0.00000E+00 (Gravity x component )
GY = 0.00000E+00 (Gravity y component )
GZ = 0.00000E+00 (Gravity z component )

OMEGAX = 0.00000E+00 (Rotation vector x component )
OMEGAY = 0.00000E+00 (Rotation vector y component )
OMEGAZ = 0.00000E+00 (Rotation vector z component )
ICORIO = 0 (Coriolis source terms )

-- Phase: 1

RO0 = 0.11786E+01 (Reference density )
VISCL0 = 0.18300E-04 (Ref. molecular dyn. visc. )
CP0 = 0.10172E+04 (Ref. specific heat )
ICP = 0 (> 0: variable CP (usphyv) )
P0 = 0.10132E+06 (Ref. total pressure )
PRED0 = 0.00000E+00 (Ref. reduced pressure )
T0 = 0.00000E+00 (Ref. temperature )

IROVAR = 0 (Density constant(0) or not(1)
IVIVAR = 0 (Molec. visc cst.(0) or not(1)

Initial reference point for pressure
XYZP0 = 0.00000E+00 0.00000E+00 0.00000E+00

---

** TURBULENCE

---

--- Phase: 1

- Commons
ITURB = 0 (Turbulence model )
IDEUCH = 0 (0: one-scale model )
(1: two-scale model )
(2: invariant wall function )
YPLULI = 0.23810E+01 (Limit Y+ )
ILOGPO = 1 (0: power law (forbidden for
k-epsilon) )
(1: one-scale log law )
IGRHOK = 0 (1: computed Grad(rho k) )
ISCALT = -1 (Temperature salar number )

--- Constants

- Commons
XKAPPA = 0.42000E+00 (Von Karman constant )
CSTLOG = 0.52000E+01 (U+=Log(y+)/kappa +CSTLOG )
APOW = 0.83000E+01 (U+=APOW (y+)**BPOW (W&W law))
BPOW = 0.14286E+00 (U+=APOW (y+)**BPOW (W&W law))

---

** SECONDARY VISCOSITY

---

--- Phase: 1
IVISSE = 1 (1: accounted for )

---

** TIME STEPPING

---

STEADY ALGORITHM

--- Global parameters

IDTVAR = -1 (-1: steady algorithm )
RELXST = 0.90000E+00 (Reference relaxation coeff. )


--- Frozen velocity field

ICCVFG = 0 (1: frozen velocity field )

--- Per variable relaxation coefficient

---

RELAXV variable

---

Pressure 0.1000E+00
Velocity 0.9000E+00
Velocity 0.9000E+00
Velocity 0.9000E+00

---

RELAXV = [0.,1.] (relaxation coefficient )

---

** CONVECTION - DIFFUSION

---

---

Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV

---

Pressure 0 1 1 1 0 0.00E+00 0.10E+01
Velocity 1 1 1 1 0 0.10E+01 0.10E+01
Velocity 1 1 1 1 0 0.10E+01 0.10E+01
Velocity 1 1 1 1 0 0.10E+01 0.10E+01

---

ICONV = 0 ou 1 (1 for convection active )
IDIFF = 0 ou 1 (1 for total diffusion active)
IDIFFT = 0 ou 1 (1 for turbulent diff. active)
ISCHCV = 0 ou 1 (SOLU or CD )
ISSTPC = 0 ou 1 (1: no slope test )
BLENCV = [0.;1.] (1-upwind proportion )
THETAV = [0.;1.] (0.5 Crank-Nicolson/AB )
(theta for convection- )
(diffusion terms uses )
((1-theta).old+theta.new )

---

** STOKES

---

IPHYDR = 0 (1: account for explicit
balance between pressure
gradient, gravity source
terms, and head losses )
ICALHY = 0 (1: compute hydrastatic
pressure for Dirichlet
conditions for pressure
on outlet )
IPRCO = 1 (1: pressure-continuity )
IPUCOU = 0 (1: reinforced U-P coupling )
NTERUP = 1 (n: n sweeps on navsto for
velocity/pressure coupling )

-- Phase: 1

IREVMC = 0 (Velocity reconstruction mode)

ARAK = 0.10000E+01 (Arakawa factor )

ISTMPF = 1 (time scheme for flow
(0: explicit (THETFL = 0 )
(1: std scheme (Saturne 1.0 )
(2: 2nd-order (THETFL = 0.5 )
THETFL = -0.99900E+03 (theta for mass flow )
IROEXT = 0 (density extrapolation
(0: explicit
(1: n+thetro with thetro=1/2
(2: n+thetro with thetro=1
THETRO = 0.00000E+00 (theta for density
((1+theta).new-theta.old
IVIEXT = 0 (total viscosity extrapolation
(0: explicit
(1: n+thetvi with thetro=1/2
(2: n+thetvi with thetro=1
THETVI = 0.00000E+00 (theta for total viscosity
((1+theta).new-theta.old
ICPEXT = 0 (specific heat extrapolation
(0: explicit
(1: n+thetcp with thetro=1/2
(2: n+thetcp with thetro=1
THETCP = 0.00000E+00 (specific heat theta-scheme
((1+theta).new-theta.old
THETSN = 0.00000E+00 (Nav-Stokes S.T. theta scheme)
((1+theta).new-theta.old
THETST = 0.00000E+00 (Turbulence S.T. theta-scheme)
((1+theta).new-theta.old
EPSUP = 0.10000E-04 (Velocity/pressure coupling
stop test )

---

** GRADIENTS CALCULATION

---

IMRGRA = 0 (Reconstruction mode )
ANOMAX = 0.78540E+00 (Non-ortho angle: limit for )
(least squares ext. neighbors)

---

Variable NSWRGR NSWRSM EPSRGR EPSRSM EXTRAG

---

Pressure 100 2 0.1000E-04 0.1000E-07 0.0000E+00
Velocity 100 1 0.1000E-04 0.1000E-07 0.0000E+00
Velocity 100 1 0.1000E-04 0.1000E-07 0.0000E+00
Velocity 100 1 0.1000E-04 0.1000E-07 0.0000E+00

---

---

Variable IRCFLU IMLIGR CLIMGR

---

Pressure 1 -1 0.1500E+01
Velocity 1 -1 0.1500E+01
Velocity 1 -1 0.1500E+01
Velocity 1 -1 0.1500E+01

---

NSWRGR = (nb sweep gradient reconstr. )
NSWRSM = (nb sweep rhs reconstrcution )
EPSRGR = (grad. reconstruction prec. )
EPSRSM = (rhs reconstruction prec. )
EXTRAG = [0.;1.] (gradients extrapolation )
IRCFLU = 0 ou 1 (flow reconstruction )
IMLIGR = < 0, 0 ou 1 (gradient limitation method )
CLIMGR = > 1 ou 1 (gradient limitation coeff. )

---

** FACE INTERPOLATION

---

IMVISF = 0 (0 arithmetic )

---

** CALCULATION OF TEMPORAL MEANS (MOMENTS)

---

NBMOMT = 0 (Number of moments )

---

** BASE ITERATIVE SOLVERS

---

---

Variable IRESOL NITMAX EPSILO IDIRCL

---

Pressure -1 10000 0.1000E-07 1
Velocity -1 10000 0.1000E-07 1
Velocity -1 10000 0.1000E-07 1
Velocity -1 10000 0.1000E-07 1

---

IRESOL = -1 (automatic solver choice )
IPOL*1000 + 0 (p conjuguate gradient )
1 (Jacobi )
IPOL*1000 + 2 (bicgstab )
avec IPOL (preconditioning degree )
NITMAX = (max number of iterations )
EPSILO = (resolution precision )
IDIRCL = 0 ou 1 (shift diagonal if
ISTAT=0 and no Dirichlet )

---

** MULTIGRID

---

NCEGRM = 30 (Max nb cells coarsest grid )
NGRMAX = 25 (Max numeber of levels )

---

Variable IMGR NCYMAX NITMGF

---

Pressure 1 100 10
Velocity 0 100 10
Velocity 0 100 10
Velocity 0 100 10

---

IMGR = 0 ou 1 (1: multigrid activated )
NCYMAX = (Max number of cycles )
NITMGF = (Max nb iter on coarsest grid)

---

** CALCULATION MANAGEMENT

---

--- Restarted calculation
ISUITE = 0 (1: restarted calculuation )
ILEAUX = 1 (1: also read suiamx )
IECAUX = 1 (1: also write suiavx )

suiamx and suiavx are the auxiliary restart files.

--- Calculation time
The numbering of time steps and the measure of simulated
physical time are absolute values, and not values
relative to the current calculation.

INPDT0 = 0 (1: 0 time step calcuation )
NTMABS = 10 (Final time step required )

--- CPU time margin
TMARUS = -0.10000E+01 (CPU time margin before stop )

---

** INPUT-OUTPUT

---

--- Restart file
NTSUIT = 0 (Checkpoint frequency )

--- EnSight, MED, or CGNS output
NTCHR = -1 (Output frequency )

Number Name
2 Pressure
3 VelocityX
4 VelocityY
5 VelocityZ
10 total_pressure
-- --

--- Probe history files
NTHIST = -1 (Output frequency )
NCAPT = 0 (Number of probes )
NTHSAV = 0 (Checkpoint frequency )

Number Name Nb. probes (-1: all)
2 Pressure -1
3 VelocityX -1
4 VelocityY -1
5 VelocityZ -1
6 Density -1
7 LamVisc -1
10 total_pressure -1
-- -- --

--- Log files
NTLIST = 1 (Output frequency )

Number Name IWARNI verbosity level
(-999: not applicable)

2 Pressure 0
3 VelocityX 0
4 VelocityY 0
5 VelocityZ 0
10 total_pressure -999
-- -- --

--- Additional post-processing variables
IPSTDV = 210 (Product of the following
values based on activation )
IPSTYP = 2 (Yplus on boundary )
IPSTCL = 3 (Variables on boundary )
IPSTFT = 5 (Thermal flow on boundary )
IPSTFO = 7 (Force exerted on boundary )
and IPSTDV = 1 (No additional output )

---

** FILES

---

--- Standard files

GEOMETRY MAIN RESTART AUX RESTART STOP
Name geomet suiamo suiamx ficstp
Unit 10 12

MAIN CHECKPOINT AUX CHECKPOINT
Name suiava suiavx

--- Radiative files

DATA RESTART JANAF CHECKPOINT
Name dp_tch rayamo JANAF rayava
Unit 25 25

--- Lagrangian files

--- RESTART --- -- CHECKPOINT - ---- POST ----
CALCUL. STAT. CALCUL. STAT. LOG HIST.
Name lagamo lasamo lagava lasava listla
Unit 80 81 82

--- Other files for Lagrangian module

Unit 50 51 52 53 54
55 56 57 58 59
60 61 62 63 64
65 66 67 68

--- Developper files for probe history
IMPHIS(1) IMPHIS(2)
Unit 30 31

--- User files for probe history
16 files
Name and unit
ush001 ush002 ush003 ush004 ush005 ush006
ush007 ush008 ush009 ush010 ush011 ush012
ush013 ush014 ush015 ush016
33 34 35 36 37 38
39 40 41 42 43 44
45 46 47 48

--- Free user files
10 files
Name and unit
usrf01 usrf02 usrf03 usrf04 usrf05 usrf06
usrf07 usrf08 usrf09 usrf10
70 71 72 73 74 75
76 77 78 79

--- Preprocessor files
Format: 1=Read preprocessor_output file
0=Standalone solver

Format 1

---

** ALE METHOD (MOVING MESH)

---

IALE = 0 (1: activated )
NALINF = -999 (Fluid initialization
iterations)

---

Finished reading: preprocessor_output

Elapsed time for reading preprocessor_output: 0.0275 s

Number of cells: 230374

---

Number of cells + halo cells: 230374

---

Number of interior faces: 452729

---

Number of boundary faces: 16038

---

Histogram of the number of interior faces per cell:

minimum value = 1
maximum value = 4

1 : [ 1 ; 2 [ = 1
2 : [ 2 ; 3 [ = 387
3 : [ 3 ; 4 ] = 229986

---

Renumbering mesh:
Mesh
Number of cells: 230374
Number of interior faces: 452729
Number of boundary faces: 16038
Number of vertices: 43489

Computing geometric quantities (0.136 s)

Mesh coordinates: minimum and maximum
X : -1.0000000e+01 1.0000000e+01
Y : -1.0000000e+01 1.0000000e+01
Z : -1.0000000e+01 1.0000000e+01

--- Main Fortran work arrays:
longia = 6911220 (Number of integers)
longra = 27644880 (Number of reals)
(4 bytes/integer, 8 bytes/real)

--- User Fortran work arrays:
nituse = 0 (Number of integers)
nrtuse = 0 (Number of reals)

--- Information on the volumes
Minimum control volume = 0.627480898E-06
Maximum control volume = 0.127208284E+00
Total volume for the domain = 0.799994193E+04

--- Phase: 1
ALMAX = 0.19940E+02 (Characteristic length )
ALMAX is the cubic root of the domain volume.

---

** VARIABLES INITIALIZATION

---

---

Variable Min. value Max. value

---

Pressure 0.0000E+00 0.0000E+00
Velocity 0.0000E+00 0.0000E+00
Velocity 0.0000E+00 0.0000E+00
Velocity 0.0000E+00 0.0000E+00

---

---

ALL PHASES : 1D-WALL THERMAL MODULE NOT ACTIVATED
NFPT1D = 0

---

PHASE 1 : HEAD LOSS TREATMENT NOT ACTIVATED
NCEPDC = 0

---

PHASE 1 : MASS SOURCE TERMS TREATMENT NOT ACTIVATED
NCETSM = 0

---

===============================================================



MAIN CALCULATION
================


===============================================================




===============================================================


INSTANT 0.100000000E+01 TIME STEP NUMBER 1
=============================================================


--- Phase: 1

---

Property Min. value Max. value

---

Density 0.1179E+01 0.1179E+01
LamVisc 0.1830E-04 0.1830E-04

---

First face with boundary condition definition error
(out of 618)
has boundary condition type 0, center (-0.0208047, 0.373085, 0.114424)



/opt/saturne-2.0.2/installer/ncs-2.0.2/src/base/cs_boundary_conditions.c:429: Fatal error.


Some boundary condition definitions are incomplete or incorrect.

For details, read the end of the calculation log,
or visualize the error postprocessing output.


Call stack:
1: 0x7f685b08e452 <cs_boundary_conditions_error+0x852> (libsaturne.so.0)
2: 0x7f685b24f106 <typecl_+0x236> (libsaturne.so.0)
3: 0x7f685b07d081 <condli_+0x1301> (libsaturne.so.0)
4: 0x7f685b2354a1 <tridim_+0x6931> (libsaturne.so.0)
5: 0x7f685b06c875 <caltri_+0x5085> (libsaturne.so.0)
6: 0x7f685b0479e3 <cs_run+0x843> (libsaturne.so.0)
7: 0x7f685b047cc5 <main+0x1f5> (libsaturne.so.0)
8: 0x7f68577e2c8d <__libc_start_main+0xfd> (libc.so.6)
9: 0x40a359 <> (cs_solver)
End of stack

Hello,

Making your mesh finer should not change anything.

Either you forgot to specify boundary conditions on some face groups, or some faces are not in the expected groups.

Did you visualize the boundary faces that have errors as recommended in the "listing" file ?
To do this, just open the error.ensight or error.med with ParaView or ParaVis.

It will show you on which faces you have issues. You may also check the groups in "check mesh" mode, and see if any faces are missing from a group.

Regards,

Yvan