From cdafb7e42be754c18ce8d6d157e1d3dba4aab042 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Fri, 14 Jun 2024 17:38:02 +0100
Subject: [PATCH] multiphaseEuler: Prevent zero blending factors when alpha is
 zero

---
 .../blendingMethods/blendingMethod/blendingMethod.C           | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/modules/multiphaseEuler/phaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C b/applications/modules/multiphaseEuler/phaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C
index b160685b3d..9e69edc527 100644
--- a/applications/modules/multiphaseEuler/phaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C
+++ b/applications/modules/multiphaseEuler/phaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2014-2023 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2014-2024 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -176,7 +176,7 @@ Foam::tmp<Foam::volScalarField> Foam::blendingMethod::x
     return
         systemSet == 0b00
       ? alpha(alphas, phaseSet, false)
-      : alpha(alphas, phaseSet, false)/alpha(alphas, systemSet, true);
+      : alpha(alphas, phaseSet, true)/alpha(alphas, systemSet, true);
 }