Code_Saturne
CFD tool
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Functions/Subroutines
inimav.f90 File Reference

This function adds $ \rho \vect{u} \cdot \vect{S}_\ij$ to the mass flux $ \dot{m}_\ij $. More...

Functions/Subroutines

subroutine inimav
 

Detailed Description

This function adds $ \rho \vect{u} \cdot \vect{S}_\ij$ to the mass flux $ \dot{m}_\ij $.

For the reconstruction, $ \gradt \left(\rho \vect{u} \right) $ is computed with the following approximated boundary conditions:

For the mass flux at the boundary we have:

\[ \dot{m}_\ib = \left[ \rho_\fib \vect{A}_u + \rho_\fib \tens{B}_u \vect{u} + \tens{B}_u \left(\gradt \vect{u} \cdot \vect{\centi \centip}\right)\right] \cdot \vect{S}_\ij \]

The last equation uses some approximations detailed in the theory guide.

Function/Subroutine Documentation

subroutine inimav ( )
Parameters
[in]nvartotal number of variables
[in]nscaltotal number of scalars
[in]ivarindex of the current variable
[in]itypflindicator (take rho into account or not)
  • 1 compute $f {u}{S} $ - 0 compute \$f \vect{u}\cdot\vect{S} $
[in]iflmb0the mass flux is set to 0 on walls and symmetries if = 1
[in]initthe mass flux is initialize to 0 if > 0
[in]incindicator
  • 0 solve an increment
  • 1 otherwise
[in]imrgraindicator
  • 0 iterative gradient
  • 1 least square gradient
[in]nswrgunumber of sweeps for the reconstruction of the gradients
[in]imliguclipping gradient method
  • < 0 no clipping
  • = 0 thank to neighbooring gradients
  • = 1 thank to the mean gradient
[in]iwarnuverbosity
[in]nfecraunit of the standard output file
[in]epsrgurelative precision for the gradient reconstruction
[in]climguclipping coeffecient for the computation of the gradient
[in]extraucoefficient for extrapolation of the gradient
[in]romcell density
[in]rombborder face density
[in]velvector variable
[in]coefavboundary condition array for the variable (Explicit part - vector array )
[in]coefbvboundary condition array for the variable (Impplicit part - 3x3 tensor array)
[in,out]flumasinterior mass flux $ \dot{m}_\fij $
[in,out]flumabborder mass flux $ \dot{m}_\fib $

Here is the call graph for this function: