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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1994

MHD simulations of boundary layer formation along the dayside Venus ionopause due to mass loading

A two-dimensional magnetohydrodynamic (MHD) simulation of mass-loaded solar wind flow around the dayside of Venus is presented. For conditions appropriate to a low-altitude ionopause the simulations show that mass loading from the pickup of oxygen ions produces a boundary layer of finite thickness along the ionopause. Within this layer the temperatures exhibit strong gradients normal to and away from the ionopause. Furthermore, there is a shear in the bulk flow velocity across the boundary layer, such that the (predominantly tangential) flow decreases in speed as the ionopause is approached and remains small along the ionopause, consistent with Pioneer Venus observations. The total mass density increases significantly as the flow approaches the ionopause, where the contribution of O(+) to the total number density is a few percent. Numerical simulations are carried out for various mass addition rates and demonstrate that the boundary layer develops when oxygen ion production exceeds approximately 2 x 10(exp 5)/cu m/s. For the upstream solar wind parameters and mass loading rates chosen for these simulations, the results are consistent with observations made on the dayside of Venus for average ionopause conditions near 300 km.

Mcgary, J. E.↗

NASA NTRS2021

HARM3D+NUC: A New Method for Simulating the Post-merger Phase of Binary Neutron Star Mergers with GRMHD, Tabulated EOS, and Neutrino Leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short γ-ray burst. In order to understand the radioactively powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino-leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

Ariadna Murguia-Berthier↗

DOE OSTI2023

Simulations of Precessing Jets and the Formation of X-shaped Radio Galaxies

Abstract Jet precession is sometimes invoked to explain asymmetries in radio galaxy (RG) jets and “X/S/Z-shaped” RGs, caused by the presence of a binary black hole companion to the source active galactic nucleus or by accretion instabilities. We present a series of simulations of RG jet precession to examine how these sources would evolve over time, including a passive distribution of cosmic-ray electrons so we can model radio synchrotron emissions and create synthetic radio maps of the sources. We find that a single source viewed from different angles can result in differing RG morphological classifications, confusing physical implications of these classifications. Additionally, the jet trajectories can become unstable due to their own self-interactions and lead to “reorientation events” that may look like the effects of external dynamics such as shocks, winds, or cold fronts in the medium. Finally, something akin to an “Odd Radio Circle” may be observed in the case of viewing the radio remnant of such a precessing source from a line of sight near the precession axis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2005

3-D GRMHD Simulations of Disk-Jet Coupling and Associated Variabilities and Emission

We have performed the fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3 r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star. We will report initial results of emission calculated based on simulations with Kerr metric using a new ray-tracing method developed by S. Feurst and K. Wu. Images include free-free emission, absorption, and electron scattering. Fluorescent iron line emission and its variabilities from black holes will also be calculated based on simulation results.

Nishikawa, K. I.↗

DOE OSTI2023

Magnetic Reconnection in Solar Flares [Slides]

Magnetic reconnection is the underlying mechanism in solar flares. Magnetic reconnection is an astrophysical process in an electrically conducting plasma where neighboring magnetic field lines running anti-parallel break and reconnection therefore exchanging the direction of their paths. The result of this process is the magnetic field energy converting to kinetic and thermal energy which is explosively released. The scale of a solar is extremely large compared to Earth so kinetic simulations are unable to portray an accurate model. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with the reconnection problem generator file created by Li et al. 2018. To further study this simulation, we began testing a heating/cooling term. The results of the original file will be compared with the modified file with an implemented heating/cooling term to better understand how solar flares function.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Here, we use the public code EBHLIGHT to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1988

An MHD simulation of plasma flow past Io - Alfven and slow mode perturbations

The flow of plasma past Io is studied using a time dependent, three-dimensional magnetohydrodynamic (MHD) simulation. In addition to observing the 'Alfven wing', a standing Alfven wave perturbation expected from analytic theory, it is found that the other MHD modes contribute important perturbations. Standing slow mode perturbations also are present in the flow.

Linker, Jon A.↗

NASA NTRS2019

Predicting the X-Ray Spectra of Stellar-mass Black Holes from Simulations

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a harm3d GRMHD simulation are post-processed by an improved version of the Monte Carlo radiation transport code pandurata (in the corona) and the Feautrier solver PTRANSX (in the disk), with xstar subroutines. The codes are run in a sequential, iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer, including features, such as the Fe Kα emission line and corresponding K-edge absorption trough, due to disk reprocessing of coronal power. For the example cases we consider—a non-rotating 10 M(sub ⊙) black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington—we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1–3 keV), steep power laws extending to high energy (Γ = 2.7–4.5), and prominent, asymmetric Fe Kα emission lines with equivalent widths in the range 40–400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry—no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., is needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

Kinch, Brooks E.↗

DOE OSTI2022

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

Abstract The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm , ipole , and igrmonty , paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Global Magnetohydrodynamic Modeling of the Solar Corona

