New Possibilities for Constructing Heuristic Solutions to Problems of Electromagnetic Diffraction

The paper formulates the foundations of a recently developed approach, named the method of fundamental components, intended for constructing heuristic solutions in problems of electromagnetic diffraction, for the first time. The difference between the new method and the known heuristic approaches lies in the application of an adjustment procedure that increases the accuracy. The possibility of the mentioned method for obtaining new results is illustrated with the help of the author’s previously published works. The advantages of the new method in constructing high-speed solvers and in the physical interpretation of numerical solutions are shown.


Introduction
The problems of electromagnetic diffraction can be solved in different ways. Mathematically rigorous formulations of the boundary value problems include those for which the existence and uniqueness of the solution have been proven. Unlike mathematically rigorous solutions (analytical or numerical), heuristic approaches are based on physical ideas about the structure and features of the solution. Heuristic solutions play an important role for the construction of high-performance electromagnetic diffraction solvers used in urgent practical problems. At the stage of deriving, all heuristic approaches need verification, i.e., revision using a reliable solution, usually a numerical one. At the stage of application, verification is no longer needed.

Relevance of Heuristic Solutions
Historically, the theory of diffraction was developed in line with the heuristic [1][2][3] approaches. In this regard, the works of Huygens, Jung, Fresnel, Helmholtz, Kirchhoff and others can be noted. In the works of Poincaré and Sommerfeld, it was first shown [1] that the solution of problems in the theory of diffraction can be found using the theory of boundary value problems of mathematical physics. This fruitful approach sometimes leads to very difficult solutions. Thus, despite the obvious importance and interest of researchers in a solution to the problem of diffraction by a plane perfectly conducting angular sector, the corresponding solution was obtained more than 100 years later than the solution of the problem of diffraction by a perfectly conducting half-plane [3,4]. Therefore, the development of heuristic methods is still an urgent task.
The need to use heuristic solutions is due to a number of practical problems, such as, for example, diffraction on radar objects, including objects with low radar signature, radio wave propagation in urban environments, etc.
A number of problems can be successfully solved using heuristic approaches, as follows: -Investigation of relatively large objects that do not have solutions obtained by more rigorous methods. Such objects include, for example, targets with low radar signature, Eng 2022, 3 28 non-axisymmetric satellite communication antennas, etc. The propagation of radio waves in urban environments may also be studied using heuristic methodologies. -Improving performance in calculating urgent problems, including those problems that have solutions obtained using other approaches. -Physical interpretation of numerical solutions, which is especially important in cases where they differ from the experimental results.

Known Heuristic Approaches
The following heuristic approaches are known: the method of geometrical optics (GO) [5], the method of physical optics (PO) [5], the geometrical theory of diffraction (GTD) [6][7][8] and the method of edge waves (MEW) [9][10][11]. The term physical theory of diffraction (PTD) can be used for the MEW and, in an extended interpretation, for all heuristic methods.
In heuristic approaches, a balance should be obtained between complexity and accuracy. Rigorous analytical decisions are accurate, but often very complex. The PO approximation can be obtained for any object, but the accuracy may not be enough. However, sometimes, the accuracy of the PO approximation is quite sufficient, as in the case of calculating the amplitude of the main beam and the first side lobe of a reflector antenna.
For relatively large objects, strict calculation is impossible due to limited computer resources. While calculating the scattering of signals on conventional radar targets, it used to be enough to take into account only the contribution of specular points; now, in problems of diffraction on targets with a reduced radar signature, they mainly use the MEW. The same applies to the calculation of propagation in an urban environment. The success of this or that approximation depends on the accuracy required in a particular practical problem.
Heuristic solutions are based on hypotheses, or "postulates" [8]. For example, the GTD postulate says that when scattered by an object of finite dimensions, the field of a finite size edge at the observation point is considered to be the same as if the incident ray field was reflected from an infinite edge. On the infinite edge, you can always find the point of the stationary phase at which the reflection occurs. At the edge of finite dimensions, the stationary phase point does not always exist; therefore, using the GTD, no solution at all can be obtained. This circumstance was the reason for the choice of the MEW (instead of the GTD) for calculating of low radar signature objects made using the "Stealth" technology [10].
The MEW postulate is that for a scatterer of finite dimensions, the field at the finite size edge is considered the same as at the edge of a semi-infinite scatterer with an edge of the same profile. In contrast to the GTD, a solution for the MEW (as well as for the PO method) can always be obtained (although sometimes with an inaccurate amplitude). The amplitude inaccuracies are caused by the fact that the field at the edge of finite dimensions is considered constant and the additional disturbance at the ends of the edge is neglected. In order to take this perturbation into account, it is necessary to investigate the field on a plane angular sector, for example, by numerical methods [12].

Reference Problems for Constructing Heuristic Solutions
When solving practical problems, both for semi-infinite objects and for objects of finite size, 1D, 2D and 3D reference problems are used.
Types of diffraction reference problems are the following: • Linear integral (3D): This is a solution in the PO approximation to the diffraction problem on a plane perfectly conducting scatterer. In this case, the far zone condition (FZC) is fulfilled both for the source and for the observation point. The integral over the area of the scatterer is reduced to an integral over its contour [13][14][15]. If the FZC is not fulfilled, the technique described in [16] can be applied.

•
Vertex waves (3D): This is a solution to the diffraction problem for a plane angular sector. The analytical solution is described in [17], the heuristic ones are described in [18][19][20][21].
Reference problems can be applied in the form of strict analytical formulas, heuristic formulas, engineering formulas (which may not have physical meaning), within the framework of hybrid approaches (combining elements of numerical and heuristic methods) or in the form of databases (pre-calculated or experimental) [3].

Motivation to Develop New Heuristic Approaches
The impossibility of using one of the well-known heuristic approaches in solving a practical problem is determined by one of the following factors: • Insufficient accuracy of known heuristic approaches. • Insufficient performance of the solver.

•
Lack of the required number of rigorous analytical solutions, on the basis of which, a solver can be built for a specific practical problem.

•
The need for a physical interpretation of the numerical solution and the identification of physical phenomena that were not taken into account by the old heuristic approaches.
These problems can be overcome by developing a new heuristic approach and constructing refined formulas.

New MFC Approach and Its Differences from Traditional Approaches
Suppose there is a specific practical problem for which one needs to construct an efficient (accurate and high-performance) solver. Suppose there are no analytical formulas for the corresponding reference problems, but there are results of numerical calculation of solutions of reference problems or experimental data. This is a common situation. With the help of a recently developed approach-the method of fundamental components (MFC) [3,4]-this situation can be used to construct efficient (compact and accurate) heuristic solutions.
At the first stage in the MFC, primary heuristic formulas are built. They are then subjected to an adjustment procedure. By this term, we mean the correction of primary formulas by functions (which we call fundamental components), based on the comparison of primary formulas with a verification solution.
In the process of adjustment, the contributions of known physical factors are consistently taken into account. Based on the difference between a strict and a heuristic solution, one can find the influence of unaccounted phenomena. In this case, we apply the adjustment procedure. One can act in two ways, as follows: • first, determine the unaccounted phenomenon and find an expression for it, then add it to the primary heuristic formula and carry out verification (as in the MEW); • or, first, find the difference between the strict solution and the primary heuristic formula, and then find a physical phenomenon that determines this difference (as in the MFC).
This difference significantly increases the number of problems for which heuristic formulas can be constructed. One can obtain a predetermined accuracy, as well as work with formulas for which there are no analytical solutions. So, the MFC has the property of universality.
As the primary heuristic formulas for any practical problem, one can choose the solutions of the simplest diffraction problems, or the heuristic formulas constructed earlier.
The fundamental components can be selected for both physical and mathematical reasons. Earlier, attempts have been also made to refine the primary heuristic formulas. It is possible to apply the multiple re-reflection of the GTD solution at the edges [18,19] or the PTD [10], and also apply the reciprocity principle [31][32][33][34]. Sometimes these approaches work well, sometimes they do not. In addition, solutions based on multiple reflections can be cumbersome. It is not always possible to reduce them into compact formulas.

Partitioning the Formula for a 2D Edge
The book [3] contains a general formula for the field V(ϕ) at the observation point (r,ϕ), excited by a source located at the point (r 0 ,ϕ 0 ) and scattered by a perfectly conducting wedge with an external aperture angle πn at normal incidence on the TH-or TE-polarized electromagnetic wave: where: Here, P(w sm ) is the source field at the saddle point w sm (from the geometric point of view, the saddle point corresponds to the situation when the edge is on the straight line connecting the observation point with the source, with ϕ ± ϕ 0 = π), and the eikonal S(ψ) and the eikonal at the saddle point S(w sm ) are equal, as follows: The area external to the wedge occupies the space of angles 0 < ϕ < πn. The geometry of the problem is shown in Figure 1.
formulas can be constructed. One can obtain a predetermined accuracy, as well as work with formulas for which there are no analytical solutions. So, the MFC has the property of universality.
As the primary heuristic formulas for any practical problem, one can choose the solutions of the simplest diffraction problems, or the heuristic formulas constructed earlier. The fundamental components can be selected for both physical and mathematical reasons.
Earlier, attempts have been also made to refine the primary heuristic formulas. It is possible to apply the multiple re-reflection of the GTD solution at the edges [18,19] or the PTD [10], and also apply the reciprocity principle [31][32][33][34]. Sometimes these approaches work well, sometimes they do not. In addition, solutions based on multiple reflections can be cumbersome. It is not always possible to reduce them into compact formulas.

