Preprint
Article

The Wave Equation for the Moving Source and the Moving Receiver

Altmetrics

Downloads

107

Views

69

Comments

0

This version is not peer-reviewed

Submitted:

16 March 2024

Posted:

18 March 2024

You are already at the latest version

Alerts
Abstract
The ordinary 3D wave equation for nondissipative, homogeneous, isotropic media admits solutions where the point sources are permitted to move, but it does not admit solutions where the receiver is allowed to move. To overcome this limitation, the new wave equation is derived, permitting both the receiver and the source to move. This new wave equation is a generalization of the standard wave equation, and it reduces to the standard wave equation when the receiver is at rest. To derive this new wave equation, we first mathematically define a diverging spherical wave caused by the stationary point source. From this purely mathematical definition, the wave equation for the stationary source and the moving receiver is derived together with the corresponding free-space Green's function. Utilizing the derived Green's function, it was shown that this new wave equation also permits the solutions where both the receiver and the source are permitted to move. In conclusion, this paper demonstrates that, instead of an ordinary wave equation, the wave equation for the moving source and the moving receiver governs the waves emitted by the moving point sources and received by the moving receivers.
Keywords: 
Subject: Physical Sciences  -   Mathematical Physics

1. Introduction

One of the most important partial differential equations in physics is the ordinary three-dimensional wave equation [1,2]. It was first derived by Leonhard Euler in the 18th century to describe the propagation of three-dimensional acoustic waves [3,4]. Although originally conceived to describe the propagation of acoustic waves without losses in homogeneous isotropic media, it was later found that under the same conditions,1 the ordinary wave equation describes a variety of wave phenomena such as propagation of the waves in solids [5], the propagation of mechanical stress waves [6], the propagation of seismic waves [7], and as shown by Maxwell in the 19th century [8], the ordinary wave equation describes the propagation of electromagnetic waves.
Regardless of whether the wave equation describes a mechanical, acoustic, or electromagnetic wave, the ordinary three-dimensional wave equation reads:
2 u 1 c 2 2 u t 2 = 0
where u is the displacement, c is the velocity of the wave propagation, t is the time, and 2 is three dimensional Laplace operator. In general, the wave equation (1) describes the propagation of a displacement, or a disturbance, propagating from the source to the receiver with finite velocity c. In the case of an acoustic wave, the function u represents the pressure or the velocity potential [9,10,11]. For an electromagnetic wave, the function u can be any Cartesian component of an electric or magnetic field [12,13].
A particularly important class of solutions to the ordinary wave equation are those that involve moving point sources. In the field of electrodynamics, the Liénard-Wiechert scalar and vector potentials are one example of such solutions [14,15]. These time retarded potentials are widely used for advanced electromagnetic calculations such as synchrotron, undulator, and radar radiation [16,17,18]. For example, the scalar Liénard-Wiechert potential takes the following mathematical form:
u = 1 4 π ϵ q R ( 1 n s · β s )
where R = r r s ( t r ) , r is the position vector of a stationary receiver, r s is the position vector of the moving source, and t r is the retarded time that satisfies the following causal relationship:
t = t r + r r s ( t r ) c
Furthermore, the vectors n s and β s are given by:
n s = r r s ( t r ) r r s ( t r )
β s = 1 c d d t r r s ( t r )
Solutions to the ordinary wave equation for moving sources are also found in the field of acoustics. For example, Schmidt and Kuperman derived the solution to the ordinary wave equation for the harmonic point source moving with constant velocity along a straight line [19]. Obrezanova and Rabinovich derived a two-dimensional solution for harmonic point source moving along arbitrary path [20]. Lim and Ozard derived the solution caused by harmonic acoustic source moving along an arbitrary trajectory in 3D [21] that has a mathematical form similar to Liénard-Wiechert scalar potential.
Modifying any of the point source solutions mentioned earlier to accommodate a moving receiver results in a violation of the ordinary wave equation. For example, if the receiver starts to move, then its position indicated by the vector r becomes the function of time, i.e. r = r ( t ) . As it will be demonstrated later in this paper, if one substitutes r by r ( t ) in the equation (2), then u is no longer the solution to the ordinary wave equation. The same applies to the acoustic point source solutions mentioned in the references [19,20,21].
Thus, the ordinary wave equation has the ability to describe waves caused by moving sources, but it cannot describe the waves seen by the moving receiver. It is shown in this paper that the wave equation that admits the solutions where both the receiver and the source are permitted to move takes the following mathematical form:
2 u 1 c 2 t v · 2 u = 0
where v is the velocity vector of the moving receiver, and where the operator t v · 2 represents the operator t v · applied twice. If the receiver is stationary, then v = 0 , and the equation (6) reduces to an ordinary wave equation.
The proof that for the moving receivers the wave equation (6) is valid instead of the ordinary wave equation is organized as follows: in the Section 2 we first mathematically define a wave function caused by stationary point source, and in the Section 3 we use this wave function to derive the wave equation (6) for the specific case of the moving receiver and the stationary source. Then in the Section 4 Green’s function for the equation (6) is derived. Finally, in the Section 5, by using the derived Green’s function, the solutions to the wave equation (6) are found where both the receiver and the source are permitted to move.
In conclusion, it is shown that instead of an ordinary wave equation, the equation (6) is the governing wave equation when both the receiver and the source are permitted to move. If the receiver becomes stationary, then the wave equation (6) reduces to an ordinary wave equation. This means that the ordinary wave equation can be considered a particular case of the wave equation for the moving receiver and the moving source. These findings are summarized in Table 1, where the wave equation for the moving source and the moving receiver and its point-source solution is compared to the ordinary wave equation and respective point-source solution.

2. Free Space Wave Function for the Stationary Source and Moving Receiver

In this section, we mathematically describe the wave caused by a stationary point source. The wave is simply a disturbance emitted by the source at retarded time t r that travels towards the receiver with finite velocity c. Because the disturbance travels with finite velocity, this disturbance arrives at the moving receiver at the time t > t r . This causal relationship can be represented by the following equation:
t = t r + R ( t ) c
where we declare the time t as an independent variable.
In the expression (7), the function R ( t ) represents the distance between the location of the source at the time of emission t r and the location of the moving receiver at the time t when the wave emitted at the time t r intercepts the receiver. This distance R ( t ) can be written as:
R ( t ) = r ( t ) r s
where the vector r ( t ) = { x ( t ) , y ( t ) , z ( t ) } is the position vector of the moving receiver at the time of wave reception, and the vector r s = { x s , y s , z s } is the position vector of a stationary source.
A simple way to define the wave function u that represents the disturbance traveling in time with velocity c is by declaring:
u = F t R ( t ) c
where F is any function twice differentiable. By substituting the equation (7) into the equation (9) we can equivalently write:
u = F t r
The equations (9) and (10) represent the disturbance that has the value F ( t r ) at the time of emission t r . This disturbance travels towards the receiver with velocity c, and it intercepts the moving receiver at the time t. As indicated in Figure 1, at the time of interception, the disturbance still has the value F ( t r ) .
From the equations (9) and (10), it is also clear that the wave function u represents the spherical wave spreading equally in all directions away from the source. If something emits spherical waves from the position r s , this also means that there is a point source located at the position r s .
Furthermore, the equations (9) and (10) imply that the amplitude of the function u does not change as the wave travels away from the source. In nature, the point sources usually generate diverging spherical waves whose amplitude is inversely proportional to the distance from the source. For example, the acoustic point source generates the diverging spherical waves [22,23], and the point source of electric potential waves and magnetic vector potential waves generates the diverging spherical waves [24,25]. Therefore, if the function u is to represent the wave generated by a physical point source, the wave function u must take the following mathematical form:
u = F t r R ( t )
If the constant is added to the argument of the function F, the function u still satisfies the wave equation. Thus, for the sake of completeness, and for the reasons that will become apparent in the Section 4 of this paper, we rewrite the wave function u as follows:
u = F a t r R ( t )
where a is any real constant. As shown in the Section 4, this particular mathematical form of the wave function allows us to derive Green’s function for the wave equation for the moving receiver in a relatively straightforward manner.

3. The Wave Equation for the Moving Receiver and the Stationary Source

In this section, we derive the wave equation for the stationary point source and the moving receiver. Here, the position of the receiver is considered the function of time t, i.e. r = r ( t ) , and the position vector r s of the stationary source is not considered a function of time:
r s = { x s , y s , z s } = c o n s t
r = r ( t ) = x ( t ) , y ( t ) , z ( t )
The distance R ( t ) between the source at the time of wave emission and the receiver at the time of wave interception is the function of time:
R ( t ) = r ( t ) r s
The derivative of the wave function u given by the equation (11) with respect to variable x i { x ( t ) , y ( t ) , z ( t ) } can be calculated using the chain rule:
u x i = F R t r x i + F x i 1 R
where R is a shorthand notation for R ( t ) , and where F = F ( a + t r ) and F = F ( a + t r ) . Note that because of the equation (7), the derivative of the retarded time with respect to variable x i is non-zero. This follows from the equation (7) where we declared the variable t as an independent variable, therefore, the retarded time t r is the function of the variable t, constant vector r s , and the vector function r ( t ) .
By differentiating the equation (16) with respect to variable x i , the second order derivative of the wave function u is obtained as:
d u 2 d x i 2 = F R t r x i 2 F R 2 t r x i R x i + F R 2 t r x i 2 + F t r x i x i 1 R + F 2 x i 2 1 R
where F = F ( a t r ) . The Laplacian of the function u can now be written as the sum:
2 u = i u 2 x i 2 = F R i t r x i 2 F R 2 i t r x i R x i + F R i 2 t r x i 2 + F i t r x i x i 1 R + F i 2 x i 2 1 R
By differentiating the equation (7) with respect to variable x i it is obtained:
t r x i = 1 c R x i
The derivative R x i can be obtained by differentiating the equation (15) with respect to variable x i as:
R x i = n i
where n i is ith Cartesian component of the unit vector n given as:
n = r ( t ) r s r ( t ) r s
Clearly, because the vector n is an unit vector, it can be written:
n · n = i n i 2 = 1
Using the equations (19), (20) and the equation (22), one can calculate the sums found in the equation (18) as:
i t r x i 2 = i 1 c 2 R x i 2 = 1 c 2 i n i 2 = 1 c 2
i t r x i R x i = i 1 c R x i n i = 1 c i n i 2 = 1 c
i 2 t r x i 2 = i x i 1 c R x i = 1 c i n i x i = 1 c · n
i t r x i x i 1 R = i 1 c R x i 1 R 2 R x i = 1 c 1 R 2 i n i 2 = 1 c 1 R 2
i 2 x i 2 1 R = 2 1 R = 4 π δ r ( t ) r s
where the symbol δ represents the 3D Dirac’s delta function. The divergence of the unit vector n can be calculated as:
· n = i x i x i x s i R = 3 R i n i 2 R = 2 R
where x s i are Cartesian components of the vector r s , i.e. x s i { x s , y s , z s } . Substituting the equations (23) - (28) into the equation (18) yields:
2 u = 1 c 2 F R + 1 c F R 2 2 c F R 2 + 1 c F R 2 4 π F δ r ( t ) r s
After cancellation of the appropriate terms on the right-hand side of the equation (29), the equation (29) simplifies to:
2 u = 1 c 2 F R 4 π F δ r ( t ) r s
We now introduce the operator t v · , where v is the velocity of the receiver:
v = d r ( t ) d t
By applying the operator t v · to the function u it is obtained:
t v · u = F R t r t F R 2 R t F R i t r x i v i + F R 2 i R x i v i
where v i are Cartesian components of the vector v . Substituting the equations (19) and (20) into the equation (32) it is obtained:
t v · u = F R t r t F R 2 R t + 1 c F R i n i v i + F R 2 i n i v i
By differentiating the equation (15) the time derivative R t is obtained as:
R t = n · v = i n i v i
Substituting the equation (34) into the equation (33) yields:
t v · u = F R t r t F R 2 R t + 1 c F R R t + F R 2 R t
The second and the fourth right-hand side terms of the equation above cancel, thus, we may write:
t v · u = F R 1 c R t + t r t
By differentiating the equation (7) with respect to time one obtains:
1 c R t + t r t = 1
Substituting the equation (37) into the equation (36) yields:
t v · u = F R
Next we will apply the operator t v · to the equation (38) to obtain:
t v · 2 u = t v · F R = = F R d t r d t F R 2 d R d t F R i t r x i v i + F R 2 i R x i v i
where the operator t v · 2 simply means the operator t v · applied twice. Substituting the equations (19) and (20) into the equation (39) it is obtained:
t v · 2 u = F R t r t F R 2 R t + 1 c F R i n i v i + F R 2 i n i v i
Furthermore, substituting the equation (34) into the equation (40) yields:
t v · 2 u = F R t r t F R 2 R t + 1 c F R R t + F R 2 R t
The second and the fourth right-hand side terms cancel, therefore:
t v · 2 u = F R 1 c R t + t r t
By substituting the equation (37) into the equation (42) it is obtained that:
t v · 2 u = F R
Subtracting the equation (43) divided by c 2 from the equation (30) yields the wave equation for the stationary source and the moving receiver:
2 u 1 c 2 t v · 2 u = 4 π F δ r ( t ) r s
where F is any function such that F = F ( a + t r ) . This means that the function F can even be a constant, which has applications in electrodynamics [24]. It is also interesting to note that if the receiver becomes stationary, then r ( t ) = r = c o n s t and the velocity of the receiver v = 0 . In that case, the wave equation for the stationary source and the moving receiver becomes the wave equation (1).

4. The Green’s Function for the Wave Equation for the Moving Receiver and the Stationary Source

The free space Green’s function for the three-dimensional version of the wave equation is typically derived by the application of Fourier or Laplace transform [26,27], which is somewhat mathematically involving. For the wave equation given by the equation (1), the Green’s function G satisfies the following inhomogeneous wave equation:
2 G 1 c 2 2 G t 2 = δ ( r r ) δ ( t t )
Using Fourier’s or Laplace transform, one can solve the equation (45) for the function G to obtain Green’s function for the wave equation in an unbounded domain [28,29]:
G = G ( r , r ; t , t ) = δ t t r r c 4 π r r
which represents the effect of an impulse propagating away from the source located at the position r .
We are now interested in finding Green’s function for the wave equation for the moving receiver and the stationary source expressed by the equation (44). Therefore, we need to find the function G that satisfies the following inhomogeneous wave equation:
2 G 1 c 2 t v · 2 G = δ r ( t ) r δ ( t t )
Instead of using Fourier or Laplace transform to find the function G, we here take the simpler approach that takes advantage of the right-hand side of the equation (44). If in the equation (8) we replace the constant vector r s by the constant vector r the function u becomes:
u = F ( a + t r ) r ( t ) r = F a + t r ( t ) r c r ( t ) r
Because r is a constant vector, the function u still satisfies the wave equation for the moving observer and the stationary source given by the equation (44). Now by letting a = t , where t is an independent real variable, the function u becomes:
u = F t t r ( t ) r c r ( t ) r
Note that because t is an independent variable, the function u still satisfies the equation (44). The function F can be any differentiable function, thus, we choose F to be the Dirac’s delta function divided by 4 π :
F t t r ( t ) r c = δ t t r ( t ) r c 4 π
With this choice, the function u becomes:
u = δ t t r ( t ) r c 4 π r ( t ) r
Substituting the function u into the equation (44) yields:
2 u 1 c 2 t v · 2 u = δ r ( t ) r δ t t r ( t ) r c
The term δ r ( t ) r is non-zero only when r ( t ) r . In that case the Dirac delta function δ t t r ( t ) r c becomes:
δ t t r ( t ) r c r ( t ) r δ ( t t )
Substituting the equation (53) into the equation (52) yields:
2 u 1 c 2 t v · 2 u = δ r ( t ) r δ t t
Thus, it is evident that the function u as defined by the equation (51) is indeed Green’s function for the wave equation for the moving receiver and the stationary source:
G = δ t t r ( t ) r c 4 π r ( t ) r
Note the difference between the Green’s functions given by the equations (46) and (55): with the equation (55), the receiver is allowed to move because its position is described by the position vector r ( t ) . With the equation (46), the receiver is not allowed to move because its position is described by the constant position vector r .

5. The Wave Equation for the Moving Source and the Moving Receiver

Thus far, it was assumed that the wave equation for the moving source and the stationary receiver given by the equation (44) admits only the solutions in which the receiver is allowed to move, but the source is stationary. In this section, however, we are attempting to find the solution to the wave equation (44) in which both the source and the receiver are allowed to move. If such solution exists, it satisfies the following inhomogeneous wave equation:
2 u 1 c 2 t v · 2 u = 4 π F ( a + t r ) δ r ( t ) r s ( t r )
Furthermore, if the function u is the solution to the equation (56) then the distance R between the location of the source at the time t r when the wave was emitted and the location of the moving receiver at the time t when this wave arrives at the receiver is the function of both the retarded time t r and the time t:
R = R ( t , t r ) = r ( t ) r s ( t r )
Substituting the equation (57) into causal relationship given by the equation (7) yields:
t = t r + r ( t ) r s ( t r ) c
Note the difference between the equation (44) and the equation (56). In the equation (44), the motion of the source is not allowed while the receiver can move along an arbitrary trajectory. In contrast, the equation (56) allows both the source and the receiver to move along arbitrary trajectory.
Before we start the mathematical procedure of finding such solution, note that the right-hand side term of the equation (56) is non-zero only when r ( t ) = r ( t r ) . This can only happen if the source and the receiver occupy the same point in 3D space, that is, if the source and the receiver collide. In that case, the retarded time t r is equal to present time t, and we can write:
2 u 1 c 2 t v · 2 u = 4 π F ( a + t ) δ r ( t ) r s ( t )
The Green’s function for the wave equation for the moving receiver and the stationary source allows us to write the solution to the equation (59) as:
u = 4 π d t d V F ( a + t ) δ r ( t ) r s ( t ) G
Substituting the equation (55) into the equation (60) yields:
u = 4 π d t d V F ( a + t ) δ r ( t ) r s ( t ) δ t t r ( t ) r c 4 π r ( t ) r
The sifting property of Dirac’s delta function δ r ( t ) r s ( t ) picks out point r = r ( t ) = r s ( t ) , thus, the equation (61) becomes:
u = 4 π d t F ( a + t ) δ t t r ( t ) r s ( t ) c 4 π r ( t ) r s ( t )
Using the following identity involving the Dirac’s delta function [30,31]:
δ ( H ( x ) ) = δ ( x x 0 ) d H ( x ) d x x = x 0
where H ( x ) is a differentiable function, and x 0 is the solution to the equation H ( x 0 ) = 0 , one can write the Dirac’s delta function under the integral on the right-hand side of the equation (62) as:
δ t t r ( t ) r s ( t ) c = δ ( t t 0 ) d d t t t r ( t ) r s ( t ) c t = t 0
where t 0 is the solution to the following equation:
t t 0 r ( t ) r s ( t 0 ) c = 0
By comparing the equations (58) and (65) it follows that the time t 0 is:
t 0 = t r
Differentiating the expression in denominator of the equation (64) while keeping in mind that t is an independent variable yields:
d d t t t r ( t ) r s ( t ) c = 1 + 1 c r ( t ) r s ( t ) r ( t ) r s ( t ) · d r s ( t ) d t
Substituting the equation (67) into the equation (64) yields:
δ t t r ( t ) r s ( t ) c = δ ( t t 0 ) 1 + 1 c r ( t ) r s ( t ) r ( t ) r s ( t ) · d r s ( t ) d t t = t 0
Using the equation (66), i.e. t 0 = t r , yields:
δ t t r ( t ) r s ( t ) c = δ ( t t r ) 1 + 1 c r ( t ) r s ( t r ) r ( t ) r s ( t r ) · d r s ( t r ) d t r
If we assume that d r s ( t r ) d t r < c the equation (70) can be written as:
δ t t r ( t ) r s ( t ) c = δ ( t t r ) 1 1 c r ( t ) r s ( t r ) r ( t ) r s ( t r ) · d r s ( t r ) d t r
The equation (70) can be written in more concise form as:
δ t t r ( t ) r s ( t ) c = δ ( t t r ) 1 β s · n
where the vectors β s and n are given as:
β s = β s ( t r ) = 1 c d r s ( t r ) d t r
n = n ( t , t r ) = r ( t ) r s ( t r ) r ( t ) r s ( t r )
Furthermore, substituting the equation (71) into the equation (62) yields:
u = 4 π d t F ( a + t ) δ t t r 4 π r ( t ) r s ( t ) 1 β s ( t r ) · n ( t , t r )
Using the sifting property of the Dirac delta function, we can rid ourselves of the integral in the equation (74) to obtain:
u = F ( a + t r ) R 1 β s · n
where R = r ( t ) r s ( t r ) , where the retarded time t r satisfies the equation (58), and where the vector n = n ( t , t r ) is the function of both time t and the retarded time t r .
Put in other words, the function u is the wave function when both the source and the receiver can move along arbitrary trajectories, and the function u satisfies the wave equation (56). This means that the equation (44) derived in the Section 3 of this paper is not only the wave equation for the moving receiver and the stationary source, but it is also the wave equation for the moving receiver and the moving source.
Note that the equation (75) does not satisfy the ordinary wave equation given by the equation (1). However, if the receiver becomes stationary, then r ( t ) = r = c o n s t , and consequently, the velocity of the receiver v = 0 . In such a scenario, the wave equation for the moving source and the moving receiver (56) simplifies to the ordinary wave equation, and the solution (75) also simplifies to the solution to the ordinary wave equation. Thus, it is evident that the ordinary wave equation is a special case of the wave equation (56).
From the point of view of electrodynamics, it is worth noting that if the function F is a constant, such that F ( a + t r ) = q / 4 π ϵ , and if the velocity of the receiver v = 0 , then the equation (75) becomes the scalar Liénard-Wiechert potential caused by the moving point charge. This potential is also the solution to the ordinary wave equation because the receiver is stationary.

6. Summary and Discussion

In Table 1, the results of the investigations conducted in this paper are concisely summarized. The table is divided into two parts: the first deals with the ordinary wave equation, and the second deals with the wave equation for the moving source and the moving receiver. The moving point-source solution to the ordinary wave equation and the corresponding parameters are obtained by setting v = 0 (and consequently r ( t ) = r = c o n s t ) in the equation (59). As indicated in the Table 1, the ordinary wave equation cannot have solutions where the receiver is in motion.
Although the point source solutions to ordinary wave equation are obtained from the equation (59), these solutions can be verified in a different manner by using the differentiation techniques described in the Appendix A. However, verifying the equations given in the Table 1 using these differentiation techniques requires a significant amount of mathematical manipulation, which would significantly and unnecessarily lengthen this paper. Alternatively, advanced symbolic computation software such as Mathematica, together with the techniques given in the Appendix A, can be employed to verify the equations in Table 1 in a different manner than the one presented in this paper.
The second part of the Table 1 presents the wave equation for the moving source and the moving receiver and its general solution for the point sources. The equations given in the second part of this table represent the summary of the Section 5 of this paper. As indicated in the second part of the Table 1, the wave equation (59) is more general than an ordinary wave equation because it admits the solutions for the moving receiver. In fact, it can be concluded that the ordinary wave equation is a particular case of the wave equation for the moving source and the moving receiver when v = 0 .
Table 1. Comparison of ordinary wave equation and the wave equation for the moving source and the moving receiver. The ordinary wave equation does not allow the motion of the receiver, while the wave equation for the moving source and the moving receiver allows the motion of both the source and the receiver.
Table 1. Comparison of ordinary wave equation and the wave equation for the moving source and the moving receiver. The ordinary wave equation does not allow the motion of the receiver, while the wave equation for the moving source and the moving receiver allows the motion of both the source and the receiver.
Name of the equation Expression(s) Source motion allowed Receiver motion allowed
Ordinary wave equation 2 u 1 c 2 2 u t 2 = 4 π F ( a + t r ) δ r r s ( t r )
Point source solution:
u = F ( a + t r ) R ( t r ) 1 β s ( t r ) · n ( t r ) X
Parameters:
R ( t r ) = r r s ( t r )     β s ( t r ) = 1 c d r s ( t r ) d t r     n s ( t r ) = r r s ( t r ) r r s ( t r )
t = t r + R ( t r ) c
F - any function differentiable twice
a - any constant
the wave equation for the moving source and the moving receiver 2 u 1 c 2 t v · 2 u = 4 π F ( a + t r ) δ r ( t ) r s ( t r )
Point source solution:
u = F ( a + t r ) R ( t , t r ) 1 β s ( t r ) · n ( t , t r )
Parameters:
R ( t , t r ) = r ( t ) r s ( t r ) β s ( t r ) = 1 c d r s ( t r ) d t r n s ( t , t r ) = r ( t ) r s ( t r ) r ( t ) r s ( t r )
v = d r ( t ) d t - the velocity of the receiver
t = t r + R ( t , t r ) c
F - any function differentiable twice
a - any constant

7. Conclusions

There are several principal results of this study, the foremost being a derivation of a three-dimensional wave equation that allows the motion of both the source and the receiver. Secondly, this study has demonstrated the inadequacy of the ordinary wave equation to admit the solutions where the receiver is allowed to move. Moreover, it was demonstrated that the ordinary wave equation represents a particular case of the wave equation derived in this paper. Therefore, it is concluded that the wave equation for the moving source and the moving receiver is a more general equation than the ordinary wave equation.
The third, but not less important result, is that the general point source solutions are derived for both the ordinary wave equation and for the wave equation for the moving source and the moving receiver. These solutions can be adapted for various physical problems involving wave propagation without losses in homogeneous isotropic media. Additionally, the free space Green’s function for the wave equation for the moving source and the moving receiver has also been derived.
Lastly, the wave equation for the moving source and the moving receiver has potential application in physics. In particular, the limitations of Liénard-Wiechert potentials in electrodynamics, namely their inability to represent moving receivers, can be resolved by altering these potentials. This modification could potentially aid with some modern challenges in the field of electrodynamics.

CRediT authorship contribution statement

Hrvoje Dodig: Conceptualization, Investigation, Methodology, Visualization, Writing - original draft.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data availability

No data was used for the research described in the article.

Acknowledgments

The author is thankful to Šimun Anđelinović, the former rector of the University of Split, Croatia for moral and otherwise support and for his insightful remarks regarding the origin of Newton’s physics.

Appendix A. Differentiation of Causal Expressions

In this appendix, we derive some mathematical identities that can be used to verify the equations given in the Table 1 in a different manner than the one presented in this paper. By differentiating the equation (58) with respect to variable x, where x = x ( t ) , one obtains:
0 = t r x + 1 c x ( t ) x s ( t r ) R ( t , t r ) 1 d x s ( t r ) d t r t r x 1 c y ( t ) y s ( t r ) R ( t , t r ) d y s ( t r ) d t r t r x 1 c z ( t ) z s ( t r ) R ( t , t r ) d z s ( t r ) d t r t r x
The equation above can be written in a simpler notation:
0 = t r x + n x c β s · n t r x
where n is given by the equation (), the vector β s is given by the equation (72), and n x is Cartesian x component of the vector n . From the equation above, it follows that:
t r x = 1 c n x 1 β s · n = 1 c n x g
where the function g is:
g = 1 β s · n
We can now generalize the result provided by the equation (A3) as:
d t r d x i = 1 c n i g
where n i is the Cartesian component of the vector n and x i { x ( t ) , y ( t ) , z ( t ) } . Using the expression (A5) one can now calculate the gradient t r as:
t r = 1 c n g
Furthermore, by differentiating the equation (58) with respect to time t one obtains:
1 = t r t + 1 c x ( t ) x s ( t r ) R ( t , t r ) d x ( t ) d t d x s ( t r ) d t r t r t + 1 c y ( t ) y s ( t r ) R ( t , t r ) d y ( t ) d t d y s ( t r ) d t r t r t + 1 c z ( t ) z s ( t r ) R ( t , t r ) d z ( t ) d t d z s ( t r ) d t r t r t
The equation above can be written in simpler notation using the vectors n and β s as:
1 = d t r d t + 1 c v · n β s · n d t r d t
where the velocity v = d r ( t ) d t is the velocity of the receiver. From the equation above, it follows that the derivative t r t is:
t r t = 1 β · n 1 β s · n = f g
where the vector β is:
β = 1 c d r ( t ) d t
and where the function f is:
f = 1 β · n
One can now use the equations (A6) and (A9) to calculate the gradient, divergence, Laplacian, and the second order time derivative of the function u given by the equation (75) to verify the equations given in the second part of the Table 1. To verify the equations in the first part of the Table 1, we set v = 0 , R = r r s ( t r ) , β = 0 . Note that verifying the equations in the Table 1 using expressions (A6) and (A9) requires a significant amount of mathematical manipulation, however, advanced symbolic computation software such as Mathematica can significantly aid that task.

References

  1. J. Kirkwood, Mathematical Physics with Partial Differential Equations, Academic Press, Amsterdam, 2013.
  2. L. Evans, Partial Differential Equations, Vol. 19., American Mathematical Society, Providence, USA, 1998.
  3. Euler, L. De la propagation du son. Mém. Acad. Sci. Berlin 1766, 15, 185–209. [Google Scholar]
  4. M. Kline, Mathematical Thought from Ancient to Modern Times, Oxford University Press, Oxford, UK, 1972.
  5. D. Royer, E. Dieulesaint, Elastic Waves in Solids I, Springer, Berlin, 2000.
  6. L. Wang, Foundations of Stress Waves, Elsevier, Amsterdam, 2007.
  7. E. Krebes, Seismic Wave Theory, Cambridge University Press, Cambridge, UK, 2019.
  8. Maxwell, J. A dynamical theory of the electromagnetic field. Phil. Trans. R. Soc. 1865, 155, 459–512. [Google Scholar]
  9. L. Beranek, Acoustics, American Institute of Physics, Inc., Woodbury, New York, 1986.
  10. J. Ginsberg, Acoustics-A Textbook for Engineers and Physicists, Springer International Publishing AG, Cham, Switzerland, 2018.
  11. A. Pierce, Acoustics - An Introduction to Its Physical Principles and Applications, 3rd Edition, Springer Nature Switzerland AG, Cham, Switzerland, 2019.
  12. C. Holt, Introduction to Electromagnetic fields and Waves, John Wiley & Sons, Inc., New York, 1963.
  13. L. Felsen, N. Marcuvitz, Radiation and Scattering of Waves, John Wiley & Sons, Inc., Hoboken, New Jersey, 1994.
  14. Wiechert, E. Elektrodynamische elementargesetze. Ann. Phys. 1901, 309, 667–689. [Google Scholar] [CrossRef]
  15. Liénard, A. Champ électrique et magnétique produit par une charge électrique concentrée en un point et animée d’un mouvement quelconque. L’Éclairage Électrque 1898, 16, 5–14. [Google Scholar]
  16. Arakelyan, S.; Grigoryan, G.; Grigoryan, R. Lienard-Wiechert potentials and synchrotron radiation from a relativistic spinning particle in pseudoclassical theory. Phys. At. Nucl. 2000, 63, 2115–2122. [Google Scholar] [CrossRef]
  17. Steinbach, A. Electromagnetic theory of range-doppler imaging in laser radar. i: Scattering theory. J. Opt. Soc. Am. A 1991, 8, 1287–1295. [Google Scholar] [CrossRef]
  18. Salehi, E.; Hosaka, M.; Katoh, M. Time structure of undulator radiation. J. Adv. Simulat. Sci. Eng. 2023, 10, 164–171. [Google Scholar] [CrossRef]
  19. Schmidt, H.; Kuperman, W. Spectral and modal representations of the doppler-shifted field in ocean waveguides. J. Acoust. Soc. Am. 1994, 96, 386–395. [Google Scholar] [CrossRef]
  20. Obrezanova, O.; Rabinovich, V. Acoustic field generated by moving sources in stratified waveguides. Wave Motion 1998, 27, 155–167. [Google Scholar] [CrossRef]
  21. Lim, P.; Ozard, J. On the underwater acoustic field of a moving point source. i. range-independent environment. J. Acoust. Soc. Am. 1994, 5, 131–137. [Google Scholar] [CrossRef]
  22. L. Kinsler, A. Frey, A. Coppens, J. Sanders, Fundamentals of acoustics, John Wiley & Sons, New York, 2000.
  23. P. Filippi, D. Habault, J. Lefebvre, A. Bergassoli, Acoustics, Academic Press, London, 1999.
  24. Dodig, H. Direct derivation of Lienard–Wiechert potentials, Maxwell’s equations and Lorentz force from Coulomb’s law. Mathematics 2021, 9, 237. [Google Scholar] [CrossRef]
  25. J. Jackson, Classical Electrodynamics, 3rd Edition, John Wiley & Sons, London, 1999.
  26. D. Duffy, Green’s Functions with Applications, Chapman & Hall/CRC, Boca Raton, 2001.
  27. K. Watanabe, Integral Transform Techniques for Green’s Function, Springer International Publishing Switzerland, Cham, Switzerland, 2014.
  28. P. Dennery, A. Krzywicki, Mathematics for Physicists, Dover Publications Inc., Mineola, New York, 1996.
  29. Hoop, A.T.D. Fields and waves excited by impulsive point sources in motion—the general 3d time-domain doppler effect. Wave Motion 2005, 43, 116–122. [Google Scholar] [CrossRef]
  30. G. Arfken, H. Weber, Mathematical Methods for Physicists, 6th Edition, Elsevier Academic Press, Amsterdam, 2006.
  31. B. Kusse, E. Westwig, Mathematical Physics, 2nd Edition, WILEY-VCH Verlag Gmbh & KGaA, Winheim, Germany, 2006.
1
These conditions assume the homogeneity and isotropy of media and the wave propagation without losses.
Figure 1. The stationary point source emits the disturbance that has the value F ( t r ) at the time t r . This disturbance travels with velocity c and spreads equally in all directions. At the moment t = t r + R ( t ) / c , the disturbance carrying the value F ( t r ) intercepts the receiver moving with the velocity v .
Figure 1. The stationary point source emits the disturbance that has the value F ( t r ) at the time t r . This disturbance travels with velocity c and spreads equally in all directions. At the moment t = t r + R ( t ) / c , the disturbance carrying the value F ( t r ) intercepts the receiver moving with the velocity v .
Preprints 101532 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated