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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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↗
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↗
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↗
DOE OSTI2021

Azimuthal magnetic field distribution in gas-puff Z -pinch implosions with and without external magnetic stabilization

An experimental study of the magnetic field distribution in gas-puff Z pinches with and without a preembedded axial magnetic field (B z0 ) is presented. Spatially resolved, time-gated spectroscopic measurements were made at the Weizmann Institute of Science on a 300 kA, 1.6 μs rise time pulsed-power driver. The radial distribution of the azimuthal magnetic field, B θ , during the implosion, with and without a preembedded axial magnetic field of B z0 = 0.26 T, was measured using Zeeman polarization spectroscopy. The spectroscopic measurements of B θ were consistent with the corresponding values of B θ inferred from current measurements made with a B-dot probe. One-dimensional magnetohydrodynamic simulations, performed with the code TRAC-II, showed agreement with the experimentally measured implosion trajectory, and qualitatively reproduced the experimentally measured radial B θ profiles during the implosion when B z0 = 0.26 T was applied. Simulation results of the radial profile of B θ without a preembedded axial magnetic field did not qualitatively match experimental results due to magneto-Rayleigh-Taylor (MRT) instabilities. Furthermore, our analysis emphasizes the importance of MRT instability mitigation when studying the magnetic field and current distributions in Z pinches. Discrepancies of the simulation results with experiment are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Three-dimensional feedback processes in current-driven metal

Using three-dimensional (3D) magnetohydrodynamic simulations, we study how a pit on a metal surface evolves when driven by intense electrical current density j. Redistribution of j around the pit initiates a feedback loop: j both reacts to and alters the electrical conductivity σ, through Joule heating and hydrodynamic expansion, so that j and σ are constantly in flux. Thus, the pit transforms into larger striation and filament structures predicted by the electrothermal instability theory. Both structures are important in applications of current-driven metal: Here, the striation constitutes a density perturbation that can seed the magneto-Rayleigh-Taylor instability, while the filament provides a more rapid path to plasma formation, through 3D j redistribution. Simulations predict distinctive self-emission patterns, thus allowing for experimental observation and comparison.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Observation of a magneto-Rayleigh-Taylor instability in magnetically collimated plasma jets

We present the direct experimental observation of the formation of a diamagnetic cavity and magneto-Rayleigh-Taylor (MRT) instability in a β ≈ 1 high energy density plasma. Proton radiography is used to measure the two dimensional path-integrated magnetic field in a laser-produced plasma propagating parallel to a preimposed magnetic field. Flutelike structures, associated with the MRT instability, are observed to grow at the surface of the cavity, with a measured wavelength of 1.2 mm and growth time of 4 ns. These measurements are in good agreement with predictions of three dimensional magnetohydrodynamic simulations using the GORGON code. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetized ablative Rayleigh-Taylor instability in three dimensions

Three-dimensional extended-magnetohydrodynamics simulations of the magnetized ablative Rayleigh-Taylor instability are presented. Previous two-dimensional (2D) simulations claiming perturbation suppression by magnetic tension are shown to be misleading, as they do not include the most unstable dimension. For perturbation modes along the applied field direction, the magnetic field simultaneously reduces ablative stabilization and adds magnetic tension stabilization; the stabilizing term is found to dominate for applied fields > 5 T, with both effects increasing in importance at short wavelengths. Additionally, for modes perpendicular to the applied field, magnetic tension does not directly stabilize the perturbation but can result in moderately slower growth due to the perturbation appearing to be 2D (albeit in a different orientation to 2D inertial confinement fusion simulations). In cases where thermal ablative stabilization is dominant the applied field increases the peak bubble-spike height. Resistive diffusion is shown to be important for short wavelengths and long timescales, reducing the effectiveness of tension stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Taylor State Merging at SSX: Experiment and Simulation

We describe experiments and simulations of dynamical merging with two Taylor state plasmas in a Swarthmore Spheromak Experiment (SSX) device. Taylor states are formed by magnetized plasma guns at opposite ends of the device. We performed experiments with Taylor states of both senses of magnetic helicity (right-handed twist or left-handed twist). We present results of both counter-helicity merging (one side left-handed, the other right-handed) and co-helicity merging (both sides left-handed). Experiments show significant ion heating, consistent with magnetic reconnection. We suggest that the merged, warm state could be a suitable target for future magneto-inertial fusion experiments. Magnetohydrodynamic simulations of these experiments reveal the structure of the final relaxed, merged state.

Brown, Michael↗