Skip to content
Snippets Groups Projects
Commit e156d196 authored by Henry's avatar Henry
Browse files

rhoSimpleFoam family: updated to use e rather than h

parent f0bb8c2c
Branches
Tags
No related merge requests found
......@@ -20,7 +20,7 @@
);
volScalarField& p = thermo.p();
volScalarField& h = thermo.h();
volScalarField& e = thermo.e();
const volScalarField& psi = thermo.psi();
Info<< "Reading field U\n" << endl;
......
{
volScalarField K("K", 0.5*magSqr(U));
fvScalarMatrix hEqn
(
fvm::div(phi, h)
- fvm::Sp(fvc::div(phi), h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
fvc::div(phi)*K - fvc::div(phi, K)
);
hEqn.relax();
hEqn.solve();
thermo.correct();
}
{
fvScalarMatrix hEqn
(
fvm::div(phi, h)
- fvm::Sp(fvc::div(phi), h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
- fvc::div(phi, 0.5*magSqr(U), "div(phi,K)")
);
pZones.addEnthalpySource(thermo, rho, hEqn);
hEqn.relax();
hEqn.solve();
thermo.correct();
}
......@@ -63,7 +63,7 @@ int main(int argc, char *argv[])
// Pressure-velocity SIMPLE corrector
{
#include "UEqn.H"
#include "hEqn.H"
#include "eEqn.H"
#include "pEqn.H"
}
......
......@@ -59,7 +59,7 @@ int main(int argc, char *argv[])
// Pressure-velocity SIMPLE corrector
{
#include "UEqn.H"
#include "hEqn.H"
#include "eEqn.H"
#include "pEqn.H"
}
......
......@@ -61,8 +61,8 @@ int main(int argc, char *argv[])
// Velocity-pressure-enthalpy SIMPLEC corrector
{
#include "UEqn.H"
#include "eEqn.H"
#include "pEqn.H"
#include "hEqn.H"
}
turbulence->correct();
......
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