Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 433 records · Page 24

Dynamic modelling and control strategy of a temperature-driven metal hydride cooling system for buildings

A temperature-driven coupled metal hydride (MH) based thermal energy storage (TES) system can allow to shave and shift the peak energy demand in buildings. The high energy density and long-term (seasonal) energy storage capability are its major advantages over other energy storage methods. The dynamic nature of the MH operation, however, requires controlled hydrogen transfer between the coupled MHs at a rate needed to meet the building's transient load. While temperature-driven MH systems are studied in the literature, their application in buildings and control are scarcely reported. Here, this paper presents a control-based dynamic modeling of the temperature-driven coupled MH-TES system for building cooling applications. The dynamic model is developed in MATLAB(R) Simulink environment, considering the thermodynamic and kinetic behaviors of the MH systems. Based on a preliminary analysis of a property database of over 337 hydrides, we select around 1600 MH pairs suitable for building cooling applications. Each of these MH pairs is studied for their performance using the dynamic model, and among all, Zr 0.76 Ti 0.24 Ni 1.16 Mn 0.63 V 0.14 Fe 0.18 -Ti 0.85 Zr 0.15 Cr 1.2 Mn 0.8 MH pair showed fast dynamics along with high coefficient of performance (COP) of 0.71. A parametric investigation is performed on this MH pair to understand the effect of operating temperatures. Finally, three proportional-integral (PI) feedback controllers are investigated to regulate the temperature, pressure and mass exchange between the coupled MH pairs. The developed PI controller is sufficiently capable of rejecting the signal noise from the hydrogen flow and internal heat exchange processes with root mean square error of 5.78 W between reference and actual cooling load.

08 HYDROGEN↗

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Georgia Tech Accelerated, Compressed, and Regularized Compute of Kinetic-based PDEs (Final Report)

This report summarizes the collaborative effort between Lawrence Livermore National Laboratory and Georgia Tech to enhance the BoBa library for tensor train computation in PDE solvers, with a target on kinetic equations and their continuum limits. We aimed to reduce computational cost and memory usage by replacing traditional array-based computations with tensor trains. We examined the compressibility of time-evolving solutions to the Euler equations with discontinuities. We also explored using the first invsicid and linear regularization of the compressible flow equations via the information geometric regularization (IGR). We explored this in a tensor train formulation. To identify that inverse terms in the IGR equations pose problems for tensor train formulations and investigate efficient methods for batched inversion of tensor trains.

97 MATHEMATICS AND COMPUTING↗

Visualizing Millisecond Atomic Dynamics of Nanocrystals in Liquid

Atomic structures of nanomaterials are inherently dynamic and continuously reshaped through interactions with chemical species and external stimuli. Such dynamics are further amplified as the size and dimensionality of nanomaterials decrease. Despite advances in analytical methods, it remains challenging to capture the structural dynamics of nanomaterials in reactive environments with both atomic spatial resolution and commensurate temporal resolution. Here, in this study, we directly visualize atomic-scale dynamics of gold (Au) nanocrystals in reactive liquid environments with millisecond-speed liquid-cell electron microscopy (EM) and deep-learning denoising. We uncover reversible fluctuations in the local crystallinity of Au nanocrystals dependent on the surrounding chemical environment. These transient fluctuations, driven by interactions at nanocrystal–liquid interfaces, critically influence the dissolution kinetics and grain boundary relaxation. By overcoming the spatiotemporal limitations in conventional liquid-cell EM, our findings provide insights into how transient nanoscale structures dictate the stability and reactivity of nanomaterials.

Kang, Sungsu [University of Chicago, IL (United St↗

Planar laser-induced fluorescence system for spatiotemporal ion velocity distribution function measurements

In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.

E×B plasma↗

Hot-Carrier Injection and Millisecond Charge Separation from a Robust Heteroleptic Iron(II) Chromophore Immobilized on TiO 2

The synthesis, spectroscopic characterization, computational analysis, and photoelectrochemical behavior of a new iron-based chromophore, [(Cpy) 2 Fe- (deeb)](PF 6 ) 2 (Fe(Cpy) 2 (deeb)), where Cpy is 1-methyl-3-(2-pyridyl)imidazole and deeb is 4,4’-(CO 2 CH 2 CH 3 ) 2 -2,2’-bipyridine, is reported. Electrochemically reversible waves assigned to a metal-centered E o (Fe III/II ) = +0.48 and a ligand-centered E o (Fe 2+/+ ) = −1.47 V vs Fc +/0 reduction were evident in cyclic voltammetry measurements. The combination of a strong σ-donor and a π-acceptor lowered the energy of the metal-toligand charge-transfer (MLCT) excited state relative to the metal-centered state. Two MLCT transitions appear in the visible region at 424 and 580 nm. TDDFT calculations revealed that the lower-energy band was well formulated as Fe(II)→deeb, and the higher-energy transition was charge transfer to both the deeb and Cpy ligands. Resonance Raman spectroscopy supports these findings showing enhanced deeb vibrational modes with 532 nm excitation, both deeb and Cpy modes with 473 nm excitation, and exclusively Cpy with 405 nm excitation. Ultrafast spectroscopy reveals a short-lived (∼2 ps) MLCT excited state and a longer-lived (∼20 ps) metal-centered state. Efficient methods to deprotect the ester groups and anchor the complex to mesoporous TiO 2 (anatase) thin films in high surface coverages, Fe(Cpy) 2 (dcb)|TiO 2 σ = 3 × 10 −8 mol/cm 2 , were established. Pulsed light excitation of Fe(Cpy) 2 (dcb)|TiO 2 resulted in rapid excited state injection (k inj > 10 8 s −1 ) and formation of a charge-separated state, Fe III (Cpy) 2 (dcb)|TiO 2 (e), which persists on the millisecond time scale before returning cleanly to the ground state with secondorder kinetics. Injection yields measured 50 ns after light excitation were found to double from Φ = 0.15 with green (532 nm) light to 0.30 with blue (457 nm) light excitation. Incident photon-to-current efficiency (% IPCE) measurements as a function of excitation wavelength in a 0.5 M LiI/I 2 /CH 3 CN electrolyte provide clear evidence for band-selective “hot carrier” injection from the remote Cpy-localized excited state. Collectively, the spectroscopic and photoelectrochemical data indicate that a semiconductor can intercept hot electrons from iron chromophores even when the excited-state dipole is oriented away from the surface-anchoring ligand.

Charge transfer↗

A Business Case Evaluation of Gas Switching Reforming (GSR) Technology: A Promising Technology for Natural Gas Reforming with Integrated CO2 Capture

Hydrogen is essential in the transition to sustainable energy, and developing low-carbon production methods is a key research focus. Traditional steam methane reforming (SMR) dominates the hydrogen industry but contributes substantially to CO2 emissions. In response, Gas Switching Reforming (GSR) has emerged as a novel process that integrates carbon capture and utilizes process heat more efficiently. Unlike other reforming methods, GSR consolidates oxidation and reduction reactions within a single reactor, which minimizes external energy inputs and simplifies scaling. Like conventional steam methane reforming (SMR), GSR can be integrated with water-gas shift and pressure swing adsorption units for pure hydrogen production. This work presents a comprehensive business case analysis of GSR technology based on experimental results in Technology Readiness Level 3, Life Cycle Assessment (LCA) and Techno-Economic (TEA) evaluation incorporating ASPEN Plus process modeling considering different configurations and energy scenarios. The TEA incorporates data from kinetic experiments from various catalysts to evaluate the GSR process under various conditions. The goal of this work is to evaluate GSR’s potential to serve as a low-carbon alternative to SMR, focusing on global warming potential and additional impact categories to evaluate a wide spectrum of environmental impacts. Comparative assessments were conducted with SMR, chemical loop reforming (CLR), and proton exchange membrane (PEM) electrolysis to explore trade-offs across environmental metrics. The environmental impact assessment of this work encompasses the entire hydrogen production lifecycle from raw material extraction to plant decommissioning, using a cradle-to-gate boundary. Preliminary findings highlight that GSR, when integrated with low-carbon energy sources, could significantly reduce environmental impacts, making it a promising candidate for low-carbon hydrogen infrastructure. The insights from this business case evaluation aim to guide industry in scale-up and commercialization of this promising clean energy technology.

03 NATURAL GAS↗

Biochemical properties of glycerol kinase from the hypersaline-adapted archaeon Haloferax volcanii

ABSTRACT Extremophilic microorganisms are promising candidates for industrial and analytical biocatalysis.Haloferax volcanii, a halophilic archaeon that prefers glycerol over glucose, channels this substrate into central metabolism through glycerol kinase (GK). Here, we report the biochemical properties ofH. volcaniiGK and its potential for biotechnological applications. An N-terminal His-tagged GK was functionalin vivoand yielded 3 mg/L culture—4.5 times more enzyme than a C-terminal StrepII-tagged version. Size exclusion chromatography revealed a glycerol-induced oligomeric shift from homodimer to a dimer-dominant state with detectable tetramer. The purified enzyme showed robust activity across broad pH and salinity ranges, with optimal activity at 100 mM NaCl and 50°C–60°C. It retained catalytic activity in 5%–10% dimethyl sulfoxide (DMSO) and crude glycerol containing methanol. His-GK was freeze-thaw stable and thermotolerant in 2 M NaCl buffers. In the absence of ligands, the enzyme’s melting temperature (T m ) was 80°C. Glycerol increased the T m to 85°C, and combinations with MgCl₂ (84°C) or ATP (88°C) provided further stabilization. The highest T m (89°C) occurred with all three ligands, suggesting a cumulative stabilizing effect. Kinetic analyses revealed positive cooperativity for glycerol, ATP, and Mg² + ; Mn² + and Co² + also supported the activity.H. volcaniiGK is the first known GK to exhibit positive cooperativity with glycerol and ATP. Its high stability and substrate flexibility support its use in biodiesel waste valorization,in vitrobiocatalysis, and biosensor development—applications demanding robust, specific, and stable enzymes. IMPORTANCE This study reveals thatH. volcaniiGK exhibits positive cooperativity for glycerol, ATP, and Mg² + , a kinetic feature not previously reported for glycerol kinases. This behavior enables steep, switch-like responses to small substrate changes, offering unique advantages for biosensor design. Importantly,H. volcaniiGK also maintains high activity under extreme salinity, temperature, broad pH, and solvent conditions that typically limit enzyme use in industrial and environmental applications. These traits make this GK an ideal candidate for enzyme-based biosensors, which often suffer from poor tolerance to pH, solvent, and thermal stress. Its robustness supports its use in cross-linked enzyme crystals, an immobilization method that enhances enzyme stability and reusability under harsh conditions. Moreover, GKs are already employed in Mg² + detection kits; however,H. volcaniiGK’s ability to tolerate and respond to diverse divalent cations (e.g., Co² + , Mn² + ) broadens their potential for pollutant detection and environmental monitoring. These features collectively positionH. volcaniiGK as a valuable biocatalyst for biosensing,in vitrodiagnostics, and biotechnological applications requiring both precision and durability.

Biotechnology & Applied Microbiology↗

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗

Polaron catastrophe within quantum acoustics

The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.

Science & Technology - Other Topics↗

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Influence of Carboxymethyl Cellulose as a Thickening Agent for Glauber’s Salt-Based Low Temperature PCM

This work is focused on a novel, promising low temperature phase change material (PCM), based on the eutectic Glauber’s salt composition. To allow phase transition within the refrigeration range of temperatures of +5 °C to +12 °C, combined with a high repeatability of melting–freezing processes, and minimized subcooling, the application of three variants of sodium carboxymethyl cellulose (Na-CMC) with distinct molecular weights (700,000, 250,000, and 90,000) is considered. The primary objective is to optimize the stabilization of this eutectic PCM formulation, while maintaining the desired enthalpy level. Preparation methods are refined to ensure repeatability in mixing components, thereby optimizing performance and stability. Additionally, the influence of Na-CMC molecular weight on stabilization is examined through differential scanning calorimetry (DSC), T-history, and rheology tests. The PCM formulation of interest builds upon prior research in which borax, ammonium chloride, and potassium chloride were used as additives to sodium sulfate decahydrate (Glauber’s salt), prioritizing environmentally responsible materials. The results reveal that CMC with molecular weights of 250 kg/mol and 90 kg/mol effectively stabilize the PCM without phase separation issues, slowing crystallization kinetics. Conversely, CMC of 700 kg/mol proved ineffective due to the disruption of gel formation at its low gel point, hindering higher concentrations. Calculations of ionic concentration indicate higher Na ion content in PCM stabilized with 90 kg/mol CMC, suggesting increased ionic interactions and gel strength. A tradeoff is discovered between the faster crystallization in lower molecular weight CMC and the higher concentration required, which increases the amount of inert material that does not participate in the phase transition. After thermal cycling, the best formulation had a latent heat of 130 J/g with no supercooling, demonstrating excellent performance. This work advances PCM’s reliability as a thermal energy storage solution for diverse applications and highlights the complex relationship between Na-CMC molecular weight and PCM stabilization.

25 ENERGY STORAGE↗

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗

Reference solutions for linear radiation transport: the Hohlraum and Lattice Benchmarks

Radiation transport describes the propagation of energetic particles through space as they interact with a surrounding material medium. In a kinetic description, radiation transport is modeled by a radiation transport equation (RTE) that prescribes the density of the radiation in position-momentum phase space. The purpose of this dataset is to provide highly resolved solutions to two benchmark problems. These two benchmarks do not possess exact solutions; moreover, the construction of a manufactured solution may require a non-physical source that is not desirable, especially if it spoils the physical nature of the solution. Thus the goal of this computational study is to provide a highly resolved reference solution for testing newer, more cost efficient methods that are currently being developed in the research community.

97 MATHEMATICS AND COMPUTING↗