Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir simulation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Potential for monitoring soil erosion features and soil erosion modeling components from remotely sensed data

Accurate estimates of soil erosion and its effects on soil productivity are essential in agricultural decision making and planning from the field scale to the national level. Erosion models have been primarily developed for designing erosion control systems, predicting sediment yield for reservoir design, predicting sediment transport, and simulating water quality. New models proposed are more comprehensive in that the necessary components (hydrology, erosion-sedimentation, nutrient cycling, tillage, etc.) are linked in a model appropriate for studying the erosion-productivity problem. Recent developments in remote sensing systems, such as Landsat Thematic Mapper, Shuttle Imaging Radar (SIR-B), etc., can contribute significantly to the future development and operational use of these models.

Langran, K. J.↗
NASA NTRS1995

Aerothermodynamic Facilities And Measurement: Flow Characterization in Shock Tunnels

This presentation will examine the key performance aspects of shock tunnels as they relate to their use as aerothermodynamic flow simulation facilities. Assessment of shock tube reservoir conditions and flow contaminants generated in the shock tube will be presented along with their limiting impact on viable test envelopes, Facility nozzle performance as it pertains to test time assessment and nozzle exit flow quality (survey of pressure, temperature, and species) will be addressed. Also included will be a discussion of free stream flow diagnostics, both intrusive and nonintrusive, for measurement of critical flow properties not directly inferred from surface mounted transducers. The use of computational fluid dynamics for purposes of validating experimental measurements as well as predicting performance in regimes where measurements are not feasible or possible will be discussed. The use of CFD for facility research and design will also be presented.

Cavolowsky, John A.↗
NASA NTRS2023

Arctic Tropospheric Ozone: Assessment of Current Knowledge and Model Performance

As the third most important greenhouse gas (GHG) after carbon dioxide (CO2) and methane (CH4), tropospheric ozone (O3) is also an air pollutant causing damage to human health and ecosystems. This study brings together recent research on observations and modeling of tropospheric O3 in the Arctic, a rapidly warming and sensitive environment. At different locations in the Arctic, the observed surface O3 seasonal cycles are quite different. Coastal Arctic locations, for example, have a minimum in the springtime due to O3 depletion events resulting from surface bromine chemistry. In contrast, other Arctic locations have a maximum in the spring. The 12 state-of-the-art models used in this study lack the surface halogen chemistry needed to simulate coastal Arctic surface O3 depletion in the springtime; however, the multi-model median (MMM) has accurate seasonal cycles at non-coastal Arctic locations. There is a large amount of variability among models, which has been previously reported, and we show that there continues to be no convergence among models or improved accuracy in simulating tropospheric O3 and its precursor species. The MMM underestimates Arctic surface O3 by 5 % to 15 % depending on the location. The vertical distribution of tropospheric O3 is studied from recent ozonesonde measurements and the models. The models are highly variable, simulating free-tropospheric O3 within a range of ±50 % depending on the model and the altitude. The MMM performs best, within ±8 % for most locations and seasons. However, nearly all models overestimate O3 near the tropopause (∼300 hPa or ∼8 km), likely due to ongoing issues with underestimating the altitude of the tropopause and excessive downward transport of stratospheric O3 at high latitudes. For example, the MMM is biased high by about 20 % at Eureka. Observed and simulated O3 precursors (CO, NOx, and reservoir PAN) are evaluated throughout the troposphere. Models underestimate wintertime CO everywhere, likely due to a combination of underestimating CO emissions and possibly overestimating OH. Throughout the vertical profile (compared to aircraft measurements), the MMM underestimates both CO and NOx but overestimates PAN. Perhaps as a result of competing deficiencies, the MMM O3 matches the observed O3 reasonably well. Our findings suggest that despite model updates over the last decade, model results are as highly variable as ever and have not increased in accuracy for representing Arctic tropospheric O3.

Arctic↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2008

Variable Conductance Heat Pipes for Radioisotope Stirling Systems

In a Stirling radioisotope system, heat must continually be removed from the GPHS modules, to maintain the GPHS modules and surrounding insulation at acceptable temperatures. Normally, the Stirling convertor provides this cooling. If the Stirling engine stops in the current system, the insulation is designed to spoil, preventing damage to the GPHS, but also ending the mission. An alkali-metal Variable Conductance Heat Pipe (VCHP) was designed to allow multiple stops and restarts of the Stirling engine. A VCHP turns on with a delta T of 30 C, which is high enough to not risk standard ASRG operation but low enough to save most heater head life. This VCHP has a low mass, and low thermal losses for normal operation. In addition to the design, a proof-of-concept NaK VCHP was fabricated and tested. While NaK is normally not used in heat pipes, it has an advantage in that it is liquid at the reservoir operating temperature, while Na or K alone would freeze. The VCHP had two condensers, one simulating the heater head, and the other simulating the radiator. The experiments successfully demonstrated operation with the simulated heater head condenser off and on, while allowing the reservoir temperature to vary over 40 to 120 C, the maximum range expected. In agreement with previous NaK heat pipe tests, the evaporator delta T was roughly 70 C, due to distillation of the NaK in the evaporator.

Anderson, William G.↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗
DOE OSTI2023

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗
DOE OSTI2023

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Field Experiments and Model Validation: The EGS Collab Project

Implementing enhanced geothermal systems (EGS) will require improvements in understanding stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The US Department of Energy (DOE) is addressing these and other challenges at multiple scales. The EGS Collab project, addressed here, is performing tests and modeling at the 10 m scale. The FORGE project is performing tests at the full reservoir scale. The EGS Collab team created an underground testbed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently building a second testbed at SURF at a depth of about 1.25 km aimed at investigating shear stimulation (Experiment 2). In Experiment 1, we drilled eight boreholes in a well-characterized phyllite and installed geophysical sensors in six of them to create a well-instrumented testbed to allow careful monitoring of stimulation events and flow tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results both in near-real-time and in detailed process studies. Stimulations performed in this testbed allowed quantification of processes occurring during stimulation and the examination of dynamic flow occurrences. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Our second testbed, targeted at shear stimulation, is currently being built at the SURF Facility at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1.

artificial intelligence↗
DOE OSTI2020

The reservoir network: A new network topology for district heating and cooling

Thermal district networks are effective solutions to substitute fossil fuels with renewable energy sources for heating and cooling. Moreover, thermal networking of buildings allows energy efficiency to be further increased. The waste heat from cooling can be reused for heating in thermal district systems. Because of bidirectional energy flows between prosumers, thermal networks require new hydraulic concepts. In this work, we present a novel network topology for simultaneous heating and cooling: the reservoir network. The reservoir network is robust in operation due to hydraulic decoupling of transfer stations, integrates heat sources and heat sinks at various temperature levels and is flexible in terms of network expansion. Furthermore, we used Modelica simulations to compare the new single-pipe reservoir network to a basecase double-pipe network, taking yearly demand profiles of different building types for heating and cooling. The electric energy consumed by the heat pumps and circulations pumps differs between the reservoir and base case networks by less than 1%. However, if the reservoir network is operated with constant instead of variable mass flow rate, the total electrical consumption can increase by 48% compared to the base case. As with any other network topology, the design and control of such networks is crucial to achieving energy efficient operation. Investment costs for piping and trenching depend on the district layout and dimensioning of the network. If a ring layout is applied in a district, the reservoir network with its single-pipe configuration is more economical than other topologies. For a linear layout, the piping costs are slightly higher for the reservoir network than for the base case because of larger pipe diameters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1993

A 3D simulation of the early winter distribution of reactive chlorine in the north polar vortex

Early in December 1991, high values of ClO are seen by the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite at latitudes south of areas of temperatures cold enough to form polar stratospheric clouds (PSCs). A 3D simulation shows that the heterogeneous conversion of chlorine reservoirs to reactive chlorine on the surfaces of PSCs (processing) takes place at high latitudes. Often the processed air must be transported to lower latitudes, where the reactive chlorine is photochemically converted to ClO, to be observed by MLS. In this simulation, one incidence of cold temperatures is associated with an anticyclone, and a second with a cyclone. The transport of processed air associated with the anticyclone is marked by shearing; a decrease in the maximum of the processed air is accompanied by growth of the area influenced by the processing. In contrast, the air processed in the cyclonic event spreads more slowly. This shows that transport and shearing is a crucial element to the evolution of reactive chlorine associated with a processing event. In particular, transport and shearing, as well as photochemical processes, can cause variations in observed ClO.

Douglass, A.↗
NASA NTRS1983

Numerically simulated two-dimensional auroral double layers

A magnetized 2 1/2-dimensional particle-in-cell system which is periodic in one direction and bounded by reservoirs of Maxwellian plasma in the other is used to numerically simulate electrostatic plasma double layers. For the cases of both oblique and two-dimensional double layers, the present results indicate periodic instability, Debye length rather than gyroradii scaling, and low frequency electrostatic turbulence together with electron beam-excited electrostatatic electron-cyclotron waves. Estimates are given for the thickness of auroral doule layers, as well as the separations within multiple auroral arcs. Attention is given to the temporal modulation of accelerated beams, and the possibilities for ion precipitation and ion conic production by the double layer are hypothesized. Simulations which include the atmospheric backscattering of electrons imply the action of an ionospheric sheath which accelerates ionospheric ions upward.

Borovsky, J. E.↗
DOE OSTI2025

Identifying Potential Geochemical and Microbial Impacts of Hydrogen Storage in a Deep Saline Aquifer

Hydrogen is valuable commodity and a promising energy carrier for variable energy production. Storage of hydrogen may occur through injection of hydrogen or a hydrogen/methane gas blend in subsurface reservoirs. However, the geochemical and biological reactions that may impact the stored hydrogen are not yet understood. Therefore, we collected samples from a deep storage aquifer located in the St. Peter Formation in southern Illinois. The reservoir material was primarily quartz with sulphur and iron deposits, while the major constituents of the fluid were chloride and sulphate. 16S rRNA gene amplicon sequencing revealed a low biomass microbial community that contained no obvious hydrogen-consuming bacteria. Next, we enriched a field sample to increase the biomass and completed a metagenomic analysis, finding a low number of genes present that are associated with hydrogen consumption. Then, we completed a series of reactor experiments under reservoir conditions with 15% H2/85% CH4 gas simulating a short-term hydrogen storage, high withdrawal scenario. We found minimal changes in the geochemistry or microbiology for the reactor experiments. This work suggests that short-term storage may be highly successful, although significant additional work needs to be completed in order to accurately evaluate the risks associated with long-term hydrogen storage scenarios. It is essential we continue to expand our understanding of the dynamics present in saline aquifers and provide new insights into how hydrogen storage may impact underground geological storage environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗
NASA NTRS1977

Multiple-compartment venting

Program simulates venting characteristics of multireservoir multiconductor configurations. Reservoirs may be finite or infinite volume.

Donehoo, L. K.↗
NASA NTRS1999

A Comparison of Observations and Model Simulations of NO(x)/NO(y) in the Lower Stratosphere

Extensive airborne measurements of the reactive nitrogen reservoir (NOY) and its component nitric oxide (NO) have been made in the lower stratosphere. Box model. simulations that are constrained by observations of radical and longlived species and which include heterogeneous chemistry systematically underpredict the NO(x) (= NO + NO2) to NOY ratio. The model agreement is substantially improved if newly measured rate coefficients for the OH + NO2 and OH + HNO3 reactions are used. When included in 2-D models, the new rate coefficients significantly increase the calculated ozone loss due to NO(x) and modestly change the calculated ozone abundances in the lower stratosphere. Ozone changes associated with the emissions of a fleet of supersonic aircraft are also altered. phase chemistry linking NO(x) and nitric acid (HNO3), which is generally the most abundant NOY species. Because the continuous daylight present at summer high latitudes limits the heterogeneous production of HNO3 by N2O5 hydrolysis, gas-phase reactions primarily control the balance between NO(x) and NO(y). Outside summer polar regions, the N2O5 hydrolysis reaction occurring on stratospheric sulfate aerosols is a more important sink of NO(x), particularly in the lower stratosphere during winter when heterogeneous pathways account for most of the HNO3 production.

Gao, R. S.↗