_{1}

^{*}

This article provides a closed form solution to the telegrapher’s equation with three space variables defined on a subset of a sphere within two radii, two azimuthal angles and one polar angle. The Dirichlet problem for general boundary conditions is solved in detail, on the basis of which Neumann and Robin conditions are easily handled. The solution to the simpler problem in cylindrical coordinates is also provided. Ways to efficiently implement the formulae are explained. Minor adjustments result in solutions to the wave equation and to the heat equation on the same domain as well, since the latter are particular cases of the more general telegrapher’s equation.

Hyperbolic partial differential equations (PDEs) are among the three main classes of PDEs in mathematical physics. The most extensively studied hyperbolic PDE is the wave equation. Another major hyperbolic PDE is the telegrapher’s equation. In rectangular coordinates, its canonical form is the following:

∂ 2 u ( x , t ) ∂ t 2 + κ 1 ∂ u ( x , t ) ∂ t = κ 2 Δ n u ( x , t ) − κ 3 u ( x , t ) + Λ ( x , t ) (1)

where x is a vector of n space variables, Δ n is the n-dimensional Laplacian operator, Λ is a “sufficiently” smooth function and κ 1 , κ 2 , κ 3 are three positive constants.

A good introduction to the fundamental properties and main applications of this PDE can be found in Zauderer [

Among all the major second order linear PDEs of mathematical physics, the telegrapher’s equation is noticeable for its generality. Indeed, all the most important non-stationary linear equations can be retrieved from (1). In particular, the wave equation is a special case of the telegrapher’s equation obtained by taking κ 1 = κ 3 = Λ ( x , t ) = 0 . The Klein-Gordon equation for a function w ( x , t ) is obtained through the following simple relation:

w ( x , t ) = exp ( − κ 1 t 2 ) u ( x , t ) (2)

An obvious parameterization of the second-order time derivative in (1) would also enable to recover the parabolic heat equation. For more general background on PDEs related to Equation (1), one can refer to [

A number of initial boundary value problems (IBVPs) involving PDE (1) have already been analytically solved, especially in rectangular coordinates. Certainly the most comprehensive compilation of known analytical formulae is the handbook by Polyanin and Nazaikinskii [

{ 0 < r 1 ≤ r ≤ r 2 < ∞ , 0 < ϕ ≤ 2 π , 0 < θ ≤ π , t ≥ 0 } (3)

Such a restriction rules out problems that require functional dependency on the boundaries of the variables ϕ or θ as well as on the boundary of r. Of course, one can always resort to the usual numerical techniques of approximation, but, in this case, it is not easy to design efficient and accurate finite difference schemes, due to the dimension of the problem and to the complexity of the telegrapher’s equation in a spherical domain. The main contribution of this paper is to provide a closed form solution for a class of IBVPs involving the linear telegrapher’s equation with three space variables and general boundary conditions on the following domain:

D = { 0 < r 1 ≤ r ≤ r 2 < ∞ , 0 < θ ≤ θ 1 ≤ π , 0 < ϕ 1 ≤ ϕ ≤ ϕ 2 ≤ 2 π , t ≥ 0 } (4)

The domain in (4) allows for partial ranges in both angular variables. As such, it includes subsets of the sphere with many more shapes than the domain in (3), so that more flexible models can be considered, with a greater number of physical specifications.

The linear and separable nature of the PDE under consideration enables the application of classical methods of separation of variables and eigenfunction expansions. One of the main advantages of this approach is that it easily allows to handle non-homogeneity in the equation or in the boundary conditions, once the homogeneous solution has been obtained. However, the solution process involves a non-standard generalized Fourier expansion of the first order initial condition, as well as computational difficulties associated with the integration of associated Legendre functions with non-integral parameters, which may explain why there is currently no available formula in the literature.

It must be emphasized that the telegrapher’s equation is the focus of this article due to its high degree of generality, not because a specific physical application is studied. Thus, it is expected that the analytical results in this article can be fruitfully applied to parabolic problems—after adjustment of the initial conditions, as well as to hyperbolic problems that do not model the same physical phenomena as those initially described by the telegrapher’s equation.

This paper is organized as follows. In Section 2, the main result is stated under the title “Result 1”, i.e. the solution of the three-dimensional telegrapher’s equation on a spherical domain with boundary conditions on all the space variables. The rest of Section 2 is dedicated to the proof of Result 1, which derives from the combination of two lemmas, along with comments on the numerical implementation of the formula. Section 3 presents two corollaries that handle non-homogeneity of the equation and of the boundary conditions, and then states, under the title “Result 2”, the solution of the three-dimensional telegrapher’s equation on a cylindrical domain with boundary conditions on all the space variables.

The main result of this article is now stated as “Result 1”.

Result 1

Solution to the Dirichlet problem for the telegrapher’s equation with three space variables on a subset of a sphere located within two radii, two azimuthal angles and one polar angle.

Let u ( r , ϕ , θ , t ) be the solution to the following partial differential equation:

∂ 2 u ∂ t 2 + κ 1 ∂ u ∂ t = κ 2 ( ∂ 2 u ∂ r 2 + 2 r ∂ u ∂ r + 1 r 2 ( ∂ 2 u ∂ θ 2 + cot ( θ ) ∂ u ∂ θ + csc 2 ( θ ) ∂ 2 u ∂ ϕ 2 ) ) − κ 3 u (5)

on the domain:

D = { 0 < r 1 ≤ r ≤ r 2 < ∞ , 0 < θ ≤ θ 1 ≤ π , 0 < ϕ 1 ≤ ϕ ≤ ϕ 2 ≤ 2 π , t ≥ 0 } (6)

with the following boundary conditions:

u ( r , θ , ϕ , 0 ) = 0 (7)

∂ u ( r , θ , ϕ , t ) ∂ t | t = 0 = f ( r , θ , ϕ ) (8)

u ( r 1 , θ , ϕ , t ) = u ( r 2 , θ , ϕ , t ) = 0 (9)

u ( r , θ 1 , ϕ , t ) = 0 (10)

u ( r , θ , ϕ 1 , t ) = u ( r , θ , ϕ 2 , t ) = 0 (11)

where κ 1 , κ 2 , κ 3 are three positive constants and f is a continuous function from D to ℝ .

Let us denote:

· P ν μ ( x ) as the associated Legendre function of the first kind with real variable x, order μ and degree ν

· J ν ( x ) and Y ν ( x ) as the Bessel functions of the first kind and second kind, respectively, with real order ν

· j ν ( x ) and y ν ( x ) as the spherical Bessel functions of the first kind and second kind, respectively, with real order ν

· Γ ( x ) as the Gamma function with real variable x

Then, the function u ( r , ϕ , θ , t ) is given by:

u ( r , θ , ϕ , t ) = ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ b m n p T ( t ) Φ m ( ϕ ) Θ n ( θ ) R p ( r ) (12)

where the functions T ( t ) , Φ m ( ϕ ) , Θ n ( θ ) and R p ( r ) are defined as follows:

T ( t ) = exp ( − κ 1 2 t ) sin ( t κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) (13)

Φ m ( ϕ ) = sin ( β m ( ϕ − ϕ 1 ) ) (14)

Θ n ( θ ) = P ν n β m ( cos θ ) (15)

R p ( r ) = π 2 λ p r 1 J ν n + 1 / 2 ( λ p r 1 ) π 2 λ p r Y ν n + 1 / 2 ( λ p r ) − π 2 λ p r 1 Y ν n + 1 / 2 ( λ p r 1 ) π 2 λ p r J ν n + 1 / 2 ( λ p r ) (16)

The numbers β m , ν n and λ p are defined as follows:

· ∀ m ∈ ℕ , β m is the m-th root of the equation

sin ( β ϕ 2 ) cos ( β ϕ 1 ) − sin ( β ϕ 1 ) cos ( β ϕ 2 ) = 0 (17)

· ∀ n ∈ ℕ , ν n is the n-th root of the equation

P ν β m ( cos θ 1 ) = 0 (18)

· ∀ p ∈ ℕ , λ p is the p-th root of the equation

y ν n ( λ r 2 ) j ν n ( λ r 1 ) − y ν n ( λ r 1 ) j ν n ( λ r 2 ) = 0 (19)

The coefficient b m n p is given by:

b m n p = 4 ( κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) − 1 / 2 β m r 1 λ p 4 J ν n + 1 / 2 2 ( λ p r 2 ) ( β m ( ϕ 2 − ϕ 1 ) − sin ( 2 β m ( ϕ 2 − ϕ 1 ) ) ) ( J ν n + 1 / 2 2 ( λ p r 1 ) − J ν n + 1 / 2 2 ( λ p r 2 ) ) ‖ P ν n β m ‖ 2 × ∫ r 1 r 2 ∫ 0 θ 1 ∫ ϕ 1 ϕ 2 f ( φ , ϑ , ρ ) sin ( λ p ( φ − ϕ 1 ) ) P ν n β m ( cos ϑ ) R p ( ρ ) ρ 2 sin ( ϑ ) d φ d ϑ d ρ (20)

where:

‖ P ν n β m ‖ 2 = ∫ 0 π ( cot ( ϑ 2 ) ) 2 β m ( ∑ k = 0 ∞ ( − 1 ) k Γ ( ν n + k + 1 ) Γ ( ν n − k + 1 ) sin ( ϑ 2 ) 2 k Γ ( k + 1 ) Γ ( k − β m + 1 ) ) 2 sin ( ϑ ) d ϑ (21)

which completes the statement of Result 1.

Proof of Result 1:

A solution to the boundary value problem under consideration can be found by the techniques of separation of variables and eigenfunction expansions. The proof consists of two lemmas.

Lemma 1

A solution which verifies both PDE (5) and the boundary conditions associated with the space variables (9)-(11) is given by:

u ( r , θ , ϕ , t ) = exp ( − κ 1 2 t ) ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ ( a m n p cos ( t κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) + b m n p sin ( t κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) ) Φ m ( ϕ ) Θ n ( θ ) R p ( r ) (22)

where a m n p and b m n p are coefficients to be determined, and the functions T ( t ) , Φ m ( ϕ ) , Θ n ( θ ) and R p ( r ) are given in Result 1.

Proof of Lemma 1

The substitution u ( r , ϕ , θ , t ) = R ( r ) Φ ( ϕ ) Θ ( θ ) T ( t ) carried out into PDE (5) yields:

1 κ 2 T ″ T + κ 1 κ 2 T ′ T = R ″ R + 2 r R ′ R + 1 r 2 ( Θ ″ Θ + ( cot θ ) Θ ′ Θ + ( csc 2 θ ) Φ ″ Φ ) − κ 3 κ 2 (23)

Equating both sides of Equation (23) to − λ 2 , where λ 2 is a first separation constant, leads to the following two equations:

1 κ 2 T ″ T + κ 1 κ 2 T ′ T + λ 2 + κ 3 κ 2 = 0 (24)

r 2 R ″ R + 2 r R ′ R + r 2 λ 2 = − ( Θ ″ Θ + ( cot θ ) Θ ′ Θ + ( csc 2 θ ) Φ ″ Φ ) (25)

The complex roots of the characteristic equation associated with the differential Equation (24) easily yield:

T ( t ) = a 1 exp ( − κ 1 2 t ) cos ( t κ 1 2 − 4 ( κ 2 λ 2 − κ 3 ) ) + a 2 exp ( − κ 1 2 t ) sin ( t κ 1 2 − 4 ( κ 2 λ 2 − κ 3 ) ) (26)

where a 1 and a 2 are arbitrary constants.

Equating both sides of Equation (25) with a second separation constant μ 2 leads to the following two differential equations:

r 2 R ″ + 2 r R ′ + ( r 2 λ 2 − μ 2 ) R = 0 (27)

Θ ″ Θ + ( cot θ ) Θ ′ Θ + ( csc 2 θ ) Φ ″ Φ + μ 2 = 0 (28)

Equation (27) is a spherical Bessel equation, while Equation (28) is a spherical harmonics equation.

Rearranging Equation (28) and equating both sides with − β 2 , where β 2 is a third separation constant, results in the following two differential equations:

Φ ″ + β 2 Φ = 0 (29)

Θ ″ + ( cot θ ) Θ ′ + ( μ 2 − β 2 csc 2 θ ) Θ = 0 (30)

The general solution to Equation (29) is given by:

Φ ( ϕ ) = b 1 cos ( β ϕ ) + b 2 sin ( β ϕ ) (31)

where b 1 and b 2 are arbitrary constants.

The boundary condition (11) yields:

{ b 1 cos ( β ϕ 1 ) + b 2 sin ( β ϕ 1 ) = 0 ( 32 ) b 1 cos ( β ϕ 2 ) + b 2 sin ( β ϕ 2 ) = 0 ( 33 )

Substituting b 1 from Equation (32) into Equation (33), we get:

b 2 ( sin ( β ϕ 2 ) − sin ( β ϕ 1 ) cos ( β ϕ 2 ) cos ( β ϕ 1 ) ) = 0 (34)

Since b 2 = 0 entails the trivial solution, β must solve:

sin ( β ϕ 2 ) cos ( β ϕ 1 ) − sin ( β ϕ 1 ) cos ( β ϕ 2 ) = 0 (35)

Classical Sturm-Liouville theory [

Eigenfunctions Φ m ( ϕ ) that accommodate the boundary conditions at ϕ 1 and ϕ 2 can thus be written as follows:

Φ m ( ϕ ) = cos ( β m ϕ 1 ) sin ( β m ϕ ) − sin ( β m ϕ 1 ) cos ( β m ϕ ) = sin ( β m ( ϕ − ϕ 1 ) ) (36)

Equation (30) is the polar form of the associated Legendre differential equation. It is well-known [

Θ ( θ ) = c 1 P ν β m ( cos θ ) + c 2 Q ν β m ( cos θ ) (37)

where Q ν β m is the associated Legendre function of the second kind.

The latter solution has to be eliminated if we want bounded solutions, since Q ν β m ( cos θ ) goes to infinity as θ approaches 0 or π . The boundary condition (10) yields:

P ν β m ( cos θ 1 ) = 0 (38)

Since θ 1 is fixed and β m has already been determined, Equation (38) has to be solved for ν . A convenient way of doing this is to use the Mehler-Dirichlet representation for the associated Legendre function P ν μ ( cos θ ) , which is valid for any θ between zero and π :

P ν μ ( cos θ ) = ( sin θ ) μ 2 π Γ ( 1 2 − μ ) ∫ 0 θ cos ( ( ν + 1 2 ) ϑ ) ( cos ϑ − cos θ ) μ + 1 2 d ϑ (39)

Some background on the formula (39) can be found in [

P − ν − 1 μ ( z ) = P ν μ ( z ) (40)

Other interesting properties of the roots of associated Legendre functions with non-integral parameters can be found in [

Each eigenfunction Θ n ( θ ) can thus be written as:

Θ n ( θ ) = P ν n β m ( cos θ ) (41)

The general solution to Equation (27) is given by:

R ( r ) = d 1 j ν n ( λ r ) + d 2 y ν n ( λ r ) (42)

There are no singularity points associated with the functions in the right hand side of (42), as the variable r is strictly positive. It is then easily seen, applying classical Sturm-Liouville theory [

y ν n ( λ r 2 ) j ν n ( λ r 1 ) − y ν n ( λ r 1 ) j ν n ( λ r 2 ) = 0 (43)

The zeros of Bessel functions have been extensively studied and this knowledge is embedded in the built-in procedures of standard scientific computing packages such as the BesselJZeros and BesselYZeros commands in Maple.

Eigenfunctions R p ( r ) that accommodate the boundary conditions at r 1 and r 2 can thus be written as follows:

R p ( r ) = j ν n ( λ p r 1 ) y ν n ( λ p r ) − y ν n ( λ p r 1 ) j ν n ( λ p r ) (44)

= π 2 λ p r 1 J ν n + 1 / 2 ( λ p r 1 ) π 2 λ p r Y ν n + 1 / 2 ( λ p r ) − π 2 λ p r 1 Y ν n + 1 / 2 ( λ p r 1 ) π 2 λ p r J ν n + 1 / 2 ( λ p r ) (45)

which completes the proof of Lemma 1.

The next step consists in finding a way to satisfy the initial conditions (7) and (8).

First, the initial condition (7) immediately yields a m n p = 0 , so that we obtain:

u ( r , θ , ϕ , t ) = ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ b m n p T ( t ) Φ m ( ϕ ) Θ n ( θ ) R p ( r ) (46)

where T ( t ) is given by (13), Φ m ( ϕ ) is given by (36), Θ n ( θ ) is given by (41), R p ( r ) is given by (45).

Next, the substitution of the initial condition (8) into (46) yields:

f ( r , θ , ϕ ) = ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ b m n p κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) Φ m ( ϕ ) Θ n ( θ ) R p ( r ) (47)

Lemma 2 shows how to expand f ( r , θ , ϕ ) into an appropriate Fourier-Legendre-Bessel triple series, thus allowing to determine the coefficients b m n p .

Lemma 2

Let f ( r , θ , ϕ ) be a continuous real function defined on D. Then, f ( r , θ , ϕ ) admits the following series expansion:

f ( r , θ , ϕ ) = ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ d m n p Φ m ( ϕ ) Θ n ( θ ) R p ( r ) (48)

where:

d m n p = 4 β m r 1 λ p 4 J ν n + 1 / 2 2 ( λ p r 2 ) ( β m ( ϕ 2 − ϕ 1 ) − sin ( 2 β m ( ϕ 2 − ϕ 1 ) ) ) ( J ν n + 1 / 2 2 ( λ p r 1 ) − J ν n + 1 / 2 2 ( λ p r 2 ) ) ‖ P ν n β m ‖ 2 (49)

∫ r 1 r 2 ∫ 0 θ 1 ∫ ϕ 1 ϕ 2 f ( φ , ϑ , ρ ) sin ( λ p ( φ − ϕ 1 ) ) P ν n β m ( cos ϑ ) R p ( ρ ) ρ 2 sin ( ϑ ) d φ d ϑ d ρ (50)

and ‖ P ν n β m ‖ 2 is given in (21).

Proof of Lemma 2:

Based on the separated form of the product Φ ( ϕ ) Θ ( θ ) R ( r ) , an appropriate approach is to break down the coefficient d m n p into the following product of three coefficients:

d m n p = d m d n d p (51)

and to solve successively for d m , d n and d p .

It is easy to verify that the set of functions { Φ m ( ϕ ) , m ∈ ℕ } is orthonormal in [ ϕ 1 , ϕ 2 ] with respect to the following norm:

‖ Φ m ‖ 2 ≜ ∫ ϕ 1 ϕ 2 [ cos ( β m ϕ 1 ) sin ( β m φ ) − sin ( β m ϕ 1 ) cos ( β m φ ) ] 2 d φ (52)

= 1 2 ( ϕ 2 − ϕ 1 ) − 1 2 β m sin ( 2 β m ( ϕ 2 − ϕ 1 ) ) (53)

Therefore, considering r and θ as constants, (48) can be rewritten as:

f ( r , θ , ϕ ) = ∑ m = 1 ∞ d m Φ m ( ϕ ) (54)

where:

d m = ∑ n = 1 ∞ ∑ p = 1 ∞ d n p Θ n ( θ ) R p ( r ) (55)

The coefficients d m in (54) are the Fourier coefficients of the expansion of f ( r , θ , ϕ ) into a series of sine functions Φ m ( ϕ ) for variable ϕ and fixed θ and r. Classical Fourier analysis yields:

d m = 1 ‖ Φ m ‖ 2 ∫ ϕ 1 ϕ 2 f ( φ , θ , r ) sin ( β m ( φ − ϕ 1 ) ) d φ (55)

where ‖ Φ m ‖ 2 is given by (53), and (54) becomes:

f ( r , θ , ϕ ) = ∑ m = 1 ∞ ( ∑ n = 1 ∞ ∑ p = 1 ∞ d n p Θ n ( θ ) R p ( r ) ) d m Φ m ( ϕ ) (56)

where d m is given by (55) and the coefficient d n p remains to be determined.

Given d m and considering r as a constant, we have:

∑ n = 1 ∞ ∑ p = 1 ∞ d n p Θ n ( θ ) R p ( r ) = ∑ n = 1 ∞ d n Θ n ( θ ) (57)

where:

d n = ∑ p = 1 ∞ d p R p ( r ) (58)

For a given m, the set of functions { Θ n ( θ ) , n ∈ ℕ } is orthogonal in [ 0 , π ] with respect to a weight given by the sine function, i.e.:

∫ 0 π P ν n 1 β m ( cos ϑ ) P ν n 2 β m ( cos ϑ ) sin ( ϑ ) d ϑ = 0 , n 1 ≠ n 2 (59)

The set of functions { Θ n ( θ ) , n ∈ ℕ } is orthonormal in [ 0 , π ] with respect to the following norm:

‖ P ν n β m ‖ 2 ≜ ∫ 0 π ( P ν n β m ( cos ϑ ) ) 2 sin ( ϑ ) d ϑ (60)

The coefficients d n are the Fourier coefficients of the expansion of f ( r , θ , ϕ ) into a series of functions Θ n ( θ ) for variable θ and fixed r and ϕ . In practice, one needs to be able to perform the numerical evaluation of the coefficients d n . The function P ν n β m ( . ) is an associated Legendre function of the first kind with non-integral order β m and non-integral degree ν n . Its moment-generating function is undefined and it is therefore less easy to evaluate the norm ‖ P ν n β m ‖ 2 than if we were dealing with a standard associated Legendre function with integral parameters. However, any function P ν n β m ( z ) is a solution to a second order, linear ordinary differential equation, with three regular singular points and, as such, it can be recast as into a hypergeometric function F 2 1 ( a , b ; c ; z ) as follows [

P ν n β m ( z ) = 1 Γ ( 1 − β m ) ( z + 1 z − 1 ) β m 2 F 2 1 ( − ν n , ν n + 1 ; 1 − β m ; 1 − z 2 ) (61)

The function F 2 1 ( a , b ; c ; z ) , also known as the Gauss hypergeometric function, admits the following series expansion [

F 2 1 ( a , b ; c ; z ) = ∑ k = 0 ∞ ( a ) k ( b ) k ( c ) k z k k ! = Γ ( c ) Γ ( a ) Γ ( b ) ∑ k = 0 ∞ Γ ( a + k ) Γ ( b + k ) Γ ( c + k ) z k k ! (62)

This series representation is well-suited to our problem, since it converges for | z | < 1 .

Letting z = cos ( ϑ ) in (61) and performing a little algebra, one can obtain the following formula for the coefficients d n :

d n = 1 ‖ P ν n β m ‖ 2 ∫ 0 π f ( ϕ , ϑ , r ) P ν n β m ( cos ϑ ) sin ( ϑ ) d ϑ (63)

where:

P ν n β m ( cos ϑ ) = ( cot ( ϑ 2 ) ) β m ∑ k = 0 ∞ ( − 1 ) k Γ ( ν n + k + 1 ) Γ ( ν n − k + 1 ) sin ( ϑ 2 ) 2 k Γ ( k + 1 ) Γ ( k − β m + 1 ) (64)

and ‖ P ν n β m ‖ 2 is given in (21).

A similar trigonometric series is derived by Bremer [

Therefore, (56) becomes:

f ( r , θ , ϕ ) = ∑ m = 1 ∞ ( ∑ n = 1 ∞ ( ∑ p = 1 ∞ d p R p ( r ) ) d n Θ n ( θ ) ) d m Φ m ( ϕ ) (65)

where d m is given by (55), d n is given by (63), and the coefficient d p is yet to be determined.

For a given ν n , one can verify that the set of functions { R p ( r ) , p ∈ ℕ } is orthogonal in [ r 1 , r 2 ] with respect to the weight function w ( r ) = r 2 , i.e. we have:

∫ r 1 r 2 [ j ν n ( λ p 1 r 1 ) y ν n ( λ p 1 ρ ) − y ν n ( λ p 1 r 1 ) j ν n ( λ p 1 ρ ) ] p 1 ≠ p 2 [ j ν n ( λ p 2 r 1 ) y ν n ( λ p 2 ρ ) − y ν n ( λ p 2 r 1 ) j ν n ( λ p 2 ρ ) ] ρ 2 d ρ = 0 (66)

The set of functions { R p ( r ) , p ∈ ℕ } is orthonormal with respect to the following norm:

‖ R p ‖ 2 = ∫ r 1 r 2 [ j ν n ( λ p r 1 ) y ν n ( λ p ρ ) − y ν n ( λ p r 1 ) j ν n ( λ p ρ ) ] 2 ρ 2 d ρ (67)

Therefore, the coefficients d p are the Fourier coefficients of the expansion of f ( r , θ , ϕ ) into a series of the combination of spherical Bessel functions defined by the function R p ( r ) , for variable r and fixed ϕ and θ , with each d p given by:

d p = 1 ‖ R p ‖ 2 ∫ r 1 r 2 f ( ϕ , θ , ρ ) [ j ν n ( λ p r 1 ) y ν n ( λ p ρ ) − y ν n ( λ p r 1 ) j ν n ( λ p ρ ) ] ρ 2 d ρ (68)

Using classical identities verified by Bessel functions, one can obtain the following formula for ‖ R p ‖ 2 :

‖ R p ‖ 2 = J ν n + 1 / 2 2 ( λ p r 1 ) − J ν n + 1 / 2 2 ( λ p r 2 ) 2 r 1 λ p 4 J ν n + 1 / 2 2 ( λ p r 2 ) (69)

Hence, the coefficients of the expansion of f ( r , θ , ϕ ) into the Fourier-Legendre-Bessel triple series generated by the product of functions Φ m ( ϕ ) Θ n ( θ ) R p ( r ) are given by d m n p = d m d n d p , where d m , d n and d p are given by (55), (63) and (68), respectively, and Lemma 2 ensues.

The combination of Lemma 1 and Lemma 2 entails Result 1.

Remark: using standard arguments, it can be proven that the continuity of the function f ( r , θ , ϕ ) ensures uniform convergence of the triple series in (48), just as it ensures the existence of a unique solution in the classical sense of the IBVP (5)-(11).

Several useful additional results can be inferred from Result 1. First, once can derive the solution to the non-homogeneous telegrapher’s equation, as given in Corollary 1.

Corollary 1

Let h ( r , ϕ , θ , t ) be a continuous function from the domain D in (6) to ℝ .

The solution of the nonhomogeneous equation:

∂ 2 u ∂ t 2 + κ 1 ∂ u ∂ t = κ 2 ( ∂ 2 u ∂ r 2 + 2 r ∂ u ∂ r + 1 r 2 ( ∂ 2 u ∂ θ 2 + cot θ ∂ u ∂ θ + csc 2 θ ∂ 2 u ∂ ϕ 2 ) ) − κ 3 u + h ( r , ϕ , θ , t ) (70)

that verifies the homogeneous boundary conditions stated in (7)-(11), is given by:

u ( r , θ , ϕ , t ) = ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ r 1 r 2 f ( ρ , ϑ , φ ) G ( r , ρ , ϕ , φ , θ , ϑ , t ) ρ 2 sin ( ϑ ) d ρ d ϑ d φ + ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ r 1 r 2 ∫ 0 t h ( ρ , ϑ , φ , s ) G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ρ 2 sin ( ϑ ) d ρ d ϑ d φ d s (71)

where G ( r , ρ , ϕ , φ , θ , ϑ , t ) is the Green’s function, which is immediately derived from Result 1:

G ( r , ρ , ϕ , φ , θ , ϑ , t ) = exp ( − κ 1 2 t ) ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ sin ( t κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) × 4 β m r 1 λ p 4 J ν n + 1 / 2 2 ( λ p r 2 ) ( κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) 1 / 2 ( β m ( ϕ 2 − ϕ 1 ) − sin ( 2 β m ( ϕ 2 − ϕ 1 ) ) ) ( J ν n + 1 / 2 2 ( λ p r 1 ) − J ν n + 1 / 2 2 ( λ p r 2 ) ) ‖ P ν n β m ‖ 2 × sin ( β m ( ϕ − ϕ 1 ) ) sin ( β m ( φ − ϕ 1 ) ) P ν n β m ( cos θ ) P ν n β m ( cos ϑ ) R p ( r ) R p ( ρ ) (72)

which completes the statement of Corollary 1.

Endowed with the Green’s function, one can also obtain a formula for the Dirichlet problem with non-homogeneous boundary conditions.

Corollary 2

Let u ( r , θ , ϕ , t ) be defined on the domain D in (6) and satisfy PDE (70) as well as the following nonhomogeneous boundary conditions:

u ( r , θ , ϕ , 0 ) = f 1 ( r , ϕ , θ ) (73)

∂ u ( r , θ , ϕ , t ) ∂ t | t = 0 = f 2 ( r , θ , ϕ ) (74)

u ( r 1 , θ , ϕ , t ) = f 3 ( θ , ϕ , t ) (75)

u ( r 2 , θ , ϕ , t ) = f 4 ( θ , ϕ , t ) (76)

u ( r , θ 1 , ϕ , t ) = f 5 ( r , ϕ , t ) (77)

u ( r , θ , ϕ 1 , t ) = f 6 ( r , θ , t ) (78)

u ( r , θ , ϕ 2 , t ) = f 7 ( r , θ , t ) (79)

where f 1 , ⋯ , f 7 are all continuous real functions such that:

f 1 ( r , θ , ϕ ) = f 3 ( θ , ϕ , 0 ) = f 4 ( θ , ϕ , 0 ) = f 5 ( r , ϕ , 0 ) = f 6 ( r , θ , 0 ) = f 7 ( r , θ , 0 ) (80)

f 3 ( θ 1 , ϕ , t ) = f 5 ( r 1 , ϕ , t ) , f 4 ( θ 1 , ϕ , t ) = f 5 ( r 2 , ϕ , t ) (81)

f 3 ( θ , ϕ 1 , t ) = f 6 ( r 1 , θ , t ) , f 4 ( θ , ϕ 1 , t ) = f 6 ( r 2 , θ , t ) (82)

f 3 ( θ , ϕ 2 , t ) = f 7 ( r 1 , θ , t ) , f 4 ( θ , ϕ 2 , t ) = f 7 ( r 2 , θ , t ) (83)

f 5 ( r 1 , ϕ 1 , t ) = f 6 ( r 1 , θ 1 , t ) , f 5 ( r 1 , ϕ 2 , t ) = f 7 ( r 1 , θ 1 , t ) (84)

f 5 ( r 2 , ϕ 1 , t ) = f 6 ( r 2 , θ 1 , t ) , f 5 ( r 2 , ϕ 2 , t ) = f 7 ( r 2 , θ 1 , t ) (85)

Then, the function u ( r , θ , ϕ , t ) is given by:

u ( r , θ , ϕ , t ) = ∂ ∂ t ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ r 1 r 2 f 1 ( ρ , ϑ , φ ) G ( r , ρ , ϕ , φ , θ , ϑ , t ) ρ 2 sin ( ϑ ) d ρ d ϑ d φ (86)

+ ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ r 1 r 2 ( f 2 ( ρ , ϑ , φ ) + κ 1 f 1 ( ρ , ϑ , φ ) ) G ( r , ρ , ϕ , φ , θ , ϑ , t ) ρ 2 sin ( ϑ ) d ρ d ϑ d φ (87)

+ κ 2 2 r 1 2 ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ 0 t f 3 ( ϑ , φ , s ) ∂ G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ∂ ρ | ρ = r 1 sin ( ϑ ) d ϑ d φ d s (88)

− κ 2 2 r 2 2 ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ 0 t f 4 ( ϑ , φ , s ) ∂ G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ∂ ρ | ρ = r 2 sin ( ϑ ) d ϑ d φ d s (89)

− κ 2 2 sin ( θ 1 ) ∫ 0 t ∫ ϕ 1 ϕ 2 ∫ r 1 r 2 f 5 ( s , φ , ρ ) ∂ G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ∂ ϑ | ϑ = θ 1 ρ 2 d ϑ d φ d s (90)

+ κ 2 2 ∫ 0 t ∫ 0 θ 1 ∫ r 1 r 2 f 6 ( s , ϑ , ρ ) ∂ G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ∂ φ | φ = ϕ 1 ρ 2 sin ( ϑ ) d ϑ d φ d s (91)

− κ 2 2 ∫ 0 t ∫ 0 θ 1 ∫ r 1 r 2 f 7 ( s , ϑ , ρ ) ∂ G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ∂ φ | φ = ϕ 2 ρ 2 sin ( ϑ ) d ϑ d φ d s (92)

+ ∫ ϕ 1 ϕ 2 ∫ 0 θ 1 ∫ r 1 r 2 ∫ 0 t h ( ρ , ϑ , φ , s ) G ( r , ρ , ϕ , φ , θ , ϑ , t − s ) ρ 2 sin ( ϑ ) d ρ d ϑ d φ d s (93)

which ends the statement of Corollary 2.

The theory underpinning the derivation of Corollary 1 and Corollary 2 is explained, e.g., in [

It should also be pointed out that Neumann or Robin boundary conditions can be easily handled in the same framework as the one used to solve the Dirichlet problem in Section 2. For example, if the radial boundary condition in (9) is replaced with the following:

∂ ∂ r u ( r , θ , ϕ , t ) | r = r 1 = ∂ ∂ r u ( r , θ , ϕ , t ) | r = r 2 = 0 (94)

Then, each λ p becomes the p-th root of the following equation:

y ′ n ( λ r 2 ) j ′ n ( λ r 1 ) − y ′ n ( λ r 1 ) j ′ n ( λ r 2 ) = 0 (95)

and the eigenfunctions of the Sturm-Liouville problem associated with R ( r ) become:

R p ( r ) = j ′ n ( λ p r 1 ) y n ( λ p r ) − y ′ n ( λ p r 1 ) j n ( λ p r ) (96)

Likewise, if the azimuthal angular boundary condition in (11) is replaced with the following:

∂ u ( r , θ , ϕ , t ) ∂ ϕ | ϕ = ϕ 1 = ∂ u ( r , θ , ϕ , t ) ∂ ϕ | ϕ = ϕ 2 = 0 (97)

Then, each β m becomes the m-th root of the following equation:

cos ( β ϕ 2 ) sin ( β ϕ 1 ) − cos ( β ϕ 1 ) sin ( β ϕ 2 ) = 0 (98)

and the eigenfunctions of the Sturm-Liouville problem associated with Φ ( ϕ ) become:

Φ m ( ϕ ) = β m ( b 1 sin ( − β m ϕ 1 ) sin ( β m ϕ ) − b 2 cos ( β m ϕ 1 ) cos ( β m ϕ ) ) (99)

Finally, the solution to the three-dimensional telegrapher’s problem in cylindrical coordinates can also be computed and is given in the following Result 2.

Result 2

Solution to the Dirichlet problem for the telegrapher’s equation with three space variables on a subset of a cylinder located within two radii, two angles and two heights.

Let u ( r , ϕ , z , t ) be the solution of the following partial differential equation:

∂ 2 u ∂ t 2 + κ 1 ∂ u ∂ t = κ 2 ( ∂ 2 u ∂ r 2 + 1 r ∂ u ∂ r + 1 r 2 ∂ 2 u ∂ ϕ 2 + ∂ 2 u ∂ z 2 ) − κ 3 u (100)

on the domain:

D = { 0 ≤ r 1 ≤ r ≤ r 2 < ∞ , 0 ≤ ϕ 1 ≤ ϕ ≤ ϕ 2 ≤ 2 π , 0 ≤ z 1 ≤ z ≤ z 2 < ∞ , t ≥ 0 } (101)

with the following boundary conditions:

u ( r , ϕ , z , 0 ) = 0 (102)

∂ u ( r , ϕ , z , t ) ∂ t | t = 0 = f ( r , ϕ , z ) (103)

u ( r 1 , ϕ , z , t ) = u ( r 2 , ϕ , z , t ) = 0 (104)

u ( r , ϕ 1 , z , t ) = u ( r , ϕ 2 , z , t ) = 0 (105)

u ( r , ϕ , z 1 , t ) = u ( r , ϕ , z 2 , t ) = 0 (106)

where κ 1 , κ 2 , κ 3 are three positive constants and f is a continuous function from D to ℝ .

Then, the function u ( r , ϕ , z , t ) is given by:

u ( r , ϕ , z , t ) = ∑ m = 1 ∞ ∑ n = 1 ∞ ∑ p = 1 ∞ b m n p T ( t ) Φ m ( ϕ ) R n ( r ) Z p ( z ) (107)

The functions T ( t ) , Φ m ( ϕ ) , R n ( r ) and Z p ( z ) are defined as follows:

T ( t ) = exp ( − κ 1 2 t ) sin ( t κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) ) (108)

Φ m ( ϕ ) = sin ( β m ( ϕ − ϕ 1 ) ) (109)

R n ( r ) = Y β m ( μ n r ) J β m ( μ n r 1 ) − Y β m ( μ n r 1 ) J β m ( μ n r ) (110)

Z p ( z ) = sin ( λ p 2 − μ n 2 ( z − z 1 ) ) (111)

The numbers β m , μ n and λ p are defined as follows:

· ∀ m ∈ ℕ , β m is the m-th root of the following equation

sin ( β ϕ 2 ) cos ( β ϕ 1 ) − sin ( β ϕ 1 ) cos ( β ϕ 2 ) = 0 (112)

· ∀ n ∈ ℕ , μ n is the n-th root of the following equation

Y β m ( μ r 2 ) J β m ( μ r 1 ) − Y β m ( μ r 1 ) J β m ( μ r 2 ) = 0 (113)

· ∀ p ∈ ℕ , λ p is the p-th root of the following equation

sin ( λ 2 − μ n 2 z 2 ) cos ( λ 2 − μ n 2 z 1 ) − sin ( λ 2 − μ n 2 z 1 ) cos ( λ 2 − μ n 2 z 2 ) = 0 (114)

The coefficient b m n p is given by:

b m n p = 8 β m λ p 2 − μ n 2 ( β m ( ϕ 2 − ϕ 1 ) − sin ( 2 β m ( ϕ 2 − ϕ 1 ) ) ) ( λ p 2 − μ n 2 ( z 2 − z 1 ) − sin ( 2 λ p 2 − μ n 2 ( z 2 − z 1 ) ) ) × J β m 2 ( μ n r 1 ) − J β m 2 ( μ n r 2 ) π 2 μ n 2 J β m 2 ( μ n r 2 ) κ 1 2 − 4 ( κ 2 λ p 2 − κ 3 ) × ∫ ϕ 1 ϕ 2 ∫ r 1 r 2 ∫ z 1 z 2 f ( φ , ρ , ξ ) Φ m ( φ ) R n ( ρ ) Z p ( ξ ) ρ d ξ d ρ d φ (115)

Proof of Result 2 is omitted for the sake of brevity, as the solution process is quite similar to the one leading to Result 1, while being simpler. In cylindrical coordinates, non-homogeneous equations, non-homogeneous boundary conditions, as well as Neumann and Robin conditions, can be handled in the same way as in spherical coordinates.

This article has shown how to obtain a closed form solution to the non-stationary version of one of the most general second order linear PDEs of mathematical physics, namely the telegrapher’s equation, in a three-dimensional spherical domain, under general boundary conditions and flexible specifications of the ranges of the various parameters. It is quite remarkable that such a sophisticated computational problem admits an exact analytical solution as a uniformly convergent series of triple integrals involving functions whose numerical evaluation is relatively easy. Although most of the problems encountered in the physical sciences are in three dimensions in space, it is worth noticing that it would be a simple extension of this work to derive an exact solution to the same problem with four space variables instead of three, i.e. in a spherical domain defined as follows:

D = { 0 < r 1 ≤ r ≤ r 2 < ∞ , 0 < θ ≤ θ 1 ≤ π , 0 < θ ′ ≤ θ 2 ≤ π , 0 < ϕ 1 ≤ ϕ ≤ ϕ 2 ≤ 2 π , t ≥ 0 } (116)

where θ ′ is an additional angle, i.e. the relations between the spherical coordinates and the rectangular ones, denoted by x 1 , x 2 , x 3 , x 4 are now given by:

x 1 = r sin θ ′ sin θ cos ϕ , x 2 = r sin θ ′ sin θ sin ϕ , x 3 = r sin θ ′ cos θ , x 4 = r cos θ ′ (117)

However, from a practical perspective, one would have to weigh the benefit of having an analytical formula, relative to the cost of increased difficulty of numerical evaluation due to the greater dimension of the integration process.

The author declares no conflicts of interest regarding the publication of this paper.

Guillaume, T. (2020) On the Telegrapher’s Equation with Three Space Variables in Non-Rectangular Coordinates. Journal of Applied Mathematics and Physics, 8, 910-926. https://doi.org/10.4236/jamp.2020.85070