广义牛顿流体

尽管水的可预测流动可以被牛顿物理学优雅地描述，但我们日常生活和工业中的许多流体——从番茄酱、油漆到血液和岩浆——其行为方式要复杂得多。它们对流动的阻力，即黏度，不是一个固定的常数，而是一个随运动变化的动态属性。这种对简单定律的违背提出了一个重大挑战，揭示了经典牛顿模型在描述一大类常见材料时的局限性。

本文旨在通过介绍​​广义牛顿流体 (GNF)​​ 的概念来填补这一知识空白，这是一个强大的框架，它扩展了经典流体动力学以捕捉这些引人入胜的行为。在接下来的章节中，您将对这一重要主题获得全面的理解。首先，在“原理与机制”一节中，我们将探讨表观黏度的基本概念，深入研究剪切稀化、剪切增稠和黏塑性流体的独特行为，并考察用于描述它们的数学模型。随后，“应用与跨学科联系”一节将展示这些原理对于理解从我们动脉中的血液流动到工业泵的设计，再到地球深处岩浆的运动等各种事物为何至关重要。

想象一下搅动一杯水。你搅得越快，感受到的阻力就越大。将勺子的速度加倍，移动它所需的力量也会加倍。这种优美简洁的线性关系，是物理学家所称的​​牛顿流体​​的标志，该名称源于 Isaac Newton，他对力学的深刻见解为描述这种流体奠定了基础。对于这些流体——如水、空气或汽油——其内部流动阻力，即我们称之为​​黏度​​的属性，是一个常数。在给定温度下，一个单一的数值 $\mu$ 就能告诉你关于其“稠度”的一切信息。剪切应力 $\tau$ 是你为使流体流动而施加的单位面积力，它与剪切速率 $\dot{\gamma}$ 成正比，剪切速率衡量流体变形的速度：$\tau = \mu \dot{\gamma}$。

几个世纪以来，这幅优雅的图景已经足够。但自然界在其无限的创造力中，很少如此简单。我们的日常生活中充满了违背这一定律的物质。想想番茄酱。它顽固地待在瓶子里，像一团浓稠的、近乎固体的物质。即使你把瓶子倒过来，它也可能一动不动。但只要用力摇晃或在瓶底拍一下，它就会突然涌出，几乎像水一样流动。一旦流动停止，它又会变稠。显然，它的黏度不是一个常数；它随运动而改变。现在考虑一种玉米淀粉和水的浓稠混合物，这是科学演示中的常客，有时被称为“欧不裂 (oobleck)”。你可以像在液体中一样慢慢地让手指穿过它，但如果你试图快速捶打或搅拌它，它会瞬间变硬，像固体一样抵抗。

这些就是​​非牛顿流体​​，它们无处不在：油漆、血液、液体肥皂、牙膏、酸奶以及工业中使用的许多聚合物溶液。它们的行为远比其牛顿流体“表亲”更丰富、更有趣。要理解它们，我们必须摒弃黏度恒定的观念，进入一个流体“稠度”是动态属性的世界，它是对其被移动方式的一种响应。

我们需要的核心概念是​​表观黏度​​，通常用 $\eta$ 表示。这是一个直接的概念：在任何时刻，我们都可以将流体的黏度定义为所施加的剪切应力与产生的剪切速率之比，即 $\eta = \tau / \dot{\gamma}$。对于牛顿流体，这个比值总是一个相同的常数 $\mu$。对于非牛顿流体，这个比值会随着我们改变剪切速率而变化。表观黏度不是一种基本的材料常数，而是流动本身的函数，即 $\eta(\dot{\gamma})$。这个看似微小的修改，开启了一个广阔而迷人的流体行为领域。

大多数非牛顿流体可分为几个关键类别：

• ​​剪切稀化（假塑性）：​​ 这类流体的表观黏度随着剪切速率的增加而降低。你搅得越快，它们就变得越“稀”。油漆就是一个绝佳的例子；它在罐中很稠，所以不会从刷子上滴下来，但当你用快速的笔触（高剪切速率）涂抹它时，它会变稀并容易铺开。血液是另一个关键例子。在极低的流速下，在你身体的微小毛细血管中，红细胞倾向于聚集在一起，形成称为​​缗钱状串列 (rouleaux)​​ 的堆叠。这些团块产生很大的阻力，导致表观黏度很高。当血液流动加快时，流体的水动力会打散这些缗钱状串列，单个细胞会变形并与流动方向对齐。这种细胞层面的重组显著降低了流动阻力，导致表观黏度下降。

• ​​剪切增稠（胀流性）：​​ 这类流体则相反；它们的表观黏度随着剪切速率的增加而增加。玉米淀粉和水的浆料是典型例子。它由悬浮在液体中的精细固体颗粒组成。在静止或缓慢变形时，水润滑着颗粒，使它们能够相互滑过。当你施加一个突然的高剪切速率（如拳击）时，颗粒没有时间移开。它们会堵塞在一起，形成一个临时的、类似固体的结构来抵抗运动。这种黏度的迅速增加使得流体感觉像固体。

• ​​黏塑性（屈服应力）：​​ 当施加的应力较低时，这些材料表现得像固体，只有当应力超过一个称为​​屈服应力​​ $\tau_y$ 的临界值时，它们才开始流动。牙膏就是一个完美的例子。它能在牙刷上保持形状（处于“未屈服”状态），因为重力产生的小应力小于其屈服应力。只有当你用力挤压牙膏管，施加大于 $\tau_y$ 的应力时，它才会像流体一样流动（进入“已屈服”状态）。许多黏塑性流体，如番茄酱，一旦开始流动，也具有剪切稀化性。

我们如何在一个数学框架中捕捉这些五花八门的行为？我们需要一个看起来像牛顿定律，但允许黏度变化的定律。这就是​​广义牛顿流体 (GNF)​​ 模型的核心。

在其完整、华丽的三维形式中，GNF 模型指出，偏应力张量 $\boldsymbol{\tau}$（代表流体中的剪切力、非压力力）与变形速率张量 $\boldsymbol{D}$（描述流体如何被拉伸和剪切）成正比： $\boldsymbol{\tau} = 2\eta(\dot{\gamma})\boldsymbol{D}$ 这个方程是物理建模的杰作。它看起来与牛顿定律几乎完全相同，保留了其优美的结构。然而，黏度 $\eta$ 不再是一个常数，而是局部剪切速率 $\dot{\gamma}$ 的函数。该剪切速率的定义方式与观察者的参考系无关，它使用了变形张量的第二不变量来定义，即 $\dot{\gamma} = \sqrt{2\boldsymbol{D}:\boldsymbol{D}}$。这种表述方式巧妙地确保了模型遵守基本的物理原理，如物质坐标系无关性（物理规律不应仅仅因为你在一个旋转平台上观察而改变）和热力学一致性（它确保流动的流体总是耗散能量，绝不会自发产生能量）。

GNF 模型的威力来自于为函数 $\eta(\dot{\gamma})$ 选择不同的数学形式：

• ​​幂律模型：​​ 一个简单但非常有效的选择是 $\eta(\dot{\gamma}) = K\dot{\gamma}^{n-1}$，其中 $K$ 是“稠度系数”，$n$ 是“流动行为指数”。如果 $n=1$，剪切速率项消失，我们便得到了黏度为 $\mu=K$ 的牛顿流体。如果 $n 1$，黏度随着 $\dot{\gamma}$ 的增加而减小，完美地描述了剪切稀化。如果 $n>1$，黏度随着 $\dot{\gamma}$ 的增加而增加，捕捉了剪切增稠的特性。该模型揭示了一个关键的行为差异：在简单剪切流中，牛顿流体的应力与速度梯度 ($U/H$) 成正比，而幂律流体的应力则与 $(U/H)^n$ 成正比，这是其非线性性质的直接结果。

• ​​黏塑性模型：​​ 为了描述具有屈服应力的流体，我们可以使用 ​​Herschel-Bulkley 模型​​，该模型将屈服应力与幂律行为结合起来：要发生流动，应力必须超过屈服应力，即 $\tau > \tau_y$，此时关系式变为 $\tau = \tau_y + K\dot{\gamma}^n$,。这个优雅的方程解释了牙膏的初始阻力及其随后的流动。如果 $n=1$，该模型简化为 ​​Bingham 模型​​。

• ​​更真实的模型：​​ 幂律模型尽管实用，但存在一些不符合物理实际的怪异之处。对于剪切稀化流体 ($n 1$)，它预测在静止时（$\dot{\gamma} \to 0$）黏度为无穷大，在无限剪切速率下黏度为零。真实流体并非如此。相反，它们的黏度通常在一个高但有限的​​零剪切黏度​​ ($\eta_0$) 和一个低但非零的​​无穷剪切黏度​​ ($\eta_\infty$) 处趋于平稳。为了捕捉这一特性，人们开发了更复杂的模型，如 ​​Carreau 模型​​或 ​​Carreau-Yasuda 模型​​。这些模型使用更复杂的数学形式，能够在低剪切速率下的 $\eta_0$ 平台和高剪切速率下的 $\eta_\infty$ 平台之间平滑过渡，从而为许多真实流体（从聚合物溶液到我们关节中的滑液等生物流体）提供了更精确的描述,。

广义牛顿流体模型是一个强大而直观的工具。它为思考一大类材料提供了一个框架。但就像任何优秀的科学模型一样，其力量也存在于其局限性之中，即它告诉我们它无法解释什么。

GNF 模型是纯黏性的。这意味着流体在任何给定时刻的应力仅取决于该确切时刻的变形速率。它没有对其过去的“记忆”。如果你停止剪切一个 GNF，其内应力会立即消失。然而，许多真实的非牛顿流体，如聚合物熔体或面包面团，是​​黏弹性​​的。它们同时表现出黏性（类液体）和弹性（类固体）的特性。如果你使这种流体变形然后突然停止，应力不会立即消失；它需要时间来松弛回零。这是因为其底层的微观结构，如长聚合物链，被拉伸并储存了能量，就像微小的橡皮筋一样。它们需要时间才能回缩到原始状态。根据其定义，GNF 模型无法捕捉这种流体记忆。

此外，GNF 模型预测应力总是与变形方向对齐（$\boldsymbol{\tau}$ 平行于 $\boldsymbol{D}$）。这带来一个至关重要的推论：它预测在简单剪切流中，不会产生垂直于流动方向的额外力。实验证明，对于许多聚合物流体而言，这一推论是错误的。如果你将一根旋转杆放入一杯黏弹性流体中，流体会神秘地沿杆向上爬升——这种效应被称为​​魏森贝格效应 (Weissenberg effect)​​。这种爬升是由称为​​法向应力差​​的力驱动的，这些力作用于与主剪切方向垂直的方向上。GNF 模型由于其自身结构，预测这些法向应力差始终为零，因此无法解释此类引人入胜的现象,。

这些局限性并非 GNF 模型的失败。恰恰相反，它们是模型的成功之处。它们精确地告诉我们该模型适用于何处，以及在何处必须由更丰富的理论——黏弹性理论——来接替。从简单的牛顿流体到广义牛顿流体的历程，是科学过程的完美例证：观察世界的丰富性，扩展我们的理论以涵盖它，并在此过程中，发现指向更深层次理解的边界。

在我们迄今的探索中，我们揭示了广义牛顿流体的秘密：其黏度不是一个固定不变、神赐的数值，而是一个响应流动本身的动态属性。这看似只是对我们简洁的牛顿世界的一个细微调整，但正是这种微妙之处，既为自然所钟爱，也为工程师所必须掌握。当黏度可以随时随地变化时，流体运动的整个特性都会被改变。流动方程变得更加丰富，现象更加出人意料，而看似遥远的科学领域之间的联系也变得更加深刻。

现在，让我们来探索这个响应性流体的世界。我们将看到，这单一的原理——一个随剪切节奏起舞的黏度——如何主导我们身体内部安静高效的运动，如何在工业过程中提出严峻的挑战，并如何重塑我们对热学、地质学乃至湍流混沌漩涡的理解。

非牛顿流最亲密、最优雅的应用或许就在我们身体内部。我们的身体并非充满了简单的水，而是充满了复杂的、有生命的流体，这些流体经过进化，能够执行非凡的任务。

思考一下滑膜关节的奇迹，比如你的膝盖或肩膀。关节表面覆盖着软骨，并由一层薄薄的滑液隔开。这种液体是一种卓越的“智能润滑剂”。当你静止或缓慢移动时，滑液浓稠且黏滞，起到坚固缓冲垫的作用。但当你需要快速移动时——挥动球棒、踢球——滑液会神奇地变稀。随着剪切速率的增加，黏度下降，从而实现了快速、低阻力的运动。原因何在？滑液中的透明质酸 (hyaluronan) 等长链聚合物会与流动方向对齐，解开缠结，从而更容易地相互滑过。这种剪切稀化行为意味着，黏性阻力扭矩 $M$ 不会像简单油液那样随关节角速度 $\omega$ 线性增加。相反，它以次线性方式增长（$M \propto \omega^n$，其中 $n 1$）。这是自然工程的奇迹，一种能够自我调节的润滑剂，总能完美地适应手头的任务。

然后是我们生命的河流：血液。血液不是一种均匀的红色液体，而是红细胞、白细胞和血小板在血浆中的密集悬浮液。这种“颗粒性”使其成为一种典型的剪切稀化流体。在低剪切速率下，在微小的毛细血管中，红细胞会聚集在一起形成称为缗钱状串列 (rouleaux) 的堆叠，从而增加黏度。在高流速的大动脉中，这些堆叠会散开，细胞与流动方向对齐，从而降低黏度。

这对于我们如何描述流动本身具有深远的影响。对于像水这样的简单牛顿流体，动量方程中的黏性力项可以很好地简化。但对于血液，黏度 $\mu(\dot{\gamma})$ 依赖于剪切速率，而剪切速率本身又随位置变化。这意味着黏度变成了一个场，一个逐点变化的量，它再也不能在运动控制方程中被从导数中提出来。问题的基本数学结构被改变了。

这不仅仅是理论家头疼的问题；它具有实际的生物学重要性。流动的血液施加在动脉壁上的力——即壁面剪切应力——是排列在我们血管内壁的内皮细胞的一个关键信号。这些细胞对该应力作出反应，控制动脉的直径并标记炎症。如果我们错误地将血液建模为具有某个平均黏度的简单牛顿流体，我们对这个至关重要的壁面剪切应力的计算可能会大错特错。将血液精确地建模为广义牛顿流体，通常是通过将幂律模型拟合到实验数据来实现的，这对于理解血管健康与疾病至关重要。

走出人体，进入工厂，我们发现工程师们几个世纪以来一直在与非牛顿流体作斗争，无论他们是否意识到这一点。油漆、番茄酱、牙膏、钻井泥浆和纸浆浆料都是广义牛顿流体。将它们当作水来处理是灾难的根源。

以通过管道泵送流体这一简单任务为例。对于水或油，工程师有一个可靠的工具：Moody 图，它根据管道的粗糙度和雷诺数给出摩擦系数。但如果你试图对纸浆浆料使用此图，它会给你错误的答案。该图建立在一个不可动摇的假设之上：黏度是恒定的牛顿黏度。对于浆料，“黏度”是流动本身的函数。该图的根基就是无效的。

其后果更为深远。对于管道中的牛顿流，驱动流动的压降 $\Delta p$ 与产生的体积流量 $Q$ 之间的关系是优美的线性关系——这就是著名的 Hagen-Poiseuille 定律。对于剪切稀化流体，情况不再如此。由于流体流速越快变得越稀，给定的压力增加会产生大于比例的流量增加。这种关系变为非线性，通常遵循幂律 $Q \propto (\Delta p)^{1/n}$，其中 $n 1$。此外，经典的抛物线速度剖面不复存在，取而代之的是一个更平钝、更像柱塞流的剖面。现代流体力学认识到这一点，“Poiseuille 流”一词也已演变为描述充分发展的管流的这种一般运动学状态，而不管压力-流量关系是否为线性。

