Changed the constructor for pd to MUST_READ.
Calculate p from pd (note p is still read in by the termo package) and correct thermo.
This commit is contained in:
parent
31b83d862b
commit
3eea88674a
@ -67,10 +67,14 @@
|
|||||||
(
|
(
|
||||||
"pd",
|
"pd",
|
||||||
runTime.timeName(),
|
runTime.timeName(),
|
||||||
mesh
|
mesh,
|
||||||
|
IOobject::MUST_READ,
|
||||||
|
IOobject::AUTO_WRITE
|
||||||
),
|
),
|
||||||
p - rho*gh - pRef,
|
mesh
|
||||||
p.boundaryField().types()
|
|
||||||
);
|
);
|
||||||
|
|
||||||
|
p = pd + rho*gh + pRef;
|
||||||
|
thermo->correct();
|
||||||
|
|
||||||
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
|
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
|
||||||
|
@ -64,12 +64,16 @@
|
|||||||
(
|
(
|
||||||
"pd",
|
"pd",
|
||||||
runTime.timeName(),
|
runTime.timeName(),
|
||||||
mesh
|
mesh,
|
||||||
|
IOobject::MUST_READ,
|
||||||
|
IOobject::AUTO_WRITE
|
||||||
),
|
),
|
||||||
p - rho*gh - pRef,
|
mesh
|
||||||
p.boundaryField().types()
|
|
||||||
);
|
);
|
||||||
|
|
||||||
|
p = pd + rho*gh + pRef;
|
||||||
|
thermo->correct();
|
||||||
|
|
||||||
|
|
||||||
label pdRefCell = 0;
|
label pdRefCell = 0;
|
||||||
scalar pdRefValue = 0.0;
|
scalar pdRefValue = 0.0;
|
||||||
|
@ -65,12 +65,16 @@
|
|||||||
(
|
(
|
||||||
"pd",
|
"pd",
|
||||||
runTime.timeName(),
|
runTime.timeName(),
|
||||||
mesh
|
mesh,
|
||||||
|
IOobject::MUST_READ,
|
||||||
|
IOobject::AUTO_WRITE
|
||||||
),
|
),
|
||||||
p - rho*gh - pRef,
|
mesh
|
||||||
p.boundaryField().types()
|
|
||||||
);
|
);
|
||||||
|
|
||||||
|
p = pd + rho*gh + pRef;
|
||||||
|
thermo->correct();
|
||||||
|
|
||||||
|
|
||||||
label pdRefCell = 0;
|
label pdRefCell = 0;
|
||||||
scalar pdRefValue = 0.0;
|
scalar pdRefValue = 0.0;
|
||||||
|
@ -173,7 +173,7 @@
|
|||||||
|
|
||||||
Info<< " Updating p from pd\n" << endl;
|
Info<< " Updating p from pd\n" << endl;
|
||||||
thermof[i].p() == pdf[i] + rhof[i]*ghf[i] + pRef;
|
thermof[i].p() == pdf[i] + rhof[i]*ghf[i] + pRef;
|
||||||
|
thermof[i].correct();
|
||||||
|
|
||||||
initialMassf[i] = fvc::domainIntegrate(rhof[i]).value();
|
initialMassf[i] = fvc::domainIntegrate(rhof[i]).value();
|
||||||
}
|
}
|
||||||
|
@ -37,7 +37,7 @@
|
|||||||
mesh
|
mesh
|
||||||
);
|
);
|
||||||
|
|
||||||
# include "compressibleCreatePhi.H"
|
#include "compressibleCreatePhi.H"
|
||||||
|
|
||||||
|
|
||||||
Info<< "Creating turbulence model\n" << endl;
|
Info<< "Creating turbulence model\n" << endl;
|
||||||
@ -69,7 +69,7 @@
|
|||||||
mesh
|
mesh
|
||||||
);
|
);
|
||||||
|
|
||||||
//pd = p - rho*gh - pRef;
|
|
||||||
p = pd + rho*gh + pRef;
|
p = pd + rho*gh + pRef;
|
||||||
|
thermo->correct();
|
||||||
|
|
||||||
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
|
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
|
||||||
|
Loading…
Reference in New Issue
Block a user