Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structured model reduction”

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,045 records · Page 58

Cloud model-based simulation of spaceborne radar observations

Simulations of observations from potential spaceborne radars are made based on storm structure generated from the three-dimensional (3D) Goddard cumulus ensemble model simulation of an intense overland convective system. Five frequencies of 3, 10, 14, 35, and 95 GHz are discussed, but the Tropical Rainfall Measuring Mission precipitation radar sensor frequency (14 GHz) is the focus of this study. Radar reflectives and their attenuation in various atmospheric conditions are studied in this simulation. With the attenuation from cloud and precipitation in the estimation of reflectivity factor (dBZ), the reflectivities in the lower atmosphere in the convective cores are significantly reduced. With spatial resolution of 4 km X 4 km, attenuation at 14 GHz may cause as large as a 20-dBZ difference between the simulated measurements of the peak, Z(sub mp) and near-surface reflectivity, Z(sub ms) in the most intense convective region. The Z(sub mp) occurs at various altitudes depending on the hydrometeor concentrations and their vertical distribution. Despite the significant attenuation in the intense cores, the presence of the rain maximum is easily detected by using information of Z(sub mp). In the stratiform region, the attenuation is quite limited (usually less than 5 dBZ), and the reduction of reflectivity is mostly related to the actual vertical structure of cloud distribution. Since Z(sub ms) suffers severe attenuation and tends to underestimate surface rainfall intensity in convective regions. Z(sub mp) can be more representative for rainfall retrieval in the lower atmosphere in these regions. In the stratiform region where attenuation is negligible, however, Z(sub mp) tends to overestimate surface rainfall and Z(sub ms) is more appropriate for rainfall retrieval. A hybrid technique using a weight between the two rain intensities is tested and found potentially usefull for future applications. The estimated surface rain-rate map based on this hybrid approach captures many of the details of the cloud model rain field but still slightly underestimates the rain-rate maximum.

Yeh, H.-Y. M.↗

Thermoplastic Matrix Composite Design for Cryotanks Using Multiscale Modeling and Bayesian Optimization

Designing lightweight, robust cryogenic storage tanks is critical for future launch vehicles, in-space propellant storage, and hydrogen powered aircraft. This work presents a multiscale modeling and Bayesian optimization framework for the design of thermoplastic matrix composite cryotanks. Molecular dynamics simulations are first used to determine temperature-dependent constituent properties for candidate thermoplastic matrices, which are homogenized to the lamina scale using NASA’s Multiscale Analysis Tool (NASMAT). These lamina properties, in combination with laminate family generation rules, are evaluated in HyperX structural optimization software to identify stacking sequences that meet all cryogenic load requirements. A Bayesian optimization framework is applied, with HyperX in the loop (via the HyperX API) to efficiently search across material and laminate design variables, yielding an optimized cryotank configuration with significant reductions in design cycle time compared to exhaustive search approaches.

thermoplastics↗

Structural and magnetic characterization of La{sub 0.8}Sr{sub 0.2}MnO{sub 3} nanoparticles prepared via a facile microwave-assisted method

Nanoparticles of La{sub 0.8}Sr{sub 0.2}MnO{sub 3} (LSMO) with different particle sizes are synthesized by a very fast, inexpensive, reproducible, and environmentally friendly method: the microwave irradiation of the corresponding mixture of nitrates. The structural and magnetic properties of the samples are investigated by the X-Ray diffraction (XRD), Fourier transform infra-red (FT-IR) spectroscopy, field-emission scanning electron microscopy (FE-SEM), transmission electron microscopy (TEM), and magnetic (DC magnetization and AC susceptibility) measurements. The XRD study coupled with the Rietveld refinement show that all samples crystallize in a rhombohedral structure with the space group of R−3C. The FT-IR spectroscopy and FE-SEM images indicate formation of the perovskite structure of LSMO. The DC magnetization measurements confirm the decrease in the particle size effects on the magnetic properties, e.g. reduction in the ferromagnetic (FM) moment and increase in the surface spin disorder. Magnetic dynamics of the samples studied by AC magnetic susceptibility shows that the magnetic behavior of the nanometer-sized samples is well-described by the Vogel-Fulcher and critical slowing down laws. Strong interaction between magnetic nanoparticles of LSMO was detected by fitting the experimental data with the mentioned models. - Graphical abstract: Temperature dependence of the magnetization M(T) was measured in the zero-field-cooling (ZFC) and field-cooling (FC) modes at the applied magnetic field of 100 Oe for the La{sub 0.8}Sr{sub 0.2}MnO{sub 3} with different size prepared via a facile microwave-assisted method. - Highlights: • Nanoparticles of La{sub 0.8}Sr{sub 0.2}MnO{sub 3} were synthesized by the microwave irradiation process. • The structural studies show that all samples crystallize in a rhombohedral structure with space group of R−3C. • The DC magnetic studies confirm tuning of the magnetic properties due to the particle size effects. • Magnetic dynamic studied by AC magnetic susceptibility indicate strong interaction between magnetic nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗

Prediction of grain structure after thermomechanical processing of U-10Mo alloy using sensitivity analysis and machine learning surrogate model

Abstract Hot rolling and annealing are critical intermediate steps for controlling microstructures and thickness variations when fabricating uranium alloyed with 10% molybdenum (U-10Mo), which is highly relevant to worldwide nuclear non-proliferation efforts. This work proposes a machine-learning surrogate model combined with sensitivity analysis to identify and predict U-10Mo microstructure development during thermomechanical processing. Over 200 simulations were collected using physics-based microstructure models covering a wide range of thermomechanical processing routes and initial alloy grain features. Based on the sensitivity analysis, we determined that an increase in rolling reduction percentage at each processing pass has the strongest effect in reducing the grain size. Multi-pass rolling and annealing can significantly improve recrystallization regardless of the reduction percentage. With a volume fraction below 2%, uranium carbide particles were found to have marginal effects on the average grain size and distribution. The proposed stratified stacking ensemble surrogate predicts the U-10Mo grain size with a mean square error four times smaller than a standard single deep neural network. At the same time, with a significant speedup (1000×) compared to the physics-based model, the machine learning surrogate shows good potential for U-10Mo fabrication process optimization.

36 MATERIALS SCIENCE↗

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three-Dimensional Electron Optics Model Developed for Traveling-Wave Tubes

A three-dimensional traveling-wave tube (TWT) electron beam optics model including periodic permanent magnet (PPM) focusing has been developed at the NASA Glenn Research Center at Lewis Field. This accurate model allows a TWT designer to develop a focusing structure while reducing the expensive and time-consuming task of building the TWT and hot-testing it (with the electron beam). In addition, the model allows, for the first time, an investigation of the effect on TWT operation of the important azimuthally asymmetric features of the focusing stack. The TWT is a vacuum device that amplifies signals by transferring energy from an electron beam to a radiofrequency (RF) signal. A critically important component is the focusing structure, which keeps the electron beam from diverging and intercepting the RF slow wave circuit. Such an interception can result in excessive circuit heating and decreased efficiency, whereas excessive growth in the beam diameter can lead to backward wave oscillations and premature saturation, indicating a serious reduction in tube performance. The most commonly used focusing structure is the PPM stack, which consists of a sequence of cylindrical iron pole pieces and opposite-polarity magnets. Typically, two-dimensional electron optics codes are used in the design of magnetic focusing devices. In general, these codes track the beam from the gun downstream by solving equations of motion for the electron beam in static-electric and magnetic fields in an azimuthally symmetric structure. Because these two-dimensional codes cannot adequately simulate a number of important effects, the simulation code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm) was used at Glenn to develop a three-dimensional electron optics model. First, a PPM stack was modeled in three dimensions. Then, the fields obtained using the magnetostatic solver were loaded into a particle-in-cell solver where the fully three-dimensional behavior of the beam was simulated in the magnetic focusing field. For the first time, the effects of azimuthally asymmetric designs and critical azimuthally asymmetric characteristics of the focusing stack (such as shunts, C-magnets, or magnet misalignment) on electron beam behavior have been investigated. A cutaway portion of a simulated electron beam focused by a PPM stack is illustrated.

Kory, Carol L.↗

NO THERMAL INVERSION AND A SOLAR WATER ABUNDANCE FOR THE HOT JUPITER HD 209458B FROM HST /WFC3 SPECTROSCOPY

The nature of the thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the past decade. There have been claims that many hot Jupiters exhibit atmospheric thermal inversions. However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high-precision Hubble Space Telescope WFC3 observations of the dayside thermal emission spectrum of the hot Jupiter HD 209458b, which was the first exoplanet suggested to have a thermal inversion. In contrast to previous results for this planet, our observations detect water in absorption at 6.2 σ confidence. When combined with Spitzer photometry, the data are indicative of a monotonically decreasing temperature with pressure over the range of 1–0.001 bars at 7.7 σ confidence. We test the robustness of our results by exploring a variety of model assumptions, including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06–10 × solar) and carbon-to-oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extrasolar planet atmospheres and to perform comparative exoplanetology.

79 ASTRONOMY AND ASTROPHYSICS↗

Lattice Anisotropy-Driven Reduction of Phonon Velocities in Black Phosphorus

Phonon dynamics and transport determine how heat is utilized and dissipated in materials. In 2D systems for optoelectronics and thermoelectrics, the impact of nanoscale material structure on phonon propagation is central to controlling thermal conduction. Here, we directly observe in-plane coherent acoustic phonon propagation in black phosphorus (BP) using ultrafast electron microscopy. We identify a significant reduction of the group velocities in directions intermediate to the armchair and zigzag lattice directions. Using a machine learning-based model with an >8000 atom supercell, we find that this slowing results from the mixing of in-plane transverse and longitudinal acoustic phonons and is independent of broken symmetries of edge reconstructions. In conclusion, this work demonstrates how coherent phonon transport is sensitive to propagation direction in the lattice plane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exact solutions for the static bending of Euler-Bernoulli beams using Eringen’s two-phase local/nonlocal model

Though widely used in modelling nano- and micro- structures, Eringen’s differential model shows some inconsistencies and recent study has demonstrated its differences between the integral model, which then implies the necessity of using the latter model. In this paper, an analytical study is taken to analyze static bending of nonlocal Euler-Bernoulli beams using Eringen’s two-phase local/nonlocal model. Firstly, a reduction method is proved rigorously, with which the integral equation in consideration can be reduced to a differential equation with mixed boundary value conditions. Then, the static bending problem is formulated and four types of boundary conditions with various loadings are considered. By solving the corresponding differential equations, exact solutions are obtained explicitly in all of the cases, especially for the paradoxical cantilever beam problem. Finally, asymptotic analysis of the exact solutions reveals clearly that, unlike the differential model, the integral model adopted herein has a consistent softening effect. Comparisons are also made with existing analytical and numerical results, which further shows the advantages of the analytical results obtained. Additionally, it seems that the once controversial nonlocal bar problem in the literature is well resolved by the reduction method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SCALE Code System

The SCALE Code System is a widely-used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.SCALE 6.2 provides many new capabilities and significant improvements of existing features.New capabilities include:• ENDF/B-VII.1 nuclear data libraries CE and MG with enhanced group structures,• Neutron covariance data based on ENDF/B-VII.1 and supplemented with ORNL data,• Covariance data for fission product yields and decay constants,• Stochastic uncertainty and correlation quantification for any SCALE sequence with Sampler,• Parallel calculations with KENO,• Problem-dependent temperature corrections for CE calculations,• CE shielding and criticality accident alarm system analysis with MAVRIC,• CE depletion with TRITON (T5-DEPL/T6-DEPL),• CE sensitivity/uncertainty analysis with TSUNAMI-3D,• Simplified and efficient LWR lattice physics with Polaris,• Large scale detailed spent fuel characterization with ORIGAMI and ORIGAMI Automator,• Advanced fission source convergence acceleration capabilities with Sourcerer,• Nuclear data library generation with AMPX, and• Integrated user interface with Fulcrum.Enhanced capabilities include:• Accurate and efficient CE Monte Carlo methods for eigenvalue and fixed source calculations,• Improved MG resonance self-shielding methodologies and data,• Resonance self-shielding with modernized and efficient XSProc integrated into most sequences,• Accelerated calculations with TRITON/NEWT (generally 4x faster than SCALE 6.1),• Spent fuel characterization with 1470 new reactor-specific libraries for ORIGEN,• Modernization of ORIGEN (Chebyshev Rational Approximation Method [CRAM] solver, API for high-performance depletion, new keyword input format)• Extension of the maximum mixture number to values well beyond the previous limit of 2147 to ~2 billion,• Nuclear data formats enabling the use of more than 999 energy groups,• Updated standard composition library to provide more accurate use of natural abundances, andvi• Numerous other enhancements for improved usability and stability.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Leveraging Intrinsic Principal Directions for Multifidelity Uncertainty Quantification

