Commit 4bb7f08e authored by Henry Weller's avatar Henry Weller
Browse files

constrainPressure: Updated to use the more efficient patch-based MRF::relative function

parent 618bf48a
......@@ -43,50 +43,37 @@ void Foam::constrainPressure
)
{
const fvMesh& mesh = p.mesh();
volScalarField::GeometricBoundaryField& pBf = p.boundaryField();
bool hasFixedFluxBCs = false;
forAll(pBf, patchi)
{
if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi]))
{
hasFixedFluxBCs = true;
break;
}
}
volScalarField::GeometricBoundaryField& pBf = p.boundaryField();
if (hasFixedFluxBCs)
{
const volVectorField::GeometricBoundaryField& UBf = U.boundaryField();
const surfaceScalarField::GeometricBoundaryField& phiHbyABf =
phiHbyA.boundaryField();
const typename RAUType::GeometricBoundaryField& rhorAUBf =
rhorAU.boundaryField();
const surfaceVectorField::GeometricBoundaryField& SfBf =
mesh.Sf().boundaryField();
const surfaceScalarField::GeometricBoundaryField& magSfBf =
mesh.magSf().boundaryField();
// Pre-compute tho relative flux for all patches: currently MRFZone does
// not support computing the relative flux for individual patches
FieldField<fvsPatchField, scalar> phiRelBf
(
MRF.relative(mesh.Sf().boundaryField() & U.boundaryField())
);
forAll(pBf, patchi)
{
if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi]))
{
refCast<fixedFluxPressureFvPatchScalarField>(pBf[patchi])
.updateCoeffs
refCast<fixedFluxPressureFvPatchScalarField>
(
pBf[patchi]
).updateCoeffs
(
(
phiHbyABf[patchi]
- rho.boundaryField()[patchi]*phiRelBf[patchi]
- rho.boundaryField()[patchi]
*MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
)
/(magSfBf[patchi]*rhorAUBf[patchi])
);
}
}
}
}
......
......@@ -68,6 +68,13 @@ public:
{
return U;
}
//- Return the argument unchanged
template<class Type>
inline const Type& relative(const Type& U, const label patchi) const
{
return U;
}
};
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment