Files
openfoam/applications/solvers/combustion/PDRFoam/UEqn.H
2008-04-15 18:56:58 +01:00

17 lines
351 B
C

fvVectorMatrix UEqn
(
betav*fvm::ddt(rho, U)
+ fvm::div(phi, U)
+ turbulence->divDevRhoReff(U)
==
betav*rho*g
);
volSymmTensorField invA = inv(I*UEqn.A() + drag->Dcu());
if (momentumPredictor)
{
U = invA & (UEqn.H() - betav*fvc::grad(p));
U.correctBoundaryConditions();
}