[论文解读] A General Approach for Producing Hamiltonian Numerical Schemes for Fluid Equations

Colin J. Cotter|ArXiv.org|Jan 26, 2005

Numerical methods for differential equations参考文献 34被引用 6

一句话总结

本论文提出一种通用的粒子-网格方法，通过结合拉格朗日粒子动力学与欧拉网格基哈密顿量评估，为流体方程构建半离散哈密顿格式。该方法确保了辛时间积分所需的典范结构，在长时间模拟中保持哈密顿量在 O(Δt) 范围内，已在二维 EP-Diff 方程上得到验证，展现出优异的守恒性质。

ABSTRACT

Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical Hamiltonian system; this system may then be integrated in time using a symplectic integrator. The method is Lagrangian, with the variables being a set of Lagrangian particle positions and their associated momenta. The canonical equations obtained yield a discrete form of Euler-Poincaré equations for $l$ when projected onto the grid, with a new form of discrete calculus to represent the gradient and divergence operators. Practical symplectic time integrators are suggested for a large family of equations which include the shallow-water equations, the EP-Diff equations and the 3D compressible Euler equations, and we illustrate the technique by showing results from a numerical experiment for the EP-Diff equations.

研究动机与目标

• 开发一种系统化方法，用于构建流体方程的典范哈密顿半离散化格式，保证其具备辛结构。
• 通过使用拉格朗日粒子恢复典范形式，克服在欧拉坐标系下对非典范李-泊松结构进行离散化的困难。
• 通过结合基于粒子的动力学与基于网格的哈密顿量评估，实现在流体模拟中长期保持哈密顿量。
• 提供一个适用于广泛流体模型类别的框架，包括浅水方程、EP-Diff 方程和可压缩欧拉方程。
• 通过二维 EP-Diff 方程的数值实验，证明该方法在计算上的可行性与稳定性。

提出的方法

• 使用拉格朗日粒子位置与动量作为典范变量，确保离散系统中具有典范辛结构。
• 将基于粒子的动量与密度插值到欧拉网格上，以网格点上的求和形式评估哈密顿量。
• 采用离散微分形式框架，表示离散欧拉-泊因卡雷方程中的梯度与散度算子。
• 通过约化拉格朗日量的勒让德变换推导离散哈密顿系统，并利用粒子-网格插值将其投影到网格上。
• 在半离散系统上应用辛时间积分格式（如辛欧拉法），以在长时间模拟中保持哈密顿量。
• 确保离散方程满足离散欧拉-泊因卡雷方程的形式，从而保持变分结构。

实验结果

研究问题

• RQ1 能否开发一种通用方法，用于在欧拉坐标系下对具有约化拉格朗日量的流体方程生成典范哈密顿半离散化格式？
• RQ2 如何将欧拉坐标系下流体方程的非典范结构转化为适合辛时间积分的典范形式？
• RQ3 粒子-网格插值在保持哈密顿结构的同时，如何实现基于网格的哈密顿量评估？
• RQ4 辛时间积分格式在长时间半离散系统中对哈密顿量的保持效果如何？
• RQ5 该框架能否应用于多种流体模型，如浅水-α 方程与三维可压缩欧拉方程？

主要发现

• 在长时间积分中，离散哈密顿量始终保持在初始值的 O(Δt) 范围内，证实了辛时间积分的有效性。
• 对二维 EP-Diff 方程的数值实验表明，解行为良好，无发散或非物理振荡。
• 该方法允许在 H¹ 空间中处理弱解，从而能够处理如 δ 函数等奇异动量分布。
• 粒子-网格格式支持使用有限元方法处理速度场，实现稳健的空间离散化。
• 该方法使得在欧拉形式中包含复杂微分算子成为可能，而不会使拉格朗日变量相互纠缠。
• 该框架可扩展至包含各种平流物理量（如标量、向量与张量）的半直积系统。