exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

fixing yPlusRAS for two phase

Browse source
This commit is contained in:
Henrik Rusche 2012-01-12 16:09:37 +01:00
parent 568e341348
commit 83e48c1bd7
1 changed files with 19 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

19
applications/utilities/postProcessing/wall/yPlusRAS/yPlusRAS.C
View file

@ -199,6 +199,24 @@ void calcTwoPhaseYPlus
# include "createPhi.H"
IOobject rhoHeader
(
"gamma",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
);
if (!rhoHeader.headerOk())
{
Info<< " no gamma field" << endl;
return;
}
Info << "Reading field gamma\n" << endl;
volScalarField gamma(rhoHeader, mesh);
Info<< "Reading transportProperties\n" << endl;
twoPhaseMixture twoPhaseProperties(U, phi, "gamma");
@ -246,6 +264,7 @@ int main(int argc, char *argv[])
#include "addRegionOption.H"
argList::validOptions.insert("compressible","");
argList::validOptions.insert("twoPhase","");
#include "setRootCase.H"
#include "createTime.H"