BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xueshang Feng

Xueshang Feng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRphysics.space-phphysics.plasm-phastro-ph.IMphysics.comp-ph

Catalog footprint

What is connected

26works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

Homologous Coronal Mass Ejections Caused by Recurring Formation and Disruption of Current Sheet within a Sheared Magnetic Arcade

The Sun often produces coronal mass ejections with similar structure repeatedly from the same source region, and how these homologous eruptions are initiated remains an open question. Here, by using a new magnetohydrodynamic simulation, we show that homologous solar eruptions can be efficiently produced by recurring formation and disruption of coronal current sheet as driven by continuously shearing of the same polarity inversion line within a single bipolar configuration. These eruptions are initiated by the same mechanism, in which an internal current sheet forms slowly in a gradually sheared bipolar field and reconnection of the current sheet triggers and drives the eruption. Each of the eruptions does not release all the free energy but with a large amount left in the post-flare arcade below the erupting flux rope. Thus, a new current sheet can be more easily formed by further shearing of the post-flare arcade than by shearing a potential field arcade, and this is favorable for producing the next eruption. Furthermore, it is found that the new eruption is stronger since the newly formed current sheet has a larger current density and a lower height. In addition, our results also indicate the existence of a magnetic energy threshold for a given flux distribution, and eruption occurs once this threshold is approached.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Simulation of Solar Magnetic Flux Emergence Using the AMR--CESE--MHD Code

Magnetic flux emergence from the solar interior to the atmosphere is believed to be a key process of formation of solar active regions and driving solar eruptions. Due to the limited capability of observation, the flux emergence process is commonly studied using numerical simulations. In this paper, we developed a numerical model to simulate the emergence of a twisted magnetic flux tube from the convection zone to the corona using the AMR--CESE--MHD code, which is based on the conservation-element solution-element method with adaptive mesh refinement. The result of our simulation agrees with that of many previous ones with similar initial conditions but using different numerical codes. In the early stage, the flux tube rises from the convection zone as driven by the magnetic buoyancy until it reaches close to the photosphere. The emergence is decelerated there and with piling-up of the magnetic flux, the magnetic buoyancy instability is triggered, which allows the magnetic field to partially enter into the atmosphere. Meanwhile, two gradually separated polarity concentration zones appear in the photospheric layer, transporting the magnetic field and energy into the atmosphere through their vortical and shearing motions. Correspondingly, the coronal magnetic field has also been reshaped to a sigmoid configuration containing a thin current layer, which resembles the typical pre-eruptive magnetic configuration of an active region. Such a numerical framework of magnetic flux emergence as established will be applied in future investigations of how solar eruptions are initiated in flux emergence active regions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The role of photospheric converging motion in initiation of solar eruptions

It is well known that major solar eruptions are often produced by active regions with continual photospheric shearing and converging motions. Here, through high accuracy magnetohydrodynamics simulation, we show how solar eruption is initiated in a single bipolar configuration as driven by first shearing and then converging motions at the bottom surface. Different from many previous simulations, we applied the converging motion without magnetic diffusion, thus it only increases the magnetic gradient across the polarity inversion line but without magnetic flux cancellation. The converging motion at the footpoints of the sheared arcade creates a current sheet in a quasi-static way, and the eruption is triggered by magnetic reconnection of the current sheet, which supports the same scenario as shown in our previous simulation with only shearing motion. With the converging motion, the current sheet is formed at a lower height and has a higher current density than with shearing motion alone, which makes reconnection more effective and eruption stronger. Moreover, the converging motion renders a fast decay rate of the overlying field with height and thus favorable for an eruption. This demonstrate that the converging flow is more efficient to create the current sheet and more favorable for eruption than by solely the shearing flow.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Simulation of a Fundamental Mechanism of Solar Eruption with Different Magnetic Flux Distributions

Solar eruptions are explosive release of coronal magnetic field energy as manifested in solar flares and coronal mass ejection. Observations have shown that the core of eruption-productive regions are often a sheared magnetic arcade, i.e., a single bipolar configuration, and, particularly, the corresponding magnetic polarities at the photosphere are elongated along a strong-gradient polarity inversion line (PIL). It remains unclear what mechanism triggers the eruption in a single bipolar field and why the one with a strong PIL is eruption-productive. Recently, using high accuracy simulations, we have established a fundamental mechanism of solar eruption initiation that a bipolar field as driven by quasi-static shearing motion at the photosphere can form an internal current sheet, and then fast magnetic reconnection triggers and drives the eruption. Here we investigate the behavior of the fundamental mechanism with different photospheric magnetic flux distributions, i.e., magnetograms, by combining theoretical analysis and numerical simulation. Our study shows that the bipolar fields of different magnetograms, as sheared continually, all exhibit similar evolutions from the slow storage to fast release of magnetic energy in accordance with the fundamental mechanism, which demonstrates the robustness of the mechanism. We further found that the magnetograms with stronger PIL produce larger eruptions, and the key reason is that the sheared bipolar fields with stronger PIL can achieve more non-potentiality, and their internal current sheet can form at a lower height and with a larger current density, by which the reconnection can be more efficient. This also provides a viable trigger mechanism for the observed eruptions in active region with strong PIL.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Continuous Null-Point Magnetic Reconnection Builds Up a Torus Unstable Magnetic Flux Rope Triggering the X9.3 Flare in Solar Active Region~12673

Two X-class solar flares occurred on 2017 September 6 from active region NOAA 12673: the first one is a confined X2.2 flare, and it is followed only $\sim 3$ hours later by the second one, which is the strongest flare in solar cycle 24, reaching X9.3 class and accompanied with a coronal mass ejection. Why these two X-class flares occurred in the same position with similar magnetic configurations, but one is eruptive while the other is not? Here we track the coronal magnetic field evolution via nonlinear force-free field extrapolations from a time sequence of vector magnetograms with high cadence. A detailed analysis of the magnetic field shows that a magnetic flux rope (MFR) forms and grows gradually before the first flare, and shortly afterwards, the MFR's growth is significantly enhanced with a much faster rise in height, from far below the threshold of torus instability to above it, while the magnetic twist only increases mildly. Combining EUV observations and the magnetic field extrapolation, we found that overlying the MFR is a null-point magnetic topology, where recurrent brightening is seen after the first flare. We thus suggest a scenario to interpret the occurrence of the two flares. The first flare occurred since the MFR reached a high enough height to activate the null point, and its continuous expansion forces the null-point reconnection recurrently. Such reconnection weakens the overlying field, allowing the MFR to rise faster, which eventually crosses the threshold of torus instability and triggers the second, eruptive flare.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-driven MHD Simulation of the Formation and Initiation of a Large-scale Pre-flare Magnetic Flux Rope in Solar Active Region 12371

Solar eruptions are the most powerful drivers of space weather. To understand their cause and nature, it is crucial to know how the coronal magnetic field evolves before eruption. Here we study the formation process of a relatively large-scale magnetic flux rope (MFR) in active region NOAA~12371 that erupts with a major flare and coronal mass ejection on 2015 June 21. A data-driven numerical magnetohydrodynamic model is employed to simulate three-dimensional coronal magnetic field evolution of one-day duration before the eruption. Comparison between the observed features and our modeled magnetic field discloses how the pre-eruption MFR forms. Initially, the magnetic field lines were weakly twisted as being simple sheared arcades. Then a long MFR was formed along the polarity inversion line due to the complex photospheric motion, which is mainly shearing rather than twisting. The presence of the MFR is evidenced by a coherent set of magnetic field lines with twist number above unity. Below the MFR a current sheet is shown in the model, suggesting that tether-cutting reconnection plays a key role in the MFR formation. The MFR's flux grows as more and more field lines are twisted due to continuous injection of magnetic helicity by the photospheric motions. Meanwhile, the height of the MFR's axis increases monotonely from its formation. By an analysis of the decay index of its overlying field, we suggest that it is because the MFR runs into the torus instability regime and becomes unstable that finally triggers the eruption.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Numerical Investigation of the Recurrent High-speed Jets as a Possibility of Solar Wind Origin

In the solar atmosphere, jets are prevalent and they are significant for the mass and energy transport. Here we conduct numerical simulations to investigate the mass and energy contributions of the recently observed high-speed jets to the solar wind. With a one-dimensional hydrodynamic solar wind model, the time-dependent pulses are imposed at the bottom to simulate the jets. The simulation results show that without other energy source, the injected plasmas are accelerated effectively to be a transonic wind with a substantial mass flux. The rapid acceleration occurs close to the Sun, and the resulting asymptotic speed, number density at 0.3 AU, as well as mass flux normalized to 1 AU are compatible with in situ observations. As a result of the high speed, the imposed pulses generate a train of shocks traveling upward. By tracing the motions of the injected plasma, it is found that these shocks heat and accelerate the injected plasmas successively step by step to push them upward and eventually allow them to escape. The parametric studies show that increasing the speed of the imposed pulses or their temperature gives a considerably faster, and hotter solar wind, while increasing their number density or decreasing their recurring period only bring a denser solar wind. These studies provide a possibility that the ubiquitous high-speed jets are a substantial mass and energy contributions to the solar wind.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparison of CME/shock propagation models with heliospheric imaging and in situ observations

The prediction of the arrival time for fast coronal mass ejections (CMEs) and their associated shocks is highly desirable in space weather studies. In this paper, we use two shock propagation models, i.e. Data Guided Shock Time Of Arrival (DGSTOA) and Data Guided Shock Propagation Model (DGSPM), to predict the kinematical evolution of interplanetary shocks associated with fast CMEs. DGSTOA is based on the similarity theory of shock waves in the solar wind reference frame, and DGSPM on the non-similarity theory in the stationary reference frame. The inputs are the kinematics of the CME front at the maximum speed moment obtained from the geometric triangulation method applied to STEREO imaging observations together with the Harmonic Mean approximation. The outputs provide the subsequent propagation of the associated shock. We apply these models to the CMEs on 2012 January 19, January 23, and March 7. We find that the shock models predict reasonably well the shock's propagation after the impulsive acceleration. The shock's arrival time and local propagation speed at Earth predicted by these models are consistent with in situ measurements of WIND. We also employ the Drag-Based Model (DBM) as a comparison, and find that it predicts a steeper deceleration than the shock models after the rapid deceleration phase. The predictions of DBM at 1 AU agree with the following ICME or sheath structure, not the preceding shock. These results demonstrate the applicability of the shock models used here for future arrival time prediction of interplanetary shocks associated with fast CMEs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dark Post-Flare Loops Observed by Solar Dynamics Observatory

Solar post-flare loops (PFLs) are arcade-like loop systems that appear during the gradual phases of eruptive flares. The extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) allow us to investigate the fine structures in PFLs. In this work, we focus on studying the dark post-flare loops (DPFLs) during X-class flares, which are more evident in SDO/AIA data than in previous EUV data. We identify and analyze the DPFLs observed by SDO and find that: (1) the DPFLs of an X5.4 flare have an average lifetime of 10.0 $\pm$ 5.5 minutes, an average width of 1022 $\pm$ 339 km, and an average maximum length of 33 $\pm$ 10 Mm, (2) blob-like falling features with a size close to the resolution of SDO/AIA are identified in the DPFLs and they have an average velocity of 76 $\pm$ 19 km s$^{-1}$, and (3) the average widths of the DPFLs slightly increase with the characteristic temperatures in AIA 304, 171, 193, and 211 Å channels. Our investigation shows that DPFLs are found in all of the 20 cases within this study, which suggests that they are a common phenomenon in X-class flares and are probably produced by the same mechanism that creates coronal rain.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How did a Major Confined Flare Occur in Super Solar Active Region 12192?

