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



  • 模型方程如下:
    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
    

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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2