Commit 82e59fa1 authored by Mark Olesen's avatar Mark Olesen
Browse files

bugfix for extra RASModel::yPlus() parameters

parent 5692a06a
......@@ -94,23 +94,26 @@ int main(int argc, char *argv[])
singlePhaseTransportModel laminarTransport(U, phi);
autoPtr<incompressible::RASModel> RASModel
autoPtr<incompressible::RASModel> rasModel
(
incompressible::RASModel::New(U, phi, laminarTransport)
);
const scalar Cmu =
rasModel->coeffDict().lookupOrDefault<scalar>("Cmu", 0.09);
const fvPatchList& patches = mesh.boundary();
forAll(patches, patchi)
forAll(patches, patchI)
{
const fvPatch& currPatch = patches[patchi];
const fvPatch& currPatch = patches[patchI];
if (typeid(currPatch) == typeid(wallFvPatch))
{
yPlus.boundaryField()[patchi] = RASModel->yPlus(patchi);
const scalarField& Yp = yPlus.boundaryField()[patchi];
yPlus.boundaryField()[patchI] = rasModel->yPlus(patchI, Cmu);
const scalarField& Yp = yPlus.boundaryField()[patchI];
Info<< "Patch " << patchi
Info<< "Patch " << patchI
<< " named " << currPatch.name()
<< " y+ : min: " << min(Yp) << " max: " << max(Yp)
<< " average: " << average(Yp) << nl << endl;
......
......@@ -184,7 +184,8 @@ tmp<scalarField> RASModel::yPlus(const label patchNo, const scalar Cmu) const
{
WarningIn
(
"tmp<scalarField> RASModel::yPlus(const label patchNo) const"
"tmp<scalarField> RASModel::yPlus"
"(const label patchNo, const scalar Cmu) const"
) << "Patch " << patchNo << " is not a wall. Returning null field"
<< nl << endl;
......
......@@ -176,7 +176,8 @@ tmp<scalarField> RASModel::yPlus(const label patchNo, const scalar Cmu) const
{
WarningIn
(
"tmp<scalarField> RASModel::yPlus(const label patchNo) const"
"tmp<scalarField> RASModel::yPlus"
"(const label patchNo, const scalar Cmu) const"
) << "Patch " << patchNo << " is not a wall. Returning null field"
<< nl << endl;
......
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