Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HYDRODYNAMIC STABILITY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
DOE OSTI2024

Applications of a Rayleigh-Taylor model to direct-drive laser fusion

Here, this paper presents a simple physics-based model for the interpretation of key metrics in laser direct drive. The only input parameters required are target scale, in-flight aspect ratio, and beam-to-target radius, and the importance of each has been quantified with a tailored set of cryogenic implosion experiments. These analyses lead to compact and accurate predictions of the fusion yield and areal density as a function of hydrodynamic stability, and suggest new ways to take advantage of direct drive. To provide examples, we will discuss how the inferred mix width behaves relative to theory, then show how it could be exploited to perform a direct drive implosion with a Lawson metric or Χ no α of 0.24 ± 0.02—using a novel parameter space at high velocities and beam radii on the OMEGA laser—that projects to ignition at a laser energy ≤2.0 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Performance enhancement of direct-drive shock-augmented ignition inertial fusion implosions through shock timing optimization

Shock-augmented ignition (SAI) [R. H. H. Scott et al., Phys. Rev. Lett. 129, 195001 (2022)] is an alternative inertial confinement fusion concept that is designed to achieve high energy gain by combining improved resilience to instabilities with enhanced fuel compression. In SAI, lower implosion velocities can improve hydrodynamic stability and limited laser intensities reduce the excitation of detrimental laser-plasma instabilities, enabling greater areal density accumulation. Here we report an experimental investigation of SAI using a series of warm D 2 implosions. By timing the augmenting shock to arrive at the implosion center immediately prior to peak compression, the areal density and hot-spot pressure increase by 50% and 112%, respectively, relative to unoptimized designs. Ignition-scale simulations reproduce the observed timing dependence and demonstrate the potential of SAI for inertial fusion energy.

Direct drive↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Correlating Solvation Free Energy to Electrolyte Properties for Lithium Metal Batteries

The electrolyte plays a critical role in lithium metal batteries. In particular, ion solvation profoundly impacts key electrolyte properties and battery performance. Here, in this study, we systematically investigate solvation-property relationships in a series of electrolytes with different solvent-diluent ratios. Through potentiometric techniques that measure the relative solvation free energies of electrolytes, we find that weaker solvation correlates with larger ion clusters, lower ionic conductivity and diffusion coefficient, and superior electrochemical stability. Weaker solvation leads to the formation of a small number of Li clusters with large hydrodynamic radii, which lowers the Li + diffusivity and ionic conductivity of the electrolyte. Concurrently, weaker solvation leads to improved electrochemical stability at both the cathode and anode interfaces. Understanding these solvation-property relationships and trade-offs is important to designing electrolytes for optimized lithium metal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2026

Hydrodynamic Modelling for Tidal Energy Kites: Cooperative Research and Development Final Report, CRADA Number CRD-24-30123

The primary goal of the work is to improve accuracy of hydrodynamic loads acting on the lifting surface of the kite turbine due to complex air-water-lifting body interaction in a low-order numerical model using high-fidelity numerical input. The result of this study will help the Participant better design control algorithms for the kite turbine so that the turbine can harvest more energy while maintaining stability of the platform as well as its reliability. The use of hydrodynamic coefficients/load components from computational fluid dynamic simulation will inform better inputs for a medium-fidelity dynamic model (specifically an OrcaFlex model) which is intended for use in rapid design iterations in the design process.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗
DOE OSTI2025

Stochastic fluctuations and the relaxation time in transient relativistic fluids

We argue that the ratio between the shear viscosity and the shear relaxation time, η/τ π , should be defined as a thermodynamic quantity obtained from the equal-time symmetric correlator of the shear-stress tensor. In kinetic theory, we show that this ratio does not depend on the type of interaction. Similarly, an exact expression for this ratio is obtained for holographic gauge theories. We also determine how stochastic fluctuations change η/τ π in transient relativistic hydrodynamics and show that thermal fluctuations do not spoil causality and stability.

