我没用过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中查看即可

   

@cccrrryyy 暴露年龄之
楼主好人，一生平安

课程材料原地址是：
https://www.tfd.chalmers.se/~hani/kurser/OS_CFD/

个人认为这里面提供的一些材料很有帮助。

前两天有一次不知道怎么回事，这个网站登不进去了，于是整个爬了下来，可以解压缩后，打开index.html即可。

链接: https://pan.baidu.com/s/1qvtJRMlIj-V4cAcAX6rlaA 提取码: j9gf

@jasper-0 这就得你自己翻code了，得看你用在算什么，dense的话需要考虑体积比什么的了，看起来这玩意就是和dense的有关

你指的SU应该是动量方程里面的源项，就是2-way coupling中的颗粒对流场的影响，那边是包涵了所有的颗粒对流场的影响的，具体你得去code里面查到底是怎么算的。

如果你考虑的是颗粒受力的话，你就得去submodel里面看了。

植入颗粒的话，下面这个教程可能说的更清晰
https://www.foamacademy.com/wp-content/uploads/2018/03/particles_slides.pdf

至于你说的phicForces 因为没做DPMFoam 所以不知道具体是什么，不过看起来这个好像用在很dense的颗粒模拟里面吧

@dxl 就是foamtovtk以后，在paraview中读取lagrangefield，然后加一个glyphy

会不会是质量有问题？sum变成了0.5就不对了吧

@vbcwl

https://www.cfd-online.com/Forums/openfoam-programming-development/241034-number-particles-inside-cell.html

const polyMesh& pmesh = this->owner().mesh();
List<DynamicList<typename CloudType::parcelType*>>& cellOccupancy =
    this->owner().cellOccupancy();

// check all cells containing particles
forAll(pmesh.cells(), celli)
{
    if ( cellOccupancy[celli].size() == 0 )
    {
        // do something
    }
    else
    {
        // do something else
    }
}

@小狗狗 DPM中半径不可以为0

@小狗狗 你这么说的话 就是质点， 不考虑半径

@小狗狗 我没处理过这样的问题，无法给你准确答案。但是我认为是不可以的，就类似你计算drag的时候同样也要考虑半径，当你考虑空间关系的话，应该也需要将其考虑为真正的小球。比如用拉格朗日模拟气泡壁面的话，肯定也要考虑气泡与避免的距离一样。