太阳低层大气小尺度活动数值模拟研究

发布时间：2017-12-15 18:05

  本文关键词：太阳低层大气小尺度活动数值模拟研究

  更多相关文章： 太阳耀斑 磁重联 磁流体 数值模拟

【摘要】：在太阳耀斑、日冕物质抛射(CME)和其他太阳活动现象中,磁重联一直都扮演着至关重要的作用。然而越来越多的作者把注意力集中在了太阳小尺度的活动之上,比如Ellerman炸弹、微耀斑、色球喷流等等。其原因在于这些小尺度活动有着相对简单的结构,比较容易研究。研究这些活动有助于我们理解太阳活动的基本物理机制,并为我们研究大尺度活动现象提供线索和依据。 微耀斑是一种小尺度、持续时间短的太阳活动现象。关于微耀斑的研究已经有几十年的历史。从统计的结果来看,微耀斑的大小从几个角秒到20个角秒不等,持续时间在10-30分钟之间,释放的能量在1025-1029尔格之内。微耀斑在各个波段都有观测,比如Ha、软X射线、硬X射线以及微波辐射等等。但是并不是所有的微耀斑同时在所有波段都有观测。微耀斑的观测特性,比如加热现象、与磁场的联系、持续时间和高能粒子的存在等等,都说明了微耀斑是磁重联的结果,与大耀斑类似。 我们主要使用磁流体(MHD)数值模拟来研究微耀斑的形成和演化。使用的模拟程序主要有两个代码。其一是基于CIP-MOCCT格式的有限差分格式,CIP格式主要处理流体部分,而MOCCT方法主要用于处理磁场部分,保持在演化过程中磁场的散度始终不变。另外一个代码是MAP程序。MAP是基于FORTRAN程序语言用于计算MHD方程组的天体物理程序。其中使用了自适应网络(AMR)的高级算法和完全MPI并行化。MAP程序中集成了三种可选的数值格式,分别为MMC (modified Mac Cormack scheme)、LF (Lax-Friedrichs scheme)和WENO (weighted essentially non-oscillatory scheme)格式。三种格式都是二阶、二步、分量形式(component wise)的用于计算双曲守恒率方程组的数值格式。TVD (total variation diminishing)限制器和近似Riemann解法也包含在内。虽然格式是二阶的,但是AMR算法可以使数值计算在达到尽可能高的分辨率的同时尽量缩短计算时间。除此之外,MAP程序通过使用EGLM (extendedgeneralized Lagrange multiplier)-MHD方程组来减少磁场散度误差。处理其他源项所使用的数值计算方法也有详细的介绍。 在考虑了重力、电离和辐射的情况下,我们使用CIP-MOCCT格式求解2.5维可压缩电阻MHD方程组来研究色球微耀斑现象。模拟中微耀斑的能量主要来自于通过磁重联过程释放的磁场能量。宁静太阳大气VALC模型和10%的氦丰度也考虑到模拟之中。我们的2.5维模拟可以定量的解释微耀斑在色球的温度增加情况。模拟还发现温度的增加量△T主要受到背景磁场强度的影响,而微耀斑的演化时间At与模拟所假设的反常电阻的大小有关系。在一系列模拟中,我们找到一个温度增加与磁场、重联高度和反常电阻值之间的标度率,即AT/At-nH-1.5B2.1η0.88。另外,模拟中向上的喷流速度远大于向下的喷流速度,X重联点的位置也会随模拟的演化有所上升和下降。 虽然用模拟得到的微耀斑可以定量解释色球微耀斑的温度增加现象,但是这个模拟的最大问题还在于其使用了十分简单的磁位形,模拟中并不考虑这种反向平行的磁场位形是如何形成的,我们也只关注了磁重联发生之后的微耀斑演化过程。此外,微耀斑在EUV或者软x射线等高温波段的辐射也不能通过色球磁重联的模拟得到。为了更加自洽的解释微耀斑的形成和演化过程,我们又进行了更加深入的MHD模拟。在新的模拟中,磁重联由新浮磁流和背景磁场之间的相互作用产生。背景磁场则是使用了更加贴近实际的伞盖磁位形,每个伞盖磁场都扎根于超米粒的边界。通过实时地改变底边界条件,新浮磁流从超米粒的中心浮出并与已经存在的伞盖磁场相互作用,从而发生磁重联。我们成功模拟出了日冕起源和色球起源的微耀斑。这两种微耀斑的区别在于其X重联点的高度分别处于日冕中和色球中。具有日冕起源的微耀斑有着较大的尺寸和更高的温度,其热的喷流的温度可以达到-1.8×106K,对应于观测中的EUV/SXR喷流,其冷的喷流的温度为-104K,对应于观测中的Hα/Ca日浪。而具有色球起源的微耀斑由于其重联高度很低,模拟中只有色球层有温度增加现象,对应于Hα/Ca亮点。参数依赖研究表明新浮磁流的尺度和磁场强度是决定重联高度的关键性参数。 最近,色球层中无处不在的喷流被Hinode卫星的SOT望远镜发现,比如色球海葵状喷流、黑子半影微喷流、针状物和黑子本影亮桥喷流等等。一般认为,这些微小的频繁出现的喷流是低层大气磁重联的结果。但是这些喷流的三维磁场位形依然不是很清楚。在这里我们提出一种可用于解释这些喷流现象的可能模型。在太阳低层大气,剪切的磁场之间的磁重联会产生一种扇型喷流,这种喷流具有一个十分有趣的特性,那就是其能沿着引导磁场的方向运动。而通常情况下,引导磁场又是比较垂直于太阳表面的,因此我们的扇型喷流可用于解释观测中发现的色球喷流现象,同时此模拟也能为天体物理其他领域的喷流现象提供了观测依据。
【学位授予单位】：南京大学
【学位级别】：博士
【学位授予年份】：2012
【分类号】：P182

【共引文献】

中国期刊全文数据库 前6条

1 阎坤;天体运行轨道的一般性Binet方程形式[J];地球物理学进展;2005年02期

2 陈寿元;相对静止论[J];山东大学学报(工学版);2002年04期

3 王云飞;;宇宙负能量的存在及作用[J];石家庄职业技术学院学报;2006年02期

4 刘俊红,钮蒸,东艳晖;白矮星的结构特点[J];石家庄师范专科学校学报;2004年03期

5 叶兴浩;真空涨落与相对论量子力学的共同基础[J];绍兴文理学院学报;2004年09期

6 费保俊;孙维瑾;肖昱;季诚响;;X射线脉冲星自主导航的基本测量原理[J];装甲兵工程学院学报;2006年03期

中国博士学位论文全文数据库 前6条

1 吴开谡;天体相关的辐射俘获反应和s—过程网络方程及中子源研究[D];中国原子能科学研究院;2003年

2 郭广海;弯曲时空中的标量场[D];大连理工大学;2005年

3 孙祖尧;暗能量相关问题的研究[D];中国科学院研究生院（上海天文台）;2007年

4 马任意;磁化黑洞吸积盘的X射线辐射[D];华中科技大学;2006年

5 周博颜;活动星系核的喷流进动[D];厦门大学;2007年

6 崔延美;太阳光球磁场特性与耀斑相关性研究[D];中国科学院研究生院（国家天文台）;2007年

中国硕士学位论文全文数据库 前10条

1 许云;γ射线暴火球模型的研究[D];华中师范大学;2003年

2 李磊;星系NGC4565天区附近的变星搜寻[D];辽宁师范大学;2003年

3 胡绍明;Blazar天体的CCD测光研究[D];云南师范大学;2004年

4 徐岩;δScuti型变星脉动特性的研究[D];天津师范大学;2003年

5 俞云伟;带壳奇异星的热演化[D];华中师范大学;2006年

6 李冬生;射电天文学历史成果中的技术作用和理论背景及其思考[D];华东师范大学;2006年

7 尤莉莎;活动星系核的演化及其中心黑洞质量[D];云南师范大学;2006年

8 张曦;数字相关频谱仪基带变换处理器的设计[D];中国科学院研究生院（上海天文台）;2007年

9 张逢春;IRAS红外选Seyfert1和Seyfert1.5型星系的光谱研究[D];河北师范大学;2007年

10 高全归;成熟γ-射线脉冲星对宇宙线正电子的贡献[D];中国科学院研究生院（云南天文台）;2007年

，

本文编号：1292903