Partitioning the Formula for a 2D Edge
The book [3] contains a general formula for the field V(φ) at the observation point (r,φ), excited by a source located at the point (r0,φ0) and scattered by a perfectly conducting wedge with an external aperture angle πn at normal incidence on the ТН-or TE-polarized electromagnetic wave: where: Here, P(wsm) is the source field at the saddle point wsm (from the geometric point of view, the saddle point corresponds to the situation when the edge is on the straight line connecting the observation point with the source, with φ ± φ0 = π), and the eikonal S(ψ) and the eikonal at the saddle point S(wsm) are equal, as follows: ( The area external to the wedge occupies the space of angles 0 < φ < πn. The geometry of the problem is shown in Figure 1.  Let us consider the integral representation for the scattered field in the case of two saddle points [3] and its particular case for diffraction by a wedge-Expression (2).
Each of the two terms on the right-hand side of Expression (2) consists of four factors. The first factor is as follows: Eng 2022, 3 31 which does not depend on the angular variable and is the product of the field value at the saddle point (i.e., on the shadow boundary) by the factor that determines the dependence of the solution on the distance to the source and the observation point.
The second factor, as follows: represents half of the diffraction coefficient. The total diffraction coefficient for a certain type of polarization is obtained by adding or subtracting the values of this factor at the observation point, as follows: Product of the third and fourth factors, as follows: is the quotient of dividing the Fresnel integral by its asymptotics and characterizes the dependence of the field on the angular distance to the "light-shadow" boundary (or, that which is the same, the "shadow" boundary). Far from the shadow boundary, this product is equal to 1. At the shadow boundary, this factor is zero and compensates for the singularity of half of the diffraction coefficient (5). This compensation results in a field equal to half of the field of geometrical optics at the shadow boundary of a semi-infinite scatterer. By changing the type of factors (4)- (7), or the parameters included in them, it is possible to construct a variety of heuristic solutions: for two-dimensional or three-dimensional scatterers, finite or infinite size, with or without fulfillment of the far-field condition, for different types of boundary conditions and edge profile.
More information on this topic can be found in [3].

Linear Integral in the PO Approximation
The linear integral in the PO approximation is used as the basis for the 3D solution of the diffraction problem.
Consider a perfectly conducting plane (with a normal vector → n ) scatterer excited by a plane wave (time dependence exp ik → n , → r ), as follows: In the PO approximation, the vector potential is defined as follows [3]: ρ is the radius vector of the point on the surface of the scatterer, R is the distance from the center of coordinates to the observation point, S is the square of the scatterer. Without loss of generality, we will assume that the plane scatterer is located in the XOY plane. In addition, let the far-field condition r ∼ = R − → n , → ρ [3,5,21] be satisfied, and then: where  depends on the selected type of polarization. The geometry of the problem is shown in Figure 2.
In the PO approximation, the vector potential is defined as follows [3]: where ρ   − ′ ′ = R n r is the distance between the point on the scatterer and the observation point, ρ  is the radius vector of the point on the surface of the scatterer, R is the distance from the center of coordinates to the observation point, S is the square of the scatterer. Without loss of generality, we will assume that the plane scatterer is located in the XOY plane. In addition, let the far-field condition [3,5,21] be satisfied, and then: where is the projection onto the scatterer surface (in our case, onto the XOY plane) of the difference ( ) between the direction vectors of the incident wave and the observation point. Applying Stokes' theorem, one can obtain the following [3,14]: where i n  is the unit internal normal to the contour C surrounding the scatterer, ρ ′  is the unit vector tangential to the contour and t is the coordinate measured along the contour. Let us denote the phase of the integration point in the exponent, as follows: Applying Stokes' theorem, one can obtain the following [3,14]: where → n i is the unit internal normal to the contour C surrounding the scatterer, → ρ is the unit vector tangential to the contour and t is the coordinate measured along the contour. Let us denote the phase of the integration point in the exponent, as follows: If the contour is a polygon with N vertices, then: Eng 2022, 3 33 where Φ j = k → ∆, → ρ j is the phase of the signal of the j-th vertex with the direction vector, which is the length of the j-th side of the polygon (located between the j − 1 th and j-th vertices).
If the observation point is located on the diffraction cone, when → ∆, → ρ = 0, then: where (ϕ 0 + π) and ϕ are the angles between the projections of the direction vectors where Φ j−1 = Φ j . Formula (13) are used as a basis for constructing heuristic formulas both for plane three-dimensional objects of finite size and semi-infinite objects.
More information on this topic can be found in [3].

Diffraction Coefficients for a Perfectly Conducting Half-Plane F and G with a Geometry Pattern
For a perfectly conducting half-plane, the expressions in (6) take the well-known [9] form of singular (at ϕ ± ϕ 0 = π) diffraction coefficients, as follows: Taking into account the relation [1], as follows: where ν(ψ) is determined by Expression (2), one can construct singular diffraction coefficients for the PO approximation, as follows:

Modifying Functions for a Perfectly Conducting Plane Angular Sector
An analytical solution to the problem of perfectly conducting quarter-plane diffraction appeared only a few years ago [17]. The geometry of the problem is shown in Figure 3. The polarization transmission coefficients D are determined by the expressions in [17], as follows: In [3,20,21], a heuristic solution to the problem of diffraction by a perfectly conducting plane angular sector is presented. In this case, the principle of "imaginary edge" was used. The imaginary edge is the direction in space in relation to which the directions to the source and to the observation point are disposed on the diffraction cone. A diffraction cone is a set of directions that make up the same angles with the edge [3,9].

Modifying Functions for a Perfectly Conducting Plane Angular Sector
An analytical solution to the problem of perfectly conducting quarter-plane diffraction appeared only a few years ago [17]. The geometry of the problem is shown in Figure  3.
The polarization transmission coefficients D are determined by the expressions in [17], as follows: In [3,20,21], a heuristic solution to the problem of diffraction by a perfectly conducting plane angular sector is presented. In this case, the principle of "imaginary edge" was used. The imaginary edge is the direction in space in relation to which the directions to the source and to the observation point are disposed on the diffraction cone. A diffraction cone is a set of directions that make up the same angles with the edge [3,9].  In Figure 4, the dotted line shows an imaginary edge making an angle γ/2 with the −z direction, while the real edge makes an angle β/2 with this direction. The currents of physical optics are designated as j ϕ (a) and j ϑ (b), depending on the types of polarization corresponding to Figure 3a,b. Components of currents directed perpendicularly and parallel to the imaginary edge are designated cr and cp, respectively.   Figure 5 illustrates the application of an imaginary edge. Figure 5a shows a scheme for solving 2D problem for diffraction of a plane wave at an infinite edge. Waves are scattered only in the directions of the diffraction cones. On the diffraction cones, the GTD and the MEW solutions coincide.
The wave vector of a wave incident on an edge can be divided into two components: one going along the edge, the other is perpendicular to it. This decomposition is shown in Figure 5c. Figure 5b shows an infinite edge, source position and viewpoint that are outside the diffraction cone, relative to the source. Diffraction cones are also shown, corresponding to situations where the source and observation point are swapped. If the positions of the   The wave vector of a wave incident on an edge can be divided into two components: one going along the edge, the other is perpendicular to it. This decomposition is shown in Figure 5c. Figure 5b shows an infinite edge, source position and viewpoint that are outside the diffraction cone, relative to the source. Diffraction cones are also shown, corresponding to situations where the source and observation point are swapped. If the positions of the source and observation points are interchanged, then decomposition becomes different, which leads to a violation of the principle of reciprocity. This violation has to be fought in the classical version of the MEW [9].
In Figure 5d a lilac dotted line shows an imaginary edge. This is the direction on the surface of the scatterer, determined so that the source and receiver are disposed on the diffraction cone in relation to the imaginary edge. Thus, for a solution using an imaginary edge, the reciprocity principle is fulfilled.  Figure 5 illustrates the application of an imaginary edge. Figure 5a shows a scheme for solving 2D problem for diffraction of a plane wave at an infinite edge. Waves are scattered only in the directions of the diffraction cones. On the diffraction cones, the GTD and the MEW solutions coincide.
The wave vector of a wave incident on an edge can be divided into two components: one going along the edge, the other is perpendicular to it. This decomposition is shown in Figure 5c. Figure 5b shows an infinite edge, source position and viewpoint that are outside the diffraction cone, relative to the source. Diffraction cones are also shown, corresponding to situations where the source and observation point are swapped. If the positions of the source and observation points are interchanged, then decomposition becomes different, which leads to a violation of the principle of reciprocity. This violation has to be fought in the classical version of the MEW [9].
In Figure 5d a lilac dotted line shows an imaginary edge. This is the direction on the surface of the scatterer, determined so that the source and receiver are disposed on the diffraction cone in relation to the imaginary edge. Thus, for a solution using an imaginary edge, the reciprocity principle is fulfilled.  The concept of an imaginary edge unites the approaches of the GTD and the MEW, since the GTD solution on an imaginary edge is used in order to obtain the MEW solution in directions outside the diffraction cones.
A more complete set of formulas describing the heuristic solution of the diffraction problem on a plane angular sector, as well as the calculation results, are given in [3,20,21] and in [35][36][37]. Calculated data from [17] were used as a verification solution. In this case, we used the linear integral in the PO approximation (13) and the diffraction coefficients (16) and (17), which replaced the diffraction coefficients (19) and (20), contained in the PO expression for the vector potential (10) [3]. Additionally, (in accordance with the principles of the MFC), additional multipliers (modifying functions) obtained as a result of the adjustment were applied, as follows: Eng 2022, 3 36 Formula (22) contain the following functions: The index "γ" means that the diffraction coefficients were taken on an imaginary edge located at an angle γ/2 to the coordinate axis (Figure 4).
The use of modifying functions clarifies the heuristic solution of the MFC in comparison with the MEW. The calculation results are shown in Figure 6.
In Figure 6, the worst coincidence with the strict solution [17] (solid gray line) gives the PO approximation, since the D coefficients for the cross polarizations Dϕϑ and Dϑϕ vanish. The best match with the strict solution is given by the modified formula in (22). The curves for the D coefficients calculated using the indicated heuristic approaches are shown by red dots and blue dashes next to the strict solution curves.

Diffraction Coefficients GDC (Generalized Diffraction Coefficients) and PODC (Physical Optics Diffraction Coefficients) for a Semitransparent Half-Plane
The boundary conditions for the impedance half-plane were considered in [38][39][40][41][42][43]. In papers [4,44], using the MFC, a heuristic solution to the problem of diffraction of an electromagnetic TH-a polarized wave on a half-plane with impedance boundary conditions, such as a thin layer-was obtained. In [45], a heuristic solution of the same problem was obtained for the case of TE polarization. The concept of an imaginary edge unites the approaches of the GTD and the MEW, since the GTD solution on an imaginary edge is used in order to obtain the MEW solution in directions outside the diffraction cones.
A more complete set of formulas describing the heuristic solution of the diffraction problem on a plane angular sector, as well as the calculation results, are given in [3,20,21] and in [35][36][37]. Calculated data from [17] were used as a verification solution. In this case, we used the linear integral in the PO approximation (13) and the diffraction coefficients (16) and (17), which replaced the diffraction coefficients (19) and (20), contained in the PO expression for the vector potential (10) [3]. Additionally, (in accordance with the principles of the MFC), additional multipliers (modifying functions) obtained as a result of the adjustment were applied, as follows: The index "γ" means that the diffraction coefficients were taken on an imaginary edge located at an angle γ/2 to the coordinate axis ( Figure 4).
The use of modifying functions clarifies the heuristic solution of the MFC in comparison with the MEW. The calculation results are shown in Figure 6.  (23) for a plane angular sector. Grey lines correspond to solution [17]. The angle φ is plotted along the horizontal axis (see Figure 1). The vertical axis shows the amplitudes of the D coefficients characterizing the relationship between the polarizations of the incident and scattered fields (21)- (23). Dφϑ: . transmission coefficients D from [17] are determined in (21).,. Horizontal lines designate the angle φ as in Figure 3.
In Figure 6, the worst coincidence with the strict solution [17] (solid gray line) gives the PO approximation, since the D coefficients for the cross polarizations Dφϑ and Dϑφ vanish. The best match with the strict solution is given by the modified formula in (22). The curves for the D coefficients calculated using the indicated heuristic approaches are shown by red dots and blue dashes next to the strict solution curves.

Diffraction Coefficients GDC (Generalized Diffraction Coefficients.) and PODC (Physical Optics Diffraction Coefficients) for a Semitransparent Half-Plane
The boundary conditions for the impedance half-plane were considered in [38][39][40][41][42][43]. In papers [4,44], using the MFC, a heuristic solution to the problem of diffraction of an electromagnetic TH-a polarized wave on a half-plane with impedance boundary conditions, such as a thin layer-was obtained. In [45], a heuristic solution of the same problem was obtained for the case of TE polarization.
Primary heuristic formulas are constructed on the basis of diffraction coefficients (16), (17), (19) and (20), adding to these expressions the reflection, R, and transmission, T, coefficients for an unbounded plane surface. The subscripts "TH" and "TE" indicate the type of polarization when the corresponding vectors are perpendicular to the edge. As the degree of translucency of the half-plane changes from opaque to completely transparent, the moduli of the reflection coefficients R change from 1 to 0.
We will call the resulting diffraction coefficients the GDC, as follow: (25) and the PODC, as follow:  TH  TH  TH  TH   R  T  T (23) for a plane angular sector. Grey lines correspond to solution [17]. The angle ϕ is plotted along the horizontal axis (see Figure 1). The vertical axis shows the amplitudes of the D coefficients characterizing the relationship between the polarizations of the incident and scattered fields (21)- (23). Dϕϑ: transmission coefficients D from [17] are determined in (21).,. Horizontal lines designate the angle ϕ as in Figure 3.
Primary heuristic formulas are constructed on the basis of diffraction coefficients (16), (17), (19) and (20), adding to these expressions the reflection, R, and transmission, T, coefficients for an unbounded plane surface. The subscripts "TH" and "TE" indicate the type of polarization when the corresponding vectors are perpendicular to the edge. As the degree of translucency of the half-plane changes from opaque to completely transparent, the moduli of the reflection coefficients R change from 1 to 0.
We will call the resulting diffraction coefficients the GDC, as follow: (25) and the PODC, as follow: Verification of these primary formulas showed that when the degree of transparency changes, formulas (24) and (25) transform to formulas (26) and (27) [44,45].

List of Fundamental Components
At the moment, we propose to include the following formulas in the list of fundamental components: 1.

2.
Linear integral with emphasis on the contributions of individual edges (13).

3.
R and T coefficients for an unbounded flat surface.

4.
Diffraction coefficients with R and T GDC (24), (25) and PODC (26), (27). Specific values of R and T are determined by a given type of boundary conditions. 5.
As new tasks are solved, the list of fundamental components will increase. The resulting heuristic formulas can be selected as primary ones for new tasks.

Prospects for the Application of New Heuristic Approaches
The difficulties in solving practical problems of diffraction by polygons and polyhedra are as follows: • Numerical solutions require large computer resources, software and programmers. Modern computers cannot cope with some tasks, so one has to use heuristic approaches. • Rigorous analytical solutions are lacking for most of the vast variety of edges that do not have analytical solutions.

•
The need to improve the performance of solvers for the study of specific types of tasks, such as, for instance, inverse problems.
In connection with the above, the MFC provides new opportunities for solving practical problems. New compact and at the same time exact (i.e., effective) heuristic formulas can be used for physical interpretation of numerical solutions and for constructing highspeed solvers.

Discussion
With the help of the MFC, on the basis of numerical solutions, it is possible to create heuristic analytical formulas with different properties. All heuristic formulas are fast. More compact and physically clear formulas can be used for a qualitative description of a numerical solution and its physical interpretation, more complex and accurate formulas are better suited for creating computers designed to solve resource-intensive problems.
One can set the accuracy depending on the purposes of the calculation. In the case of a study of propagation in urban environment, a less accurate approximation can be more appropriate. On the contrary, in the case of radar problems, it is better that the approximate (heuristic) solution is as accurate as possible, and this can be achieved using the adjustment procedure.
The novelty of the MFC in comparison with the PO, GTD and MEW lies in the way in which heuristic formulas are constructed. In the MFC, on the basis of the primary heuristic formulas GDC and PODC, it is possible to isolate and exclude from further consideration the determining factors, such as the influence of singularities at the "lightshadow" boundaries and the coefficient of reflection from an unbounded surface. Then, they focus on the adjustment procedure. As a result, within a specific formulation of the problem, it is possible to identify and describe the subtle physical features of the diffraction process, which have not yet been included in the initial heuristic formulas. For example, using the MFC during the tuning procedure, it was possible to detect and describe a characteristic of the diffraction solution, such as the semitransparency function.
Eng 2022, 3 39 The physical clarity of the formulas constructed using the MFC, taking into account their speed and accuracy, gives the method an advantage over alternative approaches: both over engineering formulas and over hybrid numerical methods, as well as over algorithms based on working with databases obtained for pre-calculated reference tasks.
Heuristic formulas are obtained on the basis of a numerical solution of the diffraction problem in a specific formulation, but the same technique can be extended to other problems. In addition, if the numerical results differ from the experiment due to the inaccurate formulation of the problem, then using the developed technique, appropriate corrections can be made so that the final heuristic formulas better correspond to the experiment. In this case, it is possible to identify physical characteristics that are absent in the numerical solution.
The unification of the MFC approach also lies in the fact that for different formulations of problems, for all types of scatterers and all types of geometry, we use the same technique for constructing heuristic solutions by adjusting the primary heuristic formulas. In this work, we construct a heuristic solution for a scatterer with boundary conditions of a specific type (a half-plane with impedance boundary conditions, such as a thin layer), but this technique for studying the diffraction process on two-dimensional structures based on a numerical solution can be extended to edges with a different profile and other types of boundary conditions.
Unlike known heuristic approaches, the MFC formulas can have a predetermined accuracy. Unlike rigorous analytical solutions, the MFC formulas have a uniform structure, can be guaranteed to be created in a predetermined time frame and are well subject to algorithmization. These qualities help improve the efficiency of application packages.
The sequence of steps in determining the scattered field is as follows: 1. conducting of ray tracing; 2.
recognition and classification of the type of obstacle in the path of the beam; 3. account according to the selected formula [47].
For the development of this approach, numerical solutions of two-dimensional and three-dimensional problems of diffraction on standard semi-infinite scatterers are required. Unfortunately, such solutions are not yet available in most well-known commercial packages.
The MFC approach can be used to calculate diffraction for large 3D objects, for which it is impossible to obtain a numerical solution. The procedure for adjusting the heuristic solution is carried out on a smaller scatterer. Having achieved the required accuracy, one can further apply the heuristic formula for a larger scatterer, which can no longer be calculated by rigorous numerical methods. In this case, the performance of heuristic formulas will not change, and the accuracy will even increase compared with the smaller size scatterer.

Conclusions
Heuristic formulas, constructed on the basis of the newly developed MFC, provide new opportunities for solving practical problems.
On the one hand, new opportunities allow an increase in the performance of solvers due to the accuracy and compactness of heuristic formulas constructed using the MFC.
On the other hand, new opportunities allow one to carry out a physical interpretation of numerical solutions with the aim of a deeper study of the physical features of the diffraction process. This becomes possible due to the decomposition of heuristic formulas into physically clear basic components.
Funding: This research received no external funding.