1. Introduction
Highly directive electromagnetic fields involved in point-to-point (P2P) communications represent a desired aspect of information/energy transmission, in those applications where not being intercepted is important. This aspect is not only required for securing a private communication channel, but also to increase the efficiency in energy transfer, thereby limiting the release of a great amount of wasted energy in neighboring media. Typical commercial solutions involve the beaming of microwaves through parabolic reflectors, laser beams carrying encoded data and short distance Wi-Fi communication systems. Free space optical (FSO) communications are used for short distance, high bandwidth transmissions. They typically rely on coherent sources, but suffer from atmospheric fluctuations in the moisture content (particularly fog), severely affecting the link. Therefore, more promising applications of FSOs are found inside chips [
1]. Radio Frequency (RF) links are preferred to secure higher bandwidths over longer distances and for extended periods of time. Nevertheless, their energy efficiency is lower in comparison to FSOs, and this aspect justifies our research.
Nowadays we may cite the ThoR consortium, active in Europe and Japan, working with state-of-the-art chip sets and modems operating in the standardized 60 and 70–80 GHz bands. Very recent applications (see [
2]) allow one to build RF wireless P2P links featuring >100 Gbps over 1 km at 300 GHz. The quest for high bandwidth links has increased the center frequency, so currently radar bands [
3] and millimeter waves (THz) [
4] are studied as possible extensions of 5G and beyond technologies, further strengthening the importance of establishing directive links. It has to be noted that, as the energy increases with center frequency, so does the concern for impacts on human health—controversial results have been found [
5]. Another aspect of relevance is wireless power transmission (WPT), fueled by the escalation of Internet-of-Things (IoT) (see [
6,
7]). Here, a crowd of sensor nodes dispersed in a physical environment requires either a distributed energy harvesting capability, or the wireless transmission of energy from a central harvesting unit. In WPT applications, three approaches have been developed: (i) near-field coupling of inductive or capacitive loads in resonant or non-resonant mode; (ii) far-field directive power beaming; (iii) far-field non-directive power transfer. Item (i) is used for extremely small scale environments; (ii), which is our aim, is used for intentional power transfer from nodes in a broader network and implements radio frequencies or microwave means; (iii) is used for electromagnetic energy harvesting of sources nonintentionally conceived for the scope (see, e.g., [
8]). The benefits of using directional antennas instead of omni-directional ones in wireless sensor networks (WSNs) have been pointed out, e.g., in [
9], in order to significantly reduce interference between network components. Several models of directional antenna in the framework of WSNs are compared in [
10]. There are currently two categories of commercial directive radio antennas: On one hand are the traditional ones, including helix, log periodic array antennas, aperture horn antennas, reflectors and patch antennas. On the other hand are smart antennas, each consisting of an array of elements that can offer adaptive transmission. Depending on the specific application, shapes, sizes and designs of directional antennas can be quite different. The signal pattern from a directional antenna consists of a cigar-shaped great lobe, usually surrounded by smaller side lobes. For a survey on radiation patterns of these devices and an application on wireless networks, we mention, for instance, [
11,
12].
4. An Extended Formulation for Electrodynamics
In order to fix the notation, we start by writing the set of Maxwell equations in SI unit convention:
where, as usual,
is the displacement field,
represents the charge sources in the medium,
is the magnetic strength field,
is the current source in the medium,
is the electric field and
is the magnetic induction field. In a vacuum and in the absence of field sources, the right hand sides of Equations (1) and (2) are set to zero. In many relevant cases, these equations are coupled with some constitutive relations—for instance:
(
being the polarization density) and
(
being the magnetization). As is customary,
and
denote the dielectric constant and the magnetic permeability in a vacuum, respectively. To complete the set, it is customary to add the equation expressing the Lorentz’s force on a given charge
q traveling in an electromagnetic field with velocity
; i.e.,
We now introduce the extensions of the equations as proposed in [
14,
15], wherein the possibility to have nonzero divergence (
) for the electric field in a vacuum, even in the absence of sources, was postulated. This assumption produces an extra current term in the Ampère’s law. As a consequence of this revision, two new variables and two new unknowns were added to the standard system. The main motivation for this extension was the possibility to expand the set of possible solutions. The reviewed model has the advantage of incorporating solitary signal-packets with compact support (solitons) within the framework of a differential theory. The term "soliton" commonly refers to solutions of certain nonlinear equations that are initially localized in space and evolve while maintaining their shape (see, e.g., [
16]). In our case, these will be smooth electromagnetic waves with compact support traveling (in the simplest situation) along a straight trajectory at the speed of light. Note that these entities are not available in the solution space of the classical Maxwell equations. In suitable circumstances, these solitons may be assimilated to “photons”, as those emanated or absorbed by matter upon proper solicitations. With this interpretation, the photon is not the “carrier” of the electromagnetic field anymore, but becomes a pure wave, carrying along electric and magnetic components. In the present paper, we do not want to stress the reader with these aspects. More insights on these issues can be found in [
15].
For the time being, let us suppose we are working in a vacuum. According to [
14], by letting the divergence of the electric field be different from zero, we arrive at the following formulation:
where, by definition,
is the divergence of the electric field. This nonlinear system counts ten unknowns, namely, the three components of the electric field
, the three components of the magnetic induction
, the three components of the wave velocity field
and the scalar function
p. Basically, here we are only considering the case of a linear homogeneous medium, corresponding to
and
, so that there is room for further extensions that for the moment are not taken into account.
The first Equation (6) is Ampère’s law where the additional term
has to be interpreted as a free flow of immaterial current with density
, associated with the evolution of an electromagnetic phenomenon. In principle, this latter term can be summed up to the classical one (not present in (6)), due to the flow of current along a conductive medium that may interact with the wave. The complete version of such a revised equation is reported in [
15], Section 2. By taking in to account these additional contributions, it is possible to model the mutual interactions between wave and guides, as studied later in this paper.
Equation (10) is basically the Euler’s equation for non-viscous fluids. On the left-hand side we find the divergence of the electric field, multiplied by the sum of two terms: the first one is the total derivative
; the second one,
, recalls the Lorentz’s law. The constant
appearing in (10) has the dimension of a charge divided by a mass and has been estimated to be approximately equal to
Coulomb/Kg (see [
15], Appendix G). On the right-hand side, there is the gradient of a potential
p, which analogously to the Euler’s equation, plays the role of pressure. We may actually treat the potential
p as a kind of pressure term that can assume both positive and negative values, and has dimensions equal to velocity over a time squared. Indeed, from a dimensional analysis it turns out that
is a force per unity of surface. Equation (9) follows from energy conservation arguments. This suggests the possibility of inducing pressure effects as a consequence of a lack of orthogonality between
and
(see [
17], where the model equations have been used to simulate some electromagnetic tweezer effects). It is crucial to observe that, when
,
and
p are identically zero, one recovers the usual Maxwell Equations (1)–(4) in a vacuum. Therefore, the revised model is actually an extension of the classical one.
Equations of the type of (10), coupled with Maxwell’s ones, are usually found in some plasma physics models—for instance, those related to magnetohydrodynamics (see, for instance, [
18] or [
19]). In such examples, a number of charged particles evolve as a real fluid. In this fashion,
takes the meaning of a fluid flow velocity and
p is a proper pressure. Moreover, a mass density
and a charge density
appear in the equations. In the present model, massive charged particles are not necessarily involved. As a consequence,
,
p and
assume different meanings.
As in the case of Maxwell’s equations, by multiplying the first two equations by
and
respectively, and summing up, we get an energy balance relation:
Finally, by taking the divergence of Equation (6), one finds the continuity equation:
which appears as an automatic consequence of the model equations.
In view of operating in cylindrical coordinates, we propose a suitable simplification based on symmetry arguments. This reduces the starting problem from three to two dimensions. We use the referring system
, where the corresponding functions do not depend on
. We set
,
and
. In this way, the divergence operator in cylindrical coordinates, applied to the electric field, takes the form:
The other meaningful equations are:
5. Free-Waves
An interesting subset of solutions is obtained when
and
p are vanishing, though
will be allowed to remain different from zero. These will be called free-waves, and in general they have the property that the triplet
is orthogonal with the additional conditions
and
. In this case, Equation (10) reduces to:
from which one could easily deduce the orthogonality of the couple
and
, and that of the couple
and
. As it is shown in Chapter 2 in [
14], the reduced system of equations is invariant under Lorentz’s transformations and contains several interesting solutions not available in the Maxwellian context. We are going to introduce those that are more relevant in this discussion.
In the three dimensional cylindrical coordinate system (
), we set:
Here
g and
f are arbitrary functions (
g is only required to tend to zero as
r goes to zero). Waves of this kind solve the whole set of Equations (6)–(10). The solution travels in the direction of the
z-axis at speed
c (see
Figure 1, left). By choosing both
f and
g with compact support, we are able to describe solitonic solutions similar in their behavior to “photons” of any size. Some plane waves of infinite extension trivially belong to this typology. The expressions in (21) are at the foundation of the research carried out in this paper. As a matter of fact, the device proposed in
Section 9 is finalized to generate this type of signal, starting from the assumption that the revised model actually predicts the possibility of such focused compact waves.
A more general situation (see again
Figure 1) is obtained by defining, for a given angle
:
where
. The choice
, corresponds to the case examined above. In particular, later we will be interested in the case where
. Here below, we just mention the case
, where the three Equations (14)–(16) reduce to:
One of the advantages of the extended model is the possibility of handling situations where
and
actually depend on
z (on the contrary, the first equation in (23) would be trivially satisfied). Nevertheless, even by neglecting the dependence on
z, the construction of exact solutions is rather complicated. The problem has been approached by various authors. In [
20,
21], efforts were made to come out with families of solutions that are successively generalized in view of applications involving waves in nonlinear media. The expressions found are nonlinear and couple the unknowns
and
in such a way that each one is a function of the other. Alternatively, expansions in terms of Bessel’s functions are available. The consequence of those studies is that the setting proposed in (22) is just a simplification. Indeed, we guess that
itself nonlinearly evolves in time depending on the behavior of the electromagnetic fields. We did not pursue this analysis further. We instead performed some numerical computations based on (22), since we believed that such a setting would not alter significantly our outcomes, at least if we were suitably far away from
.
We end this section by providing some exact solutions in the spherical context. Thus, we work with a spherical system of coordinates
. An electromagnetic wave displaying perfect axial-symmetric spherical fronts can be expressed by the following triplet:
where
f and
g are arbitrary smooth functions. These are genuine free waves, according to our definition. Note that the electric and magnetic fields have no component along the radial direction. Moreover:
. Since
g does not depend on
, the behavior is similar to that generated by an infinitesimal oscillating dipole lined up with the
z-axis. This family of solutions moves along the radial direction with velocity
c, and they are modulated by the function
f. The wave-fronts are the surfaces of spheres centered at the origin. By explicitly computing the divergence of the electric field, we get:
It turns out that there are no possibilities to impose the free divergence condition required from Equation (1) in a vacuum, unless one chooses a
g value proportional to
(which implies the existence of singularities at the poles). Therefore, waves with spherical fronts cannot be of Maxwellian type. Note that the well-known spherical Hertz’s solution for an infinitesimal dipole solves the entire set of Maxwell’s equations, including
. Nevertheless, its wave-fronts (i.e., the surfaces enveloping the electromagnetic field) are not spherical and
is not orthogonal to the propagation direction. More comments on this delicate subject are provided in Appendix D in [
15].
Let us also observe that the field in (24) has a nonvanishing divergence; i.e., . Note that this is also true for the field in (22). This aspect reveals that the model equation are not directly related to classical fluid dynamics, where the condition implies conservation of mass. No masses are involved in our case—instead, purely electromagnetic propagation in a vacuum. We are assuming however that there may be regions where is different from zero.
We can further observe that, for R tending to infinity, decreases faster than . This shows that the non-linearity of the equations becomes weaker as the fields are evaluated far from the source. Note that the superposition principle does not hold. Nevertheless, it may be recovered for sufficiently large values of R. This agrees with the fact that the superposition principle is true in many applications involving wave propagation—for instance, in constructive and destructive interference.
6. Numerical Methods
This section is devoted to brief descriptions of the numerical methods used to approximate the system (14)–(19). Other techniques for similar problems were employed in [
17,
22]. Preliminary computations based on this specific cylindrical setting were tried in [
23]. Our approach sticks to the framework of low-order finite difference schemes. More accurate numerical methods may be applied in future developments. For the reasons we are going to describe, we soon discarded the idea of applying the popular Yee scheme (see, e.g., [
24,
25]).
Consider the two Maxwell equations and relative to a region of the void space in absence of sources of any kind. Let us place ourselves in a Cartesian reference frame . Yee’s algorithm does not explicitly implement the equations and . Indeed, the space grids are structured in such a way to reproduce only the action of the curl operator. As a consequence of this choice, the divergence-free conditions are implicitly satisfied.
In a simplified context, we may assume that the fields are of the form
,
and do not depend on the variable
z. The scheme uses a system of two staggered grids of step
h. The fields are approximated on such grids, in order to ensure that every
component is surrounded by four
components, and vice versa. Time discretization is performed by fractional steps, by alternating the evaluation of
and
. In Cartesian coordinates with grid points
, the procedure reads as follows:
where the index
k refers to time. The numerical divergence of the electric field is calculated as follows:
It is easy to check that
, for any
, which means that such a quantity is preserved during the time iterations. Many approximation schemes for the electrodynamics equations trust the fact that the divergences will not be modified during the evolution. Unfortunately, this may be true only for isolated regions of space. When any wave interacts with other objects, for example, incompatibilities grow by putting together the evolution equations, the divergence-free conditions and the boundary constraints. This fact has been made evident in [
26,
27], where it is shown that these restrictions are too many, so that they cannot hold together at the same time. By allowing
to be different from zero we can come out from this dead end, thereby justifying the adoption of the consistent model we are examining in this paper. At this point, it should be clear that Yee’s scheme and similar ones are not appropriate for our study.
We then adopt a classical explicit (first-order) Euler method for advancing in time. The space derivatives in (14)–(16) are approximated by centered (second-order) differences. In cylindrical coordinates, we end up with the scheme at the grid point
:
For the computation of the divergence of the electric field, we used the following approximation of (13), valid in cylindrical coordinates:
To ensure stability we must require a Courant–Friedrichs–Lewy (CFL) type condition—i.e., that the time-step is proportional to the space discretization parameter h. Such a condition is not particularly restrictive for our purposes.
Equations (9) and (10) contain the time derivatives
,
. It is quite simple to come out with appropriate numerical schemes for their discretization. However, in this paper we are mainly concerned with the case of free-waves (see
Section 5), so that
p and the total derivative
will be set to zero at any time. Since, in practice,
and
will stay orthogonal, Equation (9) becomes trivial, whereas (10) is substituted by (20). In the particular setting we are interested in, the vectors of the triplet
will constitute an orthogonal frame and satisfy
and
. What we just wrote is always true with the exception of some transition regions that will be described later on.
Let us notice that the approximation of (10) in its general form is far from being trivial. First of all, there is no viscous term, a property that makes the theoretical analysis and the numerical implementation quite hard. Secondly, as already pointed out in the examples of
Section 5, we also have
. This means that many numerical schemes available for the discretization of the equations of fluid dynamics may turn out to be unsuitable. Finally, we observe that there may be points in the computational domain where
is zero (or relatively small). This is usually true in correspondence with points where
is also zero. As a consequence, the unreliability of the ratio
may give origins to severe instability problems in the treatment of (10). As mentioned above, here we will stay away from these situations; nevertheless, such problems exist and should be appropriately pondered within the general context.
Later in this paper, we will study the transmission of an electromagnetic wave across the interface separating media displaying different conductive properties. We describe what happens inside each medium with the help of the equations relative to free-waves. In the instant of transition, the wave is not of free type, so that suitable corrections on are taken into account. For simplicity, we will continue however to neglect the role of the potential p.
7. Design of the Directional Antenna
The geometry we are willing to test is shown in
Figure 2. The antenna device is constituted by two conductive cones separated by a little gap. These are supplied by an alternating current via a coaxial cable (not shown in figure), carrying a frequency whose wave-length, in most applications, is of the order of the magnitude of the conductive support size. Each cone features an angle of
at the vertex. The need for a gap, possibly very small, between the cones is dictated by the size of the wire connecting to the electrical source. We presume that a gap of approximately 5% of the entire size should not produce significant negative effects, but, of course, this also depends on the degree of directivity we would like to achieve.
We consider here the case where the entire antenna is immersed in a region occupied by a linear homogeneous medium
, having phase velocity of the electromagnetic radiation lower than the speed of light in vacuum. Such a medium is characterized by a permittivity coefficient
and features non-magnetic properties, so that the magnetic permeability coefficient will remain
. In
Figure 2, the shape is designed depending on some refraction properties that we are going to describe here below.
Our aim is to have a signal emerging from the antenna perfectly horizontally. In this way, the electromagnetic field follows the behavior of
Figure 1 (central picture), while it remains in
. After the transition to the outside medium (vacuum or air, in this case), the message follows instead the behavior of the right-hand side picture of
Figure 1. In terms of the velocity vector
, we have the following situation. Inside
, the magnitude of
is:
. In the open space, we have as usual:
. At the interface between the two media, we enforce the rules of geometrical optics. We denote by
the border of
, i.e., the interface between the dielectric and the outside vacuum. The determination of the surface
follows from the solution of a differential problem as explained hereafter.
As shown in
Figure 3, in the Cartesian framework
,
corresponds to the graph of a function
defined through the following constraints. First of all, the curve passes through the point
. We have
when the cones have the vertex angle equal to
(here we are neglecting the gap between the cones). A light ray comes from the origin with an angle
. As we said, the curve
s must be designed in such a way that the emerging ray is horizontal. To proceed with the determination of
s we need to recall Snell’s law (see [
28], p. 38). In
Figure 3,
denotes the difference between
and the angle determined by the direction normal to
. Similarly,
denotes the angle between the normal direction and the abscissa. In this particular situation, Snell’s law says that:
We can translate in formulas the above settings by writing:
where the prime denotes the derivative with respect to the variable
r. Using the trigonometric identity
, we get:
By recovering
we finally arrive at the nonlinear differential problem:
The above differential equation may be easily solved numerically, thereby obtaining, for different values of the parameter
, the curves of
Figure 4. One has just to pay attention to the fact that the derivative of
s at the intersection with the
r-axis is equal to infinity.
8. Numerical Simulations
This section is dedicated to the presentation of some numerical simulations, obtained by applying the schemes described in
Section 6. The computational domain is the square
, where
is a parameter. We chose by default
. The first variable is related to
r and the second one to
z, so that we are actually operating on a cylinder with a vertical axis. In
Q, we consider a squared grid, spaced by
, where
N is an integer. The total number of nodes is then proportional to
. The speed of light in a vacuum is normalized to 1. The development of the process is studied in the time interval
, for a certain
. In order to guarantee stability, the time-step
will be chosen to be suitably small.
We set up in advance the velocity field
in order to have
inside the dielectric, where we suppose that the vectors are emanating from the origin. We have
on the remaining part of the domain where the vectors are horizontal (perpendicular to the
z-axis). The deployment of the velocity field is displayed in
Figure 5, for the cases when the dielectric is missing (left) and when it is present (right). In a more rigorous framework,
should not be fixed in advance, but computed as time evolves according to Equations (18) and (19). As we specified in
Section 6, this approach is not easily implementable from the numerical viewpoint, especially when the time advancing procedure is the trivial Euler method. Numerical solutions, in this circumstance, are too much affected by viscosity, so that the results are not reliable. Since the development of more sophisticated techniques is outside of the goals of this paper, we do not insist further on these generalizations. For the above reasons, we decided to assume that
is known a priori.
As we are interested in the evolution within the region determined by the two cones of the antenna, suitable constraints are placed in the computational domain. For simplicity, we assume that these guides are perfect conductors, though the situation might be not realistic. Additionally, we should distinguish between the nodes inside the dielectric and those outside. In the transition between the two media the fields are updated by enforcing that
remains continuous, whereas
has a jump (these are consequences of the Fresnel’s conditions; see, e.g., [
28], p. 40). In practice, we truncate to zero the longitudinal component
at the first node just outside the dielectric.
As an initial guess, we take a portion of spherical wave. Based on (24), we set
for
. This modulates the displacement along the transverse direction. Concerning the longitudinal direction, we take
with
and
. The resulting electric field is shown in
Figure 6 (left). The corresponding shapes of
can be seen in the very first plot of
Figure 7 and
Figure 8. Note that this displacement has non-zero divergence (
). The corresponding initial magnetic field is orthogonal to the plane of the figure and we have
. Successively, the wave evolves according to the model equations. The single initial bump displays a wave-length comparable to the dimensions of the antenna. Of course, one may also consider a train of pulses lined up one after the other. This option provides however results perfectly similar to those pertaining to the single pulse.
The boundaries are placed at a 45 degrees angle; thus, the imposition of the constraints comes very naturally. Neumann boundary conditions for the electric field are imposed on the antenna guide. This is done by respecting the perfect conductivity relation. Again referring to
Figure 6 (right), at each time-step, the (red) vector of components
in proximity of the boundary is projected onto the direction normal to the surface, giving birth to the (green) vector
. This last vector is successively shifted up to the corresponding boundary grid point.
The parameters relative to the numerical tests are provided in the figure captions. As far as
Figure 7 is concerned, the simulations rely on the field
given in the first plot of
Figure 5. At the grid points where
is not defined, we are actually solving the system (6) and (7) with
. At the end of the guide the solution turns out to be discontinuous, thereby producing sensible numerical instabilities. The general behavior of the wave is however reasonably well reproduced.
Concerning
Figure 8, this is obtained from the field
as given in the second plot of
Figure 5. Again we set
at the remaining grid points. The solution is clearly focused in the horizontal direction. The approximation suffers from the discontinuity across the dielectric boundary. Improvements could be obtained through the adoption of a more accurate discretization scheme. Anyway, our purpose here was to show that numerical simulations are feasible and produce a consistent prediction.
Finally, in
Figure 9, we show the results when
is totally eliminated. We are not, however, enforcing the condition
. Since the initial guess itself does not satisfy
, we cannot claim that we are actually solving the full set of Maxwell’s equations in a vacuum. We are not even solving the classical wave equation, since we can just deduce that:
. On the other hand, there are very few initial data recoverable from (24) such that
. This is true, for instance, for
, which does not correspond with a consistent choice for a signal emerging from an antenna (see the discussion in [
14], Chapter 1).
The behavior of the solution shown in
Figure 9, especially regarding the second plot (to be compared with the second plot of
Figure 7), is very inefficient. As far as the successive plots are concerned, some of the energy remains concentrated at the endpoints of the guide. Noticeably, in most of the computational code available on the market, the checking of the zero divergence condition in vacuum is disregarded. This means that, if
at a certain time starts assuming nonvanishing values (because of imposing some boundary constraints deriving for example from a scattering object), the problem is not approached in the proper way. This is true, for instance, in the case of Yee’s scheme, mentioned in
Section 6, whose implementation strictly trusts the fact that
must remain zero for all times. The real trouble is that Maxwell’s equations are too restrictive. In short, the solution of the vector wave equation, in a given domain with generic boundary conditions, has no chance to have its divergence be zero. That is why we adopted the extended model in our simulations.
We end this section with a few remarks. We placed our initial guess at a distance from the real source (the gap between the two cones). In truth, we do not exactly know how the signal is actually created (by going backwards the initial guess degenerates into a singularity). This uncertainty is a consequence of the fact that the mechanism of antenna emission is still not completely clear. First of all, the guides should not be assumed to be perfect conductors, a property that involves a nontrivial discussions about boundary conditions and the way they should be approximated (see, e.g., [
29,
30,
31]). Such more complicated situations activate the whole set of modeling equations. In fact, we start with missing the orthogonality between
and
. This generates values of the potential
p different from zero, influencing in this way all the other variables. Secondly, there is a resonance process (depending on the size and shape of the device) that allows for the propagation of waves only if they are above a certain frequency. These issues are just sketched in this paper and are often underestimated in the design of a device. A transmitting antenna is an object transforming a certain signal by conferring to it the capability to travel in free space without further help from the source (see the final comments in the preliminary paper [
32]). We believe that applied physics necessitates a better theoretical understanding of the mechanism of antenna emissions. Indeed, more accurate computations cannot neglect this crucial aspect.
9. Highly Directional Emissions
The computations of
Section 8 have been realized for a relative permittivity coefficient
, corresponding, for instance, with polyethylene. This is in fact the material that we would like to use in laboratory experiments, due to the possibility of shaping the dielectric through a 3D printing process.
A further conical metal reflector can be added to the device of
Figure 2, obtaining the final antenna schematically represented in
Figure 10. The angle at the vertex of the reflector is again
. In this way, the horizontally spreading signal is driven in the vertical direction. Referring to
Figure 1, during the whole process, an annular wave starts inside the dielectric following a path like that shown in the central picture. When exiting the dielectric, the new Poynting vectors are lined up to the vectors
, as depicted in the picture on the right-hand side of
Figure 1. Finally, after the last reflection, the signal is converted into a highly directional beam, such as that of the picture on the left-hand side of
Figure 1.
During all these transitions there is no breaking of the stream lines of the magnetic field, which remain closed curves encircling the
z-axis. After projection on the plane orthogonal to the z-axis, the evolution of the wave displays a central symmetry. This peculiar property is not commonly shared by other devices—for example, parabolic reflectors (when the source is represented by a classical dipole) or horn antennas. These devices actually need to break their initial symmetries in order to radiate [
33].
Preliminary experiments without the presence of the dielectric showed that the signal emerging from the reflector has an annular shape. The fields are thus polarized in a sort of a circular fashion. The signal is absent nearby the z-axis. The global shape of the emitted wave is that of a torus. The directivity is relatively good, but, due to the lack of the dielectric lens correction, the toroidal wave grows with distance from the source, following a conical pattern. In those tests, the cones were replaced by two conductors having a complicated shape, in order to skip the construction of the dielectric medium . The outcome was interesting but still far from the directional signal we would like to achieve. The above experiments are however quite crucial in confirming some theoretical predictions. Waves as those just described do not exist in the Maxwellian framework. Indeed, the request that the divergence of the electric field must be always zero in a vacuum does not allow the inclusion of such special waves in the solution space. This means that the adoption of the revised model is rather important in the description of antennas already existing, opening the path to the design of a new generation of sustainable, efficient and performant devices.
It finally has to be noted that, if perfect directional signals exist, where the rays are all parallel, the natural receiver for this kind of antenna is the antenna itself. Indeed, as in a mirror image, the parallel rays are converted back and concentrated at the gap of the cones. There, they can be finally transformed in potentials.
10. Conclusions
We guessed the possibility of generating electromagnetic waves with extremely high directivity. This prediction is suggested by an extension of Maxwell’s equations in a vacuum aimed to simulate compact signals, traveling unperturbed at the speed of light. The enlargement of the solutions space relies on the fact that the divergence of the electric field may assume values different from zero also in the absence of charges associated with physical particles—for instance, electrons. This requirement is indeed able to explain many electromagnetic phenomena, as documented in [
15]. In addition, the revised model provides the exact link between electrodynamics equations and the rules of geometrical optics, a property that the standard Maxwellian approach is not able to reproduce, unless one accepts a sequence of rough simplifications (see, e.g., [
28], Section 3.1).
The waves belonging to the category we studied in this paper display a toroidal topology, with the magnetic field distributed along closed lines following the toroidal axis. With the help of simple devices, it is actually possible to create waves having a similar topology and developing along conical patterns. These emanations cannot be described by the classical Maxwell equation, so that the adoption of the extended model becomes a necessity dictated by reality.
The theory predicts that antennas with infinite directivity could possibly be built. Reality is a different issue, so reasonable quantitative confirmations will only be made as the results of the first experiments are made available. Once it is demonstrated that such a project is feasible, generating these waves becomes only matter of technical effort. The reader can of course understand the importance of this achievement. In this paper, we indicated a way to fabricate antennas devoted to this purpose. These are obtained by filling the space delimited by a biconic antenna with a dielectric medium placed in contact with the cones, and such that the boundary between the medium and the external vacuum is designed following a suitable curvature. The specific signals produced are expected to reach a prescribed target without dissipating their energy along the path. Moreover, since most of the theoretical passages are based on the rules of geometrical optics, we argue that the role of frequency might not be so crucial, thereby allowing transmissions ranging within a rather broad band. All these properties should guarantee a number of positive outcomes in terms of sustainable operation, safe energy transfer and secure information beaming.