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

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

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

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    594k 浏览
    H

    我写 固我在,嘿嘿

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

    这个要是openfaom的网格,就好处理了。2亿网格确实太多了,icem我也只画过最多几百万的网格。卡的都不行。你这2亿网格太疯狂了,机器内存应该很大。

  • interphaseChangeFoam

    OpenFOAM
    3
    0 赞同
    3 帖子
    4k 浏览
    Z

    谢谢...

  • 请教下离心泵的OF问题

    OpenFOAM
    3
    0 赞同
    3 帖子
    6k 浏览
    张晓宇

    @xpqiu 十分感谢!我会仔细根据您提的意见学习的,谢谢!

  • 0 赞同
    11 帖子
    11k 浏览
    ZY-GONGZ

    @李东岳
    已经看到了,非常感谢您:xiexie:

  • 0 赞同
    4 帖子
    6k 浏览
    C

    @李东岳 李老师有没有兴趣看下风工程这块的有个商业软件WT,里面是速度-压力耦合求解器,真的非常快非常快。我看了下,软件底层是从PHOENICS来的,然后搞了一套叫MIGAL的算法。

    现在CFD用的SIMPLE这种segregated solver确实太慢了,城市这种尺度下的流动本身大概率又比较简单,不可压,也没有化学反应之类的,可以暂时不考虑温度?可能会有颗粒流的问题,但是不是也可以先只考虑流动对颗粒的影响。这种情况下感觉流场本身的求解其实就看求解器的效率了。

    我搜了下,OpenFOAM这边其实已经有一些coupled solver 的工作,不知道成不成熟。

  • 如何处理奇点

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

    如果a/b 是一条曲线, b有的时候为0

    要根治的话,就是让b不为0

  • snappy画船体网格

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

    是的

  • 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; // ************************************************************************* //
  • CFD和LBM的优劣势

    CFD彩虹条
    6
    0 赞同
    6 帖子
    7k 浏览
    C

    @leo_chou 大批量的工业算例(本身可能不复杂,但就是工况多),LBM是不是也需要大规模网格和较长时间?工业应用CFD对时间的要求基本上是小时级别的。

    另外好奇问下,复杂几何下的LBM如何,需要湍流的LBM呢?我觉得大部分案例,复杂几何+湍流模型是标配,没有湍流的话可能有多相或者化学反应啥的,这些LBM是不是会更强一些?

  • 关于toalpressure边界条件的一点疑问

    OpenFOAM
    2
    0 赞同
    2 帖子
    5k 浏览
    李东岳

    @张晓宇 在 关于toalpressure边界条件的一点疑问 中说:

    在C++ documentation里面看的toalpressure的使用,有一个参数gama是气体热力学的性质参数,而且是必须要给的数值。这个意思是不是就是说toalpressure只能用于气体呢?

    总压方程,例如$p_0 = p + \frac{1}{2}u^2$来自伯努利方程。对于马赫数比较小的情况,适应于上述方程。当马赫数较大的时候,接近音速或者超音速的时候。则变为可压伯努利方程。这些在空气动力学书里面一般都可以找到。一般气体可以看错可压的,大部分情况下液体一般看做不可压缩的。并且,载推到可压缩伯努利方程的时候用到了理想气体状态方程。所以,totalPressure的可压缩部分用于气体比较合适,但是不能说不能用于液体。

    这个p0是相对压力还是绝对压力呢?

    公式中的p0为总压

    离心泵进口条件有全压的设置,这个全压能用toalpressure么?

    什么是全压?总压?可以把进口设置为总压。

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

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

  • snappy加边界层

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

    这个基本没有一针见血的办法,shm生成边界层已经是老大难了

  • M1制作免费ubuntu虚拟机

    CFD彩虹条
    8
    0 赞同
    8 帖子
    9k 浏览
    R

    感谢分享。不过虚拟机的话用UTM(免费)很方便,直接加载官方镜像安装即可。我在UTM里的Ubuntu20.04(Focal)里成功编译了OF8,9和OF v2012。

  • UDF宏自变量命名问题

    Fluent
    5
    0 赞同
    5 帖子
    6k 浏览
    E

    @深蓝 嗯嗯,现在已经明白了,谢谢 :xiezuoye:

  • openfoam中多孔介质多相渗流求解器

    OpenFOAM
    8
    0 赞同
    8 帖子
    13k 浏览
    lushiferL

    @Wayne 多谢多谢,研究一下。有啥不懂的再向大神请教吧。

  • 0 赞同
    1 帖子
    2k 浏览

    例如冬天拎起一桶水旋转,在旋转过程中水慢慢结冰,想要采用fluent中Solidification&Melting模型和Moving Reference Frames模型来模拟这一过程。
    跟仅开启Solidification&Melting模型的模拟不同,但似乎收敛性很有问题。凝固/融化模型与MRF是否可以联合使用?有谁做过类似的模拟吗?

  • 有谁用过blockMesh的mergePatchPairs功能吗

    OpenFOAM
    6
    0 赞同
    6 帖子
    11k 浏览
    张晓宇

    @xpqiu 又学到了 非常感谢!

  • lookupObject问题..

    OpenFOAM
    4
    0 赞同
    4 帖子
    6k 浏览
    W

    @mark alpha1_ 代表VOF模型中的体积分数,是一个volScalarField类型的对象。