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. rhoCentralFoam使用总温总压入口时遇到的问题

rhoCentralFoam使用总温总压入口时遇到的问题

已定时 已固定 已锁定 已移动 OpenFOAM
7 帖子 3 发布者 6.3k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 宝 离线
    宝 离线
    宝丁
    写于2019年4月20日 07:24 最后由 编辑
    #1

    上一个帖子是关于sonicFoam使用总温总压入口时遇到了一些问题,问题似乎已经解决,再次感谢李东岳老师的帮助。

    现在换一个求解器,用rhoCentralFoam又有些别的问题。

    还是给出我要用到的边界条件
    b2400160-500c-4198-b4ac-105009ee515d-1b8c06ed0b5baf6a2a6590f4d76c9af.png

    计算的问题在于进口的第二层网格会出现数值的异常,如下
    9e1302cb-ab57-4be2-b4ac-226901db22b3-图片.png

    我将文件放到这里,大家有兴趣的可以看一下:420_rhoCentralFoam.zip

    这样的问题在CFD_online上面也是有人问到过,但是并没有解决,这是链接:
    https://www.cfd-online.com/Forums/openfoam-solving/81423-rhocentralfoam-totalpressure-totaltemperature-inlet-subsonic-channel.html

    还是说rhoCentralFoam不支持总温总压的入口边界条件呢?

    还有,我将总温总压的入口边界条件换成了速度入口,整个计算收敛的非常好。

    有没有人同样遇到这样的问题呢?虽然以后用不到这个求解器,但是还是想把这个问题解决。

    1 条回复 最后回复
  • 李 在线
    李 在线
    李东岳 管理员
    写于2019年4月20日 08:30 最后由 编辑
    #2

    rhoCentralFoam使用的格式跟音速有关,你得看看你是否满足高超音速的假定,同时库郎数不只满足小于1,而是以波的传播速度来计算,小于波速,基本上要远远小于1,在满足高超音速的时候,你试试降低时间布长看看

    理论:http://dyfluid.com/rhocentralfoam.html

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

    宝 3 条回复 最后回复 2019年4月20日 08:37
  • 宝 离线
    宝 离线
    宝丁
    在 2019年4月20日 08:37 中回复了 李东岳 最后由 编辑
    #3

    @东岳 我的边界条件远没有超音,进口马赫数约0.55, 但是rhoCentralFoam只能用于高超声速吗

    1 条回复 最后回复
  • 宝 离线
    宝 离线
    宝丁
    在 2019年4月20日 08:41 中回复了 李东岳 最后由 编辑
    #4

    @东岳 好像理解错了,降低时间步长试一下

    1 条回复 最后回复
  • 宝 离线
    宝 离线
    宝丁
    在 2019年4月20日 09:21 中回复了 李东岳 最后由 编辑
    #5

    @东岳 我把库朗数降到了0.006,计算发现边界还是发散。rhoCentralFoam本身是一个跨声速求解器,为什么要满足高超音速的假定?

    1 条回复 最后回复
  • S 离线
    S 离线
    starDust
    写于2019年5月1日 06:14 最后由 编辑
    #6

    没有超音不能用rhoCentralFoam吗?
    我在算一个临界状态的喷管,喷管进口也是给的总温总压,计算结果也不是很好。

    1 条回复 最后回复
  • 李 在线
    李 在线
    李东岳 管理员
    写于2019年5月1日 08:11 最后由 编辑
    #7

    并不是这么绝对。就像是密度基求解器和压力基求解器一样。

    理论上,rhoCentralFoam是一个密度基求解器,原本用来求解欧拉方程,方程具有双曲型,不同的特征值。问题在于求解亚音速的时候,方程特征会变化,波动速度传输趋于无穷大。rhoCentralFoam调用的格式还取决于波动速读会导致时间布长太小。可以看看中心迎风格式的原始资料,讨论都是基于双曲型欧拉方程的(具备不同的特征值,非不可压缩) https://www.jianguoyun.com/p/De-T2z4Q9s3ZBhik7bkB

    不过,我并没有测试过使用密度基中心迎风格式求解亚音速问题,从理论上来看,这种收敛性很差。不过,我觉得你们用一个一维案例测试下算法,比如试试用rhoCentralFoam求解1D喷管,逐渐降低速度看看收敛性,再应用你们的实际问题更好,欢迎反馈测试结果!This is scientific research :ok2:

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

    1 条回复 最后回复
2019年4月20日 07:24

2/7

2019年4月20日 08:30

未读 5
2019年5月1日 08:11
  • 登录

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