新编译的湍流模型求解算例的时候epsilon方程发散

strive_DUT

模型方程如下：
QQ截图20200920203405.png
我写的代码如下：

volScalarField G("RASModel::G", mut_*(gradU && dev(twoSymm(gradU))));

    volScalarField E = Ceps3_*(scalar(1)-fw_)*mu()*(mut_/rho_)*fvc::magSqrGradGrad(U_);//新增
    const volScalarField T(Ts());//新增

    // Update epsilon and G at the wall
    epsilon_.boundaryField().updateCoeffs();

    // Dissipation equation
    tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(rho_, epsilon_)
      + fvm::div(phi_, epsilon_)
      - fvm::laplacian(DepsilonEff(), epsilon_)
     ==
        Ceps1_*G/T//C1_*G*epsilon_/k_
      - fvm::SuSp((2.0/3.0)*Ceps1_*rho_*divU/T, k_)//- fvm::SuSp(((2.0/3.0)*C1_ + C3_)*rho_*divU, epsilon_)
      - Ceps2_*rho_*f2()*epsilon_/T//fvm::Sp(Ceps2_*rho_*f2()/T, epsilon_)
      - E
    );

    epsEqn().relax();

    // No longer needed: matrix completes at the point of solution
    // HJ, 17/Apr/2012
//     epsEqn().completeAssembly();

    solve(epsEqn);
    bound(epsilon_, epsilon0_);

报错如下：

pEqn.H is ok
min(p) [1 -1 -2 0 0 0 0] 3.99999e+06	max(p) [1 -1 -2 0 0 0 0] 4.15888e+06
piso is ok
BiCGStab:  Solving for fw, Initial residual = 1, Final residual = 0.000964975, No Iterations 69
new cannot satisfy memory request.
This does not necessarily mean you have run out of virtual memory.
It could be due to a stack violation caused by e.g. bad use of pointers or an out of date shared library
已放弃 (核心已转储)

麻烦大家看一下epsilon方程或者别的什么地方出现了问题

李东岳

new cannot satisfy memory request.
This does not necessarily mean you have run out of virtual memory.
It could be due to a stack violation caused by e.g. bad use of pointers or an out of date shared library

重新编译一下求解器试试？

strive_DUT

@东岳我重新编译试过了，还是这个问题呢