Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 217 records · Page 12

Dynamical Unification of Tidal Disruption Events

The ~100 tidal disruption events (TDEs) observed so far exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs radiate predominantly at X-ray energies, while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. In particular, for optical TDEs, while their UV and optical emissions decline somewhat following the fallback pattern, some events can greatly rebrighten in X-rays at late time. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates in the regime of super-Eddington accretion. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al., which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation. Furthermore, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect that more X-rays will leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiative Reconnection-powered TeV Flares from the Black Hole Magnetosphere in M87

Abstract Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.

79 ASTRONOMY AND ASTROPHYSICS↗

First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization

Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here we present the first analysis of circular polarization using EHT data, acquired in 2017, which can potentially provide additional insights into the magnetic fields and plasma composition near the black hole. Interferometric closure quantities provide convincing evidence for the presence of circularly polarized emission on event-horizon scales. We produce images of the circular polarization using both traditional and newly developed methods. All methods find a moderate level of resolved circular polarization across the image ( $\langle$| v |$\rangle$ < 3.7%), consistent with the low image-integrated circular polarization fraction measured by the Atacama Large Millimeter/submillimeter Array (| v int |< 1%). Despite this broad agreement, the methods show substantial variation in the morphology of the circularly polarized emission, indicating that our conclusions are strongly dependent on the imaging assumptions because of the limited baseline coverage, uncertain telescope gain calibration, and weakly polarized signal. We include this upper limit in an updated comparison to general relativistic magnetohydrodynamic simulation models. This analysis reinforces the previously reported preference for magnetically arrested accretion flow models. We find that most simulations naturally produce a low level of circular polarization consistent with our upper limit and that Faraday conversion is likely the dominant production mechanism for circular polarization at 230 GHz in M87*.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetized Outflows from Short-lived Neutron Star Merger Remnants Can Produce a Blue Kilonova

Abstract We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino–matter interactions and is well suited to studying the resulting nucleosynthesis and kilonova emission. A magnetized wind is driven from the remnant and ejects neutron-rich material at a quasi-steady-state rate of 0.8 × 10 −1 M ⊙ s −1 . We find that the ejecta in our simulations underproduce r -process abundances beyond the second r -process peak. For sufficiently long-lived remnants, these outflows alone can produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.

Astronomy & Astrophysics↗

Subphotospheric Emission from Short Gamma-Ray Bursts: Protons Mold the Multimessenger Signals

Abstract The origin of the observed Band-like photon spectrum in short gamma-ray bursts (sGRBs) is a long-standing mystery. We carry out the first general relativistic magnetohydrodynamic simulation of an sGRB jet with initial magnetizationσ 0 = 150 in dynamical ejecta from a binary merger. From this simulation, we identify regions along the jet of efficient energy dissipation due to magnetic reconnection and collisionless subshocks. Taking into account electron and proton acceleration processes, we solve for the first time the coupled transport equations for photons, electrons, protons, neutrinos, and intermediate particle species up to close to the photosphere (i.e., up to 1 × 10 12 cm), accounting for all relevant radiative and cooling processes. We find that the subphotospheric multimessenger signals carry strong signatures of the hadronic interactions and their resulting particle cascades. Importantly, the spectral energy distribution of photons is significantly distorted with respect to the Wien one, commonly assumed below the photosphere. Our findings suggest that the bulk of the nonthermal photon spectrum observed in sGRBs can stem from hadronic processes occurring below the photosphere and previously neglected, with an accompanying energy flux of neutrinos peaking in the GeV energy range.

Astronomy & Astrophysics↗

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

Abstract We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

Astronomy & Astrophysics↗

Magnetically Arrested Circumbinary Accretion Flows

Abstract Binary systems with comparable masses and a surrounding accretion disk can accrete gas through spiral accretion streams penetrating the central cavity formed by tidal interactions. Using three-dimensional Newtonian magnetohydrodynamics simulations, we investigate the possibility of a magnetically arrested accretion flow through the cavity. Rather than solely continuously feeding the binary through spiral accretion streams, the accretion is regulated by the strong magnetic field inside the cavity. Transport of mass and angular momentum onto the binary then proceeds largely periodically in magnetic flux eruption episodes. The ejected flux tubes carry angular momentum outward and away from the binary, inject hot plasma into the disk, and can launch flares. This likely intermittent scenario could have potential implications for the emission signatures of supermassive black hole binaries and shed light onto the role magnetic fields play in the binary’s orbital evolution.

Astronomy & Astrophysics↗

Measuring Black Hole Light Echoes with Very Long Baseline Interferometry

