Parallel Code_Saturne with partitioning in 2 sub-domains Code_Saturne is running *********************** Working directory (to be periodically cleaned) : /home/caelinux/tmp_Saturne/100708.EADSIW_1ERS_CALCULS.07081106 Kernel version: /opt/code_saturne/Noyau/ncs-1.4.b Preprocessor version: /opt/code_saturne/Enveloppe/ecs-1.4.a *************************************************************** Compilation of user subroutines and linking of Code_Saturne *************************************************************** /home/caelinux/tmp_Saturne/100708.EADSIW_1ERS_CALCULS.07081106/compil.log *** "make" in a user directory ******************************************** Preparing calculation ******************************************** ******************************************** Calculation starting ******************************************** command: /home/caelinux/tmp_Saturne/100708.EADSIW_1ERS_CALCULS.07081106/cs14.exe -mpi --log 0 -param 100708_EADS-IW_1ers_calculs_laminaire_noT.xml *************************************************************** (R) Code_Saturne Version 1.4.b Copyright (C) 1998-2008 EDF S.A., France build Sun 17 May 2009 01:01:07 AM CEST MPI version 2.0 (Open MPI) 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: Thu 08 Jul 2010 11:06:40 AM CEST System: Linux 2.6.24-24-generic Machine: tsl-ordi6 Processor: Intel(R) Core(TM)2 Duo CPU E8400 @ 3.00GHz Memory: 6118336 User: caelinux (caelinux) Directory: /home/caelinux/tmp_Saturne/100708.EADSIW_1ERS_CALCULS.07081106 Reading file: preprocessor_output *************************************************************** CALCULATION PREPARATION ======================= =========================================================== No error detected during the data verification (interface, usini1 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Geometry NPRFML = 2 (Nb max. family properties ) NFML = 8 (Nb families ) --- Physics NPHAS = 1 (Nb phases ) 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.98100E+01 (Gravity y component ) GZ = 0.00000E+00 (Gravity z component ) -- Phase: 1 RO0 = 0.12506E+01 (Reference density ) VISCL0 = 0.17500E-04 (Ref. molecular dyn. visc. ) CP0 = 0.10400E+04 (Ref. specific heat ) ICP = 0 (> 0: variable CP (usphyv) ) P0 = 0.10130E+06 (Ref. total pressure ) PRED0 = 0.00000E+00 (Ref. reduced pressure ) T0 = 0.00000E+00 (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 ---------- --- Phase: 1 - Commons ITURB = 0 (Turbulence model ) IDEUCH = 0 (0: one-scale model ) (1: two-scale model ) (2: invariant wall function ) 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 ------------------- --- Phase: 1 IVISSE = 1 (1: accounted for ) ----------------------------------------------------------- ** TIME STEPPING ------------- STEADY ALGORITHM --- Global parameters IDTVAR = -1 (-1: steady algorithm ) RELXST = 0.90000E+00 (Reference relaxation coeff. ) --- Frozen velocity field ICCVFG = 0 (1: frozen velocity field ) --- Per variable relaxation coefficient --------------------- RELAXV variable --------------------- Pressure 0.1000E+00 VelocitU 0.9000E+00 VelocitV 0.9000E+00 VelocitW 0.9000E+00 ---------------------------- RELAXV = [0.,1.] (relaxation coefficient ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- ------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV ------------------------------------------------------------- Pressure 0 1 1 1 0 0.00E+00 0.10E+01 VelocitU 1 1 1 1 0 0.10E+01 0.10E+01 VelocitV 1 1 1 1 0 0.10E+01 0.10E+01 VelocitW 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 ------ 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 ) -- Phase: 1 IREVMC = 0 (Velocity reconstruction mode) 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 = 3 (Reconstruction mode ) ANOMAX = 0.78540E+00 (Non-ortho angle: limit for ) (least squares ext. neighbors) ----------------------------------------------------------- Variable NSWRGR NSWRSM EPSRGR EPSRSM EXTRAG ----------------------------------------------------------- Pressure 100 2 0.1000E-04 0.1000E-07 0.0000E+00 VelocitU 100 1 0.1000E-04 0.1000E-07 0.0000E+00 VelocitV 100 1 0.1000E-04 0.1000E-07 0.0000E+00 VelocitW 100 1 0.1000E-04 0.1000E-07 0.0000E+00 ----------------------------------------------------------- ----------------------------------- Variable IRCFLU IMLIGR CLIMGR ----------------------------------- Pressure 1 1 0.1500E+01 VelocitU 1 1 0.1500E+01 VelocitV 1 1 0.1500E+01 VelocitW 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 -1 10000 0.1000E-02 1 VelocitU -1 10000 0.1000E-02 1 VelocitV -1 10000 0.1000E-02 1 VelocitW -1 10000 0.1000E-02 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 numeber of levels ) ------------------------------ Variable IMGR NCYMAX NITMGF ------------------------------ Pressure 0 100 10 VelocitU 0 100 10 VelocitV 0 100 10 VelocitW 0 100 10 ------------------------------ IMGR = 0 ou 1 (1: multigrid activated ) NCYMAX = (Max number of cycles ) NITMGF = (Max nb iter on coarsest grid) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 1 (1: restarted calculuation ) ILEAUX = 1 (1: also read suiamx ) IECAUX = 1 (1: also write suiavx ) suiamx and suiavx are the auxiliary restart files. --- 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 = 150 (Final time step required ) --- CPU time margin TMARUS = -0.10000E+01 (CPU time margin before stop ) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- EnSight, MED, or CGNS output NTCHR = 25 (Output frequency ) Number Name 2 Pressure 3 VelocitU 4 VelocitV 5 VelocitW 10 total_pressure -- -- --- Probe history files NTHIST = 5 (Output frequency ) NCAPT = 5 (Number of probes ) NTHSAV = 0 (Checkpoint frequency ) Number Name Nb. probes (-1: all) 2 Pressure -1 3 VelocitU -1 4 VelocitV -1 5 VelocitW -1 6 Density -1 7 Lam. vis -1 10 total_pressure -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Pressure 0 3 VelocitU 0 4 VelocitV 0 5 VelocitW 0 10 total_pressure -999 -- -- -- --- Additional post-processing variables IPSTDV = 210 (Product of the following values based on activation ) IPSTYP = 2 (Yplus on boundary ) IPSTCL = 3 (Variables on boundary ) IPSTFT = 5 (Thermal flow on boundary ) IPSTFO = 7 (Force exerted on boundary ) and IPSTDV = 1 (No additional output ) ----------------------------------------------------------- ** FILES ----- --- Standard files GEOMETRY MAIN RESTART AUX RESTART STOP Name geomet suiamo suiamx ficstp Unit 10 12 MAIN CHECKPOINT AUX CHECKPOINT Name suiava suiavx --- Radiative files DATA RESTART JANAF CHECKPOINT Name dp_tch rayamo JANAF rayava Unit 25 25 --- Lagrangian files --- RESTART --- -- CHECKPOINT - ---- POST ---- CALCUL. STAT. CALCUL. STAT. LOG HIST. Name lagamo lasamo lagava lasava listla Unit 80 81 82 --- Other files for Lagrangian module Unit 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 --- Developper files for probe history IMPHIS(1) IMPHIS(2) Unit 30 31 --- User files for probe history 16 files Name and unit ush001 ush002 ush003 ush004 ush005 ush006 ush007 ush008 ush009 ush010 ush011 ush012 ush013 ush014 ush015 ush016 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 --- Free user files 10 files Name and unit usrf01 usrf02 usrf03 usrf04 usrf05 usrf06 usrf07 usrf08 usrf09 usrf10 70 71 72 73 74 75 76 77 78 79 --- Preprocessor files Format: 1=Read preprocessor_output file 0=Standalone solver Format 1 ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated ) NALINF = -999 (Fluid initialization iterations) ----------------------------------------------------------- Finished reading: preprocessor_output Reading file: domain_number_2 Finished reading: domain_number_2 ---------------------------------------------------------- Halo construction with extended neighborhood ============================================ Interface creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Extended neighborhood structures definition Halo creation times summary Creating interface: 0.0304 s Creating extended connectivity: 5.91 s Creating halo: 0.157 s Total time for halo creation: 6.09 s ---------------------------------------------------------- Number of standard cells: 708145 Histogram of the number of cells per rank: minimum value = 708145 maximum value = 708145 1 : [ 708145 ; 708145 ] = 2 ---------------------------------------------------------- Number of cells + halo cells: 718675 Histogram of number of standard + halo cells (NCELET) per rank: minimum value = 718596 maximum value = 718675 1 : [ 718596 ; 718611 [ = 1 2 : [ 718611 ; 718627 [ = 0 3 : [ 718627 ; 718643 [ = 0 4 : [ 718643 ; 718659 [ = 0 5 : [ 718659 ; 718675 ] = 1 ---------------------------------------------------------- Local number of ghost cells: 10530 in the standard neighborhood: 3389 in the extended neighborhood: 7141 Histogram of the number of ghost cells per rank: minimum value = 10451 maximum value = 10530 1 : [ 10451 ; 10466 [ = 1 2 : [ 10466 ; 10482 [ = 0 3 : [ 10482 ; 10498 [ = 0 4 : [ 10498 ; 10514 [ = 0 5 : [ 10514 ; 10530 ] = 1 ---------------------------------------------------------- Number of neighboring domains: 1 Histogram of the number of neighboring domains per rank: minimum value = 1 maximum value = 1 1 : [ 1 ; 1 ] = 2 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: ----------------------------------------------------------- ** VECTORIZATION (1: POSSIBLE, 0: IMPOSSIBLE) ------------- Interior faces Boundary faces IVECTI 0 IVECTB 0 ----------------------------------------------------------- Mesh Number of cells: 1416290 Number of interior faces: 2753241 Number of boundary faces: 158678 Number of vertices: 274812 Computing geometric quantities (0.816 s) Mesh coordinates: minimum and maximum X : -1.3415047e-14 2.0000000e+00 Y : -8.1000000e-01 8.1000000e-01 Z : -6.9860562e-13 8.0937241e-01 --- Main Fortran work arrays: LONGIA = 21560250 (Number of integers) LONGRA = 86241000 (Number of reals) (4 bytes/integer, 8 bytes/real) --- User Fortran work arrays: NITUSE = 0 (Number of integers) NRTUSE = 0 (Number of reals) --- Information on the volumes Minimum control volume = 0.819915820E-11 Total volume for the domain = 0.173234468E+01 Extended neighborhood reduced by non-orthogonality -------------------------------------------------- Size of complete cell-cell connectivity: 88707068 Size of filtered cell-cell conectivity: 369230 88337838 cells removed, for a ratio of 0.42 % used --- Phase: 1 ALMAX = 0.12008E+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 Reading dimensions complete Reading the previous time step number (restarting computation) NTPABS = 50 Reading the previous time step number (restarting computation) TTPABS = 0.5000E+02 Reading options complete Reading the pressure, velocity, turbulent variables complete Reading scalars complete Reading complete READING THE AUXILIARY FILE SUITE Start reading Fin de la lecture des options @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNNING : WHEN READING THE AUXILIARY FILE SUITE @ ********* @ LECTURE DES PROPRIETES PHYSIQUES @ @ It was not possible to read some values from the @ auxiliary file suite. @ They will be initialised by the defaults values. @ @ This situation can occur when the file @ suite is issued from a run realised with different @ options or when the file is dammaged. @ @ The run will continue... @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ 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 End reading ----------------------------------------------------------- ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ --------------------------------- Variable Min. value Max. value --------------------------------- Pressure -0.2985E+01 0.2028E+02 VelocitU -0.5213E+01 0.2184E+01 VelocitV -0.6648E+01 0.3742E+01 VelocitW -0.3766E+01 0.3021E+01 --------------------------------- ------------------------------------------------------------- ALL PHASES : 1D-WALL THERMAL MODULE NOT ACTIVATED NFPT1D = 0 ------------------------------------------------------------- PHASE 1 : HEAD LOSS TREATMENT NOT ACTIVATED NCEPDC = 0 ------------------------------------------------------------- PHASE 1 : MASS SOURCE TERMS TREATMENT NOT ACTIVATED NCETSM = 0 ------------------------------------------------------------- *************************************************************** MAIN CALCULATION ================ *************************************************************** *************************************************************** INSTANT 0.510000000E+02 TIME STEP NUMBER 51 ============================================================= --- Phase: 1 --------------------------------- Property Min. value Max. value --------------------------------- Density 0.1251E+01 0.1251E+01 Lam. vis 0.1750E-04 0.1750E-04 --------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- Phase : 1 ------------------------------------------------------------------------- Boundary type Code Nb faces (tot/loc) Min index Max index ------------------------------------------------------------------------- Entree 2 99 99 1 99 Paroi lisse 5 127144 57712 18145 75856 Paroi rugueuse 6 0 0 75857 75856 Symetrie 4 31346 18045 100 18144 Sortie libre 3 89 0 100 99 Indefini 1 0 0 1 0 ------------------------------------------------------------------------- PHASE 1 : Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = 0.00000E+00 0.99337E-01 0.10297E-01 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 -0.226323723E-17 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 0.367186179E+00 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.29230E+01 Friction velocity uet : 0.67888E-04 0.36574E+00 Friction velocity uk : 0.00000E+00 0.32274E+00 Dimensionless distance yplus : 0.30272E-01 0.13079E+03 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 34776 Total number of wall faces : 127144 ------------------------------------------------------------ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING: MESH NOT ENOUGH REFINED AT THE WALL @ ******** @ PHASE 1 @ The mesh does not seem to be enough refined at the wall @ to be able to run a laminar simulation. @ @ The last time step at which too large values for the @ dimensionless distance to the wall (yplus) have been @ observed is the time step 51 @ @ The minimum value for yplus must be lower than the @ limit value YPLULI = 0.23810E+01 @ @ Have a look at the distribution of yplus at the wall @ (with EnSight for example) to conclude on the way @ the results quality might be affected. @ @ This warning is only printed at the first two @ occurences of the problem and at the last time step @ of the calculation. The vanishing of the message does @ not necessarily mean the vanishing of the problem. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.82018E+02 1297 0.98407E-03 0.18892E+01 c VelocitU 0.24802E+00 93 0.84186E-03 0.55459E+07 c VelocitV 0.19292E+00 97 0.97069E-03 0.48254E+07 c VelocitW 0.19610E+00 93 0.94136E-03 0.35619E+07 --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.22228E+04 0.50462E+04 -- -- v VelocitU -0.10722E+05 0.12207E+05 -- -- v VelocitV -0.15528E+05 0.14838E+05 -- -- v VelocitW -0.12566E+05 0.11524E+05 -- -- v total_pressu 0.99078E+05 0.10635E+06 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 51: 0.18181E+03 *************************************************************** INSTANT 0.520000000E+02 TIME STEP NUMBER 52 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 -0.282539664E-17 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 0.158146599E-01 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.75572E+04 Friction velocity uet : 0.68908E-02 0.30439E+03 Friction velocity uk : 0.00000E+00 0.30439E+03 Dimensionless distance yplus : 0.14603E+01 0.77561E+05 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 28 Total number of wall faces : 127144 ------------------------------------------------------------ @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING: MESH NOT ENOUGH REFINED AT THE WALL @ ******** @ PHASE 1 @ The mesh does not seem to be enough refined at the wall @ to be able to run a laminar simulation. @ @ The last time step at which too large values for the @ dimensionless distance to the wall (yplus) have been @ observed is the time step 52 @ @ The minimum value for yplus must be lower than the @ limit value YPLULI = 0.23810E+01 @ @ Have a look at the distribution of yplus at the wall @ (with EnSight for example) to conclude on the way @ the results quality might be affected. @ @ This warning is only printed at the first two @ occurences of the problem and at the last time step @ of the calculation. The vanishing of the message does @ not necessarily mean the vanishing of the problem. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.31438E+08 47 0.49659E-03 0.16541E+01 c VelocitU 0.20616E+05 33 0.79002E-03 0.93062E+09 c VelocitV 0.15305E+05 32 0.88425E-03 0.85812E+09 c VelocitW 0.19919E+05 32 0.81585E-03 0.79083E+09 --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.31005E+07 0.12606E+07 -- -- v VelocitU -0.73294E+04 0.94961E+08 -- -- v VelocitV -0.10431E+08 0.12439E+08 -- -- v VelocitW -0.58350E+08 0.33906E+05 -- -- v total_pressu-0.29992E+07 0.13619E+07 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 52: 0.38262E+02 *************************************************************** INSTANT 0.530000000E+02 TIME STEP NUMBER 53 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 -0.546352508E-16 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 -0.648023418E+04 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.13658E+05 Friction velocity uet : 0.24158E-01 0.44646E+03 Friction velocity uk : 0.00000E+00 0.44646E+03 Dimensionless distance yplus : 0.29407E+01 0.19196E+06 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 1 Total number of wall faces : 127144 ------------------------------------------------------------ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.20467E+05 2490 0.97517E-03 0.97130E+00 c VelocitU 0.11698E+12 20 0.71246E-03 0.16860E+20 c VelocitV 0.17129E+11 21 0.70787E-03 0.32802E+18 c VelocitW 0.15385E+11 18 0.77020E-03 0.93179E+18 --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.43872E+09 0.53927E+11 -- -- v VelocitU -0.28713E+08 0.23849E+08 -- -- v VelocitV -0.70714E+07 0.53133E+08 -- -- v VelocitW -0.11380E+08 0.77771E+07 -- -- v total_pressu-0.43862E+09 0.53927E+11 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 53: 0.28026E+03 *************************************************************** INSTANT 0.540000000E+02 TIME STEP NUMBER 54 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 -0.244072229E-14 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 0.878800148E+01 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.32469E+00 0.10441E+08 Friction velocity uet : 0.64417E-01 0.25402E+06 Friction velocity uk : 0.64417E-01 0.25402E+06 Dimensionless distance yplus : 0.41912E+01 0.39784E+08 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 0 Total number of wall faces : 127144 ------------------------------------------------------------ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.54235E+13 34 0.68376E-03 0.11058E+02 c VelocitU 0.39677E+12 21 0.89836E-03 0.13903E+20 c VelocitV 0.70303E+12 20 0.99767E-03 0.12901E+20 c VelocitW 0.45545E+11 24 0.96963E-03 0.15459E+19 --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.32462E+14 0.26841E+13 -- -- v VelocitU -0.51800E+12 0.34489E+11 -- -- v VelocitV -0.36862E+10 0.18535E+11 -- -- v VelocitW -0.85249E+09 0.34780E+11 -- -- v total_pressu-0.32462E+14 0.26841E+13 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 54: 0.33990E+02 *************************************************************** INSTANT 0.550000000E+02 TIME STEP NUMBER 55 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 -0.138698681E-11 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 0.402928835E+09 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.40077E+00 0.33383E+11 Friction velocity uet : 0.68816E-01 0.56593E+09 Friction velocity uk : 0.68816E-01 0.56593E+09 Dimensionless distance yplus : 0.58242E+01 0.29002E+11 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 0 Total number of wall faces : 127144 ------------------------------------------------------------ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure 0.13684E+17 252 0.83150E-03 0.14195E+01 c VelocitU 0.23172E+18 81 0.94706E-03 0.76793E+36 c VelocitV 0.40487E+17 91 0.91961E-03 0.16449E+37 c VelocitW 0.83817E+17 85 0.91945E-03 0.48192E+35 --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure -0.12174E+23 0.33538E+23 -- -- v VelocitU -0.55627E+19 0.12060E+25 -- -- v VelocitV -0.64242E+18 0.32005E+24 -- -- v VelocitW -0.10319E+24 0.40456E+19 -- -- v total_pressu-0.12174E+23 0.33538E+23 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 55: 0.73657E+02 *************************************************************** INSTANT 0.560000000E+02 TIME STEP NUMBER 56 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 -0.367269666E+00 Paroi lisse 5 127144 0.783178594E-11 Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 0.000000000E+00 Sortie libre 3 89 -0.753734022E+13 Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.12480E+06 0.15846E+19 Friction velocity uet : 0.39161E+04 0.15060E+17 Friction velocity uk : 0.39161E+04 0.15060E+17 Dimensionless distance yplus : 0.32803E+06 0.78524E+19 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 0 Total number of wall faces : 127144 ------------------------------------------------------------ ** INFORMATIONS ON THE CONVERGENCE ------------------------------- --------------------------------------------------------------- Variable Rhs norm N_iter Norm. residual derive --------------------------------------------------------------- c Pressure NaN 2 0.00000E+00 NaN c VelocitU 0.13519E+37 1 NaN NaN c VelocitV 0.35875E+36 1 NaN NaN c VelocitW 0.11566E+36 1 NaN NaN --------------------------------------------------------------- ** INFORMATIONS ON THE VARIABLES ----------------------------- --------------------------------------------------------------- Variable Min. value Max. value Min clip Max clip --------------------------------------------------------------- v Pressure 0.10000E+13 -0.10000E+13 -- -- v VelocitU 0.10000E+13 -0.10000E+13 -- -- v VelocitV 0.10000E+13 -0.10000E+13 -- -- v VelocitW 0.10000E+13 -0.10000E+13 -- -- v total_pressu 0.10000E+13 -0.10000E+13 0 0 --------------------------------------------------------------- ** INFORMATIONS ON THE CLIPPINGS ----------------------------- --------------------------------------------------------------- Variable Min wo clips Max wo clips Min clip Max clip --------------------------------------------------------------- --------------------------------------------------------------- CPU TIME FOR THE TIME STEP 56: 0.26962E+02 *************************************************************** INSTANT 0.570000000E+02 TIME STEP NUMBER 57 ============================================================= ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Entree 2 99 NaN Paroi lisse 5 127144 NaN Paroi rugueuse 6 0 0.000000000E+00 Symetrie 4 31346 NaN Sortie libre 3 89 NaN Indefini 1 0 0.000000000E+00 --------------------------------------------------------------- CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 100 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 101 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 102 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 103 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 104 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 105 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 106 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 107 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 108 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 109 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 110 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 111 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 112 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 113 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 114 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 115 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 116 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 117 ............. ............. MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 75854 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 75855 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 CAUSTA SUBROUTINE CALLED FOR FACE IFAC= 75856 MAXIMUM NUMBER OF ITERATIONS REACHED FOR THE COMPUTATION OF UET: NITM= 100, DESIRED PRECISION EPS= 0.100E-02 ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : NaN NaN Friction velocity uk : 0.00000E+00 0.00000E+00 Dimensionless distance yplus : NaN NaN ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 127144 Total number of wall faces : 127144 ------------------------------------------------------------ mpirun: killing job... mpirun noticed that job rank 0 with PID 28596 on node tsl-ordi6 exited on signal 15 (Terminated). 1 additional process aborted (not shown) Error running the calculation. Check Kernel log (listing) and error* files for details ******************************************** Error in calculation stage. ********************************************