mus_physics_module Module

length (m) = dx

time (s) = dt

velocity(m/s) = dx/dt

kinematic viscosity(m^2/s) = dx^2/dt

Dynamic viscosity (Pa s) = kg/m/s

acceleration(m/s^2) = dx/dt^2

Force(N)(kg m/s^2) = rho0*dx^4/dt^2

Force per unit volume (N/m^3)(kg/s^2/m^2) = rho0*dx/dt^2

Pressure(N/m^2)(kg/m/s^2) = rho0*dx^2/dt^2

Strain Rate (1/s) = 1/dt

Energy (N-m) (kgm^2/s^2) = rho0dx^5/dt^2

mole density(mol/m^3) = mole0/dx^3 Charge density (C/m^3) = Coulomb0/dx^3

Current density (C/s/m^2) = Coulomb0/dt/dx^2

mole flux(mol/m^2/s) = moleDen0*dx/dt

mass flux(kg/m^2/s) = rho0*dx/dt

diffusivity(m^2/s) = dx^2/dt

faraday (C/mol) = coulomb0/moleDens0/dx^3

gas constant (J/mol/K) (N m/mol/K) (kg m^2/s^2/mol/K) = rho0*dx^5/dt^2/mole0/temp0

Potential (V) (kg m^2/(C*S^2))

needed to check if physics table is defined

reference length - discretization size of the coarsest level SI unit - meter

reference time - time discretization for discretization size of the coarsest level SI unit - seconds

reference physical mass density SI unit - kg/m^3

reference physical mole density SI unit - mol/m^3

reference molecular weight SI unit - kg/mol

reference temperature SI unit - Kelvin

reference fundamental electrical charge SI unit - Coulomb

mole is defined by inverse of Avogadro Constant Avogadro Constant = 6.02214129e23 [1/mol]

reference mass in kg derived from density or moleweight SI unit :: kg

Level-wise conversion factor for derived variables size: minLevel:maxLevel allocated in mus_load_physics \todo KM: conversion factor should not be level-dependent. it should be same for all levels, the lattice dx and dt for each level must be considered to scale variables in multilevel. Implemented force, visc, etc using dtL according to formula Introduced lattice speed variable: dx/dt for each level it should be same for all level for acoustic scaling and different for diffusive scaling

Pressure (strain rate) over level scale factor. This factor is meant to convert pressure in LB unit on source level to the required pressure on target level. It is mainly used in interpolation routine. It is allocated as: allocate(pFac( minLevel:maxLevel, minLevel:maxLevel)) It is allocated and initialized in routine: mus_set_scaleFac How to use it in the code: pTargetLevel = pSourceLevel * pFac( sourceLevel, targetLevel )

Velocity over level scale factor Its usage is the same as pressure scale factor

Strain rate over level scale factor Its usage is the same as pressure scale factor

physics type

flu state

global treelm mesh

scaling factor: diffusive -> 4; acoustic -> 2

reference time step if none

reference spacestep if none

This string contains lua functions to compute dt from visocosity or velocity