2025-11-22T00:37:16.330882

Structure-preserving finite element approximations of a hybrid relativistic cold fluid-particle model

Mukhamet, Kormann

We derive mixed finite element discretizations of a cold relativistics fluid model from approximations of the Poisson bracket that preserve mass, energy and the divergence constraints. For time-discretization we derive an implicit energy-conserving average-vector field method or apply an explicit strong-stability preserving Runge-Kutta scheme. We also consider a coupling of the fluid model to relativistic particles. We perform a numerical study of the scheme which shows convergence and conservation properties of the proposed methods and apply the new scheme to a plasma wake field simulation.

academic

Structure-preserving finite element approximations of a hybrid relativistic cold fluid-particle model

基本信息

论文ID: 2510.11500
标题: Structure-preserving finite element approximations of a hybrid relativistic cold fluid-particle model
作者: Tileuzhan Mukhamet, Katharina Kormann (Ruhr-University Bochum, Germany)
分类: math.NA cs.NA physics.comp-ph
发表时间: October 14, 2025 (arXiv预印本)
论文链接: https://arxiv.org/abs/2510.11500

摘要

本文从Poisson括号的近似出发，推导了相对论性冷流体模型的混合有限元离散化方法，该方法能够保持质量、能量和散度约束。在时间离散化方面，作者推导了隐式能量守恒的平均向量场方法，并应用了显式强稳定性保持Runge-Kutta格式。文章还考虑了流体模型与相对论粒子的耦合。通过数值研究验证了所提出方法的收敛性和守恒性质，并将新格式应用于等离子体尾场模拟。

研究背景与动机

要解决的问题：本文旨在构建相对论性冷等离子体流体-粒子混合模型的结构保持数值格式，该模型在激光-等离子体相互作用和粒子加速研究中具有重要应用。
问题的重要性：
- 冷流体模型为等离子体提供了集体的流体描述，在可忽略微观动力学效应时非常有用
- 在研究色散关系、波现象和平衡态附近的等离子体动力学中发挥关键作用
- 作为昂贵的粒子模拟方法的更经济替代方案
现有方法的局限性：
- 相对论方程具有更强的非线性，使得在完全离散级别上的结构保持更具挑战性
- 之前的工作主要集中在非相对论线性化冷模型上
- 缺乏对混合流体-粒子模型(1-9)的结构保持离散化研究
研究动机：构建能够在空间和时间离散级别上精确保持总质量、总能量、弱Gauss定律和divB约束的数值方法。

核心贡献

推导了基于Poisson结构的混合有限元空间离散化方法，提供了两种不同的基函数选择
开发了隐式能量守恒的平均向量场(AVF)时间积分格式，在完全离散级别保持所有原始不变量
分析了显式强稳定性保持Runge-Kutta方法的守恒性质，并提出了Gauss定律清洗程序
实现了流体和粒子的一致耦合处理，保持了混合系统的物理守恒律
通过数值实验验证了方法的收敛性和守恒性，并成功应用于等离子体尾场模拟

方法详解

任务定义

本文研究相对论性冷流体-粒子混合模型的数值求解，该模型包含：

输入：初始流体密度ρ、流体动量M、电场E、磁场B、粒子位置X_k、粒子动量U_k
输出：时间演化的物理场和粒子状态
约束条件：质量守恒、能量守恒、Gauss定律、divB=0约束

模型架构

1. 连续模型

控制方程组包括：

流体密度和动量方程(1-2)： $\partial_t\rho = -\nabla \cdot \left(\frac{M}{\gamma(mM/\rho)}\right)$ $\partial_tM = -\nabla \cdot \left(\frac{M \otimes M}{\rho\gamma(mM/\rho)}\right) + \rho\frac{e}{m}\left[E + \frac{M}{\rho c \gamma(mM/\rho)} \times B\right]$
粒子运动方程(4-5)
Maxwell方程组(6-9)

其中相对论因子： $\gamma(u) = \sqrt{1 + \frac{u \cdot u}{m^2c^2}}$

2. Hamiltonian结构

系统具有Poisson括号结构，Hamiltonian为总能量： $H = \int_\Omega \rho(\gamma(mM/\rho)-1)c^2dx + \sum_k w_k(\gamma(U_k)-1)mc^2 + \frac{1}{8\pi}\int_\Omega (E^2+B^2)dx$

3. 空间离散化

提供两种混合有限元格式：

格式一（无通量形式）：

$\rho_h \in Q_{k+1}$ , $M_h \in Q_{k+1}$ , $E_h \in \tilde{N}_k$ , $B_h \in \tilde{RT}_k$

格式二（含通量形式）：

$\rho_h \in DG_{k+1}$ , $M_h \in RT_k$ , $E_h \in \tilde{N}_k$ , $B_h \in \tilde{RT}_k$
需要处理单元间的数值通量

4. 时间离散化

隐式AVF方法：使用平均向量场梯度处理非线性项： $\frac{u^{k+1}-u^k}{\Delta t} = J(u^k,u^{k+1}) \cdot \int_0^1 \frac{\partial H}{\partial u}(u^\xi) d\xi$

显式SSP-RK方法：三阶强稳定性保持Runge-Kutta格式，配合Gauss定律清洗程序。

技术创新点

Poisson结构保持：通过离散Poisson括号确保反对称性，从而保持能量守恒
非线性处理：使用AVF梯度处理相对论非线性项，保持完全离散级别的守恒性
混合耦合：统一处理流体和粒子贡献，保持电荷守恒和Gauss定律
结构保持投影：使用L²投影确保离散解在正确的函数空间中

实验设置

数据集

使用制造解进行收敛性测试：

三维域 $\Omega = [-1,1]^3$
均匀网格细化： $h = 2^{-i}h_0$ ， $i = 2,3,4,5$
有限元次数： $k = 0,1,2$

评价指标

收敛性：L²范数误差
守恒性：
- 质量误差： $\frac{\text{mass}(t)-\text{mass}(0)}{\text{mass}(0)}$
- 能量误差： $\frac{H(t)-H(0)}{H(0)}$
- Gauss定律误差： $\max_i \int_\Omega [\nabla \cdot E - 4\pi e(\rho/m + \sum_k w_k\delta(x-X_k))] \phi_i dx$
- divB约束误差： $\|\nabla \cdot B_h\|_{L^2}$