共查询到20条相似文献，搜索用时 312 毫秒
1.
Vectorized simulation of groundwater flow and contaminant transport using analytic element method and random walk particle tracking 
下载免费PDF全文
下载免费PDF全文 Groundwater contaminant transport processes are usually simulated by the finite difference (FDM) or finite element methods (FEM). However, they are susceptible to numerical dispersion for advection‐dominated transport. In this study, a numerical dispersion‐free coupled flow and transport model is developed by combining the analytic element method (AEM) with random walk particle tracking (RWPT). As AEM produces continuous velocity distribution over the entire aquifer domain, it is more suitable for RWPT than FDM/finite element methods. Using the AEM solutions, RWPT tracks all the particles in a vectorized manner, thereby improving the computational efficiency. The present model performs a convolution integral of the response of an impulse contaminant injection to generate concentration distributions due to a permanent contaminant source. The RWPT model is validated with an available analytical solution and compared to an FDM solution, the RWPT model more accurately replicates the analytical solution. Further, the coupled AEM‐RWPT model has been applied to simulate the flow and transport in hypothetical and field aquifer problems. The results are compared with the FDM solutions and found to be satisfactory. The results demonstrate the efficacy of the proposed method. 相似文献
2.
Under certain conditions the concentration of a substance moving in a stochastic flow field is described by the stochastic convection equation. A numerical method yielding the mean solution and variance of the two-dimensional problem is described here. First, the differential operator is replaced by a discrete linear operator based on finite differences. The resulting system of stochastic equations is then replaced by a system of equations whose solution is the mean concentration. The variance of the concentration can then be calculated. In addition, and example is given for which an approximate analytical solution and its variance is known. The numerical method is applied to the example and results compared to the approximate analytical solution and variance. 相似文献
3.
本文采用传播矩阵技术研究并建立了层状正交各向异性地层中多分量感应测井响应的有效算法.首先通过Fourier变换将频率空间域中的Maxwell方程组求解问题转化为频率波数域中关于电磁场水平分量常微分方程组的定解问题.利用该方程组系数矩阵的本征值和归一化本征向量将电磁场分解成上行波和下行波模式的组合,推导出均匀正交各向异性介质中由任意方向磁偶极子产生的电磁波模式解析表达式;在此基础上,利用叠加原理和边界条件研究了电磁波在层状正交各向异性地层中的反射和透射,给出各个界面上的广义反射系数和不同地层中电磁波振幅的递推公式,进而得到电磁波模式的解析解.为了有效确定频率空间域中的电磁场,采用二维Patterson自适应求积算法结合有限连分式展开技术计算傅氏逆变换.最后通过数值模拟结果证明了该算法的有效性,考察了不同各向异性系数、不同井眼倾角以及仪器长度和工作频率变化等情况下的多分量感应测井响应特征. 相似文献
4.
Dionyssia-Pinelopi N. Kontoni Dimitri E. Beskos 《Soil Dynamics and Earthquake Engineering》1987,6(4):227-238
In this paper the applicability of an approximate Boundary Element Method to uniform half-plane elastodynamic problems is investigated. This method employs the concept of images to construct approximate fundamental solutions for the half-plane and does not require any half-plane surface discretization. The method is formulated in the frequency domain for the case of harmonic disturbances or the Laplace transform domain for the case of transient disturbances. In the latter case a numerical inversion of the transformed solution is necessary to obtain the time domain response. The proposed method can be used as an alternative to boundary element methods that either utilize the infinite plane fundamental solution and thus require a half-plane surface discretization, or employ the exact half-plane fundamental solution, which even though leads to no surface discretization, is of a very lengthy and complicated form. Two characteristics numerical examples are used to illustrate the proposed method and study its advantages and disadvantages. 相似文献
5.
6.
近些年来提出的数值方法虽克服了应用复阻尼模型求解动力响应时可能出现的发散现象,但其计算过程繁杂且不能表达出结构动力特性及响应随结构参数的变化规律。论文参照三对角Toeplitz矩阵特征值问题的数学解法,推导出水平剪切型结构各阶自振频率、振型函数的解析形式。通过对运动方程进行Fourier变换,得到复阻尼理论下结构的传递函数解析式,直观地表达出结构动力特性及响应随结构参数的变化规律。最后通过算例对比分析,结果表明:提出的复阻尼模型的解析计算方法克服了时域发散问题,并与时域数值计算方法的位移响应时程对比发现,两种方法得到的时程曲线吻合程度较高,位移峰值也基本一致。同时,对比两种方法在底层刚度变化时的动力响应,解析方法的计算曲线更光滑,避免了数值方法的离散误差问题。 相似文献
7.
Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model. 相似文献
8.
Under certain conditions the concentration and flux of a substance moving in a stochastic flow field are described by the stochastic convection equation. A numerical method for solving the one-dimensional problem is studied here. The differential operator is replaced by a discrete linear operator based on finite differences. The resulting system of stochastic equations is then replaced by a system of equations whose solution is the mean concentration or mean flux. This final system is analysed and conditions for a stable numerical solution are obtained. Finally, numerical examples are given and are compared to an approximate analytical solution to the stochastic convection equation. 相似文献
9.
An approximate solution of the classical eigenvalue problem governing the vibrations of a relatively stiff structure on a soft elastic soil is derived through the application of a perturbation analysis. The full solution is obtained as the sum of the solution for an unconstrained elastic structure and small perturbing terms related to the ratio of the stiffness of the soil to that of the superstructure. The procedure leads to approximate analytical expressions for the system frequencies, modal damping ratios and participation factors for all system modes that generalize those presented earlier for the case of stiff soils. The resulting approximate expressions for the system modal properties are validated by comparison with the corresponding quantities obtained by numerical solution of the eigenvalue problem for a nine-story building. The accuracy of the proposed approach and of the classical normal mode approach is assessed through comparison with the exact frequency response of the test structure. 相似文献
10.
An approximate solution of the classical eigenvalue problem governing the vibrations of a structure on an elastic soil is derived through the application of a perturbation analysis. For stiff soils, the full solution is obtained as the sum of the solution for a rigid-soil and small perturbing terms related to the inverse of the soil shear modulus. The procedure leads to approximate analytical expressions for the system frequencies, modal damping ratios and participation factors for all system modes that generalize those presented by other authors for the fundamental mode. The resulting approximate expressions for the system modal properties are validated by comparison with the corresponding quantities obtained by numerical solution of the eigenvalue problem for a nine-story building. The accuracy of the proposed approach and of the classical normal mode approach is assessed through comparison with the exact frequency response of the test structure. 相似文献
11.
This study develops a lattice Boltzmann method (LBM) with a two-relaxation-time collision operator (LTRT) to solve saltwater intrusion problems. A directional-speed-of-sound (DSS) technique is introduced to take into account the hydraulic conductivity heterogeneity and discontinuity, as well as the velocity-dependent dispersion coefficient. The forcing terms in the LTRT model are customized in order to recover the density-dependent groundwater flow and mass transport equations. Using the LTRT with the squared DSS achieves at least second-order accuracy. The LTRT results are verified with Henry’s analytical solution as well as compared with several numerical examples and modified Henry problems that consider heterogeneous hydraulic conductivity and velocity-dependent dispersion. The numerical results show good agreement with the Henry analytical solution and with the numerical solutions obtained by other numerical methods. 相似文献
12.
采用波函数展开法给出了板的横截面处孔洞在平面P波入射下动应力集中问题的解。本文使用大圆弧假定法来满足板中的自由表面的边界条件。数值结果表明,上表面和孔洞间的距离与下表面和孔洞间的距离都对结果有影响。如果孔洞半径相对板的厚度很小,解近似等于孔洞全空间解。所给出的解是半解析解。因此能被用来验证诸如BEM,FEM,FDM等数值法。 相似文献
13.
This paper presents a new procedure to transform an SSI system into an equivalent SDOF system using twice equivalence. A pushover analysis procedure based on the capacity spectrum method for buildings with SSI effects （PASSI） is then established based on the equivalent SDOF system, and the modified response spectrum and equivalent capacity spectrum are obtained. Furthermore, the approximate formulas to obtain the dynamic stiffness of foundations are suggested. Three steel buildings with different story heights （3, 9 and 20） including SSI effects are analyzed under two far-field and two near-field historical records and an artificial seismic time history using the two PASSI procedures and the nonlinear response history analysis （NLhRHA） method. The results are compared and discussed. Finally, combined with seismic design response spectrum, the nonlinear seismic response of a 9-story building with SSI effects is analyzed using the PASSI procedures, and its seismic performance is evaluated according to the Chinese ＇Code for Seismic Design of Buildings. The feasibility of the proposed procedure is verified. 相似文献
14.
The first part of this paper presents an extensive validation of four analytical solutions for the seismic design of circular tunnels. The validation is performed with a quasi-static finite element (FE) model which conforms to the assumptions of the analytical solutions. Analyses are performed for a wide range of flexibility ratios, slippage conditions at soil–lining interface, assuming both drained and undrained behaviour. Based on the numerical predictions the relative merits of the considered analytical solutions are discussed and recommendations are given for their use in design. The second part of this paper explores the use of equivalent linear soil properties in analytical solutions as an approximate way of simulating nonlinearity. The results of equivalent linear site response analyses are used as an input for the analytical solutions. The comparison of the analytical predictions with nonlinear numerical analysis results is very satisfactory. The results of this study suggest that analytical solutions can be used for preliminary design using equivalent linear properties and the corresponding compatible strain as an approximate way of accounting for nonlinear soil response. 相似文献
15.
在地震动数值模拟方法中,谱元法和有限元法是应用较广泛的两种方法。基于经典的Lamb问题模型,首先推导给出地表竖向位移的解析解答。然后分别利用常用的四阶谱元法和线性有限元法,模拟了地表脉冲力源作用下模型的位移响应。考虑有意义的最短波长内的采样点个数及单元高宽比的变化,对比了两种方法的模拟精度;结果表明:对于谱元法,观测点与波源之间需至少包含两个网格,在此条件下,最短波长内包含一个网格（最短波长内5个采样点）时,数值解与解析解的误差小于1%,已达很高的精度;对于有限元法,最短波长内需包含10个网格时才能达到这一精度。此外,在满足网格尺寸要求的前提下,单元水平向与垂直向尺寸的比值在1∶1到5∶1的范围内时,谱元法和有限元法的模拟精度均变化不大。因此,单位波长内采样点个数相同时,谱元法的模拟精度比有限元法高的多,同时,在一定范围内两种方法的模拟结果对于宽高比的变化不敏感。 相似文献
16.
The advection–dispersion equation with spatially variable coefficients does not have an exact analytical solution and is therefore solved numerically. However, solutions obtained with several of the traditional finite difference or finite element techniques typically exhibit spurious oscillation or numerical dispersion when advection is dominant. The mixing cell and semi-analytical solution methods proposed in this study avoid such oscillation or numerical dispersion when advection dominates. Both the mixing cell and semi-analytical solution methods calculate the spatial step size by equating numerical dispersion to physical dispersion. Because of the spatial variability of the coefficients the spatial step size varies in space. When the time step size Δt→0, the mixing cell method reduces to the semi-analytical solution method. The results of application to two cases show that the mixing cell and semi-analytical solution methods are better than a finite difference method used in the study. © 1998 John Wiley & Sons, Ltd. 相似文献
17.
An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed. 相似文献
18.
Concurrent existence of confined and unconfined zones of an aquifer can arise owing to ground water withdrawal by pumping. Using Girinskii's potential function, Chen (1974, 1983) developed an approximate analytical solution to analyze transient ground water flow to a pumping well in an aquifer that changes from an initially confined system to a system with both unconfined and confined regimes. This article presents the details of the Chen model and then compares it with the analytical model developed by Moench and Prickett (1972) for the same problem. Hypothetical pumping test examples in which the aquifer undergoes conversion from confined to water table conditions are solved by the two analytical models and also a numerical model based on MODFLOW. Comparison of the results suggests that the solutions of the Chen model give better results than the Moench and Prickett model except when the radial distance is very large or aquifer thickness is large compared with drawdown. 相似文献
19.
本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用. 相似文献
20.
动荷载作用下欧拉梁动响应的计算是一个初边值问题,通常很难得到解析解,传统数值方法一般是把空间和时间分别离散进行求解,计算相对复杂,效率也不高.针对分布动荷载作用下欧拉梁的振动偏微分方程,采用传统微分求积法,在空间和时间上同时进行离散;对于所有非0阶的初/边值条件,采用嵌入法在权系数计算中予以考虑.算例的数值结果与精确解的对比证明采用传统微分求积法处理此问题是可行的,而且是高效的.对于实际工程中的其他类似问题,该方法同样适用. 相似文献