Lewis, F.M.; Voss, C.I.; Rubin, Jacob
1986-01-01
A model was developed that can simulate the effect of certain chemical and sorption reactions simultaneously among solutes involved in advective-dispersive transport through porous media. The model is based on a methodology that utilizes physical-chemical relationships in the development of the basic solute mass-balance equations; however, the form of these equations allows their solution to be obtained by methods that do not depend on the chemical processes. The chemical environment is governed by the condition of local chemical equilibrium, and may be defined either by the linear sorption of a single species and two soluble complexation reactions which also involve that species, or binary ion exchange and one complexation reaction involving a common ion. Partial differential equations that describe solute mass balance entirely in the liquid phase are developed for each tenad (a chemical entity whose total mass is independent of the reaction process) in terms of their total dissolved concentration. These equations are solved numerically in two dimensions through the modification of an existing groundwater flow/transport computer code. (Author 's abstract)
Modeling solute advection coupled with sorption kinetics in heterogeneous formations
NASA Astrophysics Data System (ADS)
Selroos, Jan-Olof; Cvetkovic, Vladimir
1992-05-01
A method for coupling sorption kinetics and solute advection in particle-tracking models is proposed; this method is efficient for the case where sorption rate coefficients can be assumed constant field scale parameters. A simulation example of reactive solute advection in two-dimensional heterogeneous porous media is presented. The effect of sorption kinetics on solute advection is investigated. Nonequilibrium effects are exhibited as enhanced tailing in the solute breakthrough. Because high variability in the hydraulic conductivity also yields enhanced tailing, the nonequilibrium effect is more pronounced for the case of low variability. Moreover, it may be difficult to distinguish cases of low variability with nonequilibrium sorption from cases of high variability with equilibrium sorption. A comparison of Monte Carlo ensemble results is made with an analytical model for the mass arrival of kinetically sorbing solute in heterogeneous porous media obtained using first-order perturbation. The comparison indicates that the analytical model provides reasonable approximations of the expected solute breakthrough if the variance of the natural logarithm of the hydraulic conductivity is smaller than 1.
Backward fractional advection dispersion model for contaminant source prediction
NASA Astrophysics Data System (ADS)
Zhang, Yong; Meerschaert, Mark M.; Neupauer, Roseanna M.
2016-04-01
The forward Fractional Advection Dispersion Equation (FADE) provides a useful model for non-Fickian transport in heterogeneous porous media. The space FADE captures the long leading tail, skewness, and fast spreading typically seen in concentration profiles from field data. This paper develops the corresponding backward FADE model, to identify source location and release time. The backward method is developed from the theory of inverse problems, and then explained from a stochastic point of view. The resultant backward FADE differs significantly from the traditional backward Advection Dispersion Equation (ADE) because the fractional derivative is not self-adjoint and the probability density function for backward locations is highly skewed. Finally, the method is validated using tracer data from a well-known field experiment, where the peak of the backward FADE curve predicts source release time, while the median or a range of percentiles can be used to determine the most likely source location for the observed plume. The backward ADE cannot reliably identify the source in this application, since the forward ADE does not provide an adequate fit to the concentration data.
Purely Lagrangian Simulation of Advection, Dispersion, Precipitation, and Dissolution
NASA Astrophysics Data System (ADS)
Benson, D.; Zhang, Y.; Reeves, D. M.
2008-05-01
We extend the advantages of Lagrangian random walk particle tracking (RWPT) methods that have long been used to simulate advection and dispersion in highly heterogeneous media. By formulating dissolution as a random, independent decay process, the classical continuum rate law is recovered. Formulating the random precipitation process requires a consideration of the probability that two nearby particles will coincide in a given time period. This depends on local mixing (as by diffusion) and the total domain particle number density, which are fixed and therefore easy to calculate. The result is that the classical law of mass action for equilibrium reactions can be reproduced in an ensemble sense. The same number of parameters for A+B ⇌ C are needed in a probabilistic versus continuum reaction simulation-- —one each for forward and backward probabilities that correspond to rates. The random nature of the simulations allows for significant disequilibrium in any given region at any time that is independent of the numerical details such as time stepping or particle density. This is exemplified by nearby or intermingled groups of reactants and little or no product--—a result that is often noted in the field that is difficult to reconcile with continuum methods or coarse-grained Eulerian models. Our results support recent results of perturbed advection-dispersion-reaction continuum models (Luo et al., WRR 44, 2008), and suggest that many different kinds of reactions can be easily added to existing RWPT codes.
Fractional Advective-Dispersive Equation as a Model of Solute Transport in Porous Media
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. The common model is the advective-dispersive equation (ADE) describing the superposition of the advective transport and the Brownian motion in water-filled pore space. Deviations from...
Analytical solution for the advection-dispersion transport equation in layered media
Technology Transfer Automated Retrieval System (TEKTRAN)
The advection-dispersion transport equation with first-order decay was solved analytically for multi-layered media using the classic integral transform technique (CITT). The solution procedure used an associated non-self-adjoint advection-diffusion eigenvalue problem that had the same form and coef...
An exact peak capturing and essentially oscillation-free (EPCOF) algorithm, consisting of advection-dispersion decoupling, backward method of characteristics, forward node tracking, and adaptive local grid refinement, is developed to solve transport equations. This algorithm repr...
NASA Astrophysics Data System (ADS)
Maghrebi, M.; Jankovic, I.; Rabideau, A. J.; Allen-King, R. M.; Weissmann, G. S.
2011-12-01
Effects of three key transport mechanisms (advection, diffusion and sorption) on transport and contaminant tailing of chlorinated solvents have been investigated using a numerical model. Thousands of model simulations have been conducted for various combinations of transport parameters that govern three key mechanisms in order to quantify tailing and relative importance of each mechanism. Hydraulic conductivity model contains a single inclusion of constant conductivity K embedded in a homogeneous anisotropic background of conductivity Kh,Kv. The inclusion is shaped as an oblate ellipsoid and subject to uniform flow. The background represents "average" conductivity of a heterogeneous formation while inclusion is used to represent geologic units that are more or less conductive than the background. The ratio of long to short semi-axis of the inclusion (a/b) models the ratio of horizontal to vertical integral scales (Ih/Iv) of different geologic units, where integral scales can be obtained, for example, using indicator variograms. The flow solution for present problem is obtained analytically as a closed form solution with exact expressions for Darcy velocity valid both inside and outside the inclusion. Sorption is modeled as an equilibrium process governed by a linear isotherm. The effects on transport and tailing are accounted for using retardation factors. Sorption heterogeneity is created by allowing different values of retardation factor for the interior (Ri) and the exterior of the inclusion (Rb). Diffusive displacements have been added to retarded advective displacements using random walk method. Peclet number, defined as Pe=U Ih/D (U is the groundwater velocity, D is the molecular diffusion coefficient for chlorinated solvents), is used to quantify the diffusion process. Very large numbers of particles (hundreds of thousands) have been tracked using very small time steps (as small as a/10,000) to provide sufficient resolution to breakthrough curves and to
Garges, J.A.; Baehr, A.L.
1998-01-01
The relative importance of advection and dispersion for both solute and vapor transport can be determined from type curves or concentration, flux, or cumulative flux. The dimensionless form of the type curves provides a means to directly evaluate the importance of mass transport by advection relative to that of mass transport by diffusion and dispersion. Type curves based on an analytical solution to the advection-dispersion equation are plotted in terms of dimensionless time and Peclet number. Flux and cumulative flux type curves provide additional rationale for transport regime determination in addition to the traditional concentration type curves. The extension of type curves to include vapor transport with phase partitioning in the unsaturated zone is a new development. Type curves for negative Peclet numbers also are presented. A negative Peclet number characterizes a problem in which one direction of flow is toward the contamination source, and thereby diffusion and advection can act in opposite directions. Examples are the diffusion of solutes away from the downgradient edge of a pump-and-treat capture zone, the upward diffusion of vapors through the unsaturated zone with recharge, and the diffusion of solutes through a low hydraulic conductivity cutoff wall with an inward advective gradient.
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion equation and related models are indispensable for predicting or analyzing contaminant transport processes in streams and rivers, as well as in other surface water bodies. Many useful analytical solutions originated in disciplines other than surface-w...
A Lagrangian-Eulerian finite element method with adaptive gridding for advection-dispersion problems
Ijiri, Y.; Karasaki, K.
1994-02-01
In the present paper, a Lagrangian-Eulerian finite element method with adaptive gridding for solving advection-dispersion equations is described. The code creates new grid points in the vicinity of sharp fronts at every time step in order to reduce numerical dispersion. The code yields quite accurate solutions for a wide range of mesh Peclet numbers and for mesh Courant numbers well in excess of 1.
Technology Transfer Automated Retrieval System (TEKTRAN)
Analytical solutions of the advection-dispersion solute transport equation remain useful for a large number of applications in science and engineering. In this paper we extend the Duhamel theorem, originally established for diffusion type problems, to the case of advective-dispersive transport subj...
NASA Astrophysics Data System (ADS)
Sofiev, M.; Vira, J.; Kouznetsov, R.; Prank, M.; Soares, J.; Genikhovich, E.
2015-11-01
The paper presents the transport module of the System for Integrated modeLling of Atmospheric coMposition SILAM v.5 based on the advection algorithm of Michael Galperin. This advection routine, so far weakly presented in the international literature, is positively defined, stable at any Courant number, and efficient computationally. We present the rigorous description of its original version, along with several updates that improve its monotonicity and shape preservation, allowing for applications to long-living species in conditions of complex atmospheric flows. The scheme is connected with other parts of the model in a way that preserves the sub-grid mass distribution information that is a cornerstone of the advection algorithm. The other parts include the previously developed vertical diffusion algorithm combined with dry deposition, a meteorological pre-processor, and chemical transformation modules. The quality of the advection routine is evaluated using a large set of tests. The original approach has been previously compared with several classic algorithms widely used in operational dispersion models. The basic tests were repeated for the updated scheme and extended with real-wind simulations and demanding global 2-D tests recently suggested in the literature, which allowed one to position the scheme with regard to sophisticated state-of-the-art approaches. The advection scheme performance was fully comparable with other algorithms, with a modest computational cost. This work was the last project of Dr. Sci. Michael Galperin, who passed away on 18 March 2008.
İbiş, Birol
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
Ibiş, Birol; Bayram, Mustafa
2014-01-01
This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs. PMID:24578662
NASA Astrophysics Data System (ADS)
Bernabé, Y.; Wang, Y.; Qi, T.; Li, M.
2016-02-01
The main purpose of this work is to investigate the relationship between passive advection-dispersion and permeability in porous materials presumed to be statistically homogeneous at scales larger than the pore scale but smaller than the reservoir scale. We simulated fluid flow through pipe network realizations with different pipe radius distributions and different levels of connectivity. The flow simulations used periodic boundary conditions, allowing monitoring of the advective motion of solute particles in a large periodic array of identical network realizations. In order to simulate dispersion, we assumed that the solute particles obeyed Taylor dispersion in individual pipes. When a particle entered a pipe, a residence time consistent with local Taylor dispersion was randomly assigned to it. When exiting the pipe, the particle randomly proceeded into one of the pipes connected to the original one according to probabilities proportional to the outgoing volumetric flow in each pipe. For each simulation we tracked the motion of at least 6000 solute particles. The mean fluid velocity was 10-3 ms-1, and the distance traveled was on the order of 10 m. Macroscopic dispersion was quantified using the method of moments. Despite differences arising from using different types of lattices (simple cubic, body-centered cubic, and face-centered cubic), a number of general observations were made. Longitudinal dispersion was at least 1 order of magnitude greater than transverse dispersion, and both strongly increased with decreasing pore connectivity and/or pore size variability. In conditions of variable hydraulic radius and fixed pore connectivity and pore size variability, the simulated dispersivities increased as power laws of the hydraulic radius and, consequently, of permeability, in agreement with previously published experimental results. Based on these observations, we were able to resolve some of the complexity of the relationship between dispersivity and permeability.
It is well known that the fate and transport of contaminants in the subsurface are controlled by complex processes including advection, dispersion-diffusion, and chemical reactions. However, the interplay between the physical transport processes and chemical reactions, and their...
Zhao, Xiao; Liu, Wen; Fu, Jie; Cai, Zhengqing; O'Reilly, S E; Zhao, Dongye
2016-08-15
This work examined effects of model oil dispersants on dispersion, sorption and photodegradation of petroleum hydrocarbons in simulated marine systems. Three dispersants (Corexit 9500A, Corexit 9527A and SPC 1000) were used to prepare dispersed water accommodated oil (DWAO). While higher doses of dispersants dispersed more n-alkanes and PAHs, Corexit 9500A preferentially dispersed C11-C20 n-alkanes, whereas Corexit 9527A was more favorable for smaller alkanes (C10-C16), and SPC 1000 for C12-C28 n-alkanes. Sorption of petroleum hydrocarbons on sediment was proportional to TPH types/fractions in the DWAOs. Addition of 18mg/L of Corexit 9500A increased sediment uptake of 2-3 ring PAHs, while higher dispersant doses reduced the uptake, due to micelle-enhanced solubilization effects. Both dispersed n-alkanes and PAHs were susceptible to photodegradation under simulated sunlight. For PAHs, both photodegradation and photo-facilitated alkylation were concurrently taking place. The information can facilitate sounder assessment of fate and distribution of dispersed oil hydrocarbons in marine systems. PMID:27318763
NASA Astrophysics Data System (ADS)
Cartwright, Ian
Advection-dispersion fluid flow models implicitly assume that the infiltrating fluid flows through an already fluid-saturated medium. However, whether rocks contain a fluid depends on their reaction history, and whether any initial fluid escapes. The behaviour of different rocks may be illustrated using hypothetical marble compositions. Marbles with diverse chemistries (e.g. calcite + dolomite + quartz) are relatively reactive, and will generally produce a fluid during heating. By contrast, marbles with more restricted chemistries (e.g. calcite + quartz or calcite-only) may not. If the rock is not fluid bearing when fluid infiltration commences, mineralogical reactions may produce a reaction-enhanced permeability in calcite + dolomite + quartz or calcite + quartz, but not in calcite-only marbles. The permeability production controls the pattern of mineralogical, isotopic, and geochemical resetting during fluid flow. Tracers retarded behind the mineralogical fronts will probably be reset as predicted by the advection-dispersion models; however, tracers that are expected to be reset ahead of the mineralogical fronts cannot progress beyond the permeability generating reaction. In the case of very unreactive lithologies (e.g. pure calcite marbles, cherts, and quartzites), the first reaction to affect the rocks may be a metasomatic one ahead of which there is little pervasive resetting of any tracer. Centimetre-scale layering may lead to the formation of self-perpetuating fluid channels in rocks that are not fluid saturated due to the juxtaposition of reactants. Such layered rocks may show patterns of mineralogical resetting that are not predicted by advection-dispersion models. Patterns of mineralogical and isotopic resetting in marbles from a number of terrains, for example: Chillagoe, Marulan South, Reynolds Range (Australia); Adirondack Mountains, Old Woman Mountains, Notch Peak (USA); and Stephen Cross Quarry (Canada) vary as predicted by these models.
Dynamic typology of hydrothermal systems: competing effects of advection, dispersion and reactivity
NASA Astrophysics Data System (ADS)
Dolejs, David
2016-04-01
Genetic interpretation hydrothermal systems relies on recognition of (i) hydrothermal fluid source, (ii) fluid migration pathways, and (iii) deposition site identified by hydrothermal alteration and/or mineralization. Frequently, only the last object is of interest or accessible to direct observation, but constraints on the fluid source (volume) and pathways can be obtained from evaluation of the time-integrated fluid flux during hydrothermal event. Successful interpretation of the petrological record, that is, progress of alteration reactions, relies on identification of individual contributions arising from solute advection (to the deposition site), its lateral dispersion, and reaction efficiency. Although these terms are all applicable in a mass-conservation relationship within the framework of the transport theory, they are rarely considered simultaneously and their relative magnitudes evaluated. These phenomena operate on variable length and time scales, and may in turn provide insight into the system dynamics such as flow, diffusion and reaction rates, or continuous vs. episodic behavior of hydrothermal events. In addition, here we demonstrate that they also affect estimate of the net fluid flux, frequently by several orders of magnitude. The extent of alteration and mineralization reactions between the hydrothermal fluid and the host environment is determined by: (i) temperature, pressure or any other gradients across the mineralization site, (ii) magnitude of disequilibrium at inflow to the mineralization site, which is related to physico-chemical gradient between the fluid source and the mineralization site, and (iii) chemical redistribution (dispersion) within the mineralization site. We introduce quantitative mass-transport descriptors - Péclet and Damköhler II numbers - to introduce division into dispersion-dominated, advection-dominated and reaction-constrained systems. Dispersive systems are characterized by lateral solute redistribution, driven by
Huang, Y.H.; Saiers, J.E.; Harvey, J.W.; Noe, G.B.; Mylon, S.
2008-01-01
The movement of particulate matter within wetland surface waters affects nutrient cycling, contaminant mobility, and the evolution of the wetland landscape. Despite the importance of particle transport in influencing wetland form and function, there are few data sets that illuminate, in a quantitative way, the transport behavior of particulate matter within surface waters containing emergent vegetation. We report observations from experiments on the transport of 1 ??m latex microspheres at a wetland field site located in Water Conservation Area 3A of the Florida Everglades. The experiments involved line source injections of particles inside two 4.8-m-long surface water flumes constructed within a transition zone between an Eleocharis slough and Cladium jamaicense ridge and within a Cladium jamaicense ridge. We compared the measurements of particle transport to calculations of two-dimensional advection-dispersion model that accounted for a linear increase in water velocities with elevation above the ground surface. The results of this analysis revealed that particle spreading by longitudinal and vertical dispersion was substantially greater in the ridge than within the transition zone and that particle capture by aquatic vegetation lowered surface water particle concentrations and, at least for the timescale of our experiments, could be represented as an irreversible, first-order kinetics process. We found generally good agreement between our field-based estimates of particle dispersion and water velocity and estimates determined from published theory, suggesting that the advective-dispersive transport of particulate matter within complex wetland environments can be approximated on the basis of measurable properties of the flow and aquatic vegetation. Copyright 2008 by the American Geophysical Union.
Technology Transfer Automated Retrieval System (TEKTRAN)
The classical model to describe solute transport in soil is based on the advective-dispersive equation where Fick’s law is used to explain dispersion. From the microscopic point of view this is equivalent to consider that the motion of the particles of solute may be simulated by the Brownian motion....
Probabilistic exposure risk assessment with advective-dispersive well vulnerability criteria
NASA Astrophysics Data System (ADS)
Enzenhoefer, Rainer; Nowak, Wolfgang; Helmig, Rainer
2012-02-01
Time-related advection-based well-head protection zones are commonly used to manage the contamination risk of drinking water wells. According to current water safety plans advanced risk management schemes are needed to better control and monitor all possible hazards within catchments. The goal of this work is to cast the four advective-dispersive intrinsic well vulnerability criteria by Frind et al. [1] into a framework of probabilistic risk assessment framework. These criteria are: (i) arrival time, (ii) level of peak concentration, (iii) time until first arrival of critical concentrations and (iv) exposure time. Our probabilistic framework yields catchment-wide maps of probabilities to not comply with these criteria. This provides indispensable information for catchment managers to perform probabilistic exposure risk assessment and thus improves the basis for risk-informed well-head management. We resolve heterogeneity with high-resolution Monte Carlo simulations and use a new reverse formulation of temporal moment transport equations to keep computational costs low. Our method is independent of dimensionality and boundary conditions, and can account for arbitrary sources of uncertainty. It can be coupled with any method for conditioning on available data. For simplicity, we demonstrate the concept on a 2D example that includes conditioning on synthetic data.
Solution of the advection-dispersion equation: Continuous load of finite duration
Runkel, R.L.
1996-01-01
Field studies of solute fate and transport in streams and rivers often involve an. experimental release of solutes at an upstream boundary for a finite period of time. A review of several standard references on surface-water-quality modeling indicates that the analytical solution to the constant-parameter advection-dispersion equation for this type of boundary condition has been generally overlooked. Here an exact analytical solution that considers a continuous load of unite duration is compared to an approximate analytical solution presented elsewhere. Results indicate that the exact analytical solution should be used for verification of numerical solutions and other solute-transport problems wherein a high level of accuracy is required. ?? ASCE.
Bad behavior of Godunov mixed methods for strongly anisotropic advection-dispersion equations
NASA Astrophysics Data System (ADS)
Mazzia, Annamaria; Manzini, Gianmarco; Putti, Mario
2011-09-01
We study the performance of Godunov mixed methods, which combine a mixed-hybrid finite element solver and a Godunov-like shock-capturing solver, for the numerical treatment of the advection-dispersion equation with strong anisotropic tensor coefficients. It turns out that a mesh locking phenomenon may cause ill-conditioning and reduce the accuracy of the numerical approximation especially on coarse meshes. This problem may be partially alleviated by substituting the mixed-hybrid finite element solver used in the discretization of the dispersive (diffusive) term with a linear Galerkin finite element solver, which does not display such a strong ill conditioning. To illustrate the different mechanisms that come into play, we investigate the spectral properties of such numerical discretizations when applied to a strongly anisotropic diffusive term on a small regular mesh. A thorough comparison of the stiffness matrix eigenvalues reveals that the accuracy loss of the Godunov mixed method is a structural feature of the mixed-hybrid method. In fact, the varied response of the two methods is due to the different way the smallest and largest eigenvalues of the dispersion (diffusion) tensor influence the diagonal and off-diagonal terms of the final stiffness matrix. One and two dimensional test cases support our findings.
Horizontal advection and dispersion in a stratified shelf sea: The role of inertial oscillations
NASA Astrophysics Data System (ADS)
Inall, Mark E.; Aleynik, Dmitry; Neil, Clare
2013-10-01
The role played by inertial motions in horizontal dispersion within the thermocline of a broad, mid-latitude shelf sea is examined through the analysis of a deliberately released dye tracer. Our analysis is of the horizontal and vertical evolution over 40 h of a dye tracer injected into the seasonally stratified thermocline of the Celtic Sea on the NW European Shelf. The inferred diapycnal diffusivity was 1.3-1.5 × 10-5 m2 s-1, and the radial horizontal diffusivities of the depth integrated dye patch ranged from 1.9 to 4.0 m2 s-1. The inferred vertical diffusivity is in agreement with microstructure based estimates, and the depth integrated horizontal diffusivity is broadly in agreement with previous dye release derived estimates made over similar scales and time periods. Asymmetry in the horizontal evolution of the dye patch was evident. We argue that mean shear dispersion was responsible for lateral elongation of the dye patch, particularly between hours 23 and 35 after release, during which time horizontal diffusivity along the major axis, Ka, exceeded that along the minor axis, Kb, by more than a factor of 10. We further show that along-patch shear was predominantly a result of differential advection between a deep residual flow to the south-east and an oscillating wind-driven surface Ekman layer. In this region of strong low frequency (inertial) shear a time dependent model of shear dispersion (Young et al., 1982) was able to account for the observed rate of horizontal dispersion calculated on the target isopycnal surface.
Xu, Bruce S; Lollar, Barbara Sherwood; Passeport, Elodie; Sleep, Brent E
2016-04-15
Aqueous phase diffusion-related isotope fractionation (DRIF) for carbon isotopes was investigated for common groundwater contaminants in systems in which transport could be considered to be one-dimensional. This paper focuses not only on theoretically observable DRIF effects in these systems but introduces the important concept of constraining "observable" DRIF based on constraints imposed by the scale of measurements in the field, and on standard limits of detection and analytical uncertainty. Specifically, constraints for the detection of DRIF were determined in terms of the diffusive fractionation factor, the initial concentration of contaminants (C0), the method detection limit (MDL) for isotopic analysis, the transport time, and the ratio of the longitudinal mechanical dispersion coefficient to effective molecular diffusion coefficient (Dmech/Deff). The results allow a determination of field conditions under which DRIF may be an important factor in the use of stable carbon isotope measurements for evaluation of contaminant transport and transformation for one-dimensional advective-dispersive transport. This study demonstrates that for diffusion-dominated transport of BTEX, MTBE, and chlorinated ethenes, DRIF effects are only detectable for the smaller molar mass compounds such as vinyl chloride for C0/MDL ratios of 50 or higher. Much larger C0/MDL ratios, corresponding to higher source concentrations or lower detection limits, are necessary for DRIF to be detectable for the higher molar mass compounds. The distance over which DRIF is observable for VC is small (less than 1m) for a relatively young diffusive plume (<100years), and DRIF will not easily be detected by using the conventional sampling approach with "typical" well spacing (at least several meters). With contaminant transport by advection, mechanical dispersion, and molecular diffusion this study suggests that in field sites where Dmech/Deff is larger than 10, DRIF effects will likely not be
Technology Transfer Automated Retrieval System (TEKTRAN)
Understanding and modeling transport of solutes in porous media is a critical issue in the environmental protection. Contaminants from various industrial and agricultural sources can travel in soil and ground water and eventually affect human and animal health. The parabolic advective-dispersive equ...
NASA Astrophysics Data System (ADS)
Benson, D. A.; Zhang, Y.
2006-12-01
Conservative solute transport through natural media is typically "anomalous" or non-Fickian. The anomalous transport may be characterized by faster than linear growth of the centered second moment, or non-Gaussian leading or trailing edges of a plume emanating from a point source. These characteristics develop because of non-local dependence on either past (time) or far upstream (space) concentrations. Non-local equations developed to describe anomalous dispersion usually focus on constant transport parameters and/or independence of the transport on space dimension. These simplifications have been useful for fitting simple transport processes, such as laboratory column tests or 1-D projections of field data. However, they may be insufficient for real field settings, where direction-dependent depositional processes and nonstationary heterogeneity can occur. We develop a generalized, multi-dimensional, spatiotemporal fractional advection- dispersion equation (fADE) with variable parameters to characterize regional-scale anomalous dispersion processes including trapping in immobile zones and/or super-Fickian rapid transport. A Lagrangian numerical model of the space-time fractional transport equation is developed in which solute particles can disperse in both space and time, depending on the medium heterogeneity properties, such as the connectivity and statistical distributions of high versus low-permeability deposits. In the generalized fADE, the range of the order of fractional time derivative is (0 2], representing a wide range of possible trapping behavior. The extension of the order to the range (1 2] is novel to transport theory. We apply the numerical model in 1-D and 2-D to the MADE site tritium plumes, and results indicate that this method can capture the main behaviors of realistic plumes, including local variations of spreading, direction-dependent scaling rates, and arbitrary rapid transport along preferential flow paths. Since the governing equation
CELL DENSITY AND NON-EQUILIBRIUM SORPTION EFFECTS ON BACTERIAL DISPERSAL IN GROUNDWATER MICROCOSMS
The relative importance of dispersion, physical straining, non-equilibrium sorption, and cell density on the dispersal of bacteria was examined in saturated, flow-dynamic sand columns. The bacterial breakthrough as a was followed by measuring the effluent concentration of 3H-aden...
EVALUATION OF SORPTION MODELS IN THE SIMULATION OF NAPHTHALENE TRANSPORT THROUGH SATURATED SOILS
To show the effect of sorption model selection. he one-dimensional transport of naphthalene in two saturated soils was simulated using numerical models, which included the processes of dispersion, advection and sorption. hree different models describing the sorption process: 1) a...
NASA Astrophysics Data System (ADS)
Dean, A. M.; Benson, D. A.; Major, E.
2010-12-01
By adding a fractional-in-time term to the traditional advection dispersion equation, a model is able to simulate a late-time heavy-tailed contaminant breakthrough curve. This heavy-tailed breakthrough curve is observed in data collected during a conservative tracer “push-pull” test at the Macrodispersion Experiment (MADE) site. A time fractional advection dispersion equation (fADE) is able to predict power law tailing of conservative solutes by accounting for solutes transferring between the mobile and relatively immobile phases. Solutes can become trapped in a low permeability zone where the transport is controlled by diffusion instead of advection. It has been observed that the late-time heavy-tailed breakthrough curve may follow a power law due to the movement into these low flow zones. By solving the time fADE in a particle tracking program (SLIM-FAST) the model accounts for mass transfer between various phases and produces the same power law tail as observed in field data. For the implementation of the time fADE, in SLIM-FAST, the particles move based on a random-walk motion but have the ability to transition into a relatively immobile phase after (exponentially) random mobile times. Following a period in the immobile phase, the particle re-enters the mobile phase to be moved by advection and Fickian dispersion. To test the fADE approach, a recent single-well push-pull tracer test at the MADE site is reproduced using a groundwater flow code (ParFlow) and a particle tracking code (SLIM-FAST) using various immobile residence-time distributions.
Brine heterogeneity and dispersed interstitial advective flow underneath the sea of galilee, Israel
NASA Astrophysics Data System (ADS)
Weinstein, Y.; Katz, A.; Kastner, M.; Nishri, A.; Jannasch, H.
2003-04-01
Saline groundwater from submerged sources is today the main source of salts to the Sea of Galilee, supplying annually 72,000 tons chloride to the lake. MOSQUITO flux-meters were deployed during April to September 2001 at seven shallow Kinneret sites in order to study the dispersed interstitial flow. Each instrument carried 3-5 Osmo-Samplers that continuously sampled pore fluids at various depths in the sediment (0-45 cm). Samples were analyzed for their chemistry and for concentration of Na-fluorescein that was previously injected into the sediment. In general, oncentrations of conservative elements (e.g. Cl, Na, B) increase with depth into the sediment at all sites. However, concentrations vary significantly from one site to another. For instance, in sub-lacustrine brines next to the Tiberias Hot Springs Cl concentration reaches more than 14,500 mg/l at 35 cm below lake floor (not much less than the 18,000 mg/l in the nearby on-shore springs), while at other stations, brines are diluted by meteoric water to less than 1,500 mg/l Cl. Ion ratios in pore water indicate that the shallow parts of Lake Kinneret are underlain by several, separate, brine pockets, that are sometimes located very close to each other and discharge to the same area. Pore water at the east and northwest of the lake have Na/Cl ionic ratios between 0.7 and 0.8, similar to that of the overlying lake water, while at the west and south, ratios are significantly lower (<0.6 and <0.5, respectively), indicating larger degree of evaporation of the original brine end-member. Brines next to Tiberias Hot Springs have significantly higher Br/Cl and lower Mg/Cl ratios than pore water from other sites. Eastern shore sub-lacusrine brines (next to Gofra) are distinguished by their very high Sr/Cl and B/Cl ratios. Advective flow rates were derived from temporal patterns of Na-fluorescein concentration in pore water. Flow rates were between 0 and 80 cm/yr. Annual fluxes through the shallow part of the lake are
Qi, Yong; Zhang, Tian C
2016-04-01
Sorption of hormones to soil particles of different size fractions (DSFs) has been studied to understand their fate and transport (F/T) in soils. Conventional studies fractionated the soil particles into DSFs by using the high speed stirring method and/or adding surfactants to fully disperse the bulk soil. However, the natural processes (e.g., soil erosion, irrigation) often are relatively mild, and many soil particles may be still in the aggregate form. In this study, a method was developed for conducting the sorption test of a representative hormone (i.e., testosterone) to bulk soils first and then analyzing the results against DSFs. Results indicated the particle size distribution (PSD) of the two representative soils tested with partially-dispersed and fully-dispersed methods was significantly different due to the attachment of clay particles on sand and silt. Testosterone was sorbed mainly by the dominant aggregates even though they might have relatively lower sorption affinity than that of clays. However, the small particles (<2000 nm), even with ∼5% mass of the bulk soil, contributed more than 30% of sorbed testosterone in the "whole" soils. The partially-dispersed soil particles of DSFs should be used to understand the transport of hormone in runoff, because using the fully-dispersed soil particles will overestimate while the whole soil method will underestimate the transport potential. With the methodology developed in this study, the sorption tests will not compromise soil's original properties (e.g., aggregates) or the competition (e.g., sorption) among soil particles, and the contribution of DSFs (particularly the partially-dispersed aggregates) to the sorption of the "whole" soil can be determined. PMID:26826645
Downes, Barbara J; Lancaster, Jill
2010-01-01
1. In advection-dominated systems (both freshwater and marine), population dynamics are usually presumed to be dominated by the effects of migrants dispersing by advection, especially over the small spatial scales at which populations can be studied, but few studies have tested this presumption. We tested the hypothesis that benthic densities are controlled by densities of dispersers for two aquatic insects in upland streams. 2. Our study animals were two species of caddisflies (Hydropsychidae), which become sedentary filter-feeders following settlement onto substrata. Densities of dispersers in the drift (advective dispersal) were quantified using nets placed along the upstream edges of riffles, where the latter abruptly abutted a slower, upstream run. Settlement was estimated at each site using brick pavers, half of which had been fenced to prevent colonization of their top surfaces by walking hydropsychids, thus allowing us to distinguish also the mode of movement during settlement. 3. First through fifth instars of two species, Smicrophylax sp. AV2 and Asmicridea sp. AV1, were abundant and showed disparate results. Drift and settlement were relatively strongly related for Smicrophylax. The best fit lines were shown by second and third instars settling on plain bricks, suggesting that drift played a strong role in settlement, but that some drifters dropped to the bottom and located substrata by walking. Quantile regression suggested that drift sets limits to settlement in this species and that settlement success was highly variable. In contrast, settlement by Asmicridea was poorly related to drift; settlers were mainly individuals re-dispersing within sites. 4. Smicrophylax densities appear to be controlled by dispersal from upstream, but benthic density of Asmicridea is more likely linked to local demography. Our data demonstrate the dangers of assuming that supposedly drift-prone species can all be modelled in the same way. Alternative models emphasizing
NASA Astrophysics Data System (ADS)
Pelosi, A.; Schumer, R.; Parker, G.; Ferguson, R. I.
2016-03-01
Tracer pebbles are often used to study bed load transport processes in gravel bed rivers. Models have been proposed for their downstream dispersion, and also for vertical dispersion, but not for the combined effects of downstream and vertical movement. Here we use the Exner-Based Master Equation to characterize the transient coevolution of streamwise and vertical advection-diffusion of tracer pebbles under equilibrium transport conditions (no net aggradation or degradation). The coevolution of streamwise and vertical dispersion gives rise to behavior that can differ markedly from that associated with purely streamwise processes with no vertical exchange. One example is streamwise advective slowdown. Particles that are advected downward into zones where the probability of reentrainment becomes asymptotically small are essentially trapped and can no longer participate in streamwise advection. As a result, the mean streamwise velocity of the tracer plume declines in time. Qualitative and quantitative comparisons with two field experiments show encouraging agreement despite the simplified boundary conditions in the model.
NASA Astrophysics Data System (ADS)
Cornaton, F.; Perrochet, P.
2006-09-01
Groundwater age and life expectancy probability density functions (pdf) have been defined, and solved in a general three-dimensional context by means of forward and backward advection-dispersion equations [Cornaton F, Perrochet P. Groundwater age, life expectancy and transit time distributions in advective-dispersive systems; 1. Generalized reservoir theory. Adv Water Res (xxxx)]. The discharge and recharge zones transit time pdfs were then derived by applying the reservoir theory (RT) to the global system, thus considering as ensemble the union of all inlet boundaries on one hand, and the union of all outlet boundaries on the other hand. The main advantages in using the RT to calculate the transit time pdf is that the outlet boundary geometry does not represent a computational limiting factor (e.g. outlets of small sizes), since the methodology is based on the integration over the entire domain of each age, or life expectancy, occurrence. In the present paper, we extend the applicability of the RT to sub-drainage basins of groundwater reservoirs by treating the reservoir flow systems as compartments which transfer the water fluxes to a particular discharge zone, and inside which mixing and dispersion processes can take place. Drainage basins are defined by the field of probability of exit at outlet. In this way, we make the RT applicable to each sub-drainage system of an aquifer of arbitrary complexity and configuration. The case of the well-head protection problem is taken as illustrative example, and sensitivity analysis of the effect of pore velocity variations on the simulated ages is carried out.
Knopman, Debra S.; Voss, Clifford I.
1987-01-01
The spatial and temporal variability of sensitivities has a significant impact on parameter estimation and sampling design for studies of solute transport in porous media. Physical insight into the behavior of sensitivities is offered through an analysis of analytically derived sensitivities for the one-dimensional form of the advection-dispersion equation. When parameters are estimated in regression models of one-dimensional transport, the spatial and temporal variability in sensitivities influences variance and covariance of parameter estimates. Several principles account for the observed influence of sensitivities on parameter uncertainty. (1) Information about a physical parameter may be most accurately gained at points in space and time. (2) As the distance of observation points from the upstream boundary increases, maximum sensitivity to velocity during passage of the solute front increases. (3) The frequency of sampling must be 'in phase' with the S shape of the dispersion sensitivity curve to yield the most information on dispersion. (4) The sensitivity to the dispersion coefficient is usually at least an order of magnitude less than the sensitivity to velocity. (5) The assumed probability distribution of random error in observations of solute concentration determines the form of the sensitivities. (6) If variance in random error in observations is large, trends in sensitivities of observation points may be obscured by noise. (7) Designs that minimize the variance of one parameter may not necessarily minimize the variance of other parameters.
Sorption and desorption of atrazine and diuron onto water dispersible soil primary size fractions.
Wang, Peng; Keller, Arturo A
2009-03-01
In this study, a low energy separation method was employed to separate water dispersible clay-, silt-, and sand-sized fractions. The batch equilibrium method was used to conduct atrazine and diuron sorption/desorption experiments with the bulk soils and their size fractions separately. A Freundlich sorption model provided the best fit for all sorption and desorption data. A mass balance calculation, taking into account the pesticide concentration differences in the size fraction and bulk soil, showed that pesticide sorption onto the different size fractions reproduces well the total amount of the pesticide sorbed onto the bulk soils. Due to their higher soil organic carbon content, the clay fractions were much more effective sorbents for the pesticides than the bulk soils, silt, and sand fractions. For all soils, the amount of the pesticide sorbed onto the clay fractions was more than 20% of the total amount of the pesticide sorbed by the bulk soils even though the clay fractions in these soils were only 5.3-14.0% (by weight). The clay fractions had the highest desorption hysteresis among all size fractions and the bulk soils, followed by silt fractions, implying the clay fractions had the strongest bound and least desorbable pesticide molecules. Our results suggest that attention should be paid to the pesticide sorbed to the smallest colloids, the water dispersible fraction, which can be potentially mobilized under field conditions, leading to wide spreading of contamination. PMID:19147172
Healy, R.W.; Russell, T.F.
1992-01-01
A finite-volume Eulerian-Lagrangian local adjoint method for solution of the advection-dispersion equation is developed and discussed. The method is mass conservative and can solve advection-dominated ground-water solute-transport problems accurately and efficiently. An integrated finite-difference approach is used in the method. A key component of the method is that the integral representing the mass-storage term is evaluated numerically at the current time level. Integration points, and the mass associated with these points, are then forward tracked up to the next time level. The number of integration points required to reach a specified level of accuracy is problem dependent and increases as the sharpness of the simulated solute front increases. Integration points are generally equally spaced within each grid cell. For problems involving variable coefficients it has been found to be advantageous to include additional integration points at strategic locations in each well. These locations are determined by backtracking. Forward tracking of boundary fluxes by the method alleviates problems that are encountered in the backtracking approaches of most characteristic methods. A test problem is used to illustrate that the new method offers substantial advantages over other numerical methods for a wide range of problems.
Lancaster, Jill; Downes, Barbara J
2014-12-01
Many communities comprise species that select resources that are patchily distributed in an environment that is otherwise unsuitable or suboptimal. Effects of this patchiness can depend on the characteristics of patch arrays and animal movements, and produce non-intuitive outcomes in which population densities are unrelated to resource abundance. Resource mosaics are predicted to have only weak effects, however, where patches are ephemeral or organisms are transported advectively. The running waters of streams and benthic invertebrates epitomize such systems, but empirical tests of resource mosaics are scarce. We sampled 15 common macroinvertebrates inhabiting distinct detritus patches at four sites within a sand-bed stream, where detritus formed a major resource of food and living space. At each site, environmental variables were measured for 100 leaf packs; invertebrates were counted in 50 leaf packs. Sites differed in total abundance of detritus, leaf pack sizes and invertebrate densities. Multivariate analysis indicated that patch size was the dominant environmental variable, but invertebrate densities differed significantly between sites even after accounting for patch size. Leaf specialists showed positive and strong density-area relationships, except where the patch size range was small and patches were aggregated. In contrast, generalist species had weaker and variable responses to patch sizes. Population densities were not associated with total resource abundance, with the highest densities of leaf specialists in sites with the least detritus. Our results demonstrate that patchy resources can affect species even in communities where species are mobile, have advective dispersal, and patches are relatively ephemeral. PMID:25190216
Cox, T.J.; Runkel, R.L.
2008-01-01
Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme's limitations. ?? 2008 ASCE.
NASA Technical Reports Server (NTRS)
Hung, R. J.; Liaw, G. S.
1980-01-01
The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.
Lichtner, P.C.; Helgeson, H.C.
1986-06-20
A general formulation of multi-phase fluid flow coupled to chemical reactions was developed based on a continuum description of porous media. A preliminary version of the computer code MCCTM was constructed which implemented the general equations for a single phase fluid. The computer code MCCTM incorporates mass transport by advection-diffusion/dispersion in a one-dimensional porous medium coupled to reversible and irreversible, homogeneous and heterogeneous chemical reactions. These reactions include aqueous complexing, oxidation/reduction reactions, ion exchange, and hydrolysis reactions of stoichiometric minerals. The code MCCTM uses a fully implicit finite difference algorithm. The code was tested against analytical calculations. Applications of the code included investigation of the propagation of sharp chemical reaction fronts, metasomatic alteration of microcline at elevated temperatures and pressures, and ion-exchange in a porous column. Finally numerical calculations describing fluid flow in crystalline rock in the presence of a temperature gradient were compared with experimental results for quartzite.
NASA Astrophysics Data System (ADS)
Zaramella, M.; Marion, A.; Lewandowski, J.; Nützmann, G.
2016-07-01
Solute transport in rivers is controlled by surface flow hydrodynamics and by transient storage in dead zones, pockets of vegetation and hyporheic sediments where mass exchange and retention are governed by complex mechanisms. The physics of these processes are generally investigated by optimization of transient storage models (TSMs) to experimental data often yielding inconsistent and equifinal parameter sets. Uncertainty on parameters estimation is found to depend not only on the rates of exchange between the stream and storage zones, the stream-water velocity and the stream reach length according to the experimental Damkohler number (DaI), but also on the relative significance between transient storage and longitudinal dispersion on breakthrough curves (BTCs). An optimization strategy was developed and applied to an experimental dataset obtained from tracer tests in a small lowland river, analyzing BTCs generated through tracer injections under different conditions. The method supplies a tool to estimate model parameters from observed data through the analysis of the relative parameter significance. To analyze model performance a double compartment TSM was optimized by a regular fit procedure based on simple root mean square error minimization and by a fit based on a relative significance analysis of mechanism signatures. As a result consistent longitudinal dispersion and transient storage parameters were obtained when the signature targeted optimization was used.
NASA Astrophysics Data System (ADS)
Ho, D. T.; Schlosser, P.; Schlosser, P.; Schlosser, P.; Caplow, T.; Garrison, M. R.
2001-12-01
In recent years, deliberate tracer release experiments have been used in the ocean, rivers, lakes, and groundwater flow systems to study advection, mixing, air-water gas exchange, and exchanges between subsystems. Here we report results from a recent deliberate tracer release experiment conducted in the tidal Hudson River. On July 25, 2001, ca. 3.3 moles of the inert gas SF6 were injected into the Hudson River near Newburgh, NY at a depth of about 6 m. Subsequently, the SF6 was monitored from a boat (Riverkeeper) by pumping water (from 2 m depth) via a submersed pump mounted on the front of the boat, through a gas extraction unit, followed by measurement using an onboard gas chromatograph. The measurement interval was about 2 minutes and the maximum speed of the boat was about 15 km h-1. This allowed us to obtain detailed surveys of the temporal evolution of the tracer plume for 14 days. Initial results from the experiment show that during July/August 2001, there was virtually no net downward advection of the water body originally tagged with SF6. Instead, we observed rapid mixing of the tracer-tagged water up- and down-river. After one week, the tracer-tagged water could be detected over a stretch of 70 km along the axis of the river channel. At this time, the stretch of river labeled with concentrations >50% of the peak value was about 14 km. After two weeks, the tracer-tagged water had extended to over 90 km, while 28 km had SF6 concentrations >50% of the peak value. Vertical mixing into depressions on the bottom of the river reaching more than 175 feet seemed to be rapid. Dispersion coefficients and vertical turbulent exchange coefficients will be discussed.
Spatial Moment Equations for a Groundwater Plume with Degradation and Rate-Limited Sorption
In this note, we analytically derive the solution for the spatial moments of groundwater solute concentration distributions simulated by a one-dimensional model that assumes advective-dispersive transport with first-order degradation and rate-limited sorption. Sorption kinetics...
Embry, Irucka; Roland, Victor; Agbaje, Oluropo; Watson, Valetta; Martin, Marquan; Painter, Roger; Byl, Tom; Sharpe, Lonnie
2013-01-01
A new residence-time distribution (RTD) function has been developed and applied to quantitative dye studies as an alternative to the traditional advection-dispersion equation (AdDE). The new method is based on a jointly combined four-parameter gamma probability density function (PDF). The gamma residence-time distribution (RTD) function and its first and second moments are derived from the individual two-parameter gamma distributions of randomly distributed variables, tracer travel distance, and linear velocity, which are based on their relationship with time. The gamma RTD function was used on a steady-state, nonideal system modeled as a plug-flow reactor (PFR) in the laboratory to validate themore » effectiveness of the model. The normalized forms of the gamma RTD and the advection-dispersion equation RTD were compared with the normalized tracer RTD. The normalized gamma RTD had a lower mean-absolute deviation (MAD) (0.16) than the normalized form of the advection-dispersion equation (0.26) when compared to the normalized tracer RTD. The gamma RTD function is tied back to the actual physical site due to its randomly distributed variables. The results validate using the gamma RTD as a suitable alternative to the advection-dispersion equation for quantitative tracer studies of non-ideal flow systems.« less
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.; Leroy, D.; Laloui, L.; Descoeudres, F.
2001-10-01
A model was developed, to describe miscible grout propagation in a saturated deformable porous medium, based on Bear's statistical model with spatial volume averaging. In a previous paper, the model was first successfully confronted to one-dimensional laboratory experiments.In the present paper, the numerical model is used to simulate practical grouting operation in a cylindrical injection model. The cylindrical injection model lends itself to study main flow and propagation character istics for a dispersed suspension-type grout, under axisymmetric conditions close to real scale conditions.Comparison between numerical solutions and experimental results is essential to confirm the validity and accuracy of the proposed model from a phenomenological standpoint. The numerical model performances show that the underlying mathematical model constitutes a realistic predictive model reproducing most prominent features during injection of a suspension-type grout into a deformable porous medium. The basic mechanism by which injected miscible grout permeates a soil mass is discussed in detail. Such a tool leads to quality control criteria for grouting on a theoretical basis, which complements existing criteria acquired through engineering practice.
NASA Astrophysics Data System (ADS)
Bouchelaghem, F.; Vulliet, L.
2001-10-01
The development of a predictive model of behaviour of porous media during injection of miscible grout, taking into account convection, dilution and filtration of grout solution with interstitial water, as well as consolidation aspects, is presented. Model assumptions are reviewed and discussed first. During the establishment of the model, we insist on surface terms and their physical relevance in expressing adsorption effects. Constitutive laws such as Fick's law for diffusive mass transport, hydrodynamic dispersion tensor dealing with miscibility, are modified by taking into account filtration effects. A new surface term appears in mass balance equations as a consequence of filtration. According to the filtration laws used, an initial filtration rate is estimated on the basis of a one-dimensional experimental campaign. The field equations are discretized by using Galerkin finite element and -scheme standard method. For transport equation, Streamline Upwind Petrov Galerkin method is employed to prevent numerical oscillations. Lastly, confrontation of numerical results with laboratory experiments constitutes a first step to validate the model on a realistic basis.
Webb, S.W.
1996-05-01
Two models for gas-phase diffusion and advection in porous media, the Advective-Dispersive Model (ADM) and the Dusty-Gas Model (DGM), are reviewed. The ADM, which is more widely used, is based on a linear addition of advection calculated by Darcy`s Law and ordinary diffusion using Fick`s Law. Knudsen diffusion is often included through the use of a Klinkenberg factor for advection, while the effect of a porous medium on the diffusion process is through a porosity-tortuosity-gas saturation multiplier. Another, more comprehensive approach for gas-phase transport in porous media has been formulated by Evans and Mason, and is referred to as the Dusty- Gas Model (DGM). This model applies the kinetic theory of gases to the gaseous components and the porous media (or ``dust``) to develop an approach for combined transport due to ordinary and Knudsen diffusion and advection including porous medium effects. While these two models both consider advection and diffusion, the formulations are considerably different, especially for ordinary diffusion. The various components of flow (advection and diffusion) are compared for both models. Results from these two models are compared to isothermal experimental data for He-Ar gas diffusion in a low-permeability graphite. Air-water vapor comparisons have also been performed, although data are not available, for the low-permeability graphite system used for the helium-argon data. Radial and linear air-water heat pipes involving heat, advection, capillary transport, and diffusion under nonisothermal conditions have also been considered.
Alcaraz, Mar; García-Gil, Alejandro; Vázquez-Suñé, Enric; Velasco, Violeta
2016-02-01
Borehole Heat Exchangers (BHEs) are increasingly being used to exploit shallow geothermal energy. This paper presents a new methodology to provide a response to the need for a regional quantification of the geothermal potential that can be extracted by BHEs and the associated environmental impacts. A set of analytical solutions facilitates accurate calculation of the heat exchange of BHEs with the ground and its environmental impacts. For the first time, advection and dispersion heat transport mechanisms and the temporal evolution from the start of operation of the BHE are taken into account in the regional estimation of shallow geothermal resources. This methodology is integrated in a GIS environment, which facilitates the management of input and output data at a regional scale. An example of the methodology's application is presented for Barcelona, in Spain. As a result of the application, it is possible to show the strengths and improvements of this methodology in the development of potential maps of low temperature geothermal energy as well as maps of environmental impacts. The minimum and maximum energy potential values for the study site are 50 and 1800 W/m(2) for a drilled depth of 100 m, proportionally to Darcy velocity. Regarding to thermal impacts, the higher the groundwater velocity and the energy potential, the higher the size of the thermal plume after 6 months of exploitation, whose length ranges from 10 to 27 m long. A sensitivity analysis was carried out in the calculation of heat exchange rate and its impacts for different scenarios and for a wide range of Darcy velocities. The results of this analysis lead to the conclusion that the consideration of dispersion effects and temporal evolution of the exploitation prevent significant differences up to a factor 2.5 in the heat exchange rate accuracy and up to several orders of magnitude in the impacts generated. PMID:26605833
Romero-González, J; Walton, J C; Peralta-Videa, J R; Rodríguez, E; Romero, J; Gardea-Torresdey, J L
2009-01-15
The biosorption of Cr(III) onto packed columns of Agave lechuguilla was analyzed using an advective-dispersive (AD) model and its analytical solution. Characteristic parameters such as axial dispersion coefficients, retardation factors, and distribution coefficients were predicted as functions of inlet ion metal concentration, time, flow rate, bed density, cross-sectional column area, and bed length. The root-mean-square-error (RMSE) values 0.122, 0.232, and 0.285 corresponding to the flow rates of 1, 2, and 3 (10(-3))dm3min(-1), respectively, indicated that the AD model provides an excellent approximation of the simulation of lumped breakthrough curves for the adsorption of Cr(III) by lechuguilla biomass. Therefore, the model can be used for design purposes to predict the effect of varying operational conditions. PMID:18462882
Merritt, M.L.
1993-01-01
The simulation of the transport of injected freshwater in a thin brackish aquifer, overlain and underlain by confining layers containing more saline water, is shown to be influenced by the choice of the finite-difference approximation method, the algorithm for representing vertical advective and dispersive fluxes, and the values assigned to parametric coefficients that specify the degree of vertical dispersion and molecular diffusion that occurs. Computed potable water recovery efficiencies will differ depending upon the choice of algorithm and approximation method, as will dispersion coefficients estimated based on the calibration of simulations to match measured data. A comparison of centered and backward finite-difference approximation methods shows that substantially different transition zones between injected and native waters are depicted by the different methods, and computed recovery efficiencies vary greatly. Standard and experimental algorithms and a variety of values for molecular diffusivity, transverse dispersivity, and vertical scaling factor were compared in simulations of freshwater storage in a thin brackish aquifer. Computed recovery efficiencies vary considerably, and appreciable differences are observed in the distribution of injected freshwater in the various cases tested. The results demonstrate both a qualitatively different description of transport using the experimental algorithms and the interrelated influences of molecular diffusion and transverse dispersion on simulated recovery efficiency. When simulating natural aquifer flow in cross-section, flushing of the aquifer occurred for all tested coefficient choices using both standard and experimental algorithms. ?? 1993.
Parashar, R.; Cushman, J.H.
2008-06-20
Microbial motility is often characterized by 'run and tumble' behavior which consists of bacteria making sequences of runs followed by tumbles (random changes in direction). As a superset of Brownian motion, Levy motion seems to describe such a motility pattern. The Eulerian (Fokker-Planck) equation describing these motions is similar to the classical advection-diffusion equation except that the order of highest derivative is fractional, {alpha} element of (0, 2]. The Lagrangian equation, driven by a Levy measure with drift, is stochastic and employed to numerically explore the dynamics of microbes in a flow cell with sticky boundaries. The Eulerian equation is used to non-dimensionalize parameters. The amount of sorbed time on the boundaries is modeled as a random variable that can vary over a wide range of values. Salient features of first passage time are studied with respect to scaled parameters.
Healy, R.W.; Russell, T.F.
1993-01-01
Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods for solute transport problems that are dominated by advection. FVELLAM systematically conserves mass globally with all types of boundary conditions. Integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking of characteristic lines intersecting inflow boundaries. FVELLAM extends previous results by obtaining mass conservation locally on Lagrangian space-time elements. -from Authors
Healy, R.W.; Russell, T.F.
1998-01-01
We extend the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) for solution of the advection-dispersion equation to two dimensions. The method can conserve mass globally and is not limited by restrictions on the size of the grid Peclet or Courant number. Therefore, it is well suited for solution of advection-dominated ground-water solute transport problems. In test problem comparisons with standard finite differences, FVELLAM is able to attain accurate solutions on much coarser space and time grids. On fine grids, the accuracy of the two methods is comparable. A critical aspect of FVELLAM (and all other ELLAMs) is evaluation of the mass storage integral from the preceding time level. In FVELLAM this may be accomplished with either a forward or backtracking approach. The forward tracking approach conserves mass globally and is the preferred approach. The backtracking approach is less computationally intensive, but not globally mass conservative. Boundary terms are systematically represented as integrals in space and time which are evaluated by a common integration scheme in conjunction with forward tracking through time. Unlike the one-dimensional case, local mass conservation cannot be guaranteed, so slight oscillations in concentration can develop, particularly in the vicinity of inflow or outflow boundaries. Published by Elsevier Science Ltd.
HIGH TEMPERATURE SORPTION OF CESIUM AND STRONTIUM ON DISPERSED KAOLINITE POWDERS
Sorption of cesium and strontium on kaolinite powders was investigated as a means to minimize the emissions of these metals during certain high-temperature processes currently being developed to isolate and dispose of radiological and mixed wastes. In this work, nonradioactive aq...
Wagner, Brian J.; Gorelick, Steven M.
1986-01-01
A simulation nonlinear multiple-regression methodology for estimating parameters that characterize the transport of contaminants is developed and demonstrated. Finite difference containment transport simulation is combined with a nonlinear weighted least squares multiple-regression procedure. The technique provides optimal parameter estimates and gives statistics for assessing the reliability of these estimates under certain general assumptions about the distributions of the random measurement errors. Monte Carlo analysis is used to estimate parameter reliability for a hypothetical homogeneous soil column for which concentration data contain large random measurement errors. The value of data collected spatially versus data collected temporally was investigated for estimation of velocity, dispersion coefficient, effective porosity, first-order decay rate, and zero-order production. The use of spatial data gave estimates that were 2-3 times more reliable than estimates based on temporal data for all parameters except velocity. (Estimated author abstract) Refs.
NASA Astrophysics Data System (ADS)
Furbish, David Jon; Childs, Elise M.; Haff, Peter K.; Schmeeckle, Mark W.
2009-09-01
We formulate soil grain transport by rain splash as a stochastic advection-dispersion process. By taking into account the intermittency of grain motions activated by raindrop impacts, the formulation indicates that gradients in raindrop intensity, and thus grain activity (the volume of grains in motion per unit area) can be as important as gradients in grain concentration and surface slope in effecting transport. This idea is confirmed by rain splash experiments and manifest in topographic roughening via mound growth beneath desert shrubs. The formulation provides a framework for describing transport and dispersal of any soil material moveable by rain splash, including soil grains, soil-borne pathogens and nutrients, seeds, or debitage. As such it shows how classic models of topographic "diffusion" reflect effects of slope-dependent grain drift, not diffusion, and it highlights the role of rain splash in the ecological behavior of desert shrubs as "resource islands." Specifically, the growth of mounds beneath shrub canopies, where differential rain splash initially causes more grains to be splashed inward beneath the protective canopy than outward, involves the "harvesting" of nearby soil material, including nutrients. Mounds thus represent temporary storage of soil derived from areas surrounding the shrubs. As the inward grain flux associated with differential rain splash is sustained over the shrub lifetime, mound material is effectively sequestered from erosional processes that might otherwise move this material downslope. With shrub death and loss of the protective canopy, differential rain splash vanishes and the mound material is dispersed to the surrounding area, again subject to downslope movement.
Actinide and technetium sorption on iron-silicate and dispersed clay colloids
Shade, J.W.; Ames, L.L.; McGarrah, J.E.
1984-01-01
Two different colloidal suspensions, representative of those found in waste package interaction tests, were prepared from iron metal and silica powders or sodium-bentonite at 90/sup 0/C. Aliquots were spiked with /sup 233/U, /sup 235/Np, /sup 237/Pu, or /sup 95m/Tc at pH ranges from 2 to 12, then shaken for 24 hours followed by a 15 A filtration. Zeta potential measurements were made on unspiked samples. Similar sorptive properties were observed for both colloids. At 25/sup 0/C both /sup 233/U and /sup 237/Pu exhibit maximum sorption (50-90%) near pH 6. Sorption drops by about a factor of 5 at pH > 8. Slight sorption of /sup 235/Np occurs at pH 11 and decreases to zero at lower pH values. /sup 95m/Tc does not sorb on Fe-silicates and is only slightly sorbed (10%) on smectites. 11 references, 4 figures, 2 tables.
NASA Astrophysics Data System (ADS)
Parker, Jack C.; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Parker, Jack C; Kim, Ungtae
2015-11-01
The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of
Dontsova, Katerina M; Hayes, Charolett; Pennington, Judith C; Porter, Beth
2009-01-01
Explosives in soils can present environmental problems for military installations. Fine, mobile particles represent the most reactive fraction of the soil and, therefore, are expected to adsorb explosives and potentially facilitate their transport. The objective of this study was to determine the relative significance of phyllosilicate clay, organic matter, and two forms of extractable iron in adsorption of 2,4,6-trinitrotoluene (TNT) and hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) by the colloidal water-dispersible clay (WDC) fraction of the soil. The WDC fraction of two mineral and one organic soil was separated and then treated to remove organic carbon (OC) and several forms of iron (Fe(o), oxalate extractable, and Fe(d), dithionite-citrate extractable). Adsorption coefficients were determined for whole soils, untreated, and treated WDC. For mineral soils, adsorption of TNT and RDX on the WDC was greater than on the whole soil. The presence of OC increased explosives sorption by WDC. When OC was removed, iron interfered with TNT sorption. In the presence of OC, removal of Fe(o) decreased RDX adsorption and increased TNT adsorption indicating different adsorption mechanisms. Organic carbon was a more significant indicator of explosives adsorption by WDC than clays or iron oxides and hydroxides. Therefore, OC is the most likely medium for facilitated transport of TNT and RDX. PMID:19465721
Guyonvarch, Estelle; Ramin, Elham; Kulahci, Murat; Plósz, Benedek Gy
2015-10-15
The present study aims at using statistically designed computational fluid dynamics (CFD) simulations as numerical experiments for the identification of one-dimensional (1-D) advection-dispersion models - computationally light tools, used e.g., as sub-models in systems analysis. The objective is to develop a new 1-D framework, referred to as interpreted CFD (iCFD) models, in which statistical meta-models are used to calculate the pseudo-dispersion coefficient (D) as a function of design and flow boundary conditions. The method - presented in a straightforward and transparent way - is illustrated using the example of a circular secondary settling tank (SST). First, the significant design and flow factors are screened out by applying the statistical method of two-level fractional factorial design of experiments. Second, based on the number of significant factors identified through the factor screening study and system understanding, 50 different sets of design and flow conditions are selected using Latin Hypercube Sampling (LHS). The boundary condition sets are imposed on a 2-D axi-symmetrical CFD simulation model of the SST. In the framework, to degenerate the 2-D model structure, CFD model outputs are approximated by the 1-D model through the calibration of three different model structures for D. Correlation equations for the D parameter then are identified as a function of the selected design and flow boundary conditions (meta-models), and their accuracy is evaluated against D values estimated in each numerical experiment. The evaluation and validation of the iCFD model structure is carried out using scenario simulation results obtained with parameters sampled from the corners of the LHS experimental region. For the studied SST, additional iCFD model development was carried out in terms of (i) assessing different density current sub-models; (ii) implementation of a combined flocculation, hindered, transient and compression settling velocity function; and (iii
CONTAMINANT TRANSPORT IN SEDIMENT UNDER THE INFLUENCE OF ADVECTIVE FLUX
Chemical flux across the sediment/water interface is controlled by a combination of diffusive, dispersive and advective processes. The advective process is a function of submarine groundwater discharge and tidal effects. In areas where surface water interacts with groundwater, ...
Fast strontium transport induced by hydrodynamic dispersion and pH-dependent sorption
NASA Astrophysics Data System (ADS)
Prigiobbe, Valentina; Hesse, Marc A.; Bryant, Steven L.
2012-09-01
As a fluid carries solutes through a porous material, species that sorb onto the surface of the material travel more slowly than the fluid. Stronger adsorption results in slower solute migration, or increased solute retardation. The adsorption of strontium (Sr2+) onto iron-oxides is strongly pH-dependent and becomes significant at high pH. Radioactive Sr2+ isotopes are, therefore, commonly stored in alkaline solutions to maximize their retardation. Field observations and numerical simulations of the leakage of such solutions into low-pH soils, however, show that even Sr2+ stored in alkaline solutions can migrate without retardation. Migration occurs because hydrodynamic dispersion allows mixing of Sr2+ with the low-pH fluid forming an acidic Sr2+-rich plume which can travel without retardation. Here we report the first experimental observations confirming this dispersion-induced fast Sr2+ transport. We report column-flood experiments where a high-pH solution containing Sr2+ was injected into a low-pH porous medium of iron-oxide-coated beads. We observe both a strongly retarded Sr2+ front and an isolated fast pulse of Sr2+ traveling at the average fluid velocity. This dispersion-induced fast pulse of strontium must be taken into account when considering the safety of radionuclide storage in alkaline solutions.
This manual describes the next generation of the modular three-dimensional transport model, MT3D, with significantly expanded capabilities, including the addition of (a) a third-order total-variation-diminishing (TVD) scheme for solving the advection term that is mass conservativ...
Olin, T.J.; Myers, T.E.; Townsend, D.M.
1996-09-01
The sorption and transformation behavior of 2,4,6-trinitrotoluene (TNT) is important to modeling and remediation efforts at military installations where subsurface contamination exists in connection with munitions production. Processes potentially affecting the fate and transport of TNT in soils and groundwater include biotic and abiotic transformation, sorption, advection, hydrodynamic dispersion, dissolution, diffusion, and facilitated transport by organic and inorganic colloids (McGrath 1995). TNT breakthrough curves may provide indications of the type of processes occurring. The transformation rate of TNT is of particular interest in determining the long-term risk associated with TNT contamination in a soil.
Monger, Gregg R.; Duncan, Candice Morrison; Brusseau, Mark L.
2015-01-01
A gas-phase tracer test (GTT) was conducted at a landfill in Tucson, AZ, to help elucidate the impact of landfill gas generation on the transport and fate of chlorinated aliphatic volatile organic contaminants (VOCs). Sulfur hexafluoride (SF6) was used as the non-reactive gas tracer. Gas samples were collected from a multiport monitoring well located 15.2 m from the injection well, and analyzed for SF6, CH4, CO2, and VOCs. The travel times determined for SF6 from the tracer test are approximately two to ten times smaller than estimated travel times that incorporate transport by only gas-phase diffusion. In addition, significant concentrations of CH4 and CO2 were measured, indicating production of landfill gas. Based on these results, it is hypothesized that the enhanced rates of transport observed for SF6 are caused by advective transport associated with landfill gas generation. The rates of transport varied vertically, which is attributed to multiple factors including spatial variability of water content, refuse mass, refuse permeability, and gas generation. PMID:26380532
NASA Astrophysics Data System (ADS)
Pang, Liping; Goltz, Mark; Close, Murray
2003-01-01
In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data.
Pang, Liping; Goltz, Mark; Close, Murray
2003-01-01
In this note, we applied the temporal moment solutions of [Das and Kluitenberg, 1996. Soil Sci. Am. J. 60, 1724] for one-dimensional advective-dispersive solute transport with linear equilibrium sorption and first-order degradation for time pulse sources to analyse soil column experimental data. Unlike most other moment solutions, these solutions consider the interplay of degradation and sorption. This permits estimation of a first-order degradation rate constant using the zeroth moment of column breakthrough data, as well as estimation of the retardation factor or sorption distribution coefficient of a degrading solute using the first moment. The method of temporal moment (MOM) formulae was applied to analyse breakthrough data from a laboratory column study of atrazine, hexazinone and rhodamine WT transport in volcanic pumice sand, as well as experimental data from the literature. Transport and degradation parameters obtained using the MOM were compared to parameters obtained by fitting breakthrough data from an advective-dispersive transport model with equilibrium sorption and first-order degradation, using the nonlinear least-square curve-fitting program CXTFIT. The results derived from using the literature data were also compared with estimates reported in the literature using different equilibrium models. The good agreement suggests that the MOM could provide an additional useful means of parameter estimation for transport involving equilibrium sorption and first-order degradation. We found that the MOM fitted breakthrough curves with tailing better than curve fitting. However, the MOM analysis requires complete breakthrough curves and relatively frequent data collection to ensure the accuracy of the moments obtained from the breakthrough data. PMID:12498577
NASA Astrophysics Data System (ADS)
Tóth, Gábor; Keppens, Rony
2012-07-01
The Versatile Advection Code (VAC) is a freely available general hydrodynamic and magnetohydrodynamic simulation software that works in 1, 2 or 3 dimensions on Cartesian and logically Cartesian grids. VAC runs on any Unix/Linux system with a Fortran 90 (or 77) compiler and Perl interpreter. VAC can run on parallel machines using either the Message Passing Interface (MPI) library or a High Performance Fortran (HPF) compiler.
2001-01-01
The ability of species to migrate and disperse is a trait that has interested ecologists for many years. Now that so many species and ecosystems face major environmental threats from habitat fragmentation and global climate change, the ability of species to adapt to these changes by dispersing, migrating, or moving between patches of habitat can be crucial to ensuring their survival. This book provides a timely and wide-ranging overview of the study of dispersal and incorporates much of the latest research. The causes, mechanisms, and consequences of dispersal at the individual, population, species and community levels are considered. The potential of new techniques and models for studying dispersal, drawn from molecular biology and demography, is also explored. Perspectives and insights are offered from the fields of evolution, conservation biology and genetics. Throughout the book, theoretical approaches are combined with empirical data, and care has been taken to include examples from as wide a range of species as possible.
Concentration through large advection
NASA Astrophysics Data System (ADS)
Aleja, D.; López-Gómez, J.
2014-11-01
In this paper we extend the elegant results of Chen, Lam and Lou [6, Section 2], where a concentration phenomenon was established as the advection blows up, to a general class of adventive-diffusive generalized logistic equations of degenerate type. Our improvements are really sharp as we allow the carrying capacity of the species to vanish in some subdomain with non-empty interior. The main technical devices used in the derivation of the concentration phenomenon are Proposition 3.2 of Cano-Casanova and López-Gómez [5], Theorem 2.4 of Amann and López-Gómez [1] and the classical Harnack inequality. By the relevance of these results in spatial ecology, complete technical details seem imperative, because the proof of Theorem 2.2 of [6] contains some gaps originated by an “optimistic” use of Proposition 3.2 of [5]. Some of the general assumptions of [6] are substantially relaxed.
LAYER DEPENDENT ADVECTION IN CMAQ
The advection methods used in CMAQ require that the Courant-Friedrichs-Lewy (CFL) condition be satisfied for numerical stability and accuracy. In CMAQ prior to version 4.3, the ADVSTEP algorithm established CFL-safe synchronization and advection timesteps that were uniform throu...
Competitive sorption of cis-DCE and TCE in silica gel as a model porous mineral solid.
Avila, Manuel Alejandro Salaices; Breiter, Roman
2008-08-01
The competitive sorption of 1,2-cis-dichloroethene (cis-DCE) and trichloroethene (TCE) was investigated by means of column experiments using a model porous mineral solid represented by silica gel. The experimental isotherms were obtained by employing a chromatographic method. The competitive sorption isotherms were modelled with the extended Freundlich and extended Langmuir isotherms, using the parameters from single-solute experiments. The breakthrough curves were modelled with the advection-dispersion transport equation coupled with the lumped pore diffusion model. The best results were obtained when the extended Freundlich isotherm was employed. The competitive sorption was revealed with the presence of an overshoot in the breakthrough curve of cis-DCE and a decrease in the degree of sorption of cis-DCE (20%) and TCE (12%). A linear dependency of the overshoot with an increase in the concentration of cis-DCE at a fixed concentration of TCE was observed, between 16% and 20%, and at least at concentrations <6 mg L(-1) in the liquid phase. The displaced molecules of cis-DCE by TCE were accumulated through the column causing its overshoot; thus short columns may hinder its observation. Thermodynamic analysis shows an exothermic adsorption process of -34 to -41 kJ mol(-1), which is enhanced by sorption in micropores. The Gibbs free energy is positive for cis-DCE in the multi-component case, due to its displacement by TCE. PMID:18541287
Pfeifer,H.; Gueye-Girardet, A.; Reymond, D.; Schlegel, C.; Temgoua, E.; Hesterberg, D.; Chou, J.
2004-01-01
In recent years, elevated arsenic concentrations have been found in waters and soils of many countries, often resulting in a health threat for the local population. Switzerland is not an exception and this paper deals with the release and subsequent fate of arsenic in a 200-km{sup 2} mountainous watershed, characterized by crystalline silicate rocks (gneisses, schists, amphibolites) that contain abundant As-bearing sulfide ore deposits, some of which have been mined for iron and gold in the past. Using analytical methods common for mineralogical, ground water and soil studies (XRD, XRF, XAS-XANES and -EXAFS, electron microprobe, extraction, ICP, AAS with hydride generator, ion chromatography), seven different field situations and related dispersion processes of natural arsenic have been studied: (1) release by rock weathering; (2) transport and deposition by water and ice; (3) release of As to the ground and surface water due to increasing pH; (4) accumulation in humic soil horizons; (5) remobilization by reduction in water-saturated soils and stagnant ground waters; (6) remobilization by using P-rich fertilizers or dung and (7) oxidation, precipitation and dilution in surface waters. Comparison of the results with experimental adsorption studies and speciation diagrams from the literature allows us to reconstruct and identify the typical behavior of arsenic in a natural environment under temperate climatic conditions. The main parameters identified are: (a) once liberated from the primary minerals, sorption processes on Fe-oxy-hydroxides dominate over Al-phases, such as Al-hydroxides or clay minerals and limit the As concentrations in the spring and well waters between 20 and 300 {mu}g/l. (b) Precipitation as secondary minerals is limited to the weathering domain, where the As concentrations are still high and not yet too diluted by rain and soils waters. (c) Although neutral and alkaline pH conditions clearly increase the mobility of As, the main factor to
Experiments in Advective and Turbulent Hyporheic Pumping
NASA Astrophysics Data System (ADS)
Mccluskey, A. H.; Grant, S.; Stewardson, M. J.
2014-12-01
Hyporheic exchange (HE) is the mixing of stream and subsurface waters beneath the sediment-water interface (SWI). At the patch and reach scales, HE is dominated by periodic upwelling and downwelling zones, induced by pressure variation and processes within the turbulent boundary layer (TBL). This can be caused by (1) the geometry of the stream, imposing a stationary wave at the SWI or (2) by a travelling wave associated with the propagation of turbulent pressure waves generated from the TBL. Case (1) has generally been the favoured model of hyporheic exchange and has been referred to as hyporheic 'pumping' by Elliott and Brooks, and subsequently others. Case (2) can be termed turbulent pumping, and has been proposed as a mechanism to model the combined effects of turbulent dispersion alongside steady-state advection. While this has been represented numerically and analytically, conjecture remains about the physical representation of these combined processes. We present initial results from experiments undertaken to classify the spatial and temporal characteristics of pressure variation at and beneath the SWI, with a periodic sinusoidal geometry of wavelength 0.28m and height 0.02m. As an initial characterisation, the advective flow profile has been examined using time-lapse photography of dyes released across the span of a periodic downwelling zone. These tracer tests confirmed delineation of isolated upwelling and downwelling cells as noted by previous authors in modelling studies. However, their distribution deviates from the typical pumping pattern with increased discharge and stream gradient. Empirical orthogonal function (EOF) analysis of high frequency (250Hz) pressure measurements, sampled at an array along the centroid of the flume underneath one wavelength gave further insight into the spatial distribution of turbulent signatures arising from roughness-generated turbulence. A turbulent frequency of 6-10Hz dominates, however the penetration depth appears to
Radionuclide migration laboratory studies for validation of batch sorption data
Triay, I.R.; Mitchell, A.J.; Ott, M.A.
1991-12-31
Advective and diffusive migration experiments (within the Dynamic Transport Column Experiments and Diffusion Studies of the Yucca Mountain Site Characterization Project) involve utilizing crushed material, intact, and fractured tuff in order to test and improve (if necessary) transport models by experimentally observing the migration of sorbing and non-sorbing radionuclides on a laboratory scale. Performing a validation of the sorption data obtained with batch techniques (within the Batch Sorption Study) is an integral part of the mission of the Dynamic Transport Column Experiments and Diffusion Studies. In this paper the work scope of the radionuclide migration laboratory experiments (as they apply to validation of batch sorption data) is reviewed.
NASA Astrophysics Data System (ADS)
Glynn, P. D.
2002-12-01
One-dimensional (1D) geochemical transport modeling is used to demonstrate the effects of speciation and sorption reactions on the ground-water transport of Np and Pu, two redox-sensitive elements. Earlier 1D simulations (Reardon, 1981) considered the kinetically-limited dissolution of calcite and its effect on ion-exchange reactions (involving 90Sr, Ca, Na, Mg and K), and documented the spatial variation of a 90Sr partition coefficient under both transient and steady-state chemical conditions. In contrast, the simulations presented here assume local equilibrium for all reactions, and consider sorption on constant potential, rather than constant charge, surfaces. Reardon's (1981) findings documenting the spatial and temporal variability of 90Sr partitioning are reexamined and found partially caused by his assumption of a kinetically-limited reaction. In the present simulations, sorption is assumed the only retardation process controlling Pu and Np transport, and is modeled using a diffuse-double-layer-surface-complexation model. Transport simulations consider the inflow of Np- and Pu-contaminated waters into an initially uncontaminated environment, followed by the cleanup of the resultant contamination with uncontaminated water. Simulations are conducted using different spatial distributions of sorption capacities (with the same total potential sorption capacity, i.e. the same total number of sorption sites, but with different variances and spatial correlation structures). A case with a spatially uniform distribution of sorption capacities was also simulated. Results obtained differ markedly from those that would be obtained in transport simulations using constant Kd, Langmuir, or Freundlich sorption models. When possible, simulation results (breakthrough curves) are fitted to a constant Kd advection-dispersion transport model and compared to each other. Functional differences are often great enough that they prevent a meaningful fit of the simulation results with
Goode, D.J.; Konikow, L.F.
1989-01-01
The U.S. Geological Survey computer model of two-dimensional solute transport and dispersion in ground water (Konikow and Bredehoeft, 1978) has been modified to incorporate the following types of chemical reactions: (1) first-order irreversible rate-reaction, such as radioactive decay; (2) reversible equilibrium-controlled sorption with linear, Freundlich, or Langmuir isotherms; and (3) reversible equilibrium-controlled ion exchange for monovalent or divalent ions. Numerical procedures are developed to incorporate these processes in the general solution scheme that uses method-of- characteristics with particle tracking for advection and finite-difference methods for dispersion. The first type of reaction is accounted for by an exponential decay term applied directly to the particle concentration. The second and third types of reactions are incorporated through a retardation factor, which is a function of concentration for nonlinear cases. The model is evaluated and verified by comparison with analytical solutions for linear sorption and decay, and by comparison with other numerical solutions for nonlinear sorption and ion exchange.
NASA Astrophysics Data System (ADS)
Glynn, Pierre D.
2003-04-01
One-dimensional (1D) geochemical transport modeling is used to demonstrate the effects of speciation and sorption reactions on the ground-water transport of Np and Pu, two redox-sensitive elements. Earlier 1D simulations (Reardon, 1981) considered the kinetically limited dissolution of calcite and its effect on ion-exchange reactions (involving 90Sr, Ca, Na, Mg and K), and documented the spatial variation of a 90Sr partition coefficient under both transient and steady-state chemical conditions. In contrast, the simulations presented here assume local equilibrium for all reactions, and consider sorption on constant potential, rather than constant charge, surfaces. Reardon's (1981) seminal findings on the spatial and temporal variability of partitioning (of 90Sr) are reexamined and found partially caused by his assumption of a kinetically limited reaction. In the present work, sorption is assumed the predominant retardation process controlling Pu and Np transport, and is simulated using a diffuse-double-layer-surface-complexation (DDLSC) model. Transport simulations consider the infiltration of Np- and Pu-contaminated waters into an initially uncontaminated environment, followed by the cleanup of the resultant contamination with uncontaminated water. Simulations are conducted using different spatial distributions of sorption capacities (with the same total potential sorption capacity, but with different variances and spatial correlation structures). Results obtained differ markedly from those that would be obtained in transport simulations using constant Kd, Langmuir or Freundlich sorption models. When possible, simulation results (breakthrough curves) are fitted to a constant Kd advection-dispersion transport model and compared. Functional differences often are great enough that they prevent a meaningful fit of the simulation results with a constant Kd (or even a Langmuir or Freundlich) model, even in the case of Np, a weakly sorbed radionuclide under the simulation
Glynn, P.D.
2003-01-01
One-dimensional (1D) geochemical transport modeling is used to demonstrate the effects of speciation and sorption reactions on the ground-water transport of Np and Pu, two redox-sensitive elements. Earlier 1D simulations (Reardon, 1981) considered the kinetically limited dissolution of calcite and its effect on ion-exchange reactions (involving 90Sr, Ca, Na, Mg and K), and documented the spatial variation of a 90Sr partition coefficient under both transient and steady-state chemical conditions. In contrast, the simulations presented here assume local equilibrium for all reactions, and consider sorption on constant potential, rather than constant charge, surfaces. Reardon's (1981) seminal findings on the spatial and temporal variability of partitioning (of 90Sr) are reexamined and found partially caused by his assumption of a kinetically limited reaction. In the present work, sorption is assumed the predominant retardation process controlling Pu and Np transport, and is simulated using a diffuse-double-layer-surface-complexation (DDLSC) model. Transport simulations consider the infiltration of Np- and Pu-contaminated waters into an initially uncontaminated environment, followed by the cleanup of the resultant contamination with uncontaminated water. Simulations are conducted using different spatial distributions of sorption capacities (with the same total potential sorption capacity, but with different variances and spatial correlation structures). Results obtained differ markedly from those that would be obtained in transport simulations using constant Kd, Langmuir or Freundlich sorption models. When possible, simulation results (breakthrough curves) are fitted to a constant K d advection-dispersion transport model and compared. Functional differences often are great enough that they prevent a meaningful fit of the simulation results with a constant K d (or even a Langmuir or Freundlich) model, even in the case of Np, a weakly sorbed radionuclide under the
Super-diffusion versus competitive advection: a simulation
NASA Astrophysics Data System (ADS)
Del Moro, D.; Giannattasio, F.; Berrilli, F.; Consolini, G.; Lepreti, F.; Gošić, M.
2015-04-01
Context. Magnetic element tracking is often used to study the transport and diffusion of the magnetic field on the solar photosphere. From the analysis of the displacement spectrum of these tracers, it has recently been agreed that a regime of super-diffusivity dominates the solar surface. Quite habitually this result is discussed in the framework of fully developed turbulence. Aims: However, the debate whether the super-diffusivity is generated by a turbulent dispersion process, by the advection due to the convective pattern, or even by another process is still open, as is the question of the amount of diffusivity at the scales relevant to the local dynamo process. Methods: To understand how such peculiar diffusion in the solar atmosphere takes place, we compared the results from two different data sets (ground-based and space-borne) and developed a simulation of passive tracers advection by the deformation of a Voronoi network. Results: The displacement spectra of the magnetic elements obtained by the data sets are consistent in retrieving a super-diffusive regime for the solar photosphere, but the simulation also shows a super-diffusive displacement spectrum: its competitive advection process can reproduce the signature of super-diffusion. Conclusions: Therefore, it is not necessary to hypothesize a totally developed turbulence regime to explain the motion of the magnetic elements on the solar surface.
NASA Astrophysics Data System (ADS)
Ancey, Christophe; Bohorquez, Patricio; Heyman, Joris
2016-04-01
The advection-diffusion equation arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Stochastic models can also be used to derive this equation, with the significant advantage that they provide information on the statistical properties of particle activity. Stochastic models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. We develop an approach based on birth-death Markov processes, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received little attention. We show that particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due to velocity fluctuations), with the important consequence that local measurements depend on both the intrinsic properties of particle displacement and the dimensions of the measurement system.
Regenerative Sorption Refrigerator
NASA Technical Reports Server (NTRS)
Jones, Jack A.; Wen, Liang-Chi; Bard, Steven
1991-01-01
Two-stage sorption refrigerator achieves increased efficiency via regenerative-heating concept in which waste heat from praseodymium/cerium oxide (PCO) chemisorption compressor runs charcoal/krypton (C/Kr) sorption compressor. Waste heat from each PCO sorption compressor used to power surrounding C/Kr sorption compressor. Flows of heat in two compressor modules controlled by gas-gap thermal switches. Has no wearing moving parts other than extremely long life, room-temperature check valves operating about twice per hour. Virtually no measurable vibration, and has potential operating life of at least ten years.
Comparison of neptunium sorption results using batch and column techniques
Triay, I.R.; Furlano, A.C.; Weaver, S.C.; Chipera, S.J.; Bish, D.L.
1996-08-01
We used crushed-rock columns to study the sorption retardation of neptunium by zeolitic, devitrified, and vitric tuffs typical of those at the site of the potential high-level nuclear waste repository at Yucca Mountain, Nevada. We used two sodium bicarbonate waters (groundwater from Well J-13 at the site and water prepared to simulate groundwater from Well UE-25p No. 1) under oxidizing conditions. It was found that values of the sorption distribution coefficient, Kd, obtained from these column experiments under flowing conditions, regardless of the water or the water velocity used, agreed well with those obtained earlier from batch sorption experiments under static conditions. The batch sorption distribution coefficient can be used to predict the arrival time for neptunium eluted through the columns. On the other hand, the elution curves showed dispersivity, which implies that neptunium sorption in these tuffs may be nonlinear, irreversible, or noninstantaneous. As a result, use of a batch sorption distribution coefficient to calculate neptunium transport through Yucca Mountain tuffs would yield conservative values for neptunium release from the site. We also noted that neptunium (present as the anionic neptunyl carbonate complex) never eluted prior to tritiated water, which implies that charge exclusion does not appear to exclude neptunium from the tuff pores. The column experiments corroborated the trends observed in batch sorption experiments: neptunium sorption onto devitrified and vitric tuffs is minimal and sorption onto zeolitic tuffs decreases as the amount of sodium and bicarbonate/carbonate in the water increases.
Surfzone alongshore advective accelerations: observations and modeling
NASA Astrophysics Data System (ADS)
Hansen, J.; Raubenheimer, B.; Elgar, S.
2014-12-01
The sources, magnitudes, and impacts of non-linear advective accelerations on alongshore surfzone currents are investigated with observations and a numerical model. Previous numerical modeling results have indicated that advective accelerations are an important contribution to the alongshore force balance, and are required to understand spatial variations in alongshore currents (which may result in spatially variable morphological change). However, most prior observational studies have neglected advective accelerations in the alongshore force balance. Using a numerical model (Delft3D) to predict optimal sensor locations, a dense array of 26 colocated current meters and pressure sensors was deployed between the shoreline and 3-m water depth over a 200 by 115 m region near Duck, NC in fall 2013. The array included 7 cross- and 3 alongshore transects. Here, observational and numerical estimates of the dominant forcing terms in the alongshore balance (pressure and radiation-stress gradients) and the advective acceleration terms will be compared with each other. In addition, the numerical model will be used to examine the force balance, including sources of velocity gradients, at a higher spatial resolution than possible with the instrument array. Preliminary numerical results indicate that at O(10-100 m) alongshore scales, bathymetric variations and the ensuing alongshore variations in the wave field and subsequent forcing are the dominant sources of the modeled velocity gradients and advective accelerations. Additional simulations and analysis of the observations will be presented. Funded by NSF and ASDR&E.
Advection fog formation in a polluted atmosphere
Hung, R.J.; Liaw, G.S.
1981-01-01
Large quantities of atmospheric aerosols with composition SO/sub 4//sup 2 -/, NO/sub 3//sup -/ and NH/sub 4//sup +/ have been detected in highly industrialized areas. The major portions of aerosol products are the results of energy related fuel combustion. Both microphysical and macrophysical processes are considered in investigating the time dependent evolution of the saturation spectra of condensation nuclei associated with both polluted and clean atmospheres during the time periods of advection fog formation. The results show that the condensation nuclei associated with a polluted atmosphere provide more favorable conditions than condensation nuclei associated with a clean atmosphere to produce dense advection fog, and that attaining a certain degree of supersaturation is not necessarily required for the formation of advection fog with condensation nuclei associated with a polluted atmosphere for monodisperse distribution.
Antidiffusive velocities for multipass donor cell advection
Margolin, L.; Smolarkiewicz, P.K.
1999-01-01
Multidimensional positive definite advection transport algorithm (MPDATA) is an iterative process for approximating the advection equation, which uses a donor cell approximation to compensate for the truncation error of the originally specified donor cell scheme. This step may be repeated an arbitrary number of times, leading to successfully more accurate solutions to the advection equation. In this paper, the authors show how to sum the successive approximations analytically to find a single antidiffusive velocity that represents the effects of an arbitrary number of passes. The analysis is first done in one dimension to illustrate the method and then is repeated in two dimensions. The existence of cross terms in the truncation analysis of the two-dimensional equations introduces an extra complication into the calculation. The authors discuss the implementation of the antidiffusive velocities and provide some examples of applications, including a third-order accurate scheme.
MAGNETIC ADVECTION DUE TO DIFFUSIVITY GRADIENTS
NASA Astrophysics Data System (ADS)
Zita, E. J.
2009-12-01
We derive and discuss an important source of advection of magnetic fields in plasmas, for a completely general case. Magnetic diffusivity is proportional to electrical resistivity: where the value this parameter is high, it is well known that magnetic fields can leak (or diffuse) rapidly into (or out) of the plasma. Magnetohydrodynamic lore has it that where gradients, or changes in space, of the value of the diffusivity are high, magnetic fields can have enhanced flow (or advection). We derive this phenomenon rigorously, compare our results to other work in the literature, and discuss its implications, especially for kinematic dynamos. As an extra mathematical bonus, we find that the magnetic advection due to diffusivity gradients can be expressed in terms of a diffusion equation within the induction equation, making its computational implementation especially simple.
Efficient mass transport by optical advection
NASA Astrophysics Data System (ADS)
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-10-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms.
Kinetics of the sorption of triterpene saponin by hypercrosslinked polystyrene
NASA Astrophysics Data System (ADS)
Mironenko, N. V.; Brezhneva, T. A.; Selemenev, V. F.
2013-03-01
The kinetics of sorption of triterpene saponin by the polymer sorbent NM-200 is considered. The influence of the surface activity of glycoside on the rate of formation and structure of the adsorption layer on the sorbent's surface is established. The rate-determining step of sorption is found to be diffusion into the sorbent grain. The value of the activation energy demonstrates the determining role of dispersion forces in the interaction between triterpene saponin and the polymer sorbent MN-200.
SEPs Dropout Events Associated with Advected Interplanetary Magnetic Structures
NASA Astrophysics Data System (ADS)
Bruno, R.; Trenchi, L.; Telloni, D.; D'Amicis, R.; Marcucci, F.; Zurbuchen, T.; Weberg, M. J.
2013-05-01
The intensity profile of energetic particles from impulsive solar flares (SEP) often shows abrupt dropouts affecting all energies simultaneously, without time-dispersion. Part of the community thinks that these modulations are directly related to the presence of magnetic structures with a different magnetic topology advected by the wind, a sort of magnetic flux tubes. During the expansion, following the dynamical interaction between plasma regions travelling at different speed, these structures would be partially tangled up in a sort of spaghetti-like bundle. These flux tubes would be alternatively connected or not connected with the flare site and, consequently, they would be filled or devoid of SEPs. When the observer passes through them, he would observe clear particles dropout signatures. We will report about results from a detailed analysis of SEP events which showed several signatures in the local magnetic field and/or plasma parameters associated with SEP modulations. These findings corroborate the idea of a possible link between these particles events observed at the Earth's orbit and magnetic connection or disconnection of the ambient magnetic field with the flare region at the Sun. We will also discuss the advantages represented by future Solar Orbiter in-situ observations. As a matter of fact, Solar Orbiter, from its orbital vantage point during the quasi corotation phase, will be a priviledged observer of this kind of phenomenon since it will observe the advected structure of the solar wind not yet reprocessed by dynamical interaction due to wind expansion.
3D Flow Visualization Using Texture Advection
NASA Technical Reports Server (NTRS)
Kao, David; Zhang, Bing; Kim, Kwansik; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
Texture advection is an effective tool for animating and investigating 2D flows. In this paper, we discuss how this technique can be extended to 3D flows. In particular, we examine the use of 3D and 4D textures on 3D synthetic and computational fluid dynamics flow fields.
Evolution of Advection Upstream Splitting Method Schemes
NASA Technical Reports Server (NTRS)
Liou, Meng-Sing
2010-01-01
This paper focuses on the evolution of advection upstream splitting method(AUSM) schemes. The main ingredients that have led to the development of modern computational fluid dynamics (CFD) methods have been reviewed, thus the ideas behind AUSM. First and foremost is the concept of upwinding. Second, the use of Riemann problem in constructing the numerical flux in the finite-volume setting. Third, the necessity of including all physical processes, as characterised by the linear (convection) and nonlinear (acoustic) fields. Fourth, the realisation of separating the flux into convection and pressure fluxes. The rest of this review briefly outlines the technical evolution of AUSM and more details can be found in the cited references. Keywords: Computational fluid dynamics methods, hyperbolic systems, advection upstream splitting method, conservation laws, upwinding, CFD
Passive advection in a collisionless plasma
NASA Astrophysics Data System (ADS)
Kanekar, Anjor; Schekochihin, Alexander; Hammett, Greg; Dorland, William; Loureiro, Nuno
2014-10-01
We consider a simple kinetic model for the evolution of the particle distribution function in a magnetized turbulent plasma that includes both phase mixing (Landau damping) and advection by a stochastic velocity field: a ``kinetic passive scalar'' in the Batchelor regime. The advection due to stochastic velocity field allows for a stochastic version of the plasma echo by coupling the ``phase-mixing'' and the ``un-phase-mixing'' components of the free energy. We have developed a new analytical framework to diagnose the efficiency of such coupling. We have also developed a new GPU code named Gandalf that solves this kinetic model numerically. In this poster, we shall present numerical and analytical results related to this work.
Antidiffusive velocities for multipass donor cell advection
Margolin, L.G. ); Smolarkiewicz, P.K. )
1989-12-01
Smolarkiewicz describes an iterative process for approximating the advection equation. Basically, he uses a donor cell approximation to correct for the truncation error of the originally specified donor cell scheme. This step may be repeated an arbitrary number of times leading to successively more accurate solutions to the advection equation. In this report, we show how to sum the successive approximations analytically to find a single antidiffusive velocity that represents the effects of an arbitrary number of passes. The analysis is first done dimension to illustrate the method. The analysis is then repeated in two dimensions. The existence of cross terms in the truncation analysis of the two-dimensional equations introduces an extra complication into the calculation. We discuss the implementation of our new antidiffusive velocities and provide some examples of applications. 6 refs., 5 figs., 4 tabs.
Distributed Parallel Particle Advection using Work Requesting
Muller, Cornelius; Camp, David; Hentschel, Bernd; Garth, Christoph
2013-09-30
Particle advection is an important vector field visualization technique that is difficult to apply to very large data sets in a distributed setting due to scalability limitations in existing algorithms. In this paper, we report on several experiments using work requesting dynamic scheduling which achieves balanced work distribution on arbitrary problems with minimal communication overhead. We present a corresponding prototype implementation, provide and analyze benchmark results, and compare our results to an existing algorithm.
High Order Semi-Lagrangian Advection Scheme
NASA Astrophysics Data System (ADS)
Malaga, Carlos; Mandujano, Francisco; Becerra, Julian
2014-11-01
In most fluid phenomena, advection plays an important roll. A numerical scheme capable of making quantitative predictions and simulations must compute correctly the advection terms appearing in the equations governing fluid flow. Here we present a high order forward semi-Lagrangian numerical scheme specifically tailored to compute material derivatives. The scheme relies on the geometrical interpretation of material derivatives to compute the time evolution of fields on grids that deform with the material fluid domain, an interpolating procedure of arbitrary order that preserves the moments of the interpolated distributions, and a nonlinear mapping strategy to perform interpolations between undeformed and deformed grids. Additionally, a discontinuity criterion was implemented to deal with discontinuous fields and shocks. Tests of pure advection, shock formation and nonlinear phenomena are presented to show performance and convergence of the scheme. The high computational cost is considerably reduced when implemented on massively parallel architectures found in graphic cards. The authors acknowledge funding from Fondo Sectorial CONACYT-SENER Grant Number 42536 (DGAJ-SPI-34-170412-217).
Efficient mass transport by optical advection
Kajorndejnukul, Veerachart; Sukhov, Sergey; Dogariu, Aristide
2015-01-01
Advection is critical for efficient mass transport. For instance, bare diffusion cannot explain the spatial and temporal scales of some of the cellular processes. The regulation of intracellular functions is strongly influenced by the transport of mass at low Reynolds numbers where viscous drag dominates inertia. Mimicking the efficacy and specificity of the cellular machinery has been a long time pursuit and, due to inherent flexibility, optical manipulation is of particular interest. However, optical forces are relatively small and cannot significantly modify diffusion properties. Here we show that the effectiveness of microparticle transport can be dramatically enhanced by recycling the optical energy through an effective optical advection process. We demonstrate theoretically and experimentally that this new advection mechanism permits an efficient control of collective and directional mass transport in colloidal systems. The cooperative long-range interaction between large numbers of particles can be optically manipulated to create complex flow patterns, enabling efficient and tunable transport in microfluidic lab-on-chip platforms. PMID:26440069
NASA Astrophysics Data System (ADS)
Mallén, G.; Maloszewski, P.; Flynn, R.; Rossi, P.; Engel, M.; Seiler, K.-P.
2005-05-01
The bacteria Escherichia coli and Pseudomonas putida, and the bacteriophage virus H40/1 are examined both for their transport behaviour relative to inert solute tracers and for their modelability under natural flow conditions in a gravel aquifer. The microbes are attenuated in the following sequence: H40/1≥ P. putida≫ E. coli. The latter is desorbed almost completely within a few days. Breakthrough and recovery curves of the simultaneously injected non-reactive tracers are simulated with the 2D and 1D dispersion equation, in order to ascertain longitudinal dispersivity ( αL) and mean flow time ( T0). Mathematical modelling is difficult due to the aquifer heterogeneity, which results in preferential flow paths between injection and observation wells. Therefore, any attempt of fitting the dispersion model (DM) to the entire inert-tracer breakthrough curve (BTC) fails. Adequate fitting of the model to measured data only succeeds using a DM consisting of a superposition of several BTCs, each representing another set of flow paths. This gives rise to a multimodal, rather than a Gaussian groundwater velocity distribution. Only hydraulic parameters derived from the fastest partial curve, which is fitted to the rising part of the Uranine BTC, are suitable to model microbial breakthroughs. The hydraulic parameters found using 2D and 1D models were nearly identical. Their values were put into an analytical solution of 1D advective-dispersive transport combined with two-site reaction model introduced by Cameron and Klute [Cameron, D.R., Klute, A., 1977. Convective-dispersive solute transport with a combined equilibrium and kinetic adsorption model. Water Resour. Res. 13, 183-189], in order to identify reactive transport parameters (sorption/desorption) and attenuation mechanisms for the microbes migration. This shows that the microbes are almost entirely transported through preferential flow paths, which are represented by the first partial curve. Inert tracers, however
Westrich, H.R.; Brady, P.V.; Cygan, R.T.; Nagy, K.L.; Anderson, H.L.
1997-03-01
A key issue in performance assessment of low-level radioactive waste sites is predicting the transport and retardation of radionuclides through local soils under a variety of hydrologic and geochemical conditions. Improved transport codes should include a mechanistic model of radionuclide retardation. The authors have been investigating metal sorption (Cs{sup +}, Sr{sup 2+}, and Ba{sup 2+}) on a simple clay mineral (kaolinite) to better understand the geochemical interactions of common soil minerals with contaminated groundwaters. These studies include detailed characterizations of kaolinite surfaces, experimental adsorption measurements, surface complexation modeling, and theoretical simulations of cation sorption. The aluminol edge (010) site has been identified as the most likely site for metal sorption on kaolinite in natural solutions. Relative metal binding strengths decrease from Ba{sup 2+} to Sr{sup 2+} to Cs{sup +}, with some portion sorbed on both kaolinite edges and basal surfaces. Some Cs{sup +} also appears to be irreversibly sorbed on both sites. Molecular dynamics simulations suggest that Cs{sup +} is sorbed at aluminol (010) edge sites as an inner-sphere complex and weakly sorbed as an outer-sphere complex on (001) basal surfaces. These results provide the basis to understand and predict metal sorption onto kaolinite, and a framework to characterize sorption processes on more complex clay minerals.
NASA Astrophysics Data System (ADS)
Zhuravkov, S. P.; Lobanova, G. L.; Martemiyanov, D. V.; Nadeina, L. V.
2015-04-01
Experimental results of physicochemical and sorption properties of material samples obtained by electro-spark dispersion in water and by conductor electric explosion in argon are shown in the paper. Due to comparison of investigated samples sorption activity under static conditions, it was able to determine the most effective samples in the process of extraction of Cu2+ ions from aqueous solutions.
Hydrophobic conjugated microporous polymers for sorption of human serum albumin
NASA Astrophysics Data System (ADS)
Zheng, Chunli; Du, Miaomiao; Feng, Shanshan; Sun, Hanxue; Li, An; He, Chi; Zhang, TianCheng; Wang, Qiaorui; Wei, Wei
2016-02-01
This paper investigated the sorption of human serum albumin (HSA) from water by three kinds of conjugated microporous polymers (CMPs) with surface hydrophobicity and intrinsic porosity. It was found that the three CMPs captured HSA with fast sorption kinetics and good working capacity. Equilibrium was obtained at 80 min for all the tests, and the maximum sorption quantity (qm) ranged from 0.07 to 0.14 mg/mg. With the increase in the particle external surface area of the CMPs, a greater extent of HSA sorption was achieved. Moreover, promoting the dispersion of CMPs in HSA aqueous solution was also beneficial to the extraction. Attenuated Total Reflection Fourier Transform Infrared spectroscopy verified the interactions between the CMPs and the Nsbnd H, Cdbnd O, and Csbnd N groups of HSA. This paper might provide fundamental guidance for the practical application of CMPs to proteins separation and recovery.
Status Of Sorption Cryogenic Refrigeration
NASA Technical Reports Server (NTRS)
Jones, Jack A.
1988-01-01
Report reviews sorption refrigeration. Developed for cooling infrared detectors, cryogenic research, and other advanced applications, sorption refrigerators have few moving parts, little vibration, and lifetimes of 10 years or more. Describes types of sorption stages, multistage and hybrid refrigeration systems, power requirements, cooling capacities, and advantages and disadvantages of various stages and systems.
Advection and diffusion in shoreline change prediction
NASA Astrophysics Data System (ADS)
Anderson, T. R.; Frazer, L. N.
2010-12-01
We added longshore advection and diffusion to the simple cross-shore rate calculation method, as used widely by the USGS and others, to model historic shorelines and to predict future shoreline positions; and applied this to Hawaiian Island beach data. Aerial photographs, sporadically taken throughout the past century, yield usable, albeit limited, historic shoreline data. These photographs provide excellent spatial coverage, but poor temporal resolution, of the shoreline. Due to the sparse historic shoreline data, and the many natural and anthropogenic events influencing coastlines, we constructed a simplistic shoreline change model that can identify long-term behavior of a beach. Our new, two-dimensional model combines the simple rate method to accommodate for cross-shore sediment transport with the classic Pelnard-Considère model for diffusion, as well as a longshore advection speed term. Inverse methods identify cross-shore rate, longshore advection speed, and longshore diffusivity down a sandy coastline. A spatial averaging technique then identifies shoreline segments where one parameter can reasonably account for the cross-shore and longshore transport rates in that area. This produces model results with spatial resolution more appropriate to the temporal spacing of the data. Because changes in historic data can be accounted for by varying degrees of cross-shore and longshore sediment transport - for example, beach erosion can equally be explained by sand moving either off-shore or laterally - we tested several different model scenarios on the data: allowing only cross-shore sediment movement, only longshore movement, and a combination of the two. We used statistical information criteria to determine both the optimal spatial resolution and best-fitting scenario. Finally, we employed a voting method predicting the relaxed shoreline position over time.
NASA Technical Reports Server (NTRS)
Barrington, A. E.; Caruso, A. J.
1970-01-01
Modified sorption trap for use in high vacuum systems contains provisions for online regeneration of sorbent material. Trap is so constructed that it has a number of encapsulated resistance heaters and a valving and pumping device for removing gases from heated sorbing material. Excessive downtime is eliminated with this trap.
Kinetics of sorption of uranium(VI) compounds with zirconium-silica nanosorbents
NASA Astrophysics Data System (ADS)
Perlova, O. V.; Sazonova, V. F.; Perlova, N. A.; Yaroshenko, N. A.
2014-06-01
The kinetics of sorption of uranium(VI) compounds from sulfate and carbonate solutions using four samples of mesoporous zirconium-silica nanosorbents obtained by bitemplate (solubilization) synthesis was studied. The sorption equilibrium set-in time and the kinetic characteristics of sorption were shown to depend on the sorbent (its composition, specific surface area, dispersity, and pore size), the temperature, and the composition and pH of the solution from which uranium compounds are sorbed. The sorption kinetics was described by a first-order equation. The limiting stage of the process was found to be the external diffusion of uranium-containing particles to the sorbent surface.
Coupled ensemble flow line advection and analysis.
Guo, Hanqi; Yuan, Xiaoru; Huang, Jian; Zhu, Xiaomin
2013-12-01
Ensemble run simulations are becoming increasingly widespread. In this work, we couple particle advection with pathline analysis to visualize and reveal the differences among the flow fields of ensemble runs. Our method first constructs a variation field using a Lagrangian-based distance metric. The variation field characterizes the variation between vector fields of the ensemble runs, by extracting and visualizing the variation of pathlines within ensemble. Parallelism in a MapReduce style is leveraged to handle data processing and computing at scale. Using our prototype system, we demonstrate how scientists can effectively explore and investigate differences within ensemble simulations. PMID:24051840
Advective turbulent transport in the fluid plasma
NASA Astrophysics Data System (ADS)
Min, Byung-Hoon; An, Chan-Yong; Kim, Chang-Bae
2013-10-01
The Hasegawa-Wakatani model (HWM) has been employed in pedagogical analyses of the physics behind the behavior of the tokamak plasmas. In addition to the geometric simplicity HWM has an appealing feature of sustaining autonomous quasi-steady state, unstable modes providing the power that is being transported by the nonlinear interactions and is eventually dissipated by the collisional damping at small scales. Emergence of the zonal flow out of the turbulence is a main candidate to cause the transition from the low plasma confinement to the high mode. In the study of such LH transition with the HWM, the adiabaticity parameter has been shown to play an important role in forcing the zonal flow that results in the regulation of the drift-wave turbulence. Instead of concentrating on the physics of the feedback loop between the turbulence and the zonal flow the present study focuses on the presence of the advective transport of the energy. Numerical simulations of HWM are performed and the connections between the advective transport and the zonal flow will be presented. This work was supported by the Supercpmputing Center/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC-2013-C1-009).
Advection, diffusion and delivery over a network
Heaton, Luke L.M.; López, Eduardo; Maini, Philip K.; Fricker, Mark D.; Jones, Nick S.
2014-01-01
Many biological, geophysical and technological systems involve the transport of resource over a network. In this paper we present an algorithm for calculating the exact concentration of resource at any point in space or time, given that the resource in the network is lost or delivered out of the network at a given rate, while being subject to advection and diffusion. We consider the implications of advection, diffusion and delivery for simple models of glucose delivery through a vascular network, and conclude that in certain circumstances, increasing the volume of blood and the number of glucose transporters can actually decrease the total rate of glucose delivery. We also consider the case of empirically determined fungal networks, and analyze the distribution of resource that emerges as such networks grow over time. Fungal growth involves the expansion of fluid filled vessels, which necessarily involves the movement of fluid. In three empirically determined fungal networks we found that the minimum currents consistent with the observed growth would effectively transport resource throughout the network over the time-scale of growth. This suggests that in foraging fungi, the active transport mechanisms observed in the growing tips may not be required for long range transport. PMID:23005783
Waves, advection, and cloud patterns on Venus
NASA Technical Reports Server (NTRS)
Schinder, Paul J.; Gierasch, Peter J.; Leroy, Stephen S.; Smith, Michael D.
1990-01-01
The stable layers adjacent to the nearly neutral layer within the Venus clouds are found to be capable of supporting vertically trapped, horizontally propagating waves with horizontal wavelengths of about 10 km and speeds of a few meters per second relative to the mean wind in the neutral layer. These waves may possibly be excited by turbulence within the neutral layer. Here, the properties of the waves, and the patterns which they might produce within the visible clouds if excited near the subsolar point are examined. The patterns can be in agreement with many features in images. The waves are capable of transferring momentum latitudinally to help maintain the general atmospheric spin, but at present we are not able to evaluate wave amplitudes. We also examine an alternative possibility that the cloud patterns are produced by advection and shearing by the mean zonal and meridional flow of blobs formed near the equator. It is concluded that advection and shearing by the mean flow is the most likely explanation for the general pattern of small scale striations.
Kukačka, Libor; Nosek, Štĕpán; Kellnerová, Radka; Jurčáková, Klára; Jaňour, Zbyněk
2012-01-01
The objective of this study is to determine processes of pollution ventilation in the X-shaped street intersection in an idealized symmetric urban area for the changing approach flow direction. A unique experimental setup for simultaneous wind tunnel measurement of the flow velocity and the tracer gas concentration in a high temporal resolution is assembled. Advective horizontal and vertical scalar fluxes are computed from averaged measured velocity and concentration data within the street intersection. Vertical advective and turbulent scalar fluxes are computed from synchronized velocity and concentration signals measured in the plane above the intersection. All the results are obtained for five approach flow directions. The influence of the approach flow on the advective and turbulent fluxes is determined. The contribution of the advective and turbulent flux to the ventilation is discussed. Wind direction with the best dispersive conditions in the area is found. The quadrant analysis is applied to the synchronized signals of velocity and concentration fluctuation to determine events with the dominant contribution to the momentum flux and turbulent scalar flux. PMID:22649290
Sorption of strontium on bentonite
Ali Khan, S.; Riaz-ur-Rehman; Ali Khan, M.
1995-12-31
Sorption of Sr on bentonite was studied using the batch technique. Distribution coefficients (K{sub d}) were determined as a function of contact time, pH, sorbent and sorbate concentration and temperature. The data were interpreted in terms of Freundlich, Langmuir and Dubinin-Radushkevich isotherms. Thermodynamic parameters for the sorption system were determined at three different temperatures. The positive value of the heat of sorption, {Delta}H{degree} = 30.62 kJ/mol at 295 K, shows that the sorption of strontium on bentonite is endothermic. The negative value of the free energy of sorption, {Delta}G{degree} = {minus}10.69 kJ/mol at 298 K, shows the spontaneity of the reaction. {Delta}G{degree} becomes more negative with increasing temperature, which shows that the sorption process is more favorable at higher temperatures. The mean free energy for sorption, E {approximately} 9 kJ/mol, suggests that ion exchange is the predominant mode of sorption in the Sr concentration range studies, i.e., 0.01--0.3 mol/dm{sup 3}. The presence of complementary cations depresses the sorption of strontium on bentonite in the order Ca{sup 2+}>Mg{sup 2+}>K{sup +}>Na{sup +}. Some organic complexing agents and natural ligands also affect the sorption of strontium. The desorption studies with ground water at low strontium loads on bentonite show that about 90% of Sr is irreversibly sorbed on the bentonite.
TECHNETIUM SORPTION MEDIA REVIEW
DUNCAN JB; KELLY SE; ROBBINS RA; ADAMS RD; THORSON MA; HAASS CC
2011-08-25
This report presents information and references to aid in the selection of 99Tc sorption media for feasibility studies regarding the removal of 99Tc from Hanford's low activity waste. The report contains literature search material for sorption media (including ion exchange media) for the most tested media to date, including SuperLig 639, Reillex HPQ, TAM (Kruion), Purolite A520E and A530E, and Dowex 1X8. The U.S. Department of Energy (DOE), Office of River Protection (ORP) is responsible for management and completion of the River Protection Project (RPP) mission, which comprises both the Hanford Site tank farms and the Waste Treatment and Immobilization Plant (WTP). The RPP mission is to store, retrieve and treat Hanford's tank waste; store and dispose of treated wastes; and close the tank farm waste management areas and treatment facilities in a safe, environmentally compliant, cost-effective and energy-effective manner.
BUOYANT ADVECTION OF GASES IN UNSATURATED SOIL
Seely, Gregory E.; Falta, Ronald W.; Hunt, James R.
2010-01-01
In unsaturated soil, methane and volatile organic compounds can significantly alter the density of soil gas and induce buoyant gas flow. A series of laboratory experiments was conducted in a two-dimensional, homogeneous sand pack with gas permeabilities ranging from 110 to 3,000 darcy. Pure methane gas was injected horizontally into the sand and steady-state methane profiles were measured. Experimental results are in close agreement with a numerical model that represents the advective and diffusive components of methane transport. Comparison of simulations with and without gravitational acceleration permits identification of conditions where buoyancy dominates methane transport. Significant buoyant flow requires a Rayleigh number greater than 10 and an injected gas velocity sufficient to overcome dilution by molecular diffusion near the source. These criteria allow the extension of laboratory results to idealized field conditions for methane as well as denser-than-air vapors produced by volatilizing nonaqueous phase liquids trapped in unsaturated soil. PMID:20396624
NASA Astrophysics Data System (ADS)
Comolli, Alessandro; Moussey, Charlie; Dentz, Marco
2016-04-01
Transport processes in groundwater systems are strongly affected by the presence of heterogeneity. The heterogeneity leads to non-Fickian features, that manifest themselves in the heavy-tailed breakthrough curves, as well as in the non-linear growth of the mean squared displacement and in the non-Gaussian plumes of solute particles. The causes of non-Fickian transport can be the heterogeneity in the flow fields and the processes of mass exchange between mobile and immobile phases, such as sorption/desorption reactions and diffusive mass transfer. Here, we present a Continuous Time Random Walk (CTRW) model that describes the transport of solutes in d-dimensional systems by taking into account both heterogeneous advection and mobile-immobile mass transfer. In order to account for these processes in the CTRW, the heterogeneities are mapped onto a distribution of transition times, which can be decomposed into advective transition times and trapping times, the latter being treated as a compound Poisson process. While advective transition times are related to the Eulerian flow velocities and, thus, to the conductivity distribution, trapping times depend on the sorption/desorption time scale, in case of reactive problems, or on the distribution of diffusion times in the immobile zones. Since the trapping time scale is typically much larger than the advective time scale, we observe the existence of two temporal regimes. The pre-asymptotic regime is defined by a characteristic time scale at which the properties of transport are fully determined by the heterogeneity of the advective field. On the other hand, in the asymptotic regime both the heterogeneity and the mass exchange processes play a role in conditioning the behaviour of transport. We consider different scenarios to discuss the relative importance of the advective heterogeneity and the mass transfer for the occurrence of non-Fickian transport. For each case we calculate analytically the scalings of the breakthrough
Multi-moment advection scheme in three dimension for Vlasov simulations of magnetized plasma
Minoshima, Takashi; Matsumoto, Yosuke; Amano, Takanobu
2013-03-01
We present an extension of the multi-moment advection scheme [T. Minoshima, Y. Matsumoto, T. Amano, Multi-moment advection scheme for Vlasov simulations, Journal of Computational Physics 230 (2011) 6800–6823] to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.
Simulation of the advective methane transport and AOM in Shenhu area, the Northern South China Sea
NASA Astrophysics Data System (ADS)
Liu, L.; Wu, N.
2012-04-01
Anaerobic Oxidation of Methane (AOM) occurs in the transition zone between the presence of sulfate and methane. This reaction is an important process for methane and the global carbon cycle. Methane gas hydrates bearing sediments were recovered in Shenhu Area, the Northern South China Sea, and methane advective transport was detected in this area as well. A one dimension numerical simulation tool was implemented to study the AOM process combined with the advective methane transport in Shenhu Area according to the local drilling data and geochemical information. The modeled results suggest that local methane flux will be consumed in the sediment column via dissolution, sorption and AOM reaction. A portion of methane will enter water column and possibly atmosphere if the methane flux was one order of magnitude higher than current level. Furthermore, the calculated rates of AOM in Shenhu area range similar to that of gas hydrate mounds in Mexico Golf. However, AOM is ability to consume more methane than that in Golf of Mexico due to the lower permeable sediment associated with a deeper sulfate methane transition layer.
Relative dispersion in the atmosphere
NASA Astrophysics Data System (ADS)
LaCasce, Joe; Graff, Lise; Guttu, Sigmund
2014-05-01
The relative dispersion of pairs of particles in flows is of central importance when describing environmental dispersion, for example of volcanic ash. Atmospheric relative dispersion was examined previously in two balloon experiments in the Southern Hemisphere (the EOLE and TWERLE experiments). In both cases, the dispersion at scales below 1000 km grew exponentially in time, indicating the kinetic energy spectrum is steep. Subsequent analyses suggested though that the dispersion had a power law dependence on time, implying a shallower kinetic energy spectrum. The results from studies employing synthetic particles advected by reanalysis winds are similarly inconsistent, with indications of exponential growth in some cases and power law growth in others. Here we use a different statistic---the probability density function (PDF) of pair displacements---to study dispersion the dispersion of large numbers of synthetic particles, advected by ERA-Interim reanalysis winds. The particles were deployed in the troposphere and stratosphere, both in the tropics and the extra-tropics. We examine the PDFs for the different deployments and compare them to analytical expressions derived for different turbulent inertial ranges. In line with the earlier balloon experiments, the results indicate exponential growth at the sub-deformation (1000 km) scales. At larger scales, the dispersion is anisotropic (predominantly zonal) and pair motion becomes decorrelated. Structure functions calculated from the wind data are in line with these conclusions.
NASA Astrophysics Data System (ADS)
Appadu, A. R.
2016-06-01
An unconditionally positive definite scheme has been derived in [1] to approximate a linear advection-diffusion-reaction equation which models exponential travelling waves and the coefficients of advective, diffusive and reactive terms have been chosen as one. The scheme has been baptised as Unconditionally Positive Finite Difference (UPFD). In this work, we use the UPFD scheme to solve the advection-diffusion-reaction problem in [1] and we also extend our study to three other important regimes involved in this model. The temporal step size is varied while fixing the spatial step size. We compute some errors namely; L1 error, dispersion, dissipation errors. We also study the variation of the modulus of the exact amplification factor, modulus of amplification factor of the scheme and relative phase error, all vs the phase angle for the four different regimes.
Hourly and daily evapotranspiration of alfalfa under regional advection
Technology Transfer Automated Retrieval System (TEKTRAN)
Regional advection often affects the evapotranspiration rates of irrigated crops in the Southern High Plains. In 1998, during a 10-day period (13-22 June) of unusually strong advection, high evapotranspiration (ET) rates for unstressed, irrigated alfalfa (Medicago sativa) were measured with two prec...
Polzer, W.L.; Fuentes, H.R.; Essington, E.H.; Roensch, F.R.
1985-01-01
Sorption isotherms are derived from batch equilibrium data for cobalt, cesium and strontium on Bandelier Tuff. Experiments were conducted at an average temperature of 23/sup 0/C and equilibrium was defined at 48 hours. The solute concentrations ranged from 0 to 500 mg/L. The radioactive isotopes /sup 60/Co, /sup 137/Cs, and /sup 85/Sr were used to trace the sorption of the stable solutes. The Linear, Langmuir, Freundlich and a Modified Freundlich isotherm equations are evaluated. The Modified Freundlich isotherm equation is validated as a preferred general mathematical tool for representing the sorption of the three solutes. The empirical constants derived from the Modified Freundlich isotherm equation indicate that under dynamic flow conditions strontium will move most rapidly and cobalt least rapidly. On the other hand, chemical dispersion will be greatest for cesium and least for strontium. Hill Plots of the sorption data suggest that in the region of low saturation sorption of all three solutes is impeded by interactions among sorption sites; cobalt exhibits the greatest effect of interactions and strontium shows only a minimal effect. In the saturation region of 50% or more, sorption of cobalt is enhanced slightly by interactions among sorption sites whereas sorption of cesium and strontium appears to be independent of site interactions. 9 references, 4 figures, 2 tables.
Advection, diffusion, and delivery over a network
NASA Astrophysics Data System (ADS)
Heaton, Luke L. M.; López, Eduardo; Maini, Philip K.; Fricker, Mark D.; Jones, Nick S.
2012-08-01
Many biological, geophysical, and technological systems involve the transport of a resource over a network. In this paper, we present an efficient method for calculating the exact quantity of the resource in each part of an arbitrary network, where the resource is lost or delivered out of the network at a given rate, while being subject to advection and diffusion. The key conceptual step is to partition the resource into material that does or does not reach a node over a given time step. As an example application, we consider resource allocation within fungal networks, and analyze the spatial distribution of the resource that emerges as such networks grow over time. Fungal growth involves the expansion of fluid filled vessels, and such growth necessarily involves the movement of fluid. We develop a model of delivery in growing fungal networks, and find good empirical agreement between our model and experimental data gathered using radio-labeled tracers. Our results lead us to suggest that in foraging fungi, growth-induced mass flow is sufficient to account for long-distance transport, if the system is well insulated. We conclude that active transport mechanisms may only be required at the very end of the transport pathway, near the growing tips.
Advective Mechanisms in Tree Island Formation
NASA Astrophysics Data System (ADS)
Stothoff, S.
2002-05-01
Tree islands are important landscape features in the Florida Everglades. Tres islands are formed of peat deposited on the shallow limestone bedrock, and have been stressed as the system has changed in response to anthropogenic activities due to the sensitivity of organic soils to hydrologic cycles. The plume shape aligned with flow direction for typical tree islands is characteristic of advective transport, despite the rather low flow velocities in the system. Hypothesized mechanisms for the plume shape include sediment transport downstream from the head of the island (often anchored by a bedrock rise), or nutrient transport downstream allowing plants to produce more sediments in situ. Understanding mechanisms controlling tree island shape will aid in understanding the response of tree islands to hydrologic management. An integrated system of field, laboratory, and modeling studies is underway, with the first effort aimed at bounding the importance of the simpler sediment transport processes before tackling more-complex nutrient transport processes. The numerical model integrating the field and laboratory efforts is a 3D finite volume model considering water flow in the shallow groundwater/surface-water system together with sediment transport. The model can account for variable vegetative resistance through the flow column, including the important case where a dense mat forms at the surface. Model components specific for this system and associated data requirements are presented.
Vertical barriers with increased sorption capacities
Bradl, H.B.
1997-12-31
Vertical barriers are commonly used for the containment of contaminated areas. Due to the very small permeability of the barrier material which is usually in the order of magnitude of 10-10 m/s or less the advective contaminant transport can be more or less neglected. Nevertheless, there will always be a diffusive contaminant transport through the barrier which is caused by the concentration gradient. Investigations have been made to increase the sorption capacity of the barrier material by adding substances such as organoclays, zeolites, inorganic oxides and fly ashes. The contaminants taken into account where heavy metals (Pb) and for organic contaminants Toluole and Phenantrene. The paper presents results of model calculations and experiments. As a result, barrier materials can be designed {open_quotes}tailor-made{close_quotes} depending on the individual contaminant range of each site (e.g. landfills, gasworks etc.). The parameters relevant for construction such as rheological properties, compressive strength and permeability are not affected by the addition of the sorbents.
Haeri, Seyed Ammar
2016-08-15
In this study, biosorption based dispersive liquid liquid microextraction (Bio-DLLME) has been developed as a new method for the extraction of bisphenol A (BPA) from water samples. In this technique, the BPA is extracted into a stable cloudy phase. The colloidal phase is composed of micro-particles made from rhaminolipid biosurfactant and methanol, which dispersed in the water samples and facilitated the breakdown of analyte matrix bonds and provided high extraction yields. Rhaminolipid biosurfactants form a thin molecular interfacial film. This layer is formed between donor and recipient phase. This molecular layer, lowers the interfacial tension between immiscible phases (aqueous solution: colloidal particles) and allow dissimilar phases to mix and interact more easily. So the equilibrium state is achieved quickly and, therefore, the extraction time is very short. The attraction of the proposed method is that the extraction is fast, simple and can be done without toxic organic solvents. Also bioaggregates have several advantages such as higher environmental compatibility and biodegradability. Experimental parameters affecting the extraction efficiency were studied and optimized. Under the optimum conditions, relative recoveries of BPA were in the ranges of 98-103.3%. The calibration plot is linear in the range between 1 and 1000μgL(-1) (R(2)=0.998), and the relative standard deviation (RSD, for n=6) is 3.24%. PMID:27362996
NASA Astrophysics Data System (ADS)
Künzli, Pierre; Tsunematsu, Kae; Albuquerque, Paul; Falcone, Jean-Luc; Chopard, Bastien; Bonadonna, Costanza
2016-04-01
Volcanic ash transport and dispersal models typically describe particle motion via a turbulent velocity field. Particles are advected inside this field from the moment they leave the vent of the volcano until they deposit on the ground. Several techniques exist to simulate particles in an advection field such as finite difference Eulerian, Lagrangian-puff or pure Lagrangian techniques. In this paper, we present a new flexible simulation tool called TETRAS (TEphra TRAnsport Simulator) based on a hybrid Eulerian-Lagrangian model. This scheme offers the advantages of being numerically stable with no numerical diffusion and easily parallelizable. It also allows us to output particle atmospheric concentration or ground mass load at any given time. The model is validated using the advection-diffusion analytical equation. We also obtained a good agreement with field observations of the tephra deposit associated with the 2450 BP Pululagua (Ecuador) and the 1996 Ruapehu (New Zealand) eruptions. As this kind of model can lead to computationally intensive simulations, a parallelization on a distributed memory architecture was developed. A related performance model, taking into account load imbalance, is proposed and its accuracy tested.
Anisotropic Turbulent Advection of a Passive Vector Field: Effects of the Finite Correlation Time
NASA Astrophysics Data System (ADS)
Antonov, N. V.; Gulitskiy, N. M.
2016-02-01
The turbulent passive advection under the environment (velocity) field with finite correlation time is studied. Inertial-range asymptotic behavior of a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow, is investigated by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and prescribed pair correlation function. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to nontrivial fixed points of the RG equations and depend on the relation between the exponents in the energy energy spectrum ɛ ∝ k⊥1-ξ and the dispersion law ω ∝ k⊥2-η . The corresponding anomalous exponents are associated with the critical dimensions of tensor composite operators built solely of the passive vector field itself. In contrast to the well-known isotropic Kraichnan model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the dependence on the integral turbulence scale L has a logarithmic behavior: instead of power-like corrections to ordinary scaling, determined by naive (canonical) dimensions, the anomalies manifest themselves as polynomials of logarithms of L. Due to the presence of the anisotropy in the model, all multiloop diagrams are equal to zero, thus this result is exact.
Chaotic advection, diffusion, and reactions in open flows
Tel, Tamas; Karolyi, Gyoergy; Pentek, Aron; Scheuring, Istvan; Toroczkai, Zoltan; Grebogi, Celso; Kadtke, James
2000-03-01
We review and generalize recent results on advection of particles in open time-periodic hydrodynamical flows. First, the problem of passive advection is considered, and its fractal and chaotic nature is pointed out. Next, we study the effect of weak molecular diffusion or randomness of the flow. Finally, we investigate the influence of passive advection on chemical or biological activity superimposed on open flows. The nondiffusive approach is shown to carry some features of a weak diffusion, due to the finiteness of the reaction range or reaction velocity. (c) 2000 American Institute of Physics.
Burgers turbulence and passive random advection
NASA Astrophysics Data System (ADS)
Boldyrev, Stanislav Anatolievich
1999-10-01
, and the diffusivity is neglected. These considerations illustrate that even with simple statistics of the velocity field, the statistics of advected quantities are nontrivial due to nonlinear interactions of different spatial directions. The last Chapter 5 summarizes the results and discusses future directions of research.
Clay with Desiccation Cracks is an Advection Dominated Environment
NASA Astrophysics Data System (ADS)
Baram, S.; Kurtzman, D.; Sher, Y.; Ronen, Z.; Dahan, O.
2012-04-01
, indicating deep soil evaporation. Daily fluctuation of the air temperature in the desiccation cracks supported thermally induced air convection within the cracks void and could explain the deep soil salinization process. Combination of all the abovementioned observations demonstrated that the formation of desiccation cracks network in dispersive clay sediments generates a bulk advection dominated environment for both air and water flow, and that the reference to clay sediments as "hydrologically safe" should to be reconsidered.
Anomalous scaling of a scalar field advected by turbulence
Kraichnan, R.H.
1995-12-31
Recent work leading to deduction of anomalous scaling exponents for the inertial range of an advected passive field from the equations of motion is reviewed. Implications for other turbulence problems are discussed.
Magnetic helicity transport in the advective gauge family
Candelaresi, Simon; Brandenburg, Axel; Hubbard, Alexander; Mitra, Dhrubaditya
2011-01-15
Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e., when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.
Sorption isotherms and isosteric heats of sorption of Malaysian paddy.
Mousa, Wael; Ghazali, Farinazleen Mohamad; Jinap, S; Ghazali, Hasanah Mohd; Radu, Son
2014-10-01
Understanding the water sorption characteristics of cereal is extremely essential for optimizing the drying process and ensuring storage stability. Water relation of rough rice was studied at 20, 30, 40 and 50 °C over relative humidity (RH.) between 0.113 and 0.976 using the gravimetric technique. The isotherms displayed the general sigmoid, Type II pattern and exhibited the phenomenon of hysteresis where it was more pronounced at lower temperatures. The sorption characteristics were temperature dependence where the sorption capacity of the paddy increased as the temperature was decreased at fixed (RH). Among the models assessed for their ability to fit the sorption data, Oswin equation was the best followed by the third order polynomial, GAB, Smith, Chung-Pfost, and Henderson models. The monolayer moisture content was higher for desorption than adsorption and tend to decrease with the increase in temperature. Given the temperature dependence of the sorption isotherms the isosteric heats of sorption were calculated using Claussius-Clapeyron equation. The net isosteric heats decreased as the moisture content was increased and heats of desorption were greater than that of adsorption. PMID:25328208
Sorption of aniline derivatives on carbon fabric
NASA Astrophysics Data System (ADS)
Fazylova, G. F.; Valinurova, E. R.; Khamitov, E. M.; Kudasheva, F. Kh.
2015-06-01
The statistical adsorption of the chloro- and nitroderivatives of aniline on carbon fabric UVIS-AK was studied. The sorption isotherms of anilines on carbon fabric and coals were constructed at room temperature and a comparative analysis of their sorption activity was performed. The sorption isotherms were linearized in the coordinates of the Langmuir and Dubinin-Radushkevich equations. The main sorption parameters were calculated: the characteristic adsorption energies of anilines, changes in the Gibbs energies, monolayer capacities, and sorption equilibrium constants.
Advection around ventilated U-shaped burrows: A model study
NASA Astrophysics Data System (ADS)
Brand, Andreas; Lewandowski, JöRg; Hamann, Enrico; Nützmann, Gunnar
2013-05-01
Advective transport in the porous matrix of sediments surrounding burrows formed by fauna such as Chironomus plumosus has been generally neglected. A positron emission tomography study recently revealed that the pumping activity of the midge larvae can indeed induce fluid flow in the sediment. We present a numerical model study which explores the conditions at which advective transport in the sediment becomes relevant. A 0.15 m deep U-shaped burrow with a diameter of 0.002 m within the sediment was represented in a 3-D domain. Fluid flow in the burrow was calculated using the Navier-Stokes equation for incompressible laminar flow in the burrow, and flow in the sediment was described by Darcy's law. Nonreactive and reactive transport scenarios were simulated considering diffusion and advection. The pumping activity of the model larva results in considerable advective flow in the sediment at reasonable high permeabilities with flow velocities of up to 7.0 × 10-6 m s-1 close to the larva for a permeability of 3 × 10-12 m2. At permeabilities below 7 × 10-13 m2 advection is negligible compared to diffusion. Reactive transport simulations using first-order kinetics for oxygen revealed that advective flux into the sediment downstream of the pumping larva enhances sedimentary uptake, while the advective flux into the burrow upstream of the larvae inhibits diffusive sedimentary uptake. Despite the fact that both effects cancel each other with respect to total solute uptake, the advection-induced asymmetry in concentration distribution can lead to a heterogeneous solute and redox distribution in the sediment relevant to complex reaction networks.
Adaptive domain decomposition methods for advection-diffusion problems
Carlenzoli, C.; Quarteroni, A.
1995-12-31
Domain decomposition methods can perform poorly on advection-diffusion equations if diffusion is dominated by advection. Indeed, the hyperpolic part of the equations could affect the behavior of iterative schemes among subdomains slowing down dramatically their rate of convergence. Taking into account the direction of the characteristic lines we introduce suitable adaptive algorithms which are stable with respect to the magnitude of the convective field in the equations and very effective on bear boundary value problems.
Silica Precipitation and Lithium Sorption
Jay Renew
2015-09-20
This file contains silica precipitation and lithium sorption data from the project. The silica removal data is corrected from the previous submission. The previous submission did not take into account the limit of detection of the ICP-MS procedure.
A spatial SIS model in advective heterogeneous environments
NASA Astrophysics Data System (ADS)
Cui, Renhao; Lou, Yuan
2016-09-01
We study the effects of diffusion and advection for a susceptible-infected-susceptible epidemic reaction-diffusion model in heterogeneous environments. The definition of the basic reproduction number R0 is given. If R0 < 1, the unique disease-free equilibrium (DFE) is globally asymptotically stable. Asymptotic behaviors of R0 for advection rate and mobility of the infected individuals (denoted by dI) are established, and the existence of the endemic equilibrium when R0 > 1 is studied. The effects of diffusion and advection rates on the stability of the DFE are further investigated. Among other things, we find that if the habitat is a low-risk domain, there may exist one critical value for the advection rate, under which the DFE changes its stability at least twice as dI varies from zero to infinity, while the DFE is unstable for any dI when the advection rate is larger than the critical value. These results are in strong contrast with the case of no advection, where the DFE changes its stability at most once as dI varies from zero to infinity.
Can soil drying affect the sorption of pesticides in soil?
NASA Astrophysics Data System (ADS)
Chaplain, Véronique; Saint, Philippe; Mamy, Laure; Barriuso, Enrique
2010-05-01
The sorption of pesticides in soils mainly controls their further dispersion into the environment. Sorption is usually related to the physico-chemical properties of molecules but it also depends on the hydrophobic features of soils. However, the hydrophobicity of soils changes with wetting and drying cycles and this can be enhanced with climate change. The objective of this study was to measure by using controlled artificial soils the influence of the hydrophobic characteristic of soils on the retention of a model pesticide. Artificial soils consisted in silica particles covered by synthetic cationic polymers. Polymers were characterized by the molar ratio of monomers bearing an alkyl chain of 12C. Two polymers were used, with 20 and 80 % ratios, and the same degree of polymerization. In addition, porous and non-porous particles were used to study the accessibility notion and to measure the influence of diffusion on pesticide sorption kinetics. Lindane was chosen as model molecule because its adsorption is supposed mainly due to hydrophobic interactions. Results on polymers adsorption on silica showed that it was governed by electrostatic interactions, without any dependency of the hydrophobic ratio. Polymers covered the entire surface of porous particles. Kinetic measurements showed that lindane sorption was slowed in porous particles due to the molecular diffusion inside the microporosity. The adsorption of lindane on covered silica particles corresponded to a partition mechanism described by linear isotherms. The slope was determined by the hydrophobic ratio of polymers: the sorption of lindane was highest in the most hydrophobic artificial soil. As a result, modification in soil hydrophobicity, that can happen with climate change, might affect the sorption and the fate of pesticides. However additional experiments are needed to confirm these first results. Such artificial soils should be used as reference materials to compare the reactivity of pesticides, to
Linearity of iodine sorption and sorption capacities for seven soils
Sheppard, M.I.; Hawkins, J.L.; Smith, P.A.
1996-11-01
Iodine, a soluble and prevalent element in spent nuclear fuel and a pivotal element in the assessment of Canada`s nuclear fuel waste disposal option, sorbs to soils rich in organics and hydrous oxides. Biotic factors, such as microbes, enzymes and plant exudates, have been implicated in the retention of I to soils. Anion exchange of I{sup {minus}} or IO{sub 3}{sup {minus}} and chemical or biological oxidation to I{sub 2} followed by reactions with the soil organic matter are possible retention mechanisms. We have carried out sorption and desorption studies across a wide range of soil solution concentrations (10{sup {minus}7} to 10{sup 5} mg I/L, 10{sup {minus}12} to 1 M) for seven soils typical of upland and lowland soils of the Canadian Precambrian Shield. Soil solid-liquid partition values (K{sub d}), required for impact assessments, varied from 60 to 1800 L/kg and were significantly correlated with extractable Al oxide content, and background I and organic matter content. Freundlich isotherm fits show that sorption of I across our intentionally large concentration range is nonlinear; however, sorption of I across our intentionally large concentration range is nonlinear; however, sorption of I at environmental concentrations (<0.1 mg I/L soil solution) is linear and can be described by the K{sub d} model. Sorption of I was not related to peroxidase enzyme activity. Desorption percentages were small implying sorption was not easily reversed, even with a strong electrolyte, KNO{sub 3}. Desorption results and simple correlations of I sorption to soil properties suggest that the oxidation of I{sup {minus}} to I{sub 2} and complexation to organic functional groups or oxides are the major processes for I retention in Shield soils. 52 refs., 2 figs., 2 tabs.
NASA Astrophysics Data System (ADS)
Ancey, C.; Bohorquez, P.; Heyman, J.
2015-12-01
The advection-diffusion equation is one of the most widespread equations in physics. It arises quite often in the context of sediment transport, e.g., for describing time and space variations in the particle activity (the solid volume of particles in motion per unit streambed area). Phenomenological laws are usually sufficient to derive this equation and interpret its terms. Stochastic models can also be used to derive it, with the significant advantage that they provide information on the statistical properties of particle activity. These models are quite useful when sediment transport exhibits large fluctuations (typically at low transport rates), making the measurement of mean values difficult. Among these stochastic models, the most common approach consists of random walk models. For instance, they have been used to model the random displacement of tracers in rivers. Here we explore an alternative approach, which involves monitoring the evolution of the number of particles moving within an array of cells of finite length. Birth-death Markov processes are well suited to this objective. While the topic has been explored in detail for diffusion-reaction systems, the treatment of advection has received no attention. We therefore look into the possibility of deriving the advection-diffusion equation (with a source term) within the framework of birth-death Markov processes. We show that in the continuum limit (when the cell size becomes vanishingly small), we can derive an advection-diffusion equation for particle activity. Yet while this derivation is formally valid in the continuum limit, it runs into difficulty in practical applications involving cells or meshes of finite length. Indeed, within our stochastic framework, particle advection produces nonlocal effects, which are more or less significant depending on the cell size and particle velocity. Albeit nonlocal, these effects look like (local) diffusion and add to the intrinsic particle diffusion (dispersal due
Dynamic micro-mapping of CO2 sorption in coal
Radlinski, Andrzej Pawell; Melnichenko, Yuri B; Cheng, Gang; Mastalerz, Maria
2009-01-01
We have applied X-ray and neutron small-angle scattering techniques (SAXS, SANS and USANS) to study the interaction between fluids and porous media in the particular case of sub- and super-critical CO2 sorption in coal. These techniques are demonstrated to give unique, pore-size-specific insights into the kinetics of CO2 sorption in a wide range of coal pores (nano to meso), and to provide data that may be used to determine the density of the sorbed CO2. We observed densification of the adsorbed CO2 by a factor up to five compared to the free fluid at the same (p,T) conditions. Our results indicate that details of CO2 sorption into coal pores differ greatly between different coals and depend on the amount of mineral matter dispersed in the coal matrix: a purely organic matrix absorbs more CO2 per unit volume than one contaminated with mineral matter, but mineral matter markedly accelerates the sorption kinetics.
Iodine sorption by highly weathered soils
NASA Astrophysics Data System (ADS)
Bruns, W. A.; Bruns, M. R.; Kaplan, D. I.; Murphy, T.; Roberts, K.; Thompson, A.; Seaman, J. C.
2009-12-01
Radioactive iodine (I) isotopes such as I-129 are widely distributed throughout Department of Energy (DOE) and other nuclear sites as a result of past enrichment activities. While natural-abundance I serves an important role in thyroid function in human beings, anthropogenic I can cause thyroid cancer. Predominant I species are anionic and generally assumed to move conservatively through typical negatively-charged subsurface material with the mobile water phase. However, many of the 129I contaminated soils and subsurface materials on the DOE Savannah River Site (SRS) near Aiken, SC, contain significant positive charge, which could retain I-129. A series of batch and column experiments were performed in order to determine sorption of three different I species (iodide, iodate, and an organic I compound, 4-iodaniline) over a range of aqueous concentrations (10^-4 to 10^-8 M) onto surface and subsurface soils from SRS. Total aqueous I in solution after equilibration was measured using ICP-MS and used to calculate the partitioning coefficient (Kd). Kd values for all three species were much higher at lower dissolved I concentrations (i.e., non-linear partitioning), and correlated with CDB-extractable Fe and Al and well as soil organic matter content. For the species and concentrations that exhibited high Kd values (>5), we conducted column experiments to evaluate partitioning under dynamic advective systems, using tritium for comparison as a conservative hydrologic tracer. We tracked pH, EC, and hydraulic head in real-time, and the effluent was collected for I analysis as described above. The results indicate a complex pattern of I movement in the natural environment affected by isotope concentration, speciation, and soil characteristics.
Blum, A.E.; Eberl, D.D.
2004-01-01
A new method has been developed for quantifying smectite abundance by sorbing polyvinylpyrrolidone (PVP) on smectite particles dispersed in aqueous solution. The sorption density of PVP-55K on a wide range of smectites, illites and kaolinites is ???0.99 mg/m2, which corresponds to ???0.72 g of PVP-55K per gram of montmorillonite. Polyvinylpyrrolidone sorption on smectites is independent of layer charge and solution pH. PVP sorption on Si02, Fe 2O3 and ZnO normalized to the BET surface area is similar to the sorption densities on smectites. ??-Al 2O3, amorphous Al(OH)3 and gibbsite have no PVP sorption over a wide range of pH, and sorption of PVP by organics is minimal. The insensitivity of PVP sorption densities to mineral layer charge, solution pH and mineral surface charge indicates that PVP sorption is not localized at charged sites, but is controlled by more broadly distributed sorption mechanisms such as Van der Waals' interactions and/or hydrogen bonding. Smectites have very large surface areas when dispersed as single unit-cell-thick particles (???725 m2/g) and usually dominate the total surface areas of natural samples in which smectites are present. In this case, smectite abundance is directly proportional to PVP sorption. In some cases, however, the accurate quantification of smectite abundance by PVP sorption may require minor corrections for PVP uptake by other phases, principally illite and kaolinite. Quantitative XRD can be combined with PVP uptake measurements to uniquely determine the smectite concentration in such sample. ?? 2004, The Clay Minerals Society.
Modeling of sorption characteristics of backfill materials
Chitra, S.; Sasidhar, P.; Lal, K.B.; Ahmed, J.
1998-06-01
Sorption data analysis was carried out using the Freundlich, Langmuir, and Modified Freundlich isotherms for the uptake of sodium and potassium in an initial concentration range of 10--100 mg/L on backfill materials, viz., bentonite, vermiculite, and soil samples. The soil samples were collected from a shallow land disposal facility at Kalpakkam. The Freundlich isotherm equation is validated as a preferred general mathematical tool for representing the sorption of K{sup +} by all the selected backfill materials. The Modified Freundlich isotherm equation is validated as a preferred mathematical tool for representing the sorption of Na{sup +} by the soil samples. Since a negative sorption was observed for the uptake of Na{sup +} by commercial clay minerals (vermiculite and bentonite clay in the laboratory experiments), sorption analysis could not be carried out using the above-mentioned isotherm equations. Hill plots of the sorption data suggest that in the region of low saturation (10--40 mg/L), sorption of K{sup +} by vermiculite is impeded by interaction among sorption sites. In the region of higher saturation (60--100 mg/L), sorption of K{sup +} by all three backfill materials is enhanced by interaction among sorption sites. The Hill plot of the sorption data for Na{sup +} by soil suggests that irrespective of Na{sup +} concentration, sorption of Na{sup +} at one exchange size enhances sorption at other exchange sites.
Sorption of heavy metals by the soil fungi 'Aspergillus niger' and Mucor rouxii
Mullen, M.D.; Wolf, D.C.; Beveridge, T.J.; Bailey, G.W.
1992-01-01
Sorption of the nitrate salts of cadmium(II), copper(II), lanthanum(III) and silver(I) by two fungi, Aspergillus niger and Mucor rouxii, was evaluated using Freundlich adsorption isotherms and energy dispersive X-ray electron microscopy. The linearized Freundlich isotherm described the metal sorption data well for metal concentrations of 5 microM-1 mM metal. Differences in metal binding were observed among metals, as well as between fungal species. Calculated Freundlich K values indicated that metal binding decreased in the order La(3+) > or = Ag(+) > Cu(2+) > Cd(2+). However, sorption of Ag(+) was greater than that of La(3+) from solutions of 0.1 and 1 mM metal and likely due to precipitation at the cell wall surface. At the 1 mM initial concentration, there were no significant differences between the two fungi in metal sorption, except for Ag(+) binding. At the 5 microM concentration, there was no difference between the fungi in their sorption capacities for the four metals. Electron microscopy-energy dispersive X-ray analysis indicated that silver precipitated onto cells as colloidal silver. The results indicate that Freundlich isotherms may be useful for describing short-term metal sorption by fungal biomass and for comparison with other soil constituents in standardized systems. (Copyright (c) 1992 Pergamon Press plc.)
Advecting Procedural Textures for 2D Flow Animation
NASA Technical Reports Server (NTRS)
Kao, David; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
This paper proposes the use of specially generated 3D procedural textures for visualizing steady state 2D flow fields. We use the flow field to advect and animate the texture over time. However, using standard texture advection techniques and arbitrary textures will introduce some undesirable effects such as: (a) expanding texture from a critical source point, (b) streaking pattern from the boundary of the flowfield, (c) crowding of advected textures near an attracting spiral or sink, and (d) absent or lack of textures in some regions of the flow. This paper proposes a number of strategies to solve these problems. We demonstrate how the technique works using both synthetic data and computational fluid dynamics data.
Concentration polarization, surface currents, and bulk advection in a microchannel
NASA Astrophysics Data System (ADS)
Nielsen, Christoffer P.; Bruus, Henrik
2014-10-01
We present a comprehensive analysis of salt transport and overlimiting currents in a microchannel during concentration polarization. We have carried out full numerical simulations of the coupled Poisson-Nernst-Planck-Stokes problem governing the transport and rationalized the behavior of the system. A remarkable outcome of the investigations is the discovery of strong couplings between bulk advection and the surface current; without a surface current, bulk advection is strongly suppressed. The numerical simulations are supplemented by analytical models valid in the long channel limit as well as in the limit of negligible surface charge. By including the effects of diffusion and advection in the diffuse part of the electric double layers, we extend a recently published analytical model of overlimiting current due to surface conduction.
Fast multigrid solution of the advection problem with closed characteristics
Yavneh, I.; Venner, C.H.; Brandt, A.
1996-12-31
The numerical solution of the advection-diffusion problem in the inviscid limit with closed characteristics is studied as a prelude to an efficient high Reynolds-number flow solver. It is demonstrated by a heuristic analysis and numerical calculations that using upstream discretization with downstream relaxation-ordering and appropriate residual weighting in a simple multigrid V cycle produces an efficient solution process. We also derive upstream finite-difference approximations to the advection operator, whose truncation terms approximate {open_quotes}physical{close_quotes} (Laplacian) viscosity, thus avoiding spurious solutions to the homogeneous problem when the artificial diffusivity dominates the physical viscosity.
Hydrodynamic dispersion within porous biofilms.
Davit, Y; Byrne, H; Osborne, J; Pitt-Francis, J; Gavaghan, D; Quintard, M
2013-01-01
Many microorganisms live within surface-associated consortia, termed biofilms, that can form intricate porous structures interspersed with a network of fluid channels. In such systems, transport phenomena, including flow and advection, regulate various aspects of cell behavior by controlling nutrient supply, evacuation of waste products, and permeation of antimicrobial agents. This study presents multiscale analysis of solute transport in these porous biofilms. We start our analysis with a channel-scale description of mass transport and use the method of volume averaging to derive a set of homogenized equations at the biofilm-scale in the case where the width of the channels is significantly smaller than the thickness of the biofilm. We show that solute transport may be described via two coupled partial differential equations or telegrapher's equations for the averaged concentrations. These models are particularly relevant for chemicals, such as some antimicrobial agents, that penetrate cell clusters very slowly. In most cases, especially for nutrients, solute penetration is faster, and transport can be described via an advection-dispersion equation. In this simpler case, the effective diffusion is characterized by a second-order tensor whose components depend on (1) the topology of the channels' network; (2) the solute's diffusion coefficients in the fluid and the cell clusters; (3) hydrodynamic dispersion effects; and (4) an additional dispersion term intrinsic to the two-phase configuration. Although solute transport in biofilms is commonly thought to be diffusion dominated, this analysis shows that hydrodynamic dispersion effects may significantly contribute to transport. PMID:23410370
Tidal Mixing and Buoyancy Advection: Joint Influences on Lobster Distribution in Coastal Maine
NASA Astrophysics Data System (ADS)
Brooks, D. A.
2004-12-01
The Eastern Maine Coastal Current (EMCC) flows southwestward from the mouth of the Bay of Fundy to Penobscot Bay on the central Maine coast. Maximum non-tidal surface speeds reach 20-30 cm/s about 20 km offshore during April-May when the outflow from the Saint John River is strongest. Vigorous tides cause strong vertical and horizontal mixing, so that dispersal of neutral particles is influenced both by advection and tidal mixing. To survive, planktonic lobsters carried southwestward in the surface flow must settle to a nearshore cobble substrate. Larvae hatched near the Bay of Fundy can be advected to the central coast in 2-3 weeks, roughly the time needed to reach settlement stage. Over the same period, transverse tidal mixing is sufficient to raise nearshore larval concentrations to about half that offshore in the axis of the EMCC. Both processes may be necessary to explain the observed lobster distribution, which exhibits a distinct maximum in the central coastal region. The seasonal development of the EMCC is also influenced by winds and the larger circulation of the Gulf of Maine. This work is part of a multidisciplinary synthesis study funded by the NOAA Coastal Ocean Program.
Theory of advection-driven long range biotic transport
Technology Transfer Automated Retrieval System (TEKTRAN)
We propose a simple mechanistic model to examine the effects of advective flow on the spread of fungal diseases spread by wind-blown spores. The model is defined by a set of two coupled non-linear partial differential equations for spore densities. One equation describes the long-distance advectiv...
Topology preserving advection of implicit interfaces on Cartesian grids
NASA Astrophysics Data System (ADS)
Qin, Zhipeng; Delaney, Keegan; Riaz, Amir; Balaras, Elias
2015-06-01
Accurate representation of implicit interface topology is important for the numerical computation of two phase flow on Cartesian grids. A new method is proposed for the construction of signed distance function by geometrically projecting interface topology onto the Cartesian grid using a multi-level projection framework. The method involves a stepwise improvement in the approximation to the signed distance function based on pointwise, piecewise and locally smooth reconstructions of the interface. We show that this approach provides accurate representation of the projected interface and its topology on the Cartesian grid, including the distance from the interface and the interface normal and curvature. The projected interface can be in the form of either a connected set of marker particles that evolve with Lagrangian advection, or a discrete set of points associated with an implicit interface that evolves with the advection of a scalar function. The signed distance function obtained with geometric projection is independent of the details of the scaler field, in contrast to the conventional approach where advection and reinitialization cannot be decoupled. As a result, errors introduced by reinitialization do not amplify advection errors, which leads to substantial improvement in both volume conservation and topology representation.
Consistency problem with tracer advection in the Atmospheric Model GAMIL
NASA Astrophysics Data System (ADS)
Zhang, Kai; Wan, Hui; Wang, Bin; Zhang, Meigen
2008-03-01
The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG (GAMIL). Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polar-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense.
ADVECTION INFLUENCES ON EVAPOTRANSPIRATION OF ALFALFA IN A SEMIARID ENVIRONMENT
Technology Transfer Automated Retrieval System (TEKTRAN)
Advective enhancement of crop evapotranspiration (ET) occurs when drier, hotter air is transported into the crop by wind and can be an important factor in the water balance of irrigated crops in a semiarid climate. Thirteen days of moderate to extremely high ET rates of irrigated alfalfa (Medicago ...
The 'optimum' upwind advection on a triangular mesh
NASA Technical Reports Server (NTRS)
Roe, P. L.
1990-01-01
For advection schemes based on fluctuation splitting, a design criterion of optimizing the time step leads to linear schemes that coincide with those designed for least truncation error. A further stage of optimizing the time step using a nonlinear positivity criterion, leads to considerable further gains in resolution.
Tidal variability of lateral advection in a coastal plain estuary
NASA Astrophysics Data System (ADS)
Basdurak, N. B.; Valle-Levinson, A.
2013-07-01
Tidal variability of lateral advection of momentum (vuy, where u and v are along-estuary and lateral flows, respectively, and the subindex indicates differentiation with respect to the cross-estuary direction) was investigated in a coastal plain estuary with observations at Hampton Roads, which is the transition between the James River and Chesapeake Bay. Towed current velocity profiles and hydrographic profiles were captured during 9 expeditions in 2004 and 2005, to determine the intratidal and spatial changes in lateral advection of momentum and its contribution to along-channel flow. Curvature effects and lateral density gradients were important in driving lateral circulation and in modifying intratidal lateral advection of momentum. Lateral advection had the same order of magnitude as the baroclinic pressure gradient. Its contribution to the along-channel momentum balance was greatest during or just after peak flood and weakest at the end of ebb. During peak flood and peak ebb, the spatial distribution of vuy was seaward at the southern (left) side near surface and at the northern side (right) near bed (looking up-estuary), and landward in the rest of the channel. During slack periods the vuy structures were mostly landward. Observations were in good agreement with analytical model results during peak ebb and flood, but inconsistent during slack periods. The discrepancies between model results and field measurements can be attributed to bathymetry-density gradient interactions, which enhanced ebb-to-flood asymmetries in the along-channel and lateral flow.
BACTERIAL SORPTION OF HEAVY METALS
Four bacteria, Bacillus cereus, B. subtilis, Escherichia coli, and Pseudomonas aeruginosa, were examined for the ability to remove Ag+, Cd2+, Cu2+, and La3+ from solution by batch equilibration methods. d and Cu sorption over the concentration range 0.001 to 1 mM was described by...
Flaherty, L.M.
1989-01-01
This book contains papers presented at a symposium of the American Society for Testing and Materials. The topics covered include: The effect of elastomers on the efficiency of oil spill dispersants; planning for dispersant use; field experience with dispersants for oil spills on land; and measurements on natural dispersion.
NASA Astrophysics Data System (ADS)
Liao, Zijie; Lemke, Dennis; Osenbrück, Karsten; Cirpka, Olaf A.
2013-06-01
Performing stream-tracer experiments is an accepted technique to assess transport characteristics of streams undergoing hyporheic exchange. Recently, combining conservative and reactive tracers, in which the latter presumably undergoes degradation exclusively within the hyporheic zone, has been suggested to study in-stream transport, hyporheic exchange, and the metabolic activity of the hyporheic zone. The combined quantitative analysis to adequately describe such tests, however, has been missing. In this paper, we present mathematical methods to jointly analyze breakthrough curves of a conservative tracer (fluorescein), a linearly degrading tracer (resazurin), and its daughter compound (resorufin), which are synchronously introduced into the stream as pulses. In-stream transport is described by the one-dimensional advection-dispersion equation, amended with a convolution term to account for transient storage within the hyporheic zone over a distribution of travel times, transformation of the reactive tracer in the hyporheic zone, and two-site sorption of the parent and daughter compounds therein. We use a shape-free approach of describing the hyporheic travel-time distribution, overcoming the difficulty of identifying the best functional parameterization for transient storage. We discuss how this model can be fitted to the breakthrough curves of all three compounds and demonstrate the method by an application to a tracer test in the third-order stream Goldersbach in Southern Germany. The entire river water passes once through the hyporheic zone over a travel distance of about 200 m with mean hyporheic residence times ranging between 16 and 23 min. We also observed a secondary peak in the transfer functions at about 1 h indicating a second hyporheic flow path. We could jointly fit the breakthrough curves of all compounds in three monitoring stations and evaluated the parameter uncertainty of the individual and joint fits by a method based on conditional
SORPTION OF HEAVY METALS BY THE SOIL FUNGI ASPERGILLUS NIGER AND MUCOR ROUXII
Sorption of the nitrate salts of cadmium(II), copper (II), lanthanum(III) and silver (I) by two fungi, Aspergillus niger and Mucor rouxii, was evaluated using Fruendlich adsorption isotherms and energy dispersive X-ray electron microscopy. The linearized Freundlich isotherm descr...
NASA Astrophysics Data System (ADS)
Eslami, Parastou; Seo, Jung-Hee; Rahsepar, Amirali; George, Richard; Lardo, Albert; Mittal, Rajat
2014-11-01
Coronary computed tomography angiography (CTA) is a promising tool for assessment of coronary stenosis and plaque burden. Recent studies have shown the presence of axial contrast concentration gradients in obstructed arteries, but the mechanism responsible for this phenomenon is not well understood. We use computational fluid dynamics to study intracoronary contrast dispersion and the correlation of concentration gradients with intracoronary blood flow and stenotic severity. Data from our CFD patient-specific simulations reveals that contrast dispersions are generated by intracoronary advection effects, and therefore, encode the coronary flow velocity. This novel method- Transluminal Attenuation Flow Encoding (TAFE) - is used to estimate the flowrate in phantom studies as well as preclinical experiments. Our results indicate a strong correlation between the values estimated from TAFE and the values measured in these experiments. The flow physics of contrast dispersion associated with TAFE will be discussed. This work is funded by grants from Coulter Foundation and Maryland Innovation Initiative. The authors have pending patents in this technology and RM and ACL have other financial interests associated with TAFE.
NASA Astrophysics Data System (ADS)
Breus, I. P.; Mishchenko, A. A.; Shinkarev, A. A.; Neklyudov, S. A.; Breus, V. A.
2014-12-01
The diverse effect of the organic matter content on the sorption of vapors of aromatic and aliphatic hydrocarbons in soils under low moisture (<10.5%) has been revealed in sorption experiments using profile samples from two virgin heavy loamy dark gray forest soils characterized by relatively stable contents of finely dispersed mineral components. The decrease of the hydrocarbon sorption with increasing the content of organic matter under dry conditions (in the moisture range from 0 to 5-6%) indicates its lower sorption activity than that of the clay components and the blocking of the sorption sites on soil minerals by organic matter. At moisture contents above 5-6%, the effect of the soil composition on the sorption activity changes radically: it increases with increasing the content of organic matter. This is due to the inversion of the ratio between the activities of the soil components because of the hydrophilization of the surface of the mineral soil component. As a result, the sorption of water on the minerals reduces the mineral sorption activity to hydrocarbons to a lower level than the activity of organic matter. The maximum manifestation of the revealed blocking effect has been observed for the low-humus soils and this effect decreased with the accumulation of soil organic matter.
Sorption Characteristics of Organic Powder Sorption Material in Fluidized Bed with a Cooling Pipe
NASA Astrophysics Data System (ADS)
Horibe, Akihiko; Husain, Syahrul; Inaba, Hideo; Haruki, Naoto
The dynamic sorption characteristics of organic sorbent materials have been studied by using fluidized bed with a cooling pipe. The organic powder type sorbent made from a bridged complex of sodium polyacrylate which is one of the sorption polymers is adopted in this study. Sorption rate of water vapor and the variation of temperature in the sorbent bed with time were measured under various conditions. As a result, sorption ratio increases and the completion time for the sorption process decreases by using a cooling pipe. Furthermore, the non-dimensional correlation equations were obtained for water vapor mass transfer under sorption process in terms of relevant non-dimensional parameters.
Mixed refrigerant Joule-Thomson sorption cryocoolers
NASA Astrophysics Data System (ADS)
Tzabar, Nir; Grossman, Gershon
2014-01-01
Joule-Thomson (JT) sorption cryocooling is the most mature technology for cooling from a normal Room-Temperature (RT) down to temperatures below 100 K in the absence of moving parts. Therefore, high reliability and no vibrations are attainable, in comparison with other cryocoolers. Cooling to 80 - 100 K with JT cryocoolers is often implemented with pure nitrogen. Alternatively, mixed refrigerants have been suggested for reducing the operating pressures to enable closed cycle cryocooling. There is a variety of publications describing nitrogen sorption cryocoolers with different configurations of sorption compressors. In the present research we suggest a novel sorption JT cryocooler that operates with a mixed refrigerant. Merging of sorption cryocooling and a mixed refrigerant enables the use of a simple, single stage compressor for cooling to 80 - 100 K, lower operating temperatures of the sorption cycle, and thus - reduced power consumption. In previous studies we have analyzed sorption compressors for mixed gases and mixed refrigerants for JT cryocoolers, separately. In this paper the option of mixed refrigerant sorption JT cryocoolers is explored. The considerations for developing mixed refrigerants to be driven by sorption compressors and to be utilized with JT cryocoolers are provided. It appears that, unlike with pure nitrogen, mixed refrigerants can be suitable for JT cryocooling with a single stage sorption compressor.
Sorption of perfluoroalkyl substances in sewage sludge.
Milinovic, Jelena; Lacorte, Silvia; Rigol, Anna; Vidal, Miquel
2016-05-01
The sorption behaviour of three perfluoroalkyl substances (PFASs) (perfluorooctanesulfonic acid (PFOS), perfluorooctanoic acid (PFOA) and perfluorobutanesulfonic acid (PFBS)) was studied in sewage sludge samples. Sorption isotherms were obtained by varying initial concentrations of PFOS, PFOA and PFBS. The maximum values of the sorption solid-liquid distribution coefficients (Kd,max) varied by almost two orders of magnitude among the target PFASs: 140-281 mL g(-1) for PFOS, 30-54 mL g(-1) for PFOA and 9-18 mL g(-1) for PFBS. Freundlich and linear fittings were appropriate for describing the sorption behaviour of PFASs in the sludge samples, and the derived KF and Kd,linear parameters correlated well. The hydrophobicity of the PFASs was the key parameter that influenced their sorption in sewage sludge. Sorption parameters and log(KOW) were correlated, and for PFOS (the most hydrophobic compound), pH and Ca + Mg status of the sludge controlled the variation in the sorption parameter values. Sorption reversibility was also tested from desorption isotherms, which were also linear. Desorption parameters were systematically higher than the corresponding sorption parameters (up to sixfold higher), thus indicating a significant degree of irreversible sorption, which decreased in the sequence PFOS > PFOA > PFBS. PMID:26780052
Technology Transfer Automated Retrieval System (TEKTRAN)
It has been reported that this model cannot take into account several important features of solute movement through soil. Recently, a new model has been suggested that results in a solute transport equation with fractional spatial derivatives, or FADE. We have assembled a database on published solu...
Sorption of non-polar organic compounds by micro-sized plastic particles in aqueous solution.
Hüffer, Thorsten; Hofmann, Thilo
2016-07-01
The presence of microscale polymer particles (i.e., microplastics) in the environment has become a major concern in recent years. Sorption of organic compounds by microplastics may affect the phase distribution within both sediments and aqueous phases. To investigate this process, isotherms were determined for the sorption of seven aliphatic and aromatic organic probe sorbates by four polymers with different physico-chemical properties. Sorption increased in the order polyamide < polyethylene < polyvinylchloride < polystyrene. This order does not reflect the particle sizes of the investigated microplastics within the aqueous dispersions, indicating the influence of additional factors (e.g., π-π-interactions) on the sorption of aromatic compounds by polystyrene. Linear isotherms by polyethylene suggested that sorbate uptake was due to absorption into the bulk polymer. In contrast, non-linear isotherms for sorption by PS, PA, and PVC suggest a predominance of adsorption onto the polymer surface, which is supported by the best fit of these isotherms using the Polanyi-Manes model. A strong relationship between the sorption coefficients of the microplastics and the hydrophobicity of the sorbates suggests that hydrophobic interactions are of major importance. PMID:27086075
NASA Technical Reports Server (NTRS)
Wang, Xiao-Yen; Chow, Chuen-Yen; Chang, Sin-Chung
1999-01-01
Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.
Aerosol particles and the formation of advection fog
NASA Technical Reports Server (NTRS)
Hung, R. J.; Liaw, G. S.; Vaughan, O. H., Jr.
1979-01-01
A study of numerical simulation of the effects of concentration, particle size, mass of nuclei, and chemical composition on the dynamics of warm fog formation, particularly the formation of advection fog, is presented. This formation is associated with the aerosol particle characteristics, and both macrophysical and microphysical processes are considered. In the macrophysical model, the evolution of wind components, water vapor content, liquid water content, and potential temperature under the influences of vertical turbulent diffusion, turbulent momentum, and turbulent energy transfers are taken into account. In the microphysical model, the supersaturation effect is incorporated with the surface tension and hygroscopic material solution. It is shown that the aerosol particles with the higher number density, larger size nuclei, the heavier nuclei mass, and the higher ratio of the Van't Hoff factor to the molecular weight favor the formation of the lower visibility advection fogs with stronger vertical energy transfer during the nucleation and condensation time period.
Spiral defect chaos in an advection-reaction-diffusion system
NASA Astrophysics Data System (ADS)
Affan, H.; Friedrich, R.
2014-06-01
This paper comprises numerical and theoretical studies of spatiotemporal patterns in advection-reaction-diffusion systems in which the chemical species interact with the hydrodynamic fluid. Due to the interplay between the two, we obtained the spiral defect chaos in the activator-inhibitor-type model. We formulated the generalized Swift-Hohenberg-type model for this system. Then the evolution of fractal boundaries due to the effect of the strong nonlinearity at the interface of the two chemical species is studied numerically. The purpose of the present paper is to point out that spiral defect chaos, observed in model equations of the extended Swift-Hohenberg equation for low Prandtl number convection, may actually be obtained also in certain advection-reaction-diffusion systems.
Phase Segregation of Passive Advective Particles in an Active Medium.
Das, Amit; Polley, Anirban; Rao, Madan
2016-02-12
Localized contractile configurations or asters spontaneously appear and disappear as emergent structures in the collective stochastic dynamics of active polar actomyosin filaments. Passive particles which (un)bind to the active filaments get advected into the asters, forming transient clusters. We study the phase segregation of such passive advective scalars in a medium of dynamic asters, as a function of the aster density and the ratio of the rates of aster remodeling to particle diffusion. The dynamics of coarsening shows a violation of Porod behavior; the growing domains have diffuse interfaces and low interfacial tension. The phase-segregated steady state shows strong macroscopic fluctuations characterized by multiscaling and intermittency, signifying rapid reorganization of macroscopic structures. We expect these unique nonequilibrium features to manifest in the actin-dependent molecular clustering at the cell surface. PMID:26919022
Phase Segregation of Passive Advective Particles in an Active Medium
NASA Astrophysics Data System (ADS)
Das, Amit; Polley, Anirban; Rao, Madan
2016-02-01
Localized contractile configurations or asters spontaneously appear and disappear as emergent structures in the collective stochastic dynamics of active polar actomyosin filaments. Passive particles which (un)bind to the active filaments get advected into the asters, forming transient clusters. We study the phase segregation of such passive advective scalars in a medium of dynamic asters, as a function of the aster density and the ratio of the rates of aster remodeling to particle diffusion. The dynamics of coarsening shows a violation of Porod behavior; the growing domains have diffuse interfaces and low interfacial tension. The phase-segregated steady state shows strong macroscopic fluctuations characterized by multiscaling and intermittency, signifying rapid reorganization of macroscopic structures. We expect these unique nonequilibrium features to manifest in the actin-dependent molecular clustering at the cell surface.
Indium Sorption to Iron Oxides
NASA Astrophysics Data System (ADS)
White, S. J.; Sacco, S. A.; Hemond, H.; Hussain, F. A.; Runkel, R. L.; Walton-Day, K. E.; Kimball, B. A.; Shine, J. P.
2014-12-01
Indium is an increasingly important metal in semiconductors and electronics, and its use is growing rapidly as a semiconductive coating (as indium tin oxide) for liquid crystal displays (LCDs) and flat panel displays. It also has uses in important energy technologies such as light emitting diodes (LEDs) and photovoltaic cells. Despite its rapid increase in use, very little is known about the environmental behavior of indium, and concerns are being raised over the potential health effects of this emerging metal contaminant. One source of indium to the environment is acid mine drainage from the mining of lead, zinc, and copper sulfides. In our previous studies of a stream in Colorado influenced by acid mine drainage from lead and zinc mining activities, indium concentrations were found to be 10,000 times those found in uncontaminated rivers. However, the speciation and mobility of indium could not be reliably modeled because sorption constants to environmental sorbents have not been determined. In this study, we generate sorption constants for indium to ferrihydrite in the laboratory over a range of pHs, sorbent to sorbate ratios, and ionic strengths. Ferrihydrite is one of the most important sorbents in natural systems, and sorption to amorphous iron oxides such as ferrihydrite is thought to be one of the main removal mechanisms of metals from the dissolved phase in aqueous environments. Because of its relatively low solubility, we also find that indium hydroxide precipitation can dominate indium's partitioning at micromolar concentrations of indium. This precipitation may be important in describing indium's behavior in our study stream in Colorado, where modeling sorption to iron-oxides does not explain the complete removal of indium from the dissolved phase when the pH of the system is artificially raised to above 8. This study contributes much-needed data about indium's aqueous behavior, in order to better understand its fate, transport, and impacts in the
Lattice Boltzmann method for the fractional advection-diffusion equation
NASA Astrophysics Data System (ADS)
Zhou, J. G.; Haygarth, P. M.; Withers, P. J. A.; Macleod, C. J. A.; Falloon, P. D.; Beven, K. J.; Ockenden, M. C.; Forber, K. J.; Hollaway, M. J.; Evans, R.; Collins, A. L.; Hiscock, K. M.; Wearing, C.; Kahana, R.; Villamizar Velez, M. L.
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β , the fractional order α , and the single relaxation time τ , the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering.
Advective and diffusive cosmic ray transport in galactic haloes
NASA Astrophysics Data System (ADS)
Heesen, Volker; Dettmar, Ralf-Jürgen; Krause, Marita; Beck, Rainer; Stein, Yelena
2016-05-01
We present 1D cosmic ray transport models, numerically solving equations of pure advection and diffusion for the electrons and calculating synchrotron emission spectra. We find that for exponential halo magnetic field distributions advection leads to approximately exponential radio continuum intensity profiles, whereas diffusion leads to profiles that can be better approximated by a Gaussian function. Accordingly, the vertical radio spectral profiles for advection are approximately linear, whereas for diffusion they are of `parabolic' shape. We compare our models with deep Australia Telescope Compact Array observations of two edge-on galaxies, NGC 7090 and 7462, at λλ 22 and 6 cm. Our result is that the cosmic ray transport in NGC 7090 is advection dominated with V=150^{+80}_{-30} km s^{-1}, and that the one in NGC 7462 is diffusion dominated with D=3.0± 1.0 × 10^{28}E_GeV^{0.5} cm^2 s^{-1}. NGC 7090 has both a thin and thick radio disc with respective magnetic field scale heights of hB1 = 0.8 ± 0.1 kpc and hB2 = 4.7 ± 1.0 kpc. NGC 7462 has only a thick radio disc with hB2 = 3.8 ± 1.0 kpc. In both galaxies, the magnetic field scale heights are significantly smaller than what estimates from energy equipartition would suggest. A non-negligible fraction of cosmic ray electrons can escape from NGC 7090, so that this galaxy is not an electron calorimeter.
Lattice Boltzmann method for the fractional advection-diffusion equation.
Zhou, J G; Haygarth, P M; Withers, P J A; Macleod, C J A; Falloon, P D; Beven, K J; Ockenden, M C; Forber, K J; Hollaway, M J; Evans, R; Collins, A L; Hiscock, K M; Wearing, C; Kahana, R; Villamizar Velez, M L
2016-04-01
Mass transport, such as movement of phosphorus in soils and solutes in rivers, is a natural phenomenon and its study plays an important role in science and engineering. It is found that there are numerous practical diffusion phenomena that do not obey the classical advection-diffusion equation (ADE). Such diffusion is called abnormal or superdiffusion, and it is well described using a fractional advection-diffusion equation (FADE). The FADE finds a wide range of applications in various areas with great potential for studying complex mass transport in real hydrological systems. However, solution to the FADE is difficult, and the existing numerical methods are complicated and inefficient. In this study, a fresh lattice Boltzmann method is developed for solving the fractional advection-diffusion equation (LabFADE). The FADE is transformed into an equation similar to an advection-diffusion equation and solved using the lattice Boltzmann method. The LabFADE has all the advantages of the conventional lattice Boltzmann method and avoids a complex solution procedure, unlike other existing numerical methods. The method has been validated through simulations of several benchmark tests: a point-source diffusion, a boundary-value problem of steady diffusion, and an initial-boundary-value problem of unsteady diffusion with the coexistence of source and sink terms. In addition, by including the effects of the skewness β, the fractional order α, and the single relaxation time τ, the accuracy and convergence of the method have been assessed. The numerical predictions are compared with the analytical solutions, and they indicate that the method is second-order accurate. The method presented will allow the FADE to be more widely applied to complex mass transport problems in science and engineering. PMID:27176431
Anomalous diffusion of a tracer advected by wave turbulence
NASA Astrophysics Data System (ADS)
Balk, Alexander M.
2001-02-01
We consider the advection of a passive tracer when the velocity field is a superposition of random waves. Green's function for the turbulent transport (turbulent diffusion and turbulent drift) is derived. This Green's function is shown to imply sub-diffusive or super-diffusive behavior of the tracer. For the analysis we introduce the statistical near-identity transformation. The results are confirmed by numerical simulations.
Cellwise conservative unsplit advection for the volume of fluid method
NASA Astrophysics Data System (ADS)
Comminal, Raphaël; Spangenberg, Jon; Hattel, Jesper Henri
2015-02-01
We present a cellwise conservative unsplit (CCU) advection scheme for the volume of fluid method (VOF) in 2D. Contrary to other schemes based on explicit calculations of the flux balances, the CCU advection adopts a cellwise approach where the pre-images of the control volumes are traced backwards through the flow map. The donating regions of the fluxes are calculated via the streaklines of the grid intersections, represented as polygonal chains whose vertices are determined by backward tracing of particles injected in the flow at different times. High order accuracy is obtained from the fourth-order Runge-Kutta method, where intermediate velocities along pathlines are determined with quadratic temporal and bicubic spatial interpolations. The volumes of the donating regions are corrected in order to fulfill the discrete continuity of incompressible flows. Consequently, the calculation produces non-overlapping donating regions and pre-images with conforming edges to their neighbors, resulting in the conservativeness and the boundedness (liquid volume fraction inside the interval [ 0 , 1 ]) of the CCU advection scheme. Finally, the update of the liquid volume fractions is computed from the intersections of the pre-image polygons with the reconstructed interfaces. The CCU scheme is tested on several benchmark tests for the VOF advection, together with the standard piecewise linear interface calculation (PLIC). The geometrical errors of the CCU compare favorably with other unsplit VOF-PLIC schemes. Finally, potential improvements of the VOF method with the use of more precise interface representation techniques and the future extension of the CCU scheme to 3D are discussed.
Vertical Structure of Advection-dominated Accretion Flows
NASA Astrophysics Data System (ADS)
Zahra Zeraatgari, Fateme; Abbassi, Shahram
2015-08-01
We solve the set of hydrodynamic equations for optically thin advection-dominated accretion flows by assuming a radially self-similar spherical coordinate system (r,θ ,φ ). The disk is considered to be in steady state and axisymmetric. We define the boundary conditions at the pole and the equator of the disk and, to avoid singularity at the rotation axis, the disk is taken to be symmetric with respect to this axis. Moreover, only the {τ }rφ component of the viscous stress tensor is assumed, and we have set {v}θ =0. The main purpose of this study is to investigate the variation of dynamical quantities of the flow in the vertical direction by finding an analytical solution. As a consequence, we found that the advection parameter, {f}{adv}, varies along the θ direction and reaches its maximum near the rotation axis. Our results also show that, in terms of the no-outflow solution, thermal equilibrium still exists and consequently advection cooling can balance viscous heating.
Laser speckle contrast imaging is sensitive to advective flux
NASA Astrophysics Data System (ADS)
Khaksari, Kosar; Kirkpatrick, Sean J.
2016-07-01
Unlike laser Doppler flowmetry, there has yet to be presented a clear description of the physical variables that laser speckle contrast imaging (LSCI) is sensitive to. Herein, we present a theoretical basis for demonstrating that LSCI is sensitive to total flux and, in particular, the summation of diffusive flux and advective flux. We view LSCI from the perspective of mass transport and briefly derive the diffusion with drift equation in terms of an LSCI experiment. This equation reveals the relative sensitivity of LSCI to both diffusive flux and advective flux and, thereby, to both concentration and the ordered velocity of the scattering particles. We demonstrate this dependence through a short series of flow experiments that yield relationships between the calculated speckle contrast and the concentration of the scatterers (manifesting as changes in scattering coefficient), between speckle contrast and the velocity of the scattering fluid, and ultimately between speckle contrast and advective flux. Finally, we argue that the diffusion with drift equation can be used to support both Lorentzian and Gaussian correlation models that relate observed contrast to the movement of the scattering particles and that a weighted linear combination of these two models is likely the most appropriate model for relating speckle contrast to particle motion.
Sorption Energy Maps of Clay Mineral Surfaces
Cygan, Randall T.; Kirkpatrick, R. James
1999-07-19
A molecular-level understanding of mineral-water interactions is critical for the evaluation and prediction of the sorption properties of clay minerals that may be used in various chemical and radioactive waste disposal methods. Molecular models of metal sorption incorporate empirical energy force fields, based on molecular orbital calculations and spectroscopic data, that account for Coulombic, van der Waals attractive, and short-range repulsive energies. The summation of the non-bonded energy terms at equally-spaced grid points surrounding a mineral substrate provides a three dimensional potential energy grid. The energy map can be used to determine the optimal sorption sites of metal ions on the exposed surfaces of the mineral. By using this approach, we have evaluated the crystallographic and compositional control of metal sorption on the surfaces of kaolinite and illite. Estimates of the relative sorption energy and most stable sorption sites are derived based on a rigid ion approximation.
Sorption of polychlorinated biphenyls (PCBs) to sediment is a key process in determining their mobility, bioavailability, and chemical decomposition in aquatic environments. In order to examine the validity of currently used interpretation approaches for PCBs sorption, comparati...
Distinguishing resuspension and advection signals in a hypertidal estuary
NASA Astrophysics Data System (ADS)
Todd, David; Souza, Alex; Jago, Colin
2015-04-01
Terrestrial material is supplied to an estuary system by the river, while marine material is supplied by the sea. Whether the estuary acts as a trap or a bypass zone for SPM (suspended particulate matter) depends upon the properties and dynamics of both the estuary, including the tidal and residual behaviour of the currents, and the SPM, including particle sizes and settling velocities and concentration gradients, which together control the dynamics, such as the trapping efficiency, of the estuary. Whether an SPM signal is regarded as being one of resuspension or advection depends upon the area of interest, and therefore distinguishing between resuspension and advection can be complex. Material that is resuspended within the area of study is regarded as resuspension, while that which is resuspended outside, but passes through, the area of interest, is regarded as advection. The results of a measurement campaign undertaken in a hypertidal UK estuary during the pre-spring bloom February-March and post-spring bloom May-June are presented utilising a combination of acoustic and optical instruments, moorings, and CTD stations. A characteristic asymmetric "twin peak" signal is present during both time periods, implying the presence of both resuspension and advection. This is confirmed through the use of harmonic analysis. A seasonal variation in the relative importance of the resuspension and advection components is seen between the two observation periods, with the small (<122µm) and large (>122µm) particles displaying different behaviours and providing a strong indication of the presence of flocculation. Approximate point flux calculations showed a reduction in the horizontal gradient of concentration, and subsequently the flood dominance of sediment transport, between May-June and February-March. This has been attributed to changes in biological activity and atmospheric forcing between the two observational periods. Ebb-dominant concentrations brought about by the
Manganese Nitride Sorption Joule-Thomson Refrigerator
NASA Technical Reports Server (NTRS)
Jones, Jack A.; Phillips, Wayne M.
1992-01-01
Proposed sorption refrigeration system of increased power efficiency combines MnxNy sorption refrigeration stage with systems described in "Regenerative Sorption Refrigerator" (NPO-17630). Measured pressure-vs-composition isotherms for reversible chemisorption of N2 in MnxNy suggest feasibility to incorporate MnxNy chemisorption stage in Joule-Thomson cryogenic system. Discovery represents first known reversible nitrogen chemisorption compression system. Has potential in nitrogen-isotope separation, nitrogen purification, or contamination-free nitrogen compression.
Surface complexation model of uranyl sorption on Georgia kaolinite
Payne, T.E.; Davis, J.A.; Lumpkin, G.R.; Chisari, R.; Waite, T.D.
2004-01-01
The adsorption of uranyl on standard Georgia kaolinites (KGa-1 and KGa-1B) was studied as a function of pH (3-10), total U (1 and 10 ??mol/l), and mass loading of clay (4 and 40 g/l). The uptake of uranyl in air-equilibrated systems increased with pH and reached a maximum in the near-neutral pH range. At higher pH values, the sorption decreased due to the presence of aqueous uranyl carbonate complexes. One kaolinite sample was examined after the uranyl uptake experiments by transmission electron microscopy (TEM), using energy dispersive X-ray spectroscopy (EDS) to determine the U content. It was found that uranium was preferentially adsorbed by Ti-rich impurity phases (predominantly anatase), which are present in the kaolinite samples. Uranyl sorption on the Georgia kaolinites was simulated with U sorption reactions on both titanol and aluminol sites, using a simple non-electrostatic surface complexation model (SCM). The relative amounts of U-binding >TiOH and >AlOH sites were estimated from the TEM/EDS results. A ternary uranyl carbonate complex on the titanol site improved the fit to the experimental data in the higher pH range. The final model contained only three optimised log K values, and was able to simulate adsorption data across a wide range of experimental conditions. The >TiOH (anatase) sites appear to play an important role in retaining U at low uranyl concentrations. As kaolinite often contains trace TiO2, its presence may need to be taken into account when modelling the results of sorption experiments with radionuclides or trace metals on kaolinite. ?? 2004 Elsevier B.V. All rights reserved.
Spatially varying dispersion to model breakthrough curves.
Li, Guangquan
2011-01-01
Often the water flowing in a karst conduit is a combination of contaminated water entering at a sinkhole and cleaner water released from the limestone matrix. Transport processes in the conduit are controlled by advection, mixing (dilution and dispersion), and retention-release. In this article, a karst transport model considering advection, spatially varying dispersion, and dilution (from matrix seepage) is developed. Two approximate Green's functions are obtained using transformation of variables, respectively, for the initial-value problem and for the boundary-value problem. A numerical example illustrates that mixing associated with strong spatially varying conduit dispersion can cause strong skewness and long tailing in spring breakthrough curves. Comparison of the predicted breakthrough curve against that measured from a dye-tracing experiment between Ames Sink and Indian Spring, Northwest Florida, shows that the conduit dispersivity can be as large as 400 m. Such a large number is believed to imply strong solute interaction between the conduit and the matrix and/or multiple flow paths in a conduit network. It is concluded that Taylor dispersion is not dominant in transport in a karst conduit, and the complicated retention-release process between mobile- and immobile waters may be described by strong spatially varying conduit dispersion. PMID:21143474
Plutonium sorption and desorption behavior on bentonite.
Begg, James D; Zavarin, Mavrik; Tumey, Scott J; Kersting, Annie B
2015-03-01
Understanding plutonium (Pu) sorption to, and desorption from, mineral phases is key to understanding its subsurface transport. In this work we study Pu(IV) sorption to industrial grade FEBEX bentonite over the concentration range 10(-7)-10(-16) M to determine if sorption at typical environmental concentrations (≤10(-12) M) is the same as sorption at Pu concentrations used in most laboratory experiments (10(-7)-10(-11) M). Pu(IV) sorption was broadly linear over the 10(-7)-10(-16) M concentration range during the 120 d experimental period; however, it took up to 100 d to reach sorption equilibrium. At concentrations ≥10(-8) M, sorption was likely affected by additional Pu(IV) precipitation/polymerization reactions. The extent of sorption was similar to that previously reported for Pu(IV) sorption to SWy-1 Na-montmorillonite over a narrower range of Pu concentrations (10(-11)-10(-7) M). Sorption experiments with FEBEX bentonite and Pu(V) were also performed across a concentration range of 10(-11)-10(-7) M and over a 10 month period which allowed us to estimate the slow apparent rates of Pu(V) reduction on a smectite-rich clay. Finally, a flow cell experiment with Pu(IV) loaded on FEBEX bentonite demonstrated continued desorption of Pu over a 12 day flow period. Comparison with a desorption experiment performed with SWy-1 montmorillonite showed a strong similarity and suggested the importance of montorillonite phases in controlling Pu sorption/desorption reactions on FEBEX bentonite. PMID:25574607
NASA Astrophysics Data System (ADS)
Wilderbuer, Thomas; Duffy-Anderson, Janet T.; Stabeno, Phyllis; Hermann, Albert
2016-05-01
In an effort to better understand the physics of the eastern Bering Sea shelf current as it relates to flatfish advection to favorable near-shore areas, sets of multiple, satellite-tracked, oceanic drifters were released in 2010, 2012 and 2013. The release sites and dates were chosen to coincide with known spawning locations for northern rock sole (Lepidopsetta polyxystra) and known time of larval emergence. The drifters were drogued 5-each at 20 and 40 m in 2010 and 2012, and 4 at 40 m and 2 at 20 m in 2013. The locations of drifters were used to calculate divergence over a 90-day period that corresponds to the larval pelagic duration of Bering Sea shelf northern rock sole. Results indicate that there are alternating periods of positive and negative divergence with an overall trend toward drifter separation after 90 days, roughly the end of the rock sole planktonic larval period. Examination of the drifter behavior at the hourly scale indicates that semi-daily tidal forcing is the primary mechanism of drifter divergence and convergence. Field observations of early-stage northern rock sole larval distributions over the same period indicate that predominant oceanographic advection is northerly over the continental shelf among preflexion stages, though juveniles are predominantly found in nursery areas located ~ 400 km eastward and inshore. Evidence from drifter deployments suggests that behavioral movements during the postflexion and early juvenile larval phases that optimize eastward periodicity of tidal cycles is a viable mechanism to enhance eastward movement of northern rock sole larvae to favorable nursery grounds. A regional ocean modeling system (ROMS) was implemented to track the different rates of dispersion in simulations both with and without tidal forcing, and was used to estimate effective horizontal eddy diffusion in the case of both isobaric (fixed-depth) and Lagrangian (neutrally buoyant) particles. The addition of tidal forcing had a pronounced
Update on Advection-Diffusion Purge Flow Model
NASA Technical Reports Server (NTRS)
Brieda, Lubos
2015-01-01
Gaseous purge is commonly used in sensitive spacecraft optical or electronic instruments to prevent infiltration of contaminants and/or water vapor. Typically, purge is sized using simplistic zero-dimensional models that do not take into account instrument geometry, surface effects, and the dependence of diffusive flux on the concentration gradient. For this reason, an axisymmetric computational fluid dynamics (CFD) simulation was recently developed to model contaminant infiltration and removal by purge. The solver uses a combined Navier-Stokes and Advection-Diffusion approach. In this talk, we report on updates in the model, namely inclusion of a particulate transport model.
Numerical simulation of life cycles of advection warm fog
NASA Technical Reports Server (NTRS)
Hung, R. J.; Vaughan, O. H.
1977-01-01
The formation, development and dissipation of advection warm fog is investigated. The equations employed in the model include the equation of continuity, momentum and energy for the descriptions of density, wind component and potential temperature, respectively, together with two diffusion equations for the modification of water-vapor mixing ratio and liquid-water mixing ratios. A description of the vertical turbulent transfer of heat, moisture and momentum has been taken into consideration. The turbulent exchange coefficients adopted in the model are based on empirical flux-gradient relations.
Subsurface barrier design alternatives for confinement and controlled advection flow
Phillips, S.J.; Stewart, W.E.; Alexander, R.G.; Cantrell, K.J.; McLaughlin, T.J.
1994-02-01
Various technologies and designs are being considered to serve as subsurface barriers to confine or control contaminant migration from underground waste storage or disposal structures containing radioactive and hazardous wastes. Alternatives including direct-coupled flood and controlled advection designs are described as preconceptual examples. Prototype geotechnical equipment for testing and demonstration of these alternative designs tested at the Hanford Geotechnical Development and Test Facility and the Hanford Small-Tube Lysimeter Facility include mobile high-pressure injectors and pumps, mobile transport and pumping units, vibratory and impact pile drivers, and mobile batching systems. Preliminary laboratory testing of barrier materials and additive sequestering agents have been completed and are described.
Visualizing Vector Fields Using Line Integral Convolution and Dye Advection
NASA Technical Reports Server (NTRS)
Shen, Han-Wei; Johnson, Christopher R.; Ma, Kwan-Liu
1996-01-01
We present local and global techniques to visualize three-dimensional vector field data. Using the Line Integral Convolution (LIC) method to image the global vector field, our new algorithm allows the user to introduce colored 'dye' into the vector field to highlight local flow features. A fast algorithm is proposed that quickly recomputes the dyed LIC images. In addition, we introduce volume rendering methods that can map the LIC texture on any contour surface and/or translucent region defined by additional scalar quantities, and can follow the advection of colored dye throughout the volume.
Alginate and Algal-Based Beads for the Sorption of Metal Cations: Cu(II) and Pb(II).
Wang, Shengye; Vincent, Thierry; Faur, Catherine; Guibal, Eric
2016-01-01
Alginate and algal-biomass (Laminaria digitata) beads were prepared by homogeneous Ca ionotropic gelation. In addition, glutaraldehyde-crosslinked poly (ethyleneimine) (PEI) was incorporated into algal beads. The three sorbents were characterized by scanning electron microscopy (SEM) coupled with energy dispersive X-ray analysis (EDX): the sorption occurs in the whole mass of the sorbents. Sorption experiments were conducted to evaluate the impact of pH, sorption isotherms, and uptake kinetics. A special attention was paid to the effect of drying (air-drying vs. freeze-drying) on the mass transfer properties. For alginate, freeze drying is required for maintaining the porosity of the hydrogel, while for algal-based sorbents the swelling of the material minimizes the impact of the drying procedure. The maximum sorption capacities observed from experiments were 415, 296 and 218 mg Pb g(-1) and 112, 77 and 67 mg Cu g(-1) for alginate, algal and algal/PEI beads respectively. Though the sorption capacities of algal-beads decreased slightly (compared to alginate beads), the greener and cheaper one-pot synthesis of algal beads makes this sorbent more competitive for environmental applications. PEI in algal beads decreases the sorption properties in the case of the sorption of metal cations under selected experimental conditions. PMID:27598128
NASA Astrophysics Data System (ADS)
Moura, R. C.; Sherwin, S. J.; Peiró, J.
2016-02-01
This study addresses linear dispersion-diffusion analysis for the spectral/hp continuous Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion curves are obtained for the advection-diffusion problem and the role of multiple eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves' behaviour, we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion characteristics for under-resolved simulations of problems strongly dominated by either convection or diffusion. Subsequently, the linear advection equation augmented with spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of CG with SVV-based stabilization are verified to display similar non-smooth features in flow regions where convection is much stronger than dissipation or vice-versa, owing to a dependency of the standard SVV operator on a local Péclet number. First a modification is proposed to the traditional SVV scaling that enforces a globally constant Péclet number so as to avoid the previous issues. In addition, a new SVV kernel function is suggested and shown to provide a more regular behaviour for the eigencurves along with a consistent increase in resolution power for higher-order discretizations, as measured by the extent of the wavenumber range where numerical errors are negligible. The dissipation characteristics of CG with the SVV modifications suggested are then verified to be broadly equivalent to those obtained through upwinding in the discontinuous Galerkin (DG) scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is found to have a slightly higher resolution power for the same dissipation levels. These results show that improved CG-SVV characteristics can be pursued via different kernel functions with the aid of optimization algorithms.
Regenerative sorption compressors for cryogenic refrigeration
NASA Technical Reports Server (NTRS)
Bard, Steven; Jones, Jack A.
1990-01-01
Dramatic efficiency improvements for sorption coolers appear possible with use of compressor heat regeneration techniques. The general theory of sorption compressor heat regeneration is discussed in this paper, and several design concepts are presented. These designs result in long-life, low-vibration cryocoolers that potentially have efficiencies comparable to Stirling refrigerators for 65 to 90 K spacecraft instrument cooling applications.
Sorption of paracetamol onto biomaterials.
Ferchichi, Maroua; Dhaouadi, Hatem
2016-01-01
Pharmaceutical residues released into the environment are posing more and more public health problems. It is worthwhile to study the retention of pharmaceuticals residues by adsorption on solid supports. Batch sorption experiments are intended to identify the adsorption isotherms of the pharmaceutically active ingredient on the biomaterials. The results obtained in this study have shown that the retention possibilities of these compounds by bio-adsorbents (clay and sand) are not significant. The negligible sorption for these media is explained by the low hydrophobicity of paracetamol (Log K(ow) = 0.46). The retention of paracetamol on the dehydrated sewage sludge and on Posidonia oceanica showed a relatively significant adsorption with a maximal quantity of 0.956 mg g(-1) and 1.638 mg g(-1) for the dehydrate sludge and P. oceanica, respectively. On the other hand, the study of paracetamol retention on the powdered activated carbon showed a high adsorption capacity of about 515.27 mg g(-1). Isotherm data show a good fit with Langmuir's model. An infrared analysis is carried out. It shows identical bands before and after adsorption, with some modifications. PMID:27387007
Numerical Modelling of Mesoscale Atmospheric Dispersion.
NASA Astrophysics Data System (ADS)
Moran, Michael D.
Mesoscale atmospheric dispersion is more complicated than smaller-scale dispersion because the mean wind field can no longer be considered steady or horizontally homogeneous over mesoscale time and space scales. Wind shear also plays a more important role on the mesoscale, and horizontal dispersion can be enhanced and even dominated by vertical wind shear through either the simultaneous or delayed interaction of horizontal differential advection and vertical mixing over one or two diurnal periods. The CSU mesoscale atmospheric dispersion modelling system has been used in this study to simulate the transport and diffusion of a perfluorocarbon gas for episodic releases made during two North American mesoscale dispersion field experiments, the 1980 Great Plains tracer experiment and the 1983 Cross-Appalachian Tracer Experiment (CAPTEX). Ground -level and elevated tracer concentrations were measured out to distances of 600 km from the source in the first experiment and 1100 km in the second. The physiography of the two experimental domains was quite different, permitting isolation and examination of the roles of terrain forcing and differential advection in mesoscale atmospheric dispersion. Suites of numerical experiments of increasing complexity were carried out for both case studies. The experiments differed in the realism of their representation of both the synoptic-scale flow and the underlying terrain. The Great Plains nocturnal low-level jet played an important role in the first case while temporal changes in the synoptic -scale flow were very significant in the second case. The contributions of differential advection and mesoscale deformation to mesoscale dispersion dominated those of small-scale turbulent diffusion for both cases, and Pasquill's (1962) delayed-shear-enhancement mechanism for lateral dispersion was found to be particularly important. This study was also the first quantitative evaluation of the CSU mesoscale dispersion modelling system with
Ferrous iron sorption by hydrous metal oxides.
Nano, Genevieve Villaseñor; Strathmann, Timothy J
2006-05-15
Ferrous iron is critical to a number of biogeochemical processes that occur in heterogeneous aquatic environments, including the abiotic reductive transformation of subsurface contaminants. The sorption of Fe(II) to ubiquitous soil minerals, particularly iron-free mineral phases, is not well understood. Colloidal TiO2, gamma-AlOOH, and gamma-Al2O2 were used as model hydrous oxides to investigate Fe(II) sorption to iron-free mineral surfaces. Rapid Fe(II) sorption during the first few hours is followed by a much slower uptake process that continues for extended periods (at least 30 days). For equivalent solution conditions, the extent of Fe(II) sorption decreases in the order TiO2 >gamma-Al2O3 >gamma-AlOOH. Short-term equilibrium sorption data measured over a wide range of conditions (pH, ionic strength, Fe(II)-to-sorbent ratio) are well described by the diffuse double layer model. Fe(II) sorption to TiO2 is best described by a single-site model that considers formation of two surface complexes, SOFe+ and SOFeOH0. For gamma-AlOOH and gamma-Al2O3, sorption data are best described by a two-site model that considers formation of SOFe+ complexes at weak- and strong-binding surface sites. Accurate description of sorption data for higher Fe(II) concentrations at alkaline pH conditions requires the inclusion of a Fe(II) surface precipitation reaction in the model formulation. The presence of common groundwater constituents (calcium, sulfate, bicarbonate, or fulvic acid) had no significant effect on Fe(II) sorption. These results demonstrate that iron-free soil minerals can exert a significant influence on Fe(II) sorption and speciation in heterogeneous aquatic systems. PMID:16337955
Sorption kinetics of hexadecyltrimethylammonium on natural clinoptilolite
Li, Z.
1999-09-14
Sorption kinetics of hexadecyltrimethylammonium (HDTMA) chloride on a natural clinoptilolite was studied in this research. The amount of HDTMA sorbed is a function of the initial HDTMA input and the sorption time. When the initial HDTMA input is less than the external cation-exchange capacity of the clinoptilolite, the HDTMA sorption is fast and equilibrium can be established in 1 h. As the initial HDTMA input is greater than the external cation-exchange capacity of clinoptilolite, which will result in more than a monolayer HDTMA surface coverage, the time for HDTMA sorption to reach equilibrium increases exponentially. The HDTMA sorption maximum on clinoptilolite increases logarithmically with mixing time. The counterion solution concentration data suggest that at the initial stage HDTMA molecules sorb on the zeolite via micelle forms, which is manifested by a decrease in chloride solution concentration with time. When HDTMA solution concentration is depleted to less than its critical micelle concentration, the adsorbed micelles (admicelles) rearrange themselves to a more stable monolayer or bilayer configuration, which is reflected by an increase in counterion solution concentration due to the desorption of chloride from admicelles. The time required for the surface rearrangement increases exponentially as the HDTMA input increases. The data of HDTMA sorption kinetics were fitted to different kinetic models, and the parabolic diffusion model fits the data best for the HDTMA sorption, counterion sorption at the initial stage and counterion desorption at the rearrangement stage. Thus, the sorption of HDTMA on clinoptilolite surfaces is diffusion controlled. The results also indicate that it is incomplete to discuss surfactant sorption without counterion concentration data.
A cryogenic circulating advective multi-pass absorption cell
NASA Astrophysics Data System (ADS)
Stockett, M. H.; Lawler, J. E.
2012-03-01
A novel absorption cell has been developed to enable a spectroscopic survey of a broad range of polycyclic aromatic hydrocarbons (PAH) under astrophysically relevant conditions and utilizing a synchrotron radiation continuum to test the still controversial hypothesis that these molecules or their ions could be carriers of the diffuse interstellar bands. The cryogenic circulating advective multi-pass absorption cell resembles a wind tunnel; molecules evaporated from a crucible or injected using a custom gas feedthrough are entrained in a laminar flow of cryogenically cooled buffer gas and advected into the path of the synchrotron beam. This system includes a multi-pass optical White cell enabling absorption path lengths of hundreds of meters and a detection sensitivity to molecular densities on the order of 107 cm-3. A capacitively coupled radio frequency dielectric barrier discharge provides ionized and metastable buffer gas atoms for ionizing the candidate molecules via charge exchange and the Penning effect. Stronger than expected clustering of PAH molecules has slowed efforts to record gas phase PAH spectra at cryogenic temperatures, though such clusters may play a role in other interstellar phenomena.
Advection-Based Sparse Data Management for Visualizing Unsteady Flow.
Guo, Hanqi; Zhang, Jiang; Liu, Richen; Liu, Lu; Yuan, Xiaoru; Huang, Jian; Meng, Xiangfei; Pan, Jingshan
2014-12-01
When computing integral curves and integral surfaces for large-scale unsteady flow fields, a major bottleneck is the widening gap between data access demands and the available bandwidth (both I/O and in-memory). In this work, we explore a novel advection-based scheme to manage flow field data for both efficiency and scalability. The key is to first partition flow field into blocklets (e.g. cells or very fine-grained blocks of cells), and then (pre)fetch and manage blocklets on-demand using a parallel key-value store. The benefits are (1) greatly increasing the scale of local-range analysis (e.g. source-destination queries, streak surface generation) that can fit within any given limit of hardware resources; (2) improving memory and I/O bandwidth-efficiencies as well as the scalability of naive task-parallel particle advection. We demonstrate our method using a prototype system that works on workstation and also in supercomputing environments. Results show significantly reduced I/O overhead compared to accessing raw flow data, and also high scalability on a supercomputer for a variety of applications. PMID:26356969
THE ADVECTION OF SUPERGRANULES BY THE SUN'S AXISYMMETRIC FLOWS
Hathaway, David H.; Williams, Peter E.; Rosa, Kevin Dela; Cuntz, Manfred E-mail: peter.williams@nasa.go
2010-12-10
We show that the motions of supergranules are consistent with a model in which they are simply advected by the axisymmetric flows in the Sun's surface shear layer. We produce a 10 day series of simulated Doppler images at a 15 minute cadence that reproduces most spatial and temporal characteristics seen in the SOHO/MDI Doppler data. Our simulated data have a spectrum of cellular flows with just two components-a granule component that peaks at spherical wavenumbers of about 4000 and a supergranule component that peaks at wavenumbers of about 110. We include the advection of these cellular components by the axisymmetric flows-differential rotation and meridional flow-whose variations with latitude and depth (wavenumber) are consistent with observations. We mimic the evolution of the cellular pattern by introducing random variations to the phases of the spectral components at rates that reproduce the levels of cross-correlation as functions of time and latitude. Our simulated data do not include any wave-like characteristics for the supergranules yet can reproduce the rotation characteristics previously attributed to wave-like behavior. We find rotation rates which appear faster than the actual rotation rates and attribute this to projection effects. We find that the measured meridional flow does accurately represent the actual flow and that the observations indicate poleward flow to 65{sup 0}-70{sup 0} latitude with equatorward countercells in the polar regions.
Moisture advection to the Arctic : forecasted, analysed and observed
NASA Astrophysics Data System (ADS)
Dufour, Ambroise; Zolina, Olga
2015-04-01
Besides its contribution to the Arctic hydrological budget, moisture imports from mid-latitudes are also influential on shorter time scales since water vapour advection tends to occur together with extratropical cyclones. Influx of moisture to the Arctic cause the formation of clouds that have an immediate impact on the surface energy budget especially in winter. In the long run, inaccuracies in the description of cloud cover and phase lead to temperature biases in CMIP5 models. The ECMWF workshop on polar prediction has highlighted moisture advection as one of the problematic physical processes limiting the quality of forecasts. Verifying the accuracy of medium-term forecasts is of interest beyond weather prediction : it points to the ability of models to bring adequate quantities of moisture to the Arctic when they are less constrained by observations than in analyses. In this study, we have compared forecasted moisture flux fields with analyses and observations over the period 2000-2010. ECMWF's ERA-Interim provided the forecasts, extending to ten days. For the analyses, in addition to ERA-Interim, we used the Arctic System Reanalysis whose forecast model is optimized for the polar regions and runs at high resolution (30 km). Finally, the Integrated Global Radiosonde Archive data over the Arctic allowed a validation by observations.
Horizontal advection, diffusion and plankton spectra at the sea surface.
NASA Astrophysics Data System (ADS)
Bracco, A.; Clayton, S.; Pasquero, C.
2009-04-01
Plankton patchiness is ubiquitous in the oceans, and various physical and biological processes have been proposed as its generating mechanisms. However, a coherent statement on the problem is missing, due to both a small number of suitable observations and to an incomplete understanding of the properties of reactive tracers in turbulent media. Abraham (1998) suggested that horizontal advection may be the dominant process behind the observed distributions of phytoplankton and zooplankton, acting to mix tracers with longer reaction times (Rt) down to smaller scales. Conversely, Mahadevan and Campbell (2002) attributed the relative distributions of sea surface temperature and phytoplankton to small scale upwelling, where tracers with longer Rt are able to homogenize more than those with shorter reaction times. Neither of the above mechanisms can explain simultaneously the (relative) spectral slopes of temperature, phytoplankton and zooplankton. Here, with a simple advection model and a large suite of numerical experiments, we concentrate on some of the physical processes influencing the relative distributions of tracers at the ocean surface, and we investigate: 1) the impact of the spatial scale of tracer supply; 2) the role played by coherent eddies on the distribution of tracers with different Rt; 3) the role of diffusion (so far neglected). We show that diffusion determines the distribution of temperature, regardless of the nature of the forcing. We also find that coherent structures together with differential diffusion of tracers with different Rt impact the tracer distributions. This may help in understanding the highly variable nature of observed plankton spectra.
Toward enhanced subsurface intervention methods using chaotic advection.
Trefry, Michael G; Lester, Daniel R; Metcalfe, Guy; Ord, Alison; Regenauer-Lieb, Klaus
2012-01-01
Many intervention activities in the terrestrial subsurface involve the need to recover/emplace distributions of scalar quantities (e.g. dissolved phase concentrations or heat) from/in volumes of saturated porous media. These scalars can be targeted by pump-and-treat methods or by amendment technologies. Application examples include in-situ leaching for metals, recovery of dissolved contaminant plumes, or utilizing heat energy in geothermal reservoirs. While conventional pumping methods work reasonably well, costs associated with maintaining pumping schedules are high and improvements in efficiency would be welcome. In this paper we discuss how transient switching of the pressure at different wells can intimately control subsurface flow, generating a range of "programmed" flows with various beneficial characteristics. Some programs produce chaotic flows which accelerate mixing, while others create encapsulating flows which can isolate fluid zones for lengthy periods. In a simplified model of an aquifer subject to balanced pumping, chaotic flow topologies have been predicted theoretically and verified experimentally using Hele-Shaw cells. Here, a survey of the key characteristics of chaotic advection is presented. Mathematical methods are used to show how these characteristics may translate into practical situations involving regional flows and heterogeneity. The results are robust to perturbations, and withstand significant aquifer heterogeneity. It is proposed that chaotic advection may form the basis of new efficient technologies for groundwater interventions. PMID:21600670
Chaotic advection in 2D anisotropic porous media
NASA Astrophysics Data System (ADS)
Varghese, Stephen; Speetjens, Michel; Trieling, Ruben; Toschi, Federico
2015-11-01
Traditional methods for heat recovery from underground geothermal reservoirs employ a static system of injector-producer wells. Recent studies in literature have shown that using a well-devised pumping scheme, through actuation of multiple injector-producer wells, can dramatically enhance production rates due to the increased scalar / heat transport by means of chaotic advection. However the effect of reservoir anisotropy on kinematic mixing and heat transport is unknown and has to be incorporated and studied for practical deployment in the field. As a first step, we numerically investigate the effect of anisotropy (both magnitude and direction) on (chaotic) advection of passive tracers in a time-periodic Darcy flow within a 2D circular domain driven by periodically reoriented diametrically opposite source-sink pairs. Preliminary results indicate that anisotropy has a significant impact on the location, shape and size of coherent structures in the Poincare sections. This implies that the optimal operating parameters (well spacing, time period of well actuation) may vary strongly and must be carefully chosen so as to enhance subsurface transport. This work is part of the research program of the Foundation for Fundamental Research on Matter (FOM), which is part of Netherlands Organisation for Scientific Research (NWO). This research program is co-financed by Shell Global Solutions International B.V.
Multiple anisotropic collisions for advection-diffusion Lattice Boltzmann schemes
NASA Astrophysics Data System (ADS)
Ginzburg, Irina
2013-01-01
This paper develops a symmetrized framework for the analysis of the anisotropic advection-diffusion Lattice Boltzmann schemes. Two main approaches build the anisotropic diffusion coefficients either from the anisotropic anti-symmetric collision matrix or from the anisotropic symmetric equilibrium distribution. We combine and extend existing approaches for all commonly used velocity sets, prescribe most general equilibrium and build the diffusion and numerical-diffusion forms, then derive and compare solvability conditions, examine available anisotropy and stable velocity magnitudes in the presence of advection. Besides the deterioration of accuracy, the numerical diffusion dictates the stable velocity range. Three techniques are proposed for its elimination: (i) velocity-dependent relaxation entries; (ii) their combination with the coordinate-link equilibrium correction; and (iii) equilibrium correction for all links. Two first techniques are also available for the minimal (coordinate) velocity sets. Even then, the two-relaxation-times model with the isotropic rates often gains in effective stability and accuracy. The key point is that the symmetric collision mode does not modify the modeled diffusion tensor but it controls the effective accuracy and stability, via eigenvalue combinations of the opposite parity eigenmodes. We propose to reduce the eigenvalue spectrum by properly combining different anisotropic collision elements. The stability role of the symmetric, multiple-relaxation-times component, is further investigated with the exact von Neumann stability analysis developed in diffusion-dominant limit.
A cryogenic circulating advective multi-pass absorption cell
Stockett, M. H.; Lawler, J. E.
2012-03-15
A novel absorption cell has been developed to enable a spectroscopic survey of a broad range of polycyclic aromatic hydrocarbons (PAH) under astrophysically relevant conditions and utilizing a synchrotron radiation continuum to test the still controversial hypothesis that these molecules or their ions could be carriers of the diffuse interstellar bands. The cryogenic circulating advective multi-pass absorption cell resembles a wind tunnel; molecules evaporated from a crucible or injected using a custom gas feedthrough are entrained in a laminar flow of cryogenically cooled buffer gas and advected into the path of the synchrotron beam. This system includes a multi-pass optical White cell enabling absorption path lengths of hundreds of meters and a detection sensitivity to molecular densities on the order of 10{sup 7} cm{sup -3}. A capacitively coupled radio frequency dielectric barrier discharge provides ionized and metastable buffer gas atoms for ionizing the candidate molecules via charge exchange and the Penning effect. Stronger than expected clustering of PAH molecules has slowed efforts to record gas phase PAH spectra at cryogenic temperatures, though such clusters may play a role in other interstellar phenomena.
A cryogenic circulating advective multi-pass absorption cell.
Stockett, M H; Lawler, J E
2012-03-01
A novel absorption cell has been developed to enable a spectroscopic survey of a broad range of polycyclic aromatic hydrocarbons (PAH) under astrophysically relevant conditions and utilizing a synchrotron radiation continuum to test the still controversial hypothesis that these molecules or their ions could be carriers of the diffuse interstellar bands. The cryogenic circulating advective multi-pass absorption cell resembles a wind tunnel; molecules evaporated from a crucible or injected using a custom gas feedthrough are entrained in a laminar flow of cryogenically cooled buffer gas and advected into the path of the synchrotron beam. This system includes a multi-pass optical White cell enabling absorption path lengths of hundreds of meters and a detection sensitivity to molecular densities on the order of 10(7) cm(-3). A capacitively coupled radio frequency dielectric barrier discharge provides ionized and metastable buffer gas atoms for ionizing the candidate molecules via charge exchange and the Penning effect. Stronger than expected clustering of PAH molecules has slowed efforts to record gas phase PAH spectra at cryogenic temperatures, though such clusters may play a role in other interstellar phenomena. PMID:22462957
Chaotic Advection, Fluid Spreading, and Groundwater Contaminant Plumes
NASA Astrophysics Data System (ADS)
Neupauer, R. M.; Mays, D. C.
2011-12-01
In situ remediation of contaminated groundwater requires degradation reactions at the interface between the contaminant plume and an injected treatment solution containing chemical or biological amendments. Therefore a promising approach to accelerate in situ remediation is to elongate the interface between the contaminant plume and treatment solution through fluid spreading. The literature on chaotic advection describes how to accomplish spreading in laminar flows, which lack the turbulent eddies that provide spreading in streams or engineered reactors. A key result from the literature on chaotic advection is that spreading is inherent in the vicinity of certain periodic points, which are points to which fluid particles return in successive iterations of chaotic flows. Specifically, spreading is enhanced near the stable and unstable manifolds associated with hyperbolic periodic points. We investigate the transient flow created with a four-well system in which wells are operated sequentially as either injection wells or extraction wells. In particular, we identify the periodic points and demonstrate that fluid spreading occurs nearby. For appropriately designed injection and extraction sequences, the periodic points are located near the interface between the contaminant plume and treatment solution, leading to elongation of the interface, with expected benefits of enhanced reaction and accelerated remediation.
Positivity-preserving numerical schemes for multidimensional advection
NASA Technical Reports Server (NTRS)
Leonard, B. P.; Macvean, M. K.; Lock, A. P.
1993-01-01
This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.
OBSERVATION OF MAGNETIC RECONNECTION DRIVEN BY GRANULAR SCALE ADVECTION
Zeng Zhicheng; Cao Wenda; Ji Haisheng
2013-06-01
We report the first evidence of magnetic reconnection driven by advection in a rapidly developing large granule using high spatial resolution observations of a small surge event (base size {approx} 4'' Multiplication-Sign 4'') with the 1.6 m aperture New Solar Telescope at the Big Bear Solar Observatory. The observations were carried out in narrowband (0.5 A) He I 10830 A and broadband (10 A) TiO 7057 A. Since He I 10830 A triplet has a very high excitation level and is optically thin, its filtergrams enable us to investigate the surge from the photosphere through the chromosphere into the lower corona. Simultaneous space data from the Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory were used in the analysis. It is shown that the surge is spatio-temporally associated with magnetic flux emergence in the rapidly developing large granule. During the development of the granule, its advecting flow ({approx}2 km s{sup -1}) squeezed the magnetic flux into an intergranular lane area, where a magnetic flux concentration was formed and the neighboring flux with opposite magnetic polarity was canceled. During the cancellation, the surge was produced as absorption in He I 10830 A filtergrams while simultaneous EUV brightening occurred at its base. The observations clearly indicate evidence of a finest-scale reconnection process driven by the granule's motion.
An unusual H2 sorption mechanism in PCN-14: insights from molecular simulation.
Pham, Tony; Forrest, Katherine A; Space, Brian
2016-08-01
Simulations of H2 sorption were performed in PCN-14, a metal-organic framework (MOF) that consists of Cu(2+) ions coordinated to 5,5'-(9,10-anthracenediyl)diisophthalate (adip) linkers. This MOF displays an excess H2 uptake of 2.70 wt% at 77 K and 1.0 atm and an initial H2Qst value of 8.6 kJ mol(-1) according to previous experimental measurements. The experimental H2 sorption isotherms and Qst values in PCN-14 were reproduced in simulations using well-known H2 potentials that have been widely used for MOF-H2 theoretical studies. H2 sorption in PCN-14 was dominated by repulsion/dispersion energetics; this allowed the experimental observables to be reproduced by a model that includes only Lennard-Jones parameters. The most energetically favorable H2 sorption site in PCN-14 corresponds to sorption within a small cage that is enclosed by three [Cu2(O2CR)4] units and three adip linkers. The anthracenyl rings of the adip linkers represent the secondary sorption sites within the MOF. In contrast to expectations, sorption of H2 onto the Cu(2+) ions of the copper paddlewheels was not observed within the simulations at low loading. The simulations revealed that the open-metal sites in PCN-14 were occupied at high loading. Control simulations of H2 sorption in PCN-14 for different cases in which the partial positive charge of one of the paddlewheel Cu(2+) ions was increased relative to the other revealed that sorption onto the open-metal sites can be captured if there is a very high positive charge on the metal. Otherwise, the calculated partial charge for the Cu(2+) ions in PCN-14 in this work was not high enough in magnitude to facilitate strong H2-metal interactions in simulation. This study shows the power of using molecular simulations to elucidate an unusual H2 sorption behavior in a MOF. PMID:27426916
Technology Transfer Automated Retrieval System (TEKTRAN)
The convective-dispersive, or advective-dispersive, or CDE, equation has long been the model of choice for solute transport in soils. Using the total mass of soluble salts in soil profile to evaluate changes in salinity due to irrigation can be beneficial when the spatial variability of soil salini...
Characterization of lead sorption by the natural and Fe(III)-modified zeolite
NASA Astrophysics Data System (ADS)
Kragović, Milan; Daković, Aleksandra; Marković, Marija; Krstić, Jugoslav; Gatta, G. Diego; Rotiroti, Nicola
2013-10-01
The influence of contact time, temperature and particle size on lead sorption by the natural and Fe(III)-modified zeolites was investigated. Characterization of the natural and Fe(III)-modified zeolite before and after lead sorption was performed by determination of textural properties, by scanning electron microscopy and X-ray spectroscopy in energy-dispersive mode (SEM-EDS), transmission electron microscopy (TEM) and X-ray powder diffraction (XRPD) analysis. Lead sorption kinetics at 303-333 K, best represented by the pseudo-second order model and activation energy (13.5 and 8.5 kJ/mol for the natural and Fe(III)-modified zeolite respectively) confirmed an activated chemical sorption. Desorption experiments indicated that lead was irreversibly sorbed on both zeolites. XRPD, TEM and SEM results showed that modification of the natural zeolite with Fe(III) ions did not change its crystal structure and iron is mainly located at the zeolite surface, likely in form of amorphous iron oxy-hydroxides. Specific surface area significantly increases after modification of the natural zeolite with Fe(III) ions (from 30.2 for the natural to 52.5 m2/g for Fe(III)-modified zeolite). Characterization of both lead saturated sorbents suggested that besides ion exchange, lead is both chemisorbed and precipitated at their surfaces, and presence of amorphous iron in Fe(III)-modified zeolite favors sorption of lead.
Numerical modelling of mesoscale atmospheric dispersion. (Volumes I and II)
Moran, M.D.
1992-01-01
Mesoscale atmospheric dispersion is more complicated than smaller-scale dispersion because the mean wind field can no longer be considered steady or horizontally homogeneous over mesoscale time and space scales. Wind shear also plays an important role on the mesoscale, and horizontal dispersion can be enhanced and even dominated by vertical wind shear through either the simultaneous or delayed interaction of horizontal differential advection and vertical mixing. The CSU mesoscale atmospheric dispersion modelling system has been used in this study to simulate the transport and diffusion of a perfluorocarbon gas for episodic releases made during two mesoscale dispersion field experiments. The physiography of the two experimental domains was quite different, permitting isolation and examination of the roles of terrain forcing and differential advection in mesoscale atmospheric dispersion. Suites of numerical experiments of increasing complexity were carried out for both case studies. The experiments differed in the realism of their representation of both the synoptic-scale flow and the underlying terrain. The contributions of differential advection and mesoscale deformation to mesoscale dispersion dominated those of small-scale turbulent diffusion for both cases, and Pasquill's (1962) delayed-shear-enhancement mechanism for lateral dispersion was found to be particularly important. This study was also the first quantitative evaluation of the CSU mesoscale dispersion modelling system with episodic mesoscale dispersion field data. The modelling system showed considerable skill in predicting quantitative tracer-cloud characteristics such as peak concentration, maximum cloud width, arrival time, transit time, and crosswind integrated exposure. Model predictions also compared favorably with predictions made by a number of other mesoscale dispersion models for the same two case studies.
NASA Astrophysics Data System (ADS)
Hammer, Daniel X.; Noojin, Gary D.; Thomas, Robert J.; Stolarski, David J.; Rockwell, Benjamin A.; Welch, Ashley J.
1999-06-01
Spectrally resolved white-light interferometry (SRWLI) was used to measure the wavelength dependence of refractive index (i.e., dispersion) for various ocular components. The accuracy of the technique was assessed by measurement of fused silica and water, the refractive indices of which have been measured at several different wavelengths. The dispersion of bovine and rabbit aqueous and vitreous humor was measured from 400 to 1100 nm. Also, the dispersion was measured from 400 to 700 nm for aqueous and vitreous humor extracted from goat and rhesus monkey eyes. For the humors, the dispersion did not deviate significantly from water. In an additional experiment, the dispersion of aqueous and vitreous humor that had aged up to a month was compared to freshly harvested material. No difference was found between the fresh and aged media. An unsuccessful attempt was also made to use the technique for dispersion measurement of bovine cornea and lens. Future refinement may allow measurement of the dispersion of cornea and lens across the entire visible and near-infrared wavelength band. The principles of white- light interferometry including image analysis, measurement accuracy, and limitations of the technique, are discussed. In addition, alternate techniques and previous measurements of ocular dispersion are reviewed.
Fe2+ sorption onto nontronite (NAu-2)
Jaisi, Deb P.; Liu, Chongxuan; Dong, Hailiang; Blake, Ruth E.; Fein, Jeremy B.
2008-11-15
The sorption of ferrous iron to a clay mineral, nontronite (NAu-2, a ferruginous smectite), was investigated under strictly anoxic conditions as a function of pH (3-10), Fe2+ concentration (0.01-50 mM), equilibration time (1 to 35 days), and ionic strength (0.01 to 0.5 M NaClO4). The surface properties of NAu-2 were independently characterized to determine its fixed charge and amphoteric site properties for the interpretation of Fe2+ sorption. Fe2+ sorption to NAu-2 was strongly dependent on pH and ionic strength, reflecting the coupled effects of Fe2+ sorption through ion exchange and surface complexation reactions. Fe2+ sorption to NAu-2 increased with increasing pH from pH 2.5 to 4.5, reached the first plateau from pH 4.5 to 7.0, increased again with increasing pH from pH 7.0 to 8.5, and reached the maximum (the second plateau) above pH 8.5. The Fe2+ sorption below pH 7.0 increased with decreasing ionic strength. The differences of Fe2+ sorption at different ionic strength, however, diminished with increasing equilibration time. The Fe2+ sorption from pH 4.5 to 7.0 increased with increasing equilibration time up to 35 days and showed stronger kinetic behavior in higher ionic strength solutions. Increasing Fe2+ concentration nonlinearly decreased Fe2+ sorption, resulting from the mass action and site saturation. An equilibrium model that integrated ion exchange, surface complexation and aqueous speciation reactions reasonably well described the Fe2+ sorption data as a function of pH, ionic strength, and Fe2+ concentration measured at 24 hours of equilibration. Model calculations showed that species Fe(OH)+ was required to describe Fe2+ sorption above pH 8.0. Model calculation also implied that the ion exchange reactions was responsible for the kinetic behavior of Fe2+ sorption. Overall, this study demonstrated that Fe2+ sorption to NAu-2 was complexly affected by equilibrium and kinetic processes.
Colloid particle size-dependent dispersivity
NASA Astrophysics Data System (ADS)
Chrysikopoulos, C. V.; Katzourakis, V. E.
2014-12-01
Laboratory and field studies have demonstrated that dispersion coefficients evaluated by fitting advection-dispersion transport models to nonreactive tracer breakthrough curves do not adequately describe colloid transport under the same flow field conditions. Here an extensive laboratory study was undertaken to assess whether the dispersivity, which traditionally has been considered to be a property of the porous medium, is dependent on colloid particle size and interstitial velocity. A total of 49 colloid transport experiments were performed in columns packed with glass beads under chemically unfavorable colloid attachment conditions. Nine different colloid diameters, and various flow velocities were examined. The breakthrough curves were successfully simulated with a mathematical model describing colloid transport in homogeneous, water saturated porous media. The results demonstrated that the dispersivity is positively correlated with colloid particle size, and increases with increasing velocity.
Sorption of radiocesium by active silica.
McCulloch, C E; Crawford, R W; Angus, M J; Glasser, F P; Rahman, A A
1984-05-01
Cement and cement components have previously been shown to exhibit negligible sorption for cesium. Pyrogenic silica has been examined as an additive to cement materials for its ability to reduce the leachability of cesium and to provide a host material with permanent sorption sites. The incorporation of silica into cement composites can also improve the physical characteristics and strength of these materials as long ages. At neutral pH values, there is significant sorption of cesium by silica, but in high pH regimes, such as occur in cement environments, initial sorption is enhanced but this high level of sorption is followed by a gradual release of Cs. This apparent desorption is due to the consumption of SiO2 by Ca(OH)2 to form products which have little sorption potential for cesium. If, however, sufficient SiO2 is added to the system initially such that an excess remains after satisfying the demands of the Ca(OH)2 reaction, permanent sorption sites for cesium may be created. PMID:6327573
Thermally driven advection for radioxenon transport from an underground nuclear explosion
NASA Astrophysics Data System (ADS)
Sun, Yunwei; Carrigan, Charles R.
2016-05-01
Barometric pumping is a ubiquitous process resulting in migration of gases in the subsurface that has been studied as the primary mechanism for noble gas transport from an underground nuclear explosion (UNE). However, at early times following a UNE, advection driven by explosion residual heat is relevant to noble gas transport. A rigorous measure is needed for demonstrating how, when, and where advection is important. In this paper three physical processes of uncertain magnitude (oscillatory advection, matrix diffusion, and thermally driven advection) are parameterized by using boundary conditions, system properties, and source term strength. Sobol' sensitivity analysis is conducted to evaluate the importance of all physical processes influencing the xenon signals. This study indicates that thermally driven advection plays a more important role in producing xenon signals than oscillatory advection and matrix diffusion at early times following a UNE, and xenon isotopic ratios are observed to have both time and spatial dependence.
Bacterial sorption of heavy metals.
Mullen, M D; Wolf, D C; Ferris, F G; Beveridge, T J; Flemming, C A; Bailey, G W
1989-01-01
Four bacteria, Bacillus cereus, B. subtilis, Escherichia coli, and Pseudomonas aeruginosa, were examined for the ability to remove Ag+, Cd2+, Cu2+, and La3+ from solution by batch equilibration methods. Cd and Cu sorption over the concentration range 0.001 to 1 mM was described by Freundlich isotherms. At 1 mM concentrations of both Cd2+ and Cu2+, P. aeruginosa and B. cereus were the most and least efficient at metal removal, respectively. Freundlich K constants indicated that E. coli was most efficient at Cd2+ removal and B. subtilis removed the most Cu2+. Removal of Ag+ from solution by bacteria was very efficient; an average of 89% of the total Ag+ was removed from the 1 mM solution, while only 12, 29, and 27% of the total Cd2+, Cu2+, and La3+, respectively, were sorbed from 1 mM solutions. Electron microscopy indicated that La3+ accumulated at the cell surface as needlelike, crystalline precipitates. Silver precipitated as discrete colloidal aggregates at the cell surface and occasionally in the cytoplasm. Neither Cd2+ nor Cu2+ provided enough electron scattering to identify the location of sorption. The affinity series for bacterial removal of these metals decreased in the order Ag greater than La greater than Cu greater than Cd. The results indicate that bacterial cells are capable of binding large quantities of different metals. Adsorption equations may be useful for describing bacterium-metal interactions with metals such as Cd and Cu; however, this approach may not be adequate when precipitation of metals occurs. Images PMID:2515800
Sorption of methylxanthines by different sorbents
NASA Astrophysics Data System (ADS)
Dmitrienko, S. G.; Andreeva, E. Yu.; Tolmacheva, V. V.; Terent'eva, E. A.
2013-05-01
Sorption of caffeine, theophylline, theobromine, diprophylline, and pentoxyphylline on different sorbents (supercross-linked polystyrene, surface-modified copolymer of styrene and divinylbenzene Strata-X, and carbon nanomaterials Taunit and Diasorb-100-C16T) was studied in a static mode in an effort to find new sorbents suitable for sorption isolation and concentration of methylxanthines. The peculiarities of sorption of methylxanthines were explained in relation to the solution acidity, the nature of the sorbates and their concentration, the nature of the solvent, and the structural characteristics of the sorbents.
Mobile scintillometry to study heat advection over heterogeneous surfaces
NASA Astrophysics Data System (ADS)
Kleissl, J.
2007-12-01
Large Aperture Scintillometer (LAS) receivers measure the structure parameter of the refractive index from intensity fluctuations of the transmitter beam. Due to the spatial averaging over 1-4 km employed by this emerging technique the constraints for long temporal averaging (15-30 min) and associated uncertainties that have to be met by other flux measurement techniques do not apply for LASs. In this paper the constraints for temporal averaging of LASs will be examined as a function of environmental conditions and transect geometry. Moreover, analysis of data from a mobile LAS measurement across a surface gradient from rough and dry to smoother and wet will be presented. In this experiment the LAS was mounted on a pickup truck, allowing for quick redeployment of the transect after meaurement. The potential for the use of LAS to study local advection of heat in riparian or irrigated areas in the semi-arid southwest will be evaluated.
Deconvolved spectra of Two Component Advective Flow including jet
NASA Astrophysics Data System (ADS)
Mondal, Santanu; Debnath, Dipak; Chakrabarti, Sandip Kumar
2016-07-01
Outflows and winds are produced when the accretion flows have positive specific energy. Two Component Advective Flow (TCAF) model suggests that the centrifugal pressure supported region of the flow outside the black hole horizon, acts as the base of this outflow. We study the spectral properties of the TCAF which includes a jet component. We consider the jet as a conical in shape which also up-scatters the soft photons from the Keplerian disc. We see that due to the presence of jet component, spectrum become harder as the jet itself behaves like an another Compton cloud above the inner hot corona. We also see how the jet spectra depends on the flow rates. This gives the direct link in timing properties of the X-rays in CENBOL component and the radiation emitted in the jet component.
Carbon dioxide degassing by advective flow from Usu volcano, Japan.
Hernández, P A; Notsu, K; Salazar, J M; Mori, T; Natale, G; Okada, H; Virgili, G; Shimoike, Y; Sato, M; Pérez, N M
2001-04-01
Magmatic carbon dioxide (CO2) degassing has been documented before the 31 March 2000 eruption of Usu volcano, Hokkaido, Japan. Six months before the eruption, an increase in CO2 flux was detected on the summit caldera, from 120 (September 1998) to 340 metric tons per day (September 1999), followed by a sudden decrease to 39 metric tons per day in June 2000, 3 months after the eruption. The change in CO2 flux and seismic observations suggests that before the eruption, advective processes controlled gas migration toward the surface. The decrease in flux after the eruption at the summit caldera could be due to a rapid release of CO2 during the eruption from ascending dacitic dikes spreading away from the magma chamber beneath the caldera. PMID:11292867
Unification of some advection schemes in two dimensions
NASA Technical Reports Server (NTRS)
Sidilkover, D.; Roe, P. L.
1995-01-01
The relationship between two approaches towards construction of genuinely two-dimensional upwind advection schemes is established. One of these approaches is of the control volume type applicable on structured cartesian meshes. It resulted in the compact high resolution schemes capable of maintaining second order accuracy in both homogeneous and inhomogeneous cases. Another one is the fluctuation splitting approach, which is well suited for triangular (and possibly) unstructured meshes. Understanding the relationship between these two approaches allows us to formulate here a new fluctuation splitting high resolution (i.e. possible use of artificial compression, while maintaining positivity property) scheme. This scheme is shown to be linearity preserving in inhomogeneous as well as homogeneous cases.
How Hydrate Saturation Anomalies are Diffusively Constructed and Advectively Smoothed
NASA Astrophysics Data System (ADS)
Rempel, A. W.; Irizarry, J. T.; VanderBeek, B. P.; Handwerger, A. L.
2015-12-01
The physical processes that control the bulk characteristics of hydrate reservoirs are captured reasonably well by long-established model formulations that are rooted in laboratory-verified phase equilibrium parameterizations and field-based estimates of in situ conditions. More detailed assessments of hydrate distribution, especially involving the occurrence of high-saturation hydrate anomalies have been more difficult to obtain. Spatial variations in sediment properties are of central importance for modifying the phase behavior and promoting focussed fluid flow. However, quantitative predictions of hydrate anomaly development cannot be made rigorously without also addressing the changes in phase behavior and mechanical balances that accompany changes in hydrate saturation level. We demonstrate how pore-scale geometrical controls on hydrate phase stability can be parameterized for incorporation in simulations of hydrate anomaly development along dipping coarse-grained layers embedded in a more fine-grained background that is less amenable to fluid transport. Model simulations demonstrate how hydrate anomaly growth along coarse-layer boundaries is promoted by diffusive gas transport from the adjacent fine-grained matrix, while advective transport favors more distributed growth within the coarse-grained material and so effectively limits the difference between saturation peaks and background levels. Further analysis demonstrates how sediment contacts are unloaded once hydrate saturation reaches sufficient levels to form a load-bearing skeleton that can evolve to produce segregated nodules and lenses. Decomposition of such growth forms poses a significant geohazard that is expected to be particularly sensitive to perturbations induced by gas extraction. The figure illustrates the predicted evolution of hydrate saturation Sh in a coarse-grained dipping layer showing how prominent bounding hydrate anomalies (spikes) supplied by diffusive gas transport at early times
Frictional drag reduction by wavy advection of deformable bubbles
NASA Astrophysics Data System (ADS)
Oishi, Yoshihiko; Murai, Yuichi; Tasaka, Yuji; Yasushi, Takeda
2009-02-01
Bubbles can reduce frictional drag in wall turbulence, and its effect is expected to use for ships and pipelines to save their power consumptions. A number of basic experiments have been carried out to date for finding out the best condition for enhancing the drag reduction. One issue that remains at present is the difference of the performance between steady and unsteady status in terms of bubble concentration. All the experiments in the past deal with the steady effect, i.e., the drag reduction is evaluated as a function of mean void fraction or given gas flow rate of continuous injection. Despite to this, the actual phenomena highly depend on local interaction between two phases upon unsteady manner. We focus on this point and elucidate the influence of time-fluctuating void fraction on the total response to the drag reduction. This view is in fact important to estimate the persistency of the bubble-based drag reduction in the flow direction since bubbles formulate wavy advection during their migration. Our experiments are designed to measure the above-mentioned effect from laminar, transitional, and turbulent flows in a horizontal channel. For avoiding the contamination effect that worsens the reproducibility of the experiment, Silicone oil is used as carrier fluid. The oil also simulates the high Weber number bubble condition because of low surface tension. The unsteady interaction between the wavy advection of bubbles and the local skin friction, a synchronized system is constructed to connect the high-speed camera with the shear transducer, which can evaluate the interaction at 1000 fps. From the results, we confirm that the drag reduction is provided at Re>3000 in the turbulent flow regime, and also the total drag reduction is enhanced by the presence of the waves.
Multirate Runge-Kutta schemes for advection equations
NASA Astrophysics Data System (ADS)
Schlegel, Martin; Knoth, Oswald; Arnold, Martin; Wolke, Ralf
2009-04-01
Explicit time integration methods can be employed to simulate a broad spectrum of physical phenomena. The wide range of scales encountered lead to the problem that the fastest cell of the simulation dictates the global time step. Multirate time integration methods can be employed to alter the time step locally so that slower components take longer and fewer time steps, resulting in a moderate to substantial reduction of the computational cost, depending on the scenario to simulate [S. Osher, R. Sanders, Numerical approximations to nonlinear conservation laws with locally varying time and space grids, Math. Comput. 41 (1983) 321-336; H. Tang, G. Warnecke, A class of high resolution schemes for hyperbolic conservation laws and convection-diffusion equations with varying time and pace grids, SIAM J. Sci. Comput. 26 (4) (2005) 1415-1431; E. Constantinescu, A. Sandu, Multirate timestepping methods for hyperbolic conservation laws, SIAM J. Sci. Comput. 33 (3) (2007) 239-278]. In air pollution modeling the advection part is usually integrated explicitly in time, where the time step is constrained by a locally varying Courant-Friedrichs-Lewy (CFL) number. Multirate schemes are a useful tool to decouple different physical regions so that this constraint becomes a local instead of a global restriction. Therefore it is of major interest to apply multirate schemes to the advection equation. We introduce a generic recursive multirate Runge-Kutta scheme that can be easily adapted to an arbitrary number of refinement levels. It preserves the linear invariants of the system and is of third order accuracy when applied to certain explicit Runge-Kutta methods as base method.
Examination of the evolution of radiation and advection fogs. Final report
Orgill, M.M.
1993-01-01
A literature study was done on radiation and advection fog evolution. For radiation fog, six stages of fog evolution have been identified -- (1) precursor, (2) sunset, (3) conditioning, (4) mature, (5) sunrise, and (6) dissipation. The evolution of advection fog models has been in parallel with radiation fog models, but no identified stages in the evolution of advection fog have been proposed: (1) precursor, (2) initiation, (3) mature, and (4) dissipation. Radiation and advection fog models will require greater sophistication in order to study fog spatial and temporal variability. Physical aspects that require further study are discussed.
The role of phase dynamics in a stochastic model of a passively advected scalar
NASA Astrophysics Data System (ADS)
Moradi, Sara; Anderson, Johan
2016-05-01
Collective synchronous motion of the phases is introduced in a model for the stochastic passive advection-diffusion of a scalar with external forcing. The model for the phase coupling dynamics follows the well known Kuramoto model paradigm of limit-cycle oscillators. The natural frequencies in the Kuramoto model are assumed to obey a given scale dependence through a dispersion relation of the drift-wave form -βk/1 +k2 , where β is a constant representing the typical strength of the gradient. The present aim is to study the importance of collective phase dynamics on the characteristic time evolution of the fluctuation energy and the formation of coherent structures. Our results show that the assumption of a fully stochastic phase state of turbulence is more relevant for high values of β, where we find that the energy spectrum follows a k-7 /2 scaling. Whereas for lower β there is a significant difference between a-synchronised and synchronised phase states, one could expect the formation of coherent modulations in the latter case.
Shear dispersion and turbulence decorrelation by differential rotation
Garcia, O.E.; Bian, N.H.
2005-01-01
The shear enhanced dispersion of a passive scalar field subject to differential rotation is investigated analytically and interpretations are given in terms of turbulence shear decorrelation. Using the method of advected coordinates, the enhanced dispersion caused by steady and oscillatory flows with uniform shear is derived and the well-known turbulence shear decorrelation theory is recovered. The additional role of kinetic energy transfer due to differential advection of vorticity is also pointed out. Finally, the shear enhanced dispersion due to flows with periodic variations in space as well as time is given. It is found that radially alternating flows may significantly reduce the turbulence decorrelation time provided the root mean square flow shear is larger than the flow oscillation frequency. In the opposite limit of fast flow oscillations there is no turbulence decorrelation.
COSOLVENT EFFECTS ON SORPTION ISOTHERM LINEARITY
Sorption-desorption hysteresis, slow desorption kinetics, and other nonideal phenomena have been attributed to the differing sorptive characteristics of the natural organic polymers associated with soils and sediments. In this study, aqueous and mixed solvent systems were used t...
SORPTION OF HYDROPHOBIC POLLUTANTS ON NATURAL SEDIMENTS
The sorption of hydrophobic compounds (aromatic hydrocarbons and chlorinated hydrocarbons) spanning a concentration range in water solubility from 500 parts per trillion (ppt) to 1800 parts per million (ppm) on local (North Georgia) pond and river sediments was investigated. The ...
Advection and Diffusion of Substances in Biological Tissues With Complex Vascular Networks
Beard, Daniel A.; Bassingthwaighte, James B.
2010-01-01
For highly diffusive solutes the kinetics of blood–tissue exchange is only poorly represented by a model consisting of sets of independent parallel capillary–tissue units. We constructed a more realistic multicapillary network model conforming statistically to morphometric data. Flows through the tortuous paths in the network were calculated based on constant resistance per unit length throughout the network and the resulting advective intracapillary velocity field was used as a framework for describing the extravascular diffusion of a substance for which there is no barrier or permeability limitation. Simulated impulse responses from the system, analogous to tracer water outflow dilution curves, showed flow-limited behavior over a range of flows from about 2 to 5 ml min−1 g−1, as is observed for water in the heart in vivo. The present model serves as a reference standard against which to evaluate computationally simpler, less physically realistic models. The simulated outflow curves from the network model, like experimental water curves, were matched to outflow curves from the commonly used axially distributed models only by setting the capillary wall permeability–surface area (PS) to a value so artifactually low that it is incompatible with the experimental observations that transport is flow limited. However, simple axially distributed models with appropriately high PSs will fit water outflow dilution curves if axial diffusion coefficients are set at high enough values to account for enhanced dispersion due to the complex geometry of the capillary network. Without incorporating this enhanced dispersion, when applied to experimental curves over a range of flows, the simpler models give a false inference that there is recruitment of capillary surface area with increasing flow. Thus distributed models must account for diffusional as well as permeation processes to provide physiologically appropriate parameter estimates. PMID:10784090
Reversibility of strontium sorption on fracture fillings
Cui, D.; Eriksen, T.E.
1995-12-31
Granite has been chosen by several countries as a major candidate for deep geologic disposal of radioactive waste. The authors have carried out a comparative study of sorption and desorption of strontium in groundwater on separated size and magnetic fractions of fracture fillings from deep granite. Complete reversibility of the sorption process was demonstrated by identical Freundlich isotherms, isotopic exchangeability and pH dependence of the distribution coefficient R{sub d}.
NASA Astrophysics Data System (ADS)
Yoon, H.; Lee, J.; An, J.
2013-12-01
In recent times, iodate (IO3-) which can be generated under highly oxidized conditions such as the ozonation process in a water treatment plant has been receiving increasing attention due to its high toxicity to human and environment. In this respect, sorption behavior of iodide (I-) and IO3- on an activated carbon (AC) disk as a solid sorbent was investigated for the further development of efficient removal and analytical techniques. To this end, batch sorption tests were performed as a function of pH, sample volume, and initial concentration. Sorption of I- occurred preferentially on the surface of AC disk, regardless of pH levels (i.e., 4, 6, and 8). However, IO3- was quite sensitive to pH levels and the sorption capability of IO3- on the AC disk was much smaller than that of I-. Maximum sorption capacities of I- and IO3- in the different matrices (i.e., deionized water and seawater) were also assessed. In addition, the analytical strategy for the iodine speciation was also introduced. Iodine species were separated and pre-concentrated onto the AC disk based on their selective sorption properties according to the pH levels. Then, the AC disk pre-concentrated was directly analyzed using wavelength dispersive X-ray fluorescence spectrometry. Acknowledgement This research was supported by a grant from the Korea Basic Science Institute (project No. E33300).
Ali, Norizan; El-Harbawi, Mohanad; Jabal, Ayman Abo; Yin, Chun-Yang
2012-01-01
The characteristics and water/oil sorption effectiveness ofkapok fibre, sugarcane bagasse and rice husks have been compared. The three biomass types were subjected to field emission scanning electron microscopy-energy dispersive X-ray spectroscopy and surface tension analyses for liquid-air and oil-water systems were conducted. Both kapok fibre and sugarcane bagasse exhibit excellent oil sorption capabilities for diesel, crude, new engine and used engine oils as their oil sorption capacities all exceed 10 g/g. The synthetic sorbent exhibits oil sorption capacities comparable with sugarcane bagasse, while rice husks exhibit the lowest oil sorption capacities among all the sorbents. Kapok fibre shows overwhelmingly high oil-to-water sorption (O/W) ratios ranging from 19.35 to 201.53 while sugarcane bagasse, rice husks and synthetic sorbent have significantly lower O/W ratios (0.76-2.69). This suggests that kapok fibre is a highly effective oil sorbent even in well-mixed oil-water media. An oil sorbent suitability matrix is proposed to aid stakeholders in evaluating customized oil removal usage of the natural sorbents. PMID:22629620
Looney, B.B.; Scott, M.T.
1988-12-31
Recent field and laboratory data have confirmed that apparent dispersivity is a function of the flow distance of the measurement. This scale effect is not consistent with classical advection dispersion modeling often used to describe the transport of solutes in saturated porous media. Many investigators attribute this anomalous behavior to the fact that the spreading of solute is actually the result of the heterogeneity of subsurface materials and the wide distribution of flow paths and velocities available in such systems. An analysis using straightforward analytical equations confirms this hypothesis. An analytical equation based on a flow variance approach matches available field data when a variance description of approximately 0.4 is employed. Also, current field data provide a basis for statistical selection of the variance parameter based on the level of concern related to the resulting calculated concentration. While the advection dispersion approach often yielded reasonable predictions, continued development of statistical and stochastic techniques will provide more defendable and mechanistically descriptive models.
Antonov, N V; Gulitskiy, N M
2015-10-01
In this work we study the generalization of the problem considered in [Phys. Rev. E 91, 013002 (2015)] to the case of finite correlation time of the environment (velocity) field. The model describes a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial-range asymptotic behavior is studied by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and preassigned pair correlation function. Due to the presence of distinguished direction n, all the multiloop diagrams in this model vanish, so that the results obtained are exact. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their stability depends on the relation between the exponents in the energy spectrum E∝k(⊥)(1-ξ) and the dispersion law ω∝k(⊥)(2-η). In contrast to the well-known isotropic Kraichnan's model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence scale L. PMID:26565343
NASA Astrophysics Data System (ADS)
Antonov, N. V.; Gulitskiy, N. M.
2015-10-01
In this work we study the generalization of the problem considered in [Phys. Rev. E 91, 013002 (2015), 10.1103/PhysRevE.91.013002] to the case of finite correlation time of the environment (velocity) field. The model describes a vector (e.g., magnetic) field, passively advected by a strongly anisotropic turbulent flow. Inertial-range asymptotic behavior is studied by means of the field theoretic renormalization group and the operator product expansion. The advecting velocity field is Gaussian, with finite correlation time and preassigned pair correlation function. Due to the presence of distinguished direction n , all the multiloop diagrams in this model vanish, so that the results obtained are exact. The inertial-range behavior of the model is described by two regimes (the limits of vanishing or infinite correlation time) that correspond to the two nontrivial fixed points of the RG equations. Their stability depends on the relation between the exponents in the energy spectrum E ∝k⊥1 -ξ and the dispersion law ω ∝k⊥2 -η . In contrast to the well-known isotropic Kraichnan's model, where various correlation functions exhibit anomalous scaling behavior with infinite sets of anomalous exponents, here the corrections to ordinary scaling are polynomials of logarithms of the integral turbulence scale L .
Estimation of soil sorption coefficients using QSARs
Doucette, W.J.
1994-12-31
Sorption coefficients quantitatively describe the extent to which an organic chemical distributes itself between an environmental solid and the aqueous phase that it is contact with at equilibrium. Because of the difficulty and expense associated with measuring sorption coefficients, estimated values are often used in place of site specific, experimental values for fate modeling. Most reported methods for estimating the sorption of organic chemicals onto environmental solids are based on observation that for many organic chemicals, and in particular neutral hydrophobic organics, sorption is directly proportional to the quantity of organic matter associated with the solid. Normalizing soil or sediment specific sorption coefficients to the organic carbon content of the sorbent yields a new ``constant``, Koc, that is considered unique property of the organic chemical being sorbed. Values of Koc are then typically estimated from correlations between Koc and various descriptors of hydrophobicity such as octanol/water partition coefficients (Kow), aqueous solubility (S), molecular connectivity indices (MCIs) and retention times or capacity factors generated by reverse phase high performance liquid chromatography (RP-HPLC). Group contribution methods have also been described. While the so-called ``Koc approach`` for estimating sorption coefficients is most appropriate for neutral, hydrophobic organic chemicals on environmental solids containing a significant amount of organic matter, it has been applied to a wide variety of chemical and soil types. This presentation will focus on a discussion of the Koc approach, its applicability and limitations. A comparison of several widely used methods for estimating Koc will be presented.
Evaluating phenanthrene sorption on various wood chars
James, G.; Sabatini, D.A.; Chiou, C.T.; Rutherford, D.; Scott, A.C.; Karapanagioti, H.K.
2005-01-01
A certain amount of wood char or soot in a soil or sediment sample may cause the sorption of organic compounds to deviate significantly from the linear partitioning commonly observed with soil organic matter (SOM). Laboratory produced and field wood chars have been obtained and analyzed for their sorption isotherms of a model solute (phenanthrene) from water solution. The uptake capacities and nonlinear sorption effects with the laboratory wood chars are similar to those with the field wood chars. For phenanthrene aqueous concentrations of 1 ??gl-1, the organic carbon-normalized sorption coefficients (log Koc) ranging from 5.0 to 6.4 for field chars and 5.4-7.3 for laboratory wood chars, which is consistent with literature values (5.6-7.1). Data with artificial chars suggest that the variation in sorption potential can be attributed to heating temperature and starting material, and both the quantity and heterogeneity of surface-area impacts the sorption capacity. These results thus help to corroborate and explain the range of log Koc values reported in previous research for aquifer materials containing wood chars. ?? 2004 Elsevier Ltd. All rights reserved.
NASA Astrophysics Data System (ADS)
Hansel, Colleen M.; Benner, Shawn G.; Neiss, Jim; Dohnalkova, Alice; Kukkadapu, Ravi K.; Fendorf, Scott
2003-08-01
Iron (hydr)oxides not only serve as potent sorbents and repositories for nutrients and contaminants but also provide a terminal electron acceptor for microbial respiration. The microbial reduction of Fe (hydr)oxides and the subsequent secondary solid-phase transformations will, therefore, have a profound influence on the biogeochemical cycling of Fe as well as associated metals. Here we elucidate the pathways and mechanisms of secondary mineralization during dissimilatory iron reduction by a common iron-reducing bacterium, Shewanella putrefaciens (strain CN32), of 2-line ferrihydrite under advective flow conditions. Secondary mineralization of ferrihydrite occurs via a coupled, biotic-abiotic pathway primarily resulting in the production of magnetite and goethite with minor amounts of green rust. Operating mineralization pathways are driven by competing abiotic reactions of bacterially generated ferrous iron with the ferrihydrite surface. Subsequent to the initial sorption of ferrous iron on ferrihydrite, goethite (via dissolution/reprecipitation) and/or magnetite (via solid-state conversion) precipitation ensues resulting in the spatial coupling of both goethite and magnetite with the ferrihydrite surface. The distribution of goethite and magnetite within the column is dictated, in large part, by flow-induced ferrous Fe profiles. While goethite precipitation occurs over a large Fe(II) concentration range, magnetite accumulation is only observed at concentrations exceeding 0.3 mmol/L (equivalent to 0.5 mmol Fe[II]/g ferrihydrite) following 16 d of reaction. Consequently, transport-regulated ferrous Fe profiles result in a progression of magnetite levels downgradient within the column. Declining microbial reduction over time results in lower Fe(II) concentrations and a subsequent shift in magnetite precipitation mechanisms from nucleation to crystal growth. While the initial precipitation rate of goethite exceeds that of magnetite, continued growth is inhibited by