User:Tohline/MathProjects/EigenvalueProblemN1

Find Analytic Solutions to an Eigenvalue Problem

Note from J. E. Tohline to Students with Good Mathematical Skills: This is one of a set of well-defined research problems that are being posed, in the context of this online H_Book, as challenges to young, applied mathematicians. The astronomy community's understanding of the Structure, Stability, and Dynamics of stars and galaxies would be strengthened if we had, in hand, a closed-form analytic solution to the problem being posed here. A solution can be obtained numerically with relative ease, but here the challenge is to find a closed-form analytic solution. As is true with most meaningful scientific research projects, it is not at all clear whether this problem has such a solution. In my judgment, however, it seems plausible that a closed-form solution can be discovered and such a solution would be of sufficient interest to the astronomical community that it would likely be publishable in a professional astronomy or physics journal. At the very least, this project offers an opportunity for a graduate student, an undergraduate, or even a talented high-school student (perhaps in connection with a mathematics science fair project?) to hone her/his research skills in applied mathematics. Also, I would be thrilled to include a solution to this problem — along with full credit to the solution's author — as a chapter in this online H_Book. Having retired from LSU, I am not in a position to financially support or formally advise students who are in pursuit of a higher-education degree. I would nevertheless be interested in sharing my expertise — and, perhaps, developing a collaborative relationship — with any individual who are interested in pursuing an answer to the mathematical research problem that is being posed here.

The Challenge

Formally, this is an eigenvalue problem. Find one or more analytic expression(s) for the (eigen)function, <math>~\mathcal{G}_\sigma(x)</math> — and, simultaneously, the unknown value of the (eigen)frequency, <math>~\sigma</math> — that satisfies the following 2^nd-order, ordinary differential equation:

<math> (x^2\sin x ) \frac{d^2\mathcal{G}_\sigma}{dx^2} + 2 \biggl[ x \sin x + x^2 \cos x \biggr] \frac{d\mathcal{G}_\sigma}{dx} + \biggl[ \sigma^2 x^3 - 2\alpha ( \sin x - x\cos x ) \biggr] \mathcal{G}_\sigma = 0 \, , </math>

where, <math>~\alpha</math> is a known constant. The desired functional solution is subject to the following two boundary conditions: <math>~\mathcal{G}_\sigma = 0</math> and <math>~d\mathcal{G}_\sigma/dx = 0</math> at <math>~x = 0</math>. Note that, in the context of astrophysical discussions, the interval of <math>~x</math> that is of particular interest is <math>0 \le x \le \pi</math>.

Analogous Problem with Known Analytic Solutions

Here is an analogous problem whose analytic solution is known. Anyone interested in tackling the challenge, provided above, should study — and even extend — the known set of solutions of this analogous problem. This exercise should provide at least partial preparation for addressing the above challenge.

Statement of the Problem

As above, the task here is to find one or more analytic expression(s) for the (eigen)function, <math>~\mathcal{F}_\sigma(x)</math> — and, simultaneously, the unknown value of the (eigen)frequency, <math>~\sigma</math> — that satisfies the following 2^nd-order, ordinary differential equation:

<math> (1 - x^2) \frac{d^2 \mathcal{F}_\sigma}{dx^2} + \frac{4}{x}\biggl[1 - \frac{3}{2}x^2 \biggr] \frac{d\mathcal{F}_\sigma}{dx} + \biggl[3\sigma^2 - 2 \alpha \biggr] \mathcal{F}_\sigma = 0 . </math>

Try a Polynomial Expression for the Eigenfunction

Let's guess that the proper eigenfunction is a polynomial expression in <math>~x</math>. Specifically, let's try a solution of the form,

<math>\mathcal{F}_\sigma = a + bx + cx^2 + dx^3 + fx^4 + gx^5 + \cdots </math>

truncated at progressively higher- and higher-order terms.

Lowest-order mode (Mode 0)

Try,

<math>\mathcal{F} = a \, ,</math>

in which case,

<math>~\frac{d\mathcal{F}}{dx}</math>

<math>~=</math>

<math>~0 \, ,</math>

and,

<math>~\frac{d^2\mathcal{F}}{dx^2}</math>

<math>~=</math>

<math>~0 \, .</math>

So, the governing 2^nd-order ODE reduces to,

<math>~(3\sigma^2 - 2 \alpha ) a</math>

<math>~=</math>

<math>~0 \, ,</math>

which will be satisfied as long as, <math>~\sigma = (2\alpha/3)^{1/2} \, .</math> We conclude, therefore, that the eigenvector defining the lowest-order (the simplest) solution to the governing ODE has an eigenfunction given by,

<math>~\mathcal{F}_0</math>

<math>~=</math>

<math>~a = \mathrm{constant} \, ,</math>

with a corresponding eigenfrequency whose value is,

<math>~\sigma_0</math>

<math>~=</math>

<math>~\biggl(\frac{2\alpha}{3} \biggr)^{1/2} \, .</math>

Second Guess

Try,

<math>\mathcal{F} = a + bx \, ,</math>

in which case,

<math>~\frac{d\mathcal{F}}{dx}</math>

<math>~=</math>

<math>~b \, ,</math>

and,

<math>~\frac{d^2\mathcal{F}}{dx^2}</math>

<math>~=</math>

<math>~0 \, .</math>

Plugging this trial eigenfunction into the governing 2^nd-order ODE gives,

<math>~0</math>	<math>~=</math>	<math>~\frac{4}{x}\biggl[1 - \frac{3}{2}x^2 \biggr] b + \biggl[3\sigma^2 - 2 \alpha \biggr] (a + bx)</math>
	<math>~=</math>	<math>~\frac{4b}{x} + (3\sigma^2 - 2 \alpha -6) bx + (3\sigma^2 - 2 \alpha ) a \, .</math>

But this expression can be satisfied for all values of <math>~x</math> only if <math>~b = 0</math>, in which case the trial eigenfunction reduces to the earlier solution, <math>~\mathcal{F}_0</math>. We conclude, therefore, that our "second guess" does not generate a new solution to this eigenfunction problem.

Third Guess

Try,

<math>\mathcal{F} = a + cx^2\, ,</math>

in which case,

<math>~\frac{d\mathcal{F}}{dx}</math>

<math>~=</math>

<math>~2cx \, ,</math>

and,

<math>~\frac{d^2\mathcal{F}}{dx^2}</math>

<math>~=</math>

<math>~2c \, .</math>

Plugging this trial eigenfunction into the governing 2^nd-order ODE therefore gives,

<math>~0</math>	<math>~=</math>	<math>~2c(1 - x^2) + 8c(1 - \frac{3}{2}x^2 ) + (3\sigma^2 - 2 \alpha ) (a + cx^2)</math>
	<math>~=</math>	<math>~ \biggl[10c + (3\sigma^2 - 2 \alpha )a\biggr] + \biggl[ (3\sigma^2 - 2 \alpha ) -14 \biggr]cx^2 \, .</math>

This relation will be satisfied for all values of <math>~x</math> if both expressions inside the square brackets are simultaneously zero, that is, if,

<math>~10c + (3\sigma^2 - 2 \alpha )a</math>

<math>~=</math>

<math>~0 \, ,</math>

and, simultaneously,

<math>~(3\sigma^2 - 2 \alpha ) -14 </math>

<math>~=</math>

<math>~0 \, .</math>

• Mode 1:

<math>x_1 = a + b\chi_0^2</math>, in which case,

<math> \frac{dx}{d\chi_0} = 2b\chi_0; ~~~~ \frac{d^2 x}{d\chi_0^2} = 2b; </math>

<math> \frac{1}{(1 - \chi_0^2)} \biggl\{ 2b (1 - \chi_0^2) + 8b \biggl[1 - \frac{3}{2}\chi_0^2 \biggr] + A_1 \biggl(1 + \frac{b}{a}\chi_0^2 \biggr) \biggr\} = 0 , </math>

where,

<math> A_1 \equiv \frac{a}{\gamma_\mathrm{g}}\biggl[ \biggl( \frac{3}{2\pi G\rho_c} \biggr) \omega_1^2+ 2(4 - 3\gamma_\mathrm{g}) \biggr] . </math>

Therefore,

<math> (A_1 + 10b) + \biggl[ \biggl(\frac{b}{a}\biggr) A_1 - 14b \biggr] \chi_0^2 = 0 , </math>

<math> \Rightarrow ~~~~~ A_1 = - 10b ~~~~~\mathrm{and} ~~~~~ A_1 = 14a </math>

<math> \Rightarrow ~~~~~ \frac{b}{a} = -\frac{7}{5} ~~~~~\mathrm{and} ~~~~~ \frac{A_1}{a} = 14 = \frac{1}{\gamma_\mathrm{g}}\biggl[ \biggl( \frac{3}{2\pi G\rho_c} \biggr) \omega_1^2+ 2(4 - 3\gamma_\mathrm{g}) \biggr] .

</math>

Hence,

<math> \biggl( \frac{3}{2\pi G\rho_c} \biggr) \omega_1^2 = 20\gamma_\mathrm{g} -8 </math>

<math> \Rightarrow ~~~~~ \omega_1^2 = \frac{2}{3}\biggl( 4\pi G\rho_c \biggr) (5\gamma_\mathrm{g} -2) </math>

and, to within an arbitrary normalization factor,

<math> x_1 = 1 - \frac{7}{5}\chi_0^2 . </math>

Astrophysical Context

A star that is composed entirely of a gas for which the pressure varies as the square of the gas density will have an equilibrium structure that is defined by, what astronomers refer to as, an <math>~n=1</math> polytrope. Inside such an equilibrium structure, the density of the gas will vary with radial position according to the expression,

<math>~\frac{\rho}{\rho_c} = \frac{\sin x}{x} \, ,</math>

where, <math>~\rho_c</math> is the density at the center of the star, and,

<math>~x</math>

<math>~\equiv</math>

<math>~\pi\biggl(\frac{r}{R}\biggr) \, ,</math>

where, <math>~R</math> is the radius of the equilibrium star. Notice that, according to this expression, the density will drop to zero when <math>~r = R</math>, in which case, <math>~x = \pi</math>. If a star of this type is nudged out of equilibrium — for example, squeezed slightly — in such a way that it maintains its spherical symmetry, the star will begin to undergo periodic, radial oscillations about its original equilibrium radius. The 2^nd-order ODE whose solution is being sought in the above challenge is the equation that describes the behavior of these oscillations. In particular, the function,

<math>~\mathcal{G}_\sigma(x) \equiv \frac{\delta x}{x}</math>

describes the relative amplitude of the oscillation as a function of position, <math>~x</math>, within the star, and <math>~\sigma</math> gives the frequency of the oscillation.

Related Discussions

• Linear Stability Analysis
• Radial Pulsation Modes of Uniform-Density Spheres

© 2014 - 2021 by Joel E. Tohline |   H_Book Home   |   YouTube   | Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS | Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation