Lab Home | Phone | Search
Center for Nonlinear Studies  Center for Nonlinear Studies
 Home 
 People 
 Current 
 Affiliates 
 Visitors 
 Students 
 Research 
 ICAM-LANL 
 Publications 
 Conferences 
 Workshops 
 Sponsorship 
 Talks 
 Colloquia 
 Colloquia Archive 
 Seminars 
 Postdoc Seminars Archive 
 Quantum Lunch 
 CMS Colloquia 
 Q-Mat Seminars 
 Q-Mat Seminars Archive 
 Archive 
 Kac Lectures 
 Dist. Quant. Lecture 
 Ulam Scholar 
 Colloquia 
 
 Jobs 
 Students 
 Summer Research 
 Visitors 
 Description 
 Past Visitors 
 Services 
 General 
 
 History of CNLS 
 
 Maps, Directions 
 CNLS Office 
 T-Division 
 LANL 
 
Monday, August 10, 2009
10:00 AM - 11:00 AM
CNLS Conference Room (TA-3, Bldg 1690)

Seminar

Godunov-type Algorithms for MHD

Thomas Gardiner
Sandia National Laboratories

Godunov type integration algorithms have ideal dissipation characteristics for simulating moderate to high Mach number hydrodynamic systems. Unfortunately, the extension of these methods to handle ideal magnetohydrodynamics (MHD) has been traditionally hampered by stability issues related to the magnetic charge constraint, div B = 0. Adopting the method of Constrained Transport (CT) ameliorates these stability issues, but raises additional stability concerns resulting from the CT electric field algorithm which relates the Godunov fluxes to the CT electric fields. In this talk I will review the approximations (and exact relationships) which underlie the development of finite volume, Godunov-type algorithms for MHD using CT and in particular present the essential elements which went into developing the MHD algorithms in the Athena simulation code. These include the development of a new set of CT electric field algorithms, the extension of the PPM interface state algorithm for multidimensional MHD, and the incorporation of these algorithms into the unsplit corner transport upwind (CTU) integration algorithm. Along the way, I will present a variety of test problems to help clarify the discussion and highlight the characteristics of the integration algorithms. As time permits, I will close the talk with a discussion of current algorithm development efforts including adaptive mesh refinement and non-ideal extensions.

Host: Mikhail Shashkov