From d1fcf79dee8adc32fae95cd220acd9f8de6da082 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Fri, 18 Aug 2017 17:25:23 +0900
Subject: [PATCH] ignore increase in oxygen/fuel mass fraction

---
 solvers_post/LagrangianCMCFoam/CMCintegration.H | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

diff --git a/solvers_post/LagrangianCMCFoam/CMCintegration.H b/solvers_post/LagrangianCMCFoam/CMCintegration.H
index cc12438..d1e7201 100644
--- a/solvers_post/LagrangianCMCFoam/CMCintegration.H
+++ b/solvers_post/LagrangianCMCFoam/CMCintegration.H
@@ -5,7 +5,7 @@
             {
                 scalarField Ytemp(Ysize, 0.0);
 
-                for(label i=0 ; i<Ysize ; i++) //mass fraction integration over an eta-space for species i
+                for(label i = 0 ; i < Ysize ; i++) //mass fraction integration over an eta-space for species i
                 {
                     for(label k=0 ; k<=group ; k++)
                     {
@@ -18,9 +18,19 @@
                     }
                 }
 
-                for(label i=0 ; i<Ysize ; i++)
+                Switch ignoreOFIncrease = true;
+                Switch oxygenIncrease =  Ytemp[o2index] > Y[o2index][celli];
+                Switch fuelIncrease = Ytemp[fuelindex] > Y[fuelindex][celli];
+                if( ignoreOFIncrease && (oxygenIncrease || fuelIncrease))
                 {
-                    Y[i][celli] = Ytemp[i];
+                    // Do not update Yi
+                }
+                else
+                {
+                    for(label i = 0 ; i < Ysize ; i++)
+                    {
+                        Y[i][celli] = Ytemp[i];
+                    }
                 }
             }
         }