programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
resopv.f90 File Reference

This subroutine performs 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, ncesmp, icetsm, isostd, dt, rtp, rtpa, vel, propce, coefav, coefbv, coefa_dp, coefb_dp, smacel, frcxt, dfrcxt, tpucou, trav, viscf, viscb, drtp, tslagr, trava)
 

Detailed Description

This subroutine performs 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

subroutine resopv ( integer  nvar,
integer  ncesmp,
integer, dimension(ncesmp)  icetsm,
integer, dimension(nfabor+1)  isostd,
double precision, dimension(ncelet)  dt,
double precision, dimension(ncelet,*)  rtp,
double precision, dimension(ncelet,*)  rtpa,
double precision, dimension (3 ,ncelet)  vel,
double precision, dimension(ncelet,*)  propce,
double precision, dimension(3 ,ndimfb)  coefav,
double precision, dimension(3,3,ndimfb)  coefbv,
double precision, dimension(ndimfb)  coefa_dp,
double precision, dimension(ndimfb)  coefb_dp,
double precision, dimension(ncesmp,nvar)  smacel,
double precision, dimension(3,ncelet)  frcxt,
double precision, dimension(3,ncelet)  dfrcxt,
double precision, dimension(6, ncelet)  tpucou,
double precision, dimension(3,ncelet)  trav,
double precision, dimension(nfac)  viscf,
double precision, dimension(ndimfb)  viscb,
double precision, dimension(ncelet)  drtp,
double precision, dimension(ncelet,*)  tslagr,
double precision, dimension(ndim,ncelet)  trava 
)
Parameters
[in]nvartotal number of variables
[in]ncesmpnumber of cells with mass source term
[in]icetsmindex of cells with mass source term
[in]isostdindicator of standard outlet and index of the reference outlet face
[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]coefa_dpboundary conditions for the pressure increment
[in]coefb_dpboundary conditions for the pressure increment
[in]smacelvariable value associated to the mass source term (for ivar=ipr, smacel is the mass flux $ \Gamma^n $)
[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]drtptableau de travail pour increment
[in]tslagrcoupling term for the Lagrangian module
[in]travatableau de travail pour couplage