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. Algorithm
  3. 2D欧拉拉格朗日计算

2D欧拉拉格朗日计算

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

    个人经验来看,我对2D欧拉拉格朗日模拟理解存在一些疑惑,考虑相间动量耦合,在欧拉框架下,可以写为:
    \begin{equation}
    A=\frac{\rho}{\tau}(\bfU_\rd-\bfU_\rc)
    \end{equation}
    方程中不存在网格体积,因此个人认为2D欧拉欧拉方法是适用的。但是对于2D欧拉拉格朗日,曳力导致的动量交换计算大体如:
    \begin{equation}
    A=\frac{1}{V_\mathrm{cell}}\sum F_i
    \end{equation}
    其中$F_i$为每个粒子的曳力,$V_\mathrm{cell}$为网格体积。对于2D网格,如果还是按照通常的计算标准,网格厚度应该不对结果有影响,但是很明显,网格越厚,力越小,导致2D欧拉拉格朗日的模拟结果依赖于网格厚度,看起来不是很合理。

    在Kadyrov et al. Impacts of solid stress model on MP-PIC simulation of a CFB riser with EMMS drag的这篇文章中,用的就是2D欧拉拉格朗日模拟,同时,文中还专门讨论了网格厚度对结果的影响,结论是没什么影响。

    Summing up we can say that, change of domain thickness value does not significantly impact the accuracy of 2D simulations

    但是我个人很疑惑,不知道哪位老师可以讨论下

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

    1 条回复 最后回复
  • bestucanB 离线
    bestucanB 离线
    bestucan 版主 大神
    写于 最后由 编辑
    #2

    我猜一下。

    如果我写伪2D的求解器,我就把密度设成体积相关。这样就保证质量守恒了。这样体积怎么变,单位网格的质量不变。质量不变,力不变,单位时间内冲量、动量效果就一样。

    比如我认为读入的密度是当厚度为1的密度。具体计算中设定厚度是多少,我再换算成特定厚度时的密度。

    而且这种变化只影响与坐标相关的,比如扩散,但是这种变化只在z轴上,z轴上的变化又是不必关心的。

    兴许扣扣那个求解器的密度当2D的时候怎么定义的能找出答案 :shangxue:

    滚来滚去……~(~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

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 李东岳 编辑
    #3

    如果我写伪2D的求解器,我就把密度设成体积相关。这样就保证质量守恒了。这样体积怎么变,单位网格的质量不变。质量不变,力不变,单位时间内冲量、动量效果就一样。

    就像你说的。文章做了厚度方向的验证。在增加厚度的时候,同时增加了相分数(类似你说的密度)。但我觉得,应该仅调整一个参数做分析。

    CFD界这几天要推送一个2D/3D模拟的文章,可以关注下

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

    aiweimoA 1 条回复 最后回复
  • bestucanB 离线
    bestucanB 离线
    bestucan 版主 大神
    写于 最后由 编辑
    #4

    本来第一个想猜的是程序里的“体积”“面积”都是广义的,就是“网格尺寸”和“网格边界”,

    当三维时“网格尺寸”和“网格边界”分别是体积和面积,
    当二维时“网格尺寸”和“网格边界”分别是面积和长度。

    但是太没有依据了
    对着 empty 类型边界条件翻了翻找依据,翻到这翻不下去了 153行读取边界条件类型 ,以后 C++ 知识补上来了再翻。

    如果全是薄3D假装2D,那 empty 这个类型有点多余了。估摸抠出来 empty 会影响什么就知道它2D是怎么被假装出来的了。

    算法上不太懂,代码上挖大概可以这样挖。


    公众号里说 axisymmetric 会丢失信息,我记得 fluent 里也有,同时也有 sxisymmetric swirl。后者好像能捕获前者丢失的信息。OpenFOAM 里也有,不知道一样不。

    :xiezuoye:

    滚来滚去……~(~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

    1 条回复 最后回复
  • aiweimoA 离线
    aiweimoA 离线
    aiweimo
    在 中回复了 李东岳 最后由 编辑
    #5

    @李东岳 在 2D欧拉拉格朗日计算 中说:

    几天要推送一个2D/3D模拟的文章,可以关注下

    @李东岳 李老师,公众号找了一圈,好像没找到,这个公众号文章链接可以发一下么?最近同学投了一篇文章,是基于欧拉-拉格朗日方法的,做的沙尘障碍物绕流2D研究。一个审稿人意见是他不接受2D Model,并且明确表明如果下一稿还是2D Model的话不用给他审直接拒绝。我就想到之前公众号好像有一篇推文。但是搜不到了。

    aiweimoA 1 条回复 最后回复
  • aiweimoA 离线
    aiweimoA 离线
    aiweimo
    在 中回复了 aiweimo 最后由 李东岳 编辑
    #6

    打扰了,我自己翻到了。

    疏影横斜水清浅疏 李东岳李 2 条回复 最后回复
  • 疏影横斜水清浅疏 离线
    疏影横斜水清浅疏 离线
    疏影横斜水清浅
    在 中回复了 aiweimo 最后由 编辑
    #7

    @aiweimo 你好,请问你同学的那篇2D的障碍物沙尘绕流的文章题目是什么。

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    在 中回复了 aiweimo 最后由 编辑
    #8

    @aiweimo 不知道后来这个文章怎么样了?

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

    aiweimoA 1 条回复 最后回复
  • aiweimoA 离线
    aiweimoA 离线
    aiweimo
    在 中回复了 李东岳 最后由 编辑
    #9

    @李东岳 他发的开源期刊求毕业了:quwan:

    李东岳李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    在 中回复了 aiweimo 最后由 编辑
    #10

    @aiweimo 也行,有链接不

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

    1 条回复 最后回复

  • 登录

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