
我们如何用数学语言来描述水或蜂蜜等流体那种熟悉的“黏性”?固体抵抗的是形变,而流体抵抗的是其形变的速率。牛顿本构关系为此提供了一个既优美简洁又极其强大的答案,构成了现代流体动力学的基石。它在流体内部的力(应力)与流体的运动和形态变化方式(应变率)之间建立了一种线性联系。本文旨在探讨如何以一种严谨而普适的方式,为这种内摩擦力,即黏度,建立模型。读者将深入理解这一定律背后的物理原理及其惊人而广泛的影响。我们将首先探讨定义该关系的核心“原理与机制”,解构应力、应变以及模型背后的物理逻辑。随后,文章将通过考察其“应用与跨学科联系”,展示该定律的巨大效用,揭示这一个方程如何帮助我们描述世界,从地球深处的岩浆流动到生物体的发育过程。
让我们从一个简单的观察开始。拿一根橡皮筋,拉伸它,然后保持住。它会回弹,储存你施加给它的能量。它记得自己最初的形状,如果你松手,它会弹回去。现在,将你的手指伸进一罐蜂蜜中搅拌。蜂蜜会抵抗你的运动;它感觉很稠,很“黏”。但你一停止搅拌,蜂蜜也停止了流动。它完全不会试图撤销你造成的漩涡。它已经彻底忘记了自己最初的平静状态。
这个简单的对比揭示了理想固体和理想流体之间的根本区别。固体抵抗持续的形变(应变),像弹簧一样储存能量。而流体则抵抗形变率(应变率)。
想象一种材料被夹在两个大的平板之间。如果我们将上方的平板滑动一小段距离并保持不动,会发生什么?如果材料是固体,比如一块橡胶,它会变形并对上方的平板持续施加一个稳定的反作用力。它能承受应力。但如果材料是流体,比如水或空气,它会流动并重新排列。一旦上方的平板停止移动,流体也随之静止,维持平板在新位置所需的力便降为零。静止的流体无法承受剪切应力。它的决定性特征就是能够流动并忘记其形状。它所提供的阻力仅在其运动时才存在。这种阻力便是黏度的本质。
如果流体只在运动时才产生阻力,我们该如何描述这种阻力呢?让我们回到那两个平板,但这次,我们让上方的平板以一个恒定的慢速移动。为此,我们必须持续对它施加一个力,以抵消来自流体的拖拽力。这仿佛是一层层的流体相互黏附,产生了一种内摩擦。这种“黏性”就是我们所说的黏度。
对于许多常见的流体——空气、水、油、蜂蜜——艾萨克·牛顿发现了一个异常简洁的定律。他提出,使一层流体滑过另一层所需的单位面积上的力(剪切应力,记作 )与其被剪切的速度(剪切率)成正比。在我们的平板例子中,剪切率就是上方平板的速度 除以两板间的距离 。于是,我们有:
这个比例常数 就是流体的动力黏度。它是衡量流体对这种剪切运动内在阻力的度量。蜂蜜的 很高;空气的则非常低。这种简单的线性关系是牛顿流体的决定性特征。
我们用来维持平板运动所消耗的能量去哪儿了?它不像在弹簧中那样被储存起来。相反,它被持续地转化为热量,使流体升温。克服黏性拖拽力所需的功率直接衡量了这种能量耗散。因此,黏度本质上是一种耗散现象;它是一个将有序运动转变为我们称之为热的无序、随机分子运动的过程。
当然,现实世界中的流体流动很少像两平板之间那样简单。在河流中或飞机机翼周围的流体单元,同时在各个方向上被拉伸、挤压和剪切。为了描述这种复杂的状态,我们需要一个比单一力值更精密的工具:柯西应力张量 。
想象一个无穷小的流体立方体。该点的“应力”状态是对作用在该微小立方体所有面上的力的完整描述。这个不仅要描述力的大小和方向,还要描述力作用的表面朝向的物理量,是一个张量。
这个张量的美妙之处在于,我们可以将其分解为具有明确物理意义的几个部分。首先,我们可以求出作用在每个面上的推力的平均值。这部分应力在所有方向上均等作用,试图改变流体单元的体积而非其形状。我们称之为静水压力 。这就像你潜入泳池时耳朵感受到的压力。在数学上,它由应力张量的迹(对角元素之和)定义:。
从总应力中减去压力后剩下的部分称为偏应力张量 。这是应力中改变形状的部分。它代表了由流体运动引起的剪切和拉伸力。因此,流体中任意点的总应力总可以写成这两部分之和:
其中 是单位张量。正是这个偏应力 与黏度直接相关。
那么,黏性应力 是由流体运动引起的。但究竟是运动的哪个方面导致了它?关键不在于速度本身,而在于速度如何从一个点到其邻近点变化。这个信息被包含在速度梯度张量 中。
一个绝妙的数学结论揭示,任何由 描述的局部流体运动都可以分解为两个截然不同的部分。一部分是流体单元的纯旋转,仿佛它是一个微小的、刚性旋转的陀螺。这由速度梯度的反对称部分描述,与涡量 有关。另一部分是纯形变——一种改变单元形状的拉伸和剪切运动。这由速度梯度的对称部分描述,称为形变率张量 。
一个深刻的问题随之而来:黏性应力是取决于旋转、形变,还是两者兼有?由物理学基本原理决定的答案是,它只能取决于形变 。这背后有两个漂亮的理由。
首先,考虑能量。黏度导致耗散。机械能转化为热能的速率由黏性应力做功给出。一点张量代数运算表明,一个对称张量(黏性应力,为了保证角动量守恒它必须是对称的)对一个反对称张量(运动的旋转部分)做功为零。因此,纯旋转不耗散能量。只有运动的形变部分 才能与黏性应力相互作用产生热量。
其次,考虑不同观察者眼中的景象。物理定律对于任何人来说都必须是相同的,无论他们是静止站立还是在旋转木马上旋转。这就是物质标架无关性原理。现在,想象一杯以恒定速率旋转的咖啡。咖啡作为一个刚体在旋转。它有大量的涡量,但其各部分之间没有相对运动。它没有在形变,因此不应承受任何黏性应力。对于这种刚体旋转,形变率张量 为零,但涡量不为零。如果应力取决于涡量,它将不为零,这在物理上讲不通。因此,黏性应力必须只取决于 ,这个量正确地识别出没有形变正在发生。
现在我们可以组装我们的杰作了。我们正在寻找一个本构关系——一个连接原因(形变率 )和结果(黏性应力 )的方程。牛顿流体的定义正是指这种关系是能想象到的最简单的形式:线性关系。
对于不可压缩流体(其体积不能改变的流体),黏性应力与形变率张量成正比。将此与我们之前对总应力的分解相结合,我们便得到了不可压缩牛顿流体的完整本构关系:
这个异常简洁的方程是牛顿原始发现的三维推广。压力 是强制实现不可压缩性的那部分应力,而黏性项 描述了由流体运动和形变引起并抵抗这种运动和形变的应力。这一个方程,与牛顿第二运动定律()相结合,便得到了著名的纳维-斯托克斯方程,它们构成了现代流体动力学的基石。
到目前为止,我们的故事都假设流体是不可压缩的,这对于像水这样的液体在许多情况下是一个非常好的近似。但对于像空气这样容易被挤压的气体呢?当流体的体积可以改变时,我们必须考虑一种新型的形变:均匀膨胀或压缩,其速率由 给出。
这引入了第二种黏性阻力的可能性:对体积变化率的阻力。这种效应由第二个黏度系数控制,称为体积黏度 (它通过 与另一个系数 相关)。它代表了一种内摩擦,只在流体被快速压缩或膨胀时才出现,例如在声波或激波中。
在许多应用中,通常会做出斯托克斯假设,即假设这个体积黏度为零()。对于像氦这样的简单单原子气体,分子运动论证实这是一个极好的近似。其物理原因在于,在这种气体中,能量只储存在原子的平动中,而压缩产生的能量几乎瞬间就在它们之间分配完毕。
然而,斯托克斯假设只是一个近似,它的失效揭示了更深层次的物理学。在像空气(氮气和氧气)这样的多原子气体中,来自压缩的能量不仅要分配到平动中,还要分配到分子的转动和振动中。这个过程需要时间。如果你压缩气体的速度快于这个内部的“弛豫时间”,就会发生一种热力学摩擦,表现为非零的体积黏度。这就是为什么体积黏度对于理解空气中高频声的吸收以及激波的结构至关重要,。同样,对于许多液体和复杂流体,斯托克斯假设也会失效,其体积黏度可能相当大。这一细微差别显示了我们简单模型的局限性,并指向了物质丰富的内禀动力学。
我们将黏度视为一种宏观性质,但其起源在于分子的混乱微观世界。黏度不是单个水分子的属性;它是数万亿个水分子协同作用的涌现属性。
想象一下流体分层流动,每一层都比下面一层移动得稍快一些。流体中的分子并不局限于各自的层内;它们处于持续的、随机的热运动中,不停地晃动和碰撞。一个来自较高、较快层面的快速移动分子偶尔会晃动到下面较慢的层面。在随后的碰撞中,它将其多余的动量传递给它新的、较慢的邻居,给它们一个轻微的向前推动。对称地,一个来自较低层面的较慢分子可能会晃动到较快的层面,并通过碰撞起到拖拽作用,减慢那一层的速度。
相邻流体层之间这种持续的微观动量交换,是黏性剪切应力的真正物理起源。黏度系数 无非是衡量这种分子动量传递效率的尺度。气体的分子运动论表明,黏度随温度升高而增加,且基本上与压力无关,这是一个令人惊讶的结果,源于分子密度与分子在两次碰撞间平均行进距离之间的相互作用。这是连接单个粒子混乱之舞与连续介质世界平滑、可预测流动的惊人桥梁,让我们能够理解从蜂蜜的流动到木星大气壮丽漩涡的一切现象。
在探索了牛顿本构关系背后的原理与机制之后,你可能会有一种感觉,认为它虽然简洁,但或许有些枯燥,缺乏现实意义。这种应力与应变率之间的简单线性关系,真的能对我们所生活的复杂、混乱的世界说些什么吗?答案是肯定的,而且其应用的故事是一场跨越学科的旅程,揭示了自然设计中惊人的一致性。支配管道中水流动的基本定律,同样也塑造着生物胚胎的心脏,决定着地球深处熔岩的行为,甚至为驯服湍流的混沌提供了一种巧妙的技巧。
让我们从一个熟悉的场景开始:一根简单的管道。每当我们泵送流体时,无论是家中的自来水还是工业生物反应器中的培养液,我们都在与一股无形的阻力作斗争。这种力,这种拖拽力,其根源就在于牛顿本构关系的核心。当流体流动时,与管壁接触的那一层是静止的——这便是著名的无滑移条件。离管壁稍远一点,流体缓慢移动,再远一点,移动得更快。这种速度随位置的变化,即速度梯度 ,就是剪切应变率。
根据我们的定律,这个应变率会产生一个剪切应力 。这是流体试图拖动静止的管壁,反之,也是管壁试图阻碍流体。这种作用于管道整个内表面的微观应力,是所有黏性阻力的来源。当你将这个应力沿管道全长进行累加,就得到了水泵必须克服的总拖拽力。
但这场斗争是有代价的,这是由热力学定律决定的后果。必须做功来推动流体抵抗黏性摩擦,而这些能量并不会凭空消失。它被耗散掉,转化为分子的随机运动——换句话说,就是热量。决定应力的同一个项,即速度梯度,也告诉我们有多少能量正在损失。单位体积的能量耗散率与黏度和速度梯度的平方成正比,即 。所以,流体越黏,或者你试图剪切它的速度越快,你给它加热的速度就越快。每当你搅拌蜂蜜时,你都在将手臂的机械能转化为热能,使蜂蜜的温度有微乎其微的升高,这正是其牛顿性质的直接体现。
基础物理定律的美妙之处在于它们对具体情境的全然不顾。它们同样适用于熔岩和水。模拟火山通道中岩浆流动的地球物理学家,与设计工厂管道的化学工程师,使用的是完全相同的本构关系。岩浆,一种黏度极高的硅酸盐液体,向上蠕动,其速度剖面由压力和黏性力的相同平衡所塑造。它对通道壁施加的剪切应力,由 计算得出,可能巨大无比,足以使周围的岩石破裂,并影响火山喷发的动态。唯一的区别是数值的尺度——岩浆的黏度 比水高出数十亿倍,但原理是相同的。
故事变得更加引人注目。固体可以是牛顿流体吗?在人类的时间尺度上,我们认为玻璃是典型的固体。但在地质时间尺度上,或在制造时的高温下,它会流动。非晶固体中的原子缺乏固定的晶格,可以在应力下缓慢地重新排列。这种缓慢的蠕变流动通常可以被描述为一种高黏度的牛顿流体。
考虑烧结过程,即细粉末被加热形成固体物件。两个刚刚接触的微小非晶球体并未处于其最低能量状态。它们巨大的表面积代表了多余的能量,我们称之为表面张力。这种表面张力拉动颗粒,产生一种应力,试图通过将它们融合在一起来最小化表面积。这种应力驱动了缓慢的黏性流动,导致颗粒之间生长出一个物质“颈”。令人惊讶的是,我们可以通过材料的黏度,将来自表面张力的驱动应力与应变率(颈部几何形状的变化)联系起来,从而对此进行建模。牛顿本构关系使我们能够预测粉末融合的速度,这是陶瓷、冶金以及某些材料3D打印的基础过程。这里的“流体”是固体,而“泵”则是表面张力本身。
牛顿关系的微妙力量在生命领域表现得最为淋漓尽致。流动流体施加的力不仅仅是要克服的障碍;它们是引导生物体发育和功能的关键信号。
让我们进入一个发育中的斑马鱼胚胎的心脏,这个微小生物的透明身体让我们得以观察生命的展开。在其初生的循环系统中,血液——在这个尺度下表现为一种简单的牛顿流体——流过发育中的血管。当它流动时,它对构成血管壁的内皮细胞层施加剪切应力。这个应力的大小,也许只有区区 帕斯卡,看似微不足道。但对细胞而言,这却是一个响亮的信号。
细胞表面的特殊蛋白质作为力学传感器,探测到这种物理拖拽。细胞“感觉”到了流动。这个物理信号随后被转化——即转导——为细胞内部的生化级联反应,最终改变了哪些基因被开启或关闭。例如,一个名为Klf2a的基因的表达对剪切应力高度敏感。这个基因反过来又调节Notch信号通路,一个发育的主控制器。结果呢?血流的物理力量告诉细胞是该成为动脉还是静脉,并指导它们如何将自己塑造成心脏瓣膜的精巧叶片。如果血流太弱(可能由于血液黏度较低)或不存在,信号就会丢失,心脏便无法正确形成。这就是力学生物学:在非常真实的意义上,牛顿本构关系帮助书写了生命之书。
同样的定律也可以描述微观尺度上的战斗。考虑生物膜,一种可能导致感染或堵塞管道的黏滑细菌菌落。这个菌落由胞外聚合物基质聚合在一起,这赋予了它一定的内聚强度。当水流过它时,会施加一个剪切应力。在一段时间内,生物膜保持稳固。但如果水流变得足够强,施加的剪切应力最终会超过生物膜的内聚强度,大块的生物膜便会在一个称为脱落的过程中被撕裂。关键在于,真实的生物膜并非一个完美光滑的表面。它有凸起和塔状结构。这些粗糙部分伸入流中,导致流体在它们周围加速。这导致了局部剪切应力高得多的“热点”,这种应力集中现象使得即使在平均应力较低的情况下,水流也能找到一个薄弱点并引发破坏。
在21世纪,我们许多最伟大的科学和工程壮举都发生在计算机内部。牛顿本构关系是计算建模这座宏伟大厦中不可或缺的基石。
想象一下像酒杯内壁形成的“酒泪”这样精巧的现象。这种被称为马兰戈尼效应的现象是由表面张力的梯度驱动的。在气液界面,酒精蒸发得更快,增加了局部的含水量,从而提高了表面张力。这种表面张力梯度将液体沿杯壁向上拉。被拉上去的液体必须与来自下方流体的黏性剪切应力相平衡。模拟的边界条件直接陈述了这种平衡:表面的剪切应力 必须等于表面张力梯度 。没有本构关系,我们就无法模拟这种美丽的效果。
当我们模拟更复杂的场景时,比如风与桥梁的相互作用或血液流经人造心脏瓣膜,我们正在解决一个流固耦合(FSI)问题。计算机需要知道流体如何推动固体。这由动态界面条件决定,该条件是牛顿第三定律的数学表达:流体对固体施加的力矢量(面力)必须等于固体对流体施加的力矢量。流体的应力张量 包含了我们的牛顿关系,正是它使我们能够计算流体施加的力,从而使这些拯救生命的模拟成为可能。
最后,对于那些绝非简单和牛顿流动的流,比如湍射流或船后的尾流,该怎么办?在这里,本构关系激发了一种绝妙的智力建模。湍流中混乱的涡旋产生的应力远比层流复杂。对这种混沌进行完全模拟在计算上往往是不可能的。Boussinesq涡黏性假设提供了一个巧妙的变通方法。它建议我们将所有湍流涡旋的平均效应建模,就好像它是一个额外的、巨大的黏度,即“涡黏度” 。然后我们使用一个类似牛顿的本构关系来模拟平均流动,但使用的是这个新的、大得多的黏度。这承认了底层的流动并非牛顿流动,但它证明了原始概念的力量,即我们可以通过简单地保留其数学形式并创造一个新的有效参数来建立一个强大的预测模型。
从平凡到生命,再到虚拟世界,应力与应变率之间线性关系的简单前提,已被证明是一个具有非凡力量和广泛影响力的思想。它是一条连接不同世界的金线,是一个美丽的例子,说明一个简单的物理学片段如何能够阐明宇宙在各种尺度上的运作方式。