Search NASA⌕ Search

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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

A Superposed Metric for Spinning Black Hole Binaries Approaching Merger

We construct an approximate metric that represents the spacetime of spinning binary black holes (BBH) approaching merger. We build the metric as an analytical superposition of two Kerr metrics in harmonic coordinates, where we transform each black hole term with time-dependent boosts describing an inspiral trajectory. The velocities and trajectories of the boost are obtained by solving the post-Newtonian (PN) equations of motion at 3.5 PN order. We analyze the spacetime scalars of the new metric and we show that it is an accurate approximation of Einstein’s field equations in vacuum for a BBH system in the inspiral regime. Furthermore, to prove the effectiveness of our approach, we test the metric in the context of a 3D general relativistic magnetohydrodynamical (GRMHD) simulation of accreting minidisks around the black holes. We compare our results with a previous well-tested spacetime construction based on the asymptotic matching method. We conclude that our new spacetime is well-suited for long-term GRMHD simulations of spinning binary black holes on their way to the merger.

Luciano Combi↗
NASA NTRS2022

Electromagnetic signatures from supermassive binary black holes approaching merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D General Relativistic Magnetohydrodynamic (GRMHD) simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller ISCO angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

Eduardo M. Gutiérrez↗
NASA NTRS2022

Spectral Power-law Formation by Sequential Particle Acceleration in Multiple Flare Magnetic Islands

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

magnetic reconnection↗
NASA NTRS2022

The Role of Reconnection in the Onset of Solar Eruptions

Solar eruptive events such as coronal mass ejections and eruptive flares are frequently associated with the emergence of magnetic flux from the convection zone into the corona. We use three-dimensional magnetohydrodynamic numerical simulations to study the interaction of coronal magnetic fields with emerging flux and determine the conditions that lead to eruptive activity. A simple parameter study is performed, varying the relative angle between emerging magnetic flux and a preexisting coronal dipole field. We find that in all cases the emergence results in a sheared magnetic arcade that transitions to a twisted coronal flux rope via low-lying magnetic reconnection. This structure, however, is constrained by its own outer field and so is noneruptive in the absence of reconnection with the overlying coronal field. The amount of this overlying reconnection is determined by the relative angle between the emerged and preexisting fields. The reconnection between emerging and preexisting fields is necessary to generate sufficient expansion of the emerging structure so that flare-like reconnection below the coronal flux rope becomes strong enough to trigger its release. Our results imply that the relative angle is the key parameter in determining whether the resultant active regions exhibit eruptive behavior and is thus a potentially useful candidate for predicting eruptions in newly emerging active regions. More generally, our results demonstrate that the detailed interaction between the convection zone/photosphere and the corona must be calculated self-consistently in order to model solar eruptions accurately.

James E Leake↗
NASA NTRS2023

Interchange Reconnection Dynamics in A Solar Coronal Pseudo-Streamer

Context. The generation of the slow solar wind remains an open problem in heliophysics. One of the current theories among those aimed at explaining the injection of coronal plasma in the interplanetary medium is based on interchange reconnection. It assumes that the exchange of magnetic connectivity between closed and open fields allows the injection of coronal plasma in the interplanetary medium to travel along the newly reconnected open field. However, the exact mechanism underlying this effect is still poorly understood. Aims. Our objective is to study this scenario in a particular magnetic structure of the solar corona: a pseudo-streamer. This topological structure lies at the interface between open and closed magnetic field and is thought to be involved in the generation of the slow solar wind. Methods. We performed innovative 3D magnetohydrodynamic (MHD) simulations of the solar corona with a pseudo-streamer, using the Adaptively Refined MHD Solver (ARMS). By perturbing the quasi-steady ambient state with a simple photospheric, large-scale velocity flow, we were able to generate a complex dynamics of the open-and-closed boundary of the pseudo-streamer. We studied the evolution of the connectivity of numerous field lines to understand its precise dynamics. Results. We witnessed different scenarios of opening of the magnetic field initially closed under the pseudo-streamer: one-step interchange reconnection dynamics, along with more complex scenarios, including a coupling between pseudo-streamer and helmet streamer, as well as back-and-forth reconnections between open and closed connectivity domains. Finally, our analysis revealed large-scale motions of a newly opened magnetic field high in the corona that may be explained by slipping reconnection. Conclusions. By introducing a new analysis method for the magnetic connectivity evolution based on distinct closed-field domains, this study provides an understanding of the precise dynamics underway during the opening of a closed field, which enables the injection of closed-field, coronal plasma in the interplanetary medium. Further studies shall provide synthetic observations for these diverse outgoing flows, which could be measured by Parker Solar Probe and Solar Orbiter.

magnetic fields↗
NASA NTRS2024

Disentangling the Entangled Linkages of Magnetic Helicity

Magnetic helicity, H, measures magnetic linkages in a volume. The early theoretical development of helicity focused on magnetically closed systems in V bounded by ∂V. For magnetically closed systems, V ∈ R -3 = V + V∗, no magnetic flux threads the boundary, n̂ · B|∂V = 0. Berger & Field (1984) and Finn & Antonsen (1985) extended the definition of helicity to relative helicity, H, for magnetically open systems where magnetic flux may thread the boundary. Berger (1999, 2003) expressed this relative helicity as two gauge invariant terms that describe the self helicity of magnetic field that closes inside V and the mutual helicity between the magnetic field that threads the boundary ∂V and the magnetic field that closes inside V. The total magnetic field that permeates V entangles magnetic fields that are produced by current sources J in V with magnetic fields that are produced by current sources J∗ in V∗. Building on this fact, we extend Berger’s expressions for relative magnetic helicity to eight gauge invariant quantities that simultaneously characterize both of these self and mutual helicities and attribute their origins to currents J in V and/or J∗ in V∗, thereby disentangling the domain of origin for these entangled linkages. We arrange these eight terms into novel expressions for internal and external helicity (self) and internal-external helicity (mutual) based on their domain of origin. The implications of these linkages for interpreting magnetic energy is discussed and new boundary observables are proposed for tracking the evolution of the field that threads the boundary.

