Wavelet Solution of Plane Elasticity Problem in the Upper Half-plane

The plane elasticity problem includes plane strain problem and plane stress problem which are widely applied in mechanics and engineering. In this article, we first reduce the plane elasticity problem in the upper half-plane into natural boundary integral equation and then apply wavelet-Galerkin method to deal with the numerical solution of the natural boundary integral equation. The test and trial functions used are the scaling basis functions of Shannon wavelet. In our fast algorithm, the computational formulae of entries of the stiffness matrix yield simple close-form and only 3 K entries need to be computed for one 4 K ‐ 4 K stiffness matrix.     

Numerical solution of hypersingular equation using recursive wavelet on invariant set

In this paper, we construct the Chebyshev recursive wavelets on a unit interval of the first kind, the second kind and their corresponding weight functions. We apply wavelet collocation method to solve the natural boundary integral equation of the harmonic equation on the lower half-plane numerically. It is convenient and accurate to generate the stiffness matrix. Two numerical examples are presented. It is shown that the stiffness matrix is highly sparse when the order of the stiffness matrix becomes large. Current method allows choosing an appropriate weight function to increase the convergence rate and accuracy of the numerical results.     

Compatible wavelet-Galerkin method to solve stokes interior problem with circular boundary

This article addresses Neumann boundary value interior problem of Stokes equations with circular boundary. By using natural boundary element method, the Stokes interior problem is reduced into an equivalent natural integral equation with a hyper-singular kernel, which is viewed as Hadamard finite part. Based on trigonometric wavelet functions, the compatible wavelet space is constructed so that it can serve as Galerkin trial function space. In proposed compatible wavelet-Galerkin method, the simple and accurate computational formulae of the entries in stiffness matrix are obtained by singularity removal technique. It is also proved that the stiffness matrix is almost a block diagonal matrix, and its diagonal sub-blocks all are both symmetric and circulant submatrices. These good properties indicate that a 2 ^J+3 × 2 ^J+3 stiffness matrix can be determined only by its 2 ^J + 3J + 1 entries. It greatly decreases the computational complexity. Some error estimates for the compatible wavelet-Galerkin projection solutions are established. Finally, several numerical examples are given to demonstrate the validity of the proposed approach.     

Collocation method for the natural boundary integral equation

In this work, we construct the Legendre wavelet and apply it to investigate the numerical solution of the natural boundary integral equation of the Laplace equation in the upper half-plane by the collocation method. In our algorithm the coefficient matrix of the linear algebraic system is sparse when the order of the matrix is large. Two test examples show that our algorithm yields very accurate results at less computational cost.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

Hybrid boundary element formulation for acoustic wave propagation

The so‐called hybrid stress boundary element method (HSBEM) is introduced in a frequency domain formulation for the computation of acoustic radiation and scattering in closed and in finite domains. Different from other boundary element formulations, the HSBEM is based on an extended Hellinger‐Reissner variational principle and leads to a Hermitian, frequency‐dependent stiffness equation. Due to this, the method is very well suited for treating fluid structure interaction problems since the effort for the coupling the structure, discretized by a finite elements, and the fluid, discretized by the HSBEM is strongly reduced. To arrive at a boundary integral formulation, the field variables are separated into boundary variables, which are approximated by piecewise polynomial functions, and domain variables, which are approximated by a superposition of singular fundamental solutions weighed by source strength. This approximation cancels the domain integral over the equation of motion in the hybrid principle and leads to a boundary integral formulation, incorporating singular integrals. Comparing to previous results published by the authors, new considerations concerning the interpretation of singular contributions in the stiffness matrix for exterior domain problems are communicated here.     

Circulant Wavelet Preconditioners for Solving Elliptic Differential Equations and Boundary Integral Equations

In this paper is discussed solving an elliptic equation and a boundary integral equation of the second kind by representation of compactly supported wavelets. By using wavelet bases and the Galerkin method for these equations, we obtain a stiff sparse matrix that can be ill-conditioned. Therefore, we have to introduce an operator which maps every sparse matrix to a circulant sparse matrix. This class of circulant matrices is a class of preconditioners in a Banach space. Based on having some properties in the spectral theory for this class of matrices, we conclude that the circulant matrices are a good class of preconditioners for solving these equations. We called them circulant wavelet preconditioners (CWP). Therefore, a class of algorithms is introduced for rapid numerical application.     

A wavelets method for plane elasticity problem

A Neumann boundary value problem of plane elasticity problem in the exterior circular domain is reduced into an equivalent natural boundary integral equation and a Poisson integral formula with the DtN method. Using the trigonometric wavelets and Galerkin method, we obtain a fast numerical method for the natural boundary integral equation which has an unique solution in the quotient space. We decompose the stiffness matrix in our numerical method into four circulant and symmetrical or antisymmetrical submatrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT) instead of the inverse matrix. Examples are given for demonstrating our method has good accuracy of our method even though the exact solution is almost singular.     

On a Sparse Representation of an n-Dimensional Laplacian in Wavelet Coordinates

Important parts of adaptive wavelet methods are well-conditioned wavelet stiffness matrices and an efficient approximate multiplication of quasi-sparse stiffness matrices with vectors in wavelet coordinates. Therefore it is useful to develop a well-conditioned wavelet basis with respect to which both the mass and stiffness matrices are sparse in the sense that the number of nonzero elements in each column is bounded by a constant. Consequently, the stiffness matrix corresponding to the n-dimensional Laplacian in the tensor product wavelet basis is also sparse. Then a matrix–vector multiplication can be performed exactly with linear complexity. In this paper, we construct a wavelet basis based on Hermite cubic splines with respect to which both the mass matrix and the stiffness matrix corresponding to a one-dimensional Poisson equation are sparse. Moreover, a proposed basis is well-conditioned on low decomposition levels. Small condition numbers for low decomposition levels and a sparse structure of stiffness matrices are kept for any well-conditioned second order partial differential equations with constant coefficients; furthermore, they are independent of the space dimension.     

求解对流扩散方程的Haar小波方法

本文用Haar小波求解对流扩散方程，将满足初始和边界条件的常系数偏微分方程简化为较简单的代数方程组进行求解．实例说明了这种方法具有收敛速度快和计算容易的特点，同时又避免了用Daubechies小波求解微分方程需要计算相关系数的麻烦．本文所使用的方法可以求解一般的微（积）分方程．     

Stable three-point wavelet bases on general meshes

Summary. Recently, we introduced a wavelet basis on general, possibly locally refined linear finite element spaces. Each wavelet is a linear combination of three nodal basis functions, independently of the number of space dimensions. In the present paper, we show -stability of this basis for a range of , that in any case includes , which means that the corresponding additive Schwarz preconditioner is optimal for second order problems. Furthermore, we generalize the construction of the wavelet basis to manifolds. We show that the wavelets have at least one-, and in areas where the manifold is smooth and the mesh is uniform even two vanishing moments. Because of these vanishing moments, apart from preconditioning, the basis can be used for compression purposes: For a class of integral equation problems, the stiffness matrix with respect to the wavelet basis will be close to a sparse one, in the sense that, a priori, it can be compressed to a sparse matrix without the order of convergence being reduced. Received November 6, 1996 / Revised version received June 30, 1997     

Biorthogonal wavelet approximation for the coupling of FEM-BEM

Summary. We apply multiscale methods to the coupling of finite and boundary element methods to solve an exterior Dirichlet boundary value problem for the two dimensional Poisson equation. Adopting biorthogonal wavelet matrix compression to the boundary terms with N degrees of freedom, we show that the resulting compression strategy fits the optimal convergence rate of the coupling Galerkin methods, while the number of nonzero entries in the corresponding stiffness matrices is considerably smaller than . Received December 3, 1999 / Revised version received September 22, 2000 / Published online December 18, 2001     

On a quadrature algorithm for the piecewise linear wavelet collocation applied to boundary integral equations

In this paper, we consider a piecewise linear collocation method for the solution of a pseudo‐differential equation of order r=0, ?1 over a closed and smooth boundary manifold. The trial space is the space of all continuous and piecewise linear functions defined over a uniform triangular grid and the collocation points are the grid points. For the wavelet basis in the trial space we choose the three‐point hierarchical basis together with a slight modification near the boundary points of the global patches of parametrization. We choose linear combinations of Dirac delta functionals as wavelet basis in the space of test functionals. For the corresponding wavelet algorithm, we show that the parametrization can be approximated by low‐order piecewise polynomial interpolation and that the integrals in the stiffness matrix can be computed by quadrature, where the quadrature rules are composite rules of simple low‐order quadratures. The whole algorithm for the assembling of the matrix requires no more than O(N [logN]³) arithmetic operations, and the error of the collocation approximation, including the compression, the approximative parametrization, and the quadratures, is less than O(N^?(2?r)/2). Note that, in contrast to well‐known algorithms by Petersdorff, Schwab, and Schneider, only a finite degree of smoothness is required. In contrast to an algorithm of Ehrich and Rathsfeld, no multiplicative splitting of the kernel function is required. Beside the usual mapping properties of the integral operator in low order Sobolev spaces, estimates of Calderón–Zygmund type are the only assumptions on the kernel function. Copyright © 2003 John Wiley & Sons, Ltd.     

A preconditioning strategy for boundary element Galerkin methods

The Dirichlet and Neumann problems for the Laplacian are reformulated in the usual way as boundary integral equations of the first kind with symmetric kernels. These integral equations are solved using Galerkin's method with piecewise-constant and piecewise-linear boundary elements, respectively. In both cases, the stiffness matrix is symmetric and positive-definite, and has a condition number of order N, the number of degrees of freedom. By contrast, the condition number of the product of the two stiffness matrices is bounded independently of N. Hence, we can use the Neumann stiffness matrix to precondition the Dirichlet stiffness matrix, and vice versa. © 1997 John Wiley & Sons, Inc.     

Boussinesq indentation of a transversely isotropic half-space reinforced by a buried inextensible membrane

The normal indentation of a rigid circular disk into the surface of a transversely isotropic half-space reinforced by a buried inextensible thin film is addressed. By virtue of a displacement potential function and the Hankel transform, the governing equations of this axisymmetric mixed boundary value problem are represented as a dual integral equation, which is subsequently reduced to a Fredholm integral equation of the second kind. Two important results of the contact stress distribution beneath the disk region as well as the equivalent stiffness of the system are expressed in terms of the solution of the Fredholm integral equation. When the membrane is located on the surface or at the remote boundary, exact closed-form solutions are presented. For the limiting case of an isotropic half-space the results are verified with those available in the literature. As a special case, the elastic fields of a reinforced transversely isotropic half-space under the action of surface axisymmetric patch loads are also given. The effects of anisotropy, embedment depth of the membrane, and material incompressibility on both the contact stress and the normal stiffness factor are depicted in some plots.     

Boundary integral equation method of higher computational accuracy

The boundary integral equation method presented in the paper features the following: (1) no singular kernels, strong or weak, are involved, and computationally no local “element” approximations are needed; (2) the integral equations are well conditioned, including the cases of bounded and multiply connected regions, and no iterative approximations are involved; (3) no post-solution differentiation is involved. These features provide for a higher computational efficiency. The method solves in full a number of engineering problems, and can be used for the stiffness matrix formulation in more complex situations.     

Numerical Solution of Parabolic Problems with Nonlocal Boundary Conditions

In this article, the one-dimensional parabolic equation with three types of integral nonlocal boundary conditions is approximated by the implicit Euler finite difference scheme. Stability analysis is done in the maximum norm and it is proved that the radius of the stability region and the stiffness of the discrete scheme depends on the signs of coefficients in the nonlocal boundary condition. The known stability results are improved. In the case of a plain integral boundary condition, the conditional convergence rate is proved and the regularization relation between discrete time and space steps is proposed. The accuracy of the obtained estimates is illustrated by results of numerical experiments.     

边界曲线积分方程的小波解法

The boundary measure method is applied to transfer the form of the integral equation in order to use the collocation method or Galerkin method. A simple way to computer the coefficients of the wavelet series is also introduced. The way presented in this paper can be used to solve PDE problem in the two dimension region with any form of boundary.     

Prewavelet analysis of the heat equation

Summary. We consider the heat equation in a smooth domain of R with Dirichlet and Neumann boundary conditions. It is solved by using its integral formulation with double-layer potentials, where the unknown , the jump of the solution through the boundary, belongs to an anisotropic Sobolev space. We approximate by the Galerkin method and use a prewavelet basis on , which characterizes the anisotropic space. The use of prewavelets allows to compress the stiffness matrix from to when N is the size of the matrix, and the condition number of the compressed matrix is uniformly bounded as the initial one in the prewavelet basis. Finally we show that the compressed scheme converges as fast as the Galerkin one, even for the Dirichlet problem which does not admit a coercive variational formulation. Received April 16, 1999 / Published online August 2, 2000     

Wavelet approximations for first kind boundary integral equations on polygons

Summary. An elliptic boundary value problem in the interior or exterior of a polygon is transformed into an equivalent first kind boundary integral equation. Its Galerkin discretization with degrees of freedom on the boundary with spline wavelets as basis functions is analyzed. A truncation strategy is presented which allows to reduce the number of nonzero elements in the stiffness matrix from to entries. The condition numbers are bounded independently of the meshwidth. It is proved that the compressed scheme thus obtained yields in operations approximate solutions with the same asymptotic convergence rates as the full Galerkin scheme in the boundary energy norm as well as in interior points. Numerical examples show the asymptotic error analysis to be valid already for moderate values of . Received March 12, 1994 / Revised version received January 9, 1995