Add kinetic energy field K correction (v0.4.2-alpha)

staneuski · Apr 14, 2020 · 61c2137 · 61c2137
1 parent c22ef54
commit 61c2137
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 2 deletions.
diff --git a/solvers/multiCompressionFoam/createFieldRefs.H b/solvers/multiCompressionFoam/createFieldRefs.H
@@ -1,3 +1,4 @@
 volScalarField& T = thermo.T();
+volScalarField& e = thermo.he();
 const volScalarField& psi = thermo.psi(); // 1/R/T - compressibility [s^2/m^2]
 const volScalarField& mu = thermo.mu();
diff --git a/solvers/multiCompressionFoam/createFields.H b/solvers/multiCompressionFoam/createFields.H
@@ -7,8 +7,6 @@ autoPtr<fluidThermo> pThermo
 fluidThermo& thermo = pThermo();
 thermo.validate(args.executable(), "h", "e");
 
-volScalarField& e = thermo.he(); /*FIXME thermo.he() creates negative enengy*/
-
 volScalarField& p = thermo.p(); // Ref
 
 Info<< "Reading velocity field U" << endl;

diff --git a/solvers/multiCompressionFoam/multiCompressionFoam.C b/solvers/multiCompressionFoam/multiCompressionFoam.C
@@ -124,7 +124,9 @@ int main(int argc, char *argv[])
 
             UEqn.relax();
             UEqn.solve();
+
             fvOptions.correct(U);
+            K = 0.5*magSqr(U);
 
             fvScalarMatrix EEqn
             (