My Project
programmer's documentation
|
This subroutine performs 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, nfbpcd, ncmast, icepdc, icetsm, ifbpcd, ltmast, itypsm, dt, vel, vela, velk, tslagr, coefav, coefbv, cofafv, cofbfv, ckupdc, smacel, spcond, frcxt, grdphd, trava, dfrcxt, tpucou, trav, viscf, viscb, viscfi, viscbi, secvif, secvib, w1) |
This subroutine performs the velocity prediction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver.
This subroutine solves a linear system with the gauss method.
Please refer to the predvv section of the theory guide for more informations.
subroutine predvv | ( | integer | iappel, |
integer | nvar, | ||
integer | nscal, | ||
integer | iterns, | ||
integer | ncepdp, | ||
integer | ncesmp, | ||
integer | nfbpcd, | ||
integer | ncmast, | ||
integer, dimension(ncepdp) | icepdc, | ||
integer, dimension(ncesmp) | icetsm, | ||
integer, dimension(nfbpcd) | ifbpcd, | ||
integer, dimension(ncelet) | ltmast, | ||
integer, dimension(ncesmp,nvar) | itypsm, | ||
double precision, dimension(ncelet) | dt, | ||
double precision, dimension (3, ncelet) | vel, | ||
double precision, dimension (3, ncelet) | vela, | ||
double precision, dimension (3, ncelet) | velk, | ||
double precision, dimension(ncelet,*) | tslagr, | ||
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(6,ncepdp) | ckupdc, | ||
double precision, dimension(ncesmp,nvar) | smacel, | ||
double precision, dimension(nfbpcd,nvar) | spcond, | ||
double precision, dimension(3,ncelet) | frcxt, | ||
double precision, dimension(3, ncelet) | grdphd, | ||
double precision, dimension(ndim,ncelet) | trava, | ||
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 | ||
) |
[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] | nfbpcd | number of faces with condensation source terms |
[in] | ncmast | number of cells with condensation source terms |
[in] | icepdc | index of cells with head loss |
[in] | icetsm | index of cells with mass source term |
[in] | ifbpcd | index of faces with condensation source terms |
[in] | ltmast | index of cells with condensation source terms |
[in] | itypsm | type of mass source term for the variables |
[in] | dt | time step (per cell) |
[in] | vel | velocity |
[in] | vela | velocity at the previous time step |
[in] | velk | velocity at the previous sub iteration (or vela) |
[in] | tslagr | coupling term for the Lagrangian module |
[in] | coefav | boundary condition array for the variable (explicit part) |
[in] | coefbv | boundary condition array for the variable (implicit 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] | 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] | spcond | variable value associated to the condensation source term (for ivar=ipr, spcond is the flow rate ) |
[in] | frcxt | external forces making hydrostatic pressure |
[in] | trava | working array for the velocity-pressure coupling |
[in] | dfrcxt | variation of the external forces |
[in] | grdphd | hydrostatic pressure gradient to handle the imbalance 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 | same as viscf for increments |
[in] | viscbi | same as viscb for increments |
[in] | secvif | secondary viscosity at interior faces |
[in] | secvib | secondary viscosity at boundary faces |
[in] | w1 | working array |