跳到主要内容
Sandbox Physics
从磁应力走向磁重联

L2 · 磁重联与太阳耀斑信号

薄电流片怎样决定交换磁通的速度

一张稳态薄电流片,吸入多少质量,两端就要流出多少。把这个几何要求与磁扩散结合,便能推得Sweet–Parker估计。片长、片厚与速度的关系,也会告诉我们哪些情形已经超出这套估计。

先画矩形,说清两种半尺度

让 xx 沿喷流方向,yy 横穿电流片。半片长为 ℓ\ell,半片厚为 δ\delta。两侧大小为 BuB_u 的磁场沿片指向相反方向;用正值 vinv_{\rm in}、voutv_{\rm out} 表示入流和出流速度。

假设薄层近似矩形、保持稳态,进出密度相同,也没有质量在里面积起来。按单位深度计,两张长面吸入气体,两端短面把气体排走。面数与完整边长中的二都能约掉。

出口窄,出流就要比入流快,才能排掉同样多的气体。到这一步,我们只用了几何和质量守恒,还没有用到电阻。

4ρℓvin=4ρδvout⟹vinℓ=voutδ4\rho\ell v_{\rm in}=4\rho\delta v_{\rm out}\quad\Longrightarrow\quad v_{\rm in}\ell=v_{\rm out}\delta
两边取相同密度与不变方向深度。若气体明显压缩,这个关系也要相应改变。

来源:Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4

让层内和层外的电场接上

二维磁场保持稳态时,垂直平面的电场在薄层内外有相同数值。上游主要由流体带着磁场移动,电场大小约为入流速度乘场强;靠近片心,则由电阻与电流给出。

电流密度大小约为 Bu/(μ0δ)B_u/(\mu_0\delta)。乘上电阻率,再用 ηm=ηe/μ0\eta_{\rm m}=\eta_{\rm e}/\mu_0 换成磁扩散率,两边的场强就约掉了。所得关系把横穿薄层的入流与磁扩散配在一起。

这里比的是大小。实际电场的正负,要由磁场方向和坐标方向决定;电流的符号不能在计算途中被热功率替代。

∣Ez∣∼vinBu∼ηe∣Jz∣∼ηmBuδ|E_z|\sim v_{\rm in}B_u\sim\eta_{\rm e}|J_z|\sim\eta_{\rm m}\frac{B_u}{\delta}
这一步采用稳态薄层和均匀标量磁扩散率。
vin∼ηmδv_{\rm in}\sim\frac{\eta_{\rm m}}{\delta}
在其余条件保持时,更薄的扩散区能够对应更快的入流。

来源:Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4 · Zhou et al. (2019), Magnetic island merger as a mechanism for inverse magnetic energy transfer, Physical Review Research 1, 012004(R)

磁张力给出出流的速度尺度

松弛的磁场沿薄层拉动气体。比较上游单位体积磁能和流体动能,可得到出流与阿尔芬速度同量级。压力、几何与耗散分配还会影响数值系数,这个量级推导没有把系数定死。

把出流尺度代回质量守恒,再与扩散关系联立,片厚和入流速度就一起出现。定义Lundquist数时继续使用刚才的半片长,便不会在长度约定上丢掉因子。

vout∼vA=Buμ0ρ,Sℓ=ℓvAηmv_{\rm out}\sim v_A=\frac{B_u}{\sqrt{\mu_0\rho}},\qquad S_\ell=\frac{\ell v_A}{\eta_{\rm m}}
这个Lundquist数比较同一半片长上的扩散时间与阿尔芬传播时间。
δℓ∼Sℓ−1/2,vinvA∼Sℓ−1/2\frac{\delta}{\ell}\sim S_\ell^{-1/2},\qquad\frac{v_{\rm in}}{v_A}\sim S_\ell^{-1/2}
这是稳态、薄片、弱黏性条件下的Sweet–Parker尺度关系。

来源:Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4

算一个数,再回头检查假设

