Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Breakthrough Curves”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Direct Breakthrough Curve Prediction From Statistics of Heterogeneous Conductivity Fields

In this paper, we present a methodology to predict the shape of solute breakthrough curves in heterogeneous aquifers at early times and/or under high degrees of heterogeneity, both cases in which the classical macrodispersion theory may not be applicable. The methodology relies on the observation that breakthrough curves in heterogeneous media are generally well described by lognormal distributions, and mean breakthrough times can be predicted analytically. The log-variance of solute arrival is thus sufficient to completely specify the breakthrough curves, and this is calibrated as a function of aquifer heterogeneity and dimensionless distance from a source plane by means of Monte Carlo analysis and statistical regression. Using the ensemble of simulated groundwater flow and solute transport realizations employed to calibrate the predictive regression, reliability estimates for the prediction are also developed. Additional theoretical contributions include heuristics for the time until an effective macrodispersion coefficient becomes applicable, and also an expression for its magnitude that applies in highly heterogeneous systems. Finally, it is seen that the results here represent a way to derive continuous time random walk transition distributions from physical considerations rather than from empirical field calibration.

58 GEOSCIENCES↗
DOE OSTI2019

Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves

Here, we develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ~ time -(1+α) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly (α<1), the late-time decay rate of the breakthrough curve is -(1+α/2) rather than the classical -3/2. Yet, if α>1, the -3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early-time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late-time decay rate. The theoretical predictions are validated against particle tracking simulations in the three-dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.

58 GEOSCIENCES↗
DOE OSTI2022

Machine-Learning of Nonlocal Kernels for Anomalous Subsurface Transport from Breakthrough Curves

Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Improved Method for Estimating Reaction Rates During Push‐Pull Tests

Abstract The breakthrough curve obtained from a single‐well push‐pull test can be adjusted to account for dilution of the injection fluid in the aquifer fluid. The dilution‐adjusted breakthrough curve can be analyzed to estimate the reaction rate of a solute. The conventional dilution‐adjusted method assumes that the ratios of the concentrations of the nonreactive and reactive solutes in the injection fluid vs. the aquifer fluid are equal. If this assumption is invalid, the conventional method will generate inaccurate breakthrough curves and may lead to erroneous conclusions regarding the reactivity of a solute. In this study, a new method that generates a dilution‐adjusted breakthrough curve was theoretically developed to account for any possible combination of nonreactive and reactive solute concentrations in the injection and aquifer fluids. The newly developed method was applied to a field‐based data set and was shown to generate more accurate dilution‐adjusted breakthrough curves. The improved dilution‐adjusted method presented here is simple, makes no assumptions regarding the concentrations of the nonreactive and reactive solutes in the injection and aquifer fluids, and easily allows for estimating reaction rates during push‐pull tests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Chromatographic Modeling of 1-Foot Detritiation Test Column Results for the Optimization of Regeneration Time

This report documents the methods and results of a modeling study on a series of data generated by 1-foot test columns for the Water Detritiation project. In order to compliment modeling work already in progress on this project, a purely empirical curve-fitting approach is used in this study to determine chromatographic performance indicators such as retention time, curve width, and the number of chromatographic stages. The data used in this study comes from four sets of chromatographic test separations performed in March through May of 2017 using 1-foot columns. All are using the same flow rates and temperature (room temperature, nominally 25 C°). The first set of data (referred to in this report as the Column 1 data) is a series of twelve breakthrough curve measurements. Prior to each of these twelve runs, the column was regenerated by flowing hydrogen gas through the column, for increasing amounts of time, ranging from 2 hours to 25 hours. The second set of data (referred to as the Column 2 data) is comprised of five measurements of breakthrough curves measured in similar fashion as the Column 1 data, across a range of regeneration times from 2 to 21 hours. The third set of data (referred to as the Purge data) is of three breakthrough curves, measured after a regeneration time of 8 hours, which takes place after a nitrogen purge of 15, 30, or 45 minutes. The fourth data set is a series of five breakthrough curve measurements using a 2 nd batch of stationary phase, taken after regeneration times ranging from 2 to 20 hours. The overall modeling approach is as follows. First, the data is pre-processed for formatting, artifact removal, and is normalized. The data is then fit by integrated versions of four common chromatographic models. The best model is chosen based on the predictive performance of the model, and chromatographic performance values are calculated based on the best model’s fit parameters. Each of these processes are detailed below.

