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, dimension(nscamx), save	isso2t
	for each scalar, isso2t specifies the time scheme activated for the source terms of the equation for the scalar, apart from convection and diffusion (for instance: variance production, user-specified terms, ...). 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...

integer, dimension(nscamx), save	initvs
	initvs : =1 if scalar diffusivity read from checkpoint file More...

double precision, save	thetsn
	$\theta_S$ -scheme for the source terms 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 = 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 (see thetst in cs_time_scheme_t). More...

double precision, dimension(nscamx), save	thetss
	$\theta$ -scheme for the extrapolation of the nonlinear explicit source term of the scalar transport equation when the source term extrapolation has been activated (see isso2t), following the formula $(S_e)^{n+\theta}=(1+\theta)S_e^n-\theta S_e^{n-1}$ . The value of $\theta$ = thetss is deduced from the value chosen for isso2t. Generally, only the value 0.5 is used. 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...

double precision, dimension(nscamx), save	thetvs
	$\theta$ -scheme for the extrapolation of the physical property $\phi$ "diffusivity" 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$ = thetvs is deduced from the value chosen for time_extrapolated key word. 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

◆ initvs

integer, dimension(nscamx), save initvs

initvs : =1 if scalar diffusivity 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 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 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.

◆ isso2t

integer, dimension(nscamx), save isso2t

for each scalar, isso2t specifies the time scheme activated for the source terms of the equation for the scalar, apart from convection and diffusion (for instance: variance production, user-specified terms, ...).

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 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 thetss (iscal) = 0.5)
2: the linear terms $S_i\phi$ are treated in the same way as when isso2t = 1; the other terms are extrapolated according to the same formula as when isso2t = 1, but with $\theta$ = thetss (iscal) = 1.
By default, isso2t (iscal) 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 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

◆ thetss

double precision, dimension(nscamx), save thetss

$\theta$ -scheme for the extrapolation of the nonlinear explicit source term $S_e$ of the scalar transport equation when the source term extrapolation has been activated (see isso2t), following the formula $(S_e)^{n+\theta}=(1+\theta)S_e^n-\theta S_e^{n-1}$ .
The value of $\theta$ = thetss is deduced from the value chosen for isso2t. Generally, only the value 0.5 is used.

0 : explicit
1/2: extrapolated in n+1/2
1 : 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

◆ thetvs

double precision, dimension(nscamx), save thetvs

$\theta$ -scheme for the extrapolation of the physical property $\phi$ "diffusivity" 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$ = thetvs is deduced from the value chosen for time_extrapolated key word. 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)

Variables