Skip to content

adam_turbulence_kernels ​

ADAM, CPU kernels of the dynamic eddy viscosity (issue #119, T1).

mu_t = rho nu_t on the interior cells and on 2 ghost layers per active direction: the 4th-order dissipative fluxes read it at the cells -1 ... 2 along their direction. 2nd-order central velocity gradients, so ngc >= 3. The filter width is the geometric mean of the active steps of the block, and so it jumps at a 2:1 face. One kernel per model (the same body), selected on the host; the app passes the auxiliary indexes of the density and of the first velocity component.

Source: src/lib/common/adam_turbulence_kernels.F90

Dependencies

Contents ​

Subroutines ​

compute_mu_t_vreman ​

Vreman.

fortran
subroutine compute_mu_t_vreman(ir, iu, ni, nj, nk, ngc, blocks_number, c, dxyz, is_null, q_aux, mu_t)

Arguments

NameTypeIntentAttributesDescription
irinteger(kind=I4P)inAuxiliary indexes of the density and of u.
iuinteger(kind=I4P)inAuxiliary indexes of the density and of u.
niinteger(kind=I4P)inGrid dimensions.
njinteger(kind=I4P)inGrid dimensions.
nkinteger(kind=I4P)inGrid dimensions.
ngcinteger(kind=I4P)inGrid dimensions.
blocks_numberinteger(kind=I4P)inActual blocks number.
creal(kind=R8P)inModel constant.
dxyzreal(kind=R8P)inBlocks space steps [3, nb].
is_nulllogicalinNull directions.
q_auxreal(kind=R8P)inAuxiliary variables.
mu_treal(kind=R8P)inoutEddy viscosity.

Call graph

compute_mu_t_wale ​

WALE.

fortran
subroutine compute_mu_t_wale(ir, iu, ni, nj, nk, ngc, blocks_number, c, dxyz, is_null, q_aux, mu_t)

Arguments

NameTypeIntentAttributesDescription
irinteger(kind=I4P)inAuxiliary indexes of the density and of u.
iuinteger(kind=I4P)inAuxiliary indexes of the density and of u.
niinteger(kind=I4P)inGrid dimensions.
njinteger(kind=I4P)inGrid dimensions.
nkinteger(kind=I4P)inGrid dimensions.
ngcinteger(kind=I4P)inGrid dimensions.
blocks_numberinteger(kind=I4P)inActual blocks number.
creal(kind=R8P)inModel constant.
dxyzreal(kind=R8P)inBlocks space steps [3, nb].
is_nulllogicalinNull directions.
q_auxreal(kind=R8P)inAuxiliary variables.
mu_treal(kind=R8P)inoutEddy viscosity.

Call graph

compute_mu_t_sigma ​

Sigma.

fortran
subroutine compute_mu_t_sigma(ir, iu, ni, nj, nk, ngc, blocks_number, c, dxyz, is_null, q_aux, mu_t)

Arguments

NameTypeIntentAttributesDescription
irinteger(kind=I4P)inAuxiliary indexes of the density and of u.
iuinteger(kind=I4P)inAuxiliary indexes of the density and of u.
niinteger(kind=I4P)inGrid dimensions.
njinteger(kind=I4P)inGrid dimensions.
nkinteger(kind=I4P)inGrid dimensions.
ngcinteger(kind=I4P)inGrid dimensions.
blocks_numberinteger(kind=I4P)inActual blocks number.
creal(kind=R8P)inModel constant.
dxyzreal(kind=R8P)inBlocks space steps [3, nb].
is_nulllogicalinNull directions.
q_auxreal(kind=R8P)inAuxiliary variables.
mu_treal(kind=R8P)inoutEddy viscosity.

Call graph