42 ENGINEERING↗
DOE OSTI2013

Fluid bed adsorption of carbon dioxide on immobilized polyethyenimine (PEI): kinetic analysis and breakthrough behavior

The adsorption of carbon dioxide (CO{sub 2}) by immobilized polyethylenimine (PEI) on mesoporous silica was investigated in a fluid bed. The tests were performed to determine breakthrough behavior with varying bed temperature, flow rates and feed concentrations. Experimental breakthrough curves were analyzed using a theoretical 1D model developed by Bohart and Adams. The results showed that Bohart-Adams model was suitable for the normal description of breakthrough curve for the temperature ranges of 40-90{degree}C. The maximum capacity increased with temperature up to 70{degree}C and then decreased. The adsorption rate constant exhibited a negative temperature dependence decreasing as the temperature increased. Parameters characteristic of a fluid bed adsorber were inferred from these breakthrough curves including the breakthrough time, saturation time, critical reactor length, and length of mass transfer zone LMTZ. These parameters can be used to design fluid bed adsorption system without resolving the mechanistic contributions of dispersion, mixing, and intraparticle diffusion.

Breakthrough↗
DOE OSTI2023

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Implications of Grain-Scale Mineralogical Heterogeneity for Radionuclide Transport in Fractured Media

The geological disposal of nuclear waste is based on the multi-barrier concept, comprising various engineered and natural barriers, to confine the radioactive waste and isolate it from the biosphere. Some of the planned repositories for high-level nuclear waste will be hosted in fractured crystalline rock formations. The potential of these formations to act as natural transport barriers is related to two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces available in the rock matrix. Different in situ and laboratory experiments have pointed out the ubiquitous heterogeneous nature of the rock matrix: mineral surfaces and pore space are distributed in complex microstructures and their distribution is far from being homogeneous (as typically assumed by Darcy-scale coarse reactive transport models). We use a synthetically generated fracture–matrix system to assess the implications of grain-scale physical and mineralogical heterogeneity on cesium transport and retention. The resulting grain-scale reactive transport model is solved using high-performance computing technologies, and the results are compared with those derived from two alternative models, denoted as upscaled models, where mineral abundance is averaged over the matrix volume. In the grain-scale model, the penetration of cesium into the matrix is faster and the penetration front is uneven and finger-shaped. The analysis of the cesium breakthrough curves computed at two different points in the fracture shows that the upscaled models provide later first-arrival time estimates compared to the grain-scale model. The breakthrough curves computed with the three models converge at late times. These results suggest that spatially averaged upscaled parameters of sorption site distribution can be used to predict the late-time behavior of breakthrough curves but could be inadequate to simulate the early behavior.

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

INL Intern Poster Session

Characterization of off-gas constituents from nuclear fuel reprocessing can be measured using a residual gas analyzer. Characterization of this off-gas is important in a closed fuel cycle, where gasses such as Nox, Kr, I, and C-14 need to be removed from reprocessing plants to meet regulatory requirements. Both Kr and Xe have low chemical reactivity, making capture and separation difficult. Solid sorbents can be used to capture these gases as an alternative to cryogenic distillation, which is energy intensive and hazardous. AgZ-PAN and HZ-PAN are two solid sorbents that have demonstrated selective capture of Xe and Kr. A residual gas analyzer can be used to obtain breakthrough curves that demonstrate how well this adsorption and separation occur. The breakthrough curve is dependent on packing density, column size, and flow rate. Data obtained from the RGA was verified with a gas chromatographer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Pore-Scale Simulation of Intragranular Diffusion: Effects of Incomplete Mixing on Macroscopic Manifestations

Diffusive mass transfer into and out of intragranular micropores ("intragranular diffusion") plays an important role in the large-scale transport of some groundwater contaminants. We are interested in understanding the combined effect of pore-scale advection and intragranular diffusion on solute transport at the effective porous medium scale. To study this problem, we have developed a three-dimensional pore-scale model of fluid flow and solute transport that incorporates diffusion into and out of intragranular pore spaces. Our model is based on the Smoothed Particle Hydrodynamics (SPH) simulation method, which represents fluid and solid phases by a mesh-free particle discretization. In the pore spaces, fluid flow is simulated by discretizing the Navier-Stokes equations using the SPH approach. Solute transport is represented by advection, diffusion within the fluid phase, and diffusion between the fluid and solid phases. Our model is implemented on large-scale parallel computing hardware, allowing us to simulate millions of computational particles and represent fully three-dimensional systems of pores and grains with arbitrarily complex physical geometry. We have used this model system to perform numerical experiments using various model porous media systems, which allows us to draw comparisons between macroscopic measures computed from the pore-scale simulations (such as breakthrough curves) and those predicted by macroscopic formulations that assume complete mixing over the representative volume. In this paper we present results of 3D simulations of pore-scale flow and transport, including cases with and without intragranular diffusion, in two model porous media, one with randomly-packed uniform spherical grains and a second with randomly-packed spheres drawn from a binary grain size distribution. Breakthrough curves were computed from the 3D simulations at various transport distances. Comparable breakthrough curves were computed using 1D macroscopic models with parameters determined independently to the degree possible. Based on comparisons of the pore-scale and macroscopic model results, we draw two primary conclusions. First, non-Fickian behavior is persistent and ubiquitous at the scales considered, and most cases are better represented by a multi-rate mass transfer model even when there is no distinct secondary porosity (i.e., no intragranular diffusion). This suggests that diffusive mass transfer processes between preferential flow paths and relatively immobile zones within the primary porosity may have significant impact on transport, particular in low-concentration tails. Second, the application of mass transfer rate parameters based on an assumption of well-mixed concentrations at the pore-scale tend to overestimate the amount of mass transfer that occurs in heterogeneous pore geometries in which preferential flow leads to incomplete pore-scale lateral mixing.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Push-pull tracer tests: Their information content and use for characterizing non-Fickian, mobile-immobile behavior: INFORMATION CONTENT OF PUSH-PULL TESTS

Path reversibility and radial symmetry are often assumed in push-pull tracer test analysis. In reality, heterogeneous flow fields mean that both assumptions are idealizations. In this paper, to understand their impact, we perform a parametric study which quantifies the scattering effects of ambient flow, local-scale dispersion, and velocity field heterogeneity on push-pull breakthrough curves and compares them to the effects of mobile-immobile mass transfer (MIMT) processes including sorption and diffusion into secondary porosity. We identify specific circumstances in which MIMT overwhelmingly determines the breakthrough curve, which may then be considered uninformative about drift and local-scale dispersion. Assuming path reversibility, we develop a continuous-time-random-walk-based interpretation framework which is flow-field-agnostic and well suited to quantifying MIMT. Adopting this perspective, we show that the radial flow assumption is often harmless: to the extent that solute paths are reversible, the breakthrough curve is uninformative about velocity field heterogeneity. Our interpretation method determines a mapping function (i.e., subordinator) from travel time in the absence of MIMT to travel time in its presence. A mathematical theory allowing this function to be directly “plugged into” an existing Laplace-domain transport model to incorporate MIMT is presented and demonstrated. Algorithms implementing the calibration are presented and applied to interpretation of data from a push-pull test performed in a heterogeneous environment. A successful four-parameter fit is obtained, of comparable fidelity to one obtained using a million-node 3-D numerical model. In conclusion, we demonstrate analytically and numerically how push-pull tests quantifying MIMT are sensitive to remobilization, but not immobilization, kinetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Transport of perfluorocarbon tracers in the Cranfield Geological Carbon Sequestration Project

A field–scale carbon dioxide (CO 2 ) injection pilot project was conducted by the Southeast Regional Carbon Sequestration Partnership (SECARB) at Cranfield, Mississippi. Two associated campaigns in 2009 and 2010 were carried out to co–inject perfluorocarbon tracers (PFTs) and sulfur hexafluoride (SF 6 ) with CO 2 . Tracers in gas samples from two observation wells were analyzed to construct breakthrough curves. In this work, we present the field data and numerical modeling of the flow and transport of CO 2 , brine, and tracers. A high–resolution static model of the formation geology in the detailed area study (DAS) was used to capture the impact of connected flow pathways created by fluvial channels on breakthrough curves and breakthrough times of PFTs and SF6 tracers. We use the cubic–plus–association (CPA) equation of state, which takes into account the polar nature of water molecules, to describe the phase behavior of CO 2 –brine–tracer mixtures. Our simulated results show good agreement for the 2009 tracer campaign in Cranfield, while a larger discrepancy emerges by 2010. Here, the combination of multiple tracer injection pulses with detailed numerical simulations proves to be a powerful tool in constraining both formation properties and how complex flow paths develop over time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Field experiments of surface water to groundwater recharge to characterize the mobility of uranium and vanadium at a former mill tailing site

Characterizing the mobility of uranium and vanadium in groundwater with a hydraulic connection to surface water is important to inform the best management practices of former mill tailing sites. In this study, the recharge of river water to the unsaturated and saturated zones of a uranium-contaminated alluvial aquifer was simulated in a series of forced-gradient single- and multi-well injection-extraction tests. The injection fluid (river water) was traced with natural and artificial tracers that included halides, fluorobenzoates, lithium, and naphthalene sulfonate to characterize the potential mass transport mechanisms of uranium and vanadium. The extraction fluid (river water/groundwater mixture) was analyzed for the tracers, uranium, and vanadium. The results from the tracers indicated that matrix diffusion was likely negligible over the spatiotemporal scales of the tests as evident by nearly identical breakthrough curves of the halides and fluorobenzoates. In contrast, the breakthrough curves of lithium and naphthalene sulfonate indicated that sorption by cation exchange and sorption to organic matter, respectively, were potential mass transport mechanisms of uranium and vanadium. Uranium was mobilized in the saturated zone containing gypsum (gypsum-rich zone), the vadose zone (vadose-rich zone), and the saturated zone containing organic carbon (organic-rich zone) whereas vanadium was mobilized only in the saturated gypsum-rich zone. The mechanisms responsible for the mobilization of uranium and vanadium were likely dissolution of uranium- and vanadium-bearing minerals and/or desorption from the gypsum-rich zone, flushing of uranium from the vadose-rich zone, and desorption of uranium from the organic-rich zone due to the natural contrast in the geochemistry between the river water and groundwater. The experimental design of this study was unique in that it employed the use of multiple natural and artificial tracers coupled with a direct injection of native river water to groundwater. Here, these results demonstrated that natural recharge and flooding events at former mill tailing sites can mobilize uranium, and possibly vanadium, and contribute to persistent levels of groundwater contamination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Using the Diffusion Equation to Predict CST Resin Breakthroughs - 19051

The crystalline silicotitanate, CST, resin manufactured by UOP Honeywell is well-suited for single-use decontamination of caustic liquid wastes at Hanford and Savannah River Sites. There are many measurements of column loading and breakthrough for CST resin decontamination of Cs-137 in various Hanford and SRS waste simulants. Such measurements provide a very good basis for estimating CST resin demand and these more refined predictions make the CST column loading as efficient as possible given the very expensive CST resin. Although there are a variety of tools for predicting chromatographic resin column performance, [1-4], there is a particularly useful method that is suitable as a simple spreadsheet model [4]. This spreadsheet model is an analytical error function for the mass transfer zone (mtz) diffusion width in resin columns and is well-suited for resin breakthrough and column loading. This paper describes just such an error function spreadsheet model and how well the error function agrees with more traditional models and fits the observed small scale breakthrough curves for CST that can then better predict large scale column loading efficiency. The very high sorption of Cs onto CST tends to be slower than that of other lower sorption resins and as a result of the slow CST kinetics, there is also a marked asymmetry in many CST breakthrough curves. The CST breakthrough asymmetry is largely due to the distribution of CST particle sizes because smaller particles sorb much faster than larger particles. It is then the smaller particles that determine the leading mtz at a decontamination factor of 1,000 breakthrough and therefore the decontamination efficiency for a given loading. The larger CST particles correspondingly determine the tailing mtz width and therefore the overall loading efficiency and it is important to model both leading and tailing mtz widths to best optimize CST for a particular Cs-137 decontamination application. (authors)

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

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