共查询到20条相似文献，搜索用时 31 毫秒
1.
Chi Hsin Lin Vincent W. Lee Maria I. Todorovska Mihailo D. Trifunac 《Soil Dynamics and Earthquake Engineering》2010
The zero-stress boundary conditions at the surface of the half-space in the presence of surface and sub-surface cavities for in-plane, incident cylindrical P- and SV-waves have always posed challenging problems. The outgoing cylindrical P- and SV-waves can be represented by Hankel functions of radial distance coupled with the sine and cosine functions of angle. Together, at the half-space surface the P- and SV-wave functions are not orthogonal over the semi-infinite radial distance from 0 to infinity. Thus, to simultaneously satisfy the zero in-plane, normal, and shear stresses, an approximation of the geometry is often made. This paper presents an analytical formulation of the boundary-valued problem, where the Hankel wave functions are expressed in integral form, changing the representation from cylindrical to rectangular coordinates, so that the zero-stress boundary conditions at the half-space surface can be applied in a more straightforward way. 相似文献
2.
Vladimir Priklonsky Anvar Shukurov Dmitry Sokoloff Andrew Soward 《地球物理与天体物理流体动力学》2013,107(1-2):97-114
Abstract We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion. 相似文献
3.
Oversampling techniques are often used in porous media simulations to achieve high accuracy in multiscale simulations. These methods reduce the effect of artificial boundary conditions that are imposed in computing local quantities, such as upscaled permeabilities or basis functions. In the problems without scale separation and strong non-local effects, the oversampling region is taken to be the entire domain. The basis functions are computed using single-phase flow solutions which are further used in dynamic two-phase simulations. The standard oversampling approaches employ generic global boundary conditions which are not associated with actual flow boundary conditions. In this paper, we propose a flow based oversampling method where the actual two-phase flow boundary conditions are used in constructing oversampling auxiliary functions. Our numerical results show that the flow based oversampling approach is several times more accurate than the standard oversampling method. We provide partial theoretical explanation for these numerical observations. 相似文献
4.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度. 相似文献
5.
This paper is concerned with the convergence of the Galerkin finite element method applied to a groundwater flow problem containing a borehole, with special reference to quadrature effects and the accuracy of the solution. It is shown that there exists an optimal quadrature rule for every choice of piecewise polynomial basis functions. Another interesting result proved here is that, in a direct application of the method the accuracy is very nearly independent of the degree of the polynomial basis functions, but strongly dependent on the distance of the borehole from the boundary if this is small. 相似文献
6.
T. Senjuntichai S. Mani R.K.N.D. Rajapakse 《Soil Dynamics and Earthquake Engineering》2006,26(6-7):626-636
This paper considers time-harmonic vertical vibration of an axisymmetric rigid foundation embedded in a homogeneous poroelastic soil. The soil domain is represented by a homogeneous poroelastic half space that is governed by Biot's theory of poroelastodynamics. The foundation is subjected to a time-harmonic vertical load and is perfectly bonded to the surrounding half space. The contact surface can be either fully permeable or impermeable. The dynamic interaction problem is solved by employing an indirect boundary integral equation method. The kernel functions of the integral equation are the influence functions corresponding to vertical and radial ring loads, and a ring fluid source applied in the interior of a homogeneous poroelastic half space. Analytical techniques are used to derive the solution for influence functions. The indirect boundary integral equation is solved by using numerical quadrature. Selected numerical results for vertical impedance of rigid foundations are presented to demonstrate the influence of poroelastic effect, foundation geometry, hydraulic boundary condition along the contact surface and frequency of excitation. 相似文献
7.
The (simplified) Backus’ Problem (BP) consists in finding a harmonic function u on the domain exterior to the three dimensional
unit sphere S, such that u tends to zero at infinity and the norm of the gradient of u takes prescribed values g on S. Except
for a change of sign, the solution is not unique in general. However, there is uniqueness of solutions in the class of functions
with the additional property that the radial component of the gradient of u on S is nonpositive. This is the geodetically
relevant case. If a solution u with this property exists, then u is the maximal solution of the problem (and -u the minimal
one). In this paper we propose a method of successive approximations to get this particular solution of BP and prove the convergence
for functions g close to a constant function. 相似文献
8.
2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method 
下载免费PDF全文
下载免费PDF全文 Yuan Yuan Jianke Qiang Jingtian Tang Zhengyong Ren Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779
To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion. 相似文献
9.
The magnetospheric magnetic field may be conveniently described by two scalar functions (α, β), known as the Euler potentials. They are not uniquely defined, and they may be difficult to derive for configuration more complex than a simple dipole. We propose here a simple numerical method to compute one possible pair (α, β). In magnetospheric regions of closed field lines, α can be chosen as a function of the tube volume of unit magnetic flux. The method can be applied to a wide class of magnetic fields which describe the magnetospheric domain of closed field lines and the conjugated ionosphere. Here, it is used with the T87 Tsyganenko model. The results coincide with the dipolar potentials at close distances from the Earth. At larger distances, they display an increasing distortion with the radial distance (or the invariant latitude in the ionosphere) and the magnetic activity. In the magnetosphere, the contours of α and β are stretched towards the nightside. In the ionosphere, they also extend towards the nightside and present major distortions in a narrow ring at the polar cap boundary, which mapsdistant boundary layers in the magnetosphere. 相似文献
10.
In this paper we derive some similarity solutions of a nonlinear equation associated with a free boundary problem arising
in the shallow-water approximation in glaciology. In addition we present a classical potential symmetry analysis of this second-order
nonlinear degenerate parabolic equation related to non-Newtonian ice sheet dynamics in the isothermal case. After obtaining
a general result connecting the thickness function of the ice sheet and the solution of the nonlinear equation (without any
unilateral formulation), a particular example of a similarity solution to a problem formulated with Cauchy boundary conditions
is described. This allows us to obtain several qualitative properties on the free moving boundary in the presence of an accumulation-ablation
function with realistic physical properties. 相似文献
11.
We discuss the inverse medium problem associated with the reconstruction of the heterogeneous material profile of a semi-infinite (layered) soil medium, directly in the time domain, based on the complete waveform response of the medium to interrogating waves. To tackle the inversion process, we use a partial-differential-equation-constrained optimization approach, supplemented with a time-dependent regularization scheme. We introduce an absorbing boundary to truncate the semi-infinite extent of the physical domain, and propose two schemes to refine the reconstructed profiles: the first is based on iteratively re-positioning the truncation boundary until convergence, and the second is based on optimizing the observation period, so as to exclude records with information beyond the truncation boundary. We present numerical results that attest to the efficacy of the proposed schemes in reconstructing sharp profiles of semi-infinite soil domains using both noise-free and noisy data, while in the presence of absorbing boundaries. 相似文献
12.
13.
贴体网格在地质数值模拟中具有广阔的应用前景,为解决贴体网格生成时边界离散问题,提出了最大长度准则和最大面积准则,把曲线逼近和曲面网格优化问题转化为数学优化问题,为求解该问题,提出了改进的单粒子优化算法.试验表明,最大长度准则和最大面积准则的优化效果好于常规方法;以改进的单粒子优化算法求解该问题时,计算效率是智能单粒子优化算法的30倍左右（节点量为200）,从而实现最大长度准则和最大面积准则在贴体网格生成中的应用.针对最大面积准则优化曲面网格不能控制网格步长的情况,提出了限定步长的网格优化算法,使网格步长合理化,并通过实例验证了该算法的有效性.研究成果提供了生成贴体网格时边界优化准则和求解方法,对今后复杂边界的贴体网格生成具有重要意义. 相似文献
14.
We introduce a nonlinear orthogonal matching pursuit (NOMP) for sparse calibration of subsurface flow models. Sparse calibration is a challenging problem as the unknowns are both the non-zero components of the solution and their associated weights. NOMP is a greedy algorithm that discovers at each iteration the most correlated basis function with the residual from a large pool of basis functions. The discovered basis (aka support) is augmented across the nonlinear iterations. Once a set of basis functions are selected, the solution is obtained by applying Tikhonov regularization. The proposed algorithm relies on stochastically approximated gradient using an iterative stochastic ensemble method (ISEM). In the current study, the search space is parameterized using an overcomplete dictionary of basis functions built using the K-SVD algorithm. The proposed algorithm is the first ensemble based algorithm that tackels the sparse nonlinear parameter estimation problem. 相似文献
15.
George Helffrich Bruno Faria João F.B.D. Fonseca Alexandra Lodge Satoshi Kaneshima 《Earth and Planetary Science Letters》2010,289(1-2):156-161
We report on a two-year seismic deployment in the Cape Verde Islands, one goal of which was to study the upper mantle to determine its structure under a hot spot that is stationary in the hot spot reference frame. We find from analysis of P-to-S receiver functions estimated from broadband seismic recordings that, within uncertainty, the time separation between the 410 and 660 km discontinuities is normal compared to radial earth models. Thus, to exist, even stationary hot spots do not require vertical thermal anomalies from deep melting sources anchored in the lower mantle or at the core–mantle boundary or their anomalies are narrower than ～ 250 km in the upper mantle. 相似文献
16.
Martinec and Grafarend (1997) have shown how the construction of Green's function in the Stokes boundary-value problem with gravity data distributed on an ellipsoid of revolution is approached in the O(e
0
2
)-approximation. They have also expressed the ellipsoidal Stokes function describing the effect of ellipticity of the boundary as a finite sum of elementary functions. We present an effective method of avoiding the singularity of spherical and the ellipsoidal Stokes functions, and also an analytical expression for the ellipsoidal Stokes integral around the computational point suitable for numerical solution. We give the numerical results of solving the ellipsoidal Stokes boundary-value problem and their difference with respect to the spherical Stoke boundary-value problem. 相似文献
17.
Bohuslav R??ek Jaroslava Plomerová Vladislav Babu?ka 《Studia Geophysica et Geodaetica》2012,56(1):107-140
Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise
has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms
of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and
in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term
recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed
from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying
deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The
joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and
group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives
nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle.
We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using
statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological
structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively
low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly
felsic. 相似文献
18.
A new methodology is proposed for the development of parameter-independent reduced models for transient groundwater flow models. The model reduction technique is based on Galerkin projection of a highly discretized model onto a subspace spanned by a small number of optimally chosen basis functions. We propose two greedy algorithms that iteratively select optimal parameter sets and snapshot times between the parameter space and the time domain in order to generate snapshots. The snapshots are used to build the Galerkin projection matrix, which covers the entire parameter space in the full model. We then apply the reduced subspace model to solve two inverse problems: a deterministic inverse problem and a Bayesian inverse problem with a Markov Chain Monte Carlo (MCMC) method. The proposed methodology is validated with a conceptual one-dimensional groundwater flow model. We then apply the methodology to a basin-scale, conceptual aquifer in the Oristano plain of Sardinia, Italy. Using the methodology, the full model governed by 29,197 ordinary differential equations is reduced by two to three orders of magnitude, resulting in a drastic reduction in computational requirements. 相似文献
19.
Eulerian–Lagrangian localized adjoint methods (ELLAMs) provide a general approach to the solution of advection-dominated advection–diffusion equations allowing large time steps while maintaining good accuracy. Moreover, the methods can treat systematically any type of boundary condition and are mass conservative. However, all ELLAMs developed so far suffer from non-physical oscillations and are usually implemented on structured grids. In this paper, we propose a finite volume ELLAM which incorporates a novel correction step rendering the method monotone while maintaining conservation of mass. The method has been implemented on fully unstructured meshes in two space dimensions. Numerical results demonstrate the applicability of the method for problems with highly non-uniform flow fields arising from heterogeneous porous media. 相似文献
20.
M. Rahman 《Advances in water resources》1981,4(1):9-19
The forces and overturning moments exerted by second order waves on large vertical circular cylinders are analysed. The mathematical equations governing the physical system are the three-dimensional Laplace's equation satisfied by the velocity potential ?(x,y,z,t) and the boundary conditions, namely the dynamic boundary condition which is obtained from the Bernoulli's equation, kinematic boundary condition, radiation condition, bottom boundary condition and the zero radial velocity condition on the surface of the cylinder. The non-linearity of the mathematical problem is evidenced in the free surface boundary conditions viz. dynamic and kinematic boundary conditions. Analytical solutions are obtained using the perturbation technique. These solutions are compared with various experimental data. The comparison shows favourable agreement between the theory and the experimental results. 相似文献