CFD中文网

A Swin-transformer-based model for efficient compression of turbulent flow data
Image and video compression of fluid flow data

我其实想研究研究这个水几篇论文

PINN比较有意思

上图PINN，下图inr视频压缩。是不是很像？都是输入坐标和时间，用神经网络来学习一种隐式表示。感兴趣的老哥带带我

v我50， 我来跑

看了下PINN，其实和现在视频编码领域流行的隐式表征方法是一样的。不过隐式表征能直接拿偏微分方程做loss，学会流场，还是感觉挺震撼的。看起来PINN做反问题、实时仿真和流体控制是一把好手

最近流体+AI很热门啊，天气预报的工作频繁登上Nature

之前在这里学习过一段时间OpenFOAM，现在专门做AI了。

我主页：
https://air.tsinghua.edu.cn/info/1046/1555.htm
感兴趣的欢迎交流合作。

我觉得这个方向主要是CFD部分比较困难，深度学习和PyTorch很容易学。搞AI的人不愿意花时间学CFD，还是比较适合我们做CFD出身的

带我个nature子刊

mark

@hurricane007 试试，说不定能发个science

http://wiki.jikexueyuan.com/project/tensorflow-zh/tutorials/pdes.html

流不流弊

@程迪 厉害。确实有一些计算机方面的方法可以搞，不过用vtk解决了，多谢！

@李东岳 生成的随机坐标与前一个点的值有关

准备试试vtk，好像没有信息损失。。

发现paraview导出csv然后处理比较简单，虽然这样没有网格信息了，只能基于网格点上的值对任意点的T进行插值

目前做一个后处理算法，需要按照随机的位置对计算结果进行采样。后处理中需要用到复杂的算法，目前是基于python的类库实现。
最简单的想法是在python中调用probe，但是这样每采样一次就需要重新从硬盘读取整个网格的结果，太慢。
另一个想法是写一个of的后处理程序，首先读入整个网格的结果进内存，然后调用python代码，有个问题是此时python怎么访问读进来的volScalarField？即：

int main(int argc, char *argv[])
{
    #include "setRootCase.H"

    #include "createTime.H"
    #include "createMesh.H"

    Info<< "Reading field T\n" << endl;
    volScalarField T
    (
        IOobject
        (
            "T",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );

    调用一段python，里面有一个函数的功能是：对于python随机生成的坐标(x,y,z)，能够返回T[mesh.findCell(point(x,y,z))]
}

如果能在python中实现读入mesh和field，以及mesh.findCell功能，问题就简单了。

各位有何想法？多谢！

@piteqiu

    Info<< "Reading field S\n" << endl;
    volScalarField S  //1.48e-5 kg/s
    (
        IOobject
        (
            "S",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedScalar("S",dimensionSet(1 ,-3, -1 ,0 ,0 ,0 ,0),0.0)
    );
    forAll(S,i)
    {
        S[i]=0;
    }
    S[mesh.findCell(point(0,0,0))]=Q.value()/mesh.V()[mesh.findCell(point(0,0,0))];

@piteqiu 把这些物性参数都弄成场，然后随便咋赋值。。

擦擦，我是先算多孔内的达西流动，然后将结果作为边界条件传给用rans的湍流区。
参考文章：土壤大气耦合的燃气泄漏扩散数值模拟

完全实现层流区和湍流区耦合挺麻烦的（除了直接以附加源项的方式模拟多孔区，这种比较简单，fluent是这样搞的），可能也有收敛性的问题。不过我在一些其他代码中见过，如dumux。楼主要是在of里实现了这种耦合边界请借我一用！

弄好了。。

@金石为开 求个外流场文件配置。我看cfmesh的例子几何都是凹在流场里面，现在我是有几何的CAD，用了surfaceGenerateBoundingBox构建流场后不会生成网格了＝＝

@cfd-china 厉害了我的哥