Next Article in Journal
Drone Detection Using YOLOv5
Previous Article in Journal
Experimental and Numerical Analysis of a Low-Cost Solar Still
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Refinement of the Heuristic Solution for the Problem of TE-Polarized Electromagnetic Wave Diffraction on a Half-Plane with Two-Sided Impedance Boundary Conditions

by
Michael V. Vesnik
Kotelnikov Institute of Radioengineering and Electronics, Russian Academy of Sciences, 125009 Moscow, Russia
Eng 2023, 4(1), 404-415; https://doi.org/10.3390/eng4010024
Submission received: 24 December 2022 / Revised: 20 January 2023 / Accepted: 28 January 2023 / Published: 1 February 2023
(This article belongs to the Section Electrical and Electronic Engineering)

Abstract

:
A technique for obtaining refined heuristic solutions for problems with electromagnetic wave diffraction by scatterers with non-ideal boundary conditions is proposed. Based on the method of fundamental components for the problem of TE-polarized wave diffraction on a “semitransparent” half-plane with two-sided impedance boundary conditions, an improved heuristic solution for the scattered field is obtained by replacing the primary heuristic formula. A quantitative and qualitative assessment of the accuracy of the heuristic relationship is carried out as well as a physical interpretation of the rigorous solution based upon them. A classification of heuristic solutions is proposed according to the purpose of their application in practical problems.

1. Introduction

The construction of heuristic formulas for various diffraction problems is a topical task. The problem is that for many diffraction problems, obtaining rigorous analytical solutions is a long and laborious task, the ultimate success of which is not guaranteed.
Unlike rigorous analytical formulas, heuristic formulas are based not on the mathematical formulation of the problem, but on knowledge of the physical features of the diffraction process. Heuristic formulas are needed to increase the performance of solvers designed to study the resource-intensive problems of diffraction, the physical interpretation of numerical solutions, the refinement of numerical solutions that differ from experimental data as a result of the inaccurate formulation of the problem, etc. Traditional heuristic methods such as the method of edge waves (MEW) and the geometric theory of diffraction (GTD) are widely known [1,2,3].
All heuristic formulas are verified at the stage of creation for all values of their input parameters. By verification, we mean a comparison with a reliable solution, usually a numerical one. In the process of verification, the limits of the application of heuristic formulas and their accuracy are established. Later on, heuristic formulas are employed without verification.
Traditional heuristic methods are based on field behavior postulates and have inherent methodological errors [4,5] that can lead to a lack of accuracy.
Recently, the method of fundamental components (MFC) has been developed [4,5]. Along with verification, in contrast to the traditional heuristic methods of GTD and MEW, MFC involves an adjustment procedure (correction of primary heuristic formulas based on the results of verification). The adjustment procedure minimizes the inherent methodological error and allows the provision of a predetermined accuracy needed for solving a given practical problem in the entire range of parameters.
MFC consists of three stages:
  • Obtaining a verification solution.
  • Construction of primary heuristic formulas.
  • Construction of tuning formulas.
With the possibility of performing two stages (1) and (2) (see above), the MFC approach is suitable for any type of boundary conditions and allows one to obtain heuristic formulas that can be used in solvers instead of numerical solutions.
The physical meaning of the technique used in MFC is that we single out the factors that have the greatest influence on the solution (singularities and the coefficient of reflection from an unbounded surface R), and introduce them into the primary heuristic formulas. Then, during the adjustment process, we obtain additional expressions that affect the accuracy and add them to the final heuristic formulas. Additional expressions may have a physical meaning and describe the features of the scatterer that are not taken into account in the primary heuristic formulas, or can simply be a mathematical correction that increases the accuracy of the solution. As a result, we obtain a refined heuristic analytical solution that can be used in the solver of a practical problem instead of a numerical one.
If, as in GTD and MEW, unaccounted physical factors are determined in advance, formalized, and limited to this, then the resulting expressions, as a rule, contain an inherent methodological error leading to inaccurate solutions and may need further clarification if this is required by the interests of a practical problem.
Heuristic formulas for a half-plane with nonideal boundary conditions can be built based on a combination of well-known GTD and MEW solutions for a perfectly conducting half-plane with reflection and transmission coefficients for an unbounded surface. The use of heuristic formulas instead of numerical solutions is a topical line of research. Performance improvement is especially effective in multiple calculations. Works [6,7,8,9,10] are based on the use of ready-made numerical–analytical solutions; refinement can be obtained by reducing the solution to a form that satisfies the reciprocity principle. However, the calculations made by the final heuristic formulas still differ from the verification solution.
The purpose of this work is to improve the technique for working with numerical solutions of reference problems, in particular, to find a way of improving accuracy, and not to obtain a specific engineering formula for specific boundary conditions and a specific practical problem. This article refines the solution at the edges of the range of angles by switching to another primary heuristic formula.
The statement of the problem in this work is to find heuristic expressions for diffraction on semitransparent objects that could be used in solvers instead of a numerical solution to increase performance. Such a problem is relevant when a numerical solution requires significant computer resources. The technique for obtaining a heuristic solution can be applied to translucent objects with different properties. In this work, in order to demonstrate the effectiveness of the applied technique, a half-plane with two-sided impedance boundary conditions is chosen as a semitransparent object [11,12,13]. To solve the problem of diffraction on a semitransparent half-plane, S.E. Bankov [14,15] created a numerical solution based on the formulation of the diffraction problem, formulated independently of the earlier works of T.B.A. Senior [16].

2. Boundary Conditions for the TE Polarization

Two-sided impedance boundary conditions for the case of TE-polarized electromagnetic wave diffraction on a half-plane (sometimes such a half-plane is called “semitransparent”) have the form [14,15,17]:
Z H z + H z = E x E x + = E x = E x w h e n y = 0 ,
where impedance Z = iX, X is a variable parameter on which the reflection coefficient depends, and i is an imaginary unit. Here, the signs “+” and “−” correspond to the regions y > 0 and y < 0, respectively (see Figure 1).
Figure 1 shows the geometry of a two-dimensional diffraction problem for the normal incidence of a TE-polarized electromagnetic wave on the edge of a semitransparent half-plane.
Two-sided impedance boundary conditions have existed in diffraction theory for quite a long time [11,12,13,17]. They should not be confused with the Schukin–Leontovich impedance boundary conditions at the interface between two media. Two-sided impedance boundary conditions are written for a half-plane that has no physical thickness but is characterized by impedance—a proportionality coefficient between the amplitudes of electric and magnetic fields, which is included in the mathematical expressions for these boundary conditions.
Real physical objects for which two-sided impedance boundary conditions can be applied include metal gratings and thin films, including those made of resistive material. Such structures are used to obtain antennas with desired characteristics [17,18].

3. Numerical Solution and Analytical Formulas

A rigorous approach [14] based on the Wiener–Hopf method (WHM) [19,20] makes it possible to find a verification solution (“WHM solution”):
g r X , ϕ , ϕ 0 .
Introducing the parameters (X, φ, φ0) into the calculation program for the WHM, we obtain the values of the solution gr(X, φ, φ0) (2), which we will further use to verify the heuristic solutions. In the first stage, the verification solution (2) does not have a compact analytical form, and our task is to find an analytical heuristic formula that could be employed instead of (2) in practical problems.
If the boundary conditions (1) are satisfied on the unbounded half-plane, then the reflection and transmission coefficients R and T can be found:
R X , ϕ 0 = W 0 W 0 + 2 i sin ϕ 0 1 X ,   W 0 = 120 π ,
T X , ϕ 0 = 1 R X , ϕ 0 .

4. Primary Heuristic Formulas

Primary heuristic formulas for describing the solution of the diffraction problem for a TE-polarized field on a “semitransparent” half-plane with two-sided impedance boundary conditions (1) have the form [5,15,21,22,23]:
g f R , T , ϕ , ϕ 0 = 1 2 1 T cos ϕ ϕ 0 2 + R cos ϕ + ϕ 0 2 ,
g f 0 R , T , ϕ , ϕ 0 = 1 2 1 T sin ϕ ϕ 0 2 cos ϕ ϕ 0 2 + R sin ϕ + ϕ 0 2 cos ϕ + ϕ 0 2 .
Here, R and T are the reflection and transmission coefficients for a flat unbounded surface (3) and (4). A similar form of heuristic solutions was used in [6,7,8,9,10,15,21,22,23,24,25].
Primary heuristic formulas (5) and (6) are combinations of reflection and transmission coefficients for a flat unbounded surface R and T and singular diffraction coefficients for a perfectly conducting half-plane [1,4,5,22]. Primary heuristic formulas describing solutions to diffraction problems on a semitransparent and nontransparent half-plane for a TH-polarized field are given in [5,21,26]. It should be noted that the primary heuristic formulas for a semitransparent or nontransparent half-plane are significantly different. The adjustment functions and the physical interpretation of the results differ in the same way. Primary heuristic formulas can be obtained both from a rigorous solution of the two-dimensional problem of diffraction using a perfectly conducting half-plane and from the solution of the same problem in the approximation of physical optics. Let us call them, respectively, the generalized diffraction coefficient (GDC) (5) and the generalized diffraction coefficient in the approximation of physical optics (PODC) (6) [4,5,15,21,22,23]. The ODC approximation has been known in diffraction theory for a long time [24,25]. However, its accuracy for a specific statement of the problem can be established only in the process of verification.
The verification solution gr(X, φ, φ0) (2) is constructed in such a way that it can be compared with the singular diffraction coefficients GDC (5) and PODC (6). From relationships (3) and (4) it follows that at the point X = 0: R = 1, T = 0, and the relationship gf(R, T, φ, φ0) = gr(X, φ, φ0) holds for all angles φ, φ0. Therefore, the WHM solution coincides with the corresponding solution (5), which in this case gives a rigorous solution to the problem.
For the implementation of the MFC algorithm, it is not critical whether we have analytical expressions for R and T. If necessary, they can be taken in the form of numerical solutions, especially when the edge profile and the reflecting surface have a complex structure and there are no analytical solutions. The important thing is that we apply a unified approach to all tasks in which we use R and T.

