mus_mixture_module Module

This module contains all information about fluid mixture


Uses

Used by

  • module~~mus_mixture_module~~UsedByGraph module~mus_mixture_module mus_mixture_module module~mus_bc_fluid_experimental_module mus_bc_fluid_experimental_module module~mus_bc_fluid_experimental_module->module~mus_mixture_module module~mus_bc_fluid_module mus_bc_fluid_module module~mus_bc_fluid_module->module~mus_mixture_module module~mus_bc_fluid_noneqexpol_module mus_bc_fluid_nonEqExpol_module module~mus_bc_fluid_noneqexpol_module->module~mus_mixture_module module~mus_bc_fluid_turbulent_module mus_bc_fluid_turbulent_module module~mus_bc_fluid_turbulent_module->module~mus_mixture_module module~mus_bc_fluid_wall_module mus_bc_fluid_wall_module module~mus_bc_fluid_wall_module->module~mus_mixture_module module~mus_bc_general_module mus_bc_general_module module~mus_bc_general_module->module~mus_mixture_module module~mus_bc_header_module mus_bc_header_module module~mus_bc_header_module->module~mus_mixture_module module~mus_bc_nernstplanck_module mus_bc_nernstPlanck_module module~mus_bc_nernstplanck_module->module~mus_mixture_module module~mus_bc_passivescalar_module mus_bc_passiveScalar_module module~mus_bc_passivescalar_module->module~mus_mixture_module module~mus_bc_poisson_module mus_bc_poisson_module module~mus_bc_poisson_module->module~mus_mixture_module module~mus_bc_species_module mus_bc_species_module module~mus_bc_species_module->module~mus_mixture_module module~mus_field_module mus_field_module module~mus_field_module->module~mus_mixture_module module~mus_flow_module mus_flow_module module~mus_flow_module->module~mus_mixture_module module~mus_msliquid_module mus_MSLiquid_module module~mus_msliquid_module->module~mus_mixture_module module~mus_scheme_module mus_scheme_module module~mus_scheme_module->module~mus_mixture_module module~mus_scheme_type_module mus_scheme_type_module module~mus_scheme_type_module->module~mus_mixture_module

Derived Types

type, public ::  mus_mixture_type

This type contains mixture information

Read more…

Components

Type Visibility Attributes Name Initial
type(tem_ini_condition_type), public :: ic

initialization case, initial condition of the mixture
real(kind=rk), public :: rho0

mass density of the mixture SI unit: kg/m^3. Physics to lattice conversion: rho0/physics%rho0
real(kind=rk), public :: rho0LB

lattice mass density of the mixture
real(kind=rk), public :: moleDens0

number density of the mixture or total mixture molar density SI unit: mol/m^3. physics to lattice conversion: nT0/physics%mol*physics%fac(minlevel)%length^3. mixture molar density is required if initial condition are defined by molar fraction. set only if initial molefraction is space independent
real(kind=rk), public :: moleDens0LB

lattice number density of the mixture or total mixture molar density physics to lattice conversion: moleDens0/physics%mol*physics%fac(minlevel)%length^3.
real(kind=rk), public :: kine_viscosityLB

lattice kinematic shear viscosity of the mixture
real(kind=rk), public :: kine_viscosity

physical kinematic shear viscosity of the mixture
real(kind=rk), public :: bulk_viscosityLB

lattice bulk viscosity
real(kind=rk), public :: bulk_viscosity

physical bulk viscosity
real(kind=rk), public :: bulk_modulusLB

lattice bulk modulus of the liquid mixture , - speed of sound (in lattice unit: \f$ 1/\sqrt{3} \f$
real(kind=rk), public :: omega_diff

relaxation parameter, , B - free parameter unit same as resistivity
real(kind=rk), public :: omega_kine
real(kind=rk), public :: omega_hom

omega for higher order moments
type(mixRelaxation_type), public :: relaxLvl(globalMaxLevels)

relaxation parameters for each level
type(tem_temporal_type), public :: omega_ramping

temporal omega for ramping etc.
real(kind=rk), public :: paramB

free parameter \f$ B = \omega*\rho/K \f$
real(kind=rk), public :: temp0

temperature
real(kind=rk), public :: temp0LB

temperature
real(kind=rk), public :: theta_eq

equilibrium theta to choose between mixture velocity and equilibrium species velocity in the quadratic term equilibrium function. theta = 0 -> mixture velocity theta = 1 -> equilibrium species velocity \todo KM: remove theta_eq and use mixture velocity in quadratic term of equilibrium function
type(tem_spatial_type), public :: viscSpatial

spatial omega definition, e.g. for sponge layers
real(kind=rk), public :: electricField(3)

external electrical force
real(kind=rk), public :: gravityField(3)

gravitational force
real(kind=rk), public :: faraday

faraday constant(C/mol)
real(kind=rk), public :: faradayLB

faraday constant in lattice
real(kind=rk), public :: gasConst_R

gas constant R (Nm/(mol *K))
real(kind=rk), public :: gasConst_R_LB

gas constant R in lattice
character(kind=c_char, len=labelLen), public :: prop_file

eNRTL file with species properties
real(kind=rk), public :: atm_press

atmospheric pressure
real(kind=rk), public :: atm_pressLB

atmospheric pressure

type, private ::  mixRelaxation_type

relaxation paraemters multispecies required for each level

Components

Type Visibility Attributes Name Initial
real(kind=rk), public :: visc

lattice shear viscosity of the mixture for different level
real(kind=rk), public :: paramB

lattice free parameter B for different level
real(kind=rk), public :: bulkVisc

lattice bulk viscosity for different level
real(kind=rk), public :: omega_diff

relaxation parameter, , B - free parameter unit same as resistivity
real(kind=rk), public :: omega_kine

Subroutines

public subroutine mus_load_mixture(me, conf, parent, minLevel, maxLevel, physics, schemeHeader, nFields)

This routine load mixture table from scheme table. Define either mass density or number density. If mass density is specified, number density can be computed at runtime or vice versa. @note Todo Currently, the simulation is initialized by density, extend it to initialize from mixture number density/volume fraction and mole fraction \verbatim mixture = { rho0 = 1.0, omega } \endverbatim

Read more…

Arguments

Type IntentOptional Attributes Name
type(mus_mixture_type), intent(out) :: me

contains mixture information
type(flu_State) :: conf
integer, intent(in), optional :: parent
integer, intent(in) :: minLevel
integer, intent(in) :: maxLevel
type(mus_physics_type), intent(in) :: physics

physics type to convert physics to lattice unit or vice versa
type(mus_scheme_header_type), intent(in) :: schemeHeader

identifier of the scheme
integer, intent(in) :: nFields

number of fields defined in lua file

public subroutine mus_mixture_out(me, conf, schemeHeader)

This routine write mixture properties into lua file

Arguments

Type IntentOptional Attributes Name
type(mus_mixture_type), intent(in) :: me

mixture info
type(aot_out_type) :: conf
type(mus_scheme_header_type), intent(in) :: schemeHeader

identifier of the scheme

private subroutine set_omegasLvl(mixture, minLevel, maxLevel, physics)

Set the omegas according to the time step setting

Read more…

Arguments

Type IntentOptional Attributes Name
type(mus_mixture_type), intent(inout) :: mixture
integer, intent(in) :: minLevel
integer, intent(in) :: maxLevel
type(mus_physics_type), intent(in) :: physics

physics type to convert physics to lattice unit or vice versa