CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    2
    3
    712
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • strive_DUT
      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方程或者别的什么地方出现了问题

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        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
        

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

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        strive_DUT 1 条回复 最后回复 回复 引用
        • strive_DUT
          strive_DUT @李东岳 最后由 编辑

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

          1 条回复 最后回复 回复 引用
          • First post
            Last post