CFD中文网

@李东岳 目前刚接触这个领域，还在学习哈哈，我们老板主打接地气，宗旨是如果不能解决实际问题，不做也罢。还是希望用CFD和人工智能的结合做一些有意义的技术。

@李东岳 老师，我已经博后啦，博士期间弄得CFD，最近人工智能不是很火嘛，想看看能不能学科交叉一下，主要想解决一些实际问题，比如，使用神经网络替换掉CFD中的计算稠密部分，或者使用PINN解决一些已有模型的泛化性问题。当然高水平论文也是需要的，所以想了解一下这个东西的前景

@李东岳 东岳老师，发SCI的话也需要说明创新点和拟解决的问题，那么这个数据驱动CFD到底解决了什么问题，创新点就是别人没弄过我弄出来了吗？有点疑惑

最近想搞一下数据驱动CFD，就先仿照东岳老师的方法，把pisoFoam中的LES模型用神经网络替换掉了，然后算了几个案例对比了一下，结果是还可以的，但是发现神经网络的计算速度要远远慢于原始方法（LES需要10s，神经网络需要80s）。那么我有两个问题想请教一下：

1. 计算速度的差异正常的吗？还是说我植入的有问题？

2. 如果神经网络的计算就是会慢很多，那么植入它的意义是什么？

@李东岳 上面的blcokMesh确实可以，但是求解器不行，后面发现是系统没装gmp插件，有些内容编译不了，安装之后就可以了

@Amadeus 找到问题了，系统没装gmp插件，有些内容编译不了，安装之后就可以了

gcc版本也换过，结果都是一样

zf@zf-Lenovo-XiaoXinAir-14-ACN-2021:~/OpenFOAM/OpenFOAM-v2406$ blockMesh 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2406                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _630d60de3b-20240620 OPENFOAM=2406 version=v2406
Arch   : "LSB;label=32;scalar=64"
Exec   : blockMesh
Date   : Dec 25 2024
Time   : 18:21:00
Host   : zf-Lenovo-XiaoXinAir-14-ACN-2021
PID    : 372228
I/O    : uncollated
Case   : /home/zf/OpenFOAM/OpenFOAM-v2406
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time



--> FOAM FATAL ERROR: (openfoam-2406)
cannot find file "/home/zf/OpenFOAM/OpenFOAM-v2406/system/controlDict"

    From virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::readStream(Foam::regIOobject&, const Foam::fileName&, const Foam::word&, bool) const
    in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 629.

FOAM exiting

zf@zf-Lenovo-XiaoXinAir-14-ACN-2021:~/OpenFOAM/OpenFOAM-v2406$ icoFoam
找不到命令 “icoFoam”，但可以通过以下软件包安装它：
sudo apt install openfoam

最近在一台Ubuntu 22.04机子上想装一个openfoam-v2406，编译完之后也没报错，且blockMesh可以正常调用，但是求解器一个都不能用，看了一下编译后的platforms文件夹，bin里面一个求解器都没有，这是什么情况，重新编译了两次了，都是这个问题。以前在很多台电脑都装过，但是第一次遇到这种问题，望大神指教！

@派大星huangxiaojie

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/

forces1-gas
{
    type          forces;

    libs          ("libforces.so");

    writeControl  timeStep;
    timeInterval  1;

    log           yes;
    patches     (impeller_66 impeller_67);
    rho           rhoInf;     // Indicates incompressible
    log           true;
    rhoInf        700;          // Redundant for incompressible
    phase    gas;
     p p;
    CofR          (0 0.35 0);    // Rotation around centre line of propeller
    pitchAxis     (0 5 0);
}

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

把每个相的扭矩值都分别处理，然后相加即可，跟实验值较为吻合

@李东岳 好的，谢谢老师

@李东岳 另外老师还想请教个问题，scalarTransportFoam中库朗数还需要小于1吗，还是说时间步长可以大一些让库朗数大一些

@李东岳 我的意思是把multiphaseEulerFoam算出来的U.particles和phi.particles一起复制过来作为scalarTransportFoam的初值，然后把原先phi.particles中的边界条件修改一下就可以了

@李东岳 老师，问题解决了，要把固体的通量场也作为初值给scalarTransportFoam，并且要把phi里面得fixedValue类型改为calculated就可以了，谢谢老师

@李东岳 或者说是不是多相流中的单相速度场并不能直接拿到scalarTransportFoam中使用，因为单相速度场没有相含率信息，例如一个网格的固含率为0，但是其还是会有一个固相速度，然而scalarTransportFoam仅读取速度场无法知道这一点，会不会导致质量不守恒的问题。

@李东岳 老师，我继续计算了4个小时也还是没有变化，出口一直没有浓度，而且T的迭代次数已经都是0了，就是说已经稳态了，可是这样如何得到停留时间分布？

@李东岳 老师，因为我想固体颗粒在液体中应该没有扩散，所以我只考虑了对流，扩散系数设了0，是这个问题吗

@李东岳 进口固定值1，其他wall和出口都是零梯度

并且算着算着迭代次数逐渐变为0了

首先使用multiphaseEulerFoam模拟了一个气液固三相连续搅拌釜，想获得固体相在釜中的停留时间分布，所以想将固体的速度场和通量场给到scalarTransportFoam中，从入口进入示踪剂以获得固体停留时间分布。这时遇到一个问题，整个釜的空时大概4个小时，但是计算大概1000s之后整个釜内的示踪剂浓度分布就不变化了，出口的浓度也达到稳定并且极低，请问这是什么问题？附上我的设置

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,T)     Gauss linearUpwind grad(T);
}

laplacianSchemes
{
    default         none;
    laplacian(DT,T) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

solvers
{
    T
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 1;
}
relaxationFactors
{

    equations
    {
        ".*"         1;
    }
}