7.3
general documentation
Variables
Time stepping
Module for calculation options
+ Collaboration diagram for Time stepping:

Variables

integer(c_int), pointer, save ischtp
 time order of time stepping More...
 
integer, save istmpf
 time order of the mass flux scheme The chosen value for istmpf will automatically determine the value given to the variable thetfl. More...
 
integer(c_int), pointer, save itpcol
 Time scheme option: More...
 
integer(c_int), pointer, save nterup
 number of iterations on the velocity-pressure coupling on Navier-Stokes (for the U/P inner iterations scheme) More...
 
integer, save isno2t
 isno2t specifies the time scheme activated for the source terms of the momentum equation, apart from convection and diffusion (for instance: head loss, transposed gradient, ...). More...
 
integer(c_int), pointer, save isto2t
 Time scheme for source terms of turbulence equations (see isto2t in cs_time_scheme_t). More...
 
integer, save initvi
 initvi : =1 if total viscosity read from checkpoint file More...
 
integer, save initro
 initro : =1 if density read from checkpoint file More...
 
integer, save initcp
 initcp : =1 if specific heat read from checkpoint file More...
 
double precision, save thetsn
 $ \theta_S $-scheme for the source terms $S_e$ in the Navier-Stokes equations when the source term extrapolation has been activated (see isno2t), following the formula $(S_e)^{n+\theta}=(1+\theta)S_e^n-\theta S_e^{n-1}$.
The value of $theta$ = thetsn is deduced from the value chosen for isno2t. Generally only the value 0.5 is used. More...
 
real(c_double), pointer, save thetst
 The value of $theta$ (see thetst in cs_time_scheme_t). More...
 
double precision, save thetvi
 $ \theta $-scheme for the extrapolation of the physical property $\phi$ "total viscosity" when the extrapolation has been activated (see time_extrapolated key word), according to the formula $\phi^{n+\theta}=(1+\theta)\phi^n-\theta \phi^{n-1}$.
The value of $\theta$ = thetvi is deduced from the value chosen for time_extrapolated key word for the viscosity. Generally, only the value 0.5 is used. More...
 
double precision, save thetcp
 $ \theta $-scheme for the extrapolation of the physical property $\phi$ "specific heat" when the extrapolation has been activated (see time_extrapolated field key int), according to the formula $\phi^{n+\theta}=(1+\theta)\phi^n-\theta \phi^{n-1}$.
The value of $\theta$ = thetcp is deduced from the value chosen for the specific heat. Generally, only the value 0.5 is used. More...
 
real(c_double), pointer, save epsup
 relative precision for the convergence test of the iterative process on velocity-pressure coupling (inner iterations) More...
 
real(c_double), pointer, save xnrmu
 norm of the increment $ \vect{u}^{k+1} - \vect{u}^k $ of the iterative process on velocity-pressure coupling (inner iterations) More...
 
real(c_double), pointer, save xnrmu0
 norm of $ \vect{u}^0 $ (used by velocity-pressure inner iterations) More...
 

Detailed Description

Variable Documentation

◆ epsup

real(c_double), pointer, save epsup

relative precision for the convergence test of the iterative process on velocity-pressure coupling (inner iterations)

◆ initcp

integer, save initcp

initcp : =1 if specific heat read from checkpoint file

◆ initro

integer, save initro

initro : =1 if density read from checkpoint file

◆ initvi

integer, save initvi

initvi : =1 if total viscosity read from checkpoint file

◆ ischtp

integer(c_int), pointer, save ischtp

time order of time stepping

  • 2: 2nd order
  • 1: 1st order (default)

◆ isno2t

integer, save isno2t

isno2t specifies the time scheme activated for the source terms of the momentum equation, apart from convection and diffusion (for instance: head loss, transposed gradient, ...).

  • 0: "standard" first-order: the terms which are linear functions of the solved variable are implicit and the others are explicit
  • 1: second-order: the terms of the form $S_i\phi$ which are linear functions of the solved variable $\phi$ are expressed as second-order terms by interpolation (according to the formula $(S_i\phi)^{n+\theta}=S_i^n[(1-\theta)\phi^n+\theta\phi^{n+1}]$, $\theta$ being given by the value of thetav associated with the variable $\phi$); the other terms $S_e$ are expressed as second-order terms by extrapolation (according to the formula $(S_e)^{n+\theta}=[(1+\theta)S_e^n-\theta S_e^{n-1}]$, $\theta$ being given by the value of thetsn = 0.5).
  • 2: the linear terms $S_i\phi$ are treated in the same way as when isno2t = 1; the other terms $S_e$ are extrapolated according to the same formula as when isno2t = 1, but with $\theta$= thetsn = 1. By default, isno2t is initialised to 1 (second-order) when the selected time scheme is second-order (ischtp = 2), otherwise to 0.

