command: /mnt/fastdata/mep10yd/scppipe/MxCC_AIR_SEMINI_WSS/RESU/20141110-1605/cs_solver --mpi *************************************************************** (R) Code_Saturne Version 3.0.1 Copyright (C) 1998-2013 EDF S.A., France build Fri 02 Aug 2013 12:42:07 BST MPI version 2.1 (Open MPI 1.6.4) 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: Mon 10 Nov 2014 20:31:32 GMT System: Linux 2.6.32-431.23.3.el6.x86_64 Machine: node084 Processor: model name : Intel(R) Xeon(R) CPU X5650 @ 2.67GHz Memory: 24019 MB User: mep10yd (Y. Duan) Directory: /mnt/fastdata/mep10yd/scppipe/MxCC_AIR_SEMINI_WSS/RESU/20141110-1605 MPI ranks: 36 (appnum attribute: 0) I/O read method: standard input and output, parallel access I/O write method: standard input and output, serial access I/O rank step: 1 Reading metadata from file: "mesh_input" =============================================================== CALCULATION PREPARATION ======================= =========================================================== @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING : WHEN READING INPUT DATA @ ========= @ CHOICE OF TIME-SCHEME @ @ TIME-SCHEME FOR VELOCITY IS SECOND ORDER @ (THETAV = 0.50E+00) @ CERTAIN TERMES ARE HOWEVER FIRST ORDER IN TIME WITH @ THE FOLLOWING SETTINGS: @ @ parameters ISTMPF ISNO2T ISTO2T IROEXT IVIEXT ICPEXT @ Values entered 2 1 0 0 0 0 @ @ computation will go on. @ @ Check the input data given through the User Interface @ or in cs_user_parameters.f90. @ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Geometry NPRFML = 1 (Nb max. family properties ) NFML = 21 (Nb families ) --- Physics NVAR = 5 (Nb variables ) NSCAL = 1 (Nb scalars ) NSCAUS = 1 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) NPROCE = 40 (Nb cell properties ) NPROFA = 2 (Nb internal face properties ) NPROFB = 3 (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 ) -- Continuous phase: RO0 = 0.12250E+01 (Reference density ) VISCL0 = 0.17894E-04 (Ref. molecular dyn. visc. ) CP0 = 0.10064E+04 (Ref. specific heat ) ICP = 0 (> 0: variable CP (usphyv) ) P0 = 0.25331E+07 (Ref. total pressure ) PRED0 = 0.00000E+00 (Ref. reduced pressure ) T0 = 0.63300E+03 (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 ---------- --- Continuous phase: - Commons ITURB = 42 (Turbulence model ) IDEUCH = 0 (0: one-scale model ) (1: two-scale model ) (2: invariant wall function ) IWALLT = 0 (Exch. coeff. correlation ) (0: not activated ) (1: activated ) YPLULI = 0.10880E+02 (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 ) - LES (ITURB = 40, 41, 42) (Sub-grid scale model ) (40 Smagorinsky model ) (41 Dynamic model ) (42 WALE model ) CSMAGO = 0.65000E-01 (Smagorinsky constant ) CWALE = 0.25000E+00 (WALE model constant ) XLESFL = 0.20000E+01 (Filter with in a cell is ) ALES = 0.10000E+01 (written as ) BLES = 0.33333E+00 (XLESFL*(ALES*VOLUME)**(BLES)) IDRIES = 0 (=1 Van Driest damping ) CDRIES = 0.26000E+02 (Van Driest constant ) XLESFD = 0.15000E+01 (Ratio between the explicit ) (filter and LES filter ) (recommended value: 1.5 ) SMAGMX = 0.65000E+00 (Max Smagonsky in the ) (dynamic model case ) IVRTEX = 0 (Use the vortex method ) --- 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 ------------------- --- Continuous phase: IVISSE = 1 (1: accounted for ) ----------------------------------------------------------- ** TIME STEPPING ------------- UNSTEADY ALGORITHM --- Time step parameters IDTVAR = 0 (0 cst; 1,2 var (t, t-space ) IPTLRO = 0 (1: rho-related DT clipping ) COUMAX = 0.10000E+01 (Maximum target CFL ) FOUMAX = 0.10000E+02 (Maximum target Fourier ) VARRDT = 0.10000E+00 (For var. DT, max. increase ) DTMIN = 0.25000E-04 (Minimum time step ) DTMAX = 0.25000E+00 (Maximum time step ) DTREF = 0.25000E-03 (Reference time step ) With a non-constant time step (IDTVAR = 1 or 2), when the value of COUMAX or FOUMAX is negative or zero, the associated time step limitation (for CFL and Fourier respectively) is ignored. --- Frozen velocity field ICCVFG = 0 (1: frozen velocity field ) --- Per-variable properties ------------------------------------ Variable ISTAT CDTVAR ------------------------------------ Pressure 0 0.1000E+01 VelocityX 1 0.1000E+01 VelocityY 1 0.1000E+01 VelocityZ 1 0.1000E+01 Scalar 1 1 0.1000E+01 ---------------------------- ISTAT = 0 ou 1 (1 for unsteady ) CDTVAR > 0 (time step multiplier ) --- Order of base time stepping scheme ISCHTP = 2 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- --------------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV --------------------------------------------------------------------- Pressure 0 1 1 1 0 0.00E+00 0.10E+01 VelocityX 1 1 1 1 1 0.10E+01 0.50E+00 VelocityY 1 1 1 1 1 0.10E+01 0.50E+00 VelocityZ 1 1 1 1 1 0.10E+01 0.50E+00 Scalar 1 1 1 1 1 0 0.10E+01 0.50E+00 ------------------------------------------------------------- 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 ------ IVELCO = 1 (0: segregated solving of the velocity components during the prediction step 1: coupled solving of the velocity components ) 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 ) -- Continuous phase: IREVMC = 0 (Velocity reconstruction mode) RELAXV = 0.10000E+01 for pressure (relaxation) ARAK = 0.10000E+01 (Arakawa factor ) ISTMPF = 2 (time scheme for flow (0: explicit (THETFL = 0 ) (1: std scheme (Saturne 1.0 ) (2: 2nd-order (THETFL = 0.5 ) THETFL = 0.50000E+00 (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.50000E+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 5 0.1000E-04 0.1000E-03 0.0000E+00 VelocityX 100 10 0.1000E-04 0.1000E-03 0.0000E+00 VelocityY 100 10 0.1000E-04 0.1000E-03 0.0000E+00 VelocityZ 100 10 0.1000E-04 0.1000E-03 0.0000E+00 Scalar 1 100 10 0.1000E-04 0.1000E-03 0.0000E+00 ----------------------------------------------------------- ------------------------------------------- Variable IRCFLU IMLIGR CLIMGR ------------------------------------------- Pressure 1 -1 0.1500E+01 VelocityX 1 -1 0.1500E+01 VelocityY 1 -1 0.1500E+01 VelocityZ 1 -1 0.1500E+01 Scalar 1 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 = 30 (Number of moments ) ------------------------------------------------------ IMOM IMOOLD NDTMOM IDFMOM ------------------------------------------------------ 1 -1 3600 0.00000E+00 -1 0 0 0 0 2 -1 3600 0.00000E+00 -1 1 0 0 0 3 -1 3600 0.00000E+00 -1 2 0 0 0 4 -1 3600 0.00000E+00 -1 3 0 0 0 5 -1 3600 0.00000E+00 -1 4 0 0 0 6 -1 3600 0.00000E+00 -1 5 0 0 0 7 -1 3600 0.00000E+00 -1 1 2 0 0 8 -1 3600 0.00000E+00 -1 1 3 0 0 9 -1 3600 0.00000E+00 -1 1 4 0 0 10 -1 3600 0.00000E+00 -1 2 2 0 0 11 -1 3600 0.00000E+00 -1 2 3 0 0 12 -1 3600 0.00000E+00 -1 2 4 0 0 13 -1 3600 0.00000E+00 -1 3 3 0 0 14 -1 3600 0.00000E+00 -1 3 4 0 0 15 -1 3600 0.00000E+00 -1 4 4 0 0 16 -1 3600 0.00000E+00 -1 5 2 0 0 17 -1 3600 0.00000E+00 -1 5 3 0 0 18 -1 3600 0.00000E+00 -1 5 4 0 0 19 -1 3600 0.00000E+00 -1 2 2 2 0 20 -1 3600 0.00000E+00 -1 2 2 3 0 21 -1 3600 0.00000E+00 -1 2 2 4 0 22 -1 3600 0.00000E+00 -1 3 3 2 0 23 -1 3600 0.00000E+00 -1 3 3 3 0 24 -1 3600 0.00000E+00 -1 3 3 4 0 25 -1 3600 0.00000E+00 -1 4 4 2 0 26 -1 3600 0.00000E+00 -1 4 4 3 0 27 -1 3600 0.00000E+00 -1 4 4 4 0 28 -1 3600 0.00000E+00 -1 2 3 4 0 29 -1 3600 0.00000E+00 -1 -2 0 0 0 30 -1 3600 0.00000E+00 -1 -3 0 0 0 ------------------------------------------------------ IMOM = 0 ou > 0 (moment number ) IMOOLD =-1 ou > 0 (old moment corresponding ) ( to calculation restart or ) (-1 if the moment is ) ( reinitialized ) NTDMOM = 0 ou > 0 (moment calculation starting ) (time step number ) TTDMOM = 0 ou > 0 (moment calculation starting ) (time ) IDFMOM = 0 ou > 0 (number of variables of ) (which the moment is composed) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- -------------------------------------------------- Variable IRESOL NITMAX EPSILO IDIRCL -------------------------------------------------- Pressure -1 10000 0.1000E-04 1 VelocityX -1 10000 0.1000E-04 1 VelocityY -1 10000 0.1000E-04 1 VelocityZ -1 10000 0.1000E-04 1 Scalar 1 -1 10000 0.1000E-04 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 number of levels ) -------------------------------------- Variable IMGR NCYMAX NITMGF -------------------------------------- Pressure 0 100 10 VelocityX 0 100 10 VelocityY 0 100 10 VelocityZ 0 100 10 Scalar 1 0 100 10 ------------------------------ IMGR = 0 ou 1 (1: multigrid activated ) NCYMAX = (Max number of cycles ) NITMGF = (Max nb iter on coarsest grid) Multigrid rank merge parameters: mean coarse cells merge threshold: 300 total coarse cells merge threshold: 500 minimum ranks merge threshold: 1 merge stride: 1 Multigrid coarsening type: algebraic, face traveral by criteria ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) --------------------------------------------------------------------- Variable Number ISCSTH IVISLS ITURT VISLS0 SIGMAS --------------------------------------------------------------------- Scalar 1 1 1 0 0 0.2420E-01 0.1000E+01 --------------------------------------------------------------------- ------------------------------------------- Variable Number ISCAVR RVARFL ------------------------------------------- Scalar 1 1 0 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- Scalar 1 1 -1 -0.1000E+13 0.1000E+13 ------------------------------------------------------- ------------------------------------------------------------- For each scalar, the number indicates it's rank in the list of all scalars. User scalars are placed first, from 1 to NSCAUS. Specific physics scalars are placed at the end, from NSCAUS+1 to NSCAPP+NSCAUS=NSCAL. ISCSTH = -1,0, 1 ou 2 (T (C), Passive, T (K) or H ) IVISLS = 0 ou >0 (Viscosity: constant or not ) VISLS0 = >0 (Reference viscosity ) SIGMAS = >0 (Schmidt ) ISCAVR = 0 ou >0 (Associat. scalar if variance) RVARFL = >0 (Rf, cf variance dissipation ) ICLVFL = 0, 1 ou 2 (Variance clipping mode ) SCAMIN = (Min authorized value ) SCAMAX = (Max authorized value ) For variances, SCAMIN is ignored and SCAMAX is used only if ICLVFL = 2 ------------------------------------------------------ Scalar THETSS IVSEXT THETVS ------------------------------------------------------ 1 0.5000E+00 0 0.0000E+00 ------------------------------------------------------ THETSS = (theta for source terms ) ((1+theta).new-theta.old ) IVSEXT = (extrap. total viscosity ) (0: explicit ) (1: n+thetvs with thetvs=1/2 ) (2: n+thetvs with thetvs=1 ) THETVS = (theta for scalar diffusivity ((1+theta).new-theta.old ) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 0 (1: restarted calculation ) ILEAUX = 1 (1: read restart/auxiliary ) IECAUX = 1 (1: write checkpoint/auxiliary) --- 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 = 9600 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Post-processed variables Number Name 2 Pressure 3 VelocityX 4 VelocityY 5 VelocityZ 6 Scalar 1 9 Turb Viscosity 10 CFL 11 Fourier Number 12 Total Pressure 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 MoyTps29 42 MoyTps30 -- -- --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) NCAPT = 10 (Number of probes ) NTHSAV = -1 (Checkpoint frequency ) Number Name Nb. probes (-1: all) 2 Pressure -1 3 VelocityX -1 4 VelocityY -1 5 VelocityZ -1 6 Scalar 1 -1 7 Density -1 9 Turb Viscosity -1 12 Total Pressure -1 41 MoyTps29 -1 42 MoyTps30 -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Pressure 2 3 VelocityX 2 4 VelocityY 2 5 VelocityZ 2 6 Scalar 1 2 7 Density -999 9 Turb Viscosity -999 12 Total Pressure -999 41 MoyTps29 -999 42 MoyTps30 -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 0 (Force exerted at boundary) ipstyp = 0 (Yplus at boundary) ipstft = 0 (Thermal flow at boundary) ipsttb = 0 (Temperature at boundary) ipstnu = 0 (Nusselt at boundary) ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated ) NALINF = 0 (Fluid initialization iterations) ----------------------------------------------------------- Reading file: mesh_input Finished reading: mesh_input No "partition_input/domain_number_36" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Hilbert (in bounding box). Number of cells per domain (histogramm): [ 26350 ; 26350 ] = 36 Writing file: partition_output/domain_number_36 Finished writing: partition_output/domain_number_36 Partitioning finished (13.2 s) Postprocessing output writers: ------------------------------ -1: name: results_us directory: postprocessing format: EnSight Gold options: time dependency: fixed mesh output: every 400 time steps ---------------------------------------------------------- Halo construction with standard neighborhood ============================================ Face interfaces creation Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Histogram of the number of cells per rank: minimum value = 26350 maximum value = 26350 1 : [ 26350 ; 26350 ] = 36 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 31718 maximum value = 35107 1 : [ 31718 ; 32395 [ = 8 2 : [ 32395 ; 33073 [ = 10 3 : [ 33073 ; 33751 [ = 12 4 : [ 33751 ; 34429 [ = 2 5 : [ 34429 ; 35107 ] = 4 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 5368 maximum value = 8757 1 : [ 5368 ; 6045 [ = 8 2 : [ 6045 ; 6723 [ = 10 3 : [ 6723 ; 7401 [ = 12 4 : [ 7401 ; 8079 [ = 2 5 : [ 8079 ; 8757 ] = 4 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 81664 maximum value = 83866 1 : [ 81664 ; 82104 [ = 14 2 : [ 82104 ; 82544 [ = 10 3 : [ 82544 ; 82985 [ = 6 4 : [ 82985 ; 83425 [ = 0 5 : [ 83425 ; 83866 ] = 6 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 150 maximum value = 1471 1 : [ 150 ; 414 [ = 2 2 : [ 414 ; 678 [ = 6 3 : [ 678 ; 942 [ = 12 4 : [ 942 ; 1206 [ = 6 5 : [ 1206 ; 1471 ] = 10 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 4 maximum value = 6 1 : [ 4 ; 5 [ = 240 2 : [ 5 ; 6 ] = 948360 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 7 maximum value = 22 1 : [ 7 ; 10 [ = 14 2 : [ 10 ; 13 [ = 12 3 : [ 13 ; 16 [ = 6 4 : [ 16 ; 19 [ = 2 5 : [ 19 ; 22 ] = 2 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: Checking mesh renumbering for threads: ------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 27826 maximum value = 34805 1 : [ 27826 ; 29221 [ = 6 2 : [ 29221 ; 30617 [ = 5 3 : [ 30617 ; 32013 [ = 12 4 : [ 32013 ; 33409 [ = 10 5 : [ 33409 ; 34805 ] = 3 Histogram of volume mesh matrix profile/lines per rank: minimum value = 4432 maximum value = 10352 1 : [ 4432 ; 5616 [ = 16 2 : [ 5616 ; 6800 [ = 10 3 : [ 6800 ; 7984 [ = 6 4 : [ 7984 ; 9168 [ = 2 5 : [ 9168 ; 10352 ] = 2 Mesh coordinates: minimum and maximum X : -1.5699999e-02 1.5699999e-02 Y : -1.5699999e-02 1.5699999e-02 Z : -3.1400001e-01 0.0000000e+00 Mesh Number of cells: 948600 Number of interior faces: 2829540 Number of boundary faces: 32520 Number of vertices: 965107 Groups: "INLET1" boundary faces: 5100 "INLET1A" interior faces: 5100 "OUTLET1" boundary faces: 5100 "PART_1_1" interior faces: 7905 "PART_1_13" interior faces: 7905 "PART_1_14" interior faces: 7905 "PART_1_15" interior faces: 7905 "PART_1_16" interior faces: 7905 "PART_1_17" interior faces: 4185 "PART_1_18" interior faces: 2790 "PART_1_19" interior faces: 2790 "PART_1_2" interior faces: 7905 "PART_1_20" interior faces: 1395 "PART_1_3" interior faces: 7905 "PART_1_4" interior faces: 7905 "PART_1_5" interior faces: 4185 "PART_1_6" interior faces: 2790 "PART_1_7" interior faces: 2790 "PART_1_8" interior faces: 1395 "WALL1" boundary faces: 22320 --- Information on the volumes Minimum control volume = 9.4711485e-11 Maximum control volume = 7.0640867e-10 Total volume for the domain = 2.4304147e-04 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 0 --> 0 % Criterion 3: Least-Squares Gradient Quality: Number of bad cells detected: 15300 --> 2 % Criterion 4: Cells Volume Ratio: Number of bad cells detected: 0 --> 0 % Criterion 5: Guilt by Association: Number of bad cells detected: 0 --> 0 % Warning: -------- Mesh quality issue has been detected The mesh should be re-considered using the listed criteria. The calculation will run but the solution quality may be degraded... Computing geometric quantities (1.56 s) ALMAX = 0.62579E-01 (Characteristic length ) ALMAX is the cubic root of the domain volume. inisem: volume box = 1.28673036993585477E-003 inisem: [xmin, xmax] = -3.06551061179079452E-002 3.06551061179079452E-002 inisem: [ymin, ymax] = -3.06551061179079452E-002 3.06551061179079452E-002 inisem: [zmin, zmax] = -0.32815591931343091 1.41559140034951264E-002 ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Pressure 0.0000E+00 0.0000E+00 VelocityX -0.3706E+01 0.2597E+01 VelocityY -0.3173E+01 0.3075E+01 VelocityZ -0.6678E+01 0.6542E+00 Scalar 1 0.0000E+00 0.0000E+00 --------------------------------- 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 MoyTps29 0.0000E+00 0.0000E+00 MoyTps30 0.0000E+00 0.0000E+00 --------------------------------- Local Time Step 0.2500E-03 0.2500E-03 --------------------------------- Time averages (sum over the time-steps) ------------------------------------------------------------ Average Min. value Max. value Uniform/Variable in space ------------------------------------------------------------ 1 0.0000E+00 0.0000E+00 Uniform 2 0.0000E+00 0.0000E+00 Uniform 3 0.0000E+00 0.0000E+00 Uniform 4 0.0000E+00 0.0000E+00 Uniform 5 0.0000E+00 0.0000E+00 Uniform 6 0.0000E+00 0.0000E+00 Uniform 7 0.0000E+00 0.0000E+00 Uniform 8 0.0000E+00 0.0000E+00 Uniform 9 0.0000E+00 0.0000E+00 Uniform 10 0.0000E+00 0.0000E+00 Uniform 11 0.0000E+00 0.0000E+00 Uniform 12 0.0000E+00 0.0000E+00 Uniform 13 0.0000E+00 0.0000E+00 Uniform 14 0.0000E+00 0.0000E+00 Uniform 15 0.0000E+00 0.0000E+00 Uniform 16 0.0000E+00 0.0000E+00 Uniform 17 0.0000E+00 0.0000E+00 Uniform 18 0.0000E+00 0.0000E+00 Uniform 19 0.0000E+00 0.0000E+00 Uniform 20 0.0000E+00 0.0000E+00 Uniform 21 0.0000E+00 0.0000E+00 Uniform 22 0.0000E+00 0.0000E+00 Uniform 23 0.0000E+00 0.0000E+00 Uniform 24 0.0000E+00 0.0000E+00 Uniform 25 0.0000E+00 0.0000E+00 Uniform 26 0.0000E+00 0.0000E+00 Uniform 27 0.0000E+00 0.0000E+00 Uniform 28 0.0000E+00 0.0000E+00 Uniform 29 0.0000E+00 0.0000E+00 Uniform 30 0.0000E+00 0.0000E+00 Uniform ------------------------------------------------------------ ------------------------------------------------------------- ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.250000000E-03 TIME STEP NUMBER 1 ============================================================= ------------------------------------------------------------ INITIALISATIONS =============== REINITIALISATION DE LA PRESSION A L'ITERATION 1 ------------------------------------------------------------ COMPUTATION OF PHYSICAL QUANTITIES ================================== _iterative_vector_gradient: isweep = 6, residue norm: 2.398524e-06, norm: 3.289265e+05, var_num = 2 ----------------------------------------- Property Min. value Max. value ----------------------------------------- Density 0.1225E+01 0.1225E+01 Laminar Viscosit 0.1789E-04 0.1789E-04 Turb Viscosity 0.5841E-07 0.8904E-04 ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- Scalar 1 1 0.2420E-01 0.2420E-01 ----------------------------------------------- ------------------------------------------------------------ COMPUTATION OF CFL, FOURIER AND VARIABLE DT =========================================== COURANT MAX= 0.0000E+00 IN -0.1136E-01 -0.1078E-01 -0.1578E+00 COURANT MIN= 0.0000E+00 IN -0.1136E-01 -0.1078E-01 -0.1578E+00 FOURIER MAX= 0.3935E+01 IN 0.1559E-01 -0.4083E-03 -0.2912E+00 FOURIER MIN= 0.4580E-01 IN -0.3220E-03 -0.3220E-03 -0.2895E+00 COU/FOU MAX= 0.3935E+01 IN 0.1559E-01 -0.4083E-03 -0.2912E+00 COU/FOU MIN= 0.4580E-01 IN -0.3220E-03 -0.3220E-03 -0.2895E+00 ------------------------------------------------------------ SETTING UP THE BOUNDARY CONDITIONS ================================== ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 5100 Smooth wall 5 22320 Rough wall 6 0 Symmetry 4 0 Free outlet 3 5100 Undefined 1 0 ------------------------------------------------------------------------- Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = 0.32239E-03 -0.32239E-03 -0.31400E+00 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 5100 0.000000000E+00 Smooth wall 5 22320 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 5100 0.000000000E+00 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- SIGSEGV signal (forbidden memory area access) intercepted! Call stack: 1: 0x40749b (cs_solver) 2: 0x2b144b56b973 (libsaturne.so.0) 3: 0x2b144b670961 (libsaturne.so.0) 4: 0x2b144b5553fd (libsaturne.so.0) 5: 0x2b144b536925 (libsaturne.so.0) 6: 0x2b144b536a85 (libsaturne.so.0) 7: 0x328b61ed5d <__libc_start_main+0xfd> (libc.so.6) 8: 0x4054d9 <> (cs_solver) End of stack