Update the convective mass flux before the Navier Stokes equations (prediction and correction steps). More...

Functions/Subroutines
subroutine	predfl (nvar, nscal, ncesmp, icetsm, dt, rtp, rtpa, propce, propfa, propfb, smacel)

Detailed Description

Update the convective mass flux before the Navier Stokes equations (prediction and correction steps).

This function computes a potential $\varia$ solving the equation:

$D \left( \Delta t, \varia \right) = \divs \left( \rho \vect{u}^n\right) - \Gamma^n + \dfrac{\rho^n - \rho^{n-1}}{\Delta t}$

This potential is then used to update the mass flux as follows:

$\dot{m}^{n+\frac{1}{2}}_\ij = \dot{m}^{n}_\ij - \Delta t \grad_\fij \varia \cdot \vect{S}_\ij$

Function/Subroutine Documentation

subroutine predfl	(	integer	nvar,
		integer	nscal,
		integer	ncesmp,
		integer, dimension(ncesmp)	icetsm,
		double precision, dimension(ncelet)	dt,
		double precision, dimension(ncelet,*)	rtp,
		double precision, dimension(ncelet,*)	rtpa,
		double precision, dimension(ncelet,*)	propce,
		double precision, dimension(nfac,*)	propfa,
		double precision, dimension(ndimfb,*)	propfb,
		double precision, dimension(ncesmp,nvar)	smacel
	)

Parameters

[in]	nvar	total number of variables
[in]	nscal	total number of scalars
[in]	ncesmp	number of cells with mass source term
[in]	icetsm	index of cells with mass source term
[in]	dt	time step (per cell)
[in]	rtp,rtpa	calculated variables at cell centers (at current and previous time steps)
[in]	propce	physical properties at cell centers
[in,out]	propfa	physical properties at interior face centers
[in,out]	propfb	physical properties at boundary face centers
[in]	smacel	variable value associated to the mass source term (for ivar=ipr, smacel is the mass flux $\Gamma^n$ )

Here is the call graph for this function: