90 lines
2.5 KiB
C++
90 lines
2.5 KiB
C++
|
{
|
||
|
|
labelList cellBoundaryFaceCount(epsilon.size(), 0);
|
||
|
|
|
||
|
|
scalar Cmu25 = ::pow(Cmu, 0.25);
|
||
|
|
scalar Cmu75 = ::pow(Cmu, 0.75);
|
||
|
|
|
||
|
|
const fvPatchList& patches = mesh.boundary();
|
||
|
|
|
||
|
|
//- Initialise the near-wall P field to zero
|
||
|
|
forAll(patches, patchi)
|
||
|
|
{
|
||
|
|
const fvPatch& curPatch = patches[patchi];
|
||
|
|
|
||
|
|
if (isType<wallFvPatch>(curPatch))
|
||
|
|
{
|
||
|
|
forAll(curPatch, facei)
|
||
|
|
{
|
||
|
|
label faceCelli = curPatch.faceCells()[facei];
|
||
|
|
|
||
|
|
epsilon[faceCelli] = 0.0;
|
||
|
|
G[faceCelli] = 0.0;
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
//- Accumulate the wall face contributions to epsilon and G
|
||
|
|
// Increment cellBoundaryFaceCount for each face for averaging
|
||
|
|
forAll(patches, patchi)
|
||
|
|
{
|
||
|
|
const fvPatch& curPatch = patches[patchi];
|
||
|
|
|
||
|
|
if (isType<wallFvPatch>(curPatch))
|
||
|
|
{
|
||
|
|
const scalarField& rhow = rho.boundaryField()[patchi];
|
||
|
|
|
||
|
|
const scalarField muw = mul.boundaryField()[patchi];
|
||
|
|
const scalarField& mutw = mut.boundaryField()[patchi];
|
||
|
|
|
||
|
|
scalarField magFaceGradU =
|
||
|
|
mag(U.boundaryField()[patchi].snGrad());
|
||
|
|
|
||
|
|
forAll(curPatch, facei)
|
||
|
|
{
|
||
|
|
label faceCelli = curPatch.faceCells()[facei];
|
||
|
|
|
||
|
|
scalar yPlus =
|
||
|
|
Cmu25*y[patchi][facei]*::sqrt(k[faceCelli])
|
||
|
|
/(muw[facei]/rhow[facei]);
|
||
|
|
|
||
|
|
// For corner cells (with two boundary or more faces),
|
||
|
|
// epsilon and G in the near-wall cell are calculated
|
||
|
|
// as an average
|
||
|
|
|
||
|
|
cellBoundaryFaceCount[faceCelli]++;
|
||
|
|
|
||
|
|
epsilon[faceCelli] +=
|
||
|
|
Cmu75*rho[faceCelli]*::pow(k[faceCelli], 1.5)
|
||
|
|
/(kappa*y[patchi][facei]);
|
||
|
|
|
||
|
|
if (yPlus > 11.6)
|
||
|
|
{
|
||
|
|
G[faceCelli] +=
|
||
|
|
mutw[facei]*magFaceGradU[facei]
|
||
|
|
*Cmu25*::sqrt(k[faceCelli])
|
||
|
|
/(kappa*y[patchi][facei]);
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
|
||
|
|
// perform the averaging
|
||
|
|
|
||
|
|
forAll(patches, patchi)
|
||
|
|
{
|
||
|
|
const fvPatch& curPatch = patches[patchi];
|
||
|
|
|
||
|
|
if (isType<wallFvPatch>(curPatch))
|
||
|
|
{
|
||
|
|
forAll(curPatch, facei)
|
||
|
|
{
|
||
|
|
label faceCelli = curPatch.faceCells()[facei];
|
||
|
|
|
||
|
|
epsilon[faceCelli] /= cellBoundaryFaceCount[faceCelli];
|
||
|
|
G[faceCelli] /= cellBoundaryFaceCount[faceCelli];
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|