无振荡有限紧致格式在理想磁流体方程中的应用研究.pdfVIP

LHD 2012 年度夏季学术研讨会 无振荡有限紧致格式在 理想磁流体方程中的应用研究 刘利，申义庆 中国科学院力学研究所，高温气体动力学国家重点实验室（筹），北京海淀区 100190 摘要 尽管理想磁流体力学方程不是严格的双曲型方程，但由于有与流体力学方程类似的波特征结 构，因而大量的流体力学数值方法已被推广应用于磁流体方程的数值求解。有限紧致格式通过无参 数判断解在格式模板上的光滑性，将流场划分为激波区域与光滑区域，在激波区域使用激波捕捉格 式，而在光滑区域使用紧致格式，激波区域的数值通量自动作为紧致格式的边界通量，从而既能满 足无振荡捕捉激波，而在光滑区又能够获得较高的数值精度。本文发展磁流体流场的无振荡有限紧 致格式求解方法，并通过数值算例验证方法的精度及有效性。 关键词 理想磁流体，WENO 格式，有限紧致格式 理想 MHD 流体控制方程以向量形式可表 引 言 示为[5] ： 流动中的磁流体力学(MHD) 问题在近年 U F 0 (1) 来受到越来越多的关注，特别是 MHD 在热防 t 护、增加进气道流量捕获、提高燃烧室燃烧效 其中： 率等的研究方面。尽管理想的 MHD 控制方程  V  不是严格的双曲型方程，但由于有与流体力学       方程类似的波特征结构，因而大量的流体力学 V VV+p t I BB  U=  ，F=   VB BV  数值方法已被推广应用于磁流体方程的数值 B 求解[1-8]。     e e p V B V B     t      为了能准确捕捉流场中的激波以及复杂 在笛卡尔坐标系下表示为： 区域的流动小尺度结构，申义庆等[9-11]提出 U E F G    0 （2 ） 了有限紧致格式的概念并发展了系列的高分     t x y z 辨率有限紧致格式。其数值通量的构造方法主 其中 要是利用有效激波判别，将传统紧致格式在计