Data-constrained Simulations Research Articles

Data-constrained Magnetohydrodynamic Simulation of a Long-duration Eruptive Flare

Abstract We perform a zero-β magnetohydrodynamic simulation for the C7.7 class flare initiated at 01:18 UT on 2011 June 21 using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The initial condition for the simulation involves a flux rope, which we realize through the regularized Biot–Savart laws, whose parameters are constrained by observations from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) and the Extreme Ultraviolet Imager (EUVI) on the twin Solar Terrestrial Relations Observatory (STEREO). This data-constrained initial state is then relaxed to a force-free state by the magnetofrictional module in MPI-AMRVAC. The further time-evolving simulation results reproduce the eruption characteristics obtained by SDO/AIA 94 Å, 304 Å, and STEREO/EUVI 304 Å observations fairly well. The simulated flux rope possesses similar eruption direction, height range, and velocity to the observations. In particular, the two phases of slow evolution and fast eruption are reproduced by varying the density distribution in the light of the draining process of the filament material. Our data-constrained simulations also show other advantages, such as a large field of view (about 0.76 R ⊙). We study the twist of the magnetic flux rope and the decay index of the overlying field, and find that in this event, both the magnetic strapping force and the magnetic tension force are sufficiently weaker than the magnetic hoop force, thus allowing the successful eruption of the flux rope. We also find that the anomalous resistivity is necessary to keep the correct morphology of the erupting flux rope.

Radial transport modeling of high density deuterium plasmas in proto-MPEX with the B2.5-Eirene code

High density (≥6 × 1019 m−3), low temperature (2–6 eV) helicon discharges in the Prototype Material Plasma Exposure eXperiment (Proto-MPEX) are analyzed with the coupled multifluid plasma, kinetic neutrals code B2.5-Eirene. The interpretative analyses are constrained by data from multiple diagnostics, including Langmuir probes, Mach probes, filterscopes, infrared TV system, Thomson scattering, and baratrons. The objectives of the transport simulations include: investigation of the effects of heating, fueling, and plasma production; pumping, and assumed radial transport models on the calculated density and temperature distributions; plasma flow profiles and power balance. The primary objective in this report is to investigate the effects of the radial transport model in full plasma (the entire length of the plasma column in Proto-MPEX) data-constrained simulations. Results from three assumed forms of the radial transport coefficients are presented, including spatially constant, radially decreasing, and Bohm (D,χ ∼ Te/|B|). The results from each of the three transport coefficient sets agree qualitatively with the core (near axis) data. With the implicit Te dependence, the Bohm coefficients tend to decrease as functions of radius, although not as strongly as the centrally peaked set. The axial variation in the Bohm coefficients is largely due to the axial structure of the magnetic field. The agreement of the simulations and the diagnostic data with the Bohm set indicates that transport in the plasma column of Proto-MPEX is dominated by Bohm diffusion.High density (≥6 × 1019 m−3), low temperature (2–6 eV) helicon discharges in the Prototype Material Plasma Exposure eXperiment (Proto-MPEX) are analyzed with the coupled multifluid plasma, kinetic neutrals code B2.5-Eirene. The interpretative analyses are constrained by data from multiple diagnostics, including Langmuir probes, Mach probes, filterscopes, infrared TV system, Thomson scattering, and baratrons. The objectives of the transport simulations include: investigation of the effects of heating, fueling, and plasma production; pumping, and assumed radial transport models on the calculated density and temperature distributions; plasma flow profiles and power balance. The primary objective in this report is to investigate the effects of the radial transport model in full plasma (the entire length of the plasma column in Proto-MPEX) data-constrained simulations. Results from three assumed forms of the radial transport coefficients are presented, including spatially constant, radially decreasing, and Boh...

Magnetohydrodynamic Simulations for Studying Solar Flare Trigger Mechanism

Abstract In order to understand the flare trigger mechanism, we conduct three-dimensional magnetohydrodynamic simulations using a coronal magnetic field model derived from data observed by the Hinode satellite. Several types of magnetic bipoles are imposed into the photospheric boundary of the Nonlinear Force-free Field model of Active Region (AR) NOAA 10930 on 2006 December 13, to investigate what kind of magnetic disturbance may trigger the flare. As a result, we confirm that certain small bipole fields, which emerge into the highly sheared global magnetic field of an AR, can effectively trigger a flare. These bipole fields can be classified into two groups based on their orientation relative to the polarity inversion line: the so-called opposite polarity, and reversed shear structures, as suggested by Kusano et al. We also investigate the structure of the footpoints of reconnected field lines. By comparing the distribution of reconstructed field lines and observed flare ribbons, the trigger structure of the flare can be inferred. Our simulation suggests that the data-constrained simulation, taking into account both the large-scale magnetic structure and small-scale magnetic disturbance (such as emerging fluxes), is a good way to discover a flare-producing AR, which can be applied to space weather prediction.