Skip to content
  • 尖角附近的网格过度

    Meshy
    11
    0 赞同
    11 帖子
    12k 浏览
    I

    @HITSC30 icem不是很熟练,我去试一下

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

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

  • 最近固态硬盘空间不够

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

    把ubuntu装在了机械硬盘。

    用了大约4年的固态硬盘,忽然发现难以忍受机械硬盘。用了1星期之后,又装在了固态硬盘。只有40G大小。

    现在固态硬盘便宜了,250G的640元。

  • 并行计算入门求解

    Algorithm
    6
    0 赞同
    6 帖子
    11k 浏览
    Q

    @Kingkong 搜如下教材

    高性能计算并行编程技术_MPI并行程序设计_[都志辉]:非常适合入门
    Using MPI :结合上一个教材一起看
    Iterative Methods for Sparse Linear Systems [ Yousef Saad ] :介绍 deomain decomposition 的数学理论

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

    欧拉拉格朗日模型中离散型用ODE方程描述,追踪每个粒子,这和PBM模型是冲突的,PBM是数量密度函数的模型。

    目前还没见到过文献讨论这个 :cheeky:

  • 0 赞同
    6 帖子
    6k 浏览
    Z

    @李东岳
    当边界不做判断直接叠加之后是这样的
    fc3c9e32-1d39-4a41-b44a-85c2594da99b-image.png

    bottom
    {
    type codedFixedValue;
    value uniform 101325;
    name face1P;
    code #{
    scalarField& field = *this;
    const label id = this->patch().boundaryMesh().findPatchID("bottom");//在网格中找到inlet
    const volScalarField& pp = this->db().objectRegistry::lookupObject<volScalarField>("p");
    const scalarField& ppp=pp.boundaryField()[id].patchInternalField();
    forAll(field, faceI)
    {

    field[faceI]=ppp[faceI]+100; } #}; }
  • windows 10下可以跑linux!

    CFD彩虹条
    50
    0 赞同
    50 帖子
    99k 浏览
    S

    我用子系统上的OpenFOAM,总感觉并行计算没用,和我不开并行计算时间差不多。。。而且还decomposePar不了,好像是安装的时候少了啥东西。。。

  • piso和pimple

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

    @李东岳 在 piso和pimple 中说:

    数值稳定性

    李老师您好,我在用buoyantPimple求解器的PSIO模式计算流场压力、密度和速度出现明显的数值振荡,就是随时间步长周期性振荡发展,但是振幅确实是一直在减小的,可否理解为这就是数值振荡问题,通过调整PIMPLE求解器外部动量预测修正次数转化为PIMPLE增强稳定性?
    计算时序图如下:
    8055ef71-35b7-40e6-b692-ac6aa27a4f2c-image.png
    58fdb9e2-331b-4f6b-9a4a-39cb7b94f1bf-image.png
    064a77b7-1f57-4077-9480-e6395c761979-image.png

  • 每日一问:有关icem网格检查失败问题

    Meshy
    1
    0 赞同
    1 帖子
    2k 浏览
    2

    icem网格划分导入fluent时检查失败,显示Error: interface zone 3 has two adjacent cell zones.
    4f68eb65-132e-4314-b0f6-6964dc466673-image.png

  • 猪肉降价了

    CFD彩虹条
    2
    0 赞同
    2 帖子
    3k 浏览
    DY大世界D

    喜大普奔,北京市民发来贺电:chigua:

  • 水流的空化问题

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

    实验室有个试验台,里面充满水,水流速度最大60m/s,无叶轮,最近要做单相流实验,想确定增压泵打多少压力水才不会空化。按照我的想法,100℃水温时,水的汽化压力为0.1MPa,要保证不空化,需要线速度最大的地方局部压力达到0.1MPa,按照伯努利方程,需要加上动压,这个时候发现系统压力居然要达到2.8MPa。仔细想想应该是不适用伯努利方程的,所以想问一下,有没有什么好的方法可以估算系统压力呢?

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

    @李东岳 老师,如果不使用refineMesh的话,结果会好很多,而且和论文中的Cd, Cl指数可以基本对上,但是在overset的边缘处还是会有一点问题。
    648bb924-7f1b-4e8c-a39d-48cc7c40906e-image.png

    我的oversetInterpolation的method使用的是inverseDistance.

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

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

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

  • flunet udf相关问题

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

    各位大佬,在编译udf出现问题,为什么会出现解码失败呢,flunet不是使用UTF-8 解码,为什么会出现gbk' code,从而导致Error at host: open_udf_library: couldn't get udf_data
    ???????????????

    UnicodeDecodeError: 'gbk' codec can't decode byte 0x8f in position 155: illegal multibyte sequence:

    注:路径检查没有中文和非字符

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

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

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

  • Ubuntu清华镜像

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

    https://mirrors.tuna.tsinghua.edu.cn/#

    下载速度很快

  • 0 赞同
    2 帖子
    4k 浏览
    Q

    @yuan_neu 建议看看《张量分析》的相关课程
    二阶 【各向同性】的 【张量函数】 ( Dij(r) )
    只可能是上述 (6.25) 的表达形式

    所以,你首先要明白“张量”,其次再搞清楚“张量函数”的定义,然后需要理解“各向同性”的概念,最后才能明白 6.25 式是怎么推导出来的

    实际上,NS 方程的 Newton 本构关系,也是 “可以” 基于 【线性】和【各向同性】的假设直接写出来的 !

    推荐看看清华黄克智的《张量分析》教材

  • 关于谱元法和有限体积需求网格的区别

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

    @XieXiaoyang 好文章,我去研究研究

  • 请教waves2Foam中的非定常风怎么使用?

    OpenFOAM
    19
    0 赞同
    19 帖子
    22k 浏览
    C

    @LGY 请问一下你这边能够在wave2Foam中生成非定常风了吗

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

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

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

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

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

    恳请大佬们指点。