zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

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

Browse Source
This commit is contained in:
Henry Weller
2016-02-15 15:53:04 +00:00
parent 630c45d5ab
commit 51442e153e
2 changed files with 29 additions and 35 deletions
Download Patch File Download Diff File Expand all files Collapse all files

57
src/finiteVolume/cfdTools/general/constrainPressure/constrainPressure.C
View File

@ -43,48 +43,35 @@ void Foam::constrainPressure
)
{
const fvMesh& mesh = p.mesh();
volScalarField::GeometricBoundaryField& pBf = p.boundaryField();
bool hasFixedFluxBCs = false;
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();
forAll(pBf, patchi)
{
if (isA<fixedFluxPressureFvPatchScalarField>(pBf[patchi]))
{
hasFixedFluxBCs = true;
break;
}
}
if (hasFixedFluxBCs)
{
const surfaceScalarField::GeometricBoundaryField& phiHbyABf =
phiHbyA.boundaryField();
const typename RAUType::GeometricBoundaryField& rhorAUBf =
rhorAU.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]
)
/(magSfBf[patchi]*rhorAUBf[patchi])
);
}
phiHbyABf[patchi]
- rho.boundaryField()[patchi]
*MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
)
/(magSfBf[patchi]*rhorAUBf[patchi])
);
}
}
}

7
src/finiteVolume/cfdTools/general/constrainPressure/constrainPressure.H
View File

@ -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;
}
};