1. Introduction
The electromagnetic scattering involving a perfect electric conducting (PEC) wedge with arbitrary apex angle and uniform coating by a double-negative metamaterial (DNG MTM) with losses has been tackled in [
1,
2] by considering an incident plane wave that is perpendicular to the edge. This two-dimensional problem has been studied at high frequencies by applying the Uniform Asymptotic Physical Optics (UAPO) method in the framework of the Uniform Geometrical Theory of Diffraction (UTD) [
3], which gives a very attractive model that combines the diffraction contribution and the Geometrical Optics (GO) ones for the scattering computation. In the UAPO method, the electric and magnetic surface currents have been used as equivalent sources on the external surfaces of the composite structure and incorporated in the scattering integral. Closed-form expressions have been obtained for parallel and perpendicular polarizations with respect to the edge of the structure by adopting useful approximations and integral representations, and performing asymptotic evaluations of the resulting integrals. The capability of the UAPO diffracted field to balance the GO field at the shadow boundaries has been demonstrated by means of numerical tests [
1,
2].
The solution for the two-dimensional plane wave diffraction in [
1,
2], as well as that which was proposed in [
4], has an application that is limited to schemes concerning normal incidences. Accordingly, the UAPO solution for the three-dimensional scenario of plane waves at a skew incidence is needed.
This article addresses the above need by accounting for the relevance of scattering problems involving coated metallic wedges, the metallic part of which can be replaced by a PEC structure with the same geometry at high frequencies, where the UAPO approach works. As it is well known, the electromagnetic scattering of dielectric and metallic structures can be modified by adding layers on their surfaces, and therefore such a practice results very attractive for civil and military applications, i.e., airplanes and ships can necessitate to lower their radar cross section by covering the metallic surfaces. Moreover, metamaterials have been recently proposed in Intelligent Reflecting Surfaces (IRS). This article refers to a particular structure that can be allocated in a smart wireless propagation scenario. The unusual electromagnetic peculiarities of DNG MTMs are related to the negative real parts of their permittivity and permeability, and offer unconventional solutions with respect to the standard materials, so opening new frontiers for research activities and applications. The reader that is fascinated by such materials can find interesting information and applications, as well as available manufacturing processes, in several books and articles (see [
5,
6,
7,
8,
9,
10,
11,
12] also for help in finding other references).
The UAPO approach has been already applied to solve three-dimensional scattering problems. A planar junction between a double-positive material and a DNG MTM with a PEC backing has been considered in [
13], whereas the scattering by a structure that is formed by a truncated DNG MTM sheet accommodated on a PEC support has been studied in [
14]. Some references [
15,
16,
17,
18,
19,
20] contain solutions for the plane wave diffraction problems involving dielectric wedges. In particular, the time domain counterpart of the UAPO solution in [
19] has been obtained in [
20] with reference to a penetrable wedge, which has an acute apex angle. The corresponding UAPO diffracted field has been always formulated in explicit and easy-to-use forms, and the effectiveness of the solutions has been demonstrated by adopting well-assessed numerical tools.
As pointed out in a previous statement, the UAPO method is applied in this article for solving the three-dimensional scattering problem of obliquely incident plane waves, which illuminate a PEC wedge with a DNG MTM planar coating. Since the method requires the evaluation of the field, which is reflected from the structure, the reflection coefficients need to be evaluated for each wedge surface that is involved in the scattering mechanism. The equivalent transmission line (ETL) model is here applied to this purpose in accordance with [
1,
2,
13,
20], and the results are used to formulate the PO currents on the surfaces. Opportune approximations, analytical manipulations and evaluations, and the application of an asymptotic procedure permit us to extract the UAPO diffracted contribution from the scattering integral part that is associated with each wedge surface if this last is lit by the incident plane wave. In other words, each surface contributes to the diffracted field only if it interacts with the propagating plane wave.
Section 2 is devoted to the formulation of the PO equivalent surface currents incorporating the reflection coefficients that are related to each wedge surface, whereas the UAPO solution for the diffracted field is proposed in
Section 3. The radio frequency (RF) unit of Comsol Multiphysics
® is applied in
Section 4 to test the efficacy of the proposed approach. The conclusions are presented in
Section 5.
2. PO Equivalent Surface Currents
A DNG MTM uniform layer with thickness d, permittivity and permeability , where are positive quantities and are related to the free space, covers an arbitrary-angled metallic wedge that is assumed to be PEC at the working frequencies and to have an apex angle .
According to
Figure 1, the external surface
of the wedge is lit by an incident plane wave with the electric field
, where
(
is the free-space propagation constant) and the vector
locates the observation point
P. If the PO approximation holds, such a field generates the following electric
and magnetic
equivalent surface currents on
and
:
In the above expressions,
, the symbols
and
refer to parallel and perpendicular polarizations, respectively,
is the unit vector normal to the involved surface,
, and
. The ETL model is applied to calculate the reflection coefficients
R:
in which
and
are the free-space ETL characteristic impedances, and the ETL input impedances
for a DNG MTM layer covering a PEC support can be evaluated according to [
1,
2,
13,
20].
Formulas (1) and (2) provide zero surface currents if the incident wave does not illuminate the corresponding surface.
3. UAPO Diffracted Field
The equivalent sources (1) and (2) are then assumed to radiate in the free space around the wedge, and the following integral formulation can be used to evaluate the scattered field
by accounting for the linearity of the integral operator:
where
and
account for the positions of
P and the source points at
. The term
is the 3 × 3 identity matrix and
is equal to 1 if the related surface is lit by the incident plane wave or its value is 0.
The incidence direction
is assumed (see
Figure 1). According to [
13,
14], the UAPO approach is able to extract the diffraction contributions from (4) and to provide the diffracted field
at
with
(
s symbolizes the distance from the diffraction point to
P and
varies from 0 to
on the arc of circumference with radius
and center on the external edge):
Accordingly, the matrix
of the diffraction coefficients is given by the superposition of the surface contributions and uses:
The function
of the UTD model [
3] is incorporated in (6) and (7). Furthermore,
4. Numerical Examples
Two sets of figures are shown in
Section 4 to prove the effectiveness of the UAPO-based method for evaluating the electromagnetic scattering from an arbitrary-angled PEC wedge that is covered by a DNG MTM uniform layer. One should note that no limitations exist for the values of permittivity and permeability. The first set (
Figure 2,
Figure 3,
Figure 4 and
Figure 5) is devoted to the ability of the UAPO diffracted field to compensate for the jumps of the GO field at the shadow boundaries, whereas the second set (
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10) includes comparisons with the Comsol Multiphysics
® data. Both sets refer to a composite wedge that is characterized by
(
n = 1.25),
,
, and
if
is the free-space wavelength. Moreover,
P moves on a circular arc with
and, according to the
α value,
ϕ varies from 0 to 225°.
The magnitudes of the
and
components of the GO field, the UAPO diffracted field, and the total field are reported in
Figure 2 and
Figure 3, respectively, where
. Both of the wedge surfaces are lit by the incident wave (
) so that two reflection boundaries exist at
, where the GO field fails and exhibits inacceptable jumps. As expected, the UAPO diffracted field assumes significant magnitude values in the angular regions around such directions, and it is able to balance the GO jumps, so offering a continuous total field from 0 to 225°. This ability is due to the behavior of the diffraction contribution at the shadow boundaries: the value on the left side is different from the value on the right side, so obtaining the needed discontinuity for the compensation.
Figure 4 and
Figure 5, which refer to
also confirm such an ability when one reflected field boundary
and one incident field boundary
exist since the incident wave impacts only one surface (
,
).
With reference to the figures of the second set, the wave propagation is considered to be orthogonal to the edge to reduce the needed computation effort when one is running the RF unit of Comsol Multiphysics®. The apex of the wedge is centered in a square propagation scenario with the side length of and the thickness of the perfectly matched layer at the box boundary is . In order to obtain reliable data, automatic adaptive meshes from to have been used.
Save for
Figure 8, the Figures from 6 to 10 contain: (a) the amplitude of the UAPO-based field component and the corresponding Comsol Multiphysics
® data; (b) a comparison between the related phases. One should note that the reported plots show large fluctuations, so providing hard test-beds.
Figure 6 and
Figure 7 are relevant to
and
, respectively, when
, and despite the fast variability of amplitude and phase of the field, the UAPO-based data fit very well those attained by Comsol Multiphysics
®. Moreover, the interested reader can surely appreciate
Figure 8a,b, where the UAPO diffracted field components are compared with the corresponding Comsol Multiphysics
® data that were obtained by extracting the GO field contributions. The agreements are very encouraging.
The last figures refer to
. Accordingly, one surface of the wedge is involved in the evaluation of the UAPO diffracted field as in
Figure 4 and
Figure 5. Good results are again achieved, but some of the lobes of the amplitude patterns are not correctly recovered and the phase values show misalignments in correspondence of the angular region centered on the incident field boundary. It is authors’ opinion that such a case suffers from the absence of the diffraction contribution from the internal PEC wedge.