From d619453287b79e11107d3cff648fd25efbfabae0 Mon Sep 17 00:00:00 2001
From: pkuLmq <1900011062@pku.edu.cn>
Date: Thu, 4 Jan 2024 10:30:54 +0800
Subject: [PATCH] fix on coefficients of RK3

---
 applications/solvers/dfHighSpeedFoam/createFieldsSave.H | 2 +-
 applications/solvers/dfHighSpeedFoam/rhoYEqn.H          | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/dfHighSpeedFoam/createFieldsSave.H b/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
index c3fd5cd4b..37723319c 100644
--- a/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
+++ b/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
@@ -41,6 +41,6 @@ else if(ddtSchemes == "RK3SSP")
 {
     rkcoe1[0]=1.0; rkcoe2[0]=0.0; rkcoe3[0]=1.0;
     rkcoe1[1]=0.75; rkcoe2[1]=0.25; rkcoe3[1]=0.25;
-    rkcoe1[2]=0.33; rkcoe2[2]=0.66; rkcoe3[2]=0.66;
+    rkcoe1[2]=1.0/3.0; rkcoe2[2]=2.0/3.0; rkcoe3[2]=2.0/3.0;
     rk=3;
 }
diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
index fac11306c..ba634503c 100644
--- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
+++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
@@ -110,7 +110,7 @@ time_monitor_corrDiff += double(end - start) / double(CLOCKS_PER_SEC);
                     rhoYi_rhs = 
                         (
                             turbName == "laminar"
-                            ?  (fvc::laplacian(DEff(), Yi) - fvc::div(phiUc,Yi,"div(phic,Yi)")) 
+                            ?  (fvc::laplacian(DEff(), Yi) - fvc::div(phiUc,Yi,"div(phi,Yi_h)")) 
                             :  fvc::laplacian(DEff(), Yi)
                         );