DAFoam
v4.0.2
Discrete Adjoint with OpenFOAM
dafoam
src
adjoint
DASolver
DATurboFoam
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))
Generated by
1.8.17