mus_moments_module Module

This module deals with the calculation of moments from pdfs


Uses

Used by

  • module~~mus_moments_module~~UsedByGraph module~mus_moments_module mus_moments_module module~mus_derivedquantities_module2 mus_derivedQuantities_module2 module~mus_derivedquantities_module2->module~mus_moments_module module~mus_interpolate_verify_module mus_interpolate_verify_module module~mus_interpolate_verify_module->module~mus_moments_module module~mus_scheme_module mus_scheme_module module~mus_scheme_module->module~mus_moments_module module~mus_config_module mus_config_module module~mus_config_module->module~mus_scheme_module module~mus_d2q9_module mus_d2q9_module module~mus_d2q9_module->module~mus_derivedquantities_module2 module~mus_d3q19_module mus_d3q19_module module~mus_d3q19_module->module~mus_derivedquantities_module2 module~mus_d3q27_module mus_d3q27_module module~mus_d3q27_module->module~mus_derivedquantities_module2 module~mus_derquan_module mus_derQuan_module module~mus_derquan_module->module~mus_derivedquantities_module2 module~mus_derquanincomp_module mus_derQuanIncomp_module module~mus_derquanincomp_module->module~mus_derivedquantities_module2 module~mus_derquanpoisson_module mus_derQuanPoisson_module module~mus_derquanpoisson_module->module~mus_derivedquantities_module2 module~mus_dynloadbal_module mus_dynLoadBal_module module~mus_dynloadbal_module->module~mus_scheme_module module~mus_flow_module mus_flow_module module~mus_flow_module->module~mus_derivedquantities_module2 module~mus_flow_module->module~mus_interpolate_verify_module module~mus_hvs_config_module mus_hvs_config_module module~mus_hvs_config_module->module~mus_scheme_module module~mus_ibm_module mus_IBM_module module~mus_ibm_module->module~mus_derivedquantities_module2 module~mus_interpolate_average_module mus_interpolate_average_module module~mus_interpolate_average_module->module~mus_derivedquantities_module2 module~mus_interpolate_linear_module mus_interpolate_linear_module module~mus_interpolate_linear_module->module~mus_derivedquantities_module2 module~mus_interpolate_quadratic_module mus_interpolate_quadratic_module module~mus_interpolate_quadratic_module->module~mus_derivedquantities_module2 module~mus_nonnewtonian_module mus_nonNewtonian_module module~mus_nonnewtonian_module->module~mus_derivedquantities_module2 module~mus_program_module mus_program_module module~mus_program_module->module~mus_scheme_module module~mus_smagorinsky_module mus_Smagorinsky_module module~mus_smagorinsky_module->module~mus_derivedquantities_module2 module~mus_tools_module mus_tools_module module~mus_tools_module->module~mus_scheme_module program~mus_harvesting mus_harvesting program~mus_harvesting->module~mus_scheme_module

Functions

public pure function get_moment(QQ, cxDir, expX, pdf) result(mom)

Calculate the moment of a centain order The moment of a distribution is defined as:\n The fucntion argument expX is array of size 3, which contains the values of \f$p, q, r\f$

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: QQ
integer, intent(in) :: cxDir(3,QQ)
integer, intent(in) :: expX(3)
real(kind=rk), intent(in) :: pdf(QQ)

distribution value

Return Value real(kind=rk)

public pure function get_momentVector(QQ, cxDir, expX) result(mom)

get the moment vector to calculate the moment from the pdf

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: QQ
integer, intent(in) :: cxDir(3,QQ)
integer, intent(in) :: expX(3)

exponents of the moments

Return Value real(kind=rk), (QQ)

moment vector

private pure function mus_iMomVector(cxDir, expX, QQ) result(iMom)

The integer moment vector for a given cxDir and order.

Read more…

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: cxDir(:,:)

discrete velocity
integer, intent(in) :: expX(3)

order in each direction
integer, intent(in) :: QQ

number of velocity channels (include rest)

Return Value integer, (QQ)


Subroutines

public subroutine mus_init_moments(me, QQ, cxDir, label, schemeHeader)

Initialize the moment space

Arguments

Type IntentOptional Attributes Name
type(mus_moment_type), intent(inout) :: me
integer, intent(in) :: QQ
integer, intent(in) :: cxDir(3,QQ)
character(len=labelLen) :: label
type(mus_scheme_header_type), intent(in) :: schemeHeader

identifier of the scheme

public subroutine set_momentIndices(nDims, iPress, iVelMin, iVelMax, iSMin, iSMax)

set indices for accessing the pressure, velocity and the shear from a 1d vector

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nDims

number of dimensions
integer, intent(out) :: iPress

index for the pressure / density
integer, intent(out) :: iVelMin

starting index for velocity
integer, intent(out) :: iVelMax

ending index for velocity
integer, intent(out) :: iSMin

starting index for shear
integer, intent(out) :: iSMax

ending index for shear

public subroutine mus_dump_moments(me, outUnit)

Dump moments matrix: toPDF and toMoment

Arguments

Type IntentOptional Attributes Name
type(mus_moment_type), intent(in) :: me
integer, intent(in) :: outUnit

private subroutine init_transformation_matrix_fluid(QQ, cxDir, label, me, toMoment, toPdf)

Initialize Moments transformation matrix for LBM compressible and incompressible fluid model. This matrix must be consistent with the relaxation matrix used in compute kernel and interpolation routines

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: QQ
integer, intent(in) :: cxDir(3,QQ)
character(len=labelLen) :: label
type(mus_moment_type), intent(inout) :: me
real(kind=rk), intent(inout) :: toMoment(me%toMoments%nEntries(1),me%toMoments%nEntries(2))
real(kind=rk), intent(inout) :: toPdf(me%toPDF%nEntries(1),me%toPDF%nEntries(2))

private subroutine init_transformation_matrix_MS(QQ, cxDir, label, me, toMoment, toPdf)

Intialize the moment transformation matrix for multispecies. This matrix must be consistent with relaxation matrix used for multispecies MRT collision routines

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: QQ
integer, intent(in) :: cxDir(3,QQ)
character(len=labelLen) :: label
type(mus_moment_type), intent(inout) :: me
real(kind=rk), intent(inout) :: toMoment(me%toMoments%nEntries(1),me%toMoments%nEntries(2))
real(kind=rk), intent(inout) :: toPdf(me%toPDF%nEntries(1),me%toPDF%nEntries(2))