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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Observational Signatures of Coronal Heating in Magnetohydrodynamic Simulations without Radiation or a Lower Atmosphere

It is extremely difficult to simulate the details of coronal heating and also make meaningful predictions of the emitted radiation. Thus, testing realistic models with observations is a major challenge. Observational signatures of coronal heating depend crucially on radiation, thermal conduction, and the exchange of mass and energy with the transition region and chromosphere below. Many magnetohydrodynamic simulation studies do not include these effects, opting instead to devote computational resources to the magnetic aspects of the problem. We have developed a simple method of accounting approximately for the missing effects. It is applied to the simulation output ex post facto and therefore may be a valuable tool for many studies. We have used it to predict the emission from a model corona that is driven by vortical boundary motions meant to represent photospheric convection. We find that individual magnetic strands experience short-term brightenings, both scattered throughout the computational volume and in localized clusters. The former may explain the diffuse component of the observed corona, while the latter may explain bright coronal loops. Several observed properties of loops are reproduced reasonably well: width, lifetime, and quasi-circular cross section (aspect ratio not high). Our results lend support to the idea that loops are multistranded structures heated by "storms" of nanoflares.

James A. Klimchuk↗
NASA NTRS2005

General Relativistic Magnetohydrodynamic Simulations of Collapsars

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. Initially, we assume that the core collapse has failed in this star. A rotating black hole of a few solar masses is inserted by hand into the calculation. The simulation results show the formation of a disklike structure and the generation of a jetlike outflow near the central black hole. The jetlike outflow propagates and accelerated mainly by the magnetic field. The total jet velocity is approximately 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is directly converted to kinetic energy of the jet rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster.

Mizuno, Yosuke↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2024

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗
NASA NTRS2005

A General Relativistic Magnetohydrodynamic Simulation of Jet Formation

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation ofjet 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 (velocity approx.0.3c) is created, as shown by previous two-dimensional axi- symmetric simulations with mirror symmetry at the equator. The three-dimensional simulation ran over 100 light crossing time units (T(sub s) = r(sub s)/c, where r(sub s = 2GM/c(sup 2), which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted owing in part to magnetic pressure from the twisting of the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r(sub s). At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface ofthe 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 outward with a wider angle than the initial jet. The widening of the jet is consistent with the outward-moving torsional Alfven waves. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk because of the iack of streaming materiai from an accompanying star.

Nishikawa, K.-I.↗
DOE OSTI2025

Assessment of dynamic-screw-pinch-driven, current-scaled MagLIF target implosion performance using 3D magnetohydrodynamic simulations

Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

A General Relativistic Magnetohydrodynamics Simulation of Jet Formation with a State Transition

We have performed the first 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 (velocity sim 0.3c) 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. 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 shock wave. This evolution of jet-disk coupling suggests that the low/hard state of the jet system may switch to the high/soft state with a wind, as the accretion rate diminishes.

Nishikawa, K. I.↗
NASA NTRS2016

Disk Emission from Magnetohydrodynamic Simulations of Spinning Black Holes

We present the results of a new series of global, three-dimensional, relativistic magnetohydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of H/R approx. 0.05 and spin parameters of a/M = 0, 0.5, 0.9, and 0.99. Using the ray-tracing code Pandurata, we generate broadband thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Finally, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well known to the continuum-fitting method of measuring black hole spin.

accretion↗
DOE OSTI2022

Long-term general relativistic magnetohydrodynamics simulations of magnetic field in isolated neutron stars

ABSTRACT Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-C 1 code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-C 1 . A C 1 coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-C 1 extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of a stellarator plasma with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

Global magnetohydrodynamic simulation of the two-dimensional magnetosphere

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

Leboeuf, J. N.↗
NASA NTRS2006

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗
NASA NTRS1991

Magnetohydrodynamic simulation of the radial evolution and stream structure of solar-wind turbulence

A unified interpretation of observations of interplanetary fluctuations is presented in terms of nearly incompressible magnetohydrodynamics. Incompressive effects explain the rapid evolution of turbulence in slow wind containing the heliospheric current sheet. The relative constancy of the spectrum of 'inward propagating' fluctuations compared to the rapid decline in 'outward' fluctuations results from incompressive spectral transfer combined with strong dissipation of the outward fluctuations. Secondary compressive effects account for nearly pressure-balanced structures and the density fluctuation levels.

Roberts, D. A.↗
DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗
NASA NTRS1996

Dynamical Evolution of a Coronal Streamer-Flux Rope System: 2. A Self-Consistent Non-Planar Magnetohydrodynamic Simulation

The dynamical response of a helmet streamer to a flux rope escape from the sub-photosphere is examined in a physically self-consistent manner within the approximation of axisymmetric three-dimensional magnetohydrodynamics (i.e., so-called '2 1/2 D'). In contrast to the previous planar analyses of Paper 1 (Wu, Guo, and Wang), the present study shows, with the inclusion of out-of-plane components of magnetic and velocity fields, that the magnetic configuration represents a helical flux rope instead of a planar bubble as shown in Paper 1. Because of this more physically-realistic configuration, we are able to examine the dynamical evolution of the helical flux rope's interaction with the helmet streamer. This process leads to the formation of two parts of the solar mass ejection: (i) the expulsion of the helmet dome due to eruption of this flux rope, and (ii) the flux rope's eruption itself. When this two-part feature propagates out to the interplanetary space, it exhibits all the physical characteristics of observed interplanetary magnetic clouds. These numerical simulations also show that the dynamical behavior of the streamer-flux rope system has three distinct states: (i) quasi-equilibrium, (ii) non-equilibrium, and (iii) eruptive state depending on the energy level of the flux rope.

Wu, S. T.↗
DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