A robust boundary element method for nearly incompressible linear elasticity

Summary. For a Dirichlet boundary value problem in linear elasticity we consider a boundary element method which is robust for nearly incompressible materials. Based on the spectral properties of the single layer potential for the Stokes problem we introduce an orthogonal splitting of the trial space. The resulting variational problem is then well conditioned and can be discretized by using standard boundary element methods. Mathematics Subject Classification (1991):65N38     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

A new integral equation approach to the Neumann problem in acoustic scattering

We suggest a new approach of reduction of the Neumann problem in acoustic scattering to a uniquely solvable Fredholm integral equation of the second kind with weakly singular kernel. To derive this equation we placed an additional boundary with an appropriate boundary condition inside the scatterer. The solution of the problem is obtained in the form of a single layer potential on the whole boundary. The density in the potential satisfies a uniquely solvable Fredholm integral equation of the second kind and can be computed by standard codes. Copyright © 2001 John Wiley & Sons, Ltd.     

Algebraic Multigrid Preconditioners for Sparse Approximations of Boundary Element Matrices

This paper presents new algebraic multigrid preconditioners for sparse representations of boundary element matrices which arise from the so‐called adaptive‐cross‐approximation to dense boundary element matrices resulting from the standard collocation, or Galerkin boundary element discretization of the single layer potential operator.     

Solutions for the problem of boundary layer formation in a pseudo-plastic fluid: The case of gradual acceleration

We consider the development of the nonstationary boundary layer about a body that gradually starts to move in a resting fluid. Under certain conditions, we construct the solutions for the problem of formation of boundary layer in a pseudo-plastic fluid. The method used here is mainly based on a transformation which reduces the boundary layer system to a boundary value problem for a single quasilinear parabolic equation.     

Stabilized boundary element methods for exterior Helmholtz problems

In this paper we describe and analyze some modified boundary element methods to solve the exterior Dirichlet boundary value problem for the Helmholtz equation. As in classical combined field integral equations also the proposed approach avoids spurious modes. Moreover, the stability of related modified boundary element methods can be shown even in the case of Lipschitz boundaries. The proposed regularization is done based on boundary integral operators which are already included in standard boundary element formulations. Numerical examples are given to compare the proposed approach with other already existing regularized formulations.     

Boundary Integral Solution of the Time-Fractional Diffusion Equation

Here we consider initial boundary value problem for the time–fractional diffusion equation by using the single layer potential representation for the solution. We derive the equivalent boundary integral equation. We will show that the single layer potential admits the usual jump relations and discuss the mapping properties of the single layer operator in the anisotropic Sobolev spaces. Our main theorem is that the single layer operator is coercive in an anisotropic Sobolev space. Based on the coercivity and continuity of the single layer operator we finally show the bijectivity of the operator in a certain range of anisotropic Sobolev spaces.      

Weighted extended b-splines for one-dimensional electromagnetic problems

This paper considers the weighted extended b-splines as basis function for finite element method in electromagnetics and compares with the standard finite element method applied to the two-point boundary value problems with different boundary conditions. This new approach, which provides more accurate results than standard finite element method, is presented to compare other numerical techniques and applied to one-dimensional electromagnetic problems. Computed results are compared with other numerical results in literature.     

Analysis of a FEM/BEM coupling method for transonic flow computations

A sensitive issue in numerical calculations for exterior flow problems, e.g.around airfoils, is the treatment of the far field boundary conditions on a computational domain which is bounded. In this paper we investigate this problem for two-dimensional transonic potential flows with subsonic far field flow around airfoil profiles. We take the artificial far field boundary in the subsonic flow region. In the far field we approximate the subsonic potential flow by the Prandtl-Glauert linearization. The latter leads via the Green representation theorem to a boundary integral equation on the far field boundary. This defines a nonlocal boundary condition for the interior ring domain. Our approach leads naturally to a coupled finite element/boundary element method for numerical calculations. It is compared with local boundary conditions. The error analysis for the method is given and we prove convergence provided the solution to the analytic transonic flow problem around the profile exists.

     

Non-convex self-dual Lagrangians and variational principles for certain PDEʼs

We study the concept and the calculus of non-convex self-dual (Nc-SD) Lagrangians and their derived vector fields which are associated to many partial differential equations and evolution systems. They yield new variational resolutions for large class of partial differential equations with variety of linear and non-linear boundary conditions including many of the standard ones. This approach seems to offer several useful advantages: It associates to a boundary value problem several potential functions which can often be used with relative ease compared to other methods such as the use of Euler–Lagrange functions. These potential functions are quite flexible, and can be adapted to easily deal with both non-linear and homogeneous boundary value problems.     

On the boundary transmission problem of thermoelastostatics in a plane domain with interface corners

This paper deals with bimetal problems of thermoelastostatics. By means of an explicit particular solution a reduction to problems of elastostatics is given. An indirect boundary integral method is applied for solving the traction boundary value problem. The solution is represented by a potential of single layer type having Green's contact tensor as the kernel. Thus, from the first the transmission conditions are satisfied. The Fredholm property of the boundary integral operator as well as the asymptotics of the potential density at an interface corner depend on the symbol of a Mellin convolution operator. The singular functions at corners can be obtained by calculating the potential for terms in the asymptotic expansion of the density.     

Neumann and Robin problems in a cracked domain with jump conditions on cracks

The boundary value problem for the Laplace equation is studied on a domain with smooth compact boundary and with smooth internal cracks. The Neumann or the Robin condition is given on the boundary of the domain. The jump of the function and the jump of its normal derivative is prescribed on the cracks. The solution is looked for in the form of the sum of a single layer potential and a double layer potential. The solvability of the corresponding integral equation is determined and the explicit solution of this equation is given in the form of the Neumann series. Estimates for the absolute value of the solution of the boundary value problem and for the absolute value of the gradient of the solution are presented.     

A boundary parametric approximation to the linearized scalar potential magnetostatic field problem

We consider the linearized scalar potential formulation of the magnetostatic field problem in this paper. Our approach involves a reformulation of the continuous problem as a parametric boundary problem. By the introduction of a spherical interface and the use of spherical harmonics, the infinite boundary conditions can also be satisfied in the parametric framework. That is the field in the exterior of a sphere is expanded in a ‘harmonic series’ of eigenfunctions for the exterior harmonic problem. The approach is essentially a finite element method coupled with a spectral method via a boundary parametric procedure. The reformulated problem is discretized by finite element techniques which leads to a discrete parametric problem which can be solved by well conditioned iteration involving only the solution of decoupled Neumann type elliptic finite element systems and L² projection onto subspaces of spherical harmonics. Error and stability estimates given show exponential convergence in the degree of the spherical harmonics and optimal order convergence with respect to the finite element approximation for the resulting fields in L².     

定常Stokes问题的边界积分方程法

1.前言 定常Stokes问题本身虽然只反映在小雷诺数情况下不可压缩粘性流体的定常流动,然而却为处理完整的Navier-Stokes方程奠定了基础. Stokes问题一般有两种公式化途径,一是通过流函数,二是利用速度-压力公式.两种公式化途径的区域类型数值方法,如有限差分法及有限单元法,已有不少工作,见[3]和[11].近年来,对这两种公式化途径的边界类型数值方法的研究,也获得一些结果.     

Non-convex self-dual Lagrangians: New variational principles of symmetric boundary value problems

We study the concept and the calculus of Non-convex self-dual (Nc-SD) Lagrangians and their derived vector fields which are associated to many partial differential equations and evolution systems. They indeed provide new representations and formulations for the superposition of convex functions and symmetric operators. They yield new variational resolutions for large class of Hamiltonian partial differential equations with variety of linear and nonlinear boundary conditions including many of the standard ones. This approach seems to offer several useful advantages: It associates to a boundary value problem several potential functions which can often be used with relative ease compared to other methods such as the use of Euler-Lagrange functions. These potential functions are quite flexible, and can be adapted to easily deal with both nonlinear and homogeneous boundary value problems. Additionally, in most cases the solutions generated using this new method have greater regularity than the solutions obtained using the standard Euler-Lagrange function. Perhaps most remarkable, however, are the permanence properties of Nc-SD Lagrangians; their calculus is relatively manageable, and their applications are quite broad.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

Integral equation methods for scattering by infinite rough surfaces

In this paper, we consider the Dirichlet and impedance boundary value problems for the Helmholtz equation in a non‐locally perturbed half‐plane. These boundary value problems arise in a study of time‐harmonic acoustic scattering of an incident field by a sound‐soft, infinite rough surface where the total field vanishes (the Dirichlet problem) or by an infinite, impedance rough surface where the total field satisfies a homogeneous impedance condition (the impedance problem). We propose a new boundary integral equation formulation for the Dirichlet problem, utilizing a combined double‐ and single‐layer potential and a Dirichlet half‐plane Green's function. For the impedance problem we propose two boundary integral equation formulations, both using a half‐plane impedance Green's function, the first derived from Green's representation theorem, and the second arising from seeking the solution as a single‐layer potential. We show that all the integral equations proposed are uniquely solvable in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including an incident plane wave, the impedance boundary value problem for the scattered field has a unique solution under certain constraints on the boundary impedance. Copyright © 2003 John Wiley & Sons, Ltd.     

The asymmetrical mixed temperature problem for a transversely isotropic elastic layer

The problem of the uniform heating of a two-layer plate is solved. The transversely isotropic elastic layer (soft plate) investigated is in ideal contact with an absolutely rigid layer, deformable only by thermal expansion. The generalized plane temperature problem reduces to determining the stress-strain state of the soft anisotropic layer investigated using the equations of the mixed problem of elasticity theory. At the ends of the boundary layer of the soft plate (a thin contact layer), no conditions are imposed. On the remaining part of the ends of the soft plate, the boundary conditions correspond to a free boundary. The problem has a bounded smooth solution. Unlike the approach described earlier [1], it is proposed to seek an accurate solution in the form of ordinary Fourier series with respect to a single longitudinal coordinate. Solutions in polynomials are also used. It is shown that the existence of these solutions in polynomials enables the convergence of the Fourier series to be improved considerably.     

L∞-boundedness of the finite element galerkin operator for parabolic problems

In this paper the heat equation with Dirichlet boundary conditions in N ≤ 3 space dimensions - serving as model problem of second order parabolic initial boundary value problems - is considered. We prove: The standard finite element method is uniformly bounded in L_∞ with respect to space and time if the underlying finite elements are at least cubics.     

On the global existence and uniqueness of solutions to the nonstationary boundary layer system

In this paper, we study the problem of boundary layer for nonstationary flows of viscous incompressible fluids. There are some open problems in the field of boundary layer. The method used here is mainly based on a transformation which reduces the boundary layer system to an initial-boundary value problem for a single quasilinear parabolic equation. We prove the existence of weak solutions to the modified nonstationary boundary layer system. Moreover, the stability and uniqueness of weak solutions are discussed.