37 template<
class RhoType,
class RAUType,
class MRFType>
44 const RAUType& rhorAU,
48 const fvMesh&
mesh =
p.mesh();
50 volScalarField::Boundary& pBf =
p.boundaryFieldRef();
52 const volVectorField::Boundary& UBf =
U.boundaryField();
53 const surfaceScalarField::Boundary& phiHbyABf =
55 const typename RAUType::Boundary& rhorAUBf = rhorAU.boundaryField();
56 const surfaceVectorField::Boundary& SfBf =
mesh.Sf().boundaryField();
57 const surfaceScalarField::Boundary& magSfBf =
58 mesh.magSf().boundaryField();
62 typedef updateablePatchTypes::updateableSnGrad snGradType;
63 const auto* snGradPtr = isA<snGradType>(pBf[patchi]);
67 const_cast<snGradType&
>(*snGradPtr).updateSnGrad
71 -
rho.boundaryField()[patchi]
72 *
MRF.relative(SfBf[patchi] & UBf[patchi], patchi)
74 /(magSfBf[patchi]*rhorAUBf[patchi])
81 template<
class RAUType>
95 template<
class RAUType,
class MRFType>
109 template<
class RAUType>
GeometricField< vector, fvPatchField, volMesh > volVectorField
tmp< volScalarField > rAU
#define forAll(list, i)
Loop across all elements in list.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
void constrainPressure(volScalarField &p, const RhoType &rho, const volVectorField &U, const surfaceScalarField &phiHbyA, const RAUType &rhorAU, const MRFType &MRF)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField