CFD中文网

@李东岳 谢谢 看来得一步步挖到底层了

想要把符合一定条件的粒子坐标信息输出，尝试用 p.x()但是无法编译成功，不知道正确的用法应该是怎么样？
另外，如何查看p到底包含有哪些子变量？
谢谢
ps：目前使用的是reactingParcelFilmFoam求解器，尝试在ThermoSurfaceFilm.C中修改

大家好，我之前尝试在unbuntu系统上安装source pack 遇到点问题，后来我按照https://openfoam.org/download/5-0-ubuntu/链接中的办法成功安装了openfoam5,之后在src/regionModel中做了些小改动编译也没有问题。
现在我尝试将reactingParcelFilmFoam中的reactingCloud改成sprayCloud生成一个新的求解器，编译的时候出现了以下错误，与之前source pack安装的错误一致：

     -lm -o /opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/sprayFilmFoam
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_comm_null'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_group_null'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_byte'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_double'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_op_min'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_comm_world'
/opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/openmpi-system/libPstream.so: undefined reference to `ompi_mpi_op_sum'
collect2: error: ld returned 1 exit status
/opt/openfoam5/wmake/makefiles/general:140: recipe for target '/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/sprayFilmFoam' failed
make: *** [/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/sprayFilmFoam] Error 1

不知道大家有没有什么建议？谢谢！

@李东岳 你好，如何在压力方程中添加源项 包含多相流的的表面张力项？有没有推荐的文档可以学习下？谢谢

@李东岳 问题已解决，直接在.H中添加了一段简单的二分法求解方程的程序，谢谢。

@李东岳 都还不知道，主要是想知道怎么植入，有没有什么建议？另外，对于类似问题 一个正确的解决思路是什么？谢谢

各位好， 我现在在尝试更改 OpenFOAM-5.x/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.H 中关于二次粒子射出方向（splashDirection）

目前对于出射角度phiSi在模型中是直接生成随机数来表征的，我想将phiSi改为如下第一个方程中的关系式

其中alpha是液滴入射角，可以通过液滴入射速度得到，但是需要通过解下面一个方程来求中中间变量gamma，初步打算用二分法来求解方程，但是现在有点不知道下一步怎么走，是直接用C++语法实现这个过程还是另外写个头文件或者还有什么更好的办法？提前感谢大家的意见和建议，谢谢！