@白猫紧张 感谢了，我学习下。感觉我们最终也得走上自己开发的道路，边界层的处理基本是几何建模和网格的最后一个难关了。:140:

@AppleKiller 处理不了。不过你可以直接画多套网格，每套网格求解不同的方程，这个很成熟，自己写都很容易。

@AppleKiller 在 porousGasificationFoam生物质热解模拟 中说：

@xuanze 这个求解器的固相会参与化学反应吗？

会参加化学反应，这个文件夹的solidReactions就是设定固相反应的

老师同学们好，使用foam-extend-4.0编程一个新的速度边界条件时怎么可以获得patch上速度的场分布?通过网上查找资料知道当编辑速度边界条件时可以通过以下代码获得alpha场的分布：
未命名图片.png
现在的问题是编辑一个新的边界条件时候怎么可以获得本身场在patch上的分布？
提前谢谢各位同学老师！！1

@AppleKiller 是的。

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

@李东岳 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