Light passing near a black hole can follow multiple paths from an emission source to an observer due to strong gravitational lensing. Photons following different paths take different amounts of time to reach the observer, which produces an echo signature in the image. The characteristic echo delay is determined primarily by the mass of the black hole, but it is also influenced by the black hole spin and inclination to the observer. In the Kerr geometry, echo images are demagnified, rotated, and sheared copies of the direct image and lie within a restricted region of the image. Echo images have exponentially suppressed flux, and temporal correlations within the flow make it challenging to directly detect light echoes from the total light curve. In this Letter, we propose a novel method to search for light echoes by correlating the total light curve with the interferometric signal at high spatial frequencies, which is a proxy for indirect emission. We explore the viability of our method using numerical general relativistic magnetohydrodynamic simulations of a near-face-on accretion system scaled to M87-like parameters. We demonstrate that our method can be used to directly infer the echo delay period in simulated data. An echo detection would be clear evidence that we have captured photons that have circled the black hole, and a high-fidelity echo measurement would provide an independent measure of fundamental black hole parameters. Our results suggest that detecting echoes may be achievable through interferometric observations with a modest space-based very long baseline interferometry mission.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-Ray Bursts and Kilonovae from the Accretion-induced Collapse of White Dwarfs

Abstract We present the first seconds-long 2D general relativistic neutrino magnetohydrodynamic simulations of accretion-induced collapse (AIC) in rapidly rotating, strongly magnetized white dwarfs (WDs), which might originate as remnants of double-WD mergers. This study examines extreme combinations of magnetic fields and rotation rates, motivated both by the need to address the limitations of 2D axisymmetric simulations and to explore the physics of AIC under rare conditions that, while yet to be observationally confirmed, may be consistent with current theoretical models and account for unusual events. Under these assumptions, our results demonstrate that, if realizable, such systems can generate relativistic jets and neutron-rich outflows with properties consistent with long gamma-ray bursts (LGRBs) accompanied by kilonovae, such as GRB 211211A and GRB 230307A. These findings highlight the potential role of AIC in heavy r -process element production and offer a framework for understanding rare LGRBs associated with kilonova emission. Longer-duration 3D simulations are needed to fully capture magnetic field amplification, resolve instabilities, and determine the fate of the energy retained by the magnetar at the end of the simulations.

Cheong 張, Patrick Chi-Kit 志杰 (ORCID:00000003144933↗

Evolution of the solar wind structure in the outer heliosphere

Shocks and interaction regions play very important roles in the evolution of large-scale solar wind structure in the outer heliosphere. This study is based on (1) plasma and magnetic field data observed from Voyager and Pioneer spacecraft, and (2) a quantitative magnetohydrodynamic simulation model. Interaction regions bounded by a forward and a reverse shock begin to form near 1 AU at the leading edges of a large-scale stream. The total pressure in the region is greater than the ambient pressure by a factor of ten or more. Large jumps in pressure remain as a prominant feature of the interplanetary structure even as the jumps in flow speed become less visible in the outer heliosphere. The propagation of the forward and reverse shocks widens the dimension of an interaction region. As a result, two interaction regions belonging to neighboring streams coalesce to form a merged interaction region (MIR). Collision and merging of shocks take place during the coalescence process. Two MIRs can themselves merge again at greater heliocentric distances. Simulation results agree well with spacecraft observations, and they explain major restructuring of the solar wind in the outer heliosphere.

Whang, Y. C.↗

The nonlinear evolution of magnetized solar filaments

Thermal instability driven by optically thin radiation is believed to initiate the formation of plasma filaments in the solar corona. The fact that filaments are observed generally to separate regions of opposite, line-of-sight, magnetic polarity in the underlying photosphere suggests that filament formation requires the presence of a highly sheared, local magnetic field. Two-dimensional, nonlinear, magnetohydrodynamic simulations of the local genesis and growth of solar filaments in a force-free, sheared, magnetic field were performed, and the evolution of generic perturbations possessing broad spatial profiles was traced. It was found that simulations of the evolution of initial random-noise perturbations produce filamentary plasma structures that exhibit densities and temperatures characteristic of observed solar filaments. Furthermore, in each of these simulations, the filament axis lies at a finite angle with respect to the local magnetic field, consistent with solar observations.

Sparks, L.↗

Magnetic field draping at the low-latitude magnetopause

Magnetohydrodynamic simulations are used to investigate the structure of the low-latitude magnetopause for interplanetary magnetic field conditions with a dominant southward component. The structure is self-consistently calculated as an initial-value problem in which the system is allowed to evolve into a quasi-steady state. All components of the 3D velocity and magnetic field as well as compressibility, resistivity, and viscosity are included in the 2D calculation. The simulation model shows that magnetic field draping can occur at the magnetopause boundary when magnetic merging takes place in the presence of a tangential shear flow. For 'normal' (positive Bx) draping, the higher-latitude portion of the field lines are curved toward the sun on the magnetospheric side of the magnetopause and away from the sun on the magnetosheath side. The thickness of the normal draping structure scales with the viscosity. The field-aligned current system that accompanies normal magnetic draping is consistent with the sense of the region 1 currents that flow into the dayside ionosphere.

Richard, R. L.↗

A parametric study of the solar wind interaction with comets

The Naval Research Laboratory's magnetohydrodynamic simulation code is used to simulate the solar wind interction with comet Halley for two different outgassing rates and several different solar wind states. The magnetic field is more strongly draped for fast solar wind conditions than slow. For higher mass loading rates, the tail becomes wider and contains more magnetic flux. The visual appearance of the comet differs for the case in which the interplanetary magnetic field lies in the plane of the sky from the case when it lies along the line of sight. The ion tail appears shorter in the latter case. Thus variation in the IMF direction can cause significant changes in the appearance of comets. The comet also creates a large momentum flux deficit in the solar wind with a narrow enhanced region within it corresponding to the ion tail.

Russell, C. T.↗

Anomalous aspects of magnetosheath flow and of the shape and oscillations of the magnetopause during an interval of strongly northward interplanetary magnetic field

On 15 Feb. 1978, the orientation of the interplanetary magnetic field (IMF) remained steadily northward for more than 12 hours. The ISEE 1 and 2 spacecraft were located near apogee on the dawn side flank of the magnetotail. IMP 8 was almost symmetrically located in the magnetosheath on the dusk flank and IMP 7 was upstream in the solar wind. Using plasma and magnetic field data, we show the following: (1) the magnetosheath flow speed on the flanks of the magnetotail steadily exceeded the solar wind speed by 20 percent; (2) surface waves with approximately a 5-min period and very non-sinusoidal waveform were persistently present on the dawn magnetopause and waves of similar period were present in the dusk magnetosheath; and (3) the magnetotail ceased to flare at an antisunward distance of 15 R(sub E). We propose that the acceleration of the magnetosheath flow is achieved by magnetic tension in the draped field configuration for northward IMF and that the reduction of tail flaring is consistent with a decreased amount of open magnetic flux and a larger standoff distance of the subsolar magnetopause. Results of a three-dimensional magnetohydrodynamic simulation support this phenomenological model.

Chen, Sheng-Hsien↗

Numerical modeling of the energy storage and release in solar flares

This paper reports on investigation of the photospheric magnetic field-line footpoint motion (usually referred to as shear motion) and magnetic flux emerging from below the surface in relation to energy storage in a solar flare. These causality relationships are demonstrated by using numerical magnetohydrodynamic simulations. From these results, one may conclude that the energy stored in solar flares is in the form of currents. The dynamic process through which these currents reach a critical value is discussed as well as how these currents lead to energy release, such as the explosive events of solar flares.

Wu, S. T.↗

Overview of magnetospheric research

During the funding period for NASA Grant NAG5-1480 which prior to December, 1, 1990 was known as NASA Grant NAGW-78, the group has made substantial progress on the various topics originally proposed. The research performed has resulted in two Ph.D. Theses, more than 50 refereed papers in various journals and conference proceedings, and 31 invited and 104 contributed talks at conferences and symposia throughout the world. The main results from this work are summarized in each of the sections outlined in the original proposal, followed by a complete list of the group publications associated with this grant, a list of all invited talks given during the last three years, and finally a listing of the contributed talks. Research topics include global magnetohydrodynamic simulations, structure of the dayside magnetopause, merging at the dayside magnetopause, polar wind, plasma waves in the distant magnetotail, slow shocks in the distant magnetotail, collisionless tearing instabilities in the magnetotail, and other problems.

Ashour-Abdalla, Maha↗

Magnetotail instability

The stability of the geomagnetic tail is investigated on the basis of three dimensional resistive magnetohydrodynamic simulations, using different dynamic constraints and different initial equilibria. Different forms of the energy equation for isotropic pressure are found to have no significant effect on the dynamic growth of a resistive tearing instability, which is responsible for near Earth reconnection, plasmoid formation and ejection, and the generation of fast plasma flows. The constraints of a modified double adiabatic approach, however, can quench the tearing instability through the development of large, mirror type, anisotropies in the boundary regions of the plasma sheet, unless isotropization occurs on fast, nearly Alfvenic, time scales. The presence of a net cross tail magnetic field component B(sub yN) can reduce the growth of the instability without complete stabilization. An increase of B(sub z) from midnight toward the tail flanks, however, by more than a factor of about 3, apparently completely stabilizes the tearing mode. Stabilization and destabilization thus may depend on properties and constraints (and their release) in regions other than the neutral sheet where reconnection is initiated.

Birn, Joachim↗

The deformation of flux tubes in the solar wind with applications to the structure of magnetic clouds and CMEs

Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Cargill, Peter J.↗