The Lie Group SU(2) Hopf Fibration and the Fourier Equation

Adelin Mulenda Mbuto¹, Lucien Zihindula Biguru², Jean Masudi Kalongama³, Joseph Cimbela Kabongo¹, Albert Kabasele Yenga-Yenga¹

¹National Pedagogical University/Kinshasa, Kinshasa/Binza, Democratic Republic of the Congo.
²Official University of Bukavu/Bukavu, Sud-Kivu/Bukavu, Democratic Republic of the Congo.
³Institut Supérieur Pédagogique/Bukavu, Sud-Kivu/Bukavu, Democratic Republic of the Congo.

DOI: 10.4236/jamp.2020.87105                

Abstract

The Fourier equation explains the dynamics of heat transfer. But bringing this phenomenon closer to the notion of fibration seems difficult to achieve. This study then aims to find the solution of the one-dimensional Fourier equation and to interpret it in terms of bundle. And then apply the results obtained at the Kankule site in Katana in South Kivu. To do this work, we resorted to geometric or topological analysis of the Hopf fibration of the unit sphere S³ (identifiable in SU(2)). We had taken the temperatures of the thermal waters and the soil of Kankule, from 2010 to 2014, in situ. And laboratory analyses had allowed us to know the physical and chemical properties of the soil and water at each of our 14 study sites in Kankule. The data of the geomagnetic field of each site, were taken in on the site NOAA, for our period of study. We then determined the integral curve (geotherm) of the Fourier equation and wrote it as a unit quaternion which is a bundle. The constants intervened in the geotherm, for each site of Kankule, we had obtained them statistically. We have found that the geotherm of each Kankule site is a bundle. We have compared this model to the bundle model of the geomagnetic field. From there we realized that to determine the energy potential of Kankule, we should consider the thermal springs separately. We were able to find a connection between the fibration of the geomagnetic field and the heat field for the Kankule site.

Keywords

Fiber, Hopf Fiber, Geotherm, Quaternion, Geomagnetic Field and Potential Energy

1. Conceptual Fame

1.1. The Quaternions and the Sphere Unit S³

1.1.1. Quaternions

1) Definitions

In ${ℝ}^4$, consider $q=\left(a,b,c,d\right)=\left(a,u\right)$ and $q^{\prime }=\left(a^{\prime },b^{\prime },c^{\prime },d^{\prime }\right)=\left(a^{\prime },v\right)$ ; for $u=\left(b,c,d\right)$ and $v=\left(b^{\prime },c^{\prime },d^{\prime }\right)$ in ${ℝ}^3$. The addition and multiplication in ${ℝ}^4$ are defined by the relations below [1] [2] [3]:

$q+q^{\prime }=\left(a+a^{\prime },u+v\right)$ (1)

with $u+v$ the vector addition in ${ℝ}^3$, and

$\begin{array}{c}q\ast q^{\prime }=(a{a}^{\prime }-b{b}^{\prime }-c{c}^{\prime }-d{d}^{\prime },a{b}^{\prime }+b{b}^{\prime }+c{d}^{\prime }-d{c}^{\prime },\\ a{c}^{\prime }+c{a}^{\prime }-b{d}^{\prime }+d{b}^{\prime },a{d}^{\prime }+d{a}^{\prime }+b{c}^{\prime }-c{b}^{\prime })\\ =\left(a{a}^{\prime }-u\cdot v,a\cdot v+a^{\prime }\cdot v+u\wedge v\right)\end{array}$ (2)

where $u\cdot v$ and $u\wedge v$ are respectively the scalar product and the vector product in ${ℝ}^3.$ We verify that $\left({ℝ}^4,+,\ast \right)$ is a commutative body, called the quaternion body. We mark it H, in honor of Hamilton who built it first.

We have $\bar{q}=\left(a,-u\right)$ as a conjugate of q, and $-q=\left(-a,-u\right)$ being the opposite of q. The norm of q is

$\left|q\right|=\sqrt{q\ast \bar{q}}=\sqrt{a^2+b^2+c^2+d^2}$ (3)

The inverse of q is

$q^{-1}=\frac{\bar{q}}{\left|q\right|}$ (4)

We denote [4] [5] [6]

$1=\left(1,0,0,0\right),i=\left(0,1,0,0\right),j=\left(0,0,1,0\right)\text{and}k=\left(0,0,0,1\right).$ (5)

and

$i^2=j^2=k^2=-1=i\ast j\ast k$ (6)

with $ij=k=-ji$ ; $jk=i=-kj$ and $ki=j=-ik$. H is a vector space on R whose basis is $\left(1,i,j,k\right)$. Thus the quaternion

$q=\left(a,b,c,d\right)=a+b\ast i+c\ast j+d\ast k$ (7)

The set of pure quaternions is

${ℍ}_{pur}=\left\{p\in H/p=\left(0;b;c;d\right);\left(b;c;d\right)\in {ℝ}^3\right\}$ (8)

We identify ${ℍ}_{pur}$ with ${ℝ}^3$. We check that if $p\in {ℍ}_{pur}$, then $p=-\bar{p}$. So for $p\in H$ which is such that $p^2=-1$, we will have $\left|p\right|=1$ and $p=-\bar{p}$, which implies that $p\in {ℍ}_{pur}.$

2) Theorem 1

Any quaternion q can be written as

$q=\left|q\right|\left(\cos \theta +p\sin \theta \right)$ (9)

with $\theta \in R$ and $p\in {ℍ}_{pur}$ [1] [3].

Indeed, $q/\left|q\right|=\left(a;b;c;d\right)$ is a quaternion of module 1, then $a\in \left[-1;1\right]\cap R$ because $\left|q/\left|\bar{q}\right|\right|=\sqrt{a^2+b^2+c^2+d^2}=1$. There then exists $\theta \in R$ such that

$a=\cos \theta$ and $\sqrt{b^2+c^2+d^2}=\sin \theta$. When θ is not a multiple of π, we set $p=\frac{\frac{q}{\left|q\right|}-\cos \theta }{\sin \theta }$, a pure quaternion of norm 1 and whose square is −1.

When θ is a multiple of π, we immediately have the result stated by the theorem, whatever the pure quaternion considered.

1.1.2. The S² and S³ Spheres

A) Definitions

The unit sphere $S^n$, of dimension n, of ${ℝ}^{n+1}$ is defined by [2] [7] [8]:

$S^n=\left\{\left(x_0;x_1;\cdots ;x_n\right)\in {ℝ}^{n+1}/{\displaystyle {\sum }_{i=0}^0{x}_i^2}=1\right\}$ (10)

Thus $S^3$ is the part of ${ℝ}^4$ defined by

$S^3=\left\{\left(a,b,c,d\right)\in {ℝ}^4/a^2+b^2+c^2+d^2=1\right\}$ (11)

and $S^2$ the part of ${ℝ}^3$ defined by

$S^2=\left\{\left(x,y,z\right)\in {ℝ}^3/x^2+y^2+z^2=1\right\}$ (12)

Note $S^{\prime }\left(0,0,-1\right)$ and $N^{\prime }\left(0,0,1\right)$ are respectively the south and north poles of S². Consider the open sets $U_{N^{\prime }}=S^2\backslash \left\{N^{\prime }\right\}$ and $U_{S^{\prime }}=S^2\backslash \left\{S^{\prime }\right\}$. Let the applications, called stereographic projections, be as follows:

${\psi }_{N^{\prime }}:U_{N^{\prime }}\to {ℝ}^2,\text{with}{\psi }_{N^{\prime }}\left(x,y,z\right)=\left(\frac{x}{1-z},\frac{y}{1-z}\right)$ ; (13)

and

${\psi }_{S^{\prime }}:U_{S^{\prime }}\to {ℝ}^2,\text{with}{\psi }_{S^{\prime }}\left(x,y,z\right)=\left(\frac{x}{1+z},\frac{y}{1+z}\right)$. (14)

We find $\left\{\left(U_{N^{\prime }},{\psi }_{N^{\prime }}\right);\left(U_{S^{\prime }},{\psi }_{S^{\prime }}\right)\right\}$ a differentiable atlas of S².

This method, of stereographic projection, makes it possible to obtain a differentiable structure for any sphere Sⁿ.

For n = 3, the stereographic projections give for S³ the atlas $\left\{\left(U_N,{\phi }_N\right);\left(U_S,{\phi }_S\right)\right\}$ ; where $U_N=S^3\backslash \left\{\left(0,0,0;1\right)\right\}$ and $U_S=S^3\backslash \left\{\left(0,0,0;-1\right)\right\}$. We have

${\phi }_N\left(a,b,c,d\right)=\left(\frac{b}{1-a},\frac{c}{1-a},\frac{d}{1-a}\right)$ (15)

and

${\phi }_S\left(a,b,c,d\right)=\left(\frac{b}{1+a},\frac{c}{1+a},\frac{d}{1+a}\right)$ (16)

B) Link between H and SU(2)

1) Definitions

$\begin{array}{l}SU\left(2\right)=\left\{M\in U\left(2\right)/\det M=1\right\}\\ =\left\{M=\left(\begin{array}{cc}\begin{array}{c}a+ib\\ -c+id\end{array}& \begin{array}{c}c+id\\ a-ib\end{array}\end{array}\right)/a^2+b^2+c^2+d^2=1;\left(a,b,c,d\right)\in {ℝ}^4\right\}\\ =\{M=\left(\begin{array}{cc}{Z}_1& Z_2\\ \overline{-Z_2}& \overline{Z_1}\end{array}\right)/\left|Z_1\right|+\left|Z_2\right|=1;\left(Z_1,Z_2\right)\in {ℂ}^2,\\ \begin{array}{c}\\ \end{array}{Z}_1=a+ib;Z_2=c+id;\left(a,b,c,d\right)\in {ℝ}^4\}\\ =\left\{M\in GL\left(2,ℂ\right)M{\bar{M}}^t=1={\bar{M}}^{t}Met\left|M\right|=1\right\}\end{array}$ (17)

[2] [4] [5].

SU(2) is then identified with group S³, quaternions of module 1. Indeed; consider the homeomorphism f defined as below:

$f:SU\left(2\right)\to S^3$ ;

with

$\left(\begin{array}{cc}{Z}_1& Z_2\\ \overline{-Z_2}& \overline{Z_1}\end{array}\right)\to f\left(\left(\begin{array}{cc}{Z}_1& Z_2\\ \overline{-Z_2}& \overline{Z_1}\end{array}\right)\right)=\left(\begin{array}{cc}1& 0\end{array}\right)\left(\begin{array}{cc}{Z}_1& Z_2\\ \overline{-Z_2}& \overline{Z_1}\end{array}\right)=\left(\begin{array}{cc}{Z}_1& Z_2\end{array}\right)$.

Thanks to f we identify SU(2) with S³. We then consider SU(2) as a sub-body S³ of H, quaternions of module 1

2) Links between SU(2) * U(1) and U(2) as well as between SO(3) and SU(2)

a) Theorem2 "Links between SU(2) * U(1) and U(2)"

$U\left(2\right)\approx \frac{SU\left(2\right)\ast U\left(1\right)}{Z_2}$ (18)

With $I_{2\times 2}=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)$, $-I_{2\times 2}=\left(\begin{array}{cc}-1& 0\\ 0& -1\end{array}\right)$ and $Z_2=\left\{\left(I_{2\times 2},1\right);\left(-I_{2\times 2},-1\right)\right\}$ [2], [4] and [5].

Proof of Theorem 2

We find that $\forall M\in SU\left(2\right)$ and $\forall q\in Z_2$ : $q\ast M=M\ast q$. We have $Z_2$ a distinguished (or normal) subgroup of SU(2) because $\forall q\in Z_2$ ; $q\ast SU\left(2\right)\ast q^{-1}=SU(\; 2\; )$

i) Let

$f:SU\left(2\right)\ast U\left(1\right)\to U\left(2\right)$ ; $\left(M,{\text{e}}^{i\theta }\right)\to f\left(M,{\text{e}}^{i\theta }\right)={\text{e}}^{i\theta }\ast M$

We check that f is a surjective homomorphism, whose kernel is $Kerf=Z_2=\left\{\left(I_{2\times 2},1\right);\left(-I_{2\times 2},-1\right)\right\}$ and Imf = U(2).

Consequently Imf is isomorphic to (SU(2) * U(1))/Kerf. I.e. U(2) ≈ (SU(2) * U(1))/(Kerf)

b) Theorem 3 “Links between SO(3) and SU(2)”

The group SO(3) is isomorphic to the quotient of SU(2) by its center. The center of SU(2) is $Z_2=\left\{I_{2\times 2},-I_{2\times 2}\right\}$, with $I_{2\times 2}=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)$ and $-I_{2\times 2}=\left(\begin{array}{cc}-1& 0\\ 0& -1\end{array}\right)$ [2] [4].

Let us demonstrate SO(3) ≈ (SU(2))/(Z₂) in the following three steps:

1st step: We will establish a homomorphism $g:SU\left(\text{2}\right)\to SO(\; 3\; )$

2nd step: We will prove that g is surjective

3rd step: We will then show from Kerg = Z₂

Proof of Theorem 3

ii)

$\begin{array}{l}{R}_{-q}:\left({ℍ}^{*},.\right)\to \left({ℍ}^{*},.\right)\text{with}{ℍ}^{*}=\left\{q\in ℍ/q\ne 0\right\};\\ \text{and}p\to R_q\left(p\right)=q\ast p\ast \bar{q}\end{array}$ (19)

We find that $R_q$ is a bijection. It is called conjugation by the element q of the multiplicative group ${ℍ}^{*}$. We find that $\left({ℍ}_{pur},.\right)$ is a subgroup of $\left({ℍ}^{*},.\right).$

We make the identification of ${ℍ}_{pur}$ ur and ${ℝ}^3$, thanks to the following bijection ζ:

$\zeta :{ℍ}_{pur}\to {ℝ}^3;\left(0;x;y;z\right)\to \zeta \left(\left(0;x;y;z\right)\right)=\left(x;y;z\right)$ (20)

Thus $\forall q=M\in SU\left(2\right)$, we verify that $R_q$ : ${ℍ}_{pur}\to {ℝ}^3$ given by

$m\to R_q\left(m\right)=q\ast m\ast q^{-1}=q\ast m\ast \bar{q}$ (21)

Let us show that this restriction of $R_q$, to ${ℍ}_{pur}$ or ${ℝ}^3$, is a rotation. Indeed:

1˚) $\forall m_1,m_2\in {ℍ}_{pur}$ ; and $\forall \alpha ,\beta \in R$ : $R_q\left(\alpha \ast m_1+\beta \ast m_2\right)=\alpha \ast R_q\left(m_1\right)+\beta \ast R_q(\; m\; 2\; )$

2˚) $\forall q\in SU\left(2\right)$ ; $\forall m\in {ℍ}_{pur}$, we find $\bar{m}=-m$ i.e. $\bar{m}+m=0$

So $R_q\left(m\right)+R_q\left(\bar{m}\right)=0$ because $\begin{array}{c}{R}_q\left(m\right)+R_q\left(\bar{m}\right)=q*m*\bar{q}+\overline{q\ast m\ast \stackrel{¯}{q}}=q\ast m\ast \bar{q}+\overline{\stackrel{¯}{q}}\ast \bar{m}\ast \bar{q}\\ =q\ast m\ast \bar{q}+q\ast \bar{m}\ast \bar{q}=q\ast \left(\bar{m}+m\right)\ast \bar{q}=q\ast 0\ast \bar{q}=0\end{array}$

The matrix of $R_q$ is therefore a regular matrix

3˚) Let on H be the standard map $\mathfrak{n}$

$n:ℍ\to {ℝ}^{+};p\to \mathfrak{n}\left(p\right)=\sqrt{p\ast \bar{p}}$ (22)

We verify that $\begin{array}{c}n\left(R_q\left(p\right)\right)=\mathfrak{n}\left(q\ast p\ast \bar{q}\right)=\mathfrak{n}\left(q\right)\ast \mathfrak{n}\left(p\right)\ast \mathfrak{n}\left(\bar{q}\right)=\mathfrak{n}\left(p\right)\ast \mathfrak{n}\left(q\right)\ast \mathfrak{n}\left(\bar{q}\right)\\ =\mathfrak{n}\left(p\right)\ast \mathfrak{n}\left(q\ast \bar{q}\right)=\mathfrak{n}\left(p\right)\ast \mathfrak{n}\left(1\right)=\mathfrak{n}\left(p\right)\ast 1=\mathfrak{n}(\; p\; )\end{array}$

Consequently $\forall q\in SU\left(2\right)$, we find that $R_q$ is an orthogonal isomorphism of H; and therefore its restriction to ${ℍ}_{pur}$, (or ${ℝ}^3$ ) is an isometry. In other words: $\forall q\in SU\left(2\right)$, we find $R_q\in O\left(3\right)$.

4˚) As $\forall q\in SU\left(2\right)$, $\det R_q=1$ ; then $R_q$ is a rotation of ${ℝ}^3$ ; i.e. $R_q\in N/SO$ (3).

Thus $g:SU\left(2\right)\to SO\left(3\right)$ ; $q\to g\left(q\right)=R_q$ ; is an application which associates each element q with its rotation matrix $R_q.$

Like $\forall q_1,q_2\in SU\left(2\right)$, we have $R_{q_1}\text{o}{R}_{q_2}=R_{q_1\ast q_2}$ and

$g\left(q_1\ast q_2\right)=R_{q_1}\text{o}{R}_{q_2}=R_{q_1}\ast R_{q_2}$ ; (23)

then g is a group morphism.

We find that $Kerg=Z_2=\left\{I_{2\times 2},-I_{2\times 2}\right\}$.

And as $Img=g\left(SU\left(2\right)\right)=SO\left(3\right)$ ; because $\forall R_q\in SO\left(3\right)$, $\exists q\in SU\left(2\right)$ such that $g\left(q\right)=R_q$ (that is to say g is surjective); then Img = (SU(2))/(Kerg). In other words $SO\left(3\right)\approx \left(SU\left(2\right)\right)/\left(\left\{I_{2\times 2},-I_{2\times 2}\right\}\right)$.

c) Theorem 4 “Links between U(1) and SU(2)”

U(1) is a multiplicative subgroup of SU(2) [4] and [5].

Proof of Theorem 4

We have $U\left(1\right)=\left\{z\in C/z={\text{e}}^{i\theta },\theta \in R\right\}$. U(1) is a group of complex numbers of modules 1. We find the following isomorphisms $f_1$ and $f_2$ :

$f_1:U\left(1\right)\to SO\left(2\right);\text{with}{\text{e}}^{i\theta },\to \left(\begin{array}{cc}\cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right)$ (24)

and;

$f_2:U\left(1\right)\to S^1;\text{with}{\text{e}}^{i\theta },\to \left(\cos \theta ,\sin \theta \right)$ (25)

These isomorphisms identify U(1) with SO(2), or the unit sphere $S^1$ of ${ℝ}^2$.

The immersion of $S^1$ in $S^3$ allows to consider U(1) as a subgroup of SU(2).

1.2. Rotation Dies and Sphere S₃

1.2.1. Geometric Interpretation of a Quaternion

1) Theorem 5

Let q be any non-zero quaternion; noted

$q=\left|q\right|\left(\cos \theta +p\sin \theta \right)$ (26)

with $\theta \in R$ and $p\in {ℍ}_{pur}$. Then the application R_q u is a rotation of angle 2θ and of vector p [1] [3] [5] [6].

Proof of Theorem 5

By demonstrating Theorem 3, we have seen that $\forall q=M\in SU\left(2\right)$, we check that $R_q$ : ${ℍ}_{pur}\to {ℝ}^3$ given by $m\to R_q\left(m\right)=q\ast m\ast q^{-1}=q\ast m\ast \bar{q}$ ; $R_q$ is a rotation, thanks to relations (21) and (22).

Let us use relations (2), (9), (21) and (22) to determine the angle of rotation $R_q$ and his axis.

According to (9), we have $q=\left|q\right|\left(\cos \theta +p\sin \theta \right)=\cos \theta +p\sin \theta$, because $\forall q\in S^3$, $\left|q\right|=1$ and that $p\in {ℍ}_{pur}$, then $p^2=p\ast p=-1$, Thus

$\begin{array}{l}{R}_q\left(m\right)=q\ast m\ast q^{-1}=q\ast m\ast \bar{q}=\left(\cos \theta +p\sin \theta \right)\ast m\ast \left(\cos \theta -p\sin \theta \right)\\ =\left({\cos }^2\theta \right)\ast m+2\left(\sin \theta \cos \theta \right)\ast \left(pлm\right)–\left({\sin }^2\theta \right)\ast \left[m-2\ast p\ast \left(m\cdot p\right)\right]\end{array}$ (27)

Since $\{\begin{array}{l}m\ast p=\left(-m\cdot p,mлp\right)=-m\cdot p+mлp\\ p\ast m-m\ast p=2pлm\end{array}$, then

$\begin{array}{c}{R}_q\left(m\right)=\left({\cos }^2\theta \right)\ast m+\left(\sin 2\theta \right)\ast \left(pлm\right)+\left(2{\sin }^2\theta \right)\ast \left[\left(p\cdot m\right)\ast p\right]-\left({\sin }^2\theta \right)\ast m\\ =\left(\cos 2\theta \right)\ast m+\left(\sin 2\theta \right)\ast \left(pлm\right)+\left(2{\sin }^2\theta \right)\ast \left[\left(p\cdot m\right)\ast p\right]\end{array}$ (28)

Which allows to deduce:

a) If $m=p$, then $R_q\left(m\right)=R_q\left(p\right)=p$ since

$\{\begin{array}{l}pлm=0\\ p\cdot m=m\cdot p=1\end{array}$ (29)

As p is invariant by $R_q$, then p is a vector of the axis of rotation of $R_q$

b) If p is perpendicular to m, then $p\cdot m=0$. It then comes that

$R_q\left(m\right)=\left(\cos 2\theta \right)\ast m+\left(\sin 2\theta \right)\ast \left(pлm\right)$ (30)

By comparing this relation to (9), we conclude that $R_q$ is a rotation of angle 2θ, and of axis p.

Theorem 5 is therefore proven.

We therefore observe that a quaternion q is also a rotation of angle θ, and of axis p.

2) Matrix of a rotation R_q of an element q of SU(2) or S³

Consider the following matrices in SU (2) [5] [6] [9] [10] [11]:

${\sigma }_0=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)$ ; ${\sigma }_1=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right)$ ; ${\sigma }_2=\left(\begin{array}{cc}0& -i\\ i& 0\end{array}\right)$ ; and ${\sigma }_3=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)$, the last 3 of which are Pauli matrices. Note $\sigma =i{\sigma }_1+j{\sigma }_2+k{\sigma }_3=\left({\sigma }_1;{\sigma }_2;{\sigma }_3\right)$. Since $u$ is a unitary matrix of SU(2), and properties of the Pauli matrices, we find ${\left(u\cdot \sigma \right)}^2={\sigma }_0={\left(u\cdot \sigma \right)}^{2n}$ and ${\left(u\cdot \sigma \right)}^{2n+1}=u\cdot \sigma$ ; $\forall n\in N$. Consequently

$\begin{array}{c}q=\exp \left\{-i\theta u\cdot \sigma \right\}\\ ={\sigma }_0-i\left(\theta \right)\left(u\cdot \sigma \right)-{\left(\theta \right)}^2\frac{{\left(u\cdot \sigma \right)}^2}{2!}+i{\left(\theta \right)}^3\frac{{\left(u\cdot \sigma \right)}^3}{3!}\\ +\cdots +i{\left(\theta \right)}^n\frac{{\left(u\cdot \sigma \right)}^n}{n!}+\cdots \\ ={\sigma }_0\cos \theta -i\left(u\cdot \sigma \right)\sin \theta =U\left(u,\theta \right)\end{array}$

We then find that a rotation q of angle θ in $S^3=SU\left(2\right)$ corresponds to a rotation of angle 2θ in $S^2$. The matrix form of this rotation q is given by the relation:

$q=U\left(u,\theta \right)=\left(\begin{array}{cc}\cos \theta -i{u}_z\sin \theta & -\left(i{u}_x+u_y\right)\sin \theta \\ \left(-i{u}_x+u_y\right)\sin \theta & \cos \theta +i{u}_z\sin \theta \end{array}\right)$ (31)

With $u$ unit vector of ${ℝ}^4$ along the axis p of the relation (29). We also have in ${ℝ}^3$ the following unit vectors $u_x=\left(1;0;0\right)$, $u_y=\left(0;1;0\right)$, $u_z=\left(0;0;1\right)$ ; with i imaginary unity. The relation (31) indicates to us that a quaternion is a rotation, in the space of 1/2-spinner. It is a rotation in space of rank 1 spinners.

1.2.2. Euler Angles of a Quaternion

Any matrix q of SU(2) or S³ can be decomposed, using Euler angles, as below [2] [4] [6]:

$\begin{array}{c}q=\left(\begin{array}{cc}{\text{e}}^{i\psi }& 0\\ 0& {\text{e}}^{-i\psi }\end{array}\right)\ast \left(\begin{array}{cc}\cos \theta & i\sin \theta \\ i\sin \theta & \cos \theta \end{array}\right)\ast \left(\begin{array}{cc}{\text{e}}^{i\phi }& 0\\ 0& {\text{e}}^{-i\phi }\end{array}\right)\\ =\left(\begin{array}{cc}\cos \theta {\text{e}}^{i\left(\psi +\phi \right)}& i\sin \theta {\text{e}}^{i\left(\psi -\phi \right)}\\ i\sin \theta {\text{e}}^{-i\left(\psi -\phi \right)}& \cos \theta {\text{e}}^{-i\left(\psi +\phi \right)}\end{array}\right)\end{array}$ (32)

with $0\le \theta \le \pi /2$ ; $0\le \phi \le \pi$ and $-\pi \le \psi \le \pi .$

The approach formulated by relations (31) and (32) is used in Quantum Mechanics, in the theory of angular momentum.

1.3. Fibration

1.3.1. Prelude

The total space P = SU(2) from the topological point of view is identifiable with S³. Consider a group U(1), and perform a class decomposition of SU(2) with respect to U(1). For Let $g\in SU\left(2\right)\backslash U\left(1\right)$. We consider the set $\bar{g}=g\ast U\left(1\right)=\left\{g\ast w/w\in U\left(1\right)\right\}$. Then we consider an element k of SU(2), which is neither in U(1), nor in $\bar{g}$. We build the set $\bar{k}=k\ast U\left(1\right)=\left\{k\ast h/h\in U\left(1\right)\right\}$. We continue the process of class construction in SU(2)\U(1). In order to account, we succeed in writing SU(2) as an infinite union of classes of the type $\bar{g}=g\ast U\left(1\right)$. The quotient set SU(2)/U(1) obtained is identified at $S^2$. The group $P=SU\left(2\right)=S^3$ is considered to be an infinite union of circles $U\left(1\right)=S^1$ ; parameterized by the sphere $U\left(2\right)=S^2=SU\left(2\right)/U\left(1\right)$ [Courses and Seminars (1)] (Figure 1).

The structural group U(1) acts by multiplication on the right, on the total space SU(2). We have chosen U(1) as a distinguished subgroup of SU(2). This fibration in circles of $S^3=SU\left(2\right)$ is called Hopf fibration of $S^3=SU\left(2\right)$ (Figure 2).

Figure 1. Fibration of SU(2) in circles.

Figure 2. Hopf’s fibration [12].

1.3.2. Hopf Fiber of a Lie G Group

1) Algebraic formulation

Let G be a Lie group, and H a Lie subgroup of G. We assume that H is closed in G, so that the topology of the quotient group G/H is separated. We consider the equivalence relation R defining the classes to the left of H. Thus two elements g and k of G are in relation R, if and only if $k\in \bar{g}=g\ast H$. Consider the application h; defined as below [2] [13]:

$h:G\to G/H;g\to h\left(g\right)=\bar{g}=g\ast H.$ (33)

The application h defines a main fibration. Any group G is thus a principal fibered space above G/H. The structural group being H. The quotient set G/H, which is the basis of the bundle, is only a group if H is distinguished in G.

We use the notation H → G → G/H, to characterize a fibration h of G above G/H. In the previous paragraph, we then introduced this fibration U(1) → SU(2) → SU(2)/U(1) (Figure 3).

2) Geometric and topological formulation

a) Definitions of this formulation

