try ai
科普
编辑
分享
反馈
  • 牛顿流体的本构关系

牛顿流体的本构关系

SciencePedia玻尔百科
核心要点
  • 牛顿流体的定义是一个简单的线性关系,即粘性应力与应变率成正比,比例常数即为流体的粘度。
  • 应力张量和应变率张量是描述连续流体介质内部力和变形的基本数学工具。
  • 克服粘性力所做的功通过一个称为粘性耗散的过程不可逆地转化为热能,这是热力学第二定律的直接体现。
  • 这一基础模型准确地描述了许多常见流体的行为,是工程学、地球物理学、生物学和医学等多个领域的关键原理。

引言

搅动液体时感受到的阻力——即其“稠厚”或“粘滞”感——是一种我们都能理解的现象。这个被称为粘度的特性,是真实世界流体的决定性特征。然而,要从这种直观概念转变为一个精确的、可预测的科学定律,需要更严谨的语言。我们如何用数学方式来描述控制空气流过机翼、血液流过动脉或地球深处岩浆流动的内摩擦力?答案在于由 Isaac Newton 首次提出的一个流体动力学基础模型。

本文旨在弥合流体摩擦的定性体验与其定量描述之间的鸿沟。文章详细阐述了牛顿流体的本构关系,这条优美的法则是我们理解大量流动类型的基石。通过探索这一原理,您将深入了解连接流体内部作用力与其运动的基本物理学。

我们将首先探讨“原理与机制”,建立应力张量和应变张量的语言,以形式化地定义表征牛顿流体的线性关系。我们还将探讨其令人惊讶的推论,从热量的产生到无滑移条件的起源。接下来,本文将展开“应用与跨学科联系”的旅程,展示这条单一、简单的定律如何为工程学、地球物理学和生物学等领域的现象提供关键见解,揭示物理学在复杂世界中的深刻统一性。

原理与机制

想象一下搅动一杯茶和一罐蜂蜜。你都能感觉到阻力,但两者差别巨大。这种对流动的阻力,即内摩擦力,就是我们所说的​​粘度​​(viscosity)。它是区分物理入门课程中理想化的、完全光滑的流体与塑造我们世界的真实流体(从流过机翼的空气到在我们血管中奔流的血液)的决定性特征。但我们如何用物理学的精确语言来捕捉这种粘滞、复杂而又优美的特性?简单的搅动行为如何导致涡流的复杂舞蹈和流体不可避免的升温?

理解这一点的过程要求我们建立一种新的语言,一种能够描述力如何传递以及流体元如何变形的语言。

力与运动的语言

当你推一个固体块时,作用力是直接的。但流体是连续介质。施加于一点的力会在整个流体中产生复杂的内力模式。要描述这一点,我们不能只谈论单个力矢量。我们需要一个更强大的概念:​​应力张量​​(stress tensor)。

想象一个位于我们运动流体中任意位置的微小虚构立方体。在其六个面上都有力作用。这些力不一定垂直于面;它们既可以是推(或拉)力,也可以是刮擦力。为了描述该点的受力状态,我们需要指明每个面上的力。

让我们使用一个坐标系 (x,y,z)(x, y, z)(x,y,z)。考虑我们立方体的一个面,其外法线指向 xxx 方向。这个面外的流体对内部的流体施加一个力。这个力有三个分量:一个在 xxx 方向(法向推/拉),一个在 yyy 方向(切向刮擦),还有一个在 zzz 方向(另一个切向刮擦)。我们将这些单位面积上的力表示为 σxx\sigma_{xx}σxx​、σxy\sigma_{xy}σxy​ 和 σxz\sigma_{xz}σxz​。第一个下标表示面的朝向(一个“朝向”xxx 的面),第二个下标表示力分量的方向。

因此,σxx\sigma_{xx}σxx​ 是​​正应力​​(normal stress),即与压力相关的那种应力。而 σxy\sigma_{xy}σxy​ 是​​剪切应力​​(shear stress):它是作用在法线为 xxx 方向的表面上,方向为 yyy 的单位面积力。这是以不同速度移动的相邻流体层之间刮擦或摩擦的数学体现。所有九个这样的分量集合构成了​​柯西应力张量​​(Cauchy stress tensor),σ\boldsymbol{\sigma}σ,它如同一个机器,能告诉你作用在你能想象的任何穿过流体的切面上的力矢量。

那么,除了无处不在的静压之外,是什么引起了这些应力呢?答案是运动——或者更准确地说,是流体不同部分之间的相对运动。像刚体一样平移或旋转的流体不会产生内摩擦力。只有当流体发生变形时,摩擦力才会出现。为了描述这种变形,我们考察​​速度梯度​​(velocity gradient)∇u\nabla \mathbf{u}∇u,这是一个张量,告诉我们速度矢量 u\mathbf{u}u 如何随点变化。

令人惊奇的是,由 ∇u\nabla \mathbf{u}∇u 描述的任何复杂的无穷小运动都可以清晰地分解为两部分:纯旋转和纯应变(变形)。旋转部分由​​自旋张量​​(spin tensor)Ω\boldsymbol{\Omega}Ω 捕捉,它与流的​​涡度​​(vorticity)(局部旋转运动)有关。变形部分由​​应变率张量​​(rate-of-strain tensor)S\boldsymbol{S}S 捕捉。流体力学的一个关键洞见是,对于大多数流体,粘性应力不依赖于局部旋转,而只依赖于应变率。流体元旋转得多快并不重要,重要的是它被拉伸或剪切了多少。

牛顿的约定:线性比例定律

所以,我们有代表内摩擦力的应力 τ\boldsymbol{\tau}τ(σ\boldsymbol{\sigma}σ 的粘性部分),和代表变形的应变率 S\boldsymbol{S}S。它们之间有什么关系呢?在17世纪晚期,Isaac Newton 提出了最简单的可能关系:它们成正比。

这就是​​牛顿流体​​(Newtonian fluid)的本质:由粘度产生的应力是应变率的线性函数。

对于​​不可压缩流体​​(incompressible fluid)——即密度不发生变化的流体,如通常条件下的水或油——这个“牛顿约定”呈现出一种优美而简单的形式:

τ=2μS\boldsymbol{\tau} = 2\mu\boldsymbol{S}τ=2μS

在这里,τ\boldsymbol{\tau}τ 是粘性应力张量,S\boldsymbol{S}S 是应变率张量。比例常数 μ\muμ 是​​动力粘度​​(dynamic viscosity),这个数值告诉我们面对的是水还是蜂蜜。因子 2 是一个惯例,它使得在简单情况下公式更为简洁。例如,在两块平板间的简单剪切流中,这个宏大的张量方程可以优雅地简化为我们熟悉的形式,将剪切应力 τ\tauτ 与剪切率 γ˙\dot{\gamma}γ˙​ 联系起来:τ=μγ˙\tau = \mu \dot{\gamma}τ=μγ˙​。

因此,不可压缩牛顿流体的完整应力张量是各向同性压力和粘性应力之和:

σ=−pI+2μS\boldsymbol{\sigma} = -p\boldsymbol{I} + 2\mu\boldsymbol{S}σ=−pI+2μS

其中 I\boldsymbol{I}I 是单位张量。这里隐藏着一个微妙而深刻的要点。对于不可压缩流体,密度是恒定的,因此压力不能再通过像理想气体定律这样的状态方程来确定。那么,压力 ppp 是由什么决定的呢?压力变成了一种“强制执行者”。它在整个流体中瞬时调整自身,达到任何需要的值,以确保流体保持不可压缩性(∇⋅u=0\nabla \cdot \mathbf{u} = 0∇⋅u=0)。用高等力学的语言来说,压力场 ppp 充当了不可压缩性约束的​​拉格朗日乘子​​(Lagrange multiplier)。它是一种力学压力,而非热力学压力。

这个简单的线性定律并非普适真理,而是一个模型。然而,它是一个惊人成功的模型,准确地描述了包括空气、水、油和许多气体在内的大量流体。对于密度可以变化的​​可压缩流体​​(compressible fluids),该定律稍微复杂一些,涉及到与体积变化阻力相关的第二粘度系数 λ\lambdaλ。完整的关系式变为 σ=−pI+2μS+λ(∇⋅u)I\boldsymbol{\sigma} = -p\boldsymbol{I} + 2\mu\boldsymbol{S} + \lambda(\nabla \cdot \mathbf{u})\boldsymbol{I}σ=−pI+2μS+λ(∇⋅u)I。请注意,如果流是不可压缩的,则 ∇⋅u=0\nabla \cdot \mathbf{u}=0∇⋅u=0,我们就回到了更简单的形式,这揭示了其背后物理学的美妙统一性。

简单定律的惊人推论

线性关系 τ=2μS\boldsymbol{\tau} = 2\mu\boldsymbol{S}τ=2μS 看起来很普通,但其推论却很丰富,且常常与直觉相悖。

首先,粘度不仅仅与剪切有关。考虑一个纯拉伸的流动,比如一个流体在 xxx 方向被拉开,在 yyy 方向被挤压。在这样的流动中,可能完全没有剪切运动。然而,由于流体正在变形(被拉伸),应变率张量 S\boldsymbol{S}S 不为零。本构关系于是预测正应力将会有所不同!对于一个速度为 u=(ϵ˙x,−ϵ˙y,0)\mathbf{u} = (\dot{\epsilon}x, -\dot{\epsilon}y, 0)u=(ϵ˙x,−ϵ˙y,0) 的平面拉伸流,该定律告诉我们 σxx−σyy=4μϵ˙\sigma_{xx} - \sigma_{yy} = 4\mu\dot{\epsilon}σxx​−σyy​=4μϵ˙。一个理想的、无粘性的流体无法支撑这样的正压力差;它只会直接被挤出去。然而,粘性流体会抵抗被拉开,在拉伸方向产生额外的张力。这在聚合物挤出和玻璃吹制等工艺中是一个至关重要的效应。

其次,搅动咖啡的能量去哪儿了?它不会无限地加速咖啡。你克服粘性力所做的功被转化为热能,使流体变暖。这个过程被称为​​粘性耗散​​(viscous dissipation)。单位体积的能量转化率由耗散函数 Φ=τ:S\Phi = \boldsymbol{\tau} : \boldsymbol{S}Φ=τ:S 给出。对于牛顿流体,这变为 Φ=2μ(S:S)\Phi = 2\mu(\boldsymbol{S}:\boldsymbol{S})Φ=2μ(S:S)。由于粘度 μ\muμ 是正的,而 S:S\boldsymbol{S}:\boldsymbol{S}S:S(平方项之和)永远不为负,耗散总是一条单行道:机械能不可逆地损失为热能。这是热力学第二定律通过流体摩擦显现自身,是熵增的无情进程。

最后,考虑一个流体流过固体边界,比如管道中的水。我们几乎总是假设紧贴壁面的流体是静止的——即​​无滑移边界条件​​(no-slip boundary condition)。这是自然界的基本定律吗?不完全是。这是一个经验观察,源于流体和固体之间强烈的分子间作用力。我们可以通过考虑有限的界面摩擦来对此建模。这揭示了滑移程度取决于一个“滑移长度”,对于常见的流固对,这个长度通常在纳米尺度。因为这个长度远小于管道或飞机的尺寸,界面处的摩擦实际上是无限大的,无滑移条件因此成为对现实极其精确的近似。

牛顿流体世界的边缘

牛顿模型以其优美的线性,为流体动力学的绝大部分内容提供了基础。但自然界充满了“叛逆者”。想想油漆,快速刷动时很容易流动,否则就很稠;或者番茄酱,只有用力摇晃瓶子才会流出来。这些都是​​非牛顿流体​​(non-Newtonian fluids)。

对于这些材料,“牛顿约定”被打破了。粘度不再是一个恒定的材料属性。取而代之的是,它们有一个随应变率变化的​​表观粘度​​(apparent viscosity)η\etaη。本构关系被推广为 τ=2η(γ˙)S\boldsymbol{\tau} = 2\eta(\dot{\gamma})\boldsymbol{S}τ=2η(γ˙​)S,其中 γ˙\dot{\gamma}γ˙​ 是应变率总大小的度量,通常定义为 γ˙=2S:S\dot{\gamma} = \sqrt{2\boldsymbol{S}:\boldsymbol{S}}γ˙​=2S:S​。对于像油漆这样的剪切稀化流体,η\etaη 随 γ˙\dot{\gamma}γ˙​ 的增加而减小。对于像玉米淀粉和水混合物这样的剪切增稠流体,η\etaη 随之增加。

通过研究这些更复杂的流体,我们对牛顿模型有了更深的理解。它代表了粘性流体最简单、最基本的行为,是线性的基石,我们对广阔而湍流的世界的理解就建立于此。它证明了一个简单、适定的物理假设如何能解锁一个充满复杂而美丽现象的宇宙。

应用与跨学科联系

在物理学的探索中,我们常常会遇到一些异常简洁却威力惊人的定律。牛顿流体的应力与应变率之间的线性关系 τ=μγ˙\tau = \mu \dot{\gamma}τ=μγ˙​ 就是一个完美的例子。乍一看,它似乎只是教科书里一个枯燥的技术定义。但这样想就错过了其中的奇妙之处。这条简单的规则实际上是一种通用语言,为各种流动物质所共有,理解它能让我们解读从微观到行星尺度的各种现象的秘密。它支配着我们血管中的血液,脚下的岩浆,以及我们日常使用的设备中的先进材料。让我们踏上一次跨学科之旅,见证这一条优雅原理的实际应用。

工程世界:从微芯片到3D打印机

我们的旅程始于人类发明的世界。考虑一下蓬勃发展的3D打印技术。设计用于挤出熔融聚合物打印头的工程师必须精确知道在给定压力下会有多少材料流出。答案隐藏在本构关系之中。在狭窄的喷嘴内,流体粘附在壁上,而在中心流速最快。这种速度差产生了剪切率。通过应用 τ=μdudr\tau = \mu \frac{du}{dr}τ=μdrdu​,可以确定任意点的剪切应力。这个应力必须与推动流体通过的压力所产生的力相平衡。从这个基本平衡出发,通过对喷嘴半径进行积分,可以建立起整个速度剖面,然后再积分一次,就能找到总体积流率。最终的结果,即著名的哈根-泊肃叶方程(Hagen-Poiseuille equation),将宏观流率与流体粘度和喷嘴几何形状直接联系起来——这是一个源于微观摩擦定律的强大设计工具。

这个原理不仅适用于管道工和打印机操作员;它也是我们最先进技术的核心。在制造微芯片的无菌洁净室中,像等离子体增强化学气相沉积(PECVD)这样的工艺被用来沉积薄得难以想象的材料薄膜。为确保硅片上数百万个晶体管获得完全均匀的涂层,工程师们使用复杂的计算机模拟。这些模拟必须求解前驱气体在整个反应室内的流动情况。该流动的控制方程——纳维-斯托克斯方程(Navier-Stokes equations)——正是建立在牛顿本构关系的基础之上,该关系定义了气体如何抵抗剪切。正确处理好这一部分,是预测热量和化学物质如何输送到晶圆表面,并最终决定最终微芯片质量和性能的第一个也是最关键的一步。从普通的管道到先进的反应堆,工程学都依赖于对流体内摩擦的理解。

运动中的地球:来自深处的低语

现在,让我们把目光从工程世界转向自然世界,将我们的尺度从毫米扩大到公里。在地球深处,熔融的岩石,即岩浆,在通往地表的途中流过管道。地球物理学家们试图理解并希望有一天能预测火山爆发,他们对这种流动进行建模。作为一阶近似,许多类型的岩浆表现得像牛顿流体,尽管其粘度极高。当岩浆被向上推动时,它会与管道的静止岩壁产生拖曳。这种拖曳力无非就是一种剪切应力。通过了解流动的速度剖面——通常接近抛物线形,就像管道中的流动一样——科学家们可以使用 τw=μ∣dudr∣wall\tau_w = \mu \left| \frac{du}{dr} \right|_{\text{wall}}τw​=μ​drdu​​wall​ 来计算施加在周围岩石上的巨大作用力。这种应力可能导致地震微震,并影响火山喷发的动力学,提醒我们,支配蜂蜜从勺子滴落的定律,同样也描述了火山的巨大威力。

生命之舞:生物学和医学中的流体

牛顿本构关系在生命领域中的作用最为密切和出人意料。生命沐浴在流体之中,也由流体构成。流体动力学原理不是外部约束;它们是生物学法则中不可或缺的一部分。

让我们从生命之河——血液开始。我们血管的内壁是一层脆弱的、单细胞厚的组织,称为内皮。这些细胞生活在一个永恒流动的世界里,它们的健康状况严重依赖于血液施加的机械力。在健康、笔直的动脉中,血液以平滑、有序的层次流动——这种状态称为层流。内皮细胞感受到的剪切应力,通常在 1−71-71−7 帕斯卡左右,是一个至关重要的、促进健康的信号。这种温和、稳定的摩擦会刺激细胞产生分子,使血管保持松弛并防止凝血。这是一种“有益的”应力,可以预防动脉粥样硬化。

