How to Turn Jets into Cylinders near Supermassive Black Holes in 3D General Relativistic Magnetohydrodynamic Simulations
Explore the source record for details and available documents.
Valeriia Rohoza↗
SEARCH · Search NASA
Results for “magnetohydrodynamical simulations”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Analytic Model and Magnetohydrodynamic Simulations of Three-dimensional Magnetic Switchbacks
Explore the source record for details and available documents.
Chen 辰 Shi 时↗
Kinetic Simulations of Current-Sheet Formation and Reconnection at a Magnetic X Line
The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.
Black, C.↗
Numerical simulation of magnetohydrodynamic shock propagation in the corona
Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail. The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfven waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.
Wu, S. T.↗
The Substorm Current Wedge: Further Insights from MHD Simulations
Using a recent magnetohydrodynamic simulation of magnetotail dynamics, we further investigate the buildup and evolution of the substorm current wedge (SCW), resulting from flow bursts generated by near-tail reconnection. Each flow burst generates an individual current wedge, which includes the reduction of cross-tail current and the diversion to region 1 (R1)-type field-aligned currents (earthward on the dawn and tailward on the duskside), connecting the tail with the ionosphere. Multiple flow bursts generate initially multiple SCW patterns, which at later times combine to a wider single SCW pattern. The standard SCWmodel is modified by the addition of several current loops, related to particular magnetic field changes: the increase of Bz in a local equatorial region (dipolarization), the decrease of |Bx| away from the equator (current disruption), and increases in |By| resulting from azimuthally deflected flows. The associated loop currents are found to be of similar magnitude, 0.1-0.3 MA. The combined effect requires the addition of region 2 (R2)-type currents closing in the near tail through dawnward currents but also connecting radially with the R1 currents. The current closure at the inner boundary, taken as a crude proxy of an idealized ionosphere, demonstrates westward currents as postulated in the original SCW picture as well as North-South currents connecting R1- and R2-type currents, which were larger than the westward currents by a factor of almost 2. However, this result should be applied with caution to the ionosphere because of our neglect of finite resistance and Hall effects.
reconnection↗
Current-Sheet Formation and Reconnection at a Magnetic X Line in Particle-in-Cell Simulations
The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.
Black, C.↗
Contribution of ULF Wave Activity to the Global Recovery of the Outer Radiation Belt During the Passage of a High‐Speed Solar Wind Stream Observed in September 2014
Energy coupling between the solar wind and the Earth's magnetosphere can affect the electron population in the outer radiation belt. However, the precise role of different internal and external mechanisms that leads to changes of the relativistic electron population is not entirely known. This paper describes how ultralow frequency (ULF) wave activity during the passage of Alfvenic solar wind streams contributes to the global recovery of the relativistic electron population in the outer radiation belt. To investigate the contribution of the ULF waves, we searched the Van Allen Probes data for a period in which we can clearly distinguish the enhancement of electron fluxes from the background. We found that the global recovery that started on 22 September 2014, which coincides with the corotating interaction region preceding a high‐speed stream and the occurrence of persistent substorm activity, provides an excellent scenario to explore the contribution of ULF waves. To support our analyses, we employed ground‐ and space‐based observational data and global magnetohydrodynamic simulations and calculated the ULF wave radial diffusion coefficients employing an empirical model. Observations show a gradual increase of electron fluxes in the outer radiation belt and a concomitant enhancement of ULF activity that spreads from higher to lower L‐shells. Magnetohydrodynamic simulation results agree with observed ULF wave activity in the magnetotail, which leads to both fast and Alfven modes in the magnetospheric nightside sector. The observations agree with the empirical model and are confirmed by phase space density calculations for this global recovery period.
Silva, L. A. Da↗
Magnetosheath Propagation Time of Solar Wind Directional Discontinuities
Observed delays in the ground response to solar wind directional discontinuities have been explained as the result of larger than expected magnetosheath propagation times. Recently, Samsonov et al. (2017, https://doi.org/10.1002/2017GL075020) showed that the typical time for a southward interplanetary magnetic field (IMF) turning to propagate across the magnetosheath is 14 min. Here by using a combination of magnetohydrodynamic simulations, spacecraft observations, and analytic calculations, we study the dependence of the propagation time on solar wind parameters and near-magnetopause cutoff speed. Increases in the solar wind speed result in greater magnetosheath plasma flow velocities, decreases in the magnetosheath thickness and, as a result, decreases in the propagation time. Increases in the IMF strength result in increases in the magnetosheath thickness and increases in the propagation time. Both magnetohydrodynamic simulations and observations suggest that propagation times are slightly smaller for northward IMF turnings. Magnetosheath flow deceleration must be taken into account when predicting the arrival times of solar wind structures at the dayside magnetopause.
Samsonov, A. A.↗
Global magnetohydrodynamic magnetosphere simulation with an adaptively embedded particle-in-cell model
Explore the source record for details and available documents.
Xiantong Wang↗
Lattice Boltzmann model for simulation of magnetohydrodynamics
A numerical method, based on a discrete Boltzmann equation, is presented for solving the equations of magnetohydrodynamics (MHD). The algorithm provides advantages similar to the cellular automaton method in that it is local and easily adapted to parallel computing environments. Because of much lower noise levels and less stringent requirements on lattice size, the method appears to be more competitive with traditional solution methods. Examples show that the model accurately reproduces both linear and nonlinear MHD phenomena.
Chen, Shiyi↗
Modeling inner boundary values at 18 solar radii during solar quiet time for global three-dimensional time-dependent magnetohydrodynamic numerical simulation
Explore the source record for details and available documents.
Chin-Chun Wu↗
Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling
We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.
Solar prominences↗
Variability of the Reconnection Guide Field in Solar Flares
Solar flares may be the best-known examples of the explosive conversion of magnetic energy into bulk motion, plasma heating, and particle acceleration via magnetic re-connection. The energy source for all flares is the highly sheared magnetic field of a filament channel above a polarity inversion line (PIL). During the flare, this shear field becomes the so-called reconnection guide field (i.e., the non-reconnecting component), which has been shown to play a major role in determining key properties of the recon-nection including the efficiency of particle acceleration. We present new high-resolution, three-dimensional, magnetohydrodynamics simulations that reveal the detailed evolution of the magnetic shear/guide field throughout an eruptive flare. The magnetic shear evolves in three distinct phases: shear first builds up in a narrow region about the PIL, then expands outward to form a thin vertical current sheet, and finally is transferred by flare reconnection into an arcade of sheared flare loops and an erupting flux rope. We demonstrate how the guide field may be inferred from observations of the sheared flare loops. Our results indicate that initially the guide field is larger by about a factor of 5 than the reconnecting component, but it weakens by more than an order of magnitude over the course of the flare. Instantaneously, the guide field also varies spatially over a similar range along the three-dimensional current sheet. We discuss the implications of the remarkable variability of the guide field for the timing and localization of efficient particle acceleration in flares.
Solar magnetic reconnection↗
Relative timing of substorm features in MHD simulations
An investigation of the temporal sequence of substorm phenomena based on three dimensional MHD (magnetohydrodynamic) simulations of magnetic reconnection and plasmoid formation is presented. The investigation utilizes a spatially localized resistivity model which leads to a significantly faster evolution than found in previous investigations. The analysis of the results concentrates on substorm features that have received considerable attention in the past. The formation of magnetic neutral lines, the occurrence of fast flows directed both earthward and tailward, and the magnetic field changes leading to the formation of the substorm current wedge, and to the depolarization of the magnetic field earthward of the reconnection region and its dependence on the spatial distribution of resistivity, are discussed. These phenomena are seen as an integral part of the nonlinear evolution of the three dimensional tearing instability.
Hesse, Michael↗
Global Simulation of Proton Precipitation Due to Field Line Curvature During Substorms
The low latitude boundary of the proton aurora (known as the Isotropy Boundary or IB) marks an important boundary between empty and full downgoing loss cones. There is significant evidence that the IB maps to a region in the magnetosphere where the ion gyroradius becomes comparable to the local field line curvature. However, the location of the IB in the magnetosphere remains in question. In this paper, we show simulated proton precipitation derived from the Field Line Curvature (FLC) model of proton scattering and a global magnetohydrodynamic simulation during two substorms. The simulated proton precipitation drifts equatorward during the growth phase, intensifies at onset and reproduces the azimuthal splitting published in previous studies. In the simulation, the pre-onset IB maps to 7-8 RE for the substorms presented and the azimuthal splitting is caused by the development of the substorm current wedge. The simulation also demonstrates that the central plasma sheet temperature can significantly influence when and where the azimuthal splitting takes place.
Field Line Curvature↗
Global Simulation of Proton Precipitation Due to Field Line Curvature During Substorms
The low latitude boundary of the proton aurora (known as the Isotropy Boundary or IB) marks an important boundary between empty and full downgoing loss cones. There is significant evidence that the IB maps to a region in the magnetosphere where the ion gyroradius becomes comparable to the local field line curvature. However, the location of the IB in the magnetosphere remains in question. In this paper, we show simulated proton precipitation derived from the Field Line Curvature (FLC) model of proton scattering and a global magnetohydrodynamic simulation during two substorms. The simulated proton precipitation drifts equatorward during the growth phase, intensifies at onset and reproduces the azimuthal splitting published in previous studies. In the simulation, the pre-onset IB maps to 7-8 RE for the substorms presented and the azimuthal splitting is caused by the development of the substorm current wedge. The simulation also demonstrates that the central plasma sheet temperature can significantly influence when and where the azimuthal splitting takes place.
Gilson, M. L.↗
Simulations of Emerging Magnetic Flux. II. The Formation of Unstable Coronal Flux Ropes and the Initiation of Coronal Mass Ejections
We present results from three-dimensional magnetohydrodynamic simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain (approximately 36 Mm above the surface).We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include surface shearing and rotational motions, quadrupolar geometry above the surface, central sheared arcades reconnecting with oppositely orientated overlying dipole fields, the formation of coronal flux ropes underlying potential coronal field, and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout," are operating during the emergence of new active regions.
corona↗
Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars
In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.
Mathews, Grant J.↗