From 40db2359b1dc843040442b9f16cd10ab33672ffe Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 12 Sep 2013 17:02:02 +0100
Subject: [PATCH] interDyMFoam: Update surface fields following mesh change
---
.../solvers/multiphase/interFoam/alphaEqnSubCycle.H | 12 +++++++++++-
.../multiphase/interFoam/interDyMFoam/interDyMFoam.C | 2 ++
2 files changed, 13 insertions(+), 1 deletion(-)
diff --git a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
index 428876fd229..ac40a9bc9e1 100644
--- a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
@@ -3,7 +3,17 @@
if (nAlphaSubCycles > 1)
{
dimensionedScalar totalDeltaT = runTime.deltaT();
- surfaceScalarField rhoPhiSum(0.0*rhoPhi);
+ surfaceScalarField rhoPhiSum
+ (
+ IOobject
+ (
+ "rhoPhiSum",
+ runTime.timeName(),
+ mesh
+ ),
+ mesh,
+ dimensionedScalar("0", dimDensity, 0)
+ );
for
(
diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
index 4c986ef2d8a..032f2db7303 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@@ -100,6 +100,8 @@ int main(int argc, char *argv[])
// Make the flux relative to the mesh motion
fvc::makeRelative(phi, U);
+
+ interface.correct();
}
if (mesh.changing() && checkMeshCourantNo)
--
GitLab