Herschel-Bulkley 模型

我们日常生活中的许多材料都难以简单地归类为固体或液体。想想番茄酱或牙膏：它们在重力作用下能保持形状，但一经挤压就很容易流动。这些材料被称为屈服应力流体，它们的行为可以通过一个强大的数学关系式——Herschel-Bulkley 模型——来捕捉。该模型填补了我们对一类物质的理解空白，这类物质在静止时具有类似固体的“顽固性”，但在足够大的力作用下会转变为类似流体的状态。本文将引导您深入了解这个引人入胜的主题。首先，我们将探讨“原理与机制”，剖析该模型的方程及其三个关键参数，并考察屈服应力的微观起源。随后，在“应用与跨学科联系”部分，我们将探索其在现实世界中的意义，从日常工程、大规模地质事件到复杂的生物系统力学。

想象一下，你正试图倒出瓶子里最后一点番茄酱。你把瓶子倒置，什么也没发生。你等着，还是什么都没有。番茄酱在自身重力作用下拒绝流动。它有一种顽固性，一种像水这样的简单流体所没有的运动阻力。你必须给它一个有力、快速的摇晃或挤压——你必须施加一个力来克服它最初的抗拒。只有这样，它才决定移动。这个日常的“斗争”是进入屈服应力流体世界的一个完美引子，而理解它们的关键是一个极其通用的优美方程，即 ​​Herschel-Bulkley 模型​​。

从本质上讲，Herschel-Bulkley 模型是你为使流体流动所施加的力（​​剪切应力​​，$\tau$）与流体因此流动的速度（​​剪切速率​​，$\dot{\gamma}$）之间的简单关系。对于大于某个阈值的应力，该关系式为：

$\tau = \tau_y + K(\dot{\gamma})^n$

这个方程可能看起来有点复杂，但它实际上就像一个有三个角色的故事，每个角色控制着流体行为的一个不同方面。

第一个也是最重要的角色是 $\tau_y$，即​​屈服应力​​。这正是番茄酱顽固性的量化。这是你必须施加的、使流体“屈服”并开始移动的最小应力。低于这个值，材料表现得像固体——它可能会有一点变形，但不会流动。牙膏能停留在牙刷上，涂料不会从墙上滴落，正是因为它们具有屈服应力。

一旦你支付了 $\tau_y$ 这个“入场费”，方程的其余部分 $K(\dot{\gamma})^n$ 开始生效，描述流体在流动时的行为。在这里我们遇到了另外两个角色。

参数 $K$ 是​​稠度指数​​。你可以把它想象成材料一旦开始运动后的“稠重度”或“厚度”。高 $K$ 值意味着流体阻力很大，就像试图搅拌浓稠的泥浆；而低 $K$ 值则意味着它流动得更容易，像浆液一样。关于 $K$ 的一个迷人之处在于，它的物理意义甚至其单位，都与第三个参数 $n$ 相关联。它不像水或蜂蜜的粘度那样简单；它的作用更加微妙，这暗示了其表面之下隐藏着复杂的物理学。

最后一个角色 $n$ 是​​流动行为指数​​，也是最有趣的一个。它描述了流动的个性。

• 如果 $\boldsymbol{n=1}$，方程简化为 $\tau = \tau_y + K\dot{\gamma}$。这描述的是 ​​Bingham 塑料​​。一旦屈服，应力随剪切速率线性增加，就像普通的牛顿流体一样，只是有一个初始的偏移量 $\tau_y$。浓稠的泥浆或粘土浆料通常表现出这种行为。

• 如果 $\boldsymbol{n 1}$，我们得到的是​​剪切稀化​​流体。这就是涂料背后的奥秘。当你更快地剪切它时（例如用刷子涂抹），它的阻力会变小——其表观粘度会下降。这就是为什么涂料可以足够稠，不会从刷子上滴落（低剪切速率），但在墙上又能顺滑地涂抹开（高剪切速率）。

• 如果 $\boldsymbol{n>1}$，我们得到的是​​剪切增稠​​流体。这就是“欧不裂”（玉米淀粉和水的混合物）的奇怪行为。如果你慢慢搅拌它，它是一种液体。但如果你猛击它或试图快速搅拌，它会突然变得几乎像固体一样。

Herschel-Bulkley 模型之所以如此强大，是因为只需调整这三个参数，它就能描述这一整族行为，从简单的 Bingham 塑料 ($n=1$) 到典型的幂律流体（通过设置 $\tau_y=0$）。它是对我们日常遇到的各种材料的统一描述。

这个宏观模型具有极佳的描述性，但作为物理学家，我们不禁要问：为什么？这种屈服应力从何而来？为了找到答案，我们必须放大并观察材料的微观结构。事实证明，屈服应力并非单一现象；它至少可以由两种不同的潜在物理机制产生。

想象一下一种稠密的颗粒悬浮液，就像一碗浓稠的微观珠子汤。

​​机制 1：粘性网络。​​ 如果颗粒之间存在某种吸引力——有点像微弱的小磁铁——它们可以形成一个相互连接、遍布整个容器的蔓延网络。这种相互锁定的颗粒结构赋予了材料初始的固态特性。它可以通过拉伸其键合来抵抗小力。屈服应力 $\tau_y$ 是全局性地破坏这个网络所需的单位面积上的力。一旦网络被破坏，颗粒就可以相互流过，但它们仍在相互作用，键合不断断裂和重组。你剪切流体的速度越快，网络被破坏得越严重，提供的阻力就越小。这种“结构弱化”是剪切稀化行为 ($n 1$) 的物理根源。这种机制——一个具有内聚力、会断裂的网络——被完整的 Herschel-Bulkley 模型完美地捕捉到。

​​机制 2：拥挤阻塞。​​ 现在，想象颗粒不具粘性，只是被紧密地挤在一起，并由外部压力（如沙袋中的沙子）固定。它们被阻塞了。要使它们流动，你必须施加足够大的推力，以克服颗粒间的静摩擦力，使它们能够相互滑动。在这种情况下，屈服应力来自于这种集体摩擦，并且与约束压力成正比（就像在地板上拖动一个重箱子比拖动一个轻箱子更难）。一旦它们开始滑动，流动的阻力主要来自于颗粒必须推开周围液体（粘性拖曳）。这导致了一种行为，即所需的应力等于屈服应力（来自摩擦）加上一个与剪切速率成正比的项（来自粘性）。这就是 Bingham 塑料 ($n=1$) 的微观图像。

这两种图像告诉我们，我们所说的“屈服应力流体”不是一个单一的类别，而是一系列具有不同起源的行为的集合。Herschel-Bulkley 模型是让我们能够谈论所有这些行为的数学语言。

这就提出了一个微妙但深刻的问题：屈服应力流体仅仅是一种非常稠的液体吗？牙膏只是非常非常粘稠的蜂蜜吗？答案是明确的“不”，而 Herschel-Bulkley 模型的数学原理向我们揭示了原因。

让我们将流体的​​表观粘度​​定义为 $\eta_{app} = \tau/\dot{\gamma}$。对于像水这样的简单流体，这是一个常数。对于没有屈服应力的剪切稀化流体（如聚合物溶液），这种粘度在低剪切速率下可能非常高，但它始终是一个有限的数值。这意味着，如果你施加任何应力，无论多小，它都会流动，尽管非常缓慢。

但对于 Herschel-Bulkley 流体，情况发生了戏剧性的变化。当剪切速率 $\dot{\gamma}$ 趋近于零时，表观粘度变为：

$\eta_{app} = \frac{\tau_y + K\dot{\gamma}^n}{\dot{\gamma}} = \frac{\tau_y}{\dot{\gamma}} + K\dot{\gamma}^{n-1}$

由于 $\tau_y$ 是一个正常数，当 $\dot{\gamma}$ 变得无穷小时，$\tau_y/\dot{\gamma}$ 这一项会激增至​​无穷大​​。一个真正的屈服应力流体，在这个理想化的模型中，对于最微弱的流动迹象都具有无穷大的阻力。它不仅仅是一种非常稠的液体，而是一种完全不同的东西。它拥有一种必须被打破的真正固态特性。这种数学上的发散是屈服应力的鲜明而清晰的标志。

有了这个模型，我们就可以开始探索这些流体在真实世界情境中的行为，并且很快就会遇到一些挑战我们直觉、加深我们理解的有趣“悖论”。

顽固的柱塞

想象一下我们的流体被泵送到一根长直管道中。压力推动它前进。我们对剪切稀化流体的直觉可能会认为，其奇特性质会影响整个流场。但在这里，物理学在这里巧妙地分为两部分：力平衡和材料响应。

首先，一个简单的力平衡（实际上是牛顿定律）告诉我们，管道内的剪切应力在管壁处最高，在正中心为零。这对任何流体都成立，无论是水、涂料还是薄饼面糊。应力分布是由力学决定的，而不是由材料决定的。

现在，我们引入材料的个性。在靠近管壁的区域，应力很高——高于 $\tau_y$——所以材料屈服并流动。但随着我们向管道中心移动，应力下降。在某个点，应力变得等于 $\tau_y$。在这个半径以内，应力低于屈服应力。当应力低于 $\tau_y$ 时，屈服应力流体会怎样？它不会变形。它表现得像一个固体。

结果是，一个类似固体的“柱塞”在管道中心向下移动，旁边是滑动的已屈服流体层。这个柱塞的大小仅取决于应力等于屈服应力的位置。它取决于压力梯度和 $\tau_y$，但——这是令人惊讶的地方——它完全独立于流动参数 $K$ 和 $n$。流体的剪切稀化特性影响已屈服部分的流动速度，但它不改变固体柱塞的边界。这是一个绝佳的例子，说明了分离原理——力平衡和本构定律——如何能够得出清晰明确但又出人意料的结论。

滑动的固体

让我们思考另一个难题。如果我们的屈服应力流体在一个管壁是“滑”的通道中会怎样？我们假设壁面上的应力与壁面处的流体速度成正比，而不是无滑移条件。这被称为​​Navier 滑移条件​​。

现在，我们只施加一个非常小的压力梯度，小到各处的应力，甚至在壁面处，都低于屈服应力 $\tau_y$。根据我们的模型，流体不能屈服。它必须表现为固体。我们的直觉会大声疾呼，如果它不能屈服，它就不能移动。

但我们的直觉是错的。本构定律 $\dot{\gamma} = 0$ 并没有说速度为零；它说的是速度梯度为零。这意味着材料可以移动，但它必须作为一个单一的、速度均匀的刚性块体在通道中移动。这种刚性平移与“固态”完全一致。但它移动得多快呢？速度由边界条件决定。柱塞上的压力由光滑壁面的摩擦力来平衡。非零的压力梯度导致非零的壁面应力，通过滑移定律，这要求速度为非零。

因此我们得到了一个非凡的情景：一个未屈服的固态块体在通道中滑下，其运动完全由边界处的滑移促成。这精美地说明了本构模型控制内部变形，而边界条件控制物体作为一个整体如何与周围环境相互作用。

Herschel-Bulkley 模型不仅仅是一个理论上的好奇心；它是工程和科学领域的主力。但要有效地使用它，需要判断力和对建模艺术的欣赏。

​​建模者的困境。​​ 假设你的实验数据清楚地表明一种流体是剪切稀化的（$n \approx 0.7$）。你需要运行一个计算机模拟。你是否应该总是使用更精确的 Herschel-Bulkley 模型？不一定。如果你的应用只涉及一个非常窄的剪切速率范围，一个更简单的 Bingham 模型 ($n=1$) 可能是一个完全足够的近似。一个更简单的模型在计算上可能更便宜、更稳健，使其成为特定任务的更明智选择。然而，如果你的模拟涉及宽范围的剪切速率，使用错误的物理原理（如假设 $n=1$）可能会导致截然不同的错误预测。教训是，没有一个单一的“最佳”模型；正确的选择是在保真度、成本和你试图回答的问题之间进行权衡。

​​实验者的陷阱。​​ 我们如何从真实数据中找出 $\tau_y$、 $K$ 和 $n$ 的值？一个处理幂律的常用技巧是对数据取对数，将其转化为一条直线。人们可能会尝试猜测一个 $\tau_y$ 值，从应力数据中减去它，然后绘制 $\ln(\tau - \tau_y)$ 对 $\ln(\dot{\gamma})$ 的图。如果你猜对了 $\tau_y$，你会得到一条斜率为 $n$ 的直线。但如果你的猜测略有偏差呢？数学表明，你的“直线”实际上会是一条微妙的曲线。对这条曲线进行线性拟合会给你一个有偏差、不正确的 $n$ 值。科学上诚实的做法是直接将完整的、非线性的 Herschel-Bulkley 方程拟合到原始数据上。这在计算上更困难，但它尊重了模型和数据的完整性，这是任何科学探索中至关重要的原则。

从番茄酱的顽固性到先进材料的设计，Herschel-Bulkley 模型提供了一种语言，来描述、理解和预测一类迷人材料的行为。它提醒我们，有时，最深刻的见解就隐藏在最熟悉的现象中，只等待一个能赋予它们声音的框架。

掌握了 Herschel-Bulkley 模型的基本原理之后，我们现在准备踏上一段旅程。我们将看到这些思想并非局限于抽象的方程世界，实际上，它们无处不在，在我们的厨房里，在我们的身体内，在宏伟而可怕的自然舞台上，以及在现代技术的最前沿。该模型真正的美在于其力量，它能统一大量看似无关的现象，揭示出一个共同的物理原理在起作用：静止时如固体、运动时如液体的迷人二元性。

我们的第一站是我们日常使用的材料世界。你是否曾想过，为什么高质量的涂料不会从刷子上滴落，却能顺滑、毫不费力地涂抹在墙上？秘诀就在于屈服应力。在静止状态下，在重力的低应力作用下，涂料的屈服应力 $\tau_y$ 未被超过，它表现得像一个软固体，附着在刷毛上。但当你施加刷涂的高剪切力时，你克服了这个屈服应力。涂料开始流动，并且因为它也是剪切稀化的（流动指数 $n \lt 1$），其表观粘度骤降，使其能够被涂抹成薄而均匀的一层。

同样的原理在高精度制造中也至关重要。考虑使用导电油墨丝网印刷电子电路的挑战。油墨必须足够液态，才能通过刮刀的高剪切作用被压过细网格，但一旦它附着在电路板上，就必须立即“固化”并以完美的保真度保持其形状，不会在重力和表面张力的低应力下渗色或模糊。具有高屈服应力和强剪切稀化行为的材料，在印刷过程中的表观粘度远低于其静止时，使其成为这项工作的完美选择。

这些材料的旅程并未就此结束。我们如何通过数英里长的管道输送像牙膏、混凝土或食品浆料这样的浓稠物质？泵送像蜂蜜这样简单的高粘度流体需要巨大的能量。在这里，Herschel-Bulkley 模型揭示了一个奇妙的反直觉现象：柱塞流。当屈服应力流体通过管道泵送时，剪切应力在管壁处最高，在中心为零。这创造了一种有趣的情形：靠近管壁的一层流体屈服并流动，而中心的一个大核心，其应力低于 $\tau_y$，则作为一个单一的固体“柱塞”移动。这个柱塞在边缘流动层的润滑下有效地滑动。这个固态核心的大小与流体的屈服应力成正比，与驱动流动的压力梯度成反比。这不仅仅是一个奇特的现象；它是有效设计输送无数定义我们现代世界的材料的基本原则。

保持涂料在刷子上的物理原理，也同样支配着地球上一些最强大的力量。让我们把思维从管道扩大到整个地貌。火山泥流（lahar）是一种由水、火山灰和岩石组成的可怕浆体。为什么它如此具有破坏性？一场普通的水灾或许能移动汽车，但一场火山泥流可以搬运巨大的石块、房屋和桥梁。原因在于其巨大的屈服应力。

一块巨石位于火山泥流中，它受到重力向下的拉力，但也被流体内部强度的支撑力向上推。这个支撑力是屈服应力的直接结果。要使巨石被悬浮并携带走，这个支撑力必须大于巨石的浸没重量。这意味着对于一个给定屈服应力 $\tau_y$ 的火山泥流，它能携带的岩石有一个最大尺寸。任何更大的岩石都会被沉积下来。地质学家正是利用这一原理，有时通过检查古代流留下的巨石大小，来估计过去和未来火山灾害的屈服应力和破坏力。