$\text{Consider}\left(P;M;\pi \right)\text{where}\{\begin{array}{l}P:\text{is a differential variety}\\ M:\text{is a differential variety}\\ \pi :\text{the projection of}P\text{onto}M\end{array}$ (34)

We say that P is a vector bundled space of base M, of projection π and of fiber-type ${ℝ}^n$ if:

i) $\forall m\in M$ ; ${\pi }^{-1}\left(m\right)$ is a vector space

ii) $\forall m\in M$ ; $\exists$ a neighborhood U of m and a diffeomorphism (35)

$\psi :U\ast {ℝ}^n\to {\pi }^{-1}\left(U\right)$ such that $\{\begin{array}{l}(\text{a})\pi \text{o}\psi \left(m;v\right)=m\\ (\text{b})\psi \left(m;v\right)\text{is linear on the}2\text{nd argument}\end{array}$ [2] [7] [13].

Consider (P; M; π; F; G), where P, M and F are differential varieties. π: the projection of P on M, a differential manifold. G: a diffeomorphism group from F to F. We will say that P is a fibered space of base M, of projection π, of fiber type F and of structural group G if:

i) $\forall m\in M$ ; ${\pi }^{-1}\left(m\right)$ is a manifold differentiable at F

ii) $\forall m\in M$ ; $\exists$ an open neighborhood U of m in M and a diffeomorphism $\psi :U\ast F\to {\pi }^{-1}(\; U\; )$

Figure 3. Fibration of G.

Such that $\{\begin{array}{l}(\text{a})\left(\pi \text{o}\psi \right)\left(m;v\right)=m\\ (\text{b})\text{The applications}{\psi }_{\beta \alpha }={\psi }_{\beta }^{-1}\text{o}{\psi }_{\alpha }^{\left(m\right)}\text{determine the elements of}G\end{array}$

(U; ψ) being called a local trivialization (Figure 4).

b) Definitions

Consider the point $P_o\left(1;0;0\right)$ of $S^2$, a point that we identify with the quaternion $i\left(0;1;0;0\right)$. The Hopf fibration h in terms of rotation is defined as below [1] [3]:

$\begin{array}{l}h:S^3\to S^2;q\to h\left(q\right)=R_q\left(i\right)=q\ast i\ast \bar{q}\\ =\left(a^2+b^2-c^2-d^2;2\left(ad+bc\right);2\left(bd-ac\right)\right)\\ =\left(x;y;z\right)=p;\text{with}q=\left(a;b;c;d\right)\end{array}$ (36)

1.3.3. Fiber over P = (x; y; z) of S²

Let us determine the fiber above a point $P\left(x;y;z\right)=p$ of $S^2$. In other words let us determine the set of points $Q=\left(a;b;c;d\right)=q$ of $S^3$ which verify h(q) = p; i.e. $h^{-1}\left(p\right)=q$ [1] and [3].

1) For $P_o\left(1;0;0\right)=i$ ; we have $h^{-1}\left(i\right)=\left(a;b;c;d\right)$ such that $\{\begin{array}{l}{a}^2+b^2-c^2-d^2=1\\ ad+bc=0\\ bd-ac=0\end{array}$.

Hence $\{\begin{array}{l}{a}^2+b^2=1\\ b=c=0\end{array}$. There are $\theta \in \left[0;2\pi \right]\cap R$ which checks $\{\begin{array}{l}a=\cos \theta \\ b=\sin \theta \end{array}$

Thus

$q=\left(a;b;c;d\right)=\left(\cos \theta ;\sin \theta ;0;0\right)\left(\text{or}q\equiv {\text{e}}^{i\theta }\right)$ (37)

describes a circle of $S^3$.

2) For $P\left(-1;0;0\right)=-i$ ; we have $h^{-1}\left(-i\right)=\left(a;b;c;d\right)$ such that $\{\begin{array}{l}{a}^2+b^2-c^2-d^2=-1\\ ad+bc=0\\ bd-ac=0\end{array}$.

Figure 4. Fiber above the open U.

Hence $\{\begin{array}{l}{c}^2+d^2=1\\ a=b=0\end{array}$. There are $\theta \in \left[0;2\pi \right]\cap R$ which checks $\{\begin{array}{l}c=\sin \theta \\ d=\cos \theta \end{array}$

Thus

$q=\left(a;b;c;d\right)=\left(0;0;\sin \theta ;\cos \theta \right)$ (38)

describes a circle $S^3$.

3) For $P\left(x;y;z\right)=p\in S^2\backslash \left\{\left(1;0;0\right);\left(-1;0;0\right)\right\}$. Let us look for q such that $R_1\left(i\right)=R_q\left(i\right)$. We then find $q=\left(\sqrt{\frac{1+x}{2}};0;\frac{-z}{\sqrt{2\left(1+x\right)}};\frac{y}{\sqrt{2\left(1+x\right)}}\right)$. Indeed; $O{P}_o=\left(1;0;0\right)=i$, and $OP=\left(x;y;z\right)$ such that $O{P}_o\cdot OP=\cos 2\theta =x$, by exploiting the relation (26); and $O{P}_o\text{л}OP=i\text{л}\left(xi+yj+zk\right)=-zj+yk$. Thus $\Vert O{P}_o\text{л}OP\Vert =\Vert O{P}_o\Vert \cdot \Vert OP\Vert \cdot \sin 2\theta =\sin 2\theta$. We then find $\Vert O{P}_o\text{л}OP\Vert =\sqrt{y^2+z^2}=\sqrt{1-x^2}$ because $x^2+y^2+z^2=1$. Thus $\cos \theta =\sqrt{\frac{1+x}{2}}$ and $\sin \theta =\sqrt{\frac{1-x}{2}}$. Consequently a following unit vector $O{P}_o\text{л}OP$ is $u=\frac{O{P}_o\text{л}OP}{\Vert O{P}_o\text{л}OP\Vert }=\frac{1}{\sqrt{1-x}}\left(-zj+yk\right)$. This unitary vector of axis of rotation, makes it possible to bring $P_o$ in P. Consequently, according to (26), we finds $q=\cos \theta +u\sin \theta$. We therefore have

$q=\sqrt{\frac{1+x}{2}}+\frac{-zj+yk}{\sqrt{2\left(1+x\right)}}=\left(\sqrt{\frac{1+x}{2}};0;\frac{-z}{\sqrt{2\left(1+x\right)}};\frac{y}{\sqrt{2\left(1+x\right)}}\right)$ (39)

The relation (37), allows us to affirm that the relation (39) is given by quaternions q describing a circle of $S^3$. Indeed; let $q_1$ and $q_2$ in $S^3$, such that $R_{q_1}\left(i\right)=R_{q_2}\left(i\right)=P\left(x;y;z\right)$.

From relation (23), it comes then that $R_{\overline{q_2}}\text{o}{R}_{q_1}\left(i\right)=i=R_{\overline{q_2}\ast q_1}\left(i\right)$. And according to the relation (37), we obtain $\overline{q_2}\ast q_1={\text{e}}^{i\theta }$ ; with $\theta \in \left[0;2\pi \right]$. We therefore deduce that

$q_1=q_2\ast {\text{e}}^{i\theta };\text{with}\theta \in \left[0;2\pi \right],$ (40)

describes a circle $S^3$.

1.3.4. Local Trivialization of Fibration h

The local trivialization of the h fibration is obtained from relations (39) and (40) below [1] [3]:

The fiber above $p\left(x;y;z\right)\in S^2$, is the circle parameterized by

$\kappa :\left[0;2\pi \right]\to S^3;\beta \to \kappa \left(\beta \right)=q\ast {\text{e}}^{i\beta };\text{with}q=\frac{1}{\sqrt{2\left(1+x\right)}}\left(\left(1+i\right)+yj+zk\right).$ (41)

Indeed; from (9) we have $q=\cos \theta +p_1\sin \theta$ with $p_1\in {ℍ}_{pur}$. And from (26), q induces a rotation of angle 2θ; around the axis $p_1$. We can then study the rotation around $p\left(x;y;z\right)\in S^2$ as two rotations $r_1$ and $r_2$ ; bringing respectively q in i, and i in p. In other words $R_q\left(i\right)=p=r_1\text{o}{r}_2\left(q\right)$. The axis of this rotation is then the vector line perpendicular to the middle of the segment $P_o\left(1;0;0\right)=i$ and $P\left(x;y;z\right)=p$. This axis has the equation:

$w=\left(i+xi+yj+zk\right)/2=\left(\left(1+x\right)i+yj+zk\right)/2$

We find that $\Vert w\Vert =\frac{1}{2}\sqrt{{\left(1+x\right)}^2+y^2+z^2}=\frac{1}{2}\sqrt{2\left(1+x\right)}$ because $x^2+y^2+z^2=1$. We have $\Vert w\Vert =\sqrt{\frac{1+x}{2}}$. Consequently the unit vector, along this axis, is $u=\frac{w}{\Vert w\Vert }=\sqrt{\frac{2}{1+x}}\left[\frac{\left(1+x\right)i+yj+zk}{2}\right]$.

Thus $u=\sqrt{\frac{1}{2\left(1+x\right)}}\left(\left(1+x\right)i+yj+zk\right)$. Consequently $q=\frac{1}{\sqrt{2\left(1+x\right)}}\left(\left(1+x\right)i+yj+zk\right)=u$, which give

$q=\cos \left(\pi /2\right)+\left[\sin \left(\pi /2\right)\right]\left(\left[\left(1+x\right)i+yj+zk\right]\right)/\left|\left(1+x\right)i+yj+zk\right|$ ;

which is a quaternion of the fiber above p. And according to (40), this fiber is of the form $\left\{q\ast {\text{e}}^{i\beta }\right\}$ with $\beta \in \left[0;2\pi \right]\cap R$.

Therefore the application

$\psi :U\ast S^1\to S^3\backslash h^{-1}\left(-i\right),\left(\left(x;y;z\right);{\text{e}}^{i\beta }\right)\to \frac{\left(1+x\right)i+yj+zk}{\sqrt{2\left(1+x\right)}}\ast {\text{e}}^{i\beta }$ ; (42)

allows us to obtain local trivialization (U; ψ).

1.4. Stereographic Projection Visualizing the Fibration of HOPF from the Lie SU(2) Group

Stereographic projection is a method of representing a sphere, deprived of a point, on a plane. [1] [3] [5] [7] [8] [12].

The relation (13) gives us

${\psi }_{N^{\prime }}:U_{N^{\prime }}\to {ℝ}^2,\text{with}{\psi }_{N^{\prime }}\left(x,y,z\right)=\left(x/\left(1-z\right),y/\left(1-z\right)\right)$ (43)

presents the stereographic projection. This sends the point $P\left(x;y;z\right)$, from $S^2\backslash \left\{\left(0;0;1\right)\right\}$, is sent in P’ from the plane Π = Γ. Straight line NP cut Π in P’ (Figure 5).

A circle of $S^2$ which passes through the north pole $N\left(0;0;1\right)$ ; being in a plane ${\Gamma }^{\prime }$, is sent by ${\psi }_{N^{\prime }}$ to the intersection line of Γ and ${\Gamma }^{\prime }$.

Figure 5. Stereographic projection of the north pole on a non-equatorial plane [1].

A circle of $S^2$ which does not pass through the north pole N, is projected into a circle of Γ by ${\psi }_{N^{\prime }}$.

The generalization of stereographic projection to all dimensions being possible, let us take again the relation (15). He comes

$\begin{array}{l}{\phi }_N:U_N=S^3\backslash \left\{\left(0,0,0;1\right)\right\}\to {ℝ}^3,\\ \text{with}{\phi }_N\left(a,b,c,d\right)=\left(b/\left(1-a\right),c/\left(1-a\right),d/\left(1-a\right)\right)\end{array}$ (44)

Let us determine the images of the fibers above the elements of $S^2$, by means of this stereographic projection ${\phi }_N$. We find

1)

$\begin{array}{c}{\phi }_N\text{o}{h}^{-1}\left(P_0\right)={\phi }_N\text{o}{h}^{-1}\left(i\right)={\phi }_N\text{o}{h}^{-1}\left(\left(1;0;0\right)\right)={\phi }_N\left[h^{-1}\left(i\right)\right]\\ ={\phi }_N\left[\left(\cos \theta ;\sin \theta ;0;0\right)\right]\left(\text{fromrelation}\left(37\right)\right)\\ =\left(\frac{\sin \theta }{1-\cos \theta },0,0\right)=\left(\frac{\cos \frac{\theta }{2}}{\sin \frac{\theta }{2}},0,0\right)=\left(\cot \frac{\theta }{2},0,0\right)\end{array}$ (45)

which is the x-axis of the plane (x, y) or (x; z) of ${ℝ}^3=\left(x;y;z\right)$ ; that is to say a circle of infinite radius.

2)

$\begin{array}{c}{\phi }_N\text{o}{h}^{-1}\left(P\right)={\phi }_N\text{o}{h}^{-1}\left(-i\right)={\phi }_N\text{o}{h}^{-1}\left(\left(-1;0;0\right)\right)={\phi }_N\left[h^{-1}\left(-i\right)\right]\\ ={\phi }_N\left[\left(0;0;\sin \theta ;\cos \theta \right)\right]\left(\text{fromrelation}\left(38\right)\right)\\ =\left(0;\sin \theta ;\cos \theta \right)\end{array}$ (46)

which is a circle of the plane (y; z) of ${ℝ}^3=\left(x;y;z\right).$

3)

$\begin{array}{l}{\phi }_N\text{o}{h}^{-1}\left(P\right)={\phi }_N\text{o}{h}^{-1}\left(p\left(x;y;z\right)\right)\\ ={\phi }_N\left(\left(\cos \theta ;0;\frac{-z\ast \sin \theta }{\sqrt{1-x^2}};\frac{y\ast \sin \theta }{\sqrt{1-x^2}}\right)\right)\left(\text{fromrelation}\left(39\right)\right)\\ =\left(0;\frac{-z\ast \sin \theta }{\left(1-\cos \theta \right)\ast \sqrt{1-x^2}};\frac{y\ast \sin \theta }{\left(1-\cos \theta \right)\ast \sqrt{1-x^2}}\right)\\ \text{with}\{\begin{array}{l}\cos \theta =\sqrt{\frac{1+x}{2}}\\ \sin \theta =\sqrt{\frac{1-x}{2}}\end{array}=\left(0;\frac{-z\ast \cot \frac{\theta }{2}}{\sqrt{1-x^2}};\frac{y\ast \cot \frac{\theta }{2}}{\sqrt{1-x^2}}\right)\end{array}$ (47)

and according to the relation (40) is a circle of space ${ℝ}^3=\left(x;y;z\right).$

These different images of the circles, by stereographic projection ${\phi }_N$, give circles which are linked in chains like those of the Olympic games (Figure 6).

1.5. Heat Transfer Equation

1.5.1. Transport Phenomena

1) What do we mean by transport phenomenon?

Figure 6. The stereographic projection transforms the circles of the sphere into circles of the projection plane [7] [8].

In thermo-statistics, we group together a set of physical phenomena of movement of molecules in fluids or solids in transport phenomena. Transport phenomena are the processes in which there is a global transfer (or transport) of matter, energy or momentum on a macroscopic scale. The essential physical aspects of these phenomena can be described by similar methods and are characterized by an equation which links the variation in time and in space of a quantity which describes the phenomenon. In the simplest case this equation is of the form [14] [15] [16] [17]:

$\frac{\partial \xi }{\partial t}=a^2\frac{{\partial }^2\xi }{\partial x^2}$ (48)

where a is a constant characteristic of each physical situation and ξ is a quantity corresponding to the particular phenomenon of transport studied like the diffusion of the molecules, the viscosity and the conduction of heat. In this work, we will focus on the last case of these types of transport phenomena whose mathematical solutions relate to a temperature distribution.

The phenomenon of energy and/or material transport is of great importance in statistical physics. This phenomenon for heat is described by the Fourier equation which we will solve for a spatial dimension before modeling this solution for the Kankule site. The solution of this equation is called the geotherm

2) One-dimensional Fourier equation