但当这种流动受到干扰时会发生什么呢?在患有心脏瓣膜关闭不全的病人中,一股高压血流可能会反向涌入低压心腔。这种湍流射流会产生剪切率极高的区域,有时比正常值高出数千倍。我们简单的定律 τw=μγ˙w\tau_w = \mu \dot{\gamma}_wτw​=μγ˙​w​ 告诉我们接下来会发生什么。在粘度约为 3.2×10−3 Pa⋅s3.2 \times 10^{-3} \text{ Pa}\cdot\text{s}3.2×10−3 Pa⋅s 的血液中,近壁剪切率达到,比如说,9,000 s−19,000 \text{ s}^{-1}9,000 s−1,会产生近 30 Pa30 \text{ Pa}30 Pa 的剪切应力。这不是温和的、刺激性的摩擦;这是一种物理上的“喷砂”。这种病理性的高应力会物理性地损伤并剥脱内皮衬里,形成一个创面。这个损伤部位成为血小板聚集和纤维蛋白沉积的病灶,形成一个无菌赘生物,在短暂感染期间,可能被细菌定植,导致危险的感染性心内膜炎。同一个物理定律既解释了健康血流的保护作用,也解释了疾病的破坏机制。

故事变得更加微妙。细胞不仅仅是在应力下生存;它们会“读取”应力,并以惊人的精确度做出反应。这个过程,即力学转导(mechanotransduction),是物理学如何转化为生物化学的方式。在发育中的斑马鱼胚胎中,最初的原始血管中的剪切应力(直接根据血液粘度和速度梯度计算得出)提供了必要的信号,指导动脉和静脉如何正确形成。在成熟的动脉中,剪切应力的大小控制着像 KLF2 这样的关键保护性基因的表达。在实验室里,科学家可以构建他们能够完全控制流动的微流控通道。通过将培养基灌注过一层内皮细胞,他们可以精确设定壁面剪切应力。他们观察到,随着应力增加,像 KLF2 这样的基因的激活遵循一条可预测的、开关般的曲线。这使他们能够建立定量模型,如希尔函数(Hill functions),将物理输入 τw\tau_wτw​ 与生物输出(基因激活水平)直接联系起来。物理学成为了生命方程中的一个输入变量。

流体剪切力的影响渗透到生物学的各个方面。在胚胎发育期间,我们滑膜关节(如膝关节、肘关节)的形成就依赖于它。胎儿的运动剪切了新生骨骼之间的原始液体,产生的应力是清除其间细胞并塑造关节腔的必要信号。根据我们的定律,增加流体粘度的基因缺陷会按比例增加这种剪切应力,从而可能扰乱这一精细的发育过程。即使在牙科诊所,这种物理学也在发挥作用。在根管治疗过程中,一个关键目标是清除附着在根管壁上的细菌生物膜。冲洗液不仅是化学消毒剂,也是一种机械工具。清洁的效果取决于流动是否能产生足够的壁面剪切应力来克服生物膜的附着力并将其冲刷掉。

最后,考虑一个常规的医疗程序:制备细胞学涂片。技术员取一滴样本,用第二张载玻片将其涂成薄膜。这个简单的动作创造了一个平面库埃特流(Couette flow)。涂抹的速度 UUU 和间隙的高度 hhh 定义了一个剪切率 γ˙≈U/h\dot{\gamma} \approx U/hγ˙​≈U/h。样本中的细胞经历的剪切应力为 τ=μ(U/h)\tau = \mu (U/h)τ=μ(U/h)。如果技术员涂抹得太快,这个应力可能会超过细胞膜破裂的临界阈值,从而摧毁他们希望检查的细胞,使样本无法用于诊断。

简单的定律,复杂的世界

我们的旅程带领我们从微芯片的核心到火山的中心,再到活细胞的心脏。在每一种情况下,牛顿流体的应力与剪切率之间的简单线性关系都不只是一个偶然特征;它是故事的核心角色。它在工程学中是设计原则,在地球物理学中是诊断线索,在医学中是生命信号和疾病机制。发现如此多样和复杂的现象都可以通过这样一个优雅而简单的物理定律来理解,这是对科学统一性与美感的深刻证明。它提醒我们,如果我们学会说自然基本法则的语言,我们就能开始理解那些既陌生又熟悉的世界。