3 MRF.correctBoundaryVelocity(
U);
5 tmp<fvVectorMatrix>
tUEqn 26 tmp<volTensorField> tTU =
tensor(
I)*UEqn.A();
27 pZones.addResistance(UEqn, tTU.ref());
29 trTU.ref().rename(
"rAU");
35 for (
int UCorr=0; UCorr<
nUCorr; UCorr++)
37 U =
trTU() & (UEqn.H() - gradp);
39 U.correctBoundaryConditions();
45 pZones.addResistance(UEqn);
54 trAU.ref().rename(
"rAU");
autoPtr< compressible::turbulenceModel > turbulence
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
solve(UEqn==-fvc::grad(p))
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
GeometricField< vector, fvPatchField, volMesh > volVectorField
tmp< volTensorField > trTU
static const Identity< scalar > I
fvMatrix< vector > fvVectorMatrix
Switch pressureImplicitPorosity(false)
Tensor< scalar > tensor
Tensor of scalars.
tmp< volScalarField > trAU
IOporosityModelList pZones(mesh)