Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Thin Film Sensors for Surface Measurements

Advanced thin film sensors that can provide accurate surface temperature, strain, and heat flux measurements have been developed at NASA Glenn Research Center. These sensors provide minimally intrusive characterization of advanced propulsion materials and components in hostile, high-temperature environments as well as validation of propulsion system design codes. The sensors are designed for applications on different material systems and engine components for testing in engine simulation facilities. Thin film thermocouples and strain gauges for the measurement of surface temperature and strain have been demonstrated on metals, ceramics and advanced ceramic-based composites of various component configurations. Test environments have included both air-breathing and space propulsion-based engine and burner rig environments at surface temperatures up to 1100 C and under high gas flow and pressure conditions. The technologies developed for these sensors as well as for a thin film heat flux gauge have been integrated into a single multifunctional gauge for the simultaneous real-time measurement of surface temperature, strain, and heat flux. This is the first step toward the development of smart sensors with integrated signal conditioning and high temperature electronics that would have the capability to provide feedback to the operating system in real-time. A description of the fabrication process for the thin film sensors and multifunctional gauge will be provided. In addition, the material systems on which the sensors have been demonstrated, the test facilities and the results of the tests to-date will be described. Finally, the results will be provided of the current effort to demonstrate the capabilities of the multifunctional gauge.

Martin, Lisa C.↗
NASA NTRS2015

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗
NASA NTRS1992

Effects of temperature-dependent molecular absorption coefficients on the thermal infrared remote sensing of the earth surface

The effect of temperature-dependent molecular absorption coefficients on thermal infrared spectral signatures measured from satellite sensors is investigated by comparing results from the atmospheric transmission and radiance codes LOWTRAN and MODTRAN and the accurate multiple scattering radiative transfer model ATRAD for different atmospheric profiles. The sensors considered include the operational NOAA AVHRR and two research instruments planned for NASA's Earth Observing System (EOS): MODIS-N (Moderate Resolution Imaging Spectrometer-Nadir-Mode) and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer). The difference in band transmittance is as large as 6 percent for some thermal bands within atmospheric windows and more than 30 percent near the edges of these atmospheric windows. The effect of temperature-dependent molecular absorption coefficients on satellite measurements of sea-surface temperature can exceed 0.6 K. Quantitative comparison and factor analysis indicate that more accurate measurements of molecular absorption coefficients and better radiative transfer simulation methods are needed to achieve SST accuracy of 0.3 K, as required for global numerical models of climate, and to develop land-surface temperature algorithms at the 1-K accuracy level.

Wan, Zhengming↗
NASA NTRSDate not supplied

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Mars Science Laboratory↗
NASA NTRS1997

A User''s Guide to the Zwikker-Kosten Transmission Line Code (ZKTL)

This user's guide documents updates to the Zwikker-Kosten Transmission Line Code (ZKTL). This code was developed for analyzing new liner concepts developed to provide increased sound absorption. Contiguous arrays of multi-degree-of-freedom (MDOF) liner elements serve as the model for these liner configurations, and Zwikker and Kosten's theory of sound propagation in channels is used to predict the surface impedance. Transmission matrices for the various liner elements incorporate both analytical and semi-empirical methods. This allows standard matrix techniques to be employed in the code to systematically calculate the composite impedance due to the individual liner elements. The ZKTL code consists of four independent subroutines: 1. Single channel impedance calculation - linear version (SCIC) 2. Single channel impedance calculation - nonlinear version (SCICNL) 3. Multi-channel, multi-segment, multi-layer impedance calculation - linear version (MCMSML) 4. Multi-channel, multi-segment, multi-layer impedance calculation - nonlinear version (MCMSMLNL) Detailed examples, comments, and explanations for each liner impedance computation module are included. Also contained in the guide are depictions of the interactive execution, input files and output files.

Kelly, J. J.↗
NASA NTRSDate not supplied

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗
NASA NTRS1998

Evaluation of 3D Inverse Code Using Rotor 67 as Test Case

