Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow instability”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurement of the properties of the dust acoustic wave in a magnetic field

We present the first experimental observations of the dust acoustic wave where the wave was observed to propagate in the directions of gravity and magnetic field when these directions were not aligned. The experiments were conducted in the Magnetized Dusty Plasma eXperiment facility using a novel electrode system that allows for the angle between gravity and the magnetic field to be varied in a controlled way. This letter reports on measurements in an rf glow discharge argon plasma environment where the angle between direction of gravity and the magnetic field is 45 $^{\circ }$ . When there was no applied magnetic field, the wave was observed to propagate in the direction of gravity. However, as the magnetic field increased and the ions transitioned from flowing in the direction of gravity to the direction of the magnetic field, a second wave emerged and two distinct waves were observed to simultaneously propagate, one in the direction of gravity and one in the direction of the magnetic field. As the magnetic field was further increased, the wave that propagated in the direction of gravity was suppressed and the wave was only observed to propagate in the direction of the applied magnetic field. We also observe that the speed and the kinetic temperature of the dust for the mode that propagated in the direction of gravity decreased with increasing magnetic field while the speed and the kinetic temperature of the dust for the mode that propagated in the direction of the magnetic field increased with increasing magnetic field. These measurements suggest that an ion-dust streaming instability is at least partially responsible for the high temperatures that have previously been observed in dusty plasmas when the dust acoustic wave is present.

Williams, J. D. (ORCID:000000025199707X)↗
DOE OSTI2024

Wall modes and the transition to bulk convection in rotating Rayleigh-Bénard convection

We investigate states of rapidly rotating Rayleigh-Bénard convection in a cylindrical cell over a range of Rayleigh numbers 3 × 10 5 ≤ Ra ≤ 5 × 10 9 and Ekman numbers 10 − 6 ≤ Ek ≤ 10 − 4 for Prandtl number Pr = 0.8 and aspect ratios 1 / 5 ≤ Γ ≤ 5 using direct numerical simulations. We characterize, for perfectly insulating sidewall boundary conditions, the first transition to convection via wall mode instability and the nonlinear growth and instability of the resulting wall mode states, including a secondary transition to time dependence. We show how the radial structure of the vertical velocity u z and the temperature T is captured well by the linear eigenfunctions of the wall mode instability where the radial width of u z is δ u z ∼ Ek 1 / 3 r / H whereas δ T ∼ e − k r ( k is the wave number of a laterally infinite wall mode state). The disparity in spatial scales for Ek = 10 − 6 means that the heat transport is dominated by the radial structure of u z since T varies slowly over the radial scale δ u z . We further describe how the transition to a state of bulk convection is influenced by the presence of the wall mode states. We use temporal and spatial scales as measures of the local state of convection and the Nusselt number Nu as representative of global transport. Our results elucidate the evolution of the wall state of rotating convection and confirm that wall modes are strongly linked with the boundary zonal flow being the robust remnant of nonlinear wall mode states. We also show how the heat transport ( Nu ) contributions of wall modes and bulk modes are related and discuss approaches to disentangling their relative contributions. Published by the American Physical Society 2024

58 GEOSCIENCES↗
DOE OSTI2025

Image registration for accurate electrode deformation analysis in operando microscopy of battery materials

Operando imaging techniques have become increasingly valuable in both battery research and manufacturing. However, the reliability of these methods can be compromised by instabilities in the imaging setup and operando cells, particularly when utilizing high-resolution imaging systems. The acquired imaging data often include features arising from both undesirable system vibrations and drift, as well as the scientifically relevant deformations occurring in the battery sample during cell operation. For meaningful analysis, it is crucial to distinguish and separately evaluate these two factors. To address these challenges, we employ a suite of advanced image-processing techniques. These include fast Fourier transform analysis in the frequency domain, power spectrum-based assessments for image quality, as well as rigid and non-rigid image-registration methods. These techniques allow us to identify and exclude blurred images, correct for displacements caused by motor vibrations and sample holder drift and, thus, prevent unwanted image artifacts from affecting subsequent analyses and interpretations. Additionally, we apply optical flow analysis to track the dynamic deformation of battery electrode materials during electrochemical cycling. This enables us to observe and quantify the evolving mechanical responses of the electrodes, offering deeper insights into battery degradation. Together, these methods ensure more accurate image analysis and enhance our understanding of the chemomechanical interplay in battery performance and longevity.

Sun, Tianxiao↗
DOE OSTI2024

Ensuring Σ s Y s = 1 in transport of species mass fractions

When transporting species mass fractions in reacting flow simulations, there are physical constraints that must be met. Unfortunately, nonlinear transport schemes such as weighted essentially non-oscillatory (WENO) schemes do not ensure that the sum of mass fractions equals 1. In detonation simulations, errors in the sum of mass fractions are observed to worsen over time when the standard WENO scheme is used. To prevent this, typically one species is forced to absorb all or most of the error in the sum of mass fractions. Here, this work presents an alternative method for correcting the WENO-interpolated mass fractions such that inert species do not change. The method is demonstrated for both argon and nitrogen-diluted hydrogen detonations in one dimension. Using the modified WENO interpolations, the error in the sum is reduced significantly. In addition, the new approach is better able to capture the physical instability expected for nitrogen-diluted detonations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Turbulent Rotating Rayleigh–Bénard Convection

