@李东岳 感谢李老师的回复，最近这些天也在想原因，感觉其中一点还是粒子数太多，意外错误较多。之前算过小尺寸的算例计算量小出现意外状况的比较少。正在考虑能不能通过模型模化来解决。

@李东岳 oh, cool!

问题解决，this->alpha_的定义在

template < class Alpha, class Rho, class BasicTurbulenceModel, class TransportModel > class TurbulenceModel : public BasicTurbulenceModel { public: typedef Alpha alphaField; typedef Rho rhoField; typedef TransportModel transportModel; protected: // Protected data const alphaField& alpha_; const transportModel& transport_; .......... }

这里的alpha_应该表示相分数。

第一，可以使用该求解器求解Euler方程

第二，rhoCentralFoam是一个可压缩问题的显式求解器，比较常见的手段为，将对流项，粘性项取t_n时刻的值，然后进行时间推进（Runge-Kutta方法等）

$$
\frac{\partial\pmb{U}}{\partial t} = (-\int_{\partial V}[F,G,H]\cdot d\pmb{S} + \int_{\partial V} [F^\nu,G^\nu,H^\nu]\cdot d\pmb{S})|_{t = t_n}
$$
但是，显式求解可能不稳定（或者说时间步长和速度，单元大小息息相关）。在rhoCentralFoam中处理\mu\nabla\cdot\nabla\pmb{U}时似乎用fvm::laplacian，这在某种程度上提高了计算的稳定性。

@李东岳

10.1615/AtomizSpr.2017019354
可能说的是这个东西

请问关于二元液体混合的求解器有哪些，并且混合后会形成混合物，可以自定义混合物rho，miu与alpha的关系。

请问可以把twoLiquidMixingFoam的详细解析发一下吗，这个是我的邮箱2693244245@qq.com，感谢！！

@heike256
都比我强。。。有的麻烦是自己人为的：如输入OFdir。。。直接 OFdir 下运行，省了输入，不是更简单。。。

https://cpp.openfoam.org/v9/timeActivatedFileUpdate_8H_source.html

fileUpdate1 { type timeActivatedFileUpdate; libs ("libutilityFunctionObjects.so"); writeControl timeStep; writeInterval 1; fileToUpdate "$FOAM_CASE/system/fvSolution"; timeVsFile ( (-1 "$FOAM_CASE/system/fvSolution.0") (0.10 "$FOAM_CASE/system/fvSolution.10") (0.20 "$FOAM_CASE/system/fvSolution.20") (0.35 "$FOAM_CASE/system/fvSolution.35") ); }

总的来说对于复杂网格的单一方向拉伸，of是不是没法通过简单的命令进行处理呢？

这个看起来会导致非常差的网格质量。如果实现的话，应该也是不能用于CFD计算。

但对于snappyHexMesh生成后发现多个方向都成倍数扩展了，这个不知道是什么问题

我估计是不是做了网格质量控制，发现网格太差了，然后全部缩放了？log里面有说么

这个是我用openfoam10跑的。这个数据还是不太好。我看了一下原始文献以及论坛里面的一些文章。这个跟湍流模型很敏感。我尝试了一下调节进口的湍流强度，也很敏感。这个我后续不会再调了。我觉得应该留给专门做这方面的人去调试一下

impingingJet10.tar.xz

具体网格如何分解这种的算法我就不知道了。scotch与metis我也只是用。不清楚怎么分解的。如果有哪个老铁知道。我愿意学习学习。

@luofq-sysu 现在普适性CFD软件完全打不开市场。国内大多走的行业化。所以一个普适性软件就变成了行业定制软件了。

@青山 我做的不涉及三相，所以只有两相。然后你这个说只能在alphat里算感觉没道理啊，那如果是层流，这个输入都没有，就不允许闪蒸了嘛，应该不会把

@vbcwl 6

@知行 在 paraview如何绘制二维流线 中说：

StreamTracerWithCustomSource

StreamTracerWithCustomSource是在surfaceVectors下创建的

@vbcwl 原来是因为硬盘里权限问题，不能出现带：的文件:136:

water { transportModel HerschelBulkley; nu0 1e-2; tau0 8.59e-3; rho 1000; k 18.24e-3; n 0.2; //Bingham }

damBreak.tar

@李东岳 code是个好东西啊，得想想怎么拟合曲线咯，谢谢李老师的参考案例啦

各位大佬们好，我想在fsifoam流固耦合求解器中添加两种固体材料属性，请问有人做过相关的算例嘛？是否可以直接添加？还是需要修改求解器呢？