Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Grid-Forming Frequency Shaping Control for Low-Inertia Power Systems

As power systems transit to a state of high renewable penetration, little or no presence of synchronous generators makes the prerequisite of well-regulated frequency for grid-following inverters unrealistic. Thus, there is a trend to resort to grid-forming inverters which set frequency directly. We propose a novel grid-forming frequency shaping control that is able to shape the aggregate system frequency dynamics into a first-order one with the desired steady-state frequency deviation and Rate of Change of Frequency (RoCoF) after a sudden power imbalance. The no overshoot property resulting from the first-order dynamics allows the system frequency to monotonically move towards its new steady-state without experiencing frequency Nadir, which largely improves frequency security. We prove that our grid-forming frequency-shaping control renders the system internally stable under mild assumptions. The performance of the proposed control is verified via numerical simulations on a modified Icelandic Power Network test case.

control systems↗
DOE OSTI2020

Grid-Forming Frequency Shaping Control for Low-Inertia Power Systems

As power systems transit to a state of high renewable penetration, little or no presence of synchronous generators makes the prerequisite of well-regulated frequency for grid-following inverters unrealistic. As such, there is a trend to resort to grid-forming inverters which set frequency directly. We propose a novel grid-forming frequency shaping control that is able to shape the aggregate system frequency dynamics into a first-order one with the desired steady-state frequency deviation and Rate of Change of Frequency (RoCoF) after a sudden power imbalance. The no overshoot property resulting from the first-order dynamics allows the system frequency to monotonically move towards its new steady-state without experiencing frequency Nadir, which largely improves frequency security. We prove that our grid-forming frequency-shaping control renders the system internally stable under mild assumptions. The performance of the proposed control is verified via numerical simulations on a modified Icelandic Power Network test case.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Investigating shear flow through continuum gyrokinetic simulations of limiter biasing in the Texas Helimak

Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To explore this, the first gyrokinetic simulations of limiter biasing in the Helimak using the plasma physics code Gkeyll have been carried out, and results are presented here. For the scenarios considered, turbulence is mostly driven by the interchange instability, which depends on gradients of steady state density profiles. An analysis of both experimental and simulation data demonstrates that shear rates are mostly less than than local linear growth rates, and not all requirements for shear stabilization are met. Yet, the mostly vertical shear flow has an important effect on bulk transport and experimental steady-state density profiles, and changes to the gradients correspond to changes in turbulence levels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Denoising and Extension of Response Functions in the Time Domain

Response functions of quantum systems, such as electron Green’s functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Further, response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Field Validation of Wake Steering Control with Wind Direction Variability: Preprint

Wake steering is a type of wind farm control strategy wherein upstream turbines are misaligned with the wind direction to redirect their wakes away from downstream turbines, increasing overall wind plant power. Wake steering is often analyzed assuming steady mean wind directions across the wind farm. However, in practice, the wind direction varies considerably over time because of large-scale weather phenomena. Wind direction variability causes wake steering controllers to increase power production less than predicted by steady state models, but more robust wake steering strategies can be designed that account for variable wind conditions. This paper compares the achieved yaw offsets and power gains from two 2-turbine wake steering experiments at a commercial wind farm with model predictions using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented model, assuming both fixed and variable wind directions. The impact of wind direction variability is modeled by including wind direction and yaw uncertainty in the FLORIS calculations. The field results match the trends predicted assuming wind direction variability. Specifically, the yaw offsets achieved in the intended control regions are lower than desired, resulting in less power gain, while a slight loss in power occurs for wind directions outside of the intended control region because of unintentional yaw misalignment. The agreement between the model and field results suggests that the wind direction variability model can be used to design wake steering controllers that are more robust to variable wind conditions present in the field.

17 WIND ENERGY↗
DOE OSTI2024

Newton-Okounkov bodies of chemical reaction systems

Despite their noted potential in polynomial-system solving, there are few concrete examples of Newton-Okounkov bodies arising from applications. Accordingly, in this paper, we introduce a new application of Newton-Okounkov body theory to the study of chemical reaction networks and compute several examples. An important invariant of a chemical reaction network is its maximum number of positive steady states Here, we introduce a new upper bound on this number, namely the ‘Newton-Okounkov body bound’ of a chemical reaction network. Further, through explicit examples, we show that the Newton-Okounkov body bound of a network gives a good upper bound on its maximum number of positive steady states. We also compare this Newton-Okounkov body bound to a related upper bound, namely the mixed volume of a chemical reaction network, and find that it often achieves better bounds.

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

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Testing the DIII-D co/counter off-axis neutral beam injected power and ability to balance injected torque

