Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2026

Evidence of topological charge polarization at active-passive interfaces in acoustically powered active liquid crystals

Theoretical and computational studies predict topological charge separation at active-passive interfaces in active nematics, but reliable experimental validation has been lacking. We utilize a synthetic acoustically energized liquid crystal to experimentally investigate the dynamics of topological defects across an active-passive interface. The interference pattern induced by the acoustic wave inside the experimental cell creates a spatial distribution of activity, resulting in the formation of active-passive interfaces with the liquid-crystalline director field aligned parallel to those interfaces. The activity gradient drives the reorientation of positively charged topological defects along its direction, causing them to migrate into the passive zone and form a layer with net topological charge. That layer “discharges” upon cessation of activity through annihilation of positive and negative topological charges, resembling the discharge of a capacitor.

Active defects↗
DOE OSTI2020

Topological surface states in strained Dirac semimetal thin films

We computationally study the Fermi arc states in a Dirac semimetal, both in a semi-infinite slab and in the thin-film limit. We use Cd 3 As 2 as a model system, and include perturbations that break the C-4 symmetry and inversion symmetry while mirror symmetries remain. The surface states are protected by the mirror symmetries present in the bulk states and thus survive these perturbations. Additionally, the Fermi arc states persist down to very thin films, thinner than presently measured experimentally, but are affected by breaking the symmetry of the Hamiltonian. Recent experimental observations of transport in Cd 3 As 2 films strongly suggest transport by surface states. Our findings are consistent with those observations in that we show that Fermi arc states can survive to rather thin films, and our results also suggest that symmetry-breaking terms that preserve the Fermi arc states nevertheless can have a profound effect in the thin film limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2022

Defect Structure in Quantum-Cutting Yb3+-Doped CsPbCl3 Perovskites Probed by X-Ray Absorption and Atomic Pair Distribution Function Analysis

Ytterbium-doping in all-inorganic lead-halide perovskites (CsPb(Cl1-xBrx)3) generates novel properties including quantum cutting and narrow line emission, making these materials attractive spectral down-converters for solar photovoltaics. The relationship between this optical efficiency and the defect structure(s) associated with Yb3+ dopants within perovskites is not well understood. Various charge-neutral doping motifs have previously been proposed and studied computationally, including clusters involving two substitutional Yb3+ ions charge-compensated by a single local Pb2+ vacancy. Near-band-edge defect states associated with such motifs are believed to play an important mechanistic role in quantum cutting itself. Here, we report the results of X-ray absorption and X-ray total-scattering measurements on ytterbium-doped CsPbCl3. XANES shows that the dopant oxidation state is exclusively Yb3+, and a combination of Yb L3 and Pb L3 EXAFS shows that this Yb3+ substitutes exclusively at Pb2+ sites, where it adopts a pseudo-octahedral [YbCl6]3- coordination environment. Shell-by-shell fits to the data show a short Yb-Cl bond distance of 2.58 Å compared to the Pb-Cl bond distance of 2.83 Å. We confirm this finding by X-ray pair distribution function analysis, which also shows evidence of additional Pb2+ vacancy formation induced by Yb3+ doping. We evaluate whether this is the primary mechanism of charge compensation using simulated EXAFS and pair distribution function data for several computed defect structures. Together, these results resolve the local dopant structures and charge-compensation mechanisms in lanthanide-doped all-inorganic lead-halide perovskites, and, as such, significantly advance the understanding of structure-function relationships in this important class of materials.

Kluherz, Kyle (ORCID:0000000279865167)↗
DOE OSTI2022

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Charge-induced phase transition in encapsulated $\mathrm{HfTe_2}$ nanoribbons

Nanotube encapsulation is a powerful technique for coaxing solids into unconventional configurations. By synthesizing materials within the interior confines of hollow nanotubes, lower dimensional morphologies, such as one-dimensional chains or nanoribbons, are favored. We have used carbon nanotube encapsulation to realize ultranarrow, atomically precise HfTe 2 nanoribbons. A local, electron-beam-stimulated transition from the metallic 1T phase to the previously unreported semiconducting 1H phase is observed. Here, we study computationally how charging can drive the phase transition and the stability of the different atomic configurations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗
DOE OSTI2026

Enhanced sensitivity to trace 238 U impurity of sapphire via coincidence neutron activation analysis

Sapphire has mechanical and electrical properties that are advantageous for the construction of internal components of radiation detectors such as time projection chambers and bolometers. However, it has proved difficult to assess its 232 Th and 238 U content down to the picogram per gram level. Here, this work reports an experimental verification of a computational study that demonstrates 𝛾⁢𝛾 coincidence counting, coupled with neutron activation analysis (NAA), can reach ppt sensitivities. Combining results from 𝛾⁢𝛾 coincidence counting with those of earlier single-𝛾 counting based NAA shows that a sample of Saint Gobain sapphire has 232 Th and 238 U concentrations of <0.26 ppt and <2.3 ppt, respectively; the best constraints on the radiopurity of sapphire.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Computational analysis of interface-driven spin-orbit coupling in molecular adsorbates on transition metal dichalcogenides

Spin-orbit coupling (SOC) lifts molecular orbital degeneracy, enabling bilevel electronic platforms suitable for next-generation digital devices. However, common light-atom molecular feedstocks exhibit weak SOC due to the absence of heavy elements. To enhance SOC without synthesizing new materials, we leverage interfacial interactions between molecules and transition-element-based solid-state materials. This computational study investigates SOC splitting in metal-phthalocyanine adsorbed on transition metal dichalcogenides (TMDs) using density functional theory (DFT). The enhanced SOC splitting is attributed to strong orbital hybridization at the molecule-substrate interface. Specifically, Zn-phthalocyanine (ZnPC) on monolayer Mo⁢S 2 achieves a notable SOC splitting of ∼8 meV. Furthermore, when ZnPC forms self-assembled chains on Mo⁢S 2 , the splitting increases to ∼20 meV, driven by the formation of hybrid bands modulated by molecular periodicity. Furthermore, these findings highlight the role of interfacial and intermolecular interactions in inducing and enhancing SOC in surface-adsorbed molecules, providing a new strategy for molecular spintronic materials without complex synthetic efforts.

2-dimensional systems↗
DOE OSTI2025

Investigating the Combustion Performance of Dual Fuel Combustion with Diesel and Port Injected Hydrogen in a Large Bore Locomotive Engine

The heavy-duty transportation sector has primarily relied on conventional diesel combustion engines given their reliability and high thermal efficiency relative to spark ignition engines, but increased focus on reducing greenhouse gas emissions has led to investigation into alternative fuels. Gaseous hydrogen fuel has garnered a great deal of recent interest in the engine community given it has zero carbon, but hydrogen is not available at the scale and cost that petroleum fuels are currently available, and this is a barrier to adoption for industries that are looking to decarbonize their operations. Because of the fuel flexibility provided, dual fuel technology offers a pathway for some industries to adopt hydrogen as a fuel source while maintaining sufficient flexibility in times and locations where the new fuel is not yet available. This computational study investigates dual fuel combustion in a large bore locomotive engine architecture using direct injected diesel and port injected gaseous hydrogen fuel. With an optimal port fuel injection configuration from previous work, simulations of varying substitution ratio, compression ratio, manifold air temperature, diesel injection timing, and diesel injection pressure were performed to understand their effect on combustion performance. Results indicated that both increased substitution ratio and higher intake air temperature accelerates hydrogen flame propagation and can result in high peak cylinder pressures. Additionally, diesel injection timing and injection pressure were demonstrated as effective methods for controlling dual fuel combustion heat release rates.

ODonnell, Patrick Christopher↗
DOE OSTI2025

Modeling diffusion and depletion in high-aspect-ratio atomic layer deposition processes: Process parameters and manufacturing impacts

Atomic layer deposition (ALD) is a powerful technique for modifying the surface chemistry and properties of substrates with complex and nonplanar topologies. However, achieving uniform and conformal deposition on ultrahigh-aspect-ratio substrates remains challenging, typically requiring large quantities of precursors and long exposure times. Furthermore, process optimization is often performed empirically and involves substantial trial and error. In this work, we perform a combined experimental and computational study of ALD Al 2 O 3 infiltration into silica aerogel monoliths (aspect ratio >10 5 ). A reaction-diffusion model is used to explore the effects of key processing parameters, namely, exposure time per dose, precursor source temperature, number of aerogels in the reactor, and reactor volume. The model is based on quasi-static mode ALD, where the dosed precursor is held in the chamber for a fixed period of time before purging. We analyze the trade-offs between process throughput and precursor utilization for each of these parameters. Furthermore, we investigate the co-optimization and interactions between multiple process parameters, demonstrating the potential for further improvements. Furthermore, this physics-based model can be used to identify a set of process parameters for high-aspect-ratio ALD that meet specific manufacturing objective functions, including throughput, cost, and sustainability.

Aerogel↗
DOE OSTI2023

CO electroreduction on single-atom copper

Electroreduction of carbon dioxide (CO 2 ) or carbon monoxide (CO) toward C 2+ hydrocarbons such as ethylene, ethanol, acetate and propanol represents a promising approach toward carbon-negative electrosynthesis of chemicals. Fundamental understanding of the carbon–carbon (C-C) coupling mechanisms in these electrocatalytic processes is the key to the design and development of electrochemical systems at high energy and carbon conversion efficiencies. Here, we report the investigation of CO electreduction on single-atom copper (Cu) electrocatalysts. Atomically dispersed Cu is coordinated on a carbon nitride substrate to form high-density copper–nitrogen moieties. Chemisorption, electrocatalytic, and computational studies are combined to probe the catalytic mechanisms. Unlike the Langmuir-Hinshelwood mechanism known for copper metal surfaces, the confinement of CO adsorption on the single-copper-atom sites enables an Eley-Rideal type of C-C coupling between adsorbed (*CO) and gaseous [CO(g)] carbon moxide molecules. The isolated Cu sites also selectively stabilize the key reaction intermediates determining the bifurcation of reaction pathways toward different C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Manufacture and testing of biomass-derivable thermosets for wind blade recycling

Wind energy is helping to decarbonize the electrical grid, but wind blades are not recyclable, and current end-of-life management strategies are not sustainable. Here, to address the material recyclability challenges in sustainable energy infrastructure, we introduce scalable biomass-derivable polyester covalent adaptable networks and corresponding fiber-reinforced composites for recyclable wind blade fabrication. Through experimental and computational studies, including vacuum-assisted resin-transfer molding of a 9-meter wind blade prototype, we demonstrate drop-in technological readiness of this material with existing manufacture techniques, superior properties relative to incumbent materials, and practical end-of-life chemical recyclability. Most notable is the counterintuitive creep suppression, outperforming industry state-of-the-art thermosets despite the dynamic cross-link topology. Overall, this report details the many facets of wind blade manufacture, encompassing chemistry, engineering, safety, mechanical analyses, weathering, and chemical recyclability, enabling a realistic path toward biomass-derivable, recyclable wind blades.

17 WIND ENERGY↗