From a6f17edeeb8a1d48bb2035a99b2a5b70acd2ede9 Mon Sep 17 00:00:00 2001 From: andy Date: Fri, 19 Jul 2013 12:00:18 +0100 Subject: [PATCH] BUG: Corrected transportModel handling in Peclet function object --- .../functionObjects/utilities/Peclet/Peclet.C | 27 ++++++++++++++----- 1 file changed, 21 insertions(+), 6 deletions(-) diff --git a/src/postProcessing/functionObjects/utilities/Peclet/Peclet.C b/src/postProcessing/functionObjects/utilities/Peclet/Peclet.C index 885e611fa0..681e36fc9e 100644 --- a/src/postProcessing/functionObjects/utilities/Peclet/Peclet.C +++ b/src/postProcessing/functionObjects/utilities/Peclet/Peclet.C @@ -129,7 +129,6 @@ void Foam::Peclet::end() // Do nothing - only valid on write } - void Foam::Peclet::timeSet() { // Do nothing - only valid on write @@ -163,16 +162,32 @@ void Foam::Peclet::write() nuEff = model.nuEff(); } - else if (mesh.foundObject("transportProperties")) + else if (mesh.foundObject("transportProperties")) { - const transportModel& model = - mesh.lookupObject("transportProperties"); + const dictionary& model = + mesh.lookupObject("transportProperties"); - nuEff = model.nu(); + nuEff = + tmp + ( + new volScalarField + ( + IOobject + ( + "nuEff", + mesh.time().timeName(), + mesh, + IOobject::NO_READ, + IOobject::NO_WRITE + ), + mesh, + dimensionedScalar(model.lookup("nu")) + ) + ); } else { - FatalErrorIn("void Foam::wallShearStress::write()") + FatalErrorIn("void Foam::Peclet::write()") << "Unable to determine the viscosity" << exit(FatalError); }