在高科技制造业中，复杂性与日俱增。用于涂覆锂离子电池电极的浆料是流变学的奇迹。许多浆料不仅具有剪切稀化性，还具有​​屈服应力​​；在施加最小应力之前，它们表现得像固体，之后则像液体一样流动。这种特性可能是一种福音，有助于防止涂层膜破裂成细流。其他浆料具有​​触变性​​：其黏度取决于剪切历史，随时间分解，在静止时恢复。有些甚至表现出弹性的“记忆”，在变形后会轻微回弹——这种特性称为​​黏弹性​​。虽然这些行为超出了简单的广义牛顿模型，但它们遵循相同的核心原则：应力与应变率之间的关系不是简单的线性关系。

即使是单个颗粒在这种流体中的运动也是一场微妙的舞蹈。要计算悬浮液中微小球体所受的曳力，人们可能会倾向于使用斯托克斯定律 (Stokes' law)。但应该使用哪种黏度呢？流体的黏度取决于剪切速率，而剪切速率是背景流剪切和颗粒相对于流体自身运动所产生的额外剪切的组合。一个恰当的模型必须考虑到这个反馈回路，其中颗粒的运动会局部稀化其周围的流体，从而反过来改变其所受的曳力。

非牛顿特性的影响像涟漪一样扩散开来，触及了几乎所有存在流体的物理科学角落。

热传递又如何呢？如果我们改变管道中的速度剖面，我们也必须改变热量被流动携带或对流输运的方式。对于剪切稀化流体，与相似条件下的牛顿流体相比，其更平钝的速度剖面可以导致更高效的混合和增强的热传递。此外，黏性力所做的功会转化为热量——这一过程称为黏性耗散。能量方程中的这个热源项本身就依赖于非牛顿黏度。由于这些原因，为牛顿流体开发的经典热传递关联式完全失效。人们必须回归第一性原理，认识到动量方程和能量方程不仅通过速度耦合，还通过作为解本身一部分的黏度耦合。

让我们将视野放大到广阔的地质学尺度。石油、岩浆或用于环境修复的聚合物溶液是如何流过岩石和土壤中错综复杂的孔隙网络的？对于水，我们有达西定律 (Darcy's Law)，这是一个简单的线性关系，指出流量与压力梯度成正比。但如果流体是广义牛顿流体，比如钻井泥浆，情况又如何？在微观孔隙尺度上，流体的黏度随局部剪切速率而变化。当我们在大体积上平均这种行为时，宏观的达西定律继承了这种非线性。多孔介质的有效阻力不再是一个常数，而是取决于流速。其结果是广义的、幂律形式的达西定律，其中微观流变学直接决定了宏观输运定律。

最后，我们来到了整个经典物理学中尚存的一大挑战：湍流。对于牛顿流体，用于模拟湍流的雷诺平均纳维-斯托克斯 (RANS) 方程已经包含了未封闭的雷诺应力项，该项代表了混沌涡流引起的动量输运。对于非牛顿流体，出现了一个新的难题。由于黏性应力定律是非线性的，对 方程进行时间平均的过程会产生与黏度自身脉动相关的第二个未封闭项。流体的表观黏度现在与湍流脉动耦合，形成了一个令人眼花缭乱的反馈循环。模拟非牛顿湍流是现代研究的前沿领域，证明了当我们放弃恒定黏度这一简单舒适区时，会产生多么深刻而复杂的挑战。

从我们自身的关节到地球的核心，再到湍流的混沌，广义牛顿流体的原理揭示了一个充满惊人复杂性和相互联系的世界。它告诉我们，在自然界中，材料的属性与其运动的动力学常常是同一首优美交响乐中不可分割的部分。