Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Refining”

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

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load

Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward large-scale simulation of residual stress and distortion in wire and arc additive manufacturing

In this study, we aim to advance the structural analysis of wire and arc additive manufacturing (WAAM) by considering the thermomechanical features inherent in direct energy deposition. Simulation approaches including the iterative substructure method (ISM), dynamic mesh refining method (DMRM), and graphics processing unit (GPU) based explicit finite element method (FEM) were developed for accelerating additive manufacturing stress analysis that is very time consuming by conventional numerical methods. ISM and DMRM take advantage of the strong nonlinearity phenomenon near the moving heat source by reducing the global iterations and refining the local mesh, respectively. In addition to the use of GPUs, the explicit FEM is accelerated by a time scaling technique based on the inherent strain concept. The residual stress and distortion of two large builds were analyzed, showing very consistent numerical results and good agreement with experiments. Compared with the commercial software Abaqus, the novel approaches reduced the computational cost substantially without compromising accuracy. The GPU code showed the highest computational efficiency (30~70×), while DMRM and ISM had acceleration factors of 9× and 3×, respectively. Finally, such high-fidelity modeling approaches will be very useful for building up a digital twin of WAAM to reduce development time and cost.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
DOE OSTI2023

An implicit-in-time DPG formulation of the 1D1V Vlasov-Poisson equations

Efficient solution of the Vlasov equation, which can be up to six-dimensional, is key to the simulation of many difficult problems in plasma physics. The discontinuous Petrov-Galerkin (DPG) finite element methodology provides a framework for the development of stable (in the sense of Ladyzhenskaya–Babuška–Brezzi conditions) finite element formulations, with built-in mechanisms for adaptivity. While DPG has been studied extensively in the context of steady-state problems and to a lesser extent with space-time discretizations of transient problems, relatively little attention has been paid to time-marching approaches. In the present work, we study a first application of time-marching DPG to the Vlasov equation, using backward Euler for a Vlasov-Poisson discretization. We demonstrate adaptive mesh refinement for two problems: the two-stream instability problem, and a cold diode problem. Furthermore, we believe the present work is novel both in its application of unstructured adaptive mesh refinement (as opposed to block-structured adaptivity, which has been studied previously) in the context of Vlasov-Poisson, as well as in its application of DPG to the Vlasov-Poisson system. We also discuss extensive additions to the Camellia library in support of both the present formulation as well as extensions to higher dimensions, Maxwell equations, and space-time formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Interlaboratory Study: Testing Reproducibility of Solid Biofuels Component Identification Using Reflected Light Microscopy

Considering global market trends and concerns about climate change and sustainability, increased biomass use for energy is expected to continue. As more diverse materials are being utilized to manufacture solid biomass fuels, it is critical to implement quality assessment methods to analyze these fuels thoroughly. One such method is reflected light microscopy (RLM), which has the potential to complement and enhance current standard testing, leading to improving fuel quality assessment and, ultimately, preventing avoidable air pollution. An interlaboratory study (ILS) was conducted to test the reproducibility of biomass fuels component identification using a reflected light microscopy technique. The exercise was conducted on thirty photomicrographs showing biomass and various undesired components (like plastics or mineral matter), which were purposely added (by the ILS organizers) to contaminate wood pellets and charcoal-based grilling fuels. Forty-six participants had various levels of difficulty identifying the marked components, and as a result, the percentage of correct answers ranged from 52.2 to 94.4%. Among the most difficult components to distinguish were petroleum products and inorganic matter. Various reasons led to the misidentification, including insufficient morphological descriptions of the components provided to participants, ambiguities of the nomenclature, limitations of the analytical and exercise method, and insufficient experience of the participants. Overall, the results indicate that RLM has the potential to enhance the quality assessment of biomass fuels. However, they also demonstrate that the petrographic classification used in this exercise requires further refinement before it can be standardized. While a new simplified classification of solid biomass fuels components was created as an outcome of this study, future research is necessary to refine the nomenclature, develop a microscopic morphological description of the components, and verify the accuracy of component identification with a follow-up ILS.

biomass↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Multifidelity methods for uncertainty quantification of a nonlocal model for phase changes in materials

This study is devoted to the construction of a multifidelity Monte Carlo (MFMC) method for the uncertainty quantification of a nonlocal, non-mass-conserving Cahn-Hilliard model for phase transitions with an obstacle potential. Here, we are interested in estimating the expected value of an output of interest (OoI) that depends on the solution of the nonlocal Cahn-Hilliard model. As opposed to its local counterpart, the nonlocal model captures sharp interfaces without the need for significant mesh refinement. However, the computational cost of the nonlocal Cahn-Hilliard model is higher than that of its local counterpart with similar mesh refinement, inhibiting its use for outer-loop applications such as uncertainty quantification. The MFMC method augments the desired high-fidelity, high-cost OoI with a set of lower-fidelity, lower-cost OoIs to alleviate the computational burden associated with nonlocality. Most of the computational budget is allocated to sampling the cheap surrogate models to achieve speedup, whereas the high-fidelity model is sparsely sampled to maintain accuracy. For the non-mass-conserving nonlocal Cahn-Hilliard model, the use of the MFMC method results in, for a given computational budget, about an order of magnitude reduction in the mean-squared error of the expected value of the OoI relative to that of the Monte Carlo method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2020

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗
DOE OSTI2021

The intrinsic value of bioblendstocks: Prenol as a case study

Decarbonizing the transportation sector is likely to require both electrification and increased incorporation of biofuels and/or bioblendstocks. While the social and environmental benefits of bioblendstocks are well understood, their real value for the fuel producers has not been established. As such, this work considers prenol as a bioblendstock case study to identify sources of intrinsic value to fuel blenders by studying the properties of binary mixtures with gasoline components. The considered refinery blendstocks were samples of full range naphthas from the distillation, fluidized catalytic cracking, isomerization, alkylation, and reforming units. Octane numbers, Reid vapor pressure, distillation curves, and sulfur content were evaluated. Our results indicate the need for adjusting the formulation of the base fuel, depending on the interplay among the properties of the bioblendstock and those of the base fuel. Prenol increased research octane number (RON) and octane sensitivity (OS) of the base fuel, by up to 25 and 10 octane numbers, respectively. Additionally, 10 vol% prenol reduced RVP up to 2.2 psi, for the more volatile blendstock. Thus, considering prenol as a low volatility, RON/OS boosting bioblendstock, the composition of the preferred base fuel was proposed as containing reduced olefins and aromatics, and increase light fractions. The potential impact of this new gasoline formulation on refining processes and products gives rise to direct sources of value to the refiners, such as exporting products to the chemicals market, increasing the value of intermediate refinery streams, decreasing operating severity of certain refinery units, and broadening of the product suite.

09 BIOMASS FUELS↗
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 OSTI2026

Influence of rolling reduction and annealing on recrystallization and grain structure in Ta-2.5W alloys

The microstructural evolution of wrought Tantalum - 2.5 wt% Tungsten (Ta-2.5W) alloys during thermomechanical processing is critical for optimizing their mechanical reliability in demanding applications such as aerospace, chemical processing, and nuclear technology. Despite the widespread use of Ta-W alloys, a comprehensive understanding of how rolling reduction, annealing temperature, and elemental inhomogeneity interact to determine recrystallization behavior and grain refinement remains incomplete. Here, in this study, we systematically investigate the effects of cold rolling and subsequent annealing on the microstructure of Ta-2.5W, with particular attention to grain orientation, stored energy, and elemental banding. Our results demonstrate that higher rolling reduction rates lower the onset and completion temperatures for recrystallization, resulting in finer and more homogeneous grain structures. Electron backscatter diffraction (EBSD) analysis reveals that grains with a 〈111〉 parallel to the plate normal possess higher stored energy and nucleate recrystallization more readily than grains with a 〈001〉 parallel to the plate normal. Elemental mapping shows that tungsten inhomogeneity leads to localized bands of accelerated recrystallization and hardness variation. These findings provide new insights into the mechanisms of microstructure refinement in Ta-2.5W alloys, offering guidance for tailoring processing routes to achieve superior performance in demanding engineering environments.

Annealing↗
DOE OSTI2021

Thermal Processing Conditions for the Synthesis of Near Theoretical Density Li 5 La 3 Ta 2 O 12 Ceramics for Ceramic Dual-Mode Detectors

A family of Li-containing garnet compounds, extensively investigated for Li-ion solid electrolyte applications, show great potential to be developed into transparent ceramic dual-mode scintillators. In particular, Li 5 La 3 Ta 2 O 12 (LLTaO) shows promise as an efficient dual-mode scintillator due to its high Li content, high density, non-hygroscopicity, and high effective-Z. Often the synthesis of LLTaO ceramics involves solid state processing techniques starting with oxide and carbonate powders, calcining, sintering for the removal of residual porosity, and post-processing. This study focuses on two objectives: (1) establishing thermal processing boundary conditions for a flexible alkoxide-based sol-gel method for synthesizing fine-grained powders; and (2), consolidating powders via hot-pressing to form high-density ceramics. In this study, near single-phase ceramics were achieved by hot-pressing sol-gel powders with varying amounts of excess Li. Rietveld refinements were used to analyze X-ray diffraction data and provide a quantitative evaluation of the phases present in the hot pressed ceramics. Some of the diffraction data showed large low angle peak asymmetry, suggesting the presence of two cubic Li 5 La 3 Ta 2 O 12 garnet phases with different lattice parameters. Depending on the synthesis and processing conditions, the refinement revealed ~ 90 wt% concentration of a phase with a ≈ 12.8 Å and ~ 10 wt% concentration of a phase with a ≈ 12.9 Å, indicating a remaining excess of Li in the hot pressed ceramics in the large lattice parameter phase. These synthesis experiments allowed for the determination of a baseline for Li loss during hot-pressing. X-ray excited radioluminescence and photoluminescence excitation and emission spectra were obtained from the hot pressed ceramics. In conclusion, the measured emission spectra featured broad tantalate emission and line emissions likely from a rare-earth impurity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase stability in cast and additively manufactured Al-rich Al-Cu-Ce alloys

Additively manufactured (AM) eutectic Al alloy systems have been studied extensively for advantageous thermal stability and mechanical properties due to their refined microstructures. Al-Cu-Ce alloys are one subset of these AM eutectic alloys. Here we studied phase stability in AM Al-Cu-Ce alloys and compared it to that of conventionally cast ones. A new phase, Al 8 Cu 3 Ce, was identified in the microstructures of both AM and cast Al-Cu-Ce alloys. This Al 8 Cu 3 Ce phase was not previously included on experimental or thermodynamically calculated phase diagrams of the Al-Cu-Ce system. Therefore, we performed additional thermodynamic modeling of the system. These models were experimentally validated with cast and subsequently heat-treated Al-Cu-Ce alloys. We found that despite the refined microstructure of the AM alloys, the phases formed were consistent with the cast alloys, suggesting that AM processing did not significantly alter the formation and stability of phases from that in the conventional alloys. In conclusion, this work has resolved inconsistent previous descriptions of the Al-Cu-Ce ternary phase diagram in the Al-rich region and resulted in the addition of the Al 8 Cu 3 Ce as an equilibrium phase above 500 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