Abstract

Two curved targets are used to explore far-field superconvergence effects arising in numerical solutions of the electric-field and magnetic-field integral equations. Three different orders of basis and testing functions are used to discretize these equations, and three different types of target models (flat facets, quadratic-curved facets, and cubic-curved facets) are employed. Ideal far-field convergence rates are only observed when the model curvature is one degree higher than the basis order.

1. Introduction

Electromagnetic analysis has often exploited “variational” or “stationary” techniques in an attempt to obtain results that exhibit smaller error or faster convergence than might otherwise be obtained for a given order of approximation. Representations in terms of degree polynomials are expected to produce error rates of as , where is the characteristic cell dimension. For instance, a linear representation should produce convergence. Superconvergence occurs when the error decreases at a faster rate. Galerkin method-of-moments (MoM) procedures, which involve the use of identical functions for expanding the unknown surface currents and enforcing the integral equation, sometimes exhibit superconvergence in the far-zone fields or scattering cross section (SCS) [14]. Theoretical SCS convergence rates were established by Warnick for the electric-field integral equation (EFIE) and magnetic-field integral equation (MFIE) [5] and are reported below. One example of superconvergence arises with the popular Rao-Wilton-Glisson (RWG) solutions of the EFIE [6], which usually yield SCS error rates of O-O in practice instead of the rate expected from basis functions that are only complete to degree 0.

When and where superconvergence occurs are not always understood. It was demonstrated in [7] that mixed-order divergence-conforming basis functions often produce superconvergent far fields with the EFIE, but polynomial-complete basis functions apparently do not. In a complementary manner, polynomial-complete basis functions seem to produce superconvergent far fields with the MFIE, while mixed-order bases do not. Thus, the MFIE requires types of basis functions different than the EFIE for optimal SCS errors. Furthermore, ideal convergence rates should not be expected if target models are not accurate enough, for instance, if flat facets are used to model curved surfaces. Convergence rates are usually slower if the problem geometry contains edges or corners where current or charge densities may be singular, unless special basis functions containing the proper behavior are employed [8]. The purpose of the following is to attempt to clarify conditions under which one can reasonably expect to achieve the ideal convergence rates of [5]. Primarily, we investigate the issue of curved-cell models and the required accuracy in the testing integrals associated with the MoM procedure for perfectly conducting targets that are spheres or toroids.

We use mixed-order divergence-conforming vector basis functions of the Nedelec variety [9] on curved triangular patches to discretize the EFIE. The RWG basis functions [6] are the lowest-order () members of this family. Hierarchical versions of these functions are described in [10] and were implemented in a Galerkin scheme. Functions of orders , , and are considered and are reported as follows.

Local Definition of the Hierarchical Divergence-Conforming Vector Basis Functions Used with the EFIE, in Simplex Coordinates on a Standard Triangle. The Edge-Based Functions Straddle Adjacent CellsEdge-based base vectors () are,,.Order includes three additional edge-based functions,,and 2 cell-based functions,.Order includes three additional edge-based functions,,and 4 cell-based functions,,,.Inner integrals were evaluated using adaptive quadrature routines in conjunction with singularity cancellation procedures to return results accurate to some number of decimal places (6 digits for the results). Testing integrals were evaluated with a 7-point Gauss rule (that can exactly integrate degree 5 functions) for and and a 12-point Gauss rule that can exactly integrate degree 7 polynomials for . These rules were selected after numerical experimentation showed that higher-order rules do not improve the accuracy or convergence rates of the results. From Warnick’s analysis, the error in the SCS obtained from the EFIE under optimal conditions should converge at a rate of , where half-integer indices are used for to denote the mixed-order functions. For instance, the SCS error associated with functions is ideally expected to converge as , while that of the functions is ideally expected to converge as . These are superconvergent rates since they are different from the rates associated with the underlying current density, namely, . It is noteworthy that rates are observed for SCS values when point-testing is used with the EFIE instead of Galerkin testing [11].

On the other hand, [7, 12, 13] suggest that the MFIE SCS is not superconvergent unless polynomial-complete bases are employed. So for the MFIE we use polynomial-complete bases of orders (constant), (linear), and (quadratic) on curved triangular patches to treat the MFIE. The and functions are of the curl-conforming variety and belong to Nedelec’s polynomial-complete spaces [14]. The specific basis functions used are interpolatory and are reported as follows.

Local Definition of the Interpolatory Curl-Conforming Vector Basis Functions Used with the MFIE, in Simplex Coordinates on a Standard Triangle. The Edge-Based Functions Straddle Adjacent CellsOrder includes two constant cell-based functions:,.Order includes six linear edge-based functions:,,,,,.Order includes twelve quadratic functions, 9 edge-based,,,,,,,,and 3 cell-based,,.As with the EFIE, adaptive quadrature and singularity cancellation are used for the inner integrals. Galerkin testing is implemented with a 7-point Gauss rule, which was demonstrated to be adequate for all three orders. Under ideal conditions, the error in the SCS obtained from the MFIE should converge at a rate of , where integer indices are used for [5]. Thus the SCS error associated with functions is expected to converge as , while that of the functions is expected to converge as . Note that these rates are one order less than the corresponding rates for the EFIE SCS, based on the extent to which the basis functions are complete to a given degree, but are still superconvergent.

In the following, we investigate the convergence rates that can be obtained for a sphere and toroid target, as a function of the type of curvature used with the model (flat patches, quadratic-curved patches, and cubic-curved patches). This work is intended to build upon and extend the previous work of the author in [7, 11, 15] in establishing baseline convergence rates for numerical solutions of the EFIE and MFIE.

