Skip to content
  • 0 赞同
    1 帖子
    2k 浏览
    李东岳

    2016年autodesk公司把CNN用于不同外形的流场预测。但是CNN这面最爆炸性的成果当属2021年google那面发表在PNAS的工作。目前被引700多次。

    在Kochkov等的工作中,作者们在4倍糙化的网格上,实现了传统CFD做直接模拟的精度。这一篇工作并不是首创。google在2021年的另外一篇文章中,将他们的思想最开始应用于一个简单的对流传输方程中。

    2021年稍早,哈佛大学与google公司合作,在Zhuang等的工作中提出了Learned Intepolation。即通过大量的样本来训练出来的CFD格式。

    这对于咱们做CFD的很好理解,传统的面插值格式,可以直接从网格体心插值过来,也可以从2个网格体心或者多个网格体心取平均过来。对于CNN的插值格式,是通过针对样本进行多次采样,实现一个类似可以理解的在糙网格下,实现细网格上同样结果的插值格式。更好的理解可以这样:

    在细网格下,传统CFD预测的面插值是1,糙网格下,CFD预测的是0.5。CNN经过多次特征与标签的采样之后,在糙网格下预测的插值可能是0.9。通过这样来实现糙网格的高精度模拟。

    在google公司的工作中,除此之外,全部用的是CFD传统的解法。这对于使用深度学习从头到尾的计算,有很多优势。最大的优势就是,可以保证传统CFD的守恒特性。

    在上图中,第1行是CFD计算的细网格的结果。第4行是CFD糙网格计算的结果。第2行,是CFD细网格的结果,通过平均技术投影到糙网格的结果。第3行,是CNN计算的结果。

    可以看出,第3行与第1行结果相当!也就是CNN的“学习来的数值格式”可以在糙网格上进行直接模拟。

    问题还是有的,这部分数据样本还是来自于CFD。如果CFD本身样本不对。那学习的样本本身就是错误的。

    需要保持警惕的是:从错误中来学习,必然是大错特错。

    google的工作是通过他们自家的机器学习代码jax来实现的。google的工程师为了弄这个,在jax上面写了个错位网格的CFD代码。。

    CFD自家亲戚都懂的,现在很少有错位网格的CFD代码了。

    因为google那面CFD层面比较初级,所以需要处理成同位网格,时间步进还没有植入高阶格式,多重网格求解器也还没有,边界条件只有周期边界条件。。。

    目前,google团队还在公开欢迎各类CFD大佬去参与合作。岳子就是没有闲工夫,要不然岳子也要躺一躺这个浑水。

  • CFD周边笔记本

    CFD彩虹条
    5
    0 赞同
    5 帖子
    5k 浏览

    @李东岳 好嘞,感谢信任,我是很喜欢咱们的笔记本,希望做得更好:146:

  • 0 赞同
    1 帖子
    2k 浏览

    工科末流985硕士毕业,想在武汉找CFD相关的工作,家人们有没有推荐的公司。

  • numeca tabgen插值方法

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

    请问有大佬知道numeca的tabgen生成物性时的bicubic插值方法使用的是什么公式吗?它里面的插值系数又是怎么排列的呢?

  • LES和RANS中,离散的对象是相同的吗?

    Algorithm
    31
    0 赞同
    31 帖子
    54k 浏览
    H

    如果使用standard KE,计算二维流场,在某个高雷诺数下把Cd, Cl, St都算对,那么说明碰对了,意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

  • 关于piso算法中的压力梯度项

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

    @卡洛 无痛苦ns方程里面有,专门讨论的flux()函数,你看下

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

    @李东岳 李老师,在写fvm::laplacian项之前已经通过体积分数alpha计算过nu了。因为我觉得直接修改turbulence::类里面的nu很麻烦。所以我就自己定义了提及标量场nu然后正常使用alpha计算之后带入laplacian,所以获得了与文献一致的结果。

  • 理想气体的壁面条件

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

    感谢@李东岳 @xpqiu 我查了下wallShearStress会默认输出所有壁面的剪切应力,是一个volVectorField。
    另外我看rhoCentralFoam里关于壁面边界有用noSlip也有用slip的。

  • 修改dpmfoam报错

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

    这是不是代表我的源代码库里面修改的有问题,尽管编译成功了,但是有代码错误无法运行?

    有可能。经常如此

  • openfoam中的重叠网格浮体位移的提取

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

    您好,请问一下,您知道有什么好的方法可以提取重叠网格的位移吗?能否分享一个算例呢?@孤岛

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

    恭喜恭喜

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

    我看了下这个文章。这个文章非常适合练手,挺有意思。通过植入3个传输方程,处理传热流动。openfoam植入非常straightforward。做个标记。或许下次可以当做课程资料。Fluent那面不清楚了。

    你为什么不直接用传热流而是一定要自己植入一个呢?

  • 0 赞同
    21 帖子
    31k 浏览
    H

    但是RC插值没有离散,因为它是显式计算源项 -- “速度插值+高阶压力导数项”。

  • fluent 角速度设定

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

    请教大家 在fluent里如何给气泡设定初始旋转角速度?

  • Fluent那面可以模拟2D的粒子流么

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

    @袁宝强 啊我知道了 谢谢 :xiexie:

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

    你所谓的高阶格式,指的哪种高阶格式?

  • LES一个方程

    Algorithm
    5
    0 赞同
    5 帖子
    6k 浏览
    李东岳

    此问题同类于湍流动能耗散率的定义。湍流动能耗散率的定义为:
    \begin{equation}
    \varepsilon=2\nu_t\overline{S_{ij}S_{ij}}
    \end{equation}
    经过各向同性假定后有:
    \begin{equation}
    \varepsilon=\nu_t\overline{\frac{\p u_i'}{\p x_k}\frac{\p u_i'}{\p x_k}}
    \end{equation}

  • custom rate exponent如何确定?

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    梦中飞翔的阿涛

    6b6b5aa0-11d9-4944-accc-20bb6eb70b36-image.png 表中的custom rate exponent是怎么确定的?感谢各位。73ea2e62-c36f-4287-9b1d-b62745bea525-image.png

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

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

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

  • ICEM里“块”和“模型”是什么关系?

    Meshy
    2
    0 赞同
    2 帖子
    5k 浏览
    Elibathe康E

    块就是将你的整个计算区域分割成几个区域,每个区域有各自的网格分布规律;
    模型应该是你计算的几何外形吧