command: ./cs_solver *************************************************************** (R) Code_Saturne Version 7.2.0 Copyright (C) 1998-2022 EDF S.A., France revision 4a1275a4d build sam. 26 nov. 2022 09:00:41 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: lun. 05 déc. 2022 12:35:15 System: Linux 5.19.0-26-generic (Ubuntu 22.10) Machine: julien-virtual-machine Processor: model name : Intel(R) Xeon(R) CPU E5-2678 v3 @ 2.50GHz Memory: 62766 MB User: julien (julien) Directory: /home/julien/Travail/flow_past_cyl_1/RESU/20221205-1235 OpenMP threads: 1 Processors/node: 24 Compilers used for build: C compiler: gcc (Ubuntu 12.2.0-3ubuntu1) 12.2.0 C++ compiler: g++ (Ubuntu 12.2.0-3ubuntu1) 12.2.0 Fortran compiler: GNU Fortran (Ubuntu 12.2.0-3ubuntu1) 12.2.0 MPI version: 3.1 (Open MPI 4.1.4) OpenMP version: 4.5 I/O read method: standard input and output, serial access I/O write method: standard input and output, serial access =============================================================== CALCULATION PREPARATION ======================= =========================================================== No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 8 (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 omega 1 0.1000E+01 TempC 1 0.1000E+01 Wall distance 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.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 ) ----------------------------------------------------------- ** 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 ----------------------------------------------------- TempC 1 0 ------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ TempC 1 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- TempC 1 -1 0.2000E+02 0.1500E+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 = 500 (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 = 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: at end -5: name: directory: monitoring format: time_plot options: DAT 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 Before insertion of user-defined boundary Number of cells: 18956 Number of interior faces: 29253 Number of boundary faces: 38282 Number of vertices: 21336 After insertion of user-defined boundary Number of cells: 18956 Number of interior faces: 0 Number of boundary faces: 96788 Number of vertices: 117752 Before insertion of user-defined boundary Number of cells: 18956 Number of interior faces: 0 Number of boundary faces: 96788 Number of vertices: 117752 After insertion of user-defined boundary Number of cells: 18956 Number of interior faces: 0 Number of boundary faces: 96788 Number of vertices: 117752 Number of cells: 18956 ---------------------------------------------------------- Number of cells + halo cells: 18956 ---------------------------------------------------------- Number of interior faces: 0 ---------------------------------------------------------- Number of boundary faces: 96788 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 0 maximum value = 0 1 : [ 0 ; 0 ] = 18956 ---------------------------------------------------------- Extrusion: 85 boundary faces selected. 8500 cells added. Extrusion: 85 boundary faces selected. 17000 cells added. 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 number of halo-adjacent elements: 44456 Numbering for interior faces: type: default number of halo-adjacent elements: 25500 Numbering for boundary faces: type: default number of halo-adjacent elements: 198788 ---------------------------------------------------------- Matrix bandwidth for volume mesh : 42395 Matrix profile/lines for volume mesh : 203 Mesh coordinates: minimum and maximum X : -1.5000000e+00 2.5000000e+00 Y : -5.0000000e-03 2.6020852e-18 Z : -3.7500000e-01 3.7500000e-01 Mesh Number of cells: 44456 Number of interior faces: 25500 Number of boundary faces: 198788 Number of vertices: 219752 Groups: "Wall" boundary faces: 800 "auto:internal_coupling_0" "auto:internal_coupling_0_fluid" "auto:internal_coupling_0_solid" "destination_fluid" boundary faces: 43595 "destination_solid" boundary faces: 861 "fluid" cells: 43595 "inlet" interior faces: 85 boundary faces: 85 "outlet" interior faces: 85 boundary faces: 85 "solid" cells: 861 "source_fluid" boundary faces: 43595 "source_solid" boundary faces: 861 --- Information on the volumes Minimum control volume = 2.8499338e-09 Maximum control volume = 4.9340932e-07 Total volume for the domain = 1.5000000e-02 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.0653 s) Code_Saturne: /home/julien/Code_Saturne/code_saturne-7.2.0/src/mesh/cs_mesh_location.c:595: Warning The group "42" in the selection criteria: "42" does not correspond to any boundary face. --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 44456 Volume = 0.015 Center of gravity = (0.5, -0.0025, 1.409e-09) Surface = -1 (not computed) Volume zone "all_cells" id = 1 Number of cells = 44456 Volume = 0.015 Center of gravity = (0.5, -0.0025, 1.409e-09) Surface = -1 (not computed) Volume zone "fluid" id = 2 Number of cells = 43595 Volume = 0.014961 Center of gravity = (0.50131, -0.0025, 3.9679e-10) Surface = -1 (not computed) Volume zone "solid" id = 3 Number of cells = 861 Volume = 3.9205e-05 Center of gravity = (1.6892e-07, -0.0025, 3.8769e-07) Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 198788 Surface = 11.244 Center of gravity = (0.41829, -0.0025, 1.5583e-05) Perimeter = -1 (not computed) Boundary zone "wall" id = 1 Number of faces = 800 Surface = 0.04 Center of gravity = (0.5, -0.0025, -2.3539e-16) Perimeter = -1 (not computed) Boundary zone "sym1f" id = 2 Number of faces = 43595 Surface = 2.9922 Center of gravity = (0.50131, -3.4082e-18, -1.1555e-16) Perimeter = -1 (not computed) Boundary zone "sym2f" id = 3 Number of faces = 861 Surface = 0.007841 Center of gravity = (-4.7262e-19, 1.6004e-18, 5.241e-19) Perimeter = -1 (not computed) Boundary zone "inlet" id = 4 Number of faces = 85 Surface = 0.00375 Center of gravity = (-1.5, -0.0025, 1.0729e-17) Perimeter = -1 (not computed) Boundary zone "outlet" id = 5 Number of faces = 85 Surface = 0.00375 Center of gravity = (2.5, -0.0025, 2.7274e-17) Perimeter = -1 (not computed) Boundary zone "sym1s" id = 6 Number of faces = 43595 Surface = 2.9922 Center of gravity = (0.50131, -0.005, -1.1554e-16) Perimeter = -1 (not computed) Boundary zone "sym2s" id = 7 Number of faces = 861 Surface = 0.007841 Center of gravity = (-6.0596e-19, -0.005, 3.2492e-19) Perimeter = -1 (not computed) Boundary zone "interface" id = 8 Number of faces = 0 Surface = 0 Center of gravity = (0, 0, 0) Perimeter = -1 (not computed) Internal coupling coupling_id = 0 Face group selection criterion: auto:internal_coupling_0 Assign interior faces group name: auto:internal_coupling_0_fluid Assign interior faces group name: auto:internal_coupling_0_solid Volume zones: solid Locator: n dist points (total coupled boundary faces) = 0 Probe set: "probes": maximum distance between cell centers and requested coordinates: 1.456e-05 ALMAX = 0.24662E+00 (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.8500E+00 0.2700E+04 molecular_viscos 0.2500E-04 0.2500E-04 turbulent_viscos 0.0000E+00 0.1027E-01 specific_heat 0.9000E+03 0.1017E+04 ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- TempC 1 0.3350E-01 0.2000E+03 ----------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 85 Smooth wall 5 109706 Rough wall 6 0 Symmetry 4 88912 Free outlet 3 85 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 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 0.0024 0.0024 0.0024 0.0024 v omega 0.19864 0.19864 0.19864 0.19864 v TempC 150 150 150 150 v TurbVisc 0 0.01027 0.010071 0.010243 v CourantNb 0 0 0 0 v FourierNb 0 0 0 0 v total_pressure 1.0132e+05 1.0132e+05 1.0132e+05 1.0132e+05 v Local Time Step 0.0001 0.0001 0.0001 0.0001 v Wall distance 0 0 0 0 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 Boundary temperature 150 150 150 150 v Boundary Density 0.85 2700 59.303 18.894 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 85 0.000000000E+00 Smooth wall 5 109706 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 88912 0.000000000E+00 Free outlet 3 85 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.30200E-12 0.35113E-08 Friction velocity uk : 0.21007E-08 0.19578E-01 Dimensionless distance yplus : 0.28479E-03 0.33415E+01 Friction thermal sca. tstar : 0.00000E+00 0.00000E+00 Dim-less thermal sca. tplus : 0.28479E-03 0.33415E+01 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 109706 Total number of wall faces : 109706 ------------------------------------------------------------ Flexible Conjugate Gradient [Wall distance]: divergence after 1 iterations: initial residual: inf; current residual: -nan Flexible Conjugate Gradient [Wall distance]: divergence fallback from Multigrid to Jacobi (diagonal) preconditioning for re-try and subsequent solves. Flexible Conjugate Gradient [Wall distance]: divergence after 1 iterations: initial residual: inf; current residual: -nan /home/julien/Code_Saturne/code_saturne-7.2.0/src/alge/cs_sles_it.c:4994: Fatal error. Flexible Conjugate Gradient: error (divergence) solving for Wall distance Call stack: 1: 0x7f62ccf7d06a (libsaturne-7.2.so) 2: 0x7f62ccf726da (libsaturne-7.2.so) 3: 0x7f62ccf73867 (libsaturne-7.2.so) 4: 0x7f62ccb431c2 (libsaturne-7.2.so) 5: 0x7f62ccc0b4d7 (libsaturne-7.2.so) 6: 0x7f62ccc67d44 (libsaturne-7.2.so) 7: 0x7f62ccac4622 (libsaturne-7.2.so) 8: 0x7f62cdaddbea (libcs_solver-7.2.so) 9: 0x7f62cc623510 <+0x23510> (libc.so.6) 10: 0x7f62cc6235c9 <__libc_start_main+0x89> (libc.so.6) 11: 0x55ecc53580a5 <_start+0x25> (cs_solver) End of stack