Skip to content
  • vim tab设置为4个空格

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

    vim保存窗口:

    :mksession ~/.session.vim vim -S ~/.session.vim
  • 计算积分式的方法

    Algorithm
    1
    0 赞同
    1 帖子
    2k 浏览
    WeakFormW

    在破碎模型中,破碎的概率和子气泡分布一般涉及到积分式,比如说经典的Luo模型270626e3-d11f-4d47-a461-cdc3cbde1bb6-image.png 这里边openfoam官方直接用大佬提出来的化积分为多项式公式替代求解:e4c4c71f-d206-4464-bb15-5ee960dcdba8-image.png 有学者就提出来了这个exp内的判据不精准,使用表面能和韦伯数多重盘踞,然后这个积分式就这样了:
    6f926600-85cf-4d78-a936-ef166ab439ab-image.png 现在这样的话积分还怎么求解呢,因为min和max函数的存在原来的积分式公式还能使用么。还是需要分段求积。原论文是植入到fluent中大佬没写怎么处理的.

  • 请教RANS模型和LES模型的对比问题!

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

    所以LES 一定是优于 RANS 的模型,这样理解正确吗!!

    好像不能这么说。尤其你在写文章的时候。或者回复审稿人的时候。具体来讲,起码RANS算的很快嘛,有些sci,你看他们模型雷诺应力预测的也挺好的,

  • DPMFoam中颗粒时间步长在何处设置的

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

    明白了,谢谢东岳老师。

  • 求流体凝固为固体的算法

    Algorithm
    16
    0 赞同
    16 帖子
    31k 浏览
    J

    @李东岳 在 求流体凝固为固体的算法 中说:

    @北斗 最近在研究多孔介质。目前一种非常简单的凝固模型,就是利用多孔介质来处理。固体区域给定一个非常小的渗透率。

    血小板聚集形成血栓

    这种看起来像是拉格朗日粒子的黏连

    请问李博士,您的意思是通过多孔介质来实现凝固吗,如何实现,这方面资料可以发出来吗

  • 关于在linux下编译UDF

    Fluent
    18
    0 赞同
    18 帖子
    21k 浏览

    @徐小双 应该是没有删除干净的原因。

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

    该书11.8节(中文版305页,英文版411页)中的例2,推导源项的表达式时,高阶格式减去一阶迎风格式,一阶格式的值取上上游的值φu,为什么这里不取φc呢?先谢谢了。
    7a735b6e-01c9-4088-971a-359144da8189-image.png

  • 胖猫要去美国了

    CFD彩虹条
    43
    0 赞同
    43 帖子
    53k 浏览
    队长别开枪

    @东岳 https://www.linkedin.com/jobs/view/postdoc-in-cfd-of-turbulent-environmental-flows-at-university-of-houston-1628249421
    我离青千条件远着呢,先在这边混着,等我老婆博士毕业了再说:chigua:

  • 一个CFDer的日常

    CFD彩虹条
    13
    0 赞同
    13 帖子
    15k 浏览
    Z

    几年来一直想做而又没时间做的事,老铁这句话好扎心啊

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

    c7a27198-c865-41aa-bc5f-ce28c85f32d3-1728486466827.jpg

    请教各位大佬专家,立方体里边,挖掉一个圆柱体,对圆柱顶面进行切割,上半部分为流体入口边界,请问又什么好的网格划分思路吗?

    单纯立方体里挖一个圆柱能够实现,但中间多切了一条线,不知道怎么操作比较好。

    另外,有没有好的网格划分方面的教程资料推荐?

    恳请大佬们指点。

  • 一些勘误

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

    谢谢,已更新

  • 关于Fluent报错

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

    @李东岳 谢谢李老师回复,前几天关了一下,没起作用。我用其他的服务器跑了一下Fluent。

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

    顺便测试一下侧边栏回复

  • 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的。

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

  • 0 赞同
    8 帖子
    13k 浏览
    M

    QQ截图20190911102057.png

    再补充一下,用不同的电脑跑了 damBreak 算例,得到的三个 ExecutionTime 和 ClockTime 。 并不总是 clockTime 大于 ExecutionTime 。@wwzhao

    但是仔细检查发现,造成右边两个有差异的主要原因是两个时间的精度问题,如下所示:

    ExecutionTime = 1.58 s ClockTime = 1 s ExecutionTime = 1.59 s ClockTime = 1 s ExecutionTime = 1.59 s ClockTime = 1 s ExecutionTime = 1.6 s ClockTime = 1 s ExecutionTime = 1.6 s ClockTime = 1 s ExecutionTime = 1.61 s ClockTime = 1 s ExecutionTime = 1.62 s ClockTime = 1 s ExecutionTime = 1.62 s ClockTime = 1 s ExecutionTime = 1.63 s ClockTime = 1 s ExecutionTime = 1.64 s ClockTime = 1 s ExecutionTime = 1.64 s ClockTime = 2 s ExecutionTime = 1.65 s ClockTime = 2 s ExecutionTime = 1.66 s ClockTime = 2 s ExecutionTime = 1.66 s ClockTime = 2 s ExecutionTime = 1.67 s ClockTime = 2 s ExecutionTime = 1.67 s ClockTime = 2 s ExecutionTime = 1.68 s ClockTime = 2 s

    因为是一个很小的单核算例,IO也不大,所以IO基本上没有花时间,导致在保留精度不同的情况下出现了上述差异。

    多说两句:
    上图右上角是在WSL2 根目录下的结果,
    右下角是在一个 Ubuntu18.0.4 服务器上的结果,
    左边的结果是用 WSL2 在 /mnt/h 盘下的结果,

    可以看到在即使是很小的算例。这个IO耗时也是相当大。这是 WSL2目前的问题,相较于WSL 它在 根目录下确实IO能力有很大提升。但是在 /mnt/* 下 真是果断弃疗。

    版本:
    OpenFOAM-v7

    The 'ExecutionTime' is the elapsed CPU time, and the ClockTime is the elapsed wall clock time for the latest time step(approximate!!)

    https://www.cfd-online.com/Forums/openfoam-solving/67695-difference-between-executiontime-clocktime.html

  • 无锡超算的单核计算效率对比

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

    无锡超算,用openfoam计算 :

    国产申威单个核组(仅1主核,无从核加速),与intel商用节点的单核计算效率,有谁比较过吗?

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

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

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

  • 0 赞同
    3 帖子
    3k 浏览

    上传一下我的snappyhexMeshDict文件,请各位老师帮忙看看```
    code_text

    ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 10 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { format ascii; class dictionary; object snappyHexMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg" castellatedMesh true; snap true; addLayers true; geometry { cylinder { type triSurfaceMesh; file "cylinder.stl"; } refinementBox { type searchableBox; min (0.0 0.0 0.08); max (1.2 0.2 0.22); } }; castellatedMeshControls { maxLocalCells 200000; maxGlobalCells 5000000; minRefinementCells 10; maxLoadUnbalance 0.10; nCellsBetweenLevels 6; features ( { file "cylinder.eMesh"; level 3; } ); refinementSurfaces { cylinder { level (3 3); patchInfo { type wall; } } } refinementRegions { refinementBox { mode inside; level 1; } } insidePoint (0.4 0.1 0.15); } snapControls { nSmoothPatch 3; tolerance 2.0; nSolveIter 30; nRelaxIter 5; nFeatureSnapIter 10; explicitFeatureSnap true; implicitFeatureSnap false; multiRegionFeatureSnap false; } addLayersControls { layers { "(cylinder_stl_face<stlunit=MM>|cylinder_stl_top<stlunit=MM>)" { nSurfaceLayers 10; } } relativeSizes true; expansionRatio 1.05; finalLayerThickness 0.8; minThickness 0.5; } meshQualityControls {} writeFlags ( // scalarLevels // layerSets // layerFields ); mergeTolerance 1e-6; // ************************************************************************* //
  • 0 赞同
    3 帖子
    5k 浏览
    J

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

  • lagrangian颗粒信息统计

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

    @东岳 controlDict里没有统计颗粒信息的设置,只有统计流场信息的设置。颗粒信息的计算和输出都在src/lagrangian/intermediate/submodels/CloudFunctionObjects/中。