Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma instabilities”

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.

137 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2025

Microstructure, mechanical, and thermal properties of compositionally complex (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 ceramics

Novel compositionally complex borides, (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 , were fabricated using spark plasma sintering process. (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 exhibits a dual-phase microstructure, in which (Hf,Zr,Nb,Ti)B 2 is a primary phase with the hexagonal structure and LaB 6 is a secondary phase with a cubic structure. The mechanical properties of both (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 are comparable, with a combination of high hardness and moderate fracture toughness. Thermal diffusivity and conductivity of (Hf,Zr,Nb,Ti)B 2 are much lower than the individual transition metal borides but are significantly increased by the addition of LaB 6 . Herein, it is implied that the thermal properties of boride ceramics can be controlled through the appropriate design of principal metal element compositions.

borides↗
DOE OSTI2023

Fracture‐Resistant and Thermally Insulating Ultrahigh‐Temperature Carbide Foams

Dense ultrahigh‐temperature ceramics (UHTCs) carbides are recognized as potential materials for thermal protection systems (TPS) owing to properties beyond existing structural materials’ capabilities. Recent advances in UHTCs have enabled the development of multiscale porous microstructures. Herein, it is highlighted that the porosity in UHTCs are no longer treated as a defect but as a functional property specifically tailored for thermal insulation. It is a promising solution to design and fabricate bulk UHTC foams via a freeze‐drying (FD) approach followed by calibrated pressureless spark plasma sintering. Herein, monolithic TaC and HfC UHTC foams and their composite show the partial solid–solution formation of (Ta, Hf)C with porosity ≥50%. TaC–HfC foam (≈80–92 N) shows an intermediate load‐bearing capability compared to monolithic TaC (≈120–135 N) and HfC (≈28–35 N) foams, with no evident cracking on the sample surface. The thermal conductivity of partial solid‐solution TaC–HfC foam increases up to fivefold compared to parent UHTC foams. In the results, solid solutions’ efficacy and pores’ unidirectionality in providing thermal insulation to TaC–HfC while maintaining its high‐load bearing capability are illustrated. In conclusion, the developed technique establishes a new paradigm shift in UHTCs, expanding their potential for TPS in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Probing and Tuning Spatiotemporal pH Evolution in Aqueous Zinc Ion Batteries

Aqueous zinc-ion batteries (AZIBs) offer a combination of safety and low cost, but they suffer from dynamic pH fluctuations that drive parasitic reactions and capacity fading. In this review, we report the current mechanistic understanding of the dynamic pH evolution across the electrolyte-electrode interface and its role in triggering interphase instability. We summarize advanced strategies to regulate pH, including electrolyte engineering, electrode surface modification, and interphase construction. We further discuss current methods for probing pH dynamics and outline potential future approaches for gaining deeper mechanistic and design insights.

Zheng, Xueli [SLAC]↗
DOE OSTI2026

Evaluating Fracture Behavior of Bioinspired Alumina-YSZ Composites through Static and Dynamic Mechanical Testing

Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.

Bioinspired materials↗
DOE OSTI2024

Making Plasticized Polymer Electrolytes Stable Against Sodium Metal for High‐Energy Solid‐State Sodium Batteries

Solid polymer electrolytes based on plastic crystals are promising for solid-state sodium metal (Na 0 ) batteries, yet their practicality has been hindered by the notorious Na 0 -electrolyte interface instability issue, the underlying cause of which remains poorly understood. Here, in this study, by leveraging a model plasticized polymer electrolyte based on conventional succinonitrile plastic crystals, we uncover its failure origin in Na 0 batteries is associated with the formation of a thick and non-uniform solid electrolyte interphase (SEI) and whiskery Na 0 nucleation/growth. Furthermore, we design a new additive-embedded plasticized polymer electrolyte to manipulate the Na 0 deposition and SEI formulation. For the first time, we demonstrate that introducing fluoroethylene carbonate (FEC) additive into the succinonitrile-plasticized polymer electrolyte can effectively protect Na 0 against interfacial corrosion by facilitating the growth of dome-like Na 0 with thin, amorphous, and fluorine-rich SEIs, thus enabling significantly improved performances of Na//Na symmetric cells (1,800 h at 0.5 mA cm −2 ) and Na//Na 3 V 2 (PO 4 ) 3 full cells (93.0 % capacity retention after 1,200 cycles at 1 C rate in coin cells and 93.1 % capacity retention after 250 cycles at C/3 in pouch cells at room temperature). Our work provides valuable insights into the interfacial failure of plasticized polymer electrolytes and offers a promising solution to resolving the interfacial instability issue.

25 ENERGY STORAGE↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2025

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Open quantum system approach to inclusive jet production in heavy-ion collisions

We derive a factorization formula for inclusive jet production in heavy-ion collisions using the tools of Effective Field Theory (EFT). We show how physics at widely separated scales in this process can be systematically separated by matching to EFTs at successively lower virtualities. Owing to a strong scale separation, we recover a vacuum-like DGLAP evolution above the jet scale, while the additional low-energy scales induced by the medium effectively probe the internal structure of the jet. As a result, the cross section can be written as a series with an increasing number of subjets characterized by perturbative matching coefficients each of which is convolved with a distinct function. These functions encode broadening, medium-induced radiations as well as quantum interference such as the Landau-Pomeranchuk-Migdal effect and color coherence dynamics to all orders in perturbation theory. As a first application of this EFT framework, we investigate the case of an unresolved jet and show how the cross section can be factorized and fully separate the jet dynamics from the universal physics of the medium. To compare to the existing literature, we explicitly compute the medium jet function at next-to-leading order in the coupling and leading order in medium opacity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Robustness of inertial fusion energy relevant implosions to low-mode asymmetries

Inertial confinement fusion implosions for power production need to be robust. Degradation sources need to be understood and considered for an economic power plant. In this work, the impact of low-mode asymmetries on yield is studied using 2D radiation hydrodynamics simulations for inertial fusion energy relevant implosions. First, a high-gain indirectly driven capsule design is developed. This design is subsequently degraded with low-mode drive asymmetries. The results show that yield is constant with increasing amounts of asymmetry until a cliff is reached where the yield quickly drops below ignition. Robustness arises because the degraded hotspot still has enough time to get to ignition before capsule decompression. This robustness occurs even though the asymmetry couples less energy to the hotspot and causes decompression to occur sooner. The robustness is linked to the generalized Lawson parameter, χ⁠. The value of χ steadily decreases with increasing asymmetry. Once χ is close to one, the hotspot is no longer able to ignite before the capsule decompresses. A simple piston model is used to predict where the robustness cliff is for any design. The trade-off between yield and robustness is also demonstrated.

Physics - Plasma physics↗
DOE OSTI2026

Interface-induced fast Li+ transport in mixed ionic–electronic conductors

Interfacial instability between lithium metal and solid-state electrolytes limits the performance of all-solid-state lithium metal batteries (ASSLBs), leading to parasitic reactions, non-uniform Li+ flux, and dendrite growth. Here, we develop a composite interlayer composed of the anti-perovskite Li2OHCl0.75Br0.25 (AP) and carbon nanotubes (CNTs) to enhance both interfacial stability and ionic transport. The AP–CNT interlayer exhibits enhanced Li+ conductivity arising from interfacial electron transfer from AP to CNTs, which generates a built-in electric field that facilitates Li+ migration. Lithium symmetric cells incorporating this interlayer achieve a high critical current density of 2.4 mA cm−2 at 55 °C. This design integrates chemical robustness with coupled ion–electron transport, offering a generalizable strategy for safe, dendrite-free, and high-performance ASSLBs.

He, Chenche↗
DOE OSTI2026

A Chemist’s Guide to Solid-State Qubits: Diamond Color Centers as Embedded Vacancy-Centered Molecules

Diamond point defects or color centers are at the crux of solid-state quantum technology applications, from computing to sensing. Modeling molecular qubits has advanced by leaps and bounds mainly due to the advanced quantum mechanical (QM) tools in the arsenal of quantum chemists for modeling (isolated) molecular systems, which should be easily translatable for analogous systems in the solid state. Here, we simulate a series of diamond defects: a carbon vacancy (V C ) and second row substitutions that are adjacent to a V C or D C V C , with D = B, N, and O, to evaluate their magnetic and optical properties. We applied molecular orbital principles, an embedded cluster method, and a multiconfigurational QM theory to computationally characterize the above-mentioned defects for a multitude of spin, charge, and excited states. Our computational approach delivered high-quality QM insights that compare well with experiments. Furthermore, with this framework, we confirm the spin-triplets B C V C – and O C V C 0 to be analogous to N C V C – , with calculated bright transitions that respectively correspond to UV (3.80 eV) and blue (2.65 eV) light absorption and emission lifetimes (τ) of 2.3 and 5.6 ns, which are higher in energy and more short-lived (and thus brighter) than that for N C V C – : green (2.26 eV), τ = 9.5 ns.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2024

Prediction of BiS2-type pnictogen dichalcogenide monolayers for optoelectronics

Abstract In this work, we introduce a 2D materials family with chemical formula MX 2 (M={As, Sb, Bi} and X={S, Se, Te}) having a rectangular 2D lattice. This materials family has been predicted by systematic ab-initio structure search calculations in two dimensions. Using density-functional theory and many-body perturbation theory, we study the structural, vibrational, electronic, optical, and excitonic properties of the predicted MX 2 family. Our calculations reveal that the predicted SbX 2 and BiX 2 monolayers are stable while the AsX 2 layers exhibit an in-plane ferroelectric instability. All materials display strong excitonic effects and good optical absorption within the infrared-to-visible range. Hence, these monolayers can harvest solar energy and serve in optoelectronics applications. Furthermore, our results indicate that exfoliation of the predicted MX 2 monolayers from their bulk counterparts is experimentally viable.

Materials Science↗
DOE OSTI2025

Divergent evolution of slip banding in CrCoNi alloys

Abstract Metallic materials under high stress often exhibit deformation localization, manifesting as slip banding. Over seven decades ago, Frank and Read introduced the well-known model of dislocation multiplication at a source, explaining slip band formation. Here, we reveal two distinct types of slip bands (confined and extended) in compressed CrCoNi alloys through multi-scale testing and modeling from microscopic to atomic scales. The confined slip band, characterized by a thin glide zone, arises from the conventional process of repetitive full dislocation emissions at Frank–Read source. Contrary to the classical model, the extended band stems from slip-induced deactivation of dislocation sources, followed by consequent generation of new sources on adjacent planes, leading to rapid band thickening. Our findings provide insights into atomic-scale collective dislocation motion and microscopic deformation instability in advanced structural materials.

Science & Technology - Other Topics↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multi-modal dynamic radiography using short-pulse laser-generated probe beams

Radiography is an important tool for the interrogation of dynamic experiments in the fields of dynamic properties of materials, and in condensed matter, high explosive, and high-energy-density physics. Multi-modal radiography advances the hypothesis that combining the information delivered by multiple radiographic modalities can lead to more constrained (improved) “reconstruction” of the scene than can be obtained from a single probe. We identify four modalities: multi-probe, time sequence, multi-view, and multi-messenger. Multi-probe radiography is a promising candidate for a next-generation dynamic radiographic facility. High-energy X-rays are the most frequently used probe for dynamic radiography, although recent developments show the utility of proton (pRad), electron (eRad), and neutron probe beams. Because each probing species interacts with material in the radiographic scene through quantitatively different mechanisms, each returns independent information about the scene, which can add extra constraints to the reconstruction process. How to conduct detailed, quantitative “co-analysis” of multiple data streams remains an area of active research. Multi-beam, short-pulse, laser-generated probes offer sufficient dose, an appropriate spectrum, and appropriate spatio-temporal resolution to produce high-quality dynamic radiographs. This paper reports on technology development to advance the state of the art of multi-modal/multi-probe radiography and the pursuit of both deterministic and inferential (AI/ML assisted) co-analysis methodologies to produce more constrained reconstructions from multi-modal data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