command: ./cs_solver --mpi *************************************************************** (R) Code_Saturne Version 7.2-alpha Copyright (C) 1998-2022 EDF S.A., France build Mon 03 Jan 2022 20:47:52 GMT 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: Mon 10 Jan 2022 23:35:45 GMT System: Linux 3.10.0-693.5.2.el7.x86_64 Machine: node02.swmgmt.eureka Processor: model name : Intel(R) Xeon(R) CPU E5-2660 v4 @ 2.00GHz Memory: 128658 MB User: rw00793 (Wang) Directory: /users/rw00793/cases/Bohn_closed_Ra1E7_2D_refine/les_restart_2/RESU/20220110-2335 MPI ranks: 80 (appnum attribute: 0) MPI ranks per node: 24 - 28 OpenMP threads: 1 Processors/node: 1 Compilers used for build: C compiler: gcc (GCC) 4.9.3 C++ compiler: g++ (GCC) 4.9.3 Fortran compiler: GNU Fortran (GCC) 4.8.5 20150623 (Red Hat 4.8.5-16) MPI version: 3.0 (MPICH 3.1.2) OpenMP version: 4.0 External libraries: ParMETIS 4.0.3 PT-SCOTCH 6.1.1 I/O read method: collective MPI-IO (explicit offsets) I/O write method: collective MPI-IO (explicit offsets) I/O rank step: 1 =============================================================== 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. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ 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 TempK 1 0.1000E+01 Wall distance 0 0.1000E+01 Wall Y+ 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 ) ITPCOL = 0 (0: staggered; 1: collocated ) ----------------------------------------------------------- ** 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 ) 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 ) ----------------------------------------------------------- ** WALL DISTANCE COMPUTATION ------------------------- ICDPAR = -1 ( 1: std, reread if restart (-1: std, recomputed if restrt ( 2: old, reread if restart (-2: old, recomputed if restrt ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) ----------------------------------------------------- Variable Number ITURT ----------------------------------------------------- TempK 1 0 ------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ TempK 1 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- TempK 1 -1 0.0000E+00 0.1000E+13 ------------------------------------------------------- ------------------------------------------------------------- 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.5000E+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 = 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 = 372000 (Initial time step) NTMABS = 432000 (Final time step required) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 1000 (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 = 7 (Force exerted by the fluid on the boundary) ipstyp = 1 (y+ at boundary) ipsttp = 1 (T+ at boundary) ipstft = 1 (Thermal flux at boundary) ipstnu = 1 (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: every 1200 time steps, at start, 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: restart/mesh_input.csm Finished reading: restart/mesh_input.csm No "partition_input/domain_number_80" file available; ---------------------------------------------------------- Partitioning 178800 cells to 80 domains on 80 ranks (SCOTCH_dgraphPart). wall-clock time: 0.300346 s Number of cells per domain (histogramm): [ 2124 ; 2143 [ = 1 [ 2143 ; 2163 [ = 4 [ 2163 ; 2183 [ = 8 [ 2183 ; 2203 [ = 2 [ 2203 ; 2223 [ = 22 [ 2223 ; 2242 [ = 7 [ 2242 ; 2262 [ = 18 [ 2262 ; 2282 [ = 11 [ 2282 ; 2302 [ = 3 [ 2302 ; 2322 ] = 4 Writing file: partition_output/domain_number_80 Finished writing: partition_output/domain_number_80 Partitioning finished (0.336 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 = 2124 maximum value = 2322 1 : [ 2124 ; 2163 [ = 5 2 : [ 2163 ; 2203 [ = 10 3 : [ 2203 ; 2242 [ = 29 4 : [ 2242 ; 2282 [ = 29 5 : [ 2282 ; 2322 ] = 7 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 2278 maximum value = 2516 1 : [ 2278 ; 2325 [ = 5 2 : [ 2325 ; 2373 [ = 17 3 : [ 2373 ; 2420 [ = 28 4 : [ 2420 ; 2468 [ = 22 5 : [ 2468 ; 2516 ] = 8 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 132 maximum value = 210 1 : [ 132 ; 147 [ = 19 2 : [ 147 ; 163 [ = 26 3 : [ 163 ; 178 [ = 2 4 : [ 178 ; 194 [ = 18 5 : [ 194 ; 210 ] = 15 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 4307 maximum value = 4735 1 : [ 4307 ; 4392 [ = 4 2 : [ 4392 ; 4478 [ = 17 3 : [ 4478 ; 4563 [ = 25 4 : [ 4563 ; 4649 [ = 26 5 : [ 4649 ; 4735 ] = 8 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 4284 maximum value = 4703 1 : [ 4284 ; 4367 [ = 5 2 : [ 4367 ; 4451 [ = 17 3 : [ 4451 ; 4535 [ = 31 4 : [ 4535 ; 4619 [ = 21 5 : [ 4619 ; 4703 ] = 6 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 4 1 : [ 3 ; 4 ] = 178800 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 3 maximum value = 9 1 : [ 3 ; 4 [ = 41 2 : [ 4 ; 5 [ = 10 3 : [ 5 ; 6 [ = 17 4 : [ 6 ; 7 [ = 7 5 : [ 7 ; 9 ] = 5 ---------------------------------------------------------- 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 = 2125 maximum value = 2455 1 : [ 2125 ; 2191 [ = 2 2 : [ 2191 ; 2257 [ = 8 3 : [ 2257 ; 2323 [ = 33 4 : [ 2323 ; 2389 [ = 27 5 : [ 2389 ; 2455 ] = 10 Histogram of volume mesh matrix profile/lines per rank: minimum value = 97 maximum value = 262 1 : [ 97 ; 130 [ = 19 2 : [ 130 ; 163 [ = 54 3 : [ 163 ; 196 [ = 4 4 : [ 196 ; 229 [ = 1 5 : [ 229 ; 262 ] = 2 Mesh coordinates: minimum and maximum X : 0.0000000e+00 5.0000002e-04 Y : -2.3999366e-01 2.3999366e-01 Z : -2.3999366e-01 2.3999366e-01 Mesh Number of cells: 178800 Number of interior faces: 356408 Number of boundary faces: 359984 Number of vertices: 359984 Groups: "Bohn_1995_configB_2D_Ra1E7_refin" cells: 178800 boundary faces: 359984 "DOWN" boundary faces: 1192 "LEFT" boundary faces: 178800 "RIGHT" boundary faces: 178800 "UP" boundary faces: 1192 --- Information on the volumes Minimum control volume = 6.5205805e-11 Maximum control volume = 5.2499783e-10 Total volume for the domain = 6.5907197e-05 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.0045 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 178800 Volume = 6.5907e-05 Center of gravity = (0.00025, 8.2883e-09, -8.2797e-09) Surface = -1 (not computed) Volume zone "all_cells" id = 1 Number of cells = 178800 Volume = 6.5907e-05 Center of gravity = (0.00025, 8.2883e-09, -8.2797e-09) Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 359984 Surface = 0.26478 Center of gravity = (0.00025, 8.7979e-09, -8.7894e-09) Perimeter = -1 (not computed) Boundary zone "RIGHT" id = 1 Number of faces = 178800 Surface = 0.13181 Center of gravity = (0.0005, 8.751e-09, -8.7424e-09) Perimeter = -1 (not computed) Boundary zone "DOWN" id = 2 Number of faces = 1192 Surface = 0.00039265 Center of gravity = (0.00025, 1.95e-09, -1.9538e-09) Perimeter = -1 (not computed) Boundary zone "UP" id = 3 Number of faces = 1192 Surface = 0.00075389 Center of gravity = (0.00025, 2.8751e-08, -2.8781e-08) Perimeter = -1 (not computed) Boundary zone "LEFT" id = 4 Number of faces = 178800 Surface = 0.13181 Center of gravity = (7.0183e-18, 8.751e-09, -8.7424e-09) Perimeter = -1 (not computed) Probe set: "probes": maximum distance between cell centers and requested coordinates: 3.259e-07 ALMAX = 0.40393E-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 = 372000 Reading the previous time step number (restarting computation) TTPABS = 0.3582E+03 Reading options complete Read variables from restart: restart/main.csc Reading complete READING THE AUXILIARY RESTART FILE Start reading Finished reading options restart/auxiliary.csc: 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.csc: section "time_moments:sizes" not present. End reading ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1021E+01 Infinity molecular_viscos 0.1810E-04 0.1810E-04 turbulent_viscos 0.1793E-09 Infinity ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- TempK 1 0.2590E-01 0.2590E-01 ----------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 0 Smooth wall 5 2384 Rough wall 6 0 Symmetry 4 357600 Free outlet 3 0 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 TempK ** Field values on cells --------------------- field minimum maximum set mean spatial mean - ---------------- -------------- -------------- -------------- -------------- v Velocity[X] -1.3814e-11 1.7239e-11 -5.1037e-16 -3.9398e-16 v Velocity[Y] -0.34354 0.31963 0.00012551 -0.00085584 v Velocity[Z] -0.33995 0.27782 -0.010218 -0.011866 v ǁVelocityǁ 0 0.35436 0.091759 0.095616 v Pressure -0.16049 0.17511 0.0076409 0.0011287 v TempK 0 345.84 163.23 164.36 v Density 1.0208 inf inf inf v TurbVisc 1.7931e-10 inf -nan -nan v CourantNb 0 0 0 0 v FourierNb 0 0 0 0 v total_pressure 1.0132e+05 1.0133e+05 1.0133e+05 1.0132e+05 v Wall distance 0 0.057692 0.013109 0.014458 v Wall Y+ 0 207.28 23.596 25.743 v volume_forces[X] 0 0 0 0 v volume_forces[Y] -7.2807 7.7347 -0.0014006 0.0054735 v volume_forces[Z] -7.8166 7.2972 0.005871 -0.0017102 v ǁvolume_forcesǁ 0.00034577 8.0707 3.0387 3.1546 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 0 345.84 163.21 164.36 v Boundary Density 1.0208 inf inf inf v Yplus 0 0 0 0 /users/rw00793/code_saturne_72alpha/code_saturne-master/src/base/cs_log_iteration.c:842: Fatal error. Invalid (not-a-number) values detected for a field. Call stack: 1: 0x2af234c3006a <+0x16006a> (libsaturne-7.2.so) 2: 0x2af234c32228 (libsaturne-7.2.so) 3: 0x2af234b8ec5d (libsaturne-7.2.so) 4: 0x2af2348d3fe9 (libcs_solver-7.2.so) 5: 0x2af237cc7c05 <__libc_start_main+0xf5> (libc.so.6) 6: 0x401209 <> (cs_solver) End of stack