Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SATURATION”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Enhancement and Saturation of Near-Field Radiative Heat Transfer in Nanogaps between Metallic Surfaces

Near-field radiative heat transfer (NFRHT) between planar metallic surfaces was computationally explored over five decades ago by Polder and van Hove [Phys. Rev. B 4, 3303 (1971)]. These studies predicted that, as the gap size (d) between the surfaces decreased, the radiative heat flux first increases by several orders of magnitude until d is ~100 nm after which the heat flux saturates. However, despite both the fundamental and practical importance of these predictions, the combined enhancement and saturation of NFRHT at small gaps in metallic surfaces remains experimentally unverified. Here, we probe NFRHT between planar metallic (Pt, Au) surfaces and show that RHT rates can exceed the far-field rate by over a thousand times when d is reduced to ~25 nm. More importantly, we show that for small values of d RHT saturates due to the dominant contributions from transverse electric evanescent modes. Finally, our results are in excellent agreement with the predictions of fluctuational electrodynamics and are expected to inform the development of technologies such as near-field thermophotovoltaics, radiative heat-assisted magnetic recording, and nanolithography.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Common-mode inductor saturation analysis and design optimization based on spectrum concept

Understanding the CM inductor core saturation mechanism and reducing core flux density is critical for CM inductor design optimization. Instead of a time domain method, this paper introduces frequency domain spectrum concept for CM inductor core saturation analysis and design optimization, which will provide designers a better understanding of CM inductor design. First, both core permeability and converter modulation index's opposite influence on DM flux density and CM flux density are identified. Then, CM flux density is further investigated based on the spectrum concept. Three components in the CM inductor which may cause large CM flux density and core saturation are summarized: (1) switching frequency related components, (2) impedance resonance frequency related components, and (3) modulation frequency related components. Each component is investigated for CM flux density reduction and filter design optimization. A connecting AC and DC side midpoint with notch filter structure is proposed to reduce modulation frequency related components. Experiment results are presented to verify the proposed concept and method.

Chen, Ruirui↗
DOE OSTI2019

Dissolution of silica component of glass network at early stage of corrosion in initially silica-saturated solution

Corrosion of glass in silica-saturated solution has been performed with the assumption that dissolution of silicate species from the glass network would not occur. Using surface-sensitive analytical techniques, we report experimental evidence suggesting the dissolution of silicate network species from a model nuclear waste glass, called international simple glass (ISG), in an aqueous solution initially saturated with soluble silica species. Results from low energy ion scattering and X-ray photoelectron spectroscopy reveal a complete depletion of mobile element species (B, Na) from the ISG surface and an enrichment of Zr on the outmost surface. In support of spectroscopic analyses, results from topographic imaging with atomic force microscopy show a stochastic dissolution of glass surface resulting in a higher surface roughness with increasing corrosion time in aqueous solution. This study shows that a true equilibrium between soluble silica species in the solution phase and silicate species in the glass network could not be warranted by performing corrosion experiments in the solution where dissolved silica species are initially equilibrated with amorphous silica in the presence of KOH. The leaching of mobile species (B, Na) could affect the saturation level of aqueous solution and induce further dissolution of the glass surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Drying of a Fully Saturated Porous Medium With Excess Water Layers: A Numerical Study

Abstract Drying of moist porous media can be very energy inefficient. For example, in the pulp and paper industry, paper drying consumes more than two-thirds of the total energy used in paper machines. Novel drying technologies can decrease the energy used for drying and lessen the manufacturing processes' carbon footprint. Developing next-generation drying technologies to dry moist porous media may require an understanding of removing moisture from a fully saturated porous material with excess water. This paper provides a fundamental understanding of heat and mass transfer in a fully saturated porous medium with excess water. This is relevant, for example, in drying tissue as well as pulp or paper for the purpose of thermal insulation where pressing is preferred to be avoided to overcome the reduction in the sheet thickness. For this purpose, a theoretical drying model is developed where the porous medium corresponds to paper and is assumed to be sandwiched between two excess-water layers (bottom and top). The conjugate model consists of energy and mass conservation equations for each layer. The model is validated with corresponding experimental data. In the model, the thickness of each water layer is calculated as a function of drying time based on local temperature and total moisture content. The numerical model is transient and one-dimensional in space (i.e., in the thickness direction). This paper demonstrates the governing equations, boundary conditions, and results when the saturated porous medium with water layers is heated from one side. Moisture and temperature profiles are estimated in the thickness direction of the porous medium as it dries.

Engineering↗
DOE OSTI2015

Orthorhombic YBaCo{sub 4}O{sub 8.4} crystals as a result of saturation of hexagonal YBaCo{sub 4}O{sub 7} crystals with oxygen

Hexagonal YBaCo{sub 4}O{sub 7} crystals (sp. gr. P6{sub 3}mc, a{sub h} = 6.3058(4) Å, c{sub h} = 10.2442(7) Å, Z = 2) are saturated with oxygen to the YBaCo{sub 4}O{sub 8.4} composition and studied by X-ray diffraction (XRD) analysis. The saturation is completed by a structural first-order phase transition to orthorhombic crystals (sp. gr. Pbc2{sub 1}, a{sub o} = 31.8419(2) Å, b{sub o} = 10.9239(5) Å, c{sub o} = 10.0960(5) Å, Z = 20). The connection of two lattices is expressed in terms of the action of matrix (500/120/001) on the hexagonal basis. Five structural fragments of the same type but with different degrees of order alternate along the long axis of the oxygen-saturated orthorhombic structure. The XRD data on single crystals differ from the results obtained by other researchers on ceramic samples; possible causes of these differences are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Increasing Saturated Electron-Drift Velocity in Donor–Acceptor Doped pHEMT Heterostructures

Field dependences of the electron-drift velocity in typical pseudomorphic high-electron-mobility transistor (pHEMT) heteroepitaxial structures (HESs) and in those with donor–acceptor doped (DApHEMT) heterostructures with quantum-well (QW) depth increased by 0.8–0.9 eV with the aid of acceptor layers have been studied by a pulsed technique. It is established that the saturated electron-drift velocity in DA-pHEMT-HESs is 1.2–1.3 times greater than that in the usual pHEMT-HESs. The electroluminescence (EL) spectra of DA-pHEMT-HESs do not contain emission bands related to the recombination in widebandgap layers (QW barriers). The EL intensity in these HESs is not saturated with increasing electric field. This is indicative of a suppressed real-space transfer of hot electrons from QW to barrier layers, which accounts for the observed increase in the saturated electron-drift velocity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Probing gluon saturation via diffractive jets in ultra-peripheral nucleus-nucleus collisions

Abstract We argue that semi-inclusive photo-production of a pair of hard jets via coherent diffraction in nucleus-nucleus ultra-peripheral collisions at high energy is a golden channel to study gluon saturation. The dominant contribution is the diffractive production of three jets in an asymmetric configuration. Two of the jets are hard and propagate at nearly central pseudo-rapidities. The third jet is semi-hard, with transverse momentum comparable to the nuclear saturation momentum, and is well separated in pseudo-rapidity from the hard dijets. The emission of the semi-hard jet allows for strong scattering, thus avoiding the “higher-twist” suppression of the exclusive dijet production due to colour transparency. We compute the trijet cross-section using the diffractive TMD factorisation which emerges from the CGC effective theory at high energy. The cross-section is controlled by gluon saturation, which leaves its imprints on the structure of the final state, notably on the rapidity distribution.

Physics↗
DOE OSTI2017

Streamline-Based Time-Lapse-Seismic-Data Integration Incorporating Pressure and Saturation Effects

In this paper, we present an efficient history-matching technique that simultaneously integrates 4D repeat seismic surveys with well-production data. This approach is particularly well-suited for the calibration of the reservoir properties of high-resolution geologic models because the seismic data are a really dense but sparse in time, whereas the production data are finely sampled in time but spatially averaged. The joint history matching is performed by use of streamline-based sensitivities derived from either finite-difference or streamline-based flow simulation. For the most part, earlier approaches have focused on the role of saturation changes, but the effects of pressure have largely been ignored. Here, we present a streamline-based semianalytic approach for computing model-parameter sensitivities, accounting for both pressure and saturation effects. The novelty of the method lies in the semianalytic sensitivity computations, making it computationally efficient for high-resolution geologic models. The approach is implemented by use of a finite-difference simulator incorporating the detailed physics. Its efficacy is demonstrated by use of both synthetic and field applications. Lastly, for both the synthetic and the field cases, the advantages of incorporating the time-lapse variations are clear, seen through the improved estimation of the permeability distribution, the pressure profile, the evolution of the fluid saturation, and the swept volumes.

58 GEOSCIENCES↗
DOE OSTI2016

Relating Quantitative Soil Structure Metrics to Saturated Hydraulic Conductivity

Core Ideas Soil structure quantified in the field using a 3D laser scanning technique (MLT). MLT data and coefficient of linear extensibility combined into a structure metric. Soil water contents recorded at four depths in a field lysimeter and used in HYDRUS‐1D. Saturated conductivity ( K s ) estimated with a Markov chain Monte Carlo approach. New metric better correlated to K s in horizons with strong soil structure. Soil structure affects saturated hydraulic conductivity ( K s ) by creating highly conductive macropores that preferentially transmit soil water. In this study, we explored the relationship between K s and macropores in an Oxyaquic Vertic Argiudoll in northeastern Kansas. Macropores were quantified from an excavation wall using multistripe laser triangulation (MLT) scanning. Soil water contents were measured at four depths within a soil lysimeter installed within 2 m of the MLT‐scanned soil profile and adjacent to an Ameriflux tower monitoring precipitation, air temperature, and solar radiation. Selected hydraulic properties of soil horizons within the lysimeter were optimized to water content data using a Markov chain Monte Carlo technique in combination with the mobile–immobile water (MIM) model in HYDRUS‐1D. Estimates of K s varied between 4198 cm d −1 in the A horizon and 0.6 cm d −1 in a 2Btss2 horizon with strongly expressed wedge structure. Approximately 87% of the variation in K s was explained by the geometric mean of the widths of pores quantified with the MLT technique and modified by the coefficient of linear extensibility (COLE). The use of the COLE allows the widths of the macropores obtained under dry conditions to be approximated at saturation. Two models that predict K s from either texture or water retention data resulted in K s estimates that were similar to each other but significantly lower than K s values predicted with MIM in horizons where structural pores dominate water flow. This technique shows a great deal of promise in better understanding and predicting the relationship of soil structure to water flow.

Eck, Dennis V.↗
DOE OSTI2018

Characterization of Grain‐Size Distribution, Thermal Conductivity, and Gas Diffusivity in Variably Saturated Binary Sand Mixtures

Core Ideas Parameterized models are relevant to applications where different porous media mixtures are used. Models include a grain‐size distribution function to describe bimodal behavior for binary mixtures. Improved model describes observed thermal conductivity–saturation relations for binary mixtures. Combined Buckingham–Penman model used to describe observed gas diffusivity–air content relations. This work is relevant for proper simulation of mixed porous media. Characterization of differently textured porous materials, as well as different volumetric porous media mixtures, in relation to mass and heat transport is vital for many engineering and research applications. Functional relations describing physical (e.g., grain‐size distribution, total porosity), thermal, and gas diffusion properties of porous media and mixtures are necessary to optimize the design of porous systems that involve heat and gas transport processes. However, only a limited number of studies provide characterization of soil physical, thermal, and gas diffusion properties and the functional relationships of these properties under varying soil water contents, especially for soil mixtures, complicating optimization efforts. To better understand how mixing controls the physical, thermal, and gas diffusion properties of porous media, a set of laboratory experiments was performed using five volumetric mixtures of coarse‐ and fine‐grained sand particles. For each mixture, the grain‐size distribution (GSD), thermal conductivity, and gas diffusivity were obtained and parameterized using existing and suggested parametric models. Results show that the extended, two‐region Rosin–Rammler particle‐size distribution model proposed in this study could successfully describe the bimodal behavior of the GSD of binary mixtures. Further, the modified Côté and Konrad thermal conductivity model adequately described the thermal conductivity–water saturation relations observed in different mixtures. The proposed simple soil‐gas diffusivity descriptive model parameterized the upper limit, average, and lower limit behavior in gas diffusivity–air content relations in apparently texture‐invariant gas diffusivity data. Results further show a close analogy between gas diffusivity and thermal conductivity and their variation with saturation across different binary mixtures. Overall, the results of the study provide useful numerical insight into the physical, thermal, and gas transport characteristics of binary mixtures, with wide implications for future engineering and research applications that involve multicomponent porous systems.

Deepagoda, T. K. K. Chamindu↗
DOE OSTI2021

Characterization of the Sulfur-Saturated Melt Versions of the LAW Phase 4 Glasses

In this report, results are provided from the analyses of a series of sulfur-saturated melt versions of simulated nuclear waste glasses, and chemical analysis of the wash solution resulting from the preparation of these glasses. The glasses were selected and fabricated by the Pacific Northwest National Laboratory as part of a broader study of the influence of glass composition on chemical durability, sulfur retention, and other properties. The resulting data will be used in the development of enhanced property/composition models for waste vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the sulfur-saturated versions of the quenched glasses to allow for comparisons with the targeted compositions, as well as the measured compositions of the quenched glasses. The measured concentrations of chlorine, and fluorine were below the targeted values for most of the study glasses, likely because of volatility during the multiple melting steps. The measured concentrations of B 2 O 3 , K 2 0, and Li 2 0 were generally low relative to the targeted values. As expected, the measured concentrations of SO 3 in most of the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses. The wash solutions contained mainly sodium and sulfur, with moderate concentrations of potassium.

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

Vadose and Saturated Zone Flow and Transport Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

The purpose of the fate and transport modeling described in this environmental calculation file (ECF) is to evaluate the impacts to groundwater associated with waste disposal operations at Mixed Waste Low-Level Burial Ground (LLBG) Trenches 31 and 34 to satisfy requirements in DOE O 435.1, Radioactive Waste Management. The model integrates the flow and transport in the vadose zone beneath the active trenches with the saturated zone downgradient of the trenches to predict the radionuclide concentration at the point of assessment (POA). DOE M 435.1-1, Radioactive Waste Management Manual, defines the POA as the point of highest projected dose or concentration beyond a 100 m (328 ft) buffer zone surrounding the disposed waste. The modeling is conducted in accordance with the DOE G 435.1, Implementation Guide for Use with DOE M 435.1-1, performance assessment (PA) guidelines. The modeling involves evaluation of the groundwater concentrations and radionuclide arrival times during the 1,000-year compliance and 10,000-year sensitivity-uncertainty periods per DOE O 435.1 and DOE M 435.1. This analysis does not consider radionuclide release during facility operations, only the post-closure impacts of the radionuclides to the environment. The evaluation of potential radiological dose to groundwater receptors caused by releases from a closed facility containing radioactive waste typically includes the following: (1) Release of radionuclides from that facility (2) Transport of those radionuclides through the environment, and (3) Exposure to humans to environmental concentration levels of those radionuclides The fate and transport three-dimensional (3D) model analysis involves the post-closure impacts to the environment of the technetium-99, iodone-129, and uranium (all isotopes in the waste). The residual inventory estimates include several radionuclides, but technetium-99 is typically responsible for almost all of the beta-gamma dose equivalent associated with groundwater (water resources) protection per 40 CFR 141, “National Primary Drinking Water Regulations” (e.g., see the results in WCH-520, Performance Assessment of Environmental Restoration Disposal Facility, Hanford Site, Washington; hereinafter referred to as the ERDF PA), and iodone-129 can also be a significant dose contributor for some waste (e.g., RPP-RPT-59958, Performance Assessment for the Integrated Disposal Facility, Hanford Site, Washington; hereinafter referred to as the IDF PA). Uranium does not typically factor significantly into the impacts to groundwater, even during the 10,000-year sensitivity-uncertainty period, but always remains of interest as a contaminant. This ECF does not address vadose and saturated zone modeling for Trench 94 of the 200 East Area LLBG. Current information confirms the validity of the low corrosion rate of the naval reactor plant carbon steel (HY-80), and the even lower corrosion rate of the nickel-iron-chromium alloy reactor vessel (Inconel Alloy 600) presented in DOE/EIS-0259, Final Environmental Impact Statement on the Disposal of Decommissioned, Defueled Cruiser, Ohio Class, and Los Angeles Class Naval Reactor Plants. Based on these low corrosion rates, the time to breach the reactor vessel to allow release of radionuclides from the activated metal of the reactor vessel internal structure is at least 10,000 years. This time to breach precludes the need to evaluate the vadose and saturated zone transport of contaminants released from the reactor compartment disposal packages in the 200 East Area LLBGs PA (CP-63826, Waste Release Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington).

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

Characterization of the Sulfur-Saturated Melt Versions of the HLW HAlG Glasses

This report provides the results from the chemical analyses of a series of sulfur-saturated melt versions of the High-Level Waste High-Aluminum Glass study glasses, a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory. These data will be used in the development, validation, and implementation of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the sulfur-saturated melt versions of the glasses to allow for comparisons with targeted compositions as well as the measured compositions of the quenched glasses. The relative differences between the targeted and measured concentrations of B 2 O 3 , Cr 2 O 3 , K 2 O, Li 2 O, Na 2 O, and P 2 O 5 for several of the glasses were greater than ±10%. The relative differences between the targeted and measured concentrations of Al 2 O 3 and ZrO 2 for one of the glasses were greater than ±10%. As expected, the measured concentrations of SO 3 in most of the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses. The wash solutions contained mainly sodium, sulfur, and sulfate ions.

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

Characterization of the Sulfur-Saturated Melt Versions of the LAW ML1 Glasses

This report provides the results from the chemical analyses of the sulfur-saturated melt versions of the Low-Activity Waste Machine Learning study glasses, a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory. These data will be used in the development, validation, and implementation of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the sulfur-saturated melt versions of the glasses to allow for comparisons with targeted compositions as well as the measured compositions of the quenched glasses. The relative differences between the targeted and measured concentrations of F- for one glass, K 2 O for one glass, Na 2 O, for one glass, P 2 O 5 for one glass, and ZrO 2 for several of the glasses were greater than ±10%. As expected, the measured concentrations of SO 3 in most of the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses.

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

Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns: Modeling Archive

This Modeling Archive is in support of a TES-SFA publication “Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns” (Craig et al., 2021). We ran and evaluated a multi-assumption soil organic carbon (SOC) model to investigate whether alternative assumptions regarding constraints on soil microbial biomass could lead to soil carbon saturation patterns. We developed this model in the Multi-Assumption Architecture and Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_Craig2021; Walker et al. 2018). Using MAAT, we embedded three alternative hypotheses in a microbially explicit three-pool SOC model: 1) the efficiency of mineral-associated SOC formation decreases as mineral-associated SOC approaches a maximum value (“Mineral saturation”), 2) the microbial biomass turnover rate increases with increasing microbial biomass (“Density-dependent turnover”), and 3) community carbon use efficiency decreases as microbial biomass increases toward an upper limit (“Density-dependent growth”). We ran a factorial combination of these hypotheses resulting in eight models for three different classes of model (linear decay, Michaelis-Menten decay, or reverse Michaelis-Menten decay), resulting in 24 models, 12 of which are presented or discussed in the related publication. Models were parameterized using values from previous studies with similar models (Wang et al. 2013, Wieder et al. 2014, Li et al. 2014, Georgiou et al. 2017, Hassink and Whitmore 1997) and ran to an approximate steady state (200 years) at six (6) different C input rates corresponding to 0.5, 1, 2, 4, 7, and 10 times the default input value. Further model details are available in the related publication. This archive contains output from three MAAT simulations, and scripts to run these simulations and process and plot the data. Simulations are labeled “lin”, “MM_highKm”, and “RMM_highKm” reflecting factorial runs for linear, Michealis-Menten, and reverse Michaelis-Menten models, respectively. This archive contains: • 3 R scripts prepended with “init_MAAT_” to initialize runs (1 for each simulation), • 1 csv file containing years over which to run simulations (“met_year.csv”), • 1 bash script (.bs) to run MAAT, • 6 XML files that are output from MAAT describing a run (2 for each simulation), • 3 model output csv files prepended by “out_” (1 for each simulation), and • 1 analysis R script for reproducing figures 3 and 4 in Craig et al. 2021. See included user guide (Craig_2021_modeling_archive_20210315.pdf) for file organization details.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗
DOE OSTI2018

Development and evaluation of a variably saturated flow model in the global E3SM Land Model (ELM) version 1.0

Abstract. Improving global-scale model representations of near-surface soil moisture and groundwater hydrology is important for accurately simulating terrestrial processes and predicting climate change effects on water resources. Most existing land surface models, including the default E3SM Land Model (ELMv0), which we modify here, routinely employ different formulations for water transport in the vadose and phreatic zones. Clark et al. (2015) identified a variably saturated Richards equation flow model as an important capability for improving simulation of coupled soil moisture and shallow groundwater dynamics. In this work, we developed the Variably Saturated Flow Model (VSFM) in ELMv1 to unify the treatment of soil hydrologic processes in the unsaturated and saturated zones. VSFM was tested on three benchmark problems and results were evaluated against observations and an existing benchmark model (PFLOTRAN). The ELMv1-VSFM's subsurface drainage parameter, fd, was calibrated to match an observationally constrained and spatially explicit global water table depth (WTD) product. Optimal spatially explicit fd values were obtained for 79 % of global 1.9∘ × 2.5∘ grid cells, while the remaining 21 % of global grid cells had predicted WTD deeper than the observationally constrained estimate. Comparison with predictions using the default fd value demonstrated that calibration significantly improved predictions, primarily by allowing much deeper WTDs. Model evaluation using the International Land Model Benchmarking package (ILAMB) showed that improvements in WTD predictions did not degrade model skill for any other metrics. We evaluated the computational performance of the VSFM model and found that the model is about 30 % more expensive than the default ELMv0 with an optimal processor layout. The modular software design of VSFM not only provides flexibility to configure the model for a range of problem setups but also allows for building the model independently of the ELM code, thus enabling straightforward testing of the model's physics against other models.

58 GEOSCIENCES↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