日冕物质抛射和冕流结构相互作用的数值研究.pdfVIP

维普资讯 中 ，lj种 (E辑) 第 32卷 第4期 SCIENCEINCHINA (SeriesE) 2O02年 8月 日冕物质抛射和冕流结构相互作用 的 数值研 究 * 叶 占银 魏奉思 冯学尚 石 勇 姚久胜 (中国科学院空间科学与应用研究 中心 ，北京 100080) 摘要 采用具有同心圆形闭磁场结构的二维模型作为 日冕物质抛射 (CME)的触发 模型，并使这种触发模型分别在偏离冕流结构对称中心 10o和45。的位置浮 出，数值 模拟这时产生CME事件的特征 ．模拟结果可 以定性解释 SOHO飞船观测到的部分 CME事件特征 ，模拟结果表明：(1)在两种情况下，二维触发模型的浮 出都可以触 发CME事件 ，和 CME事件相关的磁场结构闭合并呈现非对称性．(2)CME事件 的 闭磁场结构在 向外传播的过程中，将不断向电流片偏转，这种偏转效应主要发生在 几十个太阳半径的范围内，最终 CME事件沿着电流片传播．(3)不同位置处浮出触 发模型触发 的CME事件，将有不 同的磁场位形，在 10o位置浮 出触发模型所触发 CME事件的磁场结构近似呈圆形，而在45。浮出触发模型所触发 CME事件的内部磁 场结构近似呈月牙形． 关键词 太阳风 日冕物质抛射 Mil1)方程组的数值模拟 飞船观测到的 CME现象丰富多彩，其 中有许多 CME现象在冕流结构附近发生或者对冕 流结构产生影响，例如 ，H1lndl1ausen…发现在 YOHKOH观测的 CME事件 中，73％的事件和冕 流结构相伴随；文献 [2]确认了这一结论 ，它研究了1996年 1月至 1998年7月期间LASCO．C2 图像 中的CME事件 ，发现 63％的 CME事件和冕流结构相伴随 ．另一方面 ，统计结果表 明大 扰动事件产生的激波存在向电流片偏转效应_3’4J，也就是说即使发生 CME事件的位置偏离冕 流结构中心 ，在它向外传播的过程 中，也会受冕流结构的影响 ．冕流结构如何导致 CME产生 的激波 向电流片偏转 ，在冕流结构和 CME事件相互作用时，磁场位形发生什么样 的变化 ，研 究这个动力学过程一般采用磁流体力学方程组 (M皿 )进行数值模拟 ． 本文在研究 CME事件和冕流结构相互作用时，采用文献 [5]中具有 同心圆形闭磁场结构 的二维模型作为 CME触发模型 ，研究触发 CME在背景太阳风 中传播的过程 ．背景太阳风取 为满足 MHD方程组的稳态 自洽解 ，在计算这个稳态 自洽解时，采用时间松弛迭代法，开始迭 代计算时磁场和太阳风参数采用文献[6]的值，这样做可以使触发 CME事件模型的中心压力 降低 ． 2001-04-23收稿 ，2001-09．19收修改稿 *国家 自然科学基金 (批准号49925412，40lO40o8)和国家 “九七三 ’’重点基础研究 (G200o78405)资助项 目 维普资讯 第 4期 叶占银等：日冕物质抛射和冕流结构相互作用 的数值研究 5l1 1 控制方程和计算方法 太阳风等离子体的动力学过程普遍地采用 MHD方程组来描述 ，为了方便数值模拟，可 以 作轴对称假设 ，采用球坐标系下二维 MHD控制方程组_6J，并采用二阶精度 的 MacCormack11 差分格式对它进行离散 ，离散格式加人 Lapidus人工黏性项，计算域为 rX =(1，20)X(一1， 181)，网格划分为 200X200，内边界采用相容投影特征线边界条件 ，外边界处太阳风是超快、 慢波速和超 Alfv6n波速流动 ，采用外推边界条件 ．极 区采用延拓边界条件 ． 2 伽 的触发模型 研究 CME和冕流间相互作用，需要根据产生 CME物理机制来选择 CME的触发模型．虽 然 目前对产生 CME的物理机制并未完全研究清楚，但是 比较普遍 的观点是 当 日冕大