2D AD Magnetohydrodynamic Blast Wave - SFG-CUHK/scorpio-test-cases GitHub Wiki

Two-fluid MHD was introduced by Tilley & Balsara (2008) with isothermal two-fluid MHD, and the adiabatic version in Tilley, Balsara & Meyer (2012). It studied the behavior of the spherical blast wave triggered by the high uniform density region within a radius r and propagate radially outward. Below are the initial conditions of the simulation setting.

Parameter Setting Parameter Setting
x, y [0.1 x 0.1] Nx, Ny 400
Gammai 5/3 Gamman 1.4
rhoiin 0.0126 rhoiout 0.00126
rhonin 10 rhonout 0.1
pnin 4347.8 pnout 4.3478
piin 0.43478 piout 4.3478 x 10-4
mui 29.0 mun 2.3
r 0.1 alphaAD 3.5468 x 10 -2
B 0,sqrt8 vi, vn 0
CFLi 0.2 CFLn 0.2
Solveri HLLMHD Solvern HLLCHD
tend 0.006

Below is the output of the code. Figure 1 shows the density of ions (left) and neutrals (right). Figure 2 shows the cross-section along the x-direction of density (left) and velocity (right) of the two types of fluids. Similar result compares with Figure 5 & Figure 6 (c, f) of Tilley, Balsara & Meyer (2012).

Figure 1. Result of the 2D AD MHD blast wave, simulated with Scorpio.

Figure 2. Result of the 2D AD MHD blast wave, simulated with Scorpio.

⚠️ **GitHub.com Fallback** ⚠️