command: ./cs_solver *************************************************************** (R) Code_Saturne Version 6.0.4 Copyright (C) 1998-2019 EDF S.A., France revision 6.0.4-4743b4795 build Sun Aug 9 14:46:39 2020 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 Oct 16 14:30:50 2020 System: Linux 2.6.32-279.14.1.el6.x86_64 (CentOS release 6.3 (Final)) Machine: compute-0-10.local Processor: model name : AMD Opteron(tm) Processor 6386 SE Memory: 129190 MB User: a_rasam () Directory: /RESU/20201016-1429 OpenMP threads: 2 Processors/node: 64 External libraries for partitioning: METIS 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 ------------------------------------------------------- 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. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 228036 Vertex-Vertex intersections: 77354 Other intersections: 150682 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 152002 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 (6.99 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. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 2765470 Vertex-Vertex intersections: 2765470 Other intersections: 0 Joining operation is conforming. Edge intersections and vertex creation done. Pre-merge for 172346 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 (236 s). All joining operations successfully finished: Wall clock time: 243 ---------------------------------------------------------- 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 Number of cells: 3153668 ---------------------------------------------------------- Number of cells + halo cells: 3495632 ---------------------------------------------------------- Number of ghost cells: 341964 in the standard neighborhood: 341964 ---------------------------------------------------------- Number of interior faces: 9628042 ---------------------------------------------------------- Number of boundary faces: 45884 ---------------------------------------------------------- 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 ---------------------------------------------------------- Synchronizing cell families Writing file: mesh_output Finished writing: mesh_output 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 number of halo-adjacent elements: 3153668 Numbering for interior faces: type: threads number of threads: 2 number of exclusive groups: 7 number of elements in group 0: 9425898 number of elements in group 1: 193436 number of elements in group 2: 6847 number of elements in group 3: 509 number of elements in group 4: 809 number of elements in group 5: 488 number of elements in group 6: 55 estimated thread imbalance: 0.026 Numbering for boundary faces: type: threads number of threads: 2 number of exclusive groups: 1 number of elements in group 0: 45884 estimated thread imbalance: 0.000 ---------------------------------------------------------- Matrix bandwidth for volume mesh : 3153668 Matrix profile/lines for volume mesh : 225088 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 (6.52 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 ================