Time-dependent electromagnetic scattering from thin layers

We consider the inverse scattering problem of determining the shape of an impenetrable two-dimensional obstacle with a generalized impedance boundary condition from near-field measurements. We propose the reverse time migration method where the imaging functional is defined as the imaginary part of the cross-correlation. Our analysis shows that, on the one hand, the imaging functional has contrast at the boundary of the scatterer and decays outside the scatterer; On the other hand, the imaging functional is always positive thus may have better stability properties. These imply that the approach has powerful imaging quality which is confirmed in our numerical experiments.

A high-order nodal discontinuous Galerkin method for a linearized fractional Cahn–Hilliard equation

The Laguerre transform is applied to the convolution product of functions of a real argument (over the time axis) with values in Hilbert spaces. The main results have been obtained by establishing a relationship between the Laguerre and Laplace transforms over the time variable with respect to the elements of Lebesgue weight spaces. This relationship is built using a special generating function. The obtained dependence makes it possible to extend the known properties of the Laplace transform to the case of the Laguerre transform. In particular, this approach concerns the transform of a convolution of functions.
 The Laguerre transform is determined by a system of Laguerre functions, which forms an orthonormal basis in the weighted Lebesgue space. The inverse Laguerre transform is constructed as a Laguerre series. It is proven that the direct and the inverse Laguerre transforms are mutually inverse operators that implement an isomorphism of square-integrable functions and infinite squares-summable sequences.
 The concept of a q-convolution in spaces of sequences is introduced as a discrete analogue of the convolution products of functions. Sufficient conditions for the existence of convolutions in the weighted Lebesgue spaces and in the corresponding spaces of sequences are investigated. For this purpose, analogues of Young’s inequality for such spaces are proven. The obtained results can be used to construct solutions of evolutionary problems and time-dependent boundary integral equations.

In this paper we consider the electromagnetic scattering problem by an obstacle characterised by a Generalized Impedance Boundary Condition in the harmonic regime. These boundary conditions are well known to provide accurate models for thin layers or imperfectly conducting bodies. We give two different formulations of the scattering problem and we provide some general assumptions on the boundary condition under which the scattering problem has at most one solution. We also prove that it is well-posed for three different boundary conditions which involve second order surface differential operators under weak sign assumptions on the coefficients defining the surface operators.

The acoustic wave equation on the whole three-dimensional space is considered with initial data and inhomogeneity having support in a bounded domain, which need not be convex. We propose and study a numerical method that approximates the solution using computations only in the interior domain and on its boundary. The transmission conditions between the interior and exterior domain are imposed by a time-dependent boundary integral equation coupled to the wave equation in the interior domain. We give a full discretization by finite elements and leapfrog time-stepping in the interior, and by boundary elements and convolution quadrature on the boundary. The direct coupling becomes stable on adding a stabilization term on the boundary. The derivation of stability estimates is based on a strong positivity property of the Calderon boundary operators for the Helmholtz and wave equations and uses energy estimates both in time and frequency domain. The stability estimates together with bounds of the consistency error yield optimal-order error bounds of the full discretization.

A semi-analytical finite element numerical mode-matching (FNMM) method is developed to simulate microscopic structures of metasurfaces on the interfaces of the 3-D multilayered structure. Similar to the recently developed spectral element NMM method, the FNMM method reduces the original 3-D problem into a series of 2-D waveguide eigenvalue problems plus a 1-D layered medium problem, which can be analytically described by a recursion procedure. The 2-D eigenvalue problems in the horizontal directions are solved numerically by the mixed finite element method (MFEM) to remove the spurious modes in the conventional FEM. Therefore, it can effectively and accurately simulate electromagnetic wave interactions with the detail metasurfaces, especially for thick layers. Metasurfaces embedded on the interfaces of 3-D structures are treated by the generalized impedance boundary conditions, thus enabling their rapid yet accurate computation. This method is applied to characterize metasurfaces consisting of microscopic patterns with complex geometric construction with details much smaller than the wavelength. Numerical experiments indicate that the FNMM method is highly efficient and accurate for the microscopic metasurface models.

Printed circuits in bounded media encompass a wide range of practical structures such as discontinuities in waveguides, planar circuits embedded in shielded multilayered media or even two-dimensional printed periodic structures. The Electromagnetic (EM) modeling of printed circuits in layered bounded media is performed via an Integral Equation (IE) technique. Green's functions (GFs) are specially defined to satisfy both the Boundary Conditions (BCs) imposed by the layered media and by the transverse boundary enclosing the entire structure. Finally, a system of IEs on the equivalent sources can be solved numerically by means of the Method of Moments (MoM). Each of the problems enumerated above has already been solved by other authors using IE-MoM techniques. Nevertheless, our formulation introduces a unified approach applicable to all the aforementioned problems. Due to the symmetry presented by a bounded layered media, the GF problem can be reduced into a two-dimensional transverse boundary problem and a one-dimensional transmission line problem in the normal direction. Both problems can be treated independently. This thesis proposes and fully develops an efficient technique that encompasses different laterally bounded multilayered problems with a seamless transition between them. The method is based on a modal representation of the transverse boundary problem and on the expansion of the equivalent surface currents by zero-curl & constant-charge Basis Functions (BFs). It offers a unified and versatile approach that, on one hand eliminates redundancy in the formulation and on the other hand simplifies each particular problem to the evaluation of constant coefficients or basic line integrals. Analytical solutions can be found for the combination of linear subsectional basis functions in rectangular and circular Perfect Electric Conductor (PEC) boundaries as well as for periodic lattices. This thesis then solves the problem of transmission line model in the longitudinal direction by proposing an efficient algorithm that guarantees numerical stability under a variety of known critical conditions where other already known formulations fail. In addition, it introduces alternate equivalent expressions of this formulation that allow new interpretations of the problem. Due to its practical interest, the method is applied for the EM modeling of multilayered boxed printed circuits. This motivated the implementation of a dedicated software tool for the efficient analysis of these topologies including losses. Extensive numerical experiments have been carried out to assess the validity of the aforementioned theory and some properties of test-structures (losses, mesh, etc).

Parallel finite difference schemes for heat equation based upon overlapping domain decomposition

We study time-harmonic scattering in [Formula: see text] ([Formula: see text]) by a planar screen (a “crack” in the context of linear elasticity), assumed to be a non-empty bounded relatively open subset [Formula: see text] of the hyperplane [Formula: see text], on which impedance (Robin) boundary conditions are imposed. In contrast to previous studies, [Formula: see text] can have arbitrarily rough (possibly fractal) boundary. To obtain well-posedness for such [Formula: see text] we show how the standard impedance boundary value problem and its associated system of boundary integral equations must be supplemented with additional solution regularity conditions, which hold automatically when [Formula: see text] is smooth. We show that the associated system of boundary integral operators is compactly perturbed coercive in an appropriate function space setting, strengthening previous results. This permits the use of Mosco convergence to prove convergence of boundary element approximations on smoother “prefractal” screens to the limiting solution on a fractal screen. We present accompanying numerical results, validating our theoretical convergence results, for three-dimensional scattering by a Koch snowflake and a square snowflake.

The local projection stabilization allows us to circumvent the Babuška–Brezzi condition and to use equal order interpolation for discretizing the Stokes problem. The projection is usually done in a two-level approach by projecting the pressure gradient onto a discontinuous finite element space living on a patch of elements. We propose a new local projection stabilization method based on (possibly) enriched finite element spaces and discontinuous projection spaces defined on the same mesh. Optimal order of convergence is shown for pairs of approximation and projection spaces satisfying a certain inf-sup condition. Examples are enriched simplicial finite elements and standard quadrilateral/hexahedral elements. The new approach overcomes the problem of an increasing discretization stencil and, thus, is simple to implement in existing computer codes. Numerical tests confirm the theoretical convergence results which are robust with respect to the user-chosen stabilization parameter.

We present full discretization of the telegraph equation with nonlocal coeffecient using Rothe-nite element method. For solving the equation numerically we use the Newton Raphson method, but the nonlocal term causes diffeculties because the Jacobien matrix is full. To remedy these diffeculties we apply the technique used by Sudhakar [4]. The optimal a priori error estimates for both semi discrete and fully discrete schemes are derived in V and H1 and a numerical experiment is described to support our theoretical result.

This article introduces the hybridizable discontinuous Galerkin (HDG) approach to numerically approximate the solution of a linear hyperbolic integro-differential equation. A priori error estimates for semi-discrete and fully discrete schemes are developed. It is shown that the optimal order of convergence is achieved for both scalar and flux variables. To achieve that, an intermediate projection is introduced for the semi-discrete error analysis, and it also shows that this projection achieves convergence of order $h^{k+3/2}$ for $k ≥ 1.$ Next, superconvergence is achieved for the scalar variable using element-by-element post-processing. For the fully discrete error analysis, the central difference scheme and the mid-point rule approximate the derivative and the integral term, respectively. Hence, the second order of convergence is achieved in the temporal direction. Finally, numerical experiments have been performed to validate the theory developed in this article.

In this paper, we present a new accelerating procedure in order to speed up the convergence of Newton-type methods. In particular, we derive iterations with a high and optimal order of convergence. This technique can be applied to any iteration with a low order of convergence.
 As expected, the convergence of the proposed methods was remarkably fast. The effectiveness of this technique is illustrated by numerical experiments.

We consider the stability of an efficient Crank–Nicolson–Adams–Bashforth method in time, finite element in space, discretization of the Leray‐α model. We prove finite‐time stability of the scheme in L2, H1, and H2, as well as the long‐time L‐stability of the scheme under a Courant‐Freidrichs‐Lewy (CFL)‐type condition. Numerical experiments are given that are in agreement with the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1155–1183, 2016

The scattering by the perfectly electric conducting (PEC) half-plane and PEC zero thickness disk placed on parallel planes is considered. The fields are represented in the spectral domain, i.e. in the domain of Fourier transform. The operator equations with respect to the Fourier amplitudes of the scattered field are obtained. The kernel functions of these equations contain poles. After regularization procedure, which is connected with the elimination of the poles, operator equations are converted to the system of singular integral equations. The convergence of the solution is based on the corresponding theorems. The scattered field consists of the plane wave, reflected by the infinite part of the half-plane, cylindrical waves, which appear as a result of scattering by the edge of the half-plane, and spherical waves, which appear as a result of scattering by the disk and multiple re-scattering by the disk-half-plane. The total near-field distribution and far-field patterns of cylindrical waves are presented.