1. Introduction
Consider the electromagnetic scattering of a timeharmonic plane wave by an open cavity, which is referred to as a bounded domain embedded in the ground with its opening aligned with the ground surface. The open cavity scattering problems have significant applications in industry and military. In computational and applied electromagnetics, one of the physical parameter of interests is the radar cross section (RCS), which measures the detectability of a target by a radar system. It is crucial to have a deliberate control in the form of enhancement or reduction of the RCS of a target in the stealth technology. The cavity RCS caused by jet engine inlet ducts or cavitybacked patch or slot antennas can dominate the total RCS of an aircraft or a device. It is indispensable to have a thorough understanding of the electromagnetic scattering characteristic of a target, particularly a cavity, in order to successfully implement any desired control of its RCS.
Due to the important applications, the open cavity scattering problems have received much attention by many researchers in both of the engineering and mathematics communities. The timeharmonic problems of cavitybacked apertures with penetrable material filling the cavity interior were introduced and studied initially by researchers in the engineering community [19, 27, 21]. The mathematical analysis for the wellposedness of the variational problems can be found in [1, 2, 3]
, where the nonlocal transparent boundary conditions, based on the Fourier transform, were proposed on the open aperture of the cavity. It has been realized that the phenomena of electromagnetic scattering by cavities not only have striking physics but also give rise to many interesting mathematical problems. As more people work on this subject, there has been a rapid development of the mathematical theory and computational methods for the open cavity scattering problems. The stability estimates with explicit dependence on the wavenumber were obtained in
[9, 10]. Various analytical and numerical methods have been proposed to solve the challenging large cavity problem [6, 11, 8, 30, 22]. The overfilled cavity problems, where the filling material inside the cavity may protrude into the space above the ground surface, were investigated in [14, 15, 25, 29], where the transparent boundary conditions, based on the Fourier series, were introduced on a semicircle enclosing the cavity and filling material. The multiple cavity scattering problem was examined in [24, 33], where the cavity is assumed to be composed of finitely many disjoint components. The mathematical analysis can be found in [7, 12] on the related scattering problems in a locally perturbed halfplane. We refer to the survey [23] and the references cited therein for a comprehensive account on the modeling, analysis, and computation of the open cavity scattering problems.There are two challenges for the open cavity scattering problems: the problems are formulated in unbounded domains; the solutions may have singularities due to possible nonsmooth surfaces and discontinuous media. In this paper, we present an adaptive finite element method with transparent boundary condition to overcome the difficulties.
The first issue concerns the domain truncation. The unbounded physical domain needs to be truncated into a bounded computational domain. An appropriate boundary condition is required on the artificial boundary of the truncated domain to avoid unwanted wave reflection. Such a boundary condition is known as a transparent boundary condition (TBC). There are two different TBCs for the open cavity scattering problems. For a regular open cavity, where the filling material is inside the cavity, the Fourier transform based TBC is imposed on the open aperture of the cavity; for an overfilled cavity, where the filling material appears to protrude out of the cavity through the open aperture into the space above the ground surface, the Fourier series based TBC is imposed on the semicircle enclosing the cavity and the protruding part. The latter is adopted in this work since it can be used to handle more general open cavities. We refer to the perfectly matched layer (PML) techniques [33, 34] and the method of boundary integral equations [5] as alternative approaches for dealing with the issue of the unbounded domains of the open cavity scattering problems.
Due to the existence of corners of cavities or the discontinuity of the dielectric coefficient for the filling material, the solutions have singularites that slow down the convergence of the finite element for uniform mesh refinements. The second issue can be resolved by using the a posteriori error estimate based adaptive finite element method. The a posteriori error estimates are computable quantities from numerical solutions. They measure the solution errors of discrete problems without requiring any a priori information of exact solutions. It is known that the meshes and the associated numerical complexity are quasioptimal for appropriately designed adaptive finite element methods.
The goal of this paper is to combine the adaptive finite element method and the transparent boundary conditions to solve the open cavity scattering problems in an optimal fashion. Specifically, we consider the scattering of a timeharmonic electromagnetic plane wave by an open cavity embedded in an infinite ground plane. Throughout, the medium is assumed to be constant in the direction. The ground plane and the cavity wall are assumed to be perfect electric conductors. The cavity is filled with a nonmagnetic and possibly inhomogeneous material, which may protrude out of the cavity to the upper halfspace in a finite extend. The infinite upper halfspace above the ground plane and the protruding part of the cavity is composed of a homogeneous medium. Two fundamental polarizations, transverse magnetic (TM) and transverse electric (TE), are studied. In this setting, the threedimensional Maxwell equations may be reduced to the two dimensional Helmholtz equation and generalized Helmholtz equation for TM and TE polarizations, respectively. Based on the DirichlettoNeumann (DtN) map for each polarization, a transparent boundary condition is imposed to reduce the scattering problem equivalently into a boundary value problem in a bounded domain. The nonlocal DtN operator is defined as an infinite Fourier series which needs to be truncated into a sum of finitely many terms in actual computation. The a posteriori error estimate is derived bewteen the solution of original scattering problem and the finite element solution of the discrete problem with the truncated DtN operator. The error estimate takes account of the finite element discretization error and the truncation error of the DtN operator. Using the asymptotic properties of the solution and DtN operator, we consider a dual problem for the error and show that the truncation error of the DtN operator decays exponentially respect to the truncation parameter, which implies that the truncation number does not need to be large. Numerical experiments are presented for both polarization cases to demonstrate the effectiveness of the proposed adaptive method. The related work can be found in [16, 17, 18, 31] on the adaptive finite element DtN method for solving other scattering problems in open domains.
The paper is organized as follows. Section 2 concerns the problem formulation. The threedimensional Maxwell equations are introduced and reduced into the twodimensional Helmholtz equation under the two fundamental modes: transverse magnetic (TM) polarization and transverse electric (TE) polarization. Sections 3 and 4 are devoted to the TM and TE polarizations, respectively. In each section, the variational problem and its finite element approximation are introduced; the a posteriori error analysis is given for the discrete problem with the truncated DtN operator; the adaptive finite element algorithm is presented. In Section 5, the stiff matrix is constructed for the the TBC part of the sesquilinear form. Section 6 describes the formulas of the backscatter radar cross section (RCS). Section 7 presents some numerical examples to illustrate the advantages of the proposed method. The paper is concluded with some general remarks and directions for future research in Section 8.
2. Problem formulation
Consider the electromagnetic scattering by an open cavity, which is a bounded domain embedded in the ground with its opening aligned with the ground surface. By assuming the time dependence , the electromagnetic wave propagation is governed by the timeharmonic Maxwell equations
(2.1) 
where is the electric field, is the magnetic field, is the magnetic flux density, is the electric flux density, is the electric current density, and is the angular frequency. For a linear medium, the constitutive relations, describing the macroscopic properties of the medium, are given by
(2.2) 
where is the magnetic permeability, is the electric permittivity, and is the electrical conductivity. Throughout, the medium is assumed to be nonmagnetic, i.e., the magnetic permeability is a constant everywhere, but the electric permittivity and the electrical conductivity are allowed to be spatial variable functions. Substituting (2.2) into (2.1) leads to a coupled system for the electric and magnetic fields
(2.3) 
Eliminating the magnetic field from (2.3), we may obtain the Maxwell system for the electric field
(2.4) 
where the wavenumber . Similarly, we may eliminate the electric field and obtain the Maxwell system for the magnetic field
(2.5) 
When the cavity has a constant cross section along the axis and the plane of incidence is in the plane, as a consequence, the electromagnetic fields are independent of the variable. The threedimensional Maxwell equations can be reduced to either the twodimensional Helmholtz equation or the twodimensional generalized Helmholtz equation.
Let be the cross section of the invariant cavity with a Lipschitz continuous boundary . Here is the cavity wall and is the open aperture of the cavity, which is aligned with the infinite ground plane . The cavity is filled with an inhomogeneous medium characterized by the dielectric permittivity , the magnetic permeability , and the electric conductivity . We point out that the inhomogeneous medium filling the cavity may protrude into the space above the ground plane, which is called an overfilled cavity. Let and be upper halfdiscs with radii and , where . Denote by and the upper semicircles. The radius can be taken to be sufficiently large such that the open exterior domain is filled with a homogeneous medium with constant permittivity and zero conductivity . Let be the bounded domain where our reduced boundary value problems are formulated. The problem geometry is shown in Figure 1.
Since the structure is invariant in the axis, we consider two fundamental polarizations for the electromagnetic fields: transverse magnetic (TM) polarization and transverse electric (TE) polarization. In TM case, the magnetic field is perpendicular to the plane of incidence and does not have the component in the axis; the electric field, being perpendicular to the magnetic field and lying in the plane, is invariant in the axis and takes the form , where is a scalar function. It is easy to verify from (2.4) that satisfies the Helmholtz equation
(2.6) 
In TE case, the electromagnetic fields are characterized by its electric field being perpendicular to the plane of incidence and contain no electric field component in the axis. The magnetic field, being perpendicular to the electric field and lying in the plane, is invariant in the axis and has the form , where is also a scalar function. It follows from (2.5) that satisfies the generalized Helmholtz equation
(2.7) 
When the ground plane and the cavity wall are assumed to be perfect conductors, the following perfectly electrically conducting (PEC) boundary condition can be imposed
(2.8) 
where
is the unit normal vector to
and . In TM polarization, the PEC boundary condition (2.8) reduces to the homogeneous Dirichlet boundary condition(2.9) 
In TE polarization, the PEC boundary condition (2.8) reduces to the homogeneous Neumann boundary condition
(2.10) 
3. TM polarization
In this section, we discuss the TM polarization and study its finite element approximation. The a posteriori analysis is carried out for both the finite element discretization error and the DtN operator truncation error. An adaptive finite element DtN method is presented for the truncated discrete problem.
3.1. Variational problem
In TM polarization, the nonzero component of the electric field satisfies the boundary value problem
(3.1) 
Since the problem is imposed in the open domain, a radiation condition is required to complete the formulation.
Consider the incidence of a plane wave
which is sent from the above to impinge the cavity. Here , is the incident angle, and is the wavenumber in the free space . It is easy to verify from (2.9) that the reflected wave is
By the Jacobi–Anger identity, the incident and reflected waves admit the following expansions:
(3.2) 
and
(3.3) 
where is the Bessel function of the first kind with order and with being the observation angle. Define the reference wave . It follows from (3.2)–(3.3) that
(3.4) 
The total field consists of the reference field and the scattered field , i.e.,
(3.5) 
where the scattered field is required to satisfy the Sommerfeld radiation condition
Let . For any , it has the Fourier series expansion
Define the trace function space where the norm is given by
It is clear that the dual space of is with respect to the scalar product in given by
As discussed in [23], a DtN operator is introduced on :
(3.6) 
where is the Hankel function of the first kind with order . It is shown in [29, Lemma 3.1] that is continuous. The TBC can be imposed for the total field as follows:
where . Substituting (3.4) into (3.6) and applying the Wronskian identity, we obtain explicitly
The original cavity scattering problem (3.1) can be reduced equivalently into the boundary value problem
which has the variational formulation: find such that
(3.7) 
Here the sesquilinear form is defined as
Theorem 3.1.
The variational problem (3.7) has a unique solution , which satisfies the estimate
Hereafter, the notation stands for , where is a positive constant whose value is not required but should be clear from the context.
3.2. Finite element approximation
Let be a regular triangulation of , where denotes the maximum diameter of all the elements in . To avoid being distracted from the main focus of the a posteriori error analysis, we assume for simplicity that and are polygonal to keep from using the isoparametric finite element space and deriving the approximation error of the boundaries and . Thus any edge is a subset of if it has two boundary vertices.
Let be a conforming finite element space, i.e.,
In practice, the DtN operator (3.6) needs to be truncated into a sum of finitely many terms
(3.8) 
Taking account of the DtN operator truncation, we obtain the finite element approximation to the variational problem (3.7): find such that
(3.9) 
where the sesquilinear form is
For sufficiently small and sufficiently large , the discrete infsup condition of the sesquilinear form can be established by an argument of Schatz [28]. It follows from the general theory in [4] that the truncated variational problem (3.9) admits a unique solution. Since our focus is the a posteriori error estimate and the associated adaptive algorithm, we assume that the discrete problem (3.9) has a unique solution .
3.3. A posteriori error analysis
For any triangular element , denoted by its diameter. Let denote the set of all the edges of . For any edge , denote by its length. For any interior edge , which is the common side of triangular elements , we define the jump residual across as
where is the unit outward normal vector on the boundary of . For any boundary edge , the jump residual is defined as
For any triangle , denote by the local error estimator as follows:
where is the Helmholtz operator defined by .
Let , where and are the solutions of the variational problems (3.7) and (3.9), respectively. Introduce a dual problem: find such that
(3.10) 
It is easy to check that is the solution of the following boundary value problem:
where is the adjoint operator of and is given by
The following three lemmas are proved in [17], where the first lemma concerns the wellposedness of the dual problem, the second lemma gives the trace result in , and the third lemma shows the error representation formulas.
Lemma 3.2.
The dual problem (3.10) has a unique solution , which satisfies the estimate
Lemma 3.3.
For any , the following estimates hold:
Lemma 3.4.
The following result concerns the truncation error of the DtN operator and plays an important role in the a posteriori error estimate.
Lemma 3.5.
Let be the solution to (3.7) and be any function in . For sufficiently large , the following estimate holds:
Proof.
By (3.5), we have , where is the reference field and is the scattered field satisfying the Sommerfeld radiation condition. For sufficiently large , it is shown in [17, Lemma 4] that
A straightforward calculation yields
By [32], for sufficiently large , we have
which give
For , it is easy to verify
Hence
Combining the above estimates, we obtain
Remark 3.6.
We notice that the result and proof of Lemma 3.5 is different from those for the scattering problems in periodic structures [18, 31, 26]. For the latter problems, the DtN operators are defined on a straight line or plane surface and have only finitely many terms when acting on the incident fields. For our case, the DtN operator is defined on a semicircle and is still an infinite series when acting on the reference field, which results in an extra term in the estimate given in Lemma 3.5.
Lemma 3.7.
Let be the solution of the dual problem (3.10). Then the following estimate holds:
Proof.
Since on , it admits the Fourier series expansion in terms of the sin functions
where are the Fourier coefficients. Following the same proof as that in [17, Lemma 5], we may show the desired result. ∎
3.4. Adaptive FEM algorithm
It is shown in Theorem 3.8 that the a posteriori error consists of two parts: the finite element discretization error and the DtN operator truncation error , where
(3.11) 
In the implementation, based on (3.11), the parameters , and can be chosen appropriately such that the finite element discretization error is not contaminated by the truncation error, i.e., is required to be small compared with , for instance, . Table 1 shows the algorithm of the adaptive finite element DtN method for solving the open cavity scattering problem in the TM polarization.

Given the tolerance and the parameter .

Fix the computational domain by choosing .

Choose and such that .

Construct an initial triangulation over and compute error estimators.

While do

refine mesh according to the strategy

denote refined mesh still by , solve the discrete problem (3.9) on the new mesh ,

compute the corresponding error estimators.

End while.
4. TE polarization
In this section, we consider the TE polarization. Since the discussions are similar to the TM polarization, we briefly present the parallel results without providing the details.
In TE polarization, the total field satisfies the boundary value problem of the generalized Helmholtz equation
(4.1) 
Consider the same plane incident wave . Due to the homogeneous Neumann boundary condition on , the reflected field is
By the Jacobi–Anger identity, the reference wave admits the following expansion:
(4.2) 
Once again, the total field is assumed to be composed of the reference field
Comments
There are no comments yet.