CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    Fluent模拟水下燃气射流(类似于水下爆炸),水中出现前导激波,如何设置出口边界条件?

    Fluent
    边界条件 fluent 前导激波
    4
    9
    2270
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 甘九_Akrasia
      甘九_Akrasia 最后由 编辑

      各位CFDer们,起个小早提个问题。:duang:
      一句话描述:利用fluent进行水下燃气射流模拟时,气、水均采用可压模型的话,在有前导激波存在的情况下,如何设定出口边界条件,使得前导激波的压强和出口边界相匹配?
      -------------------------------
      详细描述:计算模型如图所示,全域分块划分结构化网格。采用压力基求解,因为密度基不支持多相流模型,所以密度基下的自由边界条件和非反射边界条件(咱也不懂,咱也不敢问)就不能用了。
      d3d0d878-c532-417c-b4da-98ef4b258833-image.png
      EF入口处利用profile文件给定高温高压高速燃气喷射(pressure-inlet),选取vof模型,爆轰燃气采用ideal-gas,液体水采用compressible-liquid,不考虑表面张力和重力,湍流模型选取标准k-e。入口参数如下:红色是管口马赫数曲线,下x坐标轴从左到右;黑色是管口静压曲线,上x坐标轴从右到左。
      5120c85e-9f53-47e3-b3ec-a93959cb53bb-image.png
      AB、BC、CD处为压力出口边界条件(pressure-outlet),指定静压为101325Pa,其余参数默认。
      205b3056-11e2-4c34-b560-9de46f94e1c4-image.png
      当考虑了水的可压缩性之后,出口的高压燃气撞击气水接触面,会在水中形成一道前导激波,如下图所示。以下均为压力云图
      875e2e90-df3a-40e9-98e8-af7993289bcd-image.png
      当这道前导激波到达边界处时,由于边界上指定静压为101325Pa,前导激波的压力和边界处指定的压力不匹配,压力骤降会形成膨胀波回传,然后整个计算域的压力场就乱!掉!了!
      c7c1befd-e239-41f8-ab1f-93e85a5823d5-image.png ad8350f0-46ce-468c-9c85-b0ec2186c8cc-image.png
      求教如何设置这种情况下的出口边界条件。:zoule:
      我18年9月份入学(C)研(F)究(D)生(keng),这个是我毕业设计的题目。但是毕业设计的时候设置的水是不可压缩模型,这样子会导致水中的声速无穷大,球形压力波瞬间传遍全场,不会出现前导激波,边界处就不会出问题。:zoule:
      大家看这个图上的时间是19年1月份,我当时想着把水设成可压模型,算出来前导激波了很开心,但是下一秒到了边界就崩了,很难过。:zoule:
      后来开学了就跟着课题组山南海北地跑项目试验,最近想着把这个问题捡起来处理处理。谢谢大家了:ok3:

      我是山间惟一的行李,除了走动、铭记,没有一样幸福可以带走:shangxue:

      1 条回复 最后回复 回复 引用
      • 甘九_Akrasia
        甘九_Akrasia 最后由 编辑

        一楼第一行给李东岳老师:xinxin:
        由于李东岳老师的论坛插件和论坛币功能还在做,我就开个感谢帖。
        回复本帖问题者,表示感谢!(鞠躬90°):xiexie:
        对本帖问题提出有用建议的,表示非常非常感谢!!(鞠躬180°):xiexie: :xiexie: 来南京玩请吃学校食堂,坐标南京孝陵卫200号,南京食堂吃起来还行大学:baobao:
        对本问题给出突破性解决方案的老哥er,表示非常非常非常感谢!!!(鞠躬360°)的同时,可以赠送一本CFD相关书籍(几千一本的国外CFD书除外…),或者赠送一本南京诗人的签名诗集,或者你有空来南京我请你吃饭(对这种大神级人物要套关系抱大腿!):chigua2:
        总之,谢谢李东岳老师!谢谢大家!:chouchou:

        1 条回复 最后回复 回复 引用
        • 甘九_Akrasia
          甘九_Akrasia 最后由 李东岳 编辑

          本问题相关帖子
          http://www.cfd-china.com/topic/2424/关于空气自由射流模拟边界条件设置的问题
          对于帖子中 @小旗子annie 同志给出选取non-reflecting的建议,我试着算了。但是算一会就出错,压力还没传到边界上就提示边界处的turbulence参数发散了。:xiezuoye:

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

            我觉得应该用non-reflecting的边界条件,具体Fluent那面怎么做我也不太清楚 :zoule:

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

            1 条回复 最后回复 回复 引用
            • K
              Kingkong 最后由 编辑

              活捉小师弟一枚,师弟,你解决了也来解答下,封个楼。

              1 条回复 最后回复 回复 引用
              • 小
                小旗子annie @甘九_Akrasia 最后由 编辑

                @甘九_Akrasia 小校友你的问题解决了吗?个人感觉计算域不够大,另外我没看到你的湍流模型,一般带喷管的射流我用瞬态计算

                甘九_Akrasia 2 条回复 最后回复 回复 引用
                • 甘九_Akrasia
                  甘九_Akrasia @小旗子annie 最后由 编辑

                  @小旗子annie 解决了解决了,下载了个2019R1,然后在压力出口选项面板上勾选Non-Reflecting就好了。

                  1 条回复 最后回复 回复 引用
                  • 甘九_Akrasia
                    甘九_Akrasia @小旗子annie 最后由 编辑

                    @小旗子annie 计算域确实不够大,我自己的电脑配置不太强,为了计算速度计算域取得比较小。湍流模型是RNG k-e两方程模型,也是瞬态计算。一开始调不好是因为vof主次相设置没处理好,后来画时间啃了啃fluent的help文档,终于弄懂了...(然而现在好像又都忘了...:zoule:
                    现在数值计算是可以完整计算下去了,就在忙实验。等实验结果出来,如果两个能对的上,发现什么有意思的事情,希望可以整篇文章:jingya:

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

                      现在数值计算是可以完整计算下去了,

                      分享一下结果啊,看看你酷炫的结果,这个是2D还是3D算例?

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

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