Electromagnetic Scattering from Surfaces with Curved Wedges Using the Method of Auxiliary Sources (MAS)

: The method of auxiliary sources (MAS) is utilized in the analysis of Transverse Magnetic (TM) plane wave scattering from inﬁnite, conducting, or dielectric cylinders, including curved wedges. The latter are deﬁned as intersections of circular arcs. The artiﬁcial surface, including the auxiliary sources, is shaped in various patterns to study the e ﬀ ect of its form on the MAS accuracy. In juxtaposition with the standard, conformal shape, several deformations are tested, where the auxiliary sources are forced to approach the tip of the wedge. It is shown that such a procedure signiﬁcantly improves the accuracy of the numerical results. Comparisons of schemes are presented, and the optimal auxiliary source location is proposed.


Introduction
Computational Electromagnetics techniques have traditionally been utilized in the mathematical modeling of problems related to radiation and the propagation of electromagnetic waves. Standard applications involve antennas, diffraction, scattering, waveguides, propagation in complex environments, etc., while boundary surfaces are typically either perfectly conducting or dielectric. Recently, important technological advances have been based on the electromagnetic properties of Near-Zero Index (NZI) materials [1,2], plasmonics [3], metasurfaces [4], graphene [5], nanoparticles [6], etc.
Analytical and asymptotic methods were originally the only mathematical tool to perform relevant calculations, naturally being restricted to canonical geometries. However, computers demonstrated explosive progress over the last few decades, thus rendering numerical techniques indispensable in extracting approximate results of controllable accuracy for arbitrary geometries. In general terms, numerical techniques discretize either integral [7,8] or differential equations [9,10] involved in the mathematical simulation of physical problems. A variant of differential equation techniques is based on the construction of equivalent circuit or transmission line models [11][12][13]. Each group of the aforementioned techniques is characterized by advantages and disadvantages, and therefore their suitability and efficiency depend on the particular configuration. For example, in scattering problems, where the Radar Cross Section (RCS) of a target must be calculated, differential equation methods should ideally discretize the entire space all the way to infinity, which is obviously impossible. Mesh truncation is a complicated task that needs to be resolved in order to obtain accurate results. On the contrary, integral equation methods are usually more complicated in terms of calculations, and the

MAS for Eye-Shaped Scatterers
We assume an infinitely long cylinder with a cross section that resembles an eye ( Figure 1). We investigate two separate cases: perfect electric conductor (PEC) or linear-homogeneous-isotropic dielectric. The geometry of the scatterer, depicted by a solid line, comprises two circular arcs with identical radii equal to ρ, but with different centers. In particular, the Cartesian coordinates of the upper arc are given by (1): whereas those of the lower arc are given by (2): x l = ρ cos ϕ, y l = ρ sin ϕ + d, where ϕ is the azimuth angle and ±d is the vertical displacement of each arc center, taken equal to the arc apothem (see Figure 2). Obviously, ϕ does not span the entire [0, 2π) interval, but it is limited by the arc width itself, which is given by (3): Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 16

Improvement of the auxiliary surface layout
As mentioned in [19,20], MAS becomes less accurate when the auxiliary surfaces are conformal to boundaries encompassing wedges. In particular, satisfaction of the boundary condition at midpoints (MPs) (see Figure 1) close to the tips is no longer sufficiently good. To overcome this complication, the auxiliary surface, defined by radius and using expressions analogous to (1-3), may be deformed so that ASsnot only approach the tips closely, but become denser in the tip neighborhood as well. ASs may approach CPs following several patterns. In this work, two basic patterns were tested. Let be the number of ASs to be moved. Let be the polar radius of the th AS ( = 1,2, … , ), and let be the maximum polar radius difference between the th AS and

