command: ./cs_solver --param setup.xml --mpi *************************************************************** (R) Code_Saturne Version 5.0.3 Copyright (C) 1998-2017 EDF S.A., France revision 5.0.3 build Mon 17 Dec 2018 03:25:56 PM CET 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: Fri 22 Nov 2019 11:53:36 AM CET System: Linux 3.10.0-514.el7.x86_64 Machine: node102 Processor: model name : Intel(R) Xeon(R) CPU X5660 @ 2.80GHz Memory: 96510 MB User: vivaldi-dan (Danielle VIVALDI - IRSN/PSN-RES/SEMIA/LIMAR) Directory: /TMPCALCUL/Etudes_CFD/vivaldi/saturne_FSI/mesh_ok/CASE1/RESU/20191122-1153 MPI ranks: 4 (appnum attribute: 0) OpenMP threads: 1 Processors/node: 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: METIS SCOTCH 6.0.4 Reading metadata from file: "mesh_input" =============================================================== CALCULATION PREPARATION ======================= =========================================================== No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 9 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) ----------------------------------------------------------- ** VOF METHOD ---------------------------------------- IVOFMT = -1 ( -1: inactive ) ( 0: active ) ** HOMOGENEOUS MIXTURE MODEL FOR CAVITATION ---------------------------------------- ICAVIT = -1 (-1: single phase flow ) ( 0: no vap./cond. model ) ( 1: Merkle's model ) ----------------------------------------------------------- ** 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 Mesh Velocity 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 ------ -- Phase continue : 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 ) IROEXT = 0 (density extrapolation (0: explicit (1: n+thetro with thetro=1/2 (2: n+thetro with thetro=1 THETRO = 0.00000E+00 (theta for density ((1+theta).new-theta.old IVIEXT = 0 (total viscosity extrapolation (0: explicit (1: n+thetvi with thetro=1/2 (2: n+thetvi with thetro=1 THETVI = 0.00000E+00 (theta for total viscosity ((1+theta).new-theta.old ICPEXT = 0 (specific heat extrapolation (0: explicit (1: n+thetcp with thetro=1/2 (2: n+thetcp with thetro=1 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 NITMAY = 10000 (Nb iter for iter resolution ) NSWRSY = 1 (Nb iter for rhs reconstr. ) NSWRGY = 100 (Nb iter for grad. reconstr. ) IMLIGY = -1 (Gradient limitation method ) IRCFLY = 1 (Conv. Diff. flow reconstr. ) ISCHCY = 1 (Convective scheme ) ISSTPY = 0 (Slope tet use ) IWARNY = 0 (Verbosity level ) NTCMXY = 1000 (Nb iter for steady convect. ) BLENCY = 0.00000E+00 (2nd order conv. scheme prop.) EPSILY = 0.10000E-07 (Iterative solver precision ) EPSRSY = 0.10000E-04 (rhs reconstruction precision) EPSRGY = 0.10000E-04 (Gradient reconstr. precision) CLIMGY = 0.15000E+01 (Coeff. for grad. limitation ) EXTRAY = 0.00000E+00 (Coeff. for grad. extrapolat.) COUMXY = 0.50000E+04 (Max CFL for convection ) EPSCVY = 0.10000E-07 (Precision for steady conv. ) YPLMXY = 0.20000E+03 (y+ max w. damping influence ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-07 1 Pressure 0.1000E-07 1 k 0.1000E-07 1 omega 0.1000E-07 1 Mesh Velocity 0.1000E-07 1 ------------------------------------ EPSILO = (resolution precision) IDIRCL = 0 ou 1 (shift diagonal if ISTAT=0 and no Dirichlet) ----------------------------------------------------------- ** 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. INPDT0 = 0 (1: 0 time step calcuation ) NTMABS = 100000 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) Velocity 0 Pressure 0 k 0 omega 0 Mesh Velocity 0 TurbVisc -999 CourantNb -999 FourierNb -999 total_pressure -999 Mesh displacemen -999 Mesh Visc -999 Yplus -999 -- -- -- --- 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 = 1 (1: activated ) 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: binary time dependency: fixed mesh output: every 2.00000e+00 s 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 Reading file: mesh_input Finished reading: mesh_input No "partition_input/domain_number_4" file available; ---------------------------------------------------------- Partitioning by space-filling curve: Morton (in bounding box). Number of cells per domain (histogramm): [ 1193197 ; 1193198 ] = 4 Partitioning finished (15.7 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 = 1193197 maximum value = 1193198 1 : [ 1193197 ; 1193198 ] = 4 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 1236350 maximum value = 1248224 1 : [ 1236350 ; 1238724 [ = 2 2 : [ 1238724 ; 1241099 [ = 0 3 : [ 1241099 ; 1243474 [ = 0 4 : [ 1243474 ; 1245849 [ = 1 5 : [ 1245849 ; 1248224 ] = 1 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 43153 maximum value = 55026 1 : [ 43153 ; 45527 [ = 2 2 : [ 45527 ; 47902 [ = 0 3 : [ 47902 ; 50276 [ = 0 4 : [ 50276 ; 52651 [ = 1 5 : [ 52651 ; 55026 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 3539346 maximum value = 3546653 1 : [ 3539346 ; 3540807 [ = 2 2 : [ 3540807 ; 3542268 [ = 0 3 : [ 3542268 ; 3543730 [ = 0 4 : [ 3543730 ; 3545191 [ = 0 5 : [ 3545191 ; 3546653 ] = 2 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 113599 maximum value = 119454 1 : [ 113599 ; 114770 [ = 1 2 : [ 114770 ; 115941 [ = 0 3 : [ 115941 ; 117112 [ = 0 4 : [ 117112 ; 118283 [ = 2 5 : [ 118283 ; 119454 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 6 1 : [ 3 ; 4 [ = 56 2 : [ 4 ; 5 [ = 10688 3 : [ 5 ; 6 ] = 4762046 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 3 maximum value = 3 1 : [ 3 ; 3 ] = 4 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: renumbering for cells: pre-numbering: none cells adjacent to ghost cells last: no numbering: none renumbering for interior faces: cell adjacency pre-ordering: lowest id first faces adjacent to ghost cells last: no numbering: none renumbering for boundary faces: numbering: none ---------------------------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 1231028 maximum value = 1240518 1 : [ 1231028 ; 1232926 [ = 2 2 : [ 1232926 ; 1234824 [ = 0 3 : [ 1234824 ; 1236722 [ = 0 4 : [ 1236722 ; 1238620 [ = 0 5 : [ 1238620 ; 1240518 ] = 2 Histogram of volume mesh matrix profile/lines per rank: minimum value = 77199 maximum value = 86148 1 : [ 77199 ; 78988 [ = 1 2 : [ 78988 ; 80778 [ = 1 3 : [ 80778 ; 82568 [ = 0 4 : [ 82568 ; 84358 [ = 0 5 : [ 84358 ; 86148 ] = 2 Mesh coordinates: minimum and maximum X : -1.1088000e-01 1.1091197e-01 Y : -2.9102868e-09 8.8000003e-02 Z : -5.4504001e-01 5.4504001e-01 Mesh Number of cells: 4772790 Number of interior faces: 14070081 Number of boundary faces: 468660 Number of vertices: 4990100 Groups: "FLUID.1" cells: 4041774 "FLUID_DEF" cells: 178596 "PERIODIC_YNEG" boundary faces: 48210 "PERIODIC_YPOS" boundary faces: 48210 "SIDE_1" boundary faces: 4059 "SIDE_2" boundary faces: 26235 "SIDE_3" boundary faces: 4059 "SIDE_4" boundary faces: 25542 "SIDE_PIPES" boundary faces: 304821 "TUBES_DEF_BL" cells: 356796 "TUBE_BL" cells: 195624 "TUBE_WALL" boundary faces: 7524 "WALL_DEF" interior faces: 36432 "WALL_TUBE_BL" interior faces: 7524 --- Information on the volumes Minimum control volume = 4.7149380e-12 Maximum control volume = 2.0080374e-07 Total volume for the domain = 1.9639272e-02 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 198 --> 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 % 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 (3.25 s) ALMAX = 0.26980E+00 (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 omega 0.9079E-01 0.9079E-01 Mesh Velocity 0.0000E+00 0.0000E+00 Mesh Velocity 0.0000E+00 0.0000E+00 Mesh Velocity 0.0000E+00 0.0000E+00 --------------------------------- dt 0.2000E-03 0.2000E-03 --------------------------------- ------------------------------------------------------------- ------------------------------------------------------------- COUPLING MODE FOR STRUCTURES ACTIVATED WITH NBSTRU = 1 STRUCTURE(S) NEWMARK COEFFICIENTS: ALPNMK = 0.0000E+00 BETNMK = 0.2500E+00 GAMNMK = 0.5000E+00 MONITORING FILES FOR STRUCTURES: IHISTR = 1 ( 1 : activated) CODE_ASTER COUPLING MODE NOT ACTIVATED NBASTE = 0 IMPLICIT SCHEME FOR COUPING ACTIVATED NB OF MAX INNER ITERATIONS : 20 CONVERGENCE THRESHOLD : 0.1000E-06 ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.000000000E+00 ALE INITIALIZATION ============================================================= ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1000E+04 0.1000E+04 molecular_viscos 0.1000E-02 0.1000E-02 turbulent_viscos 0.6609E+01 0.6609E+01 ----------------------------------------- --- Mesh viscosity (ALE method) ----------------------------------------- Property Min. value Max. value ----------------------------------------- mesh_viscosity 0.1197E+01 0.1000E+03 ----------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 4059 Smooth wall 5 460542 Rough wall 6 0 Symmetry 4 0 Free outlet 3 4059 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = -0.23098E-06 0.44444E-03 0.54504E+00 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 4059 0.000000000E+00 Smooth wall 5 460542 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 4059 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.00000E+00 0.00000E+00 Friction velocity uk : 0.10792E-02 0.11772E-01 Dimensionless distance yplus : 0.31457E-02 0.77682E+01 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 460542 Total number of wall faces : 460542 ------------------------------------------------------------ ** WALL DISTANCE ------------- Min distance = 0.29313E-05 Max distance = 0.43120E-01 (Distance calculation done in 376 iterations) ** 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 c k 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 c omega 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 c Mesh Velocity 0.00000E+00 0 0.00000E+00 0.00000E+00 0.00000E+00 c Mesh Velocity[X] 0.00000E+00 c Mesh Velocity[Y] 0.00000E+00 c Mesh Velocity[Z] 0.00000E+00 -------------------------------------------------------------------------------- ** Computed fields 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.0006 0.0006 0.0006 0.0006 v omega 0.090789 0.090789 0.090789 0.090789 v Mesh Velocity[X] 0 0 0 0 v Mesh Velocity[Y] 0 0 0 0 v Mesh Velocity[Z] 0 0 0 0 v ǁMesh Velocityǁ 0 0 0 0 v TurbVisc 6.6087 6.6087 6.6087 6.6087 v CourantNb 0 0 0 0 v FourierNb 1.6967 44861 551.47 6.9527 v total_pressure 1.0135e+05 1.12e+05 1.0669e+05 1.0667e+05 v Mesh Visc 1.1966 100 13.485 1.965 ** Computed fields on boundary_faces --------------------------------- field minimum maximum set mean spatial mean - ----- -------------- -------------- -------------- -------------- v Yplus 0 7.7682 1.9805 3.2421 ** Computed fields on vertices --------------------------- field minimum maximum set mean - -------------------- -------------- -------------- -------------- v Mesh displacement[X] 0 0 0 v Mesh displacement[Y] 0 0 0 v Mesh displacement[Z] 0 0 0 v ǁMesh displacementǁ 0 0 0 ** Computed values on cells ------------------------ criterion minimum maximum set mean spatial mean --------------- -------------- -------------- -------------- -------------- Courant/Fourier 1.6967 44861 551.47 6.9527 TIME FOR ALE INITIALIZATION: 0.82073E+02 =============================================================== INSTANT 0.200000000E-03 TIME STEP NUMBER 1 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 4059 0.000000000E+00 Smooth wall 5 460542 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 4059 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.00000E+00 0.00000E+00 Friction velocity uk : 0.10792E-02 0.11772E-01 Dimensionless distance yplus : 0.31457E-02 0.77682E+01 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 460542 Total number of wall faces : 460542 ------------------------------------------------------------ ** WALL DISTANCE ------------- Min distance = 0.29313E-05 Max distance = 0.43120E-01 (Distance calculation done in 376 iterations) SIGSEGV signal (forbidden memory area access) intercepted! Call stack: 1: 0x407df2 (cs_solver) 2: 0x2ada3e05e2e2 (libsaturne.so.5) 3: 0x2ada3e02f6bf (libsaturne.so.5) 4: 0x2ada3e02f00f (libsaturne.so.5) 5: 0x2ada42738b35 <__libc_start_main+0xf5> (libc.so.6) 6: 0x407c49 <> (cs_solver) End of stack