拉格朗日描述法简介

在流体力学中，拉格朗日描述法（Lagrangian description）是一种随流体质点运动而记录其状态的方式。与之对应的是 欧拉描述法（Eulerian description），后者在固定的空间坐标系中观察流体的瞬时分布。拉格朗日描述的核心思想是把每一个流体质点视为“观察对象”，用其随时间的位移、速度、密度等随体变量来刻画流动。

拉格朗日描述的优势在于：

• 能够自然地处理随体变形的边界（如自由表面、固体-流体界面）。
• 在求解粒子追踪、混合过程以及相场模型时更为直观。
• 对于不可压缩流体的质量守恒，可直接使用质点的体积保持性质。

拉格朗日坐标与欧拉坐标的对比

注意：在实际计算中，拉格朗日坐标 (\mathbf{X}) 与欧拉坐标 (\mathbf{x}) 之间的映射 (\mathbf{x}=\chi(\mathbf{X},t)) 必须是可逆且光滑的，否则会出现交叉或奇异。

拉格朗日方程的推导

1. 运动方程的拉格朗日形式

对单个无粘性流体质点，受力满足牛顿第二定律：

[
\rho_0 \frac{d^2 \mathbf{x}}{dt^2}= \mathbf{f}(\mathbf{x},t)+\nabla_{\mathbf{X}}\cdot\mathbf{P},
]

其中：

• (\rho_0) 为质点的参考密度（常为常数）。
• (\mathbf{f}) 为体力（如重力）。
• (\mathbf{P}) 为第一Piola‑Kirchhoff应力张量，描述内部弹性或粘性作用。

2. 连续方程（质量守恒）

在拉格朗日框架下，质量守恒可写成：

[
\rho(\mathbf{X},t) J(\mathbf{X},t) = \rho_0(\mathbf{X}),
]

其中

[
J(\mathbf{X},t)=\det!\bigg(\frac{\partial \mathbf{x}}{\partial \mathbf{X}}\bigg)
]

是 变形梯度 (\mathbf{F}=\partial \mathbf{x}/\partial \mathbf{X}) 的行列式。对不可压缩流体，(J\equiv1)。

3. 动量方程（拉格朗日形式）

将欧拉动量方程

[
\rho\bigg(\frac{\partial \mathbf{u}}{\partial t}+(\mathbf{u}\cdot\nabla)\mathbf{u}\bigg)=\nabla\cdot\boldsymbol{\sigma}+\mathbf{f}
]

通过坐标变换，可得到拉格朗日形式：

[
\rho_0 \frac{d^2 \mathbf{x}}{dt^2}= \nabla_{\mathbf{X}}\cdot\big(\mathbf{P}\big)+\rho_0\mathbf{f}.
]

这里的导数 (\displaystyle \frac{d}{dt}=\frac{\partial}{\partial t}\bigg|_{\mathbf{X}}) 为 全导数，即随质点运动的时间导数。

示例：简单流动的拉格朗日追踪

考虑二维平面内的均匀旋转流，速度场（欧拉描述）为：

[
\mathbf{u}(x,y)=\omega\begin{pmatrix}-y\x\end{pmatrix},
]

其中 (\omega) 为常数角速度。我们希望在拉格朗日框架下求出任意质点的轨迹 (\mathbf{x}(t))。

1. 建立常微分方程

[
\frac{d\mathbf{x}}{dt}= \omega\begin{pmatrix}-y\x\end{pmatrix},
\qquad \mathbf{x}(0)=\mathbf{X}=(X,Y).
]

写成分量形式：

[
\begin{cases}
\displaystyle \frac{dx}{dt}= -\omega y,\[4pt]
\displaystyle \frac{dy}{dt}= ;;\omega x.
\end{cases}
]

2. 求解

对第一个方程求二阶导数：

[
\frac{d^2x}{dt^2}= -\omega\frac{dy}{dt}= -\omega^2 x,
]

得到简谐振子方程，其通解为

[
x(t)=X\cos(\omega t)-Y\sin(\omega t).
]

同理，

[
y(t)=X\sin(\omega t)+Y\cos(\omega t).
]

3. 结果解释

• 质点在原点周围做半径为 (\sqrt{X^2+Y^2}) 的匀速圆周运动。
• 拉格朗日映射 (\chi(\mathbf{X},t)) 为旋转矩阵：

[
\chi(\mathbf{X},t)=
\begin{pmatrix}
\cos\omega t & -\sin\omega t\
\sin\omega t & ;\cos\omega t
\end{pmatrix}
\mathbf{X}.
]

此例展示了拉格朗日描述如何直接给出质点轨迹，而不必在固定空间网格上求解偏微分方程。

数值实现要点

在实际 CFD 中，拉格朗子方法常以 粒子方法（如 SPH、DEM）或 网格移动法（ALE）实现。下面列出实现时的关键步骤：

1. 初始化粒子集合

   • 依据几何域划分粒子位置 (\mathbf{X}_i)。
   • 为每个粒子分配质量 (m_i=\rho_0 V_i)（(V_i) 为参考体积）。

2. 计算变形梯度 (\mathbf{F}_i)

   • 常用最近邻粒子加权最小二乘法：
     [
     \mathbf{F}i=\sum{j} w_{ij}(\mathbf{x}_j-\mathbf{x}_i)\otimes(\mathbf{X}_j-\mathbf{X}i)^{-1},
     ]
     其中 (w{ij}) 为核函数权重。

3. 求解应力张量 (\mathbf{P}_i)

   • 对于线性弹性流体：(\mathbf{P}= \lambda , \text{tr}(\mathbf{E})\mathbf{I}+2\mu\mathbf{E})，(\mathbf{E}=\frac{1}{2}(\mathbf{F}^T\mathbf{F}-\mathbf{I}))。
   • 对于粘性流体，可加入黏性项 (\mathbf{P}= -p\mathbf{I}+ \mu(\nabla\mathbf{u}+\nabla\mathbf{u}^T))。

4. 时间积分

   • 常用显式的 Verlet 或 Runge‑Kutta 方法：
     [
     \mathbf{x}_i^{n+1}= \mathbf{x}_i^{n}+ \Delta t,\mathbf{v}_i^{n}+ \frac{\Delta t^2}{2m_i}\mathbf{F}_i^{n},
     ]
     [
     \mathbf{v}_i^{n+1}= \mathbf{v}_i^{n}+ \frac{\Delta t}{2m_i}\big(\mathbf{F}_i^{n}+\mathbf{F}_i^{n+1}\big).
     ]

5. 边界处理

   • 对于固体壁面，可采用 镜像粒子 或 反弹条件。
   • 对自由表面，常使用 粒子重采样（增删粒子）保持分辨率。

6. 后处理

   • 将拉格朗日数据映射回欧拉网格（如体积加权平均）以便可视化。
   • 计算宏观量（如动能、动量通量）时，需要考虑 Jacobian (J) 的影响。

小结

拉格朗日描述法通过随体坐标 (\mathbf{X}) 追踪流体质点的运动，为处理变形边界、颗粒混合以及自由表面提供了天然的数学框架。其核心方程包括：

• 变形梯度 (\mathbf{F}=\partial\mathbf{x}/\partial\mathbf{X}) 与 Jacobian (J) 的几何意义；
• 拉格朗日质量守恒 (\rho J=\rho_0)；
• 拉格朗日动量方程 (\rho_0 d^2\mathbf{x}/dt^2 = \nabla_{\mathbf{X}}\cdot\mathbf{P} + \rho_0\mathbf{f})。

通过上述示例可以看出，拉格朗日方法能够直接给出质点的轨迹，而不必在固定空间网格上求解偏微分方程。数值实现时，粒子方法和 ALE（Arbitrary Lagrangian‑Eulerian）技术是两大主流手段，关键在于准确计算变形梯度、应力以及保持数值稳定性。

掌握拉格朗日描述的理论与实现技巧，对从事高精度流体模拟、自由表面预测以及多相流研究的工程技术人员具有重要意义。