Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

A study of boundary layer transition on outgassing cones in hypersonic flow

Surface heat-transfer rates and pressures were measured at hypersonic speeds on sharp cones at zero angle of attack with and without gas injection. Using the non-injection results for reference data the effects on heating and transition location of surface roughness and injectant rate, distribution and composition were determined. The transition location was sensitive to the injectant distribution. The transition Reynolds numbers were significantly greater when the injectant distribution was constant than with a variable distribution. The measured heat-transfer distribution were also strongly dependent upon the injectant distribution. Transition Reynolds number results obtained during this program with a variable injectant distribution were correlated with a limited amount of data available for a degrading model tested in a different facility. Transitional data with constant injectant distribution were correlated with earlier results. An empirical correlation of heat-transfer reduction due to gas injection in turbulent flow was developed for both distributions tested. Several effects of mass addition on heating and transition, which have been earlier reported, were observed.

Stalmach, C. J., Jr.↗
DOE OSTI2020

ELM Suppression by Boron Powder Injection and Comparison with Lithium Powder Injection on EAST

Type I edge-localized modes (ELMs) in the Experimental Advanced Superconducting Tokamak (EAST) were completely suppressed via boron powder injection into the X-point region of an upper-single null configuration over a wide range of operating conditions (2.8 < P aux < 7.5 MW, 3.8 × 10 19 < n e < 6 × 10 19 m -3 , RF-only and RF + NBI heating scenarios, both grad-B drift directions, and even He ion majority plasmas) (Sun et al. in Nucl. Fusion, 2020). A window of edge B concentration for stable long pulse operation was identified: too low and ELMs return, too high and the discharge suffers radiative collapse. The injection of boron powder above the minimum for ELM suppression coincided with the occurrence of an edge harmonic oscillation detected in magnetics (both on the high-field side and low-field side), in AXUV diodes near the upper X-point, divertor Dα emission, and in a range of other diagnostics (Diallo et al. in: Proceedings of 2020 IAEA fusion energy conference, 2021). Furthermore, no harmonic oscillation was observed when ELMs were present, and stored energy was slightly increased at constant density during ELM suppression. Core tungsten emission during ELM suppression either increased or decreased relative to ELMy H-mode, but the W emission was maintained at acceptable levels. The threshold B injection rate was measured for several conditions, and found to increase with heating power. Li powder injection into comparable discharges also resulted in a short phase of ELM suppression, but density and stored energy both decreased due to the strong pumping effect of lithium; no edge harmonic oscillation was observed with Li injection, indicating that the ELM suppression mechanisms differ. The new set of B-seeded, ELM-suppressed discharges exhibited certain characteristics of quiescent H-mode (Burrell et al. in Phys Plasmas 8:2153, 2001), but did not require high shear, counter beams, etc. The wide operating window and compatibility with RF-only discharges paves the way for future experiments targeting long pulse H-mode discharges with complete ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Scaling of Turbulence and Microphysics in a Convection–Cloud Chamber of Varying Height

Abstract The convection–cloud chamber enables measurement of aerosol and cloud microphysics, as well as their interactions, within a turbulent environment under steady‐state conditions. Increasing the size of a convection–cloud chamber, while holding the imposed temperature difference constant, leads to increased Rayleigh, Reynolds and Nusselt numbers. Large–eddy simulation coupled with a bin microphysics model allows the influence of increased velocity, time, and spatial scales on cloud microphysical properties to be explored. Simulations of a convection–cloud chamber, with fixed aspect ratio and increasing heights of H = 1, 2, 4, and (for dry conditions only) 8 m are performed. The key findings are: Velocity fluctuations scale as H 1/3 , consistent with the Deardorff expression for convective velocity, and implying that the turbulence correlation time scales as H 2/3 . Temperature and other scalar fluctuations scale as H −3/7 . Droplet size distributions from chambers of different sizes can be matched by adjusting the total aerosol injection rate as the horizontal cross‐sectional area (i.e., as H 2 for constant aspect ratio). Injection of aerosols at a point versus distributed throughout the volume makes no difference for polluted conditions, but can lead to cloud droplet size distribution broadening in clean conditions. Cloud droplet growth by collision and coalescence leads to a broader right tail of the distribution compared to condensation growth alone, and this tail increases in magnitude and extent monotonically as the increase of chamber height. These results also have implications for scaling within turbulent, cloudy mixed‐layers in the atmosphere, such as fog layers.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Reductions in Multi-Component Jet Noise by Water Injection

