From 8b71babfdd3a479206d3b3f5edf9bb473e77cc4b Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Thu, 18 Feb 2016 21:26:25 +0000
Subject: [PATCH] interFoam, interMixingFoam, interPhaseChangeFoam: construct
 rho with calculated BCs

Avoids problems with inherited complex BCs for which the controlling
fields and parameters are not initialized.
---
 applications/solvers/multiphase/interFoam/createFields.H       | 3 +--
 .../multiphase/interFoam/interMixingFoam/createFields.H        | 3 +--
 .../solvers/multiphase/interPhaseChangeFoam/createFields.H     | 3 +--
 3 files changed, 3 insertions(+), 6 deletions(-)

diff --git a/applications/solvers/multiphase/interFoam/createFields.H b/applications/solvers/multiphase/interFoam/createFields.H
index c3e65f24e3..e4891edba5 100644
--- a/applications/solvers/multiphase/interFoam/createFields.H
+++ b/applications/solvers/multiphase/interFoam/createFields.H
@@ -49,8 +49,7 @@ volScalarField rho
         mesh,
         IOobject::READ_IF_PRESENT
     ),
-    alpha1*rho1 + alpha2*rho2,
-    alpha1.boundaryField().types()
+    alpha1*rho1 + alpha2*rho2
 );
 rho.oldTime();
 
diff --git a/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H b/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
index 59180b8723..4d9f0e7f01 100644
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
@@ -50,8 +50,7 @@ volScalarField rho
         mesh,
         IOobject::READ_IF_PRESENT
     ),
-    alpha1*rho1 + alpha2*rho2 + alpha3*rho3,
-    alpha1.boundaryField().types()
+    alpha1*rho1 + alpha2*rho2 + alpha3*rho3
 );
 rho.oldTime();
 
diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H b/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H
index dc0491bd03..bb336f16c0 100644
--- a/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H
@@ -51,8 +51,7 @@ volScalarField rho
         mesh,
         IOobject::READ_IF_PRESENT
     ),
-    alpha1*rho1 + alpha2*rho2,
-    alpha1.boundaryField().types()
+    alpha1*rho1 + alpha2*rho2
 );
 rho.oldTime();