Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Controlling morphology in hybrid isotropic/patchy particle assemblies

Brownian dynamics is used to study self-assembly in a hybrid system of isotropic particles (IPs), combined with anisotropic building blocks that represent special “designer particles.” Those are modeled as spherical patchy particles (PPs) with binding only allowed between their patches and IPs. In this study, two types of PPs are considered: Octahedral PPs (Oh-PPs) and Square PPs (Sq-PPs), with octahedral and square arrangements of patches, respectively. The self-assembly is additionally facilitated by the simulated annealing procedure. Here, the resultant structures are characterized by a combination of local correlations in cubatic ordering and a symmetry-specific variation of bond orientation order parameters (SymBOPs). By varying the PP/IP size ratio, we detected a sharp crossover between two distinct morphologies in both types of systems. High symmetry phases, NaCl crystal for Oh-PP and square lattice for Sq-PP, are observed for larger size ratios. For the smaller ones, the dominant morphologies are significantly different, e.g., Oh-PPs form a compact amorphous structure with predominantly face-to-face orientation of neighboring PPs. Unusually, for a morphology without a long-range order, it is still possible to identify well organized coherent clusters of this structure, thanks to the adoption of our SymBOP-based characterization.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Characteristics of the sub-micron ash aerosol generated during oxy-coal combustion at atmospheric and elevated pressures

This paper is concerned with the effect of pressure on the particle size distribution and the size-segregated composition of the sub-micron ash aerosol created during oxy-coal combustion under near practical self-sustaining combustion conditions. The problem is important because pressurized oxy-coal combustion has been proposed as one promising technology to minimize CO 2 emissions. Sub-micron ash plays a major role in ash deposition mechanisms, which, in turn, can control boiler performance. In this work, the same bituminous coal was burned at pressures of 1, 8 and 15 bar in O 2 /CO 2 environments. Here, tests employed a 100 kW (rated) oxy-fuel combustor (OFC) operated at atmospheric pressure (1 bar) and a 300 kW (rated) entrained-flow pressurized reactor (EFPR) at elevated pressures (8 and 15 bar). Although these tests were conducted under near practical combustion conditions, confounding effects of peak flame temperature variations were minimized for the 1 bar and 15 bar tests, allowing the role of elevated pressure to be isolated. For the EFPR tests, a specially designed sampling system was used to sample sub-micron ash aerosols from the pressurized combustor and is described in detail. Results showed that at the same peak temperature but higher pressure, the fractions of ash aerosol partitioned into the PM 0.6 and PM 1 size fractions were greatly diminished. Moreover, elevated pressures caused significant changes in the composition of the (size-segregated) sub-micron aerosol, especially in its alkali content, which increased significantly. Examination of fractions of each element that ended up in the sub-micron fume suggested that, at constant temperature, the effect of pressure on vaporization of semi-volatile metals was very different from that on the release into the gas phase of non-volatile metals and could not be explained by equilibrium.

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

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗
DOE OSTI2022

Modeling Vapor Transport Deposition of Metal-halide Perovskite Thin Films for Photovoltaic and Optoelectronic Devices

Over the past decade, metal halide perovskites (MHPs) have emerged as a promising materials platform for high-efficiency solar cells and low-cost optoelectronics. However, there are challenges that frustrate the large-scale manufacturing of MHP devices, including difficulty in controlling film composition, interface formation and their device instability under ambient conditions. Vapor processing offers an attractive path to manufacturability, while also opening the door to new opportunities in device design that could favorably impact ultrahigh efficiency tandem solar cells or overall stability. Emphasis here is on the use of an alternate processing methodology, vapor transport deposition (VTD), to deposit the MHP layer. In experimental work, we found that our deposition system had significant run-to-run variations in film thickness and composition. To understand and resolve these issues, we use COMSOL Multiphysics software to model the precursor deposition rate numerically, trying to guide the broad range of parameters in the system. The result can be matched with experimental data and provides insight into system hydrodynamics and molar transport effects. Finally, we show that chamber pressure is a key factor to scale up the VTD technique.

Hsu, Wan-Ju↗
DOE OSTI2025

A Review of Quantum Computing Technologies in Power System Optimization

