mus_addElectricToAuxField_MSL_WTDF Subroutine

  1. mus_auxFieldVar_module.f90
  2. mus_auxFieldVar_module
  3. mus_addElectricToAuxField_MSL_WTDF

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


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