Commit 9879aacb authored by henry's avatar henry
Browse files

Removed ddtPhiCorr as the standard approach does not work with

moving meshes
parent 975c2d3c
...@@ -57,15 +57,15 @@ int main(int argc, char *argv[]) ...@@ -57,15 +57,15 @@ int main(int argc, char *argv[])
# include "readControls.H" # include "readControls.H"
# include "CourantNo.H" # include "CourantNo.H"
// Make the fluxes absolute
fvc::makeAbsolute(phi, U);
# include "setDeltaT.H" # include "setDeltaT.H"
runTime++; runTime++;
Info<< "Time = " << runTime.timeName() << nl << endl; Info<< "Time = " << runTime.timeName() << nl << endl;
// Make the fluxes absolute
fvc::makeAbsolute(phi, U);
mesh.update(); mesh.update();
if (mesh.changing() && correctPhi) if (mesh.changing() && correctPhi)
...@@ -96,11 +96,6 @@ int main(int argc, char *argv[]) ...@@ -96,11 +96,6 @@ int main(int argc, char *argv[])
U = rAU*UEqn.H(); U = rAU*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf()); phi = (fvc::interpolate(U) & mesh.Sf());
if (ddtPhiCorr)
{
phi += fvc::ddtPhiCorr(rAU, U, phi);
}
if (p.needReference()) if (p.needReference())
{ {
fvc::makeRelative(phi, U); fvc::makeRelative(phi, U);
......
...@@ -67,15 +67,16 @@ int main(int argc, char *argv[]) ...@@ -67,15 +67,16 @@ int main(int argc, char *argv[])
{ {
#include "readControls.H" #include "readControls.H"
#include "CourantNo.H" #include "CourantNo.H"
// Make the fluxes absolute
fvc::makeAbsolute(phi, U);
#include "setDeltaT.H" #include "setDeltaT.H"
runTime++; runTime++;
Info<< "Time = " << runTime.timeName() << nl << endl; Info<< "Time = " << runTime.timeName() << nl << endl;
// Make the fluxes absolute
fvc::makeAbsolute(phi, U);
scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime(); scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
// Do any mesh changes // Do any mesh changes
...@@ -96,9 +97,6 @@ int main(int argc, char *argv[]) ...@@ -96,9 +97,6 @@ int main(int argc, char *argv[])
#include "correctPhi.H" #include "correctPhi.H"
} }
// Keep the absolute fluxes for use in ddtPhiCorr
surfaceScalarField phiAbs0("phiAbs0", phi);
// Make the fluxes relative to the mesh motion // Make the fluxes relative to the mesh motion
fvc::makeRelative(phi, U); fvc::makeRelative(phi, U);
......
...@@ -5,11 +5,6 @@ ...@@ -5,11 +5,6 @@
U = rAU*UEqn.H(); U = rAU*UEqn.H();
surfaceScalarField phiU("phiU", (fvc::interpolate(U) & mesh.Sf())); surfaceScalarField phiU("phiU", (fvc::interpolate(U) & mesh.Sf()));
if (ddtPhiCorr)
{
phiU += fvc::ddtPhiCorr(rAU, rho, U, phiAbs0);
}
phi = phiU + phi = phiU +
( (
fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma) fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
......
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