Skip to content
  • 有关五面体网格

    Algorithm
    2
    0 赞同
    2 帖子
    5k 浏览
    李东岳

    Hi,

    我遇到的使用三菱柱的情况就是四面体的边界层网格。当然需要使用三棱柱来做边界层,有利于使用三棱柱的上下面把误差抵消掉。

    计算效率是指?
    目前出了纯结构网格,其他的网格在做索引的时候编程都要比结构网格复杂。
    如果值得是精度,三棱柱在边界层肯定要比4面体好,这个在教材中都有提及。

  • 0 赞同
    2 帖子
    6k 浏览
    李东岳

    没想到图越来越多。手机看的注意了。我这个图多。

    结构网格的拥护者简单说一下为什么我喜欢结构网格。以下只针对简单几何以及不是太复杂的几何。复杂到不行的网格我选择混合网格。我自己做的算例不用纯四面体网格。

    六面体网格具有以下优点:

    非矩形网格的歪斜(skewness)严重
    这需要做面法相梯度计算修正。稍微增加计算步骤。著名的CFD非矩形修正因此而来

    截面重构不光顺
    举例,一个正方盒子,一半是水,我们用矩形网格重构的气液截面是这样的:
    alt text
    黄色为一相,白色为一相。即使网格很糙,我们的截面也是平的。下图为三角形非结构网格:
    alt text
    可见糙网格下界面不平。这导致后处理需要耗费精力。

    高度非线性物理问题影响矩阵对角占优
    这也是非线性动力学,例如金属撞击、生物力学分析中结构网格肆虐的原因。(这个我没有测试,之前在网上看到的结论,来自链接文本

    非结构网格数量更多增加计算资源
    网格单元的面越多网格越少。六面体网格并不是数量最少的。数量最少的是多面体网格(star推荐的?)如下图
    alt text
    再看这个情况对比:只有第一个是结构网格,第二个网格多。第三个第四个和均为非正交且和流线不贴合
    alt text
    alt text
    因此目前自动生成的非结构网格一般数量是即为巨大的。计算的时间会超过你花一星期做结构网格的时间。另外,万一中途发散了怎么办。

    费结构网格高度不可控
    几何稍微一点点的变化会导致全场四面体网格做位移。结构网格(ICEMblock切出来的)具有即为即为高度可控性。可以对基本所有的网格单元进行控制。俩者区别天差地别。

    我了解的所有的离散格式对结构网格适应性更强
    也就是说结构网格对流场具有强的多的顺应性。例如旋转机械,结构网格和速度流线高度贴合。

    矩形网格对边界更好的捕获
    为何非结构网格要添加边界层网格?矩形网格对壁面贴合更好。

    复杂CFD模型在结构网格上收敛更强
    个人经验。比如双流体模型,LTS局部时间步框架下的模型,激波捕获等。当然那种不可压缩单相流不算

    太多就不写了...

    引用从业CFD27年OpenFOAM创始人Henry的一段话:

    hex and polyhedral cells are generally superior to tet cells in finite-volume and the difference in both accuracy and convergence behaviour can be large. However, it also depends on the quality of the hex and tet cells and many of the problems with tet meshes stem from the poor quality of some tets produced by many automatic tet mesh generators. In order to really appreciate the difference it might be a good idea for you to try both for your problems. For VOF calculations I would expect the solution on hex/poly cell meshes to be dramatically better.

    大意为:

    在有限体积中,六面体和多面体网格相对于四面体网格优势明显。主要在精度和收敛上。对于VOF(我上文提及的第2点),结构网格的优势是大大的。

    最后,需要耗费大量时间生成的六面体网格至今没有消失,必然有存在的道理。

  • AIChE Journal非TOP期刊

    CFD彩虹条
    4
    0 赞同
    4 帖子
    7k 浏览
    幻想de罗卜

    @东岳 惭愧啊。整个实验室就自己在零零碎碎做一些流体力学模拟。大多数人都偏催化、化学实验类的。一般上论坛来也就看看帖子。都忘记cookie过期没登陆了。。
    我们这边还是很认可AIChE J的。学校认定上还是top级别的。

  • 单元区域设置

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    W

    怎么能够对计算完的case设定单元区域,取得这个区域内自己想要的值

  • OpenFOAM冷凝模型在哪

    OpenFOAM
    2
    0 赞同
    2 帖子
    3k 浏览
    bestucanB

    https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1compressible_1_1humidityTemperatureCoupledMixedFvPatchScalarField.html

    https://www.cfd-online.com/Forums/openfoam-solving/147129-condensation-using-openfoam.html

    https://www.researchgate.net/publication/279643607_Modeling_of_Direct_Contact_Condensation_With_OpenFOAM

    https://jhmtr.semnan.ac.ir/article_467.html

    https://github.com/NimaSam/phaseChangeHeatFoam

  • 0 赞同
    1 帖子
    1k 浏览
    W

    各位前辈,我在进行一个压缩理想气体的仿真,使用了fluent自带的内燃机动网格来模拟活塞的上下运动,我这个三维压缩腔模型没有进出口,边界都是壁面,这个压缩腔内含有一部分的水,并且空气并不会溶解在水中。在时间步长为10e-05时,压缩腔内的体积平均空气质量会接近不变,但当时间步长为10e-04时,压缩腔内的体积平均空气质量会大幅度降低,导致计算结果出现很大的误差。谁知道这是怎么回事吗?270c7d84-3e4a-4a3b-a0eb-f03246e0389d-8022a33ca0ae4229656cf7dfbaaf2691.png
    32391c80-98f4-4a1f-ad2b-be648153e223-b2e6b74bb524d804dfbf9eb5304f1082.png
    bb35fd93-e478-4c6f-9769-18e0995e0a10-9d0248f522d25c2be2fd28eb600910ec.png
    596b964a-bba4-41d7-8dce-e471b145cc16-cc8a6de7ae51246153c122cd626a3b9.png

  • 博士毕业那年我才发了第一篇SCI

    CFD彩虹条
    4
    1 赞同
    4 帖子
    8k 浏览
    zousiyuZ

    心态是真的棒,怪不得科研做得好👍

  • CFD青年成长支持计划(2022)

    公告
    87
    0 赞同
    87 帖子
    216k 浏览

    我是南京航空航天大学航空学院2022级硕士研究生吴祥清,跟随张老师链接文本从事CFD研究,本人的主要研究方向为高速出入水,初步学习商软入门CFD现在准备编程写算法。查找学习资料了解到CFD中文网以及老师,所以想申请2022CFD青年成长支持计划,我目前正在使用LS-DYNA软件进行垂直破冰仿真,开始着手发相关论文。在本科期间接触流体力学便想深入了解这个专业,也申请了这个方向的研究生,如果论文写作顺利,便会在研究生二年级申请硕博连读,为国家军事建设贡献自己的一份力量。目前在上课之余也开始着手C++入门。但是CFD学习之路道阻且长,希望可以得到老师资助来更好的学习CFD。

  • 0 赞同
    4 帖子
    6k 浏览
    L

    @liuxiaodong les计算圆管管道流动的,尝试了很多都最后变成层流了,这是为什么呢

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    647k 浏览
    H

    我写 固我在,嘿嘿

  • Fluent Maxwell速度滑移和温度跳跃

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    E

    请问有大佬知道如何在Fluent中利用UDF添加Maxwell速度滑移和温度跳跃边界条件吗

  • 我,想要,鼓励!

    CFD彩虹条
    5
    0 赞同
    5 帖子
    8k 浏览
    肖艺

    鼓励鼓励!:high:

  • 读博与职业规划

    CFD彩虹条
    4
    0 赞同
    4 帖子
    6k 浏览
    李东岳

    国内形势变化太快了,一年一年的水涨船高,就像我们看现在的教授,当年竟然博士毕业就留校了。以后的年青一代看我们,估计是“这家伙当年竟然博士发了5篇文章就留校了,现在我们博士发10篇都留不下来”..

  • 0 赞同
    4 帖子
    9k 浏览
    李东岳

    @搬运工

    高斯迭代法处理20阶矩阵就已经非常慢了。通常和多重网格+高斯赛德尔求解器一起处理。并且OpenFOAM这种成熟的商软里面矩阵求解器算法非常精炼。通常由主求解器+预处理+光顺+重组构成。多重网格更是复杂。另外,在编程的时候使用动态、寻址的方法寄存系数。

    所以如果你要写一个快速的求解器:

    使用外挂包; 潜心研究稀疏线性系统求解理论;我建议花大量时间研读:YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:
  • 0 赞同
    6 帖子
    9k 浏览
    Y

    @myheart 好的好的,十分感谢,我再去研究一下

  • 0 赞同
    1 帖子
    2k 浏览
    小狗狗

    OpenFOAM中的DPM粒子追踪功能,在运动网格(ALE)上的支持并行计算吗?有没有什么限制?有大佬测试过吗

  • 0 赞同
    2 帖子
    3k 浏览

    个人感觉不一定是网格划分的问题,有可能是显示的问题

  • icem网格划分问题

    Fluent
    2
    0 赞同
    2 帖子
    2k 浏览
    李东岳

    看起来没问题呢

  • icem中模型重合面网格划分问题

    Meshy
    2
    0 赞同
    2 帖子
    5k 浏览
    Y

    solidworks导入icem时候有些地方确实有点问题,我现在先用sw把大框架画出来,自己在icem里面再编辑一下小框架,最后再切分一下,有点麻烦。你现在已经解决把,有啥好办法吗

  • 0 赞同
    5 帖子
    8k 浏览
    肖艺

    @yfclark