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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High tensile alloy of copper to mitigate current collector deformation in silicon electrodes for lithium-ion batteries

The volumetric changes of silicon electrodes, along with the strong adhesive properties of certain binders, can lead to plastic deformation of the current collector and create damage in the electrode coating. Here, in this study, we report a detailed study of silicon coatings on a high-tensile alloy (HTA) foil of copper with strength over twice that of conventional copper foils. The HTA current collectors with high mechanical strength can mitigate plastic deformation upon continuous cycling. At moderate areal capacities (2.5–3 mAh cm −2 ), conventional copper foils show significant wrinkling after only a few electrochemical cycles, whereas the HTA foils remain intact. We demonstrate viability of the HTA foils in large format xx6395 pouch cells, in which the HTA current collectors remain intact even at an areal capacity of 4.5 mAh cm −2 ; in contrast, wrinkles form in conventional copper current collectors increasing the likelihood of lithium plating. Computational studies show that stresses generated during cycling of silicon electrodes are very high in the current collector and at the current collector-coating interface, explaining the wrinkling of conventional Cu foils. Our studies highlight importance of current collector to solve the electrochemical and chemo-mechanical performance challenges associated with high-loading silicon electrodes.

Chemo-mechanical degradation↗
DOE OSTI2022

Optical and photophysical properties of platinum benzoporphyrins with C 2v and D 2h symmetry

Platinum(II) mono- and di-benzoporphyrins carrying either ester or imide peripheral groups have been successfully synthesized. These benzoporphyrins possessing C 2v and D 2h symmetry were investigated using a range of techniques including UV-Vis and phosphorescence spectroscopy, differential pulse voltammetry, and transient absorption spectroscopy covering wide spatial and temporal ranges. Up to three oxidation and three reductions within the accessible potential window of dichlorobenzene were observed. Facile reductions for both mono- and dibenzoporphyrin derivatives, more so with imide peripheral substituents were witnessed. Computational studies were supportive of push-pull type interactions within these molecular systems. Here, the formation of triplet excited states within 5 ps after laser excitation was confirmed by femtosecond transient absorption studies. Nanosecond transient absorption studies coupled with phosphorescence lifetime studies revealed relatively rapid relaxation of the triplet excited states as compared to a controlled platinum porphyrin suggesting push-pull type interactions on a slower time scale.

DFT calculations↗
DOE OSTI2025

Cryo-EM Visualization of Intermolecular π-Electron Interactions within π-Conjugated Peptidic Supramolecular Polymers

The self-assembly of “π-peptides” – molecules with π-electron cores substituted with two or more oligopeptide chains – brings organic electronic function into biologically relevant nanomaterials. π-Peptides assemble into fibrillar nanomaterials as driven by enthalpic peptide-based hydrogen bonding networks and pi-core-based quadrupolar interactions. A large body of spectroscopic, morphological and computational studies informs on the nature of the self-assembly process and the resulting nanostructures, but detailed structural information has remained elusive. Here, inspired by the recent use of cryogenic electron microscopy (cryo-EM) to provide high-resolution structures for synthetic peptide nanomaterials, we present here the use of cryo-EM to offer ca. 3 Å resolution of π-peptide nanomaterial assemblies, visualizing for the first time the nature of the intermolecular π-core electronic interactions responsible for energy transport through these supramolecular materials.

Group theory↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Stress- and interface-compatible red phosphorus anode for high-energy and durable sodium-ion batteries

Sodium-ion battery is a promising candidate for energy storage application, but the absence of high-capacity and low-cost anode materials significantly limit their practical specific energy and cost. Red phosphorus (RP) possesses a high theoretical specific capacity but suffers from large volume change, low electronic conductivity and unstable solid-electrolyte interphase (SEI). Herein, a hierarchical micro/nanostructured antimony-doped RP/carbon anode was developed, which demonstrate extraordinary electrochemical performance with high initial Coulombic efficiency of ~90%, high areal (~1.7 mAh cm-2) capacity together with good cycle stability and rate capability. Combined experimental and computational studies consistently revealed that such a unique structural design can dramatically accommodate the mechanical stress and moreover effectively restrain the undesired decomposition of electrolyte solvents regardless of electrolytes formulation, resulting in superior structural integrity and thin/robust SEI formation during cycling. The present finding has offered an alternative strategy for stress management and interface engineering on high-capacity alloying-based anode materials.

