Merge pull request deepmodeling#416 from pkuLmq/master

fix on coefficients of RK3
JX278 · Jan 4, 2024 · 6fc4848 · 6fc4848
2 parents 94c068d + d619453
commit 6fc4848
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git 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
@@ -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)
                         );