泡沫的仿真

[WeakForm]

各位前辈，要研究高粘度流体表面形成的泡沫层，有一些困惑：
泡沫是典型的流变流体，看见一些文献单独研究泡沫流动，用CFD计算+自定义剪切力模型，误差很大，想问如果整体流动包含表面的泡沫层，是不是采用TFM算是不是不合理，无法解析出来表面的泡沫，且计算误差会比较大。

[李东岳]

泡沫这个目前来看，只能是单独当做一相，就像你说的，自定义剪切力模型。把泡沫看作一种流体。

如果多相流，我还没见过怎么能够预测表面的泡沫的。要是你看见过，可以分享下我研究下思路。

TFM确实解析不出来泡沫

[WeakForm]

@李东岳 谢谢李老师回复，老师我最近看到一个开源项目，他的文章里面提到了multi-vof方法能够解决多个气泡不合并的问题，这个是他们的github地址：https://github.com/cselab/aphros
这个是文章的doi:10.1126/sciadv.abm0590
文章里面的fig6看起来解析了表面泡沫层。
还没有仔细研读，希望老师有时间有兴趣的话看一眼，给出一些建议，谢谢李老师。

[WeakForm]

@weakform https://cselab.github.io/aphros/pdf/pasc2020.pdf 这个文章对算法描述更具体

[李东岳]

我研究下

[李东岳]

https://www.ljll.math.upmc.fr/~frey/papers/bubbles-coalescence/Coyajee E., Numerical simulation of drop impact on a liquid–liquid interface with a multiple marker front-capturing method.pdf

我看了一下这个文章。目前在openfoam里面的拓展，需要1）多流体、2）最主要的是openfoam目前是VOF，不带有marker。

目前不能给你提供更深层次的算法的思路，需要看一下VOF这面怎么弄。感觉不容易。

VOF这面模拟泡沫的有文章么

[WeakForm]

@李东岳 多谢老师，我还没有看到相关的，如果有再询问老师

[WeakForm]

捞一下这个帖子，今天看公众号pof年度文章第二篇用multiphaseeulerfoam算了啤酒泡沫：Experimental and numerical investigation of beer foam ，当然了就是看体积分数

[李东岳]

@WeakForm 屌爆了啊！明天我详细看看，感觉是个神文！

[WeakForm]

@李东岳 期待李老师点评☺️

[李东岳]

@WeakForm 他这个文章里面我前后挑数值的部分，看了3遍吧算。也没太细看。目前能看出这些：

1. 目前能看出来就是用多流体模型来模拟，泡沫就当做相分数比较低罢了；没有调用我以为的什么新思路来处理。不过文章至少证明了这种思路也可以（我之前以为不可以）
2. 但是没有说里面反应是什么？看方程里面存在反应，但是没写反应方程
3. 从算例里面来看，好像就是进口注入beer，然后体积分数=0.8，然后就算就行了。

目前不确定他这个里面的反应部分起的什么角色。

• 如果泡沫跟反应没关系。如果只是注入0.8的相分数，然后通过沉积来实现泡沫的消失，那好像就非常非常简单了。
• 如果泡沫跟反应有关，但是没看到反应方程在哪

[WeakForm]

@李东岳多谢李老师，确实我看的也是有点懵，感觉写的不怎么清晰😂，尤其是边界还有反应

[WeakForm]

@WeakForm 我要尝试一下multiphaseeulerfoam能不能模拟出来粘性气泡流表面的泡沫层😂

[WeakForm]

这个文章我感觉表面能看到气泡原因是入口含气率非常大80%，所以表面的泡沫与其说是泡沫不如说是一堆水雾，因为泡沫里边含水极低。我用open foam和fluent都不能模拟出来和高粘度流体实验中类似的表面泡沫层，李老师这边还有啥建议吗，边界设置能够让气泡在液面不会马上逸散的感觉。

[李东岳]

你模拟的是高粘度流体么？高粘度流体更难模拟了。

[WeakForm]

也不知算不算高粘度，1-几百mpas都有，感觉速度对的还行，液体表面泡沫层是真的没有