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

[strive_DUT]

模型方程如下：

我写的代码如下：

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]

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