DIII-D has undergone a major upgrade and successfully injected high power off-axis neutral beams (~4 MW) in both co-current and counter-current directions. This capability of high power co/counter steerable off-axis neutral beams on a major tokamak opens a unique parameter space of broad pressure and current profiles for high beta steady-state advanced tokamak (AT) scenarios, while retaining the ability to balance the injected torque for low rotation studies. This co/counter off-axis neutral beam capability is being used to validate physics-based energetic particle and thermal transport models for designing next-step facilities based on the steady-state AT approach. This paper reports on the critical evaluation of the transmitted power and energetic ion population produced by this heating and current drive system, which is assessed through visible imaging, neutron measurements and rotation profile measurements at balanced torque. Minimal losses of neutral beam power have been achieved by optimizing the strongly focused ion sources required to pass through the aperture. Tilting of the ion source has been guided by fast visible imaging and resulted in neutral beam injection along the design centerline with empirical characterization of each beam's divergence derived from the imaging data and used in the NUBEAM description of the beam injection. Through exclusive power injection of each neutral beam into MHD quiescent plasmas across a range of neutral beam voltage, perveance and plasma current we conclude that a modest reduction (~10%–15%) of transmitted power compared to on-axis, standard focus has been incurred. We report corrections that more accurately represent the injected power. Good ability to balance the neutral beam torque has been demonstrated by injecting the new off-axis counter injecting beam against the existing off-axis co-injected beam in 2.0 T, 1.0 MA, MHD quiescent L-mode plasmas. Furthermore, the torque balance studies verify the ability to operate with balanced injection, which is critical for achieving low torque and low rotation operation for physics studies and in ITER demonstration discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Experimental and Numerical Analysis of Additively Manufactured Inconel 718 Coupons With Lattice Structure

Abstract In the present paper, two lattice geometries suitable for near surface and double wall cooling were developed and tested. The first type of unit cell consisted of six ligaments of 0.5 mm diameter joined at a common vertex near the middle. The second type of unit cell was derived from the first type by adding four mutually perpendicular ligaments in the middle plane. Two lattice configurations, referred to as L1 and L2, respectively, were obtained by repeating the corresponding unit cell in streamwise and spanwise directions in an inline fashion. Test coupons consisting of these lattice geometries embedded inside rectangular cooling channel with dimensions of 2.54 mm height, 38.07 mm width, and 38.1 mm in length were fabricated using Inconel 718 powder and selective laser sintering (SLS) process. The heat transfer and pressure drop performance was then evaluated using steady-state tests with constant wall temperature boundary condition and for channel Reynolds number ranging from 2800 to 15,000. The lattices depicted a higher heat transfer compared with a smooth channel and both the heat transfer and pressure drop increased with a decrease in the porosity from L1 to L2. Steady-state conjugate numerical results revealed formation of prominent vortical structures in the inter-unit cell spaces, which diverted the flow toward the top end wall and created an asymmetric heat transfer between the two end walls. In conclusion, these lattice structures provided an augmented heat transfer while favorably redistributing the coolant within channel.

Engineering↗
DOE OSTI2023

Development of a Pre-Combustion CO 2 Capture Process Using High-Temperature PBI Hollow-Fiber Membranes

The overall objective of this project was to evaluate the advantages of transformational polybenzimidazole (PBI) polymer hollow-fiber membrane (HFM)-based, carbon dioxide (CO 2 ) capture and purification technology at bench-scale using an actual coal-derived syngas stream from a coal gasification facility. The project was carried out over two budget periods. The technical objectives in Budget Period 1 (BP1) included preparing HFs and modules and upgrading the available skid for field testing. The technical objectives for BP2 were to field-test the skid unit with actual coal-derived syngas from an oxygen-blown gasifier to obtain performance data, update the Techno-Economic Analysis (TEA) that would assist with future process scale-up, and provide information on the design of a small pilot-scale test unit. The goal was to advance the PBI-HFM CO 2 capture and gas separation system for pre-combustion applications beyond second-generation economic performance predictions and make progress toward meeting overall fossil energy performance goals of CO 2 capture with 95% CO 2 purity at a cost of electricity (COE) 30% less than baseline capture approaches. The research program was designed with progressive technical tasks leading to both dynamic and steady-state testing of the PBI-HFM skid with actual coal-derived syngas. The work plan was to: (1) fabricate sufficient Generation-2 (GEN-2) fibers for module fabrication; (2) upgrade the fiber skid to accommodate large fiber modules for bench-scale field testing; (3) conduct dynamic and steady-state testing with coal-derived syngas from an oxygen-blown gasifier and obtain system performance data; (4) perform a TEA and environmental, health, and safety (EH&S) assessment; (5) update the State-Point Data Table, Technology Gap Analysis (TGA), and Technology Maturation Plan (TMP); (6) uninstall and return the test skid to the Recipient’s facilities; and (7) submit a Final Report that describes the results and analysis of the project research effort.

