setParameters_multispecies Subroutine

public subroutine setParameters_multispecies(field, nFields, mixture, header, layout, iLevel, tNow)

Set parameters for multispecies

Relaxation parameter for each level kine_viscosity = dtc^2(1/omega - 0.5) => omega = dt / (viscosity/c^2 + dt/2) omega for each level is stored at fluid%omLvl

Arguments

TypeIntentOptionalAttributesName
type(mus_field_type), intent(inout) :: field(nFields)
integer, intent(in) :: nFields
type(mus_mixture_type), intent(inout) :: mixture
type(mus_scheme_header_type), intent(in) :: header
type(mus_scheme_layout_type), intent(in) :: layout
integer, intent(in) :: iLevel
type(tem_time_type), intent(in) :: tNow

solver general info


Calls

proc~~setparameters_multispecies~~CallsGraph proc~setparameters_multispecies setParameters_multispecies proc~invert_matrix invert_matrix proc~setparameters_multispecies->proc~invert_matrix proc~tem_temporal_for tem_temporal_for proc~setparameters_multispecies->proc~tem_temporal_for

Called by

proc~~setparameters_multispecies~~CalledByGraph proc~setparameters_multispecies setParameters_multispecies proc~mus_update_relaxparams mus_update_relaxParams proc~mus_update_relaxparams->proc~setparameters_multispecies proc~do_recursive_multilevel do_recursive_multiLevel proc~do_recursive_multilevel->proc~mus_update_relaxparams proc~do_fast_singlelevel do_fast_singleLevel proc~do_fast_singlelevel->proc~mus_update_relaxparams proc~do_benchmark do_benchmark proc~do_benchmark->proc~mus_update_relaxparams

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: omega_diff
real(kind=rk), private :: omega_kine
real(kind=rk), private :: omega_bulk
real(kind=rk), private :: fac
integer, private :: iField
integer, private :: iDir
real(kind=rk), private, dimension( layout%fStencil%QQ, layout%fStencil%QQ ):: identity
real(kind=rk), private, dimension( layout%fStencil%QQ, layout%fStencil%QQ ):: tmpMatrix
real(kind=rk), private, dimension( layout%fStencil%QQ, layout%fStencil%QQ ):: tmpMatrixInv