mus_addElectricToAuxField_MSL Subroutine

public subroutine mus_addElectricToAuxField_MSL(fun, auxField, iLevel, time, varSys, phyConvFac, derVarPos)

This routine add electric force to momentum in auxField for multispecies liquid model Refer to Appendix in PhD Thesis of K. Masilamani "Coupled Simulation Framework to Simulate Electrodialysis Process for Seawater Desalination"

Arguments

TypeIntentOptionalAttributesName
class(mus_source_op_type), intent(inout) :: fun

Description of method to update source

real(kind=rk), intent(inout) :: auxField(:)

output auxField array

integer, intent(in) :: iLevel

current level

type(tem_time_type), intent(in) :: time

current timing information

type(tem_varSys_type), intent(in) :: varSys

variable system definition

type(mus_convertFac_type), intent(in) :: phyConvFac

Physics conversion factor for current level

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

position of derived quantities in varsys


Calls

proc~~mus_addelectrictoauxfield_msl~~CallsGraph proc~mus_addelectrictoauxfield_msl mus_addElectricToAuxField_MSL species species proc~mus_addelectrictoauxfield_msl->species

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: dens_pos
integer, private :: mom_pos(3)
integer, private :: depField
real(kind=rk), private, dimension(varSys%nStateVars):: mass_dens
real(kind=rk), private, dimension(varSys%nStateVars):: massFrac
integer, private :: iElem
integer, private :: nElems
integer, private :: elemOff
integer, private :: nInputStates
integer, private :: iField
real(kind=rk), private :: electricField(fun%elemLvl(iLevel)%nElems*3)
real(kind=rk), private :: EF_elem(3)
real(kind=rk), private, dimension(3, varSys%nStateVars ):: spcForce
real(kind=rk), private :: charge_dens
real(kind=rk), private :: diffForce_cs2inv
real(kind=rk), private :: minMolWeight
real(kind=rk), private, dimension(varSys%nStateVars):: chargeTerm
type(mus_varSys_data_type), private, pointer:: fPtr