Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permeability models”

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 379 records · Page 21

Comparative Assessment of U-Net-Based Deep Learning Models for Segmenting Microfractures and Pore Spaces in Digital Rocks

Segmentation of high-resolution X-ray microcomputed tomography (µCT) images is crucial in digital rock physics (DRP), affecting the characterization and analysis of microscale phenomena in the porous media. The complexity of geological structures and nonideal scanning conditions pose significant challenges to conventional image segmentation approaches. Motivated by the recent increasing popularity of deep learning (DL) techniques in image processing, this work undertakes a comparative study of DL models, specifically U-Net and its variants, for segmenting multiple targets with distinguished features in digital rocks, including discrete fracture networks (DFNs), pore spaces, and solid rock. Particularly, DFNs have a smaller volumetric fraction over others, bringing in a substantial challenge of imbalanced segmentation. The primary focus is to evaluate the architecture and feature enhancement strategies of various DL models, including U-Net, attention U-Net, residual U-Net, U-Net++, and residual U-Net++. The models were designed as 2.5D, utilizing a central 2D image and its two adjacent upper and lower 2D images as input to provide a pseudo-3D context. In addition, because the ground truth of segmentation was unknown for real-world digital rocks, we created a benchmark data set following the inverse operations of segmentation. The data synthesis started from the label images (i.e., solid rock, pore spaces, and DFNs), followed by simulating partial volume blurring, adding random background noise, and introducing ring artifacts to mimic real raw X-ray µCT images. The data set, which included various rock types (i.e., sandstone and artificial data), scanning resolution, and magnitudes of noise and artifacts, was divided into training and testing data sets with a 90% and 10% ratio, respectively. Moreover, in addition to the conventional pixel-wise evaluation metrics, the physics-based metric of the lattice-Boltzmann method (LBM) simulated permeability provided more comprehensive assessments. The results demonstrated that the residual connections, nested architectures, and redesigned skip connections contribute to the model performance and give the residual U-Net++ the highest accuracy. The improvements were mainly on the boundaries and small targets, especially the DFNs, which dominate the interconnectivity and therefore affect the permeability greatly. This study also rigorously evaluated the efficiency and generalization of each model, demonstrating that the sophisticated architectures achieved excellent practicability and maintained robust performance on completely unseen data, ensuring their suitability for diverse and challenging DRP applications.

58 GEOSCIENCES↗

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented here. This model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. Finally, the results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗

EXPERIMENTAL AND MODELING STUDY OF MOLECULAR TRANSPORT PHENOMENA IN COMPLEX LIQUID MIXTURES USING COVALENT ORGANIC FRAMEWORK MEMBRANES

Currently, 10-15% of the world’s energy is consumed by chemical separations, and more than 80% of that is used to purify organic liquids and recover rare earth elements and critical minerals. Membrane nanofiltration presents an energy-efficient, cost-effective, and eco-friendly alternative to current separation technologies. These complex liquid environments require high- performance, chemically-resistant membrane materials. Recently discovered Covalent Organic Frameworks (COFs) possess desirable properties as membrane materials for complex liquid separations. COFs are highly crystalline, chemically and thermally stable, with tunable size and charge. COFs, especially two-dimensional highly crystalline COFs, possess narrow pore-size distribution and controlled porous structures. Two-dimensional COFs are also resistant to swelling, another feature beneficial in complex liquid separations. Molecular interactions in complex liquid systems often dictate final membrane performance, resulting in a discrepancy between designed and apparent COF properties. Understanding and resolving this discrepancy is critical in utilizing COF membranes as a viable alternative to energy-intensive separations. The primary objective of this thesis was to use a commercially available COF TpPa-1 as a platform to investigate how mixed solvent environment affects COF membrane performance via experimental and molecular modeling studies. Specifically, the effects of solution pH, solvent- solvent-solute interactions in mixed solvents, and COF chemistry were carefully examined on apparent TpPa-1 pore size and permeability and targeted solute rejection of COF membranes for neat and mixed solvents. Experimental COF filtration performance was compared and corroborated with modeling predication. Findings from this study will provide guidance for future COF design, synthesis, and desired functional COF membrane performance.

Barnes, Anastasia↗

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗

Sorafenib derivatives-functionalized gold nanoparticles confer protection against tumor angiogenesis and proliferation via suppression of EGFR and VEGFR-2

Sorafenib is a multi-kinase inhibitor that has been highlighted as a tumor suppressor due to its anti-proliferative and anti-angiogenic properties, whereas the clinical application of Sorafenib is restricted by the side effects it may cause. The past decade has witnessed the development of a series of sorafenib derivatives to improve the clinical performance of sorafenib. Gold nanoparticles (AuNPs) have been widely utilized in drug delivery systems due to their unique properties, including biocompatible nature, simple preparation, and easy surface modification. Herein, this study is aimed to investigate the anti-tumor effects of new sorafenib derivatives-capped gold nanoparticles (AuNPs-New Sor) in tumor formation and metastasis as well as the underlying mechanisms. Initially, new sorafenib derivatives were constructed and combined with AuNPs to form AuNPs-New Sor, and the properties of synthesized AuNPs-New Sor were identified in a mouse model of tumorigenesis. The effect of AuNPs-New Sor on tumor vascular normalization was investigated by assessing vascular permeability and perfusion rate. Next, we evaluated the effect of AuNPs-New Sor on migration and viability of tumor cells and human umbilical vein endothelial cells (HUVECs) as well as on HUVEC angiogenesis in vitro. A melanoma mouse model was further established for in vivo substantiation of the anti-tumor effect of AuNPs-New Sor. According to the results, AuNPs could deliver new sorafenib derivatives into tumor tissues and downregulate the expression of epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor-2 (VEGFR-2), thereby suppressing tumor migration, EMT, and angiogenesis in vitro. In addition, AuNPs-New Sor displayed competitive anti-tumor activities in vivo. Taken together, AuNPs-New Sor may attenuate tumor development and angiogenesis through downregulation of EGFR and VEGFR-2.

60 APPLIED LIFE SCIENCES↗

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗

The impact of physiologically relevant temperatures on physical properties of thylakoid membranes: a molecular dynamics study

Thylakoid membranes are energy-converting membranes with a unique lipid composition. Though the membranes are primarily composed of proteins, their photosynthetic function is strongly influenced by the lipid constituents. Here we characterize, with molecular dynamics (MD) simulations, lipid bilayers with compositions representative of plant thylakoid membranes. We determine, in a wide range of temperatures, the physical parameters of the model membranes which are relevant for the photosynthetic function. We found a marked impact of temperature on membrane permeability due to a combination of increased compressibility and curvature of the membrane at elevated temperatures. With increasing temperatures, we observed increasingly smeared transmembrane density profiles of the membrane forming lipid headgroups predicting increased membrane flexibility. The diffusion coefficient of the lipids increased with temperature without apparent specificity for lipid species. Instead of a comprehensive experimental dataset in the relevant temperature range, we quantitatively compared and validated our MD results with MD simulations on a dipalmitoylphosphatidylcholine model system.

59 BASIC BIOLOGICAL SCIENCES↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗

Geochemical narrowing of cement fracture aperture during multiphase flow of supercritical CO2 and brine

For carbon capture and storage operations, wells are used as a necessary conduit for injecting CO2 at depth, but they can also act as leakage conduits if the cement seals are compromised. The specific objective of this research was to investigate the nature of self-healing of a fracture within cement under multiphase flow of CO2 and brine, and to compare this to a newly developed model. This was accomplished through a laboratory investigation, flowing reacting supercritical CO2 and brine through a cement fracture with multiphase flow. The influent end of the fracture was degraded as evidenced by the formation of cracks across the surface. At the effluent end of the fracture, the initial aperture of 137 µm was reduced to 50 µm. This reduction by 87 µm was similar to results from a recently developed model tested in this study, which predicted an aperture reduction of 80 µm. Our results demonstrated that the model accurately predicted the reactions which occurred in this multiphase flow system.

wellbore, cement, fracture, multiphase-flow, super↗

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗

Inhomogeneous diffusion-limited aggregation

It is demonstrated here that inhomogeneous diffusion-limited aggregation (DLA) model can be used to simulate viscous fingering in a medium with inhomogeneous permeability and homogeneous porosity. The medium consists of a pipe-pore square-lattice network in which all pores have equal volume and the pipes have negligible volume. It is shown that fluctuations in a DLA-based growth process may be tuned by noise reduction, and that fluctuations in the velocity of the moving interface are multiplicative in form.

Selinger, Robin Blumberg↗

Iron abundance and magnetic permeability of the moon

A larger set of simultaneous data from the Apollo 12 lunar surface magnetometer and the Explorer 35 Ames magnetometer are used to construct a whole-moon hysteresis curve, from which a new value of global lunar permeability is determined to be mu = 1.012 + or - 0.006. The corresponding global induced dipole moment is 2.1 x 10 to the 18th power gauss-cu cm for typical inducing fields of .0001 gauss in the lunar environment. From the permeability measurement, lunar free iron abundance is determined to be 2.5 + or - 2.0 wt. %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior: a free iron/orthopyroxene lunar composition and a free iron/olivine composition. The overall lunar total iron abundance is determined to be 9.0 + or - 4.7 wt. %. Other lunar models with a small iron core and with a shallow iron-rich layer are discussed in light of the measured global permeability.

Parkin, C. W.↗

Iron abundance and magnetic permeability of the moon

A larger set of simultaneous data from the Apollo 12 lunar surface magnetometer and the Explorer 35 Ames magnetometer are used to construct a whole-moon hysteresis curve, from which a new value of global lunar permeability is determined to be mu = 1.012 + or - 0.006. The corresponding global induced dipole moment is 2.1 times ten to the eighteenth power gauss-cu cm for typical inducing fields of one ten-thousandth gauss in the lunar environment. From the permeability measurement, lunar free iron abundance is determined to be 2.5 + or - 2.0 wt %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior: a free iron/orthopyroxene lunar composition and a free iron/olivine composition. The overall lunar total iron abundance is determined to be 9.0 + or - 4.7 wt %. Other lunar models with a small iron core and with a shallow iron-rich layer are discussed in light of the measured global permeability. Effects on permeability and iron content calculations due to a possible lunar ionosphere are also considered.

Parkin, C. W.↗

Self‐gated, dynamic contrast‐enhanced magnetic resonance imaging with compressed‐sensing reconstruction for evaluating endothelial permeability in the aortic root of atherosclerotic mice

High‐risk atherosclerotic plaques are characterized by active inflammation and abundant leaky microvessels. We present a self‐gated, dynamic contrast‐enhanced magnetic resonance imaging (DCE‐MRI) acquisition with compressed sensing reconstruction and apply it to assess longitudinal changes in endothelial permeability in the aortic root of Apoe −/− atherosclerotic mice during natural disease progression. Twenty‐four, 8‐week‐old, female Apoe −/− mice were divided into four groups (n = 6 each) and imaged with self‐gated DCE‐MRI at 4, 8, 12, and 16 weeks after high‐fat diet initiation, and then euthanized for CD68 immunohistochemistry for macrophages. Eight additional mice were kept on a high‐fat diet and imaged longitudinally at the same time points. Aortic‐root pseudo‐concentration curves were analyzed using a validated piecewise linear model. Contrast agent wash‐in and washout slopes ( b 1 and b 2 ) were measured as surrogates of aortic root endothelial permeability and compared with macrophage density by immunohistochemistry. b 2 , indicating contrast agent washout, was significantly higher in mice kept on an high‐fat diet for longer periods of time ( p = 0.03). Group comparison revealed significant differences between mice on a high‐fat diet for 4 versus 16 weeks ( p = 0.03). Macrophage density also significantly increased with diet duration ( p = 0.009). Spearman correlation between b 2 from DCE‐MRI and macrophage density indicated a weak relationship between the two parameters (r = 0.28, p = 0.20). Validated piecewise linear modeling of the DCE‐MRI data showed that the aortic root contrast agent washout rate is significantly different during disease progression. Further development of this technique from a single‐slice to a 3D acquisition may enable better investigation of the relationship between in vivo imaging of endothelial permeability and atherosclerotic plaques' genetic, molecular, and cellular makeup in this important model of disease.

Calcagno, Claudia↗

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