Skip to content
  • 0 赞同
    7 帖子
    8k 浏览
    Y

    @fu 在使用laminar(Quasi-DNS)计算的时候,壁面处的网格Y+<1 ,且整个区域网格也要相对较小;当使用湍流模型的时候,壁面处理跟laminar应该不同,使用壁面函数之类的;如果题主使用一套网格应用于两种流动模型,是否会导致计算结果的偏差呢?还请大佬指点

  • 0 赞同
    6 帖子
    15k 浏览
    S

    @硫酸亚铜 在 满足不了Y+值的计算结果是否一定不准确? 中说:

    Adapt->Y+/Y*

    thanks:happy:

  • 0 赞同
    3 帖子
    4k 浏览

    @氕氘-川 批量保存的脚本老报错: ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
    vtkCSVWriter (000001F6DAFC08C0): Unable to open file: /-0.1150.csv
    ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
    vtkCSVWriter (000001F6DAFCC8C0): Unable to open file: /-0.0650.csv
    ...

    这是我的脚本:
    import numpy as np
    from paraview import numpy_support as ns
    import paraview as pv
    import paraview.simple as ps

    afoam = FindSource('a.foam')

    for y in np.arange(-0.115,0.115,0.05):

    # Make a cut of the wing to have a local 2D profile Slice1 = Slice( Input=afoam ) Slice1.SliceType = "Plane" Slice1.Triangulatetheslice = 0 Slice1.SliceOffsetValues = [0.0] Slice1.SliceType.Origin = [0.0, y, 0.0] Slice1.SliceType.Normal = [0.0, 1.0, 0.0] descriptiveStatistics1 = DescriptiveStatistics(Slice1) descriptiveStatistics1.VariablesofInterest = ['T'] descriptiveStatistics1.AttributeMode = 0 Render() writer = CreateWriter(str(y)+".csv",OutputPort(descriptiveStatistics1, 1)) writer.FieldAssociation = "Cells" writer.UpdatePipeline() del writer #SaveData(str(y)+".csv", proxy=OutputPort(descriptiveStatistics1, 1), FieldAssociation='Cells') descriptiveStatistics1.UpdatePipeline()
  • 0 赞同
    4 帖子
    5k 浏览
    李东岳

    替代文字

    你说上面这个图。这个算力对应的不是这个图。所以不一样。不过如果你增加近气速度的话,应该也是可以的。你加大10倍试试。

  • Fluent瞬态模拟后处理

    Fluent
    5
    0 赞同
    5 帖子
    9k 浏览
    Q

    @hhp886 好的,一会我试试看,谢谢啦~~

  • When I meet god

    CFD彩虹条
    2
    0 赞同
    2 帖子
    5k 浏览
    Elibathe康E

    也许再过一个世纪,湍流就不需要上帝来解答了,哈哈

  • 发了一晚上邮件

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

    头晕眼花

    sad
    :crying:

  • 速度标准差

    OpenFOAM
    3
    0 赞同
    3 帖子
    3k 浏览
    疏影横斜水清浅

    @李东岳 好的,谢谢老师!具体的公式我也在查,文献中没有直接给出

  • codedFixedValue边界条件

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

    @李东岳 好滴好滴,非常感谢,我这就试试:xinxin:

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

    b75ecf32b61ea077acecbb62002a6c8.png
    这是我后来自己的解释

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

    @李东岳 多谢东岳老师,重新划分网格后可以运行了

  • 讨论Fluent中的各种压力

    Fluent
    3
    0 赞同
    3 帖子
    8k 浏览
    S

    @cfd-china 哈哈,没有多久不见吧:cheeky:

    操作压力(operating pressure)是不是就是参考压力? 0_1466035493090_2016-06-16_080420.png 我在模拟的时候,都把操作压力设置成了零,这样得到的压力都是绝对压力?不用怎么换算了,这样没事吧:confused: 《ANSYS Fluent理论指南》已完成约90%,同步校对中什么时候能够完稿呢?
  • 引自CFD-Online.com一个大牛的邮件:

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

    此人在2005年左右淡出CFD-online.com,个人和他有一些私交,谈到了这个话题:为什么不活跃与论坛了?我们从您的帖子中学到了很多东西。他如下回答(翻译):

    活跃于论坛是非常耗费时间,许多问题都是重复的,越来越多的人想不劳而获的获取帮助,他们也不愿意为帮助他们的人提供一点好处。CFD不是我的业余爱好,CFD是我的职业,我需要靠CFD谋生。提供免费的服务不可持续。

  • ICEM三维结构网格正交性检查失败

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

    那就在做O形网格的时候,“O穿”整个区域。就是假设那个圆柱两端无限延伸,对这个无限长的圆柱做O形网格。

    如果不知道空心实心怎么选。做的时候;part里有个part,是专门放已经被删除的block;把这个part打开。就当是个全实心的圆柱。完了再关了那个part。

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

    没想到图越来越多。手机看的注意了。我这个图多。

    结构网格的拥护者简单说一下为什么我喜欢结构网格。以下只针对简单几何以及不是太复杂的几何。复杂到不行的网格我选择混合网格。我自己做的算例不用纯四面体网格。

    六面体网格具有以下优点:

    非矩形网格的歪斜(skewness)严重
    这需要做面法相梯度计算修正。稍微增加计算步骤。著名的CFD非矩形修正因此而来

    截面重构不光顺
    举例,一个正方盒子,一半是水,我们用矩形网格重构的气液截面是这样的:
    alt text
    黄色为一相,白色为一相。即使网格很糙,我们的截面也是平的。下图为三角形非结构网格:
    alt text
    可见糙网格下界面不平。这导致后处理需要耗费精力。

    高度非线性物理问题影响矩阵对角占优
    这也是非线性动力学,例如金属撞击、生物力学分析中结构网格肆虐的原因。(这个我没有测试,之前在网上看到的结论,来自链接文本

    非结构网格数量更多增加计算资源
    网格单元的面越多网格越少。六面体网格并不是数量最少的。数量最少的是多面体网格(star推荐的?)如下图
    alt text
    再看这个情况对比:只有第一个是结构网格,第二个网格多。第三个第四个和均为非正交且和流线不贴合
    alt text
    alt text
    因此目前自动生成的非结构网格一般数量是即为巨大的。计算的时间会超过你花一星期做结构网格的时间。另外,万一中途发散了怎么办。

    费结构网格高度不可控
    几何稍微一点点的变化会导致全场四面体网格做位移。结构网格(ICEMblock切出来的)具有即为即为高度可控性。可以对基本所有的网格单元进行控制。俩者区别天差地别。

    我了解的所有的离散格式对结构网格适应性更强
    也就是说结构网格对流场具有强的多的顺应性。例如旋转机械,结构网格和速度流线高度贴合。

    矩形网格对边界更好的捕获
    为何非结构网格要添加边界层网格?矩形网格对壁面贴合更好。

    复杂CFD模型在结构网格上收敛更强
    个人经验。比如双流体模型,LTS局部时间步框架下的模型,激波捕获等。当然那种不可压缩单相流不算

    太多就不写了...

    引用从业CFD27年OpenFOAM创始人Henry的一段话:

    hex and polyhedral cells are generally superior to tet cells in finite-volume and the difference in both accuracy and convergence behaviour can be large. However, it also depends on the quality of the hex and tet cells and many of the problems with tet meshes stem from the poor quality of some tets produced by many automatic tet mesh generators. In order to really appreciate the difference it might be a good idea for you to try both for your problems. For VOF calculations I would expect the solution on hex/poly cell meshes to be dramatically better.

    大意为:

    在有限体积中,六面体和多面体网格相对于四面体网格优势明显。主要在精度和收敛上。对于VOF(我上文提及的第2点),结构网格的优势是大大的。

    最后,需要耗费大量时间生成的六面体网格至今没有消失,必然有存在的道理。

  • ICEM网格加密

    Meshy
    3
    0 赞同
    3 帖子
    7k 浏览
    G

    @ycx1042341217 @huningdong 请问两位,壁面第一层网格高度可以做到什么程度?我发现最小只能设置到大约5e-5m:136:

  • 0 赞同
    11 帖子
    14k 浏览
    liujunCFDL

    @winsway_zero 前辈,请问gmesh怎么划分block

  • ICEM导入OpenFOAM遇到问题

    Meshy
    3
    0 赞同
    3 帖子
    5k 浏览
    K

    @bestucan hhhh谢谢您,我发现我哪里错了,我画网格的时候是xoz面上的,正确的应该是在xoy画个2D的网格,然后openfoam会在z轴增加一个面。

  • 前辈们俺有一个关于ICEM得问题

    Meshy
    3
    0 赞同
    3 帖子
    5k 浏览
    liujunCFDL

    @bestucan :mihu: 确实是直接删除block就行,俺发现好多操作都跟人家一样。得到得结果却不一样,可能因为我用的不是正版

  • 0 赞同
    3 帖子
    5k 浏览

    @李东岳明白啦,谢谢东岳老师。