*Plot of Bessel function of the first kind, J*_{α}(x), for integer orders α = 0, 1, 2.

Among the most important special functions is the Bessel function. In the field of solar physics, for example, it's of inestimable importance in the analysis of solar magnetic fields and their evolution. One very key equation is the

*Helmholtz, viz.*
1/ r [¶/ ¶r ( r ¶/ ¶r)]

**B**+ (a)^{2 }**B**= 0
where r is the radial coordinate,

B

**B**the magnetic field intensity, and a a quantity called the "force free parameter". Then the axially symmetric (i.e.- in cylindrical coordinates r, z, q) Bessel function solutions areB

_{z}(r) = B_{o}J_{o}(a r)
B

_{q}(r) = B_{o}J_{1}(ar)
where the axial (top) and azimuthal magnetic field components are given, respectively, and J

_{o}(a r) is a Bessel function of the first kind, order zero and J_{1}(ar) is a Bessel function of the first kind, order unity. (See graphs at the top for Bessel functions of the orders 0, 1 and 2).
The Bessel functions are mathematically defined (cf. Menzel,

*'Mathematical Physics'*, 1961, p. 204):
J

_{m}(x) = (1/ 2^{m}m!) x^{m}[1 - x^{2}/ 2^{2}1! (m + 1) + x^{4}/ 2^{4}2! (m + 1) (m + 2) - ….(-1)^{j}x^{2j}/ 2^{2j}j! (m + 1) (m + 2)…(m + j) + …]
which we terminate with second order terms.

For m = 0 and m = 1 forms one gets:

J

_{o}(x) = 1 - x^{2}/ 2^{2}(1!)^{2}+ x^{4}/ 2^{4}(2!)^{2}- x^{6}/ 2^{6}(3!)^{2}+ ......
and:

J

_{1}(x) = x/ 2 - x^{3}/ 2^{3}·1! 2! + x^{5}/ 2^{5}·2!3! - x^{7}/ 2^{7}·3!4! - .....
The equations in B

_{z}(r), B_{q}(r), with the special Bessel functions at root, are critical in describing the respective magnetic fields for a magnetic tube. For a cylindrical magnetic flux tube (such as a sunspot represents viewed in cross-section) the “twist” is defined:
T(r) = (L B

_{q}(r))/ (r B_{z}(r))
Where L denotes the length of the sunspot-flux tube dipole and r, the radius. If the twist exceeds 2p then the magnetic configuration may be approaching instability and a solar flare.

__Problems for the Math Maven__:

1) Compute the Bessel functions for J

_{o}(x) and J_{1}(x) with x = 1 and then compare with the values obtained from the graph shown at top.
2) Find the twist in a solar loop (take it to be a magnetic tube) if: B

_{q}(r) = 0. 1T and B_{z}(r) = 0.2T. Take the radius of the tube to be r = 10^{4}km and the length L = 10^{8}m. Is the tube kink unstable or not? (Kink instability is said to obtain when: T(r) > 2p)_{}

^{}

## No comments:

Post a Comment