An experimental investigation was performed in the NASA Langley Low Speed Aeroacoustics Wind Tunnel to determine the extent of jet exhaust noise reduction that can be obtained using water injection in a hot jet environment. The effects of water parameters such as mass flow rate, injection location, and spray patterns on suppression of dominant noise sources in both subsonic and supersonic jets were determined, and extrapolations to full-scale engine noise reduction were made. Water jets and sprays were injected in to the shear layers of cold and hot circular jets operating at both subsonic and supersonic exhaust conditions. Use of convergent-divergent and convergent nozzles (2.7in. D) allowed for simulations of all major jet noise sources. The experimental results show that water injection clearly disrupts shock noise sources within the jet plume, with large reductions in radiated shock noise. There are smaller reductions in jet mixing noise, resulting in only a small decrease in effective perceived noise level when projections are made to full scale. The fact that the measured noise reduction in the direction upstream of the nozzle was consistently larger than in the noisier downstream direction contributed to keeping effective perceived noise reductions small. Variations in the operation of the water injection system clearly show that injection at the nozzle exit rather than further downstream is required for the largest noise reduction. Noise reduction increased with water pressure as well as with its mass flow, although the type of injector had little effect.

Norum, Thomas D.↗
NASA NTRS2012

Effect of Surface Nonequilibrium Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using finite-rate gas/surface interaction model provides time-accurate solutions for multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal Response and Ablation Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas momentum conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of gas/surface interaction chemistry between air and carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was a Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kang↗
NASA NTRS2012

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗
NASA NTRS2006

Design Enhancements of the Two-Dimensional, Dual Throat Fluidic Thrust Vectoring Nozzle Concept

A Dual Throat Nozzle fluidic thrust vectoring technique that achieves higher thrust-vectoring efficiencies than other fluidic techniques, without sacrificing thrust efficiency has been developed at NASA Langley Research Center. The nozzle concept was designed with the aid of the structured-grid, Reynolds-averaged Navier-Stokes computational fluidic dynamics code PAB3D. This new concept combines the thrust efficiency of sonic-plane skewing with increased thrust-vectoring efficiencies obtained by maximizing pressure differentials in a separated cavity located downstream of the nozzle throat. By injecting secondary flow asymmetrically at the upstream minimum area, a new aerodynamic minimum area is formed downstream of the geometric minimum and the sonic line is skewed, thus vectoring the exhaust flow. The nozzle was tested in the NASA Langley Research Center Jet Exit Test Facility. Internal nozzle performance characteristics were defined for nozzle pressure ratios up to 10, with a range of secondary injection flow rates up to 10 percent of the primary flow rate. Most of the data included in this paper shows the effect of secondary injection rate at a nozzle pressure ratio of 4. The effects of modifying cavity divergence angle, convergence angle and cavity shape on internal nozzle performance were investigated, as were effects of injection geometry, hole or slot. In agreement with computationally predicted data, experimental data verified that decreasing cavity divergence angle had a negative impact and increasing cavity convergence angle had a positive impact on thrust vector angle and thrust efficiency. A curved cavity apex provided improved thrust ratios at some injection rates. However, overall nozzle performance suffered with no secondary injection. Injection holes were more efficient than the injection slot over the range of injection rates, but the slot generated larger thrust vector angles for injection rates less than 4 percent of the primary flow rate.

Flamm, Jeffrey D.↗
DOE OSTI2021

Facile Synthesis of Palladium-Based Nanocrystals with Different Crystal Phases and a Comparison of Their Catalytic Properties

A relatively unexplored aspect of noble-metal nanomaterials is polymorphism, or their ability to crystallize in different crystal phases. Here, in this work, a method is reported for the facile synthesis of Ru@Pd core–shell nanocrystals featuring polymorphism, with the core made of hexagonally close-packed (hcp)-Ru while the Pd shell takes either an hcp or face-centered cubic (fcc) phase. The polymorphism shows a dependence on the shell thickness, with shells thinner than ≈1.4 nm taking the hcp phase whereas the thicker ones revert to fcc. The injection rate provides an experimental knob for controlling the phase, with one-shot and drop-wise injection of the Pd precursor corresponding to fcc-Pd and hcp-Pd shells, respectively. When these nanocrystals are tested as catalysts toward formic acid oxidation, the Ru@Pd hcp nanocrystals outperform Ru@Pd fcc in terms of both specific activity and peak potential. Density functional theory calculations are also performed to elucidate the origin of this performance enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Humate Amendment Injection Viability Testing – Supplementary Batch Mixture and Soil Column Studies

Groundwater in the Lost Lake Aquifer Zone (LLAZ) in the Southern Sector of the M-Area Hazardous Waste management Facility (HWMF) is contaminated with chlorinated ethenes, including trichloroethylene (TCE) and tetrachloroethylene (PCE). Treatment of contaminated groundwater with humic acid is being evaluated as a potential corrective action for these volatile organic compounds (VOCs) in the LLAZ in Southern Sector (SRNS, 2019). Pilot scale injections of Huma-K brand humic acid for groundwater treatment were previously performed in M-Area between 2017 and 2020 (Amidon, 2023). Groundwater was extracted from the Lost Lake Aquifer Zone (LLAZ), then a solution of humate was mixed with recovered groundwater and subsequently re-injected into the same aquifer unit. Operations of the historical humate pilot testing were discontinued in 2020 due to COVID-19 restrictions, low injection rates, and recirculation fluid spillage. No further field scale or laboratory research to investigate humate amended groundwater injection was performed at the SRS until 2024. At the request of Area Completion Projects (ACP), the Savannah River National Laboratory (SRNL) performed supplementary testing to help determine the viability of further injection of humate as a remedial option for VOCs within the LLAZ in the Southern Sector of the M-Area HWMF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Simple Statistical Models for Predicting Overpressure Due to CO2 and Low-Salinity Waste-Fluid Injection into Deep Saline Formations

Deep saline aquifers have been used for waste-fluid disposal for decades and are the proposed targets for large-scale CO2 storage to mitigate CO2 concentration in the atmosphere. Due to relatively limited experience with CO2 injection in deep saline formations and given that the injection targets for CO2 sometimes are the same as waste-fluid disposal formations, it could be beneficial to model and compare both practices and learn from the waste-fluid disposal industry. In this paper, we model CO2 injection in the Patterson Field, which has been proposed as a site for storage of 50 Mt of industrial CO2 over 25 years. We propose general models that quickly screen the reservoir properties and calculate pressure changes near and far from the injection wellbore, accounting for variable reservoir properties. The reservoir properties we investigated were rock compressibility, injection rate, vertical-to-horizontal permeability ratio, average reservoir permeability and porosity, reservoir temperature and pressure, and the injectant total dissolved solids (TDS) in cases of waste-fluid injection. We used experimental design to select and perform simulation runs, performed a sensitivity analysis to identify the important variables on pressure build-up, and then fit a regression model to the simulation runs to obtain simple proxy models for changes in average reservoir pressure and bottomhole pressure. The CO2 injection created more pressure compared to saline waste-fluids, when similar mass was injected. However, we found a more significant pressure buildup at the caprock-reservoir interface and lower pressure buildup at the bottom of the reservoir when injecting CO2 compared with waste-fluid injection.

Ansari, Esmail↗
DOE OSTI2020

Reply to Comment by Peterie Et Al. on “Accelerated Fill‐Up of the Arbuckle Group Aquifer and Links to U.S. Midcontinent Seismicity”

Abstract Peterie et al. question one observation in our paper: associating pressure increases to injection volumes at distances of up to 25 km from an injection well. In this reply, we show that the comment misunderstands our analysis and the evidence that led to this conclusion. We also show that gauge‐depth‐corrected pressures, used by the authors to produce statewide pressure maps, are discrepant with the static fluid level data, provided in our original compilation and analysis. The discrepancies are a result of the pressure correction method employed, which naïvely substitutes formation pressure for bottomhole pressure to calculate wellbore fluid density. Their linearly interpolated pressure maps, based on sparse data, contain interpolation and extrapolation artifacts that contradict injection trends in the state, the Theis solution, and the superposition principle. We reiterate that pressure and static fluid level increases in Class I wells existed prior to 2013, most notably in central Kansas, where recent earthquakes are cited in the comment as evidence of a pressure plume emanating from the Kansas‐Oklahoma border, 90 km away. We show that the space‐time pattern of seismicity in this area is inconsistent with a northward propagating pressure plume and, instead, seismicity appears to be centered on and near a cluster of high‐rate injection wells, two of which are among the highest rate wells in the state. These observations, along with recent M4 + earthquakes during continued decreases in wastewater injection in southern Kansas and northern Oklahoma, question the usefulness of the comment for understanding and managing societally significant earthquakes.

Ansari, Esmail↗
DOE OSTI2024

The Potential of Stratospheric Aerosol Injection to Reduce the Climatic Risks of Explosive Volcanic Eruptions

Sulfur-rich volcanic eruptions happen sporadically. If Stratospheric Aerosol Injection (SAI) were to be deployed, it is likely that explosive volcanic eruptions would happen during such a deployment. Here we use an ensemble of Earth System Model simulations to show how changing the injection strategy post-eruption could be used to reduce the climate risks of a large volcanic eruption; the risks are also modified even without any change to the strategy. For a medium-size eruption (10 Tg-SO 2 ) comparable to the SAI injection rate, the volcanic-induced cooling would be reduced if it occurs under SAI, especially if artificial sulfur dioxide injections were immediately suspended. Alternatively, suspending injection only in the eruption hemisphere and continuing injection in the opposite would reduce shifts in precipitation in the tropical belt and thus mitigate eruption-induced drought. Finally, we show that for eruptions much larger than the SAI deployment, changes in SAI strategy would have minimal effect.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Global atmospheric chemistry of CFC-123

The compound 1,1-dichloro--2,2,2-trifluoroethane (CFC-123) has been proposed as an industrial substitute for trichlorofluoromethane (CFC-11). The chemical destruction rates of CFC-123 by various processes is calculated here using a three-dimensional global model of the atmosphere, and it is confirmed that the chief sink of CFC-123 is destruction by OH radicals below 12 km, accounting for 88 percent of its loss. The calculated destruction rate is greatest in the equatorial region below 2 km. The calculated steady-state lifetime of CFC-123 is 1.5 years, much shorter than that of CFC-11, the destruction of which is largely confined to the stratosphere. For equal rates of CFC-123 and CFC-11 emission to the atmosphere, the molar content in the atmosphere and the injection rate of chlorine into the stratosphere are, respectively, 48 and 14 times greater for CFC-11 than for CFC-123 in steady state.

Prinn, Ronald G.↗
DOE OSTI2023

Optimization of operational strategies for rich gas enhanced oil recovery based on a pilot test in the Bakken tight oil reservoir

Horizontal well drilling and multistage hydraulic fracturing have been demonstrated as effective approaches for stimulating oil production in the Bakken tight oil reservoir. However, after multiple years of production, primary oil recovery in the Bakken is generally less than 10% of the estimated original oil in place. Gas huff ‘n’ puff (HnP) has been tested in the Bakken Formation as an enhanced oil recovery (EOR) method; however, most field pilot test results showed no significant incremental oil production. One of the factors affecting HnP EOR performance is premature gas breakthrough, which is one of the most critical issues observed in the field because of the presence of interwell fractures. Consequently, injected gas rapidly reaches adjacent production wells without contacting reservoir rock and increasing oil recovery. Proper conformance control is therefore needed to avoid early gas breakthrough and improve EOR performance. In this study, a rich gas EOR pilot in the Bakken was carefully analyzed to collect the essential reservoir and operational data. A simulation model with 16 wells was then developed to reproduce the production history and predict the EOR performance with and without conformance control. EOR operational strategies, including single- and multiple-well HnP, with different gas injection constraints were investigated. The simulation results of single-well HnP without conformance control showed that a rich gas injection rate of at least 10 MMscfd was needed to yield meaningful incremental oil production. The strategy of conformance control via water injection could significantly improve oil production in the HnP well, but injecting an excessive amount of water also leads to water breakthrough and loss of oil production in the offset wells. By analyzing the production performance of the wells individually, the arrangement of wells was optimized for multiple-well HnP EOR. The multiwell results showed that rich gas EOR could improve oil production up to 7.4% by employing conformance control strategies. Furthermore, replacing rich gas with propane as the injection gas could result in 14% of incremental oil production.

02 PETROLEUM↗
DOE OSTI2026

Defining Infrastructure Feasibility for Hub-Scale Offshore Atlantic Carbon Storage in the Northeastern United States

In the Northeast U.S., deep rock formations along the Atlantic outer continental shelf may have the potential to sequester 150–1136 million metric tons of CO2. However, the design and infrastructure necessary to develop offshore carbon storage in this region is not well defined because there has been little oil and gas exploration and no commercial production. Consequently, an infrastructure feasibility design was completed for a hub-scale offshore CO2 storage system along the Northeast U.S. Atlantic. The design included development of a detailed, site-specific geological model for a location near the Great Stone Dome geological structure in the Baltimore Canyon Trough off the coast of Delaware, Maryland, and New Jersey. A field injection system topology design was completed to portray a design with eight wells in two clusters connected by central manifolds. Reservoir simulations were completed for the injection system that showed the hub may be able to inject 17 million metric tons (MMT) of CO2 per year for thirty years, but injection rates varied substantially across the eight wells. A CO2 pipeline design determined feasible routes from the east coast shoreline to the injection field. Finally, the CO2 injection system design included subsea injection trees, manifolds, and power umbilicals. This is the first study to define large-scale carbon storage design and infrastructure options for the offshore Atlantic, which can help to progress this region towards field characterization and early-mover deployment for future decarbonization in the region.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