Denicol, Gabriel S. [Universidade Federal Fluminen↗
DOE OSTI2026

Wetted-foam target stability for laser-direct-drive inertial confinement fusion

Using highly resolved three-dimensional radiation-hydrodynamic simulations, the effects of wetted-foam (WF) ablators in direct-drive inertial confinement fusion cryogenic implosions on the OMEGA laser were investigated. The performance of these implosions was compared with that of implosions using targets with solid-plastic ablators. Simulations considering WF ablators as a homogeneous material show no mitigation of the effects of laser imprint, where the absence of this mitigation can be attributed to the use of picket pulses. Heterogeneous random-structure WF ablators introduce large short-wavelength perturbations, which are seeded by WF inhomogeneities. These perturbations exceed the perturbations seeded by the imprint and cause an additional reduction in implosion performance when the average pore size in the WF exceeds ~0.2 μm. Foam-seeded perturbations are expected to be less damaging in more stable high-adiabat (⁠ ≳ 5⁠) implosions. Here, the perturbations introduced by the WF ablators can be substantially reduced by using a foam matrix with a layered structure instead of a random one. Further study of the effects of WFs in ignition designs of the MJ laser energy scale is required in the perspective of future inertial fusion energy mass production.

Igumenshchev, I. V. [University of Rochester, NY (↗
DOE OSTI2024

Single-Step Nonthermal Plasma Synthesis of Water-Soluble and Near-Infrared-Emitting Si Quantum Dots for Bioimaging Applications

Here, we present a single-step nonthermal plasma method for the synthesis of near-infrared (NIR)-emitting and water-soluble Si quantum dots (QDs) for bioimaging applications. Oxygen gas and water vapor were introduced together with acrylic acid (AA) into the afterglow region of the synthesis plasma leading to the surface functionalization of the upstream synthesized Si QDs. The simultaneous surface oxidation and ligand grafting enabled solubility and colloidal stability of the Si QDs in water, as evidenced by strongly reduced hydrodynamic diameters. Aged Si QDs in water emitted NIR photoluminescence (PL) at around 830 nm. The PL quantum yield of the Si QDs in water increased over time from initially undetectable to ~30% after 8 days. Cell viability tests showed that >70% of 3T3 cells survived for 24 h at a concentration of 200 μg/mL of oxidized AA grafted Si QDs. The water solubility, NIR emission with a high quantum yield, and cell viability make the Si QDs promising for bioimaging applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Influence of the drive design on ablation front hydrodynamic instability growth in a capsule implosion at the National Ignition Facility

A critical aspect of inertial confinement fusion experiments lies in the control of ablation front instabilities during the implosion process. The growth of Rayleigh–Taylor and Richtmyer–Meshkov instabilities, seeded by target defects, can significantly degrade the performance of the implosion. Here, this study explores the influence of drive design on ablation front stability using the hydro-growth radiography platform at the National Ignition Facility. This platform allows the measurement of the ablation front hydrodynamic instability growth. Because these measurements are restricted to early convergence, their primary role is to constrain the initial instability growth and benchmark the simulations used to infer shell integrity and performance at peak velocity and ignition. Three ignition designs, Hybrid-E (HyE), SQ-n (“S” for scaling and “Q” for quality), and High temperature, High thickness (HiT) were analyzed. The results demonstrated a high dependence of the growth factor on the choice of design, capsule scale as well as the hohlraum conditions, mainly the radiative temperature and the gold M-band emission. HyE and SQ-n display a similar growth factor range, but their evolution is different due to the design differences in hohlraum conditions. HiT has the lowest growth factor of the studied designs, which is a result of a higher radiative temperature. These results highlight the importance of design choices in controlling instability dynamics. These insights inform future drive design strategies to enhance stability and efficiency in fusion ignition experiments.

Physics↗
DOE OSTI2026

Average hydrodynamic radius analysis reveals critical solvation thresholds in high-concentration lithium electrolytes

Understanding the solvation structures of lithium salts in carbonate- and ether-based electrolytes is central to explaining the exceptional stability of high-concentration electrolytes (HCEs) and localized high-concentration electrolytes (LHCEs). Conventional techniques such as vibrational spectroscopy and one-dimensional NMR provide only limited information, typically restricted to coordination ratios and ion-pair distributions, without revealing the actual size and mass of solution complexes. Here, in this study, we introduce an average hydrodynamic radius (AHR) analysis based on internally referenced DOSY NMR, which enables direct estimation of the average volume and molecular weight of lithium–solvent complexes in solution. Using LiFSI–EMC and LiFSI–EMC–TTE electrolytes, we demonstrate that the onset of effective lithium metal stabilization and aluminum corrosion suppression coincides with the formation of very large complexes, whose average volume exceeds 300 times that of free EMC molecules. This finding supports a new “blocking mechanism,” wherein bulky solvated complexes impede direct solvent access to reactive surfaces. The AHR analysis thus not only clarifies the fundamental origin of HCE and LHCE effectiveness, but also provides a broadly applicable experimental framework for probing complex solvation structures and guiding rational electrolyte design.

Concentrated electrolytes↗
DOE OSTI2026

Self-Assembly Interactions in Magnetite-Coated Cellulose Nanocrystals: Implications for Magnetic Hyperthermia Applications

Magnetic cellulose nanocrystal (MCNC) nanocomposites are promising sustainable and biocompatible platforms for magnetic hyperthermia; however, the molecular mechanisms governing Fe 3 O 4 adsorption and deposition onto CNCs remain poorly understood. Here, sulfated (S-CNC) and TEMPO-oxidized CNCs (T-CNC) were used to prepare nanocomposites at 1:2 and 1:4 CNC:Fe 3 O 4 mass ratios, enabling a systematic evaluation of how surface chemistry and nanoparticle loading dictate interfacial interactions and magneto-colloidal behavior. Bare magnetite nanoparticles were 21 ± 5 nm by TEM but grew to 144 ± 18 in the DLS measurement at pH 7. The S-CNC nanocomposites had hydrodynamic sizes between 144 and 210 nm, not much larger than the 140 nm long CNC rods, suggesting an enhanced dispersion stability compared to Fe3O4 alone. X-ray photoelectron spectroscopy combined with density functional theory revealed that −OH and −COOH groups drive electrostatic adsorption with charge transfer from Fe 3 O 4 to the CNC surface, while T-CNCs showed more favorable adsorption energies and evidence of covalent Fe−O bonding. Vibrating sample magnetometry demonstrated superparamagnetic behavior for all samples, with S-CNC/Fe 3 O 4 1:4 and 1:2 displaying saturation magnetizations of 78 and 77 emu/g-Fe 3 O 4 , close to the 83 emu/g of bare magnetite. The T-CNC composites showed lower (60 and 66 emu/g-Fe3O4) saturation magnetizations. Zero-field-cooled/field-cooled measurements resulted in a blocking temperature of 112 K for all samples, except T-CNC/Fe 3 O 4 1:2 (100 K). Magnetic hyperthermia studies revealed that specific absorption rate (SAR) increased with field strength and Fe 3 O 4 content; however, S-CNC/Fe 3 O 4 (1:2) achieved the highest intrinsic SAR per gram of Fe 3 O 4 (649 W/g-Fe 3 O 4 ) likely due to its anisotropy and fast magnetic relaxation. Cytotoxicity assays confirmed that all nanocomposites were nontoxic toward mammalian cells. These results establish quantitative structure−property relationships between CNC surface chemistry, interfacial bonding mechanisms, and magnetic heating performance, providing a foundation for rational design of biocompatible magnetic nanocomposites for hyperthermia and related applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