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

This function adds the explicit part of the convection/diffusion terms of a transport equation of a vector field $ \vect{\varia} $. More...

Functions/Subroutines

subroutine bilsc4 (idtvar, ivar, iconvp, idiffp, nswrgp, imligp, ircflp, ischcp, isstpp, inc, imrgra, ivisep, ippu, iwarnp, blencp, epsrgp, climgp, relaxp, thetap, pvar, pvara, coefav, coefbv, cofafv, cofbfv, flumas, flumab, viscf, viscb, secvif, icvflb, icvfli, smbr)
 

Detailed Description

This function adds the explicit part of the convection/diffusion terms of a transport equation of a vector field $ \vect{\varia} $.

More precisely, the right hand side $ \vect{Rhs} $ is updated as follows:

\[ \vect{Rhs} = \vect{Rhs} - \sum_{\fij \in \Facei{\celli}} \left( \dot{m}_\ij \left( \vect{\varia}_\fij - \vect{\varia}_\celli \right) - \mu_\fij \gradt_\fij \vect{\varia} \cdot \vect{S}_\ij \right) \]

Remark: if ivisep = 1, then we also take $ \mu \transpose{\gradt\vect{\varia}} + \lambda \trace{\gradt\vect{\varia}} $, where $ \lambda $ is the secondary viscosity, i.e. usually $ -\frac{2}{3} \mu $.

Warning:

Options:

Function/Subroutine Documentation

subroutine bilsc4 ( integer  idtvar,
integer  ivar,
integer  iconvp,
integer  idiffp,
integer  nswrgp,
integer  imligp,
integer  ircflp,
integer  ischcp,
integer  isstpp,
integer  inc,
integer  imrgra,
integer  ivisep,
integer  ippu,
integer  iwarnp,
double precision  blencp,
double precision  epsrgp,
double precision  climgp,
double precision  relaxp,
double precision  thetap,
double precision, dimension (3 ,ncelet)  pvar,
double precision, dimension (3 ,ncelet)  pvara,
double precision, dimension(3 ,nfabor)  coefav,
double precision, dimension(3,3,nfabor)  coefbv,
double precision, dimension(3 ,nfabor)  cofafv,
double precision, dimension(3,3,nfabor)  cofbfv,
double precision, dimension(nfac)  flumas,
double precision, dimension(nfabor)  flumab,
double precision, dimension (nfac)  viscf,
double precision, dimension (nfabor)  viscb,
double precision, dimension(nfac)  secvif,
integer  icvflb,
integer, dimension(nfabor)  icvfli,
double precision, dimension(3,ncelet)  smbr 
)
Parameters
[in]idtvarindicator of the temporal scheme
[in]ivarindex of the current variable
[in]iconvpindicator
  • 1 convection,
  • 0 sinon
[in]idiffpindicator
  • 1 diffusion,
  • 0 sinon
[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]ischcpindicator
  • 1 centred
  • 0 2nd order
[in]isstppindicator
  • 1 without slope test
  • 0 with slope test
[in]incindicator
  • 0 when solving an increment
  • 1 otherwise
[in]imrgraindicator
  • 0 iterative gradient
  • 1 least square gradient
[in]ivisepindicator to take $ \divv \left(\mu \gradt \transpose{\vect{a}} \right) -2/3 \grad\left( \mu \dive \vect{a} \right)$
  • 1 take into account,
  • 0 otherwise
[in]ippuindex of the variable for post-processing
[in]iwarnpverbosity
[in]blencpfraction of upwinding
[in]epsrgprelative precision for the gradient reconstruction
[in]climgpclipping coeffecient for the computation of the gradient
[in]relaxpcoefficient of relaxation
[in]thetapweightening coefficient for the theta-schema,
  • thetap = 0: explicit scheme
  • thetap = 0.5: time-centred scheme (mix between Crank-Nicolson and Adams-Bashforth)
  • thetap = 1: implicit scheme
[in]pvarvitesse resolue (instant courant)
[in]pvaravitesse resolue (instant precedent)
[in]coefavboundary condition array for the variable (Explicit part)
[in]coefbvboundary condition array for the variable (Impplicit part)
[in]cofafvboundary condition array for the diffusion of the variable (Explicit part)
[in]cofbfvboundary condition array for the diffusion of the variable (Implicit part)
[in]flumasmass flux at interior faces
[in]flumabmass flux at boundary faces
[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]secvifsecondary viscosity at interior faces
[in]icvflbglobal indicator of boundary convection flux
  • 0 upwind scheme at all boundary faces
  • 1 imposed flux at some boundary faces
[in]icvfliboundary face indicator array of convection flux
  • 0 upwind scheme
  • 1 imposed flux
[in,out]smbrright hand side $ \vect{Rhs} $