Liu, Xiang↗
DOE OSTI2024

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

A study on computational fluid dynamics modeling of a refrigerated container for COVID-19 vaccine distribution with experimental validation

A key issue with the distribution of vaccines to prevent COVID-19 is the temperature level required during transport, storage, and distribution. Typical refrigerated transport containers can provide a temperature-controlled environment down to -30 °C. However, the Pfizer vaccine must be carefully transported and stored under a lower temperature between -80 °C and -60 °C. One way to provide the required temperature is to pack the vaccine vials into small packages containing dry ice. Dry ice sublimates from a solid to a gas, which limits the allowable transport duration. This can be mitigated by transporting in a -30 °C refrigerated container. Moreover, because the dry ice will sublimate and thereby release CO 2 gas into the transport container, monitoring the CO 2 concentration within the refrigerated container is also essential. In the present work, a 3D computational fluid dynamics model was developed based on a commercially available refrigerated container and validated with experimental data. The airflow, temperature distribution, and CO 2 concentration within the container were obtained from the simulations. The modeling results can provide guidance on preparing experimental setups, thus saving time and lowering cost, and also provide insight into safety precautions needed to avoid hazardous conditions associated with the release of CO 2 during vaccine distribution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Experimental and computational kinetics study of the liquid-phase hydrogenation of C=C and C=O bonds

Solvent effects on adsorption equilibrium and reaction kinetics are evaluated for hydrogenation reactions catalyzed by Pd/alumina in a series of different solvents. Three reactants – cyclohexene, benzene and benzaldehyde – and three solvents – n-heptane (HEP), methylcyclohexane (MCHA), decalin (DL) – have been investigated. Kinetic analysis of hydrogenation of cyclohexene reveals that hydrogen adsorbs on different sites from those where cyclohexene and the solvents adsorb; however, the presence of hydrogen on these separate sites affects the heats of adsorption of the hydrocarbons. When the solvent is a weakly interacting linear alkane (HEP), the rate determining step of the reaction is the first hydrogenation of adsorbed cyclohexene. Furthermore, this conclusion is supported by DFT calculations that show a higher enthalpy barrier for the first hydrogenation than for the second, while statistical thermodynamics analysis validates the physical significance of the entropy of adsorption parameters derived from the kinetic fitting of experimental data. By contrast, with solvents such as MCHA and DL, which interact more strongly with the metal surface and compete for active sites with the reactant and the surface intermediate, the rate limiting step seems to shift to the second hydrogenation step.

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

High-Throughput Computing: Case Study of Medical Image Processing Applications

HPC is designed for large-scale simulations using monolithic codes of tightly coupled processes highly optimized to deliver decreased time to solution. Medical image processing is not a traditional field of HPC. Similar to AI applications, medical image processing parses large datasets, typically multiple times, to support a variety of studies for classification, diagnosis or monitoring purposes. The convergence of AI, HPC and Big Data encouraged more fields using image processing to transition to HPC. However, not all applications benefit from the same optimizations. In this paper we focus on high throughput medical image processing applications that analyze a huge dataset of small MRI images and that require HPC systems to decrease the time of parsing the entire dataset and not individual MRIs. We show in this research the performance of running SLANT, an image processing application for a whole brain segmentation, on large-scale systems and highlight performance limitations. We present optimizations prioritizing throughput that exhibit a 3.5x speed-up on the Summit Supercomputer that can be used as a baseline for building a high-throughput execution framework for other HPC systems.

Predescu, Maria↗
DOE OSTI2020

Distributed heterogeneous compute infrastructure for the study of additive manufacturing systems

We present the current status of a large-scale computing framework to address the need of the multidisciplinary effort to study chemical dynamics. Specifically, we are enabling scientists to process and store experimental data, run large-scale computationally expensive high-fidelity physical simulation, and analyze these results using the state-of-the-art data analytics tools, machine learning, and uncertainty quantification methods using heterogeneous computing resources, such as CPU and GPU cluster. The framework can integrate or abstract out multiple domains based on roles. In order to develop this framework, we have leveraged an existing framework coupled with in-house heterogeneous computing resources. We present the results of using this framework on a single metadata triggered workflow to accelerate an additive manufacturing use case.

