CFD中文网

@李东岳 李老师您好，很感谢您开源的代码。我的课题就是在您开源的基础上一步步做的，现在遇到了一个问题，我按照和您相同的思路，对某个算例模拟、训练、预测……，所不同的是我现在采用的是二维卷积。并且我在网络上参照图像方面的知识，添加空间注意力、残差连接等等，虽然在流场预测的趋势上有了较好的拟合结果，但是所预测的始终是这种有规律的间隙，我有点困惑，想和您请教一下这可能是什么原因造成的。
这是我去预测体积分数场的一个结果，

宽窄带呈现着一个1宽7窄的规律，下面是一个网格的内容

我在这附上一个我简单的代码,这个是我最初测试时用到的，也会有明显的条状分布

                torch::nn::Conv2d(torch::nn::Conv2dOptions(2, 5, 3).padding(1)),
                torch::nn::BatchNorm2d(5),
                torch::nn::Tanh(),
                
                torch::nn::Conv2d(torch::nn::Conv2dOptions(5, 10, 3).padding(1)),
                torch::nn::BatchNorm2d(10),
                torch::nn::Tanh(),
                
                torch::nn::Conv2d(torch::nn::Conv2dOptions(10, 10, 3).padding(1)),
                torch::nn::BatchNorm2d(10),
                torch::nn::Tanh(),
                
                torch::nn::Conv2d(torch::nn::Conv2dOptions(10, 8, 3).padding(1)),
                torch::nn::BatchNorm2d(8),
                torch::nn::Tanh(),

                torch::nn::Conv2d(torch::nn::Conv2dOptions(8, 1, 1)) 

TORCH_MODULE(CNN);

自己给自己补充一下，我现在采取了同样的网格，所不同的是，当我将前后面在blockmesh，u、alpha、p中设置为empty时，便会出现液面不平整的现象，右图为设置为empty，左图为wall

各位老师大佬好，我在做一个气泡上浮的算例，想模拟二维情况下水缸中一个气泡的上浮过程，现在遇到了一点问题，我通过setifelds设置了一个球形气泡和水气的分界线，为什么在第一个时间步下液面不平整呢?但是我同样的设置采用动网格计算的液面则是平整的，如下图
这是我setFields中的设置

defaultFieldValues
(
    volScalarFieldValue alpha.water 1
);

regions
(
    sphereToCell
    {
        centre (0 0 0.005);
        radius 0.0025;
        fieldValues
        (
            volScalarFieldValue alpha.water 0
        );
    }
    
    boxToCell
    {
        box (-0.025 -0.001 0.035) (0.025 0.001 0.04);
        radius 0.01;
        fieldValues
        (
            volScalarFieldValue alpha.water 0
        );
    }
);```

您好，请问这个有没有官方算例呢

@李东岳 李老师您好，我在复现您这个例子的时候有一些细节有些糊涂？想请问您具体是怎么做的呢，能分享一下您的算例吗

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

代码越来越复杂。这是执行一系列的卷积、池化操作等输出的8个通道的流场结果：

http://dyfluid.com/cnn.html#id6

下图是自编码器进行重组的流场：

@李东岳 李老师您好，我在看您发的那个顶盖驱动流的代码的时候有个问题，我的电脑安装了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，随后再次运行出现了这样的问题。

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

@李东岳 在 openfoam11中incompressiblevof求解器 中说：

volScalarField nuML("nuML", mesh, dimensionedScalar("nu", dimensionSet(0, 2, -1, 0, 0, 0, 0), 0.0));

你有一行代码这么写的，这么写的不对

李老师，我在定义这个变量的时候是这么定义的：

volScalarField nuML
(
    IOobject
    (
        "nuML",                
        runTime.timeName(),    
        mesh,                  
        IOobject::NO_READ,   
        IOobject::NO_WRITE   
    ),
    mesh,
    dimensionedScalar("zero", dimViscosity, 0.0) 
);

但我还是没有运行成功您的教程，想请问李老师可以把您的算例在论坛共享一下吗？然后我还有一个问题是这样的，我想在of10下icofoam的基础上先修改，但是我只是在头文件引入了#include "NN.H"，并没有在求解器中修改内容，就不能编译，我看了一下报错，大致的意思是openfoam和libtorch的数据结构不兼容，但是按照您在http://dyfluid.com/上的算例来看不应该会出现这种报错，想请教李老师一下这其中的原因。这是我这部分的log文件log.wmake

@李东岳 在 openfoam11中incompressiblevof求解器 中说：

下周给你做个算例更新到 dyfluid.com

李老师您好，感谢您百忙之中更新的数据驱动LES，但我在运行的时候遇到了一些问题，我按照您在网站上的步骤这么设计了一下。

#include <torch/torch.h>
#include "fvCFD.H"
#include "viscosityModel.H"
#include "incompressibleMomentumTransportModels.H"
#include "pisoControl.H"
#include "pressureReference.H"
#include "fvModels.H"
#include "fvConstraints.H"
#include "NN.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
    // 初始化OpenFOAM运行时间
    #include "setRootCaseLists.H"
    #include "createTime.H"
    #include "createMesh.H"

    pisoControl piso(mesh);

    #include "createFields.H"
    #include "initContinuityErrs.H"


    volScalarField nuML
    (
        IOobject
        (
            "nuML",                
            runTime.timeName(),    
            mesh,                  
            IOobject::NO_READ,   
            IOobject::NO_WRITE   
        ),
        mesh,
        dimensionedScalar("zero", dimViscosity, 0.0) 
    );


    // 创建并加载神经网络
    auto net = std::make_shared<NN>();
    torch::load(net, "net.pth"); // 加载已训练的模型
    

    while (runTime.loop())
    {
        Info << "Time = " << runTime.userTimeName() << nl << endl;

        #include "CourantNo.H"
        net->eval(); // 设置网络为评估模式
        
        // 计算速度梯度
        volTensorField gradU("", fvc::grad(U));

        // 准备输入数据
        auto inBatch = torch::zeros({nuML.size(), 9}, torch::kFloat);

        forAll(nuML, celli)
        {
            inBatch[celli][0] = gradU[celli].component(tensor::XX);
            inBatch[celli][1] = gradU[celli].component(tensor::XY);
            inBatch[celli][2] = gradU[celli].component(tensor::XZ);
            inBatch[celli][3] = gradU[celli].component(tensor::YX);
            inBatch[celli][4] = gradU[celli].component(tensor::YY);
            inBatch[celli][5] = gradU[celli].component(tensor::YZ);
            inBatch[celli][6] = gradU[celli].component(tensor::ZX);
            inBatch[celli][7] = gradU[celli].component(tensor::ZY);
            inBatch[celli][8] = gradU[celli].component(tensor::ZZ);
        }

        // 执行神经网络预测
        auto yPred = net->forward(inBatch);

        // 更新粘性系数 nuML
        forAll(nuML, celli)
        {
            nuML[celli] = yPred[celli].item<float>() / 1e5 + 2e-5;
        }


        runTime.write();
    }

    Info << "End" << endl;
    return 0;
}

但是出现了一些warning，不知道哪里出现了问题，我把我的log文件附在这里，想请李老师有空的时候能帮忙看一下问题，谢谢李老师。log.wmake

谢谢李老师

您好，我最近也在研究of11的求解器，觉得应该去modules文件夹下找你需要自定义的文件

李老师您好，打扰您了，我尝试了几次混编，想在pimple循环中对动量方程的求解不采用原来的方式，但总是调用libtorch失败，不知道李老师有没有时间能写一个简单的示例代码呢

李老师您好，我看过您发的OpenFOAM libtorch tutorial step by step，我现在有个想法想请教一下您，您写的pinn的代码相当于把整个求解器替换掉，那我如果只通过网络替换其中部分方程的求解该怎么做呢？向请教您一下。举个例子来说，我如果只想把其中的速度方程用一个简单的卷积网络来替换，其他部分还保持原interfoam求解器，我应该怎么自定义这个求解器呢？

想请问一下论坛里的各位大佬，我最近新安装了of11，看到interfoam求解器变成了incompressiblevof，我想在新版本下重新自定义一下一个新的类似于incompressiblevof求解器，但是我感觉和原来的interfoam求解器差别很大，比如原来interfoam中的ueqn的计算文件我在of11中不太清除incompressiblevof是怎么调用的，或者说在of11中incomprssiblevof是如何解这套方程的呢!
(
这是原来interfoam中计算u的文件内容

小白请教一下各位大佬