Code_Saturne
CFD tool
Functions/Subroutines
resopv.f90 File Reference

This subroutine perform the pressure correction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver. More...

Functions/Subroutines

subroutine resopv (nvar, nscal, ncepdp, ncesmp, icepdc, icetsm, itypsm, isostd, idtsca, dt, rtp, rtpa, vel, vela, propce, propfa, propfb, coefa, coefb, coefav, coefbv, coefap, ckupdc, smacel, frcxt, dfrcxt, tpucou, trav, viscf, viscb, viscfi, viscbi, drtp, tslagr, frchy, dfrchy, trava)
 

Detailed Description

This subroutine perform the pressure correction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver.

This function solves the following Poisson equation on the pressure:

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

The mass flux is then updated as follows:

\[ \dot{m}^{n+1}_\ij = \dot{m}^{n}_\ij - \Delta t \grad_\fij \delta p \cdot \vect{S}_\ij \]

Remarks:

Function/Subroutine Documentation

◆ resopv()

subroutine resopv ( integer  nvar,
integer  nscal,
integer  ncepdp,
integer  ncesmp,
integer, dimension(ncepdp)  icepdc,
integer, dimension(ncesmp)  icetsm,
integer, dimension(ncesmp,nvar)  itypsm,
integer, dimension(nfabor+1)  isostd,
integer  idtsca,
double precision, dimension(ncelet)  dt,
double precision, dimension(ncelet,*)  rtp,
double precision, dimension(ncelet,*)  rtpa,
double precision, dimension (3 ,ncelet)  vel,
double precision, dimension (3 ,ncelet)  vela,
double precision, dimension(ncelet,*)  propce,
double precision, dimension(nfac,*)  propfa,
double precision, dimension(ndimfb,*)  propfb,
double precision, dimension(ndimfb,*)  coefa,
double precision, dimension(ndimfb,*)  coefb,
double precision, dimension(3 ,ndimfb)  coefav,
double precision, dimension(3,3,ndimfb)  coefbv,
double precision, dimension(nfabor)  coefap,
double precision, dimension(ncepdp,6)  ckupdc,
double precision, dimension(ncesmp,nvar)  smacel,
double precision, dimension(ncelet,3)  frcxt,
double precision, dimension(ncelet,3)  dfrcxt,
double precision, dimension(ndim,ncelet)  tpucou,
double precision, dimension(3,ncelet)  trav,
double precision, dimension(nfac)  viscf,
double precision, dimension(nfabor)  viscb,
double precision, dimension(nfac)  viscfi,
double precision, dimension(nfabor)  viscbi,
double precision, dimension(ncelet)  drtp,
double precision, dimension(ncelet,*)  tslagr,
double precision, dimension(ncelet,ndim)  frchy,
double precision, dimension(ncelet,ndim)  dfrchy,
double precision, dimension(ndim,ncelet)  trava 
)
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 terms
[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 $)
[in]nvartotal number of variables
[in]nscaltotal number of scalars
[in]ncepdpnumber of cells with head loss
[in]ncesmpnumber of cells with mass source term
[in]icepdcindex of cells with head loss
[in]icetsmindex of cells with mass source term
[in]itypsmtype of mass source term for the variables
[in]isostdindicator of standard outlet and index of the reference outlet face
[in]idtscaindicator of non scalar time step
[in]dttime step (per cell)
[in,out]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]coefa,coefbboundary conditions
[in]ckupdcwork array for the head loss
[in]smacelvalues of variables associated to the mass source (for the pressure, smacel is the mass flux)
[in]frcxtexternal forces making hydrostatic pressure
[in]dfrcxtvariation of the external forces
[in]makingthe hydrostatic pressure
[in]tpucounon scalar time step in case of velocity pressure coupling
[in]travright hand side for the normalizing the residual
[in]viscfvisc*surface/dist aux faces internes
[in]viscbvisc*surface/dist aux faces de bord
[in]viscfiidem viscf pour increments
[in]viscbiidem viscb pour increments
[in]drtptableau de travail pour increment
[in]tslagrcoupling term for teh Lagrangian module
[in]frchytableau de travail
[in]dfrchytableau de travail
[in]travatableau de travail pour couplage
Here is the call graph for this function: