使用OpenFOAM自带的porousSimpleFoam实现多孔介质仿真
达西定律
多孔介质在石油化工、多相传热乃至航空航天等领域都有重要的研究以及应用价值。流体通过多孔介质时，其最突出的特征就是压强的降低，以及多孔介质的孔隙率、穿透性、流体粘度，并且介质的长度都有着定量的联系，这就是所谓的“达西定律”（Darcy’s Law）。
porousSimpleFoam的核心公式
porousSimpleFoam 是从稳态不可压求解器simpleFoam改写过来的，其使用SIMPLE算法来处理速度压强的耦合问题。它处理多孔介质的原理是在动量方程中添加源项来实现的，具体有以下几类多孔介质模型：
No.1：Flow through a porous plug
如上图，先是生成了cellSet体元集，在sourceInfo中设置多孔区域的左下角以及右上角，这样就划分了多孔区域了。然后再把cellSet转成cellZoneSet，为下一步作准备。
02 多孔参数设置
porousSimpleFoam求解器计算得到的压强场：
No.2 : Beavers and Joseph problem
Beavers and Joseph problem同样是管道中的算例，不一样的是，多孔介质区域布置在管道的下半区：
porousSimpleFoam求解器计算得到的速度场：
沿x=0速度曲线：
如图所示，porousSimpleFoam求解器计算结果与文献（Li，2020）的数据结果吻合较好，说明该求解器可以很好的计算Beavers and Joseph problem。
【参考文献】
Li Z, Zhang H, Liu Y*, McDonough JM. (2020) Implementation of compressible porous-fluid coupling method in an aerodynamics and aeroacoustics code part I: Laminar flow. Applied Mathematics and Computation, 364: 124682. (JCR Q1)
算例下载链接：见附件
