Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volumetric”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Object-space optimization of tomographic reconstructions for additive manufacturing

Volumetric 3D printing motivated by computed axial lithography enables rapid printing of homogeneous parts but requires a high dimensionality gradient-descent optimization to calculate image sets. Here we introduce a new, simpler approach to image-computation that algebraically optimizes a model of the printed object, significantly improving print accuracy of complex parts under imperfect material and optical precision by improving optical dose contrast between the target and surrounding regions. Quality metrics for volumetric printing are defined and shown to be significantly improved by the new algorithm. The approach is extended beyond binary printing to grayscale control of conversion to enable functionally graded materials. The flexibility of the technique is digitally demonstrated with realistic projector point spread functions, printing around occluding structures, printing with restricted angular range, and incorporation of materials chemistry such as inhibition. Finally, simulations show that the method facilitates new printing modalities such as printing into flat, rather than cylindrical packages to extend the applications of volumetric printing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation

Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.

improving volumetric energy density↗
NASA NTRS1944

Dynamics of the inlet system of a four-stroke engine

Tests were run on a single-cylinder and a multicylinder four-stroke engine in order to determine the effect of the dynamics of the inlet system upon indicated mean effective pressure. Tests on the single-cylinder engine were made at various speeds, inlet valve timings, and inlet pipe lengths. These tests indicated that the indicated mean effective pressure could be raised considerably at any one speed by the use of a suitably long inlet pipe. Tests at other speeds with this length of pipe showed higher indicated mean effective pressure than with a very short pipe, although not so high as could be obtained with the pipe length adjusted for each speed. A general relation was discovered between optimum time of inlet valve closing and pipe length; namely, that longer pipes require later inlet valve closing in order to be fully effective. Tests were also made on three cylinders connected to a single pipe. With this arrangement, increased volumetric efficiency at low speed was obtainable by using a long pipe, but only with a sacrifice of volumetric efficiency at high speed. Volumetric efficiency at high speed was progressively lower as the pipe length was increased.

ENGINES-CHRYSLER↗
DOE OSTI2021

Strategies to achieve high productivity, high conversion, and high yield in yeast fermentation of algal biomass hydrolysate

The conversion of carbohydrates in biomass via fermentation is an important component of an overall strategy to decarbonize the production of fuels and chemicals. Owing to the cost and resources required to produce biomass hydrolysates, the economic and environmental sustainability of these fermentation processes requires that they operate with high yields, sugar conversion, and productivity. Immobilized-cell technology in a continuous bioprocess can achieve significantly higher volumetric productivities than is possible from standard batch fermentation using free cells. Here, we demonstrate approaches for improvement of ethanol yield from algal hydrolysates and a mock hydrolysate medium. Saccharomyces cerevisiae was immobilized in alginate and incorporated into a two-column immobilized cell reactor system. Furthermore, the yeast quorum-sensing molecule, 2-phenylethanol, was added to improve ethanol yield by restricting growth and diverting sugar to ethanol. The bioreactor system could achieve high ethanol volumetric productivity (>20 g/ Lreactor ·h) and high glucose conversion (>99%) in mock hydrolysate, while the addition of 0.2% 2-phenylethanol resulted in 4.9% higher ethanol yield. With an algal hydrolysate of <10 g/L sugar, the ethanol volumetric productivity reached 9.8 g/ Lreactor ·h, and the addition of 0.2% 2-phenylethanol increased the ethanol yield by up to 7.4%. These results demonstrate the feasibility of novel strategies to achieve sustainability goals in biomass conversions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2021

A novel lung-inspired 3D-printed desiccant-coated heat exchanger for high-performance humidity management in buildings

Offering an independent humidity management method for buildings, desiccant-coated heat exchangers (DCHXs) are deemed a promising approach to improve the overall energy efficiency of air conditioning (AC) systems. State-of-the-art DCHXs, however, are bound with conventional HX topologies either providing limited desiccant-air interfacial areas or introducing excessive pressure drop penalties. In this paper a novel 3D-printed DCHX concept inspired by the bronchi arrangement of a human lung is introduced to address the shortcomings inherent in existing DCHX designs. The proposed lung-inspired DCHX utilizes two intertwined bicontinuous flow networks enabling highly efficient heat and mass transfer characteristics for augmented adsorption and regeneration processes at low pressure drop penalties. While the first network evenly distributes an incoming air stream through the entire volume of the lung-inspired DCHX, the second network volumetrically splits a cooling water stream within and through the first network. Effects of various parameters including air flow rate, outdoor air humidity ratio, and regeneration temperature on dehumidification performance and energy efficiency of the proposed lung-inspired DCHX were investigated. Experimental results indicated the proposed lung-inspired 3D-printed DCHX outperforms existing DCHX systems by demonstrating an excellent balance between a high volumetric adsorption rate and a low pressure drop penalty. The volumetric adsorption rate of the proposed lung-inspired DCHX technology is 54.8 g/m 3 -s, a more than two-fold improvement compared with state-of-the-art DCHX systems. Additionally, the lung-inspired DCHX offers high thermal energy efficiency of 56% at a low regeneration temperature of 40 °C. Therefore, the proposed lung-inspired 3D-printed DCHX offers a new solid-desiccant-based air dehumidification pathway for next-generation high-performance AC systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Experimental investigation on the coolability of nuclear reactor debris beds using seawater

During the Fukushima nuclear reactor accident, seawater was injected into the reactor core to cool the decay heat from the heat generating porous debris bed. However, the impact of dissolved salts in water on the coolability of a debris bed is not well understood. This paper is one of the first works to investigate the cooling mechanisms for debris bed using seawater. An experimental system was built where the packed bed test section was volumetrically heated via a high frequency induction heater. The temperature along the axis of the test section was measured using a high resolution temperature sensor based on Optical Frequency Domain Reflectometry. Experimental results show that the dryout heat flux increases with an increase in concentration, with pure water having a dryout volumetric heat flux of 1.8 MW/m and salt water of concentration 7% having a dryout volumetric heat flux of 2.8 MW/m. However, this enhanced coolability using salt water is short term, as the crystallization fouling mechanism becomes significant over time and leads to gradual temperature excursion in the bed. The fouling mechanism eventually leads to a cross section of the debris bed to plug with the deposited salts, and this leads to a pressure build up in the lower regions of the debris bed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Octane Modeling of Isobutanol Blending into Gasoline

Thirty-four gasoline blendstocks for oxygenate blending were used to create finished gasoline blends with isobutanol content of 12.5 volume percent (vol. %) and 16 vol. %. The gasoline blendstocks and finished fuels were analyzed for octane number (research [RON] and motor [MON]) to determine the effect of blending isobutanol. Volumetric and molar linear blending models were developed to predict finished fuel RON and MON, starting from the properties and composition of the gasoline blendstocks and isobutanol. Results show the molar blending model provided a better fit for the experimental data than the volumetric blending model. The volumetric model was further improved by adding nonlinear terms, improving the error to within ~1 ON. Gasoline blendstock properties impacted the finished fuel RON/MON, with paraffins having a synergistic effect with isobutanol and olefins and aromatics having an antagonistic effect.

10 SYNTHETIC FUELS↗
DOE OSTI2022

High Energy Density Aqueous Flow Battery Utilizing Extremely Stable, Branching-Induced High-Solubility Anthraquinone near Neutral pH

An anthraquinone featuring a chiral carboxylate-capped methyl-branched side chain with an ether linkage, 2,2'-((9,10-dioxo-9,10-dihydroanthracene-2,6-diyl)bis(oxy))dipropionic acid (2,6-D2PEAQ), was synthesized and evaluated for use in aqueous redox flow batteries. It was found to have an extraordinary solubility of 2 M (4 M electrons), corresponding to a theoretical volumetric capacity of 107.2 Ah/L for the negative electrolyte, which is 10 times that of its unbranched counterpart. The 2,6-D2PEAQ molecule demonstrated stability against thermal 2 decomposition and was extremely stable under cell cycling conditions. Here, a capacity fade rate of 0.02% day over 14 days was demonstrated in a 1.1 M 2,6-D2PEAQ nearly capacity-balanced cell when paired with a ferro-/ferricyanide posolyte at pH 7. Compared to other aqueous redoxactive organic molecules, its demonstrated fade rate is lower than that of any molecule with a demonstrated volumetric capacity of ≥55 Ah/L, and its volumetric capacity is greater than that of any molecule with a demonstrated fade rate of ≤0.5%/day.

25 ENERGY STORAGE↗
DOE OSTI2020

Light-Driven Expansion of Spiropyran Hydrogels

The incorporation of molecular switches in organic structures is of great interest in the chemical design of stimuli-responsive materials that mimic the complex functions of living systems. Merocyanine dyes that convert to spiropyran moieties upon exposure to visible light have been extensively studied as they can be incorporated in hydrated covalent networks that will expel water when this conversion occurs and induce a volumetric shrinkage. We report here on a sulfonatebased water soluble photoswitch that in contrast to the well-known system triggers a volumetric expansion in hydrogels upon exposure to photons. Contraction is in turn observed under dark conditions in a highly reversible manner. The novel behavior of the photoswitch incorporated in the covalent network was also predicted by coarse grain simulations of the system’s chemical structure. Using pH control and polymeric structures that differ in lower critical solution temperature we were able to develop hydrogels with highly tunable volumetric expansion. Furthermore, the novel molecular function of the systems developed here led to materials with the negative phototaxis observed in plants, and can expand the potential use of hydrogels as sensors, soft robots, and actuators.

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

Mechanically Robust Bismuth-Embedded Carbon Microspheres for Ultrafast Charging and Ultrastable Sodium-Ion Batteries

Advancements in the development of fast-charging and long-lasting microstructured alloying anodes with high volumetric capacities are essential for enhancing the operational efficiency of sodium-ion batteries (SIBs). These anodes, however, face challenges such as declined cyclability and rate capability, primarily due to mechanical degradation reduced by significant volumetric changes (over 252%) and slow kinetics of sodium-ion storage. Herein, we introduce a novel anode design featuring densely packed bismuth (Bi) embedded within highly conductive carbon microspheres to overcome the aforementioned challenges. Remarkably, the high loading Bi anode within carbon microspheres with a high tap density of 2.59 g cm -3 possesses significant mechanical strength exceeding 590 MPa and limits volume swelling of only 10.9% post-sodiation. This anode demonstrates a high volumetric capacity (908.3 mAh cm -3 ), ultrafast chargeability (200 A g -1 , full charge/discharge in just 5.5 s), and outstanding cyclability over 12,000 cycles and maintains exceptional cycling stability even at -30 °C. The full cell paired with a Na 3 V 2 (PO 4 ) 3 cathode retains over 80% capacity after 600 cycles at 36 C, demonstrating a remarkable rate capability of 126 C (full charge/discharge in 28.6 s). Our comprehensive experimental evaluations and chemo-mechanical simulations shed light on the mechanisms underpinning the anode's superior performance. In conclusion, this development marks a significant advancement in the design of durable and fast-charging anodes for high-performance SIBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗
DOE OSTI2022

Electric charge build-up and dissipation at pulsed streamer corona by EFISH and probe measurements

The aim of this work is to study the dynamics of volumetric electric charge deposited by a single pin electrode streamer corona in atmospheric air at ambient conditions. Recent studies show that, at pulse-periodic operation, electric charge remains in surrounding air after the first discharge. This residual charge significantly impacts the electric field distribution and morphology of discharges at subsequent pulses. In this study, the residual electrical potential/field was measured during and after dissipation of a pulsed 80kV streamer corona of alternating polarity. This was accomplished using electrostatic probes and the EFISH method (Electric Field Induced Second Harmonic generation). The residual volumetric charge of a positive pulse polarity streamer corona occupies a zone up to 80mm from the electrode. The subsequent application of a negative polarity pulse leads to partial neutralization of the previous charge and the generation of a region with a high-amplitude electric field. Here, the propagation of following positive polarity pulses is limited by the residual positive volumetric charge. The EFISH method reveals the magnitude of the electric field at discharge propagation and electrostatic probes inform the polarity. Intensified imaging acquisition was used to evaluate the spatial distribution and dynamics of corona discharge during the alternating polarity waveform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical Simulations of Carbon Dioxide Storage Efficiency in Heterogeneous Reservoir Models

The U.S. Department of Energy’s National Energy Technology Laboratory (DOE-NETL) has been developing methods and tools (the online Carbon Dioxide Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN) tool) to estimate carbon dioxide (CO2) storage potential in subsurface reservoirs. The CO2 storage efficiency terms are input in the tool to calculate storage potential in targeted reservoirs. In this effort, two CO2 storage efficiency terms were evaluated: volumetric displacement ( E V ) and microscopic displacement ( E d ). The first term deals with efficiency of CO2 propagation into an accessible reservoir volume, while the second term evaluates effectiveness of native fluid displacement with CO2. The interpreted well logs and core sample measurements were applied to create the heterogeneous reservoir models including geostatistical realizations of porosity and intrinsic permeability fields. Supercritical CO2 was injected over the course of 30 years into brine-saturated reservoir models for clastics, limestone, and dolomite lithologies and deltaic fluvial, aeolian, shallow marine, and reef depositional environments by means of varying reservoir parameters and injection scenarios. The reservoir models providing vertically heterogeneous petrophysical properties and designated as “layered reservoir models” (with homogeneous parameters along each layer of the model) were not determined to be a transition between the homogeneous and heterogeneous models in respect to storage efficiency. Another finding shows that high-efficiency factors do not necessarily mean increased CO2 storage; they rather indicate that the available volume and pore space are more fully utilized. The CO2 storage efficiency factors were evaluated dynamically at the select time points using P 10 ‐ P 50 ‐ P 90 percentiles. The results of this study show that the P 10 ‐ P 90 distribution for volumetric efficiency is wider when compared to the microscopic efficiency. It was found that where dominant buoyancy forces drive the plume to the top of a target formation, the volumetric efficiency is low. Tighter sandstone and carbonate formations show prevalence of capillary forces and better utilization of reservoir volume.

Myshakin, Evgeniy M.↗
DOE OSTI2022

Understanding the tool influence function during sub-aperture belt-on-wheel glass polishing

The tool influence function (TIF) during sub-aperture belt-on-wheel polishing has been evaluated as a function of various process conditions (belt use/wear, dwell time, displacement, belt velocity, and wheel modulus and diameter) on fused silica glass workpieces using C e O 2 polishing media. TIF spots are circular or elliptical in shape with a largely flat bottom character. Furthermore, the volumetric removal rate varies significantly with belt use (or wear), stabilizing after ~<!-- ~ --> 15 m i n of use. A modified Preston model, where the pressure dependence is adjusted using a different scaling of the wheel modulus ( E w 0.5 ), largely predicts the volumetric removal rate over the range of process conditions evaluated. The relatively high volumetric removal rate of 30 -<!-- - --> 60 m m 3 / h using a fixed C e O 2 -in-resin-host belt offers a rapid, and hence, more economical, initial polish of aspheric and freeform optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