Abstract

The electromagnetic vector fields, which are scattered off a highly conductive spheroid that is embedded within an otherwise lossless medium, are investigated in this contribution. A time-harmonic magnetic dipolar source, located nearby and operating at low frequencies, serves as the excitation primary field, being arbitrarily orientated in the three-dimensional space. The main idea is to obtain an analytical solution of this scattering problem, using the appropriate system of spheroidal coordinates, such that a possibly fast numerical estimation of the scattered fields could be useful for real data inversion. To this end, incident and scattered as well as total fields are written in a rigorous low-frequency manner in terms of positive integral powers of the real-valued wave number of the exterior environment. Then, the Maxwell-type problem is converted to interconnected Laplace’s or Poisson’s equations, complemented by the perfectly conducting boundary conditions on the spheroidal object and the necessary radiation behavior at infinity. The static approximation and the three first dynamic contributors are sufficient for the present study, while terms of higher orders are neglected at the low-frequency regime. Henceforth, the 3D scattering boundary value problems are solved incrementally, whereas the determination of the unknown constant coefficients leads either to concrete expressions or to infinite linear algebraic systems, which can be readily solved by implementing standard cut-off techniques. The nonaxisymmetric scattered magnetic and electric fields follow and they are obtained in an analytical compact fashion via infinite series expansions in spheroidal eigenfunctions. In order to demonstrate the efficiency of our analytical approach, the results are degenerated so as to recover the spherical case, which validates this approach.

1. Introduction

The fundamental principles of classical electromagnetism [1], initially introduced by the pioneer in the field James Clerk Maxwell, represent the basic foundation of the low-frequency scattering theory [2]. Many physical applications that are associated with the response of receivers with arbitrary shape, embedded within various media, when they are excited by several primary sources, still stand in the frontline of scientific research. Indeed, let us give emphasis to some of them like two-phase composites, Earth’s subsurface probing for mineral exploration, identification of cavities or other underground detection such as UneXploded Ordnance and buried objects, and scattering by chiral material in either chiral or nonchiral environments. Obviously, the range of the applied frequencies varies according to each case. However, low-frequency interactions exhibit a particular interest for numerous real-life problems that occur in physical areas, which are related to the present contribution. By deciphering the implicated fields to each case, information about main parameters like orientations, sizes, shapes, and magnetic and electric properties of the anomalies brings insight into the field behavior. Nevertheless, this is not at all an easy task, since the situation of an inversion scheme [3] cannot be tackled in a robust-like manner, unless efficient models of the field distribution and strong effective mathematical tools [4] are available, since in practice the measurement and the identification of the scatterer are needed simultaneously.

The large amount of vector data, the electromagnetic and geometrical complexity of the different media involved, the many configurations of sources and receivers, and the uncertainty resulting from datasets containing both the contribution of the incident and scattered fields explain the interest of solving forward and inverse electromagnetic scattering problems, first at the modeling stage and later at the inversion stage. Therefore, the already ample library of scattering by simple shapes using analytical methods is open to accept new and useful semianalytical and, why not, analytical results. In the last decade, a lot of articles were published in that direction and many references are available in the literature. Collectively, we pertain to the low-frequency scattering by a perfectly electrically conducting sphere in a conductive environment, illuminated by an arbitrarily orientated and settled magnetic dipole source field [5]. There, the 3D scattered electromagnetic fields have assumed a complete analytical formalism in terms of closed-type solutions, while a numerical implementation has validated almost perfectly the results via a comparison with the exact Mie series expansion [6]. However, this work was a first introduction to scattering by metallic objects that reflect the complete isotropy of the three-dimensional space. Moving to more complicated geometries for the same physical consideration such as spheroids [7] or ellipsoids [8], one can realize that the orientation of the body, save the position or the size, plays a crucial role not only to the mathematical treatment but also to the computational elaboration of the corresponding fields. Moreover, it is obvious that the more general the shape becomes, the more difficulty is added to a possible inversion, which, by the way, is the final goal when we deal with real-life applications. It is to this end that the utilization of elaborate numerical techniques is inevitable, as it is demonstrated in [9] for the identification of an orebody of approximately infinite conductivity, which is buried in Earth, by using an equivalent perfectly conducting triaxial model ellipsoid. High-contrast cases, where the target has an extremely large conductivity with respect to the surrounding host environment and hence is considered to be impenetrable, motivated researchers to embark upon much more complex patterns, such as the two-sphere case for two closely adjacent metallic objects [10] or the circumstance of a nonpenetrable, that is, perfect conductor, ring torus [11], both articles presuming a conductive occupied medium. These reports, each one containing a sufficient reference list, are a small but representative sample of displaying the inconvenience induced in performing analytical techniques when incorporating with different geometrical models [12, 13]. As a matter of fact, the difficulty increases due to the emergence of intricate eigenfunctions [14, 15], dealing with the involved potentials.

This project is focused now on the adjustment of the low-frequency diffusive scattering theory in solving the particular problem of identifying impenetrable metallic bodies in a lossless, that is, perfect dielectric, medium, aiming to follow the similar successful path of other contributions to the field. Bearing in mind the few, but with quite good results, attempts of simulating the environment as lossless and not conductive (e.g., see [16] for the case of two metallic spheres almost touching and [17] for the toroidal case in this direction), we realize the necessity of dealing with such brute-force analytical models. Seeking for different body structures in order to attain the representation of all kind of anomalies in lossless media, as much effectively as possible, the present paper deals with the spheroidal configuration. Although research on quasi-static fields around spheroids is vast, there are still elementary properties to be discovered; thus the demand for comparable investigations stays on continuously even nowadays. Consequently, here we extend the research to a perfectly reflective spheroid, embedded in a lossless ambient, where the appropriate and best fitting system for the particular modeling purposes is the spheroidal coordinate system [12, 13]. The prolate spheroidal geometry is considered, since the oblate one is readily obtained via a simple analytical transformation, while limiting cases are also available. On the other hand, a time-harmonic magnetic dipole source, acting as the primary source, operates at the low-frequency regime and generates the three-dimensional incident electromagnetic waves, which propagate towards the scatterer. Then, the spheroidal-shaped body responds to the excitation and generates the scattered magnetic and electric fields, where their summation comprises the total fields.

The physical nature of the problem urges us to cope with the following low-frequency technique, where the incident, the scattered, and the total electromagnetic fields are expanded in terms of positive integral powers of , where denotes the real wave number of the exterior medium at the operation frequency. Here, we elaborate the 3D vector fields at each order of for every , whereas our problem is transformed into a sequence of coupled boundary value scattering problems. The model is mathematically formulated in terms of second-order Laplace’s and Poisson’s partial differential equations, accompanied by the appropriate perfectly electrically conducting boundary conditions, that is, cancellation of normal magnetic and tangential electric fields, reflecting the impenetrable character of the boundary of the spheroid-type metallic targets. The Silver–Müller radiation conditions at infinity are applied as well. We restrict our analysis to the important terms of the expansions, confined by the static Rayleigh term for and the dynamic terms for , since terms of higher orders, that is, , are small contributors due to the low frequencies of which the source operates and, consequently, they are neglected. The analytical procedure sketched above leads either to closed relationships or to infinite linear algebraic systems for the determination of the unknown constant coefficients of the implied potential functions, where the latter can be solved approximately by enforcing usual cut-off methods. The three-dimensional scattered magnetic and electric fields for each admit closed-type forms via infinite series expansions in view of spheroidal eigenfunctions [14, 15]. The key to the validation of our approach is the degeneration of the obtained formulae in order to recover the spherical case.

The rest of the paper itself is organized as follows. In Section 2, a detailed physical interpretation and the theoretical basis via an analytical mathematical formulation are sketched, which incorporate the general features, independently of any geometry. Thereafter, a brief configuration of the prolate spheroidal system and the eigenfunctions associated with it is developed in Section 3, whereas the sought solution is constructed step by step and the electromagnetic low-frequency fields assume infinite series expansions in terms of prolate spheroidal eigenfunctions. Though the transformation of these geometrical and analytical attributes to the oblate spheroidal case is presented, some interesting limiting occasions are also discussed. Section 4 is devoted to the validation of the aforementioned analytical method via a proper reduction to the corresponding spherical problem that illustrates the above results. Finally, a short conclusion and discussion, containing an outline of our work and future steps, follow in Section 5, while an updated reference list is displayed immediately after.

2. Physical and Mathematical Interpretation

Suppose a solid spheroidal body, either prolate or oblate, which by physical means has an approximate conductivity with respect to the surrounding homogeneous, isotropic, and nonmagnetic lossless medium of conductivity . The constructed body-medium system is assumed to occupy the three-dimensional space , while the interface that distinguishes the two regions is actually the smooth surface of the perfectly electrically conducting object. Hence, the area of electromagnetic scattering is spread to the unbounded domain of dielectric permittivity and magnetic permeability , approximated by the corresponding one of the vacuum, outside the spheroid, yielding exterior-type problems. By definition of the given low circular frequency , the real-valued wave number of the medium iswhere (1) secures the conjecture of low frequencies. Harmonic time-dependence is implied for the fields via multiplication of their spatial part by (imaginary unit ); consequently they admit expressions as a function of , which is expressed via the Cartesian basis , in Cartesian coordinates , where this dependence will be omitted henceforth for reasons of writing convenience.

