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.     

Analysis of a kinematic dynamo model with FEM–BEM coupling

We consider a kinematic dynamo model in a bounded interior simply connected region Ω and in an insulating exterior region . In the so‐called direct problem, the magnetic field B and the electric field E are unknown and are driven by a given incompressible flow field w . After eliminating E , a vector and a scalar potential ansatz for B in the interior and exterior domains, respectively, are applied, leading to a coupled interface problem. We apply a finite element approach in the bounded interior domain Ω, whereas a symmetric boundary element approach in the unbounded exterior domain Ω^c is used. We present results on the well‐posedness of the continuous coupled variational formulation, prove the well‐posedness and stability of the semi‐discretized and fully discretized schemes, and provide quasi‐optimal error estimates for the fully discretized scheme. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

四阶椭圆型方程奇异摄动问题的渐近解

本文考虑了四阶椭圆型偏微分方程奇异摄动边值问题,建立了解及其导数的能量估计,并用Lyuternik-Vishik方法构造了形式渐近解.最后利用能量估计我们得到了渐近展开式余项的界.     

四阶椭圆型方程奇异摄动问题的数值解

本文对一类四阶椭圆型方程奇异摄动问题建立了指数型拟合差分格式,并且证明了这种格式在能量范数意义下关于小参数ε的一致收敛性.最后,我们给出了数值结果.     

Adaptive Frame Algorithms for Elliptic Operator Equations

We discuss the use of Gelfand frames for the adaptive numerical solution of linear elliptic operator equations. After a transformation into an equivalent infinite–dimensional system Lu = f in frame coordinates, the operator equation can be solved within a prescribed target accuracy by approximate Richardson or gradient iteration schemes. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Interface Problems for Dispersive Equations

The interface problem for the linear Schrödinger equations in one‐dimensional piecewise homogeneous domains is examined by providing an explicit solution in each domain. The location of the interfaces is known and the continuity of the wave function and a jump in their derivative at the interface are the only conditions imposed. The problem of two semi‐infinite domains and that of two finite‐sized domains are examined in detail. The problem and the method considered here extend that of an earlier paper by Deconinck et al. (2014) [1]. The dispersive nature of the problem presents additional difficulties that are addressed here.     

自共轭椭圆型方程的样条解法

本文从二元样条空间的理论出发,构造了一类新的差分格式,并利用它得到了一类自共轭椭圆型方程的样条解,并证明了这样的解的唯一性和收敛性问题.最后,给出了一个数值例子,说明了本方法是可行的.     

Composite Finite Elements for Elliptic Interface Problems

We present a Composite Finite Element Method for the approximation of linear elliptic boundary value problems of Dirichlet type with discontinuous coefficients. The challenge is the discontinuity of the coefficient (interface) which is not necessarily resolved by the underlying finite element mesh. The method is non-conforming in the sense that shape functions preserve continuity across the interface only in an approximative way. However, the construction allows to balance the non-conformity and the overall discretization error. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

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.     

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上的连续函数几乎处处相等的结论.     

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

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

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.