Improvement of the auxiliary surface layout
As mentioned in [19,20], MAS becomes less accurate when the auxiliary surfaces are conformal to boundaries encompassing wedges. In particular, satisfaction of the boundary condition at midpoints (MPs) (see Figure 1) close to the tips is no longer sufficiently good. To overcome this complication, the auxiliary surface, defined by radius and using expressions analogous to (1-3), may be deformed so that ASsnot only approach the tips closely, but become denser in the tip neighborhood as well. ASs may approach CPs following several patterns. In this work, two basic patterns were tested. Let be the number of ASs to be moved. Let be the polar radius of the th AS ( = 1,2, … , ), and let be the maximum polar radius difference between the th AS and The scatterer is illuminated by a TM plane wave impinging from an azimuth angle equal to ϕ inc . Therefore, the incident electric field E inc is given by (4): where E 0 is the amplitude of the incident electric field, k 0 is the free space wavenumber, and a hat denotes a unit vector along the corresponding direction. The incident magnetic field H inc is given by (5): where ζ 0 is the free space intrinsic impedance. To solve the scattering problem via MAS, two sets of ASs are generally defined, each one of multitude N, as shown in Figure 1. In the PEC case, only the inner set is used; the outer set is necessary only in the dielectric configuration. In standard MAS formulation, both inner and outer auxiliary surfaces are conformal to the scatterer boundary. The electric field due to the nth inner AS, located at point r n and radiating in the outer space, is as follows: where E n is the corresponding unknown weight, (n = 1, 2, . . . , N), and H (2) 0 is the Hankel function of the zero order and second kind (2D Green's function). The corresponding magnetic field of the nth auxiliary source is obviously proportional to the curl of (6), given explicitly in [20]. Similar expressions hold for the outer ASs, radiating in the inner space of the dielectric scatterer, except for k 0 and ζ 0 , which have to be replaced by k and ζ, respectively, corresponding to the scatterer's dielectric properties. The total scattered E field is expressed as the superposition of the fields in (6) and the H field accordingly. By applying the boundary conditions for both fields at the N collocation points (CPs) (x m , y m ) (m = 1, 2, . . . , N) of the scattering boundary (blue solid dots in Figure 1), we cast a linear system of equations as in (7): where {I} is the column vector of the unknown weights E n . In the PEC case, [Z] is a square matrix of size N × N with elements determined by the interaction between ASs and CPs, and {V} is the column vector of the incident E fields calculated at the CPs. In the dielectric case, the [Z] interaction matrix is of size 2N × 2N and {V} is the column vector of both the incident E and H fields calculated at the CPs.

Improvement of the Auxiliary Surface Layout
As mentioned in [19,20], MAS becomes less accurate when the auxiliary surfaces are conformal to boundaries encompassing wedges. In particular, satisfaction of the boundary condition at midpoints (MPs) (see Figure 1) close to the tips is no longer sufficiently good. To overcome this complication, the auxiliary surface, defined by radius ρ aux and using expressions analogous to (1)(2)(3), may be deformed so that ASsnot only approach the tips closely, but become denser in the tip neighborhood as well. ASs may approach CPs following several patterns. In this work, two basic patterns were tested. Let M be the number of ASs to be moved. Let ρ m be the polar radius of the mth AS (m = 1, 2, . . . , M), and let g be the maximum polar radius difference between the mth AS and the mth CP. Finally, let s be the proximity factor, defined in [0, 1], so that 0 corresponds to no approach and 1 corresponds to the maximum approach (resulting in coincident ASs and CPs). Then, the schemes proposed for the auxiliary surface deformation are defined in (8) as follows: where ν = 1 for simple and ν = 2 for progressive reach. The deformation effect is graphically described in Figure 3. Furthermore, in a way similar to [20], ASs and CPs, accordingly, should become denser close to the wedge tip. Again, there is no unique way to accomplish this. In this work, the scheme implemented multiplies the polar angle ϕ m of the mth AS location by a factor D m , 0 < D start ≤ D m ≤ 1, where D start is user-defined. For example, in the first quadrant of the 'eye', D m is defined as being close to 0 for ASs near the wedge tip, and close to 1 for ASs close to the vertical axis. For progressive densification, the scheme proposed is: ϕ m = ϕ m D 2 m . Moreover, additional ASs may be superimposed to the already existing ones close to the tips if necessary. The combined effect of the tip approach and densification is depicted in Figure 4.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 5 of 16 superimposed to the already existing ones close to the tips if necessary. The combined effect of the tip approach and densification is depicted in Figure 4.

PEC scatterer
To test the efficiency of the method, a PEC scatterer is initially defined by radius = 3 , arc displacement = 0.5 = 1.5 , incidence angle = 0, auxiliary surface radius = − 0.075 , and originally 136 CPs and therefore 136 inner ASs. The solution to the problem without any deformation yields the results of Figure 5(a), which depicts the quantified error in the generic boundary condition (BC) of the field along the boundary stretch, i.e.: where in (9) is the normal unit vector on the boundary, pointing outwards, and , are the electric fields just inside and just outside the scatterer, respectively, the former obviously vanishing for the PEC case.
To improve the satisfaction of the BC, the deformation scheme proposed above was implemented. After several trials, the following parameters were finally invoked: All ASs were displaced as well as densified and the proximity factor was set equal to = 0.88, whereas = 0.9. No extra ASs were added since their presence proved to achieve only incremental improvement.
The results are displayed in Figure 5(b), which shows that the overall field error at MPs was substantially reduced. superimposed to the already existing ones close to the tips if necessary. The combined effect of the tip approach and densification is depicted in Figure 4.

PEC scatterer
To test the efficiency of the method, a PEC scatterer is initially defined by radius = 3 , arc displacement = 0.5 = 1.5 , incidence angle = 0, auxiliary surface radius = − 0.075 , and originally 136 CPs and therefore 136 inner ASs. The solution to the problem without any deformation yields the results of Figure 5(a), which depicts the quantified error in the generic boundary condition (BC) of the field along the boundary stretch, i.e.: where in (9) is the normal unit vector on the boundary, pointing outwards, and , are the electric fields just inside and just outside the scatterer, respectively, the former obviously vanishing for the PEC case.
To improve the satisfaction of the BC, the deformation scheme proposed above was implemented. After several trials, the following parameters were finally invoked: All ASs were displaced as well as densified and the proximity factor was set equal to = 0.88, whereas = 0.9. No extra ASs were added since their presence proved to achieve only incremental improvement.
The results are displayed in Figure 5(b), which shows that the overall field error at MPs was substantially reduced.

PEC Scatterer
To test the efficiency of the method, a PEC scatterer is initially defined by radius ρ = 3λ, arc displacement d = 0.5ρ = 1.5λ, incidence angle ϕ inc = 0, auxiliary surface radius ρ aux = ρ − 0.075λ, and originally 136 CPs and therefore 136 inner ASs. The solution to the problem without any deformation yields the results of Figure 5a, which depicts the quantified error in the generic boundary condition (BC) of the E field along the boundary stretch, i.e.,: where in (9)n is the normal unit vector on the boundary, pointing outwards, and E in , E out are the electric fields just inside and just outside the scatterer, respectively, the former obviously vanishing for the PEC case.
case, keeping the multitude of CPs consistent (and thus the multitude of ASs consistent) but moving the auxiliary surface further away from the scattering boundary ( = − 0.15 ) produced the results in Figure 6(a) and (b), before and after applying the deformation scheme, respectively. Severe complications are revealed in Figure 6(a), while it may be deduced from Figure 6(b) that utilizing an improved auxiliary surface layout, compared to a conformal one, may be much less prone to accuracy degradation, presumably because caustics are still retained within the auxiliary surface.

Lossless dielectric scatterer
Next, a lossless dielectric scatterer was considered, with radius = 3 , arc displacement = 0.5 = 1.5 , relative permittivity = 2.56, incidence angle = 0, inner and outer auxiliary surface radius = − 0.3 and = + 0.45 , and, originally, 160 CPs (hence 160 inner and 160 outer ASs.). Like in the PEC case, the efficiency of MAS was initially tested using a conformal auxiliary surface and the extracted results of the BC error are presented in Figure 7(a), where the upper subplot depicts the quantified error for the field ( %) and the lower one for the field ( %), both being relatively significant near the wedge tips. To improve the satisfaction of the BC, the deformation scheme proposed above was implemented. After several trials, the following parameters were finally invoked: All ASs were displaced as well as densified and the proximity factor was set equal to s = 0.88, whereas D start = 0.9. No extra ASs were added since their presence proved to achieve only incremental improvement. The results are displayed in Figure 5b, which shows that the overall E field error at MPs was substantially reduced.
It is worth noting that the proper placement of the auxiliary surface is of great importance in order to avoid caustics exclusion by the auxiliary surface or any potential resonance effects, as presented in [23]. To illustrate the possibly detrimental influence on the BC error, specifically, in this case, keeping the multitude of CPs consistent (and thus the multitude of ASs consistent) but moving the auxiliary surface further away from the scattering boundary (ρ aux = ρ − 0.15λ) produced the results in Figure 6a,b, before and after applying the deformation scheme, respectively. Severe complications are revealed in Figure 6a, while it may be deduced from Figure 6b that utilizing an improved auxiliary surface layout, compared to a conformal one, may be much less prone to accuracy degradation, presumably because caustics are still retained within the auxiliary surface.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 16 It is worth noting that the proper placement of the auxiliary surface is of great importance in order to avoid caustics exclusion by the auxiliary surface or any potential resonance effects, as presented in [23]. To illustrate the possibly detrimental influence on the BC error, specifically, in this case, keeping the multitude of CPs consistent (and thus the multitude of ASs consistent) but moving the auxiliary surface further away from the scattering boundary ( = − 0.15 ) produced the results in Figure 6(a) and (b), before and after applying the deformation scheme, respectively. Severe complications are revealed in Figure 6(a), while it may be deduced from Figure 6(b) that utilizing an improved auxiliary surface layout, compared to a conformal one, may be much less prone to accuracy degradation, presumably because caustics are still retained within the auxiliary surface.

Lossless dielectric scatterer
Next, a lossless dielectric scatterer was considered, with radius = 3 , arc displacement = 0.5 = 1.5 , relative permittivity = 2.56, incidence angle = 0, inner and outer auxiliary surface radius = − 0.3 and = + 0.45 , and, originally, 160 CPs (hence 160 inner and 160 outer ASs.). Like in the PEC case, the efficiency of MAS was initially tested using a conformal auxiliary surface and the extracted results of the BC error are presented in Figure 7(a), where the upper subplot depicts the quantified error for the field ( %) and the lower one for the field ( %), both being relatively significant near the wedge tips.

Lossless dielectric scatterer
Next, a lossless dielectric scatterer was considered, with radius ρ = 3λ, arc displacement d = 0.5ρ = 1.5λ, relative permittivity ε r = 2.56, incidence angle ϕ inc = 0, inner and outer auxiliary Appl. Sci. 2020, 10, 2309 7 of 16 surface radius ρ aux in = ρ − 0.3λ and ρ aux out = ρ + 0.45λ, and, originally, 160 CPs (hence 160 inner and 160 outer ASs). Like in the PEC case, the efficiency of MAS was initially tested using a conformal auxiliary surface and the extracted results of the BC error are presented in Figure 7a, where the upper subplot depicts the quantified error for the E field (∆ %) and the lower one for the H field (∆H %), both being relatively significant near the wedge tips. Afterwards, the improvement of the BC was examined by employing the aforementioned deformation scheme with the following parameters, obtained through multiple trials: the portion of ASs to be displaced was 1/5 for the inner and 1/8 for the outer ones, with a proximity factor of = 0.75 and = 0.65, respectively, while all of them were densified with = 0.80 . The addition of extra ASs offered no further refinement and was skipped once more. Figure 7(b) depicts the considerably decreased quantified BC error, and the corresponding radar cross section (RCS) is displayed in Figure 7(c).
(a) (b) (c) Figure 7. Results for lossless dielectric scatterer: (a) and field boundary the condition error by employing a conformal auxiliary surface; (b) and field boundary the condition error by employing a deformed auxiliary surface; (c) the Radar Cross Section (RCS) acquired by utilizing the auxiliary surface deformation scheme.

Lossy dielectric scatterer
Finally, a lossy dielectric scatterer of radius = 5 , arc displacement = 0.5 = 2.5 , relative permittivity = 2.56 − 0.102 , incidence angle = 0, inner and outer auxiliary surface radius = − 0.5 and = + 0.75 , and originally 180 CPs (hence 180 inner and 180 outer ASs) was analyzed. Similar BC tests were conducted for conformal auxiliary surface, once again resulting in a relatively sizable BC error in the neighborhood of the wedge tips (Figure 8(a)). The deformation scheme proposed above was also implemented in order to enhance the satisfaction of the boundary condition. The parameters found to produce the lowest BC error, shown in Figure 8(b), are as follows: 1/5 of the inner ASs and 1/8 of the outer ones were displaced, with a proximity factor set equal to = 0.75 and = 0.65, respectively, whereas all of them were included in the densification with = 0.80. Adding extra ASs did not prove to be beneficial and, therefore, was not employed. The maximum field error at MPs was reduced from 5.13 * 10 to 7.27 * 10 and the maximum field error from 5.11 * 10 to 4.14 * 10 . Afterwards, the improvement of the BC was examined by employing the aforementioned deformation scheme with the following parameters, obtained through multiple trials: the portion of ASs to be displaced was 1/5 for the inner and 1/8 for the outer ones, with a proximity factor of s = 0.75 and s = 0.65, respectively, while all of them were densified with D start = 0.80. The addition of extra ASs offered no further refinement and was skipped once more. Figure 7b depicts the considerably decreased quantified BC error, and the corresponding radar cross section (RCS) is displayed in Figure 7c

Lossy Dielectric Scatterer
Finally, a lossy dielectric scatterer of radius ρ = 5λ, arc displacement d = 0.5ρ = 2.5λ, relative permittivity ε r = 2.56 − 0.102 j, incidence angle ϕ inc = 0, inner and outer auxiliary surface radius ρ aux in = ρ − 0.5λ and ρ aux out = ρ + 0.75λ, and originally 180 CPs (hence 180 inner and 180 outer ASs) was analyzed. Similar BC tests were conducted for conformal auxiliary surface, once again resulting in a relatively sizable BC error in the neighborhood of the wedge tips (Figure 8a). observable improvement, e.g., from 300° to 360° in Figure 8(e), in the azimuthal direction where and obtain their greatest values. It is worth mentioning that placing ASs too close to the wedge tip is not advisable. Indeed, an AS right on the wedge tip causes severe ill-conditioning of the linear system matrix, whereas proximity factor s = 0.99 for both inner and outer ASs lead to inadequate error reduction, namely 3.16 * 10 for the field and 4.43 * 10 for the field.  The deformation scheme proposed above was also implemented in order to enhance the satisfaction of the boundary condition. The parameters found to produce the lowest BC error, shown in Figure 8b, are as follows: 1/5 of the inner ASs and 1/8 of the outer ones were displaced, with a proximity factor set equal to s = 0.75 and s = 0.65, respectively, whereas all of them were included in the densification with D start = 0.80. Adding extra ASs did not prove to be beneficial and, therefore, was not employed.
The computed RCSs, using conformal and deformed auxiliary surface layouts, are depicted in Figure 8c,d, respectively. Although not visually discernible in these plots due to their wide dynamic range, comparing their values at angles close to the scatterer's wedges exhibits an observable improvement, e.g., from 300 • to 360 • in Figure 8e, in the azimuthal direction where ∆ and ∆H obtain their greatest values.
It is worth mentioning that placing ASs too close to the wedge tip is not advisable. Indeed, an AS right on the wedge tip causes severe ill-conditioning of the linear system matrix, whereas proximity factor s = 0.99 for both inner and outer ASs lead to inadequate error reduction, namely 3.16 * 10 −2 for the E field and 4.43 * 10 −1 for the H field.
Furthermore, the symmetry of the geometry should be accompanied by symmetry in the ASs/CPs. As a counterexample, the same lossy dielectric scatterer was experimentally analyzed via unequal CP densities in the two arcs, namely 45 for the upper and 90 for the lower one, leading to excessively high errors and completely wrong RCS plots, as demonstrated in Figure 9. Furthermore, the symmetry of the geometry should be accompanied by symmetry in the ASs/CPs. As a counterexample, the same lossy dielectric scatterer was experimentally analyzed via unequal CP densities in the two arcs, namely 45 for the upper and 90 for the lower one, leading to excessively high errors and completely wrong RCS plots, as demonstrated in Figure 9.

Investigation of the solution behavior for various geometry modifications
To check the robustness of the method described above, several modifications of the geometry were carried out and the respective results were extracted. In all cases below, lossy dielectric material with relative permittivity = 2.56 − 0.102 is considered, whereas the incidence angle is always = 0 with an arc radius of = 3 .
First, a scatterer that deviates only slightly from a circle was studied, i.e., a layout with small arc displacement = 0.1 = 0.3 , followed by a scatterer with a moderate deviation from a circle, i.e., a layout with mediocre arc displacement = 0.5 = 1.5 , and, finally, a scatterer with a significant distortion with respect to a circle, which, therefore, encompasses sharp wedges, i.e., a layout with large arc displacement = 0.7 = 2.1 . All three geometries and their computed results are depicted in Figures 10, 11, and 12, respectively, while their MAS parameters and the quantified error in the generic boundary condition of the electromagnetic field are presented in Table 1.

Investigation of the Solution Behavior for Various Geometry Modifications
To check the robustness of the method described above, several modifications of the geometry were carried out and the respective results were extracted. In all cases below, lossy dielectric material with relative permittivity ε r = 2.56 − 0.102 j is considered, whereas the incidence angle is always ϕ inc = 0 with an arc radius of ρ = 3λ.
First, a scatterer that deviates only slightly from a circle was studied, i.e., a layout with small arc displacement d = 0.1ρ = 0.3λ, followed by a scatterer with a moderate deviation from a circle, i.e., a layout with mediocre arc displacement d = 0.5ρ = 1.5λ, and, finally, a scatterer with a significant distortion with respect to a circle, which, therefore, encompasses sharp wedges, i.e., a layout with large arc displacement d = 0.7ρ = 2.1λ. All three geometries and their computed results are depicted in Figures 10-12, respectively, while their MAS parameters and the quantified error in the generic boundary condition of the electromagnetic field are presented in Table 1.
As a verification of the broad applicability of the technique, an asymmetric scatterer was also analyzed. Namely, two unequal circular arcs were connected, forming a "sorrowful" eye. The geometry parameters were chosen as follows: radius = 4 and arc displacement = 0.5 = 2 for the upper arc, radius = 5 and arc displacement = 0.5 = 2.5 for the lower arc, originally 160 CPs (hence 160 inner and 160 outer ASs), 1/2.5 of the inner ASs and 1/4 of the outer ones were displaced, with proximity factor set equal to = 0.6 and = 0.6, respectively. No densification was implemented and no extra ASs were added. The maximum field error at MPs was reduced from 4.19 * 10 to 1.57 * 10 and the maximum field error from 4.36 * 10 to 1.38 * 10 . The results are depicted in Figure 13.
Judging from the outcome of the aforementioned tests, the method proposed herein is very robust, producing reliable results for blunt, acute, or non-symmetric wedges.  (e) (f) As a verification of the broad applicability of the technique, an asymmetric scatterer was also analyzed. Namely, two unequal circular arcs were connected, forming a "sorrowful" eye. The geometry parameters were chosen as follows: radius ρ = 4λ and arc displacement d = 0.5ρ = 2λ for the upper arc, radius ρ = 5λ and arc displacement d = 0.5ρ = 2.5λ for the lower arc, originally 160 CPs (hence 160 inner and 160 outer ASs), 1/2.5 of the inner ASs and 1/4 of the outer ones were displaced, with proximity factor set equal to s = 0.6 and s = 0.6, respectively. No densification was implemented and no extra ASs were added. The maximum E field error at MPs was reduced from 4.19 * 10 −2 to 1.57 * 10 −2 and the maximum H field error from 4.36 * 10 −1 to 1.38 * 10 −1 . The results are depicted in Figure 13.  Judging from the outcome of the aforementioned tests, the method proposed herein is very robust, producing reliable results for blunt, acute, or non-symmetric wedges.

Discussion
The applicability range of the approach presented above may span various aspects of computational electromagnetics. MAS has already been used in the literature to simulate practical problems, such as scattering by a raindrop [24], although that particular geometry is smooth, without wedges. Wedge treatment, as discussed in this paper, enhances MAS accuracy for applications such as jet engine inlet modeling [25][26][27], where interior blades contain sharp edges, or even further military aircraft scattering simulation, such a stealth design. Moreover, additional possible applications may include automotive modeling, for example, the functionality of a telecommunications antenna in the presence of the vehicle surface. Furthermore, the algorithm may be used in the design of absorbers and cloaking materials [28][29][30], or generally wave manipulation [31] when the geometry layout contains wedges.
Although the numerical results presented were very good, they are limited to wedges formed

Discussion
The applicability range of the approach presented above may span various aspects of computational electromagnetics. MAS has already been used in the literature to simulate practical problems, such as scattering by a raindrop [24], although that particular geometry is smooth, without wedges. Wedge treatment, as discussed in this paper, enhances MAS accuracy for applications such as jet engine inlet modeling [25][26][27], where interior blades contain sharp edges, or even further military aircraft scattering simulation, such a stealth design. Moreover, additional possible applications may include automotive modeling, for example, the functionality of a telecommunications antenna in the presence of the vehicle surface. Furthermore, the algorithm may be used in the design of absorbers and cloaking materials [28][29][30], or generally wave manipulation [31] when the geometry layout contains wedges.
Although the numerical results presented were very good, they are limited to wedges formed by intersecting circular arcs only. It is anticipated that analogous techniques would be equally efficient in more general cases (i.e., arbitrarily shaped curved wedges); however, this is yet to be proven in future research. Likewise, other issues to be investigated include TE (Transverse Electric) incidence, the three-dimensional counterpart of this problem, unconventional materials, as well as a robust methodology to choose the optimal deformation parameters.
Furthermore, by studying various geometries, proposing schemes for wedge treatment and obtaining solutions by choosing alternative configurations for the incident wave, material properties, etc., an extensive data set can be acquired. The latter, a combination of input variables and their respective exported results may be used as training data for building a mathematical model via Machine Learning algorithms. Machine Learning, the usage of computational methods in order to produce/"learn" information directly from data without relying on explicit instructions, improves its accuracy and performance proportionally to the amount of available training data. A larger data set is, therefore, required to find hidden patterns or intrinsic properties that can be utilized in the procedure of determining the optimal settings for the scattering problem and, hopefully, automating it.

Conclusions
The method of auxiliary sources (MAS) was applied to scattering from a PEC or dielectric cylinder whose geometry contains curved wedges, formed by intersections of circular arcs. The BC error in the vicinity of the wedge tips is significant for standard, conformal auxiliary surfaces due to field singularities. To overcome such a complication, the deformation of the auxiliary sources was proposed. When close to the tips, both inner and outer (when applicable) ASs approached the CPs, and their distribution was also forced to become denser. The BC error was shown to decrease significantly for both the E and the H fields, resulting in more accurate RCS results.