Ever since the magnetohydrodynamic （MHD） method for extrapolation of the solar coronal magnetic field was first developed to study the dynamic evolution of twisted magnetic flux tubes, it has proven to be efficient in the reconstruction of the solar coronal magnetic field. A recent example is the so-called data-driven simu- lation method （DDSM）, which has been demonstrated to be valid by an application to model analytic solutions such as a force-free equilibrium given by Low and Lou. We use DDSM for the observed magnetograms to reconstruct the magnetic field above an active region. To avoid an unnecessary sensitivity to boundary conditions, we use a classical total variation diminishing Lax-Friedrichs formulation to iteratively compute the full MHD equations. In order to incorporate a magnetogram consistently and sta- bly, the bottom boundary conditions are derived from the characteristic method. In our simulation, we change the tangential fields continually from an initial potential field to the vector magnetogram. In the relaxation, the initial potential field is changed to a nonlinear magnetic field until the MHD equilibrium state is reached. Such a stable equilibrium is expected to be able to represent the solar atmosphere at a specified time. By inputting the magnetograms before and after the X3.4 flare that occurred on 2006 December 13, we find a topological change after comparing the magnetic field before and after the flare. Some discussions are given regarding the change of magnetic con- figuration and current distribution. Furthermore, we compare the reconstructed field line configuration with the coronal loop observations by XRT onboard Hinode. The comparison shows a relatively good correlation.
ZEUS is a magnetohydrodynamics simulation code widely used in astrophysical research.However,it was recently found that the code may produce artificial shocks in the rarefaction region in some numerical tests since it is not upwinded in fast and slow waves.We propose a method of magnetosonic characteristics to evolve compressional waves.The tests indicate that this method cures the "rarefaction shocks" problem to a large extent and it also greatly reduces some post shock oscillations.
ZHU XiaoShuaiWANG HuaNingFAN YuLiangDU ZhanLeHE Han
