导向场对磁重联影响机制.pdfVIP

导向场对磁重联的影响机制 王晓刚 王嘉琦 杨恒磊 （大连理工大学，物理系） 摘要：在有导向场存在的情况下，采用显格式对磁场重联进行二维三分量的 MHD 数 值模拟。并分别讨论了自发重联和受迫重联两种情况下导向场的存在对重联的影响。 在导向场很弱时，剪切长度的改变对重联影响甚微，当导向场增强，剪切长度的改变 对重联的数值影响也显著增加。并定性的描绘了空间的密度分布。 关键词：导向场，剪切长度，磁场重联 一、 引 言 TOKAMAK 等离子体中的磁重联是典型的有导向场的二维磁重联问题。在空间等离子体中，三 维磁重联在某些情况下也可以归结成有有导向场的二维磁重联。 当两块带有反相磁场的等离子体相互靠近时，电流片两侧的反相磁场在某种机制的触发下会发 生重联。但在大多数实际情况中，电流片两侧的磁场并不会严格是反平行的，这时会出现有角度的 剪切导向场。这个导向场会在剪切长度内迅速衰减为零，以保证靠近电流片的磁场为反平行。也就 是说，非反平行的磁场会在剪切长度内旋转为反平行。 本文给出了不同导向场强度及剪切长度对磁场重联的定性分析。采用显格式对撕裂模不稳定性 和边界扰动触发重联进行了二维三分量的 MHD 数值模拟。主要方程和电阻撕裂模的非线性演化和 边界扰动及其边界条件的模型建立在第二部分介绍；模型的数值计算结果及图形诊断与分析将在第 三部分介绍；结论部分将在最后介绍。 二、 数 值 模 型 [1] （1）非线性撕裂模不稳定性 ： 基本公式： ∂ρ ρ ρ ， （1） =−u ⋅∇ − ∇⋅u ∂t 2 ∂u 1 B 1 2 ρ −ρu ⋅∇u −∇p − ∇ + B ⋅∇B +υ u, （2 ） ∇ ∂t 2 μ0 μ0 ∂B [η( )] ， （3 ） =−u ⋅∇B +B ⋅∇u −B∇⋅u −∇× ∇×B ∂t ∂p =−u ⋅∇p − p ∇⋅u γ ， （4 ） ∂t ∂ ∂ ∂2 ∂2 以上方程组中的微分算符∇ ∇2 ∂x ∂z ∂x 2 ∂z 2 η 及