-
Adjusting the Potential Field Source Surface Height Based on MHD Simulations
Authors:
Zhenguang Huang,
Gabor Toth,
Jia Huang,
Nishtha Sachdeva,
Bart van der Holst,
Ward B. Manchester
Abstract:
A potential field solution is widely used to extrapolate the coronal magnetic field above the Sun's surface to a certain height. This model applies the current-free approximation and assumes that the magnetic field is entirely radial beyond the source surface height, which is defined as the radial distance from the center of the Sun. Even though the source surface is commonly specified at 2.5 Rs (…
▽ More
A potential field solution is widely used to extrapolate the coronal magnetic field above the Sun's surface to a certain height. This model applies the current-free approximation and assumes that the magnetic field is entirely radial beyond the source surface height, which is defined as the radial distance from the center of the Sun. Even though the source surface is commonly specified at 2.5 Rs (solar radii), previous studies have suggested that this value is not optimal in all cases. In this study, we propose a novel approach to specify the source surface height, by comparing the areas of the open magnetic field regions from the potential field solution with predictions made by a magnetohydrodynamics model, in our case the Alfven Wave Solar atmosphere Model. We find that the adjusted source surface height is significantly less than 2.5 Rs near solar minimum, and slightly larger than 2.5 Rs near solar maximum. We also report that the adjusted source surface height can provide a better open flux agreement with the observations near the solar minimum, while the comparison near the solar maximum is slightly worse.
△ Less
Submitted 3 March, 2024;
originally announced March 2024.
-
Solar Wind Driven from GONG Magnetograms in the Last Solar Cycle
Authors:
Zhenguang Huang,
Gabor Toth,
Nishtha Sachdeva,
Bart van der Holst
Abstract:
In a previous study, Huang et al. (2023) used the Alfven Wave Solar atmosphere Model (AWSoM), one of the widely used solar wind models in the community, driven by ADAPT-GONG magnetograms to simulate the solar wind in the last solar cycle and found that the optimal Poynting flux parameter can be estimated from either the open field area or the average unsigned radial component of the magnetic field…
▽ More
In a previous study, Huang et al. (2023) used the Alfven Wave Solar atmosphere Model (AWSoM), one of the widely used solar wind models in the community, driven by ADAPT-GONG magnetograms to simulate the solar wind in the last solar cycle and found that the optimal Poynting flux parameter can be estimated from either the open field area or the average unsigned radial component of the magnetic field in the open field regions. It was also found that the average energy deposition rate (Poynting flux) in the open field regions is approximately constant. In the current study, we expand the previous work by using GONG magnetograms to simulate the solar wind for the same Carrington rotations and determine if the results are similar to the ones obtained with ADAPT-GONG magnetograms. Our results indicate that similar correlations can be obtained from the GONG maps. Moreover, we report that ADAPT-GONG magnetograms can consistently provide better comparisons with 1 AU solar wind observations than GONG magnetograms, based on the best simulations selected by the minimum of the average curve distance for the solar wind speed and density.
△ Less
Submitted 3 March, 2024;
originally announced March 2024.
-
Analytic Model and Magnetohydrodynamic Simulations of Three-dimensional Magnetic Switchbacks
Authors:
Chen Shi,
Marco Velli,
Gabor Toth,
Kun Zhang,
Anna Tenerani,
Zesen Huang,
Nikos Sioulas,
Bart van der Holst
Abstract:
Parker Solar Probe observations reveal that the near-Sun space is almost filled with magnetic switchbacks (``switchbacks'' hereinafter), which may be a major contributor to the heating and acceleration of solar wind. Here, for the first time, we develop an analytic model of an axisymmetric switchback with uniform magnetic field strength. In this model, three parameters control the geometry of the…
▽ More
Parker Solar Probe observations reveal that the near-Sun space is almost filled with magnetic switchbacks (``switchbacks'' hereinafter), which may be a major contributor to the heating and acceleration of solar wind. Here, for the first time, we develop an analytic model of an axisymmetric switchback with uniform magnetic field strength. In this model, three parameters control the geometry of the switchback: height (length along the background magnetic field), width (thickness along radial direction perpendicular to the background field), and the radial distance from the center of switchback to the central axis, which is a proxy of the size of the switchback along the third dimension. We carry out three-dimensional magnetohydrodynamic simulations to investigate the dynamic evolution of the switchback. Comparing simulations conducted with compressible and incompressible codes, we verify that compressibility, i.e. parametric decay instability, is necessary for destabilizing the switchback. Our simulations also reveal that the geometry of the switchback significantly affects how fast the switchback destabilizes. The most stable switchbacks are 2D-like (planar) structures with large aspect ratios (length to width), consistent with the observations. We show that when plasma beta ($β$) is smaller than one, the switchback is more stable as $β$ increases. However, when $β$ is greater than one, the switchback becomes very unstable as the pattern of the growing compressive fluctuations changes. Our results may explain some of the observational features of switchbacks, including the large aspect ratios and nearly constant occurrence rates in the inner heliosphere.
△ Less
Submitted 31 March, 2024; v1 submitted 20 January, 2024;
originally announced January 2024.
-
Non-equilibrium Ionization Effects on Synthetic Spectra in the AWSoM Solar Corona
Authors:
Judit Szente,
Enrico Landi,
Bart van der Holst
Abstract:
In this work we combined AWSoM's non-equilibrium ionization [NEI] calculations from (Szente:2022) with the synthetic spectral computations of SPECTRUM (Szente:2019), to predict non-equilibrium line intensities across the entire domain of the AWSoM 3D global model. We find that the resulting spectra are strongly affected by non-equilibrium effects in the fast wind regions and streamer edges and tha…
▽ More
In this work we combined AWSoM's non-equilibrium ionization [NEI] calculations from (Szente:2022) with the synthetic spectral computations of SPECTRUM (Szente:2019), to predict non-equilibrium line intensities across the entire domain of the AWSoM 3D global model. We find that the resulting spectra are strongly affected by non-equilibrium effects in the fast wind regions and streamer edges and that these effects propagate to narrowband images from SoHO/EIT, SECCHI/EUVI and SDO/AIA. The dependence shows a different nature for each line observed resulting in significant changes in line intensity, which need to be accounted for during plasma diagnostics. However, we also find that these effects depend on the local plasma properties, and that no single correction can be developed to account for non-equilibrium effects in observed spectra and images. Comparing to observational data we saw that the changes due to NEI, while significant, are not sufficient to account for the differences between Hinode/EIS spectra and AWSoM/SPECTRUM predictions.
△ Less
Submitted 21 October, 2023;
originally announced October 2023.
-
Solar Wind with Field Lines and Energetic Particles (SOFIE) Model: Application to Historical Solar Energetic Particle Events
Authors:
Lulu Zhao,
Igor Sokolov,
Tamas Gombosi,
David Lario,
Kathryn Whitman,
Zhenguang Huang,
Gabor Toth,
Ward Manchester,
Bart van der Holst,
Nishtha Sachdeva
Abstract:
In this paper, we demonstrate the applicability of the data-driven and self-consistent solar energetic particle model, Solar-wind with FIeld-lines and Energetic-particles (SOFIE), to simulate acceleration and transport processes of solar energetic particles. SOFIE model is built upon the Space Weather Modeling Framework (SWMF) developed at the University of Michigan. In SOFIE, the background solar…
▽ More
In this paper, we demonstrate the applicability of the data-driven and self-consistent solar energetic particle model, Solar-wind with FIeld-lines and Energetic-particles (SOFIE), to simulate acceleration and transport processes of solar energetic particles. SOFIE model is built upon the Space Weather Modeling Framework (SWMF) developed at the University of Michigan. In SOFIE, the background solar wind plasma in the solar corona and interplanetary space is calculated by the Aflvén Wave Solar-atmosphere Model(-Realtime) (AWSoM-R) driven by the near-real-time hourly updated Global Oscillation Network Group (GONG) solar magnetograms. In the background solar wind, coronal mass ejections (CMEs) are launched by placing an imbalanced magnetic flux rope on top of the parent active region, using the Eruptive Event Generator using Gibson-Low model (EEGGL). The acceleration and transport processes are modeled by the Multiple-Field-Line Advection Model for Particle Acceleration (M-FLAMPA). In this work, nine solar energetic particle events (Solar Heliospheric and INterplanetary Environment (SHINE) challenge/campaign events) are modeled. The three modules in SOFIE are validated and evaluated by comparing with observations, including the steady-state background solar wind properties, the white-light image of the CME, and the flux of solar energetic protons, at energies of > 10 MeV.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Formation of Magnetic Switchbacks Observed by Parker Solar Probe
Authors:
Gabor Toth,
Marco Velli,
Bart van der Holst
Abstract:
Magnetic switchbacks are rapid high amplitude reversals of the radial magnetic field in the solar wind that do not involve a heliospheric current sheet crossing. First seen sporadically in the seventies in Mariner and Helios data, switchbacks were later observed by the Ulysses spacecraft beyond 1 au and have been recently identified as a typical component of solar wind fluctuations in the inner he…
▽ More
Magnetic switchbacks are rapid high amplitude reversals of the radial magnetic field in the solar wind that do not involve a heliospheric current sheet crossing. First seen sporadically in the seventies in Mariner and Helios data, switchbacks were later observed by the Ulysses spacecraft beyond 1 au and have been recently identified as a typical component of solar wind fluctuations in the inner heliosphere by the Parker Solar Probe spacecraft. Here we provide a simple yet predictive theory for the formation of these magnetic reversals: the switchbacks are produced by the shear of circularly polarized Alfvén waves by a transversely varying radial wave propagation velocity. We provide an analytic expression for the magnetic field variation, establish the necessary and sufficient conditions and show that the mechanism works in a realistic solar wind scenario.
△ Less
Submitted 6 January, 2023;
originally announced January 2023.
-
Solar wind modeling with the Alfven Wave Solar atmosphere Model driven by HMI-based Near-Real-Time maps by the National Solar Observatory
Authors:
Nishtha Sachdeva,
Ward B Manchester,
Igor Sokolov,
Zhenguang Huang,
Alexander Pevtsov,
Luca Bertello,
Alexei A. Pevtsov,
Gabor Toth,
Bart van der Holst
Abstract:
We explore model performance for the Alfven Wave Solar atmosphere Model (AWSoM) with near-real-time (NRT) synoptic maps of the photospheric vector magnetic field. These maps, produced by assimilating data from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), use a different method developed at the National Solar Observatory (NSO) to provide a near contemporaneou…
▽ More
We explore model performance for the Alfven Wave Solar atmosphere Model (AWSoM) with near-real-time (NRT) synoptic maps of the photospheric vector magnetic field. These maps, produced by assimilating data from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), use a different method developed at the National Solar Observatory (NSO) to provide a near contemporaneous source of data to drive numerical models. Here, we apply these NSO-HMI-NRT maps to simulate three Carrington rotations (CRs): 2107-2108 (centered on 2011/03/07 20:12 CME event), 2123 (integer CR) and 2218--2219 (centered on 2019/07/2 solar eclipse), which together cover a wide range of activity level for solar cycle 24. We show simulation results, which reproduce both extreme ultraviolet emission (EUV) from the low corona while simultaneously matching in situ observations at 1 au as well as quantify the total unsigned open magnetic flux from these maps.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Modeling the Solar Wind During Different Phases of the Last Solar Cycle
Authors:
Zhenguang Huang,
Gabor Toth,
Nishtha Sachdeva,
Lulu Zhao,
Bartholomeus van der Holst,
Igor Sokolov,
Ward Manchester,
Tamas Gombosi
Abstract:
We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfvén Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameters of the model, the Poynting flux parameter, which prescribes the energy flux pass…
▽ More
We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfvén Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameters of the model, the Poynting flux parameter, which prescribes the energy flux passing through the chromospheric boundary of the model in the form of Alfvén wave turbulence. It is found that the optimal value of the Poynting flux parameter is correlated with the area of the open magnetic field regions with the Spearman's correlation coefficient of 0.96 and anti-correlated with the average unsigned radial component of the magnetic field with the Spearman's correlation coefficient of -0.91. Moreover, the Poynting flux in the open field regions is approximately constant in the last solar cycle, which needs to be validated with observations and can shed light on how Alfvén wave turbulence accelerates the solar wind during different phases of the solar cycle. Our results can also be used to set the Poynting flux parameter for real-time solar wind simulations with AWSoM.
△ Less
Submitted 24 March, 2023; v1 submitted 5 October, 2022;
originally announced October 2022.
-
One Year in the Life of Young Suns: Data Constrained Corona-Wind Model of kappa1 Ceti
Authors:
Vladimir S. Airapetian,
Meng Jin,
Theresa Lueftinger,
Sudesha Boro Saikia,
Oleg Kochukhov,
Manuel Guedel,
Bart Van Der Holst,
W. Manchester IV
Abstract:
The young magnetically active solar-like stars are efficient generators of ionizing radiation in the form of X-ray and Extreme UV (EUV) flux, stellar wind and eruptive events. These outputs are the critical factors affecting atmospheric escape and chemistry of (exo)planets around active stars. While X-ray fluxes and surface magnetic fields can be derived from observations, the EUV emission and win…
▽ More
The young magnetically active solar-like stars are efficient generators of ionizing radiation in the form of X-ray and Extreme UV (EUV) flux, stellar wind and eruptive events. These outputs are the critical factors affecting atmospheric escape and chemistry of (exo)planets around active stars. While X-ray fluxes and surface magnetic fields can be derived from observations, the EUV emission and wind mass fluxes, Coronal Mass Ejections and associated Stellar Energetic Particle events cannot be directly observed. Here, we present the results of a three-dimensional magnetohydrodynamic (MHD) model with inputs constrained by spectropolarimetric data, HST/STIS Far UV, X-ray data data and stellar magnetic maps reconstructed at two epochs separated by 11 months. The simulations show that over the course of the year, the global stellar corona had undergone a drastic transition from a simple dipole-like to a tilted dipole with multipole field components, and thus, provided favorable conditions for Corotating Interaction Events (CIRs) that drive strong shocks. The dynamic pressure exerted by CIRs are 1300 times larger than those observed from the Sun and can contribute to the atmospheric erosion of early Venus, Earth, Mars and young Earth-like exoplanets. Our data-constrained MHD model provides the framework to model coronal environments of G-M planet hosting dwarfs. The model outputs can serve as a realistic input for exoplanetary atmospheric models to evaluate the impact of stellar coronal emission, stellar winds and CIRs on their atmospheric escape and chemistry that can be tested in the upcoming JWST and ground-based observations.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
What Sustained Multi-Disciplinary Research Can Achieve: The Space Weather Modeling Framework
Authors:
Tamas I. Gombosi,
Yuxi Chen,
Alex Glocer,
Zhenguang Huang,
Xianzhe Jia,
Michael W. Liemohn,
Ward B. Manchester,
Tuija Pulkkinen,
Nishtha Sachdeva,
Qusai Al Shidi,
Igor V. Sokolov,
Judit Szente,
Valeriy Tenishev,
Gabor Toth,
Bart van der Holst,
Daniel T. Welling,
Lulu Zhao,
Shasha Zou
Abstract:
MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environ…
▽ More
MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC).
△ Less
Submitted 27 May, 2021;
originally announced May 2021.
-
Thermodynamic Structure of the Solar Corona: Tomographic Reconstructions and MHD Modeling
Authors:
D. G. Lloveras,
A. M. Vásquez,
F. A. Nuevo,
C. Mac Cormack,
N. Sachdeva,
W. Manchester IV,
B. Van der Holst,
R. A. Frazin
Abstract:
We carry out a study of the global three-dimensional (3D) structure of the electron density and temperature of the quiescent inner solar corona ($r<1.25 R_\odot$) by means of tomographic reconstructions and magnetohydrodynamic simulations. We use differential emission measure tomography (DEMT) and the Alfvén Wave Solar Model (AWSoM), in their latest versions. Two target rotations were selected fro…
▽ More
We carry out a study of the global three-dimensional (3D) structure of the electron density and temperature of the quiescent inner solar corona ($r<1.25 R_\odot$) by means of tomographic reconstructions and magnetohydrodynamic simulations. We use differential emission measure tomography (DEMT) and the Alfvén Wave Solar Model (AWSoM), in their latest versions. Two target rotations were selected from the solar minimum between solar cycles (SCs) 23 and 24 and the declining phase of SC 24. We report in quantitative detail on the 3D thermodynamic structure of the core and outer layers of the streamer belt, and of the high latitude coronal holes (CH), as revealed by the DEMT analysis. We report on the presence of two types of structures within the streamer belt, loops with temperature decreasing/increasing with height (dubbed down/up loops), as reported first in previous DEMT studies. We also estimate the heating energy flux required at the coronal base to keep these structures stable, found to be or order $10^5 erg\, cm^{-2} s^{-1}$, consistently with previous DEMT and spectroscopic studies. We discuss how these findings are consistent with coronal dissipation of Alfvén waves. We compare the 3D results of DEMT and AWSoM in distinct magnetic structures. We show that the agreement between the products of both techniques is the best so far, with an overall agreement $\lesssim 20\%$, depending on the target rotation and the specific coronal region. In its current implementation the ASWsoM model can not reproduce down loops though. Also, in the source region of the fast and slow components of the solar wind, the electron density of the AWSoM model increases with latitude, opposite to the trend observed in DEMT reconstructions.
△ Less
Submitted 14 April, 2020;
originally announced April 2020.
-
Validation of the Alfvén Wave Solar atmosphere Model (AWSoM) with Observations from the Low Corona to 1 AU
Authors:
N. Sachdeva,
B. van der Holst,
W. B. Manchester,
G. Tóth,
Y. Chen,
D. G. Lloveras,
A. M. Vásquez,
Philippe Lamy,
Julien Wojak,
B. V. Jackson,
H. -S. Yu,
C. J. Henney
Abstract:
We perform a validation study of the latest version of the Alfvén Wave Solar atmosphere Model (AWSoM) within the Space Weather Modeling Framework (SWMF). To do so, we compare the simulation results of the model with a comprehensive suite of observations for Carrington rotations representative of the solar minimum conditions extending from the solar corona to the heliosphere up to the Earth. In the…
▽ More
We perform a validation study of the latest version of the Alfvén Wave Solar atmosphere Model (AWSoM) within the Space Weather Modeling Framework (SWMF). To do so, we compare the simulation results of the model with a comprehensive suite of observations for Carrington rotations representative of the solar minimum conditions extending from the solar corona to the heliosphere up to the Earth. In the low corona ($r < 1.25$ \Rs), we compare with EUV images from both STEREO-A/EUVI and SDO/AIA and to three-dimensional (3-D) tomographic reconstructions of the electron temperature and density based on these same data. We also compare the model to tomographic reconstructions of the electron density from SOHO/LASCO observations ($2.55 < r < 6.0$\Rs). In the heliosphere, we compare model predictions of solar wind speed with velocity reconstructions from InterPlanetary Scintillation (IPS) observations. For comparison with observations near the Earth, we use OMNI data. Our results show that the improved AWSoM model performs well in quantitative agreement with the observations between the inner corona and 1 AU. The model now reproduces the fast solar wind speed in the polar regions. Near the Earth, our model shows good agreement with observations of solar wind velocity, proton temperature and density. AWSoM offers an extensive application to study the solar corona and larger heliosphere in concert with current and future solar missions as well as being well suited for space weather predictions.
△ Less
Submitted 17 October, 2019;
originally announced October 2019.
-
Predictions for the First Parker Solar Probe Encounter
Authors:
B. van der Holst,
W. B. Manchester,
K. G. Klein,
J. C. Kasper
Abstract:
We examine Alfvén Wave Solar atmosphere Model (AWSoM) predictions of the first Parker Solar Probe (PSP) encounter. We focus on the 12-day closest approach centered on the 1st perihelion. AWSoM (van der Holst et al., 2014) allows us to interpret the PSP data in the context of coronal heating via Alfvén wave turbulence. The coronal heating and acceleration is addressed via outward-propagating low-fr…
▽ More
We examine Alfvén Wave Solar atmosphere Model (AWSoM) predictions of the first Parker Solar Probe (PSP) encounter. We focus on the 12-day closest approach centered on the 1st perihelion. AWSoM (van der Holst et al., 2014) allows us to interpret the PSP data in the context of coronal heating via Alfvén wave turbulence. The coronal heating and acceleration is addressed via outward-propagating low-frequency Alfvén waves that are partially reflected by Alfvén speed gradients. The nonlinear interaction of these counter-propagating waves results in a turbulent energy cascade. To apportion the wave dissipation to the electron and anisotropic proton temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating as described by Chandran et al. (2011). We find that during the first encounter, PSP was in close proximity to the heliospheric current sheet (HCS) and in the slow wind. PSP crossed the HCS two times, namely at 2018/11/03 UT 01:02 and 2018/11/08 UT 19:09 with perihelion occuring on the south of side of the HCS. We predict the plasma state along the PSP trajectory, which shows a dominant proton parallel temperature causing the plasma to be firehose unstable.
△ Less
Submitted 8 February, 2019;
originally announced February 2019.
-
Extended MHD modeling of the steady solar corona and the solar wind
Authors:
Tamas I. Gombosi,
Bart van der Holst,
Ward B. Manchester,
Igor V. Sokolov
Abstract:
The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind mo…
▽ More
The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind models are described in detail. This historical perspective is also applied to the present suite of global corona and solar wind models. We discuss the evolution of new ideas and their implementation into numerical simulation codes. We point out the scientific and computational challenges facing these models and discuss the ways various groups tried to overcome these challenges. Next, we discuss the latest, state-of-the art models and point to the expected next steps in modeling the corona and the interplanetary medium.
△ Less
Submitted 1 July, 2018;
originally announced July 2018.
-
Laboratory Photoionization Fronts in Nitrogen Gas: A Numerical Feasibility and Parameter Study
Authors:
William J Gray,
P. A. Keiter,
H. Lefevre,
C. R. Patterson,
J. S. Davis,
B. van Der Holst,
K. G. Powell,
R. P. Drake
Abstract:
Photoionization fronts play a dominant role in many astrophysical situations, but remain difficult to achieve in a laboratory experiment. We present the results from a computational parameter study evaluating the feasibility of the photoionization experiment presented in the design paper by Drake, R. P., Hazak, G., Keiter, P. A., Davis, J. S., Patterson, C. R., Frank, A., Blackman, E. G., & Busque…
▽ More
Photoionization fronts play a dominant role in many astrophysical situations, but remain difficult to achieve in a laboratory experiment. We present the results from a computational parameter study evaluating the feasibility of the photoionization experiment presented in the design paper by Drake, R. P., Hazak, G., Keiter, P. A., Davis, J. S., Patterson, C. R., Frank, A., Blackman, E. G., & Busquet, M. 2016, ApJ, 833, 249 in which a photoionization front is generated in a nitrogen medium . The nitrogen gas density and the Planckian radiation temperature of the x-ray source define each simulation. Simulations modeled experiments in which the x-ray flux is generated by a laser-heated gold foil, suitable for experiments using many kJ of laser energy, and experiments in which the flux is generated by a "z-pinch" device, which implodes a cylindrical shell of conducting wires. The models are run using CRASH, our block-adaptive-mesh code for multi-material radiation hydrodynamics. The radiative transfer model uses multi-group, flux-limited diffusion with thirty radiation groups. In addition, electron heat conduction is modeled using a single-group, flux-limited diffusion. In the theory, a photoionization front can exist only when the ratios of the electron recombination rate to the photoionization rate and the electron impact ionization rate to the recombination rate lie in certain ranges. These ratios are computed for several ionization states of nitrogen. Photoionization fronts are found to exist for laser driven models with moderate nitrogen densities ($\sim$10$^{21}$ cm$^{-3}$) and radiation temperatures above 90 eV. For "z-pinch" driven models, lower nitrogen densities are preferred ($<$10$^{21}$ cm$^{-3}$). We conclude that the proposed experiments are likely to generate photoionization fronts.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
Life Beyond the Solar System: Space Weather and Its Impact on Habitable Worlds
Authors:
V. S. Airapetian,
W. C. Danchi,
C. F. Dong,
S. Rugheimer,
M. Mlynczak,
K. B. Stevenson,
W. G. Henning,
J. L. Grenfell,
M. Jin,
A. Glocer,
G. Gronoff,
B. Lynch,
C. Johnstone,
T. Lueftinger,
M. Guedel,
K. Kobayashi,
A. Fahrenbach,
G. Hallinan,
V. Stamenkovic,
O. Cohen,
W. Kuang,
B. van der Holst,
C. Manchester,
G. Zank,
O. Verkhoglyadova
, et al. (11 additional authors not shown)
Abstract:
The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the sear…
▽ More
The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe.
△ Less
Submitted 16 January, 2018;
originally announced January 2018.
-
Atmospheric escape from the TRAPPIST-1 planets and implications for habitability
Authors:
Chuanfei Dong,
Meng Jin,
Manasvi Lingam,
Vladimir S. Airapetian,
Yingjuan Ma,
Bart van der Holst
Abstract:
The presence of an atmosphere over sufficiently long timescales is widely perceived as one of the most prominent criteria associated with planetary surface habitability. We address the crucial question as to whether the seven Earth-sized planets transiting the recently discovered ultracool dwarf star TRAPPIST-1 are capable of retaining their atmospheres. To this effect, we carry out numerical simu…
▽ More
The presence of an atmosphere over sufficiently long timescales is widely perceived as one of the most prominent criteria associated with planetary surface habitability. We address the crucial question as to whether the seven Earth-sized planets transiting the recently discovered ultracool dwarf star TRAPPIST-1 are capable of retaining their atmospheres. To this effect, we carry out numerical simulations to characterize the stellar wind of TRAPPIST-1 and the atmospheric ion escape rates for all the seven planets. We also estimate the escape rates analytically and demonstrate that they are in good agreement with the numerical results. We conclude that the outer planets of the TRAPPIST-1 system are capable of retaining their atmospheres over billion-year timescales. The consequences arising from our results are also explored in the context of abiogenesis, biodiversity, and searches for future exoplanets. In light of the many unknowns and assumptions involved, we recommend that these conclusions must be interpreted with due caution.
△ Less
Submitted 8 December, 2017; v1 submitted 16 May, 2017;
originally announced May 2017.
-
Chromosphere to 1 AU Simulation of the 2011 March 7th Event: A Comprehensive Study of Coronal Mass Ejection Propagation
Authors:
M. Jin,
W. B. Manchester,
B. van der Holst,
I. Sokolov,
G. Toth,
A. Vourlidas,
C. A. de Koning,
T. I. Gombosi
Abstract:
We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfvén Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with th…
▽ More
We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfvén Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with the Space Weather Modeling Framework (SWMF) extends the total domain beyond the orbit of Earth. Physical processes included in the model are multi-species thermodynamics, electron heat conduction (both collisional and collisionless formulations), optically thin radiative cooling, and Alfvén-wave turbulence that accelerates and heats the solar wind. The Alfvén-wave description is physically self-consistent, including non-Wentzel-Kramers-Brillouin (WKB) reflection and physics-based apportioning of turbulent dissipative heating to both electrons and protons. Within this model, we initiate the CME by using the Gibson-Low (GL) analytical flux rope model and follow its evolution for days, in which time it propagates beyond STEREO A. A detailed comparison study is performed using remote as well as \textit{in situ} observations. Although the flux rope structure is not compared directly due to lack of relevant ejecta observation at 1 AU in this event, our results show that the new model can reproduce many of the observed features near the Sun (e.g., CME-driven extreme ultraviolet (EUV) waves, deflection of the flux rope from the coronal hole, "double-front" in the white light images) and in the heliosphere (e.g., shock propagation direction, shock properties at STEREO A).
△ Less
Submitted 27 November, 2016;
originally announced November 2016.
-
Threaded-Field-Lines Model for the Low Solar Corona Powered by the Alfven Wave Turbulence
Authors:
Igor V. Sokolov,
Bart van der Holst,
Ward B. Manchester,
Doga Can Su Ozturk,
Judit Szente,
Aleksandre Taktakishvili,
Gabor Tóth,
Meng Jin,
Tamas I. Gombosi
Abstract:
We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfvén wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In clos…
▽ More
We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfvén wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes.
The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfvén wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores.
△ Less
Submitted 18 November, 2020; v1 submitted 8 September, 2016;
originally announced September 2016.
-
Data Constrained Coronal Mass Ejections in A Global Magnetohydrodynamics Model
Authors:
M. Jin,
W. B. Manchester,
B. van der Holst,
I. Sokolov,
G. Toth,
R. E. Mullinix,
A. Taktakishvili,
A. Chulaki,
T. I. Gombosi
Abstract:
We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GO…
▽ More
We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b).
△ Less
Submitted 27 November, 2016; v1 submitted 17 May, 2016;
originally announced May 2016.
-
A Steady-State Picture of Solar Wind Acceleration and Charge State Composition Derived from a Global Wave-Driven MHD Model
Authors:
Rona Oran,
Enrico Landi,
Bart van der Holst,
Susan T. Lepri,
Alberto M. Vásquez,
Federico. A. Nuevo,
Richard Frazin,
Ward B. Manchester IV,
Igor V. Sokolov,
Tamas I. Gombosi
Abstract:
The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSo…
▽ More
The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{é}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography.
△ Less
Submitted 29 December, 2014;
originally announced December 2014.
-
Global Numerical Modeling of Energetic Proton Acceleration in a Coronal Mass Ejection Traveling through the Solar Corona
Authors:
Kamen A. Kozarev,
Rebekah M. Evans,
Nathan A. Schwadron,
Maher A. Dayeh,
Merav Opher,
Kelly E. Korreck,
Bart van der Holst
Abstract:
The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and par…
▽ More
The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.
△ Less
Submitted 9 June, 2014;
originally announced June 2014.
-
Alfvén Wave Turbulence as a Coronal Heating Mechanism: Simultaneously Predicting the Heating Rate and the Wave-Induced Emission Line Broadening
Authors:
R. Oran,
E. Landi,
B. van der Holst,
I. V. Sokolov,
T. I. Gombosi
Abstract:
In the present work, we test the predictions of the AWSoM model, a global extended-MHD model capable of calculating the propagation and turbulent dissipation of Alfvén waves in any magnetic topology, against high resolution spectra of the quiescent off-disk solar corona. Wave dissipation is the only heating mechanism assumed in this model. Combining 3D model results with the CHIANTI atomic databas…
▽ More
In the present work, we test the predictions of the AWSoM model, a global extended-MHD model capable of calculating the propagation and turbulent dissipation of Alfvén waves in any magnetic topology, against high resolution spectra of the quiescent off-disk solar corona. Wave dissipation is the only heating mechanism assumed in this model. Combining 3D model results with the CHIANTI atomic database, we were able to create synthetic line-of-sight spectra which include the effects of emission line broadening due to both thermal and wave-related non-thermal motions. To the best of our knowledge this is the first time a global model is used to obtain synthetic non-thermal line broadening. We obtained a steady-state solution driven by a synoptic magnetogram and compared the synthetic spectra with SUMER observations of a quiescent area above the solar west limb extending between 1.04 and 1.34 solar radii at the equator. Both the predicted line widths and the total line fluxes were consistent with the observations for 5 different ions. Using the 3D solution, we were able to locate the region that contributes the most to the emission used for measuring electron properties; we found that region to be a pseudo-streamer, whose modeled electron temperature and density are consistent with the measured ones. We conclude that the turbulent dissipation assumed in the AWSoM model can simultaneously account for the observed heating rate and the non-dissipated wave energy observed in this region.
△ Less
Submitted 2 January, 2014;
originally announced January 2014.
-
Alfven Wave Solar Model (AWSoM): Coronal Heating
Authors:
Bart van der Holst,
Igor V. Sokolov,
Xing Meng,
Meng Jin,
Ward B. Manchester IV,
Gabor Toth,
Tamas I. Gombosi
Abstract:
We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimen…
▽ More
We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently.
△ Less
Submitted 2 December, 2013; v1 submitted 16 November, 2013;
originally announced November 2013.
-
A Global Wave-Driven MHD Solar Model with a Unified Treatment of Open and Closed Magnetic Field Topologies
Authors:
R. Oran,
B. van der Holst,
E. Landi,
M. Jin,
I. V. Sokolov,
T. I. Gombosi
Abstract:
We describe, analyze and validate the recently developed Alfvén Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfvén waves. In this picture, heating and acceleration of the…
▽ More
We describe, analyze and validate the recently developed Alfvén Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfvén waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counter-propagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment of any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model, and no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady-state (in the co-rotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfvén waves, the perpendicular correlation length of the turbulence, and a pseudo-reflection coefficient. We compare model results for Carrington Rotation 2063 (November-December 2007) to remote observations in the EUV and X-ray ranges from STEREO, SOHO and Hinode spacecraft, as well as to in-situ measurements performed by Ulysses. The model results are in good agreement with observations. This is the first global model capable of simultaneously reproducing the multi-wavelength observations of the lower corona and the wind structure beyond Earth's orbit.
△ Less
Submitted 17 July, 2013;
originally announced July 2013.
-
MHD Waves and Coronal Heating: Unifying Empirical and MHD Turbulence Models
Authors:
Igor V. Sokolov,
Bart van der Holst,
Rona Oran,
Cooper Downs,
Ilia I. Roussev,
Meng Jin,
Ward B. Manchester IV,
Rebekah M. Evans,
Tamas I. Gombosi
Abstract:
We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only…
▽ More
We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.
△ Less
Submitted 5 December, 2012; v1 submitted 15 August, 2012;
originally announced August 2012.
-
Simulating the long-term evolution of radiative shocks in shock tubes
Authors:
B. van der Holst,
G. Toth,
I. V. Sokolov,
B. R. Torralva,
K. G. Powell,
R. P. Drake
Abstract:
We present the latest improvements in the Center for Radiative Shock Hydrodynamics (CRASH) code, a parallel block-adaptive-mesh Eulerian code for simulating high-energy-density plasmas. The implementation can solve for radiation models with either a gray or a multigroup method in the flux-limited-diffusion approximation. The electrons and ions are allowed to be out of temperature equilibrium and f…
▽ More
We present the latest improvements in the Center for Radiative Shock Hydrodynamics (CRASH) code, a parallel block-adaptive-mesh Eulerian code for simulating high-energy-density plasmas. The implementation can solve for radiation models with either a gray or a multigroup method in the flux-limited-diffusion approximation. The electrons and ions are allowed to be out of temperature equilibrium and flux-limited electron thermal heat conduction is included. We have recently implemented a CRASH laser package with 3-D ray tracing, resulting in improved energy deposition evaluation. New, more accurate opacity models are available which significantly improve radiation transport in materials like xenon. In addition, the HYPRE preconditioner has been added to improve the radiation implicit solver. With this updated version of the CRASH code we study radiative shock tube problems. In our set-up, a 1 ns, 3.8 kJ laser pulse irradiates a 20 micron beryllium disk, driving a shock into a xenon-filled plastic tube. The electrons emit radiation behind the shock. This radiation from the shocked xenon preheats the unshocked xenon. Photons traveling ahead of the shock will also interact with the plastic tube, heat it, and in turn this can drive another shock off the wall into the xenon. We are now able to simulate the long term evolution of radiative shocks.
△ Less
Submitted 6 June, 2012;
originally announced June 2012.
-
Buildup of Magnetic Shear and Free Energy During Flux Emergence and Cancellation
Authors:
Fang Fang,
Ward Manchester IV,
William P. Abbett,
Bart van der Holst
Abstract:
We examine a simulation of flux emergence and cancellation, which shows a complex sequence of processes that accumulate free magnetic energy in the solar corona essential for the eruptive events such as coronal mass ejections (CMEs), filament eruptions and flares. The flow velocity at the surface and in the corona shows a consistent shearing pattern along the polarity inversion line (PIL), which t…
▽ More
We examine a simulation of flux emergence and cancellation, which shows a complex sequence of processes that accumulate free magnetic energy in the solar corona essential for the eruptive events such as coronal mass ejections (CMEs), filament eruptions and flares. The flow velocity at the surface and in the corona shows a consistent shearing pattern along the polarity inversion line (PIL), which together with the rotation of the magnetic polarities, builds up the magnetic shear. Tether-cutting reconnection above the PIL then produces longer sheared magnetic field lines that extend higher into the corona, where a sigmoidal structure forms. Most significantly, reconnection and upward energy-flux transfer are found to occur even as magnetic flux is submerging and appears to cancel at the photosphere. A comparison of the simulated coronal field with the corresponding coronal potential field graphically shows the development of nonpotential fields during the emergence of the magnetic flux and formation of sunspots.
△ Less
Submitted 16 May, 2012;
originally announced May 2012.
-
Dynamic Coupling of Convective Flows and Magnetic Field during Flux Emergence
Authors:
Fang Fang,
Ward Manchester IV,
William P. Abbett,
Bart van der Holst
Abstract:
We simulate the buoyant rise of a magnetic flux rope from the solar convection zone into the corona to better understand the energetic coupling of the solar interior to the corona. The magnetohydrodynamic model addresses the physics of radiative cooling, coronal heating and ionization, which allow us to produce a more realistic model of the solar atmosphere. The simulation illustrates the process…
▽ More
We simulate the buoyant rise of a magnetic flux rope from the solar convection zone into the corona to better understand the energetic coupling of the solar interior to the corona. The magnetohydrodynamic model addresses the physics of radiative cooling, coronal heating and ionization, which allow us to produce a more realistic model of the solar atmosphere. The simulation illustrates the process by which magnetic flux emerges at the photosphere and coalesces to form two large concentrations of opposite polarities. We find that the large-scale convective motion in the convection zone is critical to form and maintain sunspots, while the horizontal converging flows in the near surface layer prevent the concentrated polarities from separating. The foot points of the sunspots in the convection zone exhibit a coherent rotation motion, resulting in the increasing helicity of the coronal field. Here, the local configuration of the convection causes the convergence of opposite polarities of magnetic flux with a shearing flow along the polarity inversion line. During the rising of the flux rope, the magnetic energy is first injected through the photosphere by the emergence, followed by energy transport by horizontal flows, after which the energy is subducted back to the convection zone by the submerging flows.
△ Less
Submitted 7 November, 2011;
originally announced November 2011.
-
Simulating radiative shocks in nozzle shock tubes
Authors:
B. van der Holst,
G. Toth,
I. V. Sokolov,
L. K. S. Daldorff,
K. G. Powell,
R. P. Drake
Abstract:
We use the recently developed Center for Radiative Shock Hydrodynamics (CRASH) code to numerically simulate laser-driven radiative shock experiments. These shocks are launched by an ablated beryllium disk and are driven down xenon-filled plastic tubes. The simulations are initialized by the two-dimensional version of the Lagrangian Hyades code which is used to evaluate the laser energy deposition…
▽ More
We use the recently developed Center for Radiative Shock Hydrodynamics (CRASH) code to numerically simulate laser-driven radiative shock experiments. These shocks are launched by an ablated beryllium disk and are driven down xenon-filled plastic tubes. The simulations are initialized by the two-dimensional version of the Lagrangian Hyades code which is used to evaluate the laser energy deposition during the first 1.1ns. The later times are calculated with the CRASH code. This code solves for the multi-material hydrodynamics with separate electron and ion temperatures on an Eulerian block-adaptive-mesh and includes a multi-group flux-limited radiation diffusion and electron thermal heat conduction. The goal of the present paper is to demonstrate the capability to simulate radiative shocks of essentially three-dimensional experimental configurations, such as circular and elliptical nozzles. We show that the compound shock structure of the primary and wall shock is captured and verify that the shock properties are consistent with order-of-magnitude estimates. The produced synthetic radiographs can be used for comparison with future nozzle experiments at high-energy-density laser facilities.
△ Less
Submitted 18 September, 2011;
originally announced September 2011.
-
Obtaining Potential Field Solution with Spherical Harmonics and Finite Differences
Authors:
Gabor Toth,
Bart van der Holst,
Zhenguang Huang
Abstract:
Potential magnetic field solutions can be obtained based on the synoptic magnetograms of the Sun. Traditionally, a spherical harmonics decomposition of the magnetogram is used to construct the current and divergence free magnetic field solution. This method works reasonably well when the order of spherical harmonics is limited to be small relative to the resolution of the magnetogram, although som…
▽ More
Potential magnetic field solutions can be obtained based on the synoptic magnetograms of the Sun. Traditionally, a spherical harmonics decomposition of the magnetogram is used to construct the current and divergence free magnetic field solution. This method works reasonably well when the order of spherical harmonics is limited to be small relative to the resolution of the magnetogram, although some artifacts, such as ringing, can arise around sharp features. When the number of spherical harmonics is increased, however, using the raw magnetogram data given on a grid that is uniform in the sine of the latitude coordinate can result in inaccurate and unreliable results, especially in the polar regions close to the Sun.
We discuss here two approaches that can mitigate or completely avoid these problems: i) Remeshing the magnetogram onto a grid with uniform resolution in latitude, and limiting the highest order of the spherical harmonics to the anti-alias limit; ii) Using an iterative finite difference algorithm to solve for the potential field. The naive and the improved numerical solutions are compared for actual magnetograms, and the differences are found to be rather dramatic.
We made our new Finite Difference Iterative Potential-field Solver (FDIPS) a publically available code, so that other researchers can also use it as an alternative to the spherical harmonics approach.
△ Less
Submitted 29 April, 2011;
originally announced April 2011.
-
Crash: A Block-Adaptive-Mesh Code for Radiative Shock Hydrodynamics - Implementation and Verification
Authors:
B. van der Holst,
G. Toth,
I. V. Sokolov,
K. G. Powell,
J. P. Holloway,
E. S. Myra,
Q. Stout,
M. L. Adams,
J. E. Morel,
R. P. Drake
Abstract:
We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and uses a flux limited diffusion approximation to recover the free-streaming limit. The electrons and ions are allowed to have different temperatures and we include a…
▽ More
We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and uses a flux limited diffusion approximation to recover the free-streaming limit. The electrons and ions are allowed to have different temperatures and we include a flux limited electron heat conduction. The radiation hydrodynamic equations are solved in the Eulerian frame by means of a conservative finite volume discretization in either one, two, or three-dimensional slab geometry or in two-dimensional cylindrical symmetry. An operator split method is used to solve these equations in three substeps: (1) solve the hydrodynamic equations with shock-capturing schemes, (2) a linear advection of the radiation in frequency-logarithm space, and (3) an implicit solve of the stiff radiation diffusion, heat conduction, and energy exchange. We present a suite of verification test problems to demonstrate the accuracy and performance of the algorithms. The CRASH code is an extension of the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code with this new radiation transfer and heat conduction library and equation-of-state and multigroup opacity solvers. Both CRASH and BATS-R-US are part of the publicly available Space Weather Modeling Framework (SWMF).
△ Less
Submitted 19 January, 2011;
originally announced January 2011.
-
Towards a Realistic, Data-Driven Thermodynamic MHD Model of the Global Solar Corona
Authors:
Cooper Downs,
Ilia I. Roussev,
Bart van der Holst,
Noé Lugaz,
Igor V. Sokolov,
Tamas I. Gombosi
Abstract:
In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnet…
▽ More
In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona.
△ Less
Submitted 11 December, 2009;
originally announced December 2009.
-
A multidimensional grid-adaptive relativistic magnetofluid code
Authors:
B. van der Holst,
R. Keppens,
Z. Meliani
Abstract:
A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing Lax-Friedrichs scheme in a finite volume setting and is combined with a diffusive approach for controlling magnetic monopole errors. The consistency between the…
▽ More
A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing Lax-Friedrichs scheme in a finite volume setting and is combined with a diffusive approach for controlling magnetic monopole errors. The consistency between the primitive and conservative variables is ensured at all limited reconstructions and the spatial part of the four velocity is used as a primitive variable. Demonstrative relativistic examples are shown to validate the implementation. We recover known exact solutions to relativistic MHD Riemann problems, and simulate the shock-dominated long term evolution of Lorentz factor 7 vortical flows distorting magnetic island chains.
△ Less
Submitted 4 July, 2008;
originally announced July 2008.
-
Extragalactic jets with helical magnetic fields: relativistic MHD simulations
Authors:
R. Keppens,
Z. Meliani,
B. van der Holst,
F. Casse
Abstract:
Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind…
▽ More
Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind the Mach disk, for various jet Lorentz factors and magnetic field helicities. We investigate the influence of helical magnetic fields on jet beam propagation in overdense external medium. We use the AMRVAC code, employing a novel hybrid block-based AMR strategy, to compute ideal plasma dynamics in special relativity. The helicity of the beam magnetic field is effectively transported down the beam, with compression zones in between diagonal internal cross-shocks showing stronger toroidal field regions. In comparison with equivalent low-relativistic jets which get surrounded by cocoons with vortical backflows filled by mainly toroidal field, the high speed jets demonstrate only localized, strong toroidal field zones within the backflow vortical structures. We find evidence for a more poloidal, straight field layer, compressed between jet beam and backflows. This layer decreases the destabilizing influence of the backflow on the jet beam. In all cases, the jet beam contains rich cross-shock patterns, across which part of the kinetic energy gets transferred. For the high speed reference jet considered here, significant jet deceleration only occurs beyond distances exceeding ${\cal O}(100 R_j)$, as the axial flow can reaccelerate downstream to the internal cross-shocks. This reacceleration is magnetically aided, due to field compression across the internal shocks which pinch the flow.
△ Less
Submitted 14 February, 2008;
originally announced February 2008.