We study the physical mechanism of a major X-class solar flare that occurred in the super NOAA active region (AR) 12192 using a data-driven numerical magnetohydrodynamic (MHD) modeling complemented with observations. With the evolving magnetic fields observed at the solar surface as bottom boundary input, we drive an MHD system to evolve self-consistently in correspondence with the realistic coronal evolution. During a two-day time interval, the modeled coronal field has been slowly stressed by the photospheric field evolution,which gradually created a large-scale coronal current sheet, i.e., a narrow layer with intense current, in the core of the AR. The current layer was successively enhanced until it became so thin that a tether-cutting reconnection between the sheared magnetic arcades was set in, which led to a flare. The modeled reconnecting field lines and their footpoints match well the observed hot flaring loops and the flare ribbons, respectively, suggesting that the model has successfully "reproduced" the macroscopic magnetic process of the flare. In particular, with simulation, we explained why this event is a confined eruption-the consequent of the reconnection is the shared arcade instead of a newly formed flux rope. We also found much weaker magnetic implosion effect comparing to many other X-class flares

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Origin of Both the Fast Hot Jet and the Slow Cool Jet from Magnetic Flux Emergence and Advection in the Solar Transition Region

In the solar atmosphere, the jets are ubiquitous and found to be at various spatia-temporal scales. They are significant to understand energy and mass transport in the solar atmosphere. Recently, the high-speed transition region jets are reported from the observation. Here we conduct a numerical simulation to investigate the mechanism in their formation. Driven by the supergranular convection motion, the magnetic reconnection between the magnetic loop and the background open flux occurring in the transition region is simulated with a two-dimensional magnetohydrodynamics model. The simulation results show that not only a fast hot jet, much resemble the found transition region jets, but also a adjacent slow cool jet, mostly like classical spicules, is launched. The force analysis shows that the fast hot jet is continually driven by the Lorentz force around the reconnection region, while the slow cool jet is induced by an initial kick through the Lorentz force associated with the emerging magnetic flux. Also, the features of the driven jets change with the amount of the emerging magnetic flux, giving the varieties of both jets. These results will inspire our understanding of the formation of the prevalence of both the fast hot jet and slow cool jet from the solar transition region and chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Comparison Study of a Solar Active-Region Eruptive Filament and a Neighboring Non-Eruptive Filament

Solar active region (AR) 11283 is a very magnetically complex region and it has produced many eruptions. However, there exists a non-eruptive filament in the plage region just next to an eruptive one in the AR, which gives us an opportunity to perform a comparison analysis of these two filaments. The coronal magnetic field extrapolated using a CESE-MHD-NLFFF code (Jiang & Feng 2013) reveals that two magnetic flux ropes (MFRs) exist in the same extrapolation box supporting these two filaments, respectively. Analysis of the magnetic field shows that the eruptive MFR contains a bald-patch separatrix surface (BPSS) co-spatial very well with a pre-eruptive EUV sigmoid, which is consistent with the BPSS model for coronal sigmoids. The magnetic dips of the non-eruptive MFRs match Hα observation of the non-eruptive filament strikingly well, which strongly supports the MFR-dip model for filaments. Compared with the non-eruptive MFR/filament (with a length of about 200 Mm), the eruptive MFR/filament is much smaller (with a length of about 20 Mm), but it contains most of the magnetic free energy in the extrapolation box and holds a much higher free energy density than the non-eruptive one. Both the MFRs are weakly twisted and cannot trigger kink instability. The AR eruptive MFR is unstable because its axis reaches above a critical height for torus instability, at which the overlying closed arcades can no longer confine the MFR stably. On the contrary, the quiescent MFR is very firmly held by its overlying field, as its axis apex is far below the torus-instability threshold height. Overall, this comparison investigation supports that MFR can exist prior to eruption and the ideal MHD instability can trigger MFR eruption.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evolution of the 2012 July 12 CME from the Sun to the Earth: Data-Constrained Three-Dimensional MHD Simulations

