Euler-Lagrange的一些解析，sprayFoam

[星星星星晴]

比如你的particle diameter 应该通过p.d()获得。你的rhoc 是打算使用网格中心值还是考虑是有parcel点的插值？都是有算法的。

整体来说，你这个算法不能这么直接从别的地方生搬硬套。要根据你用的cloud，parcel，根据怎么获取parcrl这个struct内的内容来处理。建议你一步一步来，一行一行写，编译，测试。不能一口吃个胖子

[youhaoyu]

@星星星星晴 比如of他自带Uc插值，是网格中心mesh的，但没有压力P插值，我人为在求解器中间写了个Pc（模仿createField.H的Uc写的），然后在源库里面使用跟着UInterp模仿写PInterp，我这个思路我感觉应该是正确的吧，包括加什么粒子新的属性，我都是模仿他们原有的p.age这些去编写，但是最后就是有问题。
google有什么更好的搜索方式吗，我之前尝试过，请问还有什么平台我可以去询问，cfdonline我打算也问问哈哈哈哈哈

[星星星星晴]

@youhaoyu 你这个想法是对的，我们也是那么植入p的interpolation的。
cfdonline一般没人回 就是多查 有空可以看看chalmers 的openfoam的课，里面有很多学生的project。很有帮助

[星星星星晴]

@youhaoyu 因为你是基于kinematic parcel做的，说不准你做的很多东西都可能在别的地方植入过。
有的时候你觉得你植入对了，但是很有可能你并没有，这边差一点，那边差一点。

而且有一些模型并没有在kinematiccloud中真的植入。有可能就是一个空的头文件。

最好你植入一点，查一点。我建议你捋顺整个of的library，在哪个地方call cloud，cloud在什么地方call 的什么方程，到底call了几次。慢慢查，整着整着就熟了。

我之前可能在某个帖子里面发过我当时对lagranian这部分的理解的流程图，你可以借鉴一下。
不要纠结于具体某个模型，先搞懂框架，然后思考你要植入的东西具体应该植入到哪边，然后通过看看相似的其他植入，具体问题具体分析吧。git有很多repo 有的论文作者会分享出来的。

[Thomas_Huster]

@星星星星晴 大佬你好，就是我按照那个gofun在OpenFOAM9中编译的时候其OpenFOAM4中basicKinematicCollidingCloud.H的文件改名为了collidingCloud.H，但是编译一直告诉我找不到collidingCloud.H

[星星星星晴]

@Thomas_Huster
你在这include了文件，但是在你的make文件里面并没有提供相关的目录和地址，所以程序自然找不到在哪。你需要调整你的make文件的路径

[Thomas_Huster]

@星星星星晴 大佬你好，我试了下，但是不知道具体提供应该在哪添加