command: ./cs_solver --trace --mpi *************************************************************** (R) Code_Saturne Version 7.0.2 Copyright (C) 1998-2021 EDF S.A., France revision d029e16a2 build Thu 30 Jun 2022 03:27:14 PM MSK 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: Tue 16 Aug 2022 06:34:00 PM MSK System: Linux 3.10.0-1160.66.1.el7.x86_64 Machine: XeonLinux Processor: model name : Intel(R) Xeon(R) CPU E5-2678 v3 @ 2.50GHz Memory: 64331 MB User: antech (Andrew) Directory: /run/media/antech/HDD Data/Storage/Media/Work/PK-151/Saturne/RESU/Test MPI ranks: 7 (appnum attribute: 0) OpenMP threads: 1 Processors/node: 12 Compilers used for build: C compiler: gcc (GCC) 4.8.5 20150623 (Red Hat 4.8.5-44) C++ compiler: g++ (GCC) 4.8.5 20150623 (Red Hat 4.8.5-44) Fortran compiler: GNU Fortran (GCC) 4.8.5 20150623 (Red Hat 4.8.5-44) MPI version 3.0 (Open MPI 1.8.4) 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: SCOTCH 5.1.12 =============================================================== 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 ) ----------------------------------------------------------- ** 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 = 1 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** STOKES ------ -- Continuous phase: ISTMPF = 1 (time scheme for flow (0: explicit (THETFL = 0 ) (1: std scheme (Saturne 1.0 ) (2: 2nd-order (THETFL = 0.5 ) 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.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 ) ----------------------------------------------------------- ** 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. NTPABS = 0 (Initial time step) NTMABS = -1 (Final time step required) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- run_solver.log files NTLIST = 1 (Output frequency ) -- -- -- --- 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,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: ------------------------------ -5: name: directory: monitoring format: time_plot options: time dependency: fixed mesh output: every time step -1: name: results directory: postprocessing format: EnSight Gold options: separate_meshes time dependency: fixed mesh output: at end -6: name: directory: profiles format: plot options: time dependency: fixed mesh output: at end -7: name: histograms directory: histograms format: histogram options: txt time dependency: fixed mesh output: at end Reading file: mesh_input.csm Finished reading: mesh_input.csm No "partition_input/domain_number_7" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Morton (in bounding box). Number of cells per domain (histogramm): [ 1043822 ; 1043822 ] = 7 Partitioning finished (12.6 s) ---------------------------------------------------------- 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 = 1043822 maximum value = 1043822 1 : [ 1043822 ; 1043822 ] = 7 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 1053172 maximum value = 1074558 1 : [ 1053172 ; 1057449 [ = 2 2 : [ 1057449 ; 1061726 [ = 1 3 : [ 1061726 ; 1066003 [ = 1 4 : [ 1066003 ; 1070280 [ = 2 5 : [ 1070280 ; 1074558 ] = 1 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 9350 maximum value = 30736 1 : [ 9350 ; 13627 [ = 2 2 : [ 13627 ; 17904 [ = 1 3 : [ 17904 ; 22181 [ = 1 4 : [ 22181 ; 26458 [ = 2 5 : [ 26458 ; 30736 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 2054219 maximum value = 2085862 1 : [ 2054219 ; 2060547 [ = 3 2 : [ 2060547 ; 2066876 [ = 0 3 : [ 2066876 ; 2073204 [ = 2 4 : [ 2073204 ; 2079533 [ = 0 5 : [ 2079533 ; 2085862 ] = 2 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 38727 maximum value = 85284 1 : [ 38727 ; 48038 [ = 2 2 : [ 48038 ; 57349 [ = 0 3 : [ 57349 ; 66661 [ = 0 4 : [ 66661 ; 75972 [ = 3 5 : [ 75972 ; 85284 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 2 maximum value = 5 1 : [ 2 ; 3 [ = 4909 2 : [ 3 ; 4 [ = 451348 3 : [ 4 ; 5 ] = 6850497 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 3 maximum value = 6 1 : [ 3 ; 4 [ = 3 2 : [ 4 ; 5 [ = 2 3 : [ 5 ; 6 ] = 2 ---------------------------------------------------------- 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: default Numbering for boundary faces: type: default ---------------------------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 1049150 maximum value = 1069892 1 : [ 1049150 ; 1053298 [ = 2 2 : [ 1053298 ; 1057446 [ = 2 3 : [ 1057446 ; 1061595 [ = 0 4 : [ 1061595 ; 1065743 [ = 1 5 : [ 1065743 ; 1069892 ] = 2 Histogram of volume mesh matrix profile/lines per rank: minimum value = 71534 maximum value = 273113 1 : [ 71534 ; 111849 [ = 2 2 : [ 111849 ; 152165 [ = 0 3 : [ 152165 ; 192481 [ = 1 4 : [ 192481 ; 232797 [ = 0 5 : [ 232797 ; 273113 ] = 4 Mesh coordinates: minimum and maximum X : -1.6867001e+01 1.4139981e+01 Y : 2.5999999e-01 6.0200000e+01 Z : -2.3899930e+00 2.3899763e+00 Mesh Number of cells: 7306754 Number of interior faces: 14382961 Number of boundary faces: 461262 Number of vertices: 1334862 Groups: "AKTOBE_1_1" interior faces: 10136 "AKTOBE_1_2" interior faces: 26426 "AKTOBE_1_3" interior faces: 26135 "AKTOBE_1_4" interior faces: 29873 "AKTOBE_1_5" boundary faces: 63596 "AKTOBE_1_6" interior faces: 1050 "BAYPAS" boundary faces: 22495 "CONFUSOR" boundary faces: 3868 "DIFFUSOR" boundary faces: 15195 "IN" boundary faces: 241 "LUK1" boundary faces: 162 "LUK2" boundary faces: 162 "LUK3" boundary faces: 144 "LUK4" boundary faces: 162 "LUK5" boundary faces: 76 "OUT" interior faces: 56 boundary faces: 56 "PN" boundary faces: 27673 "SHG" boundary faces: 285132 "STENKA" boundary faces: 1167 "TRUBA" boundary faces: 41130 --- Information on the volumes Minimum control volume = 7.4218521e-10 Maximum control volume = 8.1427399e-02 Total volume for the domain = 1.9437676e+03 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 159172 --> 2 % Criterion 3: Least-Squares Gradient Quality: Number of bad cells detected: 2012 --> 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 (2.34 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 7306754 Volume = 1943.8 Center of gravity = (4.0629, 19.712, 1.4533e-05) Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 461262 Surface = 2144.1 Center of gravity = (3.0413, 20.74, 0.00079432) Perimeter = -1 (not computed) ALMAX = 0.12480E+02 (Characteristic length ) ALMAX is the cubic root of the domain volume. -------------------------------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1179E+01 0.1179E+01 molecular_viscos 0.1833E-04 0.1833E-04 ----------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 0 Smooth wall 5 461262 Rough wall 6 0 Symmetry 4 0 Free outlet 3 0 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- ** Field values on cells --------------------- field minimum maximum set mean spatial mean - ---------------- -------------- -------------- -------------- -------------- v Velocity[X] 0 0 0 0 v Velocity[Y] 0 0 0 0 v Velocity[Z] 0 0 0 0 v ǁVelocityǁ 0 0 0 0 v Pressure 0 0 0 0 v CFL 0 0 0 0 v Fourier Number 0 0 0 0 v Total Pressure 1.0132e+05 1.0132e+05 1.0132e+05 1.0133e+05 v volume_forces[X] 0 0 0 0 v volume_forces[Y] 0 0 0 0 v volume_forces[Z] 0 0 0 0 v ǁvolume_forcesǁ 0 0 0 0 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - ----- -------------- -------------- -------------- -------------- v Yplus 0 0 0 0 =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 0 0.000000000E+00 Smooth wall 5 461262 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 0 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 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.75030E-08 0.24508E-06 Friction velocity uk : 0.75030E-08 0.24508E-06 Dimensionless distance yplus : 0.40802E-05 0.13328E-03 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 461262 Total number of wall faces : 461262 ------------------------------------------------------------ ** Writing the main restart file ----------------------------- Start writing Wrote field names and types to checkpoint at iteration 0: checkpoint/main.csc End writing the options Wrote main variables to checkpoint: checkpoint/main.csc End writing ** Writing the auxiliary restart file ---------------------------------- Start writing End writing the dimensions and options End writing the physical properties End writing the time step Wrote boundary condition coefficients to checkpoint: checkpoint/auxiliary.csc End writing ** INFORMATION ON CONVERGENCE -------------------------- ---------------------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual Drift Time residual ---------------------------------------------------------------------------- c Velocity 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 c Velocity[X] 0.00000E+00 c Velocity[Y] 0.00000E+00 c Velocity[Z] 0.00000E+00 c Pressure 0.00000E+00 0 0.00000E+00 -0.00000E+00 0.00000E+00 ---------------------------------------------------------------------------- ** Field values on cells --------------------- field minimum maximum set mean spatial mean - ---------------- -------------- -------------- -------------- -------------- v Velocity[X] 0 0 0 0 v Velocity[Y] 0 0 0 0 v Velocity[Z] 0 0 0 0 v ǁVelocityǁ 0 0 0 0 v Pressure 0 0 0 0 v CFL 0 0 0 0 v Fourier Number 5.7564e-05 34.265 0.27166 0.0015879 v Total Pressure 1.0132e+05 1.0132e+05 1.0132e+05 1.0133e+05 v volume_forces[X] 0 0 0 0 v volume_forces[Y] 0 0 0 0 v volume_forces[Z] 0 0 0 0 v ǁvolume_forcesǁ 0 0 0 0 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - ----- -------------- -------------- -------------- -------------- v Yplus 4.0802e-06 0.00013328 2.0618e-05 6.9871e-05 ** Computed values on cells ------------------------ criterion minimum maximum set mean spatial mean --------------- -------------- -------------- -------------- -------------- Courant/Fourier 5.7564e-05 34.265 0.27166 0.0015879 =============================================================== FINAL STAGE OF THE CALCULATION ============================== =========================================================== =========================================================== END OF CALCULATION ================== =============================================================== Destroying structures and ending computation