The boundary integral method for the steady rotating Navier–Stokes equations in exterior domain (I): the existence of solution

In this paper, we apply the boundary integral method to the steady rotating Navier–Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and a infinite domain, we obtain a coupled problem by the steady rotating Navier–Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence of solution in a convex set.     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

The coupling of boundary integral and finite element methods for the bidimensional exterior steady stokes problem

In this paper, we represent a new numerical method for solving the steady-state Stokes equations in an unbounded plane domain. The technique consists in coupling the boundary integral and the finite element methods. An artificial smooth boundary is introduced separating an interior inhomogeneous region from an exterior one. The solution in the exterior domain is represented by an integral equation over the artificial boundary. This integral equation is incorporated into a velocitypressure formulation for the interior region, and a finite element method is used to approximate the resulting variational problem. This is studied by means of an abstract framework, well adapted to the model problem, in which convergence results and optimal error estimates are derived. Computer results will be discussed in a forthcoming paper.     

A new formulation of regularized meshless method applied to interior and exterior anisotropic potential problems

This paper proposes a new formulation of regularized meshless method (RMM), which differs from the traditional RMM in that the traditional formulation generates the diagonal elements of influence matrix via null-field integral equations, while our new one directly employs the boundary integral equations at the domain point to evaluate the diagonal elements. We test the present RMM formulation to two-dimensional anisotropic potential problems in finite and infinite domains in comparison with the traditional RMM. Numerical results show that the present RMM sharply outperforms the traditional RMM in the solution of interior problems, while the latter is clearly superior for exterior problems. A rigorous theoretical analysis of circular domain case also corroborates such numerical experiment observations and is provided in the appendix of this paper.     

Factorization of meromorphic matrix functions

The paper considers a class of matrix-functions defined on some contour in the complex plane that have meromorphic continuations in the interior or exterior domain of that contour. These matrix-functions generally do not admit Wiener-Hopf standard factorization. The paper studies the problem of index factorization, which is a version of Wiener-Hopf factorization. Some criterions for index factorization and exact formulas for particular indices are found along with a constructive method which applies the factorization to solution of an explicit, finite system of linear algebraic equations.     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...     

A New Non-Overlapping Domain Decomposition Method for a 3D Laplace Exterior Problem

We propose a method for solving three-dimensional boundary value problems for Laplace’s equation in an unbounded domain. It is based on non-overlapping decomposition of the exterior domain into two subdomains so that the initial problem is reduced to two subproblems, namely, exterior and interior boundary value problems on a sphere. To solve the exterior boundary value problem, we propose a singularity isolation method. To match the solutions on the interface between the subdomains (the sphere), we introduce a special operator equation approximated by a system of linear algebraic equations. This system is solved by iterative methods in Krylov subspaces. The performance of the method is illustrated by solving model problems.     

Elementare Fälle des Dirichletschen Problems für elliptische Gebiete der Ebene

Summary The Dirichlet problem is solved for an elliptic domain if the boundary values are given by a polynomial defined in the interior of the ellipse. The solution is a finite expansion in harmonic polynomials.     

A fast and efficient two-grid method for solving d-dimensional poisson equations

The aim of this paper is to introduce a fast and efficient new two-grid method to solve the d-dimensional (d=1,2,3) Poisson elliptic equations. The finite difference equations at all interior grid points form a large sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The finite difference equations are based on applying a finite difference scheme of two- and four-orders (compact finite difference method) for discretizing the spatial derivative. The obtained linear systems of Poisson elliptic equations have been solved by a new two-grid (NTG) method and we also note that the NTG method which is used for solving the large sparse linear systems is faster and more effective than that of the standard two-grid method. We utilize the local Fourier analysis to show that the spectral radius of the new two-grid method for 1D and 2D models is less than that of the standard two-grid method. As well as, we expand the corresponding algorithm to the new multi-grid method. The numerical examples show the efficiency of the new algorithms for solving the d-dimensional Poisson equations.     

The difference method in a diffraction problem: The technique of interior boundary condition

The diffraction problem for a plane wave on a half-plane covered by thin layer with an interface is solved by the difference method. The system of difference equations is derived from the variational principle. A boundary solution at infinity must be imposed; this is a radiation condition, which is used in the form of the limit absorption principle. The arising infinite system of difference equations is reduced to a finite part of the boundary (the interface) by using the technique of so-called interior boundary conditions in the sense of Ryaben’kii. The real conditions are found by the Fourier method with respect to one spatial variable in the form of Fourier or Laurent series in the corresponding variable, which converge either inside, outside, or on the unit circle. Above the upper boundary of the layer, all unknowns are eliminated by using the so-called grid Green function, that is, the resolving function for the half-plane satisfying the radiation condition at infinity. For the unknowns on the upper boundary of the layer, an equation in terms of a function of a complex variable of Wiener-Hopf type is obtained, which is solved by factorization. Factorization is performed numerically: the logarithm of the function is expanded in a bi-infinite series, which is replaced by a discrete Fourier series. The closing system in a neighborhood of the interface has order proportional to the number of points on the interface. Solving this system yields all of the required characteristics of the solution.     

Interface problem in holonomic elastoplasticity

The three-dimensional interface problem with the homogeneous Lamé system in an unbounded exterior domain and holonomic material behaviour in a bounded interior Lipschitz domain is considered. Existence and uniqueness of solutions of the interface problem are obtained rewriting the exterior problem in terms of boundary integral operators following the symmetric coupling procedure. The numerical approximation of the solutions consists in coupling of the boundary element method (BEM) and the finite element method (FEM). A Céa-like error estimate is presented for the discrete solutions of the numerical procedure proving its convergence.     

Finite Element Methods for the Stokes System in Three-Dimensional Exterior Domains

The article treats the question of how to numerically solve the Dirichlet problem for the Stokes system in the exterior of a three-dimensional bounded Lipschitz domain. In a first step, the solution of this problem is approximated by functions solving the Stokes system in a truncated domain and satisfying a suitable artificial boundary condition on the outer boundary of this truncated domain. In a second step, this new problem is approximately solved in finite element spaces related to a graded mesh as introduced by Goldstein [Math. Comp., 36, 387–404 (1981)]. The difference between this finite element approximation and the exact solution of the exterior Stokes problem is estimated in the norm of suitable unweighted L²-Sobolev spaces. These estimates are analogous to corresponding results which are known for the Poisson equation. © 1997 by B.G. Teubner Stuttgart-John Wiley & Sons, Ltd.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

Adapted kussmaul formulation of acoustic scattering problems: an overview

Problems of exterior acoustic scattering may be conveniently formulated by means of boundary integral equations. The problem seeks to find a wave function which gives velocity potential profile, pressure density profile, etc. of the acoustic wave at points in space. At the background of the formulations are two theories viz. (Helmholtz) Potential theory and the Green's representation formula. Potential theory gives rise to the so-called indirect formulation and the Green's representation formula to the direct formulations. Classical boundary integral formulations fail at the eigenfrequencies of the interior domain. That is, if a solution is sought of the exterior problem by first solving a homogeneous boundary integral equation, one is inevitably led to the conclusion that these homogeneous boundary equations have nontrivial solutions at certain wave-numbers which are the eigenvalues of the corresponding interior problem. At lower wave-numbers, these eigenfrequencies are thinly distributed but the higher the wave-number, the denser it becomes. This is a well-known drawback for both time-harmonic acoustics and elastodynamics. This is not a physical difficulty but arises entirely as a result of a deficiency in the integral equation is representation. Why then use It? The use has many advantages notably in that the meshing region is reduced from the infinite domain exterior to the body to its finite surface. This created the need for some robust formulations. A proof of the Kussmaul [1] formulation is presented. The formulation has a hypersingular kernel in the integral operator, which creates a havoc in computation (e.g., ill conditioning). The hyper-singularity can be avoided [2], as a result a new formulation is proposed. This paper presents a broad overview of the Adapted Kussmaul Formulation (AKF).     

Iterative solution of a coupled mixed and standard Galerkin discretization method for elliptic problems

In this paper, we consider approximation of a second‐order elliptic problem defined on a domain in two‐dimensional Euclidean space. Partitioning the domain into two subdomains, we consider a technique proposed by Wieners and Wohlmuth [9] for coupling mixed finite element approximation on one subdomain with a standard finite element approximation on the other. In this paper, we study the iterative solution of the resulting linear system of equations. This system is symmetric and indefinite (of saddle‐point type). The stability estimates for the discretization imply that the algebraic system can be preconditioned by a block diagonal operator involving a preconditioner for H (div) (on the mixed side) and one for the discrete Laplacian (on the finite element side). Alternatively, we provide iterative techniques based on domain decomposition. Utilizing subdomain solvers, the composite problem is reduced to a problem defined only on the interface between the two subdomains. We prove that the interface problem is symmetric, positive definite and well conditioned and hence can be effectively solved by a conjugate gradient iteration. Copyright © 2001 John Wiley & Sons, Ltd.     

Semi‐analytical solution of MHD flow through boundary integrals on the pipe wall

A mathematical model is given for the magnetohydrodynamic (MHD) pipe flow as an inner Dirichlet problem in a 2D circular cross section of the pipe, coupled with an outer Dirichlet or Neumann magnetic problem. Inner Dirichlet problem is given as the coupled convection‐diffusion equations for the velocity and the induced current of the fluid coupling also to the outer problem, which is defined with the Laplace equation for the induced magnetic field of the exterior region with either Dirichlet or Neumann boundary condition. Unique solution of inner Dirichlet problem is obtained theoretically reducing it into two boundary integral equations defined on the boundary by using the corresponding fundamental solutions. Exterior solution is also given theoretically on the pipe wall with Poisson integral, and it is unique with Dirichlet boundary condition but exists with an additive constant obtained through coupled boundary and solvability conditions in Neumann wall condition. The collocation method is used to discretize these boundary integrals on the pipe wall. Thus, the proposed procedure is an improved theoretical analysis for combining the solution methods for the interior and exterior regions, which are consolidated numerically showing the flow behavior. The solution is simulated for several values of problem parameters, and the well‐known MHD characteristics are observed inside the pipe for increasing values of Hartmann number maintaining the continuity of induced currents on the pipe wall.     

A System of Plane Elasticity Canonical Integral Equations and Its Application

In this paper, we obtain a new system of canonical integral equations for the plane elasticity problem over an exterior circular domain, and give its numerical solution. Coupling with the classical finite element method, it can be used for solving general plane elasticity exterior boundary value problems. This system of highly singular equations is also an exact boundary condition on the artificial boundary. It can be approximated by a series of nonsingular integral boundary conditions.     

Haar wavelets method for solving Poisson equations with jump conditions in irregular domain

In this paper, Haar wavelets method is used to solve Poisson equations in the presence of interfaces where the solution itself may be discontinuous. The interfaces have jump conditions which need to be enforced. It is critical for the approximation of the boundaries of the irregular domain. An irregular domain can be treated by embedding the domain into a rectangular domain and Poisson equation is solved by using Haar wavelets method on the rectangle. Firstly, we demonstrate this method in the case of 1-D region, then we consider the solution of the Poisson equations in the case of 2-D region. The efficiency of the method is demonstrated by some numerical examples.     

Vibro-acoustic Investigations Using Finite and Infinite Elements

The solution of large systems of equations often involves the use of iterative solution procedures, mostly by means of Krylov subspace methods. The performance of these methods, however, strongly depends on the spectrum of the resulting system matrices, which is affected by the polynomial approximation function space. The current work shows that finite elements based on Bernstein polynomials yield particularly good performance in combination with commonly employed Krylov solvers. Adapting the basis of Astley-Leis infinite elements, the Bernstein polynomials may also be used for exterior acoustic simulations. The efficiency of such elements for interior acoustic problems as well as for sound radiation analyses is assessed. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Unified symmetric finite element and boundary integral variational coupling methods for linear fluid–structure interaction

This article concerns the development of energy-based variational formulations and their corresponding finite element–boundary element Rayleigh–Ritz approximations for solving the time-harmonic vibration and scattering problem of an inhomogeneous penetrable fluid or solid object immersed in a compressible, inviscid, homogeneous fluid. The resulting coupled finite element and boundary integral methods (FEM-BEM) have the following attractive features: (1) Separate direct and complementary variational principles lead naturally to several alternative structure variable and fluid variable methodologies. (2) The solution in the exterior region is represented by a combined single- and double-layer potential which ensures the validity of the methods for all wave numbers; even though this representation introduces hypersingular integrals, for actual computations the hypersingular operator may be rewritten in terms of single-layer potentials, which can be integrated by standard techniques. (3) Since the discretized equations for the interior region and for the boundary are derived from the first variation of bilinear functionals the resulting algebraic systems of equations are always symmetric. In addition, the transition conditions across the interface are natural. This allows one to approximate the solutions within the interior and exterior regions independently, without imposing any boundary constraints.