恒星结构和演化-学习记录1-第二章-基本方程

第一章-恒星的实测与经验关系

  （本章大致讲了恒星各种参数的测量方式，包括光度测量，质量测量，演化过程统计等，因为后面的部分太难了，所以这里就不整理了）

第二章：基本方程

流体近似

  在研究天体的过程中，由于研究的尺度远大于粒子之间的平均自由程，因此热平衡是仅限于局部的，对于不同的部位，温度都是不一样的。可以用局部的平均值来描写运动，这样的研究手法属于（辐射）流体近似

流体方程-质量守恒

  流体在运动的过程中质量保持不变（质量流入与流出相平衡）：
$${\iiint }_V\frac{\partial \rho }{\partial t}dV+\iint \rho \vec{v}\cdot d\vec{S}=0$$
  写为相应的微分形式：
$$\frac{\partial \rho }{\partial t}+\nabla \cdot (\rho \vec{v})=0$$
  在随流体元运动系中，可以写为：
$$\frac{1}{\rho }\frac{d\rho }{dt}=-\nabla \cdot \vec{v}\text{\hspace{1em}(1)}$$

流体方程-动量守恒

  考虑位于$(x,y,z)$处的一个立方体积元，其边长分别为$\Delta x,\Delta y,\Delta z$，沿着$x,y,z$轴方向，现在尝试对此体积元分析冲量与动量变化率，从$x$方向入手：

引力带来的冲量

  首先考虑引力。体积元在引力势$\Phi$中的引力加速度为$g=-\nabla \Phi$，因此引力$x$方向分量为$-\frac{\partial \Phi }{\partial x}\Delta m=-\frac{\partial \Phi }{\partial x}\rho (\vec{x})\Delta V$

压力带来的冲量

  随后考虑作用在$x$方向两个表面上的压力（以正轴方向为正）。在$x$处的压力为$P(x,y,z)\Delta y\Delta z$，在$x+dx$处的压力为$-P(x+dx,y,z)\Delta y\Delta z$，因此，$x$方向的两个表面差为：
$$P(x,y,z)\Delta y\Delta z-P(x+dx,y,z)\Delta y\Delta z=-\frac{\partial P}{\partial x}\Delta V$$
  从左边到右边，使用了偏微分的定义

流入流出的动量

  考虑完冲量的作用后，再考虑动量的变化，首先是通过各个表面流入流出的动量速率。从$x$方向上来看，运动的过程中，会有动量从$y$方向的两个表面中流出。在这里，以流入为正，流出为负，则在$y+dy$表面处沿$y$方向流出的动量速率为：
$$-\rho (x,y+\Delta y,z)v_y(x,y+\Delta y,z)\Delta x\Delta z{v}_x(x,y+\Delta y,z)$$
  $y$表面处沿$-y$方向流出的动量速率为：
$$\rho (x,y,z)v_y(x,y,z)\Delta x\Delta z{v}_x(x,y,z)$$
  这两个表面相加，即可得到通过$y$方向的动量流出速率：
$$-\frac{\partial }{\partial y}(\rho v_y{v}_x)\Delta V$$
  再考虑$x$和$z$，即可得到总的动量流出速率：
$$-\nabla \cdot (\rho \vec{v}{v}_x)\Delta V$$

自身的动量随时间变化

  还有一项是自身的动量变化：
$$\frac{\partial }{\partial t}(\rho v_x)\Delta V$$

利用动量守恒进行推导

  根据动量守恒，冲量等于动量的变化率，因此结合上述四式，有下述等式：
$$\begin{array}{rl}\frac{\partial }{\partial t}(\rho v_x)\Delta V+\nabla \cdot (\rho \vec{v}{v}_x)\Delta V& =-\frac{\partial P}{\partial x}\Delta V-\frac{\partial \Phi }{\partial x}\rho (\vec{x})\Delta V\\ \frac{\partial }{\partial t}(\rho v_x)+\nabla \cdot (\rho \vec{v}{v}_x)& =-\frac{\partial P}{\partial x}-\frac{\partial \Phi }{\partial x}\rho (\vec{x})\end{array}$$

  结合流体随动坐标系微分的等式$\frac{d}{dt}=\frac{\partial }{\partial t}+\vec{v}\cdot \nabla$，动量守恒等式简化为：
$$\rho \frac{d{v}_x}{dt}=-\frac{\partial P}{\partial x}-\frac{\partial \Phi }{\partial x}\rho (\vec{x})$$
  最后考虑$y$和$z$这两个方向，即可得到矢量式：
$$\rho \frac{d\vec{v}}{dt}=-\nabla P-\rho \nabla P$$

流体方程-能量守恒

从热力学第一定律出发

  热力学第一定律表示为：
$$dQ=dϵ+PdV$$

做功项

  对于单位质量物质而言，其内能为$ϵ$，体积为$V=\frac{1}{\rho }$，因此，做功项可以写为：
$$PdV=-\left(\frac{P}{{\rho }^2}\right)d\rho =\frac{P}{\rho }\nabla \cdot \vec{v}dt$$
  其中，第二个等号是代入了随流体元运动系中的质量守恒方程微分形式（式（1））

热能项

  热能包括两项，第一项表示产能或损失，第二项表示表面流入或流出的能量：
$$dQ=qdt-\frac{1}{\rho }\nabla \cdot \vec{F}dt=(\rho qdt-\nabla \cdot \vec{F}dt)/\rho$$
  其中，$q$为单位质量的产能率，$\vec{F}$为能流密度（我猜测这个能流密度是per体积的，因此才会乘一个$V=\frac{1}{\rho }$）

能量方程

  将上述的做功项和热能项代入到热力学第一定律中，即可得到能量方程：
$$\rho \frac{dϵ}{dt}=-P\nabla \cdot \vec{v}-\nabla \cdot \vec{F}+\rho q$$

能量方程的另一形式

  热力学第一定律还可以变换为如下的形式：
$$dQ=dϵ+PdV={\left(\frac{\partial ϵ}{\partial V}\right)}_{V}dT+\left[{\left(\frac{\partial ϵ}{\partial V}\right)}_T+P\right]dV$$
  利用各种热力学关系以及麦克斯韦关系式，有：
$$\left[{\left(\frac{\partial ϵ}{\partial V}\right)}_T+P\right]dV=-\frac{T}{\rho }{\left(\frac{\partial P}{\partial T}\right)}_V\frac{d\rho }{\rho }=-\frac{P}{\rho }\frac{\delta }{\alpha }\frac{d\rho }{\rho }$$
  其中：
$$\begin{array}{rl}\alpha & \equiv {\left(\frac{\partial \ln \rho }{\partial \ln P}\right)}_T=-\frac{P}{V}{\left(\frac{\partial V}{\partial P}\right)}_T\\ \delta & \equiv -{\left(\frac{\partial \ln \rho }{\partial \ln T}\right)}_P=\frac{T}{V}{\left(\frac{\partial V}{\partial T}\right)}_P\end{array}$$
  由于：
$$\begin{array}{rl}{c}_V& \equiv {\left(\frac{dQ}{dT}\right)}_V={\left(\frac{\partial ϵ}{\partial T}\right)}_V\\ c_P& \equiv {\left(\frac{dQ}{dT}\right)}_P={\left(\frac{\partial ϵ}{\partial T}\right)}_P+P{\left(\frac{\partial V}{\partial T}\right)}_P\end{array}$$
  因此：
$$\begin{array}{rl}dQ& ={\left(\frac{\partial Q}{\partial T}\right)}_{P}dT+{\left(\frac{\partial Q}{\partial P}\right)}_{T}dP\\ & =c_{P}dT+T{\left(\frac{\partial S}{\partial P}\right)}_{T}dP\\ & =c_{P}dT-T{\left(\frac{\partial V}{\partial T}\right)}_{P}dP\\ & =c_{P}dT-\frac{\delta }{\rho }dP\end{array}$$
  与前面那种形式相对比：
$$c_{P}dT-\frac{\delta }{\rho }dP=(\rho qdt-\nabla \cdot \vec{F}dt)/\rho$$
  定义：
$$-q_g\equiv c_P\frac{\partial T}{\partial t}-\frac{\delta }{\rho }\frac{\partial P}{\partial t}=q-\frac{\nabla \cdot \vec{F}}{\rho }$$

能流

  能流：单位时间内，通过单位面积，从一边到另一边的能量

考虑介于$z-\lambda /2$和$z+\lambda /2$之间的粒子，设粒子的单位体积密度为$n_i$，运动速度为$v_T$，则单位时间内穿越这两个界面的粒子数为$\frac{1}{3}{n}_i{v}_T\Delta x\Delta y$

  从下层和上层进入的粒子具有能量差$-e_i(z+\lambda /2)+e_i(z-\lambda /2)$（$e_i(z)$为界面$z$处粒子的平均能量），将此能量差和穿过界面的粒子数相乘，并写为矢量形式，即可得到能流的表达式：
$${\vec{F}}_i=-\frac{1}{3}{v}_T\lambda \nabla u_i$$
  其中，$u_i=n_i{e}_i={\rho }_i{ϵ}_i$

  （$\Delta x$和$\Delta y$怎么不见了？）

  对于离子和电子而言，$u=\frac{3}{2}nkT$，如果有多种粒子同时存在， 则可以直接将它们的贡献相加

Boltzmann 输运方程与流体方程

  Boltzmann输运方程是普适的，描述了相空间分布的连续性：
$$\frac{\partial f(x_i,p_i,t)}{\partial t}+\sum _{i=1}^3\left(v_i\frac{\partial f}{\partial x_i}+{\dot{p}}_i\frac{\partial f}{\partial p_i}\right)=S$$
  其中，$S$为源函数，包含了所有非平滑的过程（例如：碰撞散射，与粒子产生、煙灭）

  这里有些深邃，没看懂

位力定理

推导过程

  从动量守恒出发：
$$\rho \frac{d^2\vec{r}}{d{t}^2}=-\nabla P-\rho \nabla \Phi$$
  两边点乘$\vec{r}$，并对空间积分：
$${\int }_m\vec{r}\cdot \ddot{\vec{r}}\rho dV=-{\int }_V\vec{r}\cdot \nabla PdV-{\int }_m\vec{r}\cdot \nabla \Phi dm$$
  先计算左边的项：
$$\vec{r}\cdot \ddot{\vec{r}}=\frac{d}{dt}(\vec{r}\cdot \dot{\vec{r}})-{\dot{\vec{r}}}^2=\frac{1}{2}\frac{d^2}{d{t}^2}(\vec{r}\cdot \vec{r})-{\dot{\vec{r}}}^2$$

$${\int }_{m}dm\vec{r}\cdot \ddot{\vec{r}}=\frac{1}{2}\frac{d^{2}I}{d{t}^2}-2KE$$

  其中，$I\equiv {\int }_m{r}^{2}dm$，$KE\equiv \frac{1}{2}{\int }_m{\dot{\vec{r}}}^{2}dm$

  随后计算右边的第一项：
$$\vec{r}\cdot \nabla P=\nabla \cdot (P\vec{r})-3P$$

$$\begin{array}{rl}-{\int }_V\vec{r}\cdot \nabla PdV& =-{\int }_V\nabla \cdot (P\vec{r})dV+{\int }_{V}3PdV\\ & =-\iint P\vec{r}\cdot d\vec{S}+{\int }_{V}3PdV\end{array}$$

  取常压强面$S$为边界，则有：
$$-\iint P\vec{r}\cdot d\vec{S}=-P_s\iint \vec{r}\cdot d\vec{S}=-P_s{\int }_V\nabla \cdot (\vec{r})dV=-3P_{S}V$$
  将物态方程假设为$P=(\gamma -1)u$（$\gamma$为绝热系数），则有：
$${\int }_{V}3PdV=3(\gamma -1)U$$
  其中，$U$为内能密度

  最后计算右边的第二项。$\Phi$是引力势，因此：
$$\Phi (\vec{r})=-G\int \frac{dm}{|\vec{r}-{\vec{r}}^{\prime }|}$$
  代入进去后计算得到的正是引力势能：
$$\Omega =-\frac{G}{2}\iint \frac{dmd{m}^{\prime }}{|\vec{r}-{\vec{r}}^{\prime }|}$$

位力定理的一般形式

  最后得到位力定理的一般形式：
$$\frac{1}{2}\frac{d^{2}I}{d{t}^2}=2KE+3(\gamma -1)U+\Omega -3P_{s}V$$
  取时间平衡后，对准稳态系统而言，$⟨\frac{d^{2}I}{d{t}^2}⟩=0$

  对于恒星的表明，取$P_S\approx 0$，而在流体静力平衡的情况下，$KE=0$，由此将一般形式简化为：
$$3(\gamma -1)U+\Omega =0$$
  对于非相对论性理想气体而言，$\gamma =5/3$，此时得到：
$$2U+\Omega =0$$

球对称稳态恒星结构与演化方程

几处近似假设

  孤立：天体不受其他天体的影响

  球对称：转动能相比于引力能而言忽略不计

  形成时化学均匀

球对称恒星结构演化方程组

质量守恒

$$\frac{1}{\rho }\frac{d\rho }{dt}=-\nabla \cdot \vec{v}$$

  改写为：
$$\frac{1}{\rho }\frac{d\rho }{dt}=-\frac{1}{r^2}\frac{d}{dr}(r^2{v}_r)$$

动量方程

$$\rho \frac{d\vec{v}}{dt}=-\nabla P-\rho \nabla \Phi$$

  改写为：
$$\begin{array}{rl}\frac{{\partial }^{2}r}{\partial t^2}& =-\frac{1}{\rho }\frac{\partial P}{\partial r}-\frac{Gm}{r^2}\\ & =-4\pi r^2\frac{\partial P}{\partial m}-\frac{Gm}{r^2}\end{array}$$
  其中：
$$m(r)={\int }_0^{r}4\pi r^{\prime 2}\rho (r^{\prime })d{r}^{\prime }$$

能量方程

$$c_P\frac{\partial T}{\partial t}-\frac{\delta }{\rho }\frac{\partial P}{\partial t}=q-\frac{\nabla \cdot \vec{F}}{\rho }$$

  改写为：
$$\frac{\partial l}{\partial m}=q-c_P\frac{\partial T}{\partial t}+\frac{\delta }{\rho }\frac{\partial P}{\partial t}=q+q_g$$

能量辐射运输方程

$$\vec{F}=-\frac{1}{3}{v}_T\lambda \nabla U$$

  改写为：
$$\frac{l}{4\pi r^2}=-\frac{c}{3\kappa \rho }\frac{d}{dr}(a{T}^4)$$
  其中，$a=7.556\times 10^{-15}ergc{m}^{-3}{K}^{-4}$，$\kappa$为质量吸收系数，单位质量的截面

重要的时标

动力学时标

  偏离流体静力平衡后，若$\gamma >4/3$便能够恢复，恢复平衡的速度非常快，这个过程所需要的时间就是动力学时标

  动力学时标可以看做是自由落体的时标，因此：
$$t_{dyn}\sim \frac{R}{v}\sim \sqrt{\frac{R^3}{GM}}\sim \frac{1}{\sqrt{G⟨\rho ⟩}}\approx 40{\left(\frac{R}{R_{\odot }}\right)}^{3/2}{\left(\frac{M}{M_{\odot }}\right)}^{-1/2}min$$

热时标

  恒星冷下来所需要的时间：热能比上表面辐射光度
$$t_{KH}\sim \frac{G{M}^2}{RL}=3\times 10^7{\left(\frac{L}{L_{\odot }}\right)}^{-1}{\left(\frac{R}{R_{\odot }}\right)}^{-1}{\left(\frac{M}{M_{\odot }}\right)}^{2}yr$$

核聚变时标

  核聚变时标是恒星核反应的持续时间，记产生核反应的核比例为$f$，核合成产能率为$\eta \approx 0.7\%$，则核聚变的时标为：
$$t_{nuc}=\frac{\eta fM{c}^2}{L}\approx 10^{10}{\left(\frac{M}{M_{\odot }}\right)}^{-2.5}yr$$
  三种时标中，动力学时标最小，随后是热时标，最长的是核聚变时标

稳态球对称流体静力平衡下基本方程

  静态平衡：所有的时间微分都等于0：
$$\vec{v}=0,\frac{d^2\vec{r}}{d{t}^2}\sim 0,\frac{dϵ}{dt}\sim 0,\frac{d\rho }{dt}\sim 0$$
  欧拉形式：
$$\{\begin{array}{rl}& \frac{dm}{dr}=4\pi r^2\rho \\ & \frac{dP}{dr}=-\rho \frac{Gm(r)}{r^2}\\ & \frac{dl}{dr}=4\pi r^2\rho q\\ & \frac{d}{dr}(a{T}^4)=-\frac{3\kappa \rho }{c}\frac{l}{4\pi r^2}\end{array}$$

  其中：

  式1为球的质量，式2为压强梯度，式3为能量梯度——两个壳层之间的核合成能量，式4为辐射方程

  将分母替换为$dm$（质量为自变量），则可以改写为拉格朗日描述：
$$\{\begin{array}{rl}& \frac{dr}{dm}=\frac{1}{4\pi r^2\rho }\\ & \frac{dP}{dm}=-\frac{Gm(r)}{4\pi r^4}\\ & \frac{dl}{dm}=q\\ & \frac{dT}{dm}=-\frac{3\kappa }{64{\pi }^{2}ac}\frac{l}{r^4{T}^3}\end{array}$$