User:Tohline/SSC/Stability/MoreGeneralApproach

From VistrailsWiki
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

More General Approach to the Parabolic Eigenvalue Problem

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

The material presented in this chapter is an extension of the chapter titled, Other Analytic Models and could also be considered to be a subsection of the associated chapter titled, Other Analytic Ramblings. More specifically, in the following "Introduction," we repeat a manipulation of the LAWE that was originally developed in the subsection of that chapter titled, "Consider Parabolic Case".



Work-in-progress.png

Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
|   Go Home   |


Introduction

In the case of a parabolic density distribution, the LAWE may be written in the form,

<math>~\frac{2}{(1-x^2)(2-x^2)} \biggl[ \biggl( \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr)(5-3x^2) -\sigma^2 \biggr]</math>

<math>~=</math>

<math>~ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} </math>


Let's try,

<math>~\mathcal{G}_\sigma</math>

<math>~=</math>

<math>~(a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^m \, ,</math>

which implies,

<math>~\mathcal{G}_\sigma^'</math>

<math>~=</math>

<math>~n(a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^m +m (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-1}(2b_2x)</math>

<math>~\Rightarrow ~~~~ \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~n(a_0 + a_2x^2)^{-1}(2a_2x^2) +m (b_0 + b_2x^2)^{-1}(2b_2x^2) </math>

 

<math>~=</math>

<math>~\frac{2x^2}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ n a_2 (b_0 + b_2x^2) +mb_2 (a_0 + a_2x^2) \biggr] </math>

 

<math>~=</math>

<math>~\frac{2x^2}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ (n a_2 b_0 + mb_2 a_0) +(na_2 b_2+ mb_2 a_2)x^2\biggr] \, ,</math>

and,

<math>~\mathcal{G}_\sigma^{' '}</math>

<math>~=</math>

<math>~n m (a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^{m-1}(2b_2x) + n(a_0 + a_2x^2)^{n-1}(2a_2) \cdot (b_0 + b_2x^2)^m + n(n-1)(a_0 + a_2x^2)^{n-2}(2a_2x)^2 \cdot (b_0 + b_2x^2)^m </math>

 

 

<math>~+m n(a_0 + a_2x^2)^{n-1}(2a_2x) \cdot (b_0 + b_2x^2)^{m-1}(2b_2x) +m (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-1}(2b_2) +m(m-1) (a_0 + a_2x^2)^n \cdot (b_0 + b_2x^2)^{m-2}(2b_2x)^2</math>

<math>~\Rightarrow ~~~~ \frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~8n m a_2b_2 x^2 (a_0 + a_2x^2)^{-1}\cdot (b_0 + b_2x^2)^{-1} + n2a_2 (a_0 + a_2x^2)^{-1} + n(n-1)4a_2^2 x^2 (a_0 + a_2x^2)^{-2} +m2b_2 (b_0 + b_2x^2)^{-1} +m(m-1)4 b_2^2 x^2 (b_0 + b_2x^2)^{-2}</math>

 

<math>~=</math>

<math>~\frac{2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)} + \biggl[ \frac{4n(n-1) a_2^2 }{ (a_0 + a_2x^2)^{2}} + \frac{8n m a_2b_2}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}+ \frac{4m(m-1) b_2^2 }{(b_0 + b_2x^2)^{2}} \biggr]x^2 </math>

So, we have for the LAWE:

LHS

<math>~=</math>

<math>~ \frac{2}{(1-x^2)(2-x^2)} \biggl[ \biggl( \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr)(5-3x^2) -\sigma^2 \biggr]</math>

 

<math>~=</math>

<math>~ \frac{2}{(1-x^2)(2-x^2)(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl\{ \biggl[ \alpha(a_0 + a_2x^2) (b_0 + b_2x^2) + 2x^2(n a_2 b_0 + mb_2 a_0) + 2x^4 (na_2 b_2+ mb_2 a_2) \biggr](5-3x^2) </math>

 

 

<math>~ -\sigma^2 (a_0 + a_2x^2) (b_0 + b_2x^2) \biggr\} \, ;</math>

RHS

<math>~=</math>

<math>~ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} </math>

 

<math>~=</math>

<math>~\frac{2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)} + \biggl[ \frac{4n(n-1) a_2^2 }{ (a_0 + a_2x^2)^{2}} + \frac{8n m a_2b_2}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}+ \frac{4m(m-1) b_2^2 }{(b_0 + b_2x^2)^{2}} \biggr]x^2 </math>

 

 