The impenetrable spheroid is excited by a primary vector source field, for example, a magnetic dipole for the purposes of this project, having an arbitrary orientation that yieldsbeing located at . The fixed singular point is excluded from the domain of scattering, since the region of interest is at a reasonably long distance from the source. Thus, we adopt the explicit definition as for the area of electromagnetic activity. According to the simple notation and , the magnetic dipolar source radiates the primary fields [2]where the symbol “” denotes dyadic product. In terms of Maclaurin’s series expansions of the exponent and implying (1), we perform some extended algebraic calculations by collecting terms as integral powers of and we retain the adequate number of them that correspond to the static Rayleigh approximation and to the first three dynamic factors. Therefore, the incident magnetic field (3) is rewritten aswhile the corresponding primary incident electric field (4) renderswhere the input of the fields of order is minor, since and consequently are very small. On the other hand, the surviving electromagnetic fields for all enjoyas far as the magnetic components are concerned, whileobtained for the electric elements, where , as it is indicated by the performed analysis. Expressions (7)–(11) are provided in view of the unit dyadic , while the handy forms, depicted on the right-hand side of (7), (8), and (10), are based on straightforward computations with the aid of trivial differential identities and by virtue ofgiven that . They contain the fundamental quantity, captured under the action of the gradient operator (note that the Laplacian is ). Though these operators could also be applied on , however for reasons of technical clarity we distinguish this feature by introducing the above notation; otherwise it will be instantly designated. In the aftermath, the primary incident fields are perturbed by the solid object, producing the corresponding scattered fields and , while the total fields emerge as the summationwhere, since the spheroidal metal body is nonpenetrable, the wave propagation is not possible inside it. By definition of the simple notation , corresponding to primary (), scattered (), and total () electromagnetic fields, all pairs satisfy the well-known Maxwell’s equations [2]. Yet, as for the scattering components that interest us in the sought solution within , they becomewhich represent a deduced form for our particular circumstance and they are divergence-free; that is,implying their solenoidal attribute, such that (14) is readily sufficed.

The low-frequency interaction high-contrast cases, as dictated by the physics that we address here, are inherited to the scattered fields and they are postulated by the serieswhere their written fashion is referenced in [2]. Substituting expansions (16) into Maxwell’s equations (14), accompanied by (15), we retrievewhich stand for the low-frequency equivalent expressions. However, for all the reasons described earlier, the behavior of the incident fields (5) and (6) with (7)–(11) reflects exactly the very same physical and mathematical treatment to the scattered fields (16) with (17) and (18); that is, we hold the corresponding nonvanishing terms, attaining fields up to the third order, since the rest of them consist of a minor correction. As a consequence, the scattered magnetic field readswhile scattered electric field is assumed to beconsidering , in the absence of the relative incident components. Obviously, taking this into account, the fundamental relations (17) and (18) are immediately fulfilled by the surviving low-frequency fields , , , , and arising from (19) and (20) for , where their evaluation specifies our final goal. Aiming to understand the contribution of each low-frequency term within (19) and (20), we reinforce (1) and we separate the real from the imaginary part, concluding tothe electric field being purely imaginary, while the magnetic field that is practically measured is a complex-valued function.

In order to proceed to the solution, we are obliged to complete the aforementioned analysis by the exposition of the related to the problem boundary and limiting conditions. Under this aspect, the properly fitted perfectly electrically conducting boundary conditions on the surface of the spheroidal target are provided in terms of the corresponding outward unit normal vector and involve the total fields viawhich demand cancellation of the normal component of the total magnetic field and of the tangential component of the total electric field. The low-frequency analogy of (23) turns out to be a combination of it with (16), whereas by virtue of (13) we are led towhile relations (24) refer only to the remaining electromagnetic terms, since for any . Additionally, for the implicated functions, the Silver–Müller radiation conditions at infinity for the scattered fields, those beingmust be instantly satisfied, also for the low-frequency regime through (16), givingsince, for , it holds that and as it is verified from (17). Exterior-type solutions, like in our case, satisfy (25) automatically, as a consequence of the appropriate construction of the potentials via the corresponding eigenfunctions [14, 15].

During our analysis, it is supposed that the conductivity of the spheroid approaches infinity. Nevertheless, herein we should include a small discussion about the transition zone, that is, how to treat materials with large, or even comparable with respect to the exterior medium, but not infinite conductivity. In this case, the problem becomes slightly more complicated, since the body becomes penetrable; hence the interior electromagnetic fields exist and contribute to the final solution via the transmission conditions, which are applied on the material’s surface instead of (23); see [2] for more details. There, it is shown that conductive or even dielectric particles can resonate with an incident electromagnetic wave over a wide range of modes related to surface excitation. This is due to the excess surface charge that causes the particles of arbitrary composition obtaining metallic properties, a fact especially important for spheres and spheroids. However, this is a completely different situation from the present investigation, concerning applications of importance in other physical areas, for example, in optics.