03 NATURAL GAS↗
DOE OSTI2026

FAST-1.2.2: Integral Assessment

This integral assessment quantifies the predictive capabilities of Fuel Analysis under Steady-state and Transients (FAST), a thermal-mechanical nuclear fuel performance code designed to analyze fuel behavior from beginning of life to burnup levels allowed by the U.S. Nuclear Regulatory Commission. FAST code calculations are shown to compare satisfactorily to a pre-selected set of experimental data with both steady-state and anticipated operating occurrence conditions and design basis accident transient operating conditions. This document describes the assessment of FAST-1.2.2, the latest version of FAST, released April 2026.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Field Validation of Wake Steering Control with Wind Direction Variability

Wake steering is a wind farm control strategy wherein upstream turbines are misaligned with the wind direction to redirect their wakes away from downstream turbines, increasing overall wind plant power. Wake steering is often analyzed assuming steady mean wind directions across the wind farm. However, in practice, the wind direction varies considerably over time because of large-scale weather phenomena. Wind direction variability causes the increase in power production from wake-steering to be less than predicted by steady-state models, but more robust wake-steering strategies can be designed that account for variable wind conditions. This paper compares the achieved yaw offsets and power gains from two different 2-turbine wake-steering experiments at a commercial wind farm with model predictions using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented model, assuming both fixed and variable wind directions. The impact of wind direction variability is modeled by including wind direction and yaw uncertainty in the FLORIS calculations. The field results match the trends predicted, assuming wind direction variability. Specifically, the yaw offsets achieved in the intended control regions are lower than desired, resulting in less power gain, while a slight loss in power occurs for wind directions outside of the intended control region because of unintentional yaw misalignment. The agreement between the model and field results suggests that the wind direction variability model can be used to design wake-steering controllers that are more robust to variable wind conditions present in the field.

17 WIND ENERGY↗
DOE OSTI2026

Bounded-Confidence Models of Multidimensional Opinions with Topic-Weighted Discordance

People’s opinions on a wide range of topics often evolve over time through their interactions with others. Models of opinion dynamics primarily focus on one-dimensional opinions, which represent opinions on one topic. However, opinions on various topics are rarely isolated; instead, they can be interdependent and correlated. In a bounded-confidence model (BCM) of opinion dynamics, agents are receptive to each other only if their opinions are sufficiently similar. Here, we extend classical agent-based BCMs—namely, the Hegselmann–Krause BCM, which has synchronous interactions, and the Deffuant–Weisbuch BCM, which has asynchronous interactions—to a multidimensional setting, in which the opinions are multidimensional vectors representing opinions of different topics and opinions on different topics are interdependent. To measure opinion differences between agents, we introduce topic-weighted discordance functions that account for opinion differences in all topics. We define regions of receptiveness for our models, and we use them to characterize the steady-state opinion clusters and provide an analytical approach to compute these regions. In addition, we numerically simulate our models on various networks with initial opinions drawn from a variety of distributions. When initial opinions are correlated across different topics, our topic-weighted BCMs yield significantly different results in both transient and steady states compared to baseline models, where the dynamics of each opinion topic are independent.

Mathematics and Computing↗
DOE OSTI2023

Multiphysics Analysis of the MSRE Experiment Using Griffin-SAM Coupled Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time since the fuel salt is circulating through the whole primary loop of the reactor. This movement of the fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and corresponding redistribution in the active core region. To capture this phenomenon accurately, neutronics and thermal hydraulics computational tools need to be modified to handle the movement of the DNPs and their decay. Idaho and Argonne National Laboratories are actively working on developing neutronics and thermal hydraulics tools to model MSRs with flowing fuel. The neutronics code Griffin and the thermal hydraulics code SAM were extended to handle flowing fuel with the drift of the DNPs under the MOOSE framework. In this work, a simplified Multiphysics model of the Molten Salt Reactor Experiment (MSRE) is developed and utilized to perform steady-state and transient analyses. The MSRE experiments were designed to show the impact of the DNPs losses on core reactivity. The following section provides a description of the MSRE experiment with the core main parameters used to develop the model. Then, the developed Multiphysics model of the MSRE is presented. Finally, steady-state verification tests and unprotected loss of flow transient tests are discussed followed by a summary and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Exploring the Feasibility of INCONEL® ALLOY 740H® for Power Plant Headers: Integrating Machine Learning with Computational Fluid Dynamics (CFD)

This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Physics-based compact modeling of electro-thermal memristors: Negative differential resistance, local activity, and non-local dynamical bifurcations

Leon Chua's Local Activity theory quantitatively relates the compact model of an isolated nonlinear circuit element, such as a memristor, to its potential for desired dynamical behaviors when externally coupled to passive elements in a circuit. However, the theory's use has often been limited to potentially unphysical toy models and analyses of small-signal linear circuits containing pseudo-elements (resistors, capacitors, and inductors), which provide little insight into required physical, material, and device properties. Furthermore, the Local Activity concept relies on a local analysis and must be complemented by examining dynamical behavior far away from the steady-states of a circuit. In this work, we review and study a class of generic and extended one-dimensional electro-thermal memristors (i.e., temperature is the sole state variable), re-framing the analysis in terms of physically motivated definitions and visualizations to derive intuitive compact models and simulate their dynamical behavior in terms of experimentally measurable properties, such as electrical and thermal conductance and capacitance and their derivatives with respect to voltage and temperature. Within this unified framework, we connect steady-state phenomena, such as negative differential resistance, and dynamical behaviors, such as instability, oscillations, and bifurcations, through a set of dimensionless nonlinearity parameters. In particular, we reveal that the reactance associated with electro-thermal memristors is the result of a phase shift between oscillating current and voltage induced by the dynamical delay and coupling between the electrical and thermal variables. We thus, demonstrate both the utility and limitations of local analyses to understand non-local dynamical behavior. Critically for future experimentation, the analyses show that external coupling of a memristor to impedances within modern sourcing and measurement instruments can dominate the response of the total circuit, making it impossible to characterize the response of an uncoupled circuit element for which a compact model is desired. However, these effects can be minimized by proper understanding of the Local Activity theory to design and utilize purpose-built instruments.

Brown, Timothy D. (ORCID:0000000347068355)↗
DOE OSTI2025

Investigation of N 2 /O 2 plasma interaction with Pt-catalyst: effect of metastable adsorbates on product hysteresis

The coupling of catalysts and atmospheric-pressure plasma has the potential to improve the efficiency of certain catalytic reactions. Understanding the changes that the catalyst surface undergoes during exposure to plasma is key to improving plasma–catalytic performance. In this work, long term exposure of Pt–Al 2 O 3 powder catalyst to an Ar/N 2 /O 2 non-equilibrium atmospheric-pressure plasma-jet was investigated. Products produced by the interaction were analyzed downstream with Fourier-transform infrared spectroscopy while surface species were analyzed operandi with diffuse reflectance infrared Fourier transform spectroscopy. During exposure, the catalyst temperature was ramped cyclically between 100 °C and 350 °C to understand how substrate temperature affects the plasma–catalyst interaction. Long-lasting changes were revealed to take place on the catalyst surface during plasma exposure. At low temperatures, Pt–O and Pt–NO accumulate on the surface which react at elevated temperatures to form NO 2 . NO 2 initially appears to spill on to the Al 2 O 3 support as nitrites and nitrates instead of desorbing. Stable surface conditions are only achieved after prolonged plasma exposure, when nitrate sites on the Al 2 O 3 support are filled. By changing the catalyst temperature at various rates, the impact of total plasma species flux to the surface was analyzed. It was found that decreasing the heating rate increased the hysteresis in the pattern of NO 2 formation during thermal cycling. The variation with temperature demonstrates that plasma exposure results in a buildup of surface NO x and oxygen species which react or desorb at high temperatures. The observed changes are discussed from the generic viewpoint that a non-equilibrium plasma interacting with a catalyst at low temperature introduces metastable steady-state surface conditions. Upon heating above a threshold temperature, the introduced surface modifications can change either due to thermal effects, or, for a plasma environment, by additional interaction with the incident plasma species flux. The surface/material changes take place in a highly predictable fashion and after sufficient time above the threshold temperature reach a steady-state condition that is different from the transient behavior that is observed during initial heating. During cooling the plasma-surface interaction exhibits a different behavior than during heating, and this results in hysteresis of diverse observables. The metastability/hysteresis description appears quite generic and analogous to hysteresis behavior seen for different systems. Furthermore, it is expected to be useful for understanding the consequences of plasma–catalyst surface interactions for various systems.

36 MATERIALS SCIENCE↗