Chandrasekhar virial equations

Consider a fluid mass ${\displaystyle M}$ of volume ${\displaystyle V}$ with density ${\displaystyle \rho (\mathbf {x} ,t)}$ and an isotropic pressure ${\displaystyle p(\mathbf {x} ,t)}$ with vanishing pressure at the bounding surfaces. Here, ${\displaystyle \mathbf {x} }$ refers to a frame of reference attached to the center of mass. Before describing the virial equations, let's define some moments.

The density moments are defined as

${\displaystyle M=\int _{V}\rho \,d\mathbf {x} ,\quad I_{i}=\int _{V}\rho x_{i}\,d\mathbf {x} ,\quad I_{ij}=\int _{V}\rho x_{i}x_{j}\,d\mathbf {x} ,\quad I_{ijk}=\int _{V}\rho x_{i}x_{j}x_{k}\,d\mathbf {x} ,\quad I_{ijk\ell }=\int _{V}\rho x_{i}x_{j}x_{k}x_{\ell }\,d\mathbf {x} ,\quad {\text{etc.}}}$

the pressure moments are

${\displaystyle \Pi =\int _{V}p\,d\mathbf {x} ,\quad \Pi _{i}=\int _{V}px_{i}\,d\mathbf {x} ,\quad \Pi _{ij}=\int _{V}px_{i}x_{j}\,d\mathbf {x} ,\quad \Pi _{ijk}=\int _{V}px_{i}x_{j}x_{k}d\mathbf {x} \quad {\text{etc.}}}$

the kinetic energy moments are

${\displaystyle T_{ij}={\frac {1}{2}}\int _{V}\rho u_{i}u_{j}\,d\mathbf {x} ,\quad T_{ij;k}={\frac {1}{2}}\int _{V}\rho u_{i}u_{j}x_{k}\,d\mathbf {x} ,\quad T_{ij;k\ell }={\frac {1}{2}}\int _{V}\rho u_{i}u_{j}x_{k}x_{\ell }\,d\mathbf {x} ,\quad \mathrm {etc.} }$

and the Chandrasekhar potential energy tensor moments are

${\displaystyle W_{ij}=-{\frac {1}{2}}\int _{V}\rho \Phi _{ij}\,d\mathbf {x} ,\quad W_{ij;k}=-{\frac {1}{2}}\int _{V}\rho \Phi _{ij}x_{k}\,d\mathbf {x} ,\quad W_{ij;k\ell }=-{\frac {1}{2}}\int _{V}\rho \Phi _{ij}x_{k}x_{\ell }d\mathbf {x} ,\quad \mathrm {etc.} \quad {\text{where}}\quad \Phi _{ij}=G\int _{V}\rho (\mathbf {x'} ){\frac {(x_{i}-x_{i}')(x_{j}-x_{j}')}{|\mathbf {x} -\mathbf {x'} |^{3}}}\,d\mathbf {x'} }$

where ${\displaystyle G}$ is the gravitational constant.

All the tensors are symmetric by definition. The moment of inertia ${\displaystyle I}$, kinetic energy ${\displaystyle T}$ and the potential energy ${\displaystyle W}$ are just traces of the following tensors

${\displaystyle I=I_{ii}=\int _{V}\rho |\mathbf {x} |^{2}\,d\mathbf {x} ,\quad T=T_{ii}={\frac {1}{2}}\int _{V}\rho |\mathbf {u} |^{2}\,d\mathbf {x} ,\quad W=W_{ii}=-{\frac {1}{2}}\int _{V}\rho \Phi \,d\mathbf {x} \quad {\text{where}}\quad \Phi =\Phi _{ii}=\int _{V}{\frac {\rho (\mathbf {x'} )}{|\mathbf {x} -\mathbf {x'} |}}\,d\mathbf {x'} }$

Chandrasekhar assumed that the fluid mass is subjected to pressure force and its own gravitational force, then the Euler equations is

${\displaystyle \rho {\frac {du_{i}}{dt}}=-{\frac {\partial p}{\partial x_{i}}}+\rho {\frac {\partial \Phi }{\partial x_{i}}},\quad {\text{where}}\quad {\frac {d}{dt}}={\frac {\partial }{\partial t}}+u_{j}{\frac {\partial }{\partial x_{j}}}}$

The Euler equations in a rotating frame of reference, rotating with an angular velocity ${\displaystyle \mathbf {\Omega } }$ is given by

${\displaystyle \rho {\frac {du_{i}}{dt}}=-{\frac {\partial p}{\partial x_{i}}}+\rho {\frac {\partial \Phi }{\partial x_{i}}}+{\frac {1}{2}}\rho {\frac {\partial }{\partial x_{i}}}|\mathbf {\Omega } \times \mathbf {x} |^{2}+2\rho \varepsilon _{i\ell m}u_{\ell }\Omega _{m}}$

where ${\displaystyle \varepsilon _{i\ell m}}$ is the Levi-Civita symbol, ${\displaystyle {\frac {1}{2}}|\mathbf {\Omega } \times \mathbf {x} |^{2}}$ is the centrifugal acceleration and ${\displaystyle 2\mathbf {u} \times \mathbf {\Omega } }$ is the Coriolis acceleration.

Consider the Navier-Stokes equations instead of Euler equations,

${\displaystyle \rho {\frac {du_{i}}{dt}}=-{\frac {\partial p}{\partial x_{i}}}+\rho {\frac {\partial \Phi }{\partial x_{i}}}+{\frac {\partial \tau _{ik}}{\partial x_{k}}},\quad {\text{where}}\quad \tau _{ik}=\rho \nu \left({\frac {\partial u_{i}}{\partial x_{k}}}+{\frac {\partial u_{k}}{\partial x_{i}}}-{\frac {2}{3}}{\frac {\partial u_{l}}{\partial x_{l}}}\delta _{ik}\right)}$

and we define the shear-energy tensor as

${\displaystyle S_{ij}=\int _{V}\tau _{ij}d\mathbf {x} .}$

With the condition that the normal component of the total stress on the free surface must vanish, i.e., ${\displaystyle (-p\delta _{ik}+\tau _{ik})n_{k}=0}$, where ${\displaystyle \mathbf {n} }$ is the outward unit normal, the second order virial equation then be

${\displaystyle {\frac {1}{2}}{\frac {d^{2}I_{ij}}{dt^{2}}}=2T_{ij}+W_{ij}+\delta _{ij}\Pi -S_{ij}.}$

This can be easily extended to rotating frame of references.

• Virial theorem
• Dirichlet's ellipsoidal problem
• Chandrasekhar tensor