heterogeneous computing, hpc, workflows, reproduci↗
DOE OSTI2023

Experimental and Computational Heat Transfer Study of sCO 2 Single-Jet Impingement

This study experimentally and computationally investigates the heat transfer capability of supercritical carbon dioxide (sCO 2 ) single jet impingement. The evaluated jet Reynolds number range is between 80,000 and 600,000, with a nondimensional jet-to-target surface spacing of 2.8. CO 2 -impinging jet stagnation conditions were maintained at approximately 200 bar and a temperature of ∼400 °C for most experiments. The goal is to understand how changes in the aforementioned parameters influence heat transfer between the working fluid and the heated surface. Additionally, due to the elevated Reynolds numbers and difference in thermodynamic properties between air and CO 2 , air-derived impingement correlations may not be appropriate for CO 2 impingement; these correlations will be evaluated against experimental sCO 2 impingement data. At the time of this study, no sCO 2 impingement data was available relevant to sCO 2 power cycles. The target surface is a 1.5-in. diameter copper block centered on the 3 mm jet orifice. A mica heating element bolted to the bottom of the copper block provides a uniform heat flux. Thermocouples embedded in the copper block are used to determine the surface temperature. Nusselt numbers obtained from experimental sCO 2 test data are compared to area-averaged Nusselt numbers from air-derived correlations. The comparisons showed that air correlations drastically underpredict the heat transfer when sCO 2 is used as the working fluid. A modified sCO 2 correlation using experimental data at discussed conditions is derived based on an existing air correlation. A CFD study is also performed to further investigate sCO 2 heat transfer characteristics, and assess the numerical model applicability to this problem type.

42 ENGINEERING↗
DOE OSTI2024

Effects of Downstream Vortex Generators on Film Cooling a Flat Plate Fed by Crossflow

Counter-rotating vortices, formed by the interaction of film-cooling jets and the hot gas flow, adversely affect the performance of conventional film-cooling designs. Downstream vortex generators have been shown to improve cooling effectiveness by mitigating the effects of the counter-rotating vortices and by deflecting the cooling jet laterally. In this study, computational and experimental methods were used to examine how cylindrical film-cooling holes (D = 3.2 mm, L/D = 6, p/D = 3, α = 30 deg) with and without downstream vortex generators perform when the coolant supply channel is perpendicular to the direction of the hot gas. For this study, the hot gas had a temperature of 650 K and an average Mach number of 0.23. Here, the hot-gas-to-coolant temperature ratio was 1.9, and two blowing ratios (0.75 and 1.0) were studied. Results from the computational fluid dynamics study show how crossflow affects the interaction between the film-cooling jet and hot gas flow with and without downstream vortex generators. The experimental measurements were based on infrared thermography in a conjugate heat transfer environment. Results were obtained for film-cooling performance in terms of overall effectiveness, film effectiveness, and local heat transfer coefficients. The downstream vortex generators can increase the laterally averaged effectiveness by a factor of 1.5 relative to cylindrical holes, but this higher performance is restricted to low crossflow velocities and higher blowing ratios.

17 WIND ENERGY↗
DOE OSTI2020

Co(salen)-Catalyzed Oxidation of Lignin Models to Form Benzoquinones and Benzaldehydes: A Computational and Experimental Study

Lignin is a highly abundant polyphenolic polymer that imparts mechanical strength to plant biomass. Transition-metal complexes can catalyze lignin oxidation to produce value-added products, but low catalytic efficiency has hampered their use in industry. Identifying the chemical and structural factors that govern catalytic activity is a prerequisite to rational design of catalysts with improved activity. In this study, we combine computational and experimental approaches to investigate the mechanism of Co(salen)-catalyzed oxidation of the monomeric lignin models syringyl (S), vanillyl (G), and 4-hydroxybenzyl alcohol (H) to produce benzoquinone and benzaldehyde products. Experimentally, S oxidation to form dimethoxybenzoquinone proceeded efficiently with a Co(salen) catalyst coordinated by a pyridine ligand, but G and H did not undergo oxidation. Density functional theory calculations reveal that catalyst regeneration is energetically unfavorable in the presence of H, which prevents oxidation. In contrast, S readily facilitates catalyst regeneration. Formation of methoxybenzoquinone from G was achieved experimentally by adding bulky, noncoordinating bases. These findings provide a fundamental baseline for enhancing the activity of Co-Schiff base catalysts toward lignin-like molecules by adding sterically hindered nitrogenous bases or potentially by including a cocatalyst that promotes catalyst regeneration.

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

Emergent properties at oxide interfaces controlled by ferroelectric polarization

Abstract Ferroelectric materials are characterized by the spontaneous polarization switchable by the applied fields, which can act as a “gate” to control various properties of ferroelectric/insulator interfaces. Here we review the recent studies on the modulation of oxide hetero-/homo-interfaces by ferroelectric polarization. We discuss the potential applications of recently developed four-dimensional scanning transmission electron microscopy and how it can provide insights into the fundamental understanding of ferroelectric polarization-induced phenomena and stimulate future computational studies. Finally, we give the outlook for the potentials, the challenges, and the opportunities for the contribution of materials computation to future progress in the area.

Chemistry↗
DOE OSTI2024

Mechanistically Guided Materials Chemistry: Synthesis of Ternary Nitrides, CaZrN 2 and CaHfN 2

Recent computational studies have predicted many new ternary nitrides, revealing synthetic opportunities in this underexplored phase space. However, synthesizing new ternary nitrides is difficult, in part because intermediate and product phases often have high cohesive energies that inhibit diffusion. Here, we report the synthesis of two new phases, calcium zirconium nitride (CaZrN 2 ) and calcium hafnium nitride (CaHfN 2 ), by solid state metathesis reactions between Ca 3 N2 and MCl 4 (M = Zr, Hf). Although the reaction nominally proceeds to the target phases in a 1:1 ratio of the precursors via Ca 3 N 2 + MCl4 → CaMN 2 + 2 CaCl 2 , reactions prepared this way result in Ca-poor materials (Ca x M 2–x N2, x < 1). A small excess of Ca 3 N 2 (ca. 20 mol %) is needed to yield stoichiometric CaMN 2 , as confirmed by high-resolution synchrotron powder X-ray diffraction. In situ synchrotron X-ray diffraction studies reveal that nominally stoichiometric reactions produce Zr 3+ intermediates early in the reaction pathway, and the excess Ca 3 N 2 is needed to reoxidize Zr 3+ intermediates back to the Zr 4+ oxidation state of CaZrN 2 . Analysis of computationally derived chemical potential diagrams rationalizes this synthetic approach and its contrast from the synthesis of MgZrN 2 . These findings additionally highlight the utility of in situ diffraction studies and computational thermochemistry to provide mechanistic guidance for synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal Stability of π-Conjugated n -Ethylene-Glycol-Terminated Quaterthiophene Oligomers: A Computational and Experimental Study

This work represents a joint computational and experimental study on a series of $\textit{n}$-ethylene glycol (PEO$\textit{n}$)-terminated quaterthiophene (4T) oligomers for 1 < $\textit{n}$ < 10 to elucidate their self-assembly behavior into a smectic-like lamellar phase. This study builds on an earlier study for $\textit{n}$ = 4 that showed that our model predictions were consistent with experimental data on the melting behavior and structure of the lamellar phase, with the latter consisting of crystal-like 4T domains and liquid-like PEO4 domains. The present study aims to understand how the length of the terminal PEO$\textit{n}$ chains modulates the disordering temperature of the lamellar phase and hence the relative stability of the ordered structure. In this work, a simplified bilayer model, where the 4T domains are not explicitly described, is put forward to efficiently estimate the disordering effect of the PEO domains with increasing n; this method is first validated by correctly predicting that layers of alkyl (PE)-capped 4T oligomers (for 1 < $\textit{n}$ < 10) stay ordered at room temperature. Both 4T-domain implicit and explicit model simulations reveal that the order-disorder temperature decreases with the length of the PEO capping chains, as the associated increase in conformational entropy drives a tendency toward disorder that overtakes the cohesive energy, keeping the ordered packing of the 4T domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