什么是Brinkman方程？

在描述多孔介质流动的达西定律中加入黏性扩散项的方程。达西定律无法满足壁面或多孔介质与自由流体界面上的速度条件，而Brinkman方程可以施加壁面无滑移条件，因此能表示边界附近的流动。

达西定律与Brinkman方程何时结果不同？

当渗透率K的平方根与通道宽度相比不可忽略时。10 mm宽通道中，K=10⁻⁸ m²时流量差2%，10⁻⁶ m²时20%，10⁻⁵ m²时58%。泡沫金属和粗颗粒填充床需要用Brinkman方程。

有效黏度怎么确定？

它表示多孔介质内部黏性的作用，常取与流体黏度相同，但实验中根据孔隙率和结构报告过1到数倍的值。示例中有效黏度加倍流量约减少10%，重要时应通过实验或孔隙尺度分析确认。

Brinkman方程

分类：流体分析（CFD） | 2026-09-30 改写

CAE visualization for brinkman equation theory - technical simulation diagram

Brinkman方程式

理论：在达西定律中加入黏性项

概述

🙋

我用达西定律计算填充泡沫金属的冷却通道，预测流量比实验多。缺了什么？

🎓

多半是壁面的影响。达西定律是“压力梯度与速度成正比”的公式，速度处处相同呈平直分布，无法表示壁面速度为零。像泡沫金属这样渗透率大的材料，壁面附近流动变慢的层厚度就不能忽略了。Brinkman方程在达西定律中加入与普通流动相同的黏性项，从而能表示这个壁面层。

达西定律与Brinkman方程

$$ \nabla p = -\frac{\mu}{K}\,\mathbf{u} $$

$$ \nabla p = -\frac{\mu}{K}\,\mathbf{u} + \mu_{eff}\,\nabla^2 \mathbf{u} $$

$K$ 为渗透率（m²），$\mu$ 为流体黏度，$\mu_{eff}$ 为有效黏度，$\mathbf{u}$ 为表观流速（整个截面上的平均流速）。加入右边第二项后，就能在壁面施加速度为零的条件。$\sqrt{K}$ 具有长度量纲，是壁面影响范围的参考。

平行平板通道的解析解

$$ u(y) = \frac{G K}{\mu}\left[1 - \frac{\cosh(y/\sqrt{K})}{\cosh(h/\sqrt{K})}\right], \qquad \frac{Q}{Q_{Darcy}} = 1 - \frac{\tanh(h/\sqrt{K})}{h/\sqrt{K}} $$

这是在宽 $2h$ 的多孔介质填充通道上施加压力梯度 $G$、有效黏度等于流体黏度时的解。$h/\sqrt{K}$ 大时趋近达西定律的平直分布，小时趋近无多孔介质的抛物线形泊肃叶流动。

Coffee Break 随笔闲谈

从喷泉到地下水，再到颗粒群

达西定律是1856年法国第戎负责设计喷泉供水的工程师达西在测量水通过砂层的流量时发现的。1947年，Brinkman为计算流体流过大量颗粒聚集体时的阻力，在达西定律中加入了黏性项，试图用一个方程把颗粒周围的流动与整个多孔介质的流动联系起来。如今它广泛用于壁面和界面影响大的场合，如泡沫金属换热器、燃料电池电极、生物组织内的流动等。

算例

例1：渗透率与流量偏差

宽10 mm（$h$ = 5 mm）通道填充多孔介质时，相对达西定律和相对泊肃叶流动的流量比：

渗透率 K	√K	达西数 K/h²	与达西之比	与泊肃叶之比
10⁻¹⁰ m²	0.01 mm	4×10⁻⁶	0.998	0.00003
10⁻⁸ m²	0.10 mm	4×10⁻⁴	0.980	0.0012
10⁻⁷ m²	0.32 mm	4×10⁻³	0.937	0.011
10⁻⁶ m²	1.0 mm	0.04	0.800	0.096
10⁻⁵ m²	3.2 mm	0.4	0.419	0.50
10⁻⁴ m²	10 mm	4	0.076	0.91

🙋

渗透率越大，达西定律高估流量越多。

🎓

对。砂和土（约10⁻¹⁰ m²）的壁面层只有0.01 mm，达西定律就够了。泡沫金属和粗颗粒填充床（10⁻⁶ m²以上）会使流量高估20%以上。反过来在10⁻⁴ m²时，多孔介质阻力几乎不起作用，流动接近空通道的泊肃叶流动——只用达西定律会高估13倍。Brinkman方程能在这两个极端之间平滑过渡。

例2：壁面层厚度

K = 10⁻⁸ m²（√K = 0.1 mm）时，速度在距壁0.23 mm处达到达西值的90%，0.46 mm处达到99%。壁面层只有 √K 的几倍厚，CFD中必须在这个厚度内布置几个以上网格，才能正确表示壁面阻力。

例3：有效黏度的影响

有效黏度 / 流体黏度	壁面层厚度参考	与达西的流量比
1	1.00 mm	0.800
2	1.41 mm	0.718
5	2.24 mm	0.563

K = 10⁻⁶ m²、宽10 mm的情况。有效黏度难以测量，常取1，但加倍后流量约减少10%。在壁面影响大的设计中，要报告该假设对结果的影响。

分析方法

模型选择

情况	使用的方程
√K 小于通道宽度的1/100，流速慢	达西定律
√K 不可忽略，有壁面或与自由流体的界面	Brinkman方程
流速快，惯性阻力起作用	加入Forchheimer项

CFD注意事项

表观流速与实际流速：方程中的 $\mathbf{u}$ 为截面平均流速，孔隙内实际流速要除以孔隙率。
多孔介质与自由流体界面：速度连续条件与允许滑移的条件（Beavers–Joseph）会给出不同结果。
渗透率各向异性：纤维和层状材料不同方向的渗透率可相差数倍，应以张量给出。

Coffee Break 随笔闲谈

“流体从壁边溜走”

在填充球体的填充床中，只有壁面附近孔隙率变高，流体容易沿壁边溜过：球沿壁面排列，距壁约半个球径处空隙变大。管径小于球径约10倍时，这种壁边流动会明显影响整体流量和传热。把Brinkman方程与随距壁距离变化的孔隙率结合，就能表示这个效应。

常见错误

错误与对策

错误	后果	对策
对高渗透率材料用达西定律	高估流量	比较 √K 与通道宽度来选择
壁面层网格太粗	壁面阻力出错	在 √K 范围内布置几个网格
混淆表观与实际流速	传热或停留时间出错	用孔隙率换算
不考虑有效黏度	流量差10%以上	确认敏感度并报告
高速流中忽略惯性	低估压力损失	加入Forchheimer项

🙋

我想自己试试渗透率和摩擦的计算。

🎓

可以用渗透率与达西定律工具试算渗透率和达西定律，用达西与范宁摩擦系数试算管道摩擦。相关方程和应用详见达西流动、Forchheimer方程、填充床和地下水流动等文章。