CFD中文网

@youhaoyu 谢谢大佬你的解答，我再看看学习学习，对于您的工作气泡长大是否已经完成？很希望和你进行交流2055977489是我的qq号

@youhaoyu 你好，能学习一下您的植入方法吗，我最近也被困在这了，期待你的回复，感谢

感谢！感谢wangfei9088大佬，您说的对，果然是量纲不同，也是我昏了头脑。在0文件夹下检查许久没发现错误，我在fmm文件的量纲给错了。再次感谢您的解答！

报错如下：

有路过的大佬了解吗？感谢！

@wangfei9088 感谢大佬解惑！一开始我看在OpenFOAM-2206版本下的liquidproperties里面的物质里面不包含下面两行代码，我就给删掉了。是我画蛇添足了。再次感谢！
#include "thermodynamicConstants.H"
using namespace Foam::constant::thermodynamic;
刚刚试了您说的两种方法，都可以正常编译，使用也没有问题。感谢！祝收敛！

目前的问题如下：

报错如下：

好像是Tstd这个变量他无法识别，这个原版本是参考：https://cfd-china.com/post/35393
这个也是我原来的帖子，感谢当时的wangfei9088大佬解答，现在向转到OpenFOAM-v2206版本上又卡了，路过的大佬帮忙看看，解答一二。

@ann 你好！请问你最后的穿透长度是怎么解决的。

@李东岳 感谢李老师的回复！谢谢您！我再多读读相关信息，谢谢！

@李东岳 谢谢李老师的解答，我也在cfd-online找到相关帖子里面有这样一句话：The dynamic Smagorinsky in previous versions was averaging the coefficient over the whole domain (which defeats the whole purpose of a dynamic solution). However you can use the dynamicLagrangian model which is supposed to perform the averaging over a streamline（在以前的版本中，动态Smagorinsky是对整个域的系数求平均值(这违背了动态解决方案的全部目的)。然而，你可以使用动态拉格朗日模型，它应该在流线上执行平均）。所以dynamicSmagorinsky和dynamicLagrangian仅仅是平均的方式不同吗？它用在喷雾雾化过程可行吗？（可能这个需要看详细代码，原谅学生的这个问题不太成熟，也是希望有做过喷雾的大佬看到这个帖子，了解的帮忙解答一二！）感谢李老师。

我想在OpenFOAM-10版本上用dynamicSmagorinsky湍流模型，在网上找到一个版本可惜它老是显示无法加载动态库参考我之前的帖子：链接文本
又在网上找到一个新的dynamicSmagorinsky模型：https://github.com/AlbertoPa/dynamicSmagorinsky
想把它改成OpenFOAM-10版本下使用，无奈我才疏学浅，不知怎么下手，只知道要改写make文件下option文件（也不知其他文件需不需要修改）：

EXE_INC = \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/incompressible/lnInclude \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \

LIB_LIBS = \
    -lturbulenceModels \
    -lincompressibleTurbulenceModels \
    -lincompressibleTransportModels \
    -lcompressibleTurbulenceModels \
    -lcompressibleTransportModels \
    -lfiniteVolume \
    -lmeshTools \
    -lOpenFOAM

路过的大佬帮忙给点意见，感谢！感谢！感谢！

@李东岳 感谢李老师的回复：
编译的信息如下：

wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file makeIncompressibleMomentumTransportModel.C
g++ -std=c++14 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3  -DNoRepository -ftemplate-depth-100 -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/physicalProperties/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/MomentumTransportModels/momentumTransportModels/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/MomentumTransportModels/incompressible/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/finiteVolume/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/meshTools/lnInclude -IlnInclude -I. -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/OpenFOAM/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/OSspecific/POSIX/lnInclude   -fPIC -c makeIncompressibleMomentumTransportModel.C -o Make/linux64GccDPInt32Opt/makeIncompressibleMomentumTransportModel.o
g++ -std=c++14 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3  -DNoRepository -ftemplate-depth-100 -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/physicalProperties/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/MomentumTransportModels/momentumTransportModels/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/MomentumTransportModels/incompressible/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/finiteVolume/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/meshTools/lnInclude -IlnInclude -I. -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/OpenFOAM/lnInclude -I/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/src/OSspecific/POSIX/lnInclude   -fPIC -fuse-ld=bfd -shared -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPInt32Opt/makeIncompressibleMomentumTransportModel.o -L/home/postgraduate/OpenFOAM/OpenFOAM-10/OpenFOAM-10/platforms/linux64GccDPInt32Opt/lib \
    -lphysicalProperties -lmomentumTransportModels -lfiniteVolume -lmeshTools -lstdc++ -fopenmp   -o /home/postgraduate/OpenFOAM/postgraduate-10/platforms/linux64GccDPInt32Opt/lib/DynSmag.so
我看在/home/postgraduate/OpenFOAM/postgraduate-10/platforms/linux64GccDPInt32Opt/lib/DynSmag.so有相应的.so文件应该是编译成功了。

@风大仙 你好，我也遇到同样的问题，请问这个问题后来你是怎么解决的？感谢！

@李东岳 这个是我在网上找的，因为我用的是OpenFOAM-10，我在网上找到了在这个版本下的dynamicSmagorinsky模型，下面是DnySmag的GitHub网址：https://github.com/halimarefat/DynSmag
李老师帮忙看看！

在reactingFoam求解器下的案例中，输入blockMesh运行时就会发生下面报错：
Create time

--> FOAM Warning :
From function void* Foam::dlOpen(const Foam::fileName&, bool)
in file POSIX.C at line 1247
dlopen error : /home/postgraduate/OpenFOAM/postgraduate-10/platforms/linux64GccDPInt32Opt/lib/DynSmag.so: undefined symbol: _ZTIN4Foam36incompressibleMomentumTransportModelE
--> FOAM Warning :
From function bool Foam::dlLibraryTable::open(const Foam::fileName&, bool)
in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 106
could not load "DynSmag.so"
因为能加载这个动态库，导致后面湍流文件不能识别DynSmag模型，计算也不能继续。
用的求解器是reactingFoam，版本号是OpenFOAM-10。求解答！感谢。

@hy1112006 请问你这里：在文件SOWFA-2.4.0/src/turbulenceModels/incompressible/LES/Make/options中加入这句话即可：
-lincompressibleLESModels
添加的-lincompressibleLESModels这个是根据你什么来的？
不知道你还记得吗？打扰了

@李东岳 谢谢李老师的回答，我再找找原因！

我的计算域是50×50×75mm，使用的是OpenFOAM-10版本下的reactingFoam求解器。之前是使用的是OpenFOAM-8版本下的sprayFoam求解器，好像在OpenFOAM-10版本下，sprayFoam求解器合并到reactingFoam里面了，所以目前使用reactingFoam求解器。一开始使用blockMesh {hex (0 1 2 3 4 5 6 7) (100 100 150) simpleGrading (1 1 1)}，这样最小网格是0.5mm，此时还可以正常计算。因为网格优点粗糙，我看有的文献网格在0.15mm以下。接着我尝试细化网格，下面是我的网格情况：


可是这时候刚开始计算就会报以下错误：

中间区域的网格尺寸应该是0.2mm。没有细化网格之前的的时间步是1e-06，细化网格后我改成1e-07.这是基于单核计算的结果。喷雾速度大小为250m/s。
湍流模型的设置如下：

LES
{
    model          kEqn;

    delta          cubeRootVol;

    turbulence      on;

    printCoeffs     on;


    kEqnCoeffs
    {
        ck      0.07;
    }

    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }

    PrandtlCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Cdelta           0.158;
    }

    vanDriestCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Aplus            26;
        Cdelta           0.158;
    }

    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        maxDeltaRatio   1.1;
    }
}

麻烦有了解LES方面的大佬帮忙看看，感谢！感谢！感谢！（重要的事情说三遍！）

@coolhhh 对，我使用这个命令：alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'时，还没有把anaconda的环境给注释掉，OpenFOAM也可以正常使用，但是并行命令不能使用。
我大概懂大佬你的意思了，如果使用了alias命令，只要没有输入Anaconda3应该就不会激活anaconda环境。可能是因为我这边前人在给anaconda设置环境的时候，有点繁琐（具体我也没弄清楚），所以我可能只使用一条alias命令，不能完全把相应的anaconda环境给包含进去。所以才导致使用alias Anaconda3='PATH=/home/install/Anaconda3/bin:$PATH'命令还不能消除使用并行命令带来的报错。再次感谢您的回答，谢谢！

@李东岳 李东岳老师，我在系统.bashrc文件中把前人使用anaconda的环境给注释掉了，后续就可以正常使用OpenFOAM的并行命令了。再次感谢您的回答。谢谢！

@coolhhh 你好，大佬，你的这个方法我试过了，可以解决打开终端时出现的 x86_64-conda_cos6-linux-gnu-cc: error: unrecognized command line option '--showme:link'报错。切换OpenFOAM版本环境也没有类似报错。但是如果运行并行计算命令，程序还是会自动退出，应该是系统bashrc文件中存在anaconda环境，二者不兼容的原因。再次感谢您的解答，谢谢！