Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

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

    OpenFOAM
    2
    3
    592
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • strive_DUT
      strive_DUT last edited by

      模型方程如下:
      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 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

        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://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        strive_DUT 1 Reply Last reply Reply Quote
        • strive_DUT
          strive_DUT @李东岳 last edited by

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

          1 Reply Last reply Reply Quote
          • First post
            Last post

          CFD中文网 | 东岳流体 | 京ICP备15017992号-2