共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution. 相似文献
2.
A semi-analytical method has been developed for calculating oil recovery in two and three dimensions, and for calculating effective relative permeabilities for coarse grids. The calculations are based on the assumption that the effects of a changing mobility field can be accounted for by using fixed streamtube geometries with flowrates updated to account for the changing mobility distribution. The single-phase pressure distribution from a numerical solution of Laplace's equation is used to calculate the pressure distribution for a two-phase flow based on a mapping of the solution of the Buckley-Leverett equation onto the streamtubes derived from the single-phase solution. The displacement calculations for oil recovery are based on theory previously developed by Dykstra and Parsons, extended to include the effects of spatially varying permeability and continuously changing mobilities, as occurs in solutions of the Buckley-Leverett equation for typical values of the mobility ratio. This idea has also been extended to the calculation of effective relative permeabilities for coarse-grid simulation and finally establishes the proper rules for averaging the results of fine-grid numerical simulations of two-phase flow for the definition of effective two-phase flow properties on coarse grids. These calculations have been generalized to three-dimensional flows by the simple device of conceptually inserting a gridded plane across the flow and defining each streamtube at that location as those streamlines which pass through any one of the grid cells. When combined with time-of-flight calculations from the gridded plane to both the producer and injector, the distribution of pore volume along each streamtube can be calculated. This information, combined with a tabulation of the single-phase, steady-state pressure distribution along each streamtube, provides all of the information needed for the semi-analytical calculation of oil recovery and effective flow properties in three-dimensional flows. © 1997 Elsevier Science Ltd. All rights reserved 相似文献
3.
We developed a new semi-analytical source zone depletion model (SZDM) for multicomponent light nonaqueous phase liquids (LNAPLs) and incorporated this into an existing screening model for estimating cleanup times for chemical spills from railroad tank cars that previously considered only single-component LNAPLs. Results from the SZDM compare favorably to those from a three-dimensional numerical model, and from another semi-analytical model that does not consider source zone depletion. The model was used to evaluate groundwater contamination and cleanup times for four complex mixtures of concern in the railroad industry. Among the petroleum hydrocarbon mixtures considered, the cleanup time of diesel fuel was much longer than E95, gasoline, and crude oil. This is mainly due to the high fraction of low solubility components in diesel fuel. The results demonstrate that the updated screening model with the newly developed SZDM is computationally efficient, and provides valuable comparisons of cleanup times that can be used in assessing the health and financial risk associated with chemical mixture spills from railroad-tank-car accidents. 相似文献
4.
Suthan Suthersan John Horst Mark Klemmer Don Malone 《Ground Water Monitoring & Remediation》2012,32(3):34-40
The remediation industry has witnessed multiple innovations arising from a greater understanding of the physical, chemical, and biological processes that control the fate and transport of chemicals in the subsurface environment. In addition, increasing emphasis is being placed on remediation solutions that are greener, simpler, and more resource efficient. The positive impacts that can be derived from this emphasis include reduced energy consumption, reduced waste emissions, and lower costs. Temperature‐activated auto‐decomposition reactions represent a potentially underutilized option for the in situ remediation of certain organic contaminants, and an option that can be both highly effective and greener than other available technologies. 相似文献
5.
《Advances in water resources》1998,21(3):237-249
The transport and fate of reactive chemicals in groundwater is governed by equations which are often difficult to solve due to the nonlinear relationship between the solute concentrations for the liquid and solid phases. The nonlinearity may cause mass balance errors during the numerical simulation in addition to numerical errors for linear transport system. We have generalized the modified Picard iteration algorithm of Celia et al.5 for unsaturated flow to solve the nonlinear transport equation. Written in a ‘mixed-form’ formulation, the total solute concentration is expanded in a Taylor series with respect to the solution concentration to linearize the transport equation, which is then solved with a conventional finite element method. Numerical results of this mixed-form algorithm are compared with those obtained with the concentration-based scheme using conventional Picard iteration. In general, the new solver resulted in negligible mass balance errors (< ∥10−8∥%) and required less computational time than the conventional iteration scheme for the test examples, including transport involving highly nonlinear adsorption under steady-state as well as transient flow conditions. In contrast, mass balance errors resulting from the conventional Picard iteration method were higher than 10% for some highly nonlinear problems. Application of the modified Picard iteration scheme to solve the nonlinear transport equation may greatly reduce the mass balance errors and increase computational efficiency. 相似文献
6.
An empirical hyperbolic scale-dependent dispersion model, which predicts a linear growth of dispersivity close to the origin and the attainment of an asymptotic dispersivity at large distances, is presented for deterministic modelling of field-scale solute transport and the analysis of solute transport experiments. A simple relationship is derived between local dispersivity, which is used in numerical simulations of solute transport, and effective dispersivity, which is estimated from the analysis of tracer breakthrough curves. The scale-dependent dispersion model is used to interpret a field tracer experiment by nonlinear least-squares inversion of a numerical solution for unsaturated transport. Simultaneous inversion of concentration-time data from several sampling locations indicates a linear growth of the dispersion process over the scale of the experiment. These findings are consistent with the results of an earlier analysis based on the use of a constant dispersion coefficient model at each of the sampling depths. 相似文献
7.
M. Oostrom M.J. Truex A.K. Rice C.D. Johnson K.C. Carroll D.J. Becker M.A. Simon 《Ground Water Monitoring & Remediation》2014,34(2):71-84
Soil vapor extraction (SVE) is a prevalent remediation remedy for volatile organic compound (VOC) contaminants in the vadose zone. To support selection of an appropriate condition at which SVE may be terminated for site closure or for transition to another remedy, an evaluation is needed to determine whether vadose zone VOC contamination has been diminished sufficiently to keep groundwater concentrations below threshold values. A conceptual model for this evaluation was developed for VOC fate and transport from a vadose zone source to groundwater when vapor‐phase diffusive transport is the dominant transport process. A numerical analysis showed that, for these conditions, the groundwater concentration is controlled by a limited set of parameters, including site‐specific dimensions, vadose zone properties, and source characteristics. On the basis of these findings, a procedure was then developed for estimating groundwater concentrations using results from the three‐dimensional multiphase transport simulations for a matrix of parameter value combinations and covering a range of potential site conditions. Interpolation and scaling processes are applied to estimate groundwater concentrations at compliance (monitoring) wells for specific site conditions of interest using the data from the simulation results. The interpolation and scaling methodology using these simulation results provides a far less computationally intensive alternative to site‐specific three‐dimensional multiphase site modeling, while still allowing for parameter sensitivity and uncertainty analyses. With iterative application, the approach can be used to consider the effect of a diminishing vadose zone source over time on future groundwater concentrations. This novel approach and related simulation results have been incorporated into a user‐friendly Microsoft® Excel®‐based spreadsheet tool entitled SVEET (Soil Vapor Extraction Endstate Tool), which has been made available to the public. 相似文献
8.
An approximate analytical solution to the advection-dispersion equation was derived to describe solute transport during spherical-flow conditions in single-well push-pull tests. The spherical-flow case may be applicable to aquifer tests conducted in packed intervals or partially penetrating wells. Using results of two-dimensional numerical simulations, we briefly illustrate the applicability of the derived spherical-flow solution and provide a comparison with its cylindrical-flow counterpart. Good agreement between simulated extraction-phase breakthrough curves and the spherical-flow solution was found when the length of the injection/extraction region was small compared to both aquifer thickness and maximum solute frontal position at the end of the injection phase. On the other hand, discrepancies between simulated breakthrough curves and the spherical-flow solution increased with increasing anisotropy in hydraulic conductivities. Several inherent limitations embedded in its derivation such as assumptions of isotropy and homogeneity warrant the cautious use of the spherical-flow solution. 相似文献
9.
In this paper we describe the transport of pollution in groundwater in the neighbourhood of a well in a uniform background flow. We compute the rate at which contaminated particles reach the well as a function of the place of the source of pollution. The motion of a particle in a dispersive flow is seen as a random walk process. The Fokker-Planck equation for the random motion of a particle is transformed using the complex potential for the advective flow field. The resulting equation is solved asymptotically after a stretching transformation. Finally, the analytical solution is compared with results from Monte Carlo simulations with the random walk model. The method can be extended to arbitrary flow fields. Then by a numerical coordinate transformation the analytical results can still be employed. 相似文献
10.
When managing large-scale ground water contamination problems, it is often necessary to model flow and transport using finely discretized domains--for instance (1) to simulate flow and transport near a contamination source area or in the area where a remediation technology is being implemented; (2) to account for small-scale heterogeneities; (3) to represent ground water-surface water interactions; or (4) some combination of these scenarios. A model with a large domain and fine-grid resolution will need extensive computing resources. In this work, a domain decomposition-based assembly model implemented in a parallel computing environment is developed, which will allow efficient simulation of large-scale ground water flow and transport problems using domain-wide grid refinement. The method employs common ground water flow (MODFLOW) and transport (RT3D) simulators, enabling the solution of almost all commonly encountered ground water flow and transport problems. The basic approach partitions a large model domain into any number of subdomains. Parallel processors are used to solve the model equations within each subdomain. Schwarz iteration is applied to match the flow solution at the subdomain boundaries. For the transport model, an extended numerical array is implemented to permit the exchange of dispersive and advective flux information across subdomain boundaries. The model is verified using a conventional single-domain model. Model simulations demonstrate that the proposed model operated in a parallel computing environment can result in considerable savings in computer run times (between 50% and 80%) compared with conventional modeling approaches and may be used to simulate grid discretizations that were formerly intractable. 相似文献
11.
Baha Y. Mirghani Kumar G. Mahinthakumar Michael E. Tryby Ranji S. Ranjithan Emily M. Zechman 《Advances in water resources》2009
Groundwater characterization involves the resolution of unknown system characteristics from observation data, and is often classified as an inverse problem. Inverse problems are difficult to solve due to natural ill-posedness and computational intractability. Here we adopt the use of a simulation–optimization approach that couples a numerical pollutant-transport simulation model with evolutionary search algorithms for solution of the inverse problem. In this approach, the numerical transport model is solved iteratively during the evolutionary search. This process can be computationally intensive since several hundreds to thousands of forward model evaluations are typically required for solution. Given the potential computational intractability of such a simulation–optimization approach, parallel computation is employed to ease and enable the solution of such problems. In this paper, several variations of a groundwater source identification problem is examined in terms of solution quality and computational performance. The computational experiments were performed on the TeraGrid cluster available at the National Center for Supercomputing Applications. The results demonstrate the performance of the parallel simulation–optimization approach in terms of solution quality and computational performance. 相似文献
12.
The transport of contaminants in aquifers is usually represented by a convection-dispersion equation. There are several well-known problems of oscillation and artificial dispersion that affect the numerical solution of this equation. For example, several studies have shown that standard treatment of the cross-dispersion terms always leads to a negative concentration. It is also well known that the numerical solution of the convective term is affected by spurious oscillations or substantial numerical dispersion. These difficulties are especially significant for solute transport in nonuniform flow in heterogeneous aquifers. For the case of coupled reactive-transport models, even small negative concentration values can become amplified through nonlinear reaction source/sink terms and thus result in physically erroneous and unstable results. This paper includes a brief discussion about how nonpositive concentrations arise from numerical solution of the convection and cross-dispersion terms. We demonstrate the effectiveness of directional splitting with one-dimensional flux limiters for the convection term. Also, a new numerical scheme for the dispersion term that preserves positivity is presented. The results of the proposed convection scheme and the solution given by the new method to compute dispersion are compared with standard numerical methods as used in MT3DMS. 相似文献
13.
A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field. 相似文献
14.
Several technologies for cleaning up DNAPLs in source zones rely on solubilizing contaminants or destroying them in situ. Typically, these approaches employ an injection/withdrawal system to recirculate this treatment fluids. Our interest is in examining factors that influence delivery efficiency. Although many factors can affect this efficiency, this study looks at the combined influence of density-driven flow and porous media heterogeneities. The analysis is based on a series of numerical simulations of hypothetical chemical floods (e.g., potassium permanganate), which are highly resolved in space and time with a scale that is typical of field installations. Results indicate that the characteristics of convective mixing, i.e., natural, forced, and mixed convection, greatly affect delivery efficiency and patterns of mass transport. The ratio of the Grashof number (Gr) and Reynolds number (Re) proved useful in interpreting the patterns of flooding in a homogeneous porous medium. When higher-permeability layers are included in the domain, they act as conduits, effectively expediting the transport of treatment chemicals from injection to recovery wells. These large fluxes of chemicals in the high-permeability layers produce significant flooding inefficiencies. The problem is less severe in heterogeneous medium where the connectivity through the treatment zone is less well developed. Overall, this paper illustrates that density effects and high-permeability pathways need to be considered in designing chemical floods. 相似文献
15.
基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段. 相似文献
16.
The current study focuses on the application of a three-dimensional numerical model for the prediction of morphological bed changes. The sediment deposition in a reservoir during a 10-year-flood was investigated and the results of the simulation were validated with data derived from a physical model study. Because of the small grain sizes in the prototype, synthetic granulate was used in the physical model. The numerical computation domain was a reproduction of the physical model, including the grain sizes and the density of the particles, in order to ensure comparability. The CFD code SSIIM, which solves the RANS-equations in three-dimensions, was used for the simulations. The sediment transport in SSIIM is divided into suspended sediment transport, computed by solving the convection-diffusion equation, and bed-load transport, calculated by an empirical formula. The results of the numerical simulation correspond well to the results of the physical model study. The simulated location and the pattern of the sediment deposition in the reservoir are an accurate representation of the observed distribution in the physical model. 相似文献
17.
Hu BX 《Ground water》2006,44(2):222-233
A Lagrangian stochastic approach is applied to develop a method of moment for solute transport in a physically and chemically nonstationary medium. Stochastic governing equations for mean solute flux and solute covariance are analytically obtained in the first-order accuracy of log conductivity and/or chemical sorption variances and solved numerically using the finite-difference method. The developed method, the numerical method of moments (NMM), is used to predict radionuclide solute transport processes in the saturated zone below the Yucca Mountain project area. The mean, variance, and upper bound of the radionuclide mass flux through a control plane 5 km downstream of the footprint of the repository are calculated. According to their chemical sorption capacities, the various radionuclear chemicals are grouped as nonreactive, weakly sorbing, and strongly sorbing chemicals. The NMM method is used to study their transport processes and influence factors. To verify the method of moments, a Monte Carlo simulation is conducted for nonreactive chemical transport. Results indicate the results from the two methods are consistent, but the NMM method is computationally more efficient than the Monte Carlo method. This study adds to the ongoing debate in the literature on the effect of heterogeneity on solute transport prediction, especially on prediction uncertainty, by showing that the standard derivation of solute flux is larger than the mean solute flux even when the hydraulic conductivity within each geological layer is mild. This study provides a method that may become an efficient calculation tool for many environmental projects. 相似文献
18.
Boundary conditions are required to close the mathematical formulation of unstable density‐dependent flow systems. Proper implementation of boundary conditions, for both flow and transport equations, in numerical simulation are critical. In this paper, numerical simulations using the FEFLOW model are employed to study the influence of the different boundary conditions for unstable density‐dependent flow systems. A similar set up to the Elder problem is studied. It is well known that the numerical simulation results of the standard Elder problem are strongly dependent on spatial discretization. This work shows that for the cases where a solute mass flux boundary condition is employed instead of a specified concentration boundary condition at the solute source, the numerical simulation results do not vary between different convective solution modes (i.e., plume configurations) due to the spatial discretization. Also, the influence of various boundary condition types for nonsource boundaries was studied. It is shown that in addition to other factors such as spatial and temporal discretization, the forms of the solute transport equation such as divergent and convective forms as well as the type of boundary condition employed in the nonsource boundary conditions influence the convective solution mode in coarser meshes. On basis of the numerical experiments performed here, higher sensitivities regarding the numerical solution stability are observed for the Adams‐Bashford/Backward Trapezoidal time integration approach in comparison to the Euler‐Backward/Euler‐Forward time marching approach. The results of this study emphasize the significant consequences of boundary condition choice in the numerical modeling of unstable density‐dependent flow. 相似文献
19.
Impact of Local Groundwater Flow Model Errors on Transport and a Practical Solution for the Issue 
下载免费PDF全文
下载免费PDF全文 A groundwater flow model is typically used to provide the flow field for conducting groundwater solute transport simulations. The advection term of the mass conserved formulation for groundwater transport assumes that the flow field is perfectly balanced and that all water flowing into a numerical grid cell is exactly balanced by outflows after accounting for sources/sinks or internal storage. However, in many complicated regional or site‐scale models, there may be localized flow balance errors that may be difficult to eliminate through tighter flow convergence tolerances due to simulation time constraints or numerical limits on convergence tolerances. Thus, if water is erroneously gained or lost within a grid cell during the flow computation, the solutes within it will also be numerically affected in the associated transport simulation. Transport solutions neglect this error in groundwater flow as the transport equations that are solved assume no error in flow. This flow imbalance error can however have consequences on the transport solution ranging from unnoticeable errors in the resulting concentrations to spurious oscillations that can grow in time and hinder further solution. An approach has been suggested here, to explicitly handle these flow imbalances during mass conserved advective transport computations and report them in the corresponding transport mass balance output, as corrections that are needed to handle errors originating in the flow solution. Example problems are provided to explain the concepts and demonstrate the impacts. 相似文献
20.
《Advances in water resources》2005,28(1):33-42
A numerical solution that is significantly more general than other semi-analytical solutions is presented for governing equations describing advective–dispersive transport with multirate mass transfer between mobile and immobile domains. The new solution approach is general in the sense that it does not impose any restrictive assumption on the spatial or temporal variability of advective and dispersive processes in the mobile domain. A single integro-differential equation (IDE) is developed for the concentration in the mobile domain by separating the concentration in the immobile domain from the set of two partial differential equations. The solution to the IDE requires the evaluation of a temporal integral of the concentration in the mobile domain, which is a function of the Laplace transform of the distribution of the mass transfer rate coefficient. The Laplace transform is not limited to flow fields with known constant velocities. The solutions for one- and two-dimensional examples obtained using the new approach agree with those obtained by existing semi-analytical and numerical approaches. 相似文献