Source author record

Yuhong Fan

Yuhong Fan 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

18works
2topics
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

18 published item(s)

preprint2026arXiv

When Magnetic Field Lines Stretch, Snap, and Expand: A New Look at Solar Flares with L-maps

Understanding the three-dimensional evolution of coronal magnetic fields during solar flares remains challenging due to the lack of direct coronal field measurements. Here we combine data-driven MHD simulations of NOAA AR 11158 (Fan et al., 2024) with flare-ribbon and coronal-dimming observations to investigate realistic coronal magnetic-field evolution during an X-class flare. We introduce L-maps - maps of natural logarithm of magnetic field-line lengths - as a diagnostic tool to track the dynamics of simulated coronal magnetic structures. Variations in L-maps identify flare ribbons through field-line shortening and coronal dimmings through field-line lengthening. Comparison with SDO/AIA observations demonstrates strong morphological and temporal agreement, validating the simulated field evolution. Applying K-means clustering to the L-map temporal profiles, we distinguish three stages of coronal evolution: (1) slow pre-flare rise phase, (2) flare reconnection accompanied by CME rise, and (3) post-reconnection CME expansion. We detect a slow pre-flare rise phase of magnetic field lines rooted in ribbon footpoints and identify reconnection dimming - area of rapid expansion of active-region core magnetic field lines during flare impulsive phase due to reconnection. Our results show that L-maps provide a powerful and physically intuitive framework for bridging simulations and observations and for tracking the full three-dimensional evolution of coronal magnetic fields during flares.

preprint2023arXiv

Solaris: A Focused Solar Polar Discovery-class Mission to achieve the Highest Priority Heliophysics Science Now

Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above.

preprint2021arXiv

SunCET: A compact EUV instrument to fill a critical observational gap

The Sun Coronal Ejection Tracker (SunCET) is an extreme ultraviolet imager and spectrograph instrument concept for tracking coronal mass ejections through the region where they experience the majority of their acceleration: the difficult-to-observe middle corona. It contains a wide field of view (0-4~\Rs) imager and a 1~Å spectral-resolution-irradiance spectrograph spanning 170-340~Å. It leverages new detector technology to read out different areas of the detector with different integration times, resulting in what we call "simultaneous high dynamic range", as opposed to the traditional high dynamic range camera technique of subsequent full-frame images that are then combined in post-processing. This allows us to image the bright solar disk with short integration time, the middle corona with a long integration time, and the spectra with their own, independent integration time. Thus, SunCET does not require the use of an opaque or filtered occulter. SunCET is also compact -- $\sim$15 $\times$ 15 $\times$ 10~cm in volume -- making it an ideal instrument for a CubeSat or a small, complementary addition to a larger mission. Indeed, SunCET is presently in a NASA-funded, competitive Phase A as a CubeSat and has also been proposed to NASA as an instrument onboard a 184 kg Mission of Opportunity.

preprint2020arXiv

Simulations of prominence eruption preceded with large amplitude longitudinal oscillations and draining

We present magnetohydrodynamic (MHD) simulations of the evolution from quasi-equilibrium to eruption of a prominence-forming twisted coronal flux rope under a coronal streamer. We have compared the cases with and without the formation of prominence condensations, and the case where prominence condensations form but we artificially initiate the draining of the prominence. We find that the prominence weight has a significant effect on the stability of the flux rope, and can significantly increase the loss-of-equilibrium height. The flux rope can be made to erupt earlier by initiating draining of the prominence mass. We have also performed a simulation where large amplitude longitudinal oscillations of the prominence are excited during the quasi-static phase. We find that the gravity force along the magnetic field lines is the major restoring force for the oscillations, in accordance with the ``pendulum model'', although the oscillation periods are higher (by about 10% to 40%) than estimated from the model because of the dynamic deformation of the field line dips during the oscillations. The oscillation period is also found to be slightly smaller for the lower part of the prominence in the deeper dips compared to the upper part in the shallower dips. The oscillations are quickly damped out after about 2-3 periods and are followed by prominence draining and the eventual eruption of the prominence. However we do not find a significant enhancement of the prominence draining and earlier onset of eruption with the excitation of the prominence oscillations compared to the case without.

preprint2016arXiv

Modeling the initiation of the 2006 December 13 coronal mass ejection in AR 10930: the structure and dynamics of the erupting flux rope

We carry out a three-dimensional magneto-hydrodynamic (MHD) simulation to model the initiation of the coronal mass ejection (CME) on 13 December 2006 in the emerging δ-sunspot active region NOAA 10930. The setup of the simulation is similar to a previous simulation by Fan (2011), but with a significantly widened simulation domain to accommodate the wide CME. The simulation shows that the CME can result from the emergence of a east-west oriented twisted flux rope whose positive, following emerging pole corresponds to the observed positive rotating sunspot emerging against the southern edge of the dominant pre-existing negative sunspot. The erupting flux rope in the simulation accelerates to a terminal speed that exceeds 1500 km/s and undergoes a counter-clockwise rotation of nearly 180 degrees such that its front and flanks all exhibit southward directed magnetic fields, explaining the observed southward magnetic field in the magnetic cloud impacting the Earth. With continued driving of flux emergence, the source region coronal magnetic field also shows the reformation of a coronal flux rope underlying the flare current sheet of the erupting flux rope, ready for a second eruption. This may explain the build up for another X-class eruptive flare that occurred the following day from the same region.

preprint2016arXiv

On the Observation and Simulation of Solar Coronal Twin Jets

We present the first observation, analysis and modeling of solar coronal twin jets, which occurred after a preceding jet. Detailed analysis on the kinetics of the preceding jet reveals its blowout-jet nature, which resembles the one studied in Liu et al. 2014. However the erupting process and kinetics of the twin jets appear to be different from the preceding one. In lack of the detailed information on the magnetic fields in the twin jet region, we instead use a numerical simulation using a three-dimensional (3D) MHD model as described in Fang et al. 2014, and find that in the simulation a pair of twin jets form due to reconnection between the ambient open fields and a highly twisted sigmoidal magnetic flux which is the outcome of the further evolution of the magnetic fields following the preceding blowout jet. Based on the similarity between the synthesized and observed emission we propose this mechanism as a possible explanation for the observed twin jets. Combining our observation and simulation, we suggest that with continuous energy transport from the subsurface convection zone into the corona, solar coronal twin jets could be generated in the same fashion addressed above.

preprint2015arXiv

$δ$-Sunspot Formation in Simulation of Active-Region-Scale Flux Emergence

$δ$-sunspots, with highly complex magnetic structures, are very productive in energetic eruptive events, such as X-class flares and homologous eruptions. We here study the formation of such complex magnetic structures by numerical simulations of magnetic flux emergence from the convection zone into the corona in an active-region-scale domain. In our simulation, two pairs of bipolar sunspots form on the surface, originating from two buoyant segments of a single subsurface twisted flux rope, following the approach of Toriumi et al. (2014). Expansion and rotation of the emerging fields in the two bipoles drive the two opposite polarities into each other with apparent rotating motion, producing a compact $δ$-sunspot with a sharp polarity inversion line. The formation of the $δ$-sunspot in such a realistic-scale domain produces emerging patterns similar to those formed in observations, e.g. the inverted polarity against Hale's law, the curvilinear motion of the spot, strong transverse field with highly sheared magnetic and velocity fields at the PIL. Strong current builds up at the PIL, giving rise to reconnection, which produces a complex coronal magnetic connectivity with non-potential fields in the Delta-spot overlaid by more relaxed fields connecting the two polarities at the two ends.

preprint2015arXiv

Differential Rotation in Solar Convective Dynamo Simulations

We carry out a magneto-hydrodynamic (MHD) simulation of convective dynamo in the rotating solar convective envelope driven by the solar radiative diffusive heat flux. The simulation is similar to that reported in Fan & Fang (2014) but with further reduced viscosity and magnetic diffusion. The resulting convective dynamo produces a large scale mean field that exhibits similar irregular cyclic behavior and polarity reversals, and self-consistently maintains a solar-like differential rotation. The main driver for the solar-like differential rotation (with faster rotating equator) is a net outward transport of angular momentum away from the rotation axis by the Reynolds stress, and we found that this transport is enhanced with reduced viscosity and magnetic diffusion.

preprint2015arXiv

Effects of Radiative Diffusion on Thin Flux Tubes in Turbulent Solar-like Convection

