11月4日至5日日冕物质抛射日地传输时间的数值模拟.doc

PAGE PAGE 2 CME日地传输时间的数值模拟 1998年11月4日至5日日冕物质抛射日地传输时间的数值模拟 熊明，郑惠南，汪毓明，傅向荣，王水，窦贤康 中国科技大学地球和空间科学学院, 合肥 230026 E-mail: HYPERLINK mailto:mxiong@ mxiong@, hue@ 摘要 在一维球坐标系下模拟了1998年11月4日至5日三个连续日冕物质抛射(CME)在行星际空间的传播和相互作用并最终形成``复杂抛射的日地传输过程。首先在磁流体力学(MHD)数值模拟中应用Harten总变差减小(TVD)格式，通过调节计算模型中的引力无量纲参数、等离子体参数和气体多方指数，构造出数值计算所需的初态背景，使之在拉格朗日点处L1的太阳风速度、质子数密度及质子热压力与磁压力的比值与ACE卫星的观测数据一致。接着仅采用速度脉冲的扰动形式，其输入的幅度和持续时间由Lasco/C2、GOES、LEAR的观测数据并结合Michalek等提出的CME``锥模型来确定。数值计算结果得到的两个激波到达时间和ACE卫星观测值的时间误差分别是3小时和4小时。这表明该模型能估算续发CME在行星际空间演化后驱动激波的到达时间和大致强度，在空间天气的激波到达时间的预报方面有潜在的应用价值。 关键词：日冕物质抛射，磁流体力学，数值计算，空间天气 A NUMERICAL SIMULATION ON THE SOLAR-TERRESTRIAL TRANSIT TIME OF SUCCESSIVE CMES DURING NOVEMBER 4-5, 1998 Xiong Ming1) Zheng Hui-Nan1) Wang Yu-Ming1) Fu Xiang-Rong1) Wang Shui1) Dou Xian-Kang1) 1) School of Earth and Space Science, University of Science and Technology of China, Heifei 230026, China Abstract The solar-terrestrial transit process of three successive CMEs during November 4-5, 1998 has been investigated numerically in one-dimensional spherical geometry. These CMEs interact with each other while they are propagating in interplanetary space and finally form a Complex Ejecta. A Hartens total variation diminishing (TVD) scheme is applied to solve MHD equations numerically, starting from an ambient solar wind equilibrium, with appropriate dimensionless gravity parameter , plasma beta , and gas polytropic index . The equilibrium is consistent in solar wind speed , proton number density, and the ratio of proton thermal pressure to magnetic pressure with the observation of ACE spacecraft at Lagrange point (L1). Merely velocity pulse is introduced in the numerical computation, whose amplitude and duration are determined by observation data of Lasco/C2, GOES, LEAR combined with CMEs Cone Model proposed by Michalek et al.. The results show that the differences of two shock arrival times (SATs) between the numerical calculation and ACE observation are 3 and 4 hours respec