As modern power grids increasingly integrate variable renewable generation, distributed energy resources, and energy storage systems, classical optimization techniques are facing unprecedented challenges. This review examines the emerging application of quantum computing to overcome these challenges in power system optimization, including optimal power flow (OPF), unit commitment (UC), economic dispatch (ED), and intelligent switching and topology optimization (IS-TO). Recent research has introduced various quantum methodologies—such as gate-based, annealing-based, variational algorithms, and quantum-inspired algorithms—to address the combinatorial complexity inherent in grid reconfiguration and energy management. The review summaries the quantum algorithms, quantum devices and the power system test cases, highlighting hybrid quantum–classical strategies that leverage the complementary strengths of both paradigms. Some quantum advantages have been observed, including theoretical speedup, accurate simulation results, scalable qubit usage, efficient QUBO mapping. In particular, the review emphasizes the importance of integrating quantum optimization techniques with classical control frameworks, these hybrid approaches demonstrate the potential to improve real-time grid management and operational reliability. A significant portion of the analysis is devoted to the practical limitations of current quantum devices. Present-day quantum hardware, operating in the noisy intermediate-scale quantum (NISQ) era, remains highly sensitive to noise and limited in qubit connectivity, which constrains the scale and accuracy of implemented algorithms. The review delves into specific challenges such as the need for qubit-efficient encoding techniques and error mitigation strategies that are critical for handling real-world grid optimization problems. In addition, the work draws attention to the performance discrepancies between theoretical quantum speedups and experimental validations, underscoring the importance of rigorous benchmark studies using representative power grid test cases. In summary, this review highlights both the promise and limitations of quantum computing for power system optimization. It provides a comprehensive overview of the state-of-the-art technologies, categorizes recent advancements in algorithm design, and discusses practical considerations for implementation, and serves as an informative resource on current research. Future research directions include developing robust hybrid frameworks, advancing qubit-efficient formulations, and scaling up experimental demonstrations to confirm the theoretical advantages of quantum methods in large-scale power system operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Modeling Vapor Transport Deposition of Metal-halide Perovskite Thin Films for Photovoltaic and Optoelectronic Devices

Over the past decade, metal halide perovskites (MHPs) have emerged as a promising materials platform for high-efficiency solar cells and low-cost optoelectronics. However, there are challenges that frustrate the large-scale manufacturing of MHP devices, including difficulty in controlling film composition, interface formation and their device instability under ambient conditions. Vapor processing offers an attractive path to manufacturability, while also opening the door to new opportunities in device design that could favorably impact ultrahigh efficiency tandem solar cells or overall stability. Emphasis here is on the use of an alternate processing methodology, vapor transport deposition (VTD), to deposit the MHP layer. In experimental work, we found that our deposition system had significant run-to-run variations in film thickness and composition. To understand and resolve these issues, we use COMSOL Multiphysics software to model the precursor deposition rate numerically, trying to guide the broad range of parameters in the system. The result can be matched with experimental data and provides insight into system hydrodynamics and molar transport effects. Finally, we show that chamber pressure is a key factor to scale up the VTD technique.

Hsu, Wan-Ju↗
DOE OSTI2022

Nanoscale Mapping and Defect-Assisted Manipulation of Surface Plasmon Resonances in 2D Bi 2 Te 3 /Sb 2 Te 3 In-Plane Heterostructures

Here, the Bi 2 Te 3 /Sb 2 Te 3 in-plane heterostructure is reported as a low-dimensional tunable chalcogenide well suited as plasmonic building block for the visible-UV spectral range. Electron-driven plasmon excitations of low-dimensional Bi 2 Te 3 /Sb 2 Te 3 are investigated by monochromated electron energy loss spectroscopy spectrum imaging. To resolve the nanoscale spatial distribution of various local plasmonic resonances, singular value decomposition is used to disentangle the spectral data and identify the individual spectral contributions of various corner, edge, and face modes. Furthermore, defect-plasmon interactions are investigated both for nanoscale intrinsic and thermally induced extrinsic polygonal defects (in situ sublimation). Signature of defect-induced red shift ranging from a several hundreds of millielectronvolts to a few electronvolts, broadening of various plasmon response, together with selective enhancement and significant variations in their intensity are detected. This study highlights the presence of a heterointerface and identifies defects as physical tuning pathways to modulate the plasmonic response over a broad spectral range. Finally, the experimental observations are compared qualitatively and validated with numerical simulations using the electron-driven discrete dipole approximation. Low-dimensional Bi 2 Te 3 /Sb 2 Te 3 as a less explored plasmonic system holds great promises as emerging platform for integrated plasmonics. Furthermore, introducing controlled structural defects can open the door for nanoengineering of plasmonic properties in such systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying the Value of Technology and Policy Innovation in Water Resource Portfolios

