Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Viscosities of the Baryon-Rich Quark-Gluon Plasma from Beam Energy Scan Data

Here, this work presents the first Bayesian inference study of the (3 + 1)D dynamics of relativistic heavy-ion collisions and quark-gluon plasma viscosities using an event-by-event (3 + 1)D hydrodynamics + hadronic transport theoretical framework and data from the Relativistic Heavy Ion Collider Beam energy scan program. Robust constraints on initial state nuclear stopping and the baryon chemical potential-dependent shear viscosity of the produced quantum chromodynamic (QCD) matter are obtained. The specific bulk viscosity of the QCD matter is found to exhibit a preferred maximum around $\sqrt{s_{NN}}$ = 19.6 GeV . This result allows for the alternative interpretation of a reduction (and/or increase) of the speed of sound relative to that of the employed lattice-QCD based equation of state for net baryon chemical potential μ B ~0.2(0.4) GeV .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Wind Tunnel Aerodynamic Characteristics of a Transport-type Airfoil in a Simulated Heavy Rain Environment

The effects of simulated heavy rain on the aerodynamic characteristics of an NACA 64-210 airfoil section equipped with leading-and trailing-edge high-lift devices were investigated in the Langley 14- by 22-Foot Subsonic Tunnel. The model had a chord of 2.5 ft, a span of 8 ft, and was mounted on the tunnel centerline between two large endplates. Aerodynamic measurements in and out of the simulated rain environment were obtained for dynamic pressures of 30 and 50 psf and an angle-of-attack range of 0 to 20 degrees for the cruise configuration. The rain intensity was varied to produce liquid water contents ranging from 16 to 46 gm/cu m. The results obtained for various rain intensity levels and tunnel speeds showed significant losses in maximum lift capability and increases in drag for a given lift as the liquid water content was increased. The results obtained on the landing configuration also indicate a progressive decrease in the angle of attack at which maximum lift occurred and an increase in the slope of the pitching-moment curve as the liquid water content was increased. The sensitivity of test results to the effects of the water surface tension was also investigated. A chemical was introduced into the rain environment that reduced the surface tension of water by a factor of 2. The reduction in the surface tension of water did not significantly alter the level of performance losses for the landing configuration.

Bezos, Gaudy M.↗
NASA NTRS2000

Aeroacoustic Prediction Codes

This report describes work performed on Contract NAS3-27720AoI 13 as part of the NASA Advanced Subsonic Transport (AST) Noise Reduction Technology effort. Computer codes were developed to provide quantitative prediction, design, and analysis capability for several aircraft engine noise sources. The objective was to provide improved, physics-based tools for exploration of noise-reduction concepts and understanding of experimental results. Methods and codes focused on fan broadband and 'buzz saw' noise and on low-emissions combustor noise and compliment work done by other contractors under the NASA AST program to develop methods and codes for fan harmonic tone noise and jet noise. The methods and codes developed and reported herein employ a wide range of approaches, from the strictly empirical to the completely computational, with some being semiempirical analytical, and/or analytical/computational. Emphasis was on capturing the essential physics while still considering method or code utility as a practical design and analysis tool for everyday engineering use. Codes and prediction models were developed for: (1) an improved empirical correlation model for fan rotor exit flow mean and turbulence properties, for use in predicting broadband noise generated by rotor exit flow turbulence interaction with downstream stator vanes: (2) fan broadband noise models for rotor and stator/turbulence interaction sources including 3D effects, noncompact-source effects. directivity modeling, and extensions to the rotor supersonic tip-speed regime; (3) fan multiple-pure-tone in-duct sound pressure prediction methodology based on computational fluid dynamics (CFD) analysis; and (4) low-emissions combustor prediction methodology and computer code based on CFD and actuator disk theory. In addition. the relative importance of dipole and quadrupole source mechanisms was studied using direct CFD source computation for a simple cascadeigust interaction problem, and an empirical combustor-noise correlation model was developed from engine acoustic test results. This work provided several insights on potential approaches to reducing aircraft engine noise. Code development is described in this report, and those insights are discussed.

Gliebe, P↗
NASA NTRS2001

A Wind-Tunnel Parametric Investigation of Tiltrotor Whirl-Flutter Stability Boundaries

A wind-tunnel investigation of tiltrotor whirl-flutter stability boundaries has been conducted on a 1/5-size semispan tiltrotor model known as the Wing and Rotor Aeroelastic Test System (WRATS) in the NASA-Langley Transonic Dynamics Tunnel as part of a joint NASA/Army/Bell Helicopter Textron, Inc. (BHTI) research program. The model was first developed by BHTI as part of the JVX (V-22) research and development program in the 1980's and was recently modified to incorporate a hydraulically-actuated swashplate control system for use in active controls research. The modifications have changed the model's pylon mass properties sufficiently to warrant testing to re-establish its baseline stability boundaries. A parametric investigation of the effect of rotor design variables on stability was also conducted. The model was tested in both the on-downstop and off-downstop configurations, at cruise flight and hover rotor rotational speeds, and in both air and heavy gas (R-134a) test mediums. Heavy gas testing was conducted to quantify Mach number compressibility effects on tiltrotor stability. Experimental baseline stability boundaries in air are presented with comparisons to results from parametric variations of rotor pitch-flap coupling and control system stiffness. Increasing the rotor pitch-flap coupling (delta(sub 3) more negative) was found to have a destabilizing effect on stability, while a reduction in control system stiffness was found to have little effect on whirl-flutter stability. Results indicate that testing in R-134a, and thus matching full-scale tip Mach number, has a destabilizing effect, which demonstrates that whirl-flutter stability boundaries in air are unconservative.

Piatak, David J.↗
DOE OSTI2026

Technical, economic, and load-following capabilities assessment of grid-connected geothermal and geothermal-solar hybrid systems

The technical and economic performance as well as the load-following capabilities of grid-connected geothermal hybrid systems were assessed in this work. The analyzed geothermal hybrid configuration is composed of a binary geothermal plant integrated with a concentrating solar-thermal system and underground thermal energy storage (UTES) through a primary heat exchanger. Physics-based models for the hybrid system for plant generation capacities of 1, 25, and 50 MW were developed from validated models for each subsystem. Also, an economic model was developed that accounts for different hybrid system capabilities, solar field sizes, and thermal storage duration. The advantage of the geothermal hybrid system was assessed by comparing the performance with the baseline benchmark geothermal plant with a similar configuration and generation capacity. It was found that hybridizing geothermal plants with concentrating solar and thermal energy storage not only improves the thermal efficiency by up to 8 percentage points when additional heat from the solar-UTES loop rises the evaporator temperatures from 70 to 125 °C, but also enhances the load-following capability for the geothermal plant, which can meet a typical residential load profile with a power rate of change 0.25 kW/s with an absolute error under 13 kW for a 1 MW plant. Other benefits of hybridization include resource preservation and a potential LCOE reduction of up to 56% for a 50 MW geothermal hybrid plant having a 50% solar share, a 1.4 solar multiple, and 24-h storage capacity. The results presented in this work demonstrate that hybridizing geothermal systems transforms them into a flexible and cost-effective solution for addressing the dynamic requirements of modern electric grids.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Much of zero emissions commitment occurs before reaching net zero emissions

Abstract We explore the response of the Earth’s coupled climate and carbon system to an idealized sequential addition and removal of CO 2 to the atmosphere, following a symmetric and continuous emissions pathway, in contrast to the discontinuous emissions pathways that have largely informed our understanding of the climate response to net zero and net negative emissions to date. We find, using both an Earth system model and an ensemble of simple climate model realizations, that warming during the emissions reduction and negative emissions phases is defined by a combination of a proportionality of warming to cumulative emissions characterized by the transient climate response to emissions (TCRE), and a deviation from that proportionality that is governed by the zero emissions commitment (ZEC). About half of the ZEC is realized before reaching zero emissions, and the ZEC thus also controls the timing between peak cumulative CO 2 emissions and peak temperature, such that peak temperature may occur before peak cumulative emissions if ZEC is negative, underscoring the importance of ZEC in climate policies aimed to limit peak warming. Thus we argue that ZEC is better defined as the committed warming relative to the expected TCRE proportionality, rather than as the additional committed warming that will occur after reaching net zero CO 2 emissions. Once established, the combined TCRE and ZEC relationship holds almost to complete removal of prior cumulative CO 2 emissions. As cumulative CO 2 emissions approach zero through negative CO 2 emissions, CO 2 concentrations drop below preindustrial values, while residual long-term climate change continues, governed by multicentennial dynamical processes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Capturing Complex Multivariate Time Series Interactions to Detect High-Risk Adverse Events During Flight

The reduction of aviation safety metrics below target thresholds continue to drive down the number of aviation fatalities and accidents. To meet future safety demands, sustained efforts by aviation agencies promoting safety assurance processes and systems have prompted ongoing research on identifying and mitigating in-flight risks. With the projected increase in passenger load factor and rollout of more autonomous systems into the national airspace, the need to detect high-risk events in-time or ahead-of-time is becoming increasingly crucial. New anomaly detection and precursor identification algorithms will need to scale to different airframes, levels of autonomy, and system complexity. While the pervasiveness of deep learning has resulted in the development of performant anomaly detection methods, these sophisticated models currently suffer from low end-user interpretability. Building off our previous work on identifying adverse events in multivariate flight data during descent, we propose a data-driven approach for detecting in-flight adverse events caused by the complex interplay of flight variables. Our approach utilizes ordinal patterns of important aircraft stability variables (e.g., airspeed and descent rate) to capture multivariate flight dynamics that can be used to predict the onset of unstable approaches, a high-risk adverse event that can occur during approach. Through the use of ordinal patterns, we aim to create more interpretable detection models of in-flight adverse events that can be translated to future autonomous systems without difficulty. Our analysis shows the presence of distinct ordinal pattern distributions that can be used to predict unstable approaches 1 minute ahead of time with an accuracy of 0.69 and a recall of 0.73 and 30 seconds ahead with an accuracy of 0.70 and a recall of 0.86.

Risk detection↗
NASA NTRSDate not supplied

Capturing Multivariate Time Series Interactions to Detect High‑Risk Instability During Approach

The reduction of aviation safety metrics below target thresholds continue to drive down the number of aviation fatalities and accidents. To meet future safety demands, sustained efforts by aviation agencies promoting safety assurance processes and systems have prompted ongoing research on identifying and mitigating in-flight risks. With the projected increase in passenger load factor and rollout of more autonomous systems into the national airspace, the need to detect high-risk events in-time or ahead-of-time is becoming increasingly crucial. New anomaly detection and precursor identification algorithms will need to scale to different airframes, levels of autonomy, and system complexity. While the pervasiveness of deep learning has resulted in the development of performant anomaly detection methods, these sophisticated models currently suffer from low end-user interpretability. Building off our previous work on identifying adverse events in multivariate flight data during descent, we propose a data-driven approach for detecting in-flight adverse events caused by the complex interplay of flight variables. Our approach utilizes ordinal patterns of important aircraft stability variables (e.g., airspeed and descent rate) to capture multivariate flight dynamics that can be used to predict the onset of unstable approaches, a high-risk adverse event that can occur during approach. Through the use of ordinal patterns, we aim to create more interpretable detection models of in-flight adverse events that can be translated to future autonomous systems without difficulty. Our analysis shows the presence of distinct ordinal pattern distributions that can be used to predict unstable approaches 1 minute ahead of time with an accuracy of 0.69 and a recall of 0.73 and 30 seconds ahead with an accuracy of 0.70 and a recall of 0.86.

Risk detection↗
DOE OSTI2020

Unconventional ideas for ionization cooling of muons

Small muon beams increase the luminosity of a muon collider. Reducing the momentum and position spreads of muons reduces emittance and leads to small, cool beams. Ionization cooling has been observed at the Muon Ionization Cooling Experiment. 6D emittance reduction by a factor of 100, 000 has been achieved in simulation. Another factor of 5 in cooling would meet the basic requirements of a high luminosity muon collider. In this paper we compare, for the first time, the amount of RF needed in a cooling channel to previous linacs. Here, we also outline three methods aimed to help achieve a final factor of 5 in 6D cooling.

43 PARTICLE ACCELERATORS↗
NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗
NASA NTRS2000

The Interaction Between Dynamics and Chemistry of Ozone in the Set-Up Phase of the Northern Hemisphere Polar Vortex

Understanding stratospheric ozone loss to the point of accurately predicting ozone in the future requires correctly distinguishing chemical from transport-induced changes in ozone. For example, evaluating the impact of chlorine reduction in controlling stratospheric ozone loss requires estimating the amount of ozone lost in chemical reactions in the polar winter and spring. The Northern Hemisphere winter polar region is a particularly crucial and interesting area because it appears that the Northern vortex may currently be poised near the threshold of extreme ozone destruction such as that which now occurs annually in the Antarctic ozone "hole." In this presentation we explore the interaction of ozone transport and chemistry through the Northern late summer and fall seasons as the vortex circulation becomes established. This phase of the seasonal cycle determines the starting point for heterogeneous processes and chlorine-driven loss that take control in the winter vortex. Using a combination of profile data from POAM, HALOE, and in situ measurements, we show that relatively low ozone at high latitudes in the middle stratosphere is associated with vortex airmasses and that these ozone abundances evolve photochemically from characteristically higher values at the end of the summer. The zonal variance of ozone mixing ratio also increases greatly at this time consistent with increasing wave-driven transport. Comparison with a three-dimensional chemistry-transport model is used to generalize the findings from the limited set of observations and quantify the relative roles of transport and chemistry in determining the ozone mixing ratio distributions.

Kawa, S.R.↗
DOE OSTI2023

diffReplication - An Energy-Aware Fault Tolerance Model for Silent Error Detection and Mitigation in Heterogeneous Extreme-scale Computing Environment

At extreme scale, the frequency of silent errors – a class of errors that remain undetected by low-level error detection mechanisms – increases significantly with the computational complexity of the application and the scale of the computing infrastructure. As hardware and software advances are made to usher in the next scientific era of computing, developing new approaches to mitigate the impact of silent errors remains a challenging problem. In this work, we propose an energy-aware fault-tolerance model, referred to diffReplication to overcome silent errors. In the proposed model, the main process is associated with one replica that executes at the same rate as the main process, and one diffReplica that is executed at a fraction of the main process' execution rate. If the main and its replica reach consensus at the end of a computation phase, the state of the diffReplica is updated and computation is resumed. If the synchronization attempt results in a disagreement, however, the diffReplica increases its execution speed to complete the computation and quickly reach the synchronization barrier. Assuming a single error over any given synchronization interval, a majority voting is used to reach consensus and tolerate silent errors. To further enhance its performance, diffReplication is augmented with speculative execution, whereby the main or its fast replica is selected to continue execution without waiting for the diffReplica. The selection process is based on the previous behaviour of the main and its replica. A performance analysis study is carried out to assess the performance of diffReplication, in terms of the energy saving and time-to-completion reduction achieved by the diffReplication scheme. The experiment shows that speculative execution reduces the time to completion with additional energy, and dynamic decision-making balances the energy consumption and time to completion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗
NASA NTRS2002

Launch Vehicle/Carrier Interaction, Improving the Analytical Integration Process

A goal of the aerospace industry is to reduce the cost of space transportation by a significant within the next decade. The present cost of launching a space transportation system which includes propulsion system, vehicle, carrier, and payload integrated together to form a system, encompasses much more than the design of the propulsion system and vehicle. The total cost includes the recurring cost of the process of integrating carriers, and payloads into the vehicle for each flight. The recurring cost of the integration of carrier/payloads systems is driven by the interaction of the vehicle. If the interaction can be well characterized and made to be very predictable for a range of payloads, or if it can be minimized then the cost of integrating a payload can be reduced significantly from today's levels. The Space Shuttle is very interactive with the payload. The interaction has been well characterized through finite element modeling and is reasonably predictable for a specific payload. Experience has shown, however, that the interaction is very manifest dependent, and small changes in one portion of a payload complement can change the interaction significantly in another portion. That is the affects of one on the other are such that if one or the other is changed slightly the environment at the interfaces can change significantly. To date the Shuttle has made in excess of one hundred flights. For each of these flights several iterations of dynamic analyses have been required in the development of each vehicle/carrier/payload system. The iterative analyses are needed because of the sensitivity of the interaction of the launch vehicle to the attached carrier/payload. The Multi Purpose Logistics Module (MPLM) is a carrier designed for flight in the Space Shuttle carrying a wide variation of cargo, supplies, and experiments to and from Space Station. Its integration process provides a unique area for improvement in the template in use today for transporting items to space. Discussion of the present MPLM integration analysis requirements and possible areas for improvement of the process are provided in the subject paper. Special emphasis is placed on reduction of requirements based on similarity of items being transported, and the use of carrier hardware for protection of the Shuttle Orbiter.

Shariett, Charles A.↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2019

Acoustic Analysis of Sensor Ports

Acquiring data in rocket engines that is representative of the actual dynamic environment can often be difficult due to a multitude of influences. One source of contamination that is often not considered entirely is the response associated with the acoustic cavity created by a sensor offset. It is common to offset a sensor due to various reasons such as for mounting accessibility, thermal isolation, shock reduction, or prevention of debris impingement. While estimating the natural frequency of the acoustic cavity is straightforward, limited analysis has been described on the determination of the overall frequency response. The sensor port design approach usually attempts to ensure the port is short enough such that the acoustic response is negligible near the frequency-of-interest, but this requires knowledge of the frequency response and simple rules-of-thumb are not always guaranteed. Data correction and/or data interpretation is also often desired for an unsatisfactory response. The limited response analysis in the literature only offers approximations or neglects important contributions. A new approach is devised theoretically and computationally that captures the true acoustic response of a sensor port. This paper summarizes the acoustics background, the port response theoretical development, and provides comparisons of a port acoustic response using an analytical model and computational acoustics. The effects of nonlinear acoustics are also examined. Additionally, the paper summarizes the design of a specialized filter using the predicted sensor port response that can be applied to data for correction.

Casiano, M. J.↗
DOE OSTI2021

An out-of-distribution-aware autoencoder model for reduced chemical kinetics

While detailed chemical kinetic models have been successful in representing rates of chemical reactions in continuum scale computational fluid dynamics (CFD) simulations, applying the models in simulations for engineering device conditions is computationally prohibitive. To reduce the cost, data-driven methods, e.g., autoencoders, have been used to construct reduced chemical kinetic models for CFD simulations. Despite their success, data-driven methods rely heavily on training data sets and can be unreliable when used in out-of-distribution (OOD) regions (i.e., when extrapolating outside of the training set). In this paper, we present an enhanced autoencoder model for combustion chemical kinetics with uncertainty quantification to enable the detection of model usage in OOD regions, and thereby creating an OOD-aware autoencoder model that contributes to more robust CFD simulations of reacting flows.Here, we first demonstrate the effectiveness of the method in OOD detection in two well-known datasets, MNIST and Fashion-MNIST, in comparison with the deep ensemble method, and then present the OOD-aware autoencoder for reduced chemistry model in syngas combustion.

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

Development of a Transferable Density-Functional Tight-Binding Model for Organic Molecules at the Water/Platinum Interface

A computationally efficient and transferable approach for modeling reactions at metal/water interfaces could significantly accelerate our understanding and ultimately the development of new catalytic transformations, particularly in the context of the emerging field of biomass conversion. Here, we present a parametrization of Pt–X (X = H, O, C) density-functional tight-binding (DFTB) for addressing this need. We first constructed Pt–H, Pt–O, and Pt–C repulsive potential splines. These pairwise parameters were then augmented to include many-body interactions using the Chebyshev Interaction Model for Efficient Simulation (ChIMES). We compare the geometrical and energetic performances of both DFTB and DFTB/ChIMES methods with DFT reference data across a variety of organic molecules at the platinum surface from nanoparticles to single-crystal surfaces. DFTB shows limited transferability between extended crystal surfaces and small nanoparticles. This transferability is significantly improved through the introduction of three-body interactions with Pt in DFTB/ChIMES, which provides consistent results across various systems, with reductions in the RMSD from around 30 kcal/mol in DFTB to around 10 kcal/mol. We demonstrate the stability and reliability of the obtained parameters by performing metadynamic simulations for the adsorption of phenol on Pt(111). We observe that DFTB itself is undersolvating the surface, leading to only one or two chemisorbed water molecules in a c(4 × 6) unit cell. In contrast, DFTB/ChIMES leads to a coverage of about 0.5 ML and successfully captures the chemisorbed mode of phenol at both the solid/liquid and the solid/gas interfaces. Furthermore, in agreement with experimental measurements, the adsorption at the solid/liquid interface is significantly weaker than that at the solid/gas interface. As a result, we highlight that even with DFTB, where we can accumulate dynamics for more than 1 ns for a given system, the simulations are not fully converged.

Adsorption↗