The dynamic process of coronal mass ejections (CMEs) in the heliosphere provides us the key information for evaluating CMEs' geo-effectiveness and improving the accurate prediction of CME induced Shock Arrival Time (SAT) at the Earth. We present a data constrained three dimensional (3D) magnetohydrodynamic (MHD) simulation of the evolution of the CME in a realistic ambient solar wind for the July 12-16, 2012 event by using the 3D COIN-TVD MHD code. A detailed comparison of the kinematic evolution of the CME between the observations and the simulation is carried out, including the usage of the time-elongation maps from the perspectives of both Stereo A and Stereo B. In this case study, we find that our 3D COIN-TVD MHD model, with the magnetized plasma blob as the driver, is able to re-produce relatively well the real 3D nature of the CME in morphology and their evolution from the Sun to Earth. The simulation also provides a relatively satisfactory comparison with the in-situ plasma data from the Wind spacecraft.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of Rotational Discontinuities in Compressive three-dimensional MHD Turbulence

Measurements of solar wind turbulence reveal the ubiquity of discontinuities. In this study, we investigate how the discontinuities, especially rotational discontinuities (RDs), are formed in magnetohydrodynamic (MHD) turbulence. In a simulation of the decaying compressive three-dimensional (3-D) MHD turbulence with an imposed uniform background magnetic field, we detect RDs with sharp field rotations and little variations of magnetic field intensity as well as mass density. At the same time, in the de Hoffman-Teller (HT) frame, the plasma velocity is nearly in agreement with the Alfvén speed, and is field-aligned on both sides of the discontinuity. We take one of the identified RDs to analyze in details its 3-D structure and temporal evolution. By checking the magnetic field and plasma parameters, we find that the identified RD evolves from the steepening of the Alfvén wave with moderate amplitude, and that steepening is caused by the nonuniformity of the Alfvén speed in the ambient turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing a Solar Coronal Magnetic Field Extrapolation Code with the Titov-Demoulin Magnetic Flux Rope Model

In the solar corona, magnetic flux rope is believed to be a fundamental structure accounts for magnetic free energy storage and solar eruptions. Up to the present, the extrapolation of magnetic field from boundary data is the primary way to obtain fully three-dimensional magnetic information of the corona. As a result, the ability of reliable recovering coronal magnetic flux rope is important for coronal field extrapolation. In this paper, our coronal field extrapolation code (CESE-MHD-NLFFF, Jiang & Feng 2012) is examined with an analytical magnetic flux rope model proposed by Titov & Demoulin (1999), which consists of a bipolar magnetic configuration holding an semi-circular line-tied flux rope in force-free equilibrium. By using only the vector field in the bottom boundary as input, we test our code with the model in a representative range of parameter space and find that the model field is reconstructed with high accuracy. Especially, the magnetic topological interfaces formed between the flux rope and the surrounding arcade, i.e., the "hyperbolic flux tube" and "bald patch separatrix surface", are also reliably reproduced. By this test, we demonstrate that our CESE-MHD-NLFFF code can be applied to recovering magnetic flux rope in the solar corona as long as the vector magnetogram satisfies the force-free constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Could the collision of CMEs in the heliosphere be super-elastic? --- Validation through three-dimensional simulations

Though coronal mass ejections (CMEs) are magnetized fully-ionized gases, a recent observational study of a CME collision event in 2008 November has suggested that their behavior in the heliosphere is like elastic balls, and their collision is probably super-elastic \citep{Shen_etal_2012}. If this is true, this finding has an obvious impact on the space weather forecasting because the direction and veliocity of CMEs may change. To verify it, we numerically study the event through three-dimensional MHD simulations. The nature of CMEs' collision is examined by comparing two cases. In one case the two CMEs collide as observed, but in the other, they do not. Results show that the collision leads to extra kinetic energy gain by 3%--4% of the initial kinetic energy of the two CMEs. It firmly proves that the collision of CMEs could be super-elastic.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear Force-Free Field Extrapolation of a Coronal Magnetic Flux Rope Supporting a Large-Scale Filament from Photospheric Vector Magnetogram

