Skip to content
  • 求关于LBM方法的算例或者教程

    Algorithm
    13
    0 赞同
    13 帖子
    19k 浏览

    @FluidGao 嗯嗯,感谢。

  • 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彩虹条
    2
    0 赞同
    2 帖子
    4k 浏览
    C

    $\mathrm{TV}=\sum_{i}^{}|\phi_{i+1}-\phi_{i}|$

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

    雷诺数从物理意义上说是惯性力比上粘性力,想问为什么雷诺数的定义会涉及到惯性力这么一个虚拟力,是怎么考虑的,对于管道流动来说,如果是匀速运动,则流体的压力和粘性摩擦力平衡,对于所研究的流体单元来说是不存在惯性力的,如果管道流动存在加速度,则流体的压力和粘性力的合力导致流体的加速,可以采用一个虚拟的惯性力来等效流体的加速,则可以认为流体在压力、粘性力和惯性力三者作用下平衡。 所以,雷诺数为什么没有压力项,而用一个虚拟的惯性力来定义?

  • bubbleFoam

    OpenFOAM
    18
    0 赞同
    18 帖子
    26k 浏览
    李东岳

    bubbleFoam已经消失很久了。然而作为双流体模型的首个求解器,并引以为傲的引入了Henry的速度Phase intensive equation。但是在后续版本中被惨痛抛弃。其中的Phase intensive equation在开发之后解决了一些相分离问题,但是被一些sci文章指出在某些情况下非守恒形式的方程会引起误差。并且在模拟工业规模的鼓泡床、萃取塔的时候,bubbleFoam没有考虑可压缩性(实际上床底和床顶压差很大),最终bubbleFoam被compressibleTwoPhaseEulerFoam(附加颗粒动理学模拟气固)完全取代。

    然而bubbleFoam依然是一个经典的求解器,或许身份等同于icoFoam。没人用,但他的意义不在于此。

  • STAR-CCM+网格转换后,checkMesh报错

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

    1.我猜是symmetry这个边界条件,试一下改成wall,然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查,外部软件的网格节点编号和blockMesh的逻辑不一样。
    2.可能是前后的empty边界问题,你可以在starccm里生成3D网格后,重新在openfoam里extrude来生成前后empty边界
    类似这样
    constructFrom patch;
    sourceCase "../w3-d_hc1-3";
    sourcePatches (symFront);

    // If construct from patch: patch to use for back (can be same as sourcePatch)
    exposedPatchName symBack;

    // Flip surface normals before usage. Valid only for extrude from surface or
    // patch.
    flipNormals false;

    //- Linear extrusion in point-normal direction
    extrudeModel linearNormal;

    nLayers 1;

    expansionRatio 1.0;

    linearNormalCoeffs
    {
    thickness 1;
    }

    // Do front and back need to be merged? Usually only makes sense for 360
    // degree wedges.
    mergeFaces false; //true;

    // Merge small edges. Fraction of bounding box.
    mergeTol 0;

  • pointwise输出二维网格到openfoam无法计算

    Meshy
    4
    0 赞同
    4 帖子
    6k 浏览
    A

    @zzcfd 改fvSolution,把p改成这个solver PCG; preconditioner DIC;试试

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

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

  • 当保安也需要硕士嘛

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

    @cfd-china 门槛啊,真高!

  • Fluent中的VOF与组分运输模型

    Algorithm
    1
    0 赞同
    1 帖子
    3k 浏览

    老师们,使用VOF能模拟气液两相流,使用组分运输能模拟同相态的质量传递。
    请问, 能用VOF与组分运输,模拟气液两相流的传质吗?

  • 一个包含CFD各种研究方向的图

    Algorithm
    3
    1 赞同
    3 帖子
    6k 浏览
    李东岳

    捕获.JPG

    重写了一个,以后用这个更新

    用的Adobe AI

  • 0 赞同
    9 帖子
    13k 浏览
    W

    gcc在Windows下的表现不如Linux。

  • fluent xy plot

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    S

    请问在fluent里怎么增加x,yplot的数据点,现在是每隔1微米一个点,怎么每0.1微米一个点?或者在别的后处理软件能实现吗?现在的网格也是1微米的。
    0b602e26-be9f-406e-bbae-aa49086f601a-image.png

  • 0 赞同
    10 帖子
    16k 浏览
    S

    @cfd-china
    有钱的话:cheeky: ,可以去构建一些关系嘛。但是这个也有点难。慢慢来嘛:happy:

  • 关于可压缩高速流求解sonicFoam求解器

    OpenFOAM
    3
    0 赞同
    3 帖子
    7k 浏览
    V

    谢谢东岳大神的回复,看了您推荐的,感觉有点明白了,还在挣扎中。

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

    @东岳 谢谢版主的解答。。这个问题暂时搁置了,以后有机会把UDF分享出来。

  • front tracking 怎么翻译?

    Algorithm
    4
    0 赞同
    4 帖子
    5k 浏览
    队长别开枪

    @东岳 不敢班门弄斧,我一直是你的粉丝:ok: :ok: :ok:

  • 修改AMI或者 GGI边界条件

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

    我想可能的话,可以让slave patch 以任意角度拷贝 master patch

  • DNS网格数量

    CFD彩虹条
    3
    0 赞同
    3 帖子
    12k 浏览
    C

    @xpqiu
    .................
    25T硬盘确实花了不少钱。:sunglasses:

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

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