In summary, we have to manipulate a little bit further the partial differential equations (17) with (18) to obtain easy-to-handle elliptic equations with partial derivatives and, then, use the spheroidal coordinate geometry and analysis, whereby the conditions (24) and (26) are applied, to reach the final desired solution.

3. Spheroidal Low-Frequency Electromagnetic Fields

We intend to derive closed analytical formulae for the main nonvanishing scattered electromagnetic fields , , , , and , since . To achieve that, we have to recover from (17), provided (18), known partial differential equations. To this direction, we apply on both sides of Maxwell’s equations (17) and we utilize repeatedly the very useful identity with being any vector, where, along with (18), we arrive at a set of boundary value problems, which are coupled to one another, from the static one at to all the rest of the dynamic ones at higher orders up to . Those are written explicitly as functions of the scalar , and the vector , harmonic potentials,providing independentlywhile utilizing (28); then the right-hand side of what follows is known; that is,from which we may recoverwhich serves as the nonhomogeneous part of the equation that follows; that is,Moreover, we work similarly to the zeroth-order field to obtainwhere we recall that the scattered fields , , , , and must be determined in the confined scattering region of wave propagation . Some extended analysis reveals that, for , which corresponds to the Rayleigh approximation term, an uncoupled Laplace’s equation is solved for the magnetic field, while there exists no electric field, since the corresponding term for the incident electric field is absent. On the contrary, the dynamic term for reveals the existence of a vector harmonic electric field with a zero-term magnetic field. As for the case, it involves a Poisson’s equation for the magnetic field, coupled to the static term. Exploiting appropriately the problem, we solve Poisson’s equation by introducing a special form of a particular solution, which is actually the outcome of the application of a standard vector identity, using the fact that and . The incident electric field does not exist, inheriting the same attribute to the scattered one. Finally, at both magnetic and electric fields appear. As far as the magnetic field is concerned, we simply solve an uncoupled Laplace’s equation. However, this is not the case for the problem corresponding to the electric field for which the solution is introduced in the form of an integral representation in terms of the fundamental solution of Laplace’s equation [14]. Problems (28)–(32) with (27) are supplemented by the perfectly reflecting boundary conditions (24) and the proper limiting behavior (26). It is obvious that the low-frequency magnetic terms of order vary like , while the corresponding electric ones vary like as increases to infinity. Thus, implying to obtain , integral within (31) is of order , which converges as .

3.1. The Prolate Spheroidal Coordinate System with Geometry and Analysis

For the sake of completeness, we provide a brief analysis of the spheroidal case and especially the prolate spheroidal geometry, since the oblate one can be readily recovered as it will be seen afterwards. Given a fixed positive number , which we consider to be the semifocal distance of our system, we define the transformed prolate spheroidal coordinates for every , , and (see Figure 1) aswith position vector , where = (, , ) designate the prescribed source’s location. In the limit as the prolate spheroid degenerates to the standard focal segment , while as it approaches a sphere located at infinity. A key point in this analysis is the exploitation of the common -axis of symmetry between the applicable coordinate system and the spheroidal body under investigation. In fact, this is a quite practical technique, which simplifies the mathematical treatment, though it retains the generality unaltered and we still cope with a genuine three-dimensional problem, since the arbitrary orientation of the magnetic dipole, which is implied by (2), defends the 3D configuration until the final solution. Nevertheless, in the very special case where the magnetic dipole is set on the axis of symmetry of the proposed prolate geometry, that is, the -axis, which means and is orientated to the -direction, yielding , then the problem under consideration turns into being axisymmetric, since any kind of dependence from the azimuthal angular is omitted. This restriction is imposed in every field in addition to the vanishing of the other two components of the magnetic dipole () and of its position () to the analytical results. However, the general nonaxisymmetric case that is examined here is concerned with a source of arbitrary location and orientation.

The main differential operators, in the prolate geometry, excluding the tip for and for every , , and , are furnished byas for the gradient andconcerning the Laplacian. In (34), the outward unit normal vector coincides with throughwhereas the other two orthonormal vectors of the spheroidal system assume the Cartesian basis formulaewhile the unit dyadic in prolate spheroidal coordinates is .

Specifying the solid surface by , then the prolate spheroidal body has major (in -direction) and minor (situated on -plane) axes and , respectively, while its eccentricity is . We remark that the low-frequency approximations are secured only when the main axes, as well as the semifocal distance, are less than the wavelength of the incident field, such a demand being ensured via the inequality . In terms of spheroidal modes, the actual region of wave propagation iswhere corresponds to by means of (33), forming an exterior domain. Apparently, since the spheroidal object is impenetrable, it lacks any field inside; that is,while on and as , the boundary conditions (23) or (24) and the limiting behavior (25) or (26), respectively, must be counted.

