Commit e86258bc authored by andy's avatar andy
Browse files

ENH: Updated sources for compressible solvers

parent a0eeb3f0
......@@ -21,8 +21,12 @@
);
EEqn.relax();
sources.constrain(EEqn);
EEqn.solve();
sources.correct(he);
thermo.correct();
}
......@@ -16,5 +16,7 @@ sources.constrain(UEqn());
if (pimple.momentumPredictor())
{
solve(UEqn() == -fvc::grad(p));
sources.correct(U);
K = 0.5*magSqr(U);
}
......@@ -15,8 +15,12 @@
);
EEqn.relax();
sources.constrain(EEqn);
EEqn.solve();
sources.correct(he);
thermo.correct();
}
......@@ -10,4 +10,8 @@
UEqn().relax();
sources.constrain(UEqn());
solve(UEqn() == -fvc::grad(p));
sources.correct(U);
......@@ -15,8 +15,12 @@
);
EEqn.relax();
sources.constrain(EEqn);
EEqn.solve();
sources.correct(he);
thermo.correct();
}
......@@ -33,6 +33,8 @@
U = trTU() & ((UEqn() == sources(rho, U))().H() - gradp);
}
U.correctBoundaryConditions();
sources.correct(U);
}
else
{
......@@ -42,6 +44,8 @@
solve(UEqn() == -fvc::grad(p) + sources(rho, U));
sources.correct(U);
trAU = 1.0/UEqn().A();
trAU().rename("rAU");
}
Supports Markdown
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