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 73 records · Page 4

A Broadband Laboratory Study of the Seismic Properties of Cracked and Fluid-Saturated Synthetic Glass Media

In this paper, for better understanding of frequency dependence (dispersion) of seismic wave velocities caused by stress-induced fluid flow, broadband laboratory measurements were performed on a suite of synthetic glass media containing both equant pores and thermal cracks. Complementary forced oscillation, resonant bar, and ultrasonic techniques provided access to millihertz-hertz frequencies, ~1 kHz frequency, and ~1 MHz frequency, respectively. The wave speeds or effective elastic moduli and associated dissipation were measured on samples under dry, argon- or nitrogen-saturated, and water-saturated conditions in sequence. The elastic moduli, in situ permeability, and crack porosity inferred from in situ X-ray computed tomography all attest to strong pressure-induced crack closure for differential (confining-minus-pore) pressures <30 MPa, consistent with zero-pressure crack aspect ratios <4 × 10 -4 . The low permeabilities of these materials allow access to undrained conditions, even at subhertz frequencies. The ultrasonically measured elastic moduli reveal consistently higher shear and bulk moduli upon fluid saturation—diagnostic of the saturated-isolated regime. For a glass rod specimen, containing cracks but no pores, saturated-isolated conditions apparently persist to subhertz frequencies—requiring in situ aspect ratios (minimum/maximum dimension) <10 -5 . In marked contrast, the shear modulus measured at subhertz frequencies on a cracked glass bead specimen of 5% porosity, is insensitive to fluid saturation, consistent with the Biot-Gassmann model for the saturated-isobaric regime. Finally, the measured dispersion of the shear modulus approaches 10% over the millihertz-megahertz frequency range for the cracked and fluid-saturated media—implying that laboratory ultrasonic data should be used with care in the interpretation of field data.

58 GEOSCIENCES↗

Predicting Saturated Hydraulic Conductivity over Time for Degrading Saltstone Vault Concrete - 20376

At the Savannah River Site (SRS), low-level liquid tank waste is mixed with cementitious materials to produce a solid, monolithic waste form referred to as saltstone. The liquid saltstone mixture is pumped into massive concrete structures referred to as saltstone disposal units (SDUs) for curing and long-term storage. SDUs are vaults made of high performance concrete forming a barrier between saltstone and the surrounding environment. Predicting the long-term performance of SDU vault concrete, as well as the saltstone material itself, is critical to understanding the potential releases of radionuclides to the surrounding environment as described in the Saltstone Disposal Facility (SDF) performance assessment (PA). The hydraulic conductivity profile for SDU vault concrete was conservatively assumed to degrade linearly from the initial saturated hydraulic conductivity of the intact concrete over a long period of time. The initial hydraulic properties of saturated SDU vault concrete have been measured, and the final hydraulic properties of degraded vault concrete are assumed to be the same as the porous backfill material that will surround the vault after closure. The assumed final condition of the vault concrete is also conservative in that it provides a greater estimate of water penetration through the vault than would generally be expected. Therefore, insights that provide more realistic but conservative estimates as to how degradation processes may change the hydraulic properties of the SDU vault concrete, the time dependent relationship of the vault concrete hydraulic properties to degradation (and thus time), and the final hydraulic properties of the vault concrete provides useful additions to the SDF PA. This paper focuses on the assumption that the saltstone vault concrete saturated hydraulic conductivity would decrease linearly as a function of degradation time; this assumption resulted from a worst-case, weighted (by relative layer thickness), arithmetic average of intact and backfill material hydraulic conductivities representing flow parallel to a layered system. This research indicates that the effective saturated hydraulic conductivity of the degraded SDU vault concrete should be based on the relative thickness-weighted geometric mean of the intact and degraded layers hydraulic conductivities. This alternative, which is supported extensively in the relevant literature and evaluated in this paper, leads to a more realistic and defensible estimate of the effective saturated hydraulic conductivity in degraded vault concrete while remaining conservative. The impact of using the more realistic estimate of the effective saturated hydraulic conductivity for the degraded concrete decreases the effective saturated hydraulic conductivity by up to several orders of magnitude for early times (i.e., when the SDU vault concrete is assumed degraded via sulfate attack and carbonation). This revised estimate of the effective saturated hydraulic conductivity has been implemented in the recently revised SDF PA. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Inorganic semiconductor quantum dots as a saturated excitation (SAX) probe for subdiffraction imaging

The photoluminescence (PL) saturation of CdSe/ZnS core/shell inorganic semiconductor quantum dots (QDs) and its utility as a probe for saturated excitation (SAX) microscopy are reported. Under saturating excitation power densities, the PL signal was demodulated and recorded at harmonics of the fundamental frequency. For commercially available Qdot® 655 ITK™ QDs, the power density required to achieve saturation was dependent upon the local environment of the QDs. For QDs deposited and dried on a glass substrate, the excitation power density required for PL saturation was less than 1 kW/cm 2 . Compared to this, saturation of PL for QDs dispersed in water required an excitation power density greater than 200 kW/cm 2 . This observation is manifested as a limitation in the imaging of hydrated samples, as demonstrated for HeLa cells labelled with biotinylated-phalloidin followed by labelling with streptavidin-coated QDs. Finally as saturation affects the obtained spatial resolution in several imaging formats, including confocal imaging, the provided data will aid in obtaining the optimal spatial resolution when using QD probes to image biological samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Low-Power Absorption Saturation in Semiconductor Metasurfaces

Saturable optical elements lie at the cornerstone of many modern optical systems. Regularly patterned quasi-planar nanostructures—metasurfaces—are known to facilitate nonlinear optical processes. Such subwavelength semiconductor nanostructures can potentially serve as saturable components. In this work, we report on the intensity-dependent reflectance of femtosecond laser pulses from semiconductor metasurfaces with Mie-type modes, caused by the absorption saturation. Moreover, arrays of GaAs nanocylinders with magnetic dipole resonances in the spectral vicinity of the GaAs bandgap demonstrate a reduced saturation intensity and increased self-modulation efficiency, an order of magnitude higher than bulk GaAs or unstructured GaAs films. By contrast, the reflection modulation is shown to be negligible in the CW regime for the same average intensities, indicating that the process is not the result of temperature effects. Our work provides a novel idea for low-power saturable elements based on nonthermal nature of saturation. We conclude by devising a high-quality metasurface that can be used, in theory, to further reduce the saturation fluence below 50 nJ/cm 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical Analysis of the EPR Spectrum Saturation Effect with Account for Spectral Diffusion in a System with a Gaussian Distribution of Spin Resonance Frequencies

We have analyzed the effect of saturation of the EPR spectrum for the model situation when a set of resonance frequencies of spin packets is described by a Gaussian distribution and spectral diffusion, which is a random process without correlation, exists. We have obtained the form of the spectrum for arbitrary values of magnetic–resonant parameters: relaxation times for the longitudinal and transverse components of the magnetization vector of the system, the spectral diffusion rate, the dispersion of the Gaussian distribution of resonance frequency, and the magnetic induction of the microwave field. We have derived an analytic expression for the saturation curve of the EPR spectrum. Detailed analysis of the dependence of the spectrum saturation curve on the spin system parameters is performed for the optimal case when the microwave field frequency coincides with the average frequency of the Gaussian distribution for resonance frequencies. The value of the microwave field induction at which the peak of the saturation curve is attained has been determined for arbitrary values of the paramagnetic relaxation time, spectral diffusion rate, and dispersion of the Gaussian spin frequency distribution. We have formulated the algorithm for determining the spin–lattice relaxation time using the experimentally determined microwave field induction for which the saturation curve passes through the peak. It should be noted that numerical calculations have been performed for parameters typical of EPR spectroscopy. The general analytic expression derived for the shape of the spectrum in the saturation conditions and qualitative conclusions concerning the effect of spectral diffusion on the saturation of stationary spectra obtained in this study can also be used in other fields of spectroscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mapping Firn Saturation over Greenland using NASA’s Soil Moisture Active Passive Satellite

