双重孔隙流体饱和介质弹性波散射二维IBIEM模拟

基于平面波势函数,采用间接边界积分方程法(indirect boundary integral equation method, IBIEM)研究了双孔隙流体饱和介质中弹性波入射下二维孔洞的散射特性。推导得到了双孔隙介质中全空间二维线源动力格林函数,并给出了各散射波的位移场和应力场。在数值精度验证的基础上,以双孔隙二维饱和全空间中孔洞为例,解决了平面P、SV波入射下的地震波散射问题。数值结果表明:双重孔隙介质中的位移幅值、环向应力幅值、孔隙压力变化规律与不同入射波形,入射频率,孔隙率和边界排水条件密切相关,位移幅值在低频(无量纲频率η≤2)入射时出现峰值。环向应力幅值与干土条件相比更为复杂,基质孔压与裂缝孔压的存在增大了双重孔隙饱和介质的能量效应,总体震动趋势大于干土条件,环向应力放大可达62%。

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS WITH NON-SEPARATED TYPE INTEGRAL BOUNDARY CONDITIONS

In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.

THE BOUNDARY INTEGRAL METHOD FOR THE HELMHOLTZ EQUATION WITH CRACKS INSIDE A BOUNDED DOMAIN

We consider a kind of scattering problem by a crack F that is buried in a bounded domain D, and we put a point source inside the domain D. This leads to a mixed boundary value problem to the Helmholtz equation in the domain D with a crack Г. Both sides of the crack F are given Dirichlet-impedance boundary conditions, and different boundary condition （Dirichlet, Neumann or Impedance boundary condition） is set on the boundary of D. Applying potential theory, the problem can be reformulated as a system of boundary integral equations. We establish the existence and uniqueness of the solution to the system by using the Fredholm theory.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

BOUNDARY INTEGRAL EQUATIONS FOR THE BENDING PROBLEM OF PLATES ON TWO-PARAMETER FOUNDATION

By means of Fourier integral transformation of generalized function, the fundamental solution for the bending problem of plates on two-parameter foundation is derived in this paper, and the fundamental solution is expanded into a uniformly convergent series. On the basis of the above work, two boundary integral equations which are suitable to arbitrary shapes and arbitrary boundary conditions are established by means of the Rayleigh-Green identity. The content of the paper provides the powerful theories for the application of BEM in this problem.

A Direct Implementation of a Modified Boundary Integral Formulation for the Extended Fisher-Kolmogorov Equation

This study is concerned with the numerical approximation of the extended Fisher-Kolmogorov equation with a modified boundary integral method. A key aspect of this formulation is that it relaxes the domain-driven approach of a typical boundary element (BEM) technique. While its discretization keeps faith with the second order accurate BEM formulation, its implementation is element-based. This leads to a local solution of all integral equation and their final assembly into a slender and banded coefficient matrix which is far easier to manipulate numerically. This outcome is much better than working with BEM’s fully populated coefficient matrices resulting from a numerical encounter with the problem domain especially for nonlinear, transient, and heterogeneous problems. Faithful results of high accuracy are achieved when the results obtained herein are compared with those available in literature.

2. 5-D RESISTIVITY TOMOGRAPHY USINGBOUNDARY INTEGRAL EQUATIONS

DC Resistivity Tomography is a non-linear inversion problem. So far there are mainly two kinds of inversion methods, based on the finite-element method and alpha centers method. In this paper, the disadvantages of these two kinds of methods were analysed,and a new method of forward modeling and inversion (Tomography) based on boundary integral equations was proposed. This scheme successfuly overcomes the difficulties of the two formarly methods. It isn’t necessary to use the linearization approximation and calculate the Jacobi matrix. Numerical modeling results given in this paper showed that the computation speed of our method is fast, and there is not any special requirement for initial model, and satisfying results of tomography can be obtained in the case of great contrast of conductivity. So it has wide applications.

ON THE UNIQUENESS OF BOUNDARY INTEGRAL EQUATION FOR THE EXTERIOR HELMHOLTZ PROBLEM

From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the Sommerfeld's condition at the infinity is changed so that it is suitable not only for the radiative wave but also for the absorptive wave when we use the boundary integral equation to describe the exterior Helmholtz problem. There fore, the total energy of the system is conservative. The mathematical dealings to guarantee the uniqueness are discussed based upon this explanation.

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regularization to the direct boundary integal equations （DBIEs）. The character is elimination of singularities by using the simple solutions. However, up to now the research of regularization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value （CPV） and Hadamard-finite-part （HFP） integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.

Solution of stress intensity factors of multiple cracks in plane elasticity with eigen COD formulation of boundary integral equation

The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of boundary integral equation is proposed in this paper, together with the solution procedures for multiple crack problems in plane elasticity. With the proposed approach, the multiple crack problems can be solved with the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix as that in the numerical Green’s function （NGF） approach but without the trouble to determine the complementary solutions since the standard boundary element discretization on the crack surface is no longer required with the proposed approach. Some numerical examples computing the stress intensity factors are presented and compared with those in literature to show the accuracy and the effectiveness of the proposed approach.

RESEARCH ON THE COMPANION SOLUTION FOR A THIN PLATE IN THE MESHLESS LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.

EQUIVALENT BOUNDARY INTEGRAL EQUATIONS WITH INDIRECT UNKNOWNS FOR THIN ELASTIC PLATE BENDING THEORY

Equivalent Boundary Integral Equations (EBIE) with indirect unknowns for thin elastic plate bending theory, which is equivalent to the original boundary value problem, is established rigorously by mathematical technique of non-analytic continuation and is fully proved by means of the variational principle. The previous three kinds of boundary integral equations with indirect unknowns are discussed thoroughly and it is shown that all previous results are not EBIE.

Entropy Formulation for Triply Nonlinear Degenerate Elliptic-Parabolic-Hyperbolic Equation with Zero-Flux Boundary Condition

In this note, we investigated existence and uniqueness of entropy solution for triply nonlinear degenerate parabolic problem with zero-flux boundary condition. Accordingly to the case of doubly nonlinear degenerate parabolic hyperbolic equation, we propose a generalization of entropy formulation and prove existence and uniqueness result without any structure condition.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

Two-dimensional polynomial eigenstrain formulation of boundary integral equation with numerical verification

The low-order polynomial-distributed eigenstrain formulation of the boundary integral equation （BIE） and the corresponding definition of the Eshelby tensors are proposed for the elliptical inhomogeneities in two-dimensional elastic media. Taking the results of the traditional subdomain boundary element method （BEM） as the control, the effectiveness of the present algorithm is verified for the elastic media with a single elliptical inhomogeneity. With the present computational model and algorithm, significant improvements are achieved in terms of the efficiency as compared with the traditional BEM and the accuracy as compared with the constant eigenstrain formulation of the BIE.

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.

BOUNDARY INTEGRAL EQUATIONS OF UNIQUE SOLUTIONS IN ELASTICITY

The properties of the fundamental solution are derived in linear elastostatics. These properties are used to show that the conventional displacement and traction boundary integral equations yield non-unique displacement solutions in a traction boundary value problem. The condition for the existence of unique displacement solutions is proposed for the traction boundary value problem. The degrees of freedom of the displacement solution are removed by the condition to obtain the boundary integral equations of unique solutions for the traction boundary value problems. Numerical example is presented to demonstrate the accuracy and efficiency of the present equations.

SINGULAR INTEGRAL EQUATIONS AND BOUNDARY ELEMENT METHOD OF CRACKS IN THERMALLY STRESSED PLANAR SOLIDS

Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method （LBIEM）, singularities in the local boundary integrals need to be treated specially. In the current paper, local integral equations are adopted for the nodes inside the domain and moving least square approximation （MLSA） for the nodes on the global boundary, thus singularities will not occur in the new al- gorithm. At the same time, approximation errors of boundary integrals are reduced significantly. As applications and numerical tests, Laplace equation and Helmholtz equation problems are considered and excellent numerical results are obtained. Furthermore, when solving the Helmholtz problems, the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions. Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.

Boundary integral equations for dynamic rupture propagation on vertical complex fault system in half-space:Theory

The boundary integral equation method （BIEM） is now widely used in numerical studies on earthquake rupture dynamics, and is proved to be a powerful tool to deal with problems on complex fault system. However, since this method heavily lies on the specific forms of Green＇s function and only the Green＇s function in full-space has a closed analytic expression, it is usually limited to a full-space medium. In this study, as a first step to extend this method to an arbitrary complex fault system in half-space, the boundary integral equations （BIEs） for dynamic strike-slip on vertical complex fault system in half-space are derived based on exact Green＇s function for isotropic and homogeneous half-space. Effect of the geometry of the complex fault system are dealt with carefully. Final BIEs is composed of two parts： contribution from full-space, which has been thoroughly investigated by Aochi and his co-workers by using the Green＇s function for full-space, and that from free surface, which is studied in detail in this study.