Liouville–Bratu–Gelfand equation explained

For Liouville's equation in differential geometry, see Liouville's equation.In mathematics, Liouville–Bratu–Gelfand equation or Liouville's equation is a non-linear Poisson equation, named after the mathematicians Joseph Liouville,[1] Gheorghe Bratu[2] and Israel Gelfand.[3] The equation reads

\nabla2\psi+λe\psi=0

The equation appears in thermal runaway as Frank-Kamenetskii theory, astrophysics for example, Emden–Chandrasekhar equation. This equation also describes space charge of electricity around a glowing wire[4] and describes planetary nebula.

Liouville's solution[5]

In two dimension with Cartesian Coordinates

(x,y)

, Joseph Liouville proposed a solution in 1853 as

λe\psi(u2+v2+1)2=2\left[\left(

\partialu
\partialx

\right)2+\left(

\partialu
\partialy

\right)2\right]

where

f(z)=u+iv

is an arbitrary analytic function with

z=x+iy

. In 1915, G.W. Walker[6] found a solution by assuming a form for

f(z)

. If

r2=x2+y2

, then Walker's solution is

8e-\psi=λ\left[\left(

r
a

\right)n+\left(

a
r

\right)n\right]2

where

a

is some finite radius. This solution decays at infinity for any

n

, but becomes infinite at the origin for

n<1

, becomes finite at the origin for

n=1

and becomes zero at the origin for

n>1

. Walker also proposed two more solutions in his 1915 paper.

Radially symmetric forms

If the system to be studied is radially symmetric, then the equation in

n

dimension becomes

\psi''+

n-1
r

\psi'+λe\psi=0

where

r

is the distance from the origin. With the boundary conditions

\psi'(0)=0,\psi(1)=0

and for

λ\geq0

, a real solution exists only for

λ\in[0,λc]

, where

λc

is the critical parameter called as Frank-Kamenetskii parameter. The critical parameter is

λc=0.8785

for

n=1

,

λc=2

for

n=2

and

λc=3.32

for

n=3

. For

n=1, 2

, two solution exists and for

3\leqn\leq9

infinitely many solution exists with solutions oscillating about the point

λ=2(n-2)

. For

n\geq10

, the solution is unique and in these cases the critical parameter is given by

λc=2(n-2)

. Multiplicity of solution for

n=3

was discovered by Israel Gelfand in 1963 and in later 1973 generalized for all

n

by Daniel D. Joseph and Thomas S. Lundgren.[7]

The solution for

n=1

that is valid in the range

λ\in[0,0.8785]

is given by

\psi=-2ln

-\psim/2
\left[e

\cosh\left(

\sqrt{λ
}e^r\right)\right]

where

\psim=\psi(0)

is related to

λ

as
\psim/2
e

=\cosh\left(

\sqrt{λ
}e^\right).

The solution for

n=2

that is valid in the range

λ\in[0,2]

is given by

\psi=ln\left[

\psim
64e
(λ
\psim
e
r2+8)2

\right]

where

\psim=\psi(0)

is related to

λ

as

(λ

\psim
e

+8)2-64

\psim
e

=0.

Notes and References

  1. Liouville, J. "Sur l’équation aux différences partielles
    d2logλ
    dudv

    \pm

    λ
    2a2

    =0

    ." Journal de mathématiques pures et appliquées (1853): 71–72. http://sites.mathdoc.fr/JMPA/PDF/JMPA_1853_1_18_A3_0.pdf
  2. Bratu, G. "Sur les équations intégrales non linéaires." Bulletin de la Société Mathématique de France 42 (1914): 113–142.http://archive.numdam.org/article/BSMF_1914__42__113_0.pdf
  3. Gelfand, I. M. "Some problems in the theory of quasilinear equations." Amer. Math. Soc. Transl 29.2 (1963): 295–381. http://www.mathnet.ru/links/aa75c5d339030f17940afb64e17793d8/rm7290.pdf
  4. Richardson, Owen Willans. The emission of electricity from hot bodies. Longmans, Green and Company, 1921.
  5. Bateman, Harry. "Partial differential equations of mathematical physics." Partial Differential Equations of Mathematical Physics, by H. Bateman, Cambridge, UK: Cambridge University Press, 1932 (1932).
  6. Walker, George W. "Some problems illustrating the forms of nebulae." Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 91.631 (1915): 410-420.https://www.jstor.org/stable/pdf/93512.pdf?refreqid=excelsior%3Af4a4cc9656b8bbd9266f9d32587d02b1
  7. Joseph, D. D., and T. S. Lundgren. "Quasilinear Dirichlet problems driven by positive sources." Archive for Rational Mechanics and Analysis 49.4 (1973): 241-269.