Source author record

Chaowei Jiang

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

25 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.

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.

preprint2022arXiv

Quantitative Characterization of Magnetic Flux Rope Properties for Two Solar Eruption Events

In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR due to the lack of direct magnetic field measurement in the corona. In this study, we investigate the MFR structures, originating from two active regions (ARs), AR 11719 and AR 12158, and estimate their magnetic properties quantitatively. We perform the nonlinear force-free field extrapolations with preprocessed photospheric vector magnetograms. In addition, remote-sensing observations are employed to find indirect evidence of MFRs on the Sun and to analyze the time evolution of magnetic reconnection flux associated with the flare ribbons during the eruption. A coherent "pre-existing" MFR structure prior to the flare eruption is identified quantitatively for one event from the combined analysis of the extrapolation and observation. Then the characteristics of MFRs for two events on the Sun before and during the eruption, forming the CME-MFR, including the axial magnetic flux, field-line twist, and reconnection flux, are estimated and compared with the corresponding in situ modeling results. We find that the magnetic reconnection associated with the accompanying flares for both events injects significant amount of flux into the erupted CME-MFRs.

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.

preprint2022arXiv

The Rotation of Magnetic Flux Rope Formed during Solar Eruption

The eruptions of solar filaments often show rotational motion about their rising direction, but it remains elusive what mechanism governs such rotation and how the rotation is related to the initial morphology of the pre-eruptive filament (and co-spatial sigmoid), filament chirality, and magnetic helicity. The conventional view regarding the rotation as a result of a magnetic flux rope (MFR) under-going the ideal kink instability still has confusion in explaining these relationships. Here we proposed an alternative explanation for the rotation during eruptions, by analyzing a magnetohydrodynamic simulation in which magnetic reconnection initiates an eruption from a sheared arcade configuration and an MFR is formed during eruption through the reconnection. The simulation reproduces a reverse S-shaped MFR with dextral chirality, and the axis of this MFR rotates counterclockwise while rising, which compares favorably with a typical filament eruption observed from dual viewing angles. By calculating the twist and writhe numbers of the modeled MFR during its eruption, we found that accompanied with the rotation, the nonlocal writhe of the MFR's axis decreases while the twist of its surrounding field lines increases, and this is distinct from the kink instability, which converts magnetic twist into writhe of the MFR axis.

preprint2021arXiv

Diagnosing a Solar Flaring Core with Bidirectional Quasi-Periodic Fast Propagating Magnetoacoustic Waves

Quasi-periodic fast propagating (QFP) waves are often excited by solar flares, and could be trapped in the coronal structure with low Alfvén speed, so they could be used as a diagnosing tool for both the flaring core and magnetic waveguide. As the periodicity of a QFP wave could originate from a periodic source or be dispersively waveguided, it is a key parameter for diagnosing the flaring core and waveguide. In this paper, we study two QFP waves excited by a GOES-class C1.3 solar flare occurring at active region NOAA 12734 on 8 March 2019. Two QFP waves were guided by two oppositely oriented coronal funnel. The periods of two QFP waves were identical and were roughly equal to the period of the oscillatory signal in the X-ray and 17 GHz radio emission released by the flaring core. It is very likely that the two QFP waves could be periodically excited by the flaring core. Many features of this QFP wave event is consistent with the magnetic tuning fork model. We also investigated the seismological application with QFP waves, and found that the magnetic field inferred with magnetohydrodynamic seismology was consistent with that obtained in magnetic extrapolation model. Our study suggest that the QFP wave is a good tool for diagnosing both the flaring core and the magnetic waveguide.

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.

preprint2021arXiv

Two Quasi-periodic Fast-propagating Magnetosonic Wave Events Observed In Active Region NOAA 11167

We report a detailed observational study of two quasi-periodic fast-propagating (QFP) magnetosonic wave events occurred on 2011 March 09 and 10, respectively. Interestingly, both the two events have two wave trains (WTs): one main and strong (WT-1) whereas the second appears small and weak (WT-2). Peculiar and common characteristics of the two events are observed, namely: 1) the two QFP waves are accompanied with brightenings during the whole stage of the eruptions; 2) both the two main wave trains are nearly propagating along the same direction; 3) EUV waves are found to be associated with the two events. Investigating various aspects of the target events, we argue that: 1) the second event is accompanied with a flux rope eruption during the whole stage; 2) the second event eruption produces a new filament-like (FL) dark feature; 3) the ripples of the two WT-2 QFP waves seem to result from different triggering mechanisms. Based on the obtained observational results, we propose that the funnel-like coronal loop system is indeed playing an important role in the two WT-1 QFP waves. The development of the second WT-2 QFP wave can be explained as due to the dispersion of the main EUV front. The co-existence of the two events offer thereby a significant opportunity to reveal what driving mechanisms and structures are tightly related to the waves.

preprint2020arXiv

Comparative Study of Data-driven Solar Coronal Field Models Using a Flux Emergence Simulation as a Ground-truth Data Set

For a better understanding of magnetic field in the solar corona and dynamic activities such as flares and coronal mass ejections, it is crucial to measure the time-evolving coronal field and accurately estimate the magnetic energy. Recently, a new modeling technique called the data-driven coronal field model, in which the time evolution of magnetic field is driven by a sequence of photospheric magnetic and velocity field maps, has been developed and revealed the dynamics of flare-productive active regions. Here we report on the first qualitative and quantitative assessment of different data-driven models using a magnetic flux emergence simulation as a ground-truth (GT) data set. We compare the GT field with those reconstructed from the GT photospheric field by four data-driven algorithms. It is found that, at least, the flux rope structure is reproduced in all coronal field models. Quantitatively, however, the results show a certain degree of model dependence. In most cases, the magnetic energies and relative magnetic helicity are comparable to or at most twice of the GT values. The reproduced flux ropes have a sigmoidal shape (consistent with GT) of various sizes, a vertically-standing magnetic torus, or a packed structure. The observed discrepancies can be attributed to the highly non-force-free input photospheric field, from which the coronal field is reconstructed, and to the modeling constraints such as the treatment of background atmosphere, the bottom boundary setting, and the spatial resolution.

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.

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.

preprint2020arXiv

On the Lorentz Force and Torque of Solar Photospheric Emerging Magnetic Fields

Magnetic flux generated and intensified by the solar dynamo emerges into the solar atmosphere, forming active regions (ARs) including sunspots. Existing theories of flux emergence suggest that the magnetic flux can rise buoyantly through the convection zone but is trapped at the photosphere, while its further rising into the atmosphere resorts to the Parker buoyancy instability. To trigger such an instability, the Lorentz force in the photosphere needs to be as large as the gas pressure gradient to hold up an extra amount of mass against gravity. This naturally results in a strongly non-force-free photosphere, which is indeed shown in typical idealized numerical simulations of flux tube buoyancy from below the photosphere into the corona. Here we conduct a statistical study of the extents of normalized Lorentz forces and torques in the emerging photospheric magnetic field with a substantially large sample of SDO/HMI vector magnetograms. We found that the photospheric field has a rather small Lorentz force and torque on average, and thus is very close to a force-free state, which is not consistent with theories as well as idealized simulations of flux emergence. Furthermore, the small extents of forces and torques seem not to be influenced by the emerging AR's size, the emergence rate, or the non-potentiality of the field. This result puts an important constraint on future development of theories and simulations of flux emergence.

preprint2020arXiv

The Relationship between Chirality, Sense of Rotation, and Hemispheric Preference of Solar Eruptive Filaments

The orientation, chirality, and dynamics of solar eruptive filaments is a key to understanding the magnetic field of coronal mass ejections (CMEs) and therefore to predicting the geoeffectiveness of CMEs arriving at Earth. However, confusion and contention remain over the relationship between the filament chirality, magnetic helicity, and sense of rotation during eruption. To resolve the ambiguity in observations, in this paper, we used stereoscopic observations to determine the rotation direction of filament apex and the method proposed by Chen et al. (2014) to determine the filament chirality. Our sample of 12 eruptive active-region filaments establishes a strong one-to-one relationship, i.e., during the eruption, sinistral/dextral filaments (located in the southern/northern hemisphere) rotate clockwise/counterclockwise when viewed from above, and corroborates a weak hemispheric preference, i.e., a filament and related sigmoid both exhibit a forward (reverse) S shape in the southern (northern) hemisphere, which suggests that the sigmoidal filament is associated with a low-lying magnetic flux rope with its axis dipped in the middle. As a result of rotation, the projected S shape of a filament is anticipated to be reversed during eruption.

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

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.

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.

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.

preprint2014arXiv

Three-dimensional Magnetic Restructuring in Two Homologous Solar Flares in the Seismically Active NOAA AR 11283

We carry out a comprehensive investigation comparing the three-dimensional magnetic field restructuring, flare energy release, and the helioseismic response, of two homologous flares, the 2011 September 6 X2.1 (FL1) and September 7 X1.8 (FL2) flares in NOAA AR 11283. In our analysis, (1) a twisted flux rope (FR) collapses onto the surface at a speed of 1.5 km/s after a partial eruption in FL1. The FR then gradually grows to reach a higher altitude and collapses again at 3 km/s after a fuller eruption in FL2. Also, FL2 shows a larger decrease of the flux-weighted centroid separation of opposite magnetic polarities and a greater change of the horizontal field on the surface. These imply a more violent coronal implosion with corresponding more intense surface signatures in FL2. (2) The FR is inclined northward, and together with the ambient fields, it undergoes a southward turning after both events. This agrees with the asymmetric decay of the penumbra observed in the peripheral regions. (3) The amounts of free magnetic energy and nonthermal electron energy released during FL1 are comparable to those of FL2 within the uncertainties of the measurements. (4) No sunquake was detected in FL1; in contrast, FL2 produced two seismic emission sources S1 and S2 both lying in the penumbral regions. Interestingly, S1 and S2 are connected by magnetic loops, and the stronger source S2 has weaker vertical magnetic field. We discuss these results in relation to the implosion process in the low corona and the sunquake generation.

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.

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.

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.

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.

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.

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).

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.