中文核心期刊

中国高校优秀科技期刊

中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

外日球层激波事件的一维磁流体力学数值模拟

郭孝城 周昱成 王赤 李晖

郭孝城, 周昱成, 王赤, 李晖. 外日球层激波事件的一维磁流体力学数值模拟[J]. 深空探测学报(中英文), 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
引用本文: 郭孝城, 周昱成, 王赤, 李晖. 外日球层激波事件的一维磁流体力学数值模拟[J]. 深空探测学报(中英文), 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
GUO Xiaocheng, ZHOU Yucheng, WANG Chi, LI Hui. One-Dimentional Numerical MHD Simulation of the Shock Events in the Outer Heliosphere[J]. Journal of Deep Space Exploration, 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
Citation: GUO Xiaocheng, ZHOU Yucheng, WANG Chi, LI Hui. One-Dimentional Numerical MHD Simulation of the Shock Events in the Outer Heliosphere[J]. Journal of Deep Space Exploration, 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059

外日球层激波事件的一维磁流体力学数值模拟

doi: 10.15982/j.issn.2096-9287.2020.20200059
基金项目: 中国科学院空间科学战略先导科技专项资助项目(XDB41000000,XDA15052500);国家自然科学基金资助项目(41874171,41674146,41574159,41731070);中国科学院前沿科学重点研究计划资助项目(QYZDJ-SSW-JSC028);民用航天技术预先研究资助项目(D030202,D020301)
详细信息
    作者简介:

    郭孝城(1978– ),男,研究员,主要研究方向:日球层物理、磁层物理。 通讯地址:北京市海淀区中关村南二条1号中国科学院国家空间科学中心(100190)电话:(010)62586362 E-mail:xcguo@swl.ac.cn

  • ● It is the first time to use the solar wind data from OMNI,the spacecraft STEREO A and B orbiting the Sun conjunctly to study the evolution and propagation of the solar wind in the outer heliosphere. ● The numerical simulations have been carried out to verify the possible correspondence between the large-scale pressure pulse structure in the heliosphere and the shock event in interstellar space.
  • 中图分类号: V467.3

