From 77bb4721ae15319ee923200787ff0162e73a7b7d Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 19 Aug 2013 11:49:26 +0100
Subject: [PATCH] twoPhaseEulerFoam: Update the fluxes from the velocities for
changes in fixedValue conditions
---
applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H b/applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H
index e9211def7f4..073fe092d0a 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H
@@ -8,6 +8,10 @@
surfaceScalarField rAlphaAU1f(fvc::interpolate(alpha1*rAU1));
surfaceScalarField rAlphaAU2f(fvc::interpolate(alpha2*rAU2));
+ // Update the phi BCs from U before p BCs are updated
+ phi1.boundaryField() == (mesh.Sf().boundaryField() & U1.boundaryField());
+ phi2.boundaryField() == (mesh.Sf().boundaryField() & U2.boundaryField());
+
volVectorField HbyA1
(
IOobject::groupName("HbyA", phase1.name()),
--
GitLab