diff --git a/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C b/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C index 94e375304a..6d5c2f7eaa 100644 --- a/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C +++ b/src/postProcessing/functionObjects/utilities/DESModelRegions/DESModelRegions.C @@ -130,7 +130,7 @@ void Foam::DESModelRegions::write() Info<< type() << " output:" << nl; } - tmp result; + tmp tresult; label DESpresent = false; if (mesh.foundObject("turbulenceModel")) @@ -142,7 +142,7 @@ void Foam::DESModelRegions::write() { const icoDESModel& des = dynamic_cast(model); - result = des.LESRegion(); + tresult = des.LESRegion(); DESpresent = true; } } @@ -155,15 +155,17 @@ void Foam::DESModelRegions::write() { const cmpDESModel& des = dynamic_cast(model); - result = des.LESRegion(); + tresult = des.LESRegion(); DESpresent = true; } } if (DESpresent) { + const volScalarField& result = tresult(); + scalar prc = - gSum(result().internalField()*mesh.V())/gSum(mesh.V())*100.0; + gSum(result.internalField()*mesh.V())/gSum(mesh.V())*100.0; if (Pstream::master()) { @@ -174,9 +176,11 @@ void Foam::DESModelRegions::write() if (log_) { Info<< " LES = " << prc << " % (volume)" << nl - << " RES = " << 100.0 - prc << " % (volume)" << nl + << " RAS = " << 100.0 - prc << " % (volume)" << nl << endl; } + + result.write(); } else {