From 5b4e84c97bf272104f432406b22fc8728d751da2 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Sun, 23 Feb 2020 22:17:54 +0000
Subject: [PATCH] reactingEulerFoam::MovingPhaseModel: Generalised support for
 derived fixedValue BCs using the assignable() member function

Resolves bug-report https://bugs.openfoam.org/view.php?id=3457
---
 .../phaseModel/MovingPhaseModel/MovingPhaseModel.C    | 11 +++--------
 1 file changed, 3 insertions(+), 8 deletions(-)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MovingPhaseModel/MovingPhaseModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MovingPhaseModel/MovingPhaseModel.C
index b22551d81e..edc13703b7 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MovingPhaseModel/MovingPhaseModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MovingPhaseModel/MovingPhaseModel.C
@@ -85,16 +85,11 @@ Foam::MovingPhaseModel<BasePhaseModel>::phi(const volVectorField& U) const
             calculatedFvPatchScalarField::typeName
         );
 
-        forAll(U.boundaryField(), i)
+        forAll(U.boundaryField(), patchi)
         {
-            if
-            (
-                isA<fixedValueFvPatchVectorField>(U.boundaryField()[i])
-             || isA<slipFvPatchVectorField>(U.boundaryField()[i])
-             || isA<partialSlipFvPatchVectorField>(U.boundaryField()[i])
-            )
+            if (!U.boundaryField()[patchi].assignable())
             {
-                phiTypes[i] = fixedValueFvPatchScalarField::typeName;
+                phiTypes[patchi] = fixedValueFvPatchScalarField::typeName;
             }
         }