爆炸计算压力震荡

[xuqiming]

求助，最近基于blastReactingFoam进行气体爆炸计算，但是得到的爆炸超压震荡非常厉害，与实际显然有出入，感觉可能是离散格式设置不合理，这一块找不到什么资料进行参考修改。各位大佬可否给予指点。
基本设置如下
求解器：blastReactingFoam，湍流模型：k-omega SST
网格：二维，全正交网格，网格尺寸0.5 mm

/*--------------------------------*- 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;
}

// ************************************************************************* //