<math>~ + \frac{8}{(a_0 + a_2x^2) (b_0 + b_2x^2)} \biggl[ (n a_2 b_0 + mb_2 a_0) +(na_2 b_2+ mb_2 a_2)x^2\biggr] </math>

 

<math>~=</math>

<math>~\frac{1}{(a_0 + a_2x^2)(b_0 + b_2x^2)} \biggl\{ 2n a_2(b_0 + b_2x^2) + 2m b_2 (a_0 + a_2x^2) + 8(n a_2 b_0 + mb_2 a_0) + 8(na_2 b_2+ mb_2 a_2)x^2 </math>

 

 

<math>~ + \biggl[8n m a_2b_2+ \frac{4n(n-1) a_2^2(b_0 + b_2x^2) }{ (a_0 + a_2x^2)} + \frac{4m(m-1) b_2^2(a_0 + a_2x^2) }{(b_0 + b_2x^2)} \biggr]x^2 \biggr\} \, . </math>

Putting these together gives,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(a_0 + a_2x^2) (b_0 + b_2x^2) + 2x^2(n a_2 b_0 + mb_2 a_0) + 2x^4 (na_2 b_2+ mb_2 a_2) \biggr](5-3x^2) -\sigma^2 (a_0 + a_2x^2) (b_0 + b_2x^2) </math>

 

 

<math>~ - \biggl[ n a_2(b_0 + b_2x^2) + m b_2 (a_0 + a_2x^2) + 4(n a_2 b_0 + mb_2 a_0) + 4(na_2 b_2+ mb_2 a_2)x^2+ 4n m a_2b_2x^2 \biggr](1-x^2)(2-x^2) </math>

 

 

<math>~ - \frac{(1-x^2)(2-x^2)}{ (a_0 + a_2x^2)(b_0 + b_2x^2)}\biggl[2n(n-1) a_2^2(b_0 + b_2x^2)^2 + 2m(m-1) b_2^2(a_0 + a_2x^2)^2 \biggr]x^2 \, . </math>

Additional Setup

Benefitting from our earlier exploration of this problem, let's divide through by the product, <math>~(a_0 b_0)</math>, and introduce the new variable notations,

<math>~\lambda \equiv \frac{a_2}{a_0} \, ,</math>       and       <math>~\eta \equiv \frac{b_2}{b_0} \, .</math>

The LAWE becomes,

<math>~ 0 </math>

<math>~=</math>

<math>~ \biggl[ \alpha(1 + \lambda x^2) (1 + \eta x^2) + 2x^2(n \lambda + m\eta ) + 2x^4 (n\lambda \eta + m\eta \lambda ) \biggr](5-3x^2) -\sigma^2 (1 + \lambda x^2) (1 + \eta x^2) </math>

 

 

<math>~ - \biggl[ n \lambda (1 + \eta x^2) + m \eta (1 + \lambda x^2) + 4(n \lambda + m\eta ) + 4(n\lambda \eta + m\eta \lambda )x^2+ 4n m \lambda \eta x^2 \biggr](1-x^2)(2-x^2) </math>

 

 

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

Multiplying through by the denominator of the last term(s) — that is, multiplying through by <math>~(1 + \lambda x^2)(1 + \eta x^2)</math> — will give us a polynomial with coefficient expressions for 6 terms <math>~(x^0, x^2, x^4, x^6, x^8, x^{10})</math> expressed in terms of 5 unknowns <math>~(\sigma^2, n, m, \lambda, \eta)</math>.

Wouldn't a better strategy be to insert yet another quadratic factor — specifically, <math>~(1+\beta x^2)^\ell</math> — which will introduce two additional unknowns but only add one more term into the polynomial expression? This would bring the total number of coefficient expressions to 7 while simultaneously raising the number of unknowns to 7. It will be tedious and messy, but worth the try.

Expanding from Two to Three Quadratic Terms

Here we rearrange terms in the "parabolic" LAWE to construct the governing ODE as,

<math>~\sigma^2 </math>

<math>~=</math>