5. Obtaining Heuristic Formulas Employing the Adjustment Process

The study of the WHM solution gr(X, φ, φ0) (2) at 0 < X < ∞ (see, for example, [15,21,22,23]) showed that for small values of the “transparency” parameter X, the WHM solution (2) behaves like the GDC (5), and for large values of X, the WHM solution (2) behaves like the PODC (6). Taking into account the influence of the parameter X on the solution, we can introduce a new concept of “adjustment function” cxg(X, φ) and approximately describe it employing the phenomenological formula:
c x g X , ϕ = sin ϕ 2 x X ,   x X = 1 exp X W 0 .
Here, x(X) is a function of the parameter X (not to be confused with the x coordinate axis). As X→0 we have x(X)→0 and cxg(X, φ)→1. As X→∞, we have: x(X)→1 and cxg(X, φ)→sin(φ/2).
Employing (5) and (7), we can construct a heuristic formula gh(R, T, φ, φ0):
g h R , T , ϕ , ϕ 0 = g f R , T , ϕ , ϕ 0 c x g X , ϕ c x g X , ϕ s .
The formula for gh(R, T, φ, φ0) (8), which describes the behavior of the WHM solution better than each of the individual diffraction coefficients GDC (5) or PODC (6), is the desired heuristic approximation of the MFC and can be used in practical applications instead of the WHM solution gr(X, φ, φ0) (2). Formula (7) can be included in the list of fundamental components along with (5), (6), (3), (4), etc. [5]. It should be noted that the form of the fundamental component (7), if necessary, can be changed according to the accuracy required to solve a specific practical problem.
The physical meaning of Formula (8) is that the function gh(R, T, φ, φ0) changes when gg0 as the parameter X increases.
The results of the calculation employing the WHM solution gr(X, φ, φ0) (2) and heuristic formulas (8), (5), and (6) are shown in Figure 2.
Figure 2. The scattering diagrams (SD) modules for X = 300, x(X) = 0.549, and φ0 = 90 are plotted along the vertical axis. Circles—GDC (5), triangles—PODC (6), solid gray line—WHM solution gr(X,φ,φ0) (2), and squares—heuristic solution gh(R, T, φ, φ0) (8).
Figure 2. The scattering diagrams (SD) modules for X = 300, x(X) = 0.549, and φ0 = 90 are plotted along the vertical axis. Circles—GDC (5), triangles—PODC (6), solid gray line—WHM solution gr(X,φ,φ0) (2), and squares—heuristic solution gh(R, T, φ, φ0) (8).
Eng 04 00024 g002
In Figure 2, the modules of the scattering diagrams of the WHM solution gr(X, φ, φ0) (2), the heuristic solution gh(R, T, φ, φ0) (8), the GDC (5), and the PODC (6) are plotted along the vertical axis. The WHM solution moves between the GDC and the PODC, and the heuristic solution moves along with the WHM solution.

6. Qualitative Evaluation of Accuracy

Let us compare the scattering diagram of the WHM solution gr(X, φ, φ0) for the case of the incidence of a TE-polarized electromagnetic wave on a semitransparent half-plane at 0 < X < ∞ with the expressions (3) and (4) substituted into approximate solutions for GDC (5) and PODC (6). Let us introduce a qualitative criterion for estimating the accuracy of approximate solutions in the form of modules for the ratios of amplitudes in the WHM solution (2) and heuristic approximations GDC (5), PODC (6), and MFC (8), respectively:
a g R , T , ϕ , ϕ 0 = g r X , ϕ , ϕ 0 g f R , T , ϕ , ϕ 0 ,
a g 0 R , T , ϕ , ϕ 0 = g r X , ϕ , ϕ 0 g f 0 R , T , ϕ , ϕ 0 ,
a h R , T , ϕ , ϕ 0 = g r X , ϕ , ϕ 0 g h R , T , ϕ , ϕ 0 .
In the case of the best coincidence of one of the solutions with the WHM solution gr(X, φ, φ0) (2), the corresponding function (9)–(11) is located closer than the others with respect to the straight line on which the modulus of the amplitude ratio is equal to one.
The results of the heuristic formula calculations (9)–(11) are shown in Figure 3.
Figure 3 shows amplitude ratio modules for ag(R, T, φ, φ0) (9), ag0(R, T, φ, φ0) (10), and ah(R, T, φ, φ0) (11).
Figure 3. The ratios of the scattering diagram modules. The ratios of the scattering diagram (SD) modules for X = 300, x(X) = 0.549, and φ0 = 90 are plotted vertically. Circles—ag(R, T, φ, φ0) (9), triangles—ag0(R, T, φ, φ0) (10), squares—ah(R, T, φ, φ0) (11), and diamonds—ah0(R, T, φ, φ0) (14).
Figure 3. The ratios of the scattering diagram modules. The ratios of the scattering diagram (SD) modules for X = 300, x(X) = 0.549, and φ0 = 90 are plotted vertically. Circles—ag(R, T, φ, φ0) (9), triangles—ag0(R, T, φ, φ0) (10), squares—ah(R, T, φ, φ0) (11), and diamonds—ah0(R, T, φ, φ0) (14).
Eng 04 00024 g003

7. Adjustment Function for PODC

Formulas (1)–(11) were given in [15]. As a disadvantage of the MFC formula (8), one can note a reduced accuracy at the edges of the zone of observed angle change (Figure 3).
When analyzing Figure 3, the idea arises that, based on PODC, the heuristic formula can be tuned more accurately at the edges. The ratio ag(R, T, φ, φ0) (9) has a drop to zero at the edges of the angular range, while the ratio ag0(R, T, φ, φ0) (10) does not have such a drop. This can make it possible to more accurately select the parameters of the heuristic formula.
Let us build a new heuristic solution no longer based on the GDC (5), but on the basis of the PODC (6). To do this, we find a new “adjustment function” cxg0(X, φ). As necessary, we employ the technique of constructing engineering formulas employing phenomenological digital coefficients.
As the initial formula for constructing the “adjustment function” cxg0(X, φ), we take the expression cos ϕ / p x X . Next, we subject this expression to a series of transformations, choosing the size of the curve in height, the vertical and horizontal offsets, and the form of the dependence p[x(X)], which changes the periodicity of the original function with respect to the angle φ.
In addition, it is necessary to maintain the equality 1 for the “adjustment function” at the singularity points φs of expressions (5) and (6) (the angles φs = π ± φ0 correspond to the directions “straight ahead” and “mirror reflection”, see Figure 1), since at these points, the values of GDC and PODC are the same—they give the correct value for the scattered field, and no correction is required. The set of selected parameters should provide the best match between the heuristic formula and the WHM solution (2). As a result, we obtain a phenomenological heuristic function cxg0(X, φ), which gives the best match for the vast majority of x(X) values and other parameters:
g h 0 R , T , ϕ , ϕ 0 = g f 0 R , T , ϕ , ϕ 0 c x g 0 X , ϕ s c x g 0 X , ϕ ,
where
c x g 0 X , ϕ s c x g 0 X , ϕ = cos ϕ s p x X 12 + 12 p 0.85 x X cos ϕ p x X 12 + 12 p 0.85 x X , p x X = 0.52 + 2 1 + 10 x X .
The amplitude ratio modulus (see item 5) for formula (12):
a h 0 R , T , ϕ , ϕ 0 = g r ϕ , ϕ 0 g h 0 R , T , ϕ , ϕ 0 .
The results of the calculation from Formula (14) are also shown in Figure 3. It can be seen from the graphs that the heuristic solution ah0(R, T, φ, φ0) (14) gives the best accuracy with respect to the WHM solution since it is located closer than others relative to the unit line.

8. Classification of Heuristic Solutions

The MFC diffraction coefficient (8) can be considered a physically clear simplified form and used for the physical interpretation of the verification solution. For example, in this work, employing a relatively simple expression for the “adjustment function” (8), one can approximately trace the nature of the transformation of the singular diffraction coefficient GDC (5) into the PODC (6) as the parameter X changes from 0 to ∞.
The problem of physical interpretation can also arise when, as a result of a simplified formulation of the problem, a mathematically rigorous solution diverges from the experimental results.
There are practical tasks in which, for example, it is required to increase the performance of a resource-intensive solver. Then accuracy is more important than physical clarity, and one can rely more on phenomenological engineering approaches, complicate heuristic expressions, and apply adjustment formulas such as (13).