从地表深入地壳，我们在岩浆的流动中发现了另一个戏剧性的应用。岩浆并非卡通里那种简单的、均匀的红色液体；它是在液态硅酸盐熔体中悬浮着固体晶体的复杂混合物。这种混合物的行为类似于 Herschel-Bulkley 流体。流变参数——屈服应力 $\tau_y$、稠度 $K$ 和流动指数 $n$——不是恒定的，而是关键地取决于晶体的分数。随着岩浆冷却和结晶，其屈服应力急剧增加。这种行为决定了一切，从熔岩流的速度到岩浆侵入岩石形成岩脉的方式。通过对这种依赖于晶体的流变学进行建模，地球物理学家可以预测岩脉中固体岩浆“柱塞”的宽度，并更好地理解火山喷发的机制以及我们星球的深层管道系统。

Herschel-Bulkley 模型不仅描述了我们周围的世界，也描述了我们内在的世界。我们的身体是流变学大师，不断地操控着复杂的流体。吞咽这个简单的动作是一个深刻的生物力学挑战。为了将一口糊状食物从口腔移动到食道，舌头必须像活塞一样，产生足够的压力来克服食物的屈服应力。如果舌头产生的压力不足以使食物屈服和流动，吞咽就会失败。这一原理不仅是学术性的；在为患有吞咽困难（dysphagia）的个体设计食物时，它具有至关重要的临床意义，确保食物既安全易咽，又美味可口。

也许最著名的非牛顿流体就是我们自己的血液。乍一看，Herschel-Bulkley 模型似乎很适合描述血液。在极低的流速下，红细胞倾向于聚集成串，形成所谓的“缗钱状”，构成一个网络，使血液具有表观屈服应力。需要一定的力来打破这个网络并启动流动。一旦流动，这些缗钱状结构就会分解，单个红细胞会变形并与流动方向对齐，导致表观粘度下降——这是一种典型的剪切稀化行为。

然而，在这里我们必须学习一个关于物理模型本质的经典的费曼式教训。虽然 Herschel-Bulkley 模型提供了一个有价值的初步近似，但它并非故事的全部。与我们的简单模型不同，血液的性质是随时间变化的（触变性）——缗钱状网络的形成或破坏都需要时间。此外，该模型预测在高剪切速率下粘度应无限降低，而血液的粘度实际上会趋于一个恒定的平台值。这并不意味着模型是“错误”的；它意味着这是一个理想化的模型，捕捉了流体丰富行为的一部分，但不是全部。它提醒我们，科学通过创建模型、测试其局限性，然后构建更精细的模型来进步。

该模型在其他生物学背景中找到了更直接的应用。细菌产生的“黏液”，即生物膜，是一个被包裹在细胞外聚合物（EPS）保护性基质中的微生物城市。这种 EPS 基质是一种屈服应力流体。这赋予了生物膜类似固体的完整性，使其极难被冲走，无论是在船体上、河石上，还是在你的牙齿上（如牙菌斑）。相反，我们呼吸道内壁的粘液和润滑关节的粘液，则将其 Herschel-Bulkley 特性用于不同目的。其屈服应力使其能形成一个稳定的保护层。然而，当受到剪切时——通过关节的滑动或肺部纤毛的摆动——它能轻松流动，充当出色的润滑剂。在这种情况下，屈服应力提供了一种形式的静摩擦力，或称“粘滞性”（stiction），必须克服它才能启动滑动，这是生物摩擦学领域的一个关键原理。

我们的旅程终点是技术前沿，在这里我们对这些材料的理解不仅让我们能够描述世界，还能以新的方式建造世界。3D 生物打印领域旨在逐层构建活体组织、支架乃至器官。此过程中的“墨水”通常是负载细胞的水凝胶。

生物打印的成功取决于精妙的流变学平衡。生物墨水在打印喷嘴的高剪切作用下必须具有足够低的表观粘度，以顺畅流动，同时又不需要高到会损伤或杀死其中活细胞的压力。然而，在被挤出的瞬间，它必须具有足够高的屈服应力来固化并保持其形状，从而能够构建复杂的、悬垂的 3D 结构。目标是找到一组完美的 Herschel-Bulkley 参数 $(\tau_y, K, n)$，既能在打印过程中最大限度地减少对细胞的应力，又能最大限度地提高最终构建物的形状保真度。在这里，该模型不仅仅是一个描述性工具；它已成为制造生命本身的设计方程。

从平凡到宏伟，从地质到生物，Herschel-Bulkley 模型提供了一种统一的语言。它向我们展示了一个单一、优雅的物理概念如何能够解释种类繁多的材料的行为，揭示了我们世界复杂性背后深刻而美丽的联系。