杨维纮2009杭州会议报告.ppt.ppt

空间无碰撞激波的数值模拟 陆启明、杨维纮 2009.7.26 2009杭州会议 Plasma Physics LaboratoryUSTC 主要内容 I. 具有间断面（激波）的问题 II. 空间无碰撞激波的简介 III. 等离子体激波的MHD模拟 IV. 空间无碰撞激波的一维模拟 V. 空间无碰撞激波的二维模拟 Plasma Physics LaboratoryUSTC I. 具有间断面（激波）的问题 若方程 的解为 Plasma Physics LaboratoryUSTC I. 具有间断面（激波）的问题 普通差分格式将出现非物理振荡（色散解） Plasma Physics LaboratoryUSTC I. 具有间断面（激波）的问题 加入人为粘滞后会将激波磨光滑（耗散解） Plasma Physics LaboratoryUSTC I. 具有间断面（激波）的问题 采用WENO(Weighted Essentially Non-Oscillatory)格式的解为 (Guang-Shan Jiang,eta, J. Comp.Phys. 150, 561, (1999)) Plasma Physics LaboratoryUSTC II. 空间无碰撞激波简介 Photograph of Bow Shock in Airflow Supersonic Airflow over sphere Plasma Physics LaboratoryUSTC II. 空间无碰撞激波简介 Numerical Simulation of Bow Shock with Fluid Dynamics Plasma Physics LaboratoryUSTC II. 空间无碰撞激波简介 Earth’s Bow Shock Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 Ideal MHD Equations Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 Conservation Form of Ideal MHD Equations Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 Vector of primitive variables Transformation: Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 The flux matrix Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 With eigenvalue Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 Scheme used for spatial discretization–   5th order WENO scheme Scheme for time evolution –  4th order Runge-Kutta scheme Plasma Physics LaboratoryUSTC III. 等离子体激波的MHD模拟 Rankine – Hugoniot jumping relationship –  used to identify the shock waves   Plasma Physics LaboratoryUSTC IV. 空间无碰撞激波的一维模拟 Interaction between quasi-perpendicular collisionless shock and plasmoid with high density 1. The Riemann data for quasi-perpendicular shock: 2. With high density plasmoid moving along with upstream flux 3. The result is similar with that of particle simulation 4. A new fast rarefaction wave is found in the simulation Plasma Physics LaboratoryUSTC Collisionless Shock  Plasmoid Plasma Physics Laborat