Solar filament are commonly thought to be supported in magnetic dips, in particular, of magnetic flux ropes (FRs). In this Letter, from the observed photospheric vector magnetogram, we implement a nonlinear force-free field (NLFFF) extrapolation of a coronal magnetic FR that supports a large-scale intermediate filament between an active region and a weak polarity region. This result is the first in that current NLFFF extrapolations with presence of FRs are limited to relatively small-scale filaments that are close to sunspots and along main polarity inversion line (PIL) with strong transverse field and magnetic shear, and the existence of a FR is usually predictable. In contrast, the present filament lies along the weak-field region (photospheric field strength $\lesssim 100$ G), where the PIL is very fragmented due to small parasitic polarities on both side of the PIL and the transverse field has a low value of signal-to-noise ratio. Thus it represents a far more difficult challenge to extrapolate a large-scale FR in such case. We demonstrate that our CESE--MHD--NLFFF code is competent for the challenge. The numerically reproduced magnetic dips of the extrapolated FR match observations of the filament and its barbs very well, which supports strongly the FR-dip model for filaments. The filament is stably sustained because the FR is weakly twisted and strongly confined by the overlying closed arcades.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extrapolation of the Solar Coronal Magnetic Field from SDO/HMI Magnetogram by a CESE--MHD--NLFFF Code

Due to the absence of direct measurement, the magnetic field in the solar corona is usually extrapolated from the photosphere in numerical way. At the moment, the nonlinear force-free field (NLFFF) model dominates the physical models for field extrapolation in the low corona. Recently we have developed a new NLFFF model with MHD relaxation to reconstruct the coronal magnetic field. This method is based on CESE--MHD model with the conservation-element/solution-element (CESE) spacetime scheme. In this paper, we report the application of the CESE--MHD--NLFFF code to \SDO/HMI data with magnetograms sampled for two active regions (ARs), NOAA AR 11158 and 11283, both of which were very non-potential, producing X-class flares and eruptions. The raw magnetograms are preprocessed to remove the force and then inputted into the extrapolation code. Qualitative comparison of the results with the \SDO/AIA images shows that our code can reconstruct magnetic field lines resembling the EUV-observed coronal loops. Most important structures of the active regions are reproduced excellently, like the highly-sheared field lines that suspend filaments in AR 11158 and twisted flux rope which corresponds to a sigmoid in AR 11283. Quantitative assess of the results shows that the force-free constraint is fulfilled very well in the strong-field regions but apparently not that well in the weak-field regions because of data noise and numerical errors in the small currents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation and Eruption of an Active Region Sigmoid: I. Study by Nonlinear Force-Free Field Modeling

We present a magnetic analysis of the formation and eruption of an active region (AR) sigmoid in AR 11283 from 2011 September 4 to 6. To follow the quasi-static evolution of the coronal magnetic field, we reconstruct a time sequence of static fields using a recently developed nonlinear force-free field model constrained by the SDO/HMI vector magnetograms. A detailed analysis of the fields compared with the SDO/AIA observations suggests the following scenario for the evolution of the region. Initially, a new bipole emerges into the negative polarity of a pre-existing bipolar AR, forming a null point topology between the two flux systems. A weakly twisted flux rope (FR) is then built up slowly in the embedded core region, largely through flux-cancellation photospheric reconnections, forming a bald patch separatrix surface (BPSS) separating the FR from its ambient field. The FR grows gradually until its axis runs into a torus instability (TI) domain near the end of the third day, and the BPSS also develops a fully S-shape. Unlike in the case of standard TI, the FR does not erupt instantly since it is still attached at the photosphere along the bald patch (BP) portion of the polarity inversion line. The combined effects of the TI-driven expansion of the FR and the line-tying at the BP tear the FR into two parts with the upper portion freely expelled and the lower portion remaining behind the post-flare arcades. This process dynamically perturbs the BPSS and results in the transient enhanced brightening of the sigmoid. The accelerated expansion of the upper portion of the FR strongly pushes its envelope flux near the null point and triggers breakout reconnection at the null, as evidenced by a remarkable circular flare ribbon, which further facilitates the eruption. We discuss the important implications of these results for the formation and disruption of sigmoid region with FR.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

