command: ./cs_solver *************************************************************** (R) Code_Saturne Version 5.0.9 Copyright (C) 1998-2018 EDF S.A., France revision 5.0.9-9b87d1bcb build Tue Dec 3 19:30:53 2019 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 Oct 16 14:17:12 2020 System: Linux 2.6.32-279.14.1.el6.x86_64 (CentOS release 6.3 (Final)) Machine: compute-0-10.local Processor: model name : AMD Opteron(tm) Processor 6386 SE Memory: 129190 MB User: a_rasam () Directory: /share/users/a_rasam/Zangeneh/AIRFOIL/RESU/20201016-1416 OpenMP threads: 2 Processors/node: 64 External libraries for partitioning: METIS Reading metadata from file: "mesh_input" =============================================================== CALCULATION PREPARATION ======================= =========================================================== @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING : WHEN READING INPUT DATA @ ========= @ CHOICE OF TIME-SCHEME @ @ TIME-SCHEME FOR VELOCITY IS SECOND ORDER @ (THETAV = 0.50E+00) @ CERTAIN 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 given through the User Interface @ or in cs_user_parameters.f90. @ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 4 (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 ---------------------------- 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 ) ----------------------------------------------------------- ** STOKES ------ -- Phase continue : ISTMPF = 2 (time scheme for flow (0: explicit (THETFL = 0 ) (1: std scheme (Saturne 1.0 ) (2: 2nd-order (THETFL = 0.5 ) THETFL = 0.50000E+00 (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.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 ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-04 1 Pressure 0.1000E-04 1 ------------------------------------ EPSILO = (resolution precision) IDIRCL = 0 ou 1 (shift diagonal if ISTAT=0 and no Dirichlet) ----------------------------------------------------------- ** 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. INPDT0 = 0 (1: 0 time step calcuation ) NTMABS = 500000 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 5 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- Log files NTLIST = 10 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) Velocity 0 Pressure 0 Laminar Viscosit -999 Turb Viscosity -999 CFL -999 Fourier Number -999 Total Pressure -999 Csdyn2 -999 CF -999 CP -999 Yplus -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 7 (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 ) 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: time dependency: fixed mesh output: 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 ------------------------------------------------------- Joining number 1: Selection criteria: "INTERFACE_LAM or INTERFACE_TURB" Parameters for the joining operation: Shortest incident edge fraction: 0.01000 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: 3153668 Number of interior faces: 9242654 Number of boundary faces: 436700 Number of vertices: 3373386 Global number of boundary faces selected for joining: 48852 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 228036 Vertex-Vertex intersections: 77354 Other intersections: 150682 Joining operation is non-conforming. Global number of new vertices to create: 150682 Equivalences broken for 0 edges. Edge intersections and vertex creation done. Pre-merge for 152002 global element couples. Max size of a merge set of vertices: 2 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: 43424 Global configuration after the joining operation: Global number of border faces to add: 0 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: 3153668 Number of interior faces: 9286078 Number of boundary faces: 387848 Number of vertices: 3367954 Joining 1 completed (12.6 s). ------------------------------------------------------- Joining number 2: Periodicity type: translation Transformation matrix: 1 0 0 0 0 1 0 0 0 0 1 0.018 Selection criteria: "PER1 or PER2" 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: 3153668 Number of interior faces: 9286078 Number of boundary faces: 387848 Number of vertices: 3367954 Global number of boundary faces selected for joining: 341964 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 2765470 Vertex-Vertex intersections: 2765470 Other intersections: 0 Joining operation is conforming. Edge intersections and vertex creation done. Pre-merge for 172346 global element couples. No need to merge 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: 512946 Global configuration after the joining operation: Global number of border faces to add: 0 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: 3153668 Number of interior faces: 9628042 Number of boundary faces: 45884 Number of vertices: 3367954 Joining 2 completed (252 s). All joining operations successfully finished: Wall clock time: 265 ---------------------------------------------------------- Composing periodicities Halo construction with extended neighborhood ============================================ Face interfaces creation Definition of periodic vertices Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Extended neighborhood structures definition Number of cells: 3153668 ---------------------------------------------------------- Number of cells + halo cells: 3495632 ---------------------------------------------------------- Number of ghost cells: 341964 in the standard neighborhood: 341964 ---------------------------------------------------------- Number of interior faces: 9628042 ---------------------------------------------------------- Number of boundary faces: 45884 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 4 maximum value = 13 1 : [ 4 ; 5 [ = 45868 2 : [ 5 ; 7 [ = 3102372 3 : [ 7 ; 9 [ = 0 4 : [ 9 ; 11 [ = 0 5 : [ 11 ; 13 ] = 5428 ---------------------------------------------------------- Synchronizing cell families Writing file: mesh_output Finished writing: mesh_output 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: multipass renumbering for boundary faces: numbering: none Number of thread pools for interior faces : 2 Number of threading groups for interior faces : 7 Estimated thread imbalance for interior faces : 2.62708e-02 ---------------------------------------------------------- Matrix bandwidth for volume mesh : 3153668 Matrix profile/lines for volume mesh : 225088 Mesh coordinates: minimum and maximum X : -5.3999822e+00 7.2001574e+00 Y : -6.0000000e+00 6.0000000e+00 Z : -3.8743020e-10 1.8000000e-02 Mesh Number of cells: 3153668 Number of interior faces: 9628042 Number of boundary faces: 45884 Number of vertices: 3367954 Periodic faces (which are also interior faces): Periodicity 1: 170982 face couples Groups: "INLET" boundary faces: 5428 "INTERFACE_LAM" interior faces: 43424 "INTERFACE_TURB" interior faces: 43424 "OUTLET" boundary faces: 4648 "PER1" interior faces: 170982 "PER2" interior faces: 170982 "WALL" boundary faces: 35808 --- Information on the volumes Minimum control volume = 3.0682996e-13 Maximum control volume = 1.1589449e-04 Total volume for the domain = 2.4427596e+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 % Criteria 6: User Specific Tag: Number of bad cells detected: 3153668 --> 100 % Warning: -------- Mesh quality issue based on user criteria has been detected The mesh should be re-considered using the listed criteria. Computing geometric quantities (15.2 s) ALMAX = 0.13468E+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 Warning: the number of variables or properties has been changed relative to the restart file. current: 2 model and 0 user variables 12 model and 2 user properties previous: 14 model and 0 user variables 1 model and 0 user properties The computation continues, with a partial restart. Reading dimensions complete restart/main: section "cavitation" not present. Reading the previous time step number (restarting computation) NTPABS = 115073 Reading the previous time step number (restarting computation) TTPABS = 0.2000E+05 Reading options complete Read variables from restart: restart/main Reading complete READING THE AUXILIARY RESTART FILE Start reading restart/auxiliary: section "cavitation" not present. Fin de la lecture des options 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: section "time_moments:sizes" not present. End reading ----------------------------------------------------------- ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Velocity -0.8197E+02 0.1714E+03 Velocity -0.6621E+02 0.1136E+03 Velocity -0.6521E+02 0.5616E+02 Pressure -0.7414E+04 0.6003E+04 --------------------------------- dt 0.1200E-05 0.1200E-05 --------------------------------- ------------------------------------------------------------- ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================