关于多域求解NS方程和Darcy方程的边界耦合问题

[李东岳]

我想模拟区域A的泄压过程，且不想使用OF自带的porous模块，因为该模块直接在方程中添加源项，仍需求解NS方程，计算复杂，求解时间较长。因此想直接在区域A求解NS方程，在区域B求解Darcy方程，并在交界面处耦合速度和压力。

你看一下无痛苦ns方程笔记最新版里面的4.8.6节，可以不适用多计算域求解器，直接实现你想做的事情

[huangyuhui723]

@李东岳 李老师，非常感谢您的解答，我去看了您说的内容，请问您的意思是直接自行植入Darcy-Brinkman-Stokes 方程，然后在空腔和多孔介质区域设置不同的孔隙率么（空腔处孔隙率设置为1，多孔介质处孔隙率设为材料实际的孔隙率）？如果是的话，我可能还有一个问题，那么在程序实现中如何设置空间分布不均匀的常数（孔隙率）呢？方便的话希望您指点一下，非常感谢。

[李东岳]

那么在程序实现中如何设置空间分布不均匀的常数（孔隙率）呢？

其他的你理解的都正确。这个问题，你定义孔隙率场，然后setFields就行了。

[huangyuhui723]

@李东岳 我大概了解了李老师，我去尝试看看，非常感谢。

[Do1975]

想请教下，为什么孔隙率在(0,1)之间，方程4.253会演变成Darcy方程？当孔隙率趋于0时，方程4.253左边第二项不应该趋于无穷大嘛？

[huangyuhui723]

@Do1975 不好意思，我当时也没有仔细推导 ，这么一说确实有一点不理解这部分。我当时尝试使用DBS方程后发现不太符合我的算例，然后我还是使用了多域求解的方法。

[Do1975]

@huangyuhui723 印象中当时用DBS是可以模拟出想要的效果，但是想进一步应用时候，却无法满足其它计算了。后续也打算试试多域求解的方法。

[huangyuhui723]

@Do1975 确实，可能还是看具体需求场景，我这是一方面DBS方法速度提升好像不多，然后在AB区域交界附近流场会出现间断，我没解决这个问题，就没用了

[ShuanmingZhao]

@huangyuhui723 你好，我最近也在做类似的工作，在三维情况下向空腔中充入液体，空腔中的气体可以通过周围的多孔区域排出，我的做法也是魔改chtMultiRegionFoam，我想请问你使用的边界方法以及是否成功？

[huangyuhui723]

@ShuanmingZhao 你好，我这个工作已经公开发表了：https://pubs.aip.org/aip/pof/article-abstract/36/10/104111/3315665/Multi-region-coupling-computational-fluid-dynamics?redirectedFrom=fulltext
然后代码也开源了：https://github.com/rararahui/chtPorousMultiRegionFoam
不过我水平有限，仅供参考哈

[ShuanmingZhao]

@huangyuhui723 感谢你的分享，明天我研究一下，有问题的话可能还得麻烦你

[huangyuhui723]

@ShuanmingZhao 好的，客气了，一起进步。

[李东岳]

@huangyuhui723 屌爆了大佬！

[huangyuhui723]

@李东岳 多谢李老师鼓励，还得感谢李老师的讨论答疑和提供的交流平台 。