command: ./cs_solver *************************************************************** (R) Code_Saturne Version 6.0.2 Copyright (C) 1998-2019 EDF S.A., France revision 6.0.2-246d8963f build mar 19 mag 2020 16:26:05 CEST MPI version 3.1 (Open MPI 2.1.1) 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: mar 26 mag 2020 16:37:13 CEST System: Linux 5.3.0-53-generic (Ubuntu 18.04.4 LTS) Machine: pici-Linux64 Processor: model name : AMD Ryzen Threadripper 1900X 8-Core Processor Memory: 64237 MB User: pici (Piero Cirrottola) Directory: /D3_Ubuntu/Code_Saturne_6_test/full_domain/case1/RESU/sub_1 OpenMP threads: 4 Processors/node: 8 I/O read method: standard input and output, serial access I/O write method: standard input and output, serial access External libraries for partitioning: ParMETIS 4.0.3 PT-SCOTCH 6.0.6 Reading metadata from file: "mesh_input/mesh_01" Reading metadata from file: "mesh_input/mesh_02" Reading metadata from file: "mesh_input/mesh_03" =============================================================== 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 = 2 (Nb scalars ) NSCAUS = 1 (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 scalar1 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 ) THETFL = -0.99900E+03 (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.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 ISCACP ITURT VISLS0 SIGMAS ------------------------------------------------------------- scalar1 1 0 0 0.6205E-01 0.1000E+01 TempC 2 1 0 0.2495E-01 0.1000E+01 ------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ scalar1 1 0.8000E+00 TempC 2 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- scalar1 1 -1 0.0000E+00 0.4000E+03 TempC 2 -1 0.0000E+00 0.4000E+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. ISCACP = 0 or 1 2 (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 2 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 = 1000 (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 = 20 (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: MED 4.0.0 options: separate_meshes time dependency: fixed mesh output: every 20 time steps and 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/mesh_01 Finished reading: mesh_input/mesh_01 Reading file: mesh_input/mesh_02 Finished reading: mesh_input/mesh_02 Reading file: mesh_input/mesh_03 Finished reading: mesh_input/mesh_03 ------------------------------------------------------- Joining number 1: Selection criteria: "all[]" 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: 1650 Number of interior faces: 3135 Number of boundary faces: 3630 Number of vertices: 3636 Global number of boundary faces selected for joining: 3630 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 482 Vertex-Vertex intersections: 166 Other intersections: 316 Joining operation is non-conforming. Global number of new vertices to create: 316 Equivalences broken for 0 edges. Edge intersections and vertex creation done. Pre-merge for 288 global element couples. Max size of a merge set of vertices: 3 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: 3625 Global configuration after the joining operation: Global number of border faces to add: 3585 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: 1650 Number of interior faces: 3175 Number of boundary faces: 3585 Number of vertices: 3622 Joining 1 completed (0.0542 s). All joining operations successfully finished: Wall clock time: 0.0581 Number of cells: 1650 ---------------------------------------------------------- Number of cells + halo cells: 1650 ---------------------------------------------------------- Number of interior faces: 3175 ---------------------------------------------------------- Number of boundary faces: 3585 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 2 maximum value = 5 1 : [ 2 ; 3 [ = 8 2 : [ 3 ; 4 [ = 264 3 : [ 4 ; 5 ] = 1378 ---------------------------------------------------------- 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: 1650 Numbering for interior faces: type: threads number of threads: 4 number of exclusive groups: 2 number of elements in group 0: 3050 number of elements in group 1: 125 estimated thread imbalance: 0.158 Numbering for boundary faces: type: threads number of threads: 4 number of exclusive groups: 1 number of elements in group 0: 3585 estimated thread imbalance: 0.002 ---------------------------------------------------------- Matrix bandwidth for volume mesh : 680 Matrix profile/lines for volume mesh : 49 Mesh coordinates: minimum and maximum X : -5.0000000e-01 1.0000000e+00 Y : -1.0000000e+00 3.0000000e+00 Z : 0.0000000e+00 1.0000000e+00 Mesh Number of cells: 1650 Number of interior faces: 3175 Number of boundary faces: 3585 Number of vertices: 3622 Groups: "1" boundary faces: 10 "2" boundary faces: 10 "21" boundary faces: 20 "22" boundary faces: 40 "23" boundary faces: 20 "24" interior faces: 40 boundary faces: 20 "28" cells: 800 boundary faces: 800 "29" boundary faces: 800 "3" boundary faces: 10 "31" boundary faces: 10 "32" interior faces: 25 "33" boundary faces: 10 "34" boundary faces: 15 "38" cells: 150 boundary faces: 150 "39" boundary faces: 150 "4" boundary faces: 50 "5" interior faces: 15 "6" boundary faces: 60 "7" boundary faces: 10 "8" cells: 700 boundary faces: 700 "9" boundary faces: 700 --- Information on the volumes Minimum control volume = 4.9266504e-04 Maximum control volume = 4.9159297e-03 Total volume for the domain = 2.0500000e+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 % Computing geometric quantities (0.00179 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 1650 Volume = 2.0500000e+00 Surface = -1 (not computed) Volume zone "all_cells" id = 1 Number of cells = 1650 Volume = 2.0500000e+00 Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 3585 Surface = 1.7100000e+01 Perimeter = -1 (not computed) Boundary zone "BC_1" id = 1 Number of faces = 10 Surface = 5.0000000e-01 Perimeter = -1 (not computed) Boundary zone "BC_2" id = 2 Number of faces = 260 Surface = 1.2000000e+01 Perimeter = -1 (not computed) Boundary zone "BC_8" id = 3 Number of faces = 3300 Surface = 4.1000000e+00 Perimeter = -1 (not computed) Boundary zone "BC_19" id = 4 Number of faces = 15 Surface = 5.0000000e-01 Perimeter = -1 (not computed) Probe set: "probes": maximum distance between cell centers and requested coordinates: 2.514e-01 ALMAX = 0.12703E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. ALMAX is the length used to initialize the turbulence. ------------------------------------------------------------- ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Velocity 0.0000E+00 0.0000E+00 Velocity 0.0000E+00 0.0000E+00 Velocity 0.0000E+00 0.0000E+00 Pressure 0.0000E+00 0.0000E+00 k 0.6000E-03 0.6000E-03 epsilon 0.1041E-05 0.1041E-05 TempC 0.2000E+02 0.2000E+02 scalar1 0.1000E+02 0.1000E+02 --------------------------------- dt 0.1000E-01 0.1000E-01 --------------------------------- ------------------------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 10 Smooth wall 5 260 Rough wall 6 0 Symmetry 4 3300 Free outlet 3 15 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.100000000E-01 TIME STEP NUMBER 1 ============================================================= ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.7257E+03 0.7257E+03 molecular_viscos 0.8951E-04 0.8951E-04 turbulent_viscos 0.2258E+02 0.2258E+02 ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- scalar1 1 0.6205E-01 0.6205E-01 TempC 2 0.2495E-01 0.2495E-01 ----------------------------------------------- ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 10 0.000000000E+00 Smooth wall 5 260 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 3300 0.000000000E+00 Free outlet 3 15 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 --------------------------------------------------------------- Flexible Conjugate Gradient [Pressure]: divergence after 1 iterations: initial residual: -1.0000e+00; current residual: 2.2928e+02 Flexible Conjugate Gradient [Pressure]: divergence fallback from Multigrid to Jacobi (diagonal) preconditioning for re-try and subsequent solves. Flexible Conjugate Gradient [Pressure]: divergence after 1 iterations: initial residual: -1.0000e+00; current residual: 1.6923e+02 /D3_Ubuntu/Scaricati/CodeSaturne/code_saturne-6.0.2/src/alge/cs_sles_it.c:5722: Fatal error. Flexible Conjugate Gradient: error (divergence) solving for Pressure Call stack: 1: 0x7f8855be9c38 (libsaturne-6.0.so) 2: 0x7f8855bdd649 (libsaturne-6.0.so) 3: 0x7f8855bded61 (libsaturne-6.0.so) 4: 0x7f885600be43 <__cs_c_bindings_MOD_sles_solve_native+0x134> (libsaturne-6.0.so) 5: 0x7f8855a38c4a (libsaturne-6.0.so) 6: 0x7f8855a1ab05 (libsaturne-6.0.so) 7: 0x7f8855a48763 (libsaturne-6.0.so) 8: 0x7f88558d3303 (libsaturne-6.0.so) 9: 0x7f885701077f (libcs_solver-6.0.so) 10: 0x7f8857010048 (libcs_solver-6.0.so) 11: 0x7f8855450b97 <__libc_start_main+0xe7> (libc.so.6) 12: 0x55fe38ae389a <_start+0x2a> (cs_solver) End of stack