CentPack Logo Center for Scientific Computation and Mathematical Modeling


Two Dimensional MHD Equations


In this example we use CentPack's two-dimensional second order semi-discrete formulation, centpack_2d_SD2, to evolve the solution of Ideal MHD equations


stands for the total pressure, and the the pressure, p, is related to the conserved quantities through the equation of state

with . The solution is evolved over the computational domain , from to . The initial conditions are those of the Orszag-Tang Vortex System

with periodic boundary conditions, , using a 288 x 288 mesh and cfl number 0.75.

The images below display, from top to bottom and left to right, the initial contours of density and pressure (at t = 0.05), and the velocity and magnetic fields. Click on the individual images to see an animation from to .

Click on the images above to see an animation


The files implementing this exaple are included in the current CentPack distributionCP-1.0. Read here how to run and compile this example.

One dimensional examples

Burgers equation
Euler equations of gas dynamics
MHD equations

Two dimensional examples

Scalar example
Euler equations of gas dynamics
MHD equations


User's guide

Jorge Balbas and Eitan Tadmor