One-Dimentional Numerical MHD Simulation of the Shock Events in the Outer Heliosphere

  • 摘要: “旅行者1号”在2012年8月穿越日球层顶后连续探测到了星际磁场的局地突然变化,通常被认为对应于激波事件,且为太阳风事件向外传播后形成的大尺度结构与日球层顶相互作用的产物。但由于观测的限制,日球层内的太阳风大尺度结构与星际空间的激波事件的对应关系尚未明确。通过一个包含太阳风等离子体和星际中性原子的双流体磁流体力学数值模型检验了这种对应关系。在距日1 AU的内边界输入2010—2017年分别位于不同经度上的OMNI,STEREO A和STEREO B的太阳风观测数据,研究了这期间太阳风与星际中性原子电荷交换后在外日球层区的传播和演化,并与“旅行者1号”和“旅行者2号”的观测数据做对比。结果表明,“旅行者1号”探测的星际空间激波事件与“旅行者2号”探测到的日球层鞘区内的压力脉冲事件具有明显的关联。
    Highlights
    ● It is the first time to use the solar wind data from OMNI,the spacecraft STEREO A and B orbiting the Sun conjunctly to study the evolution and propagation of the solar wind in the outer heliosphere. ● The numerical simulations have been carried out to verify the possible correspondence between the large-scale pressure pulse structure in the heliosphere and the shock event in interstellar space.
  • 图  1  沿着径向方向的太阳风参数

    Fig.  1  The profiles of solar wind plasma along radial direction

    图  2  2011—2017年采用3种不同内边界数据源的太阳风动压对比

    Fig.  2  The comparision of solar wind dynamical pressure between the observation(red)and the simulations(blue)during the years 2011—2017

    图  3  2012—2016年用3种不同内边界数据源,模拟得到的“旅行者1号”可能探测的动压数据

    Fig.  3  The simulated dynamical pressures that could be encountered by Voyager 1 during the years 2012—2016

  • [1] PARKER E N. The stellar wind regions[J]. The Astrophysical Journal,1961,134:20-27. doi:  10.1086/147124
    [2] BURLAGA L F,NESS N F,ACUÑA M H,et al. Magnetic fields at the solar wind termination shock[J]. Nature,2008,454(7200):75-77. doi:  10.1038/nature07029
    [3] KRIMIGIS S M,DECKER R B,ROELOF E C,et al. Search for the exit:Voyager 1 at heliosphere’s border with the galaxy[J]. Science,2013,341(6242):144-147.
    [4] SARABHAI V. Some consequences of nonuniformity of solar wind velocity[J]. Journal of Geophysical Research,1963,68(5):1555-1557. doi:  10.1029/JZ068i005p01555
    [5] GOPALSWAMY N. Properties of interplanetary coronal mass ejections[J]. Space Science Reviews,2006,124:145-168.
    [6] WHANG Y C,BURLAGA L F. Evolution and interaction of interplanetary shocks[J]. Journal of Geophysical Research,1985,90(A11):10765-10778. doi:  10.1029/JA090iA11p10765
    [7] BURLAGA L F,NESS N F,BELCHER J W. Radial evolution of corotating merged interaction regions and flows between ~14 AU and ~43 AU[J]. Journal of Geophysical Research Space Physics,1997,102(A3):4661-4672. doi:  10.1029/96JA03629
    [8] WANG C,RICHARDSON J D. Development of a strong shock in the outer heliosphere[J]. Geophysical Research Letters,2002,29(8):1181-1188. doi:  10.1029/2001GL014250
    [9] BURLAGA L F,NESS N F. Observations of the interstellar magnetic field in the outer heliosheath:Voyager 1[J]. The Astrophysical Journal,2016,829(2):1-10.
    [10] GURNETT D A,KURTH W S,STONE E C,et al. Precursors to interstellar shocks of solar origin[J]. The Astrophysical Journal,2015,809(2):1-10.
    [11] LIU Y D,RICHARDSON J D,WANG C,et al. Propagation of the 2012 march coronal mass ejections from the Sun to heliopause[J]. The Astrophysical Journal Letters,2014,788(2):1-6.
    [12] FERMO R L, POGORELOV N V, BURLAGA L F. Transient shocks beyond the heliopause[C]//14th Annual International Astrophysics Conference: Linear and Nonlinear Particle Energization Throughout the Heliosphere and Beyond. Tampa, Florida, USA: IOP, 2015.
    [13] KIM T K,POGORELOV N V,BURLAGA L F. Modelling shocks detected by Voyager 1 in the local interstellar medium[J]. The Astrophysical Journal Letters,2017,843(2):1-7.
    [14] POGORELOV N V,SUESS S T,BOROVIKOV S N,et al. Three-dimensional features of the outer heliosphere due to coupling between the interstellar and interplanetary magnetic fields. IV. Solar cycle model based on Ulysses observations[J]. The Astrophysical Journal,2013,772(2):1-17.
    [15] RICHARDSON J D,WANG C,LIU Y D,et al. Pressure pulses at Voyager 2:drivers of interstellar transients?[J]. The Astrophysical Journal,2017,834(2):1-6.
    [16] PAULS H L,ZANK G P,WILLIAMS L L. Interaction of the solar wind with the local interstellar medium[J]. Journal of Geophysical Research,1995,100(A11):21595-21604. doi:  10.1029/95JA02023
    [17] GUO X C,FLORINSKI V,WANG C. A global MHD simulation of outer heliosphere including anomalous cosmic-rays[J]. The Astrophysical Journal,2019,879(2):1-10.
    [18] FLORINSKI V,GUO X,BALSARA D S,et al. Magnetohydrodynamic modeling of solar system processes on geodesic grids[J]. The Astrophysical Journal Supplement Serisa,2013,205(2):1-12.
    [19] GUO X C. An extended HLLC Riemann solver for the magneto-hydrodynamics including strong internal magnetic field[J]. Journal of Computional Physics,2015,290:352-363. doi:  10.1016/j.jcp.2015.02.048
    [20] AXFORD W 1. The interaction of the solar wind with the interstellar medium[C]//In Solar Wind: Proceedings of the Second International Conference. [S.l.]: NASA, 1972.
    [21] WANG C,RICHARDSON J D. Energy partition between solar wind protons and pickup ions in the distant heliosphere:a three-fluid approach[J]. Journal of Geophysical Research,2001,106(A12):29401-29408. doi:  10.1029/2001JA000190
    [22] GUO X C,FLORINSKI V. Galactic cosmic-ray intensity modulation by corotating interaction region stream interfaces at 1 AU[J]. The Astrophysical Journal,2016,826(1):1-13. doi:  10.3847/0004-637X/826/1/1
    [23] RICHARDSON J D,BELCHER J W,GARCIA-GALINDO P,et al. Voyager 2 plasma observations of the heliopause and interstellar medium[J]. Nature,2019(3):1019-1023.
    [24] GURNETT D A,KURTH W S. Plasma densities near and beyond the heliopause from the Voyager 1 and 2 plasma wave instruments[J]. Nature Astronomy,2019,3(11):1024-1028. doi:  10.1038/s41550-019-0918-5
    [25] BALE S D,REINER M J,BOUGERET J L,et al. The source region of an interplanetary type II radio burst[J]. Geophysical Research Letters,1999,26(11):1573-1576. doi:  10.1029/1999GL900293
    [26] GAZIS P R. Solar wind stream structure at large heliocentric distances:pioneer observations[J]. Journal of Geophysical Research,1987,92(A3):2231-2242. doi:  10.1029/JA092iA03p02231
    [27] SZABO A, SMITH C W, SKOUG R M. The transition of interplanetary shocks through the magnetosheath[C]//AIP Conference Proceedings. Melville, NY: AIP, 2003.
    [28] JIAN L K,RUSSELL C,LUHMANN J G,et al. STEREO observations of interplanetary coronal mass ejections in 2007–2016[J]. The Astrophysical Journal,2018,855(2):1-19.
    [29] STONE E C,CUMMINGS A C,MCDONALD F B,et al. Voyager 1 observes low-energy galactic cosmic rays in a region depleted of heliospheric ions[J]. Science,2013,341(6142):150-153. doi:  10.1126/science.1236408
    [30] CUMMINGS A C,STONE E C,HEIKKILA B C,et al. Galactic cosmic rays in the local interstellar medium:Voyager 1 observations and model results[J]. The Astrophysical Journal,2016,831(1):1-21. doi:  10.3847/0004-637X/831/1/1
    [31] BURLAGA L F,NESS N F. Observations of the interstellar magnetic field in the outer heliosheath:Voyager 1[J]. The Astrophysical Journal,2013,829(2):1-10.
  • [1] 谭宝林, 谭程明, 黄静, 陈林杰.  空间甚低频太阳射电Ⅲ型爆研究进展 . 深空探测学报(中英文), 2021, 8(1): 92-99. doi: 10.15982/j.issn.2096-9287.2021.20190227002
    [2] 王玲华, 宗秋刚, 任杰.  太阳系边际的能量粒子探测 . 深空探测学报(中英文), 2020, 7(6): 567-573. doi: 10.15982/j.issn.2096-9287.2020.20200061
    [3] 辛鹏飞, 吴跃民, 荣吉利, 危清清, 刘宾, 刘鑫.  月面探测器圆形薄膜太阳翼展开动力学建模与分析 . 深空探测学报(中英文), 2020, 7(3): 255-263. doi: 10.15982/j.issn.2095-7777.2020.20191128005
    [4] 何建森, 林荣, 崔博, 王玲华, 宗秋刚.  外日球层的宽能段离子及其与湍动的耦合作用 . 深空探测学报(中英文), 2020, 7(6): 574-583. doi: 10.15982/j.issn.2096-9287.2020.20200064
    [5] 宗秋刚, 任杰, 何建森, 王玲华.  从地球磁层到外日球层及以远区域探测 . 深空探测学报(中英文), 2020, 7(6): 554-559. doi: 10.15982/j.issn.2096-9287.2020.20200063
    [6] 宁晓琳, 桂明臻, 孙晓函, 刘劲, 吴伟仁.  一种基于太阳震荡时间延迟量测的自主天文导航方法 . 深空探测学报(中英文), 2019, 6(1): 88-95. doi: 10.15982/j.issn.2095-7777.2019.01.013
    [7] 杨涛, 邵志杰, 蔡明辉, 贾鑫禹, 韩建伟.  空间高能粒子与器件布线层核反应后次级粒子LET分布研究 . 深空探测学报(中英文), 2019, 6(2): 173-178. doi: 10.15982/j.issn.2095-7777.2019.02.009
    [8] 敖先志, 刘四清, 沈华, 王晶晶, 胡骏翔, 李刚.  2 AU以内的“渐进型”太阳高能粒子事件模拟 . 深空探测学报(中英文), 2019, 6(2): 156-164. doi: 10.15982/j.issn.2095-7777.2019.02.007
    [9] 杨轩, 鄢建国, 叶茂, 金炜桐, 曲春凯, 刘素艳.  对一种月球与火星探测多程微波测量链路定轨定位的数值模拟初步分析 . 深空探测学报(中英文), 2018, 5(2): 154-161. doi: 10.15982/j.issn.2095-7777.2018.02.007
    [10] 黄护林, 李林永, 李来, 刘飞标.  等离子体磁流体发电研究进展 . 深空探测学报(中英文), 2018, 5(4): 331-346. doi: 10.15982/j.issn.2095-7777.2018.04.003
    [11] 高珊, 周文艳, 梁伟光, 刘德成, 唐玉华, 杨维廉.  地月拉格朗日L2点中继星轨道分析与设计 . 深空探测学报(中英文), 2017, 4(2): 122-129. doi: 10.15982/j.issn.2095-7777.2017.02.004
    [12] 何芸, 刘祺, 田伟, 段会宗, 叶贤基, 范淑华, 李语强.  地月第二拉格朗日点卫星激光测距技术研究 . 深空探测学报(中英文), 2017, 4(2): 130-137. doi: 10.15982/j.issn.2095-7777.2017.02.005
    [13] 邬静云, 高有涛.  利用绳系太阳帆减缓小行星自转的技术研究 . 深空探测学报(中英文), 2016, 3(1): 47-50. doi: 10.15982/j.issn.2095-7777.2016.01.007
    [14] 张建琴, 徐建明, 贾巍, 邱宝贵, 肖杰.  深空探测太阳电池阵应用及关键技术分析 . 深空探测学报(中英文), 2016, 3(1): 3-9. doi: 10.15982/j.issn.2095-7777.2016.01.001
    [15] 胡海岩.  太阳帆航天器的关键技术 . 深空探测学报(中英文), 2016, 3(4): 334-344. doi: 10.15982/j.issn.2095-7777.2016.04.005
    [16] 葛丹桐, 崔平远.  地外天体着陆点选择综述与展望 . 深空探测学报(中英文), 2016, 3(3): 197-203. doi: 10.15982/j.issn.2095-7777.2016.03.001
    [17] 姜生元, 朴松杰, 张伟伟, 沈毅, 侯绪研, 全齐全, 邓宗全.  地外天体潜入式探测典型案例分析及展望 . 深空探测学报(中英文), 2016, 3(1): 68-76. doi: 10.15982/j.issn.2095-7777.2016.01.011
    [18] 曾祥远, 龚胜平, 李俊峰, 蒋方华, 宝音贺西.  应用太阳帆悬停探测哑铃形小行星 . 深空探测学报(中英文), 2015, 2(1): 48-52. doi: 10.15982/j.issn.2095-7777.2015.01.007
    [19] 贾巍, 倪家伟, 黄三玻, 宗魏, 肖杰, 王训春, 池卫英.  火星尘埃对太阳电池阵的影响与电帘除尘研究 . 深空探测学报(中英文), 2014, 1(4): 303-307. doi: 10.15982/j.issn.2095-7777.2014.04.010
    [20] 倪彦硕, 宝音贺西, 李俊峰.  考虑太阳摄动的小行星附近轨道动力学 . 深空探测学报(中英文), 2014, 1(1): 67-74.
  • 加载中
