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, iterns, ncesmp, nfbpcd, ncmast, icetsm, ifbpcd, ltmast, isostd, dt, vel, coefav, coefbv, coefa_dp, coefb_dp, smacel, spcond, svcond, frcxt, dfrcxt, tpucou, viscf, viscb, phi, tslagr)

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:

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.

Please refer to the resopv section of the theory guide for more informations.

Function/Subroutine Documentation

◆ resopv()

subroutine resopv	(	integer	nvar,
		integer	iterns,
		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(3 ,nfabor)	coefav,
		double precision, dimension(3,3,nfabor)	coefbv,
		double precision, dimension(nfabor)	coefa_dp,
		double precision, dimension(nfabor)	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(nfac)	viscf,
		double precision, dimension(nfabor)	viscb,
		double precision, dimension(ncelet)	phi,
		double precision, dimension(ncelet,*)	tslagr
	)

Parameters

[in]	nvar	total number of variables
[in]	iterns	Navier-Stokes iteration number
[in]	ncesmp	number of cells with mass source term
[in]	nfbpcd	number of faces with condensation source term
[in]	ncmast	number of cells with condensation source terms
[in]	icetsm	index of cells with mass source term
[in]	ifbpcd	index of faces with condensation source term
[in]	ltmast	index of cells with condensation source terms
[in]	isostd	indicator of standard outlet and index of the reference outlet face
[in]	dt	time step (per cell)
[in]	vel	velocity
[in]	coefav	boundary condition array for the variable (explicit part)
[in]	coefbv	boundary condition array for the variable (implicit part)
[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 $\Gamma^n$ )
[in]	spcond	variable value associated to the condensation source term (for ivar=ipr, spcond is the flow rate $\Gamma_{s,cond}^n$ )
[in]	svcond	variable value associated to the condensation source term (for ivar=ipr, svcond is the flow rate $\Gamma_{v, cond}^n$ )
[in]	frcxt	external forces making hydrostatic pressure
[in]	dfrcxt	variation of the external forces composing the hydrostatic pressure
[in]	tpucou	non scalar time step in case of velocity pressure coupling
[in]	viscf	visc*surface/dist aux faces internes
[in]	viscb	visc*surface/dist aux faces de bord
[in]	phi	potential to be solved (pressure increment)
[in]	tslagr	coupling term for the Lagrangian module

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ resopv()