search

staneuski / dualFuelEngine

OpenFOAM solver based on phenomenological compression model for dual-fuel marine engines
GNU Affero General Public License v3.0
13 stars 3 forks source link

Unstable solution for alphaX fields #15

Closed staneuski closed 5 years ago

staneuski commented 5 years ago
Environment Parameter
Compiled using OpenFOAM v7
Solver Version v0.3-alpha
Commit d21650f
Case \ File RiemannTube/
Comments Check the concentrations plot below (right bottom corner)

alphas-0

staneuski commented 5 years ago

Commit 2eed466 close this issue #15, not the #5 (there was a typo in the commit message).

Switch div(phi,alphaX) to upwind scheme to remove instability.

If upwind is set only for alphaAir:

diff a/tutorials/multiCompression/RiemannTube/system/fvSchemes b/tutorials/multiCompression/RiemannTube/system/fvSchemes

@@ -58 +58 @@     
-  div(phi,alphaAir) Gauss cubic; // alphaAirEqn
+  div(phi,alphaAir) Gauss upwind; // alphaAirEqn

alphas_alphaAir-upwind

If upwind is set only for alphaAir and alphaGas:

diff a/tutorials/multiCompression/RiemannTube/system/fvSchemes b/tutorials/multiCompression/RiemannTube/system/fvSchemes

@@ -58 +58 @@     
-  div(phi,alphaAir) Gauss cubic; // alphaAirEqn
+  div(phi,alphaAir) Gauss upwind; // alphaAirEqn

@@ -59 +59 @@     
-  div(phi,alphaGas) Gauss cubic; // alphaGasEqn
+  div(phi,alphaGas) Gauss upwind; // alphaGasEqn

alphas_all-upwind

vl-zenkin commented 5 years ago

Можно все-таки какое-нибудь словесное резюме?