Analytical mathematics↗
NASA NTRS2023

The Complexity of the Day-Side X-Line During Southward Interplanetary Magnetic Field

High-resolution global magnetohydrodynamics (MHD) simulations include both meso- and global-scale processes occurring at the magnetopause, which interact to determine the time-dependent orientation of the day-side x-line (DXL). This study demonstrates that the global orientation of the DXL in GAMERA global MHD simulations varies on a time scale of minutes during steady southward interplanetary magnetic field conditions. This behavior manifests in observational data when reconnection outflows indicate that the direction to the x-line is opposite to the prediction from a steady-state model of the reconnection location. Because steady-state models of the DXL do not capture dynamics that are independent of solar wind variations, particularly surface waves and flux transfer events, they represent a time-averaged state of the system.

Brandon Burkholder↗
NASA NTRS2024

Exploring Standing and Reflected Slow-mode Waves in Flaring Coronal Loops: A Parametric Study Using 2.5D MHD Modeling

Abstract Recent observations of reflected propagating and standing slow-mode waves in hot flaring coronal loops have spurred our investigation into their underlying excitation and damping mechanisms. To understand these processes, we conduct 2.5D magnetohydrodynamic (MHD) simulations using an arcade active region model that includes a hot and dense loop. Our simulations allow for in-depth parametric investigations complementing and expanding our previous 3D MHD modeling results. We excite these waves using a large-amplitude, flow pulse applied at one footpoint of the loop in two distinct models as motivated by observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). The first model (Model 1) incorporates classical compressive viscosity coefficient, while the second model (Model 2) adopts a 10-times enhanced viscosity coefficient. We obtain the following major results: (1) Our 2.5D MHD simulations reinforce previous conclusions derived from 1D and 3D MHD models that significantly enhanced viscosity is crucial for the rapid excitation of standing slow waves with damping times consistent with observations by Wang et al. (2015). (2)We uncover that nonlinearity in Model 1 delays the conversion of a reflected propagating wave into a standing wave. In contrast, Model 2 exhibits a much weak influence of nonlinearity on the excitation time of standing waves, thanks to the suppression of these effects by enhanced viscosity. (3) Our results reveal that the transverse temperature structure holds more influence on wave behavior than the density structure. In Model 1, increased loop temperature contrast significantly enhances wave trapping within the structure, mitigating the impact of temperature-dependent viscous damping. Conversely, in Model 2, the impact of temperature structure on wave behavior weakens in comparison to the effect of viscosity. (4) Model 1 displays evident nonlinear coupling to the fast and kink magnetoacoustic waves and pronounced wave leakage into the corona. Model 2 exhibits significantly weaker effects in this regard. Analyzing three observed wave events by SDO/AIA aligns with Model 2 predictions, providing further support for the substantial viscosity increase. Our 2.5D study unravels the complex interplay of wave-flow phenomena and nonlinear processes in coronal loops, extending our previous 1D modeling results to incorporate more realistic loop geometry. This provides insights into scenarios where 3D effects may be neglected, thereby enhancing our understanding of the intricate dynamics of the solar corona.

Tongjiang Wang↗
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗
NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗
NASA NTRS1984

Interaction of Convection and Small-scale Magnetic Fields: Influence on the Solar Luminosity

Changes in the local solar luminosity due to the presence of a small scale structured (facular) magnetic field in the photosphere are discussed. The discussion is based on three dimensional numerical simulations of the magnetohydrodynamics of the top of the convection zone, and the adjacent stable photosphere. The simulations demonstrate that practically all of the magnetic flux present is concentrated into intense magnetic flux structures, such that the magnetic field pressure is balanced by the gas pressure of the surrounding plasma. The flux concentration is caused by the convectively unstable stratification. The average luminosity of the area is influenced by three effects: (1) the brightness of the flux concentrations, (2) their filling factor, and (3) the average luminosity of the surrounding plasma.

Nordlund, A.↗
NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗
NASA NTRS1987

Comparison of MHD simulation for the February 1986 events with interplanetary observations by the spacecraft Sakigake

During the period 3-10 February 1986 a series of major solar flares occurred on the Sun and several intense geomagnetic storms took place on the Earth. To examine the causality between the solar activity and the geomagnetic activity in this period, a magnetohydrodynamic (MHD) numerical simulation was performed using a 2 1/2 -D numerical code. In that period of February 1986, the Japanese spacecraft Sakigake was at 0.84 AU, 57 deg west of the Earth. Besides the in-situ measurements of the interplanetary plasma, Sakigake also provided Doppler scintillation observations. Comparisons between the results of the MHD simulation and the measurements made by the spacecraft Sakigake are presented.

Smith, Z.↗
NASA NTRS2009

Cavities of Weak Magnetic Field Strength in the Wake of FTEs: Results from Global Magnetospheric MHD Simulations

We use the global magnetohydrodynamic (MHD) code BATS-R-US to model multipoint observations of Flux Transfer Event (FTE) signatures. Simulations with high spatial and temporal resolution predict that cavities of weak magnetic field strength protruding into the magnetosphere trail FTEs. These predictions are consistent with recently reported multi-point Cluster observations of traveling magnetopause erosion regions (TMERs).

Kuznetsova, M. M.↗