From 9879aacb5a80b44e2e01224415b8a30500610272 Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Tue, 27 May 2008 22:51:14 +0100
Subject: [PATCH] Removed ddtPhiCorr as the standard approach does not work
with moving meshes
---
.../solvers/incompressible/icoDyMFoam/icoDyMFoam.C | 11 +++--------
.../solvers/multiphase/interDyMFoam/interDyMFoam.C | 10 ++++------
applications/solvers/multiphase/interDyMFoam/pEqn.H | 5 -----
3 files changed, 7 insertions(+), 19 deletions(-)
diff --git a/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
index 814e755326f..5b8d15c544e 100644
--- a/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
+++ b/applications/solvers/incompressible/icoDyMFoam/icoDyMFoam.C
@@ -57,15 +57,15 @@ int main(int argc, char *argv[])
# include "readControls.H"
# include "CourantNo.H"
+ // Make the fluxes absolute
+ fvc::makeAbsolute(phi, U);
+
# include "setDeltaT.H"
runTime++;
Info<< "Time = " << runTime.timeName() << nl << endl;
- // Make the fluxes absolute
- fvc::makeAbsolute(phi, U);
-
mesh.update();
if (mesh.changing() && correctPhi)
@@ -96,11 +96,6 @@ int main(int argc, char *argv[])
U = rAU*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf());
- if (ddtPhiCorr)
- {
- phi += fvc::ddtPhiCorr(rAU, U, phi);
- }
-
if (p.needReference())
{
fvc::makeRelative(phi, U);
diff --git a/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
index 9d6b9c4cd9e..053a9fa813a 100644
--- a/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interDyMFoam/interDyMFoam.C
@@ -67,15 +67,16 @@ int main(int argc, char *argv[])
{
#include "readControls.H"
#include "CourantNo.H"
+
+ // Make the fluxes absolute
+ fvc::makeAbsolute(phi, U);
+
#include "setDeltaT.H"
runTime++;
Info<< "Time = " << runTime.timeName() << nl << endl;
- // Make the fluxes absolute
- fvc::makeAbsolute(phi, U);
-
scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
// Do any mesh changes
@@ -96,9 +97,6 @@ int main(int argc, char *argv[])
#include "correctPhi.H"
}
- // Keep the absolute fluxes for use in ddtPhiCorr
- surfaceScalarField phiAbs0("phiAbs0", phi);
-
// Make the fluxes relative to the mesh motion
fvc::makeRelative(phi, U);
diff --git a/applications/solvers/multiphase/interDyMFoam/pEqn.H b/applications/solvers/multiphase/interDyMFoam/pEqn.H
index 0fdba99cd69..c6961f6c22c 100644
--- a/applications/solvers/multiphase/interDyMFoam/pEqn.H
+++ b/applications/solvers/multiphase/interDyMFoam/pEqn.H
@@ -5,11 +5,6 @@
U = rAU*UEqn.H();
surfaceScalarField phiU("phiU", (fvc::interpolate(U) & mesh.Sf()));
- if (ddtPhiCorr)
- {
- phiU += fvc::ddtPhiCorr(rAU, rho, U, phiAbs0);
- }
-
phi = phiU +
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
--
GitLab