Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Study”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Trapping Ca + inside a molecular cavity: computational study of the potential energy surfaces for Ca + -[ n ]cycloparaphenylene, n = 5–12

Ion trap quantum computing utilizes electronic states of atomic ions such as Ca + to encode information on to a qubit. To explore the fundamental properties of Ca + inside molecular cavities, we describe here a computational study of Ca + bound inside neutral [n]-cycloparaphenylenes (n = 5–12), often referred to as “nanohoops”. This ab initio study characterizes optimized structures, harmonic vibrational frequencies, potential energy surfaces, and ion molecular orbital distortion as functions of increasing nanohoop size. Here the results of this work provide a first step in guiding experimental studies of the spectroscopy of these ion-molecular cavity complexes.

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

Computational Study of Variable Fuel-Air Ratio and Hydrogen Doping in a Rotating Detonation Engine

Rotating detonation engines (RDEs) hold great promise as a technology for powering aviation engines. The pressure gain that results from continuously traversing detonation waves enables high theoretical performance compared to deflagration-based combustion. However, the harsh operating conditions inside an RDE restrict the availability of suitable experimental diagnostics to study these systems in detail and the range of spatial and temporal scales render numerical simulations of these systems difficult. In this study, we have combined adaptive mesh refinement and a robust shock-capturing approach to study RDE systems computationally using the compressible, reactive solver PeleC. Multi-species transport along with compressible Navier-Stokes equations are solved in the model along with finite-rate chemistry. Complex geometries are represented using an embedded boundary method with second-order spatial accuracy and the system is evolved in time using a second-order Runge-Kutta method. We present studies of a methane-air RDE system that is doped with varying levels of hydrogen. Our analysis focusses on how fuel mixture and doping levels effect the detonation flame structure and the generation of multiple wave modes.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Computational Study of Additively Manufactured Internally Cooled Airfoils for Industrial Gas Turbine Applications

Internal cooling features such as pin-fins, impingement jets, and rib-turbulators are necessary to keep turbine components cool, but if sufficiently advanced can potentially also eliminate the need for film cooling on turbine blades particularly in industrial gas turbines where temperatures are not extreme. Furthermore, by leveraging additive manufacturing, other advanced designs such as lattice and incremental impingement configurations are possible and have recently been experimentally tested. While the performance of such configurations has been quantified through means of overall cooling effectiveness, it is not as clear why certain designs were better than others, or what the mechanisms were behind the observed external cooling patterns. The purpose of this study was to computationally analyze different advanced turbine blade internal cooling designs previously tested by the National Energy Technology Laboratory.

CFD↗
DOE OSTI2023

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

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

Formation of Zeolites Responsible for Waste Glass Rate Acceleration: An Experimental and Computational Study for Understanding Thermodynamic and Kinetic Processes

This report summarizes work on the NEUP Project 18-15496, DE-NE0008774 (Formation of Zeolites Responsible for Waste Glass Rate Acceleration: An Experimental and Computational Study for Understanding Thermodynamic and Kinetic Processes). The content of the report is from four journal articles that resulted from work on the project.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Pd(II) and Rh(I) Catalytic Precursors for Arene Alkenylation: Comparative Evaluation of Reactivity and Mechanism Based on Experimental and Computational Studies

We combine experimental and computational investigations to compare and understand catalytic arene alkenylation using the Pd(II) and Rh(I) precursors Pd(OAc) 2 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OAc)] 2 with arene, olefin and Cu(II) carboxylate at elevated temperature (> 120 °C). Under specific conditions, previous computational and experimental efforts have identified heterotrimetallic cyclic PdCu 2 (η 2 -C 2 H 4 ) 3 (μ-OPiv) 6 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OPiv) 2 ] 2 (µ-Cu) (OPiv = pivalate) species as likely active catalysts for these processes. Further studies of catalyst speciation suggest a complicated equilibrium between Cu(II)-containing complexes containing one Rh or Pd atom with complexes containing two Rh or Pd atoms. At 120 °C, Rh catalysis produces styrene > 20-fold more rapidly than Pd. Also, at 120 °C, Rh is ~98% selective for styrene formation while Pd is ~82% selective. Furthermore, our studies indicate that Pd catalysis has a higher predilection toward olefin functionalization to form undesired vinyl ester, while Rh catalysis is more selective for arene/olefin coupling. However, at elevated temperatures, Pd converts vinyl ester and arene to vinyl arene, which is proposed to occur through low valent Pd(0) clusters that are formed in situ. Regardless of arene functionality, the regioselectivity for alkenylation of mono-substituted arenes with the Rh catalyst gives an approximate 2:1 meta:para ratio with minimal ortho C–H activation. In contrast, Pd selectivity is significantly influenced by arene electronics with electron-rich arenes giving an approximate 1:2:2 ortho:meta:para ratio while the electron deficient (α,α,α)-trifluorotoluene gives a 3:1 meta:para ratio with minimal ortho functionalization. Kinetic intermolecular arene ethenylation competition experiments find that Rh reacts most rapidly with benzene, and the rate of mono-substituted arene alkenylation does not correlate with arene electronics. In contrast, with Pd catalysis, electron-rich arenes react more rapidly than benzene while electron-deficient arenes react less rapidly than benzene. These experimental findings, in combination with computational results, are consistent with the arene C–H activation step for Pd catalysis involving significant η 1 -arenium character due to Pd-mediated electrophilic aromatic substitution character. In contrast, the mechanism for Rh catalysis is not sensitive to arene substituent electronics, which we propose indicates less electrophilic aromatic substitution character for the Rh-mediated arene C–H activation.

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

Experimental and Computational Studies of the Production of 1,3-Butadiene from 2,3-Butanediol Using SiO2-Supported H3PO4 Derivatives

Silica-supported phosphoric acid and metal phosphate catalyzed 1,3-butadiene (BDE) production from 2,3-butanediol (2,3-BDO) was studied using experimental and computational techniques. The catalyst was initially tested in a continuous flow reactor using commercially available 2,3-BDO, leading to maximum BDE yields of 63C%. Quantum chemical mechanistic studies revealed 1,2-epoxybutane is a kinetically viable and thermodynamically stable intermediate, supported by experimental demonstration that this epoxide can be converted to BDE under standard reaction conditions. Newly proposed E2 and SN2' elementary steps were studied to rationalize the formation of BDE and all detected side-products. Additionally, using quantum mechanics/molecular mechanics (QM/MM) calculations, we modeled silica-supported phosphate catalysts to study the effect of the alkali metal center. Natural population analysis showed that phosphate oxygen atoms are more negatively charged in CsH2PO4/SiO2 than in H3PO4/SiO2. In combination with temperature-programmed desorption experiments using CO2, the results of this study suggest that the improved selectivity achieved when adding the metal center is related to an increase in the basicity of the catalyst.

butadiene↗
DOE OSTI2021

Experimental and computational study of the microporous layer and hydrophobic treatment in the gas diffusion layer of a proton exchange membrane fuel cell

Enabling fuel cell operation at high current density is critical for developing competitive alternative power system to replace internal combustion engine. However, liquid water management continues to be a challenge for high humidity or high current density operation. Water condensation in the porous media hinders efficient oxygen transport to the catalyst layer, which in turns, reduces fuel cell performance. To improve water management capability, the gas diffusion layer is often impregnated with Polytetrafluoroethylene (PTFE) and coated with a thin microporous layer, which have shown to improve fuel cell performance, especially under wet conditions. However, the fundamental mechanism that drives the performance enhancement is still not well understood. In this work, the effects of PTFE impregnation and MPL were studied using both experimental and computational techniques. Both limiting current and polarization tests under dry and wet operating condition are conducted to study the oxygen transport resistance and fuel cell performance. In addition, a 2-D, two-phase, multi-physics PEMFC model is developed to simulate performance and gain a fundamental understanding of local water saturation and oxygen concentration. The combined results show that 5 wt% PTFE impregnation with MPL significantly enhances liquid water management, which enables higher current density operation of a fuel cell.

25 ENERGY STORAGE↗
DOE OSTI2021

Laser-driven flyer plate impact: Computational studies guided by experiments

We present a computational approach using a multimaterial, arbitrary Lagrangian–Eulerian code termed ALE3D to model the nanosecond/micrometer dynamics of the launch of 0.5–4.5 km/s laser-driven metal flyer plates and the impact with stationary targets of Pyrex and fused silica glasses, and Lexan and Plexiglas polymers, producing pressures in the target in the range of 5–20 GPa. The simulations are compared to experimental results where the flyer velocity profile and the velocity profile imparted to the target material were measured with high-speed velocimetry. The experimental flyer launch by a high-intensity pulsed laser is modeled by depositing heat into a thin vaporizable layer under the flyer plate. This model produces a flyer plate that has not been exposed to the laser pulse, allowing us to compare the properties of the real flyer to a simulated ideal flyer. The simulations of target impact are generally in good agreement with the experiment except at the highest impact velocities where the shock release process in the experiment is slower than that in the simulation. The cause of this disagreement is attributed to an inadequate description of the shock viscosity during the nanosecond unloading process.

Stekovic, Svjetlana (ORCID:0000000163878070)↗
DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts – often in collaboration with experimental groups – to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO 2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

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

Experimental and computational studies of melting of the spinel phase in the Fe–Al–O ternary system

The melting behavior of spinel in the Fe–Al–O system at high temperatures (1500–1800 °C) was studied by a combination of experimental and computational investigations. Differential thermal analysis (DTA) at ultra-high temperatures coupled with cooling traces on CO 2 laser-heated levitated samples provided melting temperatures and the heats of fusion of (Fe,Al) 3 O 4 spinel phases. The experimental results are in fair agreement with the predictions using a published CALPHAD description and areas for modeling improvement are identified. New insights into the melting of defect spinel are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental and computational studies on high-entropy carbide MoNbTaVWC 5 under high pressures

High-entropy carbide, MoNbTaVWC 5 , was synthesized from oxide precursors of the constituent metals, mixed with graphite powder in a microwave-generated hydrogen plasma at 26.66 kPa and 2100 °C. Ambient x-ray diffraction analysis confirms the full conversion of oxide precursors into a single-phase, face-centered cubic structure with a lattice parameter a = 4.3309 Å. Nanoindentation measured a hardness of 24.5 ± 1.3 GPa and an elastic modulus of 386 ± 22 GPa. The synthesized sample, mixed with a copper pressure marker, was studied by the radial x-ray diffraction technique with beryllium gasketing in a diamond anvil cell up to 70 GPa. The experimentally measured pressure–volume curve and shear strength were compared with theoretical predictions using the special quasi-random structure technique and density functional theory. MoNbTaVWC 5 achieved a 12% volume compression at 70 GPa and exhibited a high shear strength of 6.6 GPa. The present study demonstrates that the high-entropy carbide MoNbTaVWC 5 exhibits exceptional incompressibility and high strength under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and Computational Studies of Carbon–Carbon Bond Formation via Ketonization and Aldol Condensation over Site-Isolated Zirconium Catalysts

We report here the preparation and investigation of isolated Zr centers supported on a high surface area silica for the conversion of carboxylic acids to internal ketones by ketonic decarboxylation (ketonization) and the aldol condensation of ketones to dimeric enones. Catalysts were synthesized by the grafting of Cp 2 ZrCl 2 on the surface of amorphous silica. The connectivity of Zr was characterized by XRD, UV-vis, and Raman spectroscopy. For the lowest Zr loading, Zr is present predominantly as isolated monomeric species. As the Zr loading is increased, a progressively larger fraction of Zr forms oligomeric species and ZrO nanoparticles. Measurements of catalytic activity show that the turnover frequency for carboxylic acid ketonic decarboxylation reaction and aldol condensation of ketones decreases monotonically with increasing Zr loading. An H/D kinetic isotope effect was not observed over isolated Zr catalysts, suggesting that α-H abstraction is not the rate-determining step, rather C-C bond forming may be rate limiting for both reactions. This conclusion is supported by computational modeling of the reaction mechanism. The proposed catalytic cycle for ketonization proceeds via a β-keto acid intermediate on isolated Zr sites that are always coordinatively saturated with C-C bond formation as the rate-limiting step. C-C bond formation is also rate-determining for aldol condensation, with an apparent activation energy that is in good agreement with the experiment if the resting state is a saturated ≡ZrOH site with two adsorbed ketone molecules.

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

