Spherical One-Way Wave Equation

: The coordinate-free one-way wave equation is transferred in spherical coordinates. Therefore it is necessary to achieve consistency between gradient , divergence and Laplace operators and to establish, beside the conventional radial Nabla operator ∂ Φ / ∂ r , a new variant ∂ r Φ / r ∂ r . The two Nabla operator variants differ in the near ﬁeld term Φ / r whereas in the far ﬁeld r (cid:29) 0 there is asymptotic approximation. Surprisingly, the more complicated gradient ∂ r Φ / r ∂ r results in unexpected simpliﬁcations for – and only for – spherical waves with the 1/ r amplitude decrease. Thus the calculation always remains elementary without the wattless imaginary near ﬁelds, and the spherical Bessel functions are obsolete.


Introduction
Seismics, sonar, sodar, room and machine acoustics, ultrasound diagnosis and tomography depend on the calculation of sound paths in solid, liquid or gaseous media with heterogenous physical properties and complex geometry.Calculations are based on the well-known equations of motion listed by Augustin-Louis Cauchy 200 years ago: In a continuum with stationary cartesian coordinates x x x = {x e x e x e x , y e y e y e y , z e z e z e z } and modulus of elasticity E [Pa] a vectorial, elastic deflection s s s = s s s(x x x, t) [m] causes the stress tensor T T T = E∇ ∇ ∇s s s [Pa].The local equilibrium of the tension force div T T T = div E∇ ∇ ∇s s s = E∆s s s [N/m 3 ] due to stress tensor T T T and the inertial force ρ s s s [N/m 3 ] caused by the local acceleration s s s = ∂ 2 s s s/∂t 2 [m/s 2 ] can be written as The partial differential equation (PDE) of the 2nd order (2) delivers two mutually independent solutions.From the quadratic velocity term c 2 = (+c) 2 = (−c) 2 can be seen that there are two waves travelling in opposite directions +c and −c, hence results the designation "Two-way wave equation".Regardless of this ambiguity, irregular phantom effects occur in numerical seismic FE or FD wave calculations.
To eleminate these unwanted effects a great number of auxiliary equations have been developed, but no specific approach has been able to prevail.Because of the economic importance of the "One-Way"/"Two-Way" problem there also exist numerous patent applications beside the scientific literature [1][2][3].
Due to the above mentioned difficulties the conventionally used equilibrium of forces was -initially hypothetically -replaced by an impulse equilibrium, i.e., more precisely by an impulse flow equilibrium [4,5].Force and impulse are related and differ in one level of differentiation.Therefore, the unit of the impulse Huygens [Hy = mkg/s = mass multiplied by velocity] was used.For a particle velocity ṡ s s = ∂s s s/∂t [m/s] the local specific impulse ρ ṡ s s [Hy/m 3 ] follows.If velocity ṡ s s propagates with the vectorial wave velocity c c c, the dyadic product ρ c c c ṡ s s results in a kinetic impulse flow of the dimension [Hy/sm 2 ], i.e., Huygens [Hy] per time and per area.On the other hand, an elastic deformation ∇ ∇ ∇s s s induces the potential impulse flow T T T = E∇ ∇ ∇s s s [Hy/sm 2 ].Wheras the Cauchy force equilibrium is applied to an infinitesimal volume element dV = dx dy dz, the equilibrium of kinetic and potential impulse flow in each field point x x x results in This tensor equation being scalarly multiplied by the vectorial wave velocity c c c leads together with the known relationship E = ρ c 2 to a vectorial partial differential equation (PDE) of the 1st order respectively vectorial "One-Way wave equation" Representing the equivalent to the force based 2nd order PDE wave Equation (2).Mathematically, the 1st order PDE is much easier than the 2nd order PDE with the higher differentiation level.Also the wave propagation direction is determined.Force and impulse concept are compared below in Table 1.

Two-Way Wave Equation
One-Way Wave Equation

Balance Item
Force Because of their vector character the two wave Equations ( 2) and (4) are independent from the coordinate system and have to be transcribed to the selected coordinate system when used.Such a transcription is a mathematical formality.
While transcription to cartesian coordinates is problem-free, an unexpected conflict occurs for transcription to spherical coordinates.However, establishing the hypothetical radial Nabla variant ∇ r Φ = ∂rΦ/r∂r beside the well-known conventional radial Nabla operator ∇ r Φ = ∂Φ/∂r is a solution to the conflict.Hence, it is a further task to clarify the conflicting transcription of the vectorial one-way wave equation to spherical coordinates and to provide theoretical results that can be experimentally verified later.Therefore, calculation is restricted to spherical wave propagation without the Legendre angular functions, i.e., Φ(r, θ, φ) → Φ(r).Interestingly the rather complicated hypothetical gradient ∂rΦ/r∂r leads to unexpected simplifications for -and only for -spherical waves with 1/r amplitude decrease.Thus, calculations remain elementary and a priori, neither bulky transcendental Bessel functions nor wattless imaginary near fields appear.

Method
According to the task and the restriction to central fields with purely radial functions Φ(r, θ, φ) → Φ(r) and s s s(r, θ, φ) → s s s(r), spherical coordinates r r r = r e r e r e r are sufficient and only the radial operators gradient ∇ r Φ, the divergence div r s s s = ∇ r • s s s and Laplace r Φ are relevant in this context.These Nabla operators are used in all formulas, for the sake of certainty they are cited from [6].In Table 2 all expressions are listed and connected with the ∇ symbol.
Table 2. Acc. to their signature the conventional radial operators gradient, divergence and Laplace lead to different radial Nabla variants ∇ a , ∇ b and ∇ c .Nabla variant ∇ c is taken as hypothesis.

Radial Operator
Convention Signature Radial Nabla Variants Obviously the following three different radial Nabla variants result: Variants ∇ a and ∇ b being multiplied with each other fullfill the Laplace operator condition ∆Φ = ∇ b ∇ a Φ = (∂r 2 /r 2 ∂r) ∂Φ/∂r, but due to ∇ a = ∇ b not the identity condition ∆ = ∇•∇ = ∇ 2 .The calculated square root of Laplace ∆ 1/2 = ∇ c = ∂r/r∂r leads to a 3rd version fullfilling all conditions.With this justification the radial Nabla operator is taken as hypothesis and compared with the conventional Nabla operators below in Table 3.Control of the calculation by aligns with conventional table set ∆ r Φ = Φ + 2Φ /r and can be simplified to: Thereby following derivatives are used: Table 3.Comparison of conventional and hypothetical spherical operators ∇ r Φ, div r s s s and ∆ r Φ.
For gradient and divergence operators the calculated results of the convention and the hypothesis misalign.Only for the Laplace operator both concepts provide the same result ∆ r Φ = Φ + 2Φ /r.

Radial Operator Convention [6] Hypothesis
Gradient: Punctum saliens: For the ∇ r -differentiation of the special function Φ = Ψ/r with the obligatory 1/r amplitude decrease for a spherical wave the 1/r property is retained, even if the ∇ r differentiation is n-times repeated: To demonstrate this advantageous property, Table 4 shows the specific impedance z = p/v = ρ c on the left according to the conventional table rule ∇ r Φ = ∂Φ/∂r and on the right being calculated with the new hypothetical rule ∇ r Φ = ∂rΦ/r∂r.Subsequently the conventional calculation delivers a complex impedance, whereas the calculation with ∇ r Φ = ∂rΦ/r∂r yields a consistent, purely material-dependent real impedance z = ρ c.Table 4. Calculation of particle velocity, pressure and specific impedance z = p/v of a spherical wave with velocity potential G using different Nabla operators.For the conventional Nabla operator ∇ r Φ = ∂Φ/∂r [7] follows a complex impedance, wheras for the hypothetical Nabla operator ∇ r Φ = ∂rΦ/r∂r a purely material-dependent real impedance z = ρ c results.

Spherical Conventional Hypothetical Wave Propagation
Nabla Operator Nabla Operator

Results
An omnidirectionally unlimited three-dimensional continuum, in which solely longitudinal and transverse spatial waves but no surface waves can propagate, is taken as basis.In case of linearity and thanks to the orthogonality there is no mutual influence and the individual wave types can therefore be treated separately.With the labels s s s → s L s L s L and c c c → c L c L c L the longitudinal wave equation follows By definition, the deflections s L s L s L of the longitudinal wave lie in the direction of the wave propagation given by wave velocity c L c L c L .Without loss of generality, with c L c L c L = c L e x e x e x and s L s L s L = s L e x e x e x the plane wave should run in coordinate direction e x e x e x .Thus, the vectorial wave Equation ( 11) is simplified to scalar wave Equation ( 12) and the well-known solution (13) of a planar wave is obtained (with ω [rad/s] = circular frequency, k [rad/m]=wave number, For the spherical wave propagating in the radial direction r e r e r e r , the deflections s also have the same direction and the scalar wave equation follows with the hypothetical Nabla operator ∇ r s L = ∂rs L /r∂r As can be verified by insertion, the solution is a spherical wave with 1/r distance dependency of the amplitude (S L,0 =amplitude factor [m 2 ]): In contrast, the conventional radial gradient ∇ r s L = ∂s L /∂r describes a plane wave travelling in re r e r e r direction as shown in Table 5.This is due to the fact that ∂s L /∂r has the analogue form as the cartesian gradient ∂s L /∂x.
In case of the transverse wave (T), deflection s T s T s T and direction of propagation c T c T c T are perpendicular to each another, i.e., c T c T c T • s T s T s T = 0.In this case, the well-known two-way wave equation can be factorized by means of the identity (a and the vectorial one-way wave equation follows Transversality of the Equation ( 17) is proven by scalar multiplication with c T c T c T :, i.e., Since a spar product with two equal vectors c T c T c T is always zero follows c T c T c T • ṡ s s T T T = 0, i.e., the direction of the wave c T c T c T and the deflection ṡ s s T T T are transversal.It can also be shown that s T s T s T and rot rot rot s T s T s T are perpendicular to each another.Finally, the scalar multiplication of Equation ( 17) with s T s T s T result in i.e., the three vectors s T s T s T , c T c T c T and rot rot rot s T s T s T form an orthogonal tripod.So the vectorial transverse one-way wave Equation (17) can be expressed in scalar format as Above calculations refer to the impulse related concept "one-way wave equation".Finally, a longitudinal spherical wave should be calculated according to the "two-way wave equation" principle.The force related wave Equation ( 2 Bessel functions are transcendent and the solution can only be approximated by series expansion.An exception is the elementary half-integer function J 1/2 (kr), therefore the specific labelling as spherical Bessel function j 0 (kr).
Table 5. Spherical longitudinal wave propagation in an unlimited, elastic continuum without boundary effects according to the spherical one-way wave equation.For reasons of simplicity and definiteness linearity, homogeneity, isotropy and losslessness are assumed.The table compares the wave equations and their solutions using conventional and hypothetical Nabla operators.For the conventional Nabla operator results a plane wave with constant amplitude.For the hypothetical Nabla operator follows a solution with the usual 1/r amplitude decrease of spherical waves.

Spherical Longitudinal Conventional Hypothetical One-Way Wave Equation
Nabla Operator Nabla Operator

Discussion
Following, the specific impedance z = p/v is taken to show numerical differences between the two competing Nabla operator variants.According to Table 4 for the conventional (conv) Nabla operator the complex impedance z is normalized to ρc and transferred into polar form with the amplitude function A [-] and the phase φ [rad]: A conv = kr 1 + (kr) 2 (25) For the hypothetical (hyp) Nabla operator the impedance amplitude function A hyp [-] and the phase φ hyp [rad] have the following values To illustrate the numerical differences between the conventional Nabla operator and the hypothetical Nabla operator variant Figure 1 shows the impedance amplitude functions and Figure 2 the phases depending on source distance r devided by wavelength λ [m]; the value r/λ [-] is taken because it is more easily comprehensible than kr = 2πr/λ.According to both graphs, the conventional and hypothetical Nabla operator concepts substantially differ in the near field, but show asymptotic approximation in the far field.Because in many contexts of spherical wave propagation near field influences are neglected, above given spherical one-way wave equations and the hypothetical Nabla operator lead to significant simplifications of the related calculations.For r/λ = 1 the angular difference is 0.158 rad resp.9°and for r/λ > 3 less than 0.053 rad resp.3°.

Conclusions
With the spherical one-way wave equation and the hypothetical Nabla operator being presented in this paper the calculation of spherical wave propagation does not show wattless near fields and the solutions can be derived without the need of Bessel functions.Thereby spherical wave calculations can be significantly simplified.
ρ s s s − E∆s s s = 0 0 0 (1) By merging density ρ [kg/m 3 ] and elasticity module E sound velocity c = E/ρ [m/s] results and the governing wave equation for a homogeneous medium follows: s s s − c 2 ∆s s s = 0 0 0 (2)

Table 1 .
Comparison of force based Two-Way Wave Equation and impulse based One-Way Wave Equation in an unlimited, elastic 3D contiuum without boundary effects.For simplicity and clarification linearity, homogeneity, isotropy and losslessness are assumed; s s s[m] = elastic displacement, ṡ s s = ∂s s s/∂t [m/s] = particle velocity; a b a b a b, a a a × b b b, a a a • b b b = dyadic, vectorial, scalar product of the two vectors a a a and b b b.
) serves as the starting point for the deflection s s s = s e r e r e r s s s − c 2 ∆s s s = 0 (21) With s = −ω 2 s = −k 2 c 2 s and the spherical Laplace operator ∆s = s + 2s /r the Bessel differential equation follows: r 2 s + 2rs + (kr) 2 s = 0 (22) with bessel function of half-integer order n = 1/2 as solution s

Figure 1 .
Figure 1.Comparison of the impedance amplitude functions A conv and A hyp for spherical wave propagation as a function of the normalized distance r/λ [-].At r/λ = 1 the difference is 1.2%, for r/λ > 3 it is less than 0.14%.In practice, the near field influence is neglected in calculations.