Skip to content
  • 一个关于K方程SGS模型的问题:

    Algorithm
    5
    0 赞同
    5 帖子
    11k 浏览

    谢谢,那个我也看了(但没仔细看,只是特意找了扩散项)。应该也是和Menon一样的结论。因为Menon的文章截图比较方便,才引用他的

  • postProcessing of particle calculation

    OpenFOAM
    5
    0 赞同
    5 帖子
    7k 浏览
    Y

    @星星星星晴 谢谢!!

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

    @chpjz0391
    非常感谢,已更新。:sunglasses:

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

    @搬运工

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

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

    使用外挂包; 潜心研究稀疏线性系统求解理论;我建议花大量时间研读:YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:
  • 关于divDevRhoReff(rho,U)的疑问

    OpenFOAM
    3
    0 赞同
    3 帖子
    5k 浏览

    哦好像是定义没有弄明白,miu和nu是不一样的····

  • [分享] learn-cpp 的笔记模板 (org-mode)

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

    好像是上传格式的要求, 我用的操作系统和浏览器:

    win-10 64 Chrome Version 70.0.3538.110 (Official Build) (64-bit)
    win-7 64 Chrome Version 版本 70.0.3538.110(正式版本) (64 位)

    无效文件格式,允许的格式有:.png, .jpg, .bmp, .zip, .pdf, .tar, .tar.gz, .tin, .doc, .jpeg

    修改了几次后缀名: txt; org; 无后缀;以及一个压缩包形式:
    [0_1544141948952_README](正在上传 100%)
    [0_1544142047653_README](正在上传 100%)
    [0_1544142071028_README.txt](正在上传 100%)
    0_1544142322340_README.zip

    只有 zip 可以上传。

  • Fluent气固两相流边界条件求助

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

    @金钊 fluent推荐的颗粒初始体积分数最大是0.63

  • 清网行动/CORE成立

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

    @东岳 ok 已发送

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

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

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

  • 傅里叶变换求偏微分方程组

    Algorithm
    6
    0 赞同
    6 帖子
    12k 浏览
    WayneW

    @cfd-china 写文章,做实验。。

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

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

  • 有关紧密排布圆柱的网格绘制

    Meshy
    10
    0 赞同
    10 帖子
    13k 浏览

    @李东岳 感谢李老师,我试试看能不能画出来

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

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

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

  • 土柱CT扫描切片,用of画土壤孔隙网格

    OpenFOAM
    12
    0 赞同
    12 帖子
    15k 浏览
    Jack QQJ

    @wangxiaowei hello,我可能理解错了你的问题,我先用imageJ软件从ct扫描切片中提取出土壤介质,生成了stl文件,然后用snappyHexMesh画出孔隙的网格

  • 请问今年的培训有安排了吗

    已移动 CFD彩虹条
    3
    0 赞同
    3 帖子
    4k 浏览
    S

    @东岳 我报名,一位老师两名学生,稍后把报名表发给你

  • 前处理,进口面问题

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

    icem也有这个功能,只是我没用过,

    https://cfd-china.com/topic/6069/搅拌反应器icem结构网格划分方式/7

    Screenshot from 2024-09-28 14-40-10.png

  • RANS计算结果和理论对应

    Algorithm
    3
    0 赞同
    3 帖子
    6k 浏览
    Elibathe康E

    为何无法插入图片?我试着把公式截图发上来,但是好像不行:sad:

  • PBE中Liao-breakage 模型求助

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

    @CFD_xiaobai Luo破碎模型不是也是一个大公式么?

    破碎这面我没太看中文的书,跟你的描述可能对不上。据我了解,破碎模型有俩种:

    一种是breakge frequency(破碎频率)和daughtor size distribution (DSD)分开的, 一种是破碎频率+DSD写成一个大方程的;

    比如Laakkonen是第一种,Luo是第二种,Liao也是第二种。是这样么?

  • 有人知道这两个的区别吗?

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

    上面那个是欧拉框架下的偏导数,U为关于x,y,z,t的函数
    下面这个是拉格朗日下对速度的倒数,其中U为t的函数

    http://dyfluid.com/docs/theory.html
    公式35 and 36

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

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