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) | :: | 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 | :: | ifieldDia | ||||
| integer, | private | :: | ifieldNonDia | ||||
| integer, | private | :: | icomp | ||||
| real(kind=rk), | private | :: | matrixA(nFields,nFields) | ||||
| real(kind=rk), | private | :: | invA(nFields,nFields) |