![]() |
programmer's documentation
|
This subroutine perform the velocity prediction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver. More...
Functions/Subroutines | |
subroutine | predvv (iappel, nvar, nscal, iterns, ncepdp, ncesmp, icepdc, icetsm, itypsm, dt, rtpa, vel, vela, propce, flumas, flumab, tslagr, coefav, coefbv, cofafv, cofbfv, ckupdc, smacel, frcxt, grdphd, trava, ximpa, uvwk, dfrcxt, tpucou, trav, viscf, viscb, viscfi, viscbi, secvif, secvib, w1, w7, w8, w9, xnormp) |
This subroutine perform the velocity prediction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver.
subroutine predvv | ( | integer | iappel, |
integer | nvar, | ||
integer | nscal, | ||
integer | iterns, | ||
integer | ncepdp, | ||
integer | ncesmp, | ||
integer, dimension(ncepdp) | icepdc, | ||
integer, dimension(ncesmp) | icetsm, | ||
integer, dimension(ncesmp,nvar) | itypsm, | ||
double precision, dimension(ncelet) | dt, | ||
double precision, dimension(ncelet,*) | rtpa, | ||
double precision, dimension (3 ,ncelet) | vel, | ||
double precision, dimension (3 ,ncelet) | vela, | ||
double precision, dimension(ncelet,*) | propce, | ||
double precision, dimension(nfac) | flumas, | ||
double precision, dimension(nfabor) | flumab, | ||
double precision, dimension(ncelet,*) | tslagr, | ||
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(ncepdp,6) | ckupdc, | ||
double precision, dimension(ncesmp,nvar) | smacel, | ||
double precision, dimension(3,ncelet) | frcxt, | ||
double precision, dimension(ncelet,3) | grdphd, | ||
double precision, dimension(ndim,ncelet) | trava, | ||
double precision, dimension(ndim,ndim,ncelet) | ximpa, | ||
double precision, dimension(ndim,ncelet) | uvwk, | ||
double precision, dimension(3,ncelet) | dfrcxt, | ||
double precision, dimension(6, ncelet) | tpucou, | ||
double precision, dimension(3,ncelet) | trav, | ||
double precision, dimension(*) | viscf, | ||
double precision, dimension(nfabor) | viscb, | ||
double precision, dimension(*) | viscfi, | ||
double precision, dimension(nfabor) | viscbi, | ||
double precision, dimension(nfac) | secvif, | ||
double precision, dimension(nfabor) | secvib, | ||
double precision, dimension(ncelet) | w1, | ||
double precision, dimension(ncelet) | w7, | ||
double precision, dimension(ncelet) | w8, | ||
double precision, dimension(ncelet) | w9, | ||
double precision, dimension(ncelet) | xnormp | ||
) |
[in] | iappel | call number (1 or 2) |
[in] | nvar | total number of variables |
[in] | nscal | total number of scalars |
[in] | iterns | index of the iteration on Navier-Stokes |
[in] | ncepdp | number of cells with head loss |
[in] | ncesmp | number of cells with mass source term |
[in] | icepdc | index of cells with head loss |
[in] | icetsm | index of cells with mass source term |
[in] | itypsm | type of mass source term for the variables |
[in] | dt | time step (per cell) |
[in,out] | rtpa | calculated variables at cell centers (at previous time step) |
[in] | propce | physical properties at cell centers |
[in] | flumas | internal mass flux (depending on iappel) |
[in] | flumab | boundary mass flux (depending on iappel) |
[in] | ckupdc | work array for the head loss |
[in] | smacel | variable value associated to the mass source term (for ivar=ipr, smacel is the mass flux ![]() |
[in] | frcxt | external forces making hydrostatic pressure |
[in] | trava | working array for the velocity-pressure coupling |
[in] | ximpa | idem |
[in] | uvwk | idem (stores the velocity at the previous iteration)* |
[in] | dfrcxt | variation of the external forces |
[in] | making | the hydrostatic pressure |
[in] | grdphd | hydrostatic pressure gradient to handle the imbalance |
[in] | between | the pressure gradient and gravity source term |
[in] | tpucou | non scalar time step in case of velocity pressure coupling |
[in] | trav | right hand side for the normalizing the residual |
[in] | viscf | visc*surface/dist aux faces internes |
[in] | viscb | visc*surface/dist aux faces de bord |
[in] | viscfi | idem viscf pour increments |
[in] | viscbi | idem viscb pour increments |
[in] | secvif | secondary viscosity at interior faces |
[in] | secvib | secondary viscosity at boundary faces |
[in] | w1 | working array |
[in] | w7 | working array |
[in] | w8 | working array |
[in] | w9 | working array |
[in] | xnormp | workig array for the norm of the pressure |