command: /home/jundi/Code_Saturne/libexec/code_saturne/cs_solver --param example *************************************************************** (R) Code_Saturne Version 4.0.6 Copyright (C) 1998-2016 EDF S.A., France revision 4.0.6 build Wed 07 Dec 2016 13:18:02 GMT 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: Mon 16 Jan 2017 10:45:39 GMT System: Linux 4.4.0-53-generic Machine: jundi-pc Processor: model name : Intel(R) Core(TM) i7-2600 CPU @ 3.40GHz Memory: 5949 MB User: jundi (Jundi He) Directory: /home/jundi/Desktop/Example/example/fluid/RESU/20170116-1045_2 External libraries for partitioning: SCOTCH 5.1.12 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 ---------- --- Physics NVAR = 7 (Nb variables ) NSCAL = 1 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 1 (Nb specific physics scalars ) NPROCE = 6 (Nb cell properties ) ----------------------------------------------------------- ** PHYSICAL PROPERTIES ------------------- GX = 0.00000E+00 (Gravity x component ) GY = -0.98100E+01 (Gravity y component ) GZ = 0.00000E+00 (Gravity 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 ----------------------------------------------------------- ** 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 = 1 (0: no thermal model ) (1: temperature ) (2: enthalpy ) (3: total energy ) ITPSCL = 2 (0: none ) (1: temperature in Kelvin ) (2: temperature in Celsius ) ISCALT = 1 (Thermal scalar number ) ** TURBULENCE ---------- --- Continuous phase: - Commons ITURB = 20 (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-epsilon (ITURB = 20) ALMAX = -0.99900E+03 (Characteristic length ) UREF = 0.10000E+01 (Characteristic velocity ) ICLKEP = 0 (k-epsilon clipping model ) IKECOU = 1 (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-epsilon (ITURB = 20) Ce1 = 0.14400E+01 (Cepsilon 1: production coef.) CE2 = 0.19200E+01 (Cepsilon 2: dissipat. coef.) SIGMAK = 0.10000E+01 (Prandtl relative to k ) SIGMAE = 0.13000E+01 (Prandtl relative to epsilon ) CMU = 0.90000E-01 (Cmu constant ) ----------------------------------------------------------- ** 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-01 (Minimum time step ) DTMAX = 0.10000E+03 (Maximum time step ) DTREF = 0.10000E+00 (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 ------------------------------------ Velocity 1 0.1000E+01 Pressure 0 0.1000E+01 k 1 0.1000E+01 epsilon 1 0.1000E+01 TempC 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 --------------------------------------------------------------------- Velocity 1 1 1 1 0 0.10E+01 0.10E+01 Pressure 0 1 1 1 0 0.00E+00 0.10E+01 k 1 1 1 1 0 0.00E+00 0.10E+01 epsilon 1 1 1 1 0 0.00E+00 0.10E+01 TempC 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 ------ 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) 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 ------------------------------------------------------------------- Velocity 100 1 0.1000E-04 0.1000E-06 0.0000E+00 Pressure 100 2 0.1000E-04 0.1000E-06 0.0000E+00 k 100 1 0.1000E-04 0.1000E-06 0.0000E+00 epsilon 100 1 0.1000E-04 0.1000E-06 0.0000E+00 TempC 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 k 1 -1 0.1500E+01 epsilon 1 -1 0.1500E+01 TempC 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 ) ----------------------------------------------------------- ** BASE ITERATIVE SOLVERS ---------------------- ------------------------------------ Variable EPSILO IDIRCL ------------------------------------ Velocity 0.1000E-07 1 Pressure 0.1000E-07 1 k 0.1000E-07 1 epsilon 0.1000E-07 1 TempC 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 ) ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) ------------------------------------------------------------- Variable Number ISCACP ITURT VISLS0 SIGMAS ------------------------------------------------------------- TempC 1 1 0 0.2495E-01 0.1000E+01 --------------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ TempC 1 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- TempC 1 -1 -0.2731E+03 0.1000E+13 ------------------------------------------------------- ------------------------------------------------------------- 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. ISCACP = 0 or 1 2 (use Cp or not ) VISLS0 = >0 (Reference viscosity ) SIGMAS = >0 (Schmidt ) RVARFL = >0 (Rf, cf variance dissipation ) ICLVFL = 0, 1 or 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 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 = 2000 (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 k epsilon TempC TurbVisc CourantNb FourierNb total_pressure Local Time Step -- -- --- 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 k -1 7 epsilon -1 8 TempC -1 10 LamVisc -1 11 TurbVisc -1 12 CourantNb -1 13 FourierNb -1 14 total_pressure -1 15 Local Time Step -1 -- -- -- --- Log files NTLIST = 1 (Output frequency ) Number Name IWARNI verbosity level (-999: not applicable) 2 Velocity 0 5 Pressure 0 6 k 0 7 epsilon 0 8 TempC 0 11 TurbVisc -999 12 CourantNb -999 13 FourierNb -999 14 total_pressure -999 15 Local Time Step -999 -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 1 (Force exerted by the fluid on the boundary) ipstft = 1 (Thermal flux at boundary) ipsttb = 1 (Temperature at boundary) ipstnu = 0 (Dimensionless thermal flux at boundary) ----------------------------------------------------------- ** SYRTHES COUPLING ---------------- NBCCOU = 1 (Number of couplings ) with 1 surface coupling(s) with 0 volume coupling(s) -- Coupled scalars ------------------------------- Scalar Number ICPSYR ------------------------------- TempC 1 1 ----------------------- ICPSYR = 0 or 1 (1: scalar coupled to SYRTHES) ----------------------------------------------------------- ** 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) ----------------------------------------------------------- Unmatched SYRTHES couplings: ---------------------------- SYRTHES coupling: coupling id: 0 local name: "SYRTHES" ../../../code_saturne-4.0.6/src/base/cs_syr_coupling.c:937: Fatal error. At least 1 SYRTHES coupling was defined for which no communication with a SYRTHES instance is possible. Call stack: 1: 0x7f481af53858 (libsaturne.so.0) 2: 0x7f481aea4e01 (libsaturne.so.0) 3: 0x7f481aea4bdf (libsaturne.so.0) 4: 0x7f481aa93830 <__libc_start_main+0xf0> (libc.so.6) 5: 0x400879 <_start+0x29> (cs_solver) End of stack