src/alge/bilsca.f90 File Reference

| programmer's documentation

- Main Page
- Related Pages
- Modules
- Data Structures
- Files

- File List
- Globals

All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
bilsca.f90 File Reference

Wrapper to the function which adds the explicit part of the convection/diffusion terms of a transport equation of a scalar field . More...

## Functions/Subroutines
subroutine | bilsca (idtvar, ivar, iconvp, idiffp, nswrgp, imligp, ircflp, ischcp, isstpp, inc, imrgra, iccocg, ipp, iwarnp, imucpp, idftnp, blencp, epsrgp, climgp, extrap, relaxp, thetap, pvar, pvara, coefap, coefbp, cofafp, cofbfp, flumas, flumab, viscf, viscb, viscce, xcpp, weighf, weighb, icvflb, icvfli, smbrp)

## Detailed Description

Wrapper to the function which adds the explicit part of the convection/diffusion terms of a transport equation of a scalar field .

More precisely, the right hand side is updated as follows:

Warning:

- has already been initialized before calling bilsca!
- mind the sign minus

Options for the diffusive scheme:

- idftnp = 1: scalar diffusivity
- idftnp = 6: symmetric tensor diffusivity

Options for the convective scheme:

- blencp = 0: upwind scheme for the advection
- blencp = 1: no upwind scheme except in the slope test
- ischcp = 0: second order
- ischcp = 1: centred
- imucpp = 0: do not multiply the convective part by
- imucpp = 1: multiply the convective part by

## Function/Subroutine Documentation

subroutine bilsca | ( | integer | idtvar,
|  | integer | ivar,
|  | integer | iconvp,
|  | integer | idiffp,
|  | integer | nswrgp,
|  | integer | imligp,
|  | integer | ircflp,
|  | integer | ischcp,
|  | integer | isstpp,
|  | integer | inc,
|  | integer | imrgra,
|  | integer | iccocg,
|  | integer | ipp,
|  | integer | iwarnp,
|  | integer | imucpp,
|  | integer | idftnp,
|  | double precision | blencp,
|  | double precision | epsrgp,
|  | double precision | climgp,
|  | double precision | extrap,
|  | double precision | relaxp,
|  | double precision | thetap,
|  | double precision, dimension (ncelet) | pvar,
|  | double precision, dimension(ncelet) | pvara,
|  | double precision, dimension(nfabor) | coefap,
|  | double precision, dimension(nfabor) | coefbp,
|  | double precision, dimension(nfabor) | cofafp,
|  | double precision, dimension(nfabor) | cofbfp,
|  | double precision, dimension(nfac) | flumas,
|  | double precision, dimension(nfabor) | flumab,
|  | double precision, dimension (nfac) | viscf,
|  | double precision, dimension (nfabor) | viscb,
|  | double precision, dimension(*) | viscce,
|  | double precision, dimension(ncelet) | xcpp,
|  | double precision, dimension(2,nfac) | weighf,
|  | double precision, dimension(nfabor) | weighb,
|  | integer | icvflb,
|  | integer, dimension(nfabor) | icvfli,
|  | double precision, dimension(ncelet) | smbrp
| ) |  |

Parameters

  [in] | idtvar | indicator of the temporal scheme
  [in] | ivar | index of the current variable
  [in] | iconvp | indicator 1 convection, 0 sinon
  [in] | idiffp | indicator 1 diffusion, 0 sinon
  [in] | nswrgp | number of reconstruction sweeps for the gradients
  [in] | imligp | clipping gradient method < 0 no clipping = 0 thank to neighbooring gradients = 1 thank to the mean gradient
  [in] | ircflp | indicator 1 flux reconstruction, 0 otherwise
  [in] | ischcp | indicator 1 centred 0 2nd order
  [in] | isstpp | indicator 1 without slope test 0 with slope test
  [in] | inc | indicator 0 when solving an increment 1 otherwise
  [in] | imrgra | indicator 0 iterative gradient 1 least square gradient
  [in] | iccocg | indicator 1 re-compute cocg matrix (for iterativ gradients) 0 otherwise
  [in] | ipp* | index of the variable for post-processing
  [in] | iwarnp | verbosity
  [in] | imucpp | indicator 0 do not multiply the convectiv term by Cp 1 do multiply the convectiv term by Cp
  [in] | idftnp | indicator 1 scalar diffusivity 6 symmetric tensor diffusivity
  [in] | blencp | fraction of upwinding
  [in] | epsrgp | relative precision for the gradient reconstruction
  [in] | climgp | clipping coeffecient for the computation of the gradient
  [in] | extrap | coefficient for extrapolation of the gradient
  [in] | relaxp | coefficient of relaxation
  [in] | thetap | weightening 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] | pvar | solved variable (current time step)
  [in] | pvara | solved variable (previous time step)
  [in] | coefap | boundary condition array for the variable (Explicit part)
  [in] | coefbp | boundary condition array for the variable (Impplicit part)
  [in] | cofafp | boundary condition array for the diffusion of the variable (Explicit part)
  [in] | cofbfp | boundary condition array for the diffusion of the variable (Implicit part)
  [in] | flumas | mass flux at interior faces
  [in] | flumab | mass flux at boundary faces
  [in] | viscf | at interior faces for the r.h.s.
  [in] | viscb | at border faces for the r.h.s.
  [in] | viscce | symmetric cell tensor
  [in] | xcpp | array of specific heat (Cp)
  [in] | weighf | internal face weight between cells i j in case of tensor diffusion
  [in] | weighb | boundary face weight for cells i in case of tensor diffusion
  [in] | icvflb | global indicator of boundary convection flux 0 upwind scheme at all boundary faces 1 imposed flux at some boundary faces
  [in] | icvfli | boundary face indicator array of convection flux 0 upwind scheme 1 imposed flux
  [in,out] | smbrp | right hand side