Mapping the spatial extent of recently identified englacial hydrological features (i.e., ice slabs and perennial firn aquifers) formed by meters-thick water-saturated firn layers over the percolation facies of the Greenland Ice Sheet using L-band microwave radiometry has recently been demonstrated. However, these initial maps are binary, and do not provide a parameter to estimate the spatial variability in the thickness and volumetric fraction of meltwater stored within the firn pore space. Here, we exploit enhanced-resolution vertical-polarization L-band brightness temperature (T^(B)_(V)) imagery (2015–2019) generated using observations collected over Greenland by NASA’s Soil Moisture Active Passive (SMAP) satellite and a simple two-layer L-band brightness temperature model. We map water-saturated firn layers via a “firn saturation” parameter, and interpret our results together with ice slab and perennial firn aquifer spatial extents, estimates of snow accumulation simulated via the Regional Atmospheric Climate Model (RACMOp2.3), and airborne radar surveys collected via NASA’s Operation IceBridge (OIB) campaigns. We find that variable firn saturation parameter values are mapped in lower snow accumulation ice slab areas in western, northern, and northeastern Greenland, where firn is colder and water-saturated firn layers seasonally refreeze as solid-ice. Higher firn saturation parameter values are mapped in higher snow accumulation perennial firn aquifer areas in southeastern, southern, and northwestern Greenland, where firn is near the melting point, and meters-thick water-saturated firn layers exist. Our results have implications for identifying expansive englacial reservoirs that store significant volumes of meltwater in locations that are vulnerable to meltwater-induced hydrofracturing and accelerated outlet glacier flow year-round.

Greenland ice sheet↗

Coupled complex Ginzburg–Landau systems with saturable nonlinearity and asymmetric cross-phase modulation

We formulate and study dynamics from a complex Ginzburg–Landau system with saturable nonlinearity, including asymmetric cross-phase modulation (XPM) parameters. Such equations can model phenomena described by complex Ginzburg–Landau systems under the added assumption of saturable media. When the saturation parameter is set to zero, we recover a general complex cubic Ginzburg–Landau system with XPM. We first derive conditions for the existence of bounded dynamics, approximating the absorbing set for solutions. We use this to then determine conditions for amplitude death of a single wavefunction. We also construct exact plane wave solutions, and determine conditions for their modulational instability. In a degenerate limit where dispersion and nonlinearity balance, we reduce our system to a saturable nonlinear Schrödinger system with XPM parameters, and we demonstrate the existence and behavior of spatially heterogeneous stationary solutions in this limit. Using numerical simulations we verify the aforementioned analytical results, while also demonstrating other interesting emergent features of the dynamics, such as spatiotemporal chaos in the presence of modulational instability. In other regimes, coherent patterns including uniform states or banded structures arise, corresponding to certain stable stationary states. For sufficiently large yet equal XPM parameters, we observe a segregation of wavefunctions into different regions of the spatial domain, while when XPM parameters are large and take different values, one wavefunction may decay to zero in finite time over the spatial domain (in agreement with the amplitude death predicted analytically). We also find a collection of transient features, including transient defects and what appear to be rogue waves, while in two spatial dimensions we observe highly localized pattern formation. While saturation will often regularize the dynamics, such transient dynamics can still be observed – and in some cases even prolonged – as the saturability of the media is increased, as the saturation may act to slow the timescale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Relationship between Anisotropy in Soil Hydraulic Conductivity and Saturation

Anisotropy in unsaturated hydraulic conductivity is saturation-dependent. Accurate characterization of soil anisotropy is very important in simulating flow and contaminant (e.g., radioactive nuclides in Hanford) transport. A recently developed tensorial connectivity-tortuosity (TCT) concept describes the hydraulic conductivity tensor of the unsaturated anisotropic soils as the product of a scalar variable, the symmetric connectivity tortuosity tensor, and the hydraulic conductivity tensor at saturation. In this study, the TCT model is used to quantify soil anisotropy in unsaturated hydraulic conductivity. The TCT model can describe different types of soil anisotropy; e.g., the anisotropy coefficient, C, can be monotonically increase or decrease with saturation and can vary from greater than unity to less than unity and vice versa. Soil anisotropy is independent of soil water retention properties and can be characterized by the ratio of the saturated hydraulic conductivities and the difference of the tortuosity-connectivity coefficients in two directions. ln(C) is linearly proportional to ln(Se) with Se being the effective saturation. The log-linear relationship between C and Se allows the saturation-dependent anisotropy to be determined using linear regression with the measurements of the directional hydraulic conductivities at a minimum of two water content levels, of which one may be at full saturation. The model was tested using measurements of directional hydraulic conductivities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Joint inversion of electromagnetic measurements for the determination of water saturation profiles in concrete structures

Highlights: • DC-electrical and dielectric data are combined to estimate concrete saturation degree. • A joint inversion approach of the electromagnetic measurements is proposed. • The sensitivity of the measurements to the saturation model parameters is analysed. • The new joint approach was developed and is applied to synthetic and real data. • The benefits of the joint approach over the inversion of one data type are highlighted. Water saturation profiles in concrete are essential to assess its durability and can be determined using non-destructive techniques, especially the electric and the capacitive methods. In this paper, we propose a new inversion scheme where both resistivity and permittivity measurements are inverted jointly to retrieve the saturation profile. The finite element method is used to model the measurements in 3D, the concrete having a saturation profile with depth, represented by a continuous model taking the form of a Weibull curve with four parameters. A non-linear least-squares optimization based on the Levenberg-Marquardt scheme is developed for the inversion of measurements. Results show that information gathered from both measurements enriches the reconstructed profile, leading to a more reliable estimation of saturation profiles. We believe that the joint inversion method herein developed could lead to the study of more complex phenomena, such as the coupled water-chloride ingress.

36 MATERIALS SCIENCE↗

Nonlinear saturation of the Weibel instability

The growth and saturation of magnetic fields due to the Weibel instability (WI) have important implications for laboratory and astrophysical plasmas, and this has drawn significant interest recently. Since the WI can generate a large magnetic field from no initial field, the maximum magnitudes achieved can have significant consequences for a number of applications. Hence, an understanding of the detailed dynamics driving the nonlinear saturation of the WI is important. This work considers the nonlinear saturation of the WI when counter-streaming populations of initially unmagnetized electrons are perturbed by a magnetic field oriented perpendicular to the direction of streaming. Previous works have found magnetic trapping to be important and connected electron skin depth spatial scales to the nonlinear saturation of the WI. The results presented in this work are consistent with these findings for a high-temperature case. However, using a high-order continuum kinetic simulation tool, this work demonstrates that when the electron populations are colder, a significant electrostatic potential develops that works with the magnetic field to create potential wells. The electrostatic field develops due to transverse flows induced by the WI and in some cases is strengthened by a secondary instability. This field plays a key role in saturation of the WI for colder populations. In conclusion, the role of the electrostatic potential in Weibel instability saturation has not been studied in detail previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗

Saturated fatty acid in the phospholipid monolayer contributes to the formation of large lipid droplets

The degree of saturation of fatty acid chains in the bilayer membrane structure is known to control membrane fluidity and packing density. However, the significance of fatty acid composition in the monolayers of lipid droplets (LDs) has not been elucidated. In this study, we noted a relationship between the size of LDs and the fatty acid composition of the monolayer. To obtain large LDs, we generated NIH3T3 cells overexpressing fat-specific protein 27 (FSP27). This induced the fusion of LDs, resulting in larger LDs in FSP27-overexpressing cells compared with LDs in control cells. Moreover, the lipid extracts of LDs from FSP27-overexpressing cells reconstituted large-droplet emulsions in vitro, implying that the lipid properties of LDs might affect the size of LDs. FSP27-overexpressing cells had more saturated fatty acids in the phospholipid monolayer of the LDs compared with control cells. To further investigate the effects of the degree of phospholipid unsaturation on the size of LDs, we synthesized artificial emulsions of a lipid mixed with distearoylphosphatidylcholine (DSPC, diC18:0-PC) and with dioleoylphosphatidylcholine (DOPC, diC18:1n-9-PC) and compared the sizes of the resulting LDs. The emulsions prepared from saturated PC had larger droplets than those prepared from unsaturated PC. Our results suggest that saturated fatty acid chains in phospholipid monolayers might establish the form and/or stability of large LDs. - Highlights: • The lipid extracts of larger LDs from FSP27 cells reconstructed large-droplet emulsions. • Isolated LDs from FSP27 cells had more saturated fatty acids in the phospholipid monolayer compared with the control. • Saturated fatty acids in the phospholipid monolayer are a factor in the formation of large emulsions.

60 APPLIED LIFE SCIENCES↗

Effects of grain size and small-scale bedform architecture on CO 2 saturation from buoyancy-driven flow

Small-scale (mm-dm scale) heterogeneity has been shown to significantly impact CO 2 migration and trapping. To investigate how and why different aspects of small-scale heterogeneity affect the amount of capillary trapping during buoyancy-driven upward migration of CO 2 , we conducted modified invasion percolation simulations on heterogeneous domains. Realistic simulation domains are constructed by varying two important aspects of small-scale geologic heterogeneity: sedimentary bedform architecture and grain size contrast between the matrix and the laminae facies. Buoyancy-driven flow simulation runs cover 59 bedform architecture and 40 grain size contrast cases. Simulation results show that the domain effective CO 2 saturation is strongly affected by both grain size and bedform architecture. At high grain size contrasts, bedforms with continuous ripple lamination at the cm scale tend to retain higher CO 2 saturation than bedforms with discontinuous or cross lamination. In addition, the “extremely well sorted” grain sorting cases tend to have lower CO 2 saturation than expected for cross-laminated domains. Finally, both a denser CO 2 phase and greater interfacial tension increase CO 2 saturation. Again, variation in fluid properties seems to have a greater effect on CO 2 saturation for cross-laminated domains. This result suggests that differences in bedform architecture can impact how CO 2 saturation values respond to other variables such as grain sorting and fluid properties.

58 GEOSCIENCES↗

Decoupling of redox processes from soil saturation in Arctic tundra

Permafrost thaw in warming Arctic landscapes alters hydrology and saturation-driven biogeochemical processes. Models assume that aerobic respiration occurs in drained soils while saturated soils support methanogenesis; however, saturated soils maintain redox gradients that host a range of anaerobic metabolisms. We evaluated how redox potential and redox-active solutes vary with soil moisture in the active layer of permafrost-affected acidic and non-acidic tundra hillslopes. Oxidizing conditions persisted in highly permeable organic horizons of both unsaturated tussock tundra and saturated wet sedge meadows. Redox potential decreased with depth in all soils as increasing soil bulk density restricted groundwater flow and oxygen diffusion. High concentrations of dissolved iron, phosphate, and organic carbon coincided with redox boundaries below the soil surface in acidic tundra, indicating active iron redox cycling and potential release of adsorbed phosphate during iron (oxyhydr)oxide dissolution. In non-acidic tundra, weatherable minerals affected nutrient dynamics more than redox-driven iron cycling, especially in low-lying, saturated areas where thaw reached mineral soils. The role of thaw depth and the ability of saturated soils to maintain oxidizing conditions in organic surface layers highlight the importance of soil physical properties and hydrology in predicting biogeochemical processes and greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗

Quasiclassical evaluation of gluon saturation induced helicity effects

At sub-eikonal order in the high-energy limit, helicity dependences are generally contained in the transverse components of the gluon field. In the gluon saturation regime, novel helicity effects arise from the nonlinear interaction between the eikonal-order longitudinal gluon field and the sub-eikonal-order transverse gluon field. We derive this saturation induced helicity-dependent field both by solving the classical Yang-Mills equations and through direct diagrammatic calculations. Our analysis shows that the saturation induced helicity effect is intrinsically a two-particle (or multiparticle) correlation effect, rather than a single-particle distribution effect. Furthermore, we evaluate this effect in the context of double-spin asymmetry for incoherent diffractive dijet production in longitudinally polarized electron-nucleus collisions, using a simplified helicity-extended McLerran-Venugopalan model. We find that the saturation induced helicity effect further suppresses the back-to-back peak in the dijet azimuthal angle correlation, with a magnitude comparable to that of the direct helicity effect. This gluon saturation induced helicity effect may offer a novel avenue to probe gluon saturation in polarized collisions. Published by the American Physical Society 2025

Li, Ming (ORCID:0000000195509691)↗

Real-time deep-learning inversion of seismic full waveform data for CO 2 saturation and uncertainty in geological carbon storage monitoring

Deep-learning inversion has recently drawn attention in geological carbon storage research due to its potential of imaging and monitoring carbon storage in real time, significantly improving efficiency and safety of carbon storage operations. We present a deep-learning full waveform inversion method that after the neural network has been trained can image CO 2 saturation and its uncertainty in real time. Our deep-learning inversion method is based on the U-Net architecture with the neural network trained on pairs of synthetic seismic data and CO 2 saturation models. Accordingly, our training establishes a mapping relationship between seismic data and CO 2 saturation models and once fully trained directly estimates CO 2 saturation as a function of subsurface location. We further quantify uncertainties of CO 2 saturation estimates using the Monte Carlo dropout method and a bootstrap aggregating method. For this proof-of-concept study, the CO 2 training models and data are derived from the Kimberlina 1.2 model, a hypothetical 3D geological carbon storage model that is constructed based on various geological and hydrological data from the Southern San Joaquin Basin, California. We perform deep-learning inversion experiments using noise-free and noisy training and test data sets and compare the results. Our modelling experiments show that (1) the deep-learning inversion can estimate 2D distributions of CO 2 fairly well even in the presence of Gaussian random noise and (2) both CO 2 saturation imaging and uncertainty quantification can be done in real time. Our results suggest that the deep-learning inversion method can serve as a robust real-time monitoring tool for geological carbon storage and/or other time-varying reservoir/aquifer properties that result from injection, extraction, and/or other subsurface transport phenomena.

58 GEOSCIENCES↗

Characterization of Thermal, Hydraulic, and Gas Diffusion Properties in Variably Saturated Sand Grades

Core Ideas This study provides physical, hydraulic, and thermal properties of five sand grades. An extended three‐region thermal conductivity model is proposed. A gas diffusivity‐based variable tortuosity parameter is examined. Useful parametric functions and descriptive models were tested. Detailed characterization of partially saturated porous media is important for understanding and predicting vadose zone transport processes. While basic properties (e.g., particle‐ and pore‐size distributions and soil‐water retention) are, in general, essential prerequisites for characterizing most porous media transport properties, key transport parameters such as thermal conductivity and gas diffusivity are particularly important to describe temperature‐induced heat transport and diffusion‐controlled gas transport processes, respectively. Despite many experimental and numerical studies focusing on a specific porous media characteristic, a single study presenting a wide range of important characteristics, together with the best‐performing functional relationships, can seldom be found. This study characterized five differently textured sand grades (Accusand no. 12/20, 20/30, 30/40, 40/50, and 50/70) in relation to physical properties, water retention, hydraulic conductivity, thermal conductivity, and gas diffusivity. We used measured basic properties and transport data to accurately parameterize the characteristic functions (particle‐ and pore‐size distributions and water retention) and descriptive transport models (thermal conductivity, saturated hydraulic conductivity, and gas diffusivity). An existing thermal conductivity model was improved to describe the distinct three‐region behavior in observed thermal conductivity–water saturation relations. Applying widely used parametric models for saturated hydraulic conductivity and soil‐gas diffusivity, we characterized porous media tortuosity in relation to grain size. Strong relations among average particle diameter, characteristic pore diameter from soil‐water retention measurements, and saturated hydraulic conductivity were found. Thus, the results of this work are useful toward better describing, linking, and predicting mass transfer and pore network properties in variably saturated porous media.

Chamindu Deepagoda, T. K. K.↗

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

This model package report (MPR) documents the development of the integrated vadose and saturated zone flow and transport model developed for the performance assessment of Trenches 31 and 34 in the low-level burial grounds (LLBGs). This modeling capability is intended for use in addressing the analysis requirements outlined in DOE O 435.1, Chg 1, Radioactive Waste Management1. The overall objective of the modeling effort is to provide a basis for making informed disposal decisions pertinent to Trenches 31 and 34. The purpose of the MPR is to document the development of the three-dimensional numerical vadose zone and saturated flow and transport model test case and evaluate its adequacy to support the LLBG performance assessment. The purpose is not to present results for DOE 435.1 decision making. The use of the model to perform base case and sensitivity analysis for the LLBG performance assessment, including inputs and results, is documented separately in subsequent environmental calculation files. This report discusses the development and translation of the conceptual model for flow and contaminant transport into the LLBG performance assessment three-dimensional numerical flow and transport model evaluated using the Subsurface Transport Over Multiple Phases (STOMP©) simulator. The development of representative geologic framework is described along with the implementation of waste release models used to represent contaminant release from waste disposed in the trenches. The report also provides the technical basis for specific model parameters and boundary conditions, along with description of modeling assumptions. This MPR includes certain calculations that are necessary to demonstrate the soundness of the model. Results provided by the model include vadose zone and saturated zone flow fields, and estimates of the possible future concentration in groundwater of technetium-99 and iodine-129 as example test cases. As an evaluation of the test cases, the model estimates of current vadose conditions are compared to available and analogous field and laboratory data, and compared to the results from the original performance assessment model (WHC-EP-0645, Performance Assessment for the Disposal of Low-Level Waste in the 200 West Area Burial Grounds 2 ). The features, events, and processes applicable to vadose zone and saturated zone flow and transport model are identified, and representative initial estimates for various parameters are documented. Note that the parameter estimates presented in this MPR are for illustration purposes, and may or may not reflect values selected for eventual performance assessment. Several key topical discussions (i.e., basis for recharge estimates, basis for vadose zone modeling, basis for saturated zone model development, and calibration) are included, which serve as the groundwork for confidence building for the groundwater pathway modeling and results. Numerical simulation results based on an example test case are included to illustrate the use of the combined saturated-unsaturated model in performance assessment calculations. Sensitivity and uncertainty results are not included in this MPR. Those results will be included in future environmental calculation files. The inclusion of the trapezoidal trench geometry and construction details in the finite difference grid introduces some gross simplifications regarding the trench and liner systems. The model test case presented in this MPR includes the assumption that the trench liner system does not affect flow through or around the trenches after the liner system is assumed to fail. This MPR also includes results of other test cases that involve alternate assumptions about how to incorporate the hydraulic effects of the trenches into the vadose zone of the model. These alternate cases do not attempt to account for the presence of the liner system after its assumed failure either. None of these cases is considered to be the base case at this time. Analysis and alternate cases that attempt to account for the presence of the liner system in greater detail, and its effect on flow through and around it, are to be documented in subsequent environmental calculation files. Depending on the results of those analyses and alternate cases, the eventual base case may involve more detailed inclusion of the effects of the liner system hydraulics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