首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A distributed-parameter physically-based solute transport model using a novel approach to describe surface-subsurface interactions is coupled to an existing flow model. In the integrated model the same surface routing and mass transport equations are used for both hillslope and channel processes, but with different parametrizations for these two cases. For the subsurface an advanced time-splitting procedure is used to solve the advection-dispersion equation for transport and a standard finite element scheme is used to solve Richards equation for flow. The surface-subsurface interactions are resolved using a mass balance-based surface boundary condition switching algorithm that partitions water and solute into actual fluxes across the land surface and changes in water and mass storage. The time stepping strategy allows the different time scales that characterize surface and subsurface water and solute dynamics to be efficiently and accurately captured. The model features and performance are demonstrated in a series of numerical experiments of hillslope drainage and runoff generation.  相似文献   

2.
This work examines variable density flow and corresponding solute transport in groundwater systems. Fluid dynamics of salty solutions with significant density variations are of increasing interest in many problems of subsurface hydrology. The mathematical model comprises a set of non-linear, coupled, partial differential equations to be solved for pressure/hydraulic head and mass fraction/concentration of the solute component. The governing equations and underlying assumptions are developed and discussed. The equation of solute mass conservation is formulated in terms of mass fraction and mass concentration. Different levels of the approximation of density variations in the mass balance equations are used for convection problems (e.g. the Boussinesq approximation and its extension, fully density approximation). The impact of these simplifications is studied by use of numerical modelling.Numerical models for nonlinear problems, such as density-driven convection, must be carefully verified in a particular series of tests. Standard benchmarks for proving variable density flow models are the Henry, Elder, and salt dome (HYDROCOIN level 1 case 5) problems. We studied these benchmarks using two finite element simulators - ROCKFLOW, which was developed at the Institute of Fluid Mechanics and Computer Applications in Civil Engineering and FEFLOW, which was developed at the Institute for Water Resources Planning and Systems Research Ltd. Although both simulators are based on the Galerkin finite element method, they differ in many approximation details such as temporal discretization (Crank-Nicolson vs predictor-corrector schemes), spatial discretization (triangular and quadrilateral elements), finite element basis functions (linear, bilinear, biquadratic), iteration schemes (Newton, Picard) and solvers (direct, iterative). The numerical analysis illustrates discretization effects and defects arising from the different levels of the density of approximation. We contribute new results for the salt dome problem, for which inconsistent findings exist in literature. Applications of the verified numerical models to more complex problems, such as thermohaline and three-dimensional convection systems, will be presented in the second part of this paper.  相似文献   

3.
The alternating-direction collocation method has recently been developed for general parabolic equations. In order to test the applicability of the procedure to highly nonlinear problems, an alternating-direction collocation algorithm is developed to simulate two-dimensional flow in unsaturated porous media. The algorithm employs an alternating-direction solution procedure within the framework of a modified Picard iteration scheme. Numerical behaviour of the new procedure is compared to the behaviour of a standard two-dimensional collocation formulation. The new method is also tested on several infiltration problems of practical interest, including a layered and sloping soil. Results demonstrate the method to be accurate and highly mass conservative. The algorithm also produces significant savings in both execution time and storage.  相似文献   

4.
The “modified Picard” iteration method, which offers global mass conservation, can also be described as a form of Newton's iteration with lagged nonlinear coefficients. It converges to a time step with first-order discretization error. This paper applies second- and third-order diagonally implicit Runge Kutta (DIRK) time steps to the modified Picard method in one example. It demonstrates improvements over the first-order time step in rms error and error-times-effort model quality by factors ranging from two to over two orders of magnitude, showing that the “modified Picard” and DIRK methods are compatible.  相似文献   

5.
An efficient and accurate numerical model for multicomponent compressible single-phase flow in fractured media is presented. The discrete-fracture approach is used to model the fractures where the fracture entities are described explicitly in the computational domain. We use the concept of cross flow equilibrium in the fractures. This will allow large matrix elements in the neighborhood of the fractures and considerable speed up of the algorithm. We use an implicit finite volume (FV) scheme to solve the species mass balance equation in the fractures. This step avoids the use of Courant–Freidricks–Levy (CFL) condition and contributes to significant speed up of the code. The hybrid mixed finite element method (MFE) is used to solve for the velocity in both the matrix and the fractures coupled with the discontinuous Galerkin (DG) method to solve the species transport equations in the matrix. Four numerical examples are presented to demonstrate the robustness and efficiency of the proposed model. We show that the combination of the fracture cross-flow equilibrium and the implicit composition calculation in the fractures increase the computational speed 20–130 times in 2D. In 3D, one may expect even a higher computational efficiency.  相似文献   

6.
徐世刚  刘洋 《地球物理学报》2018,61(7):2950-2968
传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果.  相似文献   

7.
Abstract. During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models'documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport.  相似文献   

8.
This paper investigates the validity of a quasi-steady approximation for sediment transport and presents a new algorithm based on this concept. The developed non-coupled algorithm interacts among hydrodynamic, sediment, and morphology modules which are based on depth-averaged Navier-Stokes equations for the flow, the three-dimensional equation of conservation of sediment, and the mass balance between the bed and sediment (Exner equation) to simulate the reservoir sedimentation process. The non-coupled algorithm solves both the short-term scale and the relatively long-term scale problems of reservoir sedimentation. The proposed algorithm is verified using field data and by comparison with other accurate algorithms. Based upon the results of this investigation, the developed algorithm can be used to simulate long-term reservoir sedimentation while considerably decreasing the computational costs and preserving computational accuracy. The computational cost of the non-coupled algorithm is about 97% less than the conventional semi-coupled approach whereas the errors (Root Mean Square Error, Average Relative Error, and Maximum Relative Error of bed level) of the developed algorithm are approximately 15% greater than those for the semi-coupled algorithm for the average value.  相似文献   

9.
A Eulerian-Langrangian scheme is used to reformulate the equation of solute transport with ground water in saturated soils. The governing equation is decomposed into advection along characteristic path lines and propagation of the residue at a fixed grid.The method was employed to simulate transport of a conservative pollutant in a hypothetical aquifer, subject to the equivalence of real conditions. Implementation was based on data involving parameters of a heterogeneous aquifer, heavy flux stresses of densed pumpage/recharge wells, precipitation and seasonally changing flow regimes. Simulation, with coarse grid and high Peclet numbers yielded minute mass balance errors.  相似文献   

10.
The χMD matrix solver package is incorporated into USGS groundwater modeling software, such as MODFLOW-NWT, MODFLOW-USG, and MT3D. The solver is used to solve matrices assembled through numerical discretization of the groundwater flow equation, and solute transport equations. χMD has demonstrated its higher robustness, faster execution speed, and more efficient memory usage compared to the existing solvers for many types of groundwater flow problems. χMD uses preconditioned iterative Krylov-subspace methods and consists of preconditioning and acceleration modules. Because the solver package uses a variety of preconditioning features including level-based incomplete lower-upper (ILU) factorization method with a drop tolerance scheme, users must choose optimal preconditioning parameters to improve execution speed and robustness. In order to examine how the preconditioning parameters, ILU factorization level, and drop tolerance values affect the overall performance of the matrix solver, we evaluated five different groundwater model applications using MODFLOW-USG that include different numerical complexities. For those five cases, the number of discretization nodes varied from 10,000 cells to 730,300 cells. From the analysis, we found that the preconditioning parameters greatly affect execution times and memory usage of the preconditioning and acceleration procedures. In addition, a combination of the ILU level between five to seven and the drop tolerance value between 10−2 and 10−3 usually resulted in shorter overall execution time. Our study suggests that the users can elicit higher performance and robustness of the χMD matrix solver using this combination of the parameters and enhance computational efficiency of solving groundwater and solute transport problems.  相似文献   

11.
ABSTRACT

Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE.  相似文献   

12.
Abstract

The impact of pollution incidents on rivers and streams may be predicted using mathematical models of solute transport. Practical applications require an analytical or numerical solution to a governing solute mass balance equation together with appropriate values of relevant transport coefficients under the flow conditions of interest. This paper considers two such models, namely those proposed by Fischer and by Singh and Beck, and compares their performances using tracer data from a small stream in Edinburgh, UK. In calibrating the models, information on the magnitudes and the flow rate dependencies of the velocity and the dispersion coefficients was generated. The dispersion coefficient in the stream ranged between 0.1 and 0.9 m2/s for a flow rate range of 13–437 L/s. During calibration it was found that the Singh and Beck model fitted the tracer data a little better than the Fischer model in the majority of cases. In a validation exercise, however, both models gave similarly good predictions of solute transport at three different flow rates.  相似文献   

13.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant.  相似文献   

14.
The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.  相似文献   

15.
Although based on exact analytical solutions, semi‐analytical solute transport models can have significant numerical error in applications with high frequency oscillatory source terms and when parameter value combinations cause series solution approximations to converge slowly. Methods for correcting these numerical errors are presented and implemented in the AT123D code, which employs Green's functions to represent point, linear, and rectangular prismatic source zones. In order to increase its computational accuracy, a Romberg numerical integration scheme was added to AT123D with prespecified error criteria, variable time stepping, and partitioning of the integral to handle rapidly changing source terms. More rapidly converging series solution approximations for the Green's functions were also incorporated to improve both accuracy and computational efficiency for finite‐depth aquifers. AT123D also has been modified to eliminate redundant calculations at points where approximate steady‐state conditions have been reached to improve computational efficiency during numerical integration. These modifications help to decrease computer run times that can be excessive for three‐dimensional problems with large numbers of computational points, small time steps, and/or long simulation time periods. Errors in the original AT123D code also were corrected in this modified version, AT123D‐AT, in order to accurately simulate finite‐duration (pulse) source releases.  相似文献   

16.
In the paper a solution of two-dimensional (2D) nonlinear diffusive wave equation in a partially dry and wet domain is considered. The splitting technique which allows to reduce 2D problem into the sequence of one-dimensional (1D) problems is applied. The obtained 1D equations with regard to x and y are spatially discretized using the modified finite element method with the linear shape functions. The applied modification referring to the procedure of spatial integration leads to a more general algorithm involving a weighting parameter. Time integration is carried out using a two-level difference scheme with the weighting parameter as well. The resulting tri-diagonal systems of nonlinear algebraic equations are solved using the Picard iterative method. For particular sets of the weighting parameters, the proposed method takes the form of a standard finite element method and various schemes of the finite difference method. On the other hand, for the linear version of the governing equation, the proper values of the weighting parameters ensure an approximation of 3rd order. Since the diffusive wave equation can be solved no matter whether the area is dry or wet, the numerical computations can be carried out over entire domain of solution without distinguishing a current position of the shoreline which is obtained as a result of solution.  相似文献   

17.
A numerical method based on the modified method of characteristics is developed for incompressible Darcy flow. Fluid elements modeled as grid cells are mapped back in time to their twisted forms and a strict equality of volumes is imposed between the two. These relations are then cast in terms of potentials using Darcy's law and a nonlinear algebraic problem is solved for potentials. Though a general technique for obtaining Darcy flow, this method is most useful when the solute advection problem also is solved with the modified method of characteristics. The combined technique (referred to as the characteristic-conservative method) using the same characteristics to obtain both velocities and concentrations is then a direct numerical approximation to the Reynolds transport theorem. The method is implemented in three dimensions and a few sample problems featuring nonuniform flow-fields are solved to demonstrate the exact mass conservation property. Inflow and outflow boundaries do not cause any problems in the implementation. In all cases, the characteristic-conservative method obtains velocities that preserve fluid volume and, concentrations that achieve exact local and global mass balance; a desirable property that usually eludes characteristics based methods for solute advection in multidimensional, nonuniform flow-fields.  相似文献   

18.
A p finite element scheme and parallel iterative solver are introduced for a modified form of the shallow water equations. The governing equations are the three-dimensional shallow water equations. After a harmonic decomposition in time and rearrangement, the resulting equations are a complex Helmholz problem for surface elevation, and a complex momentum equation for the horizontal velocity. Both equations are nonlinear and the resulting system is solved using the Picard iteration combined with a preconditioned biconjugate gradient (PBCG) method for the linearized subproblems. A subdomain-based parallel preconditioner is developed which uses incomplete LU factorization with thresholding (ILUT) methods within subdomains, overlapping ILUT factorizations for subdomain boundaries and under-relaxed iteration for the resulting block system. The method builds on techniques successfully applied to linear elements by introducing ordering and condensation techniques to handle uniform p refinement. The combined methods show good performance for a range of p (element order), h (element size), and N (number of processors). Performance and scalability results are presented for a field scale problem where up to 512 processors are used.  相似文献   

19.
The inherent heterogeneity of geological media often results in anomalous dispersion for solute transport through them, and how to model it has been an interest over the past few decades. One promising approach that has been increasingly used to simulate the anomalous transport in surface and subsurface water is the fractional advection–dispersion equation (FADE), derived as a special case of the more general continuous time random walk or the stochastic continuum model. In FADE, the dispersion is not local and the solutes have appreciable probability to move long distances, and thus reach the boundary faster than predicted by the classical advection–dispersion equation (ADE). How to deal with different boundaries associated with FADE and their consequent impact is an issue that has not been thoroughly explored. In this paper we address this by taking one-dimensional solute movement in soil columns as an example. We show that the commonly used FADE with its fractional derivatives defined by the Riemann–Liouville definition is problematic and could result in unphysical results for solute transport in bounded domains; a modified method with the fractional dispersive flux defined by the Caputo derivatives is presented to overcome this problem. A finite volume approach is given to numerically solve the modified FADE and its associated boundaries. With the numerical model, we analyse the inlet-boundary treatment in displacement experiments in soil columns, and find that, as in ADE, treating the inlet as a prescribed concentration boundary gives rise to mass-balance errors and such errors could be more significant in FADE because of its non-local dispersion. We also discuss a less-documented but important issue in hydrology: how to treat the upstream boundary in analysing the lateral movement of tracer in an aquifer when the tracer is injected as a pulse. It is shown that the use of an infinite domain, as commonly assumed in literature, leads to unphysical backward dispersion, which has a significant impact on data interpretation. To avoid this, the upstream boundary should be flux-prescribed and located at the upstream edge of the injecting point. We apply the model to simulate the movement of Cl− in a tracer experiment conducted in a saturated hillslope, and analyse in details the significance of upstream-boundary treatments in parameter estimation.  相似文献   

20.
A quasi three-dimensional (QUASI 3-D) model is presented for simulating the subsurface water flow and solute transport in the unsaturated and in the saturated zones of soil. The model is based on the assumptions of vertical flow in the unsaturated zone and essentially horizontal groundwater flow. The 1-D Richards equation for the unsaturated zone is coupled at the phreatic surface with the 2-D flow equation for the saturated zone. The latter was obtained by averaging 3-D flow equation in the saturated zone over the aquifer thickness. Unlike the Boussinesq equation for a leaky-phreatic aquifer, the developed model does not contain a storage term with specific yield and a source term for natural replenishment. Instead it includes a water flux term at the phreatic surface through which the Richards equation is linked with the groundwater flow equation. The vertical water flux in the saturated zone is evaluated on the basis of the fluid mass balance equation while the horizontal fluxes, in that equation, are prescribed by Darcy law. A 3-D transport equation is used to simulate the solute migration. A numerical algorithm to solve the problem for the general quasi 3-D case was developed. The developed methodology was exemplified for the quasi 2-D cross-sectional case (QUASI2D). Simulations for three synthetic problems demonstrate good agreement between the results obtained by QUASI2D and two fully 2-D flow and transport codes (SUTRA and 2DSOIL). Yet, simulations with the QUASI2D code were several times faster than those by the SUTRA and the 2DSOIL codes.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号