玻尔百科为了描述物理世界，我们必须捕捉变化的本质，但“变化”本身是一个难以捉摸的概念。一个属性的变化是因为时间的流逝，还是因为我们移动到了一个该属性有所不同的新位置？物质导数（substantive derivative），也称为材料导数（material derivative），是一个强大的数学算子，旨在精确回答这个问题。它提供了一种语言，用以描述一个粒子或一个材料元在时空中运动时所经历的变化率。本文旨在弥合“随波逐流”的直观概念与对其进行建模所需的严谨数学之间的鸿沟，揭示一个在现代物理学和工程学中处于核心地位的工具。
在接下来的章节中，我们将踏上一段旅程，以全面理解这一基本概念。首先，在“原理与机制”一章中，我们将把物质导数分解为其组成部分——局部变化率和对流变化率。我们将探索其与流动几何学的优美联系，并揭示一个与物理客观性相关的深层危机，该危机导致了更复杂的“客观”导数的发展。随后，在“应用与跨学科联系”一章中，我们将见证物质导数的实际应用，了解这一个单一概念如何为描述材料变形、流体运动、机器人控制乃至前沿的计算机模拟提供统一的框架。
在理解世界的旅程中，我们面临着一个视角选择。我们可以静立不动，观察世界从我们身边流过；或者，我们可以随波逐流，被水流携带着前行。在物理学中，这两种视角被称为欧拉（Eulerian）观点和拉格朗日（Lagrangian）观点。前者就像站在桥上的一名记者，记录着那一个单点随时间变化的水位和速度。后者则像瓶中的信，亲身体验着河流的曲折和转弯。物质导数，也称为材料导数，正是这瓶中信的数学语言。它告诉我们任何属性——无论是温度、压力还是更奇特的量——在粒子运动过程中所经历的变化率。
想象你是一只小虫，站在一块凉爽的大金属板上。有人在板下打开了加热器，于是整块板开始变暖。即使你站着完全不动，你脚下的温度也会升高。这种发生在固定位置的变化，就是局部或时间变化率。这是你只看手表就会注意到的变化。在数学上，它就是我们熟悉的对时间的偏导数，。
但现在，假设这块板不仅在升温，而且还是倾斜的，顶部比底部热。如果你开始向上坡走，你感觉到的温度会升高，不是因为你所在的位置在变热，而仅仅是因为你移动到了一个更热的位置。这种因你在空间变化场中移动而引起的变化，称为对流变化率。
你所经历的总变化是这两种效应的总和：板子各处都在变暖引起的变化（局部）加上你移动到不同位置引起的变化（对流）。这个总变化率正是物质导数。对于任何一个依赖于位置 和时间 的属性，我们称之为 ，其物质导数写作 ，由一个非常直观的公式给出：
此处， 是粒子的速度（你走路的速度），而 是属性 的梯度（板上温度随位置变化的陡峭程度）。点积 巧妙地提取了你移动方向上的变化率。因此，如果你沿着一条恒温路径（等温线）行走，即使你在移动，这个对流项也为零！这个简单方程中所蕴含的物理意义是深远的。它将拉格朗日经验（）与欧拉场描述（ 和 ）联系起来。你在物理课程中可能遇到的问题，常常涉及分别计算这两个部分，然后将它们相加，以求得流体粒子所经历的总变化率，从而强化了这一基本分解。
当我们把物质导数应用于描述流动本身几何形态的量时，它的威力才真正显现出来。让我们不再考虑单个粒子，而是想象一小团虚构的流体——一小群粒子。当这团流体在流动中移动和翻滚时，它可能会拉伸、剪切，以及最重要的是，改变其体积。
在连续介质力学中，这团流体当前体积与其初始体积之比，由一个称为雅可比行列式（Jacobian）的量来描述，记为 。如果 ，这团流体就膨胀了；如果 ，它就被压缩了。现在，让我们问一个自然的问题：我们的流体团在流动过程中的体积变化速度有多快？这正是物质导数的用武之地！我们想计算 。
当你进行数学推导时，一个惊人优雅的结果出现了：
让我们来解析这个公式。左边的项是从拉格朗日视角看的我们流体团的体积变化率。右边的项 是速度场的散度（divergence）。散度是一个纯粹的欧拉量；它衡量在空间中某一点，速度矢量“散开”或“发散”的程度。你可以把它想象成一个源或一个汇。如果某点的散度为正，流体就从该点向外膨胀。如果为负，流体就向该点汇聚。
这个方程在两种视角之间架起了一座金桥。它表明，一个流体元的体积分数变化率（）恰好等于该流体元所在位置速度场的散度。这是流体动力学的一块基石。对于像水这样几乎不可压缩的流体，任何给定流体团的体积都不能改变。这意味着 必须为零，这反过来又给速度场施加了一条基本定律：。任何不可压缩流体的有效流动都必须遵守这个条件。物质导数引导我们得出了关于流动本质的一个深刻的物理约束。
到目前为止，我们的旅程一帆风顺。物质导数似乎是一个完美的工具。但一个微妙而深刻的问题潜伏在表面之下。这个问题迫使我们质疑我们所说的“变化率”到底意味着什么。
从牛顿到爱因斯坦，现代物理学的基础建立在一个简单而强大的理念上：物理定律不应依赖于观察者。无论我是在伦敦的实验室里做实验，还是在一辆平稳行驶的火车上的实验室里做实验，结果都应该是一样的。这个原则也适用于相互之间存在旋转的观察者。一个能够在这些观察者之间正确变换的物理属性或定律，被称为客观的（objective）或标架无关的（frame-indifferent）。例如，矢量是客观的；处在不同旋转坐标系中的观察者会用不同的分量来描述它，但他们描述的都是同一个物理上的箭头，并且这些分量通过一个简单的旋转规则进行变换。
现在来看危机所在：物质导数是客观的吗？让我们想象一块透明的果冻，它被保持在一个静态的、略微拉伸的状态。在它的静止坐标系中，应变张量 是恒定的，所以它的物质导数 显然为零。什么都没有发生。但现在，让我们从一个旋转的坐标系来看这块果冻——想象一下从一个旋转木马上观察它。从我们旋转的视角看，这块果冻似乎在旋转。它的应变张量的分量 现在随时间变化，因为张量本身相对于我们旋转的基矢量在旋转。如果我们天真地计算这个变化张量的物质时间导数 ，我们会得到一个非零的答案！
这是一场灾难。我们为一个物理上没有任何变化的过程计算出了一个非零的变化率。物质时间导数未能通过客观性的检验。它给出的答案依赖于观察者的任意旋转。如果我们考虑应力变化率 ，也会出现同样的失败。这种失败的原因在于，简单的物质导数将一个量的“真实”物理变化与由材料元本身刚体旋转引起的平凡变化纠缠在了一起。
我们该如何解决这个问题？答案既优雅又巧妙，正如问题本身一样。如果简单的导数因为它被旋转所迷惑而失效，那么我们必须发明一个不受其影响的“更聪明”的导数。我们需要一种方法，从一个与材料一同旋转的视角来衡量张量的变化率。
这样想：要判断一个在旋转木马上的人是否在挥舞手臂，你不能只是站在地面上观察。他们的圆周运动会淹没细节。最好的方法是和他们一起登上旋转木马。从这个同旋（co-rotating）的视角看，他们整体的圆周运动消失了，你可以清楚地看到他们手臂的“真实”运动。
物理学家们正是发展了这一思想。他们定义了新的导数算子，称为客观变化率（objective rates）或同旋变化率（co-rotational rates），这些算子在数学上实现了这一点。像Jaumann率或上随体导数率（upper convected rate）这样的变化率，其构造方法是从简单的物质导数开始，然后“减去”那些仅仅由材料局部旋转引起的部分。当你这样做时，之前困扰我们的非客观项会完美抵消，最终得到一个真正客观的变化率。
这似乎像一个抽象的数学游戏，但对于描述真实世界来说，它绝对是必不可少的。为了写出描述像钢、橡胶或面团这类材料在应力下如何变形的定律（即本构定律），我们需要将应力的变化率与变形的速率联系起来。如果我们使用非客观的物质导数，我们的物理定律将会为在旋转平台上的观察者给出不同的预测，这显然是荒谬的。通过使用客观变化率，我们确保了我们对材料行为的描述是普适的，是一条真正的自然法则，不受写下它的物理学家的任意运动的影响。物质导数的旅程，从一个简单、直观的想法，到一个深刻的危机及其优雅的解决方案，揭示了我们对世界物理描述中固有的深层结构和美。
在我们之前的讨论中，我们遇到了一个非凡的数学工具：物质导数。我们看到它提供了一种方式来讲述运动的语言，去询问一个量如何变化——不是在空间中的一个固定点，而是对于一块“物质”在流动的河流中冲浪时如何变化。它是在一个点上发生的变化与被冲到新点所带来的变化之和：。现在，你可能会想，这只是微积分里一个巧妙的技巧，但它有什么用处呢？事实证明，答案几乎涵盖了所有移动、流动、弯曲或断裂的事物。物质导数不仅仅是一个公式，它是一把钥匙，解锁了对从晶体微观变形到星系旋转之舞等各种现象的统一理解。让我们踏上旅程，看看这把钥匙能打开哪些门。
让我们从最直接的问题开始：当某物流动时，它的形状是如何改变的？想象一下，你在一个平静池塘里两个相邻水分子之间画了一个微小的直箭头。现在，水流开始了。这两个分子开始移动，连接它们的箭头——这个无穷小的物质线元，我们称之为——将被携带前行。它的长度会保持不变吗？它会指向同一个方向吗？
物质导数以惊人的优雅给出了答案。我们的小箭头的变化率不是任意的；它完全由流动的速度场决定。具体来说，可以证明箭头的变化率由给出。这个紧凑的表达式意义深远。它告诉我们，连接两个粒子的矢量的变化是由于它们之间速度的差异造成的。如果各处速度都相同，线元只是平移，保持不变。但如果速度因地而异——即存在速度梯度——那么我们的小箭头就会被拉伸、剪切和旋转。
这一个思想是连续介质力学的基础。速度梯度张量成为了主角。我们可以把它看作一个“运动分解器”。它可以被分成两部分。对称部分，称为应变率张量，告诉我们关于拉伸的信息。它是一个数学机器，当你给它一个方向，它会告诉你材料在该方向上伸长的速度有多快。事实上，如果我们要求物质箭头长度平方的变化率，我们会发现它与一个包含的量成正比。流动中最大拉伸的方向原来就是这个优美张量的主轴（特征向量）。因此，物质导数揭示了任何流动内部隐藏的几何结构，精确地告诉我们它是如何变形的。
那么速度梯度的另一部分呢？反对称部分，称为自旋或涡量张量，告诉我们另一种基本运动类型：局部旋转。它衡量一个材料元旋转的速度有多快，就像一片被卷入漩涡的小叶子。流体动力学中的一个关键问题是这个“自旋”如何演变。它会变强吗？变弱吗？控制这一变化的方程，源于对涡量张量本身取物质导数，是湍流和流体动力学研究的支柱之一。
理解运动的几何学是一回事，但物理学也关乎力、能量和材料的属性。在这里，物质导数同样不可或缺。
当你拉伸一根橡皮筋时，它会储存能量。材料科学主要关注的是找到连接施加的力（应力）与产生的变形（应变）的定律。这些定律不是静态的；它们描述一个动态过程。你通过变形对一块材料做功的速率，必须等于其内部储存能量增加的速率，加上任何以热量形式损失的能量。对于弹性材料，这种平衡通过物质导数优美地表达出来：应力功率等于应变能密度的物质时间导数 。这不仅仅是一个记账原则；它是在现代力学中从能量推导应力的根本定义。
此外，当物体变形时，我们可以选择用不同的方式来衡量变形。我们可以将其当前形状与其原始未变形形状进行比较（拉格朗日观点），或者我们可以描述其当前变形状态下某一点的应变（欧拉观点）。这些不同的应变度量，如格林-拉格朗日张量或柯西-格林张量，必须是相互关联的。物质导数提供了在它们的变化率之间进行转换的运动学词典。像或其欧拉对应物的表达式，是让工程师能够为像聚合物和生物组织这样的复杂材料制定定律的关键机制，这些材料的行为取决于它们的整个变形历史。这些关系也是计算机程序中的主力，用于模拟从车祸到钢梁锻造的一切，在这些模拟中，材料中的应力必须在每个微小的时间步长上更新。
物质导数的影响远不止于实验室的工作台。想象一下你站在我们旋转的地球上。你感受到的离心力取决于你与自转轴的距离。现在，如果你开始移动，比如说，坐在一列高速火车上呢？你体验到的离心力会改变，不仅因为时间在流逝，也因为你正在移动到一个力不同的位置。要找出运动物体感受到的离心加速度的变化率，我们必须使用物质导数。计算揭示了一个美妙的对称结果：在一个以角速度旋转的参考系中，一个速度为的粒子的离心场变化率就是。这不仅仅是一个数学上的奇趣。这个“变化的化”，与科里奥利力密切相关，后者是导致飓风大规模旋转和大规模洋流的著名的“虚拟”力。物质导数正确地解释了非惯性世界中的物理现象。
一个真正伟大思想的力量在于它能超越其最初的语境。“追随一个变化”的概念并不局限于流体或固体中的粒子。它可以应用于任何随时间演变的系统的状态。
考虑一个现代机器人手臂，一个由连杆和关节组成的工程奇迹。机器人手的位置由其关节的角度决定。手部速度与关节速度之间的关系由一个称为雅可比矩阵的矩阵给出。为了让机器人平稳快速地移动，我们不仅需要控制其速度，还需要控制其加速度。为了找到手的加速度，必须对速度方程求时间导数。这不可避免地涉及求雅可比矩阵本身的变化率。但雅可比矩阵是关节角度的函数，而这些角度正在变化！因此，我们需要它的物质导数——即当我们跟随机器ンの构型穿过其状态空间时，雅可比矩阵的变化率。在这里，“流动”不是在物理空间中，而是在所有可能的机器人构型的抽象空间中。同样的数学工具同样适用。
最后，物质导数位于最先进的计算机模拟的核心。在模拟像风中飘扬的旗帜或动脉中搏动的血液等现象时，物理域本身正在变形。固定的计算网格（欧拉方法）效率低下，而完全随材料变形的网格（拉格朗 日方法）可能会变得纠缠不清。解决方法是任意拉格朗日-欧拉（ALE）方法，其中计算网格以一种“智能”的方式移动，独立于物质流动，以保持高质量的网格。但这带来了一个难题：当你的观测网格也在移动时，你如何计算一个物质粒子的真实变化率？答案是一个广义形式的物质导数。它优美地考虑了三件事：从移动网格上的一个点看到的变化，由于材料相对于网格流动而引起的变化，以及由于网格自身变形引起的梯度变换。正是这个复杂的工具，让计算科学家能够解决科学和工程中一些最复杂的流固耦合问题。
从面团的拉伸到行星的旋转，从机器人的控制到跳动心脏的模拟，物质导数证明了它不仅仅是一个符号。它是一个捕捉动态世界中变化本质的基本概念，在一个广阔的科学探究领域中揭示了令人惊讶和美丽的统一性。