From bc3a2d60c9415e729e0d0bed758f630ecea8e188 Mon Sep 17 00:00:00 2001
From: pkuLmq <1900011062@pku.edu.cn>
Date: Tue, 18 Oct 2022 14:49:55 +0800
Subject: [PATCH] fix bug for diffusion correction in dfHighSpeedFoam

---
 applications/solvers/dfHighSpeedFoam/rhoYEqn.H | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
index 2a19deec..cf9d0bc3 100644
--- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
+++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
@@ -3,6 +3,11 @@ if (!inviscid)
     hDiffCorrFlux = Zero;
     diffAlphaD = Zero;
     sumYDiffError = Zero;
+    
+    forAll(Y, i)
+    {
+        sumYDiffError += chemistry.rhoD(i)*fvc::grad(Y[i]);
+    }
 }
 
 tmp<fv::convectionScheme<scalar>> mvConvection
@@ -44,10 +49,6 @@ tmp<fv::convectionScheme<scalar>> mvConvection
 
             if (!inviscid)
             {
-                forAll(Y, i)
-                {
-                    sumYDiffError += chemistry.rhoD(i)*fvc::grad(Y[i]);
-                }
                 const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf();
 
                 hDiffCorrFlux += chemistry.hai(i)*(chemistry.rhoD(i)*fvc::grad(Yi) - Yi*sumYDiffError);