The Fourier equation for the conduction of heat in the ground, when this energy varies as a function of time t and depth z is:

$\frac{{\partial }^{2}T}{\partial z^2}-\frac{1}{\alpha }\frac{\partial T}{\partial t}=0$ (49)

[5] [14] [15] [16] [18] [19] [20] [21] [22].

With $\alpha$ the coefficient of thermal diffusivity of the soil, z the depth and T the temperature characterizing the heat transported.

With the boundary conditions of the first species (Dirichlet problem), where the surface temperature is known at all times, this equation admits a sinusoidal solution of progressive wave of the form:

$T\left(z,t\right)={\text{e}}^{i\left(kz-\omega t\right)}$ (50)

whose

$\omega$ : is the pulsation of the waves

$\frac{\omega }{k}=v$ : is the phase velocity of these waves

i: the imaginary unit (with i² = −1)

Introduce the values of (50) in (49), we find:

$-k^2+i\frac{\omega }{\alpha }=0$ (51)

So,

$\omega =i\alpha k^2$ (52)

As a result

$k=\pm \left(1+i\right)\sqrt{\frac{\omega }{2\alpha }}$ (53)

is

$k=\frac{\pm \left(1+i\right)}{d}$ (54)

by asking

$d=\sqrt{\frac{2\alpha }{\omega }}$ (55)

d: is called the depth of penetration of the waves (or also the depth of damping).

Starting from (49), according to Wu and Nofziger [22], taking into account the boundary conditions, the previous solution can be written:

$T\left(z,t\right)=\overline{T_0}+\overline{A_0}\cdot {\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right)$ (56)

Using solution of (49), we looked for the gradient, for each site, with the relation $\partial T/\partial z=f\left(z\right)$ and deduce the depth of exploration to the Kankule site studies

3) Exploitation of geothermy

The Geoscience review of March 16, 2013, shows the operating conditions of geothermal energy to produce electricity. Having active tectonic or volcanic zones in which a surface heat source makes it possible to have a geothermal fluid at a temperature sufficient to produce electricity [23]. On page 14, Figure 8 [23], already shows that with a temperature of 50˚C at 1 km deep, we can consider the production of electricity. It will then be necessary to dig for this to more than 3 km deep, in order to trap the geothermal fluid at more than 150˚C. The same article presents a world map where high energy geothermal energy is likely to be exploited. Eastern DR. Congo is there, as in all the countries of the Albertine Rift (or Valley). From the ground temperature down to 90˚C, electricity can already be produced [24].

The latest study distinguishes between surface and deep geothermal energy as follows:

­ <30˚C and <500 m: surface geothermal energy. This takes up the very low temperatures (enthalpy).

­ between 25˚C - 30˚C and 150˚C and >500 m: deep geothermal energy. This takes over the low and medium temperatures (enthalpy).

Still the same study, specifies that if the geothermal reservoir reaches a temperature above 150˚C and is encountered at shallow depths, the resource must therefore be in a region where the geothermal gradient is normal (~30˚C/km) or beyond normal. This is for the production of electricity. This classification is consistent with that suggested by Lindal for France [25].

The exploitation of high energy geothermal energy is favorable in areas which knows geodynamic contexts, showing some form of volcanism. At these points, a significant part of the internal energy of the globe is released. The history of mountain chain formations and terrestrial magmatism allows us to understand the genesis of magma chambers. The latter are foci, zones with strong thermal gradient [26].

