Fast solvers with block‐diagonal preconditioners for linear FEM–BEM coupling

The purpose of this paper is to present optimal preconditioned iterative methods to solve indefinite linear systems of equations arising from symmetric coupling of finite elements and boundary elements. This is a block‐diagonal preconditioner together with a conjugate residual method and a preconditioned inner–outer iteration. We prove the efficiency of these methods by showing that the number of iterations to preserve a given accuracy is bounded independent of the number of unknowns. Numerical examples underline the efficiency of these methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Multiscale preconditioning for the coupling of FEM–BEM

We apply multiscale methods to the coupling of finite and boundary element methods to solve an exterior two‐dimensional Laplacian. The matrices belonging to the boundary terms of the coupled FEM–BEM system are compressed by using biorthogonal wavelet bases developed from A. Cohen, I. Daubechies and J.‐C. Feauveau (Comm. Proc. Appl. Math. 1992; 45 :485). The coupling yields a linear equation system which corresponds to a saddle point problem. As favourable solver, the Bramble–Pasciak–CG (Math. Comp. 1988; 50 :1) is utilized. A suitable preconditioner is developed by combining the BPX (Math. Comp. 1990; 55 :1) with the wavelet preconditioning (Numer. Math. 1992; 63 :315). Through numerical experiments we provide results which corroborate the theory of the present paper. Copyright © 2002 John Wiley & Sons, Ltd.     

Two Dimensional Interface Problems for Elliptic Equations

We study the structure of solutions to the interface problems for second order quasi-linear elliptic partial differential equations in two dimensional space. We prove that each weak solution can be decomposed into two parts near singular points, a finite sum of functions in the form of cr^α log^m rφ(θ) and a regular one w. The coefficients c and the C^{1,α} norm of w depend on the H¹-norm and the C^{º,α}-norm of the solution, and the equation only.     

Three Dimensional Interface Problems for Elliptic Equations

The author studies the structure of solutions to the interface problems for second order linear elliptic partial differential equations in three space dimension.The set of singular points consists of some singular lines and some isolated singular points.It is proved that near a singular line or a singular point,each weak solution can be decomposed into two parts,a singular part and a regular part.The singular parts are some finite sum of particular solutions to some simpler equations,and the regular parts are bounded in some norms,which are slightly weaker than that in the Sobolev space H~2.     

Uniform Convergence of Adaptive Multigrid Methods for Elliptic Problems and Maxwell's Equations

We consider the convergence theory of adaptive multigrid methods for second-order elliptic problems and Maxwell's equations. The multigrid algorithm only performs pointwise Gauss-Seidel relaxations on new degrees of freedom and their "immediate" neighbors. In the context of lowest order conforming finite element approximations, we present a unified proof for the convergence of adaptive multigrid V-cycle algorithms. The theory applies to any hierarchical tetrahedral meshes with uniformly bounded shape-regularity measures. The convergence rates for both problems are uniform with respect to the number of mesh levels and the number of degrees of freedom. We demonstrate our convergence theory by two numerical experiments.     

Fourier–Finite-element Approximation of Elliptic Interface Problems in Axisymmetric Domains

The paper deals with the Fourier-finite-element method (FFEM), which combines the approximate Fourier method with the finite-element method, and its application to Poisson-like equations −p̂Δ₃û = f̂ in three-dimensional axisymmetric domains Ωˆ. Here, pˆ is a piecewise constant coefficient having a jump at some axisymmetric interface. Special emphasis is given to estimates of the Fourier-finite-element error in the Sobolev space H¹(Ωˆ), if the interface is smooth or if it meets the boundary of Ωˆ at some edge. In general, the solution û contains a singularity at the interface, which is described by a tensor product representation and treated numerically by appropriate mesh grading in the meridian plane of Ωˆ. The rate of convergence of the combined approximation in H¹(Ωˆ) is proved to be 𝒪(h+N⁻¹) (h, N: the parameters of the finite-element- and Fourier-approximation, with h→0, N→∞). The theoretical results are confirmed by numerical experiments.     

Interface Problems for Elliptic Systems

Interface problems for elliptic systems of second order partial differential equations are studied. The main result is that the solution in the neighborhood of the singular point can be divided into two parts one of which is a solution to the homogeneous system with constant coefficients, and the other one possesses higher regularity.     

A Fast Algorithm for Two-Dimensional Elliptic Problems

In this paper, we extend the work of Daripa et al. [14–16,7] to a larger class of elliptic problems in a variety of domains. In particular, analysis-based fast algorithms to solve inhomogeneous elliptic equations of three different types in three different two-dimensional domains are derived. Dirichlet, Neumann and mixed boundary value problems are treated in all these cases. Three different domains considered are: (i) interior of a circle, (ii) exterior of a circle, and (iii) circular annulus. Three different types of elliptic problems considered are: (i) Poisson equation, (ii) Helmholtz equation (oscillatory case), and (iii) Helmholtz equation (monotone case). These algorithms are derived from an exact formula for the solution of a large class of elliptic equations (where the coefficients of the equation do not depend on the polar angle when written in polar coordinates) based on Fourier series expansion and a one-dimensional ordinary differential equation. The performance of these algorithms is illustrated for several of these problems. Numerical results are presented.     

A Decomposition Theorem for the Solutions to the Interface Problems of Quasi-Linear Elliptic Equations

Interface problems for second order quasi-linear elliptic partial differential equations in a two-dimensional space are studied. We prove that each weak solution can be decomposed into two parts near singular points, one of which is a finite sum of functions of the form cr^a log^m rφ(θ), where the coefficients c depend on the H^1-norm of the solution, the C^(0,δ) -norm of the solution, and the equation only; and the other one of which is a regular one, the norm of which is also estimated.     

一种半线性椭圆型方程的解的可去奇点

本文处理了一种半线性椭圆型方程解的可去奇点问题,得到方程的弱解与一个定义在RN上的连续函数几乎处处相等的结论.     

Boundary-value Problems for Integro-differential Equations of Elliptic Type

In this paper we study the existence of solutions to the Dirichlet problem for a class of integro-differential equations of elliptic type by using the weakly continuous method.     

Adaptive Wavelet BEM for Boundary Integral Equations: Theory and Numerical Experiments

We are concerned with the numerical treatment of boundary integral equations by the adaptive wavelet boundary element method. In particular, we consider the second kind Fredholm integral equation for the double layer potential operator on patchwise smooth manifolds contained in ?³. The corresponding operator equations are treated by adaptive implementations that are in complete accordance with the underlying theory. The numerical experiments demonstrate that adaptive methods really pay off in this setting. The observed convergence rates fit together very well with the theoretical predictions based on the Besov regularity of the exact solution.     

椭圆型方程的内层和边界层渐近解

本文讨论了一类半线性椭圆型方程边值问题.利用微分不等式理论,研究了边值问题内层和边界层解的存在性和渐近性态.     

Resonance Problems for a Class of Singular Elliptic Equations

In this article we study the resonance problems for Hardy-Sobolev operator $ L_{\mu }:= -\Delta _{p}-\fraca {(\mu }{|x|^{p})}, 0\le \mu \lt {\fraca {(N-p}{p} )}\,^{p} $ with unbounded nonlinearity. Existence of weak solutions is established using the minimax theorems in critical point theory.     

Free Boundary Value Problems for Abstract Elliptic Equations and Applications

The free boundary value problems for elliptic differential-operator equations are studied. Several conditions for the uniform maximal regularity with respect to boundary parameters and the Fredholmness in abstract L _p -spaces are given. In application, the nonlocal free boundary problems for finite or infinite systems of elliptic and anisotropic type equations are studied.     

Solution of Functional Equations Related to Elliptic Functions

Functional equations of the form f(x + y)g(x ? y) = Σ _j=1 ⁿ α _j (x)β _j (y) as well as of the form f₁(x + z)f₂(y + z)f₃(x + y ? z) = Σ _j=1 ^m φ _j (x, y)ψ _j (z) are solved for unknown entire functions f, g,α _j , β _j : ? → ? and f₁, f₂, f₃, ψ _j : ? → ?, φ _j : ?² → ? in the cases of n = 3 and m = 4.     

On the Solution of Some Elliptic Difference Equations

Iterative methods for the solution of some nonlinear ellipticdifference systems, approximating the first boundary value problemare considered. If h > 0 is the network step in the spaceof variables x = (x₁, x₂,..., x_p) and 2m is the order of theoriginal boundary value problem, then the iterative methodsproposed give solution of accuracy with the expenditure ofO(|In | h^–(p+m–)) and O(|In | |In h| h^–p)arithmetic operations in the case of a general region and arectangular parallelepiped respectively. In the case p = 2 theestimate O(|In | h^{–[2+ (m/2)]}) is obtained if the regionis made up of rectangles with sides parallel to the co-ordinateaxes.     

一类非线性椭圆型方程奇摄动问题

A class of singularly perturbed problems for the nonlinear elliptic equations is considered. Under suitable conditions, using the theory of differential inequalities the asymptotic behavior of solution for the boundary value problems are studied, which reduced equations possess two intersecting solutions.     

Strong Solution of the Obstacle Problem for Fully Nonlinear Elliptic Equations

In this paper we obtain the existence of W^{2, ∞} solutions of the obstacle problems for fully nonlinear elliptic equations under more general structure conditions than those in [1] by using the mollifier approach, which is also extended in our discussion.     

Coupling of FEM and BEM for a nonlinear interface problem: The h–p version

This article presents some numerical examples for coupling the finite element method (FEM) and the boundary element method (BEM) as analyzed in [11]. This coupling procedure combines the advantages of boundary elements (problems in unbounded regions) and of finite elements (nonlinear problems with inhomogeneous data). In [28], experimental rates of convergence for the h version are presented, where the accuracy of the Galerkin approximation is achieved by refining the mesh. In this article we treat the h–p version, combining an increase of the degree of the piecewise polynomials with a certain mesh refinement. In our model examples, we obtain theoretically and numerically exponential convergence, which indicates a great efficiency in particular if singularities appear. © 1995 John Wiley & Sons, Inc.