From e1fb615f87e3633a830de5153b91258d2001a1ec Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Wed, 18 Sep 2019 09:22:35 -0700
Subject: [PATCH] BUG: Fixing psib to use products mixture.

Fix Gitlab issue #1423
---
 .../reactionThermo/psiuReactionThermo/heheuPsiThermo.C        | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C b/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
index 11b6bae0938..3de79763ccb 100644
--- a/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
+++ b/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2017 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2017-2019 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
                             | Copyright (C) 2011-2016 OpenFOAM Foundation
@@ -611,7 +611,7 @@ Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::psib() const
     forAll(psibCells, celli)
     {
         psibCells[celli] =
-            this->cellReactants(celli).psi(pCells[celli], TbCells[celli]);
+            this->cellProducts(celli).psi(pCells[celli], TbCells[celli]);
     }
 
     volScalarField::Boundary& psibBf = psib.boundaryFieldRef();
-- 
GitLab