command: /home/jiahuan/Code_Saturne/4.0.5/code_saturne-4.0.5/arch/Linux_x86_64/libexec/code_saturne/cs_solver --param setup.xml *************************************************************** (R) Code_Saturne Version 4.0.5 Copyright (C) 1998-2016 EDF S.A., France revision 4.0.5 build Tue 19 Jul 2016 22:35:38 BST MPI version 3.0 (MPICH 3.0.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: Sun 02 Sep 2018 19:31:58 BST System: Linux 3.19.0-42-generic Machine: jiahuan-VirtualBox Processor: model name : Intel(R) Core(TM) i7-6560U CPU @ 2.20GHz Memory: 10907 MB User: jiahuan (jiahuan) Directory: /home/jiahuan/yanyan/code_saturne/mvg/case1/RESU/20180902-1931 I/O read method: standard input and output, serial access I/O write method: standard input and output, serial access External libraries for partitioning: ParMETIS 4.0.3 PT-SCOTCH 6.0.3 Reading metadata from file: "mesh_input" =============================================================== CALCULATION PREPARATION ======================= =========================================================== No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 4 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) NPROCE = 6 (Nb cell properties ) ----------------------------------------------------------- ** PHYSICAL PROPERTIES ------------------- GX = 0.00000E+00 (Gravity x component ) GY = 0.00000E+00 (Gravity y component ) GZ = 0.00000E+00 (Gravity z component ) ICORIO = 0 (Coriolis source terms ) -- Continuous phase: RO0 = 0.12240E+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.29315E+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 ----------------------------------------------------------- ** HOMOGENEOUS MIXTURE MODEL FOR CAVITATION ---------------------------------------- ICAVIT = -1 (-1: single phase flow ) ( 0: no vap./cond. model ) ( 1: Merkle's model ) ----------------------------------------------------------- ** THERMAL MODEL ------------- --- Continuous phase: - Commons ITHERM = 0 (0: no thermal model ) (1: temperature ) (2: enthalpy ) (3: total energy ) ITPSCL = 0 (0: none ) (1: temperature in Kelvin ) (2: temperature in Celsius ) ISCALT = -1 (Thermal scalar number ) ** TURBULENCE ---------- --- Continuous phase: - Commons ITURB = 0 (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.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) ) --- 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 = 0 (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.20000E-05 (Minimum time step ) DTMAX = 0.20000E-01 (Maximum time step ) DTREF = 0.20000E-04 (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 ------------------------------------ Velocity 1 0.1000E+01 Pressure 0 0.1000E+01 ---------------------------- ISTAT = 0 ou 1 (1 for unsteady ) CDTVAR > 0 (time step multiplier ) --- Order of base time stepping scheme ISCHTP = 1 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- --------------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV --------------------------------------------------------------------- Velocity 1 1 1 1 0 0.10E+01 0.10E+01 Pressure 0 1 1 1 0 0.00E+00 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 ------ IDILAT = 1 (1 : without unsteady term in the continuity equation 2 : with unsteady term in the continuity equation) IPOROS = 0 (0 : without porous media 1 : with porous media) 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 = 2 (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 = 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 ------------------------------------------------------------------- Velocity 100 10 0.1000E-04 0.1000E-06 0.0000E+00 Pressure 100 10 0.1000E-04 0.1000E-06 0.0000E+00 ----------------------------------------------------------- ------------------------------------------- Variable IRCFLU IMLIGR CLIMGR ------------------------------------------- Velocity 1 -1 0.1500E+01 Pressure 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 ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-07 1 Pressure 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 ) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 0 (1: restarted calculation ) ILEAUX = 0 (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 = 4300000 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 500 (Checkpoint frequency ) --- Post-processed variables Velocity Pressure CourantNb FourierNb total_pressure -- -- --- Probe history files NTHIST = -1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) NCAPT = 0 (Number of probes ) NTHSAV = -1 (Checkpoint frequency ) Number Name Nb. probes (-1: all) 2 Velocity[X] -1 3 Velocity[Y] -1 4 Velocity[Z] -1 5 Pressure -1 7 LamVisc -1 8 CourantNb -1 9 FourierNb -1 10 total_pressure -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Velocity 0 5 Pressure 0 8 CourantNb -999 9 FourierNb -999 10 total_pressure -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 1 (Force exerted by the fluid on the boundary) ipstft = 1 (Thermal flux at boundary) ipsttb = 0 (Temperature at boundary) ipstnu = 0 (Dimensionless thermal flux at boundary) ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated ) NALINF = 0 (Fluid initialization iterations) IFLXMW = 0 (ALE mass flux computation 0: thanks to vertices 1: thanks to mesh velocity) ----------------------------------------------------------- Postprocessing output writers: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: binary time dependency: fixed mesh output: every 500 time steps Reading file: mesh_input Finished reading: mesh_input ------------------------------------------------------- Joining number 1: Selection criteria: "all[]" Parameters for the joining operation: Shortest incident edge fraction: 0.01000 Maximum angle between joined face planes: 25.00000 Advanced joining parameters: Verbosity level: 1 Visualization level: 1 Deepest level reachable in tree building: 30 Max boxes by leaf: 25 Max ratio of linked boxes / init. boxes: 5.00000 Max ratio of boxes for distribution: 2.00000 Merge step tolerance multiplier: 1.00000 Pre-merge factor: 0.05000 Tolerance computation mode: 1 Intersection computation mode: 1 Max. number of equiv. breaks: 500 Max. number of subfaces by face: 200 Before joining Number of cells: 3288362 Number of interior faces: 9792824 Number of boundary faces: 144524 Number of vertices: 3361044 Global number of boundary faces selected for joining: 144524 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 0 Vertex-Vertex intersections: 0 Other intersections: 0 Joining operation is null. Edge intersections and vertex creation done. Mesh cleaning done for degenerate faces. Global number of cleaned interior faces: 0 Global number of cleaned border faces: 0 After joining Number of cells: 3288362 Number of interior faces: 9792824 Number of boundary faces: 144524 Number of vertices: 3361044 Joining 1 completed (2.37 s). All joining operations successfully finished: Wall clock time: 2.37 Number of cells: 3288362 ---------------------------------------------------------- Number of cells + halo cells: 3288362 ---------------------------------------------------------- Number of interior faces: 9792824 ---------------------------------------------------------- Number of boundary faces: 144524 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 6 1 : [ 3 ; 4 [ = 8 2 : [ 4 ; 5 [ = 1954 3 : [ 5 ; 6 ] = 3286400 ---------------------------------------------------------- Number of cells: 3288362 ---------------------------------------------------------- Number of cells + halo cells: 3288362 ---------------------------------------------------------- Number of interior faces: 9792824 ---------------------------------------------------------- Number of boundary faces: 144524 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 6 1 : [ 3 ; 4 [ = 8 2 : [ 4 ; 5 [ = 1954 3 : [ 5 ; 6 ] = 3286400 ---------------------------------------------------------- Renumbering mesh: renumbering for cells: pre-numbering: none cells adjacent to ghost cells last: no numbering: none renumbering for interior faces: cell adjacency pre-ordering: lowest id first faces adjacent to ghost cells last: no numbering: none renumbering for boundary faces: numbering: none ---------------------------------------------------------- Matrix bandwidth for volume mesh : 1796557 Matrix profile/lines for volume mesh : 33902 Mesh coordinates: minimum and maximum X : 0.0000000e+00 4.0000000e+00 Y : 0.0000000e+00 6.0000002e-01 Z : 0.0000000e+00 6.0000002e-01 Mesh Number of cells: 3288362 Number of interior faces: 9792824 Number of boundary faces: 144524 Number of vertices: 3361044 Groups: "INLET-4" boundary faces: 15198 "OUTLET-5" boundary faces: 15198 "SYMMETRY-6" boundary faces: 86800 "WALL-7" boundary faces: 27328 --- Information on the volumes Minimum control volume = 1.3724682e-11 Maximum control volume = 3.3143960e-05 Total volume for the domain = 1.4399999e+00 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 2790 --> 0 % Criterion 3: Least-Squares Gradient Quality: Number of bad cells detected: 0 --> 0 % 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 (4.31 s)