A new family of first-order boundary conditions for the Maxwell system: derivation,well-posedness and long-time behavior

This work deals with the time-dependent Maxwell system in the case of TE-polarized electromagnetic waves, when associated with a family of first-order local boundary conditions. The boundary conditions are derived by using a micro-diagonalization method, actuated by the standard one of M.E. Taylor and involving pseudodifferential technics. The conditions differ from an arbitrary function and any of them leads to a well-posed mixed problem that is described by a continuous semi-group. The arbitrary function can be seen as a parameter and an asymptotic analysis in time shows that it can be chosen so that the resulting boundary condition is absorbing: the system is related to an energy functional that converges towards zero as time tends to infinity. By involving an invariant space for the Maxwell system, the limit state can be explicitly written as a solution to a boundary-value problem depending on the initial data. The long time behavior of the solution is then completely analyzed.     

A principle of images for absorbing boundary conditions

When one uses high-order finite difference schemes for the wave equation, for instance fourth order schemes, the treatment of boundary conditions poses a real difficulty since one needs several additional equations (for the nodes close to the boundary), while one single scalar boundary condition is available. In the case of perfectly reflecting boundary conditions, namely the homogeneous Neumann or Dirichlet conditions, this difficulty can be overcomed by the use of the well-known image principle, which permits the extension of the equation outside of the domain of calculation by an appropriate symmetrization of the data. We propose in this article a generalization of this principle to the absorbing boundary conditions. Through a symmetrization process, we are led to introduce a damped wave equation with a damping term supported by the boundary. The treatment of the boundary condition is then replaced by the approximation of this new damped wave equation in the whole space. The theoretical justification of our approach is based on new energy estimates for the wave equation (when high-order absorbing boundary conditions are used), and constitutes an alternative to the use of the well-known Kreiss criterion to prove the stability of the associated initial boundary value problems. © 1994 John Wiley & Sons, Inc.     

On Local Nonreflecting Boundary Conditions for Time Dependent Wave Propagation

The simulation of wave phenomena in unbounded domains generally requires an artificial boundary to truncate the unbounded exterior and limit the computation to a finite region. At the artificial boundary a boundary condition is then needed, which allows the propagating waves to exit the computational domain without spurious reflection. In 1977, Engquist and Majda proposed the first hierarchy of absorbing boundary conditions, which allows a systematic reduction of spurious reflection without moving the artificial boundary farther away from the scatterer. Their pioneering work, which initiated an entire research area, is reviewed here from a modern perspective. Recent developments such as high-order local conditions and their extension to multiple scattering are also presented. Finally, the accuracy of high-order local conditions is demonstrated through numerical experiments.     

Remarks on the Huygens absorbing boundary conditions for electromagnetics

The Huygens absorbing boundary conditions (ABCs) are promising new implementations of operator ABCs. They have certain advantageous features which are lacked in other operator ABCs. Under certain conditions and with the Huygens ABCs, the transmitted wave depends solely upon the second derivative with respect to time or upon the double integral of the incident wave. For such cases and for problems with a Dirichlet boundary condition, the overall reflection is not unique. Two new examples of the Huygens ABCs are given for such cases. For each example and with a FDTD scheme the newly derived reflection is less than that which has been studied by Bérenger [J.-P. Bérenger, On the Huygens absorbing boundary conditions for electromagnetics, J. Comput. Phys. 226 (2007) 354–378].     

On surface radiation conditions for high-frequency wave scattering

A new approximation of the logarithmic derivative of the Hankel function is derived and applied to high-frequency wave scattering. We re-derive the on surface radiation condition (OSRC) approximations that are well suited for a Dirichlet boundary in acoustics. These correspond to the Engquist–Majda absorbing boundary conditions. Inverse OSRC approximations are derived and they are used for Neumann boundary conditions. We obtain an implicit OSRC condition, where we need to solve a tridiagonal system. The OSRC approximations are well suited for moderate wave numbers. The approximation of the logarithmic derivative is also used for deriving a generalized physical optics approximation, both for Dirichlet and Neumann boundary conditions. We have obtained similar approximations in electromagnetics, for a perfect electric conductor. Numerical computations are done for different objects in 2D and 3D and for different wave numbers. The improvement over the standard physical optics is verified.     

On the Maxwell system under impedance boundary conditions with memory

We consider an initial boundary value problem for the system of the Maxwell equations in a bounded domain with smooth boundary on a finite time interval with new boundary conditions with memory. In appropriate function spaces, we define and study the nonselfadjoint operator that is generated by the Maxwell operator under a boundary condition with memory. Using the operator method, we prove an existence and uniqueness theorem for a solution to the initial boundary value problem.     

Algebraic Dirichlet-to-Neumann mapping for linear elasticity problems with extreme contrasts in the coefficients

The convergence of iterative based domain decomposition methods is linked with the absorbing boundary conditions defined on the interface between the sub-domains. For linear elasticity problems, the optimal absorbing boundary conditions are associated with non-local Dirichlet-to-Neumann maps. Most of the methods to approximate these non-local maps are based on a continuous analysis. In this paper, an original algebraic technique based on the computation of local Dirichlet-to-Neumann maps is investigated. Numerical experiments are presented for linear elasticity problems with extreme contrasts in the coefficients.     

Well-posedness and exponential stability of Maxwell-like systems coupled with strongly absorbing layers

This work deals with a family of PDEs which are derived from the Maxwell system set into an anisotropic medium. The system consists of four coupled linear equations. The first two correspond to the Maxwell system perturbed by zero-order operators which are represented by diagonal tensors with compact support. The last two equations are ODEs. Each system is written in a bounded domain and its boundary is modelled by a Silver–Müller condition. Firstly we establish that each problem is well-posed, providing the tensors are positive with bounded terms in Ω. Secondly we address the question of the long-time stability of each model. We prove that there exists a functional of energy which is exponentially decreasing if the domain is strictly star-shaped.     

Absorbing Interface Conditions for the Simulation of Wave Propagation on Non-Uniform Meshes

We proposed absorbing interface conditions for the simulation of linear wave propagation on non-uniform meshes. Based on the superposition principle of second-order linear wave equations, we decompose the interface condition problem into two subproblems around the interface: for the first one the conventional artificial absorbing boundary conditions is applied, while for the second one, the local analytic solutions can be derived. The proposed interface conditions permit a two-way transmission of low-frequency waves across mesh interfaces which can be supported by both coarse and fine meshes, and perform a one-way absorption of high-frequency waves which can only be supported by fine meshes when they travel from fine mesh regions to coarse ones. Numerical examples are presented to illustrate the efficiency of the proposed absorbing interface conditions.     

Approximation of pseudodifferential operators in absorbing boundary conditions for hyperbolic equations

Summary Engquist and Majda [3] proposed a pseudodifferential operator as asymptotically valid absorbing boundary condition for hyperbolic equations. (In the case of the wave equation this boundary condition is valid at all frequencies.) Here, least-squares approximation of the symbol of the pseudodifferential operator is proposed to obtain differential operators as boundary conditions. It is shown that for the wave equation this approach leads to Kreiss well-posed initial boundary value problems and that the expectation of the reflected energy is lower than in the case of Taylor- and Padé-approximations [3, 4]. Numerical examples indicate that this method works even more effectively for hyperbolic systems. The least-squares approach may be used to generate the boundary conditions automatically.     

带吸收边界条件的声波方程显式差分格式的稳定性分析

１．引言 在进行无界或半无界区域上各种波动方程的数值求解时,需引进入工边界以限制计算范围,在这些边界上应加相应的人工边界条件．这种边界条件应保证所求得的有界区域上的解很好地逼近原来无界区域上的解．对波动方程来说,就是在边界上人工反射应尽可能地小,使之对区域内部解的影响在允许的误差范围以内.因而它们被称为无反射边界条件或吸收边界条件．这种边界条件还应保证所形成的有界区域上的微分方程定解问题是适定的．这也是各种数值方法稳定的必要条件。 近二十多年来,已发展了声波方程的各种类型的吸收边界条件,其中以Ｃｌ…     

On the Fredholm property of the electric field equation in the vector diffraction problem for a partially screened solid

We consider a vector problem of diffraction of an electromagnetic wave on a partially screened anisotropic inhomogeneous dielectric body. The boundary conditions and the matching conditions are posed on the boundary of the inhomogeneity domain, and under passage through it, the medium parameters have jump changes. A boundary value problem for the system of Maxwell equations in unbounded space is studied in a semiclassical statement and is reduced to a system of integro-differential equations on the body domain and the screen surfaces. We show that the quadratic form of the problem operator is coercive and the operator itself is Fredholm with zero index.     

Generalized impedance boundary conditions for the maxwell equations as singular perturbations problems

In this paper the Maxwell equations in an exterior domain with generalaized impedance boundary conditions of Engquist-Nédélec are considered. The particular form of the assumed boundary conditions can be considered to be a singular perturbation of the Dirichlet boundary conditions. The convergence of the solution of the Maxwell equations with these generalized impedance boundary conditions to that of the corresponding Dirichlet problem is proven. The proof uses a new integral equations method combined with results dealing with singular perturbation problems of a class of pseudo-differential operators.     

A fourth-order energy for the three-dimensional wave equation with second-order absorbing boundary conditions

We introduce a fourth-order energy for the three-dimensional wave equation in rectangular domains with second-order absorbing boundary conditions. A decay rate of the energy in the domain with respect to time is estimated in terms of the boundary integral. Absorbing boundary conditions considered in this paper include the treatment of the case of annihilation of obliquely incident waves arriving at the boundary. Numerical results are also shown.     

声波方程吸收边界条件的稳定性分析

引言对于无界区域中波动现象的数值模拟，必需引进人工边界将计算限制在一个有界区域上．为了确定解，需要在人工边界上加适当的边界条件．对于声波和弹性波方程，这样的一组人工边界条件，也叫吸收边界条件，在［1，2］中被系统地构造出来．对于声波方程，这些吸收边界条件恰好是单程波方程的近似．如山中所指出，减少边界反射，便于在计算中应用和稳定性是构造吸收边界条件的三点关键．Ellgqllist和Maid。用模态分析方法15］证明，带有［IJ中构造的吸收边界条件的波动方程初边值问题是适定的，并且估计了人工边界所产生的误差．对于更广…     

Persistency of wellposedness of Ventcel’s boundary value problem under shape deformations

Ventcel boundary conditions are second order differential conditions that appear in asymptotic models. Like Robin boundary conditions, they lead to well-posed variational problems under a sign condition of the coefficient. This is achieved when physical situations are considered. Nevertheless, situations where this condition is violated appeared in several recent works where absorbing boundary conditions or equivalent boundary conditions on rough surfaces are sought for numerical purposes. The well-posedness of such problems was recently investigated: up to a countable set of parameters, existence and uniqueness of the solution for the Ventcel boundary value problem holds without the sign condition. However, the values to be avoided depend on the domain where the boundary value problem is set. In this work, we address the question of the persistency of the solvability of the boundary value problem under domain deformation.     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

基于透射边界条件的高阶离散型角点条件

对波动方程的数值模拟中，在有限区域建立吸收边界条件，其中对区域角点的处理是一个很重要的问题．随着吸收边界条件阶数的提高，与之匹配的角点条件也越难建立．MTF是一种离散型吸收边界条件．在此，对于二维问题，基于MTF建立离散型高阶角点条件，对计算区域角点处理时，在区域对角线方向上建立N阶MTF公式．问题也可推广到三维．数值结果证实了我们的猜测．     

The integral equation method in problems of diffraction by a finite impedance section of a medium interface

A 3D problem of reflection of a plane electromagnetic wave by a local impedance section of a wavy surface is considered. The boundary value problem for the system of Maxwell’s equations in a region with an irregular boundary is reduced to solution of systems of hypersingular integral equations. A numerical algorithm is proposed for solution of these systems. Results of numerical computations are presented.