共查询到20条相似文献，搜索用时 15 毫秒
1.
Damping solvent extraction is a finite element method for the analysis of unbounded (visco-)elastic media which was suggested by Wolf and Song in 1994. It was originally recommended that the method should be employed with a variable domain size depending on excitation frequency. Furthermore, other researchers who have utilized this method in the context of constant domain size, have often imposed strict conditions on the mesh size for the whole domain, which reduces the effectiveness of the approach. Considering the effect of artificial damping on mesh density selection, the present study introduces damping solvent as a method in which one can relax typical mesh density requirements to a large extent by utilizing a large value of artificial damping. Therefore, it makes the finite element mesh to benefit from a large domain size which improves the results for low frequency range. Moreover, good results are simultaneously obtained for high frequency range due to employing a large value of artificial damping. To illustrate the point, a rigid strip foundation with a cross section of rectangle embedded in half plane is considered. According to some comparison between the results obtained from several finite element meshes, the best one which takes full advantages of a large value of artificial damping for dynamic stiffness coefficients of strip foundation is introduced. These comparisons are carried out on domain size, mesh density and artificial damping. 相似文献
2.
有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响. 相似文献
3.
A finite element procedure based on the extension of the saturated flow domain into the partially saturated zone above the free surface is proposed. The finite element equations are derived by using a pseudo variational principle which results into a residual or correction load vector. The steady or transient free surface is corrected by using the residual load in an iterative scheme. The proposed procedure uses only one (initial) mesh and does not require modification of the mesh during iterations. It is compared (qualitatively) with other procedures such as variable mesh and variational inequalities. The procedure provides satisfactory comparisons with a number of closed-form solutions and laboratory test results; two applications involving the latter are described in the paper. 相似文献
4.
The nodal domain integration method is applied to a one-dimensional advection—diffusion mathematical model without a source term. Comparison of the resulting numerical model to the well known Galerkin finite element, subdomain, and finite difference domain models indicates that a single numerical statement can be developed which includes the Galerkin finite element, subdomain, and finite difference models as special cases. 相似文献
5.
T.V. Hromadka 《Advances in water resources》1984,7(2):76-80
The nodal domain integration method is applied to a two-dimensional advection-diffusion process in an anisotropic inhomogeneous medium. The domain is discretised into the union of irregular triangle finite elements with vertex-located nodal points and a linear trial function is used to approximate the governing flow equation's state variable in each element. Non-linear parameters are assumed quasi-constant for small durations in time in each element. The resulting numerical model represents the Galerkin and subdomain integration weighted residual methods and the integrated finite difference method as special cases. Both Dirichlet and Neumann boundary conditions are accommodated in a manner similar to the Galerkin finite element approach. 相似文献
6.
This paper presents development of a special finite difference method for the nonlinear dynamic response analysis of semi-infinite foundation soil. Semi-infinite domain is mapped into the finite domain using special mapping. For the region of engineering interest, mapping is isometrical, and for far field, shrink mapping which transforms an infinite interval into a finite interval is adopted. Using linear and nonlinear constitutive models, the responses of semi-infinite foundation soil are computed using a proposed method with a small mesh model and an extensive mesh model. Surface loadings or incident earthquake waves are applied to the models in the computations. Good agreements were obtained among the theoretical and computed results of the two models and the effectiveness of the proposed method was demonstrated. 相似文献
7.
在复杂山地和复杂海底条件下,地表和海底的剧烈起伏对地震波数值模拟提出了更高的要求.常规有限差分法采用矩形网格对模型进行网格剖分,由于矩形网格自身的限制,起伏地表或起伏海底只能由一系列阶梯状折线代替,从而引起人为虚假绕射波.此外,在模拟液-固界面的反射波时,如果界面与网格线不一致,则需要更密的网格才能得到精确的结果.为了解决上述问题,本文将自适应网格生成技术引入到起伏海底速度模型的网格剖分中,采用高阶仿真型有限差分法(MFD)对曲线坐标下的声波方程波进行了数值模拟.利用自适应网格生成技术对速度模型进行网格剖分不仅可以准确地描述模型边界,而且可以有效消除虚假绕射波.高阶仿真型有限差分法可以有效压制频散提高计算精度.模型试算结果表明,本文方法对复杂海底模型具有很好的适应性. 相似文献
8.
We propose a domain decomposition method based on the spectral element method(DDM-SEM)for elastic wave computation in frequency domain.It combines the high accuracy of the spectral element method and the high degree of parallelism of a domain decomposition technique,which makes this method suitable for accurate and efficient simulations of large scale problems in elastodynamics.In the DDM-SEM,the original large-scale problem is divided into a number of well designed subdomains.We use the spectral element method independently for each subdomain,and the neighboring subdomains are connected by a frequency-domain version of Riemann transmission condition(RTC)for elastic waves.For the proposed method,we can employ the non-conforming meshes and different interpolation orders in different subdomains to maximize the efficiency.By separating the internal and boundary unknowns of each subdomain,an efficient and naturally parallelizable block LDU direct solver is developed to solve the final system matrix.Numerical experiments verify its accuracy and efficiency,and show that the proposed DDM-SEM can be a promising numerical tool for accurately and effectively solving large and multi-scale problems of elastic waves.It is potentially valuable for the frequency domain seismic inversion where multiple source illuminations are required. 相似文献
9.
《Advances in water resources》1999,23(3):229-237
It can be very time consuming to use the conventional numerical methods, such as the finite element method, to solve convection–dispersion equations, especially for solutions of large-scale, long-term solute transport in porous media. In addition, the conventional methods are subject to artificial diffusion and oscillation when used to solve convection-dominant solute transport problems. In this paper, a hybrid method of Laplace transform and finite element method is developed to solve one- and two-dimensional convection–dispersion equations. The method is semi-analytical in time through Laplace transform. Then the transformed partial differential equations are solved numerically in the Laplace domain using the finite element method. Finally the nodal concentration values are obtained through a numerical inversion of the finite element solution, using a highly accurate inversion algorithm. The proposed method eliminates time steps in the computation and allows using relatively large grid sizes, which increases computation efficiency dramatically. Numerical results of several examples show that the hybrid method is of high efficiency and accuracy, and capable of eliminating numerical diffusion and oscillation effectively. 相似文献
10.
根据库伦规范势的定义,推导出关于磁矢量势和电标量势的偏微分方程,为了克服由电流源引起的奇异性和数值模拟计算困难,将电磁总场分解为一次场和二次场,一次场由基于Schelkunoff势函数的一维正演算法得到,二次场由有限元法计算得到,实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例,验证该算法的计算精度。然后,利用该算法对带海底地形的三层储层模型进行正演,分析了海底地形对海洋控源电磁场各分量产生的影响。 相似文献
11.
Qian-Jiang Zhang Shi-Kun Dai Long-Wei Chen Jian-Ke Qiang Kun Li Dong-Dong Zhao 《应用地球物理》2016,13(2):257-266
To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement. 相似文献
12.
Various numerical methods have been used in the literature to simulate single and multiphase flow in fractured media. A promising approach is the use of the discrete-fracture model where the fracture entities in the permeable media are described explicitly in the computational grid. In this work, we present a critical review of the main conventional methods for multiphase flow in fractured media including the finite difference (FD), finite volume (FV), and finite element (FE) methods, that are coupled with the discrete-fracture model. All the conventional methods have inherent limitations in accuracy and applications. The FD method, for example, is restricted to horizontal and vertical fractures. The accuracy of the vertex-centered FV method depends on the size of the matrix gridcells next to the fractures; for an acceptable accuracy the matrix gridcells next to the fractures should be small. The FE method cannot describe properly the saturation discontinuity at the matrix–fracture interface. In this work, we introduce a new approach that is free from the limitations of the conventional methods. Our proposed approach is applicable in 2D and 3D unstructured griddings with low mesh orientation effect; it captures the saturation discontinuity from the contrast in capillary pressure between the rock matrix and fractures. The matrix–fracture and fracture–fracture fluxes are calculated based on powerful features of the mixed finite element (MFE) method which provides, in addition to the gridcell pressures, the pressures at the gridcell interfaces and can readily model the pressure discontinuities at impermeable faults in a simple way. To reduce the numerical dispersion, we use the discontinuous Galerkin (DG) method to approximate the saturation equation. We take advantage of a hybrid time scheme to alleviate the restrictions on the size of the time step in the fracture network. Several numerical examples in 2D and 3D demonstrate the robustness of the proposed model. Results show the significance of capillary pressure and orders of magnitude increase in computational speed compared to previous works. 相似文献
13.
Meshing effects of the 3-D FEM numerical modeling in seismo-electromagnetics: An application in selectivity of seismic electric signal (SES) 
下载免费PDF全文
下载免费PDF全文 We investigated how density and quality of mesh around interest domain affect electromagnetic (EM) responses of 3D Earth layered media using finite element method (FEM). Effect of different mesh shapes was also investigated using a method of mixing structured and unstructured mesh. As a case study, we estimated the effects of meshing on selectivity phenomenon of seismic electric signal (SES). Our results suggest that the relative errors resulting from mesh effects may not be negligible, which may lead to some unconvincing explanation of the SES selectivity based on the numerical modeling results. 相似文献
14.
The nodal domain integration method is used to develop a numerical model of the linear diffusion equation. The nodal domain integration approach is shown to represent an infinity of finite element mass matrix lumping schemes including the Galerkin and subdomain integration versions of the weighted residual method and an integrated finite difference method. Neumann, Dirichlet and mixed boundary conditions are accommodated analogous to the Galerkin finite element method. In order to reduce the overall integrated approximation relative error, a mass matrix lumping formulation is developed which is based on the Crank-Nicolson time advancement approximation. The optimum mass lumping factors are found to be strongly related to the model timestep size. 相似文献
15.
A three-dimensional baroclinic finite element model with a coarse and fine (i.e. local refinement along the shelf edge) grid
is used to examine the influence of shelf edge grid refinement upon the internal tide generation and propagation off the west
coast of Scotland. Comparisons are made with observations in the region and with a published solution using a finite difference
model. The calculations show that provided that the finite element grid is refined in the internal tide generation area and
the adjacent region through which the internal tide propagates, then a numerically accurate solution is obtained. In the regions
of strong internal tide generation with a local grid refinement, internal wave energy can accumulate at small scales and must
be removed by a scale-selective filter. 相似文献
16.
The coupling upscaling finite element method is developed for solving the coupling problems of deformation and consolidation of heterogeneous saturated porous media under external loading conditions. The method couples two kinds of fully developed methodologies together, i.e., the numerical techniques developed for calculating the apparent and effective physical properties of the heterogeneous media and the upscaling techniques developed for simulating the fluid flow and mass transport properties in heterogeneous porous media. Equivalent permeability tensors and equivalent elastic modulus tensors are calculated for every coarse grid block in the coarse-scale model of the heterogeneous saturated porous media. Moreover, an oversampling technique is introduced to improve the calculation accuracy of the equivalent elastic modulus tensors. A numerical integration process is performed over the fine mesh within every coarse grid element to capture the small scale information induced by non-uniform scalar field properties such as density, compressibility, etc. Numerical experiments are carried out to examine the accuracy of the developed method. It shows that the numerical results obtained by the coupling upscaling finite element method on the coarse-scale models fit fairly well with the reference solutions obtained by traditional finite element method on the fine-scale models. Moreover, this method gets more accurate coarse-scale results than the previously developed coupling multiscale finite element method for solving this kind of coupling problems though it cannot recover the fine-scale solutions. At the same time, the method developed reduces dramatically the computing effort in both CPU time and memory for solving the transient problems, and therefore more large and computational-demanding coupling problems can be solved by computers. 相似文献
17.
本文采用平面有限元数值模拟方法,对阿尔金断裂带及邻区的构造应力场和剪应力集中的能力与阿尔金断裂破裂性质之间的关系进行了探讨。通过与实际资料对比,获得了本区应力场计算模型。通过计算得出本断裂带各部段切应力值不均一,整个断裂带活动不频繁,未来地震的危险区位于西部段。 相似文献
18.
The finite‐volume technique is used to solve the two‐dimensional shallow‐water equations on unstructured mesh consisting of quadrilateral elements. In this paper the algorithm of the finite‐volume method is discussed in detail and particular attention is paid to accurately representing the complex irregular computational domain. The lower Yellow River reach from Huayuankou to Jiahetan is a typical meandering river. The generation of the computational mesh, which is used to simulate the flood, is affected by the distribution of water works in the river channel. The spatial information about the two Yellow River levee, the protecting dykes, and those roads that are obviously higher than the ground, need to be used to generate the computational mesh. As a result these dykes and roads locate the element interfaces of the computational mesh. In the model the finite‐volume method is used to solve the shallow‐wave equations, and the Osher scheme of the empirical function is used to calculate the flux through the interface between the neighbouring elements. The finite‐volume method has the advantage of using computational domain with complex geometry, and the Osher scheme is a method based on characteristic theory and is a monotone upwind numerical scheme with high resolution. The flood event with peak discharge of 15 300 m3/s, occurring in the period from 30 July to 10 August 1982, is simulated. The estimated result indicates that the simulation method is good for routing the flood in a region with complex geometry. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
19.
有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容。近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数。本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性。 相似文献
20.
The tsunami event generated by the great Sumatra–Andaman earthquake on 26 December 2004 was simulated with the recently developed model TsunAWI. The model is based on the finite element method, which allows for a very flexible discretization of the model domain. This is demonstrated by a triangulation of the whole Indian Ocean with a resolution of about 14 km in the deep ocean but a considerably higher resolution of about 500 m in the coastal area. A special focus is put on the Banda Aceh region in the Northern tip of Sumatra. This area was heavily hit by the tsunami and the highest resolution in this area is about 40 m in order to include inundation processes in the model simulation. We compare model results to tide gauge data from all around the Indian Ocean, to satellite altimetry, and field measurements of flow depth in selected locations of the Aceh region. Furthermore, we compare the model results of TsunAWI to the results of a nested grid model (TUNAMI-N3) with the same initial conditions and identical bathymetry and topography in the Aceh region. It turns out that TsunAWI gives accurate estimates of arrival times in distant locations and in the same mesh gives good inundation results when compared to field measurements and nested grid results. 相似文献