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

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

  • 为什么叫这个名字?

    CFD彩虹条
    3
    0 赞同
    3 帖子
    6k 浏览
    S

    fvScheme文件中各种Guass linear,大概是版主遗毒太深。

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

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

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

  • ICEM划分网格

    Meshy
    13
    0 赞同
    13 帖子
    18k 浏览
    H

    @Junren-Hou 在 ICEM划分网格 中说:

    可以试试check block,有些时候其实没有问题,就它软件自己的毛病,然后block里最后有一个check block,可以让他自己检查一下。里面有好多选项,都试试,什么check/fixed block

    回忆了一下,我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域,一方面追求地面边界层,致使网格尺度过小,一方面为了平衡网格数量,过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

  • C++对象数组初始化

    OpenFOAM
    2
    0 赞同
    2 帖子
    2k 浏览
    S

    找到了一个方法
    用动态对象数组

    IB **myIB=new IB *[3]; for(int i=0;i<3;i++) { myIB[i]=new IB(i+1); }

    注意最后delete。

  • 0 赞同
    3 帖子
    5k 浏览
    J

    请问像这种狭窄处要划分边界层,怎样处理能保证网格高质量呢

  • 0 赞同
    26 帖子
    35k 浏览

    @King_RVM 你好,使用starCCM绘制网格,然后导入到openFoam里面跑算例,计算结果会出现质量不好的情况么?请教一下~

  • 0 赞同
    5 帖子
    7k 浏览
    李东岳

    这个其实是个好话题,很多大佬在sci里面也明确说了,LES之类的模拟,如果对于这种类似的圆管,施加cyclic边界条件,会对湍流产生抑制。但是我也没见过大佬们进行对比。

    另一方面,感觉边界条件,直觉感觉,单单的一个patch引起的差异,不至于那么大。LES模型参数的影响更大。如果谁有类似的文章,可以发来看看。

  • 关于OpenFOAM矩阵求解器

    Algorithm
    9
    0 赞同
    9 帖子
    18k 浏览

    @micro 网格数量越多,网格尺度越小,计算肯定更加不稳定,您可以查一下CFL数

  • 坛子里有FLUENT MHD计算相关的算例么

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

    ANSYS Fluent Magnetohydrodynamics (MHD) Module Manual

  • 0 赞同
    6 帖子
    9k 浏览
    I

    @东岳 在 船体在ICEM中划分网格时,水域和空气域怎么定义 中说:

    是的你不需要在ICEm里面定义,在CFD里面定义
    如东岳老师所讲,你绘制好网格的时候只是定义了一个名称,真正的定义是在fluent中完成的。
    建议看一下相关教程。

  • 0 赞同
    5 帖子
    6k 浏览
    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的。

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

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

    公告
    137
    0 赞同
    137 帖子
    165k 浏览
    F

    想问一下李老师这个计划还有吗!!很想要这本湍流模型:mianmo:

  • 0 赞同
    1 帖子
    2k 浏览
    袁宝强

    请问各位大佬,对于可压缩流大涡模拟的总湍动能如何定义呢?

    对于不可压缩,大涡模型采用的亚格子模型采用的雷诺时均的形式,因此可以用雷诺时均计算大涡湍动能和亚格子湍动能,两者的和就是总湍动能。
    $K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$ 但是对于可压缩流的大涡模拟,亚格子模型采用的为法夫尔时均形式,此时如何定义总湍动能呢?我的理解是
    $\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
    请问各位大佬,我的理解对么
  • 0 赞同
    6 帖子
    7k 浏览
    麦迪文

    @Dingcy 哈哈哈哈,我就瞎看看,我主要做CFDEM方面

  • 恭喜开版!

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

    都柏林大学发来贺电!

  • 0 赞同
    3 帖子
    4k 浏览
    DY大世界D

    @东岳 感谢李老师解答,我再试试

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

    reFineMesh可以加密一个方向,这个需要指定区域,你试试topoSet+refineMesh呢》

  • 比冲?

    Algorithm
    2
    0 赞同
    2 帖子
    3k 浏览
    H

    @liujuncfd 这恐怕不行吧,比冲是推力除以单位时间消耗燃料的质量(或重量)。

  • ANSYS FLUENT UDF学习

    Fluent
    12
    0 赞同
    12 帖子
    18k 浏览
    L

    @l-j刘侃 刘老师这个成本很高吧:chigua: