mus_addElectricToAuxField_MSL_WTDF Subroutine

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

This routine add electric force to momentum in auxField for multispecies liquid model with thermodynamic factor 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_wtdf~~CallsGraph proc~mus_addelectrictoauxfield_msl_wtdf mus_addElectricToAuxField_MSL_WTDF interface~mus_calc_thermfactor mus_calc_thermFactor proc~mus_addelectrictoauxfield_msl_wtdf->interface~mus_calc_thermfactor species species proc~mus_addelectrictoauxfield_msl_wtdf->species proc~invert_matrix invert_matrix proc~mus_addelectrictoauxfield_msl_wtdf->proc~invert_matrix proc~mus_calc_thermfactor_single mus_calc_thermFactor_single interface~mus_calc_thermfactor->proc~mus_calc_thermfactor_single interface~calc_therm_factor_loc calc_therm_factor_loc proc~mus_calc_thermfactor_single->interface~calc_therm_factor_loc

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
real(kind=rk), private, dimension(varSys%nStateVars):: num_dens
real(kind=rk), private, dimension(varSys%nStateVars):: moleFrac
integer, private :: iElem
integer, private :: nElems
integer, private :: elemOff
integer, private :: nInputStates
integer, private :: iField
integer, private :: iField_2
real(kind=rk), private :: electricField(fun%elemLvl(iLevel)%nElems*3)
real(kind=rk), private, dimension(varSys%nStateVars, varSys%nStateVars):: thermodynamic_fac
real(kind=rk), private, dimension(varSys%nStateVars, varSys%nStateVars):: inv_thermodyn_fac
real(kind=rk), private, dimension(3, varSys%nStateVars ):: spcForce
real(kind=rk), private, dimension(3, varSys%nStateVars ):: spcForce_WTDF
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