9. Quantification of Accuracy

A graphical comparison of amplitude ratios (Figure 3) is a qualitative way to evaluate accuracy. We introduce a quantitative criterion for estimating accuracy in the form of an integral over the area of change in the observation angle φ. The integrand will be the average deviation from the unity of the absolute values of the ratios of the amplitudes of the WHM solution (2) and qualitative evaluation criteria according to the GDC (9), PODC (10), MFC (11), and (14), respectively:
D g R , T , ϕ , ϕ 0 = 1 2 π ϕ = 0 2 π a g R , T , ϕ , ϕ 0 1 d ϕ ,
D g 0 R , T , ϕ , ϕ 0 = 1 2 π ϕ = 0 2 π a g 0 R , T , ϕ , ϕ 0 1 d ϕ ,
D h R , T , ϕ , ϕ 0 = 1 2 π ϕ = 0 2 π a h R , T , ϕ , ϕ 0 1 d ϕ ,
D h 0 R , T , ϕ , ϕ 0 = 1 2 π ϕ = 0 2 π a h 0 R , T , ϕ , ϕ 0 1 d ϕ .
The smaller the value of any of these integrals, the greater the accuracy of the corresponding approximation. Graphs of the values of the integrals for the MFC, GDC, and PODC approximations depending on the value of x(X) are shown in Figure 4.
It follows from the graphs of the average deviation that the best agreement with the WHM solution is given by the approximations MFC and (14). For small values of X (in this case x(X)~0), the GDC approximation is also close to the WHM solution, and for large values of X (with x(X)~1), the PODC approximation is closer to the WHM solution. When incident from the side of the end, when scattering occurs at small “grazing” angles to the half-plane, the accuracy of the coincidence of the MFC approximation with the WHM solution deteriorates. It follows that this case needs to be considered more carefully. It is possible that in the case of a three-dimensional structure of finite size, the MFC approximations (8) and (12) will work equally well in all directions; however, this assumption needs to be checked separately.

10. Reasons for Decreasing the Accuracy of the Solution Based on Physical Optics for Some Values of x(X)

Let us try to understand the reasons for the decrease in the accuracy of solution (12) for small x(X) (see Figure 4). Figure 5 shows the ratios of the modules of scattering diagrams (SD) for X = 30, x(X) = 0.076, and φ0 = 160°, corresponding to the situation when approximation (12) has a large error Dh0(R, T, φ, φ0) (18) (see Figure 4).
Figure 5. The ratios of the scattering diagram modules. The ratios of the scattering diagram (SD) modules for X = 30, x(X) = 0.076, and φ0 = 160º are plotted along the vertical axis. Circles—GDC (9), triangles—PODC (10), squares—ah(R, T, φ, φ0) (11), and diamonds—ah0(R, T, φ, φ0) (14).
Figure 5. The ratios of the scattering diagram modules. The ratios of the scattering diagram (SD) modules for X = 30, x(X) = 0.076, and φ0 = 160º are plotted along the vertical axis. Circles—GDC (9), triangles—PODC (10), squares—ah(R, T, φ, φ0) (11), and diamonds—ah0(R, T, φ, φ0) (14).
Eng 04 00024 g005
Observing the behavior of the graph of the function ah0(R, T, φ, φ0) (14) shown in Figure 5, it can be seen that its deviation from the unit line is caused by the imposition of the equality condition 1 at the singularity points. Therefore, to improve the accuracy, it is necessary to either complicate the expression corresponding to the formulas gh0(R, T, φ, φ0) (12) and cxg0(X, φ) (13) or simply use another function for small x(X), for example, the MFC expression corresponding to the formulas gh(R, T, φ, φ0) (8) and cxg(X, φ) (7).

11. Discussion of the Results

The purpose of this paper is not to study certain types of boundary conditions, but to develop a specific heuristic approach for the method of fundamental components [4,5]. This article demonstrates the possibility of increasing the accuracy of the heuristic solution by replacing the primary heuristic formula with another one that has a useful property (nonzero values at the ends of the angle range). Therefore, only one polarization is considered and there is no detailed comparison with other analytical approaches; this does not matter for the heuristic technique. On the other hand, a quantitative comparison with the results of calculations employing other heuristic approaches is given. The result of this paper is a heuristic solution of increased accuracy.
When constructing a heuristic formula employing the MFC, it is important to correctly choose the main points of the solution: (A) the primary heuristic formula; (B) the type of correction function (e.g., additive or multiplicative); and (C) the location of the phenomenological adjustment factors, etc. All of these steps should be performed based on solutions to the simplest problems, experience, and physical intuition.
In the process of constructing heuristic formulas, one can act in two ways. The first way is to first determine the unaccounted-for physical features, and then write an addition to the primary formula. An example of such an approach is an MEW containing an additive, namely, Ufimtsev’s edge currents. The second way is to first find the discrepancy between the primary heuristic formula and the verification solution, and then determine the physical meaning of this discrepancy. An example of such an approach is the MFC containing a factor—an adjustment function. The determination of the physical meaning of the discrepancy with the verification solution may require some time and additional research; this is a common situation for such science as physics.
Strict analytical solutions are not required to build heuristic formulas for MFC, which is a certain advantage that gives MFC flexibility and versatility compared to other approaches. Unlike well-known heuristic approaches (for example, GTD and MEW, which use analytical solutions), this method allows one to obtain heuristic formulas as a result of verification and adjustment procedures based on numerical solutions. In fairness, it should be noted that there are hybrid methods (which also do not require strict analytical formulas) that combine the MEW and numerical approaches. Hybrid methods, just like numerical ones, require an application package.
The traditional heuristic approaches GTD and MEW are based on “postulates” that require a priori assumptions and a certain sequence of actions to obtain a heuristic solution. A priori assumptions, such as “the field on a curved edge is the same as on a straight edge”, on one hand, make it easier to obtain a solution, but on the other hand, create a laid-down methodological error.
Unlike traditional heuristic methods, the MFC approach is not based on postulates but in comparison with a verification solution. Depending on the task, the heuristic formulas of the MFC can be different. More complex and accurate heuristic formulas are better suited for use in solvers instead of a numerical reference problem calculator. Simpler and more physically clear heuristic formulas are better suited for the physical interpretation of the numerical calculation results. If the numerical verification solution differs from the experimental data (for example, due to inaccuracies in the mathematical formulation of the problem), it is possible to adjust the heuristic formulas according to the experimental results.
First, the researcher must make a “strategic” decision about which type of heuristic formula is best suited for a given practical problem. Then, one can proceed to “tactical” actions: choose the most physically clear primary heuristic formula for the “adjustment function”, carry out an adjustment procedure employing an engineering formula with phenomenological coefficients (as in this work), or interpolate with splines, or construct a solution based on the primary formula, which satisfies the principle of reciprocity. A specific formula can be developed manually, as well as employing standard programs such as Wolfram Mathematica.
Primary heuristic formulas have a clear physical meaning. If you choose them correctly, as well as the type of adjustment formula, then an adjustment will be easy. Previously, the GDC approximation was chosen as the primary heuristic formula. In this article, approximation PODC is chosen and an adjustment function containing adjustment coefficients is added to it. These coefficients were found phenomenologically, employing a numerical procedure to ensure the minimum discrepancy between the verification solution and the heuristic function. As a result, the accuracy of the heuristic solution increased significantly, which was shown by quantitative evaluation.
As already mentioned, heuristic solutions are obtained based on physical considerations, which may be different. Therefore, for the same statement of the problem, several different heuristic solutions can be constructed. Usually, by saying “exact heuristic solution”, one of the following options is meant:
  • This heuristic solution is more accurate than other known heuristic solutions;
  • This heuristic solution has an accuracy that meets the requirements of the practical problem to which it is applied.
In some cases when employing the MFC, it is possible to obtain a very accurate heuristic analytical solution [27], but this does not detract from the results of this work since the described technique can be applied to a wider class of problems.
As is known, heuristic approaches are highly efficient and widely used. In particular, the issue of replacing numerical solutions with heuristic ones in solving practical problems is topical. Heuristic solutions for wedge diffraction with impedance boundary conditions [28] were obtained in [6,7,8,9]. The Malyuzhinets function is calculated for a rather long time, so obtaining heuristic solutions significantly speeds up the calculation. An exact match with the verification solution was not obtained, but the agreement with the experimental results was satisfactory [9].
Impedance boundary conditions [28] are formulated at the interface between two media. When solving the problem of diffraction on bodies with such boundary conditions, the penetration of a wave into a conducting medium is neglected and is limited to taking into account reflection; for diffraction by a scatterer with impedance boundary conditions T = 0, 0 < |R| < 1 [6,7,8,9].
Two-sided impedance boundary conditions on the half-plane lead to the fact that the field can pass through, 0 < |T| < 1. When T = 0 from (4), it follows that R = 1. Thus, the boundary conditions in this article are different from the impedance boundary conditions.

