Code_Saturne
CFD tool
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 (nvar, nscal, ivar, itypfl, iflmb0, init, inc, imrgra, nswrgu, imligu, iwarnu, nfecra, epsrgu, climgu, extrau, rom, romb, vel, coefav, coefbv, flumas, flumab)
 

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

◆ inimav()

subroutine inimav ( integer  nvar,
integer  nscal,
integer  ivar,
integer  itypfl,
integer  iflmb0,
integer  init,
integer  inc,
integer  imrgra,
integer  nswrgu,
integer  imligu,
integer  iwarnu,
integer  nfecra,
double precision  epsrgu,
double precision  climgu,
double precision  extrau,
double precision, dimension(ncelet)  rom,
double precision, dimension(nfabor)  romb,
double precision, dimension(3,ncelet)  vel,
double precision, dimension(3,ndimfb)  coefav,
double precision, dimension(3,3,nfabor)  coefbv,
double precision, dimension(nfac)  flumas,
double precision, dimension(nfabor)  flumab 
)
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: