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"


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


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



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