How Do Radial Oscillations Affect Star Stability?

Click For Summary
Radial oscillations play a critical role in the stability of stars, as indicated by the equations derived in the discussion. The equations focus on the relationship between pressure, density, and radial displacement, suggesting that oscillations can be modeled using linearized forms around an equilibrium state. The variable ##\xi_r## represents small radial displacements, and the presence of an oscillation frequency term indicates a dependence on time and radial position. The discussion highlights the complexity of deriving the governing equations, particularly in maintaining linearity while eliminating mixed derivatives. Overall, understanding these dynamics is essential for analyzing star stability under oscillatory conditions.
ergospherical
Science Advisor
Homework Helper
Education Advisor
Insights Author
Messages
1,100
Reaction score
1,387
Homework Statement
The spherically symmetric oscillations of a star satisfy the Sturm-Liouville equation\begin{align*}
\frac{d}{dr} \left[ \frac{\gamma p}{r^2} \frac{d}{dr}(r^2 \xi_r) \right] - \frac{4}{r} \frac{dp}{dr} \xi_r + \rho \omega^2 \xi_r = 0
\end{align*}
Relevant Equations
N/A
So far I have not made much meaningful progress beyond two equations; \begin{align*}
\rho \frac{D\mathbf{u}}{Dt} = - \nabla p \implies \rho \left( \frac{\partial}{\partial t} + u \frac{\partial}{\partial r} \right)u = - \frac{\partial p}{\partial r}
\end{align*}and thermal energy:\begin{align*}
\frac{Dp}{Dt} = -\gamma p \nabla \cdot \mathbf{u} \implies \left( \frac{\partial}{\partial t} + u \frac{\partial}{\partial r} \right)p = -\gamma p \frac{\partial}{\partial r}(r^2 u)
\end{align*}I'm somewhat confused why there's no time dependence in the form given. What exactly is ##\xi_r##?
 
Physics news on Phys.org
You are dealing with oscillations, and there is an \omega^2 in the ODE. That suggests variables are assumed to be functions of r alone times e^{i\omega t}. I suspect \xi_r denotes (small) radial displacement from an equilibrium position, in which case <br /> u_r = \frac{\partial}{\partial t}(\xi_r(r)e^{i\omega t}) = i\omega \xi_r(r) e^{i\omega t}. The ODE you are trying to derive is linear and the governing PDEs are not, which suggests that the system has been linearised about the static equiblibrium state.
 
  • Like
Likes ergospherical
Thank you. I applied ##\partial/\partial t## to my first equation and ##\partial/\partial r## to my second equation, and subsequently eliminated ##\partial^2 p / \partial r \partial t##. Then I substituted for ##u = i\omega \xi_r e^{i \omega t}## and kept only terms to linear order in ##\xi_r##. After some manipulation,
\begin{align*}
(\xi_r p')' + (\gamma p(r^2 \xi_r)')' + \rho \omega^2 \xi_r = 0
\end{align*}Seems likely that I made at least one slip. I will check again tomorrow with fresher eyes!
 
ergospherical said:
I will check again tomorrow with fresher eyes!
Have you cleared that with the fresher? :eek:
 
pbuk said:
Have you cleared that with the fresher? :eek:
That would imply that they have a choice!
 
At first, I derived that: $$\nabla \frac 1{\mu}=-\frac 1{{\mu}^3}\left((1-\beta^2)+\frac{\dot{\vec\beta}\cdot\vec R}c\right)\vec R$$ (dot means differentiation with respect to ##t'##). I assume this result is true because it gives valid result for magnetic field. To find electric field one should also derive partial derivative of ##\vec A## with respect to ##t##. I've used chain rule, substituted ##\vec A## and used derivative of product formula. $$\frac {\partial \vec A}{\partial t}=\frac...

Similar threads