An Introduction to Modern Astrophysics · Problem 21.2

The Milky Way Galaxy — Carroll & Ostlie

Physics — Astrophysics · polytropic-jupiter-interior-model

Problem

Analytic functions can be derived for the pressure and density structure in the interior of Jupiter if an approximate relationship between pressure and density is assumed. A reasonable choice for a composition of pure molecular hydrogen is P(r)=Kρ2(r),P(r) = K\rho^2(r), where KK is a constant. This type of analytic model is known as a polytrope (see page 334ff for a discussion of polytropic stellar models). (a) By substituting the expression for the pressure into the hydrostatic equilibrium equation (Eq. 10.6) and differentiating, show that a second-order differential equation for the density can be obtained, namely d2ρdr2+2rdρdr+(2πGK)ρ=0.\frac{d^2\rho}{dr^2} + \frac{2}{r}\frac{d\rho}{dr} + \left(\frac{2\pi G}{K}\right)\rho = 0. (b) Show that the equation is satisfied by ρ(r)=ρc(sin⁡krkr),\rho(r) = \rho_c \left(\frac{\sin kr}{kr}\right), where ρc\rho_c is the density at the center of the planet and k≡(2πGK)1/2.k \equiv \left(\frac{2\pi G}{K}\right)^{1/2}. (c) Taking the average radius of Jupiter to be RJ=6.99×107R_J = 6.99 \times 10^7 m and assuming that the density goes to zero at the surface (i.e., kRJ=πkR_J = \pi), determine the values of kk and KK. (d) Integrate Eq. (10.7) using the analytical solution for Jupiter's density as a function of radius to find an expression for the planet's interior mass, MrM_r, written in terms of rr and ρc\rho_c. Hint: ∫r(sin⁡kr) dr=1k2sin⁡kr−rkcos⁡kr.\int r(\sin kr)\,dr = \frac{1}{k^2}\sin kr - \frac{r}{k}\cos kr. (e) Using the boundary condition that Mr=MJM_r = M_J at the surface, estimate the planet's central density. (The value obtained in this problem is lower than the result found from detailed numerical calculations, 1500 kg m−3^{-3}. One major reason for the difference is that Jupiter's composition is not purely molecular hydrogen.) (f) Make separate plots of the density and interior mass as functions of radius. (g) What is the central pressure of your model of Jupiter? (One detailed model gives a value of 8×10128 \times 10^{12} N m−2^{-2}.)

Browse the problem library