Euler-Lagrange的一些解析，sprayFoam

[好学者]

@星星星星晴 @youhaoyu 谢谢两位大佬前辈的讲解，我也有了新的认识，感谢感谢！

[youhaoyu]

@星星星星晴 您好大佬，我询问一些事情不知道您曾经是否遇到过。我在加（p.d）颗粒长大时候，假如我颗粒直径是个随着时间线性增长的函数时候，能够进行源文件编译，求解器运行，形如我代码红框那个部分（其他地方是我原始的方程代码，我没有进行使用罢了，storeddt就是我的d）

但是我修改一下代码，改成pv=nrt这种就无法使用，能够通过编译，就是无法进行求解器运行，一下子就报错

**我在想难道是他无法去识别我这个Ro_ pgo_，Pc_？？？**这个Pc是我在求解器里自定义的p*rho。但是我在其他文件下面仿造openfoam自带的一些颗粒属性都对这些进行了定义啊，这个问题困扰了我好几个月了

[星星星星晴]

@youhaoyu 应该是p.d()什么的吧。如果你在kinematicparcel.h等地方定义了的话，你应该通过p.xxx()赋值。另外如果你需要流场速度是可以通过两个方式获得，p.uc()插值的值或者通过mesh获得网格中心的数据吧

你的问题很具体，我也很久没接触OF了。所以建议你多去google。

[星星星星晴]

比如你的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]

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