programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
predvv.f90 File Reference

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)
 

Detailed Description

This subroutine perform the velocity prediction step of the Navier Stokes equations for incompressible or slightly compressible flows for the coupled velocity components solver.

Function/Subroutine Documentation

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 
)
Parameters
[in]iappelcall number (1 or 2)
[in]nvartotal number of variables
[in]nscaltotal number of scalars
[in]iternsindex of the iteration on Navier-Stokes
[in]ncepdpnumber of cells with head loss
[in]ncesmpnumber of cells with mass source term
[in]icepdcindex of cells with head loss
[in]icetsmindex of cells with mass source term
[in]itypsmtype of mass source term for the variables
[in]dttime step (per cell)
[in,out]rtpacalculated variables at cell centers (at previous time step)
[in]propcephysical properties at cell centers
[in]flumasinternal mass flux (depending on iappel)
[in]flumabboundary mass flux (depending on iappel)
[in]ckupdcwork array for the head loss
[in]smacelvariable value associated to the mass source term (for ivar=ipr, smacel is the mass flux $ \Gamma^n $)
[in]frcxtexternal forces making hydrostatic pressure
[in]travaworking array for the velocity-pressure coupling
[in]ximpaidem
[in]uvwkidem (stores the velocity at the previous iteration)*
[in]dfrcxtvariation of the external forces
[in]makingthe hydrostatic pressure
[in]grdphdhydrostatic pressure gradient to handle the imbalance
[in]betweenthe pressure gradient and gravity source term
[in]tpucounon scalar time step in case of velocity pressure coupling
[in]travright hand side for the normalizing the residual
[in]viscfvisc*surface/dist aux faces internes
[in]viscbvisc*surface/dist aux faces de bord
[in]viscfiidem viscf pour increments
[in]viscbiidem viscb pour increments
[in]secvifsecondary viscosity at interior faces
[in]secvibsecondary viscosity at boundary faces
[in]w1working array
[in]w7working array
[in]w8working array
[in]w9working array
[in]xnormpworkig array for the norm of the pressure