MHD Simulation of a Sigmoid Eruption of Active Region 11283

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in a MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Preprocess the Photospheric Vector Magnetograms for NLFFF Extrapolation using a Potential Field Model and an Optimization Method

Numerical reconstruction/extrapolation of coronal nonlinear force-free magnetic field (NLFFF) usually takes the photospheric vector magnetogram as input at the bottom boundary. Magnetic field observed at the photosphere, however, contains force which is in conflict with the fundamental assumption of the force-free model and measurement noise which is unfavorable for practical computation. Preprocessing of the raw magnetogram has been proposed by Wiegelmann, Inhester, and Sakurai (2006) to remove the force and noise for providing better input for NLFFF modeling. In this paper we develop a new code of magnetogram preprocessing which is consistent with our extrapolation method CESE-MHD-NLFFF (Jiang, Feng, and Xiang, 2012; Jiang and Feng, 2012). Basing on a magnetic-splitting rule that a magnetic field can be split into a potential field part and a non-potential part, we split the magnetogram and deal with the two parts separately. Preprocessing of the magnetogram's potential part is based on a numerical potential field model, and the non-potential part is preprocessed using the similar optimization method of Wiegelmann et al (2006). The code is applied to the SDO/HMI data and results show that the method can remove efficiently the force and noise and improve the quality of extrapolation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Using Coordinated Observations in Polarised White Light and Faraday Rotation to Probe the Spatial Position and Magnetic Field of an Interplanetary Sheath

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White-light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons, and therefore depend on both the viewing geometry and the electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component $B_\parallel$, and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modelling of an Earth-directed shock and synthesise the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance $I$ decreases with heliocentric distance $r$ roughly according to the expression $I \propto r^{-3}$. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site $\mathbf{r}_{\rm sheath}$ and the mass of plasma at that position $M_{\rm sheath}$ can be inferred from the polarisation of the shock-associated enhancement in WL radiance. From the FR measurements, the local $B_{\parallel {\rm sheath}}$ at $\mathbf{r}_{\rm sheath}$ can then be estimated. Simultaneous observations in polarised WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Code for Nonlinear Force-Free Field Extrapolation of the Global Corona

Reliable measurements of the solar magnetic field are still restricted to the photosphere, and our present knowledge of the three-dimensional coronal magnetic field is largely based on extrapolation from photospheric magnetogram using physical models, e.g., the nonlinear force-free field (NLFFF) model as usually adopted. Most of the currently available NLFFF codes have been developed with computational volume like Cartesian box or spherical wedge while a global full-sphere extrapolation is still under developing. A high-performance global extrapolation code is in particular urgently needed considering that Solar Dynamics Observatory (SDO) can provide full-disk magnetogram with resolution up to $4096\times 4096$. In this work, we present a new parallelized code for global NLFFF extrapolation with the photosphere magnetogram as input. The method is based on magnetohydrodynamics relaxation approach, the CESE-MHD numerical scheme and a Yin-Yang spherical grid that is used to overcome the polar problems of the standard spherical grid. The code is validated by two full-sphere force-free solutions from Low & Lou's semi-analytic force-free field model. The code shows high accuracy and fast convergence, and can be ready for future practical application if combined with an adaptive mesh refinement technique.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Implementation of the Magnetohydrodynamics-Relaxation Method for Nonlinear Force-Free Field Extrapolation in the Solar Corona