12. Conclusions

This article presents a technique for improving the accuracy of the solution by replacing the primary heuristic formula with another that does not have features at the boundaries of the range of angles. For the “adjustment function,” a formula containing phenomenological numerical coefficients is obtained. The technique has been tested on a certain types of boundary conditions and for the corresponding verification solution, but there is every reason to believe that it will work in other cases as well.
A quantitative comparison for the accuracy of MFC solutions showed that the accuracy of formula (12) obtained employing engineering methods, as expected, turned out to be higher than the accuracy of formulas (5), (6), and (8). It can be additionally increased, if necessary, for the purposes of a specific practical task. In some areas with a significant deviation of the “improved” curves, one can switch to the old solutions.

Funding

This research received no external funding.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The author is grateful to S.E. Bankov, the author of the programs and analytical formulas, the basis upon which heuristic formulas were verified and adjusted.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ufimtsev, P.Y. Method of Edge Waves in the Physical Theory of Diffraction; Sovetskoe Radio: Moscow, Russia; US Air Force Foreign Technology Division: Wright-Patterson Air Force Base, OH, USA, 1962; pp. 1–1154. [Google Scholar]
  2. Keller, J.B. Geometrical Theory of Diffraction. J. Opt. Soc. Am. 1962, 52, 116–130. [Google Scholar] [CrossRef] [PubMed]
  3. Borovikov, V.A.; Kinber, B.E. Geometrical Theory of Diffraction; IEE: London, UK, 1994. [Google Scholar]
  4. Vesnik, M.V. The method of the generalized eikonal. In New Approaches in the Diffraction Theory; Walter de Gruyter GmbN: Berlin, Germany; Boston, MA, USA, 2015; ISBN 978-3-11-031112-9. [Google Scholar]
  5. Vesnik, M.V. Physical Interpretation of Strict Solutions of Diffraction Problems by Heuristic Relations. J. Math. Sci. 2019, 239, 751–770. [Google Scholar] [CrossRef]
  6. Luebbers, R.J. Finite conductivity uniform GTD versus knife edge diffraction in prediction of propagation path loss. IEEE Trans. Antennas Propag. 1984, AP-32, 70–76. [Google Scholar] [CrossRef]
  7. Holm, P. A new heuristic UTD diffraction coefficient for nonperfectly conducting wedges. IEEE Trans. Antennas Propag. 2000, 48, 1211–1219. [Google Scholar] [CrossRef]
  8. El-Sallabi, H.M.; Rekanos, I.T.; Vainikainen, P. A New Heuristic Diffraction Coefficient for Lossy Dielectric Wedges at Normal Incidence. IEEE Antennas Wirel. Propag. Lett. 2002, 1, 165–168. [Google Scholar] [CrossRef]
  9. Soni, S.K.; Bhattacharya, A. New heuristic diffraction coefficient for modeling of wireless channel. Prog. Electromagn. Res. C 2010, 12, 125–137. [Google Scholar] [CrossRef]
  10. Klionovski, K.; Bankov, S. Electromagnetic wave scattering by a half-plane with generalized semitransparent boundary conditions. Waves Random Complex Media 2020, 32, 1–12. [Google Scholar] [CrossRef]
  11. Trentini, G. Partially reflecting sheet arrays. IEEE Trans. Antennas Propag. 1956, 4, 666–671. [Google Scholar] [CrossRef]
  12. Senior, T.B.A. Diffraction by a semi-infinite metallic sheet. Proc. Roy. Soc. A 1952, 213, 436. [Google Scholar]
  13. Senior, T.B.A. Half plane edge diffraction. Radio Sci. 1975, 10, 645. [Google Scholar] [CrossRef]
  14. Bankov, S.E. Integrated Microwave Optics; Fizmatlit: Moscow, Russia, 2018. [Google Scholar]
  15. Vesnik, M.V.; Bankov, S.E. Heuristic solution to the problem of diffraction of a TE-polarized electromagnetic wave on a semitransparent half-plane. Waves Random Complex Media 2021, 1–12. [Google Scholar] [CrossRef]
  16. Senior, T.B.A. Diffraction tensors for imperfectly conducting edges. Radio Sci. 1975, 10, 911. [Google Scholar] [CrossRef]
  17. Voitovich, N.N.; Katsenelenbaum, B.Z.; Korshunova, E.N.; Pangonis, L.I.; Pereyaslavets, M.L.; Sivov, A.N.; Shatrov, A.D. Electromagnetics of Antennas with Semitransparent Surfaces; A Method of Constructive Synthesis; Nauka: Moscow, Russia, 1989. [Google Scholar]
  18. Tatarnikov, D.V. Semi-Transparent Ground Planes Excited by Magnetic Line Current. IEEE Trans. Antennas Propag. 2012, 60, 2843–2852. [Google Scholar] [CrossRef]
  19. Noble, B. Methods based on the Wiener—Hopf Technique for the Solution of Partial Differential Equations; London, UK, 1958.
  20. Nefedov, E.I. Diffraction of Electromagnetic Waves on Dielectric Structures; Nauka: Moscow, Russia, 1979. (In Russian) [Google Scholar]
  21. Vesnik, M.V. Using the Method of Base Components for a Heuristic Solution to the Diffraction Problem on a Half-Plane with Nonideal Boundary Conditions. J. Commun. Technol. Electron. 2019, 64, 1211–1217. [Google Scholar] [CrossRef]
  22. Vesnik, M.V. New Possibilities for Constructing Heuristic Solutions to Problems of Electromagnetic Diffraction. Eng 2022, 3, 27–41. [Google Scholar] [CrossRef]
  23. Vesnik, M. Physical Interpretation of the Solution to the Problem of Diffraction on a Half-plane with Non-Ideal Boundary Conditions. J. Eng. Res. Sci. 2022, 1, 52–58. Available online: https://www.jenrs.com/v01/i03/p006/ (accessed on 23 December 2022). [CrossRef]
  24. Raman, C.V.; Krishnan, K.S. The diffraction of light by metallic screens. Proc. R. Soc. London. Ser. A Math. Phys. Sci. 1927, 116, 254–267. [Google Scholar] [CrossRef]
  25. Shmoys, J. Diffraction by a half-plane with a special impedance variation. IEEE Trans. Antennas Propag. 1959, 7, 88–90. [Google Scholar] [CrossRef]
  26. Bankov, S.E.; Vesnik, M.V.; Kravchenko, V.F. Heuristic solution to the diffraction problem on a superconducting half-plane. J Commun. Technol. Electron. 2020, 65, 398–405. [Google Scholar] [CrossRef]
  27. Vesnik, M.V. A technique for obtaining analytical heuristic solutions in problems of diffraction on two-dimensional semi-infinite objects with non-ideal boundary conditions. Waves Random Complex Media 2022, 1–17. [Google Scholar] [CrossRef]
  28. Maliuzhinets, G.D. Excitation, reflection and emission of surface waves from a wedge with given face impedances. Sov. Phys. Dokl. 1958, 3, 752–755. [Google Scholar]
