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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Direct prediction of nonlinear tearing mode saturation using a variational principle

It is shown that the variational principle of multi-region relaxed magnetohydrodynamics (MRxMHD) can be used to predict the stability and nonlinear saturation of tearing modes in strong guide field configurations without resolving the dynamics and without explicit dependence on the plasma resistivity. While the magnetic helicity is not a good invariant for tearing modes, we show that the saturated tearing mode can be obtained as an MRxMHD state of a priori unknown helicity by appropriately constraining the current profile. The predicted saturated island width in a tearing-unstable force-free slab equilibrium is shown to reproduce the theoretical scaling at small values of Δ' and the scaling obtained from resistive magnetohydrodynamics simulations at large Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hall instability driven seeding of helical magneto-Rayleigh–Taylor instabilities in axially premagnetized thin-foil liner Z-pinch implosions

Helical magneto-Rayleigh–Taylor instability (MRTI) structures have been observed in z-pinch-driven liner implosion experiments with a pre-imposed axial magnetic field. We show that the formation of these helical structures can be described by a Hall magnetohydrodynamical (HMHD) model. We used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) [Seyler and Martin, Phys. Plasmas 18, 012703 (2011)] to study these helical instabilities and show that a Hall interchange instability in low-density coronal plasma immediately surrounding the dense liner is responsible for producing helically oriented effects in the magnetic field and current density within the coronal layer. Here, this seeds the helical pitch angle of the MRTI even when other proposed helical seeding mechanisms are either not present in the experiments or not accounted for in the simulations. For example, this mechanism does not require low-density power-feed plasmas to be swept in from large radius or the development of electrothermal instabilities. The Hall Instability is, thus, a new, independent explanation for the origin of the helical instabilities observed in axially premagnetized liner experiments. Simulation results supporting this mechanism are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Astronomy & Astrophysics↗
DOE OSTI2024

Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius, R B = 10 3 R g . The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outside R B and suppress BH accretion to 1.5% of the Bondi rate, M ̇ B . Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of M ̇ B reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% M ̇ B . Because the BAD and RAD states tangle up the jets and destroy them well inside R B , they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

Astronomy & Astrophysics↗
DOE OSTI2023

MHD simulation on magnetic compression of field reversed configurations with NIMROD

Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Time-resolved velocity and ion sound speed measurements from simultaneous bow shock imaging and inductive probe measurements

We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (“b-dot”) probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large (R M > 10), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach–Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the sound speed to infer the value of Z̄T e , where Z̄ is the average ionization and T e is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic (M S ~ 8), super-Alfvénic (M A ~ 2) aluminum plasma generated during the ablation stage of an exploding wire array on the Magpie generator (1.4 MA, 250 ns). The velocity and Z̄T e measurements agree well with the optical Thompson scattering measurements reported in the literature and with 3D resistive magnetohydrodynamic simulations in GORGON.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Multi-fidelity equations of state and transport coefficient datasets for pulsed-power applications

Reliably simulating experiments relevant to the National Nuclear Security Administration (NNSA) requires a detailed description of material properties across a wide range of conditions. Such properties include the equations of state, charged-particle transport coefficients, and optical properties like the opacity. Together, these properties make up the material models used in radiation-magnetohydrodynamic simulations of nuclear fusion experiments. Many of these models do not incorporate uncertainties in the data used to produce them. It is unknown whether these uncertainties significantly impact the interpretation of simulation results and diagnostics. The purpose of this work is to quantify how such uncertainties impact simulations of pulsed-power experiments. We accomplished this task by first assessing discrepancies between approaches used to generate the data. This included bringing together members of the high-energy-density community spanning the three NNSA laboratories and multiple universities. Then, using these data, we developed a general framework that systematically incorporates physical uncertainties within the material models suitable for uncertainty quantification analyses. The framework utilizes machine learning, Bayesian inference, and incorporates multi-fidelity datasets. We demonstrated the framework by quantifying the impact that material model uncertainties have on simulations of pulsed-power experiments underway on Z at Sandia National Laboratories. As a result of this work, we discovered that modest uncertainties in material models (roughly 20%) correspond to significant uncertainties in the outputs from simulations. Our framework has enabled rapid construction of material models through an automated procedure and allows for the generation of material models of interest to the NNSA.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Kelvin Helmholtz Instability at the Equatorial Magnetotail Boundary: MHD Simulation and Comparison with Geotail Observations

On March 24, 1995 the Geotail spacecraft observed large fluctuations of the magnetic field and plasma properties in the Low Latitude Boundary Layer (LLBL) about 15 R(sub E) tailward of the dusk meridian. Although the magnetospheric and the magnetosheath field were strongly northward, the B(sub z) component showed strong short duration fluctuations in which B(sub z) could even reach negative values. We have used two-dimensional magnetohydrodynamic simulations with magnetospheric and magnetosheath input parameters specifically chosen for this. Geotail event to identify the processes which cause the observed boundary properties. It is shown that these fluctuations can be explained by the Kelvin-Helmholtz instability if the k vector of the instability has a component along the magnetic field direction. The simulation results show many of the characteristic properties of the Geotail observations. In particular, the quasi-periodic strong fluctuations are well explained by satellite crossings through the Kelvin-Helmholtz vortices. It is illustrated how the interior structure of the Kelvin-Helmholtz vortices leads to the rapid fluctuations in the Geotail observations. Our results suggest an average Kelvin-Helmholtz wavelength of about 5 R(sub E) with a vortex size of close to 2 R(sub E) for an average repetition time of 2.5 minutes. The growth time for these waves implies a source region of about 10 to 16 R(sub E) upstream from the location of the Geotail spacecraft (i.e., near the dusk meridian). The results also indicate a considerable mass transport of magnetosheath material into the magnetosphere by magnetic reconnection in the Kelvin-Helmholtz vortices.

Fairfield, Donald H.↗
DOE OSTI2022

Disc tearing leads to low and high frequency quasi-periodic oscillations in a GRMHD simulation of a thin accretion disc

Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broad-band variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which – if properly understood – could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work, we analyse a high-resolution 3D general relativistic magnetohydrodynamic simulation of a geometrically thin accretion disc, which is tilted by 65° with respect to the black hole spin axis. We find that the Lense–Thirring torque from the rapidly spinning 10 M ⊙ black hole causes several sub-discs to tear off within ~10–20 gravitational radii. Tearing occurs in cycles on time-scales of seconds. During each tearing cycle, the inner sub-disc precesses for 1–5 periods before it falls into the black hole. We find a precession frequency of ~3 Hz, consistent with observed low-frequency QPOs. In addition, we find a high frequency QPO (HFQPO) with centroid frequency of ~55 Hz in the power spectra of the mass-weighted radius of the inner disc and the radial mass flux. This signal is caused by radial epicyclic oscillations of a dense ring of gas at the tearing radius, which suggests a corresponding modulation of the X-ray light curve and may thus explain some of the observed HFQPOs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Do chaotic field lines cause fast reconnection in coronal loops?

Over the past decade, Boozer has argued that three-dimensional (3D) magnetic reconnection fundamentally differs from two-dimensional reconnection due to the fact that the separation between any pair of neighboring field lines almost always increases exponentially over distance in a 3D magnetic field. According to Boozer, this feature makes 3D field-line mapping chaotic and exponentially sensitive to small non-ideal effects; consequently, 3D reconnection can occur without intense current sheets. We test Boozer's theory via ideal and resistive reduced magnetohydrodynamic simulations of the Boozer–Elder coronal loop model driven by sub-Alfvénic footpoint motions [A. H. Boozer and T. Elder, Phys. Plasmas 28, 062303 (2021)]. In this work, our simulation results significantly differ from their predictions. The ideal simulation shows that Boozer and Elder under-predict the intensity of current density due to missing terms in their reduced model equations. Furthermore, resistive simulations of varying Lundquist numbers show that the maximal current density scales linearly rather than logarithmically with the Lundquist number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulated annealing of reduced magnetohydrodynamic systems

Theory of simulated annealing (SA), a method for equilibrium and stability analyses for Hamiltonian systems, is reviewed. The SA explained in this review is based on a double bracket formulation that derives from Hamiltonian structure. In addition to general theoretical aspects, the explicit formulation as well as numerical applications are presented. Both finite and infinite degree-of-freedom systems are treated, in particular, the heavy top, a toy model mimicking low-beta reduced magnetohydrodynamics (MHD) and low- and high-beta reduced MHD. Furthermore, the numerical results successfully demonstrate the usefulness of SA for equilibrium and stability analyses. At the same time, the results raise some future issues that are discussed in the paper.

Poisson Bracket↗
DOE OSTI2025

Simulated thermonuclear performance of auto-magnetizing helical liner MagLIF target implosions

Magnetized liner inertial fusion (MagLIF) implosions on the Z accelerator require magnetization of the fuel to reduce thermal conduction losses of energy from the laser-preheated deuterium fusion fuel to the surrounding liner material. While external field coils traditionally used to axially magnetize MagLIF targets are limited to 10–20 T (or perhaps ∼ 30 T with technological development), calculations suggest that higher axial magnetic fields would improve thermal insulation of the fuel and improve MagLIF target performance. We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing auto-magnetizing helical liners composed of discrete metallic helical conduction paths separated by electrically insulating material. These simulated auto-magnetizing (AutoMag) liners produce axial magnetic fields >30 T inside the fusion fuel prior to implosion. Simulations indicate that higher thermonuclear yields and burn-averaged fuel ion temperatures are attained for implosions using AutoMag liners compared to standard MagLIF implosions. Comparable implosion morphology is evident in synthetic x-ray images of AutoMag–MagLIF and standard MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗
NASA NTRS1993

Current loop coalescence studied by 3-D electromagnetic particle code

Solar flare plasma data from the Yohkoh satellite is analyzed. The interactions of current loops were observed in the active regions on the Sun. This observation pointed out the importance of the idea that the solar flare is generated by the coalescence of current loops. The three dimensional electromagnetic particle simulations are to help in understanding the global interaction between two current loops including the evolution of the twist of loops due to instabilities. Associated rapid dynamics of current loop coalescence such as reconnection, shock waves and associated kinetic processes such as energy transfer, acceleration of particles, and electromagnetic emissions are to be studied by the code to complement analytical theories and magnetohydrodynamic simulations of the current loop coalescence. The simulation results show the strong interactions between two current loops, beam and whistler instabilities, and associated parallel and perpendicular particle heating.

Nishikawa, Ken-Ichi↗
DOE OSTI2025

Numerical investigation of magneto-inertial fusion targets magnetized by dynamic enforcement of helical current flow

Magnetized Liner Inertial Fusion (MagLIF) targets require premagnetization to reduce thermal conduction losses from the laser-heated fuel to the liner material during implosion. MagLIF targets are typically magnetized by external magnetic field coils which are technologically limited to providing <20 T axial magnetic fields in the fusion fuel. We present a novel target design, AutoMag-D, which employs dynamic enforcement of helical current in the liner, resulting in >20 T axial magnetic field in the fuel region prior to implosion without the need of external magnetic field coils. These self-magnetizing liners are made of helically oriented alternating electrically conductive material and electrically insulating material surrounded by a thin, conductive outer radial layer. As the liner is pulsed with a ∼ 20 MA, ∼ 100 ns rise time drive current from Z, the outer layer of conductive material is shocked and intensely Joule heated, causing it to melt, vaporize, and turn to plasma. This allows the magnetic drive field to diffuse radially inward which dynamically enforces current flow in the helical conduction paths in the liner and produces axial magnetic field in the fuel region prior to implosion of the inner liner surface. AutoMag-D liner designs do not require dielectric breakdown of electrically insulating material (as in traditional auto-magnetizing helical liners) and do not require helical return current geometries (as in dynamic screw pinches). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing AutoMag-D liner designs. Simulated AutoMag-D targets demonstrate improved fuel conditions and thermonuclear yield compared to traditional MagLIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Discriminating Accretion States via Rotational Symmetry in Simulated Polarimetric Images of M87

In 2017 April, the Event Horizon Telescope observed the shadow of the supermassive black hole at the core of the elliptical galaxy Messier 87. While the original image was constructed from measurements of the total intensity, full polarimetric data were also collected, and linear polarimetric images are expected in the near future. We propose a modal image decomposition of the linear polarization field into basis functions with varying azimuthal dependence of the electric vector position angle. We apply this decomposition to images of ray traced general relativistic magnetohydrodynamics simulations of the Messier 87 accretion disk. Furthermore, for simulated images that are physically consistent with previous observations, the magnitude of the coefficient associated with rotational symmetry, β2, is a useful discriminator between accretion states. We find that at 20 μas resolution, $| {\beta }_{2}| $ is greater than 0.2 only for models of disks with horizon-scale magnetic pressures large enough to disrupt steady accretion. We also find that images with a more radially directed electric vector position angle correspond to models with higher black hole spin. Our analysis demonstrates the utility of the proposed decomposition as a diagnostic framework to improve constraints on theoretical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Use of hydrodynamic theory to estimate electrical current redistribution in metals

Using the analogy between hydrodynamic and electrical current flow, we study herein how electrical current density j redistributes and amplifies due to two commonly encountered inhomogeneities in metals. First, we consider flow around a spherical resistive inclusion and find significant j amplification, independent of inclusion size. Hence, even μ m-scale inclusions can affect performance in applications by creating localized regions of enhanced Joule heating. Next, we investigate j redistribution due to surface roughness, idealized as a sinusoidal perturbation with amplitude A and wavelength λ . Theory predicts that j amplification is determined by the ratio A / λ , so that even “smooth” surface finishes (i.e., small A ) can generate significant amplification, if λ is correspondingly small. We compare theory with magnetohydrodynamic simulation to illustrate both the utility and limitations of the steady-state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