共查询到20条相似文献，搜索用时 15 毫秒
1.
Analytical and numerical solutions for wave propagation in water-saturated porous layered half-space
J. Tabatabaie Yazdi S. Valliappan Chongbin Zhao 《Soil Dynamics and Earthquake Engineering》1994,13(4)
A fluid-saturated one-layer continuum underlain by a rigid half-space is considered. An exact solution is developed in frequency domain for analyzing disturbance induced by a strip footing located at the surface with vertical harmonic excitation. Since the analytical solution can be used only for very simple conditions, a finite element model has been developed also and compared with the exact solution. It is shown that finite element results are in close agreement with the results which have been obtained by a transformation technique. The proposed finite element scheme can take into account the complex geometry and inhomogeneity for practical problems. Besides this, the analytical results exhibit the overall characteristic of wave propagation in porous media and will provide a representative test problem which can be used for a quantitative evaluation of the accuracy of various numerical solution methods. 相似文献
2.
The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems. 相似文献
3.
《Advances in water resources》2002,25(3):279-291
An analytical approach is presented for solving problems of steady, two-dimensional groundwater flow with inhomogeneity boundaries. A common approach for such problems is to separate the problem domain into two homogeneous domains, search for solutions in each domain, and then attempt to match conditions, either exactly or approximately, along the inhomogeneity boundary. Here, we use classical solutions to problems with inhomogeneity boundaries with simple geometries, and map conformally the entire domain onto a new one. In this way, existing solutions are used to solve problems with more complex, and more practical, boundary geometries. The approach is general, but subject to some restrictions on the mapping functions that may be used.Using this approach, we develop explicit analytical solutions for two problems of practical interest. The first problem addresses aquifer interaction across a gap in an impermeable separating layer; flow regimes are defined and the interaction is quantified. The second solution represents flow in the vertical plane to a partially clogged stream bed that is partially penetrating the aquifer; the stream bed is modeled as a thin layer of low-permeability silt. Flow regimes for groundwater surface–water interaction are quantified analytically. 相似文献
4.
The problem of free vibration of non-linear structures is considered initially. It is shown that this problem can be represented as a non-linear eigenvalue problem. Variational principles for non-linear eigenvalue problems are defined. These variational principles are implemented with finite element models to define numerical approximations for the free vibration problem. The solution of these approximate equations provides a set of non-linear modal vectors and natural frequencies which vary with the amplitude of the solution. The non-linear eigenvalue parameters can be used in modal expansion approximations for the non-linear transient or steady state response of structural systems. To demonstrate the proposed techniques the free vibration and steady state vibration characteristics of a geometrically non-linear circular plate are determined. 相似文献
5.
A three-dimensional finite difference model was developed to study the groundwater flow system in an upland area bordering a lake. For a general perspective of the groundwater flow system, a steady state three-dimensional flow was employed. Having determined the flow net by using a three-dimensional model, the flow volumes under natural conditions have been used to establish the parameter values and for the analyses of flow patterns. Further, to study the effects of human impact and precipitation on groundwater flow conditions in a small area, a transient three-dimensional simulation was performed. Environmental tritium was used to trace the regional groundwater movement to verify the three-dimensional mathematical model. Results obtained using the three-dimensional mathematical model approach and tritium concentration analyses were in close agreement. The results demonstrated that the groundwater flow system should be analysed using a three-dimensional geometric concept of groundwater movement. 相似文献
6.
A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow. 相似文献
7.
Free surface and saturated-unsaturated analyses of borehole infiltration tests above the water table
Constant head borehole infiltration tests are widely used for the in situ evaluation of saturated hydraulic conductivities of unsaturated soils above the water table. The formulae employed in analysing the results of such tests disregard the fact that some of the infiltrating water may flow under unsaturated conditions. Instead, these formulae are based on various approximations of the classical free surface theory which treats the flow region as if it were fully saturated and enclosed within a distinct envelope, the so-called ‘free surface’. A finite element model capable of solving free surface problems is used to examine the mathematical accuracy of the borehole infiltration formulae. The results show that in the hypothetical case where unsaturated flow does not exist, the approximate formulae are reasonably accurate within·a practical range of borehole conditions. To see what happens under conditions closer to those actually encountered in the field, the effect of unsaturated flow on borehole infiltration is investigated by means of two different numerical models: a mixed explicit-implicit finite element model, and a mixed explicit-implicit integrated finite difference model. Both of these models give nearly identical results; however, the integrated finite difference model is considerably faster than the finite element model. The relatively low computational efficiency of the finite element scheme is attributed to the large number of operations required in order to re-evaluate the conductivity (stiffness) matrix at each iteration in this highly non-linear saturated-unsaturated flow problem. The saturated-unsaturated analysis demonstrates that the classical free surface approach provides a distorted picture of the flow pattern in the soil. Contrary to what one would expect on the basis of this theory, only a finite region of the soil in the immediate vicinity of the borehole is saturated, whereas a significant percentage of the flow takes place under unsaturated conditions. As a consequence of disregarding unsaturated flow, the available formulae may underestimate the saturated hydraulic conductivity of fine grained soils by a factor of two, three, or more. Our saturated-unsaturated analysis leads to an improved design of borehole infiltration tests and a more accurate method for interpreting the results of such tests. The analysis also shows how one can predict the steady state rate of infiltration from data collected during the early transient period of the test. 相似文献
8.
Katerina Spanoudaki Anastasios I. Stamou Aikaterini Nanou-Giannarou 《Journal of Hydrology》2009,375(3-4):410-427
Coupled modelling of surface and subsurface systems is a valuable tool for quantifying surface water–groundwater interactions. In the present paper, the 3-D non-steady state Navier–Stokes equations, after Reynolds averaging and with the assumption of a hydrostatic pressure distribution, are for the first time coupled to the 3-D saturated groundwater flow equations in an Integrated suRface watEr–grouNdwater modEl (IRENE). A finite-difference method is used for the solution of the governing equations of IRENE. A semi-implicit scheme is used for the discretisation of the surface water flow equations and a fully implicit scheme for the discretisation of the groundwater flow equations. The two sets of equations are coupled at the common interface of the surface water and groundwater bodies, where water exchange takes place, using Darcy’s law. A new approach is proposed for the solution of the coupled surface water and groundwater equations in a simultaneous manner, in such a fashion that gives computational efficiency at low computational cost. IRENE is verified against three analytical solutions of surface water–groundwater interaction, which are chosen so that different components of the model can be tested. The model closely reproduces the results of the analytical solutions and can therefore be used for analysing and predicting surface water–groundwater interactions in real-world cases. 相似文献
9.
Jacob Zaidel 《Ground water》2013,51(6):952-959
Known analytical solutions of groundwater flow equations are routinely used for verification of computer codes. However, these analytical solutions (e.g., the Dupuit solution for the steady‐state unconfined unidirectional flow in a uniform aquifer with a flat bottom) represent smooth and continuous water table configurations, simulating which does not pose any significant problems for the numerical groundwater flow models, like MODFLOW. One of the most challenging numerical cases for MODFLOW arises from drying‐rewetting problems often associated with abrupt changes in the elevations of impervious base of a thin unconfined aquifer. Numerical solutions of groundwater flow equations cannot be rigorously verified for such cases due to the lack of corresponding exact analytical solutions. Analytical solutions of the steady‐state Boussinesq equation, associated with the discontinuous water table configurations over a stairway impervious base, are presented in this article. Conditions resulting in such configurations are analyzed and discussed. These solutions appear to be well suited for testing and verification of computer codes. Numerical solutions, obtained by the latest versions of MODFLOW (MODFLOW‐2005 and MODFLOW‐NWT), are compared with the presented discontinuous analytical solutions. It is shown that standard MODFLOW‐2005 code (as well as MODFLOW‐2000 and older versions) has significant convergence problems simulating such cases. The problems manifest themselves either in a total convergence failure or erroneous results. Alternatively, MODFLOW‐NWT, providing a good match to the presented discontinuous analytical solutions, appears to be a more reliable and appropriate code for simulating abrupt changes in water table elevations. 相似文献
10.
An analytic element approach is presented for the modeling of steady groundwater flow through multi-aquifer systems with piecewise constant aquifer and leaky layer properties. Different properties may be specified for domains bounded by closed polygons, referred to as polygonal inhomogeneities. The boundary of these inhomogeneities is modeled with two types of high-order line elements. First, a string of single-aquifer line-doublets is used; these elements cut through all aquifers and are valid both inside and outside the inhomogeneity. Second, two strings of multi-aquifer line-sinks are used, one string that is valid inside the inhomogeneity and one string that is valid outside; the comprehensive extraction of these line-sinks is zero at any point along the string. The proposed approach results in a comprehensive flow field of which the component normal to the boundary of the inhomogeneity is continuous across the boundary at any point. Within each individual aquifer, continuity of head and the component of flow normal to the boundary are met approximately across the boundary; the accuracy increases when the order of the line elements is increased and/or when shorter line elements are used. The proposed analytic element approach produces results that are virtually identical to the exact solution for a cylindrical inhomogeneity, and a high-resolution MODFLOW2000 model of two rectangular inhomogeneities with a shared boundary. The practical application of the approach is demonstrated through the solution of a problem with an irregularly shaped inhomogeneity with rivers crossing the inhomogeneity boundary. 相似文献
11.
本文提出了一种计算不规则起伏地形中SH波散射的有效方法——局域边界元法.本方法基于传统边界元法,为计算复杂地表散射问题提供了一种更加高效的解决方案.根据地震波满足的边界积分方程中牵引力格林函数的特性,我们将自由边界分解成水平部分和起伏部分.通过公式推导,可将水平部分的位移由起伏部分的位移通过格林函数线性叠加表示,因此只需对起伏部分的位移进行直接求解,从而极大地减少了待求解的未知数个数,显著提高了计算效率.通过与半圆形山谷SH波平面波入射的解析解比较,验证了方法的正确性.数值模型比较显示,局域边界元模拟结果与传统边界元数值解完全吻合,但是大幅提高了计算效率.因此,局域边界元法可以作为模拟不规则地形中地震波散射的有效工具. 相似文献
12.
考虑地球介质电导率任意各向异性且随空间位置连续变化的情况,本文实现了直接求解电磁场的可控源音频大地电磁测深（CSAMT）三维有限元数值模拟.首先给出了电导率任意各向异性介质中CSAMT二次电场满足的控制方程及其相应变分问题,然后采用任意六面体单元对研究区域进行剖分,在网格单元中对任意各向异性电导率进行线性插值,解决了实际工作中岩矿石电导率各向异性且连续变化的情况,将变分问题转化为线性代数方程组的求解.电导率各向异性且连续变化一维模型三维有限元数值模拟结果与电导率各向异性且分层均匀渐进模型解析解结果对比验证了方法的有效性;三维地电模型电导率随位置线性变化且各向同性、主轴各向异性、方位各向异性和倾斜各向异性的数值模拟结果表明,电导率各向异性且连续变化对CSAMT视电阻率和相位数据均有明显的影响. 相似文献
13.
Saltwater intrusion problems have been usually tackled through analytical models because of its simplicity, easy implementation and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, which neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. This paper provides insight into the validity of a sharp‐interface approximation defined from a steady state solution when applied to transient seawater intrusion problems. The validation tests have been performed on a 3D unconfined synthetic aquifer, which include spatial and temporal distribution of recharge and pumping wells. Using a change of variable, the governing equation of the steady state sharp‐interface problem can be written with the same structure of the steady confined groundwater flow equation as a function of a single potential variable (?). We propose to approach also the transient problem solving a single potential equation (using also the ? variable) with the same structure of the confined groundwater flow equation. It will allow solving the problem by using the classical MODFLOW code. We have used the parameter estimation model PEST to calibrate the parameters of the transient sharp‐interface equation. We show how after the calibration process, the sharp‐interface approach may provide accurate enough results when applied to transient problems and improve the steady state results, thus avoiding the need of implementing a density‐dependent model and reducing the computational cost. This has been proved by comparing results with those obtained using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. The comparison was performed in terms of piezometric heads, seawater penetration, transition zone width and critical pumping rates. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
In this work we propose upscaling method for nonlinear Forchheimer flow in heterogeneous porous media. The generalized Forchheimer law is considered for incompressible and slightly-compressible single-phase flows. We use recently developed analytical results (Aulisa et al., 2009) [1] and formulate the resulting system in terms of a degenerate nonlinear flow equation for the pressure with the nonlinearity depending on the pressure gradient. The coarse scale parameters for the steady state problem are determined so that the volumetric average of velocity of the flow in the domain on fine scale and on coarse scale are close. A flow-based coarsening approach is used, where the equivalent permeability tensor is first evaluated following streamline methods for linear cases, and modified in order to take into account the nonlinear effects. Compared to previous works (Garibotti and Peszynska, 2009) [2], (Durlofsky and Karimi-Fard) [3], this approach can be combined with rigorous mathematical upscaling theory for monotone operators, (Efendiev et al., 2004) [4], using our recent theoretical results (Aulisa et al., 2009) [1]. The developed upscaling algorithm for nonlinear steady state problems is effectively used for variety of heterogeneities in the domain of computation. Direct numerical computations for average velocity and productivity index justify the usage of the coarse scale parameters obtained for the special steady state case in the fully transient problem. For nonlinear case analytical upscaling formulas in stratified domain are obtained. Numerical results were compared to these analytical formulas and proved to be highly accurate. 相似文献
15.
A new upstream weighting finite element technique is developed for improved solution of the two-phase immiscible flow equations. Unlike the upstream weighting technique used by previous investigators, the new technique does not employ finite difference concepts to achieve the required upstream weighting of relative permeabilities or mobilities. Instead, upstream weighting is achieved by (1) representing the relative permeabilities or mobilities as continuous functions expressed in terms of the shape functions and nodal values (2) using asymmetric weighting functions to weight the spatial terms in the flow equations. These weighting functions are constructed such that they are dependent on the flow direction along each side of an element.In conjunction with the proposed technique, two solution schemes for treating the resulting set of non-linear algebraic equations are presented. These are the fully-implicit chord slope incremental solution scheme and the Newton-Raphson solution scheme. Both schemes allow the use of large time steps without being unstable.The proposed numerical technique is applied to two problems (1) the one-dimensional Buckley-Leverett problem (2) the two-dimensional five-spot well flow problem. Results indicate that this technique is superior to not only earlier finite element schemes but also five-point upstream finite difference formulae. 相似文献
16.
Thomas LUCKNER Ulrich ZANKE 《国际泥沙研究》2007,22(2):87-102
This paper presents an analytical solution for calculating the initiation of sediment motion and the risk of river bed movement. It thus deals with a fundamental problem in sediment transport, for which no complete analytical solution has yet been found. The analytical solution presented here is based on forces acting on a single grain in state of initiation of sediment motion. The previous procedures for calculating the initiation of sediment motion are complemented by an innovative combination of optical surface measurement technology for determining geometrical parameters and their statistical derivation as well as a novel approach for determining the turbulence effects of velocity fluctuations. This two aspects and the comparison of the solution functions presented here with the well known data and functions of different authors mainly differ the presented solution model for calculating the initiation of sediment motion from previous approaches. The defined values of required geometrical parameters are based on hydraulically laboratory tests with spheres. With this limitations the derivated solution functions permit the calculation of the effective critical transport parameters of a single grain, the calculation of averaged critical parameters for describing the state of initiation of sediment motion on the river bed, the calculation of the probability density of the effective critical velocity as well as the calculation of the risk of river bed movement. The main advantage of the presented model is the closed analytical solution from the equilibrium of forces on a single grain to the solution functions describing the initiation of sediment motion. 相似文献
17.
An axisymmetric infinite element and a three-dimensional infinite element are developed to solve three-dimensional elastic wave propagation problems in unbounded media. The elements are capable of transmitting Rayleigh, shear and compressional waves in the frequency domain. A scheme to integrate numerically the characteristic matrices of the elements is formulated based upon Gauss—Laguerre quadrature. Finally, the axisymmetric infinite element is used to find the compliance functions of a rigid circular plate subjected to harmonic loading on a semi-infinite medium. By using infinite elements, the size of the near field may be kept small. Consequently, the system is characterized by relatively few degrees of freedom, thus providing the analyst with an inexpensive solution. 相似文献
18.
Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems. 相似文献
19.
A numerical method has been proposed by Ross [Ross PJ. Modeling soil water and solute transport-fast, simplified numerical solutions. Agron J 2003; 95(6): 1352–1361.] to solve one-dimensional soil water movement problems. The Ross method is a noniterative numerical scheme, that can reduce computational time without sacrificing computational accuracy. The main aim of this study is to present a general form of the Ross method for two- and three-dimensional variably saturated flow. The established numerical model (R3D) is widely tested using five problems, in which the numerical solutions of R3D are compared with analytical solutions, laboratory data, and solutions from a traditional iterative numerical model. The comparison shows that R3D accommodates various hydraulic functions and boundary conditions. Results from R3D, which does not require iteration, are as accurate as results from iterative model. With the help of the primary variable switching technique, this model is unconditionally mass conservative, and computes infiltration into dry soil more efficiently. R3D is thus considered as an efficient tool for its high accuracy and efficiency for solving two- and three-dimensional variably saturated flow problems. 相似文献
20.
Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing. 相似文献