An adaptive finite element DtN method for the three-dimensional acoustic scattering problemThanks: The work of GB is supported in part by an NSFC Innovative Group Fund (No.11621101). The research of PL is supported in part by the NSF grant DMS-1912704.

Gang Bao Address: School of Mathematical Science, Zhejiang University, Hangzhou 310027, China. Email address: baog@zju.edu.cn , Mingming Zhang Address: School of Mathematical Science, Zhejiang University, Hangzhou 310027, China. Email address: mmzaip@zju.edu.cn , Bin Hu Address: School of Mathematical Science, Zhejiang University, Hangzhou 310027, China. Email address: binh@zju.edu.cn and Peijun Li Address: Department of Mathematics, Purdue University, West Lafayette, Indiana 47907, USA. Email address: lipeijun@math.purdue.edu

Abstract.

This paper is concerned with a numerical solution of the acoustic scattering by a bounded impenetrable obstacle in three dimensions. The obstacle scattering problem is formulated as a boundary value problem in a bounded domain by using a Dirichlet-to-Neumann (DtN) operator. An a posteriori error estimate is derived for the finite element method with the truncated DtN operator. The a posteriori error estimate consists of the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to decay exponentially with respect to the truncation parameter. Based on the a posteriori error estimate, an adaptive finite element method is developed for the obstacle scattering problem. The truncation parameter is determined by the truncation error of the DtN operator and the mesh elements for local refinement are marked through the finite element approximation error. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.

Key words and phrases:

acoustic scattering problem, adaptive finite element method, transparent boundary condition, a posteriori error estimates

2010 Mathematics Subject Classification

65M30, 78A45, 35Q60

1. Introduction

Wave scattering by bounded impenetrable media is usually referred to as the obstacle scattering problem. It has played an important role in many scientific areas such as radar and sonar, non-destructive testing, medical imaging, and geophysical exploration [12]. Due to the significant applications, the obstacle scattering problem has been extensively studied in the past several decades. Consequently, a variety of methods have been developed to solve the scattering problem mathematically and numerically such as the method of boundary integral equations [11, 27] and the finite element method [23, 26]. This paper concerns a numerical solution of the acoustic wave scattering by an obstacle in three dimensions.

As an exterior boundary value problem, the obstacle scattering problem is formulated in an open domain, which needs to be truncated into a bounded computational domain when applying numerical methods such as the finite element method. It is indispensable to impose a boundary condition on the boundary of the truncated domain. The ideal boundary condition is to completely avoid artificial wave reflection by mimicking the wave propagation as if the boundary did not exist [6]. Such a boundary condition is called an absorbing boundary condition [13], a nonreflecting boundary condition [16], or a transparent boundary condition (TBC) [17]. It still remains as an active research topic in computational wave propagation [18], especially for time-domain scattering problems [2]. Since Berenger proposed the perfectly matched layer (PML) technique for the time-domain Maxwell equations [7], the PML method has been extensively studied for various wave propagation problems [5, 10, 29]. As an effective approach for the domain truncation, the basic idea of the PML technique is to surround the domain of interest by a layer of finite thickness with specially designed artificial medium that would attenuate all the waves coming from inside of the domain. Combined with the PML technique, the a posteriori error estimate based adaptive finite element methods were developed for the diffraction grating problems [8, 4] and the obstacle scattering problems [9]. It was shown that the estimates consist of the finite element discretization error and the PML truncation error which has an exponential rate of convergence with respect to the PML parameters.

Recently, an alternative adaptive finite element method was developed for solving the two-dimensional acoustic obstacle scattering problem [21], where the PML was replaced by the TBC to truncate the open domain. Since the TBC is exact, it can be imposed on the boundary which could be put as close as possible to the obstacle. Hence it does not require an extra absorbing layer of artificial medium to enclose the domain of interest. Based on a nonlocal Dirichlet-to-Neumann (DtN) operator, the TBC is given as an infinite Fourier series. Practically, the series needs to be truncated into a sum of finitely many, say $N$ , terms, where $N$ is an appropriately chosen positive integer. In [21], an a posteriori error estimate was derived for the finite element discretization but it did not include the truncation error of the DtN operator. The complete a posteriori error estimate was obtained in [22]. The new estimate takes into account both the finite element discretization error and the DtN operator truncation error. It was shown that the truncation error decays exponentially with respect to the truncation parameter $N$ . The adaptive finite element DtN method has also been applied to solve the diffraction grating problems [30] as well as the elastic wave equation in periodic structures [25]. The numerical results show that the adaptive finite element DtN method is competitive with the adaptive finite element PML method.

In this work, we extend the analysis in [22] to the three-dimensional obstacle scattering problem. It is worthy to mention that the extension is nontrivial since more complex spherical Hankel functions need to be considered and the computation is more challenging in three dimensions. Specifically, we consider the acoustic wave scattering by a sound hard obstacle. Based on a TBC, the exterior problem is formulated equivalently into a boundary value problem in a bounded domain for the three-dimensional Helmholtz equation. Using a duality argument, we derive the a posteriori error estimate which includes the finite element discretization error and the DtN operator truncation error. Moreover, we show that the truncation error has an exponential rate of convergence with respect to the truncation parameter $N$ . The a posteriori estimate is used to design the adaptive finite element algorithm to choose elements for refinements and to determine the truncation parameter $N$ . In addition, we present a technique to deal with adaptive mesh refinements of the surface. Numerical experiments are included to demonstrate the effectiveness of the proposed method.

This paper is organized as follows. In Section 2, we introduce the model problem of the acoustic wave scattering by an obstacle in three dimensions. The variational formulation is given for the boundary value problem by using the DtN operator. In Section 3, we present the finite element approximation with the truncated DtN operator. Section 4 is devoted to the a posteriori error analysis by using a duality argument. In Section 5, we discuss the numerical implementation and the adaptive finite element DtN method, and present two numerical examples to demonstrate the effectiveness of the proposed method. The paper is concluded with some general remarks and directions for future work in Section 6.

2. Problem formulation

Consider a bounded sound-hard obstacle $D$ with Lipschitz continuous boundary $\partial D$ in $\mathbb{R}^{3}$ . Denote by $B_{r}=\{x\in\mathbb{R}^{3}:|x|<r\}$ the ball which is centered at the origin and has a radius $r$ . Let $R$ and $R^{\prime}$ be two positive constants such that $R>R^{\prime}>0$ and $\overline{D}\subset B_{R^{\prime}}\subset B_{R}$ . Denote $\Omega=B_{R}\backslash\overline{D}$ . The obstacle scattering problem for acoustic waves can be modeled by the following exterior boundary value problem:

\begin{cases}\Delta u+\kappa^{2}u=0\quad&{\rm in}~\mathbb{R}^{3}\setminus\overline{D},\\ \partial_{\nu}u=-g\quad&{\rm on}~\partial D,\\ \lim\limits_{r\rightarrow\infty}r(\partial_{r}u-{\rm i}\kappa u)=0,&r=|x|,\end{cases}

(2.1)

where $\kappa>0$ is the wavenumber and $\nu$ is the unit outward normal vector to $\partial D$ . Although the results are given for the sound-hard boundary condition in this paper, the method can be applied to other types of boundary conditions, such as the sound-soft and impedance boundary conditions.

Let $\hat{x}_{1}=\sin\theta\cos\varphi$ , $\hat{x}_{2}=\sin\theta\sin\varphi$ , $\hat{x}_{3}=\cos\theta$ , $\theta\in[0,\pi]$ and $\varphi\in[0,2\pi]$ . Introduce the spherical harmonic functions

Y_{n}^{m}(\hat{x})=Y_{n}^{m}(\theta,\varphi)=\sqrt{\frac{(2n+1)(n-|m|)!}{4\pi(n+|m|)!}}P_{n}^{|m|}(\cos\theta)e^{{\rm i}m\varphi},\quad m=-n,\dots,n,\,n=0,1,\dots,

where

P_{n}^{m}(t)=(1-t^{2})^{\frac{m}{2}}\frac{{\rm d}^{m}}{{\rm d}t^{m}}P_{n}(t),\quad-1\leq t\leq 1,

are called the associated Legendre functions and $P_{n}$ are the Legendre polynomials. It is known that the spherical harmonic functions $\{Y_{n}^{m}:m=-n,\dots,n,\,n=0,1,\dots\}$ form an orthonormal system in $L^{2}(\mathbb{S}^{2})$ , where $\mathbb{S}^{2}=\{x\in\mathbb{R}^{3}:|x|=1\}$ is the unit sphere in $\mathbb{R}^{3}$ . For any function $u\in L^{2}(\partial B_{R})$ , it admits the Fourier series expansion

u(x):=u(R,\hat{x})=\sum_{n=0}^{\infty}\sum_{m=-n}^{m=n}\hat{u}_{n}^{m}(R)Y_{n}^{m}(\hat{x}),\quad\hat{u}_{n}^{m}=\int_{\mathbb{S}^{2}}u(R,\hat{x})\bar{Y}_{n}^{m}(\hat{x}){\rm d}\hat{x}.

Using the Fourier coefficients, we may define an equivalent $L^{2}(\partial B_{R})$ norm of $u$ as

\|u\|_{L^{2}(\partial B_{R})}=\left(\sum_{n=0}^{\infty}\sum_{m=-n}^{n}|\hat{u}_{n}^{m}|^{2}\right)^{\frac{1}{2}}.

The trace space $H^{s}(\partial B_{R})$ is defined by

H^{s}(\partial B_{R})=\left\{u\in L^{2}(\partial B_{R})~:~\|u\|_{H^{s}(\partial B_{R})}<\infty\right\},

where the norm may be characterized by

\|u\|_{H^{s}(\partial B_{R})}^{2}=\sum_{n=0}^{\infty}\sum_{m=-n}^{m=n}\left(1+n(n+1)\right)^{s}|\hat{u}_{n}^{m}|^{2}.

(2.2)

Clearly, the dual space of $H^{-s}(\partial B_{R})$ is $H^{s}(\partial B_{R})$ with respect to the scalar product in $L^{2}(\partial B_{R})$ defined by

\langle u,v\rangle_{\partial B_{R}}=\int_{\partial B_{R}}u\bar{v}{\rm d}s.

In the exterior domain $\mathbb{R}^{3}\backslash\overline{B}_{R}$ , the solution of the Helmholtz equation in (2.1) can be written as

u(r,\hat{x})=\sum_{n=0}^{\infty}\sum_{m=-n}^{n}\hat{u}_{n}^{m}\frac{h_{n}^{(1)}(\kappa r)}{h_{n}^{(1)}(\kappa R)}Y_{n}^{m}(\hat{x}),\quad r>R,

(2.3)

where $h_{n}^{(1)}$ is the spherical Hankel function of the first kind with order $n$ and is defined as (cf. [19])

h_{n}^{(1)}(z)=\sqrt{\frac{\pi}{2z}}H_{n+\frac{1}{2}}^{(1)}(z).

Here $H_{n+\frac{1}{2}}^{(1)}(\cdot)$ is the Hankel function of the first kind with order $n+\frac{1}{2}$ .

Define the DtN operator $T:H^{\frac{1}{2}}(\partial B_{R})\rightarrow H^{-\frac{1}{2}}(\partial B_{R})$ by

(Tu)(R,\hat{x})=\frac{1}{R}\sum_{n=0}^{\infty}\Theta_{n}(\kappa R)\sum_{m=-n}^{n}\hat{u}_{n}^{m}Y_{n}^{m}(\hat{x}),

(2.4)

where

\Theta_{n}(z)=z\frac{{h_{n}^{(1)}}^{{}^{\prime}}(z)}{h_{n}^{(1)}(z)},

which satisfies (cf. [14, 19]):

\Re\Theta_{n}(z)\leq-\frac{1}{2},\quad\Im\Theta_{n}(z)>0,\quad\Theta_{n}(z)\sim n,\quad n\rightarrow\infty.

(2.5)

The DtN operator has the following properties. The proof is similar to that of [21, Lemma 1.2] and is omitted here for brevity.

Lemma 2.1.

The DtN operator $T:H^{\frac{1}{2}}(\partial B_{R})\to H^{-\frac{1}{2}}(\partial B_{R})$ is continuous, i.e.,

\|Tu\|_{H^{-\frac{1}{2}}(\partial B_{R})}\lesssim\|u\|_{H^{\frac{1}{2}}(\partial B_{R})}.

Moreover, it satisfies

-\Re\langle Tu,u\rangle\gtrsim\|u\|_{L^{2}(\partial B_{R})}^{2},\quad\Im\langle Tu,u\rangle\geq 0.

Here $a\lesssim b$ or $a\gtrsim b$ stands for $a\leq Cb$ or $a\geq Cb$ , where $C$ is a positive constant whose specific value is not required and may be different in the context.

It follows from (2.3)–(2.4) that we have the transparent boundary condition

\partial_{r}u=Tu\quad{\rm on}~\partial B_{R}.

(2.6)

The weak formulation of (2.1) is to find $u\in H^{1}(\Omega)$ such that

a(u,v)=\langle g,v\rangle_{\partial D}\quad\forall\,v\in H^{1}(\Omega),

(2.7)

where the sesquilinear form $a:~H^{1}(\Omega)\times H^{1}(\Omega)\rightarrow\mathbb{C}$ is defined by

a(u,v)=\int_{\Omega}\nabla u\cdot\nabla\bar{v}{\rm d}x-\kappa^{2}\int_{\Omega}u\bar{v}{\rm d}x-\langle Tu,v\rangle_{\partial B_{R}}

and the linear functional

\langle g,v\rangle_{\partial D}=\int_{\partial D}g\bar{v}{\rm d}s.

Theorem 2.2.

The variational problem (2.7) has at most one solution.

Proof.

It suffices to show that $u=0$ if $g=0$ . By (2.7), we have

\int_{\Omega}(|\nabla u|^{2}-\kappa^{2}|u|^{2}){\rm d}x-\langle Tu,u\rangle_{\partial B_{R}}=0.

Taking the imaginary part of the above equation yields

\Im\langle Tu,u\rangle_{\partial B_{R}}=R\sum_{n=0}^{\infty}\sum_{m=-n}^{m}\Im\Theta_{n}(\kappa R)|\hat{u}_{n}^{m}|^{2}=0,

which gives that $\hat{u}_{n}^{m}=0$ by (2.5). Thus we have from (2.3) and (2.6) that $u=0$ and $\partial_{r}u=0$ on $\partial B_{R}$ . We conclude from the Holmgren uniqueness theorem and the unique continuation [20] that $u=0$ on $\Omega$ . ∎

Theorem 2.3.

The variational problem $\rm(\ref{wf})$ admits a unique weak solution $u$ in $H^{1}(\Omega)$ . Furthermore, there is a positive constant $C$ depending on $\kappa$ and $R$ such that

\|u\|_{H^{1}(\Omega)}\leq C\|g\|_{L^{2}(\partial D)}.

Proof.

First we show that there exists $u_{0}\in H^{1}(\Omega)$ such that $\partial_{\nu}u_{0}=g$ on $\partial D$ and

\|u_{0}\|_{H^{1}(\Omega)}\leq C\|g\|_{L^{2}(\partial D)}.

For example, $u_{0}$ can be chosen as the unique weak solution of the following boundary value problem

\begin{cases}-\Delta u_{0}+u_{0}=0\quad&\text{in}~\Omega,\\ \partial_{\nu}u_{0}=-g\quad&\text{on}~\partial D,\\ u_{0}=0\quad&\text{on}~\partial B_{R}.\end{cases}

Next is consider the variational problem: Find $u\in H^{1}(\Omega)$ such that $u-u_{0}\in H^{1}(\Omega)$ and

a(u-u_{0},v)=\langle g,v\rangle_{\partial B_{R}}-a(u_{0},v)\quad\forall\,v\in H^{1}(\Omega).

Denote $(h,v)=\langle g,v\rangle_{\partial B_{R}}-a(u_{0},v)$ . The above problem is equivalent to the following variational problem: Find $w\in H^{1}(\Omega)$ such that

a(w,v)=(h,v)\quad\forall\,v\in H^{1}(\Omega).

(2.8)

Let $a(w,v)=a_{+}(w,v)-\kappa^{2}(w,v)$ , where

a_{+}(w,v)=\int_{\Omega}\nabla w\cdot\nabla\bar{v}{\rm d}x-\langle Tw,v\rangle_{\partial B_{R}}.

By Lemma 2.1, we have

|a_{+}(v,v)|\geq\|\nabla v\|_{L^{2}(\Omega)}^{2}+|\Re\langle Tv,v\rangle_{\partial B_{R}}|\gtrsim\|\nabla v\|_{L^{2}(\Omega)}^{2}+\|v\|_{L^{2}(\partial B_{R})}^{2},

which implies from the Friedrichs inequality that

|a_{+}(v,v)|\gtrsim\|v\|_{H^{1}(\Omega)}^{2}\quad\forall\,v\in H^{1}(\Omega).

Let $Q:L^{2}(\Omega)\rightarrow H^{1}(\Omega)$ be a linear map defined by

a_{+}(Qw,v)=(w,v)\quad\forall\,v\in H^{1}(\Omega).

By the Lax–Milgram lemma, we obtain that $Q$ is bounded from $L^{2}(\Omega)$ to $H^{1}(\Omega)$ , i.e., it satisfies

\|Qw\|_{H^{1}(\Omega)}\lesssim\|w\|_{L^{2}(\Omega)}.

(2.9)

It is clear to note that (2.8) is equivalent to the operator equation

(I-\kappa^{2}Q)w=Qh.

Due to the compact embedding of $H^{1}(\Omega)$ into $L^{2}(\Omega)$ , the operator $Q$ is a compact from $L^{2}(\Omega)$ to $L^{2}(\Omega)$ . It follows from the Fredholm alternative theorem and Theorem 2.2 that the operator $I-\kappa^{2}Q$ has a bounded inverse. Hence we have

\|w\|_{L^{2}(\Omega)}\lesssim\|Qh\|_{L^{2}(\Omega)}\lesssim\|h\|_{L^{2}(\Omega)}.

(2.10)

Combining (2.9) and (2.10) yields that

\|w\|_{H^{1}(\Omega)}=\|Q(\kappa^{2}w+h)\|_{H^{1}(\Omega)}\lesssim\|\kappa^{2}w+h\|_{L^{2}(\Omega)}\lesssim\|h\|_{L^{2}(\Omega)}.

Since $u=u_{0}+w$ , we have

\|u\|_{H^{1}(\Omega)}\leq\|w\|_{H^{1}(\Omega)}+\|u_{0}\|_{H^{1}(\Omega)}\lesssim\|h\|_{L^{2}(\Omega)}+\|u_{0}\|_{L^{2}(\Omega)}.

It is easy to deduce from the definition of $h$ that

\|u\|_{H^{1}(\Omega)}\lesssim\|g\|_{L^{2}(\Omega)},

which completes the proof. ∎

By the general theory in Babuska and Aziz [1], there exists a constant $\gamma>0$ depending on $\kappa$ and $R$ such that the following inf-sup condition holds:

\sup_{0\neq v\in H^{1}(\Omega)}\frac{|a(u,v)|}{\|v\|_{H^{1}(\Omega)}}\geq\gamma\|u\|_{H^{1}(\Omega)}\quad\forall\,u\in H^{1}(\Omega).

3. Finite element approximation

In this section, we introduce the finite element approximation of (2.7) and present the a posteriori error estimate, which plays an important role in the adaptive finite element method.

Let $\mathcal{M}_{h}$ be a regular tetrahedral mesh of the domain $\Omega$ , where $h$ represents the maximum diameter of all the elements in $\mathcal{M}_{h}$ . In order to avoid using the isoparametric finite element space and discussing the approximation error of the boundaries $\partial D$ and $\partial B_{R}$ , we assume for simplicity that $\partial D$ and $\partial B_{R}$ are polyhedral. Thus any face $F\in\mathcal{M}_{h}$ is a subset of $\partial\Omega$ if it has three boundary vertices.

Let $V_{h}\subset H^{1}(\Omega)$ be a conforming finite element space, i.e.,

V_{h}:=\{v_{h}\in C(\overline{\Omega}):v_{h}|_{K}\in P_{m}(K)\quad\forall~K\in\mathcal{M}_{h}\},

where $m$ is a positive integer and $P_{m}(K)$ denotes the set of all polynomials of degree no more than $m$ . The finite element approximation to (2.7) is to seek $u_{h}\in V_{h}$ satisfying

a(u_{h},v_{h})=\langle g,v_{h}\rangle_{\partial D}\quad\forall\,v_{h}\in V_{h}.

The above variational problem involves the DtN operator $T$ defined by an infinite series in (2.4). Practically, it is necessary to truncate the infinite series by taking finitely many terms of the expansion in order to apply the finite element method. Given a positive integer $N$ , we define the truncated DtN operator

(T_{N}u)(R,\hat{x})=\frac{1}{R}\sum_{n=0}^{N}\Theta_{n}(\kappa R)\sum_{m=-n}^{m=n}\hat{u}_{n}^{m}Y_{n}^{m}(\hat{x}).

Using the truncated DtN operator $T_{N}$ , we have the truncated finite element approximation to the problem (2.7): Find $u_{h}^{N}\in V_{h}$ such that

a_{N}(u_{h}^{N},v_{h})=\langle g,v_{h}^{N}\rangle_{\partial D}\quad\forall\,v_{h}\in V_{h},

(3.1)

where the sesquilinear form $a_{N}:~V_{h}\times V_{h}\rightarrow\mathbb{C}$ is defined by

a_{N}(u,v)=\int_{\Omega}\nabla u\cdot\nabla\bar{v}{\rm d}x-\kappa^{2}\int_{\Omega}u\bar{v}{\rm d}x-\langle T_{N}u,v\rangle_{\partial B_{R}}.

(3.2)

By the argument of Schatz [28], the discrete inf-sup condition of the sesquilinear form $a_{N}$ may be established for sufficiently large $N$ and sufficiently small $h$ . It follows from the general theory in [1] that the truncated variational problem (3.1) admits a unique solution. In this work, our goal is to obtain the a posteriori error estimate and develop the associated adaptive algorithm. Thus we assume that the discrete problem (3.1) has a unique solution $u_{h}^{N}\in V_{h}$ .

4. A posteriori error analysis

First, we collect some relevant results from [24] on the Hankel functions. Let $j_{n}(t)$ and $y_{n}(t)$ be the spherical Bessel functions of the first and second kind with order $n$ , respectively. The spherical Hankel functions are

h_{n}^{(j)}(t)=j_{n}(t)\pm{\rm i}y_{n}(t),\quad j=1,2.

For fixed $t$ , the spherical Bessel functions admit the asymptotic expressions (cf. [24, Theorem 2.31])

j_{n}(t)\sim\frac{t^{n}}{(2n+1)!!},\quad y_{n}(t)\sim-\frac{(2n-1)!!}{t^{n+1}},\quad n\to\infty,

which give that

h_{n}^{(j)}(t)\sim(-1)^{j}{\rm i}\frac{(2n-1)!!}{t^{n+1}},\quad n\to\infty.

(4.1)

For any $K\in\mathcal{M}_{h}$ , let $\mathcal{B}_{F}$ represent the set of all the faces of $K$ . Denote by $h_{K}$ and $h_{F}$ the sizes of element $K$ and face $F$ , respectively. For any interior face $F$ which is the common part of elements $K_{1}$ and $K_{2}$ , we define the jump residual across $F$ as

J_{F}=-(\nabla u_{h}^{N}|_{K_{1}}\cdot\nu_{1}+\nabla u_{h}^{N}|_{K_{2}}\cdot\nu_{2}),

where $\nu_{j}$ is the unit normal vector to the boundary of $K_{j},j=1,2$ . For any boundary face $F\in\partial B_{R}$ , we define the jump residual

J_{F}=2(Tu_{h}^{N}+\nabla u_{h}^{N}\cdot\nu),

where $\nu$ is the unit outward normal on $\partial B_{R}$ . For any boundary face $F\in\partial D$ , we define the jump residual

J_{F}=2(\nabla u_{h}^{N}\cdot\nu+g),

where $\nu$ is the unit outward normal on $\partial D$ pointing toward $\Omega$ . For any $K\in\mathcal{M}_{h}$ , denote by $\eta_{K}$ the local error estimator, which is defined by

\eta_{K}=h_{K}\|(\Delta+\kappa^{2})u_{h}^{N}\|_{L^{2}(K)}+\Big(\frac{1}{2}\sum_{F\in\partial K}h_{F}\|J_{F}\|_{L^{2}(F)}^{2}\Big)^{\frac{1}{2}}.

We now state the main result, which plays an important role for the numerical experiments.

Theorem 4.1.

Let $u$ and $u_{h}^{N}$ be the solutions of (2.7) and (3.1), respectively. There exists a positive integer $N_{0}$ independent of $h$ such that the following a posteriori error estimate holds for $N>N_{0}$ :

\|u-u_{h}^{N}\|_{H^{1}(\Omega)}\lesssim\bigg(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\bigg)^{\frac{1}{2}}+\bigg(\frac{R^{\prime}}{R}\bigg)^{N}\|g\|_{L^{2}(\partial D)}.

It can be seen from Theorem 4.1 that the a posteriori error consists of two parts: the first part comes from the finite element discretization error and the second part accounts for the truncation error of the DtN operator, which decays exponentially with respect to $N$ since $R^{\prime}<R$ . We point out that the constant in the estimate may depend on $\kappa,R$ and $R^{\prime}$ , but does not depend on the truncation parameter of the DtN operator $N$ or the mesh size of the triangulation $h$ .

In the rest part of this section, we prove the a posteriori error estimator in Theorem 4.1 by using a duality argument.

Denote the error $\xi:=u-u_{h}^{N}$ . Introduce a dual problem to the original scattering problem: Find $w\in H^{1}(\Omega)$ such that

a(v,w)=(v,\xi)\quad\forall\,v\in H^{1}(\Omega).

(4.2)

It is easy to verify that $w$ satisfies the following boundary value problem:

\begin{cases}\Delta w+\kappa^{2}w=-\xi\quad&\text{in}~\mathbb{R}^{3}\setminus\overline{D},\\ \partial_{\nu}w=0\quad&\text{on}~\partial D,\\ \partial_{r}w-T^{*}w=0\quad&\text{on}~\partial B_{R},\end{cases}

(4.3)

where the adjoint operator $T^{*}$ is defined by

(T^{*}u)(R,\hat{x})=\frac{1}{R}\sum_{n=0}^{\infty}\overline{\Theta}_{n}(\kappa R)\sum_{m=-n}^{m=n}\hat{u}_{n}^{m}Y_{n}^{m}(\hat{x}).

We may follow the same proof as that for the original scattering problem (2.1) and show that the dual problem (4.3) has a unique weak solution $w\in H^{1}(\Omega)$ , which satisfies

\|w\|_{H^{1}(\Omega)}\lesssim\|\xi\|_{L^{2}(\Omega)}.

The following lemma gives the error representation formulas and is the basis for the a posteriori error analysis.

Lemma 4.2.

Let $u$ , $u_{h}^{N}$ and $w$ be the solutions of the problems $(\rm\ref{wf})$ , $(\rm\ref{fem})$ and $(\rm\ref{bi})$ , respectively. The following identities hold:

	$\displaystyle\\|\xi\\|^{2}_{H^{1}(\Omega)}=\Re\left(a(\xi,\xi)+\langle(T-T_{N})\xi,\xi\rangle_{\partial B_{R}}\right)+\Re\langle T_{N}\xi,\xi\rangle_{\partial B_{R}}+(\kappa^{2}+1)\\|\xi\\|^{2}_{L^{2}(\Omega)},$		(4.4)
	$\displaystyle\\|\xi\\|^{2}_{L^{2}(\Omega)}=a(\xi,w)+\langle(T-T_{N})\xi,w\rangle_{\partial B_{R}}-\langle(T-T_{N})\xi,w\rangle_{\partial B_{R}},$		(4.5)
	$\displaystyle a(\xi,\psi)+\langle(T-T_{N})\xi,\psi\rangle_{\partial B_{R}}=\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h})$
	$\displaystyle\hskip 113.81102pt+\langle(T-T_{N})u,\psi\rangle_{\partial B_{R}}\quad\forall\,\psi\in H^{1}(\Omega),\psi_{h}\in V_{h}.$		(4.6)

Proof.

The equality (4.4) follows directly from the definition of the sesquilinear form $a$ in (2.7). The identity (4.5) can be easily deduced by taking $v=\xi$ in (4.2). It remains to prove (4.6). It follows from (2.7) and (3.1) that

$\displaystyle a(\xi,\psi)$	$\displaystyle=$	$\displaystyle a(u-u_{h}^{N},\psi-\psi_{h})+a(u-u_{h}^{N},\psi_{h})$
	$\displaystyle=$	$\displaystyle\langle g,\psi-\psi_{h}\rangle_{\partial D}-a(u_{h}^{N},\psi-\psi_{h})+a(u-u_{h}^{N},\psi_{h})$
	$\displaystyle=$	$\displaystyle\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h})$
		$\displaystyle+a_{N}(u_{h}^{N},\psi-\psi_{h})-a(u_{h}^{N},\psi-\psi_{h})+a(u,\psi_{h})-a(u_{h}^{N},\psi_{h}).$

Since $a(u,\psi_{h})=\langle g,\psi_{h}\rangle_{\partial D}=a_{N}(u_{h}^{N},\psi_{h})$ , we have

	$\displaystyle a(\xi,\psi)$	$\displaystyle=\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h})+a_{N}(u_{h}^{N},\psi)-a(u_{h}^{N},\psi)$
		$\displaystyle=\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h})+\langle(T-T_{N})u_{h}^{N},\psi\rangle_{\partial B_{R}}$
		$\displaystyle=\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h})-\langle(T-T_{N})\xi,\psi\rangle_{\partial B_{R}}$
		$\displaystyle\quad+\langle(T-T_{N})u,\psi\rangle_{\partial B_{R}},$

which implies (4.6) and completes the proof. ∎

It is necessary to estimate (4.6) and the last term in (4.5) in order to prove Theorem 4.1. We begin with a trace regularity result.

Lemma 4.3.

For any $u\in H^{1}(\Omega)$ , the following estimates hold:

\|u\|_{H^{\frac{1}{2}}(\partial B_{R})}\lesssim\|u\|_{H^{1}(\Omega)},\quad\|u\|_{H^{\frac{1}{2}}(\partial B_{R^{\prime}})}\lesssim\|u\|_{H^{1}(\Omega)}.

Proof.

Let

B_{R}\setminus\overline{B}_{R^{\prime}}=\{(r,\theta,\varphi):0<R^{\prime}<r<R,~0<\theta<\pi,~0<\varphi<2\pi\}.

It is shown in [22, Lemma 2] that

(1+n^{2})^{\frac{1}{2}}|\zeta(R)|^{2}\lesssim\int_{R^{\prime}}^{R}\left((1+n^{2})|\zeta(r)|^{2}+|\zeta^{\prime}(r)|^{2}\right){\rm d}r,

which gives after combining (2.2) that

	$\displaystyle\\|u\\|^{2}_{H^{\frac{1}{2}}(\partial B_{R})}$	$\displaystyle=$	$\displaystyle\sum_{n=0}^{\infty}(1+n^{2})^{\frac{1}{2}}\sum\limits_{m=-n}^{n}\|\hat{u}_{n}^{m}(R)\|^{2}$
		$\displaystyle\lesssim$	$\displaystyle\sum_{n=0}^{\infty}\sum_{\|m\|\leq n}\int_{R^{\prime}}^{R}\big((1+n^{2})\|\hat{u}_{n}^{m}(r)\|^{2}+\|\hat{u}_{n}^{m^{\prime}}(r)\|^{2}\big){\rm d}r.$

Noting the fact (cf. [24, Theorem 5.34])

\|u\|^{2}_{H^{1}({B}_{R}\setminus\overline{B}_{R^{\prime}})}\geq\sum_{n=0}^{\infty}\sum_{|m|\leq n}\int_{R^{\prime}}^{R}r^{2}\Big[\Big(1+\frac{n(n+1)}{r^{2}}\Big)|\hat{u}_{n}^{m}(r)|^{2}+|\hat{u}_{n}^{m^{\prime}}(r)|^{2}\Big]{\rm d}r,

we obtain

\|u\|^{2}_{H^{\frac{1}{2}}(\partial B_{R})}\lesssim\|u\|^{2}_{H^{1}(B_{R}\setminus\overline{B}_{R^{\prime}})}\leq\|u\|^{2}_{H^{1}(\Omega)},

which shows the first inequality. The second inequality can be proved similarly by observing the identity

(R-R^{\prime})|\zeta(R^{\prime})|^{2}=\int_{R^{\prime}}^{R}|\zeta(r)|^{2}{\rm d}r+\int_{R^{\prime}}^{R}\int_{r}^{R^{\prime}}\frac{{\rm d}}{{\rm d}r}|\zeta(r)|^{2}{\rm d}t{\rm d}r.

The details are omitted here. ∎

Lemma 4.4.

Let $u$ be the solution to $\rm(\ref{wf})$ . Then the following estimate holds:

|\hat{u}_{n}^{m}(R)|\lesssim\Big(\frac{R^{\prime}}{R}\Big)^{n}|\hat{u}_{n}^{m}(R^{\prime})|.

Proof.

It is known that the solution of the scattering problem (2.1) admits the series expansion

u(r,\hat{x})=\sum\limits_{n=0}^{\infty}\sum\limits_{m=-n}^{n}\frac{h_{n}^{(1)}(\kappa r)}{h_{n}^{(1)}(\kappa R^{\prime})}\hat{u}_{n}^{m}(R^{\prime})Y_{n}^{m}(\hat{x}),\quad\hat{u}_{n}^{m}(R^{\prime})=\int_{\mathbb{S}^{2}}u(R^{\prime},\hat{x})Y_{n}^{m}(\hat{x}){\rm d}\hat{x}

(4.7)

for all $r>R^{\prime}$ . Evaluating (4.7) at $r=R$ yields

u(R,\hat{x})=\sum\limits_{n=0}^{\infty}\sum\limits_{m=-n}^{n}\frac{h_{n}^{(1)}(\kappa R)}{h_{n}^{(1)}(\kappa R^{\prime})}\hat{u}_{n}^{m}(R^{\prime})Y_{n}^{m}(\hat{x}),

which implies

\hat{u}_{n}^{m}(R)=\frac{h_{n}^{(1)}(kR)}{h_{n}^{(1)}(kR^{\prime})}\hat{u}_{n}^{m}(R^{\prime}).

Using the asymptotic expression in (4.1), we obtain

|\hat{u}_{n}^{m}(R)|=\Biggl|\frac{h_{n}^{(1)}(kR)}{h_{n}^{(1)}(kR^{\prime})}\Biggr||\hat{u}_{n}^{m}(R^{\prime})|\lesssim\Big(\frac{R^{\prime}}{R}\Big)^{n}|\hat{u}_{n}^{m}(R^{\prime})|,

which completes the proof. ∎

Lemma 4.5.

For any $\psi\in H^{1}(\Omega)$ , the following estimate holds:

|a(\xi,\psi)+\langle(T-T_{N})\xi,\psi\rangle_{\partial B_{R}}|\lesssim\bigg(\Big(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}+\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}\bigg)\|\psi\|_{H^{1}(\Omega)}.

Proof.

Define

	$\displaystyle J_{1}$	$\displaystyle=$	$\displaystyle\langle g,\psi-\psi_{h}\rangle_{\partial D}-a_{N}(u_{h}^{N},\psi-\psi_{h}),$
	$\displaystyle J_{2}$	$\displaystyle=$	$\displaystyle\langle(T-T_{N})u,\psi\rangle_{\partial B_{R}},$

where $\psi_{h}\in V_{h}$ . It follows from (4.6) that

a(\xi,\psi)+\langle(T-T_{N})\xi,\psi\rangle_{\partial B_{R}}=J_{1}+J_{2}.

Using (3.2) and the integration by parts, we obtain

J_{1}=\sum_{K\in\mathcal{M}_{h}}\biggl(\int_{K}(\Delta u_{h}^{N}+\kappa^{2}u_{h}^{N})(\bar{\psi}-\bar{\psi}_{h}){\rm d}x+\sum_{F\in\partial K}\frac{1}{2}\int_{F}J_{F}(\bar{\psi}-\bar{\psi}_{h}){\rm d}s\biggl).

Now we take $\psi_{h}=\Pi_{h}\psi\in V_{h}$ , where $\Pi_{h}$ is the Scott–Zhang interpolation operator and has the approximation properties

\|v-\Pi_{h}v\|_{L^{2}(K)}\lesssim h_{K}\|\nabla v\|_{L^{2}(\tilde{K})},\quad\|v-\Pi_{h}v\|_{L^{2}(F)}\lesssim h_{F}^{\frac{1}{2}}\|\nabla v\|_{L^{2}(\tilde{K}_{F})},

Here $\tilde{K}$ and $\tilde{K}_{F}$ are the union of all the elements in $\mathcal{M}_{h}$ , which have nonempty intersection with element $K$ and the face $F$ , respectively.

By the Cauchy–Schwarz inequality, we have

$\displaystyle\|J_{1}\|$	$\displaystyle\lesssim$	$\displaystyle\sum\limits_{K\in\mathcal{M}_{h}}\bigg(h_{K}\\|(\Delta+\kappa^{2})u_{h}^{N}\\|_{L^{2}(K)}\\|\nabla\psi\\|_{L^{2}(\tilde{K})}+\sum\limits_{F\in\partial K}\frac{1}{2}h_{F}^{\frac{1}{2}}\\|J_{F}\\|_{L^{2}(F)}\\|\psi\\|_{H^{1}(\tilde{K}_{F})}\bigg)$
	$\displaystyle\lesssim$	$\displaystyle\sum_{K\in\mathcal{M}_{h}}\bigg[h_{K}\\|(\Delta+\kappa^{2})u_{h}^{N}\\|_{L^{2}(K)}+\bigg(\sum_{F\in\partial K}\frac{1}{2}h_{F}\\|J_{F}\\|_{L^{2}(F)}^{2}\bigg)^{\frac{1}{2}}\bigg]\\|\psi\\|_{H^{1}(\Omega)}$
	$\displaystyle\lesssim$	$\displaystyle\bigg(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\bigg)^{2}\\|\psi\\|_{H^{1}(\Omega)}.$

It follows from the definitions of $T,T_{N}$ and Lemma 4.4 that

$\displaystyle\|J_{2}\|=\|\langle(T-T_{N})u,\psi\rangle_{\partial B_{R}}\|$	$\displaystyle=$	$\displaystyle\Big\|R\sum_{n>N}\sum_{\|m\|\leq n}\Theta_{n}(\kappa R)\hat{u}_{n}^{m}(R)\bar{\hat{\psi}}_{n}^{m}(R)\Big\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R)\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\Big\|\Big(\frac{R^{\prime}}{R}\Big)^{N}\hat{u}_{n}^{m}(R^{\prime})\Big\|\|\bar{\hat{\psi}}_{nm}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|.$

Using (2.5), the Cauchy–Schwarz inequality, and Lemma 4.3 yields

$\displaystyle\|J_{2}\|$	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\Big(\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|^{2}\Big)^{\frac{1}{2}}\Big(\sum_{\|m\|\leq n}\|\bar{\hat{\psi}}_{n}^{m}(R)\|^{2}\Big)^{\frac{1}{2}}$
	$\displaystyle\leq$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\Big(\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|^{2}\Big)^{\frac{1}{2}}\Big(\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\bar{\hat{\psi}}_{n}^{m}(R)\|^{2}\Big)^{\frac{1}{2}}$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\\|u\\|_{H^{\frac{1}{2}}(\partial B_{R^{\prime}})}\\|\psi\\|_{H^{\frac{1}{2}}(\partial B_{R})}$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\\|u\\|_{H^{1}(\Omega)}\\|\psi\\|_{H^{\frac{1}{2}}(\partial B_{R})}.$

Combining the above estimates gives

|J_{1}|+|J_{2}|\lesssim\bigg(\Big(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}+\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}\bigg)\|\psi\|_{H^{1}(\Omega)},

which completes the proof. ∎

Lemma 4.6.

Let $w$ be the solution to the dual problem $\rm(\ref{bi})$ . Then the following estimate holds:

|\langle(T-T_{N})\xi,w\rangle_{\partial B_{R}}|\lesssim N^{-2}\|\xi\|^{2}_{H^{1}(\Omega)}.

Proof.

It follows from (2.5), Lemma 4.3 and the Cauchy–Schwarz inequality that we have

$\displaystyle\|\langle(T-T_{N})\xi,w\rangle_{\partial B_{R}}\|$	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|n\|\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle=$	$\displaystyle\sum_{n>N}((1+n^{2})^{\frac{1}{2}}n^{3})^{-\frac{1}{2}}\sum_{\|m\|\leq n}(1+n^{2})^{\frac{1}{4}}n^{\frac{5}{2}}\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle N^{-2}\sum_{n>N}\sum_{\|m\|\leq n}(1+n^{2})^{\frac{1}{4}}n^{\frac{5}{2}}\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle N^{-2}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}(1+n(n+1))^{\frac{1}{2}}\|\hat{\xi}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}$
	$\displaystyle=$	$\displaystyle N^{-2}\\|\xi\\|_{H^{\frac{1}{2}}(\partial B_{R})}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}$
	$\displaystyle\lesssim$	$\displaystyle N^{-2}\\|\xi\\|_{H^{1}(\Omega)}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}.$

To estimate $\hat{w}_{n}^{m}(R)$ , we consider the dual problem (4.3) in the annulus $B_{R}\setminus B_{R^{\prime}}$ :

\begin{cases}\Delta w+\kappa^{2}w=-\xi\quad&\text{in}~B_{R}\setminus\overline{B}_{R^{\prime}},\\ w=w(R^{\prime},\hat{x})\quad&\text{on}~\partial R^{\prime},\\ \partial_{r}w-T^{*}w=0\quad&\text{on}~\partial B_{R},\end{cases}

which reduces to the second order equation for the coefficients $\hat{w}_{n}^{m}$ in the Fourier domain

\begin{cases}\frac{{\rm d}^{2}\hat{w}_{n}^{m}(r)}{{\rm d}r^{2}}+\frac{2}{r}\frac{{\rm d}\hat{w}_{n}^{m}(r)}{{\rm d}r}+\big(\kappa^{2}-\frac{n(n+1)}{r^{2}}\big)\hat{w}_{n}^{m}(r)=-\hat{\xi}_{n}^{m}(r),&R^{\prime}<r<R,\\ \frac{{\rm d}\hat{w}_{n}^{m}(R)}{{\rm d}r}-\frac{1}{R}\overline{\Theta}_{n}(\kappa R)\hat{w}_{n}^{m}(R)=0,&r=R,\\ \hat{w}_{n}^{m}(R^{\prime})=\hat{w}_{n}^{m}(R^{\prime}),&r=R^{\prime}.\end{cases}

By the method of the variation of parameters, we obtain the solution of the above equation

\hat{w}_{n}^{m}(r)=S_{n}(r)\hat{w}_{n}^{m}(R^{\prime})+\frac{{\rm i\kappa}}{2}\int_{R^{\prime}}^{r}t^{2}W_{n}(r,t)\hat{\xi}_{n}^{m}(t){\rm d}t+\frac{{\rm i\kappa}}{2}\int_{R^{\prime}}^{R}t^{2}S_{n}(t)W_{n}(R^{\prime},r)\hat{\xi}_{n}^{m}(t){\rm d}t,

(4.8)

where

S_{n}(r)=\frac{h_{n}^{(2)}(\kappa r)}{h_{n}^{(2)}(\kappa R^{\prime})},\quad W_{n}(r,t)=\det\left[\begin{matrix}h_{n}^{(1)}(\kappa r)&h_{n}^{(2)}(\kappa r)\\ h_{n}^{(1)}(\kappa t)&h_{n}^{(2)}(\kappa t)\end{matrix}\right].

Taking $r=R$ in (4.8), we get

\hat{w}_{n}^{m}(R)=S_{n}(R)\hat{w}_{n}^{m}(R^{\prime})+\frac{{\rm i}\kappa}{2}\int_{R^{\prime}}^{R}t^{2}S_{n}(R)W_{n}(R^{\prime},t)\hat{\xi}_{n}^{m}(t){\rm d}t.

Using the asymptotic expression (4.1) yields

S_{n}(R)\sim\Big(\frac{R^{\prime}}{R}\Big)^{n},\quad n\rightarrow\infty

and

	$\displaystyle W_{n}(R^{\prime},t)$	$\displaystyle=$	$\displaystyle 2{\rm i}j_{n}(\kappa R^{\prime})y_{n}(\kappa R^{\prime})\bigg(\frac{j_{n}(\kappa t)}{j_{n}(\kappa R^{\prime})}-\frac{y_{n}(\kappa t)}{y_{n}(\kappa R^{\prime})}\bigg)$
		$\displaystyle\sim$	$\displaystyle-\frac{2{\rm i}}{(2n+1)\kappa R^{\prime}}\bigg(\Big(\frac{t}{R^{\prime}}\Big)^{n}-\Big(\frac{R^{\prime}}{t}\Big)^{n+1}\bigg),\quad n\rightarrow\infty.$

Hence

|S_{n}(R)|\lesssim\Big(\frac{R^{\prime}}{R}\Big)^{n},\quad|W_{n}(R^{\prime},t)|\lesssim n^{-1}\Big(\frac{t}{R^{\prime}}\Big)^{n}.

Combining the above estimates, we obtain

$\displaystyle\|\hat{w}_{n}^{m}(R)\|$	$\displaystyle\leq$	$\displaystyle\|S_{n}(R)\|\|\hat{w}_{n}^{m}(R^{\prime})\|+\frac{\kappa}{2}\int_{R^{\prime}}^{R}t^{2}\|S_{n}(R)\|\|W_{n}(R^{\prime},t)\|\|\hat{\xi}_{n}^{m}(t)\|{\rm d}t,$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{n}\|\hat{w}_{n}^{m}(R^{\prime})\|+n^{-1}\Big(\frac{R^{\prime}}{R}\Big)^{n}\\|\hat{\xi}_{n}^{m}(t)\\|_{L^{\infty}([R^{\prime},R])}\int_{R^{\prime}}^{R}t^{2}\Big(\frac{t}{R^{\prime}}\Big)^{n}{\rm d}t$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{n}\|\hat{w}_{n}^{m}(R^{\prime})\|+n^{-2}\\|\hat{\xi}_{n}^{m}(t)\\|_{L^{\infty}([R^{\prime},R])},$

which gives

$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}$	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\bigg(\Big(\frac{R^{\prime}}{R}\Big)^{n}\|\hat{w}_{n}^{m}(R^{\prime})\|+n^{-2}\\|\hat{\xi}_{n}^{m}(t)\\|_{L^{\infty}([R^{\prime},R])}\bigg)^{2}$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\bigg(\Big(\frac{R^{\prime}}{R}\Big)^{2n}\|\hat{w}_{n}^{m}(R^{\prime})\|^{2}+n^{-4}\\|\hat{\xi}_{n}^{m}(t)\\|^{2}_{L^{\infty}([R^{\prime},R])}\bigg)$
	$\displaystyle:=$	$\displaystyle I_{1}+I_{2}.$

Here

	$\displaystyle I_{1}=$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\Big(\frac{R^{\prime}}{R}\Big)^{2n}\|\hat{w}_{n}^{m}(R^{\prime})\|^{2},$
	$\displaystyle I_{2}=$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}n\\|\hat{\xi}_{n}^{m}(t)\\|^{2}_{L^{\infty}([R^{\prime},R])}.$

A simple calculation yields

	$\displaystyle I_{1}$	$\displaystyle\lesssim$	$\displaystyle\max_{n>N}n^{4}\Big(\frac{R^{\prime}}{R}\Big)^{2n}\sum_{n>N}\sum_{\|m\|\leq n}n~\|\hat{w}_{n}^{m}(R^{\prime})\|^{2}\lesssim\sum_{n>N}\sum_{\|m\|\leq n}n~\|\hat{w}_{n}^{m}(R^{\prime})\|^{2}$
		$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}(1+n^{2})^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\hat{w}_{n}^{m}(R^{\prime})\|^{2}\leq\\|w\\|^{2}_{H^{\frac{1}{2}}(\partial B_{R^{\prime}})}\lesssim\\|\xi\\|^{2}_{H^{1}(\Omega)}.$

By [22, Lemma 5], we have

\|\hat{\xi}_{n}^{m}(t)\|^{2}_{L^{\infty}([R^{\prime},R])}\leq\Big(\frac{2}{\delta}+n\Big)\|\hat{\xi}_{n}^{m}(t)\|^{2}_{L^{2}([R^{\prime},R])}+n^{-1}\|\hat{\xi}_{n}^{m^{\prime}}(t)\|^{2}_{L^{2}([R^{\prime},R])},

where $\delta=R-R^{\prime}$ . Following a similar proof of Lemma 4.2 yields

\begin{split}\|\xi\|^{2}_{H^{1}(B_{R}\setminus\overline{B}_{R})}&\geq\sum_{n=0}^{\infty}\sum_{|m|\leq n}\int_{R^{\prime}}^{R}\left[(r^{2}+n(n+1))|\xi_{n}^{m}(r)|^{2}+r^{2}|\xi_{n}^{m^{\prime}}(r)|^{2}\right]{\rm d}r\\ &\geq\sum_{n=0}^{\infty}\sum_{|m|\leq n}\int_{R^{\prime}}^{R}\left[({R^{\prime}}^{2}+n(n+1))|\xi_{n}^{m}(r)|^{2}+{R^{\prime}}^{2}|\xi_{n}^{m^{\prime}}(r)|^{2}\right]{\rm d}r,\end{split}

which gives

I_{2}=\sum_{n>N}\sum_{|m|\leq n}n\|\hat{\xi}_{nm}(t)\|^{2}_{L^{\infty}([R^{\prime},R])}\lesssim\|\xi\|^{2}_{H^{1}(\Omega)}.

Therefore, we obtain

\sum_{n>N}\sum_{|m|\leq n}n^{5}|\hat{w}_{n}^{m}(R)|^{2}\lesssim\|\xi\|^{2}_{H^{1}(\Omega)},

which completes the proof. ∎

Now we prove the main theorem.

Proof.

We conclude from (2.4)–(2.5) that

\Re\langle T_{N}\xi,\xi\rangle_{\partial B_{R}}=R\sum_{n>N}\sum_{|m|\leq n}\Re(\Theta(\kappa R))|\hat{\xi}_{n}^{m}|^{2}\leq 0.

It follows from (4.4) and Lemma 4.4 that there exist two positive constants $C_{1}$ and $C_{2}$ independent of $h$ and $N$ satisfying

\|\xi\|^{2}_{H^{1}(\Omega)}\leq C_{1}\bigg(\Big(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}+\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}\bigg)\|\xi\|_{H^{1}(\Omega)}+C_{2}\|\xi\|_{L^{2}(\Omega)}.

Using (4.5) and Lemmas 4.4–4.5, we obtain

\|\xi\|^{2}_{L^{2}(\Omega)}\leq C_{3}\bigg(\Big(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}+\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}\bigg)\|\xi\|_{L^{2}(\Omega)}+C_{4}N^{-2}\|\xi\|_{H^{1}(\Omega)},

where $C_{3}$ and $C_{4}$ are positive constants independent of $h$ and $N$ . Combining the above estimates yields

\|\xi\|^{2}_{H^{1}(\Omega)}\leq C_{5}\bigg(\Big(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}+\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}\bigg)\|\xi\|_{H^{1}(\Omega)}+C_{6}N^{-2}\|\xi\|_{H^{1}(\Omega)},

where $C_{5}$ and $C_{6}$ are positive constants independent of $h$ and $N$ . We may choose a sufficiently large integer $N_{0}$ such that $C_{6}N_{0}^{-2}<1/2$ , which completes the proof by taking $N>N_{0}$ . ∎

5. Numerical experiments

In this section, we discuss the implementation of the adaptive finite element algorithm with the truncated DtN boundary condition and present two numerical examples to demonstrate the competitive performance of the proposed method. There are two components which need to be designed carefully in order to efficiently implement the h-adaptive method. The first one is an effective management mechanism of the mesh grids. The another one is an effective indicator for the adaptivity. The a posteriori error estimate from Theorem 4.1 is used to generate the indicator in our algorithm.

5.1. The hierarchy geometry tree

In our algorithm, we use the hierarchy geometry tree (HGT) or the hierarchical grids to manage the data structure of the mesh grids [3]. The structure of the grids is described hierarchically. For example, the element such as a point for 0-dimension, an edge for 1-dimension, a triangle for 2-dimension, a tetrahedron for 3-dimension is called a geometry. If a triangle is one of the faces of a tetrahedron, then it belongs to this tetrahedron. Similarly, if an edge is one of the edges of a triangle, then it belongs to this triangle. Hence all geometries in the tetrahedrons have belonging-to relationship.

A tetrahedron $T_{0}$ can be uniformly divided into eight small sub-tetrahedrons $\{T_{0,0},T_{0,1}\cdots T_{0,7}\}$ . In this refinement operation, every face of the tetrahedron is divided into four smaller triangles. This procedure can be managed by the octree data structure which is given by Figure 1, which shows that the sub-tetrahedrons $T_{0,0}$ and $T_{0,6}$ are further divided into eight smaller sub-tetrahedrons. In the octree, we name $T_{0}$ as the root node and those nodes without further subdivision like $T_{0,1}$ and $T_{0,0,0}$ as the leaf nodes. Obviously, a set of root nodes $\{T_{i}\},i=0,1,\cdots$ can form a three-dimensional initial mesh for a domain $\Omega$ and a set of all the leaf nodes of the HGTs also form a mesh.

Figure 1. A schematic of octree data structure.

By using the HGT, the refinement and even the coarsening of a mesh can be done efficiently. However, it may cause the hanging points in the direct neighbors of the refined tetrahedrons. In order to remove these hanging points, two kinds of geometries may be introduced: twin-tetrahedron and four-tetrahedron. For the twin-tetrahedron geometry as shown in Figure 2 (left), it has five degrees of freedom (DoF) and consists of two standard tetrahedrons. To conform the finite element space, the following strategy is used to construct the basis function in twin-tetrahedron geometry. For each basis function, the value is 1 at the corresponding interpolation point and the value is 0 at the other interpolation points. For the common point of the two sub-tetrahedrons in the twin-tetrahedron like A, D and E, the support of the basis function is the whole twin-tetrahedron. For the points B and C, the support of their corresponding basis function is only the tetrahedron ABED and the tetrahedron AECD, respectively. For the four-tetrahedron as shown in Figure 2 (right), the similar strategy is used. With the twin-tetrahedron geometry and the four-tetrahedron geometry, the local refinement can be implemented easily.

Figure 2. Two geometries to avoid hanging points. (left) Twin-tetrahedron geometry. (right) Four-tetrahedron geometry.

In practice, we use a polyhedral surface to approximate $\partial D$ and $\partial B_{R}$ . Since the TBC operator is represented by the spherical harmonic functions whose accuracy depends on how good the approximation is. Obviously, a rough approximation could not satisfy the computational requirement. Based on the element geometry introduced above, we present a method to deal with the surface refinement. Suppose that the domain $\Omega$ has a curved boundary and the initial mesh is given by a rough polygon. The traditional surface refinement is performed by taking the midpoint of each side of the tetrahedron. Hence the shape of the boundary cannot be well approximated. To resolve this issue, a very simple method is adopted. When the boundary elements of the mesh need to be refined, we redefine the midpoint through projecting vertically to the desired curved boundary, as shown in Figure 3. Thanks to the HGTs, it does not spend much time at all to find these boundary elements. This method works efficiently in two-dimensions. But in three-dimensions, it may cause the neighbors to become non-standard twin-tetrahedron geometry or four-tetrahedrons geometry. To handle this problem, these special tetrahedrons, whose neighbors do not need refinement, should not redefine the midpoint. So the marked boundary tetrahedron will not be refined until the neighboring boundary tetrahedrons are marked in order to keep the mesh structure.

Figure 3. Mesh refinement on the surface (red points are redefined midpoints on the boundary).

5.2. The adaptive algorithm

The numerical simulations are implemented with a C++ library: Adaptive Finite Element Package (AFEPack). The initial mesh is generated by GMSH[15]. The resulting sparse linear systems are solving by the solver called Eigen. The simulations are implemented on a HP workstation and are accelerated by using OpenMP. The a posteriori error estimate from Theorem 4.1 is adopted to generate the indicators in our algorithm.

The error consists of two parts: the finite element discretization error $\varepsilon_{h}$ and the DtN operator truncation error $\varepsilon_{N}$ which depends on $N$ . Specifically,

\displaystyle\varepsilon_{h}=\Big(\sum\limits_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\Big)^{\frac{1}{2}}=\eta_{\mathcal{M}_{h}},\quad\varepsilon_{N}=\Big(\frac{R^{\prime}}{R}\Big)^{N}\|g\|_{L^{2}(\partial D)}.

(5.1)

In the implementation, we can choose $R^{\prime}$ , $R$ , and $N$ based on (5.1) such that finite element discretization error is not contaminated by the truncation error, i.e., $\varepsilon_{N}$ is required to be very small compared with $\varepsilon_{h}$ , for example, $\varepsilon_{N}\leq 10^{-8}$ . For simplicity, in the following numerical experiments, $R^{\prime}$ is chosen such that the scatterer lies exactly in the circle $B_{R^{\prime}}$ and $N$ is taken to be the smallest positive integer satisfying $\varepsilon_{N}\leq 10^{-8}$ . Table 1 shows the adaptive finite element algorithm with the DtN boundary condition for solving the scattering problem.

Table 1. The adaptive FEM-DtN algorithm.

1	Given a tolerance $\varepsilon>0$ ;
2	Choose $R$ , $R^{\prime}$ and $N$ such that $\varepsilon_{N}<10^{-8}$ ;
3	Construct an initial tetrahedral partition $\mathcal{M}_{h}$ over $\Omega$ and compute error estimators;
4	While $\eta_{K}>\varepsilon$ , do
5	mark $K$ , refine $\mathcal{M}_{h}$ , and obtain a new mesh $\hat{\mathcal{M}}_{h}$ .
6	solve the discrete problem on the $\hat{\mathcal{M}}_{h}$ .
7	compute the corresponding error estimators;
8	End while.

5.3. Numerical examples

We present two numerical examples to illustrate the performance of the proposed method. In the implementation, the wavenumber is $\kappa=\pi$ , which accounts for the wavelength $\lambda=2\pi/\kappa=2$ .

Example 1. Let the obstacle $D=B_{0.5}$ be the ball with a radius of 0.5 and $\Omega=B_{1}\setminus\overline{B}_{0.5}$ be the computational domain. The boundary condition $g$ is chosen such that the exact solution is

u(x)=\frac{e^{{\rm i}\kappa r}}{r},\quad r=|x|.

Refer to caption — Figure 4. Example 1: (left) initial mesh on the $x_{1}x_{2}$ -plane. (right) adaptive mesh on the $x_{1}x_{2}$ -plane.

Figure 5. Example 1: quasi-optimality of the a priori and a posteriori error estimates.

The initial mesh and an adaptive mesh is shown in Figure 4. Figure 5 displays the curves of $\log e_{h}$ and $\log\varepsilon_{h}$ versus $\log{\rm DoF}_{h}$ for our adaptive DtN method, where $e_{h}=\|\nabla(u-u_{h}^{N})\|_{L^{2}(\Omega)}$ is the a priori error, $\varepsilon_{h}$ is the a posteriori error given in (5.1), and ${\rm Dof}_{h}$ denotes the degree of freedom or the number of nodal points of the mesh $\mathcal{M}_{h}$ in the domain $\Omega$ . It indicates that the meshes and associated numerical complexity are quasi-optimal, i.e., $\|\nabla(u-u_{h}^{N})\|_{L^{2}(\Omega)}=\mathcal{O}({\rm DoF}_{h}^{-\frac{1}{3}})$ holds asymptotically.

Example 2. This example concerns the scattering of the plane wave $u^{\rm inc}=e^{{\rm i}\kappa x_{3}}$ by a U-shaped obstacle $D$ which is contained in the box $\{x\in\mathbb{R}^{3}:-0.25\leq x_{1},x_{2},x_{3}\leqslant 0.25\}$ . There is no analytical solution for this example and the solution contains singularity around the corners of the obstacle. The Neumann boundary condition is set by $g=\partial_{\nu}u^{\rm inc}$ on $\partial D$ . We take $R=1$ , $R^{\prime}=\frac{\sqrt{3}}{4}$ for the adaptive DtN method. Figure 6 shows the cross section of the obstacle and the adaptive mesh of 63898 elements, and the curve of $\log\varepsilon_{h}$ versus $\log{\rm DoF}_{h}$ . It implies that the decay of the a posteriori error estimate is $\mathcal{O}({\rm DoF}_{h}^{-1/3})$ , which is optimal.

6. Conclusion

In this paper, we have presented an adaptive finite element method with the transparent boundary condition for the three-dimensional acoustic obstacle scattering problem. The truncated DtN operator was considered for the discrete problem. A dual argument was developed in order to derive the a posteriori error estimate. The error consists of the finite element approximation error and the DtN operator truncation error which was shown to exponentially decay with respect to the truncation parameter $N$ . Numerical results show that the method is effective to solve the three-dimensional acoustic obstacle scattering problem. Possible future work is to extend the adaptive FEM-DtN method for solving the three-dimensional electromagnetic and elastic obstacle scattering problems, where the wave propagation is governed by the Maxwell equations and the Navier equation, respectively. We hope to report the progress on solving these problems elsewhere in the future.

References

[1] I. Babuška and A. Aziz, Survey Lectures on Mathematical Foundations of the Finite Element Method, in The Mathematical Foundations of the Finite Element Method with Application to the Partial Differential Equations, ed. by A. Aziz, Academic Press, New York, 1973, 5–359.
[2] G. Bao, Y. Gao, and P. Li, Time-domain analysis of an acoustic-elastic interaction problem, Arch. Ration. Mech. Anal., 229 (2018), 835–884.
[3] G. Bao, G. Hu, D. Liu, An h-adaptive finite element solver for the calculations of the electronic structures, J. Comput. Phys., 231 (2012), 4967–4979.
[4] G. Bao, P. Li, and H. Wu, An adaptive edge element method with perfectly matched absorbing layers for wave scattering by periodic structures, Math. Comp., 79 (2010), 1–34.
[5] G. Bao and H. Wu, Convergence analysis of the perfectly matched layer problems for time-harmonic Maxwell’s equations, SIAM J. Numer. Anal., 43 (2005), 2121–2143.
[6] A. Bayliss and E. Turkel, Radiation boundary conditions for numerical simulation of waves, Comm. Pure Appl. Math., 33 (1980), 707–725.
[7] J.-P. Berenger, A perfectly matched layer for the absorption of electromagnetic waves, J. Comput. Phys., 114 (1994), 185–200.
[8] Z. Chen and H. Wu, An adaptive finite element method with perfectly matched absorbing layers for the wave scattering by periodic structures, SIAM J. Numer. Anal., 41 (2003), 799–826.
[9] Z. Chen and X. Liu, An adaptive perfectly matched layer technique for time-harmonic scattering problems, SIAM J. Numer. Anal., 43 (2005), 645–671.
[10] F. Collino and P. Monk, The perfectly matched layer in curvilinear coordinates, SIAM J. Sci. Comput., 19 (1998), 2061–2090.
[11] D. Colton and R. Kress, Integral Equation Methods in Scattering Theory, John Wiley $\&$ Sons, New York, 1983.
[12] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory, Second Edition, Springer, Berlin, New York, 1998.
[13] B. Engquist and A. Majda, Absorbing boundary conditions for the numerical simulation of waves, Math. Comp., 31 (1977), 629–651.
[14] Q. Fang, D. Nicholls, and J. Shen, A stable, high-order method for three-dimensional, bounded-obstacle, acoustic scattering, J. Comput. Phys., 224 (2007), 1145–1169.
[15] C. Geuzaine, J.-F. Remacle, GMSH: a three-dimensional finite element mesh generator with built-in pre- and post-processing facilities, International Journal for Numerical Methods in Engineering, 79 (2009), 1309–1331.
[16] M. Grote and J. Keller, On nonreflecting boundary conditions, J. Comput. Phys., 122 (1995), 231–243.
[17] M. Grote and C. Kirsch, Dirichlet-to-Neumann boundary conditions for multiple scattering problems, J. Comput. Phys., 201 (2004), 630–650.
[18] T. Hagstrom, Radiation boundary conditions for the numerical simulation of waves, Acta Numerica, 8 (1999), 47–106.
[19] I. Harari and T. Hughes, Analysis of continuous formulations underlying the computation of time-harmonic acoustics in exterior domains, Comput. Methods Appl. Mech. Engrg., 97 (1992), 103–124.
[20] D. Jerison and C. Kenig, Unique continuation and absence of positive eigenvalues for Schrodinger operators, Ann. Math., 121 (1985), 463–488.
[21] X. Jiang, P. Li, and W. Zheng, Numerical solution of acoustic scattering by an adaptive DtN finite element method, Commun. Comput. Phys., 13 (2013), 1227–1244.
[22] X. Jiang, P. Li, J. Lv, and W. Zheng, An adaptive finite element method for the wave scattering with transparent boundary condition, J. Sci. Comput., 72 (2017), 936–956.
[23] J. Jin, The Finite Element Method in Electromagnetics, New York: Wiley, 1993.
[24] A. Kirsch and F. Hettlich, The Mathematical Theory of Time-Harmonic Maxwell’s Equations, Springer International Publishing, 2015.
[25] P. Li and X. Yuan, Convergence of an adaptive finite element DtN method for the elastic wave scattering by periodic structures, Comput. Methods Appl. Mech. Engrg., 360 (2020), 112722.
[26] P. Monk, Finite Element Methods for Maxwell’s Equations, Clarendon Press, Oxford, 2003.
[27] J.-C. Nédélec, Acoustic and Electromagnetic Equations Integral Representations for Harmonic Problems, Springer-Verlag, New York, 2001.
[28] A. H. Schatz, An observation concerning Ritz–Galerkin methods with indefinite bilinear forms, Math. Comp., 28 (1974), 959–962.
[29] E. Turkel and A. Yefet, Absorbing PML boundary layers for wave-like equations, Appl. Numer. Math., 27 (1998), 533–557.
[30] Z. Wang, G. Bao, J. Li, P. Li, and H. Wu, An adaptive finite element method for the diffraction grating problem with transparent boundary condition, SIAM J. Numer. Anal., 53 (2015), 1585–1607,

	$\displaystyle\\|u\\|^{2}_{H^{\frac{1}{2}}(\partial B_{R})}$	$\displaystyle=$	$\displaystyle\sum_{n=0}^{\infty}(1+n^{2})^{\frac{1}{2}}\sum\limits_{m=-n}^{n}\|\hat{u}_{n}^{m}(R)\|^{2}$
		$\displaystyle\lesssim$	$\displaystyle\sum_{n=0}^{\infty}\sum_{\|m\|\leq n}\int_{R^{\prime}}^{R}\big((1+n^{2})\|\hat{u}_{n}^{m}(r)\|^{2}+\|\hat{u}_{n}^{m^{\prime}}(r)\|^{2}\big){\rm d}r.$

$\displaystyle\|J_{1}\|$	$\displaystyle\lesssim$	$\displaystyle\sum\limits_{K\in\mathcal{M}_{h}}\bigg(h_{K}\\|(\Delta+\kappa^{2})u_{h}^{N}\\|_{L^{2}(K)}\\|\nabla\psi\\|_{L^{2}(\tilde{K})}+\sum\limits_{F\in\partial K}\frac{1}{2}h_{F}^{\frac{1}{2}}\\|J_{F}\\|_{L^{2}(F)}\\|\psi\\|_{H^{1}(\tilde{K}_{F})}\bigg)$
	$\displaystyle\lesssim$	$\displaystyle\sum_{K\in\mathcal{M}_{h}}\bigg[h_{K}\\|(\Delta+\kappa^{2})u_{h}^{N}\\|_{L^{2}(K)}+\bigg(\sum_{F\in\partial K}\frac{1}{2}h_{F}\\|J_{F}\\|_{L^{2}(F)}^{2}\bigg)^{\frac{1}{2}}\bigg]\\|\psi\\|_{H^{1}(\Omega)}$
	$\displaystyle\lesssim$	$\displaystyle\bigg(\sum_{K\in\mathcal{M}_{h}}\eta_{K}^{2}\bigg)^{2}\\|\psi\\|_{H^{1}(\Omega)}.$

$\displaystyle\|J_{2}\|=\|\langle(T-T_{N})u,\psi\rangle_{\partial B_{R}}\|$	$\displaystyle=$	$\displaystyle\Big\|R\sum_{n>N}\sum_{\|m\|\leq n}\Theta_{n}(\kappa R)\hat{u}_{n}^{m}(R)\bar{\hat{\psi}}_{n}^{m}(R)\Big\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R)\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\Big\|\Big(\frac{R^{\prime}}{R}\Big)^{N}\hat{u}_{n}^{m}(R^{\prime})\Big\|\|\bar{\hat{\psi}}_{nm}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|.$

$\displaystyle\|J_{2}\|$	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|\|\bar{\hat{\psi}}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\Big(\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|^{2}\Big)^{\frac{1}{2}}\Big(\sum_{\|m\|\leq n}\|\bar{\hat{\psi}}_{n}^{m}(R)\|^{2}\Big)^{\frac{1}{2}}$
	$\displaystyle\leq$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\Big(\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\hat{u}_{n}^{m}(R^{\prime})\|^{2}\Big)^{\frac{1}{2}}\Big(\sum_{n>N}(1+n(n+1))^{\frac{1}{2}}\sum_{\|m\|\leq n}\|\bar{\hat{\psi}}_{n}^{m}(R)\|^{2}\Big)^{\frac{1}{2}}$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\\|u\\|_{H^{\frac{1}{2}}(\partial B_{R^{\prime}})}\\|\psi\\|_{H^{\frac{1}{2}}(\partial B_{R})}$
	$\displaystyle\lesssim$	$\displaystyle\Big(\frac{R^{\prime}}{R}\Big)^{N}\\|u\\|_{H^{1}(\Omega)}\\|\psi\\|_{H^{\frac{1}{2}}(\partial B_{R})}.$

$\displaystyle\|\langle(T-T_{N})\xi,w\rangle_{\partial B_{R}}\|$	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|\Theta_{n}(\kappa R)\|\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\lesssim$	$\displaystyle\sum_{n>N}\sum_{\|m\|\leq n}\|n\|\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle=$	$\displaystyle\sum_{n>N}((1+n^{2})^{\frac{1}{2}}n^{3})^{-\frac{1}{2}}\sum_{\|m\|\leq n}(1+n^{2})^{\frac{1}{4}}n^{\frac{5}{2}}\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle N^{-2}\sum_{n>N}\sum_{\|m\|\leq n}(1+n^{2})^{\frac{1}{4}}n^{\frac{5}{2}}\|\hat{\xi}_{n}^{m}(R)\|\|\hat{w}_{n}^{m}(R)\|$
	$\displaystyle\leq$	$\displaystyle N^{-2}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}(1+n(n+1))^{\frac{1}{2}}\|\hat{\xi}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}$
	$\displaystyle=$	$\displaystyle N^{-2}\\|\xi\\|_{H^{\frac{1}{2}}(\partial B_{R})}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}$
	$\displaystyle\lesssim$	$\displaystyle N^{-2}\\|\xi\\|_{H^{1}(\Omega)}\bigg(\sum_{n>N}\sum_{\|m\|\leq n}n^{5}\|\hat{w}_{n}^{m}(R)\|^{2}\bigg)^{\frac{1}{2}}.$