求解一维磁流体方程的无震荡中心差分格式.docVIP

求解一维磁流体方程的无震荡中心差分格式 颜克清 魏伟平 (长安大学 理学院 陕西 西安 710064) 摘 要：本文使用求解双曲型守恒律的高分辨率，无震荡中心格式来求解一维理想磁流体力学方程的解。应用二阶和三阶中心格式求解一维激波管问题，两个数值算例表明其结果比迎风格式要好，且中心格式只需要最大速度值的估计，避免了相应雅克比矩阵及其特征信息的估计，同时消除了维数分裂。 关键词：理想磁流体方程；高分辨率中心格式；无震荡重构 【中图分类号O242；O361 文献标识码A Non-oscillatory central scheme for one-dimensional ideal MHD equations YAN Ke-qing, WEI Wei-ping (School of Science, Changan University, Xian 710064, China) Abstract: In this paper,We utilize a family of high-resolution,non-oscillatory central schemes for the approximate solution of the one-dimensional MHD equations. Solution of one-dimensional shock-tube problems is carried out using second- and third- order central schemes. A qualitative comparison reveals an excellent agreement with previous results based on upwind schemes.Central schemes,however, require little knowledge about the eigenstructure of the problem-in fact,we even aviod an explicit evalution of the corresponding Jacobians, while at the same time they eliminate the need for dimensional splitting. Keywords: Ideal magnetohydrodynamics equations; High-resolution central schemes;Non-oscillatory reconstructions 1 引言 在本文中对理想的磁流体力学方程(1)-(4)的近似解我们提出二阶和三阶无震荡中心差分格式 ， （1） ， （2） ， （3） ， （4） 其中，和是标量，分别代表质量密度和总内能，是速度场，它的-二范数表示为，和分别代表磁场及其-二范数。最后，结合压力，内能表示为，在这儿是固定的比热比。通过螺线管型的约束增强系统，也就是说，如果初始时满足条件，那么通过（3）式它将保持不变。 由于磁流体方程的内在复杂性，对于计算磁流体方程（1）-（4）的近似解，中心格式类方法将有效替代迎风格式类。本文我们使用的中心差分格式是基于交错网格上单元平均的演化。中心差分格式的交错型版本[1]-[3]将在下面第2小节中引出。中心差分格式去除了对雅克比矩阵的特征结构详细知识的需要。时间演化上使用简单的求积公式代替了迎风格式的近似黎曼求解器模块。这个方法不仅让我们节省了高成本的雅克比矩阵的特征分解，事实上让我们完全地避免了一维空间中的雅克比矩阵的高代价评估。此外，中心差分格式避免了在多维MHD系统中特别相关的维数分裂。本文用中心差分格式解决一维磁流体方程更简便，更精确。这一点，在后文的数值模拟中我们将看到证明。 在本文中我集中关注2个原型磁流体问题。我们使用二阶和三阶无震荡中心差分格式去计算一维Brio-Wu激波管模型[1]的近似解。两个不同的MHD激波管问题在3.1和3.2小节中研究。 我们的结果发现与解决同类问题的迎风格式[4][5]的模拟结果极好的一致。结果表明用中心格式去探测和解决这类模型的间断解，保持了有效性和简单性。 2 一维中心格式 我们实施交错网格上的预估-校正中心格式近似(1)-(4)的解。此格式分两步： 基于单元平均的点值，进行无震荡分片多项式重构； 基于交错单元平均，实现这