Abstract Recent water infrastructure planning models demonstrate that explicitly accounting for hydrological changes in long‐term water planning reduces costs and increases system robustness. However, current models do not consider the effects of other changes occurring in the system over long time horizons, namely technology and policy innovation. By leveraging state‐of‐the‐art dynamic control techniques, we propose a suite of tools that explicitly relate hydrologic change, technology innovation, and policy interventions to system‐level water costs. We design robust and cost‐optimal strategies for water supply expansion under hydrological variability and hundreds of plausible innovation states, for example, treatment cost and diversification improvements, deployment times, and demand reduction campaigns. Our approach circumvents the need to pre‐assign probabilities to innovation scenarios and, instead, directly calculates the effects of technology, policy, and hydrology variations on system cost and planning strategies to identify avenues of impactful innovation. These tools can support technology development hubs, urban, state, and federal water planners in anticipating cost implications of technology and policy innovation at the local level, and prioritizing high‐impact innovation goals based on quantitative targets. Results for the case study of Santa Barbara, California identify high‐value innovation attributes for the city as well as combinations of innovation attributes across technologies and policies, and quantify their potential utility cost reductions.

13 HYDRO ENERGY↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Advances in Kriging-Based Autonomous X-Ray Scattering Experiments

Autonomous experimentation is an emerging paradigm for scientific discovery, wherein measurement instruments are augmented with decision-making algorithms, allowing them to autonomously explore parameter spaces of interest. We have recently demonstrated a generalized approach to autonomous experimental control, based on generating a surrogate model to interpolate experimental data, and a corresponding uncertainty model, which are computed using a Gaussian process regression known as ordinary Kriging (OK). We demonstrated the successful application of this method to exploring materials science problems using x-ray scattering measurements at a synchrotron beamline. Here, we report several improvements to this methodology that overcome limitations of traditional Kriging methods. The variogram underlying OK is global and thus insensitive to local data variation. We augment the Kriging variance with model-based measures, for instance providing local sensitivity by including the gradient of the surrogate model. As with most statistical regression methods, OK minimizes the number of measurements required to achieve a particular model quality. However, in practice this may not be the most stringent experimental constraint; e.g. the goal may instead be to minimize experiment duration or material usage. We define an adaptive cost function, allowing the autonomous method to balance information gain against measured experimental cost. We provide synthetic and experimental demonstrations, validating that this improved algorithm yields more efficient autonomous data collection.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the dual advantage of placing observations through forward sensitivity analysis

The four-dimensional variational data assimilation methodology for assimilating noisy observations into a deterministic model has been the workhorse in forecasting centres for over three decades. While this method provides a computationally efficient framework for dynamic data assimilation, it is largely silent on the important question concerning the minimum number and placement of observations. To answer this question, we demonstrate the dual advantage of placing the observations where the square of the sensitivity of the model solution with respect to the unknown control variables, called forward sensitivities, attains its maximum. By following this approach, we can force the observability Gramian to be of full rank, which in turn guarantees efficient recovery of the optimal values of the control variables, which is the first of the two advantages of this strategy. We further show that the proposed strategy of placing observations has another inherent optimality: the square of the sensitivity of the optimal estimates of the control with respect to the observations (used to obtain these estimates) attains its minimum value, a second advantage that is a direct consequence of the above strategy for placing observations. Our analytical framework and numerical experiments on linear and nonlinear systems confirm the effectiveness of our proposed strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗
DOE OSTI2025

Optimizing Grid Regulation With Gravity Storage Systems: A Comparative Analysis With Different Motor Inertias

The integration of renewable energy sources into power grids necessitates solutions for grid support and stability during fluctuations in electricity generation and demand. Gravity energy storage systems (GESS) are emerging as a promising technology for managing the balance between energy supply and demand. However, their capacity to optimize energy flow and offer voltage and frequency regulation amid imbalances in generation and demand is less reported. This paper investigates the control of GESS for optimizing energy flow during voltage and frequency regulation. The study evaluates the regulation capabilities of GESS with different motor inertias during a Texas grid event: one with a high-speed, low-inertia motor and another with a low-speed, high-inertia motor. Results indicate that both GESS scenarios provide fast frequency response by converting potential energy into kinetic energy and vice versa, with a response time of 1.5 s from the frequency variation. This aligns with grid requirements for primary frequency response from traditional synchronous generators and motors with large inertia. Furthermore, a GESS based on a high-inertia motor may be able to operate over a broader range of frequency variations, whereas a low-inertia system may be limited by thermal constraints of the motor.

frequency regulation↗
DOE OSTI2024

Optimizing Grid Regulation With Gravity Storage Systems: A Comparative Analysis With Different Motor Inertias: Preprint

The integration of renewable energy sources into power grids necessitates solutions for grid support and stability during fluctuations in electricity generation and demand. Gravity energy storage systems (GESS) are emerging as a promising technology for managing the balance between the energy supply and demand. However, their capacity to optimize energy flow and offer voltage and frequency regulation amid generation-demand imbalances is less reported. This paper investigates the control of GESS for optimizing energy flow during voltage and frequency regulation. The study evaluates the regulation capabilities of GESS with different motor inertias during a Texas grid event: one with a high-speed, low-inertia motor and another with a low-speed, high-inertia motor. Results indicate that both GESS systems provide fast frequency response by converting potential energy into kinetic energy and vice versa, with a response time of 1.5 seconds from the frequency variation. This aligns with grid requirements for primary frequency response from traditional synchronous generators and motors with large inertia. Furthermore, a GESS with a high-inertia motor may be able operate over a broader range of frequency variations whereas a low-inertia system may be limited by thermal constraints of the motor.

frequency regulation↗
DOE OSTI2024

Optimizing Grid Regulation with Gravity Energy Storage Systems: A Comparative Analysis with Different Motor Inertias

The integration of renewable energy sources into power grids necessitates solutions for grid support and stability during fluctuations in electricity generation and demand. Gravity energy storage systems (GESS) are emerging as a promising technology for managing the balance between energy supply and demand. However, their capacity to optimize energy flow and offer voltage and frequency regulation amid imbalances in generation and demand is less reported. This paper investigates the control of GESS for optimizing energy flow during voltage and frequency regulation. The study evaluates the regulation capabilities of GESS with different motor inertias during a Texas grid event: one with a high-speed, low-inertia motor and another with a low-speed, high-inertia motor. Results indicate that both GESS scenarios provide fast frequency response by converting potential energy into kinetic energy and vice versa, with a response time of 1.5 s from the frequency variation. This aligns with grid requirements for primary frequency response from traditional synchronous generators and motors with large inertia. Furthermore, a GESS based on a highinertia motor may be able to operate over a broader range of frequency variations, whereas a low-inertia system may be limited by thermal constraints of the motor.

ENERGY STORAGE,POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Enabling Realistic Communications Evaluations for ADMS

Increasing penetration levels of inverter-based distributed energy resources (DERs) in distribution systems are creating challenges in system operation. Currently, wired or proprietary wireless communications are used by advanced distribution management systems (ADMS) to detect the state of the system and control the DERs for optimal operation. Wireless communications are an alternative solution that can enable connectivity between the dispersed assets used in distribution system applications. As demand for the seamless integration of distributed energy resources (DERs) continues to increase, proper networking and communications systems are becoming important to support the DER infrastructure. This report summarizes validation of a private wireless LTE network in a laboratory environment. The purpose of this work is to 1) demonstrate the effectiveness of private wireless LTE communication for grid applications and 2) provide a means of evaluating latency and overall performance of wireless communication in grid applications. First, private wireless LTE communication was validated in a direct transfer trip relaying scenario as an example of one-way communications using grid equipment. This study provided a baseline understanding of the variation in latency for seven different kinds of wireless scenarios shown in Figure A below. As expected, weak (attenuated) wireless signal-based communication yielded the widest ranges in latency. When higher priority was set for the trip signal communication, this resulted in lower latency in communication timing. In the final phase of this project, the private wireless LTE network was validated using two-way communications in an advanced distribution management system (ADMS) testbed which was developed in the ESIF at NREL. The purpose of this was to study the impact of wireless communication on grid performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Relating land surface, water table, and weathering fronts with a conceptual valve model for headwater catchments

Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the earth-air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re-orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes toward an “ideal hill”. In this steady state, reaction-front valves partition water volumes into shallow and deep flowpaths. Furthermore, these flows dissolve low- and high-solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.

54 ENVIRONMENTAL SCIENCES↗