Magnetic field in the solar corona is usually extrapolated from photospheric vector magnetogram using a nonlinear force-free field (NLFFF) model. NLFFF extrapolation needs a considerable effort to be devoted for its numerical realization. In this paper we present a new implementation of the magnetohydrodynamics (MHD)-relaxation method for NLFFF extrapolation. The magneto-frictional approach which is introduced for speeding the relaxation of the MHD system is novelly realized by the spacetime conservation-element and solution-element (CESE) scheme. A magnetic field splitting method is used to further improve the computational accuracy. The bottom boundary condition is prescribed by changing the transverse field incrementally to match the magnetogram, and all other artificial boundaries of the computational box are simply fixed. We examine the code by two types of NLFFF benchmark tests, the Low & Lou (1990) semi-analytic force-free solutions and a more realistic solar-like case constructed by van Ballegooijen et al. (2007). The results show that our implementation are successful and versatile for extrapolations of either the relatively simple cases or the rather complex cases which need significant rebuilding of the magnetic topology, e.g., a flux rope. We also compute a suite of metrics to quantitatively analyze the results and demonstrate that the performance of our code in extrapolation accuracy basically reaches the same level of the present best-performing code, e.g., that developed by Wiegelmann (2004).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Study of the 3D Coronal Magnetic Field of Active Region 11117 Around the Time of a Confined Flare Using a Data-Driven CESE--MHD Model

We apply a data-driven MHD model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element (CESE) scheme, is designed to focus on the magnetic-field evolution and to consider a simplified solar atomsphere with finite plasma $β$. Magnetic vector-field data derived from the observations at the photoshpere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria basing on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager (HMI) on board the {\it Solar Dynamic Observatory (SDO)} around the time of flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and the coronal loops observed by the Atmospheric Imaging Assembly (AIA), which shows that the coronal field can indeed be well characterized by the MHD equilibrium in most time. The magnetic configuration changes very limited during the studied time interval of two hours. A topological analysis reveals that the small flare is correlated with a bald patch (BP, where the magnetic field is tangent to the photoshpere), suggesting that the energy release of the flare can be understood by magnetic reconnection associated with the BP separatrices. The total magnetic flux and energy keep increasing slightly in spite of the flare, while the computed magnetic free energy drops during the flare with an amount of $\sim 10^{30}$ erg, which seems to be adequate to provide the energy budget of the minor C-class confined flare.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient magnetohydrodynamic simulations on graphics processing units with CUDA

Magnetohydrodynamic (MHD) simulations based on the ideal MHD equations have become a powerful tool for modeling phenomena in a wide range of applications including laboratory, astrophysical, and space plasmas. In general, high-resolution methods for solving the ideal MHD equations are computationally expensive and Beowulf clusters or even supercomputers are often used to run the codes that implemented these methods. With the advent of the Compute Unified Device Architecture (CUDA), modern graphics processing units (GPUs) provide an alternative approach to parallel computing for scientific simulations. In this paper we present, to the author's knowledge, the first implementation of MHD simulations entirely on GPUs with CUDA, named GPU-MHD, to accelerate the simulation process. GPU-MHD supports both single and double precision computation. A series of numerical tests have been performed to validate the correctness of our code. Accuracy evaluation by comparing single and double precision computation results is also given. Performance measurements of both single and double precision are conducted on both the NVIDIA GeForce GTX 295 (GT200 architecture) and GTX 480 (Fermi architecture) graphics cards. These measurements show that our GPU-based implementation achieves between one and two orders of magnitude depending on the used graphics card, problem size, and precision when comparing to the original serial CPU MHD implementation. In addition, we extend GPU-MHD to support the visualization of the simulation results and thus the whole MHD simulation and visualization process can be performed entirely on GPUs.

0 citations0 reviewsNo code linkedOpen access