Figure 1. The geometry of the problem.
Figure 1. The geometry of the problem.
Eng 04 00024 g001
Figure 4. Quantitative characteristics of the deviation. Graphs that quantitatively describe the deviations of heuristic approximations from a rigorous solution for various angles of incidence φ0. The x(X) values are plotted horizontally. The amplitudes of the functions describing the mean deviation are plotted along the vertical axis. Circles—GDC (15), triangles—PODC (16), squares—Dh(R, T, φ, φ0) (17), and diamonds—Dh0(R, T, φ, φ0) (18).
Figure 4. Quantitative characteristics of the deviation. Graphs that quantitatively describe the deviations of heuristic approximations from a rigorous solution for various angles of incidence φ0. The x(X) values are plotted horizontally. The amplitudes of the functions describing the mean deviation are plotted along the vertical axis. Circles—GDC (15), triangles—PODC (16), squares—Dh(R, T, φ, φ0) (17), and diamonds—Dh0(R, T, φ, φ0) (18).
Eng 04 00024 g004
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

Vesnik, M.V. Refinement of the Heuristic Solution for the Problem of TE-Polarized Electromagnetic Wave Diffraction on a Half-Plane with Two-Sided Impedance Boundary Conditions. Eng 2023, 4, 404-415. https://doi.org/10.3390/eng4010024

AMA Style

Vesnik MV. Refinement of the Heuristic Solution for the Problem of TE-Polarized Electromagnetic Wave Diffraction on a Half-Plane with Two-Sided Impedance Boundary Conditions. Eng. 2023; 4(1):404-415. https://doi.org/10.3390/eng4010024

Chicago/Turabian Style

Vesnik, Michael V. 2023. "Refinement of the Heuristic Solution for the Problem of TE-Polarized Electromagnetic Wave Diffraction on a Half-Plane with Two-Sided Impedance Boundary Conditions" Eng 4, no. 1: 404-415. https://doi.org/10.3390/eng4010024

Article Metrics

Back to TopTop