CFD中文网

@wangfei9088
谢谢大佬，我这就试试~

各位老师，大佬们好：
我目前正在学习如何基于OpenFOAM10自定义湍流模型，在尝试编写湍流模型的过程中，我按照旧版本的教程进行(实际上我没找到OpenFOAM10的教程), 湍流模型编译过程没有产生任何错误，但是在仿真运行时会出现

pimpleFoam: symbol lookup error: /home/cfd/OpenFOAM/cfd-10/platforms/linux64GccDPInt32Opt/lib/libbitm.so: undefined symbol: _ZN4Foam9RASModels10mykEpsilonINS_36incompressibleMomentumTransportModelEEC1ERKNS_17geometricOneFieldES6_RKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS7_IdNS_13fvsPatchFieldENS_11surfaceMes[files](/assets/uploads/files/1732866400685) [log.pimpleFoam](/assets/uploads/files/1732866400697-log.pimplefoam) [log.wmake](/assets/uploads/files/1732866400698-log.wmake) [myincompressibleMomentumTransportModels.C](/assets/uploads/files/1732866400701-myincompressiblemomentumtransportmodels.c) [mykEpsilon.C](/assets/uploads/files/1732866400703-mykepsilon.c) [mykEpsilon.H](/assets/uploads/files/1732866400704-mykepsilon.h) [options](/assets/uploads/files/1732866400706) hEEESJ_RKNS_9viscosityERKNS_4wordE

这一的错误。
即使我未对kepsilon做任何改动，仅仅是按照教程把所有的kepsilon替换为mykepsilon。
同时，我也通过myincompressibleMomentumTransportModels.C对湍流模型进行了实例化。
最后的结果就是OpenFOAM可以识别湍流模型，但是会产生这个报错。
请问能否有大佬指点我一下如何在openfoam10中正确的植入一个湍流模型？非常感谢！
以下是我本次使用的文件
log.pimpleFoam
myincompressibleMomentumTransportModels.C
mykEpsilon.C
mykEpsilon.H
files
options

自回一下：
问题已解决，经过多次尝试需要在film的physical properteis 中增加温度限制，粒子数据之类的信息，还需要查询具体的粘度系数和热熔系数公式，我把文件内容放在下面，希望对后来人有用

thermoType
    {
        type        heRhoThermo;
        mixture     pureMixture;
        transport   polynomial;
        thermo      janaf;
        equationOfState perfectGas;
        specie      specie;
        energy      sensibleInternalEnergy;
    }

mixture
{
    specie
    {
        molWeight       100.20557;
    }
    thermodynamics
    {
        Tlow        200;
        Thigh       3500;
        Tcommon     1000;
        highCpCoeffs    ( 20.4565203 0.0348575357 -1.09226846e-05 1.67201776e-09 -9.8102485e-14 -32555.6365 -80.4405017 );
        lowCpCoeffs     ( 11.1532994 -0.00949419773 0.000195572075 -2.49753662e-07 9.84877715e-11 -26768.8904 -15.9096837 );
    }
    transport
    {
        muCoeffs<8>     ( -10.378 0.244 -0.000054  0 0 0 0 0 ); //( 1000 -0.05 0.003 0 0 0 0 0 );
        kappaCoeffs<8>    ( 80.10 3.45542 0.00128822 -0.00000166835 0 0 0 0 );
    }
    equationOfState
    {
        rho         684.0;  // 庚烷的密度，单位 kg/m³
    }
    specie
    {
        Hf          0;
    }
}
sigma
{
    type    constant;
    sigma   0.07;
}

deltaWet    1e-4;

各位老师，师兄，师姐，大佬们好，
在OpenFOAM11中，我使用surfaceFilm模型进行一个3D仿真。仿真流程如下：

1. 先运行一个关闭surfaceFilm和spray injection 模型，让T_inlet=900k 运行0.05s 初始化流场。
2. 激活燃烧，喷雾和薄膜模型，向模型中注入正庚烷C7H16（LPT），然后进行燃烧仿真。
   在只使用燃烧和喷雾模型的时候，仿真可以正确的进行并获得有物理意义的结果。
   在激活surfaceFilm模型后，模型设置参照openFOAM自带算例的“heatBox”，仿真疑似因为在薄膜位置计算温度时浮点溢出而不能运行。
   在"CFDonline"中，我查到在之前的版本中 可以使用“fvOption” 来对温度的上下限进行限制以防温度失控，但是在OpenFOAM11版本中似乎没有该功能。我在中外各个网站查了一圈，有关于该模型的使用案例似乎特别少。
   请问这个问题出现的原因是什么？我应该如何解决呢
   下附报错图片，我认为比较有用的设置文件截图，以及我的仿真文件
   报错如下
   
   可能有用的文件截图：

film的温度

film的fvModel

关闭薄膜模型时的结果

求求大佬们救救孩子吧~~~~~~~
附上我的仿真文件下载链接:
https://www.jianguoyun.com/p/DTKpvfgQlPqQDBi01tUFIAA

你好，我做过一个类似的仿真，该仿真是可以很好的收敛的，你可以对比一下。
https://www.jianguoyun.com/p/DXJmh5YQxfn5CxjGx8EFIAA

呜呜呜可以注册了，感谢老师