programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
itrgrv.f90 File Reference

This function adds the explicit part of the divergence of the mass flux due to the pressure gradient (routine analog to diften.f90). More...

Functions/Subroutines

subroutine itrgrv (init, inc, imrgra, iccocg, nswrgp, imligp, ircflp, iphydp, iwarnp, nfecra, epsrgp, climgp, extrap, frcxt, pvar, coefap, coefbp, cofafp, cofbfp, viscf, viscb, viscel, weighf, weighb, diverg)
 

Detailed Description

This function adds the explicit part of the divergence of the mass flux due to the pressure gradient (routine analog to diften.f90).

More precisely, the divergence of the mass flux side $ \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij $ is updated as follows:

\[ \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij = \sum_{\fij \in \Facei{\celli}} \dot{m}_\fij - \sum_{\fij \in \Facei{\celli}} \left( \tens{\mu}_\fij \gradv_\fij P \cdot \vect{S}_\ij \right) \]

Function/Subroutine Documentation

subroutine itrgrv ( integer  init,
integer  inc,
integer  imrgra,
integer  iccocg,
integer  nswrgp,
integer  imligp,
integer  ircflp,
integer  iphydp,
integer  iwarnp,
integer  nfecra,
double precision  epsrgp,
double precision  climgp,
double precision  extrap,
double precision, dimension(3,ncelet)  frcxt,
double precision, dimension(ncelet)  pvar,
double precision, dimension(nfabor)  coefap,
double precision, dimension(nfabor)  coefbp,
double precision, dimension(nfabor)  cofafp,
double precision, dimension(nfabor)  cofbfp,
double precision, dimension(nfac)  viscf,
double precision, dimension(nfabor)  viscb,
double precision, dimension(6,ncelet), target  viscel,
double precision, dimension(2,nfac)  weighf,
double precision, dimension(nfabor)  weighb,
double precision, dimension(ncelet)  diverg 
)
Parameters
[in]initindicator
  • 1 initialize the mass flux to 0
  • 0 otherwise
[in]incindicator
  • 0 when solving an increment
  • 1 otherwise
[in]imrgraindicator
  • 0 iterative gradient
  • 1 least square gradient
[in]iccocgindicator
  • 1 re-compute cocg matrix (for iterativ gradients)
  • 0 otherwise
[in]nswrgpnumber of reconstruction sweeps for the gradients
[in]imligpclipping gradient method
  • < 0 no clipping
  • = 0 thank to neighbooring gradients
  • = 1 thank to the mean gradient
[in]ircflpindicator
  • 1 flux reconstruction,
  • 0 otherwise
[in]iphydpindicator
  • 1 hydrostatic pressure taken into account
  • 0 otherwise
[in]iwarnpverbosity
[in]epsrgprelative precision for the gradient reconstruction
[in]climgpclipping coeffecient for the computation of the gradient
[in]extrapcoefficient for extrapolation of the gradient
[in]frcxtbody force creating the hydrostatic pressure
[in]pvarsolved variable (pressure)
[in]coefapboundary condition array for the variable (Explicit part)
[in]coefbpboundary condition array for the variable (Impplicit part)
[in]cofafpboundary condition array for the diffusion of the variable (Explicit part)
[in]cofbfpboundary condition array for the diffusion of the variable (Implicit part)
[in]viscf$ \mu_\fij \dfrac{S_\fij}{\ipf \jpf} $ at interior faces for the r.h.s.
[in]viscb$ \mu_\fib \dfrac{S_\fib}{\ipf \centf} $ at border faces for the r.h.s.
[in]viscelsymmetric cell tensor $ \tens{\mu}_\celli $
[in]weighfinternal face weight between cells i j in case of tensor diffusion
[in]weighbboundary face weight for cells i in case of tensor diffusion
[in,out]divergdivergence of the mass flux