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
where
![](../equation_images/figures/MHD1d_total_press.png)
stands for the total pressure, and the the pressure, p, is related to the conserved quantities through the equation of state
![](../equation_images/figures/MHD1d_state.png)
with . The solution is evolved over the computational domain , from to . The initial conditions are those of the Orszag-Tang Vortex System
![](../equation_images/figures/MHD2d_ic.png)
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 .
![](../animations/mhd2d/MHD_rho_01.gif) ![](../animations/mhd2d/MHD_p_01.gif)
![](../animations/mhd2d/MHD_u_00.gif) ![](../animations/mhd2d/MHD_b_00.gif)
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.
|