
在研究流体和固体等变形材料时,准确描述物理性质如何变化是一项基础而又复杂的挑战。物理定律应在任何观察者参考系下都成立——这一概念被称为物质标架无关性原理,或称客观性原理。然而,测量变化最直观的方式——物质时间导数——却无法满足这一要求,因为其值取决于观察者的旋转运动。这种差异使得我们在为连续介质力学建立普适物理定律时面临着巨大的障碍。
本文通过引入一系列客观时间导数来解决这个问题,这些数学工具旨在提供一种与标架无关的变化量度。在接下来的章节中,您将对这些关键概念获得全面的理解。“原理与机制”一节将剖析观察者依赖性问题,并引出下随体导数,解释其数学结构以及与变形表面的深刻几何联系。随后,“应用与跨学科联系”一节将展示这些理论结构如何应用于模拟从高分子流变学到地球物理学等领域中真实材料的行为,从而说明它们在现代科学与工程中不可或缺的作用。
要理解世界,我们必须衡量它如何变化。但这看似简单的任务却出人意料地微妙。想象一下,在一张橡胶片上画一幅画,然后拉伸并旋转这张橡胶片。如果你问:“这幅画是如何变化的?”,答案完全取决于你的视角。你是一个站在橡胶片上,随之一起被拉伸和旋转的微小观察者吗?还是一个从远处观看的静止观察者?我们相信,物理定律不应依赖于观察者的任意选择。这一基本思想,被称为物质标架无关性原理,或简称为客观性,是我们探索之旅的起点。
在流体和固体等连续材料的力学中,我们经常追踪由称为张量的数学对象表示的物理性质。你可以将一个简单的矢量(如速度)看作一阶张量。但许多性质更为复杂。例如,材料内部的应力——即粒子间相互施加的内力——是一个二阶张量,因为在任何一点,都有作用在不同方向表面上的力。另一个例子是描述流体中高分子分子平均取向和拉伸的张量。
测量这样一个张量 变化率最直接的方法是跟随一小块材料,观察 如何随时间变化。这被称为物质时间导数,记作 。它告诉你,如果你像冲浪一样随着一个质点在空间中运动,你会看到的变化率。
但问题在于:这个简单的导数不是客观的。为了说明原因,考虑一个正在进行纯刚体旋转的流体,就像杯子里被搅拌的咖啡。想象一个张量 ,它代表某种固定的微观结构特征,比如微小、冻结在其中的纤维的排列。对于一个与流体一同旋转的观察者来说,这个特征是恒定的;没有任何变化。然而,对于一个观察杯子的静止观察者来说,这些纤维显然在旋转。物质导数 会捕捉到这种旋转并且不为零。事实证明,对于由自旋张量 描述的纯旋转,物质导数为 。这种变化率在物理意义上并非“真实”的;它只是我们静止视角的产物。它告诉我们的是物体的旋转,而不是其内部的任何内在变化。由于物理定律应该描述内在变化,因此物质导数 不适合用来构建物理定律。
为了构建独立于观察者运动的自然法则,我们需要发明新类型的导数——客观时间导数。这些数学工具被巧妙地设计用来减去由刚性旋转引起的“伪”变化,只留下对应于材料实际拉伸和变形的变化。它们就像特殊的相机滤镜,可以去除相机本身的运动,让我们只看到场景中真正发生的事情。
实现这一点的方法不止一种,从而产生了一族客观率,每一种都有其自身关于何为“真实”变化的观点。这一族中最重要的成员是Jaumann率、上随体导数和下随体导数。
要理解它们的结构,我们首先需要剖析运动本身。流体的局部运动由速度梯度张量 描述。这个张量包含了速度如何随点变化的全部信息。我们可以将 分解为两部分:
通过这种分解,我们可以看到不同的客观率是如何构建的。
Jaumann率(或共旋率)是最直观的修正方法。它取物质导数 ,然后简单地减去由纯旋转引起的部分:
这个变化率衡量的是一个与局部材料元一起旋转但不随之拉伸的观察者所看到的变化。
然而,随体导数则更为深刻。它们描述的是一个不仅旋转,而且还随材料一同拉伸和变形的观察者视角下的变化率。它们有两种形式:上随体和下随体。
上随体导数定义为:
而下随体导数定义为:
这些定义起初可能看起来有些随意,但它们正是确保客观性所必需的。涉及 和 的附加项恰好可以完美抵消所有因观察者改变而产生的非客观项。如果我们使用 分解为 和 的形式,这些导数的结构就变得非常清晰:
这个优美的结果表明,这三种率都共享对自旋的核心修正(即Jaumann率)。它们的不同之处仅在于如何处理拉伸 。上随体率减去拉伸效应,而下随体率则加上它们。两者之间的差异与变形直接相关;例如,对于应变率张量 本身,它们的差为 。
那么,下随体导数背后的物理意义是什么?为什么我们要加上与拉伸相关的项?答案在于我们所测量的量的几何性质。
下随体导数是测量协变张量变化率的自然方式。一种理解协变张量的简单方法是将其视为与嵌入材料中的面元相关联的东西。想象一个画在变形橡胶片上的小方块。当橡胶片变形时,这个方块也随之变形。下随体导数测量的是与这个变形方块相关的属性(如其面积或作用于单位面积上的力)从方块本身视角看是如何变化的。因为这个视角“附着”在变形的材料上,所以得到的变化率是内在地客观的。
通过几个关键例子,这种几何解释得到了优美的阐释。在力学中,单位张量 可以被认为是空间的度规——它是我们用来测量长度和角度的工具。如果我们问这个度规随着材料流动如何变化,我们必须使用一个客观导数。计算单位张量的下随体导数会得到一个惊人简洁的结果:
这个方程是一个深刻的陈述:我们空间测量尺的客观变化率恰好是材料本身拉伸和变形率的两倍。该方程的迹(对角元素之和)告诉我们,度规体积的变化率是材料元体积变化率()的两倍,直接将空间几何与可压缩性的物理特性联系起来。
下随体导数的真正威力在于它得到零结果的时候。在连续介质力学中,Finger张量 是一种变形的度量。它是一个协变张量,实质上追踪了面元从其在未变形参考状态下的原始形状被拉伸和剪切的情况。如果我们计算它的下随体导数,会发现:
这个结果远非微不足道。它意味着从一个嵌入在材料表面并随之变形的观察者视角来看,Finger张量是恒定的。对于这个视角而言,它是自然的、不变的应变度量。这就像有一把用完美弹性橡胶制成的尺子;即使尺子被拉伸到其两倍长度,生活在尺子上的蚂蚁仍然会测量其长度为“一个尺子长度”。正是这种一致性使得Finger张量和下随体导数的组合在为粘弹性流体和橡胶状固体等材料构建本构模型时如此强大。
这也阐明了为什么不同的客观率是不可互换的。对于纯刚体旋转,材料旋转但不拉伸(),三种客观率——Jaumann率、上随体率和下随体率——变得完全相同。正如预期的那样,它们都正确地报告了一个仅随物体旋转的张量的变化率为零,而非客观的物质导数则报告了一个非零的变化率。在物理定律中选择使用哪种导数取决于所讨论张量的几何特性:它是否与线相关(逆变的,使用上随体导数),与面相关(协变的,使用下随体导数),还是完全是其他类型?物理学与几何学之间这种深刻的联系揭示了连续介质力学潜在的统一性与优雅。
在深入探讨了客观导数的原理和机制之后,您可能会倾向于将它们视为一种有些深奥的数学工具。但事实远非如此。这些概念不仅仅是抽象的形式主义;它们正是我们用来描述真实材料——从洗发水瓶中的高分子链到我们星球的地幔岩石——在流动、拉伸和旋转的胁迫下如何表现的语言。正是在这里,物理学才真正鲜活起来,将我们的数学框架与可触及、可测量的现象联系起来。这是一段从张量的抽象世界到工程、地质和化学实践世界的旅程。
首先,让我们问一个基本问题:一个导数是“客观的”意味着什么?想象你是一个微小的观察者,一个麦克斯韦妖之类的存在,骑在一小团流体上。这团流体在移动时可能会剧烈地翻滚和旋转。如果流体的某个性质,比如它的内应力,在你这个翻滚的参考系中看起来是恒定的,那么你测量的任何“客观”变化率都必须为零。这是客观性的基本检验。
在纯刚体旋转中,材料的每个部分都一同旋转,没有任何拉伸或扭曲。对于任何只是被这种旋转“携带”的张量属性——意味着它在共旋参考系中是恒定的——其物质导数绝对不为零。在一个固定的实验室参考系中的观察者会看到张量的分量在变化,仅仅因为物体在转动。客观导数,如上随体、下随体和Jaumann率的魔力在于,它们被专门设计用来减去这种因纯旋转而产生的虚假变化。当应用于刚体旋转中的这样一个张量时,它们都正确地得出一个零值,通过了作为真实、与参考系无关的物理测量的石蕊试纸测试。它们捕捉的是材料变形所固有的变化率,而不是观察者的视角。
乍一看,多种客观导数(如上随体率和下随体率)的存在似乎是一种不必要的复杂化。为什么不止一种?答案揭示了连续介质力学核心中一种优美而深刻的二元性,一种反映了几何对象被变形介质输运的不同方式的对称性。
想象一条注入透明流动糖浆中的染料线。随着糖浆变形,这条线被拉伸和旋转。连接这条线上两点的矢量是一个逆变矢量,一个被流动“前推”和拉伸的对象。由这类矢量构建的结构(如高分子链的构象)的演化,最自然地由上随体导数描述。事实上,如果我们将高分子溶液中的长链分子建模为被流动仿射携带的微小哑铃,它们的平均形状,由一个构象张量 描述,其演化方式恰好使其上随体导数为零(在没有热松弛的情况下)。这种导数是为那些随材料拉伸的事物量身定做的。
现在,考虑一个不同的对象:一个嵌入在同一流动糖浆中的小面元。可以把它想象成一个小立方体的面。当材料变形时,这个面不仅移动,而且其方向和面积也会改变。这个面的法向量是一个协变矢量;它的变换方式与拉伸的线元不同。它的演化由从参考状态的“拉回”来描述。下随体导数是描述由这些协变对象构建的张量演化的自然语言。的确,如果构建一个代表随体材料表面的张量,其下随体导数被发现为零。这种二元性不是数学上的怪癖;它是变形基本几何结构的反映。
这种区别并不仅仅是学术性的。在流变学这门研究物质流动与变形的科学中,客观导数的选择具有戏剧性的、可实验验证的后果。考虑一个简单的“Maxwell”流体,这是一个用于描述像高分子熔体或Boger流体这类同时表现出粘性(类液体)和弹性(类固体)特性的材料的基本模型。
如果我们使用上随体导数构建一个本构模型(上随体Maxwell或UCM模型),并将其与使用下随体导数的模型(下随体Maxwell或LCM模型)进行比较,它们对流体在简单剪切流(如两平行板之间的流动)中的行为做出了惊人不同的预测。虽然两种模型都预测了相同的剪切应力和相同的第一法向应力差 (),后者是导致像Weissenberg效应(流体爬上旋转杆)这类现象的原因,但它们在第二法向应力差 () 上存在分歧。UCM模型预测 ,而LCM模型预测一个非零的负值 。对高分子溶液的实验表明, 确实非零,并且通常比 小得多。这告诉我们,这两种简单的模型都不完美,但它们的差异是真实且可测量的,直接源于客观率的选择。
Oldroyd-B模型是高分子流变学的基石,它优雅地结合了这些思想。它将稀高分子溶液建模为一种简单牛顿溶剂(如水)和一种聚合物组分的混合物,后者的应力由上随体Maxwell模型控制。这一选择有其物理动机:上随体导数正确地捕捉了高分子链逆变构象张量的仿射输运。
这些概念的用途远远超出了高分子溶液。在地球物理学中,模拟地球地幔的流动或构造板块的变形涉及在巨大的压力和温度下、在地质时间尺度上作用的材料。在这些条件下,岩石可以流动并经历巨大的旋转。为了建立一个关联岩石应力率和变形率的“亚弹性”定律,必须使用一个客观应力率来正确处理构造块的大旋转。地球物理学家使用了多种变化率,包括Jaumann率、Green-Naghdi率(基于材料内部坐标系的真实旋转)和Truesdell率,每种都有其自身的细微差别和适用范围。对于准确预测地震和造山等地质过程中的应力积累和释放,这种选择至关重要。
在现代,这些理论原理对价值数十亿美元的计算机辅助工程产业产生了直接影响。在编写计算流体动力学(CFD)或计算固体力学(CSM)软件时,确保材料客观性不是一个选项;它是有效模拟的必要条件。一个搅拌釜反应器或汽车轮胎的模拟,如果其材料应力没有使用客观公式,将会产生依赖于实验室坐标系任意方向的、物理上无意义的结果。客观性的数学验证是开发和验证高级计算代码的标准程序。
最后,对更好本构模型的追求将我们带到了一个引人入胜的计算和物理难题:高魏森贝格数问题(HWNP)。魏森贝格数 () 是一个无量纲数,它比较了材料的特征松弛时间与施加流动的速率。在高 下,流动速度快,材料没有时间松弛,导致高度弹性的行为。
简单而优雅的Oldroyd-B模型,凭借其上随体导数,做出了一个可怕的预测:在强拉伸流(一种拉伸流动,像拉太妃糖)中,一旦 超过一个临界值(对于平面拉伸为0.5),聚合物应力预计将无限增长!这种“应力爆增”是该模型的一个特征,反映了一个真实的物理现象,称为线圈-拉伸转变,即高分子链突然解开。然而,这种数学上的奇异性对数值模拟构成了巨大的挑战,导致它们灾难性地失败。
这并非客观导数概念的失败,而是我们物理模型过于简单的标志。自然界有更聪明的方式来处理这个问题。Giesekus模型引入了一个更现实的物理成分:各向异性阻力。它假定高分子链感受到的摩擦力取决于其相对于邻近分子的方向。这个看似微小的改进带来了深远的数学后果。它在应力的演化方程中引入了一个*非线性阻尼项*。这个项是二次的,并且作为一个强大的“汇”,随着聚合物的拉伸而迅速增长,从而阻止应力达到无穷大。它驯服了简单模型预测的灾难性爆增,导致应力饱和,这既更符合物理现实,也更易于计算管理。
这种美妙的相互作用——更深刻的物理洞察导致更稳健的数学模型,从而解决了一个棘手的计算问题——是现代连续介质物理学的标志。它表明,客观导数家族并非教科书中一个已完结的章节,而是科学家和工程师们推动我们理解和模拟能力边界的一个活跃且必不可少的工具包。