Skip to content
  • 0 赞同
    2 帖子
    3k 浏览
    N

    自己捞一下,有没有大佬知道怎么整

  • 0 赞同
    2 帖子
    6k 浏览
    W

    Navier-Stokes 方程中的对流项写作 $\nabla \cdot (\textbf{u}\phi)$,扩散项写作 $\nabla \cdot (\nu \nabla \phi)$。

    对流表示流体微团从空间某点运动到另一点的过程,而扩散则表示流体物理性质浓度变化的过程。

    对动量方程而言 $\phi=\textbf{u}$,因此对流项变为 $\nabla \cdot (\textbf{u} \textbf{u})$,为非线性项。一些算法将其中一个速度作为已知量,另一个作为未知量离散求解。可参考icoFoam解析。

  • 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点),结构网格的优势是大大的。

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

  • 0 赞同
    4 帖子
    4k 浏览
    Z

    @renao7 是电脑的运行内存不够高吗?:134:

  • 0 赞同
    1 帖子
    1k 浏览

    各位老师好,我使用paraview5.9.1进行颗粒的后处理,vtk文件大约900MB,使用Glyph Arrow时可以正常运行,但当换成Sphere后出现无响应,弹窗情况,如下
    08c7546d-d581-4698-ad66-4a35b8c9393f-屏幕截图 2024-10-10 223916.png
    无报错提示,随后闪退,查看系统事件簿出现如下错误提示
    c7b727f2-a987-42b4-8f29-39bea0225940-image.png
    请问各位老师这是什么原因呢

  • 0 赞同
    3 帖子
    4k 浏览
    2

    感谢东岳老师,我尝试着改小了,问题得到了部分解决。

  • DEM和DPM

    已移动 OpenFOAM
    103
    0 赞同
    103 帖子
    471k 浏览
    Y

    @AppleKiller 行,我再看看,研究了好久感觉都没问题,我再重新看看我的函数

  • CFD算法编程课 买家秀

    CFD彩虹条
    23
    0 赞同
    23 帖子
    80k 浏览
    WeakFormW

    好的感谢李老师

  • fluent udf中体积分数宏 C_VOF 的取值问题

    Fluent
    5
    0 赞同
    5 帖子
    10k 浏览
    W

    thread_loop_c 获取subthread层的值似乎有点问题,见https://muchong.com/html/201211/5150612_2.html,最好是用mp_thread_loop_c,配合**pt和C_VOF(c, pt[0])来进行取值。

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

    bdbd6309-2879-476a-a5fd-79f685f05403-1731635955698.png
    udf也没用编译错误,网格也是正常的

  • 0 赞同
    5 帖子
    7k 浏览
    zzkluckZ

    说一点我的看法,不一定对,仅供参考:

    non-orthogonality

    dca809a8-9cd5-4f74-9329-6d8459966264-image.png

    non-orthogonality由snappyHexMeshDict.meshQualityControls.maxNonOrtho直接控制,你的snappyHexMeshDict里这个值是45,这个值一般不需要控制的这么小,65是更常见的配置值。

    另外,你的snappyHexMeshDict里,relaxed.maxNonOrtho给的是75,这个值通常是和maxNonOrtho = 65相配合的。如果你真的要把maxNonOrtho设为45,按逻辑来说relaxed.maxNonOrtho也应该相应的减少一点。

    网格纵横比

    这里截取你的背景网格配置:

    vertices ( (-340 -302 -3) // Slightly smaller than the STL model bounds (332 -302 -3) (332 341 -3) (-340 341 -3) (-340 -302 95) // Slightly larger than the STL model bounds (332 -302 95) (332 341 95) (-340 341 95) ); blocks ( hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1) // Adjust mesh density as needed );

    如上面B老师所说,你网格的大纵横比就是这里导致的。如果不是有意要在Z方向上做加密的话没必要这么分块,我口算一下,大概(70 65 10)就可以。

    castellate

    从上面可以看到,你背景网格的量级是100 * 100 * 100 = 1,000,000;而你SHM的maxGlobalCells给了2,000,000。这看起来不太对,这几乎没给SHM进一步细分的空间,从结果来看也是这样,网格划分的不够细。

    另外,截取一段level的配置:

    features ( { file "building.eMesh"; level 3; } ...... ); refinementSurfaces { building { level (3 4); } ...... }

    我理解features的level应该不低于Surface的。

    胡乱写了一些个人看法,希望能有所帮助。

  • 0 赞同
    3 帖子
    7k 浏览
    A

    您好,请问这里的Re_thetac如果不修改,即为none时,Re_thetac的默认值是多少?

  • 李博最近不水群了,大家甚是想念

    CFD彩虹条
    2
    0 赞同
    2 帖子
    5k 浏览
    李东岳

    @J 9年过去了...

  • fluent的overset出问题

    Fluent
    2
    0 赞同
    2 帖子
    3k 浏览
    C

    楼主,请问您解决这个问题了么,我最近做重叠网格入水也遇到这个问题,求助

  • 边界条件梯度为0 具体代表什么?

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

    @ywj123450
    《数值传热学(第2版)》是我唯一力荐的一本中文CFD理论书。

  • ICEM绘制的网格导入openfoam之后出现异常

    Meshy
    4
    0 赞同
    4 帖子
    4k 浏览

    @李东岳 好的,谢谢李老师

  • 什么是averaging turbulence?

    Algorithm
    6
    0 赞同
    6 帖子
    11k 浏览
    .J..

    @showhand 多谢大神回复!

  • 0 赞同
    4 帖子
    10k 浏览
    hangszH

    SRF和MRF只是坐标系变换,将旋转问题变成stationary问题,所以湍流模型还是适用的。

    但是对于旋转问题或者流线弯曲问题,湍流模型是要进行修正的。

    参看: 《Some improvements in Menter's k-omega SST turbulence model》

  • 每日一问:有关icem网格检查失败问题

    Meshy
    1
    0 赞同
    1 帖子
    2k 浏览
    2

    icem网格划分导入fluent时检查失败,显示Error: interface zone 3 has two adjacent cell zones.
    4f68eb65-132e-4314-b0f6-6964dc466673-image.png

  • 0 赞同
    16 帖子
    11k 浏览
    Y

    @coolhhh 好的好的 真的万分感激