UEqnTurbo.H
Go to the documentation of this file.
1 // Solve the Momentum equation
2 
3 MRF.correctBoundaryVelocity(U);
4 
5 tmp<fvVectorMatrix> tUEqn(
6  fvm::div(phi, U)
7  + MRF.DDt(rho, U)
8  + turbulencePtr_->divDevRhoReff(U));
9 fvVectorMatrix& UEqn = tUEqn.ref();
10 
11 UEqn.relax();
12 
13 // get the solver performance info such as initial
14 // and final residuals
15 SolverPerformance<vector> solverU = solve(UEqn == -fvc::grad(p));
16 
17 DAUtility::primalResidualControl(solverU, printToScreen_, "U", daGlobalVarPtr_->primalMaxRes);
18 
19 // bound U
20 DAUtility::boundVar(allOptions, U, printToScreen_);
21 U.correctBoundaryConditions();
U
volVectorField & U
Definition: createRefsPimpleDyM.H:7
UEqn
fvVectorMatrix & UEqn
Definition: UEqnTurbo.H:9
solverU
SolverPerformance< vector > solverU
Definition: UEqnTurbo.H:15
solve
nuTilda1Eqn solve(solverDictNuTilda)
MRF
IOMRFZoneListDF & MRF
Definition: createRefsRhoSimple.H:18
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
tUEqn
tmp< fvVectorMatrix > tUEqn(fvm::div(phi, U)+MRF.DDt(rho, U)+turbulencePtr_->divDevRhoReff(U))