Computational study of the adsorption of bimetallic clusters on alumina substrate

We performed computational investigation of the adsorption of bimetallic Pd 3 M 2 (where M changes from Ag, Au, Co, Cu, Mn, Ni, Pt, and Ru) cluster on the hydroxylated alumina surface. Previously, it was shown that small silver cluster can control the rate of discharge at the cathode of lithium-oxygen battery. The gap near the fermi energy was shown to control the oxygen reduction, an important reaction for LiO 2 formation. Controlling the gap would ultimately control the rate of LiO 2 formation. One can vary the size of the cluster to vary the gap, however, this "knob" provides limited variance in the gap. Alloying, in combination with size variation offer a much wider control of the gap, hence the LiO 2 formation. Using Density Functional Theory (DFT), we determined the most stable geometry of the bimetallic clusters Pd 3 M 2 and calculated the binding energies of these clusters on the alumina substrate which ranges from 0.2 eV to 0.25 eV depending upon the composition of the alloy-cluster, its orientation and the adsorption site. We also find that Pd atoms bind strongly with the substrate oxygen atoms with an average short bond-length of about 2.2 angstrom. We explored how the gap at the Fermi level of the system varies as a function of elemental composition and the calculated gap ranges from 0 meV to 90 meV. Charges distribution using Bader analysis was also performed to probe how charges are transferred between the cluster and the substrate. These preliminary results will open the door for more systematic studies of alloy clusters of different size and stoichiometry for Li-O 2 battery cathode design.

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

High-Fidelity Computational Studies of Intermediate Pressure Capacitively Coupled Plasmas

This report discussed progress made in the development of a computational model for intermediate pressure capacitively coupled plasma discharge. A new full fidelity argon chemistry model has been developed that includes electron, argon ion, three lumped states for electronically excited species, and the background ground state argon. As such the mechanism has been validated over a wide range of intermediate pressure conditions from 100 mTorr to a few Torr. Uncertainties in the electron impact cross sections and the rate coefficients for the reactions are propagated through the model to establish uncertainties that emerge in the model due the reaction mechanism and its impact on the discharge properties. The model is therefore validated in face of uncertainties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mechanisms in the Synthesis of S -Alcohols with 1,4-NADH Biomimetic Co-factor N-Benzyl-1,4-dihydronicotinamide using Horse Liver Alcohol Dehydrogenase: A Hybrid Computational Study

The enantioselective reduction of prochiral ketones catalyzed by horse liver alcohol dehydrogenase (HLADH), was investigated via a hybrid computational approach, for molecular reactions involved in chiral synthesis of S-alcohols, when the natural co-factor, 1,4-dihyronicotinamide adenine dinucleotide, 1,4-NADH, was replaced with biomimetic co-factor, N-benzyl-1,4-dihydronicotinamide, 1. We surmised that different hydride and proton transfer mechanisms were involved using co-factor, 1. An alternative mechanism, where the hydride transfer step occurred, via an η 1 -keto-S-η 2 -5,6-1,4-dihydronicotinamide-Zn(II) complex, was previously investigated with a model of the HLADH-Zn(II) catalytic site (J. Organometal. Chem. 2021, 943, 121810). Presently, we studied canonical and alternative mechanisms compared to models of the entire enzyme structure. We disproved the η 2 -Zn(II) complex, and discovered a canonical hydride transfer from biomimetic 1,4-NADH, 1, to the Zn(II) bound prochiral ketone substrate, followed by a new proton relay, consisting of a water chain connecting His51 to Ser48 that accomplished the S-alkoxy anion's protonation to yield the final S-alcohol product. The HLADH catalysis, with biomimetic co-factor, 1, that replaced the ribose group, the 5'-diphosphate groups, and the adenine nucleotide with a N-benzyl group, has provided a new paradigm for the design of other structures of 1,4-NADH biomimetic co-factors, including their economic value in biocatalysis reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