关于多域求解NS方程和Darcy方程的边界耦合问题
-
各位老师好，我已使用OpenFOAM一段时间，但是对底层算法和程序了解有限，最近想要实现一个多域求解的问题，但是有些问题无法想通因此发帖求助。
该算例情况如下：
上图为二维计算域示意图，区域A为空腔，初始时刻压力较高；区域B为多孔介质区域，初始时刻压力较低，外部为出口。
我想模拟区域A的泄压过程，且不想使用OF自带的porous模块，因为该模块直接在方程中添加源项，仍需求解NS方程，计算复杂，求解时间较长。因此想直接在区域A求解NS方程，在区域B求解Darcy方程，并在交界面处耦合速度和压力。
在论坛中搜索发现 https://www.cfd-china.com/topic/779/求助-如何使用openfoam进行多孔介质内数值模拟 该算例与我的问题类似，但是仍无法解决我的问题，我目前有两种思路想询问各位老师是否可行：
1、魔改chtMultiRegionFoam求解器，将区域B的热传导方程改为Darcy方程，将耦合变量改为速度和压力；
2、新编写求解器，直接在一个求解器中分两步分别计算两个算例（区域Ａ的NS方程和区域Ｂ的Darcy方程），每一个时间步计算完后映射更新交界处的速度和压力。
想请教各位老师以上思路的可行性和实现难度，如有更多思路和建议就更加感谢了。
谢谢，麻烦各位老师不吝赐教。 -
@huangyuhui723 印象中当时用DBS是可以模拟出想要的效果，但是想进一步应用时候，却无法满足其它计算了。后续也打算试试多域求解的方法。
-
@huangyuhui723 你好，我最近也在做类似的工作，在三维情况下向空腔中充入液体，空腔中的气体可以通过周围的多孔区域排出，我的做法也是魔改chtMultiRegionFoam，我想请问你使用的边界方法以及是否成功？
-
-
@huangyuhui723 感谢你的分享，明天我研究一下，有问题的话可能还得麻烦你
-
@ShuanmingZhao 好的，客气了，一起进步。
-
@huangyuhui723 屌爆了大佬！
。