◆ istmpf

integer, save istmpf

time order of the mass flux scheme The chosen value for istmpf will automatically determine the value given to the variable thetfl.

  • 2: theta scheme with theta > 0 (theta=0.5 means 2nd order) the mass flow used in the momentum equations is extrapolated at n+ thetfl (= n+1/2) from the values at the two former time steps (Adams Bashforth); the mass flow used in the equations for turbulence and scalars is interpolated at time n+ thetfl (= n+1/2) from the values at the former time step and at the newly calculated $n+1$ time step.
  • 0: theta scheme with theta = 0 (explicit): the mass flow calculated at the previous time step is used in the convective terms of all the equations (momentum, turbulence and scalars)
  • 1: implicit scheme (default) : the mass flow calculated at the previous time step is used in the convective terms of the momentum equation, and the updated mass flow is used in the equations of turbulence and scalars. By default, istmpf=2 is used in the case of a second-order time scheme (if ischtp=2) and istmpf = 1 otherwise.

◆ isto2t

integer(c_int), pointer, save isto2t

Time scheme for source terms of turbulence equations (see isto2t in cs_time_scheme_t).

◆ itpcol

integer(c_int), pointer, save itpcol

Time scheme option:

  • 0: staggered time scheme. On the time grids, the velocity is half a time step behind the density and the buoyant scalar. (See the thesis of [Pierce:2004])
  • 1: collocated time scheme. On the time grids, the velocity is at the same location as the density and the buoyant scalar. (See [Ma:2019])

◆ nterup

integer(c_int), pointer, save nterup

number of iterations on the velocity-pressure coupling on Navier-Stokes (for the U/P inner iterations scheme)

◆ thetcp

double precision, save thetcp

$ \theta $-scheme for the extrapolation of the physical property $\phi$ "specific heat" when the extrapolation has been activated (see time_extrapolated field key int), according to the formula $\phi^{n+\theta}=(1+\theta)\phi^n-\theta \phi^{n-1}$.
The value of $\theta$ = thetcp is deduced from the value chosen for the specific heat. Generally, only the value 0.5 is used.

  • 0 : explicit
  • 1/2: extrapolated in n+1/2
  • 1 : extrapolated in n+1

◆ thetsn

double precision, save thetsn

$ \theta_S $-scheme for the source terms $S_e$ in the Navier-Stokes equations when the source term extrapolation has been activated (see isno2t), following the formula $(S_e)^{n+\theta}=(1+\theta)S_e^n-\theta S_e^{n-1}$.
The value of $theta$ = thetsn is deduced from the value chosen for isno2t. Generally only the value 0.5 is used.

  • 0 : second viscosity explicit
  • 1/2: second viscosity extrapolated in n+1/2
  • 1 : second viscosity extrapolated in n+1

◆ thetst

real(c_double), pointer, save thetst

The value of $theta$ (see thetst in cs_time_scheme_t).

◆ thetvi

double precision, save thetvi

$ \theta $-scheme for the extrapolation of the physical property $\phi$ "total viscosity" when the extrapolation has been activated (see time_extrapolated key word), according to the formula $\phi^{n+\theta}=(1+\theta)\phi^n-\theta \phi^{n-1}$.
The value of $\theta$ = thetvi is deduced from the value chosen for time_extrapolated key word for the viscosity. Generally, only the value 0.5 is used.

  • 0 : explicit
  • 1/2: extrapolated in n+1/2
  • 1 : extrapolated in n+1

◆ xnrmu

real(c_double), pointer, save xnrmu

norm of the increment $ \vect{u}^{k+1} - \vect{u}^k $ of the iterative process on velocity-pressure coupling (inner iterations)

◆ xnrmu0

real(c_double), pointer, save xnrmu0

norm of $ \vect{u}^0 $ (used by velocity-pressure inner iterations)