🔢 mathematics

Kinetic derivation of thermal viscous models for nematic liquid crystal dynamics

本文通过对具有 BGK 型碰撞的动力学模型应用 Chapman-Enskog 展开和约束最大化方法，推导出了向列型液晶的宏观热力学理论，从而将现有的无粘理论推广到包含可压缩及不可压缩流体的粘性、热效应和自旋扩散效应。

原作者： P. E. Farrell, J. Málek, O. Souček, U. Zerbinati

发布于 2026-07-01

📖 1 分钟阅读🧠 深度阅读

原作者： P. E. Farrell, J. Málek, O. Souček, U. Zerbinati

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一个装满微小、刚性棒状物（类似于火柴棒）并在流体中漂浮的盒子。有时，这些棒状物会随机向各个方向漂浮（就像在盒子里摇晃的火柴）；而有时，它们会突然排列整齐，指向同一个方向，从而形成一种“液晶”态。这篇论文是一套数学配方，用于精确预测这些棒状物如何运动、旋转以及如何与热量和压力相互作用，特别是当它们正处于开始排列的过程时。

以下是这篇论文的故事，通过简单的概念进行了拆解：

1. 起点：“微观”视角

作者从单个分子的角度观察流体。他们使用了一种名为**动力学理论（Kinetic Theory）**的数学工具，这就像是在追踪盒子里每一根火柴的位置和速度。

问题： 追踪数十亿根棒状物是不可能的。
解决方案： 他们使用了一个“碰撞算子”（描述棒状物如何相互碰撞的规则）。他们没有计算每一次微小的碰撞，而是使用了一个简化的规则，称为 BGK 算子。你可以把它看作是一个“弛豫规则”：如果棒状物的排列不一致，这个规则会将它们轻轻地推向一个舒适的、平均的状态，就像一群人在行走时自然而然地找到一种共同节奏一样。

2. 巨大的飞跃：从微观到宏观

这篇论文的主要成就，是将单个棒状物的混乱舞蹈转化为平滑的大尺度定律（例如你可能在天气预报中见过的流体动力学定律）。

方法： 他们使用了 Chapman–Enskog 展开技术。想象一下，从一张高分辨率的人群照片缩放到一个模糊的视图，其中你只能看到人群的整体流动。这种方法使他们能够从微观规则推导出关于质量、动量、能量和熵（无序度）的方程。
转折点： 之前的模型假设棒状物已经排列整齐了。而这篇论文很特别，因为它解释了它们最初是如何排列起来的。它考虑了“涌现”行为，即棒状物由于彼此的相互作用而自发决定对齐，而不是被强迫对齐。

3. “交通警察”：熵最大化

一旦有了描述流体如何运动的方程，他们就需要确定具体的“本构关系”——即告诉我们流体如何抵抗流动（粘性）或传导热量的规则。

策略： 他们使用了一个名为**熵产生最大化（Maximization of Entropy Production）**的原理。
类比： 想象一条河流顺着山坡流下。水想尽可能快地到达底部，但它必须应对岩石和摩擦力。自然界似乎会选择那条能以最有效的方式产生最多“摩擦”（熵）的路径。作者利用这个想法，在数学上“锁定”了关于液晶如何流动和传导热量的正确公式。

4. 新的成分：粘性、热量与自旋

论文将一个较旧、较简单的模型（该模型仅适用于“无粘”或无摩擦流体）进行了泛化，使其包含了现实世界的影响：

粘性（摩擦）： 流体抵抗流动。论文根据棒状物的排列情况推导出了这种阻力的特定公式。
热传导： 流体携带热量。论文展示了棒状物的排列如何影响热量在流体中的移动。
自旋扩散（Spin Diffusion）： 这是液晶的一个独特特征。因为棒状物具有方向性，它们可以相对于流动进行“自旋”或旋转。论文推导出了一个关于这种旋转运动如何扩散并减慢的新规则，类似于旋转陀螺最终如何摇晃并停止。

5. 结果：一套完整的“行动计划”

作者最终得到了一套完整的方程组（一个“封闭系统”），可以描述：

可压缩流体： 可以被挤压的流体（如气体）。
不可压缩流体： 无法被挤压的流体（如水）。

这些方程是用于模拟液晶的著名的 Leslie–Ericksen 方程的高级版本。这个新版本是“热力学一致的”，这意味着它严格遵守有关能量和热量的物理定律。

简而言之

你可以把这篇论文看作是对一种复杂流体的说明书进行的升级。

旧版说明书： “假设棒状物已经排列整齐。忽略摩擦和热量。”
新版说明书（本论文）： “这里有关于棒状物如何自主决定对齐的精确过程。这里有它们如何互相摩擦（摩擦力）、如何携带热量，以及它们的旋转运动如何减慢。我们是从底层开始，从单个碰撞的物理学出发，推导出了这些规则。”

这篇论文并不声称解决了特定的现实工程问题（如制造更好的屏幕）或医学应用。相反，它提供了关于这些复杂流体如何行为的严谨数学基础和“源代码”，确保未来基于这项工作的任何模拟或模型在物理上都是准确的。

技术摘要：向列型液晶动力学热粘性模型的动力学推导

问题陈述
近期的研究建立了一个基于动力学理论推导的、可压缩且无粘性的向列型液晶 Leslie–Ericksen 方程变体。然而，该方法依赖于围绕已对齐状态进行的矩闭合，实际上是预设了集体取向行为，而非允许其从动力学中涌现。此外，所得模型是严格无粘性的，未能捕捉到真实液晶流体中固有的耗散、热效应及自旋扩散效应。本文旨在解决这些局限性，通过推导一个考虑了对齐涌现以及热与粘性效应的可压缩粘性 Leslie–Ericksen 方程。

方法论
作者从有序流体的动力学理论出发，开发了一种宏观热力学理论。该方法通过以下阶段进行：

动力学框架： 出发点是 [7] 中提出的有序流体动力学理论，该理论通过 Vlasov 势引入了平均场对齐相互作用。系统假设存在时间尺度分离，即取向弛豫发生得比平移动量弛豫更快。分布函数 $f(x, p, \nu, \varsigma, t)$ 取决于位置、线性动量、微观结构构型（取向 $\nu$ ）以及共轭动量 $\varsigma$ 。
碰撞算子： 为了使分析具有可操作性，完整的 Boltzmann 碰撞算子被单个弛豫项的 Bhatnagar–Gross–Krook (BGK) 算子所取代。该算子保持了完整动力学理论的守恒律和熵产生结构。
平衡律： 作者通过对动力学方程求矩，推导出了质量、线性动量、角动量、能量和熵的宏观平衡方程。一个关键特征是在角动量平衡中识别出了一个平均场转矩项，该项驱动了对齐的涌现。
热力学结构： Helmholtz 自由能被确定为热力学势。作者证明了平衡分布（麦克斯韦分布）使该自由能最小化，并推导了相关的 Gibbs 分布取向轮廓。他们证明了在特定相互作用势（Kuramoto 和 Maier–Saupe 类型）下，系统随着温度降低会经历从各向同性到向列型的相变，且 Vlasov 转矩在对齐平衡态时消失。
Chapman–Enskog 展开： 在局部麦克斯韦平衡态附近进行了一阶 Chapman–Enskog 展开。与标准推导不同，展开中的时间导数是通过零阶平衡律（包含 Ericksen 应力和偶应力项）而非标准的压缩 Euler 方程进行代换的。这确保了取向弹性能够正确地传播到一阶修正项中。
本构闭合： 作者采用了 Rajagopal 和 Srinivasa 的约束极大化原理。他们将通过动力学展开得到的熵产生率与从 Helmholtz 自由能推导出的热力学结构相结合。通过在满足平衡律的条件下极大化熵产生，推导出了 Cauchy 应力、偶应力以及热通量的本构关系。

主要贡献与结果

涌现对齐： 不同于以往假设对齐的无粘性模型，此推导允许通过 Vlasov 转矩使各向同性到向列型的转变自然涌现。分析确认了平均场转矩在对齐 Gibbs 平衡态时消失，从而确保了与宏观平衡律的一致性。
粘性和热效应： 本文推导了一套完整的关于可压缩、粘性、导热及自旋扩散向列型流体的本构方程组。所得应力张量包括：
- 牛顿剪切粘度。
- 源于平移-转动能量不平衡的各向同性体积粘度。
- 与 Ericksen 张量及速度梯度耦合的各向异性粘性应力。
- 包含自旋扩散效应（一种 Eringen 型微极项）及与指向标梯度耦合的偶应力张量。
输运系数： 显式表达了由 BGK 弛豫时间 $\tau$ 和 Frank 模量 $K$ 决定的输运系数（剪切粘度、热导率、体积粘度和自旋粘度）。热导率满足双原子气体的 Eucken 关系。
不可压缩极限： 该框架通过在极大化过程中利用拉格朗日乘子强制执行散度为零的约束，自然地产生了不可压缩对应的模型，而无需单独推导。
结构新颖性： 推导出的粘性应力张量受单一粘度参数控制（由于各向同性 BGK 近似），并且无法恢复完整的五个系数的 Leslie 粘性张量。作者指出，若要恢复完整的 Leslie 张量，需要一个具有多个弛豫时间尺度的各向异性 BGK 算子，这是正在进行中的工作。

意义与主张
本文声称将 Farrell, Russo, 和 Zerbinati [1] 最近的无粘性动力学理论进行了推广，以涵盖粘性、热效应及自旋扩散效应。其主要意义在于直接从动力学理论中提供了一种热力学一致的耗散效应推导方法，弥补了微观粒子相互作用与有序流体宏观连续介质力学之间的鸿沟。

作者强调，他们的方法克服了以往模型不预设对齐的局限性；相反，对齐是通过自由能最小化以及 Vlasore 转矩在平衡态时消失而涌现的。所得系统（方程 8.21）代表了 Leslie–Ericksen 方程的一个可压缩、粘性、导热及自旋扩散变体。在无粘性极限（ $\tau \to 0$ ）下，该模型退化为 [1] 中推导出的无粘性系统。

这项工作被视为迈向更完整的液晶动力学理论的一步，作者谦虚地指出，恢复具有五个独立粘度的完整 Leslie 粘性张量仍是一个开放性问题，需要更复杂的各向异性碰撞算子。