UEqnSimple.H

Go to the documentation of this file.

// Momentum predictor
 
MRF.correctBoundaryVelocity(U);
 
if (hasFvSource_)
{
    daFvSourcePtr_->calcFvSource(fvSource);
}
 
tmp<fvVectorMatrix> tUEqn(
    fvm::div(phi, U)
    + fvm::Sp(alphaPorosity, U)
    + MRF.DDt(U)
    + turbulencePtr_->divDevReff(U)
    - fvSource
    - fvOptions(U));
fvVectorMatrix& UEqn = tUEqn.ref();
 
UEqn.relax();
 
fvOptions.constrain(UEqn);
 
if (simple.momentumPredictor())
{
    // get the solver performance info such as initial
    // and final residuals
    SolverPerformance<vector> solverU = solve(UEqn == -fvc::grad(p));
 
    DAUtility::primalResidualControl(solverU, printToScreen_, "U", daGlobalVarPtr_->primalMaxRes);
 
    fvOptions.correct(U);
}
 
// bound U
DAUtility::boundVar(allOptions, U, printToScreen_);