7.2
general documentation
predvv.f90 File Reference

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, icepdc, icetsm, itypsm, dt, vel, vela, velk, da_uu, tslagr, coefav, coefbv, cofafv, cofbfv, ckupdc, smacel, frcxt, trava, dfrcxt, tpucou, trav, viscf, viscb, viscfi, viscbi, secvif, secvib)
 

Detailed Description

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

  • at the first call, the predicted velocities are computed and also an estimator on the predicted velocity is computed.
  • at the second call, a global estimator on Navier Stokes is computed. This second call is done after the correction step (cs_pressure_correction).

Please refer to the predvv section of the theory guide for more informations.

Function/Subroutine Documentation

◆ predvv()

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 (3, ncelet)  vel,
double precision, dimension (3, ncelet)  vela,
double precision, dimension (3, ncelet)  velk,
double precision, dimension (6, ncelet)  da_uu,
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(3,ncelet)  frcxt,
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 
)
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]velvelocity
[in]velavelocity at the previous time step
[in]velkvelocity at the previous sub iteration (or vela)
[in,out]da_uuvelocity matrix
[in]tslagrcoupling term for the Lagrangian module
[in]coefavboundary condition array for the variable (explicit part)
[in]coefbvboundary condition array for the variable (implicit part)
[in]cofafvboundary condition array for the diffusion of the variable (explicit part)
[in]cofbfvboundary condition array for the diffusion of the variable (implicit part)
[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
[out]dfrcxtvariation of the external forces
[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]viscfisame as viscf for increments
[in]viscbisame as viscb for increments
[in]secvifsecondary viscosity at interior faces
[in]secvibsecondary viscosity at boundary faces