Interference in Phased Array and Conservation of Energy

Aug 24, 2022·
Zakir Hussain Shaik
Zakir Hussain Shaik
· 16 min read

Have you ever encountered a question that looks innocently subtle but still haunts your brain for justification of an apparent paradox? One such riddle that most students studying EM waves or antenna theory come across is how power behaves when coherent sources interfere, as shown in the figure below.

The law of conservation of energy is one of the fundamental laws of physics, and there is no known phenomenon to date where it has failed. This means there must be some flaw in the argument presented in the above figure. This post addresses this question from two perspectives. First, it examines what happens to energy when individual EM waves interfere constructively or destructively. Second, it explores how radiated power behaves when we have multiple phased antennas.

Problem description of an apparent paradox

First, let us formally present the apparent paradox depicted in the figure. Assume case 1: where we have a sinusoidal signal x1(t)=Asin(2πfct)x_1(t) = A\sin(2\pi f_c t) and the power of the signal is obtained by using the formula given below, which is most commonly utilized in signal theory books as

Power=1T0Tx(t)2dt{\rm Power} = \frac{1}{T}\int_{0}^{T} \vert x(t) \vert^2 dt

where AA and TT are the amplitude and time period of the signal x(t)x(t), respectively. Accordingly, the power of signal x1(t)x_1(t) is P1=A2/2P_1 = A^2/2 units. Now consider case 2: where we have a constructive combination of two exactly identical signals. By the superposition principle, we have x2(t)=2x1(t)=2Asin(2πfct)x_2(t) = 2x_1(t) = 2A\sin(2\pi f_c t). Using the power definition, the power of the combined signal is P2=4A2/2=4P1P_2 = 4A^2/2 = 4P_1 units. However, the total power of the individual signals is only 2P12P_1 units. An obvious but strange question arises: where did the extra 2P12P_1 units of power come from? This gives birth to the paradox and the apparent failure of conservation of energy.

Short Answer:

  • If the model in context is wired, then it is impossible to build such a situation.
  • If the model in context is a phased array, then the received power is always less than radiated power.
  • If the model in context is the intersection of two EM waves traveling exactly in opposite directions, even then the energy is conserved because whenever electric energy is doubled the magnetic energy is zero and vice-versa. The total energy in the wave is the sum of electric and magnetic field energy. Thus, energy is conserved.

Mathematical justification for the above statements follows next.

Equivalent circuit for the above phenomenon

The example to be presented is a first-counter argument that it is impossible to build a circuit to give the above paradox. Moreover, if something is impossible in a wired circuit, then it is impossible in wireless. Let’s consider two cases similar to the previous discussion. Case 1: consider an AC circuit with voltage source V1(t)=Asin(ωct)V_1(t) = A\sin(\omega_c t) (units are in Volt) as shown below:

Network 1

The AC current flowing in such circuit is I1(t)=Asin(ωct)I_1(t) = A\sin(\omega_c t) (units are in Ampere). Now, if we compute the average power dissipated across the resistor:

P1=1T0TV1(t)2dtP_1 = \frac{1}{T}\int_{0}^{T} V_1(t)^2 dt

this results in P1=A2/2P_1 = A^2/2 W, and the power delivered by the voltage source is

Ps1=1T0TV1(t)I1(t)dtP_{s1} = \frac{1}{T}\int_{0}^{T} V_1(t)I_1(t) dt

which also gives Ps1=A2/2P_{s1} = A^2/2 W. So far, so good. Now take the second case

Network 2

The total voltage across the resistor is V2(t)=2Asin(ωct)V_2(t) = 2A\sin(\omega_c t) V, and the total power dissipated across the resistor

P2=1T0TV2(t)2dtP_2 = \frac{1}{T}\int_{0}^{T} V_2(t)^2 dt

which gives P2=4A2/2P_2 = 4A^2/2 W. Thus, the circuit has indeed delivered 4P14P_1 W, similar to the claim in the paradox. However, in this case, the power delivered by each source, say Ps1P_{s1} and Ps2P_{s2}, is not A2/2A^2/2 W as was the case when they were operating individually. The power delivered by each source is computed as follows:

Ps1=Ps2=1T0TV1(t)I2(t)dtP_{s1} = P_{s2} = \frac{1}{T}\int_{0}^{T} V_1(t)I_2(t) dt

where I2(t)=2Asin(ωct)=2I1(t)I_2(t) = 2A\sin(\omega_c t) = 2I_1(t). That is to say, when connected in series with another voltage source with zero phase difference, the same voltage source delivers more current than when acting alone. With the above calculation, Ps1=Ps2=2A2/2P_{s1} = P_{s2} = 2A^2/2 W. Thus, each voltage source is delivering twice the power it delivered when operating alone. Hence, the law of energy conservation holds. One important lesson from this example is that the power delivered by sources is not always proportional to the square of voltage or current. Rather, it is always safer to compute the product of voltage and current whenever there are multiple sources. In the example presented above, the voltage and current amplitudes are not the same.

Disclaimer: further reading is intended for readers who are interested to know how interference behaves in wireless communication. Some fundamental background in EM and basic antenna theory is expected.

Phased Antenna Array

Now we shall see how multiple antenna elements radiating together behave. Here too, we shall start with two cases: the first one involves only one source, and the second involves two sources and then compute received and radiated powers. For the sake of explanation, we shall use only isotropic antennas, although they don’t exist but are often used as references.

Case 1: We consider here a single isotropic source. Let the electric field of an isotropic antenna in the far field at a reference distance of rr in terms of elevation angle θ\theta and azimuthal angle ϕ\phi be E1(θ,ϕ)=E0rE_1(\theta,\phi) = \frac{E_0}{r}, and let the magnetic field at that reference distance be H1H_1. Then, the time-average Poynting vector is given as

S(θ,ϕ)=12{E1×H1}=12Z0E1(θ,ϕ)2[W/m2] \begin{align} S(\theta,\phi) &= \frac{1}{2}\Re\{E_1\times H_1^{*}\}\\ &=\frac{1}{2Z_0}\vert E_1 (\theta,\phi) \vert^2\quad[{\rm W/m^2}] \end{align}

where H1H_1 is the magnetic field at the reference distance. Note that in the far field E1/H1=Z0\vert E_1\vert/\vert H_1\vert = Z_0, where Z0Z_0 is the intrinsic impedance of free space. Then the radiation intensity is U(θ,ϕ)=r2S(θ,ϕ)U(\theta,\phi) = r^2S(\theta,\phi) [W/sr][\rm W/sr], which is also called the angular distribution of time-averaged power and denoted as P(θ,ϕ)/Ω\partial P(\theta,\phi)/\partial \Omega. The total radiated power can now be computed as follows:

P1=02π0πU(θ,ϕ)sin(θ,ϕ)dθdϕ=12Z0E024π=12.56A12[W] \begin{align} P_1 &= \int_{0}^{2\pi}\int_{0}^{\pi}U(\theta,\phi)\sin(\theta,\phi)d\theta d\phi\\ &=\frac{1}{2Z_0}\vert E_0 \vert^2 4\pi\\ &=12.56 A_1^2\quad [{\rm W}] \end{align}

where A12=12Z0E02A_1^2 = \frac{1}{2Z_0}\vert E_0 \vert^2.

For an isotropic antenna, we can write S=P14πr2S = \frac{P_1}{4\pi r^2} and the effective area of receiving isotropic antenna is λ24π\frac{\lambda^2}{4\pi}, where λ\lambda is the wavelength of the wave radiated by the antenna, then the received power at any direction can also be written as

P1(θ,ϕ)=βP1\partial P_1\left(\theta,\phi\right) = \beta P_1

where β=λ2(4πr)2\beta = \frac{\lambda^2}{(4\pi r)^2} and in textbooks, it is often referred to as line-of-sight (LoS) channel gain for an isotropic receiving antenna.

Case 2: Now consider two isotropic sources and let’s assume there are placed on the y-axis at (0,Δ2,0)(0,-\frac{\Delta}{2},0) and (0,Δ2,0)(0,\frac{\Delta}{2},0), respectively. Let the observing reference is in the far field and in the direction nˉ=(sin(θ)cos(ϕ),sin(θ)sin(ϕ),cos(θ))\bar{n} = (\sin(\theta)\cos(\phi),\sin(\theta)\sin(\phi),\cos(\theta)). Then, the approximate path length difference between the wave travel from these two sources is obtained by nˉdˉ\bar{n}\cdot \bar{d}, where dˉ\bar{d} is the vector distance between two sources, i.e. dˉ=(0,Δ,0)\bar{d} = (0,\Delta,0).

With the assumptions made above, the total phase difference between the waves in the far field would be

ψ=2πλnˉdˉ+ψ0 \psi = \frac{2\pi}{\lambda} \bar{n}\cdot \bar{d} + \psi_0

where ψ0\psi_0 is the initial phase difference between the sources. Now from the principle of superposition, the total electric field will be given by

E2=E1+E1eiψ=2E1eiψ2cos(ψ2) \begin{align} E_2 &= E_1 + E_1e^{i\psi}\\ &= 2E_1e^{i\frac{\psi}{2}}\cos\left(\frac{\psi}{2}\right) \end{align}

Then the effective radiation intensity is given by (Note: if below equations are not rendering properly in mobile browser then in browser options select “desktop site”)

U2(θ,ϕ)=r22Z0E22=12Z04E02cos2(πλnˉdˉ+ψ02)=4A12[cos2(ψ02)cos2(πΔλsin(θ)sin(ϕ))+sin2(ψ02)sin2(πΔλsin(θ)sin(ϕ))12sin(ψ0)sin(2πΔλsin(θ)sin(ϕ))] \begin{align} U_2(\theta,\phi) &= \frac{r^2}{2Z_0}\vert E_2 \vert^2\\ &=\frac{1}{2Z_0}4\vert E_0 \vert^2 \cos^2\left(\frac{\pi}{\lambda}\bar{n}\cdot\bar{d} + \frac{\psi_0}{2}\right)\\ &=4A_1^2 \left[\cos^2\left(\frac{\psi_0}{2}\right)\cos^2 \left(\frac{\pi\Delta}{\lambda} \sin(\theta)\sin(\phi)\right)\right.\\ &\qquad + \sin^2\left(\frac{\psi_0}{2}\right)\sin^2 \left(\frac{\pi\Delta}{\lambda} \sin(\theta)\sin(\phi)\right)\\ &\qquad -\left.\frac{1}{2}\sin\left(\psi_0\right)\sin\left(\frac{2\pi\Delta}{\lambda}\sin(\theta)\sin(\phi)\right)\right] \end{align}

The total radiated power

P2=02π0πU2(θ,ϕ)sin(θ)dθdϕ P_2 = \int_{0}^{2\pi}\int_{0}^{\pi}U_2(\theta,\phi)\sin(\theta) d\theta d\phi

Now we present three different situations and give an interpretation behind them:

Case 2a: Δ=λ2\Delta = \frac{\lambda}{2} and ψ_0=π2\psi\_0 = \frac{\pi}{2}

\begin{align} P\_2 &= 2A\_1^2 \int\_{0}^{2\pi}\int\_{0}^{\pi} \left\[1 - \sin\left(\frac{\pi}{2}\sin(\theta)\sin(\phi)\right)\right]d\theta d\phi \\\\ & = 8\pi A\_1^2\\\\ & = 2P\_1 \end{align}

So, in this case, the total radiated power equals the sum of radiated powers of individual sources. In this case, we say there is no mutual coupling between antennas.

Then one may wonder what is the benefit of having two antennas. The answer lies in the observation that in the case of practical wireless systems, the power is radiated outwards like a balloon but is not pointed as a sharp laser beam. So first consider the total radiated power in both cases (single source and two sources) to be same, i.e. let the total power radiated in a single antenna case be 2P_12P\_1 then let us consider the direction in which the radiation intensity (or average power per solid angle) is maximum, and this happens in any direction for the case of a single source (since it is isotropic) and for the case of two sources θ=π2,ϕ=ϕ2\theta = \frac{\pi}{2}, \phi = -\frac{\phi}{2} (one of the case) and is given as

(P_2(θ,ϕ)Ω)_max=4A_12 =22(P_1(θ,ϕ)Ω)_max 0.16P_2 \begin{align} \left(\frac{\partial P\_2(\theta,\phi)}{\partial \Omega}\right)\_{\rm max} &= 4A\_1^2 \\\ & = 2\cdot 2\left(\frac{\partial P\_1(\theta,\phi)}{\partial \Omega}\right)\_{\rm max} \\\ & \approx 0.16 P\_2 \end{align}

Thus, we can see that in the most favourable direction, the power gain with two isotropic antennas is 22 compared to a single isotropic antenna source. This is called beamforming gain or array power gain. What we mean by beamforming gain of 2 for two elements isotropic antenna array is that we receive 22 times more power when compared to a single source case. In other words, to receive the same power in two cases, two antenna (isotropic) systems need half the power as that of a single source (isotropic) case. However, the fraction of power that is received is (P2(θ,ϕ)Ω)max0.16P2\left(\frac{\partial P_2(\theta,\phi)}{\partial \Omega}\right)_{\rm max} \approx 0.16 P_2 which is still very less than total radiated power. Thus, in wireless as well, there is no paradox.

Similar to the single antenna case, we can write the fraction of power received at the isotropic antenna that is placed in the most favourable direction

P2(π2,π2)=2βP2=4βP1\partial P_2\left(\frac{\pi}{2},-\frac{\pi}{2}\right) = 2\beta P_2 = 4\beta P_1

As such, this case may sometimes lead to an incorrect impression that received power is a fraction of 4P14 P_1. What is actually happening is that in the case of two isotropic antennas, the power pattern is not uniformly distributed on a sphere as was the case with a single antenna source. Rather, it is dumbbell-shaped, and power is just redistributed. In simple words, if we have 2P12 P_1 power in the case of a single antenna case (these to make fare comparison between two cases in terms of radiated power), then all the points on the sphere would receive β2P1\beta 2P_1 power and now if we remove the power from half of the points on that sphere and redistribute it uniformly in other places then we would have two regions with received powers as β2P1+β2P1=β4P1\beta 2P_1 + \beta 2P_1 = \beta 4P_1 and β2P1β2P1=0\beta 2P_1 - \beta 2P_1 =0, respectively. These respective regions are locations of constructive and destructive interference.

Case 2b: Δ=λ4\Delta = \frac{\lambda}{4} and ψ0=0\psi_0 = 0

P2=4A1202π0πcos2(π4sin(θ)sin(ϕ))dθdϕ3.27(4π)A12=3.27P1 \begin{align} P_2 &= 4A_1^2\int_{0}^{2\pi}\int_{0}^{\pi}\cos^2\left(\frac{\pi}{4}\sin(\theta)\sin(\phi)\right) d\theta d\phi \\ & \approx 3.27 (4\pi)A_1^2 \\ &= 3.27P_1 \end{align}

In this case, the total radiated power is larger than the sum of radiated powers of individual sources when acting alone. In this case, we say there is positive mutual coupling between antennas. What is happening is that the radiation of one antenna is influencing the other antenna and overall changing the impedance seen by the two sources, which in turn affects the radiated power. More specifically, in this case, the radiated power by at least one of the antenna source is larger than P1P_1, i.e. the antenna is radiating more power than it was acting alone. To know more precisely how much each antenna radiates depends on the phases of the exicitation currents of the sources and mutual impedences and these concepts are beyond the intended aim of this article.

If in case the total input drive power is just 2P12P_1 then we need to renormalize U2(θ,ϕ)U_2(\theta,\phi) by 2/3.272/3.27 and if the total drive power is P1P_1 then the renormalizing factor is 1/3.271/3.27. In this case, (P2(θ,ϕ)Ω)max1.22(P1(θ,ϕ)Ω)max\left(\frac{\partial P_2(\theta,\phi)}{\partial \Omega}\right)_{\rm max} \approx 1.22 \left(\frac{\partial P_1(\theta,\phi)}{\partial \Omega}\right)_{\rm max} and hence the beamforming gain is 1.22. Thus, it shows if there is coupling between the antennas, then there is a significant reduction in the beamforming gain [1], [5]. In the most extreme case, if Δ0\Delta \approx 0, the total radiated power is four times the single antenna case. This case resembles the same situation as that of the circuit example presented previously, and in principle, each antenna radiates twice the power than it was radiating when acting alone [3-4]. Also, to add if we introduce phase difference in the circuit example, we can produce these three cases in a circuit as well. Thus, the circuit example can be thought of as maximum positive coupling among the sources.

This case of mutual coupling may sometimes lead to the misimpression that the received power can be computed without renormalization and that more power is radiated than the total drive power.

Case 2c: Δ=λ4\Delta = \frac{\lambda}{4} and ψ0=π\psi_0 = \pi

P2=02π0πcos2(π4sin(θ)sin(ϕ)+π2),dθ,dϕ0.73(4π)A12=0.73P1 \begin{align} P_2 &= \int_{0}^{2\pi}\int_{0}^{\pi}\cos^2\left(\frac{\pi}{4}\sin(\theta)\sin(\phi)+ \frac{\pi}{2}\right),d\theta ,d\phi\\ & \approx 0.73 (4\pi)A_1^2\\ &= 0.73P_1 \end{align}

In this case, the total radiated power is much less than the source acting alone. This is said due to the fact that the effective impedance seen by the antenna array is very high, and hence radiated power is reduced. In the most extreme case, if Δ0\Delta \approx 0, then total radiated power is zero, and this can be thought of as a global destructive interference because, at each point in space, the net electric field is zero. However, this doesn’t violate any energy conservation law because the power doesn’t radiate at all and is within the system itself because the effective impedance seen by the antenna array approaches infinity.

So, in general, with two sources, total radiated power can be written as

Ptotal=P1+P2+PcouplingP_{\rm total} = P_1 + P_2 + P_{\rm coupling}

where P1P_1 and P2P_2 is the power radiated by two sources and PcouplingP_{\rm coupling} is the coupling power which can be negative, zero and positive. So in most cases (P1+P2)PcouplingP1+P2-(P_1+P_2)\leq P_{\rm coupling} \leq P_1+P_2. Negative coupling doesn’t mean power is lost but should be interpreted as the power didn’t radiate from antennas. Such examples can also be produced by adjusting the phase difference in wired circuits and studying how power drawn from each source behave.

If you are interested, then you can use the below MATLAB code snippet to play around with Δ\Delta and ψ_0\psi\_0 and see how radiated power varies:

% Antenna array spacing and initial phase difference
delta_lambda = 0.25; % Delta/Lambda
psi_0 = pi;

% Radiation intensity
U = @(x,y)4*sin(x).*(cos(psi_0/2 + pi*delta_lambda*sin(x).*sin(y))).^2;

% Radiated power
Prad = integral2(U,0,pi,0,2*pi);

% Factor of radiated power compared to single isotropic antenna
factor = Prad/(4*pi);

The channel gain and beamforming gains presented were only for identically excited isotropic antenna elements. However, these values depend on the antenna type, spacing between them, initial phase difference and also on excitation of antenna elements. So an interesting question to ponder is how large can the array power gain be, or is there any upper bound?

Constructive/Destructive interference with EM waves

There’s another perspective where the conservation of energy paradox arises, and that is the case when considering EM waves traveling in the exact opposite direction. Even in such scenarios, the energy is conserved. This is because the total energy that an EM wave carries is the sum of energies of electric and magnetic fields. The important point to note here is that whenever there is constructive interference of electric fields, there would be destructive interference of magnetic fields, and the total energy is preserved. We shall describe this mathematically next.

Let the wave be traveling in a lossless medium in the x-direction and let there be linear polarization of electric field and magnetic fields. Assume that electric field is polarized in y-direction i.e., E=E(xct)y^E = E(x-ct)\hat{y} and accordingly the magnetic field B=B(xct)z^B = B(x-ct)\hat{z}. The shapes of the field can be chosen as sinusoidal as it is widely used for representation. Here, cc is the speed of light in a lossless medium such as a vacuum. Also, c=1/ϵ_0μ_0c = 1/\sqrt{\epsilon\_0 \mu\_0} where ϵ_0\epsilon\_0 and μ_0\mu\_0 are permittivity and permeability of the vacuum, respectively.

Then the energy of the EM wave comprises both electric and magnetic components and is obtained by integrating corresponding energy densities as follows

U_E=12_all volumeϵ0E2dVolU\_E = \frac{1}{2}\int\_{{\rm all\ volume}} \epsilon_0 E^2 dVolU_M=12_all volume1μ0B2dVolU\_M = \frac{1}{2}\int\_{{\rm all\ volume}} \frac{1}{\mu_0} B^2 dVol

when a single EM wave is traveling then we have U_E=U_MU\_E = U\_M, and the total energy UU is given by

U=U_E+U_M=2U_E=2U_MU = U\_E + U\_M = 2U\_E = 2U\_M

Now consider the case of interest here, for two waves propagating in opposite directions with similar waveforms

E_1(x,t)=E(xct),B_1(x,t)=B(xct)E\_1(x,t) = E(x-ct), B\_1(x,t) = B(x-ct)

and

E_2(x,t)=±E(x+ct),B_2(x,t)=B(x+ct)E\_2(x,t) = \pm E(x+ct), B\_2(x,t) = \mp B(x+ct)

the reverse signs of magnetic fields can be obtained by using the right-hand rule for the direction of the electric field, magnetic fields and direction of wave propagation. Now total energies before overlapping

Utotal=UE,1+UM,1+UE,2+UM,2=2UE+2UM=4UE=4UMU_{\rm total} = U_{E,1} + U_{M,1} + U_{E,2} + U_{M,2} = 2U_E + 2U_M = 4U_E = 4U_M

The total electric/magnetic fields after superposition are:

Etotal=E1+E2=E(xct)±E(x+ct)E_{\rm total} = E_1 + E_2 = E(x-ct) \pm E(x+ct)Btotal=B1+B2=B(xct)B(x+ct)B_{\rm total} = B_1 + B_2 = B(x-ct) \mp B(x+ct)

Constructive interference (if the focus is electric fields) at t=0t=0

Etotal(t=0)=2E(x),Btotal(t=0)=0E_{\rm total}(t=0) = 2 E(x),\quad B_{\rm total}(t=0)=0UE,total(t=0)=12all volumeϵ0Etotal2dVol=412all volumeϵ0E2dVol=4UE \begin{align} U_{E,{\rm total}}(t=0) &= \frac{1}{2}\int_{\rm all\ volume} \epsilon_0 E_{\rm total}^2 dVol \\ &= 4 \frac{1}{2}\int_{\rm all\ volume} \epsilon_0 E^2 dVol \\ &= 4U_E \end{align}

and

UM,total(t=0)=12all volume1μ0Btotal2dVol=0U_{M,{\rm total}}(t=0) = \frac{1}{2}\int_{\rm all\ volume} \frac{1}{\mu_0} B_{\rm total}^2 dVol = 0Utotal(t=0)=UE,total(t=0)+UM,total(t=0)=4UE=UtotalU_{\rm total}(t=0) = U_{E,{\rm total}}(t=0) + U_{M,{\rm total}}(t=0) = 4U_E = U_{\rm total}

Constructive interference destroys magnetic energy and double electric energy, and as such total energy is preserved.

Destructive interference (if the focus is electric fields) at t=0t=0

Etotal(t=0)=0,Btotal(t=0)=2B(x)E_{\rm total}(t=0) = 0, \quad B_{\rm total}(t=0)=2B(x)UE,total(t=0)=12all volumeϵ0Etotal2dVol=0U_{E,{\rm total}}(t=0) = \frac{1}{2}\int_{\rm all\ volume} \epsilon_0 E_{\rm total}^2 dVol = 0

and

UM,total(t=0)=12all volume1μ0Btotal2dVol=412all volume1μ0B2dVol=4UM \begin{align} U_{M,{\rm total}}(t=0) &= \frac{1}{2}\int_{\rm all\ volume} \frac{1}{\mu_0} B_{\rm total}^2 dVol \\ &= 4 \frac{1}{2}\int_{\rm all\ volume} \frac{1}{\mu_0} B^2 dVol \\ &= 4U_M \end{align} Utotal(t=0)=UE,total(t=0)+UM,total(t=0)=4UM=UtotalU_{\rm total}(t=0) = U_{E,{\rm total}}(t=0) + U_{M,{\rm total}}(t=0) = 4U_M = U_{\rm total}

Destructive interference destroys electric energy and double magnetic energy, and as such total energy is preserved. However, in both the cases poynting vector becomes zero because you have standing wave and hence there is no net energy flow in a standing wave. This just says net flow is zero but there is equal flow of energy in opposite directions and energy is no where lost as shown above. For more on standing waves refer to [8].

Thus, this explanation shows there is no breakdown of energy conservation law.

Thank you for reading through.

Now food to your brains, what happens if equal and opposite electric charges pass through each other?

Some important references for further reading:

[1] Marlow, Daniel, and Kirk McDonald. A Phased Antenna Array. http://kirkmcd.princeton.edu/examples/endfire.pdf.

[2] Drosd, Robert, Leonid Minkin, and Alexander S. Shapovalov. “Interference and the law of energy conservation.” The Physics Teacher 52.7 (2014): 428-430.

[3] Schantz, Hans G. “On the superposition and elastic recoil of electromagnetic waves.” arXiv preprint arXiv:1407.1800 (2014).

[4] Levine, Richard C. “False paradoxes of superposition in electric and acoustic waves.” American Journal of Physics 48.1 (1980): 28-31.

[5] https://space.mit.edu/RADIO/CST_online/mergedProjects/3D/special_postpr/special_postpr_farfield_array.htm

[6] Rowland, David R. “Comment on “What happens to energy and momentum when two oppositely-moving wave pulses overlap?,” by N. Gauthier [Am. J. Phys. 71 (8), 787–790 (2003)].” American journal of physics 72.11 (2004): 1425-1429.

[7] McDonald, Kirk T. “Does Destructive Interference Destroy Energy?.” (2014).

[8] https://scholar.harvard.edu/files/david-morin/files/waves_electromagnetic.pdf