计算物理

• • 上一篇    

一种适用于磁流体力学切向间断的HLLC黎曼算子

习心悦, 郭孝城, 王赤   

  1. 1. 中国科学院国家空间科学中心,北京100190);
    2. 中国科学院大学地球和行星科学学院, 北京100049
  • 收稿日期:2021-07-20 修回日期:2021-10-09 接受日期:2021-10-09
  • 作者简介:习心悦,女,硕士研究生,研究方向为空间等离子体数值模拟 E-mail:xyxi@swl.ac.cn
  • 基金资助:
    国家自然科学基金项目(41874171,41674146,41574159,41731070),中国科学院空间科学战略先导科技专项资助项目(XDB41000000,XDA15052500),中国科学院前沿科学重点研究计划资助项目(QYZDJ-SSW-JSC028), 中国科学院国家空间科学中心大学生“科创计划”的支持

An HLLC Riemann solver designed for MHD tangential discontinuities

Xi Xinyue, Guo Xiaocheng, Wang Chi   

  1. 1. National Space Science Center, Chinese Academy of Sciences, 100190, Beijing, CHINA
    2. School of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Chinese Academy of Sciences, 100049, Beijing, CHINA
  • Received:2021-07-20 Revised:2021-10-09 Accepted:2021-10-09

摘要: 磁场的存在使得磁流体力学特征波模不同于流体力学,因此直接由流体力学HLLC黎曼算子导出的HLLC双中间态在交界的间断面会出现不守恒的问题,通常降级采用HLL单磁场中间态代替HLLC双磁场中间态以实现守恒和计算稳定,代价是切向间断的模拟精度不足。本文对此进行改进,在模拟切向间断时仍然保留原有的HLLC双磁场中间态,同时各守恒量仍然能够满足Toro相容条件;改进型HLLC算子在间断两侧的磁场分量存在差异,因此能够更精确还原切向间断面。基于不同的数值测试,包括一维激波管和切向间断的时变模拟,以及地球磁层三维数值模拟,我们将改进型HLLC算子与已发展的HLLC和HLLD算子的模拟结果进行对比;结果表明相比于已发展的HLLC算子,改进型HLLC算子对切向间断具有更好的捕捉精度,在模拟中能够达到或接近耗时更多的HLLD算子的模拟精度。

关键词: 磁流体力学, 黎曼算子, HLLC, 切向间断

Abstract: The HLLC Riemann solver for magnetohydrodynamics (MHD) numerical simulation was previously developed based on the HLLC solver for hydrodynamics (HD). The inclusion of magnetic field leads to extra characteristic waves, leading to an inconsistency problem of the jump condition across the contact discontinuity. Usually, for the magnetic field, the single HLL intermediate state is used to replace the two HLLC intermediate states to achieve conservation and computational stability, at the cost of insufficient simulation accuracy for tangential discontinuity. In this paper, the previously developed HLLC solver is modified to deal with MHD tangential discontinuities accurately and satisfy the Toro condition. Based on different numerical tests, such as the time-dependent simulation of one-dimensional shocktube, tangential discontinuities, and the global MHD simulation of Earth's magnetosphere, we compare the numerical results of the modified HLLC solver with those of the standard HLLC and HLLD solvers. The results indicate that the modified HLLC solver has better capture accuracy for tangential discontinuities than the standard HLLC solver, and reaches the accuracy of HLLD solver in some situations.

Key words: Magnetohydrodynamics, Riemann solver, HLLC, tangential discontinuities

中图分类号: