Next Article in Journal
New and Improved Criteria on Fundamental Properties of Solutions of Integro—Delay Differential Equations with Constant Delay
Next Article in Special Issue
Asymptotic Behavior of Solutions of Integral Equations with Homogeneous Kernels
Previous Article in Journal
An Oscillator without Linear Terms: Infinite Equilibria, Chaos, Realization, and Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions

by
Antonella Lupica
1,†,
Clemente Cesarano
2,*,†,
Flavio Crisanti
1,† and
Artur Ishkhanyan
3,4,†
1
Department of Economics, Engineering, Society and Business Organization (DEIm), University of Tuscia, Largo dell’Università, 01100 Viterbo, Italy
2
Section of Mathematics, International Telematic University Uninettuno, 00186 Roma, Italy
3
Department of General Physics, Russian-Armenian University, Yerevan 0051, Armenia
4
Institute for Physical Research, Ashtarak 0203, Armenia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(24), 3316; https://doi.org/10.3390/math9243316
Submission received: 22 November 2021 / Revised: 14 December 2021 / Accepted: 16 December 2021 / Published: 20 December 2021

Abstract

:
We present some solutions of the three-dimensional Laplace equation in terms of linear combinations of generalized hyperogeometric functions in prolate elliptic geometry, which simulates the current tokamak shapes. Such solutions are valid for particular parameter values. The derived solutions are compared with the solutions obtained in the standard toroidal geometry.

1. Introduction

In 2019, the European Union (27 countries) produced a total of around 617.52 Mtoe of electricity for its energy consumption needs, of which 100.63 Mtoe concern solid fossil fuels [1]. To produce energy of this magnitude, large quantities of fuel are required, at great expense. Fusion energy can provide a viable alternative energy source to fossil fuels.
One of the biggest obstacles to fusion is the question of how to hold the reactants long enough for the energy production to exceed the energy input. It is essential to understand the dynamics and confinement of plasma.
Recently, these types of experiments have been carried out with the help of Tokamak devices, which use a powerful magnetic field to confine the plasma to a toroidal shape. The tokamak is one of several types of magnetic confinement devices developed to produce controlled thermonuclear fusion energy. As of 2021, it is the leading candidate for a practical fusion reactor.
A crucial component in this task is the magnetohydrodynamic equilibrium (MHD) which defines the geometry of the confinement magnetic field. The MHD equations produce a second-order nonlinear differential equation known as the Grad–Shafranov equation [2,3], which is the equilibrium equation in MHD ideal for a two-dimensional plasma (e.g., axial-symmetric toroidal plasma in a tokamak). In axisymmetry the Grad–Shafranov equation is given by:
2 ψ R 2 1 R ψ R + 2 ψ Z 2 = μ 0 R J ϕ ,
where J ϕ is the toroidal plasma current, ψ is the flux function and ( R , θ * , Z ) are the standard cylindrical coordinates.
The literature contains a large amount of scientific articles that address the tokamak balance problem, adopting different approaches.
For example, if the reference system is oblate toroidal, the Grad–Shafranov equation in vacuum admits an analytic solution highlighted in the works [4,5]. In case the source term assumes simple forms, the existence of analytical solutions of the nonlinear Grad–Shafranov equation has been shown [6,7,8,9].
The problem of finding the solution of the Grad–Shafranov equation was also addressed from the numerical point of view, through the development of some predictive equilibrium codes [10,11,12,13,14]. Semi-analytical methods have also been adopted in [15,16].
An interesting approach to analytically solve the Grad–Shafranov equation with non-constant source terms through the technique of separable variables was used in [17,18].
The exact analytical solution of the Grad–Shafranov equation in vacuum has been found only if the reference system has a standard circular shape or has an oblate elliptical toroidal geometry [4,5]. Both of these geometries are unsuitable for current tokamak experiments, which are all based on prolate elliptical geometry. In 2019, Crisanti cite Crisanti observed that, in axisymmetry, the Grad–Shafranov equation of vacuum coincides with the Laplace equation for the toroidal component of the vector potential. Cristanti faces the problem of finding the analytical solution for the Grad–Shafranov equation in vacuum (and therefore also of the Laplace equation) when the reference system is written in toroidal prolate elliptic cap-cyclide coordinates. A detailed study of the geometric and metric properties of these coordinates allows us to elaborate the analytical solution of both equations in terms of the Wangerin functions, the analytical expression of which is, however, not known.
The work of Crisanti [19] opens up the possibility of creating a reconstructive code of equilibrium based on an elongated geometry. Wangerin functions were first evaluated in [4]. However, to use them in a budget code, an independent evaluation and a cross-check will be required [19]. Furthermore, it will be necessary to evaluate the derivatives, to derive the poloidal magnetic field and find the Green’s function for this [19] geometry.
The literature shows us a wide range of scientific articles that address the problem of the analytical solution of the Laplace equation (in its various variants), see for example [20,21,22,23,24].
In this article, starting from the recent work of Crisanti [19], we transform the Laplace equation into the Heun equation, and based on the work of A.M. Ishkhanyan [25], we provide an analytical solution of the Laplace equation for some parameter values. Since the Laplace equation and the Grad–Shafranov equation differ by one sign, the procedure presented in this manuscript can be repeated to find the analytical solution of the Grad–Shafranov equation. In Section 2, we introduce the cap-cyclide geometry and, through some transformations of variables, we pass from the Grad–Shafranov equation to the Heun equation. In Section 3, we define the standard toroidal geometry as a special case of the cap-cyclide geometry. In Section 4 we find the analytical solution of the Laplace equation valid in standard toroidal geometry. In Section 5 we find the analytical solution of the Laplace equation in cap-cyclide coordinates for some parameter values. In Section 6, we present the conclusions.

1.1. Hypergeometric Functions

Definition 1. 
A hypergeometric series is formally defined as a power series
c 0 + c 1 z + c 2 z 2 + = n 0 c n z n ,
in which the ratio of successive coefficients is a rational function of n, that is
c n + 1 c n = A ( n ) B ( n ) ,
with A ( n ) and B ( n ) polynomials in n.
It is customary to assume c 0 to be 1. The polynomials A ( n ) and B ( n ) can be factored into linear factors of the form ( a j + n ) and ( b k + n ) respectively, where the a j and b k are complex numbers. For historical reasons, it is assumed that ( 1 + n ) is a factor of B ( n ) .
The ratio between consecutive coefficients now has the form
c n + 1 c n = c ( a 1 + n ) ( a p + n ) d ( b 1 + n ) ( b q + n ) ( 1 + n ) ,
where c and d are the leading terms of A ( n ) and B ( n ) . The series then has the form
1 + a 1 a p b 1 b q · 1 c z d + a 1 a p b 1 b q · 1 ( a 1 + 1 ) ( a p + 1 ) ( b 1 + 1 ) ( b q + 1 ) · 2 c z d 2 + ,
or, by scaling z by the appropriate factor and rearranging
1 + a 1 a p b 1 b q z 1 ! + a 1 ( a 1 + 1 ) a p ( a p + 1 ) b 1 ( b 1 + 1 ) b q ( b q + 1 ) z 2 2 ! +
This series is usually denoted by
p F q ( a 1 , , a p ; b 1 , , b q ; z ) .
The series, if convergent, defines a generalized hypergeometric function, which may then be defined over a wider domain of the argument by analytic continuation.
Historically, the most important are the functions of the form 2 F 1 ( a 1 , a 2 ; b 1 ; z ) . These are sometimes called Gauss’s hypergeometric functions.

1.2. Heun Functions

Definition 2. 
The local Heun function H ( a , Q ; α , β , γ , δ ; z ) (Karl L. W. Heun 1889) is the solution of Heun’s differential equation, that is
d 2 H d z 2 + γ z + δ z 1 + ε z a d H d z + α β z Q z ( z 1 ) ( z a ) H = 0 ,
where ε is a real number such that the Fuchsian condition
γ + δ + ε = α + β + 1 ,
is satisfied. This relation is needed to ensure regularity of the point at ∞.
Function H ( a , Q ; α , β , γ , δ ; z ) is holomorphic and such that H ( a , Q ; α , β , γ , δ ; 0 ) = 1 .
The local Heun function is called a Heun function if it is regular at z = 1 , and is called a Heun polynomial if it is regular at all three finite singular points z = 0 , 1 , a .
The complex number Q is called the accessory parameter. Heun’s equation has four regular singular points: 0 , 1 , a and with exponents ( 0 , 1 γ ) ,   ( 0 , 1 δ ) ,   ( 0 , 1 ε ) , and ( α , β ) .

2. From Grad–Shafranov Equation to Heun Equation

Starting from [19], the vacuum Grad–Shafranov equation are tackled in the elliptical prolate toroidal cap-cyclide coordinates framework.
In axisymmetric cylindrical coordinates ( R , ϑ * , Z ) , ϑ * = 0 , the vacuum Grad–Shafranov equation is given by:
2 ψ R 2 1 R ψ R + 2 ψ Z 2 = 0 .
This equation is formally similar to the Laplace equation:
2 ψ R 2 + 1 R ψ R + 2 ψ Z 2 = 0 .
(Note that Equations (9) and (10) differ for a sign in the first derivative part).
Consequently, it is obvious that any analytical solution found for the Laplace equation will be correlated to the analytical solution of the Grad–Shafranov equation.
For the present toroidal elongated tokamak, the most appropriate coordinate system is the cap-cyclide one [26], that is
x = Λ a s Γ s n ( μ , k ) d n ( ν , k 1 ) cos ϕ , y = Λ a s Γ s n ( μ , k ) d n ( ν , k 1 ) sin ϕ , z = k Π 2 a s Γ ,
where ( x , y , z ) are the standard Cartesian coordinates, d n ,   c n , and s n are the Jacobi elliptic functions, a s is a dimensional scale parameter and
Λ = 1 d n 2 ( μ , k ) s n 2 ( ν , k 1 ) , Γ = s n 2 ( μ , k ) d n 2 ( ν , k 1 ) + Λ k + c n ( μ , k ) d n ( μ , k ) s n ( ν , k 1 ) c n ( ν , k 1 ) 2 , Π = Λ 2 k s n 2 ( μ , k ) d n 2 ( ν , k 1 ) + c n 2 ( μ , k ) d n 2 ( μ , k ) s n 2 ( ν , k 1 ) c n 2 ( ν , k 1 ) .
The variation range of new variables ( μ , ν , ϕ ) is defined as
0 μ K ,   0 ν K ,   0 ϕ 2 π .
Here K and i K are respectively the real and the imaginary complete elliptic integrals
K ( k ) = 0 π 2 d θ 1 k 2 sin 2 θ , i K ( k 1 ) = i 0 π 2 d θ 1 k 1 2 sin 2 θ , k 2 + k 1 2 = 1 ,
where k and k 1 are respectively the parameter and the complementary parameter of the elliptic integrals.
By varying the parameter k, the coordinate transformation (11) describes a large set of quite different geometries [19]. For μ 0 , independently of k, the geometry resembles the standard toroidal geometry; for larger values of μ the shape of the constant μ surfaces depend on the value of the k parameter. For k 0 the surfaces tends to a bean shape. For intermediate values of k the surfaces can be either D or purely elliptical prolate shaped. For k 1 all the surfaces are similar to the standard toroidal ones, independently of the value of μ [19].
As it is proved in [26] and subsequently reported by Crisanti [19], the Laplace Equation (10) in the cap-cyclide coordinates admits a quasi-separable solution of the type
ψ ( μ , ν , ϕ ) = Γ Λ M ( μ ) N ( ν ) Φ ( ϕ ) ,
where the functions M ( μ ) ,   N ( ν ) ,   Φ ( ϕ ) satisfy the following ordinary differental equations
d 2 M d μ 2 + c n ( μ , k ) d n ( μ , k ) s n ( μ , k ) d M d μ + k 2 s n 2 ( μ , k ) α 2 α 3 k 2 s n 2 ( μ , k ) + 1 s n 2 ( μ , k ) M = 0 , d 2 N d ν 2 k 2 c n ( ν , k ) s n ( ν , k ) d n ( ν , k ) d N d ν + d n 2 ( ν , k ) + α 2 k 2 + α 3 d n 2 ( ν , k ) + k 2 d n 2 ( ν , k ) N = 0 , d 2 Φ d ϕ 2 + α 3 Φ = 0 ,
with α 2 = p 2 and α 3 = q 2 , p and q being constants, not necessarily integers [26].
In the rest of the paper, we refer to the Equation (14) as the Laplace equations.
The analytic solution of (14) 3 is obtained immediately and is written as:
Φ ( ϕ ) = c 1 cos α 3 ϕ + c 2 sin α 3 ϕ ,   c 1 c 2 R .
From [19], by substituting z 1 = s n 2 ( μ , k ) and z 2 = d n 2 ( ν , k ) , ( 0 z 1 , z 2 1 ) the two equations for M ( μ ) and N ( ν ) can be written as
d 2 Z d z 2 + 1 2 1 z 1 + 1 z a + 2 z d Z d z + 1 4 A 0 + A 1 z + A 2 z 2 ( z 1 ) ( z a ) z 2 Z = 0 ,
where
( z , a , A 0 , A 1 , A 2 ) = z 1 , 1 k 2 , q 2 k 2 , p 2 k 2 , 1 q 2 ,   for   equation   in   M ( μ ) ( z , a , A 0 , A 1 , A 2 ) = z 2 , k 2 , q 2 k 2 , p 2 k 2 , 1 q 2 ,   for   equation   in   N ( ν ) ,
Following [25] and applying the transformation
Z = z σ u ( z ) , σ = ± A 0 4 a = ± q 2 ,
Equation (16) is reduced to the general Heun equation [25,27]
d 2 u d z 2 + γ z + δ z 1 + ε z a d u d z + α β z Q z ( z 1 ) ( z a ) u = 0 ,
where
δ = ε = 1 2 ,   γ = 1 + 2 σ ,   α β = σ + a A 2 A 0 4 a ,   Q = σ a + 1 2 A 0 + a A 0 + a A 1 4 a ,
and under Fuchsian condition (8)
α + β = γ .
Remark 1. 
Knowing the solution u ( z ) of Heun’s Equation (19), it is consequential to go back to the solution of Wangerin’s Equation (16) given by
Z ( z ) = z σ u ( z ) ,
so that, the solution of the Equation (14) 1 , 2 in terms of μ and ν is given by
M ( μ ) = sn 2 ( μ , k ) σ u ( sn 2 ( μ , k ) ) ,   N ( ν ) = dn 2 ( ν , k 1 ) σ u ( dn 2 ( ν , k 1 ) )

3. Standard Toroidal Geometry as a Particular Case of Cap-Cyclide Geometry

As already pointed out in [19,26], if k = 1 and μ 0 , the cap-cyclide geometry reduces to the standard toroidal one. Toroidal coordinates are a three-dimensional orthogonal coordinate system that results from rotating the two-dimensional bipolar coordinate system about the axis that separates its two foci. Thus, the two foci F 1 and F 2 in bipolar coordinates become a ring of radius a in the x y plane of the toroidal coordinate system; the y-axis is the axis of rotation. The focal ring is also known as the reference circle. The toroidal coordinates are given by the following expressions [26]:
x = a sinh μ cosh μ cos ν cos ϕ , y = a sinh μ cosh μ cos ν sin ϕ , z = a sin ν cosh μ cos ν ,
together with sign ( z ) = sign ( ° ) and a R . The ν coordinate of a point P equals the angle F 1 P F 2 and the μ coordinate equals the natural logarithm of the ratio of the distances d 1 and d 2 to opposite sides of the focal ring
μ = ln d 1 d 2 .
The coordinate ranges are π < ν π , μ 0 and 0 ϕ < 2 π . Both geometries (standard toroidal and cap-cyclide) are obtained by rotating a two-dimensional coordinate system around the y axis, being respectively:
If k = 1 , then the two-dimensional coordinate system of cap-cyclide geometry becomes
x c c k = 1 = cosh ( μ ) sinh ( μ ) a s ( cosh ( 2 μ ) + sin ( 2 ν ) ) , z c c k = 1 = cos ( 2 ν ) 2 a s ( cosh ( 2 μ ) + sin ( 2 ν ) ) .
Moreover, when μ is close to zero, then the cap-cyclide bi-dimensional system approaches standard toroidal geometry, as in the Figure 1.
We emphasize that even if in the limit for μ 0 and k = 1 , the two geometries have the same shape, the metrics are different. In fact, the metric of the standard toroidal geometry is not obtained as a particular case of that of the cap-cyclide geometry (as also can be seen in [26]). This can also be noted in the formulas in Table 1: in fact even when μ 0 , x s t x c c and z s t z c c .
In [26], the authors obtained the analytic solution of Laplace equation in standard toroidal coordinates as
ψ ( μ , ν , ϕ ) = cosh μ cos ν H ( μ ) Θ ( ν ) Φ ( ϕ ) ,
where
H ( μ ) = A P p 1 2 q ( cosh μ ) + B Q p 1 2 q ( cosh μ ) , Θ ( ν ) = A sin ( p ν ) + B cos ( p ν ) , Φ ( ϕ ) = A sin ( p ϕ ) + B cos ( p ϕ )
where P p 1 2 q ,   Q p 1 2 q are Legendre functions and A and B are constants.

4. Analytic Solution of Laplace Equation (14) 1 in Standard Toroidal Geometry

In this section, we determine the analytic form of the solution of Equation (14) 1 valid in standard toroidal geometry in terms of hypergeometric functions.
In the case of standard toroidal geometry, since k = 1 , then a = 1 and Equation (16) becomes
d 2 Z d z 2 + 1 z 1 + 1 z d Z d z + 1 4 A 0 + A 1 z + A 2 z 2 ( z 1 ) 2 z 2 Z = 0 ,
with
( z , a , A 0 , A 1 , A 2 ) = z 1 , 1 , q 2 , p 2 , 1 q 2 ,   for   equation   in   M ( μ ) , ( z , a , A 0 , A 1 , A 2 ) = z 2 , 1 , q 2 , p 2 , 1 q 2 ,   for   equation   in   N ( ν ) .
After substitution (18), Equation (25) becomes
u ( z ) + γ z + δ + ε z 1 u ( z ) + α β z Q z ( z 1 ) 2 u ( z ) = 0 ,
where
δ = ε = 1 2 ,   γ = 1 + 2 σ ,   α β = σ + 1 4 ,   Q = σ + 2 q 2 A 1 4 .
The last equation is the Heun Equation (19) for k = 1 . Let
c = 1 2 γ α β ( γ α β ) 2 4 α β + 4 Q .
Since, in our case, Fuchsian condition reduces to α + β = γ , then
c = Q α β .
Following [28], consider the change of the dependent variable u ( z ) = ( z 1 ) c G ( z ) . After substitution, from Equation (27), we get a hypergeometric equation
z ( z 1 ) G ( z ) + ( 2 c + α + β + 1 ) z γ G ( z ) + ( c + α ) ( c + β ) G ( z ) = 0 ,
whose general solution is given by
G ( z ) = c 1   2 F 1 c + α , c + β ; γ ; z + + c 2 ( 1 ) 1 γ z 1 γ   2 F 1 1 + c + α γ , 1 + c + β γ ; 2 γ ; z ,
with c 1 , c 2 C . Thus, the analytic solution of the Heun Equation (27) is given as
u ( z ) = ( z 1 ) c [ c 1   2 F 1 c + α , c + β ; γ ; z + + c 2 ( 1 ) 1 γ z 1 γ   2 F 1 1 + c + α γ , 1 + c + β γ ; 2 γ ; z ] ,
with c 1 , c 2 C . Accordingly, the analytic solution of Equation (14) 1 for k = 1 reads
M ( μ ) = sn 2 ( μ ) σ sn 2 ( μ ) 1 c [ c 1   2 F 1 c + α , c + β ; γ ; sn 2 ( μ ) + c 2 ( 1 ) 1 γ sn 2 ( μ ) 1 γ   2 F 1 1 + c + α γ , 1 + c + β γ ; 2 γ ; sn 2 ( μ ) ] ,
with c 1 , c 2 C . Similarly, the analytic solution of Equation (14) 2 for k = 1 reads
N ( ν ) = dn 2 ( ν ) σ dn 2 ( ν ) 1 c [ c 1   2 F 1 c + α , c + β ; γ ; dn 2 ( ν ) + c 2 ( 1 ) 1 γ dn 2 ( ν ) 1 γ   2 F 1 1 + c + α γ , 1 + c + β γ ; 2 γ ; dn 2 ( ν ) ] .

5. Analytic Solutions of Laplace Equation (14) 1 in Cap-Cyclide Geometry

In this section, we determine the analytical form of the solution of Equation (14) 1 valid in the cap-cyclide geometry in terms of hypergeometric functions and for particular values of the parameters.
Let’s solve the Heun Equation (19) first.
Following [25], we impose that the following series developed around the singularity z = a :
u = z a 1 a μ n = 0 + c n z a 1 a n ,   μ = 0 , 1 ε
is a solution of the Heun equation and that it reduces to a generalized hypergeometric function.
Let q be integer. This condition will allow us to construct a solution in terms of hypergeometric functions.
Let μ = 0 . The generalized hypergeometric function r F s is defined as the series [29]
r F s a 1 , , a r ; b 1 , , b s ; z a 1 a = n = 0 c n z a 1 a n ,
where the coefficients satisfy the two-term recurrence relation
c n c n 1 = 1 n k = 1 r ( a k 1 + n ) k = 1 s ( b k 1 + n ) .
Since q is integer, then γ is also integer. Assume that γ is negative, i.e., γ = N ,   N = 1 , 2 , 3 , . Following [25], we set r = N + 2 , s = N + 1 , and
a 1 , , a N , a N + 1 , a N + 2 = 1 + e 1 , , 1 + e N , α , β ,
b 1 , , b N , b N + 1 = e 1 , , e N , ε ,
with e 1 , , e N being certain constants.
It can be shown that the Heun Equation (19) admits the following solution [25]
u = N + 2 F 1 + N 1 + e 1 , , 1 + e N , α , β ; e 1 , , e N , ε ; z a 1 a .
if and only if the following polynomial Π ( n ) in an auxiliary variable n is identically zero [25]:
Π ( n ) = R n ( α 1 + n ) ( β 1 + n ) n ( ε 1 + n ) k = 1 N ( e k + n ) + Q n 1 k = 1 N ( e k 1 + n ) + + P n 2 ( n 1 ) ( ε 2 + n ) ( α 2 + n ) ( β 2 + n ) k = 1 N ( e k 2 + n ) ,
where
R n = a n ( ε 1 + n ) , Q n = n ( a 1 ) ( n + α + β ) + a ( n 1 ) + n ( δ + a ε ) + a α β Q , P n = ( 1 a ) ( α + n ) ( β + n ) .
The term of degree N + 2 in Π ( n ) is automatically canceled, while the term of degree N + 1 is zero owing to the Fuchsian condition. Then, by canceling the coefficients of the polynomial Π ( n ) , we obtain a system of N + 1 algebraic equations, N of which are used to determine the N constants e 1 , , e N and the last equetion imposes a restriction on the accessory parameter Q: this restriction is obtained by substituting the expressions of the constants e 1 , , e N in the last equation thus arriving at a polynomial equation of degree N + 1 in the variable Q.
Therefore, there are many solutions of Equation (19) in terms of a single generalized hypergeometric function r F s if γ is a negative interger and the accessory parameter Q satisfies a certain polynomial equation [25].
Similarly, another solution can be constructed in terms of a power series in the vicinity of the singularity z = 1 :
u = z 1 a 1 μ n = 0 + c n z 1 a 1 n ,   μ = 0 , 1 δ .
Let μ = 0 . Such a series can be reduced to a generalized hypergeometric series for any negative integer γ = N ,   N = 1 , 2 , 3 , [25], i.e.,
u = N + 2 F 1 + N 1 + s 1 , , 1 + s N , α , β ; s 1 , , s N , δ ; z 1 a 1 .
The accessory parameter Q and the parameters s 1 , s 2 , , s N involved in solution (39) are determined from a system of N + 1 algebraic equations constructed by equating the coefficients of the following polynomial Π ˜ ( n ) in an auxiliary variable n to zero [25]:
Π ˜ ( n ) = R ˜ n ( α 1 + n ) ( β 1 + n ) n ( δ 1 + n ) k = 1 N ( e k + n ) + Q ˜ n 1 k = 1 N ( e k 1 + n ) + + P ˜ n 2 ( n 1 ) ( δ 2 + n ) ( α 2 + n ) ( β 2 + n ) k = 1 N ( e k 2 + n ) ,
where
R ˜ n = n ( δ 1 + n ) , Q ˜ n = n ( 1 a ) ( n + α + β ) + n 1 + n ( δ + a ε ) + α β Q , P ˜ n = ( a 1 ) ( α + n ) ( β + n ) .
The restrictions imposed on Q obtained from the polynomial Π ( n ) (36) and then from the polynomial Π ˜ ( n ) (40) are identical.
The two solutions (35) and (39) are linearly independent, indeed their Wronskian turns out to be:
W ( z ) = α β ( a 1 ) e 1 e 2 e N s 1 s 2 s N δ ε [ ( 1 + e 1 ) ( 1 + e 2 ) ( 1 + e N ) s 1 s 2 s N δ · · N + 2 F N + 1 2 + e 1 , 2 + e 2 , , 2 + e N , 1 + α , 1 + β ; 1 + e 1 , 1 + e 2 , , 1 + e N , 1 + ε ; z a 1 a · · N + 2 F N + 1 1 + s 1 , 1 + s 2 , , 1 + s N , α , β ; s 1 , s 2 , , s N , δ ; z 1 a 1 + + e 1 e 2 e N ( 1 + s 1 ) ( 1 + s 2 ) ( 1 + s N ) ε · · N + 2 F N + 1 1 + e 1 , 1 + e 2 , , 1 + e N , α , β ; e 1 , e 2 , , e N , ε ; z a 1 a · · N + 2 F N + 1 2 + s 1 , 2 + s 2 , , 2 + s N , 1 + α , 1 + β ; 1 + s 1 , 1 + s 2 , , 1 + s N , 1 + δ ; z 1 a 1 ] 0
Thus, if γ = N , N Z and Q satisfies a certain polynomial equation of N + 1 degree, the general solution of the Heun Equation (19) is given by:
u ( z ) = c 1   N + 2 F N + 1 1 + e 1 , , 1 + e N , α , β ; e 1 , , e N , ε ; z a 1 a + + c 2   N + 2 F N + 1 1 + s 1 , , 1 + s N , α , β ; s 1 , , s N , δ ; z 1 a 1 ,   c 1 , c 2 C .
Similarly, the general analytic solution of Laplace equation is given by:
M ( μ ) = sn 2 σ ( μ , k ) [ c 1   N + 2 F N + 1 1 + e 1 , , 1 + e N , α , β ; e 1 , , e N , ε ; sn 2 ( μ , k ) a 1 a + + c 2   N + 2 F N + 1 1 + s 1 , , 1 + s N , α , β ; s 1 , , s N , δ ; sn 2 ( μ , k ) 1 a 1 ] ,   c 1 , c 2 C .
Now we write hypergeometric solution of the Heun equation for some values of parameter γ .

5.1. Case: N = 0     γ = 0

In this case, the Heun equation becomes
d 2 u d z 2 + δ z 1 + ε z a d u d z + α β z Q z ( z 1 ) ( z a ) u = 0 ,
and admits a solution in terms of the ordinary hypergeometric function:
u = c 1   2 F 1 α , β ; ε ; z a 1 a + c 2   2 F 1 α , β ; δ ; z 1 a 1 ,   c 1 , c 2 C .
Since the Fuchsian condition (8) reduces to α + β = 0 , from Equations (36) and (40), we obtain the restriction on Q
Q = 0 .
  • σ = q 2
    From γ = 0 , we get q = 1 and Fuchsian condition (8) α = β . Furthermore, from Equation (17) 1 we obtain
    a = 1 k 2 ,   A 0 = 1 k 2 ,   A 1 = p 2 k 2 ,   A 2 = 0 ,   for   equation   in   M ( μ ) .
    Since α β = σ + a A 2 A 0 4 a , then from Equation (47), we get
    α 1 , 2 = ± 1 2 ,   β 1 , 2 = 1 2 .
    The expression Q = σ a + 1 2 A 0 + a A 0 + a A 1 4 a , together with Equation (46), gives
    p 2 4 k 2 = 0 .
    Thus we get p 2 = 0 . Consequently, Heun equation reduces to
    u ( z ) + 1 2 ( z 1 ) + 1 2 ( z a ) u ( z ) z 4 z ( z 1 ) ( z a ) = 0 ,
    whose solution,   k is given by
    u = c 1   2 F 1 1 2 , 1 2 ; 1 2 ; z a 1 a + c 2   2 F 1 1 2 , 1 2 ; 1 2 ; z 1 a 1     z [ 0 , 1 ] .
    Hence, the solution of the corresponding Wangerin Equation (16)   k is given by
    Z ( z ) = 1 z [ c 1   2 F 1 1 2 , 1 2 ; 1 2 ; z a 1 a + + c 2   2 F 1 1 2 , 1 2 ; 1 2 ; z 1 a 1 ]     z ( 0 , 1 ] .
    The solution of the corresponding (14) 1 equation   0 < k < 1 is given by
    M ( μ ) = 1 sn 2 ( μ , k ) [ c 1   2 F 1 1 2 , 1 2 ; 1 2 ; sn 2 ( μ , k ) a 1 a + + c 2   2 F 1 1 2 , 1 2 ; 1 2 ; sn 2 ( μ , k ) 1 a 1 ] ,
      μ ( 0 , K ] .
  • σ = q 2   The same solution as for q = 1 !

5.2. Case N = 1 , γ = 1

In this case, the Heun equation becomes
d 2 u d z 2 + 1 z + δ z 1 + ε z a d u d z + α β z Q z ( z 1 ) ( z a ) u = 0 ,
and admits a solution in terms of the generalized hypergeometric function:
u = c 1   3 F 2 1 + e 1 , α , β ; e 1 , ε ; z a 1 a + c 2   3 F 2 1 + s 1 , α , β ; s 1 , δ ; z 1 a 1 ,   c 1 , c 2 C ,
with constants e 1 , s 1 . Since the Fuchsian condition (8) reduces to α + β = 1 , the polynomials (36) and (40) give
e 1 = a α β Q ,   s 1 = α β Q .
Besides, the accessory parameter Q satisfies the following polinomial equation of second degree
2 Q 2 1 + a Q + 2 a α β = 0 .
  • σ = q 2
    From γ = 1 , we get q = 2 and the Fuchsian condition (8) reads α = β 1 . Besides, from Equation (17) we obtain
    a = 1 k 2 ,   A 0 = 4 k 2 ,   A 1 = p 2 k 2 ,   A 2 = 3 ,   for   equation   in   M ( μ )
    Since α β = σ + a A 2 A 0 4 a , then from Equation (56), we get
    α 1 = 3 2 ,   β 1 = 1 2 , α 2 = 1 2 ,   β 2 = 3 2
    From Equation (55), we get the Q values
    Q 1 , 2 = 1 4 1 + 1 k 2 ± 1 + 14 k 2 + 1 k 4 .
    The expression Q = σ a + 1 2 A 0 + a A 0 + a A 1 4 a gives the relation between k and p
    1 + 1 ± 1 + 1 k 4 + 14 k 2 k 2 + p 2 4 k 2 = 0 .
    The Heun equation reduces to
    u ( z ) + 1 z + 1 2 ( z 1 ) + 1 2 ( z a ) u ( z ) 3 4 z + Q z ( z 1 ) ( z a ) = 0
    whose solution,   k is given by
    u = c 1   3 F 2 1 + e 1 , 3 2 , 1 2 ; e 1 , 1 2 ; z a 1 a + c 2   3 F 2 1 + s 1 , 3 2 , 1 2 ; s 1 1 2 ; z 1 a 1 ,
      c 1 , c 2 C and z [ 0 , 1 ] .
    Hence, the solution of the corresponding Wangerin equation for any k is given by
    Z ( z ) = 1 z [ c 1   3 F 2 1 + e 1 , 3 2 , 1 2 ; e 1 , 1 2 ; z a 1 a + + c 2   3 F 2 1 + s 1 , 3 2 , 1 2 ; s 1 1 2 ; z 1 a 1 ] ,
      c 1 , c 2 C and z ( 0 , 1 ] .
    Consequently, the solution of the corresponding (14) 1 is written as
    M ( μ ) = 1 sn 2 ( μ , k ) [ c 1   3 F 2 3 2 , 1 2 , 1 + e 1 ; 1 2 , e 1 ; sn 2 ( μ , k ) a 1 a + + c 2   3 F 2 3 2 , 1 2 , 1 + s 1 ; 1 2 , s 1 ; sn 2 ( μ , k ) 1 a 1 ] ,
    c 1 , c 2 C and μ ( 0 , K ] .
  • σ = q 2   The same solution as for q = 2 .

5.3. Positive Integer γ

Let γ now be a positive integer. Following [25], we apply a change of the dependent variable in order to obtain a Heun equation with altered parameters and with a negative characteristic exponent γ . Let u = z 1 γ w : this change transforms the Heun equation into another Heun equation with the altered parameter γ * = 2 γ , i.e.,
w ( z ) + γ * z + δ z 1 + ε z a w ( z ) + α * β * z Q * z ( z 1 ) ( z a ) w ( z ) = 0 ,
where α * = ( α + 1 γ ) ,   β * = ( β + 1 γ ) and Q * = Q + ( 1 γ ) ( ε + a δ ) .
Then, for γ   2 , we have a Heun equation with a zero or negative integer γ * = N , N = 1 , 2 , 3 , . . As a result, we obtain the solution
u = z 1 γ [ c 1   N + 2 F 1 + N e 1 + 1 , , e N + 1 , α + 1 γ , β + 1 γ ; e 1 , , e N , ε ; z a 1 a + + c 2   N + 2 F 1 + N s 1 + 1 , , s N + 1 , α + 1 γ , β + 1 γ ; s 1 , , s N , δ ; z 1 a 1 ] .
The only exception is the case γ = 1 : we do not know an r F s solution in this exceptional case [25].

γ = 2

