UEqnRhoSimple.H
Go to the documentation of this file.
1 // Solve the Momentum equation
2 
3 MRF.correctBoundaryVelocity(U);
4 
5 if (hasFvSource_)
6 {
7  daFvSourcePtr_->calcFvSource(fvSource);
8 }
9 
10 tmp<fvVectorMatrix> tUEqn(
11  fvm::div(phi, U)
12  + MRF.DDt(rho, U)
13  + turbulencePtr_->divDevRhoReff(U)
14  - fvSource
15  - fvOptions(rho, U));
16 fvVectorMatrix& UEqn = tUEqn.ref();
17 
18 UEqn.relax();
19 
20 fvOptions.constrain(UEqn);
21 
22 // get the solver performance info such as initial
23 // and final residuals
24 SolverPerformance<vector> solverU = solve(UEqn == -fvc::grad(p));
25 
26 DAUtility::primalResidualControl(solverU, printToScreen_, "U", daGlobalVarPtr_->primalMaxRes);
27 
28 fvOptions.correct(U);
29 
30 // bound U
31 DAUtility::boundVar(allOptions, U, printToScreen_);
U
volVectorField & U
Definition: createRefsPimpleDyM.H:7
solve
nuTilda1Eqn solve(solverDictNuTilda)
MRF
IOMRFZoneListDF & MRF
Definition: createRefsRhoSimple.H:18
fvSource
volScalarField & fvSource
Definition: createRefsHeatTransfer.H:7
p
volScalarField & p
Definition: createRefsPimpleDyM.H:6
phi
surfaceScalarField & phi
Definition: createRefsPimpleDyM.H:8
rho
volScalarField & rho
Definition: createRefsRhoPimple.H:8
allOptions
const dictionary & allOptions
Definition: createRefsPimpleDyM.H:14
solverU
SolverPerformance< vector > solverU
Definition: UEqnRhoSimple.H:24
UEqn
fvVectorMatrix & UEqn
Definition: UEqnRhoSimple.H:16
tUEqn
tmp< fvVectorMatrix > tUEqn(fvm::div(phi, U)+MRF.DDt(rho, U)+turbulencePtr_->divDevRhoReff(U) - fvSource - fvOptions(rho, U))