SEARCH · Engineering Papers

Results for “hydraulics”

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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Understanding Water Controls on Shale Gas/Oil Mobilization into Fractures (FY 2019-2021(FY2022, No Cost Extension), Final Report)

The large volumes of injected and produced water, and costs associated with both hydraulic fracturing fluids and produced water treatment motivate development of hydraulic fracturing strategies that use water more efficiently. The overall objective of this project is to improve understanding on how water is distributed within unconventional shales during and after hydraulic fracturing in order to provide a more rational basis for reservoir stimulation strategies that increase hydrocarbon recovery per unit of water used. The strong capillary imbibition of water into shale matrix pores and immiscibility with hydrocarbons makes it highly desirable to reduce water use during reservoir development. Given the importance, and competing influences of capillarity and gravity in controlling water in fractured porous media, this project also seeks to understand how these factors interact to obstruct or allow gas flow.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2024

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

DOE OSTI2024

FORGE STRESS annual report

The project's goal is to combine high-fidelity numerical models and true-triaxial block fracturing tests at high temperatures to understand the relationship between in situ stress, thermal effects, wellbore orientations and hydraulic fracture patterns. The numerical models are calibrated against field data, such as well pressures and microseismic data, and employed to estimate the in situ stress at the FORGE site. Laboratory experiments investigate the complex physics driving hydraulic fracture nucleation in EGS, enhancing understanding of the role of parameters like temperature, well orientation, and stress. Additionally, they are employed to validate some of the numerical tools used in the project. The project will have a significant impact by: (1) improving the characterization of the in-situ stress field at FORGE; (2) demonstrating the use of high-fidelity modeling tools for EGS; (3) providing a unique set of high-temperature hydraulic fracturing results to identify key components for in-situ stress estimation, validate current theories, and propose new ones; (4) offering a validated set of numerical tools within an open-source simulation framework, GEOS, that will be available to any future user.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Collaborations of the NEAMS Center of Excellence in 2025

The mission of the NEAMS Center of Excellence for Thermal-fluids Applications in Nuclear Energy is to advance the goals of the NEAMS program by providing leadership, best practices, research, and support and training for computational thermal hydraulics. In particular NEAMS tools are being developed to ultimately span the entire range of length- and time-scales required for a comprehensive reactor design and safety analysis. The focus of The Center is on the use, best practices, and deployment to stakeholders of the thermal hydraulic codes. The Center of Excellence, originally launched in 2018, has as its key goal to serve as a “front door” to industry. As part of this goal, The Center is in its sixth year of a program to engage in collaborative efforts between the laboratories and industry with the objective of stimulating cooperation and increasing the adoption of thermal hydraulics tools developed under NEAMS by the industry at large. With an overwhelmingly positive response from participants, previous short-term projects have led to more in-depth collaborations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Hydropower Biological Evaluation Toolset Best Practice Guide

Field studies using live fish are necessary for the evaluation of turbine biological performance, but they cannot determine the specific hydraulic conditions or physical stresses experienced by the fish, the locations where deleterious conditions occur, or the specific causes of the biological response. Using the Sensor Fish (SF) sensing technology, this deficiency can be overcome because the SF can be released independently or concurrently with live fish directly into operating infrastructure, and it takes high-frequency measurements of hydraulic conditions such as pressure, acceleration, and rotation acting on a body in situ during downstream passage. The Hydropower Biological Evaluation Tools (HBET; Hou et al. 2018) software package, developed by Pacific Northwest National Laboratory (PNNL), is designed to assemble, organize, and process data collected by the PNNL-developed SF and by live fish. HBET was developed specifically to design SF field studies, process the raw data, and analyze the processed data efficiently and scientifically. Its objectives are to facilitate SF studies focused on characterizing hydraulic conditions and to apply SF data for evaluating the impacts on fish from passage through hydro-structures. HBET allows users to design new studies, analyze data, perform statistical analyses, and evaluate predicted biological responses. It can be used by researchers, turbine designers, hydropower operators, and regulators to evaluate hydro-structures to enhance environmental sustainability in a cost-effective manner.

13 HYDRO ENERGY↗

DOE OSTI2021

Hydropower Energy Storage Capacity Dataset

The Hydropower Energy Storage Capacity (HESC) Dataset catalogues estimates of nominal energy storage capacity based on varying levels of detail. Dams and reservoirs selected were selected based on those reported in the National Inventory of Dams (NID 2019) and/or the Global Reservoir and Dam (GRanD v1.3) datasets. These data provide a foundation for understanding available resources at existing hydropower facilities and their potential to provide storage of energy and more flexible generation. Current estimates include Level 1 (based on maximum storage capacities and hydraulic head) and Level 2 (based on historical models or observations of reservoir volume and hydraulic head). For facilities where installed capacity is known, there are also estimates for discharge duration or the length of time when a facility could provide generation at a given capacity. Essential information used to calculate the energy storage capacity and discharge duration (volume, hydraulic head, and details about the sources or records used to obtain those parameters) and summaries of historical generation (for context) are also included.

13 HYDRO ENERGY↗

DOE OSTI2023

Hydropower Energy Storage Capacity (HESC) Dataset

The Hydropower Energy Storage Capacity (HESC) Dataset catalogs characteristics that are relevant to evaluating reservoir storage and estimates of energy storage capacity based on varying levels of detail. Hydropower dams and reservoirs were included based on information from the National Inventory of Dams (NID; USACE, 2021) and Global Reservoir and Dam (GRanD v1.3) and Existing Hydropower Assets datasets. These data provide a foundation for understanding available resources at existing hydropower facilities and their potential to provide storage of energy and more flexible generation. Estimates of energy storage capacity include: • Level 1 – nominal energy storage capacity based on maximum storage capacities and hydraulic head • Level 2 – nominal energy storage capacity based on historical models or observations of reservoir volume and hydraulic head. These estimates are provided based on capacity from the entire historical period as well as monthly values. • Level 3 – modeled energy generation based on volume-elevation relationships, historical storage, observed/modeled inflows, and hydraulic capacity of turbines and calculated both as overall and on a monthly basis. • Level 4 – modeled energy generation incorporating information from Level 3 and operational constraints. For facilities where installed capacity is known, there are also estimates for discharge duration (the length of time when a facility could provide generation at a given capacity).

13 HYDRO ENERGY↗

DOE OSTI2024

SPRUCE Xylem Native Embolism and Leaf Traits of Picea mariana and Larix laricina, 2019

This data set contains measurements of native embolism in branch xylem and associated branch and leaf traits from Picea mariana (Black Spruce) and Larix laricina (Tamarack) from September-October 2019 at the SPRUCE experiment (Hanson et al. 2017). Data are presented in one comma-separated (*.csv) file. Native embolism, a measurement of in-situ embolism in the xylem tissue that blocks water movement, measurements were conducted at the end of the growing season on cut branches using the hydraulic pipette method (see Peters et al. 2023 for full method). In short, branch segments were connected to hydraulic apprentice and flow rates of perfusion liquid were measured using graduated pipettes and stopwatch. After initial conductance measurements, branch segments were flushed using vacuum infiltration and hydraulic conductance was remeasured to calculate the percent loss in conductance due to embolism (PLC). Three branch segments (distal, middle, and proximal) were measured from each branch representing different diameters size classes. This dataset also contains leaf area associated with each measured branch segment, calculate leaf mass per area (LMA), sapwood specific conductivity (Ks), leaf area specific conductivity (Kl) and the sapwood area to leaf area ratio (Huber value). Measurements were made on mature trees from all ten treatment enclosures. One branch from each of five trees per species were used from each enclosure where available. Some plots do not contain five individual Larix laricina, in which case all available trees were sampled. All measurements were conducted in late September/early October 2019, at the end of growing season but before Larix laricina needle senescence. There are 10 experimental plots at SPRUCE: five temperature treatments (+0, +2.25, +4.5, +6.75, +9°C) at ambient CO2, and the same five temperature treatments at elevated CO2 (+500 ppm). These data were collected after the treatments had been running in full for three years, meaning much of the material measured was grown under treatment conditions.

Huber value↗

DOE OSTI2025

Comparison of spatial dynamics and point kinetics approaches in multiphysics modeling of the molten salt reactor experiment

In this work, we present validation test results of fully coupled neutronics and thermal-hydraulics models of the Molten Salt Reactor Experiment (MSRE) against experimental data of the zero power pump transients and the natural circulation tests at low power. To capture the strong coupling between neutronics and thermal-hydraulics due to fuel circulation, and to account for the delayed neutron precursor (DNP) distribution, the porous media thermal-hydraulics solver Pronghorn was fully coupled to the spatial neutron dynamics code Griffin, which solves the neutron diffusion equation, and to the 0-D point kinetics solver Squirrel, using a 2-D homogenized representation of the MSRE. The validation test results show very good agreement with experimental data for both point kinetics and spatial dynamics simulations, capturing the strong feedback effect and DNP losses in the MSRE. The 0-D code Squirrel accurately predicted the time-dependent behavior in the MSRE given the steady-state spatial dynamics solution of Griffin.

42 - ENGINEERING↗

DOE OSTI2025

Assessing the Heat Transfer Modeling Capabilities of CFD Software for Involute-Shaped Plate Research Reactors

The ongoing efforts to convert High-Performance Research Reactors (HPRRs) using Highly Enriched Uranium (HEU) to Low-Enriched Uranium (LEU) fuel require reliable thermal–hydraulic assessments of modified core designs. The involute-shaped fuel plates used in several major HPRRs present unique modeling challenges due to their compact core geometries and high heat flux conditions. This study evaluates the capability of three commercial CFD tools, STAR-CCM+, COMSOL, and ANSYS CFX, to predict cladding-to-coolant heat transfer using Reynolds-Averaged Navier–Stokes (RANS) methods within the thermal–hydraulic regimes of involute-shaped plate reactors. Broad sensitivity analysis was conducted across a range of reactor-relevant parameters using two turbulence models (k−ϵ and k−ω SST) and different near-wall treatment strategies. The results were benchmarked against the Sieder–Tate correlation and experimental data from historic studies. The codes produced consistent results, showing good agreement with the empirical correlation of Sieder–Tate and the experimental measurements. The findings support the use of these commercial CFD codes as effective tools for assessing the thermal–hydraulic performance of involute-shaped plate HPRRs and guide future LEU core development.

CFD↗

DOE OSTI2020

Post-Calibration Uncertainty Analysis for Travel Times at a Naval Weapons Industrial Reserve Plant

The Naval Weapons Industrial Reserve Plant (NWIRP) in McGregor, Texas began manufacturing explosives in 1980 and several hazardous chemicals were discovered in lakes and streams surrounding the site in 1998. Contaminants traveled to local lakes and streams much faster than initially predicted. This research estimated contaminant travel times and identified locations where monitoring wells should be installed to yield the greatest reductions in uncertainties in travel-time predictions. To this end, groundwater and particle-tracking models for NWIRP site were built to predict hydraulic heads and contaminant travel times. Next, parameter (hydraulic conductivities) uncertainties, parameter identifiabilities, observation (hydraulic heads) worth, and predictive (contaminant travel times) uncertainties were quantified. Parameter uncertainties were reduced by up to 92%; a total of 19 of 158 parameters were at least moderately identifiable; travel-time uncertainties were reduced up to 92%. Additionally, travel-time predictions and post-calibration parameter distributions were generated using the null-space Monte Carlo (NSMC) technique. NSMC predicted that conservative tracers exited the flow system within a year, which matches with field data. Finally, an observations-worth analysis found that additional 11 more measurements would reduce travel-time uncertainties by factors from 1.04 to 4.3 over existing data if monitoring wells were installed at the suggested locations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2022

Modeling Monthly Nitrate Concentration in a Karst Spring with and without Discrete Conduit Flow

Understanding the coupled continuum pipe-flow framework for modeling contaminant transport in karst systems is critical for protecting water resources therein. This study simulated point and non-point source nitrate concentration in a karst spring and investigated the results generated from the flow and transport model with and without discrete conduit flow. CFPv2 and CMT3D models were integrated to address the changes in nitrate concentration at a monthly scale in a karst spring, and the results were compared with that from an equivalent porous media (EPM) model with high-hydraulic conductivity (K) zones set in the main karstified area to represent conduits. The results show that the CFPv2+CMT3D model is able to describe well the recession of nitrate concentration in spring discharge, and the relatively larger deviation (slower nitrate recession) from the observed trend for the EPM model is probably a result of the limitation of utilizing high hydraulic conductivity cells to represent conduit. Moreover, simulated hydraulic heads in poorly karstified areas from the two models both show slight differences from the observations (the head RMSE values of calibration/validation for CFPv2 and MODFLOW models are 0.16 m/0.25 m and 0.26 m/0.17 m, respectively), indicating the inclusion of conduits may not affect the simulation considerably, and the lower the proportion of karstic area, the slight effects brought from the inclusion of conduits in the model. For highly karstified areas, the CFPv2+CMT3D model could provide more accurate results (head RMSE of calibration/validation for CFPv2 and MODFLOW are 0.22 m/0.06 and 0.52 m/0.47 m, respectively), showing the coupled continuum pipe-flow framework may be more appropriate for applying to highly and maturely karstified areas where the variations in the behavior of flow and contaminant transport are more affected by turbulent flow regime.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Coordination of rooting, xylem, and stomatal strategies explains the response of conifer forest stands to multi-year drought in the southern Sierra Nevada of California

Abstract. Extreme droughts are a major determinant of ecosystem disturbance that impacts plant communities and feeds back into climate change through changes in plant functioning. However, the complex relationships between aboveground and belowground plant hydraulic traits and their role in governing plant responses to drought are not fully understood. In this study, we use a model, the Functionally Assembled Terrestrial Ecosystem Simulator in a configuration that includes plant hydraulics (FATES-Hydro), to investigate ecosystem responses to the 2012–2015 California drought in comparison with observations at a site in the southern Sierra Nevada that experienced widespread tree mortality during this drought. We conduct a sensitivity analysis to explore how different plant water sourcing and hydraulic strategies lead to differential responses during normal and drought conditions. The analysis shows the following. Deep roots that sustain productivity through the dry season are needed for the model to capture observed seasonal cycles of evapotranspiration (ET) and gross primary productivity (GPP) in normal years, and deep-rooted strategies are nonetheless subject to large reductions in ET and GPP when the deep soil reservoir is depleted during extreme droughts, in agreement with observations. Risky stomatal strategies lead to greater productivity during normal years as compared to safer stomatal control, but they also lead to a high risk of xylem embolism during the 2012–2015 drought. For a given stand density, stomatal and xylem traits have a stronger impact on plant water status than on ecosystem-level fluxes. Our study highlights the significance of resolving plant water sourcing strategies to represent drought impacts on plants and consequent feedbacks in models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modeling the mechanisms of coastal vegetation dynamics and ecosystem responses to changing water levels

Coastal forests are increasingly experiencing mortality due to inundation by fresh- and seawater, leading to their replacement by marshes. These shifts alter vegetation composition, biogeochemical cycling, carbon storage, and hydrology. Using a hydraulically enabled ecosystem demography model (FATES-Hydro), we conducted numerical experiments to investigate the mechanisms behind inundation-driven forest loss and the ecosystem-scale consequences of forest-to-marsh transitions. We compared mortality processes and their effects across broadleaf and conifer trees at two coastal sites – Lake Erie (freshwater) and Chesapeake Bay (saline). Our simulations show that hydraulic failure, driven by root loss under prolonged flooding, is the primary mortality mechanism across both tree types and sites. Forest replacement by marsh reduced ecosystem-scale leaf area index (LAI), gross primary production (GPP), transpiration, and deep soil water uptake in conifer forests, while broadleaf forests experienced smaller changes due to lower initial LAI and greater marsh compensation. Marsh invasion occurred following canopy thinning driven by tree mortality. These findings suggest that, under similar root loss, hydraulic failure dominates coastal tree mortality regardless of species or water type, with denser forests experiencing stronger ecosystem impacts. Our study identifies key mortality mechanisms and offers testable hypotheses for future empirical studies on coastal vegetation change.

Ding, Junyan [Occidental College, Los Angeles, CA ↗

DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗

DOE OSTI2025

Temporal persistence of postfire flood hazards under present and future climate conditions in southern Arizona, USA

Changes to soil hydraulic properties that reduce infiltration capacity following fire can increase flash flood risks. These risks are exacerbated by rainfall intensification associated with a warming climate. However, the potential effects of climate-change-driven rainfall intensification on postfire floods remain largely unexplored. Using rainfall and runoff observations from a 49.4 km 2 watershed in southern Arizona, USA, and a hydrologic model (KINEROS2), we examined the temporal evolution following a historic fire of three crucial hydrologic parameters: soil saturated hydraulic conductivity (K sp ), net capillary drive (G p ), and hydraulic roughness (n c ). We explored how the effect of fire on these parameters may influence peak flow under future climate scenarios derived from CMIP6, specifically the medium emissions scenario (SSP2-4.5) and high emissions scenario (SSP5-8.5). Results demonstrate an increase in K sp from 11 mm h −1 in the first postfire year to 60 mm h −1 in postfire year 3. G p similarly increased from 19 mm in the first postfire year to 30 mm in the third, while nc was relatively constant. The highest simulated Q p occurred in the first postfire year. Under the SSP2-4.5 scenario, the likelihood of a 100-year flood is projected to be twice as large by the middle of the century relative to its historical magnitude. Simulations further indicate that the maximum expected discharge associated with a postfire flood, as derived from historical data, could be triggered by a 10-year rainstorm under the SSP5-8.5 scenario by the late century. Simulations also demonstrate that rainfall intensification will lead to greater persistence of elevated flood hazards following fire by the late century under both the SSP2-4.5 and the SSP5-8.5 scenarios.

54 ENVIRONMENTAL SCIENCES↗