Proceeding to the analysis, we introduce the interior (regular as ) and the exterior (regular as ) prolate spheroidal harmonic eigenfunctions of degree and of order , in terms of the associated Legendre functions of the first and of the second kind (for more information about these functions and their properties, refer to [7, 14, 15], where a detailed appendix contains such details) according torespectively, and for any , , and , which both are regular on the axis of symmetry. The azimuthal angular dependence readswhere stands for the even () or the odd () part of the eigenfunctions, while the prime denotes derivation with respect to the argument here and thereafter. Consequently, every harmonic function in prolate spheroidal geometry belongs to the kernel space of Laplace’s operator (35), that is, , and assumeswhere and are arbitrary constants, which selectively vanish once we deal with either an interior problem () or an exterior one (), like in our case, in order to obtain the mathematically accepted solution. On the other hand, we join the internal productneeded in our project, being the Kronecker delta and , while the trigonometric functions for and are orthogonal with respect towhere and are again the Kronecker deltas with and otherwise.

Furthermore, a very crucial piece of information for the forthcoming analytical subsections incorporates the position of the dipole at with the very useful Green’s function expansion for prolate geometry [14], so that in terms of (41) and for any observation region (39), we obtainwhere at we attain singularity and it holdswhich ends the well-known, but necessary, gathering of preliminaries for the sequel.

Finally, for completeness, we invoke the exact expressions of gradient action upon interior and exterior solid spheroidal harmonics (41) via the surface ones for the shake of clarity. This is feasible, since and for , , and belong to the subspace produced by the surface spheroidal harmonics. After cumbersome analysis, profiting by recurrence relations for and , defined for different arguments, that is, for and [14, 15], we arrive atas far as the gradient on the interior solid harmonics is concerned, whileregarding the gradient on the exterior solid harmonics. Both (48) and (49) are defined for any , , and , while only when ; else it is zero. It is obvious that expressions (48) and (49) are also applicable at , if we make the replacement , whereas then the prime declares derivation with respect to .

Next, we perform long and tedious calculations in order to pursuit an analytical form for the scattered fields. In what follows, we present the basic and interesting steps of them by starting from the magnetic modes to obtain , , and and proceed to the electric ones , , bearing in mind the vanishing electromagnetic terms , where and are the most cumbersome cases, which we endeavor to retrieve in a handy form.

3.2. The , , and Magnetic Fields

Beginning with the situation for the Rayleigh static term, the relative potential boundary value problem that has to be solved embodies (27) and (28) with the Neumann boundary condition (24) for , yieldingwith , where is provided by (36) for . The harmonic potential is rendered via the exterior-type part of (43); hence, by virtue of (50), we conclude towhere for , , and denote the constant coefficients to be recovered from boundary condition (50). Expansion (51) is compatible with the acquired limiting behavior (26). In the first place, we imply definition (34) upon (51) and we calculate the inner productfor every and . On the other hand, the corresponding inner product with the incident field (7) encounters a complexity, caused by the double action of the gradient operator over Green’s function quantity. Aiming to exceed this difficultness, we work as follows. Firstly, we take the unit normal projection of (34), where the outcome issince the dyadic is symmetric, and we employ (12) to rewrite (53) aswhich, upon introduction of eigenexpansion (46) with (47) for , becomeswhere that has to do with the source is defined in (47) and , while and . Apparently, the gradient is a known quantity and satisfies (2). Combining (52) and (55), condition (50) reveals thatat any and . Orthogonality arguments based on (44) and (45) lead to the obvious and compact resultat the singularity point . Substituting relationship (57) into (51), we reach the sought solution, that is,for every , whereas and are provided via the right-hand side formulae of (41), while arises from (49), the same thing for , by a simple switch of arguments, all written for , , and .

