Commit d812879f authored by Henry Weller's avatar Henry Weller
Browse files

TurbulenceModels: updated the mapping of k to R BCs

Now k BCs which do not have a symmTensor equivalent are converted to "calculated"
parent 7321a5ee
......@@ -84,6 +84,23 @@ Foam::eddyViscosity<BasicTurbulenceModel>::R() const
{
tmp<volScalarField> tk(k());
// Get list of patchField type names from k
wordList patchFieldTypes(tk().boundaryField().types());
// For k patchField types which do not have an equivalent for symmTensor
// set to calculated
forAll(patchFieldTypes, i)
{
if
(
!fvPatchField<symmTensor>::patchConstructorTablePtr_
->found(patchFieldTypes[i])
)
{
patchFieldTypes[i] = calculatedFvPatchField<symmTensor>::typeName;
}
}
return tmp<volSymmTensorField>
(
new volSymmTensorField
......@@ -98,7 +115,7 @@ Foam::eddyViscosity<BasicTurbulenceModel>::R() const
false
),
((2.0/3.0)*I)*tk() - (nut_)*dev(twoSymm(fvc::grad(this->U_))),
tk().boundaryField().types()
patchFieldTypes
)
);
}
......
Supports Markdown
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