[论文解读] GEM -- An Energy Conserving Electromagnetic Gyrofluid Model

Bruce Scott|arXiv (Cornell University)|Jan 23, 2005

Ionosphere and magnetosphere dynamics被引用 6

一句话总结

本文提出GEM模型，该模型通过严格耦合静电势、密度、温度和平行通量，利用极化方程实现总涨落自由能的守恒。该模型在所有尺度上确保能量守恒，从而实现对离子和电子回旋半径尺度下双尺度湍流的稳定模拟，其结果与gyrokinetic基准在ITG和ETG湍流中高度一致。

ABSTRACT

The details of fluctuation free energy conservation in the gyrofluid model are examined. The polarisation equation relates ExB flow and eddy energy to combinations of the potential and the density and perpendicular temperature. These determine the combinations which must appear under derivatives in the moment equations so that not only thermal free energy but its combination with the ExB energy is properly conserved by parallel and perpendicular compressional effects. The resulting system exhibits the same qualitative energy transfer properties as corresponding Braginskii or Landau fluid models. One clear result is that the numerical model built on these equations is well behaved for arbitrarily large perpendicular wavenumber, allowing exploration of two scale phenomena linking dynamics at the ion and electron gyroradii. When the numerical formulation is done in the globally consistent flux tube model, the results with adiabatic electrons are consistent with the ``Cyclone Base Case'' results of gyrokinetic models.

研究动机与目标

• 通过确保涨落自由能守恒，为gyrofluid模型建立严格的能量基础。
• 解决以往gyrofluid模型中能量守恒被忽略或近似处理所导致的不一致性。
• 实现涉及离子和电子回旋半径的双尺度湍流的长期、稳定模拟。
• 在核心湍流（包括Cyclone Base Case）中与gyrokinetic结果保持一致。
• 支持局部和全局模拟，具备几何与能量一致性。

提出的方法

• 推导出一个极化方程，将ExB流、涡旋能量与静电势、密度和垂直温度的组合联系起来。
• 构建包含对特定变量组合求导的矩方程，以保持总自由能守恒。
• 在全局一致的通量管几何中实现该模型，以维持几何保真度。
• 采用六矩形式，包括离子和电子密度、温度及平行热通量。
• 应用人工耗散（扩散和超扩散）以稳定模拟，同时不破坏能量守恒。
• 通过Cyclone Base Case和ETG湍流基准，与gyrokinetic模拟结果进行验证。

实验结果

研究问题

• RQ1 如何构建一个gyrofluid模型，使其精确守恒总涨落自由能，包括ExB运动和热涨落的贡献？
• RQ2 在存在压缩性和剪切效应的情况下，空间导数下必须出现哪些流体变量的特定组合，才能确保能量守恒？
• RQ3 该模型能否以高精度和计算效率重现ITG和ETG湍流的关键gyrokinetic结果？
• RQ4 在任意大的垂直波数下，该模型的行为如何，特别是在双尺度湍流区域？
• RQ5 在无显式耗散的情况下，极化方程在连接热能与ExB分量间能量传递中起什么作用？

主要发现

• GEM模型精确守恒总涨落自由能，包括热能、ExB能和平行压缩能的贡献。
• 对于ITG湍流，该模型在全谱范围内与Cyclone Base Case结果的偏差在20%以内，非线性阈值与统计误差范围内一致。
• 对于ETG湍流，该模型成功捕捉到径向扩展的结构，输运水平为χₑ ≈ 6ρₑ²Ve/LT，峰值出现在ky = 10，与gyrokinetic模拟结果一致。
• 即使在大垂直波数（kyρi > 10）下，该模型仍保持数值稳定且行为良好，能够解析离子和电子尺度的动力学。
• 电子矩振幅（ne, Te∥）达到约0.2，而离子矩振幅小两个数量级，证实ETG占主导。
• 在ETG模拟中，该模型成功捕捉到t ≈ 8时的非线性转变和t ≈ 12时的饱和，且持续存在径向扩展的结构。