给假想薄层指定 Sℓ=104S_\ell=10^4,半片厚便约为半片长的百分之一,入流约为阿尔芬速度的百分之一。算出来确实很薄,和推导的起点相容。

若Lundquist数再增加一百倍,两个比例各缩小十倍。反过来,把它降到一的量级,片厚就会接近片长,原先的薄片近似也就站不住了。

黏性还引入磁Prandtl数 Pm=ν/ηmP_{\rm m}=\nu/\eta_{\rm m},其中 ν\nu 为运动黏度。强黏性可以压低出流速度。交互算例会求出相应的流动和耗散热,不能拿它直接拟合弱黏性稳态薄片的通用系数。

指定的估计片厚比例入流比例
Sℓ=104S_\ell=10^4δ/ℓ∼0.01\delta/\ell\sim0.01vin/vA∼0.01v_{\rm in}/v_A\sim0.01
Sℓ=106S_\ell=10^6δ/ℓ∼0.001\delta/\ell\sim0.001vin/vA∼0.001v_{\rm in}/v_A\sim0.001

来源:Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4 · Zhou et al. (2019), Magnetic island merger as a mechanism for inverse magnetic energy transfer, Physical Review Research 1, 012004(R)

薄片本身稳不稳定,要另问一次

磁扩散率降低时,稳态估计中的薄片会越来越薄。这样的平衡本身可能不稳定。“扰动增长多快”和“稳态入流多快”是两个问题,可以出现不同趋势。

Loureiro等人采用完整片长 LCS=2ℓL_{\rm CS}=2\ell,因而有 Sfull=2SℓS_{\rm full}=2S_\ell。对应的薄度参数写在下面。把长度约定带在身边,读原文时就不会漏掉根号二。

进阶指南将重算那篇论文的线性色散关系,回答小扰动怎样增长。要得到非线性的磁岛链与并合过程,还须充分分辨的演化计算。这里的周期初始应力模型则给出一个暂态能量转换问题,并保留自己的网格检查与停止时刻。

ϵ=2δLCS∼2Sfull,Γ0=2vALCS\epsilon=\frac{2\delta}{L_{\rm CS}}\sim\sqrt{\frac{2}{S_{\rm full}}},\qquad\Gamma_0=\frac{2v_A}{L_{\rm CS}}
固定2007v1阅读与参考脚本采用这套完整片长和流动率约定。

来源:Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4 · Zhou et al. (2019), Magnetic island merger as a mechanism for inverse magnetic energy transfer, Physical Review Research 1, 012004(R)

自己推一步

保持半片长和上游阿尔芬速度不变,把磁扩散率增加九倍,Sweet–Parker估计会怎样变?

Lundquist数变为九分之一,δ/ℓ\delta/\ell 与 vin/vAv_{\rm in}/v_A 都变为三倍。还要检查片仍然够薄,以及其余稳态假设是否成立。

为什么更大的Lundquist数,能对应更慢的稳态入流,也能对应更快的不稳定性增长?

稳态入流来自质量流和扩散平衡;线性不稳定性研究的是围绕那张平衡薄片的小扰动,而薄片厚度也随Lundquist数改变。两者问的物理量不同。

能否用一个周期暂态快照,测出通用的Sweet–Parker系数?

周期快照可以用来提出暂态比较问题。估计稳态系数,还须建立矩形薄层、固定上游区和相应的弱黏性出流条件。

比较磁扩散与黏性阻尼

  1. 先预测:提高磁扩散率和提高黏性,会直接改变同一条方程吗?
  2. 在中间保存时刻比较两个算例,并把探针拖过电流集中区。
  3. 借助能量曲线,区分磁场被抹平与流体运动被抑制。

原始来源

  1. Zhou et al. (2019), Magnetic island merger as a mechanism for inverse magnetic energy transfer, Physical Review Research 1, 012004(R)
  2. Loureiro, Schekochihin & Cowley (2007), astro-ph/0703631v1, equations 3 and 18–20, Figure 4
反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策