CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    komega模型求解传热算例的报错

    OpenFOAM
    3
    5
    949
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 一
      一颗鸭蛋 最后由 编辑

      各位老师好,我是在DNS求解器的基础上引入湍流kOMega模型,并加入温度方程引入湍流热扩散系数(参照OF2.3中的传热求解器buoyanBoussinesqPimpleFoam 中的热扩散系数代码),求解器编译通过。

      参考代码:

      alphat = turbulence->nut()/Prt;
      alphat.correctBoundaryConditions();
      
      volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);
      

      报错如下:

      Reading field T
      
      Reading/calculating face flux field phi
      
      Reading field phi_0
      
      Creating turbulence model
      
      Selecting incompressible transport model Newtonian
      Selecting turbulence model type RASModel
      Selecting RAS turbulence model kOmegaMHD//这里以上都能正确读取各个声明的场和模型,接下去本应该读取热扩散系数的alpha场的,但报以下一些错误。
      bounding omega, min: 0 max: 0 average: 0
      [39] [1] #0  Foam::error::printStack(Foam::Ostream&)[2] [11] [44] #[25] #0      
      Foam::error::printStack(Foam::Ostream&)#0  
      Foam::error::printStack(Foam::Ostream&)[26] #0      
      Foam::error::printStack(Foam::Ostream&)[40] #0  
      
      [16] #1  Foam::sigFpe::sigHandler(int) at ??:?//这个错误是除零错误
      [37] #1  Foam::sigFpe::sigHandler(int)[20] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
      [2] #1  Foam::sigFpe::sigHandler(int) at ??:?
      [47] #4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
      [23] #3  Foam::divide(Foam::Field<double>&, 
      Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
      
      [kdd:28991] *** Process received signal ***
      [kdd:28991] Signal: Floating point exception (8)//浮点数例外错误
      [kdd:28991] Signal code:  (-6)
      [kdd:28991] Failing at address: 0x3e80000713f
      
      --------------------------------------------------------------------------
      mpirun noticed that process rank 30 with PID 31201 on node kdd exited on signal 8 (Floating point exception).
      --------------------------------------------------------------------------
      2 total processes killed (some possibly by mpirun during cleanup)
      

      自带的模型中的Omega方程中的热扩散系数alpha如何理解?我本来将代码中的综合热扩散系数alphaEff去掉,并把alpht改成alpha来表示湍流热扩散系数来与模型保持一致,但依然报如上错误。

      OF自带的模型中的Omega方程

      tmp<fvScalarMatrix> omegaEqn
      (
          fvm::ddt(omega_)
        + fvm::div(phi_, omega_)
        - fvm::laplacian(DomegaEff(), omega_)
       ==
          alpha_*G*omega_/k_
        - fvm::Sp(beta_*omega_, omega_)
      );
      
      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        Selecting incompressible transport model Newtonian
        Selecting turbulence model type RASModel
        Selecting RAS turbulence model kOmegaMHD//这里以上都能正确读取各个声明的场和模型,接下去本应该读取热扩散系数的alpha场的,但报以下一些错误。
        bounding omega, min: 0 max: 0 average: 0
        [39] [1] #0  Foam::error::printStack(Foam::Ostream&)[2] [11] [44] #[25] #0      
        Foam::error::printStack(Foam::Ostream&)#0  
        Foam::error::printStack(Foam::Ostream&)[26] #0      
        Foam::error::printStack(Foam::Ostream&)[40] #0  
        

        不太好判断呀,一上来就报错了,不太好判断

        2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • 一
          一颗鸭蛋 最后由 编辑

          是啊,老师。希望有人遇到过,给我解解惑,哈哈~

          1 条回复 最后回复 回复 引用
          • bestucan
            bestucan 版主 副教授 最后由 编辑

            我记得openfoam的DNS不能用并行。可以试试不用并行能不能算。

            滚来滚去……~(~o ̄▽ ̄)~o 滚来滚去都不能让大家看出来我不是老师么 O_o

            异步沟通方式(《posting style》from wiki)(下载后打开):
            https://www.jianguoyun.com/p/Dc52X2sQsLv2BRiqnKYD
            提问的智慧(github在gitee的镜像):
            https://gitee.com/bestucan/How-To-Ask-Questions-The-Smart-Way/blob/master/README-zh_CN.md

            1 条回复 最后回复 回复 引用
            • 一
              一颗鸭蛋 最后由 编辑

              问题解决了,是因为我将Omega设置为零了,所以在计算的时,湍流模型在算湍流黏性系数是nut=k/Omega 出现除零错误。

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