OpenFOAM libtorch tutorial step by step

李东岳

@chengan-wang 太老了，gcc版本啥的都太低了，我没有折腾过，你得亲自自己都尝试一下看看哪个能用

chengan.wang

@李东岳李老师，http://dyfluid.com/pinn.html ，#include "output.H"这个文件能补充上传么，想研究一下您的完整代码，谢谢

李东岳

@chengan-wang 那个我已经删了，你可以参考这个写个类似的，纯粹是C++的输出数据：

        if (iter % 500 == 0) 
        {
            cout<< iter << " loss = " << loss.item<double>() << endl;

            auto meshxy = torch::cat({mesh,x,y}, 1);
            std::ofstream filex("results");
            filex << meshxy << "\n";
            filex.close();
        }

http://dyfluid.com/pinn2.html

李东岳

我给你手写了一个output.H：

        if (iter % 500 == 0) 
        {
            cout<< iter << " loss = " << loss.item<double>() << endl;
 
            std::ofstream filex("results");
            filex << upred << "\n";
            filex.close();
        }

chengan.wang

@李东岳谢谢李老师

chengan.wang

@李东岳我看了您提供的PINN算例，感觉目前还没有真正跟of结合起来，是吧？只不过用了of的c++编译环境，可以这么理解么？

李东岳

是。PINN不需要跟FVM结合，只不过用了of的c++编译环境。想跟OpenFOAM结合，除了画网格能结合一下，其他的也结合不上。

李东岳

@chengan-wang

You can use the following code to initialize libtorch tensors from OpenFOAM mesh. It is much easier to generate any kinds of geometries. Please see the following clipped cavity geometry. I am using English since I hope it can be helpful for all the people.

微信图片_20241017082858.png

    // openfoam initialization
    vectorField ofMesh = mesh.C();
    auto meshTorch = torch::full({ofMesh.size(),1,2}, 0.0);
    forAll(ofMesh, i)
    {
        double x = ofMesh[i].x();
        double y = ofMesh[i].y();
        auto posi = torch::tensor({x,y}).reshape({1,2});
        meshTorch[i] = posi;
    }

    const surfaceVectorField& Cf = mesh.Cf(); 
    label topID = mesh.boundaryMesh().findPatchID("top");
    label leftID = mesh.boundaryMesh().findPatchID("left");
    label rightID = mesh.boundaryMesh().findPatchID("right");
    label bottomID = mesh.boundaryMesh().findPatchID("bottom");
    vectorField cfTop = Cf.boundaryField()[topID];
    vectorField cfLeft = Cf.boundaryField()[leftID];
    vectorField cfRight = Cf.boundaryField()[rightID];
    vectorField cfBottom = Cf.boundaryField()[bottomID];

    auto top = torch::full({cfTop.size(),1,2}, 0.0);
    auto left = torch::full({cfLeft.size(),1,2}, 0.0);
    auto right = torch::full({cfRight.size(),1,2}, 0.0);
    auto bottom = torch::full({cfBottom.size(),1,2}, 0.0);

    forAll(cfTop, i)
    {
        double x = cfTop[i].x();
        double y = cfTop[i].y();
        auto posi = torch::tensor({x,y}).reshape({1,2});
        top[i] = posi;
    }

    forAll(cfRight, i)
    {
        double x = cfRight[i].x();
        double y = cfRight[i].y();
        auto posi = torch::tensor({x,y}).reshape({1,2});
        right[i] = posi;
    }

    forAll(cfLeft, i)
    {
        double x = cfLeft[i].x();
        double y = cfLeft[i].y();
        auto posi = torch::tensor({x,y}).reshape({1,2});
        left[i] = posi;
    }

    forAll(cfBottom, i)
    {
        double x = cfBottom[i].x();
        double y = cfBottom[i].y();
        auto posi = torch::tensor({x,y}).reshape({1,2});
        bottom[i] = posi;
    }

    meshTorch = meshTorch.reshape({-1,2});
    top = top.reshape({-1,2});
    left = left.reshape({-1,2});
    right = right.reshape({-1,2});
    bottom = bottom.reshape({-1,2});

    auto u_top = torch::full({cfTop.size()}, vel);
    auto v_top = torch::full({cfTop.size()}, 0.0);

    std::cout<< "top: " << top << std::endl;
    std::cout<< "left: " << left << std::endl;
    std::cout<< "right: " << right << std::endl;
    std::cout<< "bottom: " << bottom << std::endl;

水润晨阳

@李东岳在 OpenFOAM libtorch tutorial step by step 中说：

net.H

#include <torch/torch.h>

class NN
:
    public torch::nn::Module 
{
    torch::nn::Sequential net_;

public:

    NN();

    torch::Tensor forward(torch::Tensor x);
};

net.C

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Copyright (C) 2022-2023 OpenFOAM Foundation
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.

Class

See also

SourceFiles

\*---------------------------------------------------------------------------*/
#include "net.H"

NN::NN()
{
    net_ = register_module
    (
        "net", 
        torch::nn::Sequential
        (
            torch::nn::Linear(2,20),
            torch::nn::ReLU(),
            torch::nn::Linear(20,30),
            torch::nn::ReLU(),
            torch::nn::Linear(30,30),
            torch::nn::ReLU(),
            torch::nn::Linear(30,20),
            torch::nn::ReLU(),
            torch::nn::Linear(20,20),
            torch::nn::ReLU(),
            torch::nn::Linear(20,1)
        )
    );
};

torch::Tensor NN::forward(torch::Tensor x) 
{
    return net_->forward(x);
}

这个要编译成库，挂到PINNFoam上面

李老师请问这个如何编译并挂载呢可以详细的指导一下吗

DJohn1211

@leaon 您好，我也和您遇到了一样的问题，请问您如何解决的呢

DJohn1211

@李东岳李老师您好，我在看您发的那个顶盖驱动流的代码的时候有个问题，我的电脑安装了openfoam10和11，同样的文件内容为什么在of10下我可以成功编译，但在of11下不可以呢。在11下，一开始提示找不到libtriSurface.so和libsurfMesh.so，然后我在bashrc中添加了，export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/home/john/OpenFOAM/OpenFOAM-11/platforms/linux64GccDPInt32Opt/lib，随后再次运行出现了这样的问题。

李东岳

of10跟of11的Make下面的options不一样，我先用of10吧。等过段时间再debug

CFD中文网

OpenFOAM libtorch tutorial step by step