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"
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| 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 |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| 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 |