MHD现代算法介绍及存在的问题.pptVIP

Thanks More … Diffusive Control(Dedner)Modify induction equation to diffuse away div(B) Hyperbolic divergence cleaning(Dedner et al.)may not work well for periodical boundary conditions Vector potential(Londrillo Del Zanna)use vector potential A instead of Bnumerically identical to CT methodreduce accuracy Better to combine diffusive control, etc. with 8-wave formulation * MHD现代算法介绍及存在的问题 马骁鹏 2009.9 Overview MHD SchemesClassical schemes: Lax-Wendroff, MacCormack, Lax-Friedrichs etc.Riemann solvers: Godunov, Roe, MUSCL, PPM, HLL( also HLLC, HLLE, etc.)High order schemes: MUSCL, PPM, WENO, etc. Problems Difficulties in MHD Simulationsdiv(B) control schemes Governing Equations Conservative form Non-conservative form has many advantages which will be discussed later. HD ? MHD Div(B) constraintprojection method, modify MHD equations, staggered grid (CT). Magnetic energy may dominate total energy especially in low-βplasma.Negative pressure density Alfvén speed may dominate maximum speed.CFL condition Other… Classical Schemes Lax-Friedrichs scheme MacCormack scheme Warming-Beam, Upwind, etc. Riemann Problem 1D shock tube Accurate solutions can be obtained but costly (iteration required), especially in MHD. Solutions consist of several waves corresponding to each eigen-value (Advection:1, Euler:3, MHD:7). x t 0 Godunov Scheme Godunov, 1959 Assume constant value in each cell Solve local Riemann problem at …, j-1/2, j+1/2, …accurately. Update cell values with the integral of accurate solution (over time). Only first order accuracy Time-consuming and hard to implement j j+1 j+2 j-1 Approximate Riemann Solvers Linear Riemann SolversSuffer from negative pressure and densityEigen-decomposition requiredhigher resolution HLL-type Riemann Solverspreserve positivityeasy to implementlower resolution Linear Riemann Solvers Reference state: Jacobian Matrix: Eigen-values and (left and right) Eigen vectors: Flux function: Flux Vector splitting requir