Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. OpenFOAM11中使用surfaceFilme 模型进行燃烧仿真,温度失控导致浮点溢出

OpenFOAM11中使用surfaceFilme 模型进行燃烧仿真,温度失控导致浮点溢出

已定时 已固定 已锁定 已移动 OpenFOAM
3 帖子 2 发布者 1.5k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • BlancheDHB 离线
    BlancheDHB 离线
    BlancheDH
    写于 最后由 编辑
    #1

    各位老师,师兄,师姐,大佬们好,
    在OpenFOAM11中,我使用surfaceFilm模型进行一个3D仿真。仿真流程如下:

    1. 先运行一个关闭surfaceFilm和spray injection 模型,让T_inlet=900k 运行0.05s 初始化流场。
    2. 激活燃烧,喷雾和薄膜模型,向模型中注入正庚烷C7H16(LPT),然后进行燃烧仿真。
      在只使用燃烧和喷雾模型的时候,仿真可以正确的进行并获得有物理意义的结果。
      在激活surfaceFilm模型后,模型设置参照openFOAM自带算例的“heatBox”,仿真疑似因为在薄膜位置计算温度时浮点溢出而不能运行。
      在"CFDonline"中,我查到在之前的版本中 可以使用“fvOption” 来对温度的上下限进行限制以防温度失控,但是在OpenFOAM11版本中似乎没有该功能。我在中外各个网站查了一圈,有关于该模型的使用案例似乎特别少。
      请问这个问题出现的原因是什么?我应该如何解决呢
      下附报错图片,我认为比较有用的设置文件截图,以及我的仿真文件
      报错如下
      923c4a3e-4fb9-4a2b-b54d-25bc8304cb0a-image.png
      可能有用的文件截图:

    film的温度

    75936490-1654-4132-84fb-120deadfa39e-image.png

    film的fvModel

    fb323613-b62c-4f7e-8895-5958d77ff541-image.png

    关闭薄膜模型时的结果

    77f671c0-e7d5-4036-8341-fd5a4466a7e8-image.png
    求求大佬们救救孩子吧~~~~~~~:qichuang:
    附上我的仿真文件下载链接:
    https://www.jianguoyun.com/p/DTKpvfgQlPqQDBi01tUFIAA

    1 条回复 最后回复
  • BlancheDHB 离线
    BlancheDHB 离线
    BlancheDH
    写于 最后由 编辑
    #2

    自回一下:
    问题已解决,经过多次尝试需要在film的physical properteis 中增加温度限制,粒子数据之类的信息,还需要查询具体的粘度系数和热熔系数公式,我把文件内容放在下面,希望对后来人有用

    thermoType
        {
            type        heRhoThermo;
            mixture     pureMixture;
            transport   polynomial;
            thermo      janaf;
            equationOfState perfectGas;
            specie      specie;
            energy      sensibleInternalEnergy;
        }
    
    mixture
    {
        specie
        {
            molWeight       100.20557;
        }
        thermodynamics
        {
            Tlow        200;
            Thigh       3500;
            Tcommon     1000;
            highCpCoeffs    ( 20.4565203 0.0348575357 -1.09226846e-05 1.67201776e-09 -9.8102485e-14 -32555.6365 -80.4405017 );
            lowCpCoeffs     ( 11.1532994 -0.00949419773 0.000195572075 -2.49753662e-07 9.84877715e-11 -26768.8904 -15.9096837 );
        }
        transport
        {
            muCoeffs<8>     ( -10.378 0.244 -0.000054  0 0 0 0 0 ); //( 1000 -0.05 0.003 0 0 0 0 0 );
            kappaCoeffs<8>    ( 80.10 3.45542 0.00128822 -0.00000166835 0 0 0 0 );
        }
        equationOfState
        {
            rho         684.0;  // 庚烷的密度,单位 kg/m³
        }
        specie
        {
            Hf          0;
        }
    }
    sigma
    {
        type    constant;
        sigma   0.07;
    }
    
    deltaWet    1e-4;
    
    
    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #3

    :146: :146: :146:

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]