In this work we propose an approach for accelerating Uncertainty Quantification (UQ) analysis in the context of Multifidelity applications. In the presence of complex multiphysics applications, which often require a prohibitive computational cost for each evaluation, multifidelity UQ techniques try to accelerate the convergence of statistics by leveraging the in- formation collected from a larger number of a lower fidelity model realizations. However, at the-state-of-the-art, the performance of virtually all the multifidelity UQ techniques is related to the correlation between the high and low-fidelity models. In this work we proposed to design a multifidelity UQ framework based on the identification of independent important directions for each model. The main idea is that if the responses of each model can be represented in a common space, this latter can be shared to enhance the correlation when the samples are drawn with respect to it instead of the original variables. There are also two main additional advantages that follow from this approach. First, the models might be correlated even if their original parametrizations are chosen independently. Second, if the shared space between models has a lower dimensionality than the original spaces, the UQ analysis might benefit from a dimension reduction standpoint. In this work we designed this general framework and we also tested it on several test problems ranging from analytical functions for verification purpose, up to more challenging application problems as an aero-thermo-structural analysis and a scramjet flow analysis.

42 ENGINEERING↗

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations over the Indian Ocean

Observational analyses and a hierarchy of ocean general circulation model (OGCM) experiments were performed to understand the influence of oceanic processes on the warm sea surface temperature anomalies (SSTAs) prior to the convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden-Julian Oscillation (MJO), in the equatorial Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool region during the November-April season of the 2001-2012 period. 17/39 ISO events initiated in the Seychelles-Chagos Thermocline Ridge (SCTR) before propagating eastward; the remaining events initiated in the southern Arabian Sea (6) or Warm Pool (16) regions. The SCTR event set was notable in that it contained more global-scale MJOs (71-76%), as defined by the RMM and OMI indices, than the WP events (25-44%). Additionally, ~24% (44%) of the SCTR (Warm Pool) events were preceded by strong oceanic process-induced SSTAs of similar magnitude to those of shortwave radiative and turbulent heat fluxes. The Arabian Sea events, however, were not associated with statistically significant SSTA signals prior to convection. Based on a mixed layer heat budget analysis, entrainment and upwelling reduction were the dominant oceanic processes contributing to the warming, in contrast with boreal summer, when horizontal advection dominated. We examined several case studies, including primary MJO events, where oceanic Rossby waves were associated with the entrainment and upwelling reduction. Two simple atmospheric boundary layer convergence models revealed that the SSTAs contributed at least half of the total convergence and suggested that the ocean dynamical effect was responsible for the majority of SSTA-forced convergence for those case studies. These results underscore the need for climate prediction models to accurately represent the ocean structure and processes to include the effects of oceanic predictors.

West, B. Jason↗

Adaptive Aerostructures for Revolutionary Civil Supersonic Transportation

To enable commercially viable civil supersonic transport (SST) aircraft, innovative solutions must be developed to meet noise and efficiency requirements for overland flight. This research effort consists of a multidisciplinary team of academic and industrial experts exploring for the first time the potential of small real-time geometric outer mold line (OML) reconfigurations to minimize sonic boom signatures and aircraft drag in response to changing ambient conditions, thereby enabling noise-compliant overland supersonic flight. The team utilizes recent advances in supersonic computational fluid dynamic (CFD) methods, new noise prediction tools, and new design approaches to consider embedded highly energy-dense shape memory alloy (SMA) actuators for local shape modifications to an SST aircraft leading to optimal low boom signature and low drag in different environments. This university-led program will provide strategic leadership toward technology convergence that advances NASA's Aerospace Research Mission Directorate's (ARMD) research objectives with regard to Thrust 2: “Innovation in Commercial Supersonic Aircraft” by exploring for the first time enabling low-boom operation across a range of flight conditions via structural adaptivity, and will promote education of the next generation of engineers. The overall research strategy is to pursue three critical areas: the design of configurations for reducing boom, SMA material development and modeling, and technology feasibility demonstration in a relevant environment. Initially, the team will identify potential applications where structure or geometry adaptivity provides a benefit in noise or drag across the entire flight envelope. For selected applications/structural locations, required OML geometry changes will be determined based on analysis of sonic boom ground signature and drag reduction using new design tools, trade studies, and atmospheric sensing techniques. Designs will be developed and evaluated against requirements on loading, stroke length, and operational temperature. New alloy formulations will be developed tailored for both autonomous and controlled actuation modes. As the SMA material development matures, integrated system-level factors will be investigated. Optimized designs for small-scale distributed adaptivity applications of maximum benefit will then be matured and tested, moving toward demonstration of the innovative technology approaches at a TRL 4-5 and showing that sonic booms can be reduced by reconfiguration on demand.

Koushik Datta↗

Hybrid RANS/LES Simulation of Jet Surface Interaction Noise

NASA’s Commercial Supersonic Technology (CST) project has formulated a technical challenge to design a quiet propulsion system for a low boom supersonic aircraft that meets Federal Aviation Authority’s airport noise regulations with sufficient margin. Several proposed configurations take advantage of shielding from the wing or other air-frame components. Development of carefully validated computational tools are necessary for critically evaluating installation concepts that are currently being proposed to meet the technical challenge. Semi-empirical models that predict the noise reduction potential of arbitrary shielding surfaces are yet to mature. Another key challenge is the systematic assessment of additional noise from the interaction between high speed jet turbulence and a surface in it’s vicinity. As a first step towards predicting noise reduction due to radical installation concepts from first principles, we simulate the noise generated by a high speed turbulent round jet near a simple planar surface. Detailed comparisons are made with a dedicated experiment conducted at NASA’s Glenn Research Center. Sensitivity of far-field noise predictions to grid resolution is systematically documented. A permeable Ffowcs Williams Hawkings (FWH) surface enclosing both the jet and the shielding surface is used to predict far-field noise from the simulated flowfield. Details of the structured overset grids, numerical discretization, and turbulence model are provided. Near-field comparisons to PIV data and far-field comparisons to microphone array measurements are discussed. Excellent agreement for an initial validation study on an isolated free round jet was obtained and the findings were utilized in the jet surface interaction study. The split between shielded and reflected side of the microphone array was captured with good agreement, as well as the peak in the noise spectra due to scattering of turbulent energy into sound by the trailing edge of the surface.

Stich, Gerrit-Daniel↗

Atomic Ordering Enhanced Electrocatalytic Activity of Nanoalloys for Oxygen Reduction Reaction

For oxygen reduction reaction (ORR) over alloy electrocatalysts, the understanding of how the atomic arrangement of the metal species in the nanocatalysts is responsible for the catalytic enhancement is challenging for achieving better design and tailoring of nanoalloy catalysts. This paper reports results of an investigation of the atomic structures and the electrocatalytic activities of ternary and binary nanoalloys, aiming at revealing a fundamental insight into the unique atomic-scale structure-electrocatalytic activity relationship. PtIrCo catalyst and its binary counterparts (PtCo and PtIr) are chosen as a model system for this study. The effect of thermochemical treatment temperature on the atomic-scale structure of the catalysts was examined as a useful probe to the structure-activity correlation. The structural characterization of the binary and ternary nanoalloy catalysts was performed by combining surface sensitive techniques such as XPS and 3D atomic ordering sensitive techniques such as high-energy X-ray diffraction (HE-XRD) coupled to atomic pair distribution function (PDF) analysis (HE-XRD/PDFs) and computer simulations. The results show that the thermal treatment temperature tunes the nanoalloy’s atomic and chemical ordering in a different way depending on the chemical composition, leading to differences in the nanoalloy’s mass and specific activities. A unique structural tunability of the atomic ordering in a platinum-iridium-cobalt nanoalloy has been revealed for enhancing greatly the electrocatalytic activity toward oxygen reduction reaction, which has significant implication for rational design and nanoengineering of advanced catalysts for electrochemical energy conversion and storage.

Environmental Molecular Sciences Laboratory↗

Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO 2 reduction reaction (CO 2 RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO 2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO 2 RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO 2 RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO 2 RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at –0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Furthermore, our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

30 DIRECT ENERGY CONVERSION↗

Probing Active-Site Relocation in Cu/SSZ-13 SCR Catalysts during Hydrothermal Aging by in situ EPR Spectroscopy, Kinetic Studies, and DFT Calculations

Cu/SSZ-13 selective catalytic reduction catalysts display activity loss after mild hydrothermal aging without support degradation or loss of isolated Cu ion content. By applying electron paramagnetic resonance spectroscopy to compare hydrated, dehydrated and NH3-saturated catalysts before and after hydrothermal aging, as well as density functional theory calculations, it is found that hydrothermal aging induces Cu relocation, leading to stronger Cu-support interactions. This is particularly evidenced by ZCuOH conversion to Z2Cu, which leads to the generation of a unique structure with two Cu(II) ions positioned in a double 6-membered ring prism. The Cu(II) ions that do not undergo such dramatic changes also appear to stay closer to framework windows as evidenced by the overall increased anisotropy during EPR measurements. In situ EPR measurements under SCR conditions demonstrate that, in comparison to the fresh catalyst, higher percentages of Cu species stay as Cu(II) in the aged catalyst. This suggests that the reduction half-cycle of the SCR redox, i.e., Cu(II) ? Cu(I), becomes slower for the hydrothermally aged catalyst.

42 ENGINEERING↗