command: /opt/saturne-3.0/libexec/code_saturne/cs_solver --param LES.xml --mpi *************************************************************** (R) Code_Saturne Version 3.0.0 Copyright (C) 1998-2013 EDF S.A., France build Thu 02 May 2013 08:06:16 AM PKT MPI version 2.1 (Open MPI 1.4.3) 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: Tue 14 Jun 2016 10:05:15 AM PKT System: Linux 3.2.0-60-generic Machine: fahad-virtual-machine Processor: model name : AMD Opteron(tm) Processor 6348 Memory: 122453 MB User: fahad (Fahad) Directory: /home/fahad/Desktop/Code_Saturne_CircularCylinder/CIRCULARCYLINDER/LES/RESU/20160614-1001 MPI ranks: 16 (appnum attribute: 0) 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: ParMETIS 3.1 PT-SCOTCH 1.12.2147483631 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 ---------- --- Geometry NPRFML = 1 (Nb max. family properties ) NFML = 7 (Nb families ) --- Physics NVAR = 4 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) NPROCE = 6 (Nb cell properties ) NPROFA = 1 (Nb internal face properties ) NPROFB = 2 (Nb boundary face properties ) ----------------------------------------------------------- ** PHYSICAL PROPERTIES ------------------- GX = 0.00000E+00 (Gravity x component ) GY = 0.00000E+00 (Gravity y component ) GZ = 0.00000E+00 (Gravity z component ) OMEGAX = 0.00000E+00 (Rotation vector x component ) OMEGAY = 0.00000E+00 (Rotation vector y component ) OMEGAZ = 0.00000E+00 (Rotation vector z component ) ICORIO = 0 (Coriolis source terms ) -- Continuous phase: RO0 = 0.11786E+01 (Reference density ) VISCL0 = 0.18300E-04 (Ref. molecular dyn. visc. ) CP0 = 0.10172E+04 (Ref. specific heat ) ICP = 0 (> 0: variable CP (usphyv) ) P0 = 0.10132E+06 (Ref. total pressure ) PRED0 = 0.00000E+00 (Ref. reduced pressure ) T0 = 0.29315E+03 (Ref. temperature ) IROVAR = 0 (Density constant(0) or not(1) IVIVAR = 0 (Molec. visc cst.(0) or not(1) Initial reference point for pressure XYZP0 = 0.00000E+00 0.00000E+00 0.00000E+00 ----------------------------------------------------------- ** TURBULENCE ---------- --- Continuous phase: - Commons ITURB = 0 (Turbulence model ) IDEUCH = 0 (0: one-scale model ) (1: two-scale model ) (2: invariant wall function ) IWALLT = 0 (Exch. coeff. correlation ) (0: not activated ) (1: activated ) YPLULI = 0.23810E+01 (Limit Y+ ) ILOGPO = 1 (0: power law (forbidden for k-epsilon) ) (1: one-scale log law ) IGRHOK = 0 (1: computed Grad(rho k) ) ISCALT = -1 (Temperature salar number ) --- Constants - Commons XKAPPA = 0.42000E+00 (Von Karman constant ) CSTLOG = 0.52000E+01 (U+=Log(y+)/kappa +CSTLOG ) APOW = 0.83000E+01 (U+=APOW (y+)**BPOW (W&W law)) BPOW = 0.14286E+00 (U+=APOW (y+)**BPOW (W&W law)) ----------------------------------------------------------- ** SECONDARY VISCOSITY ------------------- --- Continuous phase: IVISSE = 1 (1: accounted for ) ----------------------------------------------------------- ** TIME STEPPING ------------- UNSTEADY ALGORITHM --- Time step parameters IDTVAR = 2 (0 cst; 1,2 var (t, t-space ) IPTLRO = 0 (1: rho-related DT clipping ) COUMAX = 0.10000E+01 (Maximum target CFL ) FOUMAX = 0.10000E+02 (Maximum target Fourier ) VARRDT = 0.10000E+00 (For var. DT, max. increase ) DTMIN = 0.10000E-04 (Minimum time step ) DTMAX = 0.10000E+00 (Maximum time step ) DTREF = 0.10000E-03 (Reference time step ) With a non-constant time step (IDTVAR = 1 or 2), when the value of COUMAX or FOUMAX is negative or zero, the associated time step limitation (for CFL and Fourier respectively) is ignored. --- Frozen velocity field ICCVFG = 0 (1: frozen velocity field ) --- Per-variable properties ------------------------------------ Variable ISTAT CDTVAR ------------------------------------ Pressure 0 0.1000E+01 VelocityX 1 0.1000E+01 VelocityY 1 0.1000E+01 VelocityZ 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 ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- --------------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV --------------------------------------------------------------------- Pressure 0 1 1 1 0 0.00E+00 0.10E+01 VelocityX 1 1 1 1 0 0.10E+01 0.10E+01 VelocityY 1 1 1 1 0 0.10E+01 0.10E+01 VelocityZ 1 1 1 1 0 0.10E+01 0.10E+01 ------------------------------------------------------------- ICONV = 0 ou 1 (1 for convection active ) IDIFF = 0 ou 1 (1 for total diffusion active) IDIFFT = 0 ou 1 (1 for turbulent diff. active) ISCHCV = 0 ou 1 (SOLU or CD ) ISSTPC = 0 ou 1 (1: no slope test ) BLENCV = [0.;1.] (1-upwind proportion ) THETAV = [0.;1.] (0.5 Crank-Nicolson/AB ) (theta for convection- ) (diffusion terms uses ) ((1-theta).old+theta.new ) ----------------------------------------------------------- ** STOKES ------ IVELCO = 1 (0: segregated solving of the velocity components during the prediction step 1: coupled solving of the velocity components ) IPHYDR = 0 (1: account for explicit balance between pressure gradient, gravity source terms, and head losses ) ICALHY = 0 (1: compute hydrastatic pressure for Dirichlet conditions for pressure on outlet ) IPRCO = 1 (1: pressure-continuity ) IPUCOU = 0 (1: reinforced U-P coupling ) NTERUP = 1 (n: n sweeps on navsto for velocity/pressure coupling ) -- Continuous phase: IREVMC = 0 (Velocity reconstruction mode) RELAXV = 0.10000E+01 for pressure (relaxation) ARAK = 0.10000E+01 (Arakawa factor ) 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 ) ----------------------------------------------------------- ** GRADIENTS CALCULATION --------------------- IMRGRA = 0 (Reconstruction mode ) ANOMAX = 0.78540E+00 (Non-ortho angle: limit for ) (least squares ext. neighbors) ------------------------------------------------------------------- Variable NSWRGR NSWRSM EPSRGR EPSRSM EXTRAG ------------------------------------------------------------------- Pressure 100 5 0.1000E-04 0.1000E-06 0.0000E+00 VelocityX 100 10 0.1000E-04 0.1000E-06 0.0000E+00 VelocityY 100 10 0.1000E-04 0.1000E-06 0.0000E+00 VelocityZ 100 10 0.1000E-04 0.1000E-06 0.0000E+00 ----------------------------------------------------------- ------------------------------------------- Variable IRCFLU IMLIGR CLIMGR ------------------------------------------- Pressure 1 -1 0.1500E+01 VelocityX 1 -1 0.1500E+01 VelocityY 1 -1 0.1500E+01 VelocityZ 1 -1 0.1500E+01 ----------------------------------- NSWRGR = (nb sweep gradient reconstr. ) NSWRSM = (nb sweep rhs reconstrcution ) EPSRGR = (grad. reconstruction prec. ) EPSRSM = (rhs reconstruction prec. ) EXTRAG = [0.;1.] (gradients extrapolation ) IRCFLU = 0 ou 1 (flow reconstruction ) IMLIGR = < 0, 0 ou 1 (gradient limitation method ) CLIMGR = > 1 ou 1 (gradient limitation coeff. ) ----------------------------------------------------------- ** FACE INTERPOLATION ------------------ IMVISF = 0 (0 arithmetic ) ----------------------------------------------------------- ** CALCULATION OF TEMPORAL MEANS (MOMENTS) --------------------------------------- NBMOMT = 0 (Number of moments ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- -------------------------------------------------- Variable IRESOL NITMAX EPSILO IDIRCL -------------------------------------------------- Pressure 0 10000 0.1000E-07 1 VelocityX -1 10000 0.1000E-07 1 VelocityY -1 10000 0.1000E-07 1 VelocityZ -1 10000 0.1000E-07 1 ----------------------------------- IRESOL = -1 (automatic solver choice ) IPOL*1000 + 0 (p conjuguate gradient ) 1 (Jacobi ) IPOL*1000 + 2 (bicgstab ) avec IPOL (preconditioning degree ) NITMAX = (max number of iterations ) EPSILO = (resolution precision ) IDIRCL = 0 ou 1 (shift diagonal if ISTAT=0 and no Dirichlet ) ----------------------------------------------------------- ** MULTIGRID --------- NCEGRM = 30 (Max nb cells coarsest grid ) NGRMAX = 25 (Max number of levels ) -------------------------------------- Variable IMGR NCYMAX NITMGF -------------------------------------- Pressure 1 100 10 VelocityX 0 100 10 VelocityY 0 100 10 VelocityZ 0 100 10 ------------------------------ IMGR = 0 ou 1 (1: multigrid activated ) NCYMAX = (Max number of cycles ) NITMGF = (Max nb iter on coarsest grid) Multigrid rank merge parameters: mean coarse cells merge threshold: 300 total coarse cells merge threshold: 500 minimum ranks merge threshold: 1 merge stride: 1 Multigrid coarsening type: algebraic, face traveral by criteria ----------------------------------------------------------- ** 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 = 5 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Post-processed variables Number Name 2 Pressure 3 VelocityX 4 VelocityY 5 VelocityZ 8 CourantNb 9 FourierNb 10 total_pressure 11 LocalTime -- -- --- Probe history files NTHIST = -1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) NCAPT = 0 (Number of probes ) NTHSAV = -1 (Checkpoint frequency ) Number Name Nb. probes (-1: all) 2 Pressure -1 3 VelocityX -1 4 VelocityY -1 5 VelocityZ -1 6 Density -1 7 LamVisc -1 8 CourantNb -1 9 FourierNb -1 10 total_pressure -1 11 LocalTime -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Pressure 0 3 VelocityX 0 4 VelocityY 0 5 VelocityZ 0 8 CourantNb -999 9 FourierNb -999 10 total_pressure -999 11 LocalTime -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 1 (Force exerted at boundary) ipstyp = 1 (Yplus at boundary) ipstft = 1 (Thermal flow at boundary) ipsttb = 0 (Temperature at boundary) ipstnu = 0 (Nusselt at boundary) ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated ) NALINF = 0 (Fluid initialization iterations) ----------------------------------------------------------- Reading file: mesh_input Finished reading: mesh_input No "partition_input/domain_number_16" file available; ---------------------------------------------------------- Partitioning 11275730 cells to 16 domains on 16 ranks (ParMETIS_V3_PartKway). Total number of faces on parallel boundaries: 6610 wall-clock time: 8.938235 s Number of cells per domain (histogramm): [ 633739 ; 642772 [ = 1 [ 642772 ; 651806 [ = 0 [ 651806 ; 660840 [ = 0 [ 660840 ; 669874 [ = 0 [ 669874 ; 678908 [ = 0 [ 678908 ; 687941 [ = 0 [ 687941 ; 696975 [ = 2 [ 696975 ; 706009 [ = 3 [ 706009 ; 715043 [ = 5 [ 715043 ; 724077 ] = 5 Writing file: partition_output/domain_number_16 Finished writing: partition_output/domain_number_16 Partitioning finished (10.9 s) Postprocessing output writers: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: binary time dependency: fixed mesh output: every 1 time steps ---------------------------------------------------------- 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 = 633739 maximum value = 724077 1 : [ 633739 ; 651806 [ = 1 2 : [ 651806 ; 669874 [ = 0 3 : [ 669874 ; 687941 [ = 0 4 : [ 687941 ; 706009 [ = 5 5 : [ 706009 ; 724077 ] = 10 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 634302 maximum value = 724690 1 : [ 634302 ; 652379 [ = 1 2 : [ 652379 ; 670457 [ = 0 3 : [ 670457 ; 688534 [ = 0 4 : [ 688534 ; 706612 [ = 5 5 : [ 706612 ; 724690 ] = 10 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 563 maximum value = 1411 1 : [ 563 ; 732 [ = 9 2 : [ 732 ; 902 [ = 2 3 : [ 902 ; 1071 [ = 2 4 : [ 1071 ; 1241 [ = 0 5 : [ 1241 ; 1411 ] = 3 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 947884 maximum value = 1083251 1 : [ 947884 ; 974957 [ = 1 2 : [ 974957 ; 1002030 [ = 0 3 : [ 1002030 ; 1029104 [ = 0 4 : [ 1029104 ; 1056177 [ = 5 5 : [ 1056177 ; 1083251 ] = 10 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 1273490 maximum value = 1454980 1 : [ 1273490 ; 1309788 [ = 1 2 : [ 1309788 ; 1346086 [ = 0 3 : [ 1346086 ; 1382384 [ = 0 4 : [ 1382384 ; 1418682 [ = 5 5 : [ 1418682 ; 1454980 ] = 10 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 2 maximum value = 3 1 : [ 2 ; 3 ] = 11275730 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 2 maximum value = 10 1 : [ 2 ; 3 [ = 11 2 : [ 3 ; 5 [ = 3 3 : [ 5 ; 6 [ = 0 4 : [ 6 ; 8 [ = 1 5 : [ 8 ; 10 ] = 1 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: Checking mesh renumbering for threads: ------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 633690 maximum value = 724028 1 : [ 633690 ; 651757 [ = 1 2 : [ 651757 ; 669825 [ = 0 3 : [ 669825 ; 687892 [ = 0 4 : [ 687892 ; 705960 [ = 5 5 : [ 705960 ; 724028 ] = 10 Histogram of volume mesh matrix profile/lines per rank: minimum value = 1425 maximum value = 5918 1 : [ 1425 ; 2323 [ = 2 2 : [ 2323 ; 3222 [ = 5 3 : [ 3222 ; 4120 [ = 1 4 : [ 4120 ; 5019 [ = 1 5 : [ 5019 ; 5918 ] = 7 Mesh coordinates: minimum and maximum X : -7.0000000e+00 1.5000000e+01 Y : -5.0000000e+00 5.0000000e+00 Z : 0.0000000e+00 1.0000000e+00 Mesh Number of cells: 11275730 Number of interior faces: 16863169 Number of boundary faces: 22652312 Number of vertices: 11376582 Groups: "Back" boundary faces: 11275730 "Bottom" boundary faces: 275 "Front" boundary faces: 11275730 "Left" boundary faces: 125 "Top" boundary faces: 275 "right" boundary faces: 125 --- Information on the volumes Minimum control volume = 5.5578274e-11 Maximum control volume = 5.2654104e-03 Total volume for the domain = 2.1685841e+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: 2113 --> 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 (4.22 s) ALMAX = 0.59972E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Pressure 0.0000E+00 0.0000E+00 VelocityX 0.0000E+00 0.0000E+00 VelocityY 0.0000E+00 0.0000E+00 VelocityZ 0.0000E+00 0.0000E+00 --------------------------------- LocalTime 0.1000E-03 0.1000E-03 --------------------------------- ------------------------------------------------------------- ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.100000000E-03 TIME STEP NUMBER 1 ============================================================= ----------------------------------------- Property Min. value Max. value ----------------------------------------- Density 0.1179E+01 0.1179E+01 LamVisc 0.1830E-04 0.1830E-04 ----------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 125 Smooth wall 5 0 Rough wall 6 100052 Symmetry 4 22552010 Free outlet 3 125 Undefined 1 0 ------------------------------------------------------------------------- Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = 0.15000E+02 0.61542E-11 0.50000E+00 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 125 0.000000000E+00 Smooth wall 5 0 0.000000000E+00 Rough wall 6 100052 0.000000000E+00 Symmetry 4 22552010 0.000000000E+00 Free outlet 3 125 0.000000000E+00 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR ROUGH WALLS -------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.14672E-11 0.93805E-10 Friction velocity uk : 0.14672E-11 0.93805E-10 Rough dimensionless dist yplus : 0.44874E-02 0.33143E+00 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 ------------------------------------------------------------ ** INFORMATION ON CONVERGENCE -------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.10351E+02 8910 0.24802E-03 0.98187E+00 c VelocityX 0.84801E-02 5 0.49500E-12 0.94168E+06 c VelocityY 0.84801E-02 5 0.49500E-12 0.65982E+04 c VelocityZ 0.84801E-02 5 0.49500E-12 0.13542E-23 --------------------------------------------------------------- ** INFORMATION ON VARIABLES ------------------------ --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure 0.13924E+04 0.24269E+07 -- -- v VelocityX -0.12487E+00 0.29212E+02 -- -- v VelocityY -0.14914E+02 0.14642E+02 -- -- v VelocityZ -0.11451E-12 0.11085E-12 -- -- v CourantNb 0.00000E+00 0.00000E+00 0 0 v FourierNb 0.18792E-05 0.17161E+02 0 0 v total_pressu 0.10272E+06 0.25282E+07 0 0 v LocalTime 0.10000E-04 0.11000E-03 2 0 --------------------------------------------------------------- ** INFORMATION ON CLIPPINGS ------------------------ --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- ** INFORMATION ON THE TIME STEP ----------------------------- --------------------------------------------------------------- Criterion Value at xc yc zc --------------------------------------------------------------- Courant min 0.00000E+00 -0.94559E+00 0.32543E+00 0.50000E+00 Courant max 0.00000E+00 -0.94559E+00 0.32543E+00 0.50000E+00 Fourier min 0.18792E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Fourier max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 Cou_Fou min 0.18792E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Cou_Fou max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 --------------------------------------------------------------- Dt min 0.10000E-04 0.98501E+00 0.17249E+00 0.50000E+00 Dt max 0.11000E-03 -0.99416E+00 0.10796E+00 0.50000E+00 --------------------------------------------------------------- TIME FOR THE TIME STEP 1: 0.19530E+04 =============================================================== INSTANT 0.200000000E-03 TIME STEP NUMBER 2 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 125 -0.117862000E+03 Smooth wall 5 0 0.000000000E+00 Rough wall 6 100052 0.000000000E+00 Symmetry 4 22552010 0.000000000E+00 Free outlet 3 125 0.117861727E+03 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR ROUGH WALLS -------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.21458E-03 0.27062E+03 Friction velocity uk : 0.21458E-03 0.27062E+03 Rough dimensionless dist yplus : 0.44874E-02 0.33143E+00 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 ------------------------------------------------------------ ** INFORMATION ON CONVERGENCE -------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.38417E+03 625 0.27578E-04 0.99346E+00 c VelocityX 0.62610E+05 303 0.34244E-01 0.41402E+07 c VelocityY 0.62610E+05 303 0.34244E-01 0.38766E+06 c VelocityZ 0.62610E+05 303 0.34244E-01 0.21831E-25 --------------------------------------------------------------- ** INFORMATION ON VARIABLES ------------------------ --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.74799E+07 0.29342E+06 -- -- v VelocityX -0.17530E+07 0.22103E+07 -- -- v VelocityY -0.27362E+05 0.66148E+06 -- -- v VelocityZ -0.15399E-09 0.35983E-11 -- -- v CourantNb 0.64542E-03 0.19458E+02 0 0 v FourierNb 0.20672E-05 0.17161E+02 0 0 v total_pressu-0.73786E+07 0.39475E+06 0 0 v LocalTime 0.10000E-04 0.12100E-03 7540223 0 --------------------------------------------------------------- ** INFORMATION ON CLIPPINGS ------------------------ --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- ** INFORMATION ON THE TIME STEP ----------------------------- --------------------------------------------------------------- Criterion Value at xc yc zc --------------------------------------------------------------- Courant min 0.64542E-03 -0.10000E+01 -0.26305E-04 0.50000E+00 Courant max 0.19458E+02 -0.16183E+00 -0.98684E+00 0.50000E+00 Fourier min 0.20672E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Fourier max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 Cou_Fou min 0.21674E-01 0.89047E+01 0.20709E+00 0.50000E+00 Cou_Fou max 0.21552E+02 -0.16183E+00 -0.98684E+00 0.50000E+00 --------------------------------------------------------------- Dt min 0.10000E-04 -0.93315E+00 0.35949E+00 0.50000E+00 Dt max 0.12100E-03 -0.10000E+01 0.13502E-02 0.50000E+00 --------------------------------------------------------------- TIME FOR THE TIME STEP 2: 0.47864E+03 =============================================================== INSTANT 0.300000000E-03 TIME STEP NUMBER 3 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 125 -0.117862000E+03 Smooth wall 5 0 0.000000000E+00 Rough wall 6 100052 0.000000000E+00 Symmetry 4 22552010 0.000000000E+00 Free outlet 3 125 0.117862135E+03 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR ROUGH WALLS -------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.27966E-03 0.91301E+08 Friction velocity uk : 0.27966E-03 0.91301E+08 Rough dimensionless dist yplus : 0.44874E-02 0.33143E+00 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 ------------------------------------------------------------ ** INFORMATION ON CONVERGENCE -------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.24775E+03 615 0.15251E-05 0.13450E+01 c VelocityX 0.24132E+07 17600 0.14572E-03 0.69484E+07 c VelocityY 0.24132E+07 17600 0.14572E-03 0.97910E+06 c VelocityZ 0.24132E+07 17600 0.14572E-03 0.18616E-25 --------------------------------------------------------------- ** INFORMATION ON VARIABLES ------------------------ --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.58742E+07 0.10515E+07 -- -- v VelocityX -0.65909E+06 0.16464E+06 -- -- v VelocityY -0.38231E+06 0.41477E+06 -- -- v VelocityZ -0.13380E-10 0.73289E-10 -- -- v CourantNb 0.34887E-03 0.30861E+06 0 0 v FourierNb 0.22739E-05 0.17161E+02 0 0 v total_pressu-0.57729E+07 0.11528E+07 0 0 v LocalTime 0.10000E-04 0.13310E-03 7533389 0 --------------------------------------------------------------- ** INFORMATION ON CLIPPINGS ------------------------ --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- ** INFORMATION ON THE TIME STEP ----------------------------- --------------------------------------------------------------- Criterion Value at xc yc zc --------------------------------------------------------------- Courant min 0.34887E-03 0.10000E+01 -0.26327E-04 0.50000E+00 Courant max 0.30861E+06 -0.90219E+00 -0.43133E+00 0.50000E+00 Fourier min 0.22739E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Fourier max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 Cou_Fou min 0.23842E-01 0.89047E+01 0.20709E+00 0.50000E+00 Cou_Fou max 0.30863E+06 -0.90219E+00 -0.43133E+00 0.50000E+00 --------------------------------------------------------------- Dt min 0.10000E-04 -0.94189E+00 0.33594E+00 0.50000E+00 Dt max 0.13310E-03 -0.10000E+01 0.13502E-02 0.50000E+00 --------------------------------------------------------------- TIME FOR THE TIME STEP 3: 0.73796E+04 =============================================================== INSTANT 0.400000000E-03 TIME STEP NUMBER 4 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 125 -0.117862000E+03 Smooth wall 5 0 0.000000000E+00 Rough wall 6 100052 0.000000000E+00 Symmetry 4 22552010 0.000000000E+00 Free outlet 3 125 0.117862028E+03 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR ROUGH WALLS -------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.52706E-04 0.69454E+07 Friction velocity uk : 0.52706E-04 0.69454E+07 Rough dimensionless dist yplus : 0.44874E-02 0.33143E+00 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 ------------------------------------------------------------ ** INFORMATION ON CONVERGENCE -------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.82076E+02 711 0.96147E-05 0.66073E+00 c VelocityX 0.42210E+07 17078 0.32768E-04 0.74308E+05 c VelocityY 0.42210E+07 17078 0.32768E-04 0.19823E+06 c VelocityZ 0.42210E+07 17078 0.32768E-04 0.56168E-26 --------------------------------------------------------------- ** INFORMATION ON VARIABLES ------------------------ --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.43242E+07 0.11278E+07 -- -- v VelocityX -0.50646E+06 0.49320E+05 -- -- v VelocityY -0.28018E+06 0.80354E+05 -- -- v VelocityZ -0.15483E-10 0.38116E-10 -- -- v CourantNb 0.33668E-03 0.14718E+06 0 0 v FourierNb 0.25013E-05 0.17161E+02 0 0 v total_pressu-0.42229E+07 0.12291E+07 0 0 v LocalTime 0.10000E-04 0.14641E-03 7532691 0 --------------------------------------------------------------- ** INFORMATION ON CLIPPINGS ------------------------ --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- ** INFORMATION ON THE TIME STEP ----------------------------- --------------------------------------------------------------- Criterion Value at xc yc zc --------------------------------------------------------------- Courant min 0.33668E-03 -0.10000E+01 0.29129E-03 0.50000E+00 Courant max 0.14718E+06 -0.90222E+00 -0.43128E+00 0.50000E+00 Fourier min 0.25013E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Fourier max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 Cou_Fou min 0.26226E-01 0.89047E+01 0.20709E+00 0.50000E+00 Cou_Fou max 0.14720E+06 -0.90219E+00 -0.43133E+00 0.50000E+00 --------------------------------------------------------------- Dt min 0.10000E-04 -0.94881E+00 0.31587E+00 0.50000E+00 Dt max 0.14641E-03 -0.10000E+01 0.13502E-02 0.50000E+00 --------------------------------------------------------------- TIME FOR THE TIME STEP 4: 0.77039E+04 =============================================================== INSTANT 0.500000000E-03 TIME STEP NUMBER 5 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 125 -0.117862000E+03 Smooth wall 5 0 0.000000000E+00 Rough wall 6 100052 0.000000000E+00 Symmetry 4 22552010 0.000000000E+00 Free outlet 3 125 0.117862084E+03 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR ROUGH WALLS -------------------------------------- ------------------------------------------------------------ Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.39780E-04 0.10427E+08 Friction velocity uk : 0.39780E-04 0.10427E+08 Rough dimensionless dist yplus : 0.44874E-02 0.33143E+00 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 ------------------------------------------------------------ Write final restart files checkpoint at iteration 5, Physical time 0.50000E-03 ** Writing the main restart file ----------------------------- Start writing End writing the dimensions End writing the options End writing the main variables End writing ** Writing the auxilliary restart file ----------------------------------- Start writing End writing the dimensions and options End writing the physical properties End writing the time step End writing the mass fluxes End writing the boundary conditions End writing Time for restart files: 0.10358E+02 ** INFORMATION ON CONVERGENCE -------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.43058E+02 531 0.13012E-04 0.48370E+00 c VelocityX 0.15331E+07 30895 0.49287E-04 0.30086E+05 c VelocityY 0.15331E+07 30895 0.49287E-04 0.17738E+05 c VelocityZ 0.15331E+07 30895 0.49287E-04 0.41605E-26 --------------------------------------------------------------- ** INFORMATION ON VARIABLES ------------------------ --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.30826E+07 0.14696E+07 -- -- v VelocityX -0.33814E+06 0.38694E+05 -- -- v VelocityY -0.18810E+06 0.27257E+05 -- -- v VelocityZ -0.18336E-10 0.21093E-10 -- -- v CourantNb 0.91439E-04 0.15142E+06 0 0 v FourierNb 0.27514E-05 0.17161E+02 0 0 v total_pressu-0.29813E+07 0.15709E+07 0 0 v LocalTime 0.10000E-04 0.16105E-03 7525525 0 --------------------------------------------------------------- ** INFORMATION ON CLIPPINGS ------------------------ --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- ** INFORMATION ON THE TIME STEP ----------------------------- --------------------------------------------------------------- Criterion Value at xc yc zc --------------------------------------------------------------- Courant min 0.91439E-04 0.10000E+01 -0.26327E-04 0.50000E+00 Courant max 0.15142E+06 -0.90219E+00 -0.43133E+00 0.50000E+00 Fourier min 0.27514E-05 -0.48000E+01 -0.49794E+01 0.50000E+00 Fourier max 0.17161E+02 -0.90219E+00 -0.43133E+00 0.50000E+00 Cou_Fou min 0.28849E-01 0.89047E+01 0.20709E+00 0.50000E+00 Cou_Fou max 0.15144E+06 -0.90219E+00 -0.43133E+00 0.50000E+00 --------------------------------------------------------------- Dt min 0.10000E-04 -0.94693E+00 0.32147E+00 0.50000E+00 Dt max 0.16105E-03 -0.99999E+00 0.63119E-02 0.50000E+00 --------------------------------------------------------------- TIME FOR THE TIME STEP 5: 0.13633E+05 =============================================================== =============================================================== FINAL STAGE OF THE CALCULATION ============================== =========================================================== ** TIME FOR FINAL WRITING: 0.13430E+00 ----------------------- =========================================================== END OF CALCULATION ================== =============================================================== Destroying structures and ending computation Calculation time summary: User CPU time: 27816.447 s System CPU time: 3391.300 s Total CPU time: 499327.658 s Elapsed time: 31212.496 s CPU / elapsed time 1.000 Memory use summary: Total memory used: 32.532 GiB local minimum: 1.849 GiB (rank 10) local maximum: 2.085 GiB (rank 12) Theoretical instrumented dynamic memory: 10.821 GiB local minimum: 622.703 MiB (rank 10) local maximum: 711.472 MiB (rank 15)