programmer's documentation
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, nfbpcd, ncmast, icetsm, ifbpcd, ltmast, isostd, dt, vel, propce, coefav, coefbv, coefa_dp, coefb_dp, smacel, spcond, svcond, frcxt, dfrcxt, tpucou, trav, viscf, viscb, dpvar, 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  nfbpcd,
integer  ncmast,
integer, dimension(ncesmp)  icetsm,
integer, dimension(nfbpcd)  ifbpcd,
integer, dimension(ncelet)  ltmast,
integer, dimension(nfabor+1)  isostd,
double precision, dimension (1:ncelet), target  dt,
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(nfbpcd,nvar)  spcond,
double precision, dimension(ncelet,nvar)  svcond,
double precision, dimension(3,ncelet)  frcxt,
double precision, dimension(3,ncelet)  dfrcxt,
double precision, dimension (1:6,1:ncelet), target  tpucou,
double precision, dimension(3,ncelet)  trav,
double precision, dimension(nfac)  viscf,
double precision, dimension(ndimfb)  viscb,
double precision, dimension(ncelet)  dpvar,
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]nfbpcdnumber of faces with condensation source term
[in]ncmastnumber of cells with condensation source terms
[in]icetsmindex of cells with mass source term
[in]ifbpcdindex of faces with condensation source term
[in]ltmastindex of cells with condensation source terms
[in]isostdindicator of standard outlet and index of the reference outlet face
[in]dttime step (per cell)
[in]velvelocity
[in]propcephysical properties at cell centers
[in]coefavboundary condition array for the variable (explicit part)
[in]coefbvboundary condition array for the variable (implicit part)
[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]spcondvariable value associated to the condensation source term (for ivar=ipr, spcond is the flow rate $ \Gamma_{s,cond}^n $)
[in]svcondvariable value associated to the condensation source term (for ivar=ipr, svcond is the flow rate $ \Gamma_{v, cond}^n $)
[in]frcxtexternal forces making hydrostatic pressure
[in]dfrcxtvariation of the external forces composing the 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]dpvartableau de travail pour increment
[in]tslagrcoupling term for the Lagrangian module
[in]travatableau de travail pour couplage