SEARCH · Search NASA

Results for “solute transport”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Predicting Water Cycle Characteristics from Percolation Theory and Observational Data

The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth’s surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET 0 /P, with P being the precipitation and ET 0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko’s model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2006

Isotopic Dependence of GCR Fluence behind Shielding

In this paper we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR), nuclear fragmentation cross-sections, and isotopic-grid on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. For the nuclear interaction data-base and transport solution, we use the quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, respectively. The QMSFRG model is shown to accurately describe existing fragmentation data including proper description of the odd-even effects as function of the iso-spin dependence on the projectile nucleus. The principle finding of this study is that large errors (+/-100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotopic-grid (approx.170 ions) to ones that use a reduced isotopic-grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (<+/-20%) occur in the elemental-fluence spectra. Because a complete isotopic-grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

DOE OSTI2020

Layer size polydispersity in hydrated montmorillonite creates multiscale porosity networks

The aluminosilicate layers of the swelling clay mineral montmorillonite, and the saturated pores they delineate, control the mechanical properties and the transport of solutes in many natural and engineered environments. However, the structural basis of montmorillonite porosity remains poorly characterized due to the difficulty in visualizing hydrated samples in their native state. Here in this paper, we used cryogenic transmission electron microscopy (cryo-TEM) and cryo electron tomography (cryo-ET) to show that stacking defects within minimally altered, fully hydrated montmorillonite particles define multiscale porosity networks. Variations in layer lateral dimensions over tens to thousands of nanometers cause a range of topological and dynamic defects that generate pervasive curvature and introduce previously uncharacterized solute transport pathways. Observations of long-range rotational order between neighboring layers indicate that the layer-layer interactions that govern clay swelling involve three dimensional orienting forces that operate across nanoscale pores. These direct observations of the hierarchical structure of hydrated montmorillonite pore networks with nanoscale resolution reveal potentially general aspects of colloidal interactions in fluid-saturated clay minerals.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Idaho National Laboratory’s Mobile Hot Cell Transportation: Engineering Solutions for Global Disused Sealed Radioactive Sources.

Title: Idaho National Laboratory’s Mobile Hot Cell Transportation: Engineering Solutions for Global Disused Sealed Radioactive Sources. Abstract: The Mobile Hot Cell (MHC), currently under development by Idaho National Laboratory (INL) for the Off-Site Source Recovery Project (OSRP), is designed to help international partners meet the unique challenges of end-of-life radioactive material management. The MHC will provide a critical resource for countries that require assistance securing and disposing of Disused Sealed Radioactive Sources (DSRS) and orphaned sources in challenging environments, allowing these sources to be secured against misuse and nefarious activities. The MHC is a rapidly deployable system for conditioning and preparing end-of-life radioactive sources for transportation or storage. It is designed to handle sources of up to 1,000 Ci Co-60 equivalent while maintaining full radiological and biological containment. It will be deployable within 48 hours of an alert, making it ideal for emergency situations. The MHC features an operational suite for control, support racks for electronics, pneumatics, and welding systems, and a modular robust steel structure providing radiological shielding and internal robotic support. This design allows configurations for multiple device types to be conditioned and the ability to safely manage routine issues such as leaking or damaged sources. The MHC has been designed with the transportation challenges of rapid deployment to difficult environments in mind. The system weighs approximately 150,000 pounds, with individual systems breaking down into pieces not exceeding 20,000 pounds. Components are to be transportable on standard 20ft ISO containers, with shielding shells on 20ft flat racks. It is estimated that a total of eight containers and flat racks will be required. The use of 20ft containers, as opposed to 40ft containers, minimizes the impact on less developed road infrastructures, enabling the MHC to be positioned in constrained environments such as hospital parking lots. The system’s modularity also allows for deployment using smaller equipment, such as a 10-ton boom truck or forklift, which is crucial given the potential logistical challenges in different countries. This transportation strategy, evaluated in collaboration with Utah State University, ensures the MHC can be deployed via ground, rail, sea, or air, addressing the primary concern of international transport logistics.

99 - GENERAL AND MISCELLANEOUS↗

NASA NTRS1986

Rationalizing VICAR within a TAE framework: Some problems and some solutions

Transportable Applications Executive (TAE) implementation may impose a strain on centers with modest resources. This may be eased in a number of ways. The balance of a small number of expert users and a large number of computing novices at the Interactive Planetary Image-Processing System (IPIPS) imposes special constraints. Some solutions to these and other particular problems are described.

Harris, Colin P.↗

DOE OSTI2020

Stress effects on flow and transport in three-dimensional fracture networks

We investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three-dimensional fracture networks. We use well-established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low-velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

58 GEOSCIENCES↗

DOE OSTI2024

Mechanistic definition and prediction of the mass exchange coefficient between rivers and hyporheic zones: The $α$ of two $Ω$s

Solute transport in interconnected rivers and hyporheic zones is typically modeled through dual-domain models where first-order solute mass transfer between the two domains, Ω R and Ω HZ , is represented by a coefficient α. The transient storage model (TSM) is an example of such an approach. In practice, α is determined by fitting the tails of solute tracer breakthrough curves using a TSM. This approach has led to ambiguity regarding α’s physical meaning and transferability. Here, in this work, we investigated the physical basis for α and tested it with virtual experiments through the fully coupled multiphysics model hyporheicFoam for the Ω R – Ω HZ system. hyporheicFoam explicitly simulated coupled flow and solute transport over a kilometer with centimeter-scale resolution. Model results were analyzed to calculate α following its theoretical definition directly. Using the determined α within a TSM enables accurate reproduction of solute transport, underscoring α’s physical relevance and precision.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2003

Effect of Gravity Level on the Particle Shape and Size During Zeolite Crystal Growth

A microscopic diffusion model is developed to represent solute transport in the boundary layer of a growing zeolite crystal. This model is used to describe the effect of gravity on particle shape and solute distribution. Particle dynamics and crystal growth kinetics serve as the boundary conditions of flow and convection-diffusion equations. A statistical rate theory is used to obtain the rate of solute transport across the growing interface, which is expressed in terms of concentration and velocity of solute species. Microgravity can significantly decrease the solute velocity across the growing interface compared to its earth-based counterpart. The extent of this reduction highly depends on solute diffusion constant in solution. Under gravity, the flow towards the crystal enhances solute transport rate across the growing interface while the flow away from crystals reduces this rate, suggesting a non-uniform growth rate and thus an elliptic final shape. However, microgravity can significantly reduce the influence of flow and obtain a final product with perfect spherical shape. The model predictions compare favorably with the data of space experiment of zeolites grown in space.

Song, Hong-Wei↗

DOE OSTI2020

Unifying Advective and Diffusive Descriptions of Bedform Pumping in the Benthic Biolayer of Streams

Abstract Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment‐water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow‐dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.

Grant, Stanley B.↗

DOE OSTI2020

The Complexity of Nonlinear Flow and non-Fickian Transport in Fractures Driven by Three-Dimensional Recirculation Zones

Subsurface fluid flow and solute transport are typically described by first-order or linear rate laws. However, deviations from the first-order rate laws (or anomalous behavior) are typical and lead to nonlinear flow and non-Fickian transport phenomena. The shared underlying mechanisms for these “anomalous” flow and transport phenomena, which are actually typical of fractures, are seldom analyzed. We investigated the co-occurrence of nonlinear flow and non-Fickian transport through numerical simulation experiments of pore-scale flow and solute transport processes through three-dimensional natural fractures. In these fractures, 3-D recirculation zones (RZs) developed and grew with increasing flow rates. Flow nonlinearity resulted from the RZ development and growth, which shrank the main flow channel. Moreover, the same RZs also captured solutes from the main flow channel and released them back later. This retardation resulted in non-Fickian transport. Thus, nonlinear flow and non-Fickian transport are intertwined via their shared dependence on RZs.

58 GEOSCIENCES↗

NASA NTRS1987

Dendritic growth in a supercooled alloy melt

A simple model which describes the growth of an 'array' of dendrites into a supercooled, binary, alloy melt is presented. Solute diffusion is calculated by superposing the solutions given by Flemings and Zener, and also, by superposing the solutions given by Ivantsov and Flemings. A general expression for the transport solution is suggested from which all other dendrite growth models presented earlier may be obtained as special cases. It is shown that both 'free' and 'constrained' growth may be described by a single transport solution, which indicates that (1) both thermal and solutal effects will be important during 'free' growth in dilute alloys, (2) only solutal effects are predominant during 'free' growth in concentrated alloys and during 'constrained' growth. An examination of the relevant dimensionless parameters also suggests that all dendrite growth models, regardless of the assumptions used to determine the tip radius (marginal stability, minimum undercooling, maximum velocity, minimum entropy production) should predict the experimentally observed extrema in tip radius and growth velocity in dilute alloys, during 'free' dendritic growth. Experimental data in binary H2O-NaCl and succinonitrile-acetone solutions are shown to be in good agreement with the model.

Laxmanan, V.↗

DOE OSTI2022

Quantifying subsurface parameter and transport uncertainty using surrogate modelling and environmental tracers

Here, we combine physics-based groundwater reactive transport modelling with machine-learning techniques to quantify hydrogeological model and solute transport predictive uncertainties. We train an artificial neural network (ANN) on a dataset of groundwater hydraulic heads and 3 H concentrations generated using a high-fidelity groundwater reactive transport model. Using the trained ANN as a surrogate model to reproduce the input–output response of the high-fidelity reactive transport model, we quantify the posterior distributions of hydrogeological parameters and hydraulic forcing conditions using Markov chain Monte Carlo calibration against field observations of groundwater hydraulic heads and 3 H concentrations. We demonstrate the methodology with a model application that predicts Chlorofluorocarbon-12 (CFC-12) solute transport at a contaminated field site in Wyoming, United States. Our results show that including 3 H observations in the calibration dataset reduced the uncertainty in the estimated permeability field and infiltration rates, compared to calibration against hydraulic heads alone. However, predictive uncertainty quantification shows that CFC-12 transport predictions conditioned to the parameter posterior distributions cannot reproduce the field measurements. We found that calibrating the model to hydraulic head and 3 H observations results in groundwater mean ages that are too large to explain the observed CFC-12 concentrations. The coupling of the physics-based reactive transport model with the machine-learning surrogate model allows us to efficiently quantify model parameter and predictive uncertainties, which is typically computationally intractable using reactive transport models alone.

58 GEOSCIENCES↗

DOE OSTI2020

PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data

Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

Effects of Convective Transport of Solute and Impurities on Defect-Causing Kinetics Instabilities in Protein Crystallization

The objective of the proposed research is to obtain further insight into the onset and development of the defect-causing instabilities that anise due to the coupling of the bulk transport and nonlinear-interfacial kinetics during growth in the mixed regime, utilizing the reduction of the convective contribution to the bulk transport under microgravity. These studies will build upon the data on the effects of quantitative variations of the forced convection velocity on the averaged and time-dependent kinetic behavior of protein crystal growth systems that have recently been obtained in our laboratory.

Vekilov, Peter G.↗

DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

Fluid Flow and Solidification Under Combined Action of Magnetic Fields and Microgravity

Mathematical models, both 2-D and 3-D, are developed to represent g-jitter induced fluid flows and their effects on solidification under combined action of magnetic fields and microgravity. The numerical model development is based on the finite element solution of governing equations describing the transient g-jitter driven fluid flows, heat transfer and solutal transport during crystal growth with and without an applied magnetic field in space vehicles. To validate the model predictions, a ground-based g-jitter simulator is developed using the oscillating wall temperatures where timely oscillating fluid flows are measured using a laser PIV system. The measurements are compared well with numerical results obtained from the numerical models. Results show that a combined action derived from magnetic damping and microgravity can be an effective means to control the melt flow and solutal transport in space single crystal growth systems.

Li, B. Q.↗

NASA NTRS2001

Magnetic Damping of g-Jitter Induced Double-Diffusive Convection

This paper describes a numerical study of the g-jitter driven double diffusive convective flows, thermal and concentration distributions in binary alloy melt systems subject to an external magnetic field. The study is based on the finite element solution of transient magnetohydrodynamic equations governing the momentum, thermal and solutal transport in the melt pool. Numerical simulations are conducted using the synthesized single- and multi- frequency g-jitter as well as the real g-jitter data taken during space flights with or without an applied magnetic field. It is found that for the conditions studied, the main melt flow follows approximately a lineal- superposition of velocity components induced by individual g-jitter components, regardless of whether a magnetic field exists or not. The flow field is characterized by a recirculating double diffusive convection loop oscillating in time with a defined frequency equal to that of the driving g-jitter force. An applied magnetic field has little effect on the oscillating recirculating pattern, except around the moment in time when the flow reverses its direction. The field has no effect on the oscillation period, but it changes the phase angle. It is very effective in suppressing the flow intensity and produces a notable reduction of the solutal striation and time fluctuations in the melt. For a given magnetic field strength, the magnetic damping effect is more pronounced on the velocity associated with the largest g-jitter component present and/or the g-jitter spiking peaks. A stronger magnetic field is more effective in suppressing the melt convection and also is more helpful in bringing the convection in phase with the g-jitter driving force. The applied field is particularly useful in suppressing the effect of real g-jitter spikes on both flow and solutal distributions. With appropriately selected magnetic fields, the convective flows caused by g-jitter can be reduced sufficiently and diffusion dominant. solutal transport in the melt is possible.

Shu, Y.↗