In the sequel, we treat the magnetic field-case for , which exhibits a particular difficulty, mainly due to the coupling with the static term, indicated by (29). However, secondarily speaking, since is a pure vector field, needing a set of three separate boundary conditions to be determined, then, beyond the application of the scalar condition in the first part of (24) at , it is indirectly associated with the total electric field on the boundary through the second part of (24) by implying (10) and (30). All the above are mathematically summarized as follows:where, provided directly from (8) and given (46) with (47),which completes the boundary value problem that has to be solved for this situation, by the imposition of the divergence-free property of that must be always valid viaInserting, now, (57) into (51) and combining it with the fact that , then the zeroth-order potential is known and admitswhereas for , , and follow from (47). On the other hand, the vector harmonic function keeps the exterior-type solution fashion, leading towhere (63) meets the radiation criteria imposed by (26), whiledenote the vector character of the unknown constant coefficients , while, willing to be consistent with (64) and work within the frame of purely prolate spheroidal geometry, we are obliged to write the spheroidal representation of the position vector, entering (59) asThe field problem under consideration has an additional physical difficulty, which is inherited from the corresponding mathematical complexity. Thus, in what follows, we give emphasis on the main steps due to the large number of calculations to perform. To this end, we handle (60) in such a manner so as to obtain three separate sets of boundary conditions, one from the -component and two from relations, which correspond to the , -components. This is accomplished by operating from (36) into (60) and tackling with (46) for . Once done, we reach cumbersomely the following three relations for the coefficients, satisfyingcounting for the aforementioned three sets of boundary relations, where , , and . Omitting dependence on and for reasons of writing clarity, the leading functions of the unknown constant coefficients and those comprising the constant terms of (66) are furnished initially for by the formulaewhile for it is readily obtainedFinally, the case is structured by the functionsAll notations, (67)–(69), are valid for any , , , and , while the unit normal vectors , , and are defined in (36)–(38), the ad hoc function comes from (47), the dipole source is met by (2), and the position vector is taken from (65), while by definitionfor , , and . Notice that derivations are with respect to the argument. Certainly, relationships (66) accompanied by definitions (67)–(69) are very complicated to accept elaboration with respect to recurrence relations, since in the end a numerical implementation for the resulting systems is inevitable. However, in order to achieve a compact and easy-to-handle analytical formation of these conditions, we expand the already known functions (67)–(69) in terms of the priori orthonormal basis for , , and viaand for every , , and , where, by virtue of the orthogonality properties (44) and (45), we end up with the easily amenable integral representationsfor every , , and with . The last step includes supersession of expressions (71), provided (72), into (66) and processing of a handy orthogonal argumentation based on (44) and (45), so as to conclude towhile , as well as , , and , offer all existing combinations regarding (73). In doing so, we achieved conveying the arisen difficulty from the boundary conditions (66) to the simple calculation of integrals (72), which just incorporate trigonometric and associated Legendre functions, so as to become a matter of performing trivial integrations or applying a standard computational code. Thereafter, (73) represents systems of linear algebraic equations, which can be solved via usual cut-off techniques by the imposition of an indispensable common upper limit for both the degree indexes; let us assume , in order to obtain quadrature forms, such aswhere for and , as well as , we havewhich denote the squared-type invertible matrix of the coefficients of the unknowns, the vector of the unknown coefficients, and the vector of the known constants, respectively. The vector specifies the system’s solution through (74), being able to evaluate the unknown constant coefficients for every value , , and up to a certain level, where the desired accuracy is achieved. Thus, by substitution of (62) and (63) into (59), we reach the analytical form of the scattered magnetic field for this occasion, that is,for , satisfying (61), wherein, upon solution (74), the prolate spheroidal modes, indicated by (64) and (65), can be readily embodied.

The simplest case of this project concerns the third-order scattered magnetic field at due to the constant character of the corresponding incident field (9). Similarly to the Rayleigh approximation, here we are obliged to solve the potential boundary value problem (32) with the Neumann boundary condition (24) for , that is,by means of , where is given through (36) for . In order to comply with the demand of well-behaved solutions at infinity via (26), we again choosewhereas for , , and stand for the constant coefficients to be determined from (77). Consequently, invoking the source (9) and projecting from (2), then, accordingly to (78), the convenient relation (49) is applied to (77) in order to receivefor every and . Since within (79) we can write , , and , obviously standard orthogonality arguments in terms of (44) and (45) yieldwhile the remainder ones readonce we enforce the useful identityprime denoting derivation with respect to . Inserting (80) and (81) into (78), we get which introduces a compact dyadic formation, bearing in mind thatare conveniently chosen functions for every and .

3.3. The and Electric Fields

The first nonvanishing electric term appears as at the low-frequency regime and it is interconnected with via (30). Hence, we incorporate the latter with the already known field (76) and with aid of a simple differential identity we obtainsince and , where is the right-hand side expression of (41), while is taken from (49), the singular quantity follows from (47), and, finally, the constant coefficients are provided via the magnetic situation, all written in view of , , and .

