Skip to content
  • 高等数学的CFD相关问题

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

    页数按照pdf的页数,不是书上的页数

    上册:

    https://www.jianguoyun.com/p/DVcobNoQ9s3ZBhjejV8

    73页【直线运动的速度】:导数的直观描述,CFD方程这面全是导数,大部分为偏导数,后面会介绍 75页、76页【公式1-4】 138页【公式3-3】:泰勒公式,CFD很多方程都是基于泰勒公式推导出来的

    下册

    https://www.jianguoyun.com/p/DQnanmcQ9s3ZBhjijV8

    23页:流量计算公式 71-75页:偏导数,CFD这面都是偏导,重点看一下,理解一下 77页提到了个拉普拉斯方程,随便扫一眼 78页全微分,跟CFD的物质导数关联 79页3-3 81页有个例子 112页介绍了梯度,但是感觉描述的不是很简单,可以看看 141-144页 重积分,为之后理解有限体积法做准备 224-225页 曲面积分的定义,为之后理解有限体积法做准备。其中225页有个图,比较形象 230-231页,看高数中对通量的定义 243-245页,高数中明确的对速度散度的定义
  • 计算对流扩散方程没有扩散效果?

    Algorithm
    8
    0 赞同
    8 帖子
    13k 浏览
    李东岳

    不好意思,那个算例里面求解器我这面编译不了(openfoam-4.x and 5.x).. 估计后来问题已经解决了吧

  • ICEM导出inp文件出错?

    Meshy
    2
    0 赞同
    2 帖子
    4k 浏览
    J

    导出时候选择了如下这些选项
    image.png

  • 压力驱动的方腔流动

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

    最近计算一个方腔流动遇到一些问题,这个流动沿流向(x方向)是周期性边界,yz方向是壁面,流动固定x方向的平均速度,自动求解压力梯度,在采用rans计算时是正常的,可以计算出一个合理的结果

    pressure gradient = 0.179097

    但是当采用PANS(自己编译的湍流模型)计算的时候,压力梯度会变得非常小,导致计算结果很差,

    pressure gradient = 0.0903267

    自己编译的模型可以在channel395中正常计算,但是在该算例中出现压力的降低,不知道是什么原因引起的

  • 0 赞同
    2 帖子
    3k 浏览
    bestucanB

    https://www.cfd-online.com/Forums/openfoam-meshing/133850-snappyhexmesh-problems.html

  • ANSYS Fluent中QMOM的植入缺陷

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    S

    @东岳 谢谢:baobao:

  • ANSYS 19.2 Linux安装问题

    Fluent
    2
    0 赞同
    2 帖子
    4k 浏览
    bestucanB

    c运行库太老了,新的接口,软件用了,系统还不认识。更新C库可以解决。
    按说这库装哪都可以,家目录也行,只要配置好环境让软件找到就行。
    实际操作比较复杂,这种操作个性化太强,没有手把手式的教程,自己配方案,只能正着学,不能倒着学。倒着一堆坑,磨死人

  • 均匀各项同性湍流的初场设置问题

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    B

    @shifc 不知道你现在解决了没有,我最近做格栅湍流也碰到了同样的问题,我的解决方法是谱空间的复速度对三个坐标轴必须都是共轭对称的。

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

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

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

  • 使用KomegaSST模型做湍流模拟的y+值

    OpenFOAM
    1
    0 赞同
    1 帖子
    3k 浏览
    I

    对于使用KomegaSST模型作低雷诺数的湍流管流模拟,当y+值小于1时,是不是可以认为网格足够精细,无需再细化?
    另外k omega nut三个量大概在哪个范围内合理?我的结果大概是以下:

    k 10的负3次方

    omega 最大50w+,最小500左右
    nut 10的负10到负5次方

  • forAll boundary

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

    应该是Sb.boundaryField() 返回const reference
    Sb.boundaryField().ref() 返回 reference
    根据我的经验是这样,如果错误的话把编译错误贴上来
    认真看看代码很简单的,编译器的错误已经告诉你所有信息了!

  • ICEM划分网格

    Meshy
    13
    0 赞同
    13 帖子
    19k 浏览
    H

    @Junren-Hou 在 ICEM划分网格 中说:

    可以试试check block,有些时候其实没有问题,就它软件自己的毛病,然后block里最后有一个check block,可以让他自己检查一下。里面有好多选项,都试试,什么check/fixed block

    回忆了一下,我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域,一方面追求地面边界层,致使网格尺度过小,一方面为了平衡网格数量,过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

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

    @队长别开枪

    谢谢您的建议,感觉这个方案可行。 我尽早试试,给您答复。。。

    谢谢

  • 0 赞同
    2 帖子
    3k 浏览
    L

    《A Thorough Description Of How Wall Functions Are Implemented In OpenFOAM》一文2.6节描述了nutWallFunction,nutAtmRoughWallFunction中nutw公式与之相同,其他量有改变,我再看看。

  • 直升机炫图一张

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

    替代文字

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

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

  • 质量入口概念问题

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

    @a602657563 谢谢同志解答,你讲的很精炼,我同意你的观点,应该是这个笔者自己没弄清楚入口的概念问题,入口要么就是指定压力,要么就是指定速度,而且你讲的initial pressure很准确,只是给定一个初始条件值,而不是计算初始值。

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

    额,自问自答?

  • 动网格设置两块旋转域问题?

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

    如果想要设置两块旋转域(对转),dynamicMeshDict应该怎么编写呢?我像下面那样编写,结果只有sector_part2旋转。

    dynamicFvMesh dynamicMotionSolverFvMesh; motionSolverLibs ("libfvMotionSolvers.so"); motionSolver solidBody; cellZone sector_part1; solidBodyMotionFunction rotatingMotion; origin (0 0 0); axis (0 0 1); omega 112.573737; // rad/s cellZone sector_part2; solidBodyMotionFunction rotatingMotion; origin (0 0 0); axis (0 0 1); omega -112.573737; // rad/s

    正确的编写形式应该是怎样的呢?

  • 关于气液分离的采出分率问题

    Algorithm
    2
    0 赞同
    2 帖子
    4k 浏览
    C

    补充一下,不知道这些图中的这些数据点是如何提取的,谢谢大家