图(3)
计量
  • 文章访问数:  85
  • HTML全文浏览量:  19
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-19
  • 修回日期:  2020-09-25
  • 网络出版日期:  2021-01-18
  • 刊出日期:  2021-01-27

外日球层激波事件的一维磁流体力学数值模拟

doi: 10.15982/j.issn.2096-9287.2020.20200059
    基金项目:  中国科学院空间科学战略先导科技专项资助项目(XDB41000000,XDA15052500);国家自然科学基金资助项目(41874171,41674146,41574159,41731070);中国科学院前沿科学重点研究计划资助项目(QYZDJ-SSW-JSC028);民用航天技术预先研究资助项目(D030202,D020301)
    作者简介:

    郭孝城(1978– ),男,研究员,主要研究方向:日球层物理、磁层物理。 通讯地址:北京市海淀区中关村南二条1号中国科学院国家空间科学中心(100190)电话:(010)62586362 E-mail:xcguo@swl.ac.cn

  • ● It is the first time to use the solar wind data from OMNI,the spacecraft STEREO A and B orbiting the Sun conjunctly to study the evolution and propagation of the solar wind in the outer heliosphere. ● The numerical simulations have been carried out to verify the possible correspondence between the large-scale pressure pulse structure in the heliosphere and the shock event in interstellar space.
  • 中图分类号: V467.3

