command: ./cs_solver --mpi --param new.xml *************************************************************** (R) Code_Saturne Version 5.0.4 Copyright (C) 1998-2017 EDF S.A., France revision 5.0.4 build Mon Oct 30 13:31:29 2017 MPI version 3.1 (MPICH 3.2) 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: Fri Nov 10 18:40:35 2017 System: Linux 3.0.101-0.46.1_1.0502.8871-cray_ari_c Machine: nid00009 Processor: model name : Intel(R) Xeon(R) CPU E5-2697 v2 @ 2.70GHz Memory: 64523 MB Directory: /fs2/e556/e556/yanyan/LES/5.0TSR3/case1/RESU/20171110-1837 MPI ranks: 192 (appnum attribute: 0) OpenMP threads: 1 Processors/node: 1 I/O read method: collective MPI-IO (explicit offsets) I/O write method: collective MPI-IO (explicit offsets) I/O rank step: 1 Reading metadata from file: "mesh_input/mesh_01" Reading metadata from file: "mesh_input/mesh_02" =============================================================== CALCULATION PREPARATION ======================= =========================================================== @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING : WHEN READING INPUT DATA @ ========= @ CHOICE OF TIME-SCHEME @ @ TIME-SCHEME FOR VELOCITY IS SECOND ORDER @ (THETAV = 0.50E+00) @ CERTAIN TERMS 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 ---------- --- Physics NVAR = 4 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) ----------------------------------------------------------- ** VOF METHOD ---------------------------------------- IVOFMT = -1 ( -1: inactive ) ( 0: active ) ** HOMOGENEOUS MIXTURE MODEL FOR CAVITATION ---------------------------------------- ICAVIT = -1 (-1: single phase flow ) ( 0: no vap./cond. model ) ( 1: Merkle's model ) ----------------------------------------------------------- ** TIME STEPPING ------------- --- 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 = 2 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** STOKES ------ -- Phase continue : 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 ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-04 1 Pressure 0.1000E-04 1 ------------------------------------ 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 = 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 = 10000 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- Log files NTLIST = -1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) Velocity 0 Pressure 0 TurbVisc -999 CourantNb -999 FourierNb -999 total_pressure -999 Yplus -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 1 (Force exerted by the fluid on the boundary) ipstyp = 1 (y+ at boundary) ipsttp = 0 (T+ at boundary) ipstft = 1 (Thermal flux 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: transient connectivity output: every 100 time steps and at calculation end -5: name: directory: monitoring format: time_plot options: time dependency: fixed mesh output: every 1 time steps -6: name: directory: profiles format: plot options: time dependency: fixed mesh output: at calculation end Reading file: mesh_input/mesh_01 Finished reading: mesh_input/mesh_01 Reading file: mesh_input/mesh_02 Finished reading: mesh_input/mesh_02 No "partition_input/domain_number_192" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Morton (in bounding box). Number of cells per domain (histogramm): [ 1537 ; 1538 ] = 192 Partitioning finished (1.58 s) ------------------------------------------------------- Joining number 1: Selection criteria: "internal" Parameters for the joining operation: Shortest incident edge fraction: 0.10000 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: 295236 Number of interior faces: 588888 Number of boundary faces: 593640 Number of vertices: 593636 Global number of boundary faces selected for joining: 0 Element selection successfully done. Stop joining algorithm: no face selected... After joining Number of cells: 295236 Number of interior faces: 588888 Number of boundary faces: 593640 Number of vertices: 593636 Joining 1 completed (0.00491 s). All joining operations successfully finished: Wall clock time: 0.007 ---------------------------------------------------------- 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 = 1537 maximum value = 1538 1 : [ 1537 ; 1538 ] = 192 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 1685 maximum value = 1976 1 : [ 1685 ; 1743 [ = 77 2 : [ 1743 ; 1801 [ = 81 3 : [ 1801 ; 1859 [ = 23 4 : [ 1859 ; 1917 [ = 9 5 : [ 1917 ; 1976 ] = 2 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 147 maximum value = 439 1 : [ 147 ; 205 [ = 76 2 : [ 205 ; 263 [ = 81 3 : [ 263 ; 322 [ = 25 4 : [ 322 ; 380 [ = 8 5 : [ 380 ; 439 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 3094 maximum value = 3316 1 : [ 3094 ; 3138 [ = 3 2 : [ 3138 ; 3182 [ = 16 3 : [ 3182 ; 3227 [ = 111 4 : [ 3227 ; 3271 [ = 57 5 : [ 3271 ; 3316 ] = 5 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 3074 maximum value = 3255 1 : [ 3074 ; 3110 [ = 163 2 : [ 3110 ; 3146 [ = 13 3 : [ 3146 ; 3182 [ = 11 4 : [ 3182 ; 3218 [ = 1 5 : [ 3218 ; 3255 ] = 4 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 4 1 : [ 3 ; 4 ] = 295236 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 3 maximum value = 13 1 : [ 3 ; 5 [ = 9 2 : [ 5 ; 7 [ = 54 3 : [ 7 ; 9 [ = 65 4 : [ 9 ; 11 [ = 39 5 : [ 11 ; 13 ] = 25 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Writing file: mesh_output Finished writing: mesh_output 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 ---------------------------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 1537 maximum value = 1960 1 : [ 1537 ; 1621 [ = 39 2 : [ 1621 ; 1706 [ = 71 3 : [ 1706 ; 1790 [ = 66 4 : [ 1790 ; 1875 [ = 13 5 : [ 1875 ; 1960 ] = 3 Histogram of volume mesh matrix profile/lines per rank: minimum value = 105 maximum value = 328 1 : [ 105 ; 149 [ = 71 2 : [ 149 ; 194 [ = 85 3 : [ 194 ; 238 [ = 30 4 : [ 238 ; 283 [ = 4 5 : [ 283 ; 328 ] = 2 Mesh coordinates: minimum and maximum X : -8.8981520e+00 8.7594056e+00 Y : -8.8006768e+00 8.8734923e+00 Z : 0.0000000e+00 1.0000000e-02 Mesh Number of cells: 295236 Number of interior faces: 588888 Number of boundary faces: 593640 Number of vertices: 593636 Groups: "a" cells: 13246 "b" cells: 13727 "farfield" boundary faces: 590472 "group-1" cells: 45780 "group-10" cells: 9912 "group-11" cells: 1682 "group-12" cells: 11600 "group-13" cells: 841 "group-14" cells: 11600 "group-2" cells: 19293 "group-3" cells: 35721 "group-4" cells: 11600 "group-5" cells: 9912 "group-6" cells: 9381 "group-7" cells: 45780 "group-8" cells: 45780 "group-9" cells: 9381 "inlet" boundary faces: 358 "interface1" boundary faces: 729 "interface2" boundary faces: 729 "internal" interior faces: 2939 "outlet" boundary faces: 371 "wall1" boundary faces: 327 "wall2" boundary faces: 327 "wall3" boundary faces: 327 --- Information on the volumes Minimum control volume = 1.4730753e-11 Maximum control volume = 3.3560305e-04 Total volume for the domain = 2.4503063e+00 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: 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 % Computing geometric quantities (0.00389 s) ALMAX = 0.13482E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Velocity 0.0000E+00 0.0000E+00 Velocity 0.0000E+00 0.0000E+00 Velocity 0.0000E+00 0.0000E+00 Pressure 0.0000E+00 0.0000E+00 --------------------------------- dt 0.1000E-02 0.1000E-02 --------------------------------- ------------------------------------------------------------- ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.100000000E-02 TIME STEP NUMBER 1 ============================================================= ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1225E+01 0.1225E+01 molecular_viscos 0.1830E-04 0.1830E-04 turbulent_viscos 0.0000E+00 0.0000E+00 ----------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 358 Smooth wall 5 981 Rough wall 6 0 Symmetry 4 590472 Free outlet 3 371 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = 0.86098E+01 -0.15783E+01 0.50000E-02 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 358 0.000000000E+00 Smooth wall 5 981 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 590472 0.000000000E+00 Free outlet 3 371 0.000000000E+00 Free inlet 14 0 0.000000000E+00 Convective inlet 16 0 0.000000000E+00 Free surface 15 0 0.000000000E+00 Undefined 1 0 0.000000000E+00 ---------------------------------------------------------------