我没用过2.3 不过一般情况下自己改模型的标准流程是：

1. copy
2. 修改文件内模型命名，用sed -i可以批量修改，比如原来叫kOmega，你现在可以改为MykOmega，注意看看是不是所有的class都改了，是不是有的地方不应该改，然后sed就给你全局 修改了？
3. 修改Make中的files和option

files中注意你的目标文件名是否已经修改为心的 .C文件
sprayFoamTest.C//原来是sprayFoam
//是否输出到你的USER_LIBBIN，
//而且有一个心的命名，如果成功编译在你的USER_APPBIN中会出现一个新的文件，以你设置的下面的名字命名，例如我的就是sprayFoamTest
EXE = $(FOAM_USER_APPBIN)/sprayFoamTest

首行添加引用的地方
LIB_USER_SRC = $(WM_PROJECT_USER_DIR)/src

EXE_INC = \
    -I. \
    -I../reactingParcelFoam \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I${LIB_SRC}/meshTools/lnInclude \
    -I${LIB_SRC}/sampling/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
    -I$(LIB_USER_SRC)/lagrangian/basic/lnInclude \
    -I$(LIB_USER_SRC)/lagrangian/intermediate/lnInclude \
    -I$(LIB_USER_SRC)/lagrangian/spray/lnInclude \
    -I$(LIB_USER_SRC)/lagrangian/distributionModels/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
    -I$(LIB_SRC)/transportModels/compressible/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/liquidMixtureProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/solidProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/solidMixtureProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/thermophysicalFunctions/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/SLGThermo/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
    -I$(LIB_SRC)/ODE/lnInclude \
    -I$(LIB_SRC)/regionModels/regionModel/lnInclude \
    -I$(LIB_SRC)/regionModels/surfaceFilmModels/lnInclude \
    -I$(LIB_SRC)/combustionModels/lnInclude

EXE_LIBS = \
    -L$(FOAM_USER_LIBBIN) \ //这个地方也别忘了添加，有可能编译的时候找不到改位置 如果不设置的话
    -lturbulenceModels \
    -lcompressibleTurbulenceModels \
    -llagrangian \
    -llagrangianIntermediate \
    -llagrangianTurbulence \
    -llagrangianSpraytest \  //这边要引用你自己修改的模型的名称
    -lspecie \
    -lcompressibleTransportModels \
    -lfluidThermophysicalModels \
    -lliquidProperties \
    -lliquidMixtureProperties \
    -lsolidProperties \
    -lsolidMixtureProperties \
    -lthermophysicalFunctions \
    -lreactionThermophysicalModels \
    -lSLGThermo \
    -lchemistryModel \
    -lradiationModels \
    -lODE \
    -lregionModels \
    -lsurfaceFilmModels \
    -lcombustionModels \
    -lfiniteVolume \
    -lfvOptions \
    -lmeshTools \
    -lsampling

希望能帮到你

有条件的可以自行油管：https://www.youtube.com/watch?v=K-nAF3qAPTc
感谢youtuber:Mark Kimber
侵删

1. 首先你要有一个自己的进口速度场的公式

   

   或者自行拟合一个，根据实际情况来。

   

2. 在Terminal中输入

   writeCellCentres
   

   获取网格每个cell的坐标，会生成ccx，ccy，ccz三个文件，内涵每个patch中的每个cell的坐标(x,y,z)，找到inlet patch，我的是topWall。3个文件中的topwall（就是我的inlet)的2100代表的是这个patch上又2100个网格。同理对于ccy文件同样找到这个patch的所有cell坐标。分别复制这些数据，生成两个新的文本，叫inletX, inletY。根据实际网格情况，以及坐标轴自行选择用哪2个文件。

   
   

3. 进入octave（一个类似matlab的程序，youtuber是这么说的）这一步主要是根据你的速度场的公式，xy坐标，计算出进口的vector。
   

   x= load('inletX');   //读取inlet patch的x坐标
   y=load('inletY');    //读取inlet patch的y坐标
   Uz==1.1+0.18988/0.01415/sqrt(3.1415926/2)*exp(-2*(x.^2+y.^2)/0.1415.^2);  //公式
   fid= fopen('inletU','w')// inletU文件要提前新建好，不知道为什么视频中输出结果为5，我的结果一直为3，不过不影响使用
   N=length(x)
   
   for i=1:N
   fprintf(fid,'(%8.4f %8.4f %8.4f) \n',0,0,Uz(i));
   end
   
   fclose(fid)
   
   

   

   

   其实第3步完全可以用excel做。。。。。。。

4. 把得到的vector的值带入到0/U中
   

5. 进入paraview中查看即可

   

@zhe 因为最后一步你重新给p.U()速度了啊。。

这个可以写的地方有很多，比如cloudfunction，或者cloud中，在solver中也能求。
也不知道你这个写在哪了。。
我是在particlecollector这个cloudfunction基础上写的 直接输出file 比较好后处理，

你如果想再paraview 中直接看的话，我们同事是在cloud中写这个field，然后在solver中对其进行平均，比如每100 iteration 平均一次。。

@李东岳 anaconda 下这个软件 然后什么包都不用装了 anaconda里面的spyder的ide用起来很简单，但是不能旋转图了，还是terminal出来的图的那个还可以旋转～
大佬不用python么？

可以多组分的，你得在thermo那边设置 mixture properties，然后在sprayclouldproperties中设置parcel的mixture的其他东西。
建议研究一下tut的reactingparcel，sprayparcel就是在reactingparcel的基础上增加了一些其他模型，也不明白为什么of将sprayparcel单独拿出来，不是放到reactingparcel哪个文件夹中。。

你用的是哪个版本
有的时候需要添加这么几句话在intermediate的options 中

在files中要对你自己intermediate的lib重新命名

@李东岳

https://pythonpedia.com/en/knowledge-base/38698277/plot-normal-distribution-in-3d

https://wizardforcel.gitbooks.io/matplotlib-intro-tut/content/matplotlib/31.html
把x3y3dz改为你的数据就行，剩下的就是设置图的其他东西了
为啥恢复间隔300s。。。

python

@dyb 图2 max Courant = 1.8，你用的PISO 还是PIMPLE
减小deltaT试试

@李东岳
大哥，你这样我压力好大= =