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]):
where the electric and magnetic fields
are related to the electric and magnetic fluxes
and the tensors
are, respectively, the electric permittivity, the magnetic permeability, and the two magneto-electric coupling parameters.
By assuming
- (a)
Cartesian coordinates ,
- (b)
time–harmonic field dependence,
- (c)
Invariant geometry along z and stratification along y, and
- (d)
Sources constituted of plane waves with z-dependence , where depends on the skewness angle with respect to z (with at normal incidence on z),
We obtain the transverse differential equations in matrix form for layered planar media:
where
is the four-dimensional column vector (Throughout the article, we assume the notation
for vectors and not for the modulus of a vector).
Based on the nature of Maxwell’s equations,
is a second-order four-dimensional matrix differential operator of the form
where the explicit forms of the matrices
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
where
and
We introduce here an analysis of the operator
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
is semi-simple, we compute its eigenvalues
and eigenvectors
as follows:
i.e.,
where
and
(dependence on
and
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].
This yields the following four eigenvalues:
where
with
We note that the column vectors
of
provide a basis in the space
, where we define the transverse electromagnetic field
, while the column vectors
of
in the reciprocal space will be fundamental in obtaining functional equations through the characteristic Green function procedure. Each couple
is related to a single
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 (
,
,
),
assumes the following forms:
where in the presence of losses (
) we have
and
, which are respectively related to progressive (
) and regressive (
) waves with regard to
y of form
. 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,
) with a positive real part
representing progressive waves and two (say,
) with a negative real part
representing regressive waves, yielding all four
y-longitudinal propagation constants with
(
y-progressive/regressive waves
, i.e., assuming time-harmonic dependence
, we have
x,
y,
z progressive waves
, respectively, with
).
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):
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
:
The application of a Fourier transform along
reduces (
19) to
where
and
Based on the link between
and
, we have
with the same eigenvectors
of
, and the following relationship between the eigenvalues
and
can be obtained:
resulting in the following “oblique”
v-longitudinal propagation constants:
which is in agreement with the relationship between
and
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
). We note that the quantities
,
, and
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
.
The benefit of using the Laplace transform is correlated with the incorporation of boundary conditions through initial conditions with the term
. In (
26)–(
27), the condition
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.
Now, considering the asymptotic behavior in (
28) of exponential functions in
v for
, we need to have
; at the same time, the first couple of integrals are null, as
and
are related to progressive and regressive waves, respectively. For this reason, setting
, we obtain the following spectral field representation along face
o:
By weighting (
29) with the reciprocal vectors
of
, we obtain the following functional equations:
where we use the following definition of the Laplace transform:
With reference to
Figure 1a and its caption, analyzing and expanding
in (
30) using Maxwell’s equations, we rephrase the functional equations for region 1 as
where
is the
Laplace transform of the tangent-to-face-
o field components (i.e., at
) in Cartesian
and cylindrical coordinates
(omitting the
z coordinate for invariance)
where
is the
Laplace transform of the tangent-to-face-
a field components (i.e., at
) in local-to-face-
a Cartesian
coordinates and global cylindrical coordinates
(located at
)
and
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 (
) and face
a (
). Furthermore, we observe that the angle
is essential in determining the impact of anisotropies through
.
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
:
where
has the same expression of
(
27) but with a different support
, and allows the incorporation of boundary conditions along face
b (
). 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
and the source term
, 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
, we need to have
; at the same time, the second couple of integrals are null. For this reason, setting
, we obtain
By weighting (
37) with the reciprocal vectors
of
, we obtain the functional equations
where we have used the definition of the
v left Laplace transform
Note the differences and similarities between Laplace transformations (
31) and (
39), which yield the same definition of the
Laplace transform in
but are applied to different quantities. Furthermore, the regularity properties of the
Laplace transform are inherited from
(
) according to (
24)–(
25).
With reference to
Figure 1a and its caption, analyzing and expanding
in (
38), we rephrase the functional equations into
In (
40),
is again (
35); as for region 1,
is needed for
in region 2 with respect to
in region 1,
is the
Laplace transform of the tangent-to-face-
o field components reported in (
33), and
is the
Laplace transform of the tangent-to-face-
b field components (i.e., at
) in the local-to-face-
b Cartesian
coordinates and global cylindrical coordinates
of
Figure 1a.
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
where
which is the
Laplace transform of the tangent-to-face-
b field components (i.e., now at
) in local-to-face-
b Cartesian
coordinates and global cylindrical coordinates
in
Figure 1b. Note that in (
42) we assumed that region 2′ is homogeneous to region 1, yielding the same
; 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 (
) and face
b (
). In (
42), note the new dependence of
(
35) on
due to the effect of anisotropies while changing the orientation of face
b from
to
. Furthermore, in the case of symmetric media (
), we have
; 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),
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
components (
33) related to two combinations of
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
(
33) and
(
43).
We further note that in (
44) the components of face
o and face
are functions of the spectral variables
and
, respectively, which are related to one another via (
24)–(
25). We can reverse the role of the variables
and
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
). 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:
In the second pair of Equation (
45), we have used the subscript
to refer to a rotated coordinated system
with respect to
; 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
,
,
(from
) because of their dependence on the constitutive tensorial parameters
of region 1, which are redefined in the reference coordinate system
, i.e.,
, for example,
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
In the first pair of equations in (
47), we have used the subscript
to make reference to a rotated coordinated system
with respect to
; 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
,
, and
because of their dependence on the constitutive parameters
, which are redefined in the reference coordinate system
, i.e.,
, for example,
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.,
) we have
; 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
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
plane (i.e., one of the
that all depends on
) applied to an arbitrary
(i.e., the standard form is obtained simply by replacing
with
):
where
is the non-standard contribution of
in the
plane. Assuming a lossy medium, we observe that, in general,
is located with the positive imaginary part for
, 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
and
, yielding a path of
for
from right to left because of
. On the contrary, for
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
, 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:
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
and the second to the field component
). 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
where on the LHS we have plus field unknowns in
and on the RHS we have minus field unknowns in
. The apparent redundancy in (
51) after the imposition of boundary conditions is exploited to obtain integral representations only in terms of the field components
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:
recalling that all occurrences of
are functions of
, i.e.,
. 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
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:
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:
Given the expressions of
in (
24)–(
25) with
in lossy media (
Figure 4), Equation (
54) requires the computation of
which can be performed either numerically or analytically while paying attention to the branch cuts of
. Furthermore, in (
54) we also need to consider n.s. singularities related to the field.
The validity of the estimation of
extends to complex values of
as long as
does not cross the singularity line determined by the behavior of
for
, 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
, 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
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
. It is important to highlight that the system only depends on the spectral variable
, ensuring that functions do not rely on
outside of the integration sign. This property is fundamental to avoid the analysis of unknowns defined in different complex planes (
and various
) 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
From the homogeneous portion of the solution in (
56), we obtain the definitions of arbitrary coefficients in terms of the field components at
(face
o):
The particular integrals in (
56) are terms related to face
a via
. Due to the linearity of the problem, we apply the superposition principle; we can interpret (
56) as the result of an equivalent theorem, where
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
where we note that
impacts all terms of the solution, as already reported in
Section 3, i.e.,
,
,
,
, and the field components. The solution takes the form
where
in
Figure 2b is now different from
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
(face
a):
The particular integrals in (
59) are terms related to face
o via
. 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
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
; 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
and subregion B
). After the face spectra at
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. 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
is
where
z is the distinguished axis of the medium and
,
. As reported in [
15], this vector problem is separable into two equivalent scalar problems for the
(H-mode) and
(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:
where
,
,
,
, and
.
The corresponding eigenvectors
, from which we can easily compute the reciprocal vectors
through inversion, are
The problem shows simplification because of
; see for instance the impact of the anisotropies on (
35) or
However, we keep the procedure as general as possible, i.e., applicable to wedge problems (arbitrary
); from (
27) and (
31), we obtain
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:
Similarly, the procedure can be repeated for region 2. The complete set of equations highlights the decoupling of
from
polarization. Applying the PEC boundary conditions on the faces, after some manipulations we obtain the following equations for the
and
polarizations:
Now, we impose
, i.e., the angular regions are defined for the half-plane problem. In (
144), we note that
polarization behaves as a half-plane problem immersed in a classical isotropic region [
33] but with a propagation constant of
, i.e., this can be interpreted via a network representation with characteristic impedance
, confirming the findings of [
15].
With further mathematical manipulation of (
145), we obtain
and the second equation in (
146) shows the same WH kernel as that of Equation (
25) in [
15]:
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.