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. Fluent
  3. 我用fluent模拟含不凝气体水蒸气凝结的例子,自编udf,运行时存在问题,希望大神们帮忙看一下

我用fluent模拟含不凝气体水蒸气凝结的例子,自编udf,运行时存在问题,希望大神们帮忙看一下

已定时 已固定 已锁定 已移动 Fluent
19 帖子 6 发布者 13.9k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 中 离线
    中 离线
    中国石油大学研究生
    写于 最后由 编辑
    #1

    案例中有air、水蒸气和液体水,自己编写了定义物性的udf和水蒸气、液体水的源项以及能量源项所有的udf都可正常编译和加载。
    单独加载物性udf时,模拟正常运行。
    当加载水蒸气源项时,一运行就会提示error:received a fatal sognal (segmentation fault)
    当加载液体水源项时,一运行就会提示error:divergence detected in amg solve : x-momentum
    1.png
    (该案例所打开的模型,vof模型和组分输运模型,组分输运模型定义气体mixture包含air和水蒸气,vof模型:第一相是气体mixture,第二相是液体水)
    5.png
    6.png
    (所加载的液体源项)
    7.png
    (加载气体源项时,运行出现的错误)
    8.png
    (加载液体源项时,运行出现的错误)
    下面是我所修改的蒸汽源项的udf
    11.png
    12.png
    液体的源项udf跟气相udf大致相同
    QQ图片20210317193807.png
    (这个是源项的公式)
    求大神们帮忙看一下问题出在大概那个位置
    拜谢!!!

    姜 1 条回复 最后回复
  • zousiyuZ 离线
    zousiyuZ 离线
    zousiyu
    写于 最后由 编辑
    #2

    按我的经验来看,自定义源项后发散,一般都是源项太大了。你可以这么测试,

    1. 设定源项大小为较小的常数,比如 1e-6,如果不发散,基本就确定是源项太大了;
    2. 用 UDM 把源项储存下来,看看大小。
    中 1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 zousiyu 最后由 编辑
    #3

    @zousiyu 感谢大神的建议,上午的时候我将源项直接设定为一个很小的常数,fluent还是运行不了,还是udf本身存在问题我准备一块一块的试错,看看到底是那个地方出现了问题

    1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    写于 最后由 编辑
    #4

    问题已经解决了,是因为上面的程序没有对网格进行循环,只是对一个网格的所有面进行了循环,所以才存在问题

    D B 2 条回复 最后回复
  • D 离线
    D 离线
    Ding
    写于 最后由 编辑
    #5

    你好,我是大四学生,正在做毕设,换热器湿工况的模拟,也在编写冷凝的udf,想向学长咨询一下,请问你做的是二维工况么?为什么在气液共存区也会发生冷凝啊,不应该只发生在接触表面么? 能加个联系方式交流一下么,万分感谢。qq1554024438

    中 2 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 Ding 最后由 编辑
    #6

    @Ding 我做的是三维工况,二维工况也是同理。你这样想一下,冷凝液是会变厚的,原理上来说应该涉及到相平衡,气液界面就是一个气液分子热运动平衡界面,当符合相应条件,是会发生相应的冷凝或是蒸发工况的

    1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 Ding 最后由 编辑
    #7

    @Ding 我现在也在忙着自己毕业的事情,有点自顾不暇,没有太多空闲的时间对别人进行帮助。
    如果你有问题直接在这里问就可以,有大神会帮你解答,我空闲时间也会进来看一眼,如果我知道我会直接回答的。
    见谅

    D 1 条回复 最后回复
  • D 离线
    D 离线
    Ding
    在 中回复了 中国石油大学研究生 最后由 编辑
    #8

    @中国石油大学研究生 那能分享一下udf代码么?邮箱1554024438@qq.com

    中 1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 Ding 最后由 编辑
    #9

    @Ding 这个还是得你自己编写,我可以帮你看看问题,我的成品不能直接给你

    1 条回复 最后回复
  • D 离线
    D 离线
    Ding
    在 中回复了 中国石油大学研究生 最后由 Ding 编辑
    #10

    @中国石油大学研究生 对所有网格循环是怎么编写呢?

    中 1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 Ding 最后由 编辑
    #11

    @Ding 这个需要用begin_c_loop宏进行实现,即在外层套上这个宏即可

    1 条回复 最后回复
  • D 离线
    D 离线
    Ding
    写于 最后由 编辑
    #12

    QQ截图20210427163433.png

    1 条回复 最后回复
  • D 离线
    D 离线
    Ding
    写于 最后由 编辑
    #13

    QQ截图20210427163502.png

    中 1 条回复 最后回复
  • D 离线
    D 离线
    Ding
    写于 最后由 编辑
    #14

    能帮我看下问题么?初始化就卡住了,怀疑程序是不是程序是死循环QQ截图20210427165654.png 我参考的这篇文献

    1 条回复 最后回复
  • 中 离线
    中 离线
    中国石油大学研究生
    在 中回复了 Ding 最后由 编辑
    #15

    @Ding 你需要一句一句的看。
    具体方法是:从上至下把无关的先注释掉,一句一句往里加,每一次都用fluent进行编译、初始化、模拟,看是加入那一句之后才出现的问题,这样才能看出是哪里出现了问题。
    你这样单独看是很难发现问题的

    1 条回复 最后回复
  • 彼 离线
    彼 离线
    彼岸
    写于 最后由 编辑
    #16

    楼主,最后问题解决了吗

    1 条回复 最后回复
  • 姜 离线
    姜 离线
    姜磊
    在 中回复了 中国石油大学研究生 最后由 编辑
    #17

    @中国石油大学研究生
    您好,
    1.麻烦问一下您udf中的mass_transfer_coeff=vapor_densitydiffusivity/distance(C_YI(cell,first,0)- water_mass_fraction)/(1-water_mass_fraction);
    替代文字
    为什么/distance*(C_YI(cell,first,0)- water_mass_fraction)这一段能表示公式最后的偏导数?

    我的理解是您想把偏ω改为Δω,distance是否表示偏n?

    您udf中说water_mass_fraction应该和C_YI(cell,first,0)所获取的值相同,那么C_YI(cell,first,0)- water_mass_fraction是否结果为0,该如何理解?

    偏ω/偏n是否可以想办法写为相邻两网格的ω差/两网格的体心距离?

    2.我看一些文章里说上式的ρ应该是湿空气的密度,您udf里用的是水蒸气的密度,请问用哪一个密度更合理

    3.air_molecular应该是干空气气体摩尔质量,取28.9合适吗?

    4.最后乘10倍是系数吗?我看一些文章中是没有乘10的。

    1 条回复 最后回复
  • B 离线
    B 离线
    being curious
    在 中回复了 中国石油大学研究生 最后由 编辑
    #18

    @中国石油大学研究生 不是把,define source 宏会自动循环所有网格的啊?

    1 条回复 最后回复
  • B 离线
    B 离线
    being curious
    写于 最后由 编辑
    #19

    想问下楼主最后成功实现仿真了,我目前已经完成了UDF的编写,并且在fluent中也能正常计算,但是液膜厚度一直保持在第一层网格内,不知道是什么原因。希望可以交流一下。

    1 条回复 最后回复

  • 登录

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