<math>~5(1-\tfrac{3}{5}x^2) \biggl[ \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr] - (1-x^2)(1-\tfrac{1}{2}x^2)\biggl[ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} \biggr] </math>

Let's try,

<math>~\mathcal{G}_\sigma</math>

<math>~=</math>

<math>~(1 + \lambda x^2)^n \cdot (1 + \eta x^2)^m \cdot (1 + \beta x^2)^\ell \, ,</math>

or, in an effort to permit writing more compact expressions,

<math>~\mathcal{G}_\sigma</math>

<math>~=</math>

<math>~N^n \cdot M^m \cdot L^\ell \, ,</math>

where,

<math>~N \equiv (1 + \lambda x^2)\, ;</math>       <math>~M \equiv (1 + \eta x^2)\, ;</math>       and       <math>~L \equiv (1 + \beta x^2)\, .</math>      

This implies (after some whiteboard derivations),

<math>~\frac{x\mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~ \frac{2x^2}{N\cdot M \cdot L} \biggl[ \ell \beta M\cdot N + m\eta L\cdot N + n\lambda L\cdot M\biggr] \, , </math>

<math>~\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}</math>

<math>~=</math>

<math>~\frac{4 x^2}{N\cdot M \cdot L} \biggl[ \ell \beta (m\eta N + n\lambda M) + m\eta (\ell \beta N + n\lambda L) + n\lambda (\ell \beta M + m\eta L) \biggr] </math>

 

 

<math>~ + \frac{2\ell \beta}{L^2} \biggl[ 1 + x^2 \beta (2\ell -1)\biggr] + \frac{2m\eta}{M^2}\biggl[ 1+x^2 \eta(2m-1)\biggr] + \frac{2n\lambda}{N^2}\biggl[1 + x^2 \lambda(2n-1) \biggr] \, . </math>

Specific Values of Quadratic Coefficients

Now, if we assume that,

<math>~\lambda = -1 \, ;</math>       <math>~\eta = -\tfrac{1}{2} \, ;</math>       and       <math>~\beta = - \tfrac{3}{5} \, .</math>      

the "parabolic" LAWE becomes,


<math>~L \cdot \sigma^2 </math>

<math>~=</math>

<math>~5L^2 \biggl[ \alpha + \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma}\biggr] - N\cdot M\cdot L\biggl[ \biggl(\frac{\mathcal{G}_\sigma^{' '}}{\mathcal{G}_\sigma}\biggr) +\frac{4}{x^2} \cdot \frac{x \mathcal{G}_\sigma^'}{\mathcal{G}_\sigma} \biggr] \, . </math>

Then, plugging in the expressions for <math>~\mathcal{G}_\sigma</math> and its derivatives, we have,

<math>~L \cdot \sigma^2 </math>

<math>~=</math>

<math>~5L^2 \biggl\{ \alpha + \frac{2x^2}{N\cdot M \cdot L} \biggl[ \ell \beta M\cdot N + m\eta L\cdot N + n\lambda L\cdot M\biggr] \biggr\} </math>

 

 

<math>~ - \biggl\{ 4x^2\biggl[ \ell \beta (m\eta N + n\lambda M) + m\eta (\ell \beta N + n\lambda L) + n\lambda (\ell \beta M + m\eta L) \biggr] + 8\biggl[ \ell \beta M\cdot N + m\eta L\cdot N + n\lambda L\cdot M \biggr] \biggr\} </math>

 

 

<math>~ - N\cdot M\cdot L\biggl\{ \frac{2\ell \beta}{L^2} \biggl[ 1 + x^2 \beta (2\ell -1)\biggr] + \frac{2m\eta}{M^2}\biggl[ 1+x^2 \eta(2m-1)\biggr] + \frac{2n\lambda}{N^2}\biggl[1 + x^2 \lambda(2n-1) \biggr] \biggr\} </math>

 

<math>~=</math>

<math>~5L^2 \biggl\{ \alpha - \frac{2x^2}{N\cdot M \cdot L} \biggl[ \tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N + n L\cdot M\biggr] \biggr\} </math>

 

 

<math>~ + \biggl\{ - 4x^2\biggl[ \tfrac{3}{5}\ell (\tfrac{1}{2}m N + n M) + \tfrac{1}{2}m (\tfrac{3}{5}\ell N + n L) + n (\tfrac{3}{5}\ell M + \tfrac{1}{2}m L) \biggr] + 8\biggl[ \tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N + n L\cdot M \biggr] \biggr\} </math>

 

 

<math>~ + \biggl\{ \frac{6\ell N\cdot M}{5 L} \biggl[ 1 - \tfrac{3}{5} (2\ell -1)x^2 \biggr] + \frac{m N\cdot L}{M}\biggl[ 1 - \tfrac{1}{2}(2m-1)x^2 \biggr] + \frac{2nM\cdot L}{N}\biggl[1 - (2n-1) x^2\biggr] \biggr\} </math>

<math>~\Rightarrow ~~~ N\cdot M\cdot L^2 \sigma^2 </math>

<math>~=</math>

<math>~5L^2 \biggl\{ N\cdot M\cdot L\cdot \alpha - 2x^2 [ \tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N + n L\cdot M ] \biggr\} </math>

 

 

<math>~ + N\cdot M\cdot L\biggl\{ - 4x^2 [ \tfrac{3}{5}\ell (\tfrac{1}{2}m N + n M) + \tfrac{1}{2}m (\tfrac{3}{5}\ell N + n L) + n (\tfrac{3}{5}\ell M + \tfrac{1}{2}m L) ] + 8 [ \tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N + n L\cdot M ] \biggr\} </math>

 

 

<math>~ + \biggl\{ \tfrac{6}{5} \ell N^2\cdot M^2 [ 1 - \tfrac{3}{5} (2\ell -1)x^2 ] + m N^2\cdot L^2 [ 1 - \tfrac{1}{2}(2m-1)x^2 ] + 2nM^2\cdot L^2 [1 - (2n-1) x^2] \biggr\} </math>

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

<math>~=</math>

<math>~5L \biggl\{ \alpha - 2x^2 \biggl[ \frac{3}{5}\ell \biggl( \frac{1}{L}\biggr) + \frac{1}{2}m \biggl(\frac{1}{M}\biggr) + n \biggl(\frac{1}{N}\biggr) \biggr] \biggr\} </math>

 

 

<math>~ + \frac{1}{L}\biggl\{ - 4x^2 [ \tfrac{3}{5}\ell (\tfrac{1}{2}m N + n M) + \tfrac{1}{2}m (\tfrac{3}{5}\ell N + n L) + n (\tfrac{3}{5}\ell M + \tfrac{1}{2}m L) ] + 8 [ \tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N + n L\cdot M ] \biggr\} </math>

 

 

<math>~ + \frac{1}{L}\biggl\{ \frac{6\ell}{5} \biggl[ \frac{N\cdot M}{L}\biggl] \biggl[ 1 - \frac{3}{5} (2\ell -1)x^2 \biggl] + m \biggl[ \frac{N \cdot L}{M} \biggl]\biggl[ 1 - \tfrac{1}{2}(2m-1)x^2 \biggl] + 2n\biggl[ \frac{M \cdot L}{N} \biggl] \biggl[ 1 - (2n-1) x^2\biggl] \biggr\} </math>

 

<math>~=</math>

<math>~5L \biggl\{ \alpha + 2 \biggl[ (L-1)\ell \biggl( \frac{1}{L}\biggr) + (M-1)m \biggl(\frac{1}{M}\biggr) + (N-1)n \biggl(\frac{1}{N}\biggr) \biggr] \biggr\} </math>

 

 

<math>~ + \frac{4}{L}\biggl\{ \tfrac{6}{5}\ell M\cdot N + m L\cdot N +2 n L\cdot M + (L-1)\ell (\tfrac{1}{2}m N + n M) + (M-1) m (\tfrac{3}{5}\ell N + n L) + (N-1)n (\tfrac{3}{5}\ell M + \tfrac{1}{2}m L) \biggr\} </math>

 

 

<math>~ + \frac{1}{L}\biggl\{ \frac{6\ell}{5} \biggl[ \frac{N\cdot M}{L}\biggl] \biggl[ 1 + (L-1)(2\ell -1) \biggl] + m \biggl[ \frac{N \cdot L}{M} \biggl]\biggl[ 1 +(M-1)(2m-1) \biggl] + 2n\biggl[ \frac{M \cdot L}{N} \biggl] \biggl[ 1+ (N-1) (2n-1) \biggl] \biggr\} </math>

 

<math>~=</math>

<math>~5L \biggl\{ \alpha + 2 \biggl[ \biggl(1-\frac{1}{L} \biggr)\ell + \biggl(1-\frac{1}{M} \biggr)m + \biggl(1-\frac{1}{N} \biggr)n \biggr] \biggr\} </math>

 

 

<math>~ + \frac{4}{L}\biggl\{ \tfrac{6}{5}\ell M\cdot N + m L\cdot N +2 n L\cdot M + \ell (\tfrac{1}{2}m N\cdot L + n M\cdot L) -\ell (\tfrac{1}{2}m N + n M) + m (\tfrac{3}{5}\ell N\cdot M + n L\cdot M) - m (\tfrac{3}{5}\ell N + n L) + n (\tfrac{3}{5}\ell M\cdot N + \tfrac{1}{2}m L\cdot N) - n (\tfrac{3}{5}\ell M + \tfrac{1}{2}m L) \biggr\} </math>

 

 

<math>~ + \frac{1}{L}\biggl\{ \frac{6\ell}{5} \biggl[ \frac{N\cdot M}{L}\biggl] \biggl[ 2(1-\ell) + L(2\ell -1) \biggl] + m \biggl[ \frac{N \cdot L}{M} \biggl]\biggl[ 2(1-m) +M(2m-1) \biggl] + 2n\biggl[ \frac{M \cdot L}{N} \biggl] \biggl[ 2(1-n) + N (2n-1) \biggl] \biggr\} </math>

 

<math>~=</math>

<math>~5L \biggl\{ \alpha + 2(\ell + m + n) - 2 \biggl[ \biggl(\frac{\ell}{L} \biggr) + \biggl(\frac{m}{M} \biggr) + \biggl(\frac{n}{N} \biggr) \biggr] \biggr\} </math>

 

 

<math>~ + \frac{2}{L}\biggl\{ \frac{6\ell}{5} N\cdot M \biggl[ 2 + m + n \biggr] + m L\cdot N \biggl[ 2 + \ell + n \biggr] + 2nL \cdot M \biggl[ 2 + \ell + m \biggr] \biggr\} \biggr\} + \frac{1}{L}\biggl\{ \frac{6\ell}{5} N\cdot M \biggl[ (2\ell -1) \biggl] + m L \cdot N \biggl[ (2m-1) \biggl] + 2n L \cdot M \biggl[ (2n-1) \biggl] \biggr\} </math>

 

 

<math>~ + \frac{1}{L}\biggl\{ \frac{6\ell}{5} N\cdot M \biggl[ \frac{2}{L}(1-\ell) \biggl] + m L \cdot N \biggl[ \frac{2}{M} (1-m) \biggl] + 2n L \cdot M \biggl[ \frac{2}{N}(1-n) \biggl] \biggr\} - \frac{4}{L}\biggl\{ \biggl[ \ell m N (\tfrac{1}{2} + \tfrac{3}{5}) + \ell n M(1 + \tfrac{3}{5} ) + m n L (1 + \tfrac{1}{2})\biggr] </math>

 

<math>~=</math>

<math>~5L \biggl\{ \alpha + 2(\ell + m + n) - 2 \biggl[ \biggl(\frac{\ell}{L} \biggr) + \biggl(\frac{m}{M} \biggr) + \biggl(\frac{n}{N} \biggr) \biggr] \biggr\} + \frac{(3 + 2m + 2n + 2\ell)}{L}\biggl[ \frac{6\ell}{5} N\cdot M + m L \cdot N + 2n L \cdot M \biggr] </math>

 

 

<math>~ - \frac{4}{L}\biggl\{ \frac{3}{5} \biggl[ \frac{N\cdot M}{L} \biggr] \ell(\ell - 1) + \frac{1}{2}\biggl[ \frac{L \cdot N}{M}\biggr] m(m-1) + \biggl[ \frac{L \cdot M }{N} \biggr] n(n-1) \biggr\} - \frac{4}{L}\biggl\{ \biggl[ \ell m N (\tfrac{1}{2} + \tfrac{3}{5}) + \ell n M(1 + \tfrac{3}{5} ) + m n L (1 + \tfrac{1}{2})\biggr] \biggr\} </math>


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