2. Results for Scattering from a Sphere

Triangular-cell models are obtained by dividing a sphere uniformly in between its north pole and its south pole and dividing it uniformly in so that the number of cell edges along the equator is twice the number of divisions in . Each octant of the sphere surface is therefore a curved triangle that is divided into smaller triangles of fairly uniform size. Figure 1 shows one example employing 72 flat patches. For curved patches, the cells are mapped from triangles to patches defined by the usual 6 Lagrange points (quadratic) or 10 Lagrange points (cubic), with all the points on the sphere surface. The mapping procedures for the divergence-conforming (EFIE) and curl-conforming (MFIE) bases are described in detail in [16]. The surface area of each subsectional-cell model is normalized to that of the desired sphere.

Relative accuracy is determined by evaluating the error in the bistatic SCS at 5° increments in and over the entire spherical range, with the 2-norm error computed according towhere denotes the error at point , , and where denotes the value of at point . The reference is the exact eigenfunction “Mie series” expansion [1].

Consider a sphere with or a surface area of approximately 12.57. Thus a density of 100 unknowns per square wavelength of surface area is reached when the total number of unknowns equals 1257. Figure 2 shows the error in SCS versus number of unknowns for results obtained from the MFIE. Results are reported for three representations, a flat-faceted model used in conjunction with basis and test functions (piecewise constant), a quadratic-curved model used with functions (linear), and a cubic-curved model used with functions. As expected, the higher-order bases and higher-order target models produce more accurate solutions. Figure 3 shows the same data plotted versus the average edge length used in the models. It is apparent that the result closely follows an trajectory, the result follows an trajectory, and the result follows an trajectory, for small . These are the ideal rates suggested in [5]. Also shown in Figure 3 is a plot of the corresponding error for the case when quadratic-curved patches are used in the model. With the quadratic-patch models, the results follow a slower rate, near . Similar behavior was observed by the author for other targets (not shown), suggesting the conclusion that bases must be used with cubic-patch models to produce optimal convergence rates.

The SCS errors produced using the EFIE are presented in Figures 4 and 5 for the same sphere target. Data were obtained using (RWG), , and mixed-order basis and test functions. Figure 4 shows results versus the number of unknowns, with all models represented by cubic-curved patches. In common with the MFIE results, once the number of unknowns is in a reasonable range for a problem of this electrical size, higher-order bases produce more accurate results for a given number of unknowns. Figure 5 shows EFIE results plotted versus the average edge length used in the models, but for the case when the results are obtained with flat-faceted models and the results are obtained with quadratic-curved models. The results follow an trajectory, while the results in this case are somewhere between and . The results appear to follow a trajectory between and . Thus, these data are close to the ideal rates of [5], but as observed in [7, 15], the ideal rates are not always obtained in practice. The author has also investigated the use of bases with quadratic-patch sphere models (not shown) but always observed a substantially slower convergence rate for that combination. This suggests that bases must be used with cubic-patch models to approach optimal convergence rates.

3. Results for Scattering from a Torus

A torus is a simple target that, like a sphere, avoids the complications introduced by edges and corners. For the following comparisons, a torus with major radius of λ and minor radius of λ, for a total surface area of 18.0, is employed. The torus is centered at the origin, parallel to the - plane, and is illuminated by a uniform plane wave in the -direction with the electric field polarized along . For illustration, Figure 6 shows a torus model containing 196 flat triangular patches.

The torus is divided uniformly along both circumferences into curved quadrilateral cells, which are each divided into two (curved) triangles. The models employed twice as many divisions along the larger circumference as along the smaller and are scaled to the proper surface area. For comparison purposes, a reference solution was obtained using the formulation described in [17], with basis functions complete to polynomial degree 5, and a total of 9216 unknowns. Relative accuracy is determined from (1) by evaluating the error in the bistatic SCS at 5° increments in and 90° increments in . In the figures for the torus example, a density of 100 unknowns per square wavelength of surface area is reached when the total number of unknowns equals 1800.

Figure 7 shows the SCS error produced by the MFIE versus the number of unknowns for flat-faceted models with functions, quadratic-curved models with functions, and cubic-curved models with functions. Figure 8 shows the same data plotted versus the average edge length used in the models. For these results, the data exhibit a convergence rate between and , while the and data closely approximate the ideal trajectories of and .

Figures 9 and 10 show EFIE results for the SCS error for the same torus, plotted versus the average edge length. The and results converge somewhat slower than the ideal rates of and , while the data appear to closely follow the rate. Figure 10 shows that the convergence rate for the results is slower than if quadratic-patch models are used, while it is between and when cubic-patch models are used.

4. Conclusions

Far-field convergence rates are studied using two curved targets. Numerical solutions of the EFIE and MFIE are employed with higher-order vector basis functions and target models that involve flat facets, quadratic-curved patches, or cubic-curved patches. Superconvergence in the SCS at rates close to those suggested in [5] is observed in the results, provided that the models in use have high enough curvature, and the testing integration is done with quadrature rules of sufficient accuracy.

An interesting result of these observations is that it appears necessary to employ cubic-curved-patch models in order to obtain ideal convergence rates with quadratic basis functions and the MFIE. This is at odds with the classical literature on numerical methods where it is supposed that quadratic models should be adequate with quadratic expansions. It also appears necessary to employ cubic models with the mixed-order basis functions (which are only complete to degree 2) used with the EFIE to approach ideal SCS convergence rates. A similar conclusion applies to the lower-order bases; namely, the linear functions (used with the MFIE) and functions (used with the EFIE) require quadratic-patch models to yield near-ideal convergence rates.

Competing Interests

The author declares that he has no competing interests.