|
My Project
programmer's documentation
|
Update the convective mass flux before the Navier Stokes equations (prediction and correction steps). More...
Functions/Subroutines | |
| subroutine | predfl (nvar, ncesmp, icetsm, dt, smacel) |
Update the convective mass flux before the Navier Stokes equations (prediction and correction steps).
This function computes a potential 
![\[
D \left( \Delta t, \varia \right) = \divs \left( \rho \vect{u}^n\right)
- \Gamma^n
+ \dfrac{\rho^n - \rho^{n-1}}{\Delta t}
\]](form_440.png)
This potential is then used to update the mass flux as follows:
![\[
\dot{m}^{n+\frac{1}{2}}_\ij = \dot{m}^{n}_\ij
- \Delta t \grad_\fij \varia \cdot \vect{S}_\ij
\]](form_441.png)
| subroutine predfl | ( | integer | nvar, |
| integer | ncesmp, | ||
| integer, dimension(ncesmp) | icetsm, | ||
| double precision, dimension(ncelet) | dt, | ||
| double precision, dimension(ncesmp,nvar) | smacel | ||
| ) |
| [in] | nvar | total number of variables |
| [in] | ncesmp | number of cells with mass source term |
| [in] | icetsm | index of cells with mass source term |
| [in] | dt | time step (per cell) |
| [in] | smacel | variable value associated to the mass source term (for ivar=ipr, smacel is the mass flux |