Parallel Code_Saturne with partitioning in 8 sub-domains Code_Saturne is running *********************** Working directory (to be periodically cleaned) : /home/alexandreguilloux/tmp_Saturne/110927_Etuve_reel.05-RADIATION.10191229 Kernel version: /opt/code_saturne-2.0-rc2/cs-2.0-rc2 Preprocessor: /opt/code_saturne-2.0-rc2/cs-2.0-rc2/bin *************************************************************** Compilation of user subroutines and linking of Code_Saturne *************************************************************** Use the radiative transfer module ******************************************** Preparing calculation ******************************************** ******************************************** Starting calculation ******************************************** Warning: /opt/code_saturne-2.0-rc2/cs-2.0-rc2/bin/cs_partition not found. The partitioner may not have been installed (this is the case if neither METIS nor. SCOTCH are avaialable). Unoptimized partitioning will be used, so parallel performance may be degraded. command: /home/alexandreguilloux/tmp_Saturne/110927_Etuve_reel.05-RADIATION.10191229/cs_solver --mpi --log 0 --param 111018_Etuve_reel_comparaison_exp_Radiation.xml *************************************************************** (R) Code_Saturne Version 2.0.0-rc2 Copyright (C) 1998-2009 EDF S.A., France build mer. 18 août 2010 15:46:16 CEST MPI version 2.1 (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: mer. 19 oct. 2011 12:29:49 CEST System: Linux 2.6.32-24-generic Machine: TSL-ORDI7 Processor: AMD Opteron(tm) Processor 6128 Memory: 12329428 User: alexandreguilloux (Alexandre Guilloux) Directory: /home/alexandreguilloux/tmp_Saturne/110927_Etuve_reel.05-RADIATION.10191229 MPI ranks: 8 =============================================================== CALCULATION PREPARATION ======================= =========================================================== Reading file: preprocessor_output No error detected during the data verification (interface, usini1 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Geometry NPRFML = 1 (Nb max. family properties ) NFML = 18 (Nb families ) --- Physics NPHAS = 1 (Nb phases ) NVAR = 7 (Nb variables ) NSCAL = 1 (Nb scalars ) NSCAUS = 1 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) NPROCE = 15 (Nb cell properties ) NPROFA = 1 (Nb internal face properties ) NPROFB = 10 (Nb boundary face properties ) ----------------------------------------------------------- ** PHYSICAL PROPERTIES ------------------- GX = 0.00000E+00 (Gravity x component ) GY = 0.00000E+00 (Gravity y component ) GZ = -0.98100E+01 (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 ) -- Phase: 1 RO0 = 0.12930E+01 (Reference density ) VISCL0 = 0.30000E-04 (Ref. molecular dyn. visc. ) CP0 = 0.10600E+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.00000E+00 (Ref. temperature ) IROVAR = 1 (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 = 60 (Turbulence model ) IDEUCH = 1 (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 ) - k-omega SST (ITURB = 60) ALMAX = -0.10000E+14 (Characteristic length ) UREF = 0.10000E+01 (Characteristic velocity ) IKECOU = 0 (k-epsilon coupling mode ) IGRAKE = 1 (Account for gravity ) --- 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)) - k-omega SST (ITURB = 60) CKWSK1 = 0.11765E+01 (sigma_k1 constant ) CKWSK2 = 0.20000E+01 (sigma_k2 constant ) CKWSW1 = 0.20000E+01 (sigma_omega1 constant ) CKWSW2 = 0.11682E+01 (sigma_omega2 constant ) CKWBT1 = 0.75000E-01 (beta1 constant ) CKWBT2 = 0.82800E-01 (beta2 constant ) CKWGM1 = 0.53933E+00 (gamma1 constant ) CKWGM2 = 0.41667E+00 (gamma2 constant ) CKWA1 = 0.31000E+00 (a1 constant to compute mu_t ) CKWC1 = 0.10000E+02 (c1 const. for prod. limiter ) CMU = ----------------------------------------------------------- ** SECONDARY VISCOSITY ------------------- --- Phase: 1 IVISSE = 1 (1: accounted for ) ----------------------------------------------------------- ** RADIATIVE THERMAL TRANSFER -------------------------- --- Phase: 1 IIRAYO = 1 (0: no; 1: DOM; 2: P-1 ) ICSALT = 1 (Assoc. thermal scalar num. ) ISCSTH = 1 (-1: T(C); 1: T(K); 2: H ) --- Options: ISUIRD = 0 (0: no restart; 1: restart ) NFREQR = 1 (Radiation pass frequency ) NDIREC = 128 (32 or 128 directions(if DOM)) IDIVER = 2 (0 1 or 2: compute radiat. ST) IMODAK = 0 (1: modak absor coef; 0 else ) IIMPAR = 1 (0 1 or 2: print wall temp. ) IIMLUM = 0 (0 1 or 2: print solver info ) --- Graphical output: NBRVAF = Wall_temp NBRVAF = Incident_flux NBRVAF = Th_conductivity NBRVAF = Thickness NBRVAF = Emissivity NBRVAF = Net_flux NBRVAF = Flux_convectif NBRVAF = Convective_exch_coef ----------------------------------------------------------- ** TIME STEPPING ------------- STEADY ALGORITHM --- Global parameters IDTVAR = -1 (-1: steady algorithm ) RELXST = 0.60000E+00 (Reference relaxation coeff. ) --- Frozen velocity field ICCVFG = 0 (1: frozen velocity field ) --- Per variable relaxation coefficient --------------------- RELAXV variable --------------------- Pression 0.1000E+00 VitesseX 0.6000E+00 VitesseY 0.6000E+00 VitesseZ 0.6000E+00 EnerTurb 0.7000E+00 omega 0.7000E+00 TempK 0.6000E+00 ---------------------------- RELAXV = [0.,1.] (relaxation coefficient ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- ------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV ------------------------------------------------------------- Pression 0 1 1 1 0 0.00E+00 0.10E+01 VitesseX 1 1 1 1 0 0.00E+00 0.10E+01 VitesseY 1 1 1 1 0 0.00E+00 0.10E+01 VitesseZ 1 1 1 1 0 0.00E+00 0.10E+01 EnerTurb 1 1 1 1 0 0.00E+00 0.10E+01 omega 1 1 1 1 0 0.00E+00 0.10E+01 TempK 1 1 1 1 0 0.00E+00 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 ----------------------------------------------------------- Pression 100 2 0.1000E-04 0.1000E-07 0.0000E+00 VitesseX 100 1 0.1000E-04 0.1000E-07 0.0000E+00 VitesseY 100 1 0.1000E-04 0.1000E-07 0.0000E+00 VitesseZ 100 1 0.1000E-04 0.1000E-07 0.0000E+00 EnerTurb 100 1 0.1000E-04 0.1000E-07 0.0000E+00 omega 100 1 0.1000E-04 0.1000E-07 0.0000E+00 TempK 100 1 0.1000E-04 0.1000E-07 0.0000E+00 ----------------------------------------------------------- ----------------------------------- Variable IRCFLU IMLIGR CLIMGR ----------------------------------- Pression 1 1 0.1500E+01 VitesseX 1 1 0.1500E+01 VitesseY 1 1 0.1500E+01 VitesseZ 1 1 0.1500E+01 EnerTurb 1 1 0.1500E+01 omega 1 1 0.1500E+01 TempK 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 ) ----------------------------------------------------------- ** WALL DISTANCE COMPUTATION ------------------------- ICDPAR = 1 ( 1: std, reread if restart (-1: std, recomputed if restrt ( 2: old, reread if restart (-2: old, recomputed if restrt NITMAY = 10000 (Nb iter for iter resolution ) NSWRSY = 1 (Nb iter for rhs reconstr. ) NSWRGY = 100 (Nb iter for grad. reconstr. ) IMLIGY = 1 (Gradient limitation method ) IRCFLY = 1 (Conv. Diff. flow reconstr. ) ISCHCY = 1 (Convective scheme ) ISSTPY = 0 (Slope tet use ) IMGRPY = 0 (Multigrid algorithm ) IWARNY = 0 (Verbosity level ) NTCMXY = 1000 (Nb iter for steady convect. ) BLENCY = 0.00000E+00 (2nd order conv. scheme prop.) EPSILY = 0.10000E-07 (Iterative solver precision ) EPSRSY = 0.00000E+00 (rhs reconstruction precision) EPSRGY = 0.10000E-04 (Gradient reconstr. precision) CLIMGY = 0.15000E+01 (Coeff. for grad. limitation ) EXTRAY = 0.00000E+00 (Coeff. for grad. extrapolat.) COUMXY = 0.50000E+04 (Max CFL for convection ) EPSCVY = 0.10000E-07 (Precision for steady conv. ) YPLMXY = 0.20000E+03 (y+ max w. damping influence ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------------ Variable IRESOL NITMAX EPSILO IDIRCL ------------------------------------------ Pression -1 10000 0.1000E-04 1 VitesseX -1 10000 0.1000E-04 1 VitesseY -1 10000 0.1000E-04 1 VitesseZ -1 10000 0.1000E-04 1 EnerTurb -1 10000 0.1000E-04 1 omega -1 10000 0.1000E-04 1 TempK -1 10000 0.1000E-04 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 ------------------------------ Pression 1 100 10 VitesseX 0 100 10 VitesseY 0 100 10 VitesseZ 0 100 10 EnerTurb 0 100 10 omega 0 100 10 TempK 0 100 10 ------------------------------ IMGR = 0 ou 1 (1: multigrid activated ) NCYMAX = (Max number of cycles ) NITMGF = (Max nb iter on coarsest grid) ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) ------------------------------------------------------------- Variable Number IPHSCA ISCSTH IVISLS VISLS0 SIGMAS ------------------------------------------------------------- TempK 1 1 1 0 0.4312E-04 0.1000E+01 ------------------------------------------------------------- ------------------------------------------ Variable Number IPHSCA ISCAVR RVARFL ------------------------------------------ TempK 1 1 0 0.8000E+00 ------------------------------------------ ------------------------------------------------------ Variable Number IPHSCA ICLVFL SCAMIN SCAMAX ------------------------------------------------------ TempK 1 1 -1 0.0000E+00 0.8000E+03 ------------------------------------------------------ ------------------------------------------------------------- For each scalar, the number indicates it's rank in the list of all scalars. User scalars are placed first, from 1 to NSCAUS. Specific physics scalars are placed at the end, from NSCAUS+1 to NSCAPP+NSCAUS=NSCAL. IPHSCA = (Carrier phase ) ISCSTH = -1,0, 1 ou 2 (T (C), Passive, T (K) or H ) IVISLS = 0 ou >0 (Viscosity: constant or not ) VISLS0 = >0 (Reference viscosity ) SIGMAS = >0 (Schmidt ) ISCAVR = 0 ou >0 (Associat. scalar if variance) RVARFL = >0 (Rf, cf variance dissipation ) ICLVFL = 0, 1 ou 2 (Variance clipping mode ) SCAMIN = (Min authorized value ) SCAMAX = (Max authorized value ) For variances, SCAMIN is ignored and SCAMAX is used only if ICLVFL = 2 ------------------------------------------------------ Scalar THETSS IVSEXT THETVS ------------------------------------------------------ 1 0.0000E+00 0 0.0000E+00 ------------------------------------------------------ THETSS = (theta for source terms ) ((1+theta).new-theta.old ) IVSEXT = (extrap. total viscosity ) (0: explicit ) (1: n+thetvs with thetvs=1/2 ) (2: n+thetvs with thetvs=1 ) THETVS = (theta for scalar diffusivity ((1+theta).new-theta.old ) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 0 (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 = 50000 (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 = 1000 (Output frequency ) Number Name 2 Pression 3 VitesseX 4 VitesseY 5 VitesseZ 6 EnerTurb 7 omega 8 TempK 9 MasseVol 11 ViscTurb 14 total_pressure -- -- --- Probe history files NTHIST = 1 (Output frequency ) NCAPT = 15 (Number of probes ) NTHSAV = 0 (Checkpoint frequency ) Number Name Nb. probes (-1: all) 2 Pression -1 3 VitesseX -1 4 VitesseY -1 5 VitesseZ -1 6 EnerTurb -1 7 omega -1 8 TempK -1 9 MasseVol -1 10 ViscLam -1 11 ViscTurb -1 14 total_pressure -1 15 Lumin -1 16 Qxrad -1 17 Qyrad -1 18 Qzrad -1 19 Srad1 -1 21 Absorp1 -1 22 Emiss1 -1 23 CoefAb_1 -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Pression 0 3 VitesseX 0 4 VitesseY 0 5 VitesseZ 0 6 EnerTurb 0 7 omega 0 8 TempK 0 9 MasseVol -999 -- -- -- --- Additional post-processing variables IPSTDV = 4 (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 No "domain_number_8" file available; domain partitioning will use a space-filling curve. ---------------------------------------------------------- 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 Interface creation: 0.0181 s Extended connectivity creation: 0.846 s Halo creation: 0.197 s Total time for halo creation: 1.06 s ---------------------------------------------------------- Histogram of the number of cells per rank: minimum value = 183016 maximum value = 183017 1 : [ 183016 ; 183017 ] = 8 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 200216 maximum value = 210651 1 : [ 200216 ; 202303 [ = 2 2 : [ 202303 ; 204390 [ = 1 3 : [ 204390 ; 206477 [ = 1 4 : [ 206477 ; 208564 [ = 3 5 : [ 208564 ; 210651 ] = 1 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 17200 maximum value = 27634 1 : [ 17200 ; 19286 [ = 2 2 : [ 19286 ; 21373 [ = 1 3 : [ 21373 ; 23460 [ = 1 4 : [ 23460 ; 25547 [ = 3 5 : [ 25547 ; 27634 ] = 1 Histogram of the number of ghost cells in the standard neighborhood per rank: minimum value = 7547 maximum value = 12100 1 : [ 7547 ; 8457 [ = 1 2 : [ 8457 ; 9368 [ = 1 3 : [ 9368 ; 10278 [ = 1 4 : [ 10278 ; 11189 [ = 4 5 : [ 11189 ; 12100 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 412332 maximum value = 415343 1 : [ 412332 ; 412934 [ = 1 2 : [ 412934 ; 413536 [ = 1 3 : [ 413536 ; 414138 [ = 2 4 : [ 414138 ; 414740 [ = 1 5 : [ 414740 ; 415343 ] = 3 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 13544 maximum value = 15525 1 : [ 13544 ; 13940 [ = 2 2 : [ 13940 ; 14336 [ = 2 3 : [ 14336 ; 14732 [ = 1 4 : [ 14732 ; 15128 [ = 1 5 : [ 15128 ; 15525 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 1 maximum value = 5 1 : [ 1 ; 2 [ = 1 2 : [ 2 ; 3 [ = 380 3 : [ 3 ; 4 [ = 28612 4 : [ 4 ; 5 ] = 1435142 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 5 maximum value = 7 1 : [ 5 ; 6 [ = 1 2 : [ 6 ; 7 ] = 7 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families Renumbering mesh: Mesh Number of cells: 1464135 Number of interior faces: 3275896 Number of boundary faces: 115836 Number of vertices: 543500 Computing geometric quantities (0.219 s) Mesh coordinates: minimum and maximum X : -2.0800000e+00 7.2000000e-02 Y : -8.9500000e-01 8.9300000e-01 Z : -2.3116950e-11 2.1600000e+00 --- Main Fortran work arrays: longia = 6061260 (Number of integers) longra = 24245040 (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.193488578E-09 Maximum control volume = 0.144580590E-03 Total volume for the domain = 0.823937671E+01 Extended neighborhood reduced by non-orthogonality -------------------------------------------------- Size of complete cell-cell connectivity: 68234752 Size of filtered cell-cell conectivity: 2758263 65476489 connections removed, for a ratio of 4 % used --- Phase: 1 ALMAX = 0.20183E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. ALMAX is the length used to initialize the turbulence. ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ --------------------------------- Variable Min. value Max. value --------------------------------- Pression 0.0000E+00 0.0000E+00 VitesseX 0.0000E+00 0.0000E+00 VitesseY 0.0000E+00 0.0000E+00 VitesseZ 0.0000E+00 0.0000E+00 EnerTurb 0.6000E-03 0.6000E-03 omega 0.1214E-01 0.1214E-01 TempK 0.6581E+03 0.6581E+03 --------------------------------- ------------------------------------------------------------- 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.100000000E+01 TIME STEP NUMBER 1 ============================================================= --- Phase: 1 --------------------------------- Property Min. value Max. value --------------------------------- MasseVol 0.5363E+00 0.5363E+00 ViscLam 0.3000E-04 0.3000E-04 ViscTurb 0.2652E-01 0.2652E-01 --------------------------------- --- Diffusivity: --------------------------------------- Scalar Number Min. value Max. value --------------------------------------- TempK 1 0.4312E-04 0.4312E-04 --------------------------------------- *** Temperature de la face d entree : usclim Temperature TempExp = 662.14999999999998 Temperature TempEntree = 662.14999999999998 -------------Verification------------------- Temperature TempExp = 662.14999999999998 Temperature Temp01m = 0.0000000000000000 Temperature TempEntree = 662.14999999999998 ** INFORMATIONS SUR LE MODULE DE RAYONNEMENT ------------------------------------------ Initialisation de la temperature de paroi (TPAROI) avec le profil utilisateur (TINTP) et du flux incident aux parois (QINCID). ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- Phase : 1 ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 6322 Smooth wall 5 106270 Rough wall 6 0 Symmetry 4 0 Free outlet 3 3244 Undefined 1 0 ------------------------------------------------------------------------- PHASE 1 : Boundary faces with free inlet/outlet detected Update of reference point for total pressure XYZP0 = -0.87119E+00 0.89796E-01 0.19536E+01 ** BOUNDARY MASS FLOW INFORMATION ------------------------------ Phase : 1 --------------------------------------------------------------- Boundary type Code Nb faces Mass flow --------------------------------------------------------------- Inlet 2 6322 0.000000000E+00 Smooth wall 5 106270 0.000000000E+00 Rough wall 6 0 0.000000000E+00 Symmetry 4 0 0.000000000E+00 Free outlet 3 3244 0.000000000E+00 Undefined 1 0 0.000000000E+00 --------------------------------------------------------------- ** BOUNDARY CONDITIONS FOR SMOOTH WALLS --------------------------------------- ------------------------------------------------------------ Phase 1 Minimum Maximum ------------------------------------------------------------ Rel velocity at the wall uiptn : 0.00000E+00 0.00000E+00 Friction velocity uet : 0.00000E+00 0.00000E+00 Friction velocity uk : 0.00000E+00 0.13416E-01 Dimensionless distance yplus : 0.77118E-01 0.43240E+01 ------------------------------------------------------ Nb of reversal of the velocity at the wall : 0 Nb of faces within the viscous sub-layer : 93459 Total number of wall faces : 106270 ------------------------------------------------------------ ** WALL DISTANCE ------------- Min distance = 0.30857E-03 Max distance = 0.61514E+00 (Distance calculation done in 737 iterations) ** INFORMATIONS SUR LE TERME SOURCE RADIATIF ----------------------------------------- Bi-CGstab [Rayon002]: divergence after 32 iterations: initial residual: 2.2123e+02; current residual: 1.1205e+07 Writing of "chr" (EnSight Gold) summary: CPU time for meshes: 0.130 CPU time for variables: 0.040 Elapsed time for meshes: 0.130 Elapsed time for variables: 0.037 Writing of "error" (EnSight Gold) summary: CPU time for meshes: 0.140 CPU time for variables: 0.170 Elapsed time for meshes: 0.142 Elapsed time for variables: 0.182 -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- /opt/code_saturne-2.0-rc2/ncs-2.0.0-rc2/src/base/cs_sles.c:1746: Fatal error. Bi-CGstab: error (divergence) solving for Rayon002 Call stack: 1: 0x7fd25983d8e7 (libsaturne.so.0) 2: 0x7fd2599af140 (libsaturne.so.0) 3: 0x7fd2598645b3 (libsaturne.so.0) 4: 0x40dd04 (cs_solver) 5: 0x7fd259b8181c (libsaturne.so.0) 6: 0x7fd259a2525c (libsaturne.so.0) 7: 0x7fd259855ac5 (libsaturne.so.0) 8: 0x7fd259830c03 (libsaturne.so.0) 9: 0x7fd259830ee5 (libsaturne.so.0) 10: 0x7fd256acec4d <__libc_start_main+0xfd> (libc.so.6) 11: 0x40ccf9 <> (cs_solver) End of stack -------------------------------------------------------------------------- mpiexec has exited due to process rank 0 with PID 15919 on node TSL-ORDI7 exiting without calling "finalize". This may have caused other processes in the application to be terminated by signals sent by mpiexec (as reported here). -------------------------------------------------------------------------- Error running the calculation. Check Kernel log (listing) and error* files for details ******************************************** Error in calculation stage. ********************************************