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