SEARCH · Engineering Papers

Results for “Conditional Normalizing Flow”

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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2011

On the Shear Band Direction in Metallic Glasses

This paper shows that it is inappropriate to relate the angle between loading axis and shear-band (or fracture) plane in metallic glasses under uniaxial loading conditions to the coefficient of internal friction in the Mohr-Coulomb model. Shear bands in metallic glasses are a result of material instability (which can be predicted from constitutive parameters and loading conditions), but not a yield phenomenon. Specifically, the shear band directions depend on Poisson s ratio, the ratios of three deviatoric principal stresses to Mises stress, the coefficient of internal friction, and the dilatancy factor. The last parameter describes whether the plastic flow is associative or non-associative. Theoretical predictions based on the classic Rudnicki-Rice model agree well with a compilation of observations in uniaxial mechanical tests. Furthermore, we identify three (two) regimes under the two-dimensional cylindrical (three-dimensional spherical) contact where different shear-band directions may occur. When using bonded interface technique to visualize the shear bands under the contact, it should be noted that the stress component normal to the bonded interface is released, resulting into the commonly observed semi-circular shear bands whose directions are predicted to follow the larger in-plane principal stress.

Gao, Yanfei↗

DOE OSTI2023

Development, calibration, and validation of a novel gray-box energy model for residential split air conditioners

Energy models for vapor compression refrigeration in residential air conditioners have been developed through white-box, gray-box, and black-box methods in decades. However, existing white-box and gray-box models require complicated equations with detailed geometries while black-box models require substantial experimental data. Further, this paper aims to develop and validate a simple gray-box steady-state energy model without the need for detailed geometries, which can accurately predict the cooling capacity and electrical power input based on outdoor and indoor air conditions, and supply airflow rates. First three state variables, including the evaporation and condensation temperatures, and refrigerant mass flowrate, are applied to develop the energy model and are solved by three physical equations, including the energy conservations at the evaporator and condenser, and the refrigerant volume-mass flow correlation of the compressor. Secondly, seven performance property functions related to three state variables and three physical equations are identified and calibrated by simple temperature and power measurements. Finally, field experiments are conducted on a residential air conditioner to calibrate these performance property functions and validate the developed model. The validated results reveal the model can accurately predict the cooling capacity and electrical power input, with the normalized root mean square errors of 2.3% and 0.87% respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Thermal oxidation of nuclear graphite and pyrolytic carbon coatings

The oxidation of pyrolytic carbon (PyC) deposited via fluidized bed chemical vapor deposition was characterized and compared with that of standard nuclear-grade graphite. The materials were heated at 700 to 1000 °C in a thermogravimetric analysis system under 20% v/v O 2 flow, allowing for direct comparison of dynamic oxidative mass change in each material. Further, three different PyC samples fabricated under different conditions exhibited variation in total mass loss and mass loss rate, varying by as much as 709 mg/cm 2 in total mass loss and 14.2 (mg/cm 2 )/min in mass loss rate at a single temperature. These variations highlight the correlation between PyC microstructure/defect density and oxidation susceptibility. Additionally, changes in the microstructure and composition between PyC and graphite were characterized via scanning electron microscopy and correlated to the mass loss results. The results of this work have implications toward the safety of tristructural isotropic (TRISO) and other coated particle fuels, especially under off-normal conditions, given the limited information that exists about the oxidation behavior of PyC.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Rapid, Untargeted Chemical Profiling of Single Cells in Their Native Environment

We report a method that enables untargeted, high throughput, and quantitative mass spectrometric analysis of single cells from cell suspension without needing additional sample preparation procedures (e.g., molecular tagging) through the combination of single-cell printer technology and liquid vortex capture–mass spectrometry (SCP-LVC-MS). The operating principle behind the SCP-LVC-MS technology is single cell isolation via small droplet piezoelectric ejection followed by capture of the droplet into an LVC-MS sampling probe. Once exposed to an appropriate solvent, the cell is lysed, extracted, and analyzed by MS. The SCP-LVC-MS approach was validated by measuring the lipid composition of microalgae, Chlamydomonas reinhardtii (ChRe) and Euglena gracilis (EuGr), and HeLa cells in their native growth media. Numerous diacylglyceryltrimethylhomo-Ser (DGTS), phosphatidylcholine (PC), monogalactosyldiacylglycerol (MGDG), and digalactosyldiacylglycerol (DGDG) lipids were observed in single cells. Continuous solvent flow ensures that cells are analyzed rapidly, and no signal carryover between cells is observed. ChRe and EuGr microalgae mixed together in the same solution were differentiated cell-by-cell in real-time based on differences between levels of diacylglyceryltrimethylhomo-Ser (DGTS) and phosphatidylcholine (PC) lipids measured in each cell. Several DGTS lipids present in ChRe were quantified with single-cell resolution by normalizing to a DGTS(32:0) internal standard added to the LVC probe solvent during analysis. Quantitative peak areas were validated by comparing to bulk lipid extracts. Furthermore, peak area distributions comprised of hundreds of cells were compared for ChRe after 5 days of nitrogen-limited and normal growth conditions, which show clear differences and the ability to resolve cellular population differences with single-cell resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Induction Heating for Tritium Storage Beds (FY19 PDRD Project Summary Report)

Metal hydride beds used for storing hydrogen isotopes at the Savannah River Site are conductively heated and cooled to desorb and absorb gases during storage cycling. While several storage bed configurations have been utilized over the years, many achieve heating by utilizing electric cartridge heaters positioned laterally along the longest axis of the storage bed (FACE Gen 2 Beds). Others achieve heating and cooling through hot and cold nitrogen flow (Gen 1 Beds). While these beds/heating configurations are well-designed for storage cycling operations, it is unclear if they can be regenerated to operating conditions once too much helium has accumulated in them. Full regeneration of hydride beds occurs around 625 - 650 °C. It is unclear if the currently employed electric cartridge heaters in Gen 2 beds can reach these temperatures. Furthermore, since swelling can occur during routine heating operations, it is unclear if existing cartridge heaters in Gen 2 beds can be swapped for higher-power cartridges that can reach 625 - 650 °C . Similarly, nitrogen air flow (to the best of our knowledge) cannot cause Gen 1 beds to reach these temperatures. Thus, existing Gen 1 and Gen 2 beds need to be heated through an alternative method to removal helium, so the beds can either return to normal operation or be disposed of. The goal of this PDRD was to investigate if induction heating can be used for various heating processes involving hydride storage beds (e.g. regeneration, gas cycling, helium heel removal). In the original call for proposals, this PDRD fell under "Gas Processing" area of interest # 3: "Develop methods to deploy alternate ways to heat existing process equipment to minimize impact to glove box environment."

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

DOE OSTI2016

A computational study of radiation and gravity effect on temperature and soot formation in a methane air co-flow diffusion flame

An effort has been made for a quantitative assessment of the soot formed under steady state in a methane air co flow diffusion flame by a numerical simulation at normal gravity and at lower gravity levels of 0.5 G, 0.1 G and 0.0001 G (microgravity). The peak temperature at microgravity is reduced by about 50 K than that at normal gravity level. There is an augmentation of soot formation at lower gravity levels. Peak value at microgravity multiplies by a factor of ∼7 of that at normal gravity. However, if radiation is not considered, soot formation is found to be much more.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

The effect of guide-field and boundary conditions on collisionless magnetic reconnection in a stressed X-point collapse

Works of Tsiklauri and Haruki [Phys. Plasmas 15, 102902 (2008); 14, 112905 (2007)] are extended by inclusion of the out-of-plane magnetic (guide) field. In particular, magnetic reconnection during collisionless, stressed X-point collapse for varying out-of-plane guide-fields is studied using a kinetic, 2.5D, fully electromagnetic, relativistic particle-in-cell numerical code. For zero guide-field, cases for both open and closed boundary conditions are investigated, where magnetic flux and particles are lost and conserved, respectively. It is found that reconnection rates, out-of-plane currents and density in the X-point increase more rapidly and peak sooner in the closed boundary case, but higher values are reached in the open boundary case. The normalized reconnection rate is fast: 0.10-0.25. In the open boundary case it is shown that an increase of guide-field yields later onsets in the reconnection peak rates, while in the closed boundary case initial peak rates occur sooner but are suppressed. The reconnection current changes similarly with increasing guide-field; however for low guide-fields the reconnection current increases, giving an optimal value for the guide-field between 0.1 and 0.2 times the in-plane field in both cases. Also, in the open boundary case, it is found that for guide-fields of the order of the in-plane magnetic field, the generation of electron vortices occurs. Possible causes of the vortex generation, based on the flow of decoupled particles in the diffusion region and localized plasma heating, are discussed. Before peak reconnection onset, oscillations in the out-of-plane electric field at the X-point are found, ranging in frequency from approximately 1 to 2 ω{sub pe} and coinciding with oscillatory reconnection. These oscillations are found to be part of a larger wave pattern in the simulation domain. Mapping the out-of-plane electric field along the central lines of the domain over time and applying a 2D Fourier transform reveal that the waves predominantly correspond to the ordinary and the extraordinary mode and hence may correspond to observable radio waves such as solar radio burst fine structure spikes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Unsteady Land-Sea Breeze Circulations in the Presence of a Synoptic Pressure Forcing

Unsteady land-sea breezes (LSBs) that result from time-varying surface temperature contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, M g , oriented from sea to land (α = 0°) or land to sea (α = 180°). Large eddy simulations reveal the development of four distinctive regimes, depending on the joint interaction between M g , α, and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the velocity between the unsteady and the corresponding steady scenarios at Δθ = Δθ max , are found to be significant and to extend 2 hr longer for α = 0° compared to α = 180°. These diurnal dynamics result in nonequilibrium conditions that are significantly affected by the flow history, and that behave differently over the two patches for the different α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land to sea (α = 180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the turbulence-mean flow equilibrium response. Finally, the vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution and can thus be used to identify these regimes using k-means clustering.

58 GEOSCIENCES↗

DOE OSTI2015

Rayleigh-Taylor instability at spherical interfaces between viscous fluids: Fluid/vacuum interface

For a spherical interface of radius R separating two different homogeneous regions of incompressible viscous fluids under the action of a radially directed acceleration, we perform a linear stability analysis in terms of spherical surface harmonics Y n to derive the dispersion relation. The instability behavior is investigated by computing the growth rates and the most-unstable modes as a function of the spherical harmonic degree n. This general methodology is applicable to the entire parameter space spanned by the Atwood number, the viscosity ratio, and the dimensionless number B = (α R Ρ² 2 /μ² ² )¹ / ³ R (where α R , Ρ 2 and μ 2 are the local radial acceleration at the interface, and the density and viscosity of the denser overlying fluid, respectively). While the mathematical formulation here is general, this paper focuses on instability that arises at a spherical viscous fluid/vacuum interface as there is a great deal to be learned from the effects of one-fluid viscosity and sphericity alone. To quantify and understand the effect that curvature and radial accelerationhave on the Rayleigh-Taylor instability, a comparison of the growth rates, under homologous driving conditions, between the planar and spherical interfaces is performed. The derived dispersion relation for the planar interface accounts for an underlying finite fluid region of thickness L and normal acceleration α R . Under certain conditions, the development of the most-unstable modes at a spherical interface can take place via the superposition of two adjacent spherical harmonics Y n and Y n+1 . This bimodality in the evolution of disturbances in the linear regime does not have a counterpart in the planar configuration where the most-unstable modes are associated with a unique wave number.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Development of comprehensive two-dimensional liquid chromatography for investigating aging of plastic bonded explosives

The feasibility of measuring the aging and degradation of PBX 9501 via online two dimensional liquid chromatography (LC × LC) is investigated, and a preliminary instrumental setup and method is developed. Plastic-Bonded eXplosive (PBX) 9501 is nominally composed of 94.9 wt% HMX, 2.5 wt% Estane® 5703 (poly (ester urethane)), 2.5 wt% BDNPA/F (nitroplasticizer), 0.1 wt% Irganox 1010 and PBNA (N-phenyl-naphthylamine) at low concentrations. When exposed to various environmental conditions, PBX 9501 will degrade through different pathways. Because PBX 9501 is composed of both low molecular weight compounds (BDNPA/F, Irganox 1010, PBNA, and potential degradation products) and high molecular weight compounds (Estane® 5703), analysis is normally performed via two independent analyses. Additionally, the low molecular weight species are analyzed via high pressure liquid chromatography (HPLC) and the high molecular weight species via size exclusion chromatography (SEC). While these individual techniques yield information about the aging of PBX 9501, the combination of HPLC and SEC (i.e. HPLC × SEC) can simplify and streamline the analyses while also providing additional chemical information. A simplified sample preparation method is proposed for LC × LC analysis. Various SEC columns and HPLC column selection, flow rate, and gradient ramps were investigated for their application of measuring aged PBX 9501. Lastly, two LC × LC separations of a library standard of PBX 9501 and a sample of aged PBX 9501 are compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2011

Nanopatterning of Si surfaces by normal incident ion erosion: Influence of iron incorporation on surface morphology evolution

The surface morphology of Si(100) induced by 1200 eV Ar{sup +} ion bombardment at normal incidence with and without Fe incorporation is presented. The formation of nanodot patterns is observed only when the stationary Fe areal density in the surface is above a threshold value of 8x10{sup 14} cm{sup -2}. This result is interpreted in terms of an additional surface instability due to locally nonuniform sputtering in connection with the presence of a Fe rich amorphous phase at the peak of the nanodots. At Fe concentrations below the threshold, smoothing dominates and pattern formation is inhibited. The transition from a k{sup -2} to a k{sup -4} behavior in the asymptotic power spectral density function supports the conclusion that under these conditions ballistic smoothing and ion-enhanced viscous flow are the two dominant mechanisms of surface relaxation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Improvements in the reactor core isolation cooling (RCIC) pump model

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems that was observed during the Fukushima accidents and to extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine are being developed and tested with the RELAP-7 code to simulate the RCIC system. In the previous work, we developed a set of analytical models for the normal working conditions, as well as the two-phase off-design conditions, of the RCIC system, based on Sandia’s original work. An under-expanded jet model was developed to replace the CFD-based reduced-order model in the Sandia work to obtain the velocity and thermodynamic conditions of the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. For the two-phase off-design conditions, two well-established choking models were used - the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model. The new twophase Terry turbine model used the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. The two-phase nozzle expansion model was also developed by considering the vapor phase with a similar model for the single-phase case with the assumption that the liquid phase would slip along the wall with a much slower speed and will not contribute to the impulse on the rotor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Quantification of Off-Channel Inundated Habitat for Pacific Chinook Salmon ( Oncorhynchus tshawytscha ) along the Sacramento River, California, Using Remote Sensing Imagery

Off-channel areas are one of the most impacted aquatic habitats by humans globally, as extensive agricultural and urban development has limited them to roughly 10% of historical extent. This is also true for California’s Sacramento River Valley, where historically frequent widespread inundation has been reduced to a few off-channel water bodies along the mid-Sacramento River. This remaining shallow-water habitat provides crucial ecological benefits to multiple avian and fish species, but especially to floodplain-adapted species such as Chinook salmon (Oncorhynchus tshawytscha). Characterizing spatiotemporal off-channel dynamics, including inundation extent and residence time, is fundamental to better understanding the intrinsic value of such habitats and their potential to support recovery actions. Remote sensing techniques have been increasingly used to map surface water at regional and local scales, with improved resolutions. As such, this study maps off-channel inundation areas and describes their temporal dynamics by analyzing pixel-based time- series of multiple water indices, modified Normalized Difference Water Index (mNDWI) and the Automated Water Extraction Index (AWEI), generated from LandSat-8 and Sentinel-2 data between 2013–2021. Quantified off-channel area was similar with each water index and method used, but improved performance was associated with Sentinel-2 products and AWEI index to identify wetted areas under lower mainstem discharges. Results indicate an uneven distribution of off-channel habitat in the study area, with limited inundated areas in upstream reaches (<16% of total off-channel area for greater flows). In addition, much less habitat exists for flows under 400 m 3 /s, an important migration cue for endangered winter-run Chinook salmon, limiting juvenile access to areas with enhanced rearing conditions. Off-channel habitat residence times averaged between 7 and 16 days, primarily defined by the rate of receding flows, with rapid flow recession providing marginal off-channel habitat. This study shows reasonable performance of moderate resolution LandSat-8 and Sentinel-2 remote sensing imagery to characterize shallow-water inundated habitat in higher-order rivers, and as a method to inform restoration and native fish recovery efforts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Hyporheic Exchange in Sand Dunes Under a Freely Deforming River Water Surface

Bedform-driven hyporheic exchange is conditioned by the head gradients at the sediment-water interface. Local exchange phenomena between the surface and the subsurface are often driven by the dynamic forces related to the velocity distributions around sediment waves. In open channels, the static forces are represented by the water depth, yet most computational fluid dynamics models use a rigid-lid approximation. Here, we investigated whether and when the deformation of the river’s free-surface influences bedform-driven hyporheic exchange. This was done through simulations of coupled open channel and hyporheic flows with the air-water interface modeled either as a free-surface or a rigid-lid across increasing subcritical Froude numbers from ~0.05 to ~0.95. The normalized hyporheic flux was higher when considering free-surface deformation across most of the range of Froude numbers considered. When the Froude number was larger than 0.6 and smaller than 0.85, both hydrostatic and nonhydrostatic-driven fluxes increased significantly compared to the rigid-lid approach and the total flux was about five times that predicted with the rigid-lid. These results indicate that bedform-driven hyporheic fluxes are underestimated in most studies that use a rigid-lid assumption in subcritical flows.

58 GEOSCIENCES↗

DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗

DOE OSTI2021

Damage Modeling of Solid Oxide Fuel Cells Accounting for Redox Effects

This work applies a multiscale mechanistic damage model developed for brittle ceramics and implemented in commercial finite element (FE) packages via user subroutines to study progressive damage in solid oxide fuel cells (SOFC) subjected to thermomechanical loading under normal operating and shutdown conditions including redox effects. The damage model captures the micromechanics of stiffness reduction due to material porosity change and microcracking and integrates the as-obtained stiffness reduction law into a continuum damage mechanics (CDM) formulation for the evolution of microcracks up to fracture. The volumetric “swelling” that occurs during redox is treated in constitutive modeling similarly to thermal expansion, but swelling strains are irreversible. Furthermore, this damage model was first validated through predictions of strength and stress-strain response for the SOFC electrode materials. Next, it has been applied to predict the potential for degradation in a generic planar SOFC stack with large active area cells. Multicell stack models were simulated in both co-flow and counter-flow configurations. In addition, a constant temperature redox cycle was also simulated to capture overall cell electrode damage due to volumetric swelling of the nickel (Ni)-based anode in the anode-supported cells.

25 ENERGY STORAGE↗

DOE OSTI2017

Hydrological Drought in the Anthropocene: Impacts of Local Water Extraction and Reservoir Regulation in the U.S.

Hydrological drought is a substantial negative deviation from normal hydrologic conditions and is influenced by climate and human activities such as water management. By perturbing the streamflow regime, climate change and water management may significantly alter drought characteristics in the future. Here we utilize a high-resolution integrated modeling framework that represents water management in terms of both local surface water extraction and reservoir regulation, and use the Standardized Streamflow Index (SSI) to quantify hydrological drought. We explore the impacts of water management on hydrological drought over the contiguous US in a warming climate with and without emissions mitigation. Despite the uncertainty of climate change impacts, local surface water extraction consistently intensifies drought that dominates at the regional to national scale. However, reservoir regulation alleviates drought by enhancing summer flow downstream of reservoirs. The relative dominance of drought intensification or relief is largely determined by the water demand, with drought intensification dominating in regions with intense water demand such as the Great Plains and California, while drought relief dominates in regions with low water demand. Here, at the national level, water management increases the spatial extent of extreme drought despite some alleviations of moderate to severe drought. In an emissions mitigation scenario with increased irrigation demand for bioenergy production, water management intensifies drought more than the business-as-usual scenario at the national level, so the impacts of emissions mitigation must be evaluated by considering its benefit in reducing warming and evapotranspiration against its effects on increasing water demand and intensifying drought.

54 ENVIRONMENTAL SCIENCES↗