src/base/coditv.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
coditv.f90 File Reference

This function solves an advection diffusion equation with source terms for one time step for the vector variable . More...

## Functions/Subroutines
subroutine | coditv (idtvar, ivar, iconvp, idiffp, ireslp, ndircp, nitmap, imrgra, nswrsp, nswrgp, imligp, ircflp, ivisep, ischcp, isstpp, iescap, idftnp, iswdyp, imgrp, ncymxp, nitmfp, ippu, ippv, ippw, iwarnp, blencp, epsilp, epsrsp, epsrgp, climgp, relaxp, thetap, pvara, pvark, coefav, coefbv, cofafv, cofbfv, flumas, flumab, viscfm, viscbm, viscfs, viscbs, secvif, secvib, icvflb, icvfli, fimp, smbrp, pvar, eswork)

## Detailed Description

This function solves an advection diffusion equation with source terms for one time step for the vector variable .

The equation reads:

This equation is rewritten as:

where and

It is in fact solved with the following iterative process:

where

Be careful, it is forbidden to modify here!

## Function/Subroutine Documentation

subroutine coditv | ( | integer | idtvar,
|  | integer | ivar,
|  | integer | iconvp,
|  | integer | idiffp,
|  | integer | ireslp,
|  | integer | ndircp,
|  | integer | nitmap,
|  | integer | imrgra,
|  | integer | nswrsp,
|  | integer | nswrgp,
|  | integer | imligp,
|  | integer | ircflp,
|  | integer | ivisep,
|  | integer | ischcp,
|  | integer | isstpp,
|  | integer | iescap,
|  | integer | idftnp,
|  | integer | iswdyp,
|  | integer | imgrp,
|  | integer | ncymxp,
|  | integer | nitmfp,
|  | integer | ippu,
|  | integer | ippv,
|  | integer | ippw,
|  | integer | iwarnp,
|  | double precision | blencp,
|  | double precision | epsilp,
|  | double precision | epsrsp,
|  | double precision | epsrgp,
|  | double precision | climgp,
|  | double precision | relaxp,
|  | double precision | thetap,
|  | double precision, dimension(3,ncelet) | pvara,
|  | double precision, dimension(3,ncelet) | pvark,
|  | double precision, dimension(3,ndimfb) | coefav,
|  | double precision, dimension(3,3,ndimfb) | coefbv,
|  | double precision, dimension(3,ndimfb) | cofafv,
|  | double precision, dimension(3,3,ndimfb) | cofbfv,
|  | double precision, dimension(nfac) | flumas,
|  | double precision, dimension(nfabor) | flumab,
|  | double precision, dimension(*) | viscfm,
|  | double precision, dimension(nfabor) | viscbm,
|  | double precision, dimension(*) | viscfs,
|  | double precision, dimension(nfabor) | viscbs,
|  | double precision, dimension(nfac) | secvif,
|  | double precision, dimension(nfabor) | secvib,
|  | integer | icvflb,
|  | integer, dimension(nfabor) | icvfli,
|  | double precision, dimension(3,3,ncelet) | fimp,
|  | double precision, dimension(3,ncelet) | smbrp,
|  | double precision, dimension(3,ncelet) | pvar,
|  | double precision, dimension(3,ncelet) | eswork
| ) |  |

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] | ireslp | indicator 0 conjugate gradient 1 jacobi 2 bi-cgstab
  [in] | ndircp | indicator (0 if the diagonal is stepped aside)
  [in] | nitmap | maximum number of iteration to solve the iterative process
  [in] | imrgra | indicateur 0 iterative gradient 1 least square gradient
  [in] | nswrsp | number of reconstruction sweeps for the Right Hand Side
  [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] | ivisep | indicator to take 1 take into account, 0 otherwise
  [in] | ischcp | indicator 1 centred 0 2nd order
  [in] | isstpp | indicator 1 without slope test 0 with slope test
  [in] | iescap | compute the predictor indicator if 1
  [in] | idftnp | indicator 1 the diffusivity is scalar 6 the diffusivity is a symmetric tensor
  [in] | iswdyp | indicator 0 no dynamic relaxation 1 dynamic relaxation depending on 2 dynamic relaxation depending on and
  [in] | imgrp | indicator 0 no multi-grid 1 otherwise
  [in] | ncymxp | max. number of multigrid cycles
  [in] | nitmfp | number of equivalent iterations on fine mesh
  [in] | ippu | index of the variable for post-processing
  [in] | ippv | index of the variable for post-processing
  [in] | ippw | index of the variable for post-processing
  [in] | iwarnp | verbosity
  [in] | blencp | fraction of upwinding
  [in] | epsilp | precision pour resol iter
  [in] | epsrsp | relative precision for the iterative process
  [in] | epsrgp | relative precision for the gradient reconstruction
  [in] | climgp | clipping coeffecient for the computation 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] | pvara | variable at the previous time step
  [in] | pvark | variable at the previous sub-iteration . If you sub-iter on Navier-Stokes, then it allows to initialize by something else than pvara (usually pvar=pvara)
  [in] | coefav | boundary condition array for the variable (Explicit part)
  [in] | coefbv | boundary condition array for the variable (Impplicit part)
  [in] | cofafv | boundary condition array for the diffusion of the variable (Explicit part)
  [in] | cofbfv | 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] | viscfm | at interior faces for the matrix
  [in] | viscbm | at border faces for the matrix
  [in] | viscfs | at interior faces for the r.h.s.
  [in] | viscbs | at border faces for the r.h.s.
  [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] | fimp |
  [in] | smbrp | Right hand side
  [in,out] | pvar | current variable
  [out] | eswork | prediction-stage error estimator (if iescap > 0)