Code_Saturne
CFD tool
Functions/Subroutines
predfl.f90 File Reference

Update the convective mass flux before the Navier Stokes equations (prediction and correction steps). More...

Functions/Subroutines

subroutine predfl (nvar, nscal, ncesmp, icetsm, dt, rtp, rtpa, propce, propfa, propfb, smacel)
 

Detailed Description

Update the convective mass flux before the Navier Stokes equations (prediction and correction steps).

This function computes a potential $ \varia $ solving the equation:

\[ D \left( \Delta t, \varia \right) = \divs \left( \rho \vect{u}^n\right) - \Gamma^n + \dfrac{\rho^n - \rho^{n-1}}{\Delta t} \]

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 \]

Function/Subroutine Documentation

◆ predfl()

subroutine predfl ( integer  nvar,
integer  nscal,
integer  ncesmp,
integer, dimension(ncesmp)  icetsm,
double precision, dimension(ncelet)  dt,
double precision, dimension(ncelet,*)  rtp,
double precision, dimension(ncelet,*)  rtpa,
double precision, dimension(ncelet,*)  propce,
double precision, dimension(nfac,*)  propfa,
double precision, dimension(ndimfb,*)  propfb,
double precision, dimension(ncesmp,nvar)  smacel 
)
Parameters
[in]nvartotal number of variables
[in]nscaltotal number of scalars
[in]ncesmpnumber of cells with mass source term
[in]icetsmindex of cells with mass source term
[in]dttime step (per cell)
[in]rtp,rtpacalculated variables at cell centers (at current and previous time steps)
[in]propcephysical properties at cell centers
[in,out]propfaphysical properties at interior face centers
[in,out]propfbphysical properties at boundary face centers
[in]smacelvariable value associated to the mass source term (for ivar=ipr, smacel is the mass flux $ \Gamma^n $)
Here is the call graph for this function: