Next Article in Journal
MD-TransUNet: An Image Segmentation Network for Car Front Face Design
Previous Article in Journal
Combining CBAM and Iterative Shrinkage-Thresholding Algorithm for Compressive Sensing of Bird Images
Previous Article in Special Issue
Controlling the Focusing Ability of the Parabolic Graphene Reflector in Free Space at Microwave Frequencies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectral Analysis of Electromagnetic Diffraction Phenomena in Angular Regions Filled by Arbitrary Linear Media

by
Vito G. Daniele
and
Guido Lombardi
*,†
The Department of Electronics and Telecommunications (DET), Politecnico di Torino, 10129 Torino, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2024, 14(19), 8685; https://doi.org/10.3390/app14198685
Submission received: 21 July 2024 / Revised: 19 August 2024 / Accepted: 27 August 2024 / Published: 26 September 2024

Abstract

:
A general theory for solving electromagnetic diffraction problems with impenetrable/penetrable wedges immersed in/made of an arbitrary linear (bianistropic) medium is presented. This novel and general spectral theory handles complex scattering problems by using transverse equations for layered planar and angular structures, the characteristic Green function procedure, the Wiener–Hopf technique, and a new methodology for solving GWHEs. The technique has been proven effective for analyzing problems involving wedges immersed in isotropic media; in this study, we extend the theory to more general cases while providing all necessary mathematical tools and corresponding validations. We obtain generalized Wiener–Hopf equations (GWHEs) from spectral functional equations in angular regions filled by arbitrary linear media. The equations can be interpreted with a network formalism for a systematic view. We recall that spectral methods (such as the Sommerfeld–Malyuzhinets (SM) method, the Kontorovich–Lebedev (KL) transform method, and the Wiener–Hopf (WH) method) are well-consolidated, fundamental, and effective tools for the correct and precise analysis of electromagnetic diffraction problems constituted by abrupt discontinuities immersed in media with one propagation constant, although they are not immediately applicable to multiple-propagation-constant problems. To the best of our knowledge, the proposed mathematical technique is the first extension of spectral analysis to electromagnetic problems in the presence of angular regions filled by complex arbitrary linear media, thereby providing novel mathematical tools. Validation through fundamental examples is proposed.

1. Introduction

The theory of wave diffraction constitutes one of the fundamental problems in mathematical physics. Apart from its direct relevance to engineering and physics, this subject gives rise to significant methodologies in applied mathematics.
Spectral methods play a crucial role in the study of electromagnetic diffraction. Notably, the Sommerfeld–Malyuzhinets (SM) method, the Kontorovich–Lebedev (KL) transform method, and the Wiener–Hopf (WH) method are fundamental and complementary in studying diffraction problems in the presence of sharp discontinuities. These methods have been extensively and effectively applied for studying wedge diffraction in isotropic regions; see [1,2,3,4,5,6] for SM, Refs. [7,8,9,10] for KL, and Refs. [11,12] for WH, along with the references therein starting from the early 2000s literature produced by the authors of this study. Moreover, using synergy among the three methods (WH, SM, and KL), the authors obtained a complete network representation of the angular region in the presence of isotropic media [13], which helped to build a systematic methodology for analysis.
The main advantage of the aforementioned techniques (SM and KL) is also a limitation, i.e., the utilization of the spectral complex angular plane derived from the Sommerfeld integral theory [14], which, since the early 2000s, the authors of this study have also effectively used in the WH framework for Fredholm factorization [11,12] and in the definition of rotating waves in isotropic angular regions [11,12]. The definition of this complex plane is intricately connected to the physics of the problem, as it specifically requires spectral transformations associated with the propagation constant. Consequently, this methodology is applicable to problems involving a single propagation constant, such as isotropic media in electromagnetic fields, as well as other specific problem configurations with decoupling properties in propagation modalities. Different attempts have been developed to extend the spectral analysis to diffraction problems in more complex media, such as gyrotropic media and/or uniaxial media. For example, we recall the analysis of scattering by a perfect electrically conducting (PEC) half-plane immersed in such anisotropic media; see [15,16,17,18,19,20,21,22,23,24]. However, to the best of our knowledge no spectral method for scattering problems by wedges in arbitrary linear media (i.e., bianisotropic media [25,26,27]) characterized by multiple propagation constants has been developed to date. One of the most important results obtained in the presence of anisotropic media is the exact solution obtained by Felsen in the case of scattering by a PEC wedge immersed in a uniaxial medium illuminated by plane waves at normal incidence [17,18]. However, the method used for this problem is substantially that of the separation of variables after transformations in the physical domain, which does not present the powerful characteristics of spectral methods such as asymptotic evaluation of fields and physical interpretation of field components in terms of structural and source spectral singularities. Other important works have examined the behavior of the field near the edge of a wedge immersed in a complex medium [28] and diffraction by a wedge immersed in the special case of an isotropic chiral medium with the SM method [29].
Given our experience in spectral analysis of complex electromagnetic scattering problems in isotropic media [11,12,30,31,32], and with the help of the theory proposed in [33] for the analysis of structures embedded in layered media, in this work we develop a new theory in the spectral domain with proper mathematical tools that allow for the representation of scattering problems immersed in arbitrary linear media of an angular shape. In particular, these new formulations are in the spectral domain (Laplace domain) without introducing angular complex planes, and as such are not limited to one-propagation-constant problems. In [34], we developed the general theory in abstract form to model angular regions filled by arbitrary linear media and we reported its implementation only for isotropic media.
In the present work, we propose a complete theoretical package for solving diffraction problems with impenetrable wedges immersed in an arbitrary linear medium; in addition, this package is extendable to multiple penetrable angular regions. The proposed method exploits the combination and the extension of powerful mathematical tools developed in different contexts. The first tool is the Bresler–Marcuvitz (BM) transverse equation theory for layered media [33,35], the second is the characteristic Green function procedure [36,37], the third one is the Wiener–Hopf technique [33,38] in its generalized form [11,12], and the fourth one (which is a completely novel contribution) is the direct application of Fredholm factorization to generalized Wiener–Hopf equations (GWHEs).
The method starts with an extension of transverse equation theory for layered arbitrary linear media applied to the stratification of an angular shaped region with the help of abstract BM notation. We then apply the characteristic Green function procedure to obtain the solution of equations in angular geometries. The solutions defined at the faces of the angular region are spectral functional equations that relate continuous (tangential) field components of the two faces delimiting a homogeneous angular region. The application of boundary conditions yields a system of generalized Wiener–Hopf equations (GWHEs), where generalized means that the field components of each face are defined into different complex planes but related to one another. The GWHEs preserve the characteristic form of classical Wiener–Hopf equations (CWHEs), where the system of equations presents a kernel and plus and minus unknowns; however, the plus and minus unknowns are defined into different complex planes which are related to one another. The functional equations and GWHEs of angular regions can be suitably interpreted with a network formalism, as is common for classical layered regions using transmission line theory. This circuit/network modeling representation of angular regions allows for the description of the technique with systematic steps to avoid redundancy. This capability is particularly useful when dealing with complex scattering problems, where we break down the complexity of a geometry into subdomains of canonical shape. These subdomains are modeled via spectral functional equations or related integral representations that can be interpreted through a network approach (obtained once and for all) and are capable of modeling the entire complex problem through the composition of circuital relationships; see, for instance, [30,31,32].
In the presence of an isotropic medium (and further special cases of more general media), a suitable mapping reduces the GWHEs to CWHEs, which in certain cases are amenable of exact solutions; alternatively, we can resort to the semi-analytical/approximate general-purpose factorization method, that is, Fredholm factorization, which reduces the factorization problem to Fredholm Integral Equations (FIEs) of the second kind. This technique was presented in the early 2000s for CWHEs; it was first effectively applied in impenetrable and penetrable wedge problems immersed in isotropic media [11,12], and more recently in isotropic complex scattering problems [30,31,32].
The main constraint in the present work resides in the complexity of media that do not allow mappings between complex planes of GWHEs for their transformation into CWHEs. Consequently, when dealing with arbitrary linear media in particular, we propose relying on a novel version of the versatile approximate method known as Fredholm factorization. Here, for the first time, we apply the Fredholm factorization method directly to GWHEs as a regularization tool. This regularized method can also be derived before the imposition of boundary conditions, i.e., directly on spectral functional equations (that is, before obtaining the GWHEs of the problem) by reversing the classical order of imposing boundary conditions and then applying Fredholm regularization to obtain the same effectiveness. We call this new methodology direct Fredholm factorization.
We observe that the impossibility of mapping GWHEs to CWHEs in arbitrary linear media is similar to the impossibility of defining a unique angular complex plane for the SM, KL, and WH methods; however, the new WH methodology proposed in this study overcomes this obstacle by resorting to direct Fredholm factorization applied to GWHEs.
From the solution of the GWHEs inherent to the angular region problem, we obtain the spectral representation of the field components along the faces delimiting homogeneous angular regions. The complete spectral analysis of the diffraction problems is then obtained by resorting again to spectral functional equations written for an arbitrary azimuthal direction. Finally, spectral inversion yields field components in the physical domain for any point in the angular regions. An alternative method for obtaining the field is also proposed; because of the linearity, it is based on the use of superposition on spectral representations prior to spectral inversion by identifying the spectral contributions of the faces of the angular regions through the equivalence theorem.
All of the theoretical properties of the mathematical statements are fully described in the text, although completely rigorous mathematical proofs are sometimes limited. On the other hand, validation of the proposed novel theoretical package through examples is reported, starting by demonstrating the effectiveness of direct Fredholm factorization applied to GWHEs in the scattering from a PEC wedge immersed in an isotropic medium and ending with validation of functional equations of angular regions in arbitrary linear media with the analysis of a PEC half-plane immersed in particular anisotropic media.
While implementing the method, we observe that the main difficulty resides in correctly estimating the kernel functions in the GWHEs and the corresponding FIE formulations for the presence of multivalued functions that need particular attention in their definition and calculation. The following sections highlight all multivalued functions and their correct estimations and assumptions.
In summary, we highlight in brief the main novelties of this work with respect to the state of the art reported in the introduction:
  • Development of a novel spectral method capable of handling scattering in arbitrary linear media with multiple propagation constants.
  • Introduction of a novel solution procedure for GWHEs, particularly for those with multiple propagation constants, which we call direct Fredholm factorization.
  • A description of spectral functional equations and related integral representations for angular regions filled by arbitrary linear media in terms of the network interpretation.
  • Computation of the field at each point within the angular region by resorting to the equivalence theorem and using Kirchhoff representations in the spectral domain;
  • Improving the quality of the approximate spectral solutions by reimposing GWHEs, referred to as iteration.
It is important to highlight that the applicability of the proposed WH technique to arbitrary linear media resides in its formulation directly in the Laplace domain while avoiding other complex planes; notably, other techniques such as SM use complex angular planes based on Sommerfeld representations, which are applicable only to isotropic media or special cases of anisotropic media. Moreover, although SM also uses Fredholm integral equations in the complex angular plane for approximate solutions [2,6], it is limited to isotropic media; again, the proposed WH method is extended to arbitrary linear media with direct Fredholm factorization because it is directly formulated in the Laplace domain. Furthermore, another important result is that while Sommerfeld–Malyuzhinets solutions combined with asymptotic methods require analytical extension of the spectral solutions in the improper sheet to compute the far field, our application of equivalence theorem in the context of the proposed method can be directly applied to approximate WH spectral solutions in the Laplace domain. This result is due to the direct solution of the GWHEs which provides the complete spectra of the field on the two faces of an angular region, useful and sufficient for asymptotic estimations.
This article is organized into seven sections and one appendix. In Section 1, we introduce the motivation and the scope of the present work, then report the state of the art related to the spectral analysis of diffraction in complex media. Section 2 presents the main mathematical steps for obtaining spectral functional equations in an angular region filled by an arbitrary linear medium and with arbitrary boundary conditions, starting from the abstract BM notation for transverse equations in layered planar regions and extending this theory to layered angular regions filled by arbitrary linear media. Section 3 develops the theory, starting from the spectral functional equation to obtain regularized integral representations for angular regions in arbitrary linear media with the direct application of the Fredholm factorization method. If boundary conditions are applied, the representations are GWHEs. Section 4 presents the route for obtaining an asymptotic estimation of the far field inside the angular region after the face spectra on the two limiting faces have been obtained. To demonstrate the efficacy of the proposed methodology, in particular direct Fredholm factorization, Section 5 reports a validation using the simple case of a PEC wedge immersed in an isotropic medium. To further validate the method in arbitrary linear media, Section 6 presents an example of the application of functional equations in arbitrary linear media, specifically, a PEC half-plane immersed in a gyrotropic medium. Section 7 presents the conclusion of the paper. Finally, the Appendix A reports the full and explicit formulas and equations, as abstract notation is used in the main text with the dual purpose of enhancing readability and ensuring completeness.

2. Spectral Functional Equations in an Angular Region Filled by Arbitrary Linear Media

Spectral functional equations in angular regions filled by arbitrary linear media are obtained by exploiting the combination and extension of the following powerful mathematical tools developed in different contexts: first, the Bresler–Marcuvitz (BM) transverse equation theory for layered media [33,35], and second, the characteristic Green function procedure [36,37]. In this section, following [34], we first briefly revisit the BM theory for layered planar arbitrary linear media as a fundamental step in analyzing layered angular regions. We then apply the characteristic Green function procedure to derive solutions of the obtained system of differential equations. Finally, we provide the spectral functional equations by evaluating the solution at the faces of the angular region. In particular, the functional equations relate continuous (tangential) spectral field components defined at the two faces of the angular region.
We start from the application of BM theory to Maxwell’s equations in layered arbitrary nondispersive homogeneous linear media with tensorial constitutive relations (i.e., bianisotropic media [25,26,27]):
D = ε ̲ ̲ · E + ξ ̲ ̲ · H B = ζ ̲ ̲ · E + μ ̲ ̲ · H
where the electric and magnetic fields ( E , H ) are related to the electric and magnetic fluxes ( D , B ) and the tensors ( ε ̲ ̲ , μ ̲ ̲ , ξ ̲ ̲ , ζ ̲ ̲ ) are, respectively, the electric permittivity, the magnetic permeability, and the two magneto-electric coupling parameters.
By assuming
(a)
Cartesian coordinates ( z , x , y ) ,
(b)
e + j ω t time–harmonic field dependence,
(c)
Invariant geometry along z and stratification along y, and
(d)
Sources constituted of plane waves with z-dependence e j α o z , where α o depends on the skewness angle with respect to z (with α o = 0 at normal incidence on z),
We obtain the transverse differential equations in matrix form for layered planar media:
y ψ y ( x , y ) = M y ( j α o , x ) · ψ y ( x , y )
where ψ y is the four-dimensional column vector (Throughout the article, we assume the notation | | for vectors and not for the modulus of a vector).
ψ y = | E t t , H t t | t , w i t h E t = | E z , E x | t , H t = | H z , H x | t .
Based on the nature of Maxwell’s equations, M y ( j α o , x ) is a second-order four-dimensional matrix differential operator of the form
M y ( j α o , x ) = M y o + M y 1 x + M y 2 2 x 2 ,
where the explicit forms of the matrices M y o , M y 1 , M y 2 for an arbitrarily linear medium (1) are reported in Appendix A with (A2)–(A9). The application of the Fourier transform along x reduces (2) to
d d y ψ y ( η , y ) = M y ( j α o , j η ) · ψ y ( η , y ) ,
where ψ y ( x , y ) 1 2 π ψ y ( η , y ) e j η x d η and
M y ( j α o , j η ) = M y o j η M y 1 η 2 M y 2 .
We introduce here an analysis of the operator M y ( j α o , j η ) of the layered planar arbitrarily linear media; this is necessary to obtain the solution of (2) in terms of eigenvalues, eigenvectors (with the characteristic Green function procedure), and boundary conditions. The same analysis is needed to obtain a solution for layered angular arbitrarily linear media. Supposing for the general case (removing exceptions) that M y is semi-simple, we compute its eigenvalues λ i and eigenvectors u i as follows:
M y u i = λ i u i
i.e.,
M y = U y J y U y 1
where J y = d i a g { λ 1 , λ 2 , λ 3 , λ 4 } and U y = ( u 1 , u 2 , u 3 , u 4 ) (dependence on η and α o is omitted). The computation of the eigenvalues is obtained from the zeros of the characteristic equation of order four (9), the coefficients of which can be written using Bocher’s formula, as shown below [39].
d e t [ M y λ i I ] = λ i 4 + a λ i 3 + b λ i 2 + c λ i + d = 0
a = t r ( M y ) , b = a t r ( M y ) + t r ( M y 2 ) 2 , c = b t r ( M y ) + a t r ( M y 2 ) + t r ( M y 3 ) 3 , d = d e t [ M y ]
This yields the following four eigenvalues:
λ 1 = a 4 + T + M a + Q 2 , λ 2 = a 4 + T M a + Q 2 ,
λ 3 = a 4 T + M a Q 2 , λ 4 = a 4 T M a Q 2
where
T = a 2 4 + 3 a c + b 2 + 12 d 3 u 3 + u 3 2 b 3 , Q = a 3 4 a b + 8 c 4 T , M a = 3 a 2 4 2 b T
with
u = s + v 2 , v = 9 3 a 2 d a b c 8 b d + 3 c 2 + 2 b 3 , s = v 2 4 3 a c + b 2 + 12 d 3 .
We note that the column vectors u i = 1 , 2 , 3 , 4 of U y provide a basis in the space C 4 , where we define the transverse electromagnetic field ψ y , while the column vectors v i = 1 , 2 , 3 , 4 of
V y = U y 1
in the reciprocal space will be fundamental in obtaining functional equations through the characteristic Green function procedure. Each couple ( u i , v i ) is related to a single λ i for which the explicit forms are, in general, the cumbersome expressions reported in (11) and (12), which depend on η . In the most simple case, i.e., the isotropic medium ( ε ̲ ̲ = ε I ̲ ̲ , μ ̲ ̲ = μ I ̲ ̲ , ξ ̲ ̲ = ζ ̲ ̲ = 0 ̲ ̲ ), λ i assumes the following forms:
λ 1 = λ 2 = λ 3 = λ 4 = ( α o 2 + η 2 ) k 2 = j ( k 2 α o 2 ) η 2 = j ξ i s o , k = ω ε μ
where in the presence of losses ( k = k r j k i ; k r , k i > 0 ) we have Re [ λ 1 , 2 ] > 0 and Re [ λ 3 , 4 ] < 0 , which are respectively related to progressive ( i = 1 , 2 ) and regressive ( i = 3 , 4 ) waves with regard to y of form e j η x e λ i y e j α o z . In this framework, we associate the direction of propagation with attenuation phenomena; we let the phase variation be free of constraints to allow modeling of left-handed materials. In a general arbitrary (even small) lossy linear medium, we always have two eigenvalues (say, i = 1 , 2 ) with a positive real part λ i = + j ξ i representing progressive waves and two (say, i = 3 , 4 ) with a negative real part λ i = j ξ i representing regressive waves, yielding all four y-longitudinal propagation constants with Im [ ξ i ] < 0 (y-progressive/regressive waves e j ξ i y , i.e., assuming time-harmonic dependence e + j ω t , we have x, y, z progressive waves e j η x e j ξ i y e j α o z , respectively, with Im [ η , ξ i , α o ] < 0 ).
We affirm here the importance of retaining the generality of the medium; while investigating the scattering of objects immersed in arbitrary linear media, the scatterer can be arbitrarily oriented with respect to the principal axis of the (crystal) medium. However, when the problem permits the definition of a coordinate system that coincides with the principal axes of the crystal medium, we obtain tensorial constitutive relations with diagonal tensors (1). These media are called biaxial, uniaxial, and isotropic, while the three terms in the diagonal are, respectively, all different, one different, and all equal. Other special cases are gyrotropic media, which represent media in which the tensorial constitutive relations are Hermitian with respect to the coordinate system and of the following type (in (17), we limit the description to the permittivity, i.e., a gyroelectric medium):
ε ̲ ̲ = ε + j ε g 0 j ε g ε 0 0 0 ε a .
Starting from planar layered regions, we extend the theory to angular-shaped regions of aperture γ filled by arbitrary linear media, as already done in Section 3 of [34] for angular regions filled by isotropic media. With reference to region 1 in Figure 1a, we derive from (2) the oblique transverse Equation (19) using an oblique system of Cartesian axes ( z , u x , v ) :
x = u + v cos γ , y = v sin γ ,
v ψ y ( u , v ) = M γ ( j α o , u ) · ψ y ( u , v ) .
The application of a Fourier transform along u = x reduces (19) to
d d v ψ y ( η , v ) = M γ ( j α o , j η ) · ψ y ( η , v ) ,
where ψ y ( u , v ) 1 2 π ψ y ( η , v ) e j η u d η and
M γ ( j α o , j η ) = M γ o j η M γ 1 η 2 M γ 2 ,
M γ o = M yo sin γ , M γ 1 = M y 1 sin γ I t cos γ , M γ 2 = M y 2 sin γ .
Based on the link between M γ and M y , we have M γ with the same eigenvectors u i of M y , and the following relationship between the eigenvalues λ γ i ( γ ) and λ i can be obtained:
λ γ i ( γ ) = j η cos γ + λ i sin γ , i = 1 4
resulting in the following “oblique” v-longitudinal propagation constants:
m i ( γ ) = j λ γ i ( γ ) = + η cos γ + ξ i sin γ , i = 1 , 2
m i ( γ ) = + j λ γ i ( γ ) = η cos γ + ξ i sin γ , i = 3 , 4
which is in agreement with the relationship between λ i and ξ i as well as with the correlated progressive and regressive propagating interpretation along the longitudinal direction y and along the oblique “longitudinal” direction v (progressive/regressive e j m i v ). We note that the quantities M γ ( j α o , j η ) , λ γ i ( γ ) , and m i ( γ ) depend on the geometrical parameter γ and on the spectral variable η .
With reference to region 1 in Figure 1a, we obtain the functional equations with a circuital interpretation as a mathematical manipulation of the solution of the differential Equation (19) using the Laplace domain ψ ˜ y ( η , v ) 0 e j η u ψ y ( u , v ) d u .
d d v ψ ˜ y ( η , v ) = M γ ( j α o , j η ) · ψ ˜ y ( η , v ) + ψ s a ( v ) , v > 0
ψ s a ( v ) = M γ 1 · ψ y ( 0 + , v ) + j η M γ 2 · ψ y ( 0 + , v ) M γ 2 · u ψ y ( u , v ) u = 0 +
The benefit of using the Laplace transform is correlated with the incorporation of boundary conditions through initial conditions with the term ψ s a ( v ) . In (26)–(27), the condition u = 0 + , v > 0 imposes boundary conditions on the fields along face a in Figure 1a. The solution is obtained by using the characteristic Green function procedure [34] in terms of homogeneous and particular solutions, yielding the representation shown below.
ψ ˜ y ( η , v ) = i = 1 4 C i e λ γ i ( γ ) v u i i = 1 2 u i v i · 0 v e λ γ i ( γ ) ( v v ) ψ s a ( v ) d v + i = 3 4 u i v i · v e λ γ i ( γ ) ( v v ) ψ s a ( v ) d v
Now, considering the asymptotic behavior in (28) of exponential functions in v for v + , we need to have C 3 = C 4 = 0 ; at the same time, the first couple of integrals are null, as R e [ λ 1 , 2 ] > 0 and R e [ λ 3 , 4 ] < 0 are related to progressive and regressive waves, respectively. For this reason, setting v = 0 , we obtain the following spectral field representation along face o:
ψ ˜ o + ( η ) ψ ˜ y ( η , 0 ) = C 1 u 1 + C 2 u 2 + i = 3 4 u i v i · 0 e λ γ i ( γ ) ( v v ) ψ s a ( v ) d v .
By weighting (29) with the reciprocal vectors v 3 , v 4 of M γ , we obtain the following functional equations:
v i · ψ ˜ o + ( η ) = v i · ψ ˜ s a + ( m i ( γ ) ) , i = 3 , 4
where we use the following definition of the Laplace transform:
ψ ˜ s a + ( m i ( γ ) ) 0 e j m i ( γ ) v ψ s a ( v ) d v = 0 e j m i ( γ ) ρ ψ s a ( ρ ) d ρ .
With reference to Figure 1a and its caption, analyzing and expanding ψ ˜ s a + ( m i ( γ ) ) in (30) using Maxwell’s equations, we rephrase the functional equations for region 1 as
v i · ψ ˜ o + ( η ) = v i · T ( γ ) · ψ ˜ a + ( m i ( γ ) , γ ) , i = 3 , 4 ,
where ψ ˜ o + ( η ) is the η Laplace transform of the tangent-to-face-o field components (i.e., at φ = 0 ) in Cartesian ( z , x , y ) and cylindrical coordinates ( z , ρ , φ ) (omitting the z coordinate for invariance)
ψ ˜ o + ( η ) = 0 | E z ( x , 0 ) , E x ( x , 0 ) , H z ( x , 0 ) , H x ( x , 0 ) | t e j η x d x = 0 | E z ( ρ , 0 ) , E ρ ( ρ , 0 ) , H z ( ρ , 0 ) , H ρ ( ρ , 0 ) | t e j η ρ d ρ ,
where ψ ˜ a + ( m i ( γ ) , γ ) is the m i ( γ ) Laplace transform of the tangent-to-face-a field components (i.e., at φ = γ ) in local-to-face-a Cartesian ( z , X 1 , Y 1 ) coordinates and global cylindrical coordinates ( z , ρ , φ ) (located at φ = + γ )
ψ ˜ a + ( m i ( γ ) , γ ) = 0 | E z ( X 1 , 0 ) , E X 1 ( X 1 , 0 ) , H z ( X 1 , 0 ) , H X 1 ( X 1 , 0 ) | t e j m i ( γ ) X 1 d x = 0 | E z ( ρ , γ ) , E ρ ( ρ , γ ) , H z ( ρ , γ ) , H ρ ( ρ , γ ) | t e j m i ( γ ) ρ d ρ ,
and
T ( γ ) = sin ( γ ) ( α o ξ y y + ζ x y ξ y y ω μ x y ω ϵ y y ) ω ( μ y y ϵ y y ζ y y ξ y y ) + cos ( γ ) 0 sin ( γ ) ( α o μ y y + ζ x y μ y y ω ζ y y μ x y ω ) ω ( μ y y ϵ y y ζ y y ξ y y ) 0 sin ( γ ) ( ζ z y ξ y y ω + η ξ y y + μ z y ω ϵ y y ) ω ( μ y y ϵ y y ζ y y ξ y y ) 1 sin ( γ ) ( ζ y y μ z y ω ζ z y μ y y ω + η μ y y ) ω ( μ y y ϵ y y ζ y y ξ y y ) 0 sin ( γ ) ( α o ϵ y y ξ y y ω ϵ x y + ξ x y ω ϵ y y ) ω ( μ y y ϵ y y ζ y y ξ y y ) 0 cos ( γ ) sin ( γ ) ( α o ζ y y ζ y y ξ x y ω + μ y y ω ϵ x y ) ω ( μ y y ϵ y y ζ y y ξ y y ) 0 sin ( γ ) ( ξ y y ω ϵ z y ϵ y y ( η + ξ z y ω ) ) ω ( μ y y ϵ y y ζ y y ξ y y ) 0 sin ( γ ) ( μ y y ω ϵ z y ζ y y ( η + ξ z y ω ) ) ω ( μ y y ϵ yy ζ y y ξ y y ) 1
Note that (32) are functional equations that relate the Laplace transforms of combinations of field components on the boundaries of angular region 1 in Figure 1a, i.e., face o  u > 0 , v = 0 ( φ = 0 ) and face a  u = 0 , v > 0 ( φ = γ ). Furthermore, we observe that the angle γ is essential in determining the impact of anisotropies through T ( γ ) .
Repeating the same procedure for region 2 in Figure 1a, we obtain the functional equations as the solution of the differential Equation (19) in the Laplace domain using the left Laplace transform ψ ˜ y ( η , v ) 0 e j η u ψ y ( u , v ) d u :
d d v ψ ˜ y ( η , v ) = M γ ( j α o , j η ) · ψ ˜ y ( η , v ) + ψ s b ( v ) , v < 0
where ψ s b ( v ) has the same expression of ψ s a ( v ) (27) but with a different support v < 0 , and allows the incorporation of boundary conditions along face b ( u = 0 + , v < 0 ). The application of the characteristic Green function procedure yields expression (28) for region 2 in Figure 1a, which is identical to that for region 1 except for C i and the source term ψ s b ( v ) , which depend on the local constitutive parameters and boundary conditions of region 2. Now, considering the asymptotic behavior of the exponential function in v for v , we need to have C 1 = C 2 = 0 ; at the same time, the second couple of integrals are null. For this reason, setting v = 0 , we obtain
ψ ˜ o + ( η ) ψ ˜ y ( η , 0 ) = C 3 u 3 + C 4 u 4 i = 1 2 u i v i · 0 e λ i ( γ ) ( v v ) ψ s b ( v ) d v .
By weighting (37) with the reciprocal vectors v 1 , v 2 of M γ , we obtain the functional equations
v i · ψ ˜ o + ( η ) = v i · ψ ˜ s b + ( m i ( γ ) ) , i = 1 , 2 ,
where we have used the definition of the v left Laplace transform
ψ ˜ s b + ( m i ( γ ) ) 0 e j m i ( γ ) v ψ s b ( v ) d v = 0 e j m i ( γ ) ρ ψ s b ( ρ ) d ρ .
Note the differences and similarities between Laplace transformations (31) and (39), which yield the same definition of the m i ( γ ) Laplace transform in ρ but are applied to different quantities. Furthermore, the regularity properties of the m i ( γ ) Laplace transform are inherited from ξ i ( I m [ ξ i ] < 0 ) according to (24)–(25).
With reference to Figure 1a and its caption, analyzing and expanding ψ ˜ s b + ( m i ( γ ) ) in (38), we rephrase the functional equations into
v i · ψ ˜ o + ( η ) = v i · T ( γ ) · P · ψ ˜ b + ( m i ( γ ) , π + γ ) , i = 1 , 2 .
In (40), T ( γ ) is again (35); as for region 1, P = d i a g { 1 , 1 , 1 , 1 } is needed for v = X 2 in region 2 with respect to v = X 1 in region 1, ψ ˜ o + ( η ) is the η Laplace transform of the tangent-to-face-o field components reported in (33), and ψ ˜ b + ( m i ( γ ) , π + γ ) is the m i ( γ ) Laplace transform of the tangent-to-face-b field components (i.e., at φ = π + γ ) in the local-to-face-b Cartesian ( z , X 2 , Y 2 ) coordinates and global cylindrical coordinates ( z , ρ , φ ) of Figure 1a.
ψ ˜ b + ( m i ( γ ) , π + γ ) = 0 | E z ( X 2 , 0 ) , E X 2 ( X 2 , 0 ) , H z ( X 2 , 0 ) , H X 2 ( X 2 , 0 ) | t e j m i ( γ ) X 2 d x = 0 | E z ( ρ , π + γ ) , E ρ ( ρ , π + γ ) , H z ( ρ , π + γ ) , H ρ ( ρ , π + γ ) | t e j m i ( γ ) ρ d ρ
While considering the wedge scattering problem with symmetry with respect to the x axis, in combination with region 1 in Figure 1a, we need to consider region 2′ in Figure 1b, where γ π γ with respect to region 2 in Figure 1a, i.e., for the same face a, we change the orientation of face b at φ = γ from φ = π + γ to φ = γ . The functional equations of region 2′ become
v i · ψ ˜ o + ( η ) = v i · T ( π γ ) · P · ψ ˜ b + ( m i ( π γ ) , γ ) , i = 1 , 2 ,
where
ψ ˜ b + ( m i ( π γ ) , γ ) = 0 | E z ( X 2 , 0 ) , E X 2 ( X 2 , 0 ) , H z ( X 2 , 0 ) , H X 2 ( X 2 , 0 ) | t e j m i ( π γ ) X 2 d x = 0 | E z ( ρ , γ ) , E ρ ( ρ , γ ) , H z ( ρ , γ ) , H ρ ( ρ , γ ) | t e j m i ( π γ ) ρ d ρ ,
which is the m i ( π γ ) Laplace transform of the tangent-to-face-b field components (i.e., now at φ = γ ) in local-to-face-b Cartesian ( z , X 2 , Y 2 ) coordinates and global cylindrical coordinates ( z , ρ , φ ) in Figure 1b. Note that in (42) we assumed that region 2′ is homogeneous to region 1, yielding the same u i , v i ; otherwise, specific vectors would be needed. Equations (42) are functional equations that relate the Laplace transforms of combinations of field components on the boundaries of the angular region 2′ in Figure 1b, i.e., face o  u > 0 , v = 0 ( φ = 0 ) and face b  u = 0 , v < 0 ( φ = γ ). In (42), note the new dependence of T ( · ) (35) on π γ due to the effect of anisotropies while changing the orientation of face b from π + γ to γ . Furthermore, in the case of symmetric media ( λ 1 , 2 = λ 3 , 4 ), we have m 3 , 4 ( γ ) = m 1 , 2 ( π γ ) ; see (24)–(25).
In general, the system of functional equations in (32) and (42) allows for the analysis of angular regions that are symmetric with respect to the x axis; these are at the base of the analysis of scattering problems constituted by impenetrable and penetrable wedges surrounded/made by arbitrary linear media. In the following, to investigate a practical scattering problem, we impose boundary conditions at the faces of each angular region on the functional Equations (32) and (42), yielding a system of GWHEs.

3. From Functional Equations to GWHEs and Their Regularized Integral Representations with a Network Interpretation

Network representations of angular regions in isotropic media for electromagnetic scattering were extensively studied in multiple spectral domains in [13] using algebraic and integral formalisms. The proposed equations were effectively applied in several works to practical wedge scattering problems; see [11,12] and the references therein. Furthermore, network formalisms have been effectively applied for complex canonical problems containing angular and layered regions in isotropic media; see for instance the double wedge [31], flanged dielectric loaded waveguide [32], and wedge over dielectric layer [30], among others.
In arbitrary linear media, the system of functional Equations (32) and (42),
v i · ψ ˜ o + ( η ) = v i · T ( γ ) · ψ ˜ a + ( m i ( γ ) , γ ) , i = 3 , 4 v i · ψ ˜ o + ( η ) = v i · T ( π γ ) · P · ψ ˜ b + ( m i ( π γ ) , γ ) , i = 1 , 2
constitutes two systems of network relations that link the spectral field components in region 1 and region 2′ (Figure 1), respectively, via a sort of two-port transmission relations in algebraic form. Looking at the first system in (44), we have two combinations of ψ ˜ o + ( η ) components (33) related to two combinations of ψ ˜ a + ( m i ( γ ) , γ ) components (34), i.e., with reference to Figure 1a, the tangential field components of face o are related to tangential field components of face a. A similar interpretation can be repeated for the second system in (44) for region 2 with field components defined at faces o and b, respectively, in ψ ˜ o + ( η ) (33) and ψ ˜ b + ( m i ( π γ ) , γ ) (43).
We further note that in (44) the components of face o and face a , b are functions of the spectral variables η and m i ( · ) , respectively, which are related to one another via (24)–(25). We can reverse the role of the variables η and m i ( · ) in the arguments of the components of these faces. In this way, we double the equations of region 1, i.e., the first line of (44), which is again reported in (45), yielding the second line of (45), which relates the components of face a (now functions of the variable η ) to the components of face o (now functions of m i ( · ) ). The second line of (45) is obtained by defining region 1 as region 2′ (Figure 1) after a clockwise rotation of an angle + γ , yielding the following complete set of equations for region 1:
v i · ψ ˜ o + ( η ) = v i · T ( γ ) · ψ ˜ a + ( m i ( γ ) , γ ) , i = 3 , 4 , v iY 1 · ψ ˜ a + ( η ) = v iY 1 · T Y 1 ( π γ ) · P · ψ ˜ o + ( m i Y 1 ( π γ ) , γ ) , i = 1 , 2 .
In the second pair of Equation (45), we have used the subscript Y 1 to refer to a rotated coordinated system ( z , X 1 , Y 1 ) with respect to ( z , x , y ) ; see region 1 in Figure 2a and the related region 2′ in Figure 2b.
We note that the second pair of equations in (45) can be easily derived by studying the classical region 2′ (see the second pair of equations in (44)) except with modified definitions of the quantities v iY 1 , T Y 1 ( γ ) , m i Y 1 ( γ ) (from λ i Y 1 ( γ ) ) because of their dependence on the constitutive tensorial parameters ( ε ̲ ̲ , μ ̲ ̲ , ξ ̲ ̲ , ζ ̲ ̲ ) of region 1, which are redefined in the reference coordinate system ( z , X 1 , Y 1 ) , i.e., ( ε ̲ ̲ Y 1 , μ ̲ ̲ Y 1 , ξ ̲ ̲ Y 1 , ζ ̲ ̲ Y 1 ) , for example,
ε ̲ ̲ Y 1 = R ̲ ̲ Y 1 1 · ε ̲ ̲ · R ̲ ̲ Y 1 , R ̲ ̲ Y 1 = cos ( γ ) sin ( γ ) 0 sin ( γ ) cos ( γ ) 0 0 0 1 ,
due to the rotation of + γ .
The same rationale is applied to region 2 to double the equations of that region (second line of (44), as also reported in (47)) by obtaining
v iY 2 · ψ ˜ b + ( η ) = v iY 2 · T Y 2 ( γ ) · ψ ˜ o + ( m i Y 2 ( γ ) , γ ) , i = 3 , 4 , v i · ψ ˜ o + ( η ) = v i · T ( π γ ) · P · ψ ˜ b + ( m i ( π γ ) , γ ) , i = 1 , 2 .
In the first pair of equations in (47), we have used the subscript Y 2 to make reference to a rotated coordinated system ( z , X 2 , Y 2 ) with respect to ( z , x , y ) ; see region 2 in Figure 2c and the related region 1′ in Figure 2d. We note that the first pair of equations in (47) can be easily derived by studying the classical region 1 (see the first pair of equations in (44)) except with modified definitions of v iY 2 , T Y 2 ( γ ) , and m i Y 2 ( γ ) because of their dependence on the constitutive parameters ( ε ̲ ̲ , μ ̲ ̲ , ξ ̲ ̲ , ζ ̲ ̲ ) , which are redefined in the reference coordinate system ( z , X 2 , Y 2 ) , i.e., ( ε ̲ ̲ Y 2 , μ ̲ ̲ Y 2 , ξ ̲ ̲ Y 2 , ζ ̲ ̲ Y 2 ) , for example,
ε ̲ ̲ Y 2 = R ̲ ̲ Y 2 1 · ε ̲ ̲ · R ̲ ̲ Y 2 , R ̲ ̲ Y 2 = cos ( γ ) sin ( γ ) 0 sin ( γ ) cos ( γ ) 0 0 0 1 ,
due to the rotation of γ .
The sets of Equations (45) and (47) constitute a complete set of functional equations that respectively describe regions 1 and 2 in Figure 3. In the case of symmetric media (i.e., λ 1 , 2 = λ 3 , 4 ) we have m 3 , 4 ( γ ) = m 1 , 2 ( π γ ) ; see (24)–(25).
In isotropic media, it is always possible to introduce the angular complex plane w and the KL transform method [13], where functional equations become two-port admittance relations of the Norton type in integral and algebraic forms in a unique complex plane. In arbitrary linear media, the definition of such complex planes is not possible; however, a novel method that resorts to the following Cauchy decomposition formula in the m ( η ) plane is introduced. This is a fundamental tool that allows for the description of angular region problems in arbitrary linear media without introducing further complex planes other than the initial Laplace transforms. In particular, to obtain regularized integral equations from GWHEs, it is not necessary to map the GWHEs into CWHEs with suitable transformations before the application of Fredholm factorization (originally ideated and valid only for the CWHEs). This revisited novel version of the regularization procedure can be called the direct Fredholm factorization method.
At the origin of this method, we introduce the following generalized form of the standard Cauchy decomposition formula in the m ( η ) plane (i.e., one of the m i ( · ) that all depends on η ) applied to an arbitrary F + ( m ( η ) ) (i.e., the standard form is obtained simply by replacing m ( η ) with η ):
F + ( m ( η ) ) = 1 2 π j F + ( η ) η + m ( η ) d η + F + n . s . ( m ( η ) ) , η R
where F + n . s . ( m ( η ) ) is the non-standard contribution of F + ( m ( η ) ) in the m ( η ) plane. Assuming a lossy medium, we observe that, in general, m ( η ) is located with the positive imaginary part for η R , i.e., it is in the upper half-plane of complex plane η ; thus, the application of (49) to plus functions is justified. For example, see Figure 4, where we have assumed that k = 1 0.1 j and γ = 0.7 π , yielding a path of m ( t ) for t R from right to left because of γ > π / 2 . On the contrary, for γ < π / 2 we obtain a similar path located in the upper half-plane, but in the opposite direction. We anticipate that the application of (49) to GWHEs with multiple propagation constants, i.e., with more than one m i ( η ) , is fundamental for developing a solution in the η plane, as (49) transforms the GWHEs into integral equations in the unique complex plane η .
The complete sets of Equation (45) for region 1 can be represented in the following form:
A 11 E ( η ) E o z ( η ) + A 12 E ( η ) E o ρ ( η ) + A 11 ( η ) H o z ( η ) + A 12 ( η ) H o ρ ( η ) = B 11 E ( η ) E a z ( m 3 ) + B 12 E ( η ) E a ρ ( m 3 ) + B 11 ( η ) H a z ( m 3 ) + B 12 ( η ) H a ρ ( m 3 ) A 21 E ( η ) E o z ( η ) + A 22 E ( η ) E o ρ ( η ) + A 21 ( η ) H o z ( η ) + A 22 ( η ) H o ρ ( η ) = B 21 E ( η ) E a z ( m 4 ) + B 22 E ( η ) E a ρ ( m 4 ) + B 21 ( η ) H a z ( m 4 ) + B 22 ( η ) H a ρ ( m 4 ) A 31 E ( η ) E a z ( η ) + A 32 E ( η ) E a ρ ( η ) + A 31 ( η ) H a z ( η ) + A 32 ( η ) H a ρ ( η ) = B 31 E ( η ) E o z ( m 1 ) + B 32 E ( η ) E o ρ ( m 1 ) + B 31 ( η ) H o z ( m 1 ) + B 32 ( η ) H o ρ ( m 1 ) A 41 E ( η ) E a z ( η ) + A 42 E ( η ) E a ρ ( η ) + A 41 ( η ) H a z ( η ) + A 42 ( η ) H a ρ ( η ) = B 41 E ( η ) E o z ( m 2 ) + B 42 E ( η ) E o ρ ( m 2 ) + B 41 ( η ) H o z ( m 2 ) + B 42 ( η ) H o ρ ( m 2 )
where the spectral field components of face o and face a are related to one another (Throughout this study, in spectral equations we assume a notation with two subscripts for the spectral field; the first subscript is related to the considered face ( o , a , b ) and the second to the field component ( z , x , y ) ). Moreover, the complete set of equations for region 2 (47) has a similar representation. The imposition of boundary conditions makes these equations a well-posed mathematical problem, resulting in a system of GWHEs. In particular, if the region is surrounded by something modeled with impenetrable impedance boundary conditions, then we establish relations among field components on the boundary faces. On the contrary, if the region is surrounded by penetrable regions, then we establish continuity through tangent components to neighboring regions that provide further functional equations (coupled together). In any case, the types of completed functional equations and constraints with boundary conditions are always of the same form and are a well-posed mathematical problem of GWHEs.
As a simple example to illustrate the procedure, let us consider a problem constituted by only region 1 with PEC boundary conditions when filled by arbitrary linear media. In this case, we obtain
A 11 ( η ) H o z ( η ) + A 12 ( η ) H o ρ ( η ) = B 11 ( η ) H a z ( m 3 ) + B 12 ( η ) H a ρ ( m 3 ) A 21 ( η ) H o z ( η ) + A 22 ( η ) H o ρ ( η ) = B 21 ( η ) H a z ( m 4 ) + B 22 ( η ) H a ρ ( m 4 ) A 31 ( η ) H a z ( η ) + A 32 ( η ) H a ρ ( η ) = B 31 ( η ) H o z ( m 1 ) + B 32 ( η ) H o ρ ( m 1 ) A 41 ( η ) H a z ( η ) + A 42 ( η ) H a ρ ( η ) = B 41 ( η ) H o z ( m 2 ) + B 42 ( η ) H o ρ ( m 2 )
where on the LHS we have plus field unknowns in η and on the RHS we have minus field unknowns in m i ( ) . The apparent redundancy in (51) after the imposition of boundary conditions is exploited to obtain integral representations only in terms of the field components H o z ( η ) , H o ρ ( η ) , H a z ( η ) , H a ρ ( η ) in the unique complex plane η using (49). Furthermore, the application of the novel version of the Fredholm factorization method allows us to obtain regularized integral equations. We assert that this procedure is applicable to GWHEs in general, not only for the specific problem presented in this simple example. The application of (49) to the RHS of (51) yields the following:
A 11 ( η ) H o z ( η ) + A 12 ( η ) H o ρ ( η ) = B 11 ( η ) 2 π j H a z ( η ) η + m 3 d η + B 12 ( η ) 2 π j H a ρ ( η ) η + m 3 d η + H a z n . s ( m 3 ) + H a ρ n . s ( m 3 ) A 21 ( η ) H o z ( η ) + A 22 ( η ) H o ρ ( η ) = B 21 ( η ) 2 π j H a z ( η ) η + m 4 d η + B 22 ( η ) 2 π j H a ρ ( η ) η + m 4 d η + H a z n . s ( m 4 ) + H a ρ n . s ( m 4 ) A 31 ( η ) H a z ( η ) + A 32 ( η ) H a ρ ( η ) = B 31 ( η ) 2 π j H o z ( η ) η + m 1 d η + B 32 ( η ) 2 π j H o ρ ( η ) η + m 1 d η + H o z n . s ( m 1 ) + H o ρ n . s ( m 1 ) A 41 ( η ) H a z ( η ) + A 42 ( η ) H a ρ ( η ) = B 41 ( η ) 2 π j H a z ( η ) η + m 2 d η + B 42 ( η ) 2 π j H a ρ ( η ) η + m 2 d η + H o z n . s ( m 2 ) + H o ρ n . s ( m 2 )
recalling that all occurrences of m i are functions of η , i.e., m i ( η ) . The integral equations in (52) are of the singular type; for this reason, we resort to the Fredholm factorization method to obtain regularized expressions. The procedure consists of γ 1 t Cauchy smile contour integration [11] on both sides of each equation and consequent mathematical elaboration. Focusing our attention on the LHS for each term of each equation in (52), using dummy subscripts, we have the following regularized expression:
1 2 π j γ 1 t A ( t ) H + ( t ) t η d t = 1 2 π j γ 1 t ( A ( t ) A ( η ) ) H + ( t ) t η d t + A ( η ) 2 π j γ 1 t H + ( t ) t η d t = 1 2 π j ( A ( t ) A ( η ) ) H + ( t ) t η d t + A ( η ) H + ( η ) A ( η ) H + n . s ( η ) .
Focusing our attention on the RHS for each term of each equation in (52), using dummy subscripts and going back to the representation in (51), we have the following regularized expression:
1 2 π j γ 1 t B ( t ) H + ( m ( t ) ) t η d t = 1 2 π j γ 1 t ( B ( t ) B ( η ) ) H + ( m ( t ) ) t η d t + B ( η ) 2 π j γ 1 t H + ( m ( t ) ) t η d t = 1 ( 2 π j ) 2   ( B ( t ) B ( η ) ) H + ( η ) ( t η ) ( η + m ( t ) ) d t d η + B ( η ) ( 2 π j ) 2   γ 1 t 1 ( t η ) ( η + m ( t ) ) d t H + ( η ) d η + n . s . t e r m s .
Given the expressions of m i ( η ) in (24)–(25) with I m [ m i ( t ) ] < 0 in lossy media (Figure 4), Equation (54) requires the computation of
M e ( η , η ) = γ 1 t 1 ( t η ) ( η + m ( t ) ) d t ,
which can be performed either numerically or analytically while paying attention to the branch cuts of m ( t ) . Furthermore, in (54) we also need to consider n.s. singularities related to the field.
The validity of the estimation of M e ( η , η ) extends to complex values of η as long as η does not cross the singularity line determined by the behavior of m ( t ) for t R , as shown in Figure 4.
The expressions in (53) and (54) are regularized integral terms, as their kernels are compact. Moreover, they include the n.s. terms of the field components in η and m i , respectively. The detailed proof of this assertion needs to be performed for specific problems. While numerically implementing the method we observe that one of the main difficulties resides in the correct estimation of kernel functions A ( η ) , B ( η ) for the presence of multivalued functions that need particular attention in their definition and calculation.
For simplicity and compactness of discussion, we examine the properties of integral equations in the simple case of a PEC wedge immersed in an isotropic medium in Section 5.2. Equation (51) yields a 4 × 4 system of Fredholm integral equations of the second kind by utilizing (52), (53), and (54). This system is expressed in terms of H o z ( η ) , H o ρ ( η ) , H a z ( η ) , H a ρ ( η ) . It is important to highlight that the system only depends on the spectral variable η , ensuring that functions do not rely on m i outside of the integration sign. This property is fundamental to avoid the analysis of unknowns defined in different complex planes ( η and various m i ) that are correlated through cumbersome improper sheet properties.

4. Asymptotic Estimation of the Field in the Angular Region

Having obtained the spectra at the faces of the angular region, we can estimate the asymptotic behavior of the far field inside of the angular region.
Going back to the solution of (26) in Section 2 for region 1, according to (28) we have
ψ ˜ y ( η , v ) = i = 1 2 C i e λ γ i ( γ ) v u i + i = 3 4 u i v i · v e λ γ i ( γ ) ( v v ) ψ s a ( v ) d v , v > 0 .
From the homogeneous portion of the solution in (56), we obtain the definitions of arbitrary coefficients in terms of the field components at v = 0 (face o):
v i · ψ ˜ y ( η , 0 ) = C i , i = 1 , 2 .
The particular integrals in (56) are terms related to face a via ψ s a ( v ) . Due to the linearity of the problem, we apply the superposition principle; we can interpret (56) as the result of an equivalent theorem, where ψ ˜ y ( η , v ) is represented through equivalent sources at faces o and a. Similarly, the spectral field in region 1 can be considered as a result of the analysis of the rotated region 2′ from Figure 2b in Section 3, yielding
d d v ψ ˜ Y 1 ( η , v ) = M π γ ( j α o , j η ) · ψ ˜ Y 1 ( η , v ) + ψ s o ( v ) , v < 0 ,
where we note that γ π γ impacts all terms of the solution, as already reported in Section 3, i.e., u iY 1 , v iY 1 , λ i Y 1 , m i Y 1 , and the field components. The solution takes the form
ψ ˜ Y 1 ( η , v ) = i = 3 4 C i e λ γ i Y 1 ( π γ ) v u iY 1 i = 1 2 u iY 1 v iY 1 · v e λ γ i Y 1 ( π γ ) ( v v ) ψ s o ( v ) d v , v < 0 ,
where v = x in Figure 2b is now different from v = X 1 in Figure 2a. From the homogeneous portion of the solution in (59), we obtain the definitions of arbitrary coefficients in terms of the field components at v = 0 (face a):
v iY 1 · ψ ˜ Y 1 ( η , 0 ) = C i , i = 3 , 4 .
The particular integrals in (59) are terms related to face o via ψ s o ( v ) . Due to the linearity of the problem, we again apply the superposition principle and can interpret (59) as the result of an equivalent theorem, where ψ ˜ Y 1 ( η , v ) is represented through equivalent sources at faces a and o.
Using the superposition principle and considering only homogeneous portions of (56) and (59), we can represent the complete field without the particular integrals. Each contribution originating from (56) and (59) is a spectral component that can be inversely Fourier/Laplace transformed into the physical domain ( u , v ) ; these contributions represent the fields from equivalent currents distributed in half-planes (face o and face a, respectively). The application of the asymptotic representation of fields for each component in a unique global system of cylindrical coordinates provides the estimation of the field in terms of the classical GTD for angular region 1, but as a superposition of the GTD for two half-planes (face o and face a), as in Kirchhoff representations. This procedure will be examined in detail in the practical examples reported in the following sections and is a fundamental tool for estimating the GTD directly in the Fourier/Laplace domain for an angular region filled by arbitrarily linear media, where the GTD in the w plane is not available (as is common in isotropic angular region problems). Indeed, the computation of the GTD for an angular region filled by arbitrarily linear media is proposed here by effectively resorting to the computation of the GTD in two half-plane problems.
An alternative way to obtain the far field is based on the computation of the spectral field for any azimuthal direction φ by splitting the angular region into two subregions at any observation angle φ (subregion A 0 < φ < φ and subregion B φ < φ < γ ). After the face spectra at φ = 0 , γ have been obtained for the entire angular region, as proposed in the previous sections, we can then relate the spectra at φ to those of the two faces using the functional equations of the two subregions. These φ -parametric spectral representations of field spectra allow for asymptotic evaluation of the far field at any φ . We observe that the functional equations are written in terms of continuous field components at the boundary faces of the angular region; see Section 4. The properties analyzed in this section can be interpreted as a novel and original form of the electromagnetic equivalence theorem in the spectral domain in the specific context of problems involving angular regions filled by an arbitrary linear medium.

5. Validation of the Novel Regularization Procedure with a Simple Example: Direct Fredholm Factorization Applied to a PEC Wedge in an Isotropic Region

In order to validate the procedure from a mathematical point of view, let us first demonstrate its efficacy in the simple case of angular region 1 (Figure 2a,b), representing a PEC wedge filled by an isotropic medium, where a closed-form WH solution is available. From the following:
m = m i ( π γ ) = m i + 2 ( γ ) = η cos γ + ξ sin γ , i = 1 , 2 ; ξ = k 2 α o 2 η 2
u 1 = τ o 2 ω ε ξ α o η ω ε ξ 0 1 , u 2 = α o η ω ε ξ ( ξ 2 + α o 2 ) ω ε ξ 1 0 , u 3 = τ o 2 ω ε ξ α o η ω ε ξ 0 1 , u 4 = α o η ω ε ξ ( ξ 2 + α o 2 ) ω ε ξ 1 0
v 1 = ξ 2 + α o 2 2 ω μ ξ α 0 η 2 ω μ ξ 0 1 2 , v 2 = α 0 η 2 ω μ ξ k 2 α o 2 2 ω μ ξ 1 2 0 v 3 = ξ 2 + α o 2 2 ω μ ξ α 0 η 2 ω μ ξ 0 1 2 , v 4 = α 0 η 2 ω μ ξ k o 2 α o 2 2 ω μ ξ 1 2 0
we have the functional equations [34] for region 1 (the first two equations in (45)):
α o η E o ρ ( η ) + ( η 2 k 2 ) E o z ( η ) + k ξ Z o H o ρ ( η ) = α o η E a ρ ( m ) [ η ξ sin ( γ ) + cos ( γ ) ( k 2 η 2 ) ] E a z ( m ) + k ξ Z o H a ρ ( m ) sin ( γ ) α o k Z o H a z ( m )
( k 2 α o 2 ) E o ρ ( η ) + α o η E o z ( η ) + k ξ Z o H o z ( η ) = ( k 2 α o 2 ) E a ρ ( m ) + α o [ cos ( γ ) η sin ( γ ) ξ ] E a z ( m ) + k Z o [ sin ( γ ) η + cos ( γ ) ξ ] H a z ( m ) .
At normal incidence ( α o = 0 ), we obtain
ξ E o z ( η ) + k Z o H o ρ ( η ) = [ η sin ( γ ) + ξ cos ( γ ) ] E a z ( m ) + k Z o H a ρ ( m ) ,
k E o ρ ( η ) + ξ Z o H o z ( η ) = k E a ρ ( m ) + Z o [ η sin ( γ ) + ξ cos ( γ ) ] H a z ( m ) ,
where we note the decoupling of Equations (66) and (67) for E z and H z polarization, respectively. The imposition of the PEC boundary conditions on functional Equations (66) and (67) yields the GWHEs
H o ρ ( η ) = H a ρ ( m ) ,
ξ H o z ( η ) = [ η sin ( γ ) + ξ cos ( γ ) ] H a z ( m ) ,
with plus/minus field unknowns, respectively, in η , m . We note that the regularity properties of the problem depend on the multivalued function ξ = k 2 η 2 (due to physical reasons) [11], which defines proper and improper sheets of the η plane.

5.1. Classical Solution of the GWHEs of the Problem in Different Complex Planes

In order to illustrate and validate the new direct Fredholm factorization procedure of Section 3 in the following subsections, this subsection presents the classical WH solution of (68) and (69) obtained in closed form [11] with the help of a specialized mapping, factorization, and decomposition and with the extraction of source terms such as geometrical optics (GO) fields for plane wave illumination. This subsection also clarifies important properties related to different complex planes (including the angular complex plane w) where the problem and the solutions are represented.
The specialized mapping for the closed-form WH solution is
η = k cos γ π arccos η ¯ k ,
which was introduced for the first time in 2001 and has been extensively used in isotropic wedge scattering problems, as reported in [11,12]. The mapping transforms plus unknowns in the η plane and minus unknowns in the m plane (61) into plus and minus unknowns in the η ¯ plane, respectively, yielding the following classical Wiener–Hopf equations in the new complex plane η ¯ :
H o ρ + ( η ¯ ) = H a ρ + ( η ¯ )
ξ H o z + ( η ¯ ) = [ η sin ( γ ) + ξ cos ( γ ) ] H a z + ( η ¯ )
where ξ and η become functions of η ¯ and
m = k cos γ π arccos η ¯ k + γ .
From this point, the solution proceeds as it does for CWHEs, with factorization, decomposition, and the application of Liouville’s theorem while considering plane wave illumination with E z and H z polarization and incident waves.
E z i ( ρ , φ ) = E o e j k ρ cos ( φ φ o ) , H ρ i ( ρ , φ ) = 1 j ω μ ρ E z i ( ρ , φ ) φ = k ω μ sin ( φ φ o ) e j k ρ cos ( φ φ o ) E o
H z i ( ρ , φ ) = H o e j k ρ cos ( φ φ o ) , E ρ i ( ρ , φ ) = 1 j ω ε ρ H z i ( ρ , φ ) φ = k ω ε sin ( φ φ o ) e j k ρ cos ( φ φ o ) H o
Due to the PEC boundary conditions, we obtain the following GO source terms, which are tangential to faces a and o, respectively, of angular region 1.
H x G O ( ρ , 0 ) = 2 E o Z o sin φ o e j k ρ cos φ o , H ρ G O ( ρ , γ ) = 2 E o Z o sin ( γ φ o ) e j k ρ cos ( γ φ o )
H z G O ( ρ , 0 ) = 2 H o e j k ρ cos ( φ o ) , H z G O ( ρ , γ ) = 2 H o e j k ρ cos ( γ φ o )
In the spectral domain, according to the Laplace transforms (33)–(34), (76)–(77) become
H o ρ G O ( η ) = 2 j E o sin φ o Z o ( η η o ) , H a ρ G O ( m ) = 2 j E o sin ( γ φ o ) Z o ( m m o )
H o z G O ( η ) = 2 j H o η η o , H a z G O ( m ) = 2 j H o m m o
with η o = k cos φ o and m o = k cos ( γ φ o ) . In the η ¯ plane (70), pole η o is mapped into η ¯ o = k cos ( π γ φ o ) . In the following, we assume φ o < γ / 2 to locate η ¯ o in the upper half-plane of the complex plane η ¯ , yielding non-standard plus unknowns; generalization is straightforward, yielding η ¯ o in the η ¯ -lower half-plane, while γ / 2 < φ o < γ .
Focusing our attention on E z polarization, due to the simplicity of Equation (71), we observe the absence of a need for factorization; thus, we perform decomposition to highlight the non-standard contribution in the plus unknown H o ρ + ( η ¯ ) , which is constituted by H o ρ G O ( η ) = R / ( η η o ) (78) and can be mapped into the η ¯ plane (70), yielding H o ρ G O ( η ¯ ) = T / ( η η o ) . We obtain
H o ρ + ( η ¯ ) T η ¯ η ¯ o = H a ρ + ( η ¯ ) T η ¯ η ¯ o
with
T = R d η ¯ d η η o = 2 j π γ E o Z o sin π γ φ o , R = 2 j E o sin φ o Z o , d η ¯ d η η o = π γ sin π γ φ o sin φ o .
Due to the regularity and asymptotic behavior of the LHS and RHS of (80), applying Liouville’s theorem, (80) is equal to zero; thus, we obtain the following simple closed-form solutions:
H o ρ + ( η ¯ ) = T η ¯ η ¯ o , H a ρ + ( η ¯ ) = T η ¯ + η ¯ o .
The solutions in (82) can be mapped into the η plane using the inverse mapping of (70):
η ¯ = k cos π γ arccos η k .
We recall that the regularity properties of the problem in (68)–(69) in the η plane depend on the multivalued function ξ = k 2 η 2 (due to physical reasons), and now, after the application of the mapping in (70), on the multivalued function κ = k 2 η ¯ 2 in the η ¯ plane through the log representation of arccos η ¯ / k ; see Section 3.4 of [11]. Contrary to (70), the transformation in (83) requires particular attention, as it maps η ¯ into η for 0 < γ < π without covering the entire proper sheet of the η plane defined by the ξ function. For this reason, a portion of the proper sheet of the η plane falls into the improper sheet of the η ¯ plane (defined by κ ), and because the closed-form solution is obtained in the η ¯ plane, this solution must be considered correct (not offending) only in the proper sheet of η ¯ even after applying the transformation in (83). To easily control the proper/improper sheets of the η and η ¯ planes, we can resort to their visualization in the complex plane w ( η = k cos w ; thus, η ¯ = k cos π γ w and m = k cos ( w + γ ) ). The w plane shows the proper sheets of both planes ( η , η ¯ ) in a unique plane. In particular, for real w, the proper segments originating from η and η ¯ (related to ξ and κ , respectively) are π < w < 0 and γ < w < 0 , respectively; see Section 3.4 of [11]. This means that the closed-form solution obtained in the proper sheet of η ¯ is not valid in the entire proper sheet of the η plane but only in a portion, which is due to the properties of (83).
Let us now consider the CWHE of H z polarization (72)
G ( η ¯ ) H o z + ( η ¯ ) = H a z + ( η ¯ ) , G ( η ¯ ) = ξ / n ,
with n = η sin ( γ ) ξ cos ( γ ) = k 2 m 2 . According to [33], we have the factorization
G ( η ¯ ) = G ( η ¯ ) G + ( η ¯ ) , G + ( η ¯ ) = ξ ξ n + , ξ = k η ¯ , n + = k + η ¯ .
Confirming the same assumption that φ o < γ / 2 for simplicity, η ¯ o is located in the η ¯ upper half-plane, yielding a nonstandard plus unknown H o z + ( η ¯ ) that is constituted by the source nonstandard component H o z G O ( η ) = R H / ( η η o ) (79), which in the η ¯ plane becomes
H o z G O ( η ¯ ) = T H η ¯ η ¯ o , T H = R H d η ¯ d η η o = 2 j H o π γ sin π γ φ o sin φ o , R H = 2 j H o .
Applying factorization and decomposition to (84), we obtain
G + ( η ¯ ) H o z + ( η ¯ ) G + ( η ¯ o ) H o z G O ( η ¯ ) = G 1 ( η ¯ ) H a z + ( η ¯ ) G + ( η ¯ o ) H o z G O ( η ¯ ) .
Due to the regularity and asymptotic behavior of the LHS and RHS of (87), applying Liouville’s theorem, (87) is equal to zero; thus, we obtain the following simple closed-form solutions:
H o z + ( η ¯ ) = G + 1 ( η ¯ ) G + ( η ¯ o ) H o z G O ( η ¯ ) , H a z + ( η ¯ ) = G ( η ¯ ) G + ( η ¯ o ) H o z G O ( η ¯ ) .
Again, the closed-form solutions in (88) with the H z polarization obtained in the proper sheet of the η ¯ plane can be mapped into the η plane using the inverse mapping in (83); however, we need to consider these solutions valid only for η values belonging to the proper sheet of the η ¯ plane. Moreover, this property can be ascertained by checking that (68)–(69) (from which we obtained the solutions in η ¯ ) are enforced only for η values belonging to the proper sheet of the η ¯ plane.
In order to obtain solutions that are valid in the entire proper sheet of the η plane or beyond (i.e., also in the improper sheet), we need to resort to an analytical continuation technique that, in case of a unique propagation constant problems, can be implemented via representation of GWHEs in the complex w plane yielding difference equations (for instance, see the examples in [11,12]). Another option is to describe the problem with a unique propagation constant directly in the w plane, where the concept of proper and improper sheets of the η and η ¯ planes is expanded periodically into the w plane with an alternative vision of Riemann sheets. In this case, the closed-form solutions corresponding to (88) and (89) are valid in the entire w plane, as opposed to approximate solutions obtained with numerical line integration located in a particular sheet in either the η ¯ or w plane. In this last case, which originates from the classical implementation of Fredholm factorization [12], we again need to resort to difference equations for analytical continuation.
H o z + ( w ) = 2 j H o π csc w sin π w γ k γ cos π w γ + k γ cos π φ o γ , H a z + ( w ) = 2 j H o π csc w sin π w γ k γ cos π w γ + k γ cos π φ o γ

5.2. Regularized Integral Equation Method for the Direct Solution of the GWHEs in Angular Regions (Direct Fredholm Factorization)

Following the procedure in Section 3, which is simplified in the current problem because of the isotropic medium, we duplicate the equations. For E z polarization, we have
H o ρ ( η ) = H a ρ ( m ) , H a ρ ( η ) = H o ρ ( m ) ,
while for H z polarization, we have
ξ H o z ( η ) = [ η sin ( γ ) + ξ cos ( γ ) ] H a z ( m ) ξ H a z ( η ) = [ η sin ( γ ) + ξ cos ( γ ) ] H o z ( m )
where m = m ( η ) is defined in (61). Notice that when applying (70) to each of (90)–(91), the duplicated equations assume the same form of the original ones (both CWHEs), with only the replacement of η ¯ with η ¯ .
Both systems of equations can be considered a particular case of
G ( η ) F + ( η ) = H ( η ) X + ( m ) G a ( η ) X + ( η ) = H a ( η ) F + ( m )
which is suitable for describing more general cases. To describe the procedure, for simplicity, let us assume that F + ( η ) is a non-standard plus η unknown while X + ( m ) is a standard minus m unknown; generalization is possible with a little effort.
Applying the Cauchy decomposition Formula (49) to the unknowns defined in m ( η ) ,
F + ( m ) = 1 2 π j F + ( η ) η + m d η + F + n . s . ( m ) , η R X + ( m ) = 1 2 π j X + ( η ) η + m d η , η R .
From (92), we obtain the following system of integral equations:
G ( η ) F + ( η ) = 1 2 π j H ( η ) X + ( η ) η + m ( η ) d η G a ( η ) X + ( η ) = 1 2 π j H a ( η ) F + ( η ) η + m ( η ) d η + H a ( η ) F + n . s . ( m ( η ) ) .
This is not a system of Fredholm integral equations of the second kind (non-compact kernel). To regularize (94), we follow the procedure presented in Section 3. Performing a smile integration of (94), after mathematical manipulation, on the LHS, we have
1 2 π j γ 1 t G ( t ) F + ( t ) t η d t = G ( η ) F + ( η ) + 1 2 π j ( G ( t ) G ( η ) ) F + ( t ) t η d t G ( η ) F + n s ( η ) 1 2 π j γ 1 t G a ( t ) X + ( t ) t η d t = G a ( η ) X + ( η ) + 1 2 π j ( G a ( t ) G a ( η ) ) X + ( t ) t η d t
and on the RHS we have
1 2 π j γ 1 t 1 2 π j H ( t ) t η X + ( η ) η + m ( t ) d η d t = 1 ( 2 π j ) 2 M ( η , η ) X + ( η ) d η
and
1 2 π j γ 1 t 1 2 π j H a ( t ) t η F + ( η ) η + m ( t ) d η d t = 1 ( 2 π j ) 2 M a ( η , η ) F + ( η ) d η 1 2 π j γ 1 t H a ( t ) t η F + n s ( m ( t ) ) d t = 1 2 π j [ H a ( t ) H a ( η ) ] F + n s ( m ( t ) ) t η d t + H a ( η ) 2 π j γ 1 t F + n s ( m ( t ) ) t η d t
where
M ( η , η ) = γ 1 t H ( t ) ( t η ) ( η + m ( t ) ) d t = H ( t ) H ( η ) ( t η ) ( η + m ( t ) ) d t + H ( η ) γ 1 t 1 ( t η ) ( η + m ( t ) ) d t M a ( η , η ) = γ 1 t H a ( t ) ( t η ) ( η + m ( t ) ) d t = H a ( t ) H a ( η ) ( t η ) ( η + m ( t ) ) d t + H a ( η ) γ 1 η 1 ( t η ) ( η + m ( t ) ) d t .
Merging (95) and (96)–(98), we obtain the following FIEs of the second kind:
G ( η ) F + ( η ) + 1 2 π j ( G ( t ) G ( η ) ) F + ( t ) t η d t = 1 ( 2 π j ) 2 M ( η , η ) X + ( η ) d η + G ( η ) ) F + n s ( η )
G a ( η ) X + ( η ) + 1 2 π j ( G a ( t ) G a ( η ) ) X + ( t ) t η d t = 1 ( 2 π j ) 2 M a ( η , η ) F + ( η ) d η + 1 2 π j ( H a ( t ) H a ( η ) ) F + n s ( m ( t ) ) t η d t + H a ( η ) 2 π j γ 1 t F + n s ( m ( t ) ) t η d t .
It can be observed that the regularized FIEs (99) and (100) are particularly efficient from a computational point of view due to the presence of compact kernels integrated along the real axis, with the exception of the smile integration included in (98):
M e ( η , η ) = γ 1 t 1 ( t η ) ( η + m ( t ) ) d t .
The evaluation of the integral in (101) can be effectively performed by warping the smile contour γ 1 t in the lower half of the complex t plane into the integration path λ 1 wrapped around the vertical branch cut Γ 1 of m ( t ) (61), which originated in branch point + k (see Figure 5).
By collapsing λ 1 onto Γ 1 , we obtain
M e ( η , η ) = Γ 1 Δ ( 1 ( t η ) ( η + m ( t ) ) d t ,
where
Δ ( 1 ( t η ) ( η + m ( t ) ) = 4 ( k t ) ( k + t ) sin ( γ ) ( t η ) k 2 + 2 t 2 + η 2 4 t η cos ( γ ) + k 2 cos ( 2 γ ) .
Assuming that t = k j v ( v > 0 ) , the representation in (102) is quickly numerically convergent. A closed-form expression of (102) is obtainable by considering the following:
  • Selection of the branch cut Γ 1 as the line t = k u (with real u > 1 ) replacing λ 1 and use of the mapping t = k u under the integration sign.
  • Expansion of (103) with a minimal denominator.
  • Careful mathematical manipulation of multivalued functions.
We obtain
M e ( η , η ) = 2 sin ( γ ) F ( u 1 ( η ) ) ( u 1 ( η ) u 2 ( η ) ) ( u 1 ( η ) u 3 ( η ) ) F ( u 2 ( η ) ) ( u 1 ( η ) u 2 ( η ) ) ( u 2 ( η ) u 3 ( η ) ) + F ( u 3 ( η ) ) ( u 1 ( η ) u 3 ( η ) ) ( u 2 ( η ) u 3 ( η ) )
with
F ( u ) = j u log ( 2 ) 1 u 2 log ( u + j 1 u 2 )
and the poles
u 1 ( η ) = η / k , u 2 ( η ) = η cos γ k 2 η 2 sin γ k , u 3 ( η ) = η cos γ + k 2 η 2 sin γ k .
Let us now go back to particular cases and consider equations for H z polarization (91) in angular region 1 representing a PEC wedge, written in the form (92) with the following definitions:
F + ( η ) = H o z ( η ) , X + ( η ) = H a z ( η ) , G ( η ) = G a ( η ) = ξ η sin ( γ ) + ξ cos ( γ ) , H ( η ) = H a ( η ) = 1 .
The set of FIEs (99)–(100) is simplified; in particular, M ( η , η ) = M a ( η , η ) = M e ( η , η ) , as reported in (102). Equations (99)–(100) respectively become the following system of FIEs:
G ( η ) H o z ( η ) + 1 2 π j ( G ( t ) G ( η ) ) H o z ( t ) t η d t = 1 ( 2 π j ) 2 M e ( η , η ) H a z ( η ) d η + s 1 ( η )
and
G ( η ) H a z ( η ) + 1 2 π j ( G ( t ) G ( η ) ) H a z ( t ) t η d t = 1 ( 2 π j ) 2 M e ( η , η ) H o z ( η ) d η + s 2 ( η )
with
s 1 ( η ) = G ( η ) H o z n s ( η ) , s 2 ( η ) = 1 2 π j γ 1 t H o z n s ( m ( t ) ) t η d t .
Let us focus our attention on the source term (110); for simplicity, we assume that only F + ( η ) = H o z ( η ) is nonstandard:
F + n s ( η ) = H o z n s ( η ) = 2 j H o η η o
where η o = k cos ( φ o ) , 0 < φ o < π / 2 and k has small losses ( k = k r j k i , k i < < k r ). From (111), according to the properties of m ( η ) (see also Figure 4), H o z n s ( m ( η ) ) in the proper lower half of the complex η plane shows poles originating from the zeros of m ( η ) + η o (in the m plane, we have the pole m o = η o ). The poles can be related to GO waves, i.e., connected to the last couple of reflections from faces a and o that create shadow boundaries (for instance, see [40]). As an example, if φ o < π γ , then we have one reflection from face a and one reflection from face o reflected again by face a. In fact, from a mathematical point of view, in this case we have pole m o related to poles η r a = k cos ( γ φ o ) (reflection from face a) and η r a r o = k cos ( γ + φ o ) (reflection from a after o), which are associated with incoming azimuthal directions γ φ o with respect to the reference face a, i.e., the incoming directions 2 γ φ o with respect to face o. However, we also need to note that the residues of the poles in the selected test problems are always related only to the incident field. This means that the primary spectra of H o z n s ( m ( η ) ) in (110) are more similar to a replica of the incident spectrum for η r a , η r a r o , similar to what was described in the w plane in [13].
Indeed, the integrand of the source term in (110) also exhibits a branch cut of m ( η ) . Thus, we estimate (110) by warping γ 1 t into λ 1 , as follows:
s 2 ( η ) = 1 2 π j λ 1 H o z n s ( m ( t ) ) t η d t + R a η η r a + R a o η η r a r o
where R a and R o are the residues of H o z ( m ( η ) ) in η r a and η r a r o , respectively.
H o z n s ( m ( η ) ) = 2 j H o m ( η ) + η o = T m o m ( η ) + η o , T m o = 2 j H o
R a , a o = T m o d η d m η r a , η r a r o = 2 j H o cos γ + η sin γ k 2 η 2 η r a , η r a r o
Using the same mathematical procedure that yields (102)–(104) from (101), from (112) we obtain
1 2 π j λ 1 H o z n s ( m ( t ) ) t η d t = H o π λ 1 1 ( t η ) ( m ( t ) η o ) d t = H o π M e ( η , η o ) .
Thus,
s 2 ( η ) = H o π M e ( η , η o ) + R a η η r a + R a o η η r a r o .
The final set of FIEs for H z polarization when the structure is illuminated by a plane wave with 0 < φ o < π / 2 is then (108)–(109) (a specialization of (99)–(100)), with sources s 1 , 2 ( η ) being defined and calculated in (110), (111), and (116). Note that s 1 ( η ) and s 2 ( η ) are spectral components defined in the η plane of H o z ( η ) and H a z ( η ) , respectively, i.e., with the reference coordinate system of face o and face a.
We now examine the convergence properties of the FIEs in (108)–(109) to obtain accurate numerical results [41]. According to the classical Fredholm factorization method [11], the regularization procedure provides compact kernels of the type reported on the LHS of (108)–(109), i.e., it is square-integrable. The further integral operator reported on the RHS of (108)–(109) in terms of M e ( η , η ) is related to the coupling term between the spectra of delimiting faces. This kernel is again compact, as (101) shows that M e ( η , η ) is never singular for η t and η m ( t ) and (104) shows that M e ( η , η ) is square-integrable according to its asymptotic behavior in terms of (106). Similar considerations can be repeated in more complex and general cases of angular regions immersed in/made of arbitrary linear media.

5.3. Implementation of a Numerical Example and Validation of Direct Fredholm Factorization

Let us consider region 1 in Figure 1 with an aperture angle π / 2 < γ < π , which is filled by a homogeneous isotropic medium with a propagation constant k ( k = k r j k i , k i < < k r ) and terminated by PEC boundary conditions. The angular region is illuminated by an H z polarized plane wave with an incoming direction φ o ( 0 < φ o < π γ ) and intensity H o . The spectral solution ( H o z ( η ) , H a z ( η ) ) can be provided by the system of FIEs in (108)–(109). Due to the convergence properties of the kernel [41], simple sample-and-hold approximation is enforced with the truncation of integration intervals at ± A and integration step h such that A / h N . We tested our novel direct FIE solution against the classical exact closed-form solution provided in Section 5.1 in the η and w planes, respectively ((88) and (89)). Furthermore, we compared the asymptotic results in terms of GTD coefficients. We examine the case where γ = 0.7 π , k = 1 j 0.1 , H o = 1 A / m , and φ o = 0.1 π in detail. Because we have 0 < φ o < π γ , the GO field is constituted by the incident waves, waves reflected by face a, and doubly-reflected waves (from face o and then from face a); only the plus spectral unknown along face o, i.e., H o z ( η ) , is nonstandard in the WH formulation (91), as reported in the example in the previous subsection. To enhance the convergence of the approximate FIE solution provided by (108)–(109), we warp the integration line constituted by the real axis into a straight line located in the first and third quadrants of the complex plane at an angle θ with respect to the real axis (the singularities of the kernel and the sources are located in the second and fourth quadrants; see Sections 5.3 and 5.4 of [12]):
α t ( t ) = t exp j θ , t R , 0 < θ < π / 2 .
According to the physical parameters of the test problem, we have η r a = 0.309017 0.0309017 j , η r a r o = 0.809017 0.0809017 j , with both located in the lower half-plane, as considered in (116). The discretization of (108)–(109) by the sample-and-hold technique with A , h yields a linear system of equations with dimensions of 2 ( 2 ( A / h ) + 1 ) :
d ̲ ̲ G + K ̲ ̲ G M ̲ ̲ e M ̲ ̲ e d ̲ ̲ G + K ̲ ̲ G H ̲ o z H ̲ a z = s ̲ 1 s ̲ 2
where the diagonal matrix d ̲ ̲ G , the matrix K ̲ ̲ G , and the matrix M ̲ ̲ e contain samples of G ( η ) , G ( t ) G ( η ) t η , and M e ( η , η ) , respectively, while the vectors H ̲ o z , H ̲ a z and s ̲ 1 , s ̲ 2 contain samples of H o z ( η ) , H a z ( η ) and s 1 ( η ) , s 2 ( η ) , respectively. Note that M ̲ ̲ e is a coupling matrix that is much weaker than the remaining terms. The sampled solutions allow for a representation of H o z ( η ) , H a z ( η ) to be built by substituting them into the integral part of (108)–(109):
H o z , a z ( η ) = h 2 π j A / h A / h [ G 1 ( η ) G ( α t ( h i ) ) 1 ] H o z , a z ( α t ( h i ) ) α t ( h i ) η + h G 1 ( η ) ( 2 π j ) 2 A / h A / h M e ( η , α t ( h i ) ) H a z , o z ( α t ( h i ) ) + G 1 ( η ) s 1 , 2 ( η ) .
These approximate expressions of H o z ( η ) , H a z ( η ) are valid for analytic continuation in the proper sheet of the η plane, which is useful for correctly estimating fields in the physical domain through the asymptotics of half-planes, as discussed in Section 4. This property limits the requirement to know the spectra only in the proper sheet, as acquired in the procedure, which is a novelty and represents progress with respect to classical Fredholm factorization combined with spectral mapping in GWHE wedge problems.
To highlight the performance of this method, we compare the spectra along the real axis of the η plane and a segment of the η plane that is useful for asymptotics according to the steepest descent path (SDP) method, which in an isotropic medium corresponds to η = k cos w with π < w < 0 , i.e., the segment that connects k with k.
To study the convergence of the method, we select physical parameters of region 1 with an aperture angle of γ = 0.7 π and plane wave illumination at H z polarization with H o = 1 A / m , φ o = 0.1 π , and k = 1 j 0.1 . We select the quadrature parameters 5 A 40 , 0.2 h 0.25 , θ = 0.1 such that A / h N . The numerical results are provided in Figure 6 along the segment for asymptotic estimation. In the figure, we note that along the segment we have a degradation of the spectral solution near w = π , 0 , which corresponds to the branch points η = k , k of ξ (where ξ defines the proper and improper sheet of the η plane [11]). We recall that the solution of the FIEs was obtained through the simple sample-and-hold quadrature and estimation of M e ( η , η ) , which saturates the precision, particularly near η = k , k (the branch point η = k is a local offending singularity for the plus spectra that should not appear, while η = k is related to physical structural properties of the problem). An improvement would be obtained with a specialized quadrature (and method of moments) capable of taking into account non-algebraic behavior such as branch points (see [42]); however, the scope of the present method involves obtaining a very simple, fast, and convergent solution that cannot incorporate sophisticated quadratures. Furthermore, we observe that the lack of precision near η = k , k is mitigated while computing asymptotics because plus spectral unknowns are multiplied by sin w , providing locally-smoothing errors. However, while the offending η = k is a very local perturbation, the physical η = k is more present, as it is physical.
To recover the quality of the solution near w = π , 0 ( η = k , k ) , we resort to spectral considerations based on the properties of the original GWHE formulation in (91). Equation (91) can be applied to the approximate solutions obtained from the FIEs to obtain a new improved representation of the plus spectra. This application allows us to obtain an approximate spectrum near w = π , 0 ( η = k , k ) that originates from another portion of the η plane according to m ( η ) . This procedure is particularly effective and valid because m ( η ) for η in the proper sheet lies in just a portion of the proper sheet of the η plane. This property is particularly evident (in this simplified isotropic problem) by rewriting (91) in the w plane:
H o z ( k cos w ) = n ( k cos w ) ξ ( k cos w ) H a z ( k cos ( w + γ ) ) H a z ( k cos w ) = n ( k cos w ) ξ ( k cos w ) H o z ( k cos ( w + γ ) )
with ξ ( k cos w ) = k sin w and n = k sin ( w + γ ) . We note that π < w < 0 on the LHS corresponds to π + γ < w < γ on the RHS due to ( m = k cos ( w + γ ) ) , where the unknowns are correctly computed. This methodology (named iteration) reimposes GWHEs on the initial FIE-approximated spectra and shifts the lack of precision to a region where the spectral solution is good, yielding a homogenization of the error level (see Figure 7). In the figure, we report (1) the exact solution, (2) the approximate solution obtained from the quadrature of FIEs with A = 40 , h = 0.025 , (3) the approximate solution obtained from the quadrature of FIEs with A = 40 , h = 0.025 plus the application of (91), and (4) the approximate solution obtained from the application of (91) to the sources of FIEs while ignoring integral terms, i.e., using H o z , a z ( η ) = + G 1 ( η ) s 1 , 2 ( η ) .
Note that when considering (91), the map in (120) is only limited; thus, we cannot interpret this procedure as a first iteration of the application of the contraction theorem. In fact, according to our studies, successive iterations do not yield any benefit in the convergence of the solution. This is also justified by the fact that in the w plane, the multiple applications of (120) correspond to recursive equations/difference equations that further shift the spectra in the w plane, thereby navigating replicas of the proper and improper sheets of the η plane defined by ξ ; see [11,13]. Moreover, the map is also excluded from compensating for all physical behaviors of the problem starting from roughly approximate solutions. Figure 7 shows the importance of the quality of starting spectra that originate from the solution of the FIEs before the application of (120).
Finally, we observe that while the FIEs provide good spectra except near the branch cuts, the iteration in (91) enforces correct modeling of structural spectral properties such as branch cuts.
To further compare the solutions and validate the proposed procedure, we have computed the GTD diffraction coefficient as outlined in Section 4 using asymptoptics. Using superposition, we can compute the diffraction by applying asymptotics individually to the spectral solutions at faces o and a while considering only homogeneous terms in (56) and (59). This takes care of the different reference coordinates (see the discussion in Section 3 and Section 4 considering region 1 characterized by γ as region 2′ characterized by π γ ; see Figure 1 and Figure 2, (45), and (59)):
ψ ˜ y ho ( η , v ) = i = 1 2 v i · ψ ˜ y ( η , 0 ) e λ γ i ( γ ) v u i , v > 0 ,
ψ ˜ Y 1 ho ( η , v ) = i = 3 4 v iY 1 · ψ ˜ Y 1 ( η , 0 ) e λ γ i Y 1 ( π γ ) v u iY 1 , v < 0 .
Starting from the inversion of the contribution of face o (121), we have
ψ y ho ( u , v ) = 1 2 π B r ψ ˜ y ho ( η , v ) e j η u d η .
According to the coordinate mapping in (18), from (16) and (23) we have
λ γ i ( γ ) v j η u = j η cos γ v j ξ i sin γ v j η ( x v cos γ ) = j ( η x + ξ i y ) , i = 1 , 2
with ξ i = ξ , i = 1 , 2 ; thus,
ψ y ho ( x , y ) = 1 2 π B r i = 1 2 v i · ψ ˜ y ( η , 0 ) u i e j ( η x + ξ y ) d η ,
where B r is the Bromwich contour (over all singularities). The asymptotic estimation of (125) at the far field is composed of GO terms (captured poles) and GTD diffracted components (due to the saddle point with the application of the SDP method); in global cylindrical coordinates, this is
ψ y ho , gtd ( ρ , φ ) = k 2 π ρ e j ( k ρ π / 4 ) i = 1 2 v i · ψ ˜ y ( k cos φ , 0 ) u i sin | φ | .
In our test problem (region 1 with PEC faces with H z polarization), this reduces to the third component
ψ y ho , gtd ( ρ , φ ) [ 3 ] = H o z g t d ( ρ , φ ) = k 2 π ρ e j ( k ρ π / 4 ) H o z ( k cos φ ) 2 sin | φ |
according to the definition of ψ ˜ y ( η , v ) in (3) and u i , v i in (62)–(63). We obtain the following GTD diffraction coefficient component due to face o:
D H o z g t d ( φ ) = k H o z ( k cos φ ) sin | φ | j 2 H o .
Now, we repeat the procedure starting from the inversion of the contribution of face a (122) using the notation in Figure 1 and Figure 2b:
ψ Y 1 ho ( u , v ) = 1 2 π B r ψ ˜ Y 1 ho ( η , v ) e j η u d η .
According to the coordinate mapping
X 1 = u + v cos ( π γ ) , Y 1 = v sin ( π γ ) ,
from (16) and (23) we have
λ γ i Y 1 ( π γ ) v j η u = + j η cos γ v + j ξ i sin γ v j η ( X 1 + v cos γ ) = j η X 1 + j ξ i Y 1 , i = 3 , 4
where ξ i = ξ , i = 3 , 4 ; thus,
ψ Y 1 ho ( X 1 , Y 1 ) = 1 2 π B r i = 3 4 v iY 1 · ψ ˜ Y 1 ( η , 0 ) u iY 1 e j η X 1 + j ξ Y 1 d η ,
where B r is a Bromwich contour. The asymptotic estimation of (132) at the far field is composed of GO terms and GTD diffracted components in the following global cylindrical coordinates:
ψ Y 1 ho , gtd ( ρ , φ ) = k 2 π ρ e j ( k ρ π / 4 ) i = 3 4 v iY 1 · ψ ˜ Y 1 ( k cos ( φ γ ) , 0 ) u iY 1 sin | φ γ | .
For our test problem (region 1 with PEC faces with H z polarization), this reduces to the third component
ψ Y 1 ho , gtd ( ρ , φ ) [ 3 ] = H a z g t d ( ρ , φ ) = k 2 π ρ e j ( k ρ π / 4 ) H a z ( k cos ( φ γ ) ) 2 sin | φ γ |
according to the definition of ψ ˜ y ( η , v ) in (3) and u iY 1 = u i , v iY 1 = v i in (62)–(63). Note that the invariance of u iY 1 = u i , v iY 1 = v i in the rotation of the reference system is allowable only in isotropic regions; otherwise, a more complex procedure is required for their definition in arbitrary linear media (see Section 2).
Finally, we obtain the following GTD diffraction coefficient component due to face a:
D H a z g t d ( φ ) = k H a z ( k cos ( φ γ ) ) sin | φ γ | j 2 H o .
The complete GTD coefficient is just the sum for the superposition of (128) and (135):
D H z g t d ( φ ) = D H o z g t d ( φ ) + D H a z g t d ( φ ) .
Figure 8 shows the GTD diffraction coefficient for the test problem under consideration, that is, γ = 0.7 π and plane wave illumination with H z polarization with H o = 1 A / m , φ o = 0.1 π , and k = 1 j 0.1 . The figure reports the exact GTD coefficient in terms of the absolute value and phase together with those obtained with the FIE-approximated estimation of the spectra without and with the application of an iteration while selecting A = 20 , h = 0.05 , θ = 0.1 . Figure 9 shows the corresponding relative error of the GTD diffraction coefficient on a log 10 scale.
We note, as expected, that the solution with the iteration is correct, while the solution without the iteration is lacking in its estimation near the faces of the angular regions, i.e., face o for φ = 0 and face a for φ = γ , as they are related to the spectra of H o z ( η ) near η = k (128) and H a z ( η ) near η = k (135), respectively ( η = k corresponds to w = π , and is a physical branch point). In addition, note that the spectra of H o z ( η ) , H a z ( η ) near η = k ( w = 0 ) are not used for the GTD computation; thus, the lack of possible precision in the offending branch point does not impact the quality of the solution. Moreover, the change in the slope and level of the relative error in Figure 9 is obtained using the reported algorithm to improve the quality of the approximate solution provided by the direct application of Fredholm factorization. In fact, FIE+iteration involves computation of the GTD diffraction coefficient in (136) via (128) and (135), where the spectra H o z ( k cos w ) and H a z ( k cos w ) are obtained by enforcing (120) on the approximate spectra obtained through the direct application of Fredholm factorization. This procedure mixes the spectral resolution properties of the two faces, thereby improving the quality of the spectra and recovering the degradation of the spectral resolution near w = π , 0 , i.e., η = k , k .
Finally, we comment that the direct implementation of FIEs in the w plane yields high-precision results in isotropic angular region problems [12], exceeding the precision of the current procedure in terms of spectra; however, we recall that the scope of the present work involves presenting an effective procedure for computing diffraction that is implementable in problems where the w plane cannot be defined, such as in arbitrarily linear media.

6. An Example of the Application of the Functional Equations in Complex Media: Scattering from a PEC Half-Plane in a Gyrotropic Medium

The scattering of an electromagnetic plane wave at normal incidence by a perfectly conducting semi-infinite screen embedded in a homogeneous gyrotropic medium (such as plasma) is presented in this section with the goal of validating the proposed method, the functional equations, and the WH equations in nonisotropic media. As our formulation is in terms of field components, we selected [15,16,20,21] as studies for comparison; we did not consider other literature where definitions in terms of potentials were employed. In particular, we selected [15], where the distinguished axis of the electric gyrotropic medium is parallel to the edge of the half-plane, i.e., as in plasma, with a uniform magnetic field impressed along the edge direction. This medium enforces a tensorial electric permittivity, which in our reference system of coordinates ( z , x , y ) is
ε ̲ ̲ = ε 3 0 0 0 ε 1 + j ε 2 0 j ε 2 ε 1 ,
where z is the distinguished axis of the medium and μ ̲ ̲ = μ o I ̲ ̲ , ξ ̲ ̲ = ζ ̲ ̲ = 0 ̲ ̲ . As reported in [15], this vector problem is separable into two equivalent scalar problems for the E z (H-mode) and H z (E-mode) polarizations.
By applying the procedure described in Section 2 and with simplified definitions of the quantities reported in Appendix A, we obtain the following (progressive, regressive) eigenvalues:
λ 1 , 3 = ± η 2 k 1 2 = ± j ξ 1 , λ 2 , 4 = ± η 2 k 2 2 = ± j ξ 2
where k 1 2 = ω 2 μ o ε / ε 1 = k o 2 ε / ε 1 , k 2 2 = ω 2 μ o ε 3 = k o 2 ε r 3 , ε r i = ε i / ε o , ε = ε 1 2 ε 2 2 , and k o = ω ε o μ o .
The corresponding eigenvectors u i , from which we can easily compute the reciprocal vectors v i through inversion, are
u 1 = 0 j ( ε 2 η + ε 1 j ξ 1 ) ε ω 1 0 , u 2 = μ o ω ε 2 0 0 1 , u 3 = 0 j ( ε 2 η + ε 1 j ξ 1 ) ε ω 1 0 , u 4 = μ o ω ε 2 0 0 1 .
The problem shows simplification because of γ = π ; see for instance the impact of the anisotropies on (35) or
m = m i ( π γ ) = m i + 2 ( γ ) = η ; i = 1 , 2 .
However, we keep the procedure as general as possible, i.e., applicable to wedge problems (arbitrary γ ); from (27) and (31), we obtain
ψ ˜ s a + ( m i ( γ ) ) = E a z cos ( γ ) , E a ρ + η H a z sin ( γ ) ω ε 1 , H a z cos ( γ ) j H a z ε 2 sin ( γ ) ε 1 , H a ρ E a z η sin ( γ ) μ o ω t .
From here on, we omit the spectral dependence in the field components for better compactness of the formulae. Applying (32), we obtain the following two functional equations for region 1 in explicit form:
E o x ω ε + H o z ξ 1 ε 1 + j η H o z ε 2 = H a z sin ( γ ) ( η ε 1 j ξ 1 ε 2 ) + cos ( γ ) ( ξ 1 ε 1 + j η ε 2 ) + E a ρ ε ω ,
H o x μ o ω E o z ξ 2 = H a ρ μ o ω .
Similarly, the procedure can be repeated for region 2. The complete set of equations highlights the decoupling of E z from H z polarization. Applying the PEC boundary conditions on the faces, after some manipulations we obtain the following equations for the E z and H z polarizations:
H o x = H a ρ 2 + H b ρ 2 E o z ξ 2 μ o ω = H a ρ 2 H b ρ 2
E o x ω ε + H o z ξ 1 ε 1 + j η H o z ε 2 = H a z sin ( γ ) ( η ε 1 j ξ 1 ε 2 ) + cos ( γ ) ( ξ 1 ε 1 + j η ε 2 ) E o x ω ε + H o z ξ 1 ε 1 j η H o z ε 2 = H b z sin ( γ ) ( η ε 1 + j ξ 1 ε 2 ) + cos ( γ ) ( ξ 1 ε 1 j η ε 2 ) .
Now, we impose γ = π , i.e., the angular regions are defined for the half-plane problem. In (144), we note that E z polarization behaves as a half-plane problem immersed in a classical isotropic region [33] but with a propagation constant of k 2 2 = ω 2 μ o ε 3 = k o 2 ε r 3 , i.e., this can be interpreted via a network representation with characteristic impedance Z E z = ω μ o / ξ 2 , confirming the findings of [15].
With further mathematical manipulation of (145), we obtain
2 H o z + 2 i E o x η ω ε ε 2 ξ 1 2 ε 1 2 + η 2 ε 2 2 = H a z + H b z 2 E o x ξ 1 ω ε ε 1 ξ 1 2 ε 1 2 + η 2 ε 2 2 = H a z H b z
and the second equation in (146) shows the same WH kernel as that of Equation (25) in [15]:
G H z 1 = ξ 1 2 ε 1 2 + η 2 ε 2 2 2 ξ 1 ω ε ε 1 = ( k 1 2 η 2 ) ε 1 2 + η 2 ε 2 2 2 k 1 2 η 2 ω ε ε 1 = k 1 2 ε 1 2 / ε η 2 2 k 1 2 η 2 ω ε 1
with the sole irrelevant distinction of a multiplication by a scalar. Moreover, the characteristic pole of the surface wave phenomenon is easily recognizable from the numerator, as also found in [15]. The solutions to the problem can be achieved with the proposed approximate technique that was validated in previous sections or via the classical WH procedure, as in [15]; however, this is beyond the scope of the present study.

7. Conclusions

Spectral methods (such as SM, KL, and WH) are well-consolidated, fundamental, and effective tools for the correct and precise analysis of electromagnetic diffraction problems with one propagation constant, although they are not immediately applicable to problems with multiple propagation constants.
In this study, we propose a comprehensive theoretical package in the spectral domain with all necessary mathematical tools that for the first time extends the possibilities of spectral analysis to electromagnetic problems involving wedges immersed in an arbitrary linear medium, and is additionally extendable to multiple penetrable angular regions. The theory is presented in an exhaustive way, showing the theoretical background, implementation, and validation. The methodology is based on transverse equations for layered angular structures, the characteristic Green function procedure, the Wiener–Hopf technique, and our novel direct Fredholm factorization method, which reduces GWHEs with multiple propagation constants to integral representations in a unique complex plane. Validation through examples is applied, starting from demonstrating the effectiveness of direct Fredholm factorization applied to GWHEs in the scattering from a PEC wedge immersed in an isotropic medium and ending with the validation of functional equations of angular regions in arbitrary linear media with the analysis of a PEC half-plane immersed in particular anisotropic media. While numerically implementing the method, we observed that one of the main difficulties resides in the correct estimation of kernel functions for the presence of multivalued functions that need particular attention in their definition and calculation.
The proposed equations are interpreted using a network formalism in order to provide a systematic perspective, particularly for the analysis of complex scattering problems where the complexity of the geometry is broken into subdomains of canonical shapes in which the angular regions are immersed in/made of arbitrarily linear media.
This work presents significant advancements in the spectral analysis of electromagnetic problems from different mathematical, physical, and engineering perspectives: (1) the first spectral method capable of handling scattering in arbitrary linear media with multiple propagation constants; (2) direct Fredholm factorization, a novel solution procedure for GWHEs, particularly those with multiple propagation constants; (3) a network interpretation of spectral functional equations and related integral representations for angular regions filled by arbitrary linear media (4) computation of the field at each point within the angular region, avoiding spectral analytical extension; and (5) improved quality of approximate spectral solutions from FIEs by reimposing the GWHEs of the problem, which we name iteration.
This theoretical package has been validated and is ready for future applications.

Author Contributions

V.G.D. and G.L. collaborated in the conceptualization, methodology, mathematics, formal analysis, investigation, validation, writing—original draft preparation, writing—review and editing, project administration, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Universities and Research under PRIN Grant 2017NT5W7Z GREEN TAGS and Next-Generation EU-PNRR M4C2-Inv 1.4–National Centre for HPC, Big Data, and Quantum Computing (HPC)—Multiscale Modeling and Engineering App.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in this study are included in the article and the Appendix A; further inquiries can be directed to the corresponding authors. The data presented in this study were obtained by means of an in-house software code implementing the proposed method.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Explicit Form of Expressions Used in the Main Text

Appendix A contains details and supplements to the main text, in particular explicit expressions that would be cumbersome to report in the main text due to issues with readability and preserving completeness.
The explicit expressions of the 4 × 4 matrices M y o , M y 1 , M y 2 used in Section 2 for (4), i.e.,
M y ( j α o , x ) = M y o + ( x ) M y 1 + ( x ) 2 M y 2
for arbitrary linear media are presented in factorized form as follows:
M y o = j ( ω M y o ( 0 ) + α o M y o ( 1 ) + α o 2 M y o ( 2 ) / ω ) ϵ y y μ y y ζ y y ξ y y
M y o ( 0 ) = ζ x y ζ y z ξ y y ζ x z ζ y y ξ y y + ζ x z μ y y ϵ y y ζ y z μ x y ϵ y y ζ x y μ y y ϵ y z + ζ y y μ x y ϵ y z ζ x x ζ y y ξ y y + ζ x y ζ y x ξ y y ζ x y μ y y ϵ y x + ζ y y μ x y ϵ y x + ζ x x μ y y ϵ y y ζ y x μ x y ϵ y y ζ y y ζ z z ξ y y ζ y z ζ z y ξ y y + ζ y z μ z y ϵ y y ζ z z μ y y ϵ y y ζ y y μ z y ϵ y z + ζ z y μ y y ϵ y z ζ y x ζ z y ξ y y + ζ y y ζ z x ξ y y ζ y y μ z y ϵ y x + ζ z y μ y y ϵ y x + ζ y x μ z y ϵ y y ζ z x μ y y ϵ y y ζ y z ξ y y ϵ x y + μ y y ϵ x y ϵ y z + ζ y y ξ y y ϵ x z μ y y ϵ x z ϵ y y + ζ y z ξ x y ϵ y y ζ y y ξ x y ϵ y z ζ y y ξ y y ϵ x x μ y y ϵ x x ϵ y y ζ y x ξ y y ϵ x y + μ y y ϵ x y ϵ y x ζ y y ξ x y ϵ y x + ζ y x ξ x y ϵ y y ζ y z ξ z y ϵ y y + μ y y ϵ y y ϵ z z + ζ y y ξ z y ϵ y z μ y y ϵ y z ϵ z y + ζ y z ξ y y ϵ z y ζ y y ξ y y ϵ z z ζ y y ξ z y ϵ y x μ y y ϵ y x ϵ z y ζ y x ξ z y ϵ y y + μ y y ϵ y y ϵ z x ζ y y ξ y y ϵ z x + ζ y x ξ y y ϵ z y
ζ x y μ y y ξ y z + ζ x y μ y z ξ y y + ζ y y μ x y ξ y z ζ y y μ x z ξ y y μ x y μ y z ϵ y y + μ x z μ y y ϵ y y ζ x y μ y x ξ y y ζ x y μ y y ξ y x ζ y y μ x x ξ y y + ζ y y μ x y ξ y x + μ x x μ y y ϵ y y μ x y μ y x ϵ y y ζ y y μ z y ξ y z + ζ y y μ z z ξ y y + ζ z y μ y y ξ y z ζ z y μ y z ξ y y μ y y μ z z ϵ y y + μ y z μ z y ϵ y y ζ y y μ z x ξ y y ζ y y μ z y ξ y x ζ z y μ y x ξ y y + ζ z y μ y y ξ y x + μ y x μ z y ϵ y y μ y y μ z x ϵ y y ζ y y ξ x y ξ y z + ζ y y ξ x z ξ y y + μ y y ξ y z ϵ x y μ y z ξ y y ϵ x y μ y y ξ x z ϵ y y + μ y z ξ x y ϵ y y ζ y y ξ x x ξ y y ζ y y ξ x y ξ y x μ y x ξ y y ϵ x y + μ y y ξ y x ϵ x y + μ y x ξ x y ϵ y y μ y y ξ x x ϵ y y ζ y y ξ y y ξ z z + ζ y y ξ y z ξ z y + μ y y ξ z z ϵ y y μ y z ξ z y ϵ y y μ y y ξ y z ϵ z y + μ y z ξ y y ϵ z y ζ y y ξ y x ξ z y ζ y y ξ y y ξ z x μ y x ξ z y ϵ y y + μ y y ξ z x ϵ y y + μ y x ξ y y ϵ z y μ y y ξ y x ϵ z y
M y o ( 1 ) = ζ y z ξ y y μ y y ϵ y z ξ y y ( ζ y x ζ x y ) μ y y ϵ y x + μ x y ϵ y y μ y z ξ y y μ y y ξ y z μ y y ( ζ x y + ξ y x ) + ζ y y μ x y + μ y x ξ y y 0 ζ z y ξ y y μ z y ϵ y y 0 ζ z y μ y y ζ y y μ z y ζ y y ϵ y z ζ y z ϵ y y ξ y y ϵ x y + ζ y y ϵ y x ϵ y y ( ζ y x + ξ x y ) ζ y y ξ y z μ y z ϵ y y ζ y y ( ξ y x ξ x y ) + μ y y ϵ x y + μ y x ( ϵ y y ) 0 ξ z y ϵ y y ξ y y ϵ z y 0 ζ y y ξ z y μ y y ϵ z y
M y o ( 2 ) = 0 ξ y y 0 μ y y 0 0 0 0 0 ϵ y y 0 ζ y y 0 0 0 0
M y 1 = M y 1 ( 0 ) + α o M y 1 ( 1 ) / ω ϵ y y μ y y ζ y y ξ y y
M y 1 ( 0 ) = μ x y ϵ y y ζ x y ξ y y 0 ζ y y μ x y ζ x y μ y y 0 ζ y z ξ y y + ζ z y ξ y y μ z y ϵ y y + μ y y ϵ y z μ y y ϵ y x ζ y x ξ y y ζ y y μ z y + ζ z y μ y y + μ y y ξ y z μ y z ξ y y μ y y ξ y x μ y x ξ y y ξ y y ϵ x y ξ x y ϵ y y 0 μ y y ϵ x y ζ y y ξ x y 0 ζ y z ϵ y y + ξ z y ϵ y y + ζ y y ( ϵ y z ) ξ y y ϵ z y ζ y x ϵ y y ζ y y ϵ y x ζ y y ξ y z + ζ y y ξ z y + μ y z ϵ y y + μ y y ( ϵ z y ) μ y x ϵ y y ζ y y ξ y x
M y 1 ( 1 ) = ξ y y 0 μ y y 0 0 ξ y y 0 μ y y ϵ y y 0 ζ y y 0 0 ϵ y y 0 ζ y y
M y 2 = 0 0 0 0 ξ y y 0 μ y y 0 0 0 0 0 ϵ y y 0 ζ y y 0
We note that the expressions of M y o , M y 1 , M y 2 reported in (A2), (A6), and (A9), respectively, when in factorized form with respect to α o , allow for immediate evaluation in the case of normal incidence α o = 0 by nullifying the contribution of M y o ( 1 ) , M y o ( 2 ) , M y 1 ( 1 ) in (A2) and (A6).

References

  1. Malyuzhinets, G. Excitation, reflection and emission of surface waves from a wedge with given face impedances. Sov. Phys.-Dokl. 1958, 3, 752–755. [Google Scholar]
  2. Budaev, B.V. Diffraction by Wedges; Number 322 in Pitman Research Notes in Mathematics Series; Longman Scientific & Technical: Harlow, UK, 1995. [Google Scholar]
  3. Senior, T.B.A.; Volakis, J.L. Approximate Boundary Conditions in Electromagnetics; Number 41 in IEE Electromagnetic Waves Series; Institution of Electrical Engineers: London, UK, 1995. [Google Scholar]
  4. Osipov, A.V.; Norris, A.N. The Malyuzhinets theory for scattering from wedge boundaries: A review. Wave Motion 1999, 29, 313–340. [Google Scholar] [CrossRef]
  5. Babich, V.M.; Lyalinov, M.A.; Grikurov, V.E. Diffraction Theory: The Sommerfeld-Malyuzhinets Technique; Alpha Science International: Oxford, UK, 2008. [Google Scholar]
  6. Lyalinov, M.A.; Zhu, N.Y. Scattering of Wedges and Cones with Impedance Boundary Conditions; Mario Boella Series on Electromagnetism in Information and Communication Series; SciTech Publishing: Raleigh, NC, USA, 2013. [Google Scholar]
  7. Oberhettinger, F. Diffraction of waves by a wedge. Commun. Pure Appl. Math. 1954, 7, 551–563. [Google Scholar] [CrossRef]
  8. Osipov, A.V. On the method of Kontorovich–Lebedev’s integrals for the problems of diffraction in sectorial media. In Problems of Diffraction and Propagation of Waves; St Petersburg University Publications: Saint Petersburg, Russia, 1993; Volume 25, pp. 173–219. [Google Scholar]
  9. Rawlins, A.D. Diffraction by, or Diffusion into, a Penetrable Wedge. Proc. Math. Phys. Eng. Sci. 1999, 455, 2655–2686. [Google Scholar] [CrossRef]
  10. Salem, M.A.; Kamel, A.H.; Osipov, A.V. Electromagnetic fields in the presence of an infinite dielectric wedge. Proc. R. Soc. A Math. Phys. Eng. Sci. 2006, 462, 2503–2522. [Google Scholar] [CrossRef]
  11. Daniele, V.G.; Lombardi, G. Scattering and Diffraction by Wedges 1: The Wiener-Hopf Solution—Theory; John Wiley & Sons: Hoboken, NJ, USA, 2020. [Google Scholar]
  12. Daniele, V.G.; Lombardi, G. Scattering and Diffraction by Wedges 2: The Wiener-Hopf Solution Advanced Applications; John Wiley & Sons: Hoboken, NJ, USA, 2020. [Google Scholar] [CrossRef]
  13. Daniele, V.G.; Lombardi, G.; Zich, R.S. Network representations of angular regions for electromagnetic scattering. PLoS ONE 2017, 12, e0182763. [Google Scholar] [CrossRef] [PubMed]
  14. Sommerfeld, A. Mathematische Theorie der Diffraction. Math. Ann. 1896, 47, 317–374. [Google Scholar] [CrossRef]
  15. Seshadri, S.; Rajagopal, A. Diffraction by a perfectly conducting semi-infinite screen in an anisotropic plasma. IEEE Trans. Antennas Propag. 1963, 11, 497–502. [Google Scholar] [CrossRef]
  16. Jull, E.V. Diffraction by a conducting half-plane in an anisotropic plasma. Can. J. Phys. 1964, 42, 1455–1468. [Google Scholar] [CrossRef]
  17. Felsen, L.B. Propagation and diffraction in uniaxially anisotropic regions. Part 1. Theory. Proc. Inst. Electr. Eng. 1964, 111, 445–453. [Google Scholar] [CrossRef]
  18. Felsen, L.B. Propagation and diffraction in uniaxially anisotropic regions. Part 2. Applications. Proc. Inst. Electr. Eng. 1964, 111, 454–464. [Google Scholar] [CrossRef]
  19. Williams, W.E. Electromagnetic Diffraction in an Anisotropic Medium. IMA J. Appl. Math. 1966, 2, 186–196. [Google Scholar] [CrossRef]
  20. Hurd, R.A.; Przeździecki, S. Diffraction by a half-plane perpendicular to the distinguished axis of a gyrotropic medium. J. Math. Phys. 1976, 17, 1838–1847. [Google Scholar] [CrossRef]
  21. Przedziecki, S.; Hurd, R.A. Diffraction by a half-plane perpendicular to the distinguished axis of a general gyrotropic medium. Can. J. Phys. 1977, 55, 305–324. [Google Scholar] [CrossRef]
  22. Hurd, R.; Przedziecki, S. Half-plane diffraction in a gyrotropic medium. IEEE Trans. Antennas Propag. 1985, 33, 813–822. [Google Scholar] [CrossRef]
  23. Przedziecki, S. Half-plane diffraction in a chiral medium. Wave Motion 2000, 32, 157–200. [Google Scholar] [CrossRef]
  24. Daniele, V.; Graglia, R.D. Diffraction by an imperfect half plane in a bianisotropic medium. Radio Sci. 2007, 42, 1–16. [Google Scholar] [CrossRef]
  25. Tellegen, B.D. The gyrator, a new electric network element. Philips Res. Rep 1948, 3, 81–101. [Google Scholar]
  26. Cheng, D.; Kong, J.A. Covariant descriptions of bianisotropic media. Proc. IEEE 1968, 56, 248–251. [Google Scholar] [CrossRef]
  27. Lindell, I.; Sihvola, A.; Tretyakov, S.; Viitanen, A.J. Electromagnetic Waves in Chiral and Bi-Isotropic Media; Antennas and Propagation Library, Artech House: Boston, MA, USA, 1994. [Google Scholar]
  28. Olyslager, F. The behavior of electromagnetic fields at edges in bi-isotropic and bi-anisotropic materials. IEEE Trans. Antennas Propag. 1994, 42, 1392–1397. [Google Scholar] [CrossRef]
  29. Vashtalov, S.; Fisanov, V. Diffraction of a plane wave from a wedge in a chiral medium. Russ. Phys. J. 1993, 36, 982–989. [Google Scholar] [CrossRef]
  30. Daniele, V.G.; Lombardi, G.; Zich, R.S. The Electromagnetic Field for a PEC Wedge Over a Grounded Dielectric Slab: 2. Diffraction, Modal Field, Surface Waves, and Leaky Waves. Radio Sci. 2017, 52, 1492–1509. [Google Scholar] [CrossRef]
  31. Daniele, V.; Lombardi, G.; Zich, R.S. The Double PEC Wedge Problem: Diffraction and Total Far Field. IEEE Trans. Antennas Propag. 2018, 66, 6482–6499. [Google Scholar] [CrossRef]
  32. Daniele, V.; Lombardi, G.; Zich, R.S. Radiation and Scattering of an Arbitrarily Flanged Dielectric-Loaded Waveguide. IEEE Trans. Antennas Propag. 2019, 67, 7569–7584. [Google Scholar] [CrossRef]
  33. Daniele, V.G.; Zich, R.S. The Wiener-Hopf Method in Electromagnetics; Mario Boella Series on Electromagnetism in Information and Communication Series; SciTech Publishing: Raleigh, NC, USA, 2014. [Google Scholar]
  34. Daniele, V.G.; Lombardi, G. The generalized Wiener–Hopf equations for wave motion in angular regions: Electromagnetic application. Proc. R. Soc. A Math. Phys. Eng. Sci. 2021, 477, 20210040. [Google Scholar] [CrossRef]
  35. Bresler, A.D.; Marcuvitz, N. Operator Methods in Electromagnetic Field Theory; Report R-495,56, PIB425; MRI Polytechnic Institute of Brooklyn: New York, NY, USA, 1956. [Google Scholar]
  36. Felsen, L.B.; Marcuvitz, N. Radiation and Scattering of Waves; Prentice-Hall: Englewood Cliffs, NJ, USA, 1973. [Google Scholar]
  37. Daniele, V.G.; Zich, R.S. Radiation by arbitrary sources in anisotropic stratified media. Radio Sci. 1973, 8, 63–70. [Google Scholar] [CrossRef]
  38. Noble, B. Methods Based on the Wiener-Hopf Technique for the Solution of Partial Differential Equations; Pergamon Press: London, UK, 1958. [Google Scholar]
  39. DeRusso, P.; Roy, R.; Charles, M. State Variables for Engineers; Wiley: New York, NY, USA, 1965. [Google Scholar]
  40. Lombardi, G. Skew Incidence on Concave Wedge With Anisotropic Surface Impedance. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 1141–1145. [Google Scholar] [CrossRef]
  41. Kantorovich, L.V. Approximate Methods of Higher Analysis; Noordhoff: Groningen, The Netherlands, 1958. [Google Scholar]
  42. Lombardi, G.; Papapicco, D. Quadrature of functions with endpoint singular and generalised polynomial behaviour in computational physics. Comput. Phys. Commun. 2024, 299, 109124. [Google Scholar] [CrossRef]
Figure 1. Angular regions and oblique Cartesian coordinates. (a) The figure reports the Cartesian coordinates z , x , y and the oblique Cartesian coordinate system z , u x , v with reference to angular region 1 of aperture γ ( 0 < φ < γ ) with 0 < γ < π , delimited by faces a and o. In the figure, a second region is identified ( π + γ < φ < 0 ) delimited by faces b and o. The figure also reports the local-to-face Cartesian coordinate systems Z 1 z , X 1 , Y 1 and Z 2 z , X 2 , Y 2 for face a of region 1 and face b of region 2, respectively. The local-to-face Cartesian coordinate systems are obtained from the z , x , y Cartesian coordinate system through rotation for a positive γ and a negative π γ , respectively. (b) The figure shows the new framework of the space divided into two angular regions, which is useful for the study of wedge structures. The figure reports both the Cartesian coordinates z , x , y and the oblique Cartesian coordinate system z , u x , v , where γ is the aperture angle of region 2′. The figure also reports the local-to-face-b Cartesian coordinate system of region 2′ Z 2 z , X 2 , Y 2 , which is obtained from the z , x , y Cartesian coordinate system through a rotation of an angle γ . Finally, in both figures we also use cylindrical coordinates ( z , ρ , φ ) .
Figure 1. Angular regions and oblique Cartesian coordinates. (a) The figure reports the Cartesian coordinates z , x , y and the oblique Cartesian coordinate system z , u x , v with reference to angular region 1 of aperture γ ( 0 < φ < γ ) with 0 < γ < π , delimited by faces a and o. In the figure, a second region is identified ( π + γ < φ < 0 ) delimited by faces b and o. The figure also reports the local-to-face Cartesian coordinate systems Z 1 z , X 1 , Y 1 and Z 2 z , X 2 , Y 2 for face a of region 1 and face b of region 2, respectively. The local-to-face Cartesian coordinate systems are obtained from the z , x , y Cartesian coordinate system through rotation for a positive γ and a negative π γ , respectively. (b) The figure shows the new framework of the space divided into two angular regions, which is useful for the study of wedge structures. The figure reports both the Cartesian coordinates z , x , y and the oblique Cartesian coordinate system z , u x , v , where γ is the aperture angle of region 2′. The figure also reports the local-to-face-b Cartesian coordinate system of region 2′ Z 2 z , X 2 , Y 2 , which is obtained from the z , x , y Cartesian coordinate system through a rotation of an angle γ . Finally, in both figures we also use cylindrical coordinates ( z , ρ , φ ) .
Applsci 14 08685 g001
Figure 2. (a) Angular region 1 of aperture γ , which is delimited by faces a and o with the original reference Cartesian coordinate system z , x , y . The figure also reports the local-to-face-a Cartesian coordinate systems Z 1 z , X 1 , Y 1 . (b) Angular region 1 after a clockwise rotation of an angle γ , becoming region 2’. The figure shows the reference systems of region 1 after the rotation. (c) Angular region 2 of aperture γ , which is delimited by faces b and o with the original reference Cartesian coordinate system z , x , y . The figure also reports the local-to-face-b Cartesian coordinate systems Z 2 z , X 2 , Y 2 . (d) Angular region 2 after a clockwise rotation of an angle γ , becoming region 1’. The figure shows the reference systems of region 2 after the rotation.
Figure 2. (a) Angular region 1 of aperture γ , which is delimited by faces a and o with the original reference Cartesian coordinate system z , x , y . The figure also reports the local-to-face-a Cartesian coordinate systems Z 1 z , X 1 , Y 1 . (b) Angular region 1 after a clockwise rotation of an angle γ , becoming region 2’. The figure shows the reference systems of region 1 after the rotation. (c) Angular region 2 of aperture γ , which is delimited by faces b and o with the original reference Cartesian coordinate system z , x , y . The figure also reports the local-to-face-b Cartesian coordinate systems Z 2 z , X 2 , Y 2 . (d) Angular region 2 after a clockwise rotation of an angle γ , becoming region 1’. The figure shows the reference systems of region 2 after the rotation.
Applsci 14 08685 g002
Figure 3. Two angular regions that are symmetric with respect to the x axis of the aperture angle γ and represent wedge problems immersed in arbitrary linear media; they are modeled using the complete sets of Equations in (45) and (47).
Figure 3. Two angular regions that are symmetric with respect to the x axis of the aperture angle γ and represent wedge problems immersed in arbitrary linear media; they are modeled using the complete sets of Equations in (45) and (47).
Applsci 14 08685 g003
Figure 4. Cauchy smile contour integration line γ 1 t and an example of line m ( t ) for t R , k = 1 0.1 j , γ = 0.7 π (if γ < π / 2 , then the behavior of m ( t ) has a similar path but opposite direction; to intuitively understand this property in an isotropic medium, consider the definition of m in the w plane for supplementary angles γ ).
Figure 4. Cauchy smile contour integration line γ 1 t and an example of line m ( t ) for t R , k = 1 0.1 j , γ = 0.7 π (if γ < π / 2 , then the behavior of m ( t ) has a similar path but opposite direction; to intuitively understand this property in an isotropic medium, consider the definition of m in the w plane for supplementary angles γ ).
Applsci 14 08685 g004
Figure 5. Vertical branch cuts Γ 1 , 2 of m ( t ) originating in branch point ± k assuming a lossy medium (for visibility, k = 1 j ), smile contour integration line γ 1 t , and frown contour integration line γ 2 t with corresponding warped contours λ 1 and λ 2 wrapped around vertical branch cuts Γ 1 and Γ 2 . Note that γ 1 t and γ 2 t assume different observation points for indentation in the figure.
Figure 5. Vertical branch cuts Γ 1 , 2 of m ( t ) originating in branch point ± k assuming a lossy medium (for visibility, k = 1 j ), smile contour integration line γ 1 t , and frown contour integration line γ 2 t with corresponding warped contours λ 1 and λ 2 wrapped around vertical branch cuts Γ 1 and Γ 2 . Note that γ 1 t and γ 2 t assume different observation points for indentation in the figure.
Applsci 14 08685 g005
Figure 6. On top, plots of the absolute values of the spectral solutions | H o z ( k cos w ) | and | H a z ( k cos w ) | obtained as exact solutions and with the FIE approximations for different values of A , h . On the bottom, the corresponding relative errors between the exact solutions and the FIE solutions for different values of A , h on a log 10 scale. Degradation of the spectral solution can be observed near w = π , 0 , which corresponds to η = k , k . The branch point η = k is an offending singularity for the plus spectra, while η = k is related to the physical structural properties of the problem.
Figure 6. On top, plots of the absolute values of the spectral solutions | H o z ( k cos w ) | and | H a z ( k cos w ) | obtained as exact solutions and with the FIE approximations for different values of A , h . On the bottom, the corresponding relative errors between the exact solutions and the FIE solutions for different values of A , h on a log 10 scale. Degradation of the spectral solution can be observed near w = π , 0 , which corresponds to η = k , k . The branch point η = k is an offending singularity for the plus spectra, while η = k is related to the physical structural properties of the problem.
Applsci 14 08685 g006
Figure 7. On top, plots of the absolute values of the spectral solutions | sin w H o z ( k cos w ) | and | sin w H a z ( k cos w ) | obtained as exact solutions and with (1) the FIEs with A = 40 ,   h = 0.025 , (2) the FIEs plus application of one iteration of (91) (FIE+iter), and (3) application of (91) to the source terms of the FIEs (GO+iter). On the bottom, the corresponding relative errors between the exact solutions and the approximate solutions. It is possible to observe an improvement of the solutions near w = π , 0 after applying an iteration of (91) to the approximate solution from the FIEs, yielding a homogenization of the error.
Figure 7. On top, plots of the absolute values of the spectral solutions | sin w H o z ( k cos w ) | and | sin w H a z ( k cos w ) | obtained as exact solutions and with (1) the FIEs with A = 40 ,   h = 0.025 , (2) the FIEs plus application of one iteration of (91) (FIE+iter), and (3) application of (91) to the source terms of the FIEs (GO+iter). On the bottom, the corresponding relative errors between the exact solutions and the approximate solutions. It is possible to observe an improvement of the solutions near w = π , 0 after applying an iteration of (91) to the approximate solution from the FIEs, yielding a homogenization of the error.
Applsci 14 08685 g007
Figure 8. GTD diffraction coefficient (absolute value and phase) for the test problem under consideration: γ = 0.7 π and plane wave illumination with H z polarization with H o = 1 A / m , φ o = 0.1 π , and k = 1 j 0.1 . The figure reports the exact GTD coefficient together with those obtained with the FIE-approximated estimation of the spectra without and with the application of an iteration while selecting A = 20 , h = 0.05 , θ = 0.1 .
Figure 8. GTD diffraction coefficient (absolute value and phase) for the test problem under consideration: γ = 0.7 π and plane wave illumination with H z polarization with H o = 1 A / m , φ o = 0.1 π , and k = 1 j 0.1 . The figure reports the exact GTD coefficient together with those obtained with the FIE-approximated estimation of the spectra without and with the application of an iteration while selecting A = 20 , h = 0.05 , θ = 0.1 .
Applsci 14 08685 g008
Figure 9. Relative error of the GTD diffraction coefficient on a log 10 scale corresponding to the results in Figure 8.
Figure 9. Relative error of the GTD diffraction coefficient on a log 10 scale corresponding to the results in Figure 8.
Applsci 14 08685 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Daniele, V.G.; Lombardi, G. Spectral Analysis of Electromagnetic Diffraction Phenomena in Angular Regions Filled by Arbitrary Linear Media. Appl. Sci. 2024, 14, 8685. https://doi.org/10.3390/app14198685

AMA Style

Daniele VG, Lombardi G. Spectral Analysis of Electromagnetic Diffraction Phenomena in Angular Regions Filled by Arbitrary Linear Media. Applied Sciences. 2024; 14(19):8685. https://doi.org/10.3390/app14198685

Chicago/Turabian Style

Daniele, Vito G., and Guido Lombardi. 2024. "Spectral Analysis of Electromagnetic Diffraction Phenomena in Angular Regions Filled by Arbitrary Linear Media" Applied Sciences 14, no. 19: 8685. https://doi.org/10.3390/app14198685

APA Style

Daniele, V. G., & Lombardi, G. (2024). Spectral Analysis of Electromagnetic Diffraction Phenomena in Angular Regions Filled by Arbitrary Linear Media. Applied Sciences, 14(19), 8685. https://doi.org/10.3390/app14198685

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop