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

   

@chengan-wang 上面已写具体做法

@chengan-wang 1-way 就是 parcel收到流场的影响，2way的话就是流场也收到parcel的影响
newFormat 就是chemkin库，具体这部分你要去thermo查，
如何计算某个温度下的表面张力和粘度什么的 polynomial 还是什么的，具体参数都是什么，都需要你自己查的。 这玩意是有标准库的。
chemkin不是很熟，我平时都不用这个，文件里面写了

Parcel和Particle的区别，两者并不相同

在使用Parcel概念的时候一定要注意Parcel 和Particle 是有区别的
Parcel是Particle的合集，所以diameter和velocity是相同的，但是质量是所有particle的和

$m_{parcel} = \sum m_{particles}= \rho_p V * np$

所以再考虑到parcel质量的时候一定别忘了np
在考虑mass的时候也别忘了体积和直径是立方关系。
再做统计的时候一定要考虑好到底是做parcel的统计还是particle的统计。

统计某一个平面的parcel的相关信息

因为of里是单独将lagrangianfield 拿出来的，所以说，再paraview中是无法使用plot over line查看拉格朗日场的。于是解决方法通常有两种

1. 将拉格朗日场转到Euler场

template<class CloudType>
Foam::volScalarField& Foam::IVT8_StochasticCollisionModel_2022<CloudType>::PPC() {
	if (!PPCPtr_.valid() ) {
		const fvMesh& mesh = this->owner().mesh();
		PPCPtr_.reset
		(
		    new volScalarField
		    (
		        IOobject
		        (
		            this->owner().name() + ":PPC",
		            mesh.time().timeName(),
		            mesh,
		            IOobject::READ_IF_PRESENT,
		            IOobject::AUTO_WRITE
		        ),
		        mesh,
		        dimensionedScalar ("zero", dimless, 0.0),
		        zeroGradientFvPatchScalarField::typeName
		    )
		);
	}

	return PPCPtr_();
}

然后再调用的时候这么调用即可
scalar& PPC = PPCPtr_->primitiveFieldRef() [cellI];
PPC += 1.0;

2. 输出文件

然后用python或者什么后处理。
具体可以参阅ParticleCollector，再particleCollector下的postMove中添加如下code, 具体参数就根据自己的具体要求添加把。

     forAll (hitFaceIDs_, i) {
            std::ofstream outfile;
            outfile.setf (ios_base::fixed, ios_base::floatfield);
            outfile.precision (8);
            outfile.open ("output.txt", ios_base::app);
            outfile << this->owner().db().time().timeName() << tab
                    << p.origId()                       << tab
                    << p.typeId()                       << tab
                    << p.age()                          << tab
                    << p.position().component (0)       << tab
                    << p.position().component (1)       << tab
                    << p.position().component (2)       << tab
                    << p.U().component (0)              << tab
                    << p.U().component (1)              << tab
                    << p.U().component (2)              << tab
                    << p.d()                            << tab
                    << p.Rep()                          << tab
                    << p.nParticle()                    << tab
                    << p.mass() * 1e18                  << tab
                    << nl;
        }

两种方法利弊：
方法1，
pro: 简单快捷，不需要大量的算力，会吃一些内存，但是完全取决于你的cell,可以再parview中简单直接的看到随时间的变化
con: 当你需要进行一些统计学的运算的话，需要建立一些fields，而且无法精确统计每个parcel， 另外如果需要做时均的话，会损失一些信息。。（血与泪的教训）

方法2，
pro: 每一条输出就是一条数据，你可以输出任何你想要的数据，做任何统计学上的后处理
con: 需要后处理，可能会输出几十上百G的file，python或者matlab等后处理，需要一定的时间。

aachanbomb_sprayFoam设置解析.pdf

默认的Lagrangian solver

particleFoam, rhoParticleFoam 这两个是不解流场的，所以在使用这个solver的时候，需要有一个frozen的流场或者如果单纯解拉格朗日的话，就不用考虑流场的东西了

所以只有one-way coupling对应之前版本icoUncoupledKinematicParcelFoam和uncoupledKinematicParcelFoam

其他几个solver是解流场的，除了sprayFoam是用过的，其他的就不献丑了

因为要写给project partner关于我的一部分工作，所以写了这个东西，仅涉及了默认的OpenFOAM的一些东西，不能说是本着开源精神，毕竟自己做的东西并没有involve，但是希望能对论坛里做欧拉拉格朗日的小伙伴们一些帮助，尤其是新手。此贴会长期更新，一点一点来。

本人也是资历有限，在OF里一番摸爬滚打，不免会有理解错误，望见谅～
因为工作不涉及到thermo部分，所以这部分基本不会。。哈哈

This document is just for the explaination of the sprayFoam.pdf

@香柏树 另外你直接把整个PIMPLE循环给关掉了，是不是保留一下然后PISO模式运行就行呢, 瞬态

@香柏树 除以0可能是最大的问题，但是0在哪呢，请问你怎么处理parcel 极小的部分的，会不会是parcel部分，因为你点火了，就出现了这种情况，在clockTime以后就是计算Co 然后solving parcel了，也有可能是你parcel的maxCo=1， parcel在你计算的时候已经飞出了对应的cell 也有可能报错，尝试分析分析这俩地方，看看你的parcel都什么样，看看你的MaxCo减小到0.1或者0.01看看，还是turbulence爆了

@chengan-wang 请善用谷歌，这些东西在网上都有，只是零碎的，等我忙完这段时间或许会写个关于喷雾模拟的东西，但是现在请善用搜索引擎。而且你需要多读code才明白为什么这个东西是这么设置的

请见谅说的有些直，但是事实如此。

@田畔的风 兄弟，给个原生编译的教程呗？我这边只是M1的air 和mini 还是跑的虚拟机