Skip to content
  • 为什么网格质量会影响收敛?

    Algorithm
    4
    0 赞同
    4 帖子
    11k 浏览
    S

    在评估网格质量的时候,一般常用下面4种标准进行衡量。
    1 均匀性(Uniformity)
    定义:网格从最小单元过渡到最大单元的速度(反之亦然)。
    网格单元之间尺寸变化越缓慢,则均匀性越好,结果越精确。
    一般来说,当整个计算域空间的所有网格都在同一尺寸时,均匀性是最好的。然而,在一些区域,所允许的最大网格尺寸非常小,因此,如果在整个计算域都采用这种尺度的网格将会导致过多的网格数目,这在实际应用中是不现实的。
    有了这个约束,在网格划分阶段,决定从哪里开始进行网格单元增长,以多大的速率增长就非常重要了,最终目的是在保证关键位置网格解析度的同时尽可能减少总的网格数目

    2 纵横比(aspect ratio)
    定义:下面两个值中较大的一个:
    网格单元边界框各个面中最大面积和最小面积的比值
    表达式的值
    捕获.PNG
    其中,ax、ay、az是网格单元边界框各个面的面积,V是网格单元的体积
    纵横比越接近1,结果越精确。
    这类网格质量问题最常见于两种情景:
    • 边界层网格
    • 六面体网格的各向异性细化(refinement)
    前者最为常见,为了捕捉近壁面处流动,需要在垂直壁面方向上布置非常精细的网格,尤其对于Yplus≈1的情形,有可能会导致纵横比高达几千。
    后者常见于网格加密,譬如在模拟船体兴波时,对于自由液面处网格的加密,将导致垂向上的网格尺度远远小于纵向上的尺度。

    3 正交性(Orthogonality)
    定义:连接两个相邻单元中心的直线和两个单元所共面的法线之间的夹角。
    越接近于0,精度越高
    在多面体网格中,正交性经常被提及,因为在有限的空间中,网格单元的形状可以有多种选择。正交性对于模拟精度具有重要影响,因为它与变量梯度、粘性和对流通量的计算息息相关。

    4 偏斜度(Skewness)
    定义:连接相邻两个单元中心的直线与所共面的交点,距离该面中心点的长度。
    偏斜度越接近于0,精度越高
    较为复杂的几何结构加上比较粗糙的网格,是导致偏斜度大的主要罪魁祸首。同样,该因素对于结果的精确性也有很大影响,它影响着通量平衡的计算。
    针对此因素,提高网格质量的方法一般是在几何较为复杂的地方,适当减小网格尺寸和增加网格密度。

  • 想问下这个张量公式怎么推导

    Algorithm
    3
    0 赞同
    3 帖子
    4k 浏览
    S

    个人建议去看一下北大吴望一老师的流体力学第一章,看了之后推导这类方程简直小菜一碟

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

    我错了,是$S_{ii}=0$。因为满足不可压缩流体的连续性方程。尴尬。。。:shangxue:

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

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

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

  • 32天回国了

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

    @zhanghan

    科研环境感觉还是国外好。

    但是国外好山好水好无聊。没有狐朋狗友也没有烤串火锅。语言也不通。同时毕竟是亚洲人,发展中国家,很多对中国人有歧视。日本人地位比中国人要高。

  • 0 赞同
    4 帖子
    5k 浏览
    bestucanB

    @jdzhang 你的颗粒会动么?不会动,为什么颗粒间还有距离?会动,那对着“多相流”查比对着“多孔介质”查更容易找到解决方案些。

    这有一个画网格的:
    https://www.bilibili.com/video/BV1eS4y1F7Wm/?vd_source=c01e8d4fff0db460c75b32e0d0191b24

    这有些相关的:
    https://www.researchgate.net/publication/308996206_Pore-scale_numerical_simulation_of_mud_erosion_in_the_subsea_sand-mud_alternate_layer_using_lattice_Boltzmann_method

    https://ascelibrary.org/doi/10.1061/(ASCE)HY.1943-7900.0001489

  • 一维Sod激波管问题

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

    其实一楼这个问题我还是有点没想明白。就是一个x方向的一维算例,出现了y方向的速度,应该如何理解..

  • 气液两相 没有纯气体

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览
    H

    DZ44GBIONA}7`WR})4VI(2V.png
    GU0962NTFA1OXF5XX{B5.png
    上面是加热壁面,下面全都是多孔区域(孔隙率0.73),然后我产生的液相分布是这样的~ 几乎没有纯气体区域~甚至在固液界面处的液相含量反而高~我已经调过不同的毛细力模型,发现不是毛细力的影响~
    现在还可以做的是增大相变系数~但是我相变系数一增大,就会在开始相变的时候发散~

    有没有老师可以分析一下 可能是哪里导致的问题~

    solvers { "alpha.liquid.*" { cAlpha 1; nAlphaCorr 2; nAlphaSubCycles 2; MULESCorr yes; nLimiterIter 5; solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0; } rho { solver diagonal; tolerance 1e-7; relTol 0.1; } rhoFinal { $rho; tolerance 1e-7; relTol 0; } mDotSmearFinal { solver PCG; tolerance 1e-7; preconditioner DIC; relTol 0.00; smoother DIC; } "(U|T)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-7; relTol 0.; minIter 15; maxIter 50; }; "(U|T)Final" { $U; tolerance 1e-7; relTol 0; maxIter 50; } p { tolerance 1e-8; solver PCG; preconditioner DIC; relTol 0.001; smoother DIC; } pFinal { $p; tolerance 1e-07; relTol 0; } "pcorr.*" { $p_rgh; tolerance 1e-9; relTol 0; }; } PIMPLE { momentumPredictor yes; nCorrectors 5; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 1e5; } relaxationFactors { alpha.liquid 0.2; U 0.5; p 0.5; } ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { div(phi,alpha) Gauss vanLeer01; div(phirb,alpha) Gauss vanLeer01; //div(phirb,alpha) Gauss interfaceCompression; //div(rhoPhi,U) Gauss vanLeerV; //div(rhoPhi,U) Gauss limitedLinearV 1; div(rhoPhi,U) Gauss linearUpwind; div(rhoPhiByEps,U) Gauss linearUpwind; div((mu*dev2(grad(U).T()))) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(((rho*nu)*dev2(T(grad(U))))) Gauss linear; div(rhoCpPhi,T) Gauss vanLeer; div((interpolate(cp)*rhoPhi),T) Gauss vanLeer; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; }
  • interFoam求解报错

    OpenFOAM
    11
    0 赞同
    11 帖子
    8k 浏览
    A

    @李东岳 您好,我刚刚解决了之前权限的问题,可是这次键入interFoam命令之后又出现了另一种错误,看提示好像是wmake出错?请问这又是什么原因,这个算例本身应该是没问题的。
    ![)K435ZJCL}TWILHXJP`O42.png

  • 导出网格数据

    Meshy
    6
    0 赞同
    6 帖子
    8k 浏览
    李东岳

    @禾斗石开 fluent

  • 0 赞同
    3 帖子
    4k 浏览

    @bestucan 非常感谢老师指导:xinxin: 我回去试一试

  • shm网格看起来不共节点

    Meshy
    5
    0 赞同
    5 帖子
    8k 浏览
    T

    @我是河滩 您好,请问一下这种网格在snappyHexMeshDict是怎么画出来呢?我现在正想画出这种网格,也就是让笛卡尔网格沿着我给的边界进行平分加密

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

    首先是CPU 但是CPU的优化性和效果是非常复杂的问题 一般来说 目前的主流 最好最大主频不要低于2.8 同时并行核心数不要低于32 从我的经验48核心比较好
    其次是内存 内存对一些内存要求多的算例的提升效果比较明显 能够占到30%左右
    硬盘一般跟软件面板操作相关 有的时候也是一些时间
    显卡的话 如果你是Fluent可以适用CUDA的方法 提升效果比较明显 两张4W的卡 大概coupled的下 提升30%左右
    最后就是操作系统 你的多核心下 有些操作系统能够提升速度几百倍(显然有些系统对于40以上的核心效果很差) 一般Linux下比同等的windows server要快5%

  • 管道内的圆柱绕流

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

    @zhaodl 这个目前我尚不清楚

  • 0 赞同
    8 帖子
    9k 浏览
    B

    @LuoS 在 snappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说:

    教程中似乎有个建筑风的案例:
    /Tutorials/inco

    好的谢谢 我试一下

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

    @东岳 对,现在打开了

  • 关于udf载入出错原因

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

    应该是,udf.h头文件无法读取

  • 0 赞同
    6 帖子
    7k 浏览
    H

    @李东岳 老师,我看了一下我的z盘是网络驱动器,所以读取相比固态硬盘会慢很多~也就是说如果我直接放在Ubuntu根目录下运行算例的话,用时相对会少~
    但是就有另外的问题:在根目录下运行时,对电脑内存是否有要求(实际占用的是哪里的内存?);后期后处理时,安装的win版本的paraview,好像打开不到wsl下~
    1640918395(1).png
    1640918235(1).png

  • 0 赞同
    2 帖子
    3k 浏览

    这里给出了一种方法:fluent导入形成流场的外部结构

  • 关于Fluent中壁面润滑力的UDF

    Fluent
    2
    0 赞同
    2 帖子
    4k 浏览
    waltenW

    这个问题 我自己想明白了,这个方向向量不需要指定,该udf使用DEFINE_SOURCE,在指定source时同时指定udf作用的相和方向即可。