mus_species_module Module

  1. mus_species_module.f90
  2. mus_species_module

This module contains mus_species_type and routines to load species table from config file.


Uses

  • module~~mus_species_module~~UsesGraph module~mus_species_module mus_species_module module~tem_tools_module tem_tools_module module~mus_species_module->module~tem_tools_module module~mus_physics_module mus_physics_module module~mus_species_module->module~mus_physics_module module~env_module env_module module~mus_species_module->module~env_module module~aotus_module aotus_module module~mus_species_module->module~aotus_module module~tem_aux_module tem_aux_module module~mus_species_module->module~tem_aux_module module~aot_table_module aot_table_module module~mus_species_module->module~aot_table_module module~aot_out_module aot_out_module module~mus_species_module->module~aot_out_module module~tem_param_module tem_param_module module~mus_species_module->module~tem_param_module module~tem_logging_module tem_logging_module module~mus_species_module->module~tem_logging_module module~mus_physics_module->module~tem_tools_module module~mus_physics_module->module~env_module module~mus_physics_module->module~aotus_module module~mus_physics_module->module~tem_aux_module module~mus_physics_module->module~aot_table_module module~mus_physics_module->module~aot_out_module module~mus_physics_module->module~tem_logging_module module~treelmesh_module treelmesh_module module~mus_physics_module->module~treelmesh_module module~tem_geometry_module tem_geometry_module module~mus_physics_module->module~tem_geometry_module

Used by

  • module~~mus_species_module~~UsedByGraph module~mus_species_module mus_species_module module~mus_field_module mus_field_module module~mus_field_module->module~mus_species_module module~mus_bc_species_module mus_bc_species_module module~mus_bc_species_module->module~mus_species_module module~mus_field_prop_module mus_field_prop_module module~mus_field_prop_module->module~mus_species_module

Contents

Derived Types

Subroutines


Derived Types

type, public :: mus_species_type

this type contains species parameters

Read more…

Components

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: molWeight

molecular weight of the species
real(kind=rk), private :: molWeightInv

Inverse of molecular weight of the species. This parameter is required to convert mass density to mole density
real(kind=rk), private :: molWeigRatio

ratio of molecular weight \f$ \phi_\sigma = min(M)/M_\sigma i \f$
real(kind=rk), private, allocatable:: diff_coeff(:)

coefficient of diffusivity of the species (size of nspecies)
real(kind=rk), private, allocatable:: resi_coeff(:)

coefficient of resisivity of species which is reciprocal of diffusivity of the species KM:

Read more…
type(mrt_species_type), private :: mrt(globalMaxLevels)

molar fraction of this species in the mixture Volume fraction of is species in the mixture mrt relaxation for each level
real(kind=rk), private :: omBulk

bulk relaxation parameter omBulk_k = (2-phi_k)/3*bulkViscosity
real(kind=rk), private :: ombulkLvl(globalMaxLevels)

bulk relaxation parameter for each level
real(kind=rk), private :: omega

relaxation paramete for Nernst-Planck equation
real(kind=rk), private :: chargeNr

charge number of the species

type, private :: mrt_species_type

MRT species type

Components

TypeVisibilityAttributesNameInitial
real(kind=rk), private, allocatable:: s_mrt(:,:)

relaxation matrix for mrt size of this matrix is (layout%QQ, layout%QQ)
real(kind=rk), private, allocatable:: omegaMoments(:,:)

transformed relaxation matrix-moments factor omegaMoments = (Moments^-1.s_mrt.Moments) .(I+(Moments^-1.s_mrt.Moments)/2.0)^-1
real(kind=rk), private, allocatable:: omegaMomForce(:,:)

Omega factor for 2nd order force term omegaMomForce = (I+(Moments^-1.s_mrt.Moments)/2.0)^-1

Subroutines

public subroutine mus_load_species(me, conf, parent, minLevel, nFields, physics, cs_lattice)

this routines load species table from config file

Read more…

Arguments

TypeIntentOptionalAttributesName
type(mus_species_type), intent(out) :: me
type(flu_State) :: conf
integer, intent(in), optional :: parent
integer, intent(in) :: minLevel
integer, intent(in) :: nFields
type(mus_physics_type), intent(in) :: physics

physics type to convert physics to lattice unit or vice versa
real(kind=rk), intent(in) :: cs_lattice

lattice speed of sound calculated for defined stencil layout required to compute omega from potential diffusivity

public subroutine compute_molWeightRatio(molWeights, molWeigRatios)

This routine computes the molecular weight ratio for all species based asinari model

Read more…

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: molWeights(:)

molecular weight of the species
real(kind=rk), intent(out) :: molWeigRatios(:)

ratio of molecular weight \f$ \phi_\sigma = min(M)/M_\sigma \f$

public subroutine compute_bulkViscOmega(species, bulkvisc, bulkviscLvl, minLevel, maxLevel)

This routine compute bulk viscosity omega for species for all levels omega_bulk = (2-molWeigRatio_k)/(3*bulk_visc)

Arguments

TypeIntentOptionalAttributesName
type(mus_species_type), intent(inout) :: species

contains species information
real(kind=rk), intent(in) :: bulkvisc

bulk viscosity of the mixture
real(kind=rk), intent(in) :: bulkviscLvl(globalMaxLevels)
integer, intent(in) :: minLevel
integer, intent(in) :: maxLevel

public subroutine mus_species_out(me, conf)

writes species propertries into a lua file

Arguments

TypeIntentOptionalAttributesName
type(mus_species_type), intent(in) :: me
type(aot_out_type) :: conf