求助，最近基于blastReactingFoam进行气体爆炸计算，但是得到的爆炸超压震荡非常厉害，与实际显然有出入，感觉可能是离散格式设置不合理，这一块找不到什么资料进行参考修改。各位大佬可否给予指点。
基本设置如下
求解器：blastReactingFoam，湍流模型：k-omega SST
网格：二维，全正交网格，网格尺寸0.5 mm
6130e168-4c33-46c8-8b4d-a4e814099a85-image.png
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fluxScheme Kurganov;
ddtSchemes
{
default CrankNicolson 1;//Euler;
timeIntegrator RK2SSP;
}
gradSchemes
{
default Gauss linear;
grad(U) cellLimited leastSquares 1;
grad(k) cellLimited leastSquares 1;
grad(epsilon) cellLimited leastSquares 1;
grad(omega) cellLimited leastSquares 1;
}
divSchemes
{
default bounded Gauss linearUpwind limited;
div(tauMC) Gauss linear;
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(rhoPhi,Yi_h) Gauss MUSCL;
div(rhoPhid,p) Gauss linearUpwind grad(U);
div(rhoPhi,K) Gauss limitedLinear 0.5;//
div(rhoPhi,omega) Gauss linearUpwind default;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear orthogonal;
}
interpolationSchemes
{
default linear;
reconstruct(rho) vanLeer;
reconstruct(U) vanLeerV;
reconstruct(h) vanLeer;
reconstruct(e) vanLeer;
reconstruct(p) vanLeer;
reconstruct(speedOfSound) vanLeer;
}
snGradSchemes
{
default orthogonal;
}
wallDist
{
method meshWave;
}
// ************************************************************************* //
