My Project
programmer's documentation
Functions/Subroutines
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, 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)
 

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.

This subroutine solves a linear system with the gauss method.

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  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 
)
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]nfbpcdnumber of faces with condensation source terms
[in]ncmastnumber of cells with condensation source terms
[in]icepdcindex of cells with head loss
[in]icetsmindex of cells with mass source term
[in]ifbpcdindex of faces with condensation source terms
[in]ltmastindex of cells with condensation source terms
[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]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]spcondvariable value associated to the condensation source term (for ivar=ipr, spcond is the flow rate $ \Gamma_{s, cond}^n $)
[in]frcxtexternal forces making hydrostatic pressure
[in]travaworking array for the velocity-pressure coupling
[in]dfrcxtvariation of the external forces
[in]grdphdhydrostatic pressure gradient to handle the imbalance between the 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]viscfisame as viscf for increments
[in]viscbisame as viscb for increments
[in]secvifsecondary viscosity at interior faces
[in]secvibsecondary viscosity at boundary faces
[in]w1working array