Our final task is the calculation of the exact low-frequency problem at , where here the main difficulty is due to the integral implication with the electric field (85) of the offered solution; see (31) for instance. Henceforth, we deal with a boundary value problem, being the outcome of the combination of (31) and of the second part of (24) for , that is,for every and on with and , while we recall the corresponding incident electric field (11). Field (86) is solenoidal, which is secured viasince , whereas (12) can also be employed. Obviously is a vector field, retaining three components, which require three scalar boundary conditions. However, we notice that (86) offers two of them, forcing us to retrieve the third one by the application of the restriction condition (87) on the surface of the metallic prolate spheroid, meaningwhich completes the well-posed analysis. In order for it to be compatible with the radiation condition (26) for , we expand as usual, being an exterior potential,where, similarly to (64), we haveIn the prolate spheroidal geometry, we utilize the infinitesimal elementary volumeto rewrite the particular integral solution that appears into (86) for any asand , whereas we have implied a specific technique to bound the singular point, by introducing a very small positive number , while for reasons of simplicity and without loss of the generality, we assume that by definition. The first and the third integral of (92) are analytic, since they exclude the singular point; hence therein we can use collectively (46) to represent the different areas, by applying the lower branch for and the upper one for , respectively. This action conducts towhich are harmonic, where the leading constants for , , and ,respectively, are written by virtue of (47) at the relative integration interval on . On the other hand, this is not the issue with the second integral, which contains the singularity and implying (91), it provides us with the particular solution in (31) viain terms of the Dirac function in Laplace’s fundamental solution [14]. Yet, this term is treated differently, by expanding it with respect to the orthonormal basis ,for , where, taking into account (91) and the inner products (44) and (45), we havefor . Gathering all information in between (92) and (98), integral (92) is rewritten as for every , which combined with (89) yields the scattered field (86). Under the aim of expanding the incident field (11) in the same manner, we incorporate (46) and (47) for , which is in fact the surface of the body, whereas is activated and since , thenwhereat . We are ready to precede fitting of the above analysis to (86) and (88), which are the boundary conditions to be satisfied. From this point forward, we deal with the problem in a very similar way to that for obtaining , since it possess the same complexity. Hence, putting all the mathematical tools together, (86) and (88) with (89)–(101) merge towhere, for , the case, originated from reinforcement (88), constitutesBearing in mind, now, the two components of condition (77), then for we are led to the relationswith the case being set by the relationsExpressions (103)–(105) act for , , , and , where the vectors , , and are with (36)–(38), the dipole source is given in (2), and for , , and assume (48) and (49), respectively, derivations are with respect to the argument, and notations (94), (95), (98), and (101) have been implied. Those are expanded through the orthonormal basis for , , and so thatas for , , and , whereas projection in view of (44) and (45) reveals thatfor every , , and with . Including (106) and (116) with (107) and (108) within (102), we end up withwhile , as well , , and . In matrix form, (109) readswhich are linear algebraic systems, where, imposing a fictitious upper limit , then for any and , as well as , we havewhich stand for the squared-type invertible matrix of the coefficients of the unknowns, the vector of the unknown coefficients, and the vector of the known constants, respectively, of this particular problem. A standard cut-off method is applied to (110) in order to decide the truncation amplitude , which supplies with an acceptable accuracy of a possible numerical implementation. Afterwards, the coefficients are determined for , , and , while the scattered electric field, comprising this situation, is built up from (86), encountering (89) and (99), that is,being established by the prolate spheroidal configuration of in (90) for , , and , once notations (94), (95), and (98) are considered.

3.4. Recapitulation: The Oblate Spheroid and Limiting Cases

Summarizing all the above analysis, we have retrieved in a closed-type analytical fashion the surviving low-frequency terms of the scattered magnetic , , and electric , fields under consideration in the prolate spheroidal domain of interest from well-posed boundary value problems, while . Collecting the associated terms (58), (76), (83) and (85), (112), respectively, we may substitute them into (21) and (22), so as to obtain the electromagnetic fields in a compact analytical form as infinite power series expansions, being ready for further possible numerical implementation.

Recapitulating and for reasons of completeness, we mention that the corresponding results of the oblate spheroidal geometry are recovered through the simple transformation [12]where and are the new characteristic variables of the system. On the other hand, the asymptotic case of the needle of length can be reached by a prolate spheroid when , while, for the case where , the oblate spheroid takes the shape of a circular disk of radius .

4. Analytical and Numerical Validation

The spheroidal geometry yields a reliable and quite representative coordinate system when we have to deal with scattering wave problems, where the field is perturbed by a perfectly conducting body. However, this is true only when the low-frequency approximation is implied, since we are able to deal with boundary value problems of Laplace’s equation, which is fully separable in this geometry. Hence, problems similar to our case adopt this fitting system of coordinates to obtain analytical or semianalytical results for the corresponding fields.

Nevertheless, it is not always easy and feasible to pursue fully analytical solutions in closed forms without computational error, or the obtained closed formulae do not have the required accuracy, even in the low-frequency realm. Indeed, even though the scattered electromagnetic fields have been calculated for (higher order terms are not of substantial interest) in a closed analytical form of infinite series in terms of the spheroidal harmonic eigenexpansions, they are not given in fully compact fashion. This is due to the fact that the constant coefficients in some circumstances are evaluated up to a certain order of accuracy via cut-off techniques. There, accuracy is controlled by the sum of terms taken in each case, which inherits the semianalytical character to our method. To this end, we need validation methods, both analytical and numerical, in order to cross-check the credibility of our results, those being the analytical reduction to the sphere-body case and the implementation of the obtained formulae.

4.1. Reduction to the Spherical Case

The prolate spheroidal geometry degenerates to the spherical one [12] when . For the corresponding analytical reduction, the limiting process involves an appropriate combination of the semifocal distance with the coordinate variable , such aswhere is the radial component of the spherical coordinate systemwhich is readily recovered from (33), pointing out that the other two radial variables remain unaltered, because they coincide in both systems. Obviously, if we multiply nominator and dominator of (36) and (37) by and use (114), we obtain the spherical radial unit normal vectorwhere and , while the remaining vectors assume similar spherical form.

The corresponding reduction of formulae (58), (76), (83) and (85), (112) for the electromagnetic fields is mainly based on a proper limiting procedure, which is a consequence of the definition of the associated Legendre functions of the first and the second kind [14, 15], revealing thatupon usage of (114), where for any value of it isThe degeneration of the aforementioned low-frequency electric and magnetic fields, in order to recover the case of a spherical metallic body of radius , embedded within a lossless medium, is then a straightforward sequence of steps, since any kind of indeterminacies is absent.

4.2. Numerical Implementation

We adopt the appropriate prolate spheroidal geometry to describe a perfectly conducting () spheroidal body of major and minor axis and , respectively; therefore, taking the semifocal distance as , the characteristic prolate spheroidal variable is . We wish to implement the most frequently measurable field, which is the scattered magnetic field , given by expression (21), implying the low-frequency solutions (58), (76), and (83). The particular geometry we use to represent the metallic target, which is embedded in a homogeneous infinite space of permeability and permittivity (values of vacuum), is shown in Figure 1, whereas the axis of symmetry is comprised by the -axis.

In order to provide understandable and symmetric diagrams so as to validate the achieved relations for the implicated fields, we choose a magnetic dipole source of strength A·m2, which is located at and illuminates the prolate spheroid at the low frequency of , where , while the magnetic field is measured along a line parallel to the -axis, set at . We show numerical results by providing illustrations for the real and the imaginary parts of the complex approximated scattered magnetic field in units , readily shown in the sequel within Figure 2.

Additionally, under the aim of providing numerical demonstrations for the limiting occasions of the prolate spheroid, we also plot the two interesting limits of the needle with length and the sphere of radius . Furthermore, in order to compare all three cases, we insert in the same figure (Figure 2) both the needle and the sphere circumstances, along with the general case of a genuine prolate spheroidal object.

The series solution obtained for the evaluation of converges relatively fast as , while it must be pointed out that a few terms of the series expansions are enough for determining the scattered magnetic field within a very good accuracy. Moreover, the number of the adequate terms needed to converge and obtain a specific accuracy depends upon the ratio of the major axes over one of the two minor axes of the prolate spheroid, which means an increase of as this ratio increases. Certainly, the number of the utilized terms is increased, a fact that it is partly justified by the necessity of the applied cut-off techniques. Further discussion involves the numerical behavior of the analytical solution, presented in Figure 2. Obviously, the -components are absent, since we have chosen our observation line at . As expected and in accordance with the position of the dipole, the results are symmetric with respect to the line . On the other hand, all the three situations behave similarly, while in terms of magnitudes there is no significant change between the prolate spheroidal and the spherical case, which is not the fact for the case of a needle, meaning that, for extreme circumstances, there is a slight difference in the final result.

5. Conclusions and Discussion

An analytical method to the electromagnetic low-frequency scattering problem for a perfectly electrically conducting spheroidal body in a conductive surrounding, which is illuminated by a magnetic dipole with arbitrary close to the object location and orientation, is developed. The general methodology for the identification of the metal spheroid is presented, based on rigorous, yet versatile, mathematical tools. Its highly conducting morphology is taken into account via impenetrable conditions, while the limiting behavior of the fields at infinity is readily secured. This radiation approximation of a vanishing field at sufficiently long distance from the source creates an error of minor significance. Upon definition of the medium’s wave number, the three-dimensional fields assumed expansions at low frequencies, where a few terms are needed to provide an accurate description of the phenomenon. The spheroidal coordinate system was best fitted to the physical configuration of the problem and therein the 3D solution has been obtained in terms of handy spheroidal eigenfunctions. The commensurable spherical occasion is recovered via a careful limiting procedure, which demonstrated our cumbersome analysis at hand.

Suchlike analytical solutions and formulae attaining closed-type forms have quite important advantages compared with the pure numerical methods. Indeed, the validity of numerical solutions can be verified by analytical or semianalytical techniques. On the other hand, bearing in mind that very important physical laws can be derived from analytical methods, we can understand the necessity of tackling with a confident mathematical basis before incorporating with an algorithmic procedure. Therefore, even nowadays, there is always room for such kind of methods that coexist with pure numerical codes and aim to the solution of boundary value problems in physical applications of major importance.

Mathematical and computational work is currently in progress and involves research into several directions, such as the introduction of more complicated geometries for the metallic target or the accomplishment of actually difficult inversion algorithms, taking profit from the proposed framework.

Conflicts of Interest

The author declares that they have no conflicts of interest.