日冕物質拋射(CME)是災害性空間天氣的主要驅動源，如何針對特定的CME/太陽風暴事件,了解其日冕/行星際傳播及演化過程，預報是否以及何時到達地球軌道，預測近地空間太陽風狀態(tài)的變化，是空間天氣科學界一直關注的重要課題之一。

　　空間中心空間天氣學國家重點實驗室沈芳研究員和博士生劉佑生、博士后楊易以漸變圓柱殼（GCS）模型為基礎，建立了一種新的三維CME磁繩模型，在三維MHD背景下對新的CME磁繩模型的初始幾何參數和物理設定對于不同位置處行星際空間模擬結果的影響進行了細致的參數研究。第一個模擬工作主要針對初始CME傳播方向沿著日地連線的情況下展開，討論了CME初始參數，包括密度、幾何大小、質量以及磁場對地球和火星位置處數值模擬結果的影響，結果表明，改變初始磁場會影響CME到達地球時的地磁效應，但幾乎不會影響CME在行星際空間傳播的動力學過程；改變初始密度或幾何尺寸會影響CME的傳播過程，包括CME激波到達地球和火星的時間以及CME在行星際空間的傳播速度；若同時改變初始密度和幾何尺寸但保持初始質量不變，CME的傳播過程則沒有明顯改變，上述工作發(fā)表在The Astrophysical Journal Supplement Series上。

　　隨后，在第一個模擬工作的基礎之上，沈芳等人進一步研究了CME磁繩模型的初始參數對相同日心距、不同經緯度位置的模擬結果的影響，本研究共選取了不同相同質量但不同尺寸、密度及磁場的5個Case展開模擬，分別比較了位于地球附近、南北緯10度、20度、30度；東西經10度、20度和30度位置處的激波到達時間、速度、密度和磁場的峰值。模擬結果表明（圖1和圖2），與第一個模擬結果不同的是，在CME初始質量保持不變的情況下，在遠離CME傳播的方向上，其初始大小和密度會對CME在行星際空間的傳播過程有著顯著影響，并且CME模型的初始大小對不同經緯度位置的影響并不是對稱的，如，初始CME大小對CME在行星際空間傳播過程中緯度方向的影響要遠遠比經度方向顯著。另外，由于CME在行星際空間傳播過程中收到背景太陽風結構的影響，其偏轉在經度和緯度方向均有發(fā)生，具體來說在5個Case中的所有的CME都發(fā)生了向南和向東的偏轉，這個偏轉的過程是由于快速CME與背景太陽風中的CIR（共轉相互作用區(qū)）發(fā)生相互作用產生的。在未來的工作中，將結合Parker Solar Probe，Solar Orbiter，ASOS等衛(wèi)星的觀測數據，對更為豐富的初始參數，不同CME速度，不同太陽風背景，以及對不同的行星際位置處的影響將展開進一步的參數研究，希望對行星際空間的空間天氣數值預報等方面做出新的貢獻。相關工作發(fā)表在The Astrophysical Journal上。

　　日冕區(qū)連接太陽表面和行星際空間，三維日冕太陽風模擬是太陽風三維磁流體模擬中的重要組成部分，其中日冕太陽風加速加熱是日冕三維太陽風模擬中亟待解決的難題之一，而磁場散度處理也是磁流體力學模擬中需要解決的關鍵問題。近日，沈芳研究員和碩士生劉暢等人在太陽風日冕數值建模中磁場散度保持和加熱加速處理方面取得新進展。本研究利用近年來發(fā)展的 三維日冕行星際COIN-TVD MHD模型對日冕區(qū)三維太陽風進行模擬研究，首先采用體積加熱法對日冕加速加熱過程進行模擬，通過調節(jié)體積加熱項的參數，生成符合日冕區(qū)太陽風高溫高速特點的穩(wěn)態(tài)結構；進而探討不同的磁場散度處理方法對日冕太陽風結構的影響，分別采用Powell法、擴散法、擴散-Powell組合法處理磁場散度問題，模擬結果表明（如圖3所示），將Powell法和擴散法組合起來處理磁場散度可將相對磁場散度誤差控制在10-9-10-6，擴散-Powell組合法與COIN-TVD模型結合能充分發(fā)揮其控制磁場離散的能力，它不僅可以有效減小磁場散度誤差，而且可以確保MHD方程的守恒性。相關工作發(fā)表在Frontiers in Physics上。

　　Citation：

　　(1)Shen, F.*, Liu, Y. and Yang, Y., 2021. Numerical Research on the Effect of the Initial Parameters of a CME Flux-rope Model on Simulation Results. The Astrophysical Journal Supplement Series, 253,12. https://doi.org/10.3847/1538-4365/abd4d2

　　(2)Shen, F.*, Liu, Y. and Yang, Y., 2021. Numerical Research on the Effect of the Initial Parameters of a CME Flux-rope Model on Simulation Results. II. Different Locations of Observers. Astrophys. J., 915,30. https://iopscience.iop.org/article/10.3847/1538-4357/ac004e

　　(3)Liu C., Shen F.*, Liu Y., Zhang M. and Liu X. (2021) Numerical Study of Divergence Cleaning and Coronal Heating/Acceleration Methods in the 3D COIN-TVD MHD Model. Front. Phys. 9:705744. doi: 10.3389/fphy.2021.705744

圖 1不同初始幾何形狀的5組CME事件(不同顏色線條)在不同緯度位置處（a-g分別表示地球、北緯10、20、30度和南緯10、20、30度）的數值模擬結果對比（從上到下依次為：密度、速度、總磁場和Bz）

圖 2 不同初始幾何形狀的5組CME事件(不同顏色線條)在不同緯度位置處（a-g分別表示地球、西經10、20、30度和東經10、20、30度）的數值模擬結果對比（從上到下依次為：密度、速度、總磁場和Bz）

圖3 三種磁場散度處理方法對日冕太陽風模擬過程中磁場散度誤差消除作用的比較，其中紅線代表擴散法，綠線表示Powell方法，藍線代表擴散-Powell組合法

　　（供稿：天氣室）