command: ./cs_solver *************************************************************** (R) Code_Saturne Version 4.0.2 Copyright (C) 1998-2015 EDF S.A., France revision 4.0.2 build Thu 27 Aug 2015 12:32:23 BST MPI version 3.0 (Open MPI 1.8.4) 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 01 Sep 2015 13:46:07 BST System: Linux 3.16.0-46-generic Machine: Marputor Processor: model name : Intel(R) Core(TM) i7-5500U CPU @ 2.40GHz Memory: 15960 MB User: martian (Martian) Directory: /home/martian/Practice/CS/utube/serial/RESU/20150901-1344 I/O read method: standard input and output, serial access I/O write method: standard input and output, serial access External libraries for partitioning: ParMETIS 4.0.3 Reading metadata from file: "mesh_input" =============================================================== CALCULATION PREPARATION ======================= =========================================================== @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ ATTENTION : ENTREE DES DONNEES @ ========= @ REFERENCE VELOCITY UREF WAS NOT DEFINED @ or is ill defined (NEGATIVE value ? ). @ It IS EQUAL -0.10000E+14 @ @ calculation can only run if turbulence is defined @ via a restart file or the interface or the @ cs_user_initialization subroutine. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 6 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) NPROCE = 6 (Nb cell properties ) ----------------------------------------------------------- ** PHYSICAL PROPERTIES ------------------- GX = 0.00000E+00 (Gravity x component ) GY = 0.00000E+00 (Gravity y component ) GZ = 0.00000E+00 (Gravity z component ) ICORIO = 0 (Coriolis source terms ) -- Continuous phase: RO0 = 0.42827E+02 (Reference density ) VISCL0 = 0.31816E-04 (Ref. molecular dyn. visc. ) CP0 = 0.26863E+04 (Ref. specific heat ) ICP = 0 (> 0: variable CP (usphyv) ) P0 = 0.15000E+08 (Ref. total pressure ) PRED0 = 0.00000E+00 (Ref. reduced pressure ) T0 = 0.83315E+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 ----------------------------------------------------------- ** HOMOGENEOUS MIXTURE MODEL FOR CAVITATION ---------------------------------------- ICAVIT = -1 (-1: single phase flow ) ( 0: no vap./cond. model ) ( 1: Merkle's model ) ----------------------------------------------------------- ** THERMAL MODEL ------------- --- Continuous phase: - Commons ITHERM = 0 (0: no thermal model ) (1: temperature ) (2: enthalpy ) (3: total energy ) ITPSCL = 0 (0: none ) (1: temperature in Kelvin ) (2: temperature in Celsius ) ISCALT = -1 (Thermal scalar number ) ** TURBULENCE ---------- --- Continuous phase: - Commons ITURB = 60 (Turbulence model ) IDEUCH = 1 (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) ) - k-omega SST (ITURB = 60) ALMAX = -0.99900E+03 (Characteristic length ) UREF = -0.10000E+14 (Characteristic velocity ) IKECOU = 0 (k-epsilon coupling mode ) IGRAKE = 1 (Account for gravity ) - Rotation/curvature correction IRCCOR = 0 (0: desactivated ) (1: activated ) --- 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.10000E+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 = 0.90000E-01 (Cmu (or Beta*) constant for ) omega/epsilon conversion) ----------------------------------------------------------- ** SECONDARY VISCOSITY ------------------- --- Continuous phase: IVISSE = 1 (1: accounted for ) ----------------------------------------------------------- ** TIME STEPPING ------------- STEADY ALGORITHM --- Global parameters IDTVAR = -1 (-1: steady algorithm ) RELXST = 0.70000E+00 (Reference relaxation coeff. ) --- Frozen velocity field ICCVFG = 0 (1: frozen velocity field ) --- Per variable relaxation coefficient ----------------------------- Variable RELAXV ----------------------------- Velocity 0.7000E+00 Pressure 0.3000E+00 Turb Kinetic Ene 0.7000E+00 Omega 0.7000E+00 ---------------------------- RELAXV = [0.,1.] (relaxation coefficient ) ----------------------------------------------------------- ** CONVECTION - DIFFUSION ---------------------- --------------------------------------------------------------------- Variable ICONV IDIFF IDIFFT ISCHCV ISSTPC BLENCV THETAV --------------------------------------------------------------------- Velocity 1 1 1 1 0 0.10E+01 0.10E+01 Pressure 0 1 1 1 0 0.00E+00 0.10E+01 Turb Kinetic Ene 1 1 1 1 0 0.00E+00 0.10E+01 Omega 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 ------ IDILAT = 1 (1 : without unsteady term in the continuity equation 2 : with unsteady term in the continuity equation) IPOROS = 0 (0 : without porous media 1 : with porous media) 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) 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 ------------------------------------------------------------------- Velocity 100 1 0.1000E-04 0.1000E-06 0.0000E+00 Pressure 100 2 0.1000E-04 0.1000E-06 0.0000E+00 Turb Kinetic Ene 100 1 0.1000E-04 0.1000E-06 0.0000E+00 Omega 100 1 0.1000E-04 0.1000E-06 0.0000E+00 ----------------------------------------------------------- ------------------------------------------- Variable IRCFLU IMLIGR CLIMGR ------------------------------------------- Velocity 1 -1 0.1500E+01 Pressure 1 -1 0.1500E+01 Turb Kinetic Ene 1 -1 0.1500E+01 Omega 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 ) ----------------------------------------------------------- ** WALL DISTANCE COMPUTATION ------------------------- ICDPAR = 1 ( 1: std, reread if restart (-1: std, recomputed if restrt ( 2: old, reread if restart (-2: old, recomputed if restrt NITMAY = 10000 (Nb iter for iter resolution ) NSWRSY = 1 (Nb iter for rhs reconstr. ) NSWRGY = 100 (Nb iter for grad. reconstr. ) IMLIGY = -1 (Gradient limitation method ) IRCFLY = 1 (Conv. Diff. flow reconstr. ) ISCHCY = 1 (Convective scheme ) ISSTPY = 0 (Slope tet use ) IWARNY = 0 (Verbosity level ) NTCMXY = 1000 (Nb iter for steady convect. ) BLENCY = 0.00000E+00 (2nd order conv. scheme prop.) EPSILY = 0.10000E-07 (Iterative solver precision ) EPSRSY = 0.10000E-04 (rhs reconstruction precision) EPSRGY = 0.10000E-04 (Gradient reconstr. precision) CLIMGY = 0.15000E+01 (Coeff. for grad. limitation ) EXTRAY = 0.00000E+00 (Coeff. for grad. extrapolat.) COUMXY = 0.50000E+04 (Max CFL for convection ) EPSCVY = 0.10000E-07 (Precision for steady conv. ) YPLMXY = 0.20000E+03 (y+ max w. damping influence ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-07 1 Pressure 0.1000E-07 1 Turb Kinetic Ene 0.1000E-07 1 Omega 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 ) ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 0 (1: restarted calculation ) ILEAUX = 0 (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 = 400 (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 Velocity Pressure Turb Kinetic Energy Omega Laminar Viscosity Turb Viscosity Total Pressure -- -- --- 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 Velocity[X] -1 3 Velocity[Y] -1 4 Velocity[Z] -1 5 Pressure -1 6 Turb Kinetic Ene -1 7 Omega -1 10 Turb Viscosity -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Velocity 0 5 Pressure 0 6 Turb Kinetic Ene 0 7 Omega 0 9 Laminar Viscosit -999 10 Turb Viscosity -999 13 Total Pressure -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 0 (Force exerted by the fluid on the boundary) ipstft = 0 (Thermal flux at boundary) ipsttb = 0 (Temperature at boundary) ipstnu = 0 (Dimensionless thermal flux at boundary) ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated ) NALINF = 0 (Fluid initialization iterations) IFLXMW = 0 (ALE mass flux computation 0: thanks to vertices 1: thanks to mesh velocity) ----------------------------------------------------------- Postprocessing output writers: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: time dependency: fixed mesh output: at calculation end Reading file: mesh_input Finished reading: mesh_input Number of cells: 427500 ---------------------------------------------------------- Number of cells + halo cells: 427500 ---------------------------------------------------------- Number of interior faces: 1269975 ---------------------------------------------------------- Number of boundary faces: 25050 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 4 maximum value = 6 1 : [ 4 ; 5 [ = 120 2 : [ 5 ; 6 ] = 427380 ---------------------------------------------------------- Renumbering mesh: renumbering for cells: pre-numbering: none cells adjacent to ghost cells last: no numbering: none renumbering for interior faces: cell adjacency pre-ordering: lowest id first faces adjacent to ghost cells last: no numbering: none renumbering for boundary faces: numbering: none ---------------------------------------------------------- Matrix bandwidth for volume mesh : 416480 Matrix profile/lines for volume mesh : 43586 Mesh coordinates: minimum and maximum X : -4.1000000e+02 1.0000000e+01 Y : 0.0000000e+00 1.2000000e+03 Z : -1.0000000e+01 1.0000000e+01 Mesh Number of cells: 427500 Number of interior faces: 1269975 Number of boundary faces: 25050 Number of vertices: 440436 Groups: "Group 26" boundary faces: 22800 "Group 27" boundary faces: 1125 "Group 28" boundary faces: 1125 --- Information on the volumes Minimum control volume = 2.3480806e-01 Maximum control volume = 2.3941333e+00 Total volume for the domain = 4.7479290e+05 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 0 --> 0 % Criterion 3: Least-Squares Gradient Quality: Number of bad cells detected: 0 --> 0 % Criterion 4: Cells Volume Ratio: Number of bad cells detected: 0 --> 0 % Criterion 5: Guilt by Association: Number of bad cells detected: 0 --> 0 % Computing geometric quantities (0.47 s) ALMAX = 0.77674E+02 (Characteristic length ) ALMAX is the cubic root of the domain volume. ALMAX is the length used to initialize the turbulence. ----------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Variable Min. value Max. value ----------------------------------------- Velocity 0.0000E+00 0.0000E+00 Velocity 0.7430E-01 0.7430E-01 Velocity 0.0000E+00 0.0000E+00 Pressure 0.1500E+08 0.1500E+08 Turb Kinetic Ene 0.7500E+00 0.7500E+00 Omega 0.3000E-01 0.3000E-01 --------------------------------- ------------------------------------------------------------- ------------------------------------------------------------- =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.100000000E+01 TIME STEP NUMBER 1 ============================================================= ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.4283E+02 0.4283E+02 molecular_viscos 0.3182E-04 0.3182E-04 turbulent_viscos 0.1071E+04 0.1071E+04 ----------------------------------------- /home/martian/APP-INST/code_saturne-4.0.2/src/fvm/fvm_selector_postfix.c:1281: Fatal error. Error parsing expression: Group 28 ^ Expected operator instead of operand. Call stack: 1: 0x7f7ef90f0db9 <+0x1ffdb9> (libsaturne.so.0) 2: 0x7f7ef90f220e <+0x20120e> (libsaturne.so.0) 3: 0x7f7ef90f3d2b (libsaturne.so.0) 4: 0x7f7ef911770d <+0x22670d> (libsaturne.so.0) 5: 0x7f7ef9119afa (libsaturne.so.0) 6: 0x7f7ef900a111 (libsaturne.so.0) 7: 0x7f7ef900aa4e (libsaturne.so.0) 8: 0x40325b (cs_solver) 9: 0x7f7ef9095f90 (libsaturne.so.0) 10: 0x7f7ef8f73a96 (libsaturne.so.0) 11: 0x7f7ef8f5513e (libsaturne.so.0) 12: 0x7f7ef8f54c32 (libsaturne.so.0) 13: 0x7f7ef85c4ec5 <__libc_start_main+0xf5> (libc.so.6) 14: 0x4030b9 <> (cs_solver) End of stack