Under this contract, we have continued our investigations of the large scale structure of the solar corona and inner heliosphere using global magnetohydrodynamic (MHD) simulations. These computations have also formed the basis for studies of coronal mass ejections (CMES) using realistic coronal configurations. We have developed a technique for computing realistic magnetohydrodynamic (MHD) computations of the solar corona and inner heliosphere. To perform computations that can be compared with specific observations, it is necessary to incorporate solar observations into the boundary conditions. We have used the Wilcox Solar Observatory synoptic maps (collected during a solar rotation by daily measurements of the line-of-sight magnetic field at central meridian) to specify the radial magnetic field (B,) at the photosphere. For the initial condition, we use a potential magnetic field consistent with the specified distribution of B, at the lower boundary, and a wind solution consistent with the specified plasma density and temperature at the solar surface. Together this initial condition forms a (non-equilibrium) approximation of the state of the solar corona for the time-dependent MHD computation. The MHD equations are then integrated in time to steady state. Here we describe solutions relevant to a recent solar eclipse, as well as Ulysses observations. We have also developed a model configuration of solar minimum, useful for studying CME initiation and propagation.

Linker, Jon A.↗

NASA NTRS1986

An MHD simulation of the interaction of the solar wind with the outflowing plasma from a comet

The interaction between the solar wind and the outflowing plasmas from a comet has been studied by using a two-dimensional time-dependent magnetohydrodynamic (MHD) simulation. The model reproduced several features of the comet-solar wind interaction predicted by earlier theories and observed on the recent cometary probes. These include the formation of the contact surface and the cometary magnetotail. For a constant interplanetary magnetic field (IMF) the cometary plasma captures field lines which drape over the comet to form an antiparallel magnetic field configuration in the tail and a thin plasma sheet. Eventually, tail magnetic reconnection begins to occur at several points. When the IMF orientation is reversed dayside magnetic reconnection occurs at the subsolar point and a large disturbance propagation down the tail.

Ogino, T.↗

NASA NTRS2004

Theory and Simulation of Real and Ideal Magnetohydrodynamic Turbulence

Incompressible, homogeneous magnetohydrodynamic (MHD) turbulence consists of fluctuating vorticity and magnetic fields, which are represented in terms of their Fourier coefficients. Here, a set of five Fourier spectral transform method numerical simulations of two-dimensional (2-D) MHD turbulence on a 512(sup 2) grid is described. Each simulation is a numerically realized dynamical system consisting of Fourier modes associated with wave vectors k, with integer components, such that k = |k| less than or equal to k(sub max). The simulation set consists of one ideal (non-dissipative) case and four real (dissipative) cases. All five runs had equivalent initial conditions. The dimensions of the dynamical systems associated with these cases are the numbers of independent real and imaginary parts of the Fourier modes. The ideal simulation has a dimension of 366104, while each real simulation has a dimension of 411712. The real runs vary in magnetic Prandtl number P(sub M), with P(sub M) is a member of {0.1, 0.25, 1, 4}. In the results presented here, all runs have been taken to a simulation time of t = 25. Although ideal and real Fourier spectra are quite different at high k, they are similar at low values of k. Their low k behavior indicates the existence of broken symmetry and coherent structure in real MHD turbulence, similar to what exists in ideal MHD turbulence. The value of PM strongly affects the ratio of kinetic to magnetic energy and energy dissipation (which is mostly ohmic). The relevance of these results to 3-D Navier-Stokes and MHD turbulence is discussed.

Shebalin, John V.↗

NASA NTRS1984

Reconnection of Magnetic Fields

Spacecraft observations of steady and nonsteady reconnection at the magnetopause are reviewed. Computer simulations of three-dimensional reconnection in the geomagnetic tail are discussed. Theoretical aspects of the energization of particles in current sheets and of the microprocesses in the diffusion region are presented. Terrella experiments in which magnetospheric reconnection is simulated at both the magnetopause and in the tail are described. The possible role of reconnection in the evolution of solar magnetic fields and solar flares is discussed. A two-dimensional magnetohydrodynamic computer simulation of turbulent reconnection is examined. Results concerning reconnection in Tokamak devices are also presented.

Source record↗

DOE OSTI2025

Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions

Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.

Compton scattering↗

NASA NTRS1998

Large-Scale Dynamics of the Magnetospheric Boundary: Comparisons between Global MHD Simulation Results and ISTP Observations

Understanding the large-scale dynamics of the magnetospheric boundary is an important step towards achieving the ISTP mission's broad objective of assessing the global transport of plasma and energy through the geospace environment. Our approach is based on three-dimensional global magnetohydrodynamic (MHD) simulations of the solar wind-magnetosphere- ionosphere system, and consists of using interplanetary magnetic field (IMF) and plasma parameters measured by solar wind monitors upstream of the bow shock as input to the simulations for predicting the large-scale dynamics of the magnetospheric boundary. The validity of these predictions is tested by comparing local data streams with time series measured by downstream spacecraft crossing the magnetospheric boundary. In this paper, we review results from several case studies which confirm that our MHD model reproduces very well the large-scale motion of the magnetospheric boundary. The first case illustrates the complexity of the magnetic field topology that can occur at the dayside magnetospheric boundary for periods of northward IMF with strong Bx and By components. The second comparison reviewed combines dynamic and topological aspects in an investigation of the evolution of the distant tail at 200 R(sub E) from the Earth.

Berchem, J.↗

NASA NTRS2022

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

supermassive black holes↗

DOE OSTI2022

On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath

Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta's isotropic equivalent mass along the pole at the time of the BH formation is ≲10 -4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta's isotropic equivalent mass along the pole must be small (≲10 -2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ~1.5).

Jets↗