User:Tohline/Appendix/Ramblings/PPToriPt2

From VistrailsWiki
Jump to navigation Jump to search

Stability Analyses of PP Tori (Part 2)

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

This is a direct extension of our Part 1 discussion. Here we continue our effort to check the validity of the Blaes85 eigenvector. The relevant reference is:

Start From Scratch

Basic Equations from Blaes85

  Blaes85

Eq. No.

<math>~(\beta\eta)^2</math>

<math>~=</math>

<math>~x^2(1+xb) \, ;</math>

(2.6)

<math>~b</math>

<math>~\equiv</math>

<math>~3\cos\theta - \cos^3\theta \, ;</math>

(2.6)

<math>~f</math>

<math>~=</math>

<math>~1-\eta^2 \, .</math>

(2.5)

  Blaes85

Eq. No.

<math>~LHS \equiv \hat{L}W</math>

<math>~=</math>

<math> ~fx^2 \cdot \frac{\partial^2 W}{\partial x^2} + f \cdot \frac{\partial^2 W}{\partial \theta^2} + \biggl[ \frac{fx(1-2x\cos\theta)}{(1-x\cos\theta)} + nx^2\cdot \frac{\partial f}{\partial x}\biggr]\frac{\partial W}{\partial x} </math>

 

 

 

<math> + \biggl[ \frac{fx\sin\theta}{(1-x\cos\theta)} + n\cdot \frac{\partial f}{\partial \theta}\biggr]\frac{\partial W}{\partial \theta} + \biggl[ \frac{2nx^2m^2}{\beta^2(1-x\cos\theta)^4} - \frac{m^2 x^2 f}{(1-x\cos\theta)^2} \biggr]W </math>

(4.2)

<math>~RHS</math>

<math>~=</math>

<math> ~-\frac{2nm^2}{\beta^2} \cdot (\beta\eta)^2 \biggl[ M \biggl(\frac{\nu}{m}\biggr)^2 + \frac{N}{m} \biggl(\frac{\nu}{m}\biggr)\biggr] W </math>

(4.1)

 

<math>~=</math>

<math> ~-\frac{2nm^2}{\beta^2} \biggl[ x^2 \biggl(\frac{\nu}{m}\biggr)^2 + \frac{2x^2}{(1-x\cos\theta)^2} \biggl(\frac{\nu}{m}\biggr)\biggr] W </math>

(4.2)

<math>~\frac{W}{A_{00}}</math>

<math>~=</math>

<math> ~1 + \beta^2 m^2 \biggl\{ 2\eta^2\cos^2\theta - \frac{3\eta^2}{4(n+1)} - \frac{(4n+1)}{4(n+1)^2} ~\pm~i~\biggl[ \frac{2^3\cdot 3}{(n+1)}\biggr]^{1/2} \eta\cos\theta \biggr\} </math>

(4.13)

<math>~\frac{\nu}{m}</math>

<math>~=</math>

<math> ~-1 ~\pm ~ i~\biggl[ \frac{3}{2(n+1)} \biggr]^{1/2} \beta </math>

(4.14)

Our Manipulation of These Equations

<math>~\Lambda \equiv \frac{2^2(n+1)^2}{m^2}\biggl[\frac{W}{A_{00}}-1\biggr]</math>

<math>~=</math>

<math>~\beta^2 \biggl\{ 2^3(n+1)^2 \eta^2\cos^2\theta - 3\eta^2(n+1)^2 - (4n+1) ~\pm~i~[ 2^7\cdot 3(n+1)^3 ]^{1/2} \eta\cos\theta \biggr\} </math>

 

<math>~=</math>

<math>~- (4n+1)\beta^2 + (\beta\eta)^2 (n+1)^2[ 2^3 \cos^2\theta - 3] ~\pm~i~\beta [ 2^7\cdot 3(n+1)^3 ]^{1/2} (\beta\eta) \cos\theta \, ; </math>

<math>~\Rightarrow~~~~\frac{W}{A_{00}} </math>

<math>~=</math>

<math>~1+ \biggl[ \frac{m}{2(n+1)} \biggr]^2 \Lambda </math>


<math>~\frac{LHS}{A_{00}} </math>

<math>~=</math>

<math>~\biggl[ \frac{m}{2(n+1)} \biggr]^2 f ~\biggl[ x^2 \cdot \frac{\partial^2 \Lambda}{\partial x^2} + \frac{\partial^2 \Lambda}{\partial \theta^2}\biggr] + \biggl[ \frac{m}{2(n+1)} \biggr]^2\biggl[ \frac{fx(1-2x\cos\theta)}{(1-x\cos\theta)} + nx^2\cdot \frac{\partial f}{\partial x}\biggr]\frac{\partial \Lambda}{\partial x} </math>

 

 

<math> + \biggl[ \frac{m}{2(n+1)} \biggr]^2\biggl[ \frac{fx\sin\theta}{(1-x\cos\theta)} + n\cdot \frac{\partial f}{\partial \theta}\biggr]\frac{\partial \Lambda}{\partial \theta} + \biggl[ \frac{2nx^2m^2}{\beta^2(1-x\cos\theta)^4} - \frac{m^2 x^2 f}{(1-x\cos\theta)^2} \biggr]\biggl\{1+ \biggl[ \frac{m}{2(n+1)} \biggr]^2 \Lambda\biggr\} </math>

 

<math>~=</math>

<math>~\biggl[ \frac{m}{2(n+1)} \biggr]^2 f \biggl\{ ~\biggl[ x^2 \cdot \frac{\partial^2 \Lambda}{\partial x^2} + \frac{\partial^2 \Lambda}{\partial \theta^2}\biggr] + \biggl[ \frac{x(1-2x\cos\theta)}{(1-x\cos\theta)} \biggr]\frac{\partial \Lambda}{\partial x} + \biggl[ \frac{x\sin\theta}{(1-x\cos\theta)} \biggr]\frac{\partial \Lambda}{\partial \theta} - \biggl[ \frac{m^2 x^2 }{(1-x\cos\theta)^2} \biggr] \biggl[ \frac{2^2(n+1)^2}{m^2} + \Lambda\biggr]\biggr\} </math>

 

 

<math> + n\biggl[ \frac{m}{2(n+1)} \biggr]^2 \biggl\{ x^2\cdot \frac{\partial f}{\partial x}\cdot \frac{\partial \Lambda}{\partial x} ~+~ \frac{\partial f}{\partial \theta}\cdot \frac{\partial \Lambda}{\partial \theta} ~+~ \biggl[ \frac{2x^2m^2}{\beta^2(1-x\cos\theta)^4} \biggr]\biggl[ \frac{2^2(n+1)^2}{m^2} + \Lambda\biggr]\biggr\} </math>

 

<math>~=</math>

<math>~\frac{x^2 f}{(1-x\cos\theta)^2} \biggl[ \frac{m}{2(n+1)} \biggr]^2 \biggl\{ ~(1-x\cos\theta)^2\biggl[ \frac{\partial^2 \Lambda}{\partial x^2} + \frac{1}{x^2}\cdot \frac{\partial^2 \Lambda}{\partial \theta^2}\biggr] + \frac{(1-x\cos\theta)}{x} \biggl[ (1-2x\cos\theta) \frac{\partial \Lambda}{\partial x} + \sin\theta\cdot \frac{\partial \Lambda}{\partial \theta} \biggr] - [ 2^2(n+1)^2 + m^2\Lambda ]\biggr\} </math>

 

 

<math> + ~\frac{x^2 n}{\beta^2(1-x\cos\theta)^4} \biggl[ \frac{m}{2(n+1)} \biggr]^2 \biggl\{\beta^2 (1-x\cos\theta)^4\biggl[ \frac{\partial f}{\partial x}\cdot \frac{\partial \Lambda}{\partial x} ~+~ \frac{1}{x^2}\cdot \frac{\partial f}{\partial \theta}\cdot \frac{\partial \Lambda}{\partial \theta} \biggr] ~+~ [ 2^3(n+1)^2 + 2m^2\Lambda ]\biggr\} \, . </math>

Also,

<math>~\frac{RHS}{A_{00}}</math>

<math>~=</math>

<math> ~-\frac{2n x^2}{\beta^2(1-x\cos\theta)^2} \biggl[ \frac{m}{2(n+1)} \biggr]^2 \biggl[ (1-x\cos\theta)^2\biggl(\frac{\nu}{m}\biggr)^2 + 2\biggl(\frac{\nu}{m}\biggr)\biggr] [ 2^2(n+1)^2 + m^2\Lambda ] </math>

 

<math>~=</math>

<math> ~-\frac{x^2n}{\beta^2(1-x\cos\theta)^4} \biggl[ \frac{m}{2(n+1)} \biggr]^2 \biggl[ (1-x\cos\theta)^4\biggl(\frac{\nu}{m}\biggr)^2 + 2(1-x\cos\theta)^2\biggl(\frac{\nu}{m}\biggr)\biggr] [ 2^3(n+1)^2 + 2m^2\Lambda ] \, . </math>

Putting the two together implies,

<math>~0</math>

<math>~=</math>

<math>~\frac{1}{x^2}\biggl[\frac{LHS}{A_{00}} - \frac{RHS}{A_{00}}\biggr]\biggl[ \frac{2(n+1)}{m} \biggr]^2 (1-x\cos\theta)^4</math>

 

<math>~=</math>

<math>~f (1-x\cos\theta)^2 \biggl\{ ~(1-x\cos\theta)^2\biggl[ \frac{\partial^2 \Lambda}{\partial x^2} + \frac{1}{x^2}\cdot \frac{\partial^2 \Lambda}{\partial \theta^2}\biggr] + \frac{(1-x\cos\theta)}{x} \biggl[ (1-2x\cos\theta) \frac{\partial \Lambda}{\partial x} + \sin\theta\cdot \frac{\partial \Lambda}{\partial \theta} \biggr] - [ 2^2(n+1)^2 + m^2\Lambda ]\biggr\} </math>

 

 

<math> + ~\frac{n}{\beta^2} \biggl\{ (1-x\cos\theta)^4\biggl[ \frac{\partial \Lambda}{\partial x} \cdot \frac{\partial (\beta^2 f)}{\partial x} ~+~ \frac{\partial \Lambda}{\partial \theta} \cdot \frac{\partial (\beta^2 f/x^2)}{\partial \theta} \biggr] ~+~ \biggl[ (1-x\cos\theta)^4\biggl(\frac{\nu}{m}\biggr)^2 + 2(1-x\cos\theta)^2\biggl(\frac{\nu}{m}\biggr)+ 1 \biggr] [ 2^3(n+1)^2 + 2m^2\Lambda ] \biggr\} \, . </math>

The first line of this governing, two-line expression contains the function, <math>~f</math>, as a leading factor, while the leading factor in the second line is the ratio, <math>~n/\beta^2</math>. Presumably the three terms (hereafter, TERM1, TERM2, & TERM3, respectively) inside the curly brackets on the first line must cancel — to a sufficiently high order in <math>~x</math> — and, independently, the two terms (hereafter, TERM4 & Term5, respectively) inside the curly brackets on the second line must cancel. Furthermore, these cancellations must occur separately for the real parts and the imaginary parts of each bracketed expression.

Evaluating various terms using the parameter set,    <math>~(n, \theta, x/\beta) = (1, \tfrac{\pi}{3}, \tfrac{1}{4})</math>    as begun in our "Part 1" analysis, we have:

TERM4

<math>~=</math>

<math>~ -x \ell^4\biggl[ (2+3xb)\cdot \frac{\partial\Lambda}{\partial x} - 3\sin^3\theta \cdot \frac{\partial\Lambda}{\partial \theta} \biggr] </math>

 

<math>~=</math>

<math>~ -x \ell^4\biggl[ (2+3xb)\cdot [~1.515625000~\pm~i~36.23373732 ~\beta] - 3\sin^3\theta \cdot [~-2.388335684~\pm~i~(-1)15.36617018 ~\beta] \biggr] </math>

 

<math>~=</math>

<math>~ -x\ell^4 [~9.248046874~\pm~i~139.7753772~\beta] </math>

TERM5 (Case B)

<math>~=</math>

<math>~ \biggl[ \ell^4 [1-0.75\beta^2~\pm~i~(-1)\sqrt{3}\beta] +2\ell^2[ -1~\pm~i~\sqrt{0.75}\beta ] + 1 \biggr] \cdot \biggl[~2^5 + 2\cancelto{1}{m^2}[~- 5\beta + 0.167968750~\pm~i~8.031189202 ~\beta]~\biggr] </math>

 

<math>~=</math>

<math>~ \biggl[ \ell^4 [1-0.75\beta^2] - 2\ell^2 + 1 \biggr] \cdot \biggl[~[2^5 - 10\beta + (2)0.167968750]~\pm~i~[(2)8.031189202 ~\beta]~\biggr] </math>

 

 

<math>~ \pm~\sqrt{3}\beta\biggl[ \ell^2-~\ell^4 \biggr] \cdot \biggl[~i~ [2^5 - 10\beta + (2)0.167968750]~-~[(2)8.031189202 ~\beta]~\biggr] </math>

 

<math>~=</math>

<math>~ \biggl[ \ell^4 [1-0.75\beta^2] - 2\ell^2 + 1 \biggr] \cdot \biggl[~[2^5 - 10\beta + (2)0.167968750]\biggr] \pm~(-1)\sqrt{3}\beta\biggl[ \ell^2-~\ell^4 \biggr] \cdot \biggl[[(2)8.031189202 ~\beta]~\biggr] </math>

 

 

<math>~ \pm~i~\biggl\{\biggl[ \ell^4 [1-0.75\beta^2] - 2\ell^2 + 1 \biggr] \cdot \biggl[[(2)8.031189202 ~\beta]~\biggr] +~\sqrt{3}\beta\biggl[ \ell^2-~\ell^4 \biggr] \cdot \biggl[~ [2^5 - 10\beta + (2)0.167968750]~\biggr] \biggr\} \, . </math>

Evaluating this TERM5 expression for the case of <math>~\beta = 1</math>, we have,


TERM5 (Case B)

<math>~=</math>

<math>~ \biggl[ 0.25\ell^4 - 2\ell^2 + 1 \biggr] \cdot \biggl[~[2^5 - 10 + (2)0.167968750]\biggr] \pm~(-1)\sqrt{3}\biggl[ \ell^2-~\ell^4 \biggr] \cdot \biggl[[(2)8.031189202 ]~\biggr] </math>

 

 

<math>~ \pm~i~\biggl\{\biggl[ \ell^4 [1-0.75] - 2\ell^2 + 1 \biggr] \cdot \biggl[[(2)8.031189202]~\biggr] +~\sqrt{3}\biggl[ \ell^2-~\ell^4 \biggr] \cdot \biggl[~ [2^5 - 10 + (2)0.167968750]~\biggr] \biggr\} </math>

 

<math>~=</math>

<math>~ [ -0.38470459 ] \cdot [22.3359375] \pm~(-1)[ ~0.31080502 ] \cdot [~16.0623784 ~] </math>

 

 

<math>~ \pm~i~\biggl\{[ -0.38470459 ] \cdot [~16.0623784 ~] +~[ ~0.31080502 ] \cdot [22.3359375]\biggr\} </math>

 

<math>~=</math>

<math>~[~-13.58500545~] \pm~i~[~0.76285080~] \, . </math>

Testing for Expected Cancellations

Note first that, adopting the shorthand notation,

<math>~\ell</math>

<math>~\equiv</math>

<math>~(1-x\cos\theta)</math>

<math>~\Rightarrow ~~~~\ell^2</math>

<math>~=</math>

<math>~1-2\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + \beta^2 \biggl(\frac{x}{\beta}\biggr)^2\cos^2\theta + \mathcal{O}(\beta^3) \, ;</math>

<math>~\ell^3</math>

<math>~=</math>

<math>~1-3\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + 3\beta^2 \biggl(\frac{x}{\beta}\biggr)^2\cos^2\theta + \mathcal{O}(\beta^3) \, ;</math>

<math>~\ell^4</math>

<math>~=</math>

<math>~1-4\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + 6\beta^2\biggl(\frac{x}{\beta}\biggr)^2\cos^2\theta + \mathcal{O}(\beta^3) \, .</math>


Real Parts

<math>~\mathrm{Re}\biggl[\frac{\mathrm{TERM4}}{\ell^4}\biggr]</math>

<math>~=</math>

<math>~ \biggl\{ (n+1)[2^3(n+1)\cos^2\theta -3]x(2+3xb)\biggr\} \cdot \biggl[ -x(2+3xb) \biggr] </math>

 

 

<math>~ +~ (n+1)\sin\theta \biggl\{ -2^4 (n+1) (\beta\eta)^2 \cos\theta + 3x^3 \sin^2\theta \biggl[3 - 2^3(n+1)\cos^2\theta \biggr] \biggr\} \cdot \biggl[ 3x\sin^3\theta \biggr] </math>

 

<math>~=</math>

<math>~ -~(n+1)[2^3(n+1)\cos^2\theta -3]x^2(2+3xb)^2 </math>

 

 

<math>~ -~ 3x^3(n+1)\sin^4\theta \biggl\{ 2^4 (n+1) (1+xb) \cos\theta + 3x \sin^2\theta [2^3(n+1)\cos^2\theta -3] \biggr\} </math>

 

<math>~=</math>

<math>~ -~x^2 \cdot 2^2 (n+1)[2^3(n+1)\cos^2\theta -3]\biggl(1+\frac{3xb}{2}\biggr)^2 </math>

 

 

<math>~ -~ x^3 \cdot 2^4\cdot 3(n+1)^2 \cos\theta\sin^4\theta (1+xb) ~-~x^4\cdot 3^2(n+1)\sin^6\theta [2^3(n+1)\cos^2\theta -3] \, .</math>

 

<math>~=</math>

<math>~ -x\biggl\{~x[~18.37695315~] + x^2[~72.5625~] + x^3[~7.59375~]~~\biggr\} = -x[~9.24804688~]\, . </math>

<math>~\mathrm{Re}\biggl[\mathrm{TERM5}\biggr]</math>

<math>~=</math>

<math>~ \mathrm{Re}\biggl[ \ell^4\biggl(\frac{\nu}{m}\biggr)^2 + 2\ell^2\biggl(\frac{\nu}{m}\biggr)+ 1 \biggr] \cdot \mathrm{Re}[ 2^3(n+1)^2 + 2m^2\Lambda ] -\mathrm{Im}\biggl[ \ell^4\biggl(\frac{\nu}{m}\biggr)^2 + 2\ell^2\biggl(\frac{\nu}{m}\biggr)+ 1 \biggr] \cdot \mathrm{Im}[ 2^3(n+1)^2 + 2m^2\Lambda ] </math>

Case B:

<math>~=</math>

<math>~ \biggl\{ \ell^4\biggl[1-\frac{3\beta^2}{2(n+1)}\biggr] + 2\ell^2\biggl(-1\biggr)+ 1 \biggr\} \cdot \biggl\{ 2^3(n+1)^2 + 2m^2\biggl[ ~- (4n+1)\beta^2 + (n+1)^2(2^3 \cos^2\theta - 3) x^2(1+xb)\biggr] \biggr\} </math>

 

 

<math>~ -~\biggl\{ \ell^4(-1)\biggl[\frac{2\cdot 3\beta^2}{(n+1)}\biggr]^{1/2} + 2\ell^2\biggl[ \frac{3\beta^2}{2(n+1)}\biggr]^{1/2} \biggr\} \cdot 2m^2\beta [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta \cdot x(1+xb)^{1/2} </math>

 

<math>~=</math>

<math>~ \biggl\{1 - 2\ell^2 + \ell^4-\frac{3\beta^2\ell^4}{2(n+1)} \biggr\} \cdot \biggl\{ \biggl[ 2^3(n+1)^2 - 2m^2(4n+1)\beta^2\biggr] + x^2\cdot 2m^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr\} </math>

 

 

<math>~ -~x\beta^2 \cdot m^2[\ell^2 - \ell^4 ] \cdot [ 2^{10}\cdot 3^2(n+1)^2 ]^{1/2} \cos\theta (1+xb)^{1/2} \, . </math>

When added together, we obtain,

<math>~\mathrm{Re}[\mathrm{TERM4} + \mathrm{TERM5}]</math>

<math>~=</math>

<math>~ -~\beta^2 \biggl(\frac{x}{\beta}\biggr)^2 \ell^4\cdot 2^2 (n+1)[2^3(n+1)\cos^2\theta -3 ]\biggl(1+\frac{3xb}{2}\biggr)^2 </math>

 

 

<math>~ -~ \beta^3 \biggl(\frac{x}{\beta}\biggr)^3\ell^4\cdot 2^4\cdot 3(n+1)^2 \cos\theta\sin^4\theta (1+xb) ~-~\beta^4\biggl(\frac{x}{\beta}\biggr)^4 \ell^4\cdot 3^2(n+1)\sin^6\theta [2^3(n+1)\cos^2\theta-3] </math>

 

 

<math>~ +~\biggl\{1 - 2\ell^2 + \ell^4 \biggr\} \cdot \biggl\{ 2^3(n+1)^2 + 2m^2\beta^2\biggr[ - (4n+1) + \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr]\biggr\} </math>

 

 

<math>~ -~\frac{3\beta^2\ell^4}{2(n+1)} \biggl\{ 2^3(n+1)^2 + 2m^2\beta^2\biggr[ - (4n+1) + \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr]\biggr\} </math>

 

 

<math>~ -~\beta^3\biggl(\frac{x}{\beta}\biggr) \cdot m^2[\ell^2 - \ell^4 ] \cdot [ 2^{10}\cdot 3^2(n+1)^2 ]^{1/2} \cos\theta (1+xb)^{1/2} </math>

 

<math>~=</math>

<math>~ \beta^0 \cdot 2^3(n+1)^2\biggl\{1 - 2\ell^2 + \ell^4 \biggr\} </math>

 

 

<math>~ -~\beta^2 \cdot 2m^2 [ 1 - 2\ell^2 + \ell^4 ] \cdot \biggr[ (4n+1) - \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr] </math>

 

 

<math>~ -~\beta^2\ell^4 2^2\cdot 3 (n+1) + \beta^2 \biggl(\frac{x}{\beta}\biggr)^2 \ell^4\cdot 2^2 (n+1)[3 - 2^3(n+1)\cos^2\theta ]\biggl(1+\frac{3xb}{2}\biggr)^2 </math>

 

 

<math>~ -~\cancelto{0}{\beta^3}\biggl(\frac{x}{\beta}\biggr) \cdot m^2[\ell^2 - \ell^4 ] \cdot [ 2^{10}\cdot 3^2(n+1)^2 ]^{1/2} \cos\theta (1+xb)^{1/2} </math>

 

 

<math>~ -~ \cancelto{0}{\beta^3} \biggl(\frac{x}{\beta}\biggr)^3\ell^4\cdot 2^4\cdot 3(n+1)^2 \cos\theta\sin^4\theta (1+xb) ~-~\cancelto{0}{\beta^4}\biggl(\frac{x}{\beta}\biggr)^4 \ell^4\cdot 3^2(n+1)\sin^6\theta [2^3(n+1)\cos^2\theta-3] </math>

 

 

<math>~ +~\frac{3\cancelto{0}{\beta^4}\ell^4 m^2}{(n+1)} \biggr[ (4n+1) - \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr] </math>

 

<math>~\approx</math>

<math>~ \beta^0 \cdot 2^3(n+1)^2\biggl\{1 - 2\biggl[ 1-2\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + \beta^2 \biggl(\frac{x}{\beta}\biggr)^2\cos^2\theta + \cancelto{0}{\mathcal{O}(\beta^3)}\biggr] + \biggl[ 1-4\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + 6\beta^2\biggl(\frac{x}{\beta}\biggr)^2\cos^2\theta + \cancelto{0}{\mathcal{O}(\beta^3)} \biggr] \biggr\} </math>

 

 

<math>~ -~\beta^2 \cdot 2m^2 [ 1 - 2 + 1 ] \cdot \biggr[ (4n+1) - \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) (1+\cancelto{0}{x}b) \biggr] </math>

 

 

<math>~ -~\beta^2 2^2\cdot 3 (n+1) + \beta^2 \biggl(\frac{x}{\beta}\biggr)^2 2^2 (n+1)[3 - 2^3(n+1)\cos^2\theta ]\biggl(1+\frac{3\cancelto{0}{x}b}{2}\biggr)^2 </math>

 

<math>~\approx</math>

<math>~ \beta^0 \cdot 2^3(n+1)^2\biggl\{1 - 2+ 1 \biggr\} +~\beta^1 \biggl(\frac{x}{\beta}\biggr) \cdot 2^3(n+1)^2\biggl\{4\cos\theta -4\cos\theta \biggr\} </math>

 

 

<math>~ +~\beta^2 \biggl(\frac{x}{\beta}\biggr)^2 \cdot 2^5(n+1)^2 \cos^2\theta </math>

 

 

<math>~ -~\beta^2 \cdot 2m^2 [ 1 - 2 + 1 ] \cdot \biggr[ (4n+1) - \biggl(\frac{x}{\beta}\biggr)^2(n+1)^2(2^3 \cos^2\theta - 3) \biggr] </math>

 

 

<math>~ -~\beta^2 2^2\cdot 3 (n+1) \biggl[1 - \biggl(\frac{x}{\beta}\biggr)^2\biggr] - \beta^2 \biggl(\frac{x}{\beta}\biggr)^2 [2^5(n+1)^2\cos^2\theta ] </math>

 

<math>~=</math>

<math>~-~\beta^2 2^2\cdot 3 (n+1) \biggl[1 - \biggl(\frac{x}{\beta}\biggr)^2\biggr] \, .</math>

So we see that the coefficients of the lowest-order <math>(\beta^0 ~\mathrm{and} ~ \beta^1)</math> terms are zero, and the coefficient of the <math>~\beta^2</math> term is almost zero!

Imaginary Parts

<math>~\mathrm{Im}\biggl[\frac{\mathrm{TERM4}}{\ell^4}\biggr]</math>

<math>~=</math>

<math>~ \biggl\{ \beta\cos\theta [2^5\cdot 3 (n+1)^3]^{1/2} \cdot \frac{x(2+3xb)}{(\beta\eta)}\biggr\} \cdot \biggl[ -x(2+3xb) \biggr] </math>

 

 

<math>~ -~ \beta \sin\theta [2^7\cdot 3 (n+1)^3 (\beta\eta)^2]^{1/2}\biggl\{ 1 +\frac{3x^3}{2}\cdot\biggl[ \frac{\sin^2\theta \cos\theta}{(\beta\eta)^2} \biggr]\biggr\} \cdot \biggl[ 3x\sin^3\theta \biggr] </math>

 

<math>~=</math>

<math>~ -~x \cdot 2\beta\cos\theta [2^7\cdot 3 (n+1)^3]^{1/2} \cdot (1+xb)^{-1/2}\cdot \biggl(1+\frac{3xb}{2}\biggr)^2 </math>

 

 

<math>~ -~ x^2\cdot 3\beta \sin^4\theta [2^7\cdot 3 (n+1)^3 ]^{1/2} (1+xb)^{1/2} \biggl\{ 1 +\frac{3x}{2}\cdot\biggl[ \frac{\sin^2\theta \cos\theta}{(1+xb)} \biggr]\biggr\} </math>

 

<math>~=</math>

<math>~ -x\biggl\{~[~109.8335164~] + x[~2.16086768~]~\biggr\} </math>


<math>~\mathrm{Im}\biggl[\mathrm{TERM5}\biggr]</math>

<math>~=</math>

<math>~ \mathrm{Re}\biggl[ \ell^4\biggl(\frac{\nu}{m}\biggr)^2 + 2\ell^2\biggl(\frac{\nu}{m}\biggr)+ 1 \biggr] \cdot \mathrm{Im}[ 2^3(n+1)^2 + 2m^2\Lambda ] -\mathrm{Im}\biggl[ \ell^4\biggl(\frac{\nu}{m}\biggr)^2 + 2\ell^2\biggl(\frac{\nu}{m}\biggr)+ 1 \biggr] \cdot \mathrm{Re}[ 2^3(n+1)^2 + 2m^2\Lambda ] </math>

Case B:

<math>~=</math>

<math>~ x\cdot 2 \beta m^2 \biggl\{1 - 2\ell^2 + \ell^4 -\frac{3\beta^2\ell^4}{2(n+1)} \biggr\} \cdot [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta \cdot (1+xb)^{1/2} </math>

 

 

<math>~ -~\beta \biggl[ \frac{2\cdot 3}{(n+1)}\biggr]^{1/2} [\ell^2 -\ell^4] \cdot \biggl\{ \biggl[ 2^3(n+1)^2 ~- 2m^2(4n+1)\beta^2\biggr] + x^2 \cdot 2m^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr\} </math>

When added together, we obtain,

<math>~\mathrm{Im}[\mathrm{TERM4} + \mathrm{TERM5}]</math>

<math>~=</math>

<math>~ -~\beta^2 \biggl(\frac{x}{\beta}\biggr) \ell^4 \cos\theta [2^9\cdot 3 (n+1)^3]^{1/2} \cdot (1+xb)^{-1/2}\cdot \biggl(1+\frac{3xb}{2}\biggr)^2 </math>

 

 

<math>~ -~\cancelto{0}{\beta^3} \biggl(\frac{x}{\beta}\biggr)^2\cdot 3 \ell^4 \sin^4\theta [2^7\cdot 3 (n+1)^3 ]^{1/2} (1+xb)^{1/2} \biggl\{ 1 +\frac{3x}{2}\cdot\biggl[ \frac{\sin^2\theta \cos\theta}{(1+xb)} \biggr]\biggr\} </math>

 

 

<math>~+\beta^2 \biggl(\frac{x}{\beta}\biggr)\cdot 2 m^2 [1 - 2\ell^2 + \ell^4 ] \cdot [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta \cdot (1+xb)^{1/2} </math>

 

 

<math>~-\cancelto{0}{\beta^4} \biggl(\frac{x}{\beta}\biggr) \biggl[\frac{3 m^2\ell^4}{(n+1)} \biggr] \cdot [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta \cdot (1+xb)^{1/2} </math>

 

 

<math>~ -~\beta [ 2^7\cdot 3 (n+1)^3]^{1/2} [\ell^2 -\ell^4] </math>

 

 

<math>~ +~\cancelto{0}{\beta^3} \biggl[ \frac{2\cdot 3}{(n+1)}\biggr]^{1/2} [\ell^2 -\ell^4] \cdot \biggl[ 2m^2(4n+1) - \biggl(\frac{x}{\beta}\biggr)^2 2m^2(n+1)^2(2^3 \cos^2\theta - 3) (1+xb) \biggr] </math>

 

<math>~\approx</math>

<math>~ -~\beta^1 [ 2^7\cdot 3 (n+1)^3]^{1/2} \biggl\{ \biggl[ 1-2\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + \cancelto{0}{\mathcal{O}(\beta^2)} \biggr] - \biggl[ 1-4\beta \biggl(\frac{x}{\beta}\biggr)\cos\theta + \cancelto{0}{\mathcal{O}(\beta^2)} \biggr] \biggr\} </math>

 

 

<math>~ -~\beta^2 \biggl(\frac{x}{\beta}\biggr) \cos\theta [2^9\cdot 3 (n+1)^3]^{1/2} \cdot (1+\cancelto{0}{x}b)^{-1/2}\cdot \biggl(1+\frac{3\cancelto{0}{x}b}{2}\biggr)^2 </math>

 

 

<math>~+\beta^2 \biggl(\frac{x}{\beta}\biggr)\cdot 2 m^2 [1 - 2 + 1 ] \cdot [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta \cdot (1+\cancelto{0}{x}b)^{1/2} </math>

 

<math>~\approx</math>

<math>~ -~\beta^1 [ 2^7\cdot 3 (n+1)^3]^{1/2} [1 - 1] </math>

 

 

<math>~ -~\beta^2 \biggl(\frac{x}{\beta}\biggr)\cos\theta [ 2^9\cdot 3 (n+1)^3]^{1/2} -~\beta^2 \biggl(\frac{x}{\beta}\biggr) \cos\theta [2^9\cdot 3 (n+1)^3]^{1/2} </math>

 

 

<math>~+\beta^2 \biggl(\frac{x}{\beta}\biggr)\cdot 2 m^2 [1 - 2 + 1 ] \cdot [ 2^7\cdot 3(n+1)^3 ]^{1/2} \cos\theta </math>

See Also


Whitworth's (1981) Isothermal Free-Energy Surface

© 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