User:Tohline/StabilityVariationalPrincipal
From VisTrailsWiki
Contents 
FreeEnergy Stability Analysis
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
Most General Case
Consider a freeenergy function of the form,



where, and are constants, and the dimensionless configuration radius,
is defined in terms of a characteristic length, , which is likely to be different for each type of problem.
Virial Equilibrium
The first variation (first derivative) of this function with respect to the configuration's radius is,



According to the virial theorem, the radius of an equilibrium configuration is obtained by setting and identifying the roots of the resulting equation. For example, identifying roots of the polynomial expression,



Stability
Let's rewrite the first variation of the freeenergy function in terms of three coefficients which, in general, we will permit to have different values from the original three ,



The first variation (first derivative) of this function with respect to the configuration's radius — which, in effect, represents the second variation of the freeenergy function — gives,



If we evaluate this function by setting , the sign of the resulting expression should indicate stability (positive) or dynamical instability (negative); and the marginally unstable configuration is identified by the value of for which .
PressureTruncated Configurations
Expectations
For pressuretruncated polytropes, we set and let represent the chosen polytropic index. In this situation, then, we have,
Freeenergy expression: 




Virial equilibrium: 




Stability indicator: 



Hence, the (critical) equilibrium radius of the marginally unstable configuration is given by the expression,















Notice that, if , this gives,






Energies and Structural Form Factors
Old Approach
As has been developed in, for example, our accompanying review, we adopt the following normalizations:












Then, from separate summaries — both here and here — we can write,





















where the structural form factors are defined as follows:









This gives, specifically for specifically for pressuretruncated polytropic configurations,









New Approach
In order to accommodate the structural integrals required by the Ledoux variational principle, let's rederive some of these key energy and formfactor expressions. Basically, we will be repeating some earlier derivations.
Mass
Defining as the total mass of the isolated configuration, while is the truncated configuration's mass; defining as the truncated configuration's (not necessarily equilibrium) radius; and being careful to define the mean density of the truncated configuration such that,
we have,


















Acknowledging that when the upper integration limit goes to , we see that the "mass" formfactor is,



Now, from the,
we realize that,



So these last two expressions may also be written as,



and,



Modified Internal Energy
Now we want to develop the appropriately scaled integral definition of a "variational" internal energy having the form,









Hence, the coefficient, , in the freeenergy expression is,



or, if , then,



where,



When , then according to Viala & Horedt (1974), this integral over polytropic functions becomes,






which matches the expression for derived earlier.
Modified Gravitational Potential Energy
Similarly, we have,















Hence, the coefficient, , in the freeenergy expression is,



or, if , then,



where,



Now, according to Viala & Horedt (1974), when , this integral over polytropic functions becomes,






As we have detailed elsewhere, from this, we have deduced that, for polytropic configurations,






Test Virial Equilibrium Condition
If the correct normalized equilibrium radius, , is specified, our expectation regarding virial equilibrium is that,



Let's see if this expression is valid when we plug in the expressions for the equilibrium parameter pair — and — that has been given by Horedt (1970), namely,






where we have taken into account the shift in normalization factors,
Switch from Hoerdt's (1970) Normalization  







We therefore have:
First Term












Second Term





















Third Term






Combined
Combining the three terms, the virial expression becomes,















Q. E. D.
The Ledoux Variational Principle
Drawing from a separate presentation of Ledoux's variational principle, let's normalize his Lagrangian using the same normalizations that have been used, above. His expression is …



Our normalization produces,


















Given that,
and,
this expression for the Lagrangian becomes,






In an effort to help identify the various terms in this expression as well as the relationship between the entire expression and our unperturbed free energy expression, let's ignore all terms involving the radial eigenfunction, , and its derivative. In this case, we have,









See Also
© 2014  2019 by Joel E. Tohline 