共查询到20条相似文献,搜索用时 31 毫秒
1.
《Atmospheric environment(England)》1978,12(6-7):1365-1369
A numerical model is described for computing pollutant concentration distributions downwind from a source. It is based on the three-dimensional dispersion equation governing the time-dependent advective and diffusive transport of air pollutants and is solved numerically by a mixed Lagrangian-Eulerian finite-difference scheme. The model includes the vertical wind shear, the turning of the actual wind, and vertical variations of the vertical eddy diffusivity. In this paper the model is used to simulate the pollutant dispersion process in a stratified planetary boundary layer. The vertical profiles of horizontal mean wind and vertical eddy diffusivities are calculated numerically from a planetary boundary layer model. The influence of the ground roughness and the atmospheric stability on the pollutant distribution is investigated. The results indicate that both parameters essentially determine the air pollutant dispersion process. 相似文献
2.
Huali Yuan John C. Little Eva Marand Zhe Liu 《Atmospheric environment (Oxford, England : 1994)》2007,41(40):9300-9308
Structural insulated panels (SIPs) have significant environmental and energy advantages. However, the tight structure that results may cause degraded indoor air quality and the potential release of volatile organic compounds (VOCs) from these layered materials must be considered. A physically based model for predicting VOC emissions from multi-layer materials is described. Fugacity is used to eliminate the concentration discontinuities at the interface between layers. This avoids an obstacle associated with numerically simulating mass transfer in composite materials. The numerical model is verified for a double-layer system by comparing predicted concentrations to those obtained with a previously published analytical model. In addition, hexanal emissions from multi-layer SIPs are simulated to demonstrate the usefulness of the fugacity approach. Finally, the multi-layer model is used to investigate the impact that clay/polyurethane nanocomposite diffusion barriers can have on VOC emissions. Indoor gas-phase concentrations can be greatly reduced with a barrier layer on the surface, thereby minimizing the environmental impact of SIPs. 相似文献
3.
The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions. 相似文献
4.
Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases. 相似文献
5.
The paper presents a model for contaminant transport and flux through a consolidating subaqueous sediment and overlying cap. The formulation is based on the effect of consolidation and excess pore pressure dissipation on transient, nonlinear advective component of transport through sediment and the cap. The consolidation is induced by the buoyant weight of the cap when it is placed on the contaminated sediments. One equation is presented for advective-diffusive transport through the sediment that is dependent upon soil/contaminant properties and transient advective velocity, which is calculated from a second equation based on the Terzaghi consolidation theory. A third equation is provided to describe the transport of contaminants in the cap. The parameters, including advective velocity, and boundary conditions used for contaminant transport through the cap are derived from the solution of the first two equations. The finite difference method is used to solve the system of equations for consolidation and contaminant transport. A hypothetical case is analyzed to demonstrate the formulation, and the results show that advection due to consolidation can accelerate breakthrough of contaminant through the cap by orders of magnitude. The derivation and results show that consolidation should be included for cap design, and that reactive caps are essential for delaying and reducing dissolved contaminant flux. 相似文献
6.
An analytical solution for vertical transport of volatile chemicals in the vadose zone 总被引:1,自引:0,他引:1
An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability. 相似文献
7.
A two-dimensional flow and transport model was developed for simulating transient water flow and nonreactive solute transport in heterogeneous, unsaturated porous media containing air and water. The model is composed of a unique combination of robust and accurate numerical algorithms for solving the Richards', Darcy flux, and advection-dispersion equations. The mixed form of Richards' equation is solved using a finite-element formulation and a modified Picard iteration scheme. Mass lumping is employed to improve solution convergence and stability behavior. The flow algorithm accounts for hysteresis in the pressure head-water content relationship. Darcy fluxes are approximated with a Galerkin and Petrov-Galerkin finite-element method developed for random heterogeneous porous media. The transport equation is solved using an Eulerian-Lagrangian method. A multi-step, fourth-order Runge-Kutta, reverse particle tracking technique and a quadratic-linear interpolation scheme are shown to be superior for determining the advective concentration. A Galerkin finite-element method is used for approximating the dispersive flux. The unsaturated flow and transport model was applied to a variety of rigorous problems and was found to produce accurate, mass-conserving solutions when compared to analytical solutions and published numerical results. 相似文献
8.
This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials. 相似文献
9.
To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions. 相似文献
10.
Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils. 相似文献
11.
《Atmospheric environment(England)》1987,21(1):201-212
A three-dimensional, grid-based numerical air pollution model for the estimation of air pollutant concentrations in an urban area is developed. Based on the continuity equation, the modeling system incorporates the combined influences of advective transport, turbulent diffusion, chemical transformation, source emissions and surface removal of air contaminants. Recent developments in plume rise and plume penetration processes, objective wind field analysis procedures and numerical solution techniques incorporated into the model are described. 相似文献
12.
《Atmospheric environment (Oxford, England : 1994)》2001,35(28):4777-4790
Sulfur dioxide absorption dynamics in a raindrop are studied numerically by means of a fully numerical simulation method (FNSM) in which a composite orthogonal grid system consisting of both gas- and liquid-phase is adopted. When a raindrop with fixed radius falls in association with terminal velocity, a recirculation bubble always accompanies the gas-phase flow field in the aft region of the drop. With regard to the drop internal flow structure it has a drastic variation with drop size. When the drop radius is small (e.g. rs=200 μm), only a single vortex motion is seen inside the drop. Under such a situation, sulfur dioxide absorbed from the interface is mainly transported from the area in front of the aft stagnation point. In contrast, as the drop is relatively large, say, rs=500 μm, it is of interest to find that a double-vortex motion, composed of a primary and a secondary vortexes, is clearly observed. As a result, the onset of SO2 transport process occurs at the drop’s surface near where the two vortexes meet. By defining a drop mass transport number, it indicates that the mass transported via internal circulation is always much faster than that by mass diffusion and the latter is highly relevant to the drop radius. Accordingly, the SO2 transport dynamics in a raindrop is essentially determined by drop size. 相似文献
13.
David Werner Hrissi K. Karapanagioti Patrick Hhener 《Journal of contaminant hydrology》2009,108(1-2):54-63
A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers. 相似文献
14.
Josephine A. Jackson Jude T. Sommerfeld 《Journal of the Air & Waste Management Association (1995)》2013,63(11):1328-1243
The DSS/2 software system was used to numerically simulate oxygen stratification effects in a typical industrial boiler stack. The model developed for the simulation involves turbulent mass transfer in a tube with a constant, high mass flux from the walls of the tube. The concentration profile is considered to be developing, but the velocity profile is assumed to be fully developed. The model is for a gaseous system at constant temperature and pressure; the transport properties are shown to be constant. The universal velocity profile for turbulent flow and von Karman’s expressions for the turbulent eddy diffusivity are incorporated in the model. An expression for a wall concentration profile corresponding to conditions of high mass flux at the wall is developed. The model is then expressed as a system of partial differential equations which are solved using the numerical method of lines together with the DSS/2 numerical integration system. The resulting concentration profiles are presented for various flux rates which might be encountered in a typical stack. These profiles indicate that a potential for oxygen stratification in stacks does exist; they also suggest that careful placement of the oxygen cell probe can help minimize such errors. 相似文献
15.
A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required. 相似文献
16.
High performance computing has made possible the development of high resolution, multidimensional, multicomponent reactive transport models that can be used to analyze complex geochemical environments. However, as increasingly complex processes are included in these models, the accuracy of the numerical formulation coupling the nonlinear processes becomes difficult to verify. Analytical solutions are not available for realistically complex problems and benchmark solutions are not generally available for specific problems. We present an advective reactive streamtube (ARS) transport technique that efficiently provides accurate solutions of nonlinear multicomponent reactive transport in nonuniform multidimensional velocity fields. These solutions can be compared with results from Eulerian-based advection-dispersion-reaction models to evaluate the accuracy of the numerical formulation used. The ARS technique includes mixed equilibrium and kinetic complexation and precipitation-dissolution reactions subject to the following assumptions: (1) transport is purely advective (i.e., no explicit diffusion or dispersion), and (2) chemistry is described by a canonical system of reactions that evolves with time and is unaffected by position in space. Results from the ARS technique are compared with results from the massively parallel, multicomponent reactive transport model MCTRACKER on a test problem involving irreversible oxidation of organic carbon and reaction of the oxidation products with two immobile mineral phases, gypsum and calcite, and fifteen aqueous complexes. Truncation error, operator splitting error, and the nonlinear transformation of these errors in the high-resolution reactive transport model are identified for this problem. 相似文献
17.
Kosakowski G 《Journal of contaminant hydrology》2004,72(1-4):23-46
Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters. 相似文献
18.
19.
《Atmospheric environment(England)》1984,18(10):2109-2124
The mathematical model presented in this paper describes in detail the gas-phase chemistry (80 reactions), gas-phase/liquid-phase equilibrium (12 equilibria) and liquid-phase chemistry (27 reactions and equilibria) of acid formation. The model is used to simulate acid formation under various conditions:
- 1.(1) nonprecipitating clouds in the Adirondacks,
- 2.(2) raining clouds in the Ohio River valley,
- 3.(3) night-time fog in the Los Angeles basin and
- 4.(4) night-time nonprecipitating clouds in the Los Angeles basin. Model calculations are compared with observed values of sulfate and nitrate formation. The model is used to analyze the relative importance of sulfate and nitrate formation pathways under these various conditions.
20.
Numerical simulation tools have been used to study the dominating processes during transport of aromatic hydrocarbons in the unsaturated soil zone. Simulations were based on field observations at an experimental site located on a glacial delta plain with pronounced layered sedimentary structures. A numerical model for transport in the unsaturated zone, SWMS-3D, has been extended to incorporate coupled multispecies transport, microbial degradation following Monod kinetics and gas diffusive transport of oxygen and hydrocarbons. The flow field parameters were derived from previous work using nonreactive tracers. Breakthrough curves (BTC) from the hydrocarbon field experiment were used to determine sorption parameters and Monod kinetic parameters using a fitting procedure. The numerical simulations revealed that the assumption of homogeneous layers resulted in deviations from the field observations. The deviations were more pronounced with incorporation of reactive transport, compared with earlier work on nonreactive transport. To be able to model reasonable BTC, sorption had to be reduced compared to laboratory experiments. The initial biomass and the maximum utilisation rate could be adjusted to capture both the initial lag phase and the overall degradation rate. Nevertheless, local oxygen limitation is predicted by the model, which was not observed in the field experiment. Incorporation of evaporation and diffusive gas transport of the hydrocarbons did not significantly change the local oxygen demand. The main cause of the observed discrepancies between model and field are attributed to channelling as a result of small-scale heterogeneities such as biopores. 相似文献