Computer Methods in Applied Mechanics and Engineering
The discrete artificial boundary condition on a polygonal artificial boundary for the exterior problem of Poisson equation by using the direct method of lines☆
Introduction
When computing the numerical solution of a boundary value problem of partial differential equations in an unbounded domain, one often introduces an artificial boundary to cut off the unbounded part of the domain and sets up an artificial boundary condition at the artificial boundary. Then the original problem is reduced to a boundary value problem defined in a bounded computational domain. In order to limit the computational cost, the bounded computational domain must not be too large. Then the artificial boundary condition must be a good approximation of the exact boundary condition on the artificial boundary. Thus the accuracy of the artificial boundary condition and the computational cost are closely related. Therefore designing an artificial boundary condition with high accuracy on a given artificial boundary has become an effective and important method for solving partial differential equations in an unbounded domain which arises in various fields of engineering, such as fluid flow around obstacles, coupling of structures with foundation, wave propagation and so on.
There are many authors who have worked on this subject for various problems by different techniques. For instance, Engquist and Majda [1] designed absorbing boundary conditions for the wave equation. Goldstein [2] presented the exact boundary condition and a sequence of its approximations at an artificial boundary for Helmholtz-type equation in waveguides. Feng [3] proposed the asymptotic radiation conditions for the reduced wave equation by using the asymptotic approximation of Hankel functions. Han and Wu [4], [5] obtained the exact boundary conditions and a series of their approximations at an artificial boundary for the Laplace equation and the linear elastic system. The exact boundary condition at an artificial boundary for partial differential equations in an infinite cylinder was proposed by Hagstrom and Keller [6], [7]. Shortly after, they used this technique to solve nonlinear problems. A family of artificial boundary conditions for unsteady Oseen equations in the velocity pressure formulation with small viscosity was developed by Halpern and Schatzman [8], which was then applied to unsteady Navier–Stokes equations. Nataf [9] designed an open boundary condition for a steady Oseen equation in streamfunction vorticity formulation, which is applied to viscous incompressible fluid flow around a body in a flat channel with slip boundary conditions on the wall. Hagstrom [10], [11] proposed asymptotic boundary conditions at an artificial boundary for the simulation of time-dependent fluid flows. Han et al. [12] designed discrete artificial boundary conditions for incompressible viscous flows in an infinite channel by using a fast iterative method. Han and Bao [13], [14] proposed discrete artificial boundary conditions for incompressible viscous flows in a channel by using the method of lines. Han et al. [15] developed artificial boundary conditions for the problem of infinite elastic foundation. Recently Ben-Porat and Givoli [16] considered an elliptic artificial boundary for the Laplace equation. One can find more references in Ref. [17].
Because of the restriction of the techniques they used, many authors mainly consider the regular artificial boundaries, such as circumferences, straight lines or a segment, as artificial boundaries in solving two-dimensional problems. As we know, it is very easy to implement discretizing a boundary value problem in a bounded domain with a polygonal boundary by using the finite element method [18]. Thus from the engineering point of view, it is natural to introduce a polygonal boundary as an artificial boundary for the problem in an unbounded domain. Thus how to design an artificial boundary condition with high accuracy at a given polygonal artificial boundary becomes an interesting open problem. In this paper, we propose a method to design a discrete artificial boundary condition at a given polygonal artificial boundary for the exterior problem of Poisson equation by using the direct method of lines. Then the problem is reduced to a boundary value problem defined in a bounded computational domain. Finite element approximation of the reduced problem is also considered. Furthermore numerical examples show that the discrete artificial boundary condition presented in this paper is very effective.
Section snippets
The discrete artificial boundary condition at a polygonal artificial boundary
Let Γi be a bounded, simple and closed curve in and be the unbounded domain with boundary Γi. We consider the following model problem, the exterior problem of Poisson equation:where g is a given function on Γi, f is a given function in and we assume that its support is compact. This problem is defined in the unbounded domain . In Ref. [4], this problem was considered. They introduced a circumference as an artificial boundary and designed a series of
The numerical solution of the problem (2.1)–(2.3)
On the bounded computational domain , we now consider the numerical solution of the problem (2.1)–(2.3). As we know that the restriction of u, the solution of the problem (2.1)–(2.3), satisfies the boundary value problem (2.9)–(2.11). Let denote the usual Sobolov space on [19] and assume thatThen the boundary value problem (2.9)–(2.11) is equivalent to the following variational problem:with
Numerical implementation and examples
Let denote the exterior domain of the unit square, namelyWe consider the numerical solution of the original problem (2.1)–(2.3) with given f and g. We take the artificial boundary . Hence and . Since the solution of each example, u(x1,x2), is symmetric about x2 axes and antisymmetric about x1 axes, respectively, the domain of computation is taken to be the part of lying in the
References (20)
An open boundary condition for the computation of the steady incompressible Navier-Stokes equations
J. Comput. Phys.
(1989)- et al.
A discrete artificial boundary condition for steady incompressible viscous flows in a no-slip channel using a fast iterative method
J. Comput. Phys.
(1994) - et al.
Numerical simulation for the problem of infinite elastic foundation
Comput. Methods Appl. Mech. Engrg.
(1997) - et al.
Absorbing boundary conditions for the numerical simulation of waves
Math. Comp.
(1977) A finite element method for solving Helmholtz type equations in waveguides and other unbounded domains
Math. Comp.
(1982)Asymptotic radiation conditions for reduced wave equations
J. Comput. Math.
(1984)- et al.
Approximation of infinite boundary condition and its application to finite element method
J. Comput. Math.
(1985) - et al.
The approximation of exact boundary condition at an artificial boundary for linear elastic equation and its application
Math. Comp.
(1992) - et al.
Exact boundary conditions at artificial boundary for partial differential equations in cylinders
SIAM J. Math. Anal.
(1986) - et al.
Asymptotic boundary conditions and numerical methods for nonlinear elliptic problems on unbounded domains
Math. Comp.
(1987)
Cited by (14)
The direct method of lines for elliptic problems in star-shaped domains
2018, Journal of Computational and Applied MathematicsCitation Excerpt :The corresponding numerical example is Example 4.2. What is more, it can also be applied to the boundary value problem outside a bounded domain [37,38]. The corresponding numerical example is Example 4.3.
Solving two dimensional second order elliptic equations in exterior domains using the inverted finite elements method
2016, Computers and Mathematics with ApplicationsCitation Excerpt :Ciskowski and Brebbia [2], Beer et al. [3] and Colton and Kress [4]), Artificial boundary condition methods (see, e.g., [5–9]), Absorbing Boundary Conditions (ABC) methods (see, e.g., Enquist and Majda [10,11], or Bayliss and Turkel [12]) or Perfectly Matched Layers (PML) methods (see Bérenger [13]). Among non truncature methods we can mention the Infinite Element Method (IEM) (see Bettess [14], Bettess and Zienkiewicz [15], Burnett [16], Gerdes and Demkowicz [17], Shirron and Babuvska [18] and Toselli [19]) or spectral-like methods (see Shen and Wang [20] or Boulmezaoud et al. [21]).
Numerical simulations of fracture problems by coupling the FEM and the direct method of lines
2001, Computer Methods in Applied Mechanics and EngineeringEfficient numerical method for Stokes flows in unbounded domains with informative boundary condition using axial Green function method
2023, International Journal for Numerical Methods in Fluids
- ☆
This work was supported partly by the Climbing Program of National Key Project of Foundation, Doctoral Program foundation of Institution of Higher Education and the National Natural Science Foundation of China. Computation was supported by the State Key Lab. of Scientific and Engineering Computing in China.