A design modification of Rotor 67 is carried out with a full 3D inverse method. The blade camber surface is modified to produce a prescribed pressure loading distribution, with the blade tangential thickness distribution and the blade stacking line at midchord kept the same as the original Rotor 67 design. Because of the inviscid-flow assumption used in the current version of the method, Rotor 67 geometry is modified for use at a design point different from the original design value. A parametric study with the prescribed pressure loading distribution yields the following results. In the subsonic section, smooth pressure loading shapes generally produce blades with well-behaved blade surface pressure distributions. In the supersonic section, the study shows that the strength and position of the passage shock correlate with the characteristics of the blade pressure loading shape. In general, "smooth" prescribed blade pressure loading distributions generate blade designs with reverse cambers which have the effect of weakening the passage shock.

Dang, T.↗
DOE OSTI2022

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Power flow in magnetically insulated transmission lines with ion backscatter effects

Ion backscattering off of surfaces in magnetically insulated transmission lines (MITLs) is often ignored in kinetic simulations of MITL power flow. Backscattering reduces ion current losses and the surface impact heating, which dictates the rate at which surface-adsorbed contaminants are liberated into the anode–cathode gap. Backscatter probabilities are difficult to implement in a kinetic code because there are limited data for incident ion energies less than a few keV. This paper presents an analytic model based on the Rutherford scattering formula that reproduces the measured backscatter probabilities at high incident energies and transitions to the highly reflective behavior expected at low incident energies. The backscatter model is implemented in power flow simulations, which are validated with current loss experiments conducted on the 0.4 TW Mykonos accelerator at Sandia National Laboratories. This simulation setup is then used in a high-current Z-machine shot. In conclusion, backscatter effects are found to be unimportant in the low-current Mykonos regime but significantly reduce current losses at the Z-machine scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

A Boundary Condition for Simulation of Flow Over Porous Surfaces

A new boundary condition is presented.for simulating the flow over passively porous surfaces. The model builds on the prior work of R.H. Bush to eliminate the need for constructing grid within an underlying plenum, thereby simplifying the numerical modeling of passively porous flow control systems and reducing computation cost. Code experts.for two structured-grid.flow solvers, TLNS3D and CFL3D. and one unstructured solver, USM3Dns, collaborated with an experimental porosity expert to develop the model and implement it into their respective codes. Results presented,for the three codes on a slender forebody with circumferential porosity and a wing with leading-edge porosity demonstrate a good agreement with experimental data and a remarkable ability to predict the aggregate aerodynamic effects of surface porosity with a simple boundary condition.

Frink, Neal T.↗
NASA NTRS1987

The Uses and Abuses of the Acoustic Analogy in Helicopter Rotor Noise Prediction

This paper is theoretical in nature and addresses applications of the acoustic analogy in helicopter rotor noise prediction. It is argued that in many instances the acoustic analogy has not been used with care in rotor noise studies. By this it is meant that approximate or inappropriate formulations have been used. By considering various mechanisms of noise generation, such abuses are identified and the remedy is suggested. The mechanisms discussed are thickness, loading, quadrupole, and blade-vortex interaction noise. The quadrupole term of the Ffowcs Williams-Hawkings equation is written in a new form which separates the contributions of regions of high gradients such as shock surfaces. It is shown by order of magnitude studies that such regions are capable of producing noise with the same directivity as the thickness noise. The inclusion of this part of quadrupole sources in current acoustic codes is quite practical. Some of the difficulties with the use of loading noise formulations of the first author in predictions of blade-vortex interaction noise are discussed. It appears that there is a need for development of new theoretical results based on the acoustic analogy in this area. Because of the impulsive character of the blade surface pressure, a time scale of integration different from that used in loading and thickness computations must he used in a computer code for prediction of blade-vortex interaction noise.

Farassat, F.↗
NASA NTRS2021

A Practical Guide to Writing a Radiative Transfer Code

Using our decades-long experience in radiative transfer (RT) code development for Earth science, we endeavor to reduce the knowledge gap of bringing RT from theory to code quickly. Despite numerous classic and recent literature, it is still hard to develop anRT code from scratch within a few weeks. It is equally hard to understand, not to mention modify, an existing “monster” RT code, for which the developer is either located remotely or has retired. Following the format of “Numerical Recipes” by Press et al., we collocate in this paper small pieces of necessary theory with corresponding small pieces of RT code. These are arranged in an order that is natural for code development, which is often opposite of the natural order for laying out the theoretical basis. We focus on the transfer of unpolarized monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface. Both the surface and the atmosphere are homogeneous (uniform) at all directions. The multiple scattering is numerically solved using the deterministic method of Gauss-Seidel iterations. Except for the presented Python-Numba open-source RT code gsit, the paper does not report any new scientific results, but rather serves as an academic demonstration. If development time is an issue or the reader is familiar with basic concepts of RT theory, we recommend proceeding directly to Sec.3 “RT code development.

multiple light scattering↗
NASA NTRSDate not supplied

Surface and Atmospheric Shortwave Radiative Fluxes From ISCCP Data and Fu-Liou Algorithm at Higher Spatial and Vertical Resolution

NASA/GEWEX Surface Radiation Budget (SRB) has provided surface and TOA shortwave and longwave radiative fluxes using International Satellite Cloud Climatology data as the primary input. Release 4.0 came out in 2019, with data from July 1983 through June 2017, and improved validation stats relative to Release 3.0 (Stackhouse et al., 2012). Release 4.0 used the newly reprocessed ISCCP H data set. The shortwave algorithm for versions 3 and 4 has been based on the Pinker and Laszlo (1992) algorithm, which relies on the satellite-measured shortwave radiance to produce a TOA albedo. This background TOA albedo is made utilizing the background clear-sky visible radiance (called clear-sky composite radiance in the ISCCP data) and applying narrowband-to-broadband angular distribution models. Using a climatological aerosol assumption, the surface albedo is estimated using radiative transfer based (i.e., Fu/Liou radiative transfer model; Fu and Liou, 1993) look-up tables. Using this surface albedo and an instantaneous TOA albedo, the look-up tables are used to adjust both the aerosol optical depth (for clear sky) and cloud optical depth (for cloudy sky) for each scene. These properties and look up tables are then used to generate surface fluxes. While this has worked well to provide useful fluxes, it has occasionally resulted in inaccurate direct and diffuse surface SW fluxes due to the overreliance on aerosol optical depth. With increased computational power and storage now available, a more direct approach is now possible. Here we present a different approach for a future high resolution SRB Release 5.0. We replace the Pinker-Laszlo lookup table approach with a forward call to the Fu-Liou radiative transfer model as modified by the CERES team (Rose et al., 2006). In addition to being a proven radiation code, Fu-Liou allows the calculation of fluxes at different atmospheric levels and spectral bands, which will provide more insight into the surface radiation budget, its variability and attribution. It is expected that horizontal resolution will increase to at least 0.5°x0.5° from the current 1°x1°. Here we demonstrate a preliminary look at the next SRB, and show comparisons with the current release.

Earth Radiation Budget↗
NASA NTRS1988

Large-scale Advanced Prop-fan (LAP) high speed wind tunnel test report

High Speed Wind Tunnel testing of the SR-7L Large Scale Advanced Prop-Fan (LAP) is reported. The LAP is a 2.74 meter (9.0 ft) diameter, 8-bladed tractor type rated for 4475 KW (6000 SHP) at 1698 rpm. It was designated and built by Hamilton Standard under contract to the NASA Lewis Research Center. The LAP employs thin swept blades to provide efficient propulsion at flight speeds up to Mach .85. Testing was conducted in the ONERA S1-MA Atmospheric Wind Tunnel in Modane, France. The test objectives were to confirm that the LAP is free from high speed classical flutter, determine the structural and aerodynamic response to angular inflow, measure blade surface pressures (static and dynamic) and evaluate the aerodynamic performance at various blade angles, rotational speeds and Mach numbers. The measured structural and aerodynamic performance of the LAP correlated well with analytical predictions thereby providing confidence in the computer prediction codes used for the design. There were no signs of classical flutter throughout all phases of the test up to and including the 0.84 maximum Mach number achieved. Steady and unsteady blade surface pressures were successfully measured for a wide range of Mach numbers, inflow angles, rotational speeds and blade angles. No barriers were discovered that would prevent proceeding with the PTA (Prop-Fan Test Assessment) Flight Test Program scheduled for early 1987.

Campbell, William A.↗
NASA NTRS1997

Application of Chimera Navier-Stokes Code for High Speed Flows

The primary task for this year was performed in support of the "Trailblazer" project. The purpose of the task was to perform an extensive CFD study of the shock boundary-layer interaction between the engine-diverters and the primary body surfaces of the Trailblazer vehicle. Information gathered from this study would be used to determine the effectiveness of the diverters in preventing the boundary-layer coming off of the vehicle forebody from entering the main engines. The PEGSUS code was used to define the "holes" and "boundaries" for each grid. Two sets of CFD calculations were performed.Extensive post-processing of the results was performed.

Ajmani, Kumud↗
NASA NTRS1995

SeaWiFS calibration and algorithm validation

A substantial amount of field research has been completed during this period (1 October 1994 - 31 October 1995). The applicability of different in-flight calibration methods was studied. In addition, various field instruments and new equipment were tested to appraise their usefulness during such calibration activities. The purposes and (expected) results of these field experiments are detailed in separate sections of this report. One section contains a synopsis of calibrations conducted over both land and water targets. A second section contains the results of an experiment conducted to determine the spatial uniformity of surface reflectance at Lake Tahoe, the chosen water calibration site. Advances were also made in the area of radiative transfer software development. An existing exact radiative transfer code, Successive Orders (SO), was modified to provide outputs of the radiance distribution at arbitrary altitudes above a selected calibration target. These modifications were needed to utilize the SO software in radiance-based sensor calibrations. Serious thought was given to the general applicability of various calibration techniques to calibrations conducted over water surfaces. In particular, the practicality of the radiance-based and reflectance-based methods of in-flight radiometric calibration was evaluated.

Parada, Robert J., Jr.↗
NASA NTRS2000

Test Cases for the Benchmark Active Controls: Spoiler and Control Surface Oscillations and Flutter

As a portion of the Benchmark Models Program at NASA Langley, a simple generic model was developed for active controls research and was called BACT for Benchmark Active Controls Technology model. This model was based on the previously-tested Benchmark Models rectangular wing with the NACA 0012 airfoil section that was mounted on the Pitch and Plunge Apparatus (PAPA) for flutter testing. The BACT model had an upper surface spoiler, a lower surface spoiler, and a trailing edge control surface for use in flutter suppression and dynamic response excitation. Previous experience with flutter suppression indicated a need for measured control surface aerodynamics for accurate control law design. Three different types of flutter instability boundaries had also been determined for the NACA 0012/PAPA model, a classical flutter boundary, a transonic stall flutter boundary at angle of attack, and a plunge instability near M = 0.9. Therefore an extensive set of steady and control surface oscillation data was generated spanning the range of the three types of instabilities. This information was subsequently used to design control laws to suppress each flutter instability. There have been three tests of the BACT model. The objective of the first test, TDT Test 485, was to generate a data set of steady and unsteady control surface effectiveness data, and to determine the open loop dynamic characteristics of the control systems including the actuators. Unsteady pressures, loads, and transfer functions were measured. The other two tests, TDT Test 502 and TDT Test 5 18, were primarily oriented towards active controls research, but some data supplementary to the first test were obtained. Dynamic response of the flexible system to control surface excitation and open loop flutter characteristics were determined during Test 502. Loads were not measured during the last two tests. During these tests, a database of over 3000 data sets was obtained. A reasonably extensive subset of the data sets from the first two tests have been chosen for Test Cases for computational comparisons concentrating on static conditions and cases with harmonically oscillating control surfaces. Several flutter Test Cases from both tests have also been included. Some aerodynamic comparisons with the BACT data have been made using computational fluid dynamics codes at the Navier-Stokes level (and in the accompanying chapter SC). Some mechanical and active control studies have been presented. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. Cases for static angles of attack, static trailing-edge and upper-surface spoiler deflections are included for a range of conditions near those for the oscillation cases. Cases for trailing-edge control and upper-surface spoiler oscillations for a range of Mach numbers, angle of attack, and static control deflections are included. Cases for all three types of flutter instability are selected. In addition some cases are included for dynamic response measurements during forced oscillations of the controls on the flexible mount. An overview of the model and tests is given, and the standard formulary for these data is listed. Some sample data and sample results of calculations are presented. Only the static pressures and the first harmonic real and imaginary parts of the pressures are included in the data for the Test Cases, but digitized time histories have been archived. The data for the Test Cases are also available as separate electronic files.

Bennett, Robert M.↗
DOE OSTI2023

Introducing the embedded random phase approximation: H 2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. Further, we explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H 2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078–7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208–3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