Rotation with thermally induced buoyancy governs many astrophysical and geophysical processes in the atmosphere, ocean, sun, and Earth's liquid-metal outer core. Rotating Rayleigh–Bénard convection (RRBC) is an experimental system that has features of rotation and buoyancy, where a container of height H and temperature difference Δ between its bottom and top is rotated about its vertical axis with angular velocity Ω. The strength of buoyancy is reflected in the Rayleigh number (~ H 3 Δ) and that of the Coriolis force in the Ekman and Rossby numbers (~Ω -1 ). Rotation suppresses the convective onset, introduces instabilities, changes the velocity boundary layers, modifies the shape of thermal structures from plumes to vortical columns, affects the large-scale circulation, and can decrease or enhance global heat transport depending on buoyant and Coriolis forcing. RRBC is an extremely rich system, with features directly comparable to geophysical and astrophysical phenomena. Here we review RRBC studies, suggest a unifying heat transport scaling approach for the transition between rotation-dominated and buoyancy-dominated regimes in RRBC, and discuss non-Oberbeck–Boussinesq and centrifugal effects.

42 ENGINEERING↗
DOE OSTI2022

Investigations into EGR dilution tolerance in a pre-chamber ignited GDI engine

EGR is a useful means to improve fuel economy in spark-ignited gasoline engines, especially under part-load conditions by means of reducing throttling losses. There is also a significant reduction in NOx emissions due to the reduced peak cylinder temperatures. However, the ability of conventional spark ignition systems to reliably ignite dilute mixtures limits the dilution tolerance due to the onset of combustion instability. In this study, an active pre-chamber system was evaluated at part load condition with various EGR rates using a regular grade full boiling range gasoline. A premixed air-fuel mixture was supplied to the pre-chamber, which assists in scavenging the residual mass fractions to enable reliable ignition of EGR dilute mixtures in the main combustion chamber. A range of pre-chamber injection, air-fuel mixture formulation, and ignition timing strategies were evaluated during the EGR sweeps. In addition, 0-D and 1-D simulations were conducted in GT-Power to assess the thermodynamic state and composition in the pre-chamber, as well as estimate the jet momentum at the various dilution rates. The active pre-chamber extended the EGR dilution limit from 20% for conventional SI to above 30%, with a higher pre-chamber flow rate (more scavenging) resulting in increased combustion stability. A tradeoff between the jet momentum and spark timing was observed when analyzing the GT-Power simulations and the pre-chamber jet momentum was found to decrease with increasing EGR rate. At unstable high EGR operation, no relation could be found between the pre-chamber indicated data and main chamber combustion stability.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Ensuring Transient Stability With Guaranteed Region of Attraction in DC Microgrids

DC microgrids have promising applications in renewable integration due to their better energy efficiency when connecting DC components. However, they might be unstable since many loads in a DC microgrid are regulated as constant power loads (CPLs) that have a destabilizing negative impedance effect. As a result, the state trajectory displacement caused by abrupt load changes or contingencies can easily lead to instability. Many existing works have been devoted to studying the region of attraction (ROA) of a DC microgrid, in which the system is guaranteed to be asymptotically stable. Nevertheless, existing work either focuses on using numerical methods for ROA approximations that generally have no performance guarantees or cannot ensure a desired ROA for a general DC microgrid. To close this gap, this paper develops an innovative control synthesis algorithm to make a general DC microgrid have a theoretically guaranteed ROA, for example, to cover the entirety of its operating range regarding state trajectories. Here, we first study the nonlinear dynamics of a DC microgrid to derive a novel transient stability condition to rigorously certify whether a given operating range is a subset of the ROA; then, we formulate a control synthesis optimization problem to guarantee the condition’s satisfaction. This condition is a linear constraint, and the optimization problem resembles an optimal power flow problem and has a good computational behavior. Simulation case studies verify the validity of the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2022

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Empirical probability and machine learning analysis of m , n = 2, 1 tearing mode onset parameter dependence in DIII-D H-mode scenarios

m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: I. Computational approach

Theoretical computations reveal insight into parameters that govern the accelerated crucible rotation technique (ACRT) applied to the gradient freeze (VGF) growth of cadmium zinc telluride (CZT). A metric based on the classic Mullins and Sekerka instability criterion is put forth that provides a quantitative means of assessing the impact of different ACRT rotation schedules on morphological stability and inclusion formation during CZT growth. Here this metric is employed in a 2-k factorial design of experiments that identifies the maximum rotation rate as the most important ACRT parameter, followed by the acceleration period and the rest period. Subsequent optimization of rotation schedule is performed for a specific VGF growth system, and mechanistic explanations are provided for the preferred ACRT schedule, which improves interface stability via slower rotation rates and longer acceleration periods than would be applied using classical ACRT schedules.

36 MATERIALS SCIENCE↗
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↗