1.5.2. Hopf Fiber of the Heat Transfer Phenomenon

Relations (50) or (56) determine solutions of Equation (49).

Take

$q=\frac{T\left(z,t\right)}{\left|T\left(z,t\right)\right|}=\left(\cos \theta ;\sin \theta ;0;0\right);\text{with}\theta =\left(\pm z/d-\omega t\right)\in \left[0;2\pi \right]\cap R$ (57)

By comparing to relation (38), we observe that q is a fiber above the north pole. Our q being a solution of (50) or (56); we conclude that each solution of these equations is associated with a fiber in $SU\left(\text{2}\right)=S^3$.

For the three-dimensional Fourier equation

$\frac{\partial T}{\partial t}=\alpha \Delta T;\text{with}\Delta =\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial x{z}^2}$ (58)

the Laplacian; there is always a quaternion of the type (56) which is a solution; by means of a modification of x by a three-dimensional vector. In short, there is always a quaternion unit, which is a fiber. From the above, we conclude that for each solution of the Fourier heat transfer equation, there is an associated fiber in $SU\left(\text{2}\right)=S^3$. This result is analogous to what in classical electromagnetism is called magnetic field lines. What we can visualize as below (Figure 7).

Figure 7. Model of a heat fiber [7] [12].

This model is analogous to that obtained by T.M. and Dirk Bouwmeester [27], for the modeling of lights (electromagnetic waves of the visible spectrum) with Hopf fibration. Indeed, each fiber appears as a line of Earth’s magnetic field. And since two distinct lines of fields don’t join, so do two fibers in circles that don’t cut. And since two distinct lines of fields don’t join, so do two fibers in circles that don’t cut. As a result each geomagnetic field line is representative of a heat fiber. By exploiting Maxwell’s equations of the Terrestrial geomagnetic field, we find that

$B=B_0\ast {\text{e}}^{-k_{1}x}{\text{e}}^{i\left(k_{2}x-\bar{\omega }t\right)}$ [28], (59)

which is a solution of the same type as (50).

2. Methodology

To do this work, we used the following methodological approach:

2.1. Analytic

To do this we sought:

­ Let us relate it between Quaternions units and the Lie group SU(2);

­ determining the coefficients of the geotherm;

­ The links between the Geotherm of the one-dimensional Fourier equation with the Hopf fibration of the sphere $S^3$ ;

­ See how to visualize Hopf fibration using stereographic projection;

­ The significance in terms of bundle of the geotherm.

2.2. Experimental

1) From 8 h to 14 h, with the ground thermometer, we took samples from 14 sites, called Kankule I, Kankule II,…, and Kankule XIV:

We took soil temperatures up to 15 centimeters deep, every day, in 52 weeks. This in order to characterize each site by the average temperature of the surface soil;

Also took the temperatures of the thermal waters, over the same period from September 2010 to December 2014, in order to see how to also differentiate the sites from the electrolytic properties of its waters. The analysis of which was carried out at the Lwiro Natural Sciences Research Center (CNRS/Lwiro) in Katana.

2) With a Garmin brand GPS, we determined the geographic coordinates of each site. Using these coordinates we had:

­ Geo referenced each site

­ Downloaded on the NOAA site, the geomagnetic fields of each site, for the same years of our temperature samples [29].

2.3. Comparative

To do this, we compared our results from the ground geotherm with various accepted results, in order to see if each thermal source can be compared to a bundle.

2.4. Statistical

Using statistics we were able to:

­ Determine different coefficients involved in the geotherm of each site.

­ Perform the linear regression of the geomagnetic field for each site according to the temperature field of the same site.

3. Results and Analysis

1) Our study rather leads to the result:

$T\left(z,t\right)=\overline{T_0}+\overline{A_0}\cdot {\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)$ (60)

as solution of Equation (49).

Indeed, consider the depreciated waves real parts of the solution (50), which are $U_1\left(z,t\right)={\text{e}}^{-\frac{z}{d}}\cos \left(\omega t-\frac{z}{d}\right)$ and $U_2\left(z,t\right)={\text{e}}^{-\frac{z}{d}}\sin \left(-\omega t+\frac{z}{d}\right)$. Superposing the two waves, we find the solution:

$\begin{array}{c}{U}_3\left(z,t\right)=U_1\left(z,t\right)+U_2\left(z,t\right)\\ ={\text{e}}^{-\frac{z}{d}}\cos \left(\omega t-\frac{z}{d}\right)+{\text{e}}^{-\frac{z}{d}}\sin \left(-\omega t+\frac{z}{d}\right)\\ ={\text{e}}^{-\frac{z}{d}}\left(\cos \left(\omega t-\frac{z}{d}\right)+\sin \left(-\omega t+\frac{z}{d}\right)\right)\\ ={\text{e}}^{-\frac{z}{d}}\left[\cos \left(\omega t-\frac{z}{d}\right)-\sin \left(\omega t-\frac{z}{d}\right)\right]\\ ={\text{e}}^{-\frac{z}{d}}\left[\cos \left(\omega t-\frac{z}{d}\right)-\cos \left(\frac{\pi }{2}-\left(\omega t-\frac{z}{d}\right)\right)\right]\end{array}$

$\begin{array}{l}={\text{e}}^{-\frac{z}{d}}\left[-2\sin \left(\frac{\omega t-\frac{z}{d}+\frac{\pi }{2}-\omega t+\frac{z}{d}}{2}\right)\ast \sin \left(\frac{\omega t-\frac{z}{d}+\frac{\pi }{2}-\omega t+\frac{z}{d}}{2}\right)\right]\\ ={\text{e}}^{-\frac{z}{d}}\left[-2\sin \left(\frac{\pi }{4}\right)\ast \sin \left(\omega t-\frac{z}{d}-\frac{\pi }{4}\right)\right]\\ ={\text{e}}^{-\frac{z}{d}}\left[-2\frac{\sqrt{2}}{2}\ast \sin \left(\omega t-\frac{z}{d}-\frac{\pi }{4}\right)\right]\\ =\sqrt{2}{\text{e}}^{-\frac{z}{d}}\sin \left(\pi +\omega t-\frac{z}{d}-\frac{\pi }{4}\right)\\ =\sqrt{2}{\text{e}}^{-\frac{z}{d}}\sin \left(\omega t-\frac{z}{d}+\frac{3\pi }{4}\right)\end{array}$

Multiplying $U_3\left(z,t\right)$ by a real constant A, we still get a solution

$U_4\left(z,t\right)=A\ast U_3\left(z,t\right)$ of Equation (49). By posing $\overline{A_0}=A\ast \sqrt{2}$, then we write

$U_4\left(z,t\right)=\overline{A_0}\ast {\text{e}}^{-\frac{z}{d}}\sin \left(\omega t-\frac{z}{d}+\frac{3\pi }{4}\right)$. We have a particular solution of (49)

which is $U_5\left(0,t\right)=\overline{T_0}$, with $\overline{T_0}$ a real constant. Finally by supposing the waves $U_4\left(z,t\right)$ and $U_5\left(0,t\right)$, we get the solution (60). In other words, we have:

$T\left(z,t\right)=U_5\left(0,t\right)+U_4\left(z,t\right)=\overline{T_0}+\overline{A_0}\cdot {\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)$

Using statistical methods, we obtained, for each Kankule site, the coefficients of the geotherm (56) or ((56)-bis) given by Table 1.

2) By seeking the linear regression of the geomagnetic field B, expressed in micron-tesla, as a function of the temperature field T of the soil, expressed in degrees Celsius, we have obtained the relation of the type B = aT + b. The results of this regression are summarized in Table 2.

In this table, Ho is the null hypothesis. It translates the acceptance of the linear dependence between B and T. On each column, where Ho appears, the value 1 means acceptance of the null hypothesis. On the other hand the value 0, implies the rejection of the null hypothesis.

In order to better understand whether each site should be considered separately, we have analyzed which ions are found in different thermal waters, per site. Table 3 gives us a summary of the results found.

3) The geotherm (60) gives us the gradient

$\text{grad}\left(z,t\right)=-\sqrt{2}\frac{A_0}{d}{\text{e}}^{-\frac{z}{d}}\cos \left(\frac{2\pi \left(t-t_0\right)}{365}--\frac{z}{d}+\frac{\pi }{2}\right)$ (61)

Is $\begin{array}{l}g\left(z\right)=\frac{\partial T}{\partial z}=-\frac{A_0}{z}{\text{e}}^{-\frac{z}{d}}\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\\ -\frac{A_0}{z}{\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\end{array}$ a function that varies with depth z.

Table 1. The coefficients of the geotherm by Kankule thermal spring.

An array where the coefficients of the two parallel lines differ.

Table 2. Reduction of the magnetic field B as a function of the soil temperature.

Table 3. Thermal water electrolysis data and associated thermal field.

Let’s determine the maximum values of this function, using differential calculus. It comes at a given depth

$\begin{array}{l}\{\begin{array}{l}g\left(z\right)>0\\ g^{\prime }\left(z\right)=0\\ g^{″}\left(z\right)<0\end{array}\\ \iff \{\begin{array}{l}-\frac{A_0}{z}{\text{e}}^{-\frac{z}{d}}\left(\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]+{\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\right)>0\\ 2\frac{A_0}{d^2}{\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]=0\\ 2\frac{A_0}{d^3}{\text{e}}^{-\frac{z}{d}}\left(\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]-\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\right)<0\end{array}\\ \iff \{\begin{array}{l}g\left(z\right)=\frac{A_0}{z}{\text{e}}^{-\frac{z}{d}}>0\\ \frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}=\left(2k+1\right)\frac{\pi }{2}\\ \sin \left(2k+1\right)\frac{\pi }{2}<0\end{array};\text{with}k\in \mathbb{Z}\end{array}$ (62)

With $\frac{z}{d}>0$, the first value of z which makes the gradient maximum, is obtained for k = −1.

We then find the following results for the Kankule IV site:

(Abscissa of the maximum) $z=d\ast \pi \left(1.25+60/365\right)=4.537828208\text{m}$ and the maximum gradient

$g\left(z\right)=\partial T/\partial z=0.125471195˚\text{C}/\text{m}$ (63)

For all site of Kankule the average maxima gradient give

$\bar{m}=0.117\left[˚\text{C}/\text{m}\right]$ (64)

4. Discussion of the Results

4.1. Which of the Geotherms (56) and (60) Verify (49)?

To find out, we then introduce these results into Equation (49) and see more.

For Equation (49) for which Wu and Nofziger find the solution.

$T=T\left(z,t\right)=T_a+A_0{\text{e}}^{-\frac{z}{d}}\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right]$,

with $d=\sqrt{\frac{2\alpha }{\omega }}$ and $\omega =2\pi /365.$

He then comes

1)

$\begin{array}{c}{T^{\prime }}_t\left(z;t\right)=\frac{1}{\alpha }\frac{\partial T}{\partial t}=\frac{1}{\alpha }\frac{2\pi }{365}{A}_0{\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right]\\ =\frac{A_0}{\alpha }\omega {\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right]\end{array}$ (65)

2)

$\begin{array}{l}{T^{″}}_z\left(z;t\right)=\frac{{\partial }^{2}T}{\partial z^2}=\frac{\partial }{\partial z}\left(\frac{\partial T}{\partial z}\right)\\ =\frac{\partial }{\partial z}\left[-\frac{A_0}{d}{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right)-\frac{A_0}{d}{\text{e}}^{-\frac{z}{d}}\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right)\right]\\ =-\frac{A_0}{d}\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\left(\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{2}\right)-\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)\right)\right]\\ =\frac{2A_0}{d}\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)\sin \left(-\frac{\pi }{4}\right)\right]\\ =-\frac{\sqrt{2}}{d}{A}_0\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)\right]\end{array}$

$\begin{array}{l}=-\frac{\sqrt{2}}{d}{A}_0\left[-\frac{1}{d}{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)-\frac{1}{d}{\text{e}}^{-\frac{z}{d}}\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)\right]\\ =\frac{\sqrt{2}}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\left[\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)+\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)\right]\\ =\frac{\sqrt{2}}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\left[\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)+\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}+\frac{2\pi }{4}\right)\right]\\ =\frac{\sqrt{2}}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\left[\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}-\frac{\pi }{4}\right)+\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{4}\right)\right]\\ =\frac{2}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)=\frac{A_0}{\alpha }\omega {\text{e}}^{-\frac{z}{d}}\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right]\end{array}$ (66)

Comparing solutions (65) and (66), we find that the two results differ. We then conclude that solution (56), given by Wu and Nofziger is not suitable as solution of (49).

The result (60) that we found for the geotherm being

$T=T\left(z,t\right)=T_a+A_0{\text{e}}^{-\frac{z}{d}}\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]$,

with $d=\sqrt{\frac{2\alpha }{\omega }}$ and $\omega =\frac{2\pi }{365}$.

Then he comes:

1)

$\begin{array}{c}{T^{\prime }}_t\left(z;t\right)=\frac{1}{\alpha }\frac{\partial T}{\partial t}=\frac{1}{\alpha }\frac{2\pi }{365}{A}_0{\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\\ =\frac{A_0}{\alpha }\omega {\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right]\end{array}$ (67)

2)

$\begin{array}{l}{T^{″}}_z\left(z;t\right)=\frac{{\partial }^{2}T}{\partial z^2}=\frac{\partial }{\partial z}\left(\frac{\partial T}{\partial z}\right)\\ =\frac{\partial }{\partial z}\left[-\frac{A_0}{d}{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)-\frac{A_0}{d}{\text{e}}^{-\frac{z}{d}}\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)\right]\\ =-\frac{A_0}{d}\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\left(\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)+\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}+\frac{2\pi }{4}\right)\right)\right]\\ =-\frac{2A_0}{d}\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}+\frac{\pi }{4}\right)\cos \left(\frac{\pi }{4}\right)\right]\\ =-\frac{\sqrt{2}}{d}{A}_0\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\pi \right)\right]\end{array}$

$\begin{array}{l}=\frac{\sqrt{2}}{d}{A}_0\frac{\partial }{\partial z}\left[{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)\right]\\ =\frac{\sqrt{2}}{d}{A}_0\left[-\frac{1}{d}{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)-\frac{1}{d}{\text{e}}^{-\frac{z}{d}}\cos \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)\right]\\ =-\frac{\sqrt{2}}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\left[\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}\right)+\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{2}\right)\right]\\ =-\frac{\sqrt{2}}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\left[2\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{4}\right)\cos \left(-\frac{\pi }{4}\right)\right]\end{array}$

$\begin{array}{l}=-\frac{2}{d^2}{A}_0{\text{e}}^{-\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{4}\right)\\ =-\frac{A_0}{\alpha }\omega {\text{e}}^{-\frac{z}{d}}\sin \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{4}\right]\\ =\frac{A_0}{\alpha }\omega {\text{e}}^{-\frac{z}{d}}\cos \left[\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{\pi }{4}+\frac{2\pi }{4}\right]\end{array}$ (68)

Comparing solutions (67) and (68), we find that the two results are equal. We then conclude that the relation (60), is suitable as solution of (49).

This analysis of the solution used by Wu and Nofziger shows that it is not the integral of the one-dimensional Fourier Equation (49). This is why it has not been exploited by Alonso/Finn [15]. The result of the latter verifies the Fourier equation. Its only problem is that some constants of the geotheme are not well specified

4.2. Does Each Thermal Water Site Represent a Bundle?

As a result, the geotherm (50) takes the form: $T\left(z,t\right)={\text{e}}^{i\left(\pm \left(\frac{1+i}{d}\right)z-\omega t\right)}$ ; a quaternion of module $\left|T\left(z,t\right)\right|=\{\begin{array}{l}{\text{e}}^{\frac{-z}{d}}\\ {\text{e}}^{\frac{z}{d}}\end{array}$.

The solutions which have a physical meaning are those for which $\left|T\left(z,t\right)\right|={\text{e}}^{\frac{-z}{d}}$, corresponding to the propagation of the heat-absorbed waves.

Always starting from (50), we get the following real solution for (60):

$T\left(z,t\right)=\overline{T_0}+\overline{A_0}{\text{e}}^{\_\frac{z}{d}}\sin \left(\frac{2\pi \left(t-t_0\right)}{365}-\frac{z}{d}+\frac{3\pi }{4}\right)$

Result called integral curve of (49) or geothermal. It is this result that Adelin MULENDA [1] applies for the Kankule site.

If we introduce the values of Table 1, in the geotherms (60) or (67), then in (57), we find that we are dealing with a bundle, and this for each site. The coefficients obtained from the geotherm of each, are different from another site in Kankule. We can then say that each site is a different bundle than that of a site from the same environment. Paragraph 2.1.4: The visualization of bundles by means of stereographic projection, allows us to say that two different bundles in circles, do not meet. Such a situation should be verified with the physical or chemical properties of the thermal waters of the sites. For that, we will have to check if at least one same source feeds two different sites. This could help us understand why each thermal bath in the world’s largest power plant, The Geysers in the United States, exploits each source to generate electric current [25]. It is the same for the geothermal power station of Soultz for Alsace in France.

4.3. Comparison of Ion Density Values in Site Thermal Baths

Can we say that the same source of thermal water supplies two different sites? Here is the result which seems to appear from the study of Table 3. We have the same ions in all thermal waters: Cl⁻ (mg/L), SO− 4 (mg/L), HCO− 3 (mg/L), PO− 4 (μmol/L), NH+ 4 (μmole/L) and NO− 3 (μmol/L), which suggests that we have the same source everywhere. However, the analysis of the geological map of the DRC [30] or that of Kivu [31], what can we say? We are in the presence of the same recent basic lava. What our first impressions can only translate is that our Kankule thermal springs have only been in contact with the same source rock. To realize this, we divide the numerical data of the two parallel lines each time. If the results found give a value of 1 everywhere, then two different sites should be supplied by the same groundwater source. What we find as a result of the calculations, for example between Kankule I and II, are really different. We find that the ion concentrations or densities are very different between two different sites. We come to the conclusion that in depth, two different sources do not meet after being in contact with the heat coming from the magma chamber (or hot rock). In other words, deep heat islands of the same basalt, but differently heated at two different points, are sources of energy observable in Kankule through the thermal waters. This corroborates the idea of bundles discussed in the previous paragraph. In fibration language, what are we going to call the thermal water sites and the magma chambers that supply them? The visualization of fibration by stereographic projection, from paragraph 2.4.3., allows us to say that the deep magma chamber is the point above which the bundle is made. This hot rock is a source of heat, so we call it the heat core. The thermal water site, which is a bundle above the magma chamber, we call it bundle of heat.

Are there parallels between the heat bundles and the magnetic field lines?

The Hopf fibration of electromagnetic waves obtained by [27] and the one we obtained using relations (33), (36), (37) and (57), lead us to think that these two structures on S³, have links. In an electromagnetic wave there is always an interaction between a magnetic field and an electric field. And this at all times. The results in Table 3 help us to see that the heat from each Kankule site has contributed to the ionization of the thermal waters. There can then exist a certain bond between the geomagnetic field in a site with the field created by heat, represented by the temperature of the ground of the site of thermal baths. This observation leads us to discuss the results of Table 2. Is there a linear regression of the geomagnetic field B in temperature of the soil of the corresponding site? We are tempted to answer with yes when the null hypothesis Ho is accepted, and with no when this null hypothesis is rejected. In all cases, we arrive at a non-zero correlation between the temperature and the Earth’s magnetic field of the corresponding site. We have found by way of illustration for Kankule I, using the R software, and exploiting the statistics, the following results:

>summary (dt1)

Call:

lm (formula = zmatrix.temp.champ1[, 1] ~ zmatrix.temp.champ1[, 2])

Residuals:

Min 1Q Median 3Q Max

−1.9712 −0.6506 0.0449 0.5635 3.7750

Coefficients:

Estimate Std. Error t value Pr (>|t|)

(Intercept) 4.964e−15 1.399e−01 0.000 1.000

zmatrix.temp.champ1[, 2] 5.196e−02 1.412e−01 0.368 0.714

Residual standard error: 1.009 on 50 degrees of freedom

Multiple R-squared: 0.0027, Adjusted R-squared: −0.01725

F-statistic: 0.1354 on 1 and 50 DF, p-value: 0.7145

Result that we can write

$B=5.196\text{th}-02T-4.964\text{th}-15$ ; 5%. (69)

In other words, the linear dependency hypothesis between B and T is accepted up to the 5% threshold. The B magnetic field is expressed in micron-Tesla and the temperature T in degrees Celsius. Matching the site’s magnetic field with soil temperature in situ helps us understand how the geomagnetic field can model heat fiber. In this passage we bypassed the topological approach. Introducing (69) in (49), we find:

$\frac{{\partial }^{2}B}{\partial z^2}-\frac{1}{\alpha }\frac{\partial B}{\partial t}=0$ (70)

This equation is that of Fourier, applied to the terrestrial geomagnetic field. She explains that the Earth’s magnetic field is increasing from the surface of the ground to the center of the Earth, fact which is justified by the nonzero thermal gradient, or then the nonzero geomagnetic gradient, phenomenon that allows various applications of geothermal energy [16]. From Maxwell’s equations of the classical electromagnetic field, Faraday-Henry’s law gives us:

$\text{rot}E=-\frac{\partial B}{\partial t}$ (71)

Equation expresses that the variation of the magnetic field in time generates the electric field. The integral solution of (70), being of type (50) or (57), then the earth’s magnetic field is really a bundle. It is a Hopf fibration generated by Charge movements, according to Maxwell’s equations. Our approach allows us to validate the Hopf fibration model that we carried out in paragraphs 2.1.3 and 2.1.4. Thus the electromagnetic waves are then modeled as a Hopf Fibration. This is also the result of the work of William T.M. and Dirk Bouwmeester [27]. And reciprocally, the linear links between B and T, allows us to corroborate the notion of fibration of the field T, for the site of Kankule I.

For the other Kankule sites, we had obtained similar results. Even if the null hypothesis is rejected, as in Kankule IV. There is always a certain linearity, as weak as it is, between the magnetic field of the site and the geomagnetic field. This implies that there is an interaction between the magnetic field of the site and the soil of the corresponding site, however small. Indeed, the ionization of thermal waters observed with the results of Table 3, allows the creation of the electromagnetic field in place [16]. And this field in thermal waters interacts with the geomagnetic field of the place. In view of the above, we can say that the geomagnetic field allows us to detect the state of electrification at any Kankule geothermal site. There is therefore parallelism, in terms of interaction, between the B and T fields of a given Kankule site.

4.4. Which Geothermy Could Be Exploited in Kankule?

The average Kankule thermal gradient of 0.117˚C/m differs from the global average thermal gradient of 0.033˚C/m. This Kankule gradient is 3.5 times greater than the global average thermal gradient. With this average Kankule thermal gradient, we should be able to exploit different types of geothermal energy as indicated [24]. The different types of geothermal energy that can be used in Kankule are:

• <30˚C and <500 m: surface geothermal energy. This takes up the very low temperatures (enthalpy). Because at 257 m deep, we will have temperatures of at most 30˚C.

• between 25˚C - 30˚C and 150˚C and >500 m: deep geothermal energy. Because at a depth of 500 m, we will already have a temperature of 58.5˚C. In order to generate electricity, it would have to be over 100˚C. By digging up to 1500 m, or 1.5 km, deep, we will have temperatures reaching 175.5˚C. Going up to 2000 m deep, we will have a temperature of 234˚C. The production of electricity is then possible from 50˚C as already pointed out [23].

• For Kankule IV, with a gradient of 125, 471, 195˚C/km. You can reach desired temperatures by drawing water more than 2 km deep, reaching temperatures of 234˚C. A good result when we already compare them to that of Soultz in France [25] [26].

5. Conclusions

In the present work, we have found that:

­ Each thermal water site is represented by the temperatures of the soil or of its waters;

­ Each site is a heat beam taken as a fibrous equivalent to a line of the magnetic field of the place;

­ There is interaction between the bundle of the site’s geomagnetic field and the bundle of heat.

Symbols Used

q: quaternion

(a, b, c, d): unit quaternion

Sⁿ: sphere unit of dimension n

${\psi }_{N^{\prime }}$ : stereographic projection in R² from the North Pole N'

${\phi }_N$ : stereographic projection in R³ from North Pole N

i: imaginary unit

T: temperature

$\overline{T_0}$ : average annual soil temperature

$\overline{A_0}$ : the annual average thermal amplitude

d: damping depth

α: thermal diffusivity coefficient

z: depth in the ground

ω: annual pulse of thermal waves

B: module of the Earth’s geomagnetic field

t: time

σ: temperature variance

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Lyons, D.W. (2006) An Elementary Introduction to the Hopf Fibration. Lebannonn Valley College, Annville, 1-16.
[2]	Kalongama, J.M. (2009-2010) Méthodes de Géométrie différentielle appliquées en Sciences Physique. Séminaire de DEA animé à l’UPN/ISP en RDC; Inédit.
[3]	Lafuente-Gavy, L., et al. (2009) La fibration de Hopf. Université Libre de Bruxelles, Bruxelles, 26-36.
[4]	Kalongama, J.M. (2007-2008) Sur la généralisation isotopique de la théorie quantique du moment angulaire. Thèse de doctorat défendue à l’UNIKIN.
[5]	Faraut, J. (2006) Analyse sur les groupes de Lie. Calrenge &Mounet, Paris.
[6]	Hladik, J., et al. (2009) Mécanique Quantique. Atomes et Noyaux: Applications. Dunod, Paris.
[7]	Jean-François Arbour, UQAM (2012) S3 et la fibration de Hopf.
[8]	https://fr.wikipedia.org/wiki/Projection_stérréographique
[9]	Messiah, A. (1960) Mécanique Quantique, Tome 2. Dunod, Paris.
[10]	Tannoudji, C.C., et al. (1977) Mécanique Quantique, Tome 2. Herman, Paris.
[11]	Tannoudji, C.C., et al. (2007) Mécanique Quantique, Tome 2. Herman, Paris.
[12]	http://www.dimensions-math.org/Dim_CH7.html
[13]	Navez, J. (2009-2010) Théorie des groupes. Séminaire de DEA animé à l’UPN/ISP en RDC, Inédit.
[14]	Mulenda, A. (2013) Modélisation du transfert des températures sol-air-eau thermale de Kankule dans Katana au Sud-Kivu. Mémoire de D.E.A présenté à l’Université Pédagogique Nationale de Kinshasa (UPN /Kin), Faculté des Sciences, Département de Physique.
[15]	Finn, A., et al. (1986) Physique générale 1, Mécanique et thermodynamique. Inter Editions, Paris.
[16]	Finn, A., et al. (1977) Physique générale 2, Champs et Ondes. Inter Editions, Paris.
[17]	Reif, F. (1989) Physique statistique, Berkeley: Cours de Physique, Volume 5. Armand Colin, Paris.
[18]	Lavrentiev, M., et al. (1972) Méthodes de la théorie des fonctions d’une variable complexe. Mir, Moscou.
[19]	Piskounov, N. (1972) Calcul différentiel et intégral, Tome 2. Mir, Moscou.
[20]	Nofzinger, D.I., et al. (2003) Soil Temperature Change with Time and Depth: Theory.
[21]	Chatterji, S.D. (1998) Cours d’Analyse. 3 Equations différentielles ordinaires et aux dérivées partielles. Presses Polytechnique et Universitaires romandes, Bienne.
[22]	Koudriatsev, V., et al. (1986) Cours élémentaire de mathématiques supérieures. Mir, Moscou.
[23]	Géosciences (2013) La revue du BRGM pour une Terre durable, No. 16.
[24]	Ecorem, S.A. (2011) Etude des obstacles à la géothermie profonde (basse et haute énergie). Rapport final, Département de l’énergie, Wallonie.
[25]	https://fr.wikipedia.org/wiki/Géothermie
[26]	Boden, J.-P., et al. (1988) Biologie et Géologie, Première S. Bordas, Paris.
[27]	Irvine, W.T.M. and Bouwmeester, D. (2008) Linked and Knotted Beams of Light. Nature Physics, 4, 716-721. https://doi.org/10.1038/nphys1056
[28]	http://www.geologie.ens.fr
[29]	http://www.arl.noaa.gov/ready.geomagn
[30]	Musée Royal de l’Afrrique Central (2005) Carte géologique de la République Démocratique du Congo, Tervuren.
[31]	De la carte géologique du Sud-Kivu, 1939.