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: Fri 12 Aug 2022 12:23:06 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/20220812-1222 MPI ranks: 12 (appnum attribute: 0) OpenMP threads: 1 Processors/node: 1 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 = 7 (Nb variables ) NSCAL = 1 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 1 (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 k 1 0.1000E+01 epsilon 1 0.1000E+01 TempC 1 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.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 ) ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) ----------------------------------------------------- Variable Number ITURT ----------------------------------------------------- TempC 1 0 ------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ TempC 1 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- TempC 1 -1 0.0000E+00 0.6000E+03 ------------------------------------------------------- ------------------------------------------------------------- For each scalar, the number indicates it's rank in the list of all scalars. User scalars are placed first, from 1 to NSCAUS. Specific physics scalars are placed at the end, from NSCAUS+1 to NSCAPP+NSCAUS=NSCAL. is_temerature = 0 or 1 (use Cp or not ) VISLS0 = >0 (Reference viscosity ) SIGMAS = >0 (Schmidt ) RVARFL = >0 (Rf, cf variance dissipation ) ICLVFL = 0, 1 or 2 (Variance clipping mode ) SCAMIN = (Min authorized value ) SCAMAX = (Max authorized value ) For variances, SCAMIN is ignored and SCAMAX is used only if ICLVFL = 2 ------------------------------------------------------ Scalar THETSS THETVS ------------------------------------------------------ 1 0.0000E+00 0.0000E+00 ------------------------------------------------------ THETSS = (theta for source terms ) ((1+theta).new-theta.old ) THETVS = (theta for scalar diffusivity ((1+theta).new-theta.old ) ----------------------------------------------------------- ** 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 = 10000000 (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: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: separate_meshes time dependency: fixed mesh output: at end -5: name: directory: monitoring format: time_plot options: CSV time dependency: fixed mesh output: every time step -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_12" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Morton (in bounding box). Number of cells per domain (histogramm): [ 608882 ; 608883 ] = 12 Partitioning finished (9.04 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 = 608882 maximum value = 608883 1 : [ 608882 ; 608883 ] = 12 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 619442 maximum value = 643139 1 : [ 619442 ; 624181 [ = 4 2 : [ 624181 ; 628920 [ = 1 3 : [ 628920 ; 633660 [ = 5 4 : [ 633660 ; 638399 [ = 0 5 : [ 638399 ; 643139 ] = 2 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 10559 maximum value = 34257 1 : [ 10559 ; 15298 [ = 4 2 : [ 15298 ; 20038 [ = 1 3 : [ 20038 ; 24777 [ = 5 4 : [ 24777 ; 29517 [ = 0 5 : [ 29517 ; 34257 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 1201669 maximum value = 1234571 1 : [ 1201669 ; 1208249 [ = 5 2 : [ 1208249 ; 1214829 [ = 3 3 : [ 1214829 ; 1221410 [ = 3 4 : [ 1221410 ; 1227990 [ = 0 5 : [ 1227990 ; 1234571 ] = 1 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 11873 maximum value = 51536 1 : [ 11873 ; 19805 [ = 1 2 : [ 19805 ; 27738 [ = 1 3 : [ 27738 ; 35670 [ = 2 4 : [ 35670 ; 43603 [ = 3 5 : [ 43603 ; 51536 ] = 5 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 2 maximum value = 4 1 : [ 2 ; 3 [ = 4925 2 : [ 3 ; 4 ] = 7301661 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 2 maximum value = 7 1 : [ 2 ; 3 [ = 1 2 : [ 3 ; 4 [ = 2 3 : [ 4 ; 5 [ = 1 4 : [ 5 ; 6 [ = 2 5 : [ 6 ; 7 ] = 6 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Extrusion: 56 boundary faces selected. 168 cells added. ---------------------------------------------------------- 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 = 608882 maximum value = 608961 1 : [ 608882 ; 608897 [ = 9 2 : [ 608897 ; 608913 [ = 1 3 : [ 608913 ; 608929 [ = 0 4 : [ 608929 ; 608945 [ = 1 5 : [ 608945 ; 608961 ] = 1 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 619538 maximum value = 643139 1 : [ 619538 ; 624258 [ = 4 2 : [ 624258 ; 628978 [ = 1 3 : [ 628978 ; 633698 [ = 5 4 : [ 633698 ; 638418 [ = 0 5 : [ 638418 ; 643139 ] = 2 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 10577 maximum value = 34257 1 : [ 10577 ; 15313 [ = 4 2 : [ 15313 ; 20049 [ = 1 3 : [ 20049 ; 24785 [ = 5 4 : [ 24785 ; 29521 [ = 0 5 : [ 29521 ; 34257 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 1201858 maximum value = 1234571 1 : [ 1201858 ; 1208400 [ = 5 2 : [ 1208400 ; 1214943 [ = 3 3 : [ 1214943 ; 1221485 [ = 3 4 : [ 1221485 ; 1228028 [ = 0 5 : [ 1228028 ; 1234571 ] = 1 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 11873 maximum value = 51536 1 : [ 11873 ; 19805 [ = 1 2 : [ 19805 ; 27738 [ = 1 3 : [ 27738 ; 35670 [ = 2 4 : [ 35670 ; 43603 [ = 3 5 : [ 43603 ; 51536 ] = 5 ---------------------------------------------------------- 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 = 2 maximum value = 7 1 : [ 2 ; 3 [ = 1 2 : [ 3 ; 4 [ = 2 3 : [ 4 ; 5 [ = 1 4 : [ 5 ; 6 [ = 2 5 : [ 6 ; 7 ] = 6 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Writing file: mesh_output.csm Finished writing: mesh_output.csm 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 = 616885 maximum value = 636566 1 : [ 616885 ; 620821 [ = 4 2 : [ 620821 ; 624757 [ = 2 3 : [ 624757 ; 628693 [ = 1 4 : [ 628693 ; 632629 [ = 2 5 : [ 632629 ; 636566 ] = 3 Histogram of volume mesh matrix profile/lines per rank: minimum value = 44860 maximum value = 165653 1 : [ 44860 ; 69018 [ = 3 2 : [ 69018 ; 93177 [ = 1 3 : [ 93177 ; 117335 [ = 0 4 : [ 117335 ; 141494 [ = 0 5 : [ 141494 ; 165653 ] = 8 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.41 s) Code_Saturne: ../../../source/src/mesh/cs_mesh_location.c:595: Warning The group "VLM_SH" in the selection criteria: "VLM_SH" does not correspond to any boundary face. Code_Saturne: ../../../source/src/mesh/cs_mesh_location.c:595: Warning The group "VLM_VPR" in the selection criteria: "VLM_VPR" does not correspond to any boundary face. Code_Saturne: ../../../source/src/mesh/cs_mesh_location.c:595: Warning The group "VLM_ECO" in the selection criteria: "VLM_ECO" does not correspond to any boundary face. Code_Saturne: ../../../source/src/mesh/cs_mesh_location.c:595: Warning The group "VLM_GVH" in the selection criteria: "VLM_GVH" does not correspond to any boundary face. --- 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) Volume zone "all_cells" id = 1 Number of cells = 7306754 Volume = 1943.8 Center of gravity = (4.0629, 19.712, 1.4533e-05) Surface = -1 (not computed) Volume zone "SH" id = 2 Number of cells = 0 Volume = 0 Center of gravity = (0, 0, 0) Surface = -1 (not computed) Volume zone "VPR" id = 3 Number of cells = 0 Volume = 0 Center of gravity = (0, 0, 0) Surface = -1 (not computed) Volume zone "ECO" id = 4 Number of cells = 0 Volume = 0 Center of gravity = (0, 0, 0) Surface = -1 (not computed) Volume zone "GVH" id = 5 Number of cells = 0 Volume = 0 Center of gravity = (0, 0, 0) 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) Boundary zone "Outlet" id = 1 Number of faces = 56 Surface = 10.462 Center of gravity = (11.75, 60.2, -0.00013365) Perimeter = -1 (not computed) Boundary zone "Inlet" id = 2 Number of faces = 241 Surface = 7.2824 Center of gravity = (-16.867, 2.5, -5.0539e-06) Perimeter = -1 (not computed) Boundary zone "Wall" id = 3 Number of faces = 460965 Surface = 2126.4 Center of gravity = (3.0667, 20.608, 0.00080163) Perimeter = -1 (not computed) ALMAX = 0.12480E+02 (Characteristic length ) ALMAX is the cubic root of the domain volume. ALMAX is the length used to initialize the turbulence. -------------------------------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.6160E+00 0.6160E+00 molecular_viscos 0.2807E-04 0.2807E-04 turbulent_viscos 0.9416E+01 0.9416E+01 specific_heat 0.1100E+04 0.1100E+04 ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- TempC 1 0.4329E-01 0.4329E-01 ----------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 241 Smooth wall 5 460965 Rough wall 6 0 Symmetry 4 0 Free outlet 3 56 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- ** Clippings for computed fields ----------------------------- field initial min initial max clips to min clips to max - ------- -------------- -------------- ------------ ------------ a k a epsilon a TempC ** 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 k 1.5 1.5 1.5 1.5 v epsilon 0.013248 0.013248 0.013248 0.013248 v TempC 300 300 300 300 v Density 0.61604 0.61604 0.61604 0.61604 v LamVisc 2.8068e-05 2.8068e-05 2.8068e-05 2.8068e-05 v TurbVisc 9.416 9.416 9.416 9.416 v CourantNb 0 0 0 0 v FourierNb 0 0 0 0 v total_pressure 1.0064e+05 1.0132e+05 1.0102e+05 1.011e+05 v SpecHeat 1099.8 1099.8 1099.8 1099.8 v Local Time Step 0.0001 0.0001 0.0001 0.0001 v ThermalCond 0.043292 0.043292 0.043292 0.043292 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 300 300 300 300 v Boundary Density 0.61604 0.61604 0.61604 0.61604 v Yplus 0 0 0 0 =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.100000000E-03 TIME STEP NUMBER 1 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 241 0.000000000E+00 Smooth wall 5 460965 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 56 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.52285E-12 0.85707E-12 Friction velocity uet : 0.40019E-13 0.13360E-12 Friction velocity uk : 0.45936E+00 0.67082E+00 Dimensionless distance yplus : 0.26109E+01 0.40682E+04 Friction thermal sca. tstar : 0.00000E+00 0.00000E+00 Dim-less thermal sca. tplus : 0.18616E+01 0.23379E+02 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 0 Total number of wall faces : 460965 ------------------------------------------------------------ ** INFORMATION ON CONVERGENCE -------------------------- ---------------------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual Drift Time residual ---------------------------------------------------------------------------- c Velocity 0.47926E+04 14 0.10377E+00 0.33660E+07 0.10172E+05 c Velocity[X] 0.97332E+06 c Velocity[Y] 0.23734E+07 c Velocity[Z] 0.19235E+05 c Pressure 0.10185E+02 54 0.79045E+00 0.43678E+00 0.10789E+05 c k 0.54804E+05 5 0.10878E-02 0.15899E+02 0.26395E+03 c epsilon 0.23214E+06 103 0.19228E+00 0.16610E+09 0.27423E+06 c TempC 0.11883E+11 4 0.14702E-03 0.41444E+05 0.64685E+02 ---------------------------------------------------------------------------- ** Clippings for computed fields ----------------------------- field initial min initial max clips to min clips to max - --------------- -------------- -------------- ------------ ------------ a k a epsilon -0.90651 1939 a TempC a Local Time Step 2.97e-08 512827 ** Field values on cells --------------------- field minimum maximum set mean spatial mean - --------------- -------------- -------------- -------------- -------------- v Velocity[X] -29.858 82.692 0.11506 4.6102 v Velocity[Y] -38.543 116.44 11.625 9.3317 v Velocity[Z] -34.384 37.672 0.0021319 -0.0017468 v ǁVelocityǁ 1.4962e-07 116.47 12.126 12.695 v Pressure 5594.8 1.2845e+07 5.4397e+06 5.8604e+06 v k 0.8909 6.005 1.4967 1.5009 v epsilon 4.9032e-06 38243 433.56 3.7267 v TempC 300 430.97 300 300.07 v Density 0.61604 0.61604 0.61604 0.61604 v LamVisc 2.8068e-05 2.8068e-05 2.8068e-05 2.8068e-05 v TurbVisc 9.416 9.416 9.416 9.416 v CourantNb 6.1259e-11 0.69797 0.038289 0.01528 v FourierNb 0.062976 336.69 10.568 1.159 v total_pressure 1.0623e+05 1.2947e+07 5.5408e+06 5.9615e+06 v SpecHeat 1099.8 1099.8 1099.8 1099.8 v Local Time Step 1e-06 0.00011 4.4044e-05 0.00010819 v ThermalCond 0.043292 0.043292 0.043292 0.043292 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 300 552 300.16 302.09 v Boundary Density 0.61604 0.61604 0.61604 0.61604 v Yplus 0 4068.2 138.4 1247 ** Computed values on cells ------------------------ criterion minimum maximum set mean spatial mean --------------- -------------- -------------- -------------- -------------- Courant/Fourier 0.067093 336.69 10.606 1.1743 INSTANT 0.200000000E-03 TIME STEP NUMBER 2 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 241 -0.195708011E+03 Smooth wall 5 460965 -0.162438938E-16 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 56 0.195708016E+03 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.50992E+02 0.77013E+02 Friction velocity uet : 0.24984E-08 0.75738E+01 Friction velocity uk : 0.32220E+00 0.12210E+01 Dimensionless distance yplus : 0.18354E+01 0.40112E+04 Friction thermal sca. tstar : 0.00000E+00 0.00000E+00 Dim-less thermal sca. tplus : 0.13087E+01 0.23345E+02 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 84248 Nb of faces within the viscous sub-layer : 22 Total number of wall faces : 460965 ------------------------------------------------------------ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING: MESH TOO REFINED AT THE WALL @ ======== @ The mesh seems to be too refined at the wall to use @ a high-Reynolds turbulence model. @ @ The last time step at which too small values for the @ dimensionless distance to the wall (yplus) have been @ observed is the time step 2 @ @ The minimum value for yplus must be greater than the @ limit value YPLULI = 0.23810E+01 @ @ Have a look at the distribution of yplus at the wall @ (with EnSight or ParaView for example) to conclude on @ the way the results quality might be affected. @ @ This warning is only printed at the first two @ occurences of the problem and at the last time step @ of the calculation. The vanishing of the message does @ not necessarily mean the vanishing of the problem. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ ** INFORMATION ON CONVERGENCE -------------------------- ---------------------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual Drift Time residual ---------------------------------------------------------------------------- c Velocity 0.61411E+06 16 0.57991E-02 0.41563E+05 0.12095E+04 c Velocity[X] 0.11487E+05 c Velocity[Y] 0.29167E+05 c Velocity[Z] 0.90867E+03 c Pressure 0.12645E+02 48 0.10023E+00 0.66654E+00 0.98556E+04 c k 0.49323E+05 7 0.20008E-04 0.11152E+01 0.66438E+02 c epsilon 0.10499E+06 16 0.68130E-02 0.18477E+06 0.53652E+04 c TempC 0.10877E+11 8 0.45106E-04 0.11606E+05 0.33641E+02 ---------------------------------------------------------------------------- ** Clippings for computed fields ----------------------------- field initial min initial max clips to min clips to max - --------------- -------------- -------------- ------------ ------------ a k a epsilon -1.5194 1990 a TempC a Local Time Step ** Field values on cells --------------------- field minimum maximum set mean spatial mean - --------------- -------------- -------------- -------------- -------------- v Velocity[X] -25.835 73.263 0.10309 4.1048 v Velocity[Y] -29.749 96.291 10.424 8.287 v Velocity[Z] -32.198 32.998 0.0020299 -0.0009157 v ǁVelocityǁ 2.3883e-07 96.337 10.874 11.28 v Pressure -1.2853e+06 153.59 -5.457e+05 -5.8815e+05 v k 0.8591 7.2485 1.4969 1.5014 v epsilon 7.9329e-06 35314 427.52 3.6757 v TempC 299.99 513.3 300.01 300.1 v Density 0.50142 0.61604 0.61603 0.61597 v LamVisc 2.8068e-05 3.2376e-05 2.8069e-05 2.8071e-05 v TurbVisc 1.2481e-06 25442 2.2024 7.7969 v CourantNb 1.9331e-10 0.68135 0.03705 0.014729 v FourierNb 2.9906e-05 10 1.9882 0.69109 v total_pressure -1.184e+06 1.0079e+05 -4.4467e+05 -4.8706e+05 v SpecHeat 1099.8 1134.9 1099.8 1099.8 v Local Time Step 1.1e-06 0.000121 4.7774e-05 0.00011889 v ThermalCond 0.043292 0.051653 0.043292 0.043296 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 300 552 300.18 302.16 v Boundary Density 0.50142 0.61604 0.61597 0.61532 v Yplus 0 4011.2 138.31 1246.8 ** Computed values on cells ------------------------ criterion minimum maximum set mean spatial mean --------------- -------------- -------------- -------------- -------------- Courant/Fourier 3.8828e-05 10.447 2.0252 0.70582 Options set or changed by "control_file": ------------------------- max_time_step 3 (current: 3) Finished reading "control_file". INSTANT 0.300000000E-03 TIME STEP NUMBER 3 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 241 -0.171714304E+03 Smooth wall 5 460965 0.286583053E-14 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 56 0.171714594E+03 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.33087E+02 0.66046E+02 Friction velocity uet : 0.11188E-07 0.59632E+01 Friction velocity uk : 0.31524E+00 0.13451E+01 Dimensionless distance yplus : 0.17958E+01 0.40074E+04 Friction thermal sca. tstar : 0.00000E+00 0.00000E+00 Dim-less thermal sca. tplus : 0.12805E+01 0.23343E+02 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 83312 Nb of faces within the viscous sub-layer : 23 Total number of wall faces : 460965 ------------------------------------------------------------ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING: MESH TOO REFINED AT THE WALL @ ======== @ The mesh seems to be too refined at the wall to use @ a high-Reynolds turbulence model. @ @ The last time step at which too small values for the @ dimensionless distance to the wall (yplus) have been @ observed is the time step 3 @ @ The minimum value for yplus must be greater than the @ limit value YPLULI = 0.23810E+01 @ @ Have a look at the distribution of yplus at the wall @ (with EnSight or ParaView for example) to conclude on @ the way the results quality might be affected. @ @ This warning is only printed at the first two @ occurences of the problem and at the last time step @ of the calculation. The vanishing of the message does @ not necessarily mean the vanishing of the problem. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ Write final restart files checkpoint at iteration 3, Physical time 0.30000E-03 ** Writing the main restart file ----------------------------- Start writing Wrote field names and types to checkpoint at iteration 3: 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.43248E+06 24 0.14710E-01 0.21264E+04 0.27114E+03 c Velocity[X] 0.61435E+03 c Velocity[Y] 0.14199E+04 c Velocity[Z] 0.92210E+02 c Pressure 0.97231E+01 50 0.13844E+00 0.82406E+00 0.35607E+05 c k 0.45015E+05 5 0.16762E-04 0.97546E+00 0.59074E+02 c epsilon 0.92632E+05 17 0.39144E-02 0.13075E+06 0.51811E+04 c TempC 0.99233E+10 6 0.14313E-03 0.57517E+04 0.21977E+02 ---------------------------------------------------------------------------- ** Clippings for computed fields ----------------------------- field initial min initial max clips to min clips to max - --------------- -------------- -------------- ------------ ------------ a k a epsilon -0.54223 1491 a TempC a Local Time Step ** Field values on cells --------------------- field minimum maximum set mean spatial mean - --------------- -------------- -------------- -------------- -------------- v Velocity[X] -24.826 71.003 0.10025 3.9886 v Velocity[Y] -28.474 88.838 10.18 8.0512 v Velocity[Z] -28.413 30.245 0.0021487 0.00011055 v ǁVelocityǁ 4.3349e-08 88.891 10.618 10.959 v Pressure -2.626e+05 1659.4 -1.1006e+05 -1.1843e+05 v k 0.7319 8.3338 1.4964 1.5019 v epsilon 6.3596e-06 32461 421.08 3.6187 v TempC 299.95 527.54 300.01 300.12 v Density 0.44892 0.61604 0.61603 0.61595 v LamVisc 2.8068e-05 3.4836e-05 2.8069e-05 2.8072e-05 v TurbVisc 1.2411e-06 15737 2.1414 7.866 v CourantNb 3.4495e-11 0.73748 0.039717 0.015748 v FourierNb 3.2358e-05 10 2.0146 0.73153 v total_pressure -1.6131e+05 1.023e+05 -9036.1 -17328 v SpecHeat 1099.8 1157.4 1099.8 1099.8 v Local Time Step 1.21e-06 0.0001331 5.1793e-05 0.00013062 v ThermalCond 0.043292 0.056908 0.043292 0.043298 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 300 552 300.18 302.17 v Boundary Density 0.44892 0.61604 0.61596 0.61514 v Yplus 0 4007.4 138.08 1247 ** Computed values on cells ------------------------ criterion minimum maximum set mean spatial mean --------------- -------------- -------------- -------------- -------------- Courant/Fourier 4.2632e-05 10.432 2.0544 0.74728 =============================================================== FINAL STAGE OF THE CALCULATION ============================== =========================================================== =========================================================== END OF CALCULATION ================== =============================================================== Destroying structures and ending computation