We study the combined effects of convection and radiative diffusion on the evolution of thin magnetic flux tubes in the solar interior. Radiative diffusion is the primary supplier of heat to convective motions in the lower convection zone, and it results in a heat input per unit volume of magnetic flux tubes that has been ignored by many previous thin flux tube studies. We use a thin flux tube model subject to convection taken from a rotating spherical shell of turbulent, solar-like convection as described by Weber, Fan, and Miesch (2011, Astrophys. J., 741, 11; 2013, Solar Phys., 287, 239), now taking into account the influence of radiative heating on flux tubes of large-scale active regions. Our simulations show that flux tubes of less than or equal to 60 kG subject to solar-like convective flows do not anchor in the overshoot region, but rather drift upward due to the increased buoyancy of the flux tube earlier in its evolution as a result of the inclusion of radiative diffusion. Flux tubes of magnetic field strengths ranging from 15 kG to 100 kG have rise times of less than or equal to 0.2 years, and exhibit a Joy's Law tilt-angle trend. Our results suggest that radiative heating is an effective mechanism by which flux tubes can escape from the stably stratified overshoot region, and that flux tubes do not necessarily need to be anchored in the overshoot region to produce emergence properties similar to those of active regions on the Sun.

preprint2015arXiv

MHD Equilibria and Triggers for Prominence Eruption

Magneto-hydrodynamic (MHD) simulations of the emergence of twisted magnetic flux tubes from the solar interior into the corona are discussed to illustrate how twisted and sheared coronal magnetic structures (with free magnetic energy), capable of driving filament eruptions, can form in the corona in emerging active regions. Several basic mechanisms that can disrupt the quasi-equilibrium coronal structures and trigger the release of the stored free magnetic energy are discussed. These include both ideal processes such as the onset of the helical kink instability and the torus instability of a twisted coronal flux rope structure and the non-ideal process of the onset of fast magnetic reconnections in current sheets. Representative MHD simulations of the non-linear evolution involving these mechanisms are presented.

preprint2015arXiv

Numerical Study on Emergence of Kinked Flux Tube for Understanding of Possible Origin of Delta-spot Regions

We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence, the magnetic twist at the kinked apex of the tube is greatly reduced, although the other parts of the tube is still strongly twisted. This leads to the formation of a complex quadrupole structure: a pair of the coherent, strongly twisted spots and a narrow complex bipolar pair between it. The quadrupole is formed by the submergence of a portion of emerged magnetic fields. This result is relevant for understanding of the origin of the complex multipolar $δ$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws.

preprint2014arXiv

A simulation of convective dynamo in the solar convective envelope: maintenance of the solar-like differential rotation and emerging flux

We report the results of a magneto-hydrodynamic (MHD) simulation of a convective dynamo in a model solar convective envelope driven by the solar radiative diffusive heat flux. The convective dynamo produces a large-scale mean magnetic field that exhibits irregular cyclic behavior with oscillation time scales ranging from about 5 to 15 years and undergoes irregular polarity reversals. The mean axisymmetric toroidal magnetic field is of opposite signs in the two hemispheres and is concentrated at the bottom of the convection zone. The presence of the magnetic fields is found to play an important role in the self-consistent maintenance of a solar-like differential rotation in the convective dynamo model. Without the magnetic fields, the convective flows drive a differential rotation with a faster rotating polar region. In the midst of magneto-convection, we found emergence of strong super-equipartition flux bundles at the surface, exhibiting properties that are similar to emerging solar active regions.

preprint2014arXiv

Rotating Solar Jets in Simulations of Flux Emergence with Thermal Conduction

We study the formation of coronal jets through numerical simulation of the emergence of a twisted magnetic flux rope into a pre-existing open magnetic field. Reconnection inside the emerging flux rope in addition to that between the emerging and pre-existing fields give rise to the violent eruption studied. The simulated event closely resembles the coronal jets ubiquitously observed by Hinode/XRT and demonstrates that heated plasma is driven into the extended atmosphere above. Thermal conduction implemented in the model allows us to qualitatively compare simulated and observed emission from such events. We find that untwisting field lines after the reconnection drive spinning outflows of plasma in the jet column. The Poynting flux in the simulated jet is dominated by the untwisting motions of the magnetic fields loaded with high-density plasma. The simulated jet is comprised of spires of untwisting field that are loaded with a mixture of cold and hot plasma and exhibit rotational motion of order 20 km/s and match contemporary observations.

preprint2014arXiv

Simulation of homologous and cannibalistic Coronal Mass Ejections produced by the emergence of a twisted flux rope into the Solar Corona

We report the first results of a magnetohydrodynamic (MHD) simulation of the development of a homologous sequence of three coronal mass ejections (CMEs) and demonstrate their so-called cannibalistic behavior. These CMEs originate from the repeated formations and partial eruptions of kink unstable flux ropes as a result of continued emergence of a twisted flux rope across the lower boundary into a pre-existing coronal potential arcade field. The simulation shows that a CME erupting into the open magnetic field created by a preceding CME has a higher speed. The second of the three successive CMEs simulated is cannibalistic catching up and merging with the first into a single fast CME before exiting the domain. All the CMEs including the leading merged CME, attained speeds of about 1000 km/s as they exit the domain. The reformation of a twisted flux rope after each CME eruption during the sustained flux emergence can naturally explain the X-ray observations of repeated reformations of sigmoids and "sigmoid-under-cusp" configurations at a low-coronal source of homologous CMEs.

preprint2013arXiv

A Theory on the Convective Origins of Active Longitudes on Solar-like Stars

Using a thin flux tube model in a rotating spherical shell of turbulent, solar-like convective flows, we find that the distribution of emerging flux tubes in our simulation is inhomogeneous in longitude, with properties similar to those of active longitudes on the Sun and other solar-like stars. The large-scale pattern of flux emergence our simulations produce exhibits preferred longitudinal modes of low order, drift with respect to a fixed reference system, and alignment across the Equator at low latitudes between 15 degrees. We suggest that these active-longitude-like emergence patterns are the result of columnar, rotationally aligned giant cells present in our convection simulation at low latitudes. If giant convecting cells exist in the bulk of the solar convection zone, this phenomenon, along with differential rotation, could in part provide an explanation for the behavior of active longitudes.

preprint2013arXiv

Three-Dimensional MHD Simulations of Emerging Active Region Flux in a Turbulent Rotating Solar Convective Envelope: the Numerical Model and Initial Results

We describe a 3D finite-difference spherical anelastic MHD (FSAM) code for modeling the subsonic dynamic processes in the solar convective envelope. A comparison of this code with the widely used global spectral anlastic MHD code, ASH (Anelastic Spherical Harmonics), shows that FSAM produces convective flows with statistical properties and mean flows similar to the ASH results. Using FSAM, we first simulate the rotating solar convection in a partial spherical shell domain and obtain a statistically steady, giant-cell convective flow with a solar-like differential rotation. We then insert buoyant toroidal flux tubes near the bottom of the convecting envelope and simulate the rise of the flux tubes in the presence of the giant cell convection. We find that for buoyant flux tubes with an initial field strength of 100 kG, the magnetic buoyancy largely determines the rise of the tubes although strong down flows produce significant undulation and distortion to the shape of the emerging $Ω$-shaped loops. The convective flows significantly reduce the rise time it takes for the apex of the flux tube to reach the top. For the weakly twisted and untwisted cases we simulated, the apex portion is found to rise nearly radially to the top in about a month, and produce an emerging region (at a depth of about 30 Mm below the photosphere) with an overall tilt angle consistent with the active region tilts, although the emergence pattern is more complex compared to the case without convection. Near the top boundary at a depth of about 30 Mm, the emerging flux shows a retrograde zonal flow of about 345 m/s relative to the mean flow at that latitude.

preprint2012arXiv

Comparing Simulations of Rising Flux Tubes Through the Solar Convection Zone with Observations of Solar Active Regions: Constraining the Dynamo Field Strength

We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1 degree to 40 degrees in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J., 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 10^20 Mx and 10^21 Mx, and more simulations of the previously investigated case of 10^22 Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy's Law as in Article 1, and in addition the scatter of tilt angles about the Joy's Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of \geq 40 kG are good candidates for the progenitors of large (10^21 Mx to 10^22 Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10^22 Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy's Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.

preprint2011arXiv

The Rise of Active Region Flux Tubes in the Turbulent Solar Convective Envelope

We use a thin flux tube model in a rotating spherical shell of turbulent convective flows to study how active region scale flux tubes rise buoyantly from the bottom of the convection zone to near the solar surface. We investigate toroidal flux tubes at the base of the convection zone with field strengths ranging from 15 kG to 100 kG at initial latitudes ranging from 1 degree to 40 degrees with a total flux of 10^22 Mx. We find that the dynamic evolution of the flux tube changes from convection dominated to magnetic buoyancy dominated as the initial field strength increases from 15 kG to 100 kG. At 100 kG, the development of Omega-shaped rising loops is mainly controlled by the growth of the magnetic buoyancy instability. However, at low field strengths of 15 kG, the development of rising Omega-shaped loops is largely controlled by convective flows, and properties of the emerging loops are significantly changed compared to previous results in the absence of convection. With convection, rise times are drastically reduced (from years to a few months), loops are able to emerge at low latitudes, and tilt angles of emerging loops are consistent with Joy's Law for initial field strengths of greater than or equal to 40 kG. We also examine other asymmetries that develop between the leading and following legs of the emerging loops. Taking all the results together, we find that mid-range field strengths of approximately 40 - 50 kG produce emerging loops that best match the observed properties of solar active regions.