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");
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)
GeometricField< vector, fvPatchField, volMesh > volVectorField
tmp< volTensorField > trTU
static const Identity< scalar > I
Foam::fvConstraints & fvConstraints
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevTau(U)==fvModels.source(rho, U))
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
fvMatrix< vector > fvVectorMatrix
Switch pressureImplicitPorosity(false)
Tensor< scalar > tensor
Tensor of scalars.
tmp< volScalarField > trAU
IOporosityModelList pZones(mesh)