![]() |
programmer's documentation
|
This function solves an advection diffusion equation with source terms for one time step for the symmetric tensor variable .
More...
Functions/Subroutines | |
subroutine | coditts (idtvar, ivar, iconvp, idiffp, ndircp, imrgra, nswrsp, nswrgp, imligp, ircflp, ischcp, isstpp, idftnp, iswdyp, iwarnp, blencp, epsilp, epsrsp, epsrgp, climgp, relaxp, thetap, pvara, pvark, coefats, coefbts, cofafts, cofbfts, flumas, flumab, viscfm, viscbm, viscfs, viscbs, visccs, weighf, weighb, icvflb, icvfli, fimp, smbrp, pvar) |
This function solves an advection diffusion equation with source terms for one time step for the symmetric tensor 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!
subroutine coditts | ( | integer | idtvar, |
integer | ivar, | ||
integer | iconvp, | ||
integer | idiffp, | ||
integer | ndircp, | ||
integer | imrgra, | ||
integer | nswrsp, | ||
integer | nswrgp, | ||
integer | imligp, | ||
integer | ircflp, | ||
integer | ischcp, | ||
integer | isstpp, | ||
integer | idftnp, | ||
integer | iswdyp, | ||
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(6,ncelet) | pvara, | ||
double precision, dimension(6,ncelet) | pvark, | ||
double precision, dimension(6,ndimfb) | coefats, | ||
double precision, dimension(6,6,ndimfb) | coefbts, | ||
double precision, dimension(6,ndimfb) | cofafts, | ||
double precision, dimension(6,6,ndimfb) | cofbfts, | ||
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(6,ncelet) | visccs, | ||
double precision, dimension(2,nfac) | weighf, | ||
double precision, dimension(nfabor) | weighb, | ||
integer | icvflb, | ||
integer, dimension(nfabor) | icvfli, | ||
double precision, dimension(6,6,ncelet) | fimp, | ||
double precision, dimension(6,ncelet) | smbrp, | ||
double precision, dimension(6,ncelet) | pvar | ||
) |
[in] | idtvar | indicator of the temporal scheme |
[in] | ivar | index of the current variable |
[in] | iconvp | indicator
|
[in] | idiffp | indicator
|
[in] | ndircp | indicator (0 if the diagonal is stepped aside) |
[in] | imrgra | indicateur
|
[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
|
[in] | ircflp | indicator
|
[in] | ischcp | indicator
|
[in] | isstpp | indicator
|
[in] | idftnp | indicator
|
[in] | iswdyp | indicator
|
[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 coefficient for the computation of the gradient |
[in] | relaxp | coefficient of relaxation |
[in] | thetap | weighting coefficient for the theta-schema,
|
[in] | pvara | variable at the previous time step ![]() |
[in] | pvark | variable at the previous sub-iteration ![]() |
[in] | coefats | boundary condition array for the variable (Explicit part) |
[in] | coefbts | boundary condition array for the variable (Impplicit part) |
[in] | cofafts | boundary condition array for the diffusion of the variable (Explicit part) |
[in] | cofbfts | 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 | ![]() |
[in] | viscbm | ![]() |
[in] | viscfs | ![]() |
[in] | viscbs | ![]() |
[in] | visccs | symmetric cell tensor ![]() |
[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
|
[in] | icvfli | boundary face indicator array of convection flux
|
[in] | fimp | ![]() |
[in] | smbrp | Right hand side ![]() |
[in,out] | pvar | current variable |