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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Variational Path Sampling of Rare Dynamical Events

This article reviews the concepts and methods of variational path sampling. These methods allow computational studies of rare events in systems driven arbitrarily far from equilibrium. Based upon a statistical mechanics of trajectory space and leveraging the theory of large deviations, they provide a perspective from which dynamical phenomena can be studied with the same types of ensemble reweighting ideas that have been used for static equilibrium properties. Applications to chemical, material, and biophysical systems are highlighted.

Singh, Aditya N↗
DOE OSTI2020

Enhancing the efficiency benefit of thermal barrier coatings for homogeneous charge compression ignition engines through application of a low-k oxide

Prior experiments reported by the authors have proven the hypothesis that achieving a dynamic temperature swing on the combustion chamber surface will lead to improved thermal and combustion efficiencies of the homogeneous charge compression ignition engine. A thin layer of yttria-stabilized zirconia, roughly 150 μm, was plasma sprayed on the piston top. It led to markedly advanced ignition and heat release in the gasoline homogeneous charge compression ignition engine, accompanied with reduced unburnt hydrocarbon and carbon monoxide emissions, improved combustion efficiency, and a higher thermal efficiency. A related computational study highlighted the critical role of coating thermal conductivity in achieving a desired dynamic response; hence, the second phase of experimental investigations focused on introducing structured porosity in the yttria-stabilized zirconia coating, as a means of reducing effective conductivity. Indeed, additional incremental improvements were observed, as well as limitations related to adverse effects of the surface roughness and the fuel interactions with the surface roughness and open pores. Erosion can also be a problem in a direct injection engine. Therefore, the third round of investigations focused on a material with a natively low conductivity (low-k), sprayed on the top of an Al piston in a relatively dense form, and in a way that yields a smooth surface. The objective was to capitalize on the low conductivity, while avoiding the pitfalls accompanying high-porosity formulations. The heat-storage capacity was limited by keeping the thickness relatively low. The results verify the paramount importance of thermal conductivity in the context of high “temperature swing” behavior and indicate a potential to improve the homogeneous charge compression ignition engine’s combustion efficiency roughly 1.5%, with the overall indicated efficiency improvement on the order of 5%, on a relative basis. In addition, the low-k oxide thermal barrier applied to the piston extended significantly the low-load homogeneous charge compression ignition operability limit.

42 ENGINEERING↗
DOE OSTI2023

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Investigation of Thermal Radiation under Pressurized Oxy-combustion Conditions

Thermal radiation of the gaseous and particle phases in a pilot-scale pressurized oxy-combustor is computationally studied. In particular, the radiation characteristics of gases and particles are estimated by employing the statistical narrow-band model and the large-particle model. It is found that thermal radiation of the particle cloud dominates in the combustor under a furnace temperature of 1500 K and when there is no substantial loss of particles to the walls. Another important observation is that radiation from the gas and particles can be approximately treated as a graybody under these conditions. More specifically, the results on the spectral radiation intensity of a gas comprising 40% (vol) H 2 O and 60% CO 2 show that when the pressure is increased to 15 bar, and the radiation pathlength is 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by employing the large-particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of the particles. Finally, this outcome of the present study is expected to be used to validate the assumption of the gray-gas model adopted in the numerical simulations of pressurized oxy-combustion.

large-particle model↗
DOE OSTI2022

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗
DOE OSTI2019

Scalable Polymer Precursor Based LiPON Type Li-Metal Interface Enabled High Performance Li-Sulfur Battery

The Battery 500 program objective is to research, develop, and demonstrate lithium battery technologies capable of achieving 500 Wh/kg with stable lifetime of 1000 cycles. This Battery 500 seedlings project has demonstrated a new approach and a pathway toward the program’s goals in just one awarded and executed phase. Cells have been assembled, tested, and evaluated. This report details efforts to develop polymer precursor electrolytes that offer properties anticipated to be similar or superior to LiPON glasses. Such precursors offer the potential to be used to process LiPON-like thin glass/ceramic coatings for use in lithium batteries. LiPON glasses provide a design basis for the synthesis of sets of oligomers/polymers by lithiation of OP(NH 2 ) 3-x (NH) x [from OP(NH) 3 ], OP(NH 2 ) 3-x (NHSiMe 3 ) x and [P=N] 3 (NHSiMe 3 ) 6-x (NH) x . Polymer electrolytes impregnated in/on celgard exhibit Li+ conductivities up to ~1×10 -5 S/cm at room temperature and are thermally stable to ≈150 °C. A Li-S battery assembled using a Li 6 SiPON composition polymer electrolyte exhibits an initial reversible capacity of 1500 mAh/g sulfur and excellent cycle performance at 0.25 and 0.5 C rate over 120 cycles, declining to 750 mAh/g at after a 220 h test period. At the same time, computational studies were performed to understand the conduction mechanism of LIPON and its interaction with Li metal. These initial project results provide a foundation for further development of the system to achieve the program goal.

25 ENERGY STORAGE↗