In this case, the Heun equation becomes
d 2 u d z 2 + 2 z + δ z 1 + ε z a d u d z + α β z Q z ( z 1 ) ( z a ) u = 0 .
  • Let σ = q 2
    From γ = 2 , we have q = 1 and from Fuchsian condition (8) we have α = 2 β . Besides, from Equation (17) we obtain
    a = 1 k 2 ,   A 0 = 1 k 2 ,   A 1 = p 2 k 2 ,   A 2 = 0 ,   for   equation   in   M ( μ )
    Since α β = σ + a A 2 A 0 4 a , then from Equation (66), we get
    α 1 = 1 2 ,   α 2 = 3 2 , β 1 = 3 2 ,   β 2 = 1 2 ,
    and expression Q = σ a + 1 2 A 0 + a A 0 + a A 1 4 a gives
    Q = 2 + 2 k 2 + p 2 4 k 2 .
    Consider the change u = w z . We obtain another Heun equation that reads as:
    d 2 w d z 2 + γ * z + δ z 1 + ε z a d w d z + α * β * z Q * z ( z 1 ) ( z a ) w = 0 ,
    with γ * = 0 , α * = α 1 ,   β * = β 1 and Q * = Q ( a δ + ε ) .
    We already know the solution of Equation (67) which for any k is written as
    w ( z ) = c 1   2 F 1 α * , β * ; ε ; z a 1 a + c 2   2 F 1 α * , β * ; δ ; z 1 a 1 ,   α * = 1 2 ,   β * = 1 2 ,
    with c 1 , c 2 C , ( p * ) 2 = 0 and Q * = 0 . Consequently, Q = a δ + ε and we get the relation between k and p
    p 2 4 k 2 = 0 .
    Thus, the solution of the Heun Equation (65) for any k is given by
    u ( z ) = 1 z c 1   2 F 1 α * , β * ; ε ; z a 1 a + c 2   2 F 1 α * , β * ; δ ; z 1 a 1 ,
    with α * = 1 2 ,   β * = 1 2 ,   c 1 , c 2 C and z ( 0 , 1 ] .
    Hence, the solution of the corresponding Wangerin equation for any k is given by
    Z ( z ) = 1 z c 1   2 F 1 α * , β * ; ε ; z a 1 a + c 2   2 F 1 α * , β * ; δ ; z 1 a 1 ,
    with α * = 1 2 ,   β * = 1 2 ,   c 1 , c 2 C and z ( 0 , 1 ] .
    The solution of the corresponding Equation (14) 1 for any 0 < k < 1 is given by
    M ( μ ) = 1 sn 2 ( μ , k ) [ c 1   2 F 1 α * , β * ; ε ; sn 2 ( μ , k ) a 1 a + + c 2   2 F 1 α * , β * ; δ ; sn 2 ( μ , k ) 1 a 1 ] ,
    with α * = 1 2 ,   β * = 1 2 ,   c 1 , c 2 C and μ ( 0 , K ]
  • σ = q 2   The same solution as for q = 1 .

6. Conclusions and Discussion

In this paper, we have presented the analytic solution of Laplace’s Equation (14) 1 in terms of linear combinations of generalized hypergeometric functions which hold for integer values of γ and when the accessory parameter Q satisfies a particular polynomial equation, the degree of which is related to the value of the integer γ . γ being an integer is not a particularly significant limitation, since γ = 2 σ + 1 = q + 1 and q is generally an integer.
The second Laplace Equation (14) 2 still remains to be solved: it too is reduced to a Heun equation (but with different parameters) and can therefore be solved with the technique applied here.
Once this is done, the complete solution of the three-dimensional Laplace Equation is constructed.
At this respect, since the Laplace equation and the Grad–Shafranov equation differ by one sign (plus instead of minus in the term of the first derivative), one can also apply the technique presented in this manuscript to the Grad–Shafranov equation to construct its solution in terms of a linear combination of hypergeometric functions. We will investigate these latter aspects in future work.

Author Contributions

Conceptualization, A.L., C.C., F.C. and A.I.; Data curation, A.L., C.C., F.C. and A.I.; Formal analysis, A.L., C.C., F.C. and A.I.; Investigation, A.L., C.C., F.C. and A.I.; Methodology, C.C., F.C. and A.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The work by Artur Ishkhanyan was supported by the Russian-Armenian University, the Armenian Science Committee (SC Grant No. 20RF-171, No. 21SC-BRFFR-1C021 and No. 21AG-1C064) and the Armenian National Science and Education Fund (ANSEF Grant No. PS-2520). The work of Antonella Lupica is supported by GNFM-INdAM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Commission, E. Europe Energy Information Administration. 2021. Available online: https://ec.europa.eu/energy/home_en (accessed on 12 November 2021).
  2. Grad, H.; Rubin, H. Hydromagnetic equilibria and force-free fields. J. Nucl. Energy 1958, 7, 284–285. [Google Scholar] [CrossRef]
  3. Shafranov, V. On magnetohydrodynamical equilibrium configurations. Sov. Phys. JETP 1958, 6, 1013. [Google Scholar]
  4. Aikawa, K.; Takahara, M. Wangerin Functions; Report of Faculty of Engineering; Yamanashi University: Yamanashi, Japan, 1975; pp. 126–133. [Google Scholar]
  5. Honma, T.; Kito, M.; Kaji, I.; Seki, M. The Vacuum Poloidal Flux Functions Satisfying of the Grad–Shafranov Equation in the Flat-Ring Cyclide Coordinate System. Hokkaido Daigaku Kogakubu Kenkyu Hokoku 1979, 94, 123–131. [Google Scholar]
  6. Solov’ev, L. The theory of hydromagnetic stability of toroidal plasma configurations. Sov. Phys. JETP 1968, 26, 400–407. [Google Scholar]
  7. Greene, J. An analytic large aspect ratio, high beta equilibrium. Plasma Phys. Control. Fusion 1988, 30, 327. [Google Scholar] [CrossRef]
  8. Zheng, S.; Wootton, A.; Solano, E.R. Analytical tokamak equilibrium for shaped plasmas. Phys. Plasmas 1996, 3, 1176–1178. [Google Scholar] [CrossRef]
  9. Cerfon, A.J.; Freidberg, J.P. “One size fits all” analytic solutions to the Grad–Shafranov equation. Phys. Plasmas 2010, 17, 032502. [Google Scholar] [CrossRef]
  10. Feneberg, W.; Lackner, K. Multipole tokamak equilibria. Nucl. Fusion 1973, 13, 549. [Google Scholar] [CrossRef]
  11. Cenacchi, G.; Galvao, R.; Taroni, A. Numerical computation of axisymmetric MHD-equilibria without conducting shell. Nucl. Fusion 1976, 16, 457. [Google Scholar] [CrossRef]
  12. Blum, J.; Le Foll, J.; Thooris, B. The self-consistent equilibrium and diffusion code SCED. Comput. Phys. Commun. 1981, 24, 235–254. [Google Scholar] [CrossRef]
  13. Albanese, R.; Blum, J.; De Barbieri, O. Numerical studies of the Next European Torus via the PROTEUS code. In Proceedings of the 12th Conference on Numerical Simulation of Plasmas, San Francisco, CA, USA, 20–24 September 1987. [Google Scholar]
  14. Albanese, R.; Ambrosino, R.; Mattei, M. CREATE-NL+: A robust control-oriented free boundary dynamic plasma equilibrium solver. Fusion Eng. Des. 2015, 96, 664–667. [Google Scholar] [CrossRef]
  15. Alladio, F.; Crisanti, F. Analysis of MHD equilibria by toroidal multipolar expansions. Nucl. Fusion 1986, 26, 1143. [Google Scholar] [CrossRef]
  16. Alladio, F.; Crisanti, F.; Marinucci, M.; Micozzi, P.; Tanga, A. Analysis of tokamak configurations using the toroidal multipole method. Nucl. Fusion 1991, 31, 739. [Google Scholar] [CrossRef]
  17. Atanasiu, C.; Günter, S.; Lackner, K.; Miron, I. Analytical solutions to the Grad–Shafranov equation. Phys. Plasmas 2004, 11, 3510–3518. [Google Scholar] [CrossRef]
  18. Guazzotto, L.; Freidberg, J.P. A family of analytic equilibrium solutions for the Grad–Shafranov equation. Phys. Plasmas 2007, 14, 112508. [Google Scholar] [CrossRef] [Green Version]
  19. Crisanti, F. Analytical solution of the Grad Shafranov equation in an elliptical prolate geometry. J. Plasma Phys. 2019, 85, 905850210. [Google Scholar] [CrossRef]
  20. Holm, S. Dispersion analysis for wave equations with fractional Laplacian loss operators. Fract. Calc. Appl. Anal. 2019, 22, 1596–1606. [Google Scholar] [CrossRef]
  21. Guariglia, E. Riemann zeta fractional derivative—Functional equation and link with primes. Adv. Differ. Equ. 2019, 2019, 1–15. [Google Scholar] [CrossRef]
  22. Fărcăşeanu, M. On an eigenvalue problem involving the fractional (s, p)-Laplacian. Fract. Calc. Appl. Anal. 2018, 21, 94–103. [Google Scholar] [CrossRef] [Green Version]
  23. Arioua, Y.; Ma, L. On Criteria of Existence for Nonlinear Katugampola Fractional Differential Equations with p–Laplacian Operator. Fract. Differ. Calc. 2021, 1, 55–68. [Google Scholar] [CrossRef]
  24. Guariglia, E. Fractional calculus, zeta functions and Shannon entropy. Open Math. 2021, 19, 87–100. [Google Scholar] [CrossRef]
  25. Ishkhanyan, A.M. Generalized hypergeometric solutions of the Heun equation. Theor. Math. Phys. 2020, 202, 1–10. [Google Scholar] [CrossRef] [Green Version]
  26. Moon, P.; Spencer, D.E. Field Theory Handbook: Including Coordinate Systems, Differential Equations and Their Solutions; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  27. Ronveaux, A.; Arscott, F. Heun’s Differential Equations; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  28. Figueiredo, B. Solutions of Heun’s general equation and elliptic Darboux equation. arXiv 2020, arXiv:2011.07940. [Google Scholar] [CrossRef]
  29. Slater, L.J. Generalized Hypergeometric Functions. 1966. Available online: http://inis.jinr.ru/sl/vol1/UH/_Ready/Mathematics/Advanced%20calculus/Slater%20L.J.%20Generalized%20hypergeometric%20functions%20(Cambridge,%201966)(L)(143s).pdf (accessed on 20 November 2021).
Figure 1. Two-dimensional coordinate system of cap-cyclide geometry with k = 1 and μ close to 0.
Figure 1. Two-dimensional coordinate system of cap-cyclide geometry with k = 1 and μ close to 0.
Mathematics 09 03316 g001
Table 1. Two-dimensional coordinate system of standard toroidal and cap-cyclide geometries.
Table 1. Two-dimensional coordinate system of standard toroidal and cap-cyclide geometries.
Standard ToroidalCap-Cyclide
x s t = a sinh μ cosh μ cos ν , z s t = a sin ν cosh μ cos ν , x c c = Λ a s Γ s n ( μ , k ) d n ( ν , k 1 ) z c c = k Π 2 a s Γ
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lupica, A.; Cesarano, C.; Crisanti, F.; Ishkhanyan, A. Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions. Mathematics 2021, 9, 3316. https://doi.org/10.3390/math9243316

AMA Style

Lupica A, Cesarano C, Crisanti F, Ishkhanyan A. Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions. Mathematics. 2021; 9(24):3316. https://doi.org/10.3390/math9243316

Chicago/Turabian Style

Lupica, Antonella, Clemente Cesarano, Flavio Crisanti, and Artur Ishkhanyan. 2021. "Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions" Mathematics 9, no. 24: 3316. https://doi.org/10.3390/math9243316

APA Style

Lupica, A., Cesarano, C., Crisanti, F., & Ishkhanyan, A. (2021). Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions. Mathematics, 9(24), 3316. https://doi.org/10.3390/math9243316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop