src/base/resopv.f90 File Reference

| programmer's documentation

- Main Page
- Related Pages
- Modules
- Data Structures
- Files

- File List
- Globals

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:

The mass flux is then updated as follows:

Remarks:

- an iterative process is used to solve the Poisson equation.
- if the coefficient arak is set to 1, the the Rhie & Chow filter is activated.

## 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] | nvar | total number of variables
  [in] | ncesmp | number of cells with mass source term
  [in] | icetsm | index of cells with mass source term
  [in] | isostd | indicator of standard outlet and index of the reference outlet face
  [in] | dt | time step (per cell)
  [in,out] | rtp,rtpa | calculated variables at cell centers (at current and previous time steps)
  [in] | propce | physical properties at cell centers
  [in] | coefa_dp | boundary conditions for the pressure increment
  [in] | coefb_dp | boundary conditions for the pressure increment
  [in] | smacel | variable value associated to the mass source term (for ivar=ipr, smacel is the mass flux )
  [in] | frcxt | external forces making hydrostatic pressure
  [in] | dfrcxt | variation of the external forces
  [in] | making | the hydrostatic pressure
  [in] | tpucou | non scalar time step in case of velocity pressure coupling
  [in] | trav | right hand side for the normalizing the residual
  [in] | viscf | visc*surface/dist aux faces internes
  [in] | viscb | visc*surface/dist aux faces de bord
  [in] | drtp | tableau de travail pour increment
  [in] | tslagr | coupling term for the Lagrangian module
  [in] | trava | tableau de travail pour couplage