derive untransformed pdf velocity of species by solving system of equations of nSpecies

Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|

real(kind=rk), | intent(in) | :: | moleFraction(nFields) | molefraction of all species |
||

real(kind=rk), | intent(in) | :: | first_Moments(3,nFields) | momentum from transformed pdf of all species |
||

integer, | intent(in) | :: | nFields | number of species |
||

real(kind=rk), | intent(in) | :: | inv_thermodyn_fac(nFields,nFields) | inverse of thermodynamic factor |
||

real(kind=rk), | intent(in) | :: | phi(nFields) | molecular weight ratio of all species |
||

real(kind=rk), | intent(in) | :: | resi_coeff(nFields,nFields) | resistivity coefficients |
||

real(kind=rk), | intent(in) | :: | omega_fac | relaxation parameter, omega_diff*0.5_rk |
||

real(kind=rk), | intent(in) | :: | paramBInv | free parameter B |

return actual momentum

Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

integer, | private | :: | ifield | ||||

integer, | private | :: | ifield_2 | ||||

integer, | private | :: | ifield_3 | ||||

integer, | private | :: | icomp | ||||

real(kind=rk), | private | :: | matrixA(nFields,nFields) | ||||

real(kind=rk), | private | :: | invA(nFields,nFields) |