摘要: “旅行者1号”在2012年8月穿越日球层顶后连续探测到了星际磁场的局地突然变化,通常被认为对应于激波事件,且为太阳风事件向外传播后形成的大尺度结构与日球层顶相互作用的产物。但由于观测的限制,日球层内的太阳风大尺度结构与星际空间的激波事件的对应关系尚未明确。通过一个包含太阳风等离子体和星际中性原子的双流体磁流体力学数值模型检验了这种对应关系。在距日1 AU的内边界输入2010—2017年分别位于不同经度上的OMNI,STEREO A和STEREO B的太阳风观测数据,研究了这期间太阳风与星际中性原子电荷交换后在外日球层区的传播和演化,并与“旅行者1号”和“旅行者2号”的观测数据做对比。结果表明,“旅行者1号”探测的星际空间激波事件与“旅行者2号”探测到的日球层鞘区内的压力脉冲事件具有明显的关联。

注释:
1)  ● It is the first time to use the solar wind data from OMNI,the spacecraft STEREO A and B orbiting the Sun conjunctly to study the evolution and propagation of the solar wind in the outer heliosphere. ● The numerical simulations have been carried out to verify the possible correspondence between the large-scale pressure pulse structure in the heliosphere and the shock event in interstellar space.

English Abstract

郭孝城, 周昱成, 王赤, 李晖. 外日球层激波事件的一维磁流体力学数值模拟[J]. 深空探测学报(中英文), 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
引用本文: 郭孝城, 周昱成, 王赤, 李晖. 外日球层激波事件的一维磁流体力学数值模拟[J]. 深空探测学报(中英文), 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
GUO Xiaocheng, ZHOU Yucheng, WANG Chi, LI Hui. One-Dimentional Numerical MHD Simulation of the Shock Events in the Outer Heliosphere[J]. Journal of Deep Space Exploration, 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
Citation: GUO Xiaocheng, ZHOU Yucheng, WANG Chi, LI Hui. One-Dimentional Numerical MHD Simulation of the Shock Events in the Outer Heliosphere[J]. Journal of Deep Space Exploration, 2020, 7(6): 560-566. doi: 10.15982/j.issn.2096-9287.2020.20200059
    • 源于太阳表面的超声速太阳风等离子体向外径向膨胀时会与包含星际等离子体和中性原子的邻近星际介质相互作用形成一泡状结构即日球层,它们的交界面就是日球层顶[1]。在经历漫长的旅程后,“旅行者1号”( Voyager 1)和“旅行者2号”( Voyager 2)飞船分别在距太阳94 AU和84 AU左右穿越终止激波,最终在120 AU左右穿越日球层顶进入星际空间[2-3]。太阳风事件,包括共转太阳风作用区(Corotating Interaction Regions,CIRs)[4]和行星际日冕抛射事件(Interplanetary Coronal Mass Ejections,ICMEs)[5],在日球层内传输时会引起系列反应,比如飞船发现在几十AU处太阳风事件传输时会逐渐汇聚合并形成更大尺度的等离子体脉冲或激波事件[6-8]

      “旅行者1号”在2012年8月穿越日球层顶后,磁强计在2012.92(十进制年份)和2014.66年(十进制年份)分别探测到两个前向激波,在2013.35年探测到一个疑似后向激波[9]。这些激波被认为与“旅行者1号”的电场波动仪器探测到的射电爆信号相关[10]

      学者尝试对这些星际空间激波的来源做分析。Liu等[11]首先采用了一维磁流体力学数值模拟研究了太阳风事件在外日球层区的传播和演化,提出“旅行者1号”在2013.35年观测到的射电爆和激波信号可能源自Wind飞船在2012年3月探测到的一系列的行星际日冕物质抛射事件。然而,他们的模拟没有考虑星际等离子体的入流,所以结果没有终止激波和日球层鞘区的形成,太阳风始终为超声速;在计算激波传播时类比地球弓激波上下游的情况并对预定终止激波以上区域的激波速度进行人为校正。Fermo等[12]采用更复杂的全球磁流体力学数值模拟研究了太阳风事件在外日球层的传播和演化,在距日1 AU处引入OMNI(https://omniweb.gsfc.nasa.gov)的太阳风观测数据,发现产生的日球层顶和激波位置比实际观测到的数据要远约30 AU。后来,Kim等[13]扩展了该工作,在1 AU的低纬度地区采用OMNI数据,但高纬度区采用经验给出的高速太阳风数值[14]。结果显示多重的太阳风共转作用区对第一个前向激波的形成具有重要的作用,而第二个前向激波可能与太阳风共转作用区和行星际日冕抛射事件的共同作用有关。

      当“旅行者1号”在星际空间探测到这3个激波时,“旅行者2号”此时在内日球层鞘区的太阳风里。Richandson等[15]通过“旅行者2号”的观测数据发现鞘区内的太阳风大尺度压力脉冲与星际空间的激波具有一定的关联性,经过对波速的计算认为5个合并太阳风作用区可能最终驱动了相应激波在星际空间的形成。然而,由于“旅行者1号”和“旅行者2号”并不沿着相同的方向飞行,这些压力脉冲事件与星际激波的因果关联仍有待确认。本文将通过一个一维磁流体力学数值模拟对这个因果关系进行分析和讨论。

    • 首先假设太阳风在赤道面附近为沿着径向方向运动的球对称流,这里忽略了太阳风在传播过程中可能受到的来自侧面方向太阳风的影响。此时理想磁流体力学方程组在球坐标系下可以简化为沿着径向${{r}}$方向的一维问题为

      $$\left\{ \begin{array}{*{20}{l}} {\begin{array}{*{20}{l}} {\dfrac{\partial \rho }{\partial t}+\dfrac{1}{{r}^{2}}\dfrac{\partial }{\partial r}\left({r}^{2}\rho {u}_{r}\right)={Q}_{N}}\\ {\dfrac{\partial \left(\rho {u}_{r}\right)}{\partial t}+\dfrac{1}{{r}^{2}}\dfrac{\partial }{\partial r}\left[{r}^{2}\left(\rho {u}_{r}^{2}+{p}_{{\rm{T}}}-{B}_{r}^{2}\right)\right]-\left(2{p}_{{\rm{T}}}-{B}_{\phi }^{2}+\rho {u}_{\phi }^{2}\right)/r-\rho G{M}_{{\rm{s}}}/{r}^{2}={Q}_{{\rm M}r}}\\ { \dfrac{\partial \left(\rho {u}_{\phi }\right)}{\partial t}+\dfrac{1}{{r}^{2}}\dfrac{\partial }{\partial r}\left[{r}^{2}\left(\rho {u}_{r}{u}_{\phi }-{B}_{r}{B}_{\phi }\right)\right]-\left(\rho {u}_{\phi }{u}_{r}-{B}_{\phi }{B}_{r}\right)/r={Q}_{{\rm{M}}\phi }}\\ { \dfrac{\partial {B}_{\phi }}{\partial t}+\dfrac{1}{{r}^{2}}\dfrac{\partial }{\partial r}\left[{r}^{2}\left({u}_{r}{B}_{\phi }-{B}_{r}{u}_{\phi }\right)\right]+\left({u}_{\phi }{B}_{r}-{B}_{\phi }{u}_{r}\right)/r=0} \end{array}}\\ {\dfrac{\partial E}{\partial t}+\dfrac{1}{{r}^{2}}\dfrac{\partial }{\partial r}\left\{{r}^{2}\left[\left(E+{p}_{{\rm{T}}}\right){u}_{r}-{B}_{r}\left({u}\cdot {B}\right)\right]\right\}+\rho G{M}_{{\rm{s}}}{u}_{r}/{r}^{2}={Q}_{{\rm{E}}}} \end{array} \right. $$ (1)

      其中:ρuBE分别是太阳风等离子体的密度、速度、磁场强度和能量密度。

      总压$ {p}_{\rm T}^{\;} $和能量密度${{E}}$的具体形式为

      $$ E=P/({{\gamma}} {-1})+{{\rho}} {{u}}^{2}/2+{{ B}}^{2}/2 $$ (2)
      $$ p_{\rm T} = p + { B^2}/2 $$ (3)

      其中:$ p $是等离子热压。

      考虑了太阳风等离子体离子与星际中性原子的电荷交换,方程组右边的$ {Q}_{N} $${Q}_{{\rm{M}}r}$${Q}_{{\rm{M}}\phi }$${Q}_{{\rm{E}}}$分别表示因电荷交换引起的质量、径向动量、方位角动量和能量的变化[16-17]。太阳的重力效应也被考虑,其中$ G $是万有引力常数,${M}_{{\rm{s}}}$是太阳质量。采用MUSCL数值格式和有限体积方法数值求解上述方程组[18],并运用拓展型HLLC黎曼算子完成对方程通量的计算[19],在时间推进上采用龙格–库塔格式,最终数值格式在空间和时间上均具有二级精度。

      数值模拟区域范围为r方向从1~200 AU,共1万个网格点。采用非均匀网格,设N为总格点数,则格面值${r}_{i}(i=3, \cdots ,N-1)$由代数式为

      $$ {r}_{i}={{{r}}}_{\rm{min}}+ \left({{{r}}}_{\rm{max}}-{{{r}}}_{\rm{min}}\right)\left[{\rm{e}}^{{\rm{\delta }} \left(i-3\right)/(N-4)}-1\right]/\left({\rm{e}}^{\rm{\delta }}-1\right)$$ (4)

      其余四个格面值$ {r}_{1} $$ {r}_{2} $$ {r}_{N} $$ {r}_{N+1} $由插值给出;而格心值为两相邻格面值的中点;$ {\rm{\delta }} $参数表征最大网格与最小网格的差异程度,$ {\rm{\delta }}\to 0 $时网格趋向均匀网格,这里取$ {\rm{\delta }} $ = 2.0。这种网格规模能够使得模拟获得较高的分辨率,比如在内边界网格尺寸约0.006 2 AU,在终止激波(~90 AU)附近约0.024 AU,在日球层顶(~120 AU)约为0.03 AU。初始时刻,内边界处太阳风取值为:数密度5 ${\rm /cm}^{3}$,速度${{{u}}}_{{{r}}}=400\;{\rm{k}}{\rm{m}}/{\rm{s}}$,温度$ 1.5\times {10}^{5}\;{\rm{K}} $,磁场强度${B}_{r}=2.8 \;\rm nT$。星际介质方面,这里只考虑中性原子,其分布形式为[20]

      $$ {n}_{H}={n}_{H0}{\rm e}^{ \left(-\dfrac{{r}_{0}}{r}\right)} $$ (5)

      其中:${n}_{H0}=0.15\;{\rm /cm}^{3}$为临近星际空间的中性原子密度初始值;${r}_{0}=4\;\rm AU$为星际中性原子进入日球层的深度,用距离日心的距离表征。

      简化计,模拟没有考虑星际等离子体,因此太阳风与星际介质相互作用的形式主要以电荷交换形式出现。

      从初态开始,太阳风经过5年的演化达到准稳态,如图1的点线所示。注意到,由于没有星际等离子体的作用,终止激波和日球层顶均没有出现。由于太阳风离子与中性原子的电荷交换作用导致太阳风被加热,在外日球层区域其温度沿着径向并没有出现绝热冷却的现象,而是在5 AU以上逐渐升高[21]。进一步以此为初态,在1 AU的内边界处代之以2010—2017年间的OMNI、STEREO A和B(以下简写STA和STB)等飞船的观测数据。因为STA在2014年9月—2015年11月间部分数据缺失,而STB自2014年10月起已经失去联系,所以缺失的数据这里以上述恒定的初态条件代替。后面讨论的外日球层激波事件在2015年之前,假设激波速度500 km/s,那么从1~120 AU以上区域激波传播需要1.14 a左右,因此可以认为这种代替的负面影响是小的。图1中的实线表示以OMNI数据为输入,模拟数据在2013.42年时刻的物理量径向分布,可见在太阳风事件沿着径向方向向外传播,出现的激波信号将会在120 AU以上的星际空间被探测到。

      图  1  沿着径向方向的太阳风参数

      Figure 1.  The profiles of solar wind plasma along radial direction

    • 需要说明的是,在模拟中太阳风被简化为接近赤道面的球对称流。“旅行者1号”在2010—2017年间在HGI(Heliopsheric Inertial Coordinate)坐标系下的纬度变化为34.4°N~ 34.8°N,“旅行者2号”从32.3°N~28.8°N;经度变化上,“旅行者1号”是173.6°S~175.1°S,“旅行者2号”是216.8°S~218.3°S;可见两艘飞船均离赤道面有一定的距离。但由于太阳风结构在外日球层区域传播存在大区域传播的特性,在经度和纬度上覆盖面广[22]。这里假设模拟的结果能够一定程度上反映两艘飞船观测得到的等离子体大尺度结构的特征。

    • “旅行者2号”最终在2007年距日84 AU处穿越终止激波[2],之后一直在内日球层鞘区直至2018年底被确认穿越日球层顶进入星际空间[23]。给出2011—2017年间的模拟与观测的动压对比,在此期间“旅行者2号”从距日94.17 AU处运动至113.14 AU,如图2所示。绿色横线段表示发生在星际空间的等离子体振荡事件,由“旅行者1号”携带的等离子体电场波动仪器测量得到[24]。绿色的竖线对应于“旅行者1号”磁强计测量得到的磁场跳变信号,对应于星际空间传播的激波[9]

      图  2  2011—2017年采用3种不同内边界数据源的太阳风动压对比

      Figure 2.  The comparision of solar wind dynamical pressure between the observation(red)and the simulations(blue)during the years 2011—2017

      一般认为图2中的前3个激波信号与等离子体振荡时间关联得较好,后者被认为来自星际激波前方的低能电子束流,类似于II型的太阳射电爆发信号[25]。第4个等离子体振荡信号由于没有对应的磁场跳变,一般认为它可能来源于星际其他地方的激波且未被“旅行者1号”的磁强计探测得到。6个字母A~F和对应的垂直虚线表示“旅行者2号”观测到的太阳风合并作用区动压的局地最高点,这些压力脉冲被认为与“旅行者1号”探测到的等离子体爆发事件信号相关,其中C、D和E被认为对应于第1、3和4个等离子体振荡[15]

      由于没有终止激波,模拟得到的动压相对观测值较大;出于方便考虑,将观测的动压结果乘以1.5后与模拟结果做对比。图2中可见,很难将模拟结果与“旅行者2号”的观测做精确的逐一比较,一方面是模拟的有效区域在近赤道面,与“旅行者2号”的纬度存在较大差异,另外球对称流不存在与径向方向相切的流动,实际的太阳风传播要复杂得多,此外终止激波的缺失也是一方面,以及其他的物理因素如湍流、拾起离子的影响等。

      大体上,3个不同的内边界输入源会导致不同的模拟结果,但趋势大致相近,存在近似1年左右的长周期变化[26]。观测数据发现,激波穿越地球弓激波后速度为原先速度的0.7~1倍左右[27]。因此,设终止激波在84 AU左右,激波速度在~500 km/s,假设其在图2中距离(94~113 AU)上穿越终止激波后以350 km/s运动10~30 AU左右,预计图2中的模拟激波信号比实际传播快0.04~0.12年左右,这种差异在图中的刻度范围内并不显著。从图2中看,OMNI源的模拟结果对A、B、C和D动压脉冲具有一定的表现,但对后期两个压力脉冲E和F并不明显;而STA源的模拟结果对B、D和F具有一定的表现,对A、C和E表现不明显;对STB源的模拟结果来说,前期在动压的变化趋势上比较一致,但对后期的脉冲比如F表现不明显,这种情况很可能是因为STB数据在自2014年10月后的缺失引起。

      模拟结果表明,对一维模拟来说,1 AU处不同经度上的内边界输入在外日球层区虽然在大尺度结构上具有一定的一致性,但细节方面仍具有差异,原因在于OMNI,STA和STB在1 AU轨道上的存在经度差异,它们对太阳风事件的测量并不总是一致,比如存在有些太阳风事件只被某一经度上的飞船探测到的情况[28]。因此在与“旅行者2号”的观测数据对比时,需要综合考虑3个不同内边界输入的情况。以往人们在模拟外日球层太阳风传播时侧重于使用OMNI数据[13],该模拟算例表明在更高维度上如何处理好内边界问题是一个很具有挑战性的难题。

    • “旅行者1号”在2013年8月于距日121.6 AU处穿越了日球层顶[29],自那以后人们开始能够持续测量星际空间的等离子体波动[10],银河宇宙线[30]和磁场信号[9]。然而,由于“旅行者1号”的等离子体仪器自1980年开始就已经失效,无法将模拟得到的等离子体参数与观测数据直接比较。

      图3分别给出了采用OMNI,STA和STB源作为内边界输入模拟得到“旅行者1号”自2012—2016年间测量得到的等离子体动压变化。黑色竖虚线表示“旅行者1号”在2012.65时刻观测到的穿越日球层顶的位置,另外3个竖线表征时刻2012.904、2013.356和2014.644,分别对应于“旅行者1号”通过磁场强度跳变观测得到的前向激波FS1,后向激波RS和前向激波FS2[9]。可以看出,模拟得到3个激波附近均具有持续时间较长的压力脉冲结构甚至激波组合,时间范围大致包括2012.4—2013.0、2013.2—2013.8和2014.6—2015.2。

      图  3  2012—2016年用3种不同内边界数据源,模拟得到的“旅行者1号”可能探测的动压数据

      Figure 3.  The simulated dynamical pressures that could be encountered by Voyager 1 during the years 2012—2016

      由于模拟并未包含终止激波和日球层顶,实际上这些压力脉冲在与终止激波和日球层顶作用后能够进一步形成激波,进而在星际空间被探测到。以OMNI数据的第2个动压脉冲结构中的2013.36处的激波为例,注意到在模拟后的第1 086 d到达85 AU,而在第1 212 d到达120 AU左右(本文未列出),可以估算激波速度为482 km/s。如果激波在85 AU左右穿越终止激波后速度下降为70%,则在120 AU处激波会迟54 d左右到达。而在星际空间里,磁场观测显示快磁声波速大概为40 km/s[31],因此激波穿越日球层顶后速度会降至原来的1/10左右,这样越往后模拟到的动压脉冲实际对应的激波信号延迟越大。Richardson等[15]提出经过修正后,“旅行者2号”观测到的动压脉冲C、D、E分别对应第1、3和4个等离子体振荡信号,而第2个振荡信号来源仍然未知。比较图2和图3可见,图2中的C,D和E基本对应图3中的三个动压脉冲组合,即模拟得到的3个动压脉冲也同样可能对应于实际观测到的第1、3和4个等离子体振荡信号。此外,从模拟结果看,仍然无法对第2个振荡信号的来源做出判断。

      以往的三维数值模拟结果认为第1前向激波FS1来自于太阳风的共转作用区在外日球层区的合并发展,第2个前向激波是由共转作用区与行星际日冕物质抛射共同作用的结果[13]。可能的后向激波RS最初被认为与2012年3月爆发的大规模日冕物质抛射事件有关[11]

      由于模型的限制,这里没有对模拟得到的3个动压脉冲结构的驱动源做分析,后续将采用更为全面的日球层全球数值模拟开展相关研究。

    • 本文通过简化的一维太阳风−星际中性原子相互作用的磁流体力学数值模型,对太阳风在外日球层区的演化进行了数值模拟,其中内边界设在距日1 AU处,外边界在200 AU的星际空间。采用1 AU轨道上不同经度的OMNI,STA和STB这3个不同观测数据作为输入条件,时间范围从2010—2017年。模拟结果表明,采用不同内边界输入后,虽然整体上看太阳风在外日球层区的大尺度结构比较类似,但细节表现会有不同,因此在做数值模拟时需要做综合考虑。模拟了“旅行者1号”在星际空间可能探测到的大尺度动压脉冲结构,这些结构被认为能够与日球层顶相互作用并进一步产生激波,从而触发“旅行者1号”观测到的等离子电场振荡信号,因此可以认为这类结构也与观测到的星际激波信号相关。

参考文献 (31)

目录

    /

    返回文章
    返回