command: ./cs_solver --mpi *************************************************************** (R) Code_Saturne Version 6.0.5 Copyright (C) 1998-2019 EDF S.A., France revision 6.0.5-3e9e9e94c build Sun Aug 9 14:46:39 2020 MPI version 3.0 (Open MPI 1.10.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: - System: Linux 4.4.0-190-generic (Ubuntu 16.04.1 LTS) Machine: Ubuntu1604-001 Processor: model name : AMD Opteron(tm) Processor 6134 Memory: 7983 MB User: zangeneh (${Username}) Directory: RESU/20201013-0419 MPI ranks: 4 (appnum attribute: 0) OpenMP threads: 4 Processors/node: 2 I/O read method: collective MPI-IO (explicit offsets) I/O write method: collective MPI-IO (explicit offsets) I/O rank step: 1 External libraries for partitioning: PT-SCOTCH 6.0.9 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) @ SOME 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 ) ----------------------------------------------------------- ** HOMOGENEOUS MIXTURE MODEL VoF ---------------------------------------- IVOFMT = 0 ( 0 : disabled ) ( > 0: enabled ) ----------------------------------------------------------- ** 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 ------ -- Continuous phase: 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 ) THETVI = 0.00000E+00 (theta for total viscosity ((1+theta).new-theta.old THETCP = -0.99900E+03 (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 ) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 1 (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. NTPABS = 115073 (Initial time step) NTMABS = 500000 (Final time step required) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 5 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- run_solver.log files NTLIST = 10 (Output frequency ) -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 7 (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,2: CDO ) 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: separate_meshes time dependency: fixed mesh output: 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 -7: name: histograms directory: histograms format: histogram options: txt time dependency: fixed mesh output: at calculation end Reading file: mesh_input Finished reading: mesh_input No "partition_input/domain_number_4" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Morton (in bounding box). Number of cells per domain (histogramm): [ 788417 ; 788417 ] = 4 Partitioning finished (27.9 s) ------------------------------------------------------- Joining number 1: Selection criteria: "INTERFACE_LAM or INTERFACE_TURB" 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: 3153668 Number of interior faces: 9242654 Number of boundary faces: 436700 Number of vertices: 3373386 Global number of boundary faces selected for joining: 48852 Element selection successfully done. Global min/max. tolerance: Glob. Num. | Tolerance | Coordinates 82012 | 0.000003 | 7.0423021115e-02 7.2714890893e-02 0.0000000000e+00 | ORI 202161 | 0.000045 | 5.7484254713e-01 -1.3995147804e-01 4.4999998063e-03 | ORI Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 281625 Vertex-Vertex intersections: 101480 Other intersections: 180145 Joining operation is non-conforming. Global number of new vertices to create: 150682 Equivalences broken for 0 edges. Edge intersections and vertex creation done. Pre-merge for 300117 global element couples. Max size of a merge set of vertices: 2 Excessive transitivity for 0 set(s) of vertices. Edge removed for 0 faces (global). Join mesh cleaning done. Edge removed for 0 faces (global). Join mesh cleaning done. Merge vertices and mesh update done. Global number of faces after splitting: 43424 Global configuration after the joining operation: Global number of border faces to add: 0 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: 3153668 Number of interior faces: 9286078 Number of boundary faces: 387848 Number of vertices: 3367954 Joining 1 completed (4.9 s). ------------------------------------------------------- Joining number 2: Periodicity type: translation Transformation matrix: 1 0 0 0 0 1 0 0 0 0 1 0.018 Selection criteria: "PER1 or PER2" 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: 3153668 Number of interior faces: 9286078 Number of boundary faces: 387848 Number of vertices: 3367954 Global number of boundary faces selected for joining: 341964 Element selection successfully done. Global min/max. tolerance: Glob. Num. | Tolerance | Coordinates 29115 | 0.000000 | 3.5999424779e-02 -2.4882761289e-02 0.0000000000e+00 | ORI 1 | 0.000450 | -3.6840002179e+00 -4.2008739725e+00 0.0000000000e+00 | ORI Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 5198322 Vertex-Vertex intersections: 5198322 Other intersections: 0 Joining operation is conforming. Edge intersections and vertex creation done. Pre-merge for 335620 global element couples. No need to merge vertices. Edge removed for 0 faces (global). Join mesh cleaning done. Edge removed for 0 faces (global). Join mesh cleaning done. Merge vertices and mesh update done. Global number of faces after splitting: 512946 Global configuration after the joining operation: Global number of border faces to add: 0 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: 3153668 Number of interior faces: 9628042 Number of boundary faces: 45884 Number of vertices: 3367954 Joining 2 completed (134 s). All joining operations successfully finished: Wall clock time: 139 ---------------------------------------------------------- Composing periodicities Halo construction with extended neighborhood ============================================ Face interfaces creation Definition of periodic vertices Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Extended neighborhood structures definition Histogram of the number of cells per rank: minimum value = 788417 maximum value = 788417 1 : [ 788417 ; 788417 ] = 4 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 897017 maximum value = 1062861 1 : [ 897017 ; 930185 [ = 1 2 : [ 930185 ; 963354 [ = 0 3 : [ 963354 ; 996523 [ = 1 4 : [ 996523 ; 1029692 [ = 0 5 : [ 1029692 ; 1062861 ] = 2 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 108600 maximum value = 274444 1 : [ 108600 ; 141768 [ = 1 2 : [ 141768 ; 174937 [ = 0 3 : [ 174937 ; 208106 [ = 1 4 : [ 208106 ; 241275 [ = 0 5 : [ 241275 ; 274444 ] = 2 Histogram of the number of ghost cells in the standard neighborhood per rank: minimum value = 106744 maximum value = 268720 1 : [ 106744 ; 139139 [ = 1 2 : [ 139139 ; 171534 [ = 0 3 : [ 171534 ; 203929 [ = 1 4 : [ 203929 ; 236324 [ = 0 5 : [ 236324 ; 268720 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 2417807 maximum value = 2501696 1 : [ 2417807 ; 2434584 [ = 1 2 : [ 2434584 ; 2451362 [ = 0 3 : [ 2451362 ; 2468140 [ = 1 4 : [ 2468140 ; 2484918 [ = 0 5 : [ 2484918 ; 2501696 ] = 2 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 9980 maximum value = 12448 1 : [ 9980 ; 10473 [ = 1 2 : [ 10473 ; 10967 [ = 0 3 : [ 10967 ; 11460 [ = 1 4 : [ 11460 ; 11954 [ = 0 5 : [ 11954 ; 12448 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 4 maximum value = 13 1 : [ 4 ; 5 [ = 45868 2 : [ 5 ; 7 [ = 3102372 3 : [ 7 ; 9 [ = 0 4 : [ 9 ; 11 [ = 0 5 : [ 11 ; 13 ] = 5428 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 4 maximum value = 4 1 : [ 4 ; 4 ] = 4 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Writing file: mesh_output Finished writing: mesh_output ---------------------------------------------------------- Ignoring periodicity for graph-based partitioning. Partitioning 3153668 cells to 4 domains on 4 ranks (SCOTCH_dgraphPart). wall-clock time: 19.741627 s Number of cells per domain (histogramm): [ 779100 ; 780884 [ = 1 [ 780884 ; 782669 [ = 0 [ 782669 ; 784454 [ = 0 [ 784454 ; 786239 [ = 0 [ 786239 ; 788024 [ = 1 [ 788024 ; 789808 [ = 0 [ 789808 ; 791593 [ = 1 [ 791593 ; 793378 [ = 0 [ 793378 ; 795163 [ = 0 [ 795163 ; 796948 ] = 1 Writing file: partition_output/domain_number_4 Finished writing: partition_output/domain_number_4 Partitioning finished (21.8 s) ---------------------------------------------------------- Composing periodicities Halo construction with extended neighborhood ============================================ Face interfaces creation Definition of periodic vertices Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Extended neighborhood structures definition Histogram of the number of cells per rank: minimum value = 779100 maximum value = 796948 1 : [ 779100 ; 782669 [ = 1 2 : [ 782669 ; 786239 [ = 0 3 : [ 786239 ; 789808 [ = 1 4 : [ 789808 ; 793378 [ = 1 5 : [ 793378 ; 796948 ] = 1 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 874087 maximum value = 886990 1 : [ 874087 ; 876667 [ = 1 2 : [ 876667 ; 879248 [ = 0 3 : [ 879248 ; 881828 [ = 1 4 : [ 881828 ; 884409 [ = 1 5 : [ 884409 ; 886990 ] = 1 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 89726 maximum value = 94987 1 : [ 89726 ; 90778 [ = 2 2 : [ 90778 ; 91830 [ = 0 3 : [ 91830 ; 92882 [ = 0 4 : [ 92882 ; 93934 [ = 0 5 : [ 93934 ; 94987 ] = 2 Histogram of the number of ghost cells in the standard neighborhood per rank: minimum value = 87934 maximum value = 93512 1 : [ 87934 ; 89049 [ = 2 2 : [ 89049 ; 90165 [ = 0 3 : [ 90165 ; 91280 [ = 0 4 : [ 91280 ; 92396 [ = 0 5 : [ 92396 ; 93512 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 2383305 maximum value = 2434115 1 : [ 2383305 ; 2393467 [ = 1 2 : [ 2393467 ; 2403629 [ = 0 3 : [ 2403629 ; 2413791 [ = 2 4 : [ 2413791 ; 2423953 [ = 0 5 : [ 2423953 ; 2434115 ] = 1 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 10616 maximum value = 12420 1 : [ 10616 ; 10976 [ = 2 2 : [ 10976 ; 11337 [ = 0 3 : [ 11337 ; 11698 [ = 0 4 : [ 11698 ; 12059 [ = 1 5 : [ 12059 ; 12420 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 4 maximum value = 13 1 : [ 4 ; 5 [ = 45868 2 : [ 5 ; 7 [ = 3102372 3 : [ 7 ; 9 [ = 0 4 : [ 9 ; 11 [ = 0 5 : [ 11 ; 13 ] = 5428 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 3 maximum value = 3 1 : [ 3 ; 3 ] = 4 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: renumbering for cells: pre-numbering: no renumbering cells adjacent to ghost cells last: no numbering: no renumbering renumbering for interior faces: cell adjacency pre-ordering: lowest id first faces adjacent to ghost cells last: no numbering: multipass renumbering for boundary faces: numbering: no shared cell across threads renumbering for vertices: numbering: no renumbering Numbering for cells: type: default Numbering for interior faces: type: threads minimum maximum mean number of threads: 4 4 4 number of exclusive groups: 3 4 3 number of elements in group 0: 2282196 2330097 2306337 number of elements in group 1: 100100 104023 102424 number of elements in group 2: 1009 1128 1087 number of elements in group 3: 294 499 392 estimated thread imbalance: 0.093 0.100 0.097 Numbering for boundary faces: type: threads minimum maximum mean number of threads: 4 4 4 number of exclusive groups: 1 1 1 number of elements in group 0: 10616 12420 11471 estimated thread imbalance: 0.000 0.000 0.000 ---------------------------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 870853 maximum value = 883381 1 : [ 870853 ; 873358 [ = 1 2 : [ 873358 ; 875864 [ = 1 3 : [ 875864 ; 878369 [ = 1 4 : [ 878369 ; 880875 [ = 0 5 : [ 880875 ; 883381 ] = 1 Histogram of volume mesh matrix profile/lines per rank: minimum value = 58989 maximum value = 60188 1 : [ 58989 ; 59228 [ = 3 2 : [ 59228 ; 59468 [ = 0 3 : [ 59468 ; 59708 [ = 0 4 : [ 59708 ; 59948 [ = 0 5 : [ 59948 ; 60188 ] = 1 Mesh coordinates: minimum and maximum X : -5.3999822e+00 7.2001574e+00 Y : -6.0000000e+00 6.0000000e+00 Z : -3.8743020e-10 1.8000000e-02 Mesh Number of cells: 3153668 Number of interior faces: 9628042 Number of boundary faces: 45884 Number of vertices: 3367954 Periodic faces (which are also interior faces): Periodicity 1: 170982 face couples Groups: "INLET" boundary faces: 5428 "INTERFACE_LAM" interior faces: 43424 "INTERFACE_TURB" interior faces: 43424 "OUTLET" boundary faces: 4648 "PER1" interior faces: 170982 "PER2" interior faces: 170982 "WALL" boundary faces: 35808 --- Information on the volumes Minimum control volume = 3.0682996e-13 Maximum control volume = 1.1589449e-04 Total volume for the domain = 2.4427596e+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 % Criteria 6: User Specific Tag: Number of bad cells detected: 3153668 --> 100 % Warning: -------- Mesh quality issue based on user criteria has been detected The mesh should be re-considered using the listed criteria. Computing geometric quantities (3.45 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 3153668 Volume = 2.4427596e+00 Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 45884 Surface = 8.1521397e-01 Perimeter = -1 (not computed) ALMAX = 0.13468E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. ------------------------------------------------------------- ----------------------------------------------------------- ** READING MAIN AND AUXILIARY RESTART FILES ---------------------------------------- READING THE MAIN RESTART FILE Start reading Reading dimensions complete Reading the previous time step number (restarting computation) NTPABS = 115073 Reading the previous time step number (restarting computation) TTPABS = 0.2000E+05 Reading options complete Read variables from restart: restart/main Reading complete READING THE AUXILIARY RESTART FILE Start reading Fin de la lecture des options restart/auxiliary: section "ref_presstot01" not present. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING: WHEN READING THE AUXILIARY RESTART FILE @ ======= @ Lecture des proprietes physiques @ @ It was not possible to read some values from the @ auxiliary restart file. @ They will be initialized by the default values. @ @ This situation can occur when the restart file @ originates from a run realised with different @ options or when the file is damaged. @ @ The run will continue... @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ Fin de la lecture des proprietes physiques Fin de la lecture du pas de temps Fin de la lecture des flux de masse Fin de la lecture des conditions aux limites restart/auxiliary: section "time_moments:sizes" not present. End reading ----------------------------------------------------------- ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Velocity -0.8197E+02 0.1714E+03 Velocity -0.6621E+02 0.1136E+03 Velocity -0.6521E+02 0.5616E+02 Pressure -0.7414E+04 0.6003E+04 --------------------------------- dt 0.1200E-05 0.1200E-05 --------------------------------- ------------------------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 5428 Smooth wall 5 35808 Rough wall 6 0 Symmetry 4 0 Free outlet 3 4648 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- =============================================================== MAIN CALCULATION ================