跳到主要内容
Sandbox Physics
风压与气体回落

风压与气体回落 · L2

从一阵风压算到盘面的恢复力

风除了要够强,还要有时间把气体推走。变化很慢的压力,可以与盘面上方最大的恢复力比较;短促的一阵风,则更适合先算它传来了多少动量。气体自身的响应时间,决定哪种比较更有用。

先选一小片气体

取气体盘里一块面积为AA、质量为mm的小片。质量除以面积得到柱密度Σg\Sigma_g,单位是千克每平方米。迎风面积相同,背后堆着的质量越多,风就越难把它加速。

坐标采用星系静止系。zz轴垂直盘面,正方向顺着入流;RR表示盘面内离中心的距离。风以正速度VV到来,气体片层自己也有法向速度vzv_z。决定碰撞动量交换的是两者的相对速度。

当片层顺风跑得比入流还快时,阻力会转向。完整表达式保留绝对值和带符号的速度差。只有片层速度远小于风速时,才适合直接用入流压强来估算受力。

路径偏置定义为q=b/rcq=b/r_c,其中bb是最近接团心的距离,rcr_c是环境的核心尺度。我们只计算零、0.37、1.5、2.72和四这五种偏置;换一条路径,会得到另一整段风压经历,星系速度仍相同。

Σg=mA,aram=ρ∣V−vz∣(V−vz)Σg\Sigma_g=\frac mA,\qquad a_{\rm ram}=\frac{\rho\lvert V-v_z\rvert(V-v_z)}{\Sigma_g}
ρ\rho是给定的入流密度。加速度等于力除以质量;固定面积相除后,质量通过柱密度进入分母。完整轨迹还保留了造星对气体质量的消耗。这条教学片层力律,把法向动量耦合系数取为一。

对应阅读:Original fixed-slice and static-restoring-force materials

盘面上方的引力并非越高越强

在光滑、上下对称的盘面正中,朝两侧的竖直引力互相抵消。稍微抬高气体,向下的恢复力逐渐增强;离星系足够远以后,引力又会减弱。在固定半径上,向下的加速度通常有一个最大值。

一股缓慢增强的风,可以先把片层移到新的平衡高度。如果压力超过这个半径所能提供的最大恢复力,它便无法再找到同样的局部静态平衡。这就是Gunn–Gott类压力比较的物理出发点。

我们用给定的核球、暗晕和恒星盘势,逐个初始半径寻找正高度上的最大恢复加速度。图中蓝线使用初始气体柱,金线则使用按既定耗尽规律演化到名义最近接时刻的气体柱。片层面积固定,气体转成恒星以后,柱密度会下降。

Prestore(R)=Σg(R)max⁡z>0∣∂Φ(R,z)∂z∣P_{\rm restore}(R)=\Sigma_g(R)\max_{z>0}\left|\frac{\partial\Phi(R,z)}{\partial z}\right|
Φ\Phi表示每单位质量的引力势能,并令无穷远处为零。这里求最大值时固定RR;真正运动的片层还会改变自己的半径。
初始与名义最近接气体柱的固定半径静态恢复压力
在6464个初始半径,独立数值搜索给定势的恢复力峰值。蓝线用初始气体柱,金线多乘exp⁡(−500/2000)\exp(-500/2000);横线为中央入流峰值。这是局部准静态参照,有限脉冲逃逸仍要跟随运动。 查看完整图像

对应阅读:Köppen, Jáchym, Taylor & Palouš, Ram Pressure Stripping Made Easy, arXiv:1806.05887v1 (2018) · Original fixed-slice and static-restoring-force materials

短风更像给了气体一脚

若风压最强的阶段很短,短到气体还来不及明显移动,可以把它看成一次推踢。把每个小时间段里的“压强乘以时间”加起来,就得到每单位面积收到的冲量。画在压强—时间图上,它就是曲线下的面积。

再用这个冲量除以气体柱密度,便得到最简单的速度增量估计。压强加倍、持续时间减半,在这个近似下给出相同的冲量。这里假设气体柱几乎没变,并暂时忽略推踢过程中引力与相对速度的变化。

例如,压强为10−13 Pa10^{-13}\,\mathrm{Pa},持续一亿年,作用在柱密度0.002 kg m−20.002\,\mathrm{kg}\,\mathrm m^{-2}的气体上。按本课程的时间单位,短脉冲估计约增加每秒一百五十八千米的速度。若把同样的风延长许多倍,气体已有时间响应引力,就要回到完整运动过程。使用这项估计的前提,是气体的竖直响应时间远长于这段一亿年的持续时间。

Δvz≃1Σg∫Pinc(t) dt\Delta v_z\simeq\frac{1}{\Sigma_g}\int P_{\rm inc}(t)\,dt
积分符号表示沿时间累加。这里估算的是短脉冲极限;完整轨迹没有用这条近似代替逐时受力。
Δt=100 Myr=3.156×1015 sΔvz≃10−13 Pa (3.156×1015 s)2×10−3 kg m−2=1.578×105 m s−1≃158 km s−1\begin{aligned}\Delta t&=100\,\mathrm{Myr}\\&=3.156\times10^{15}\,\mathrm s\\\Delta v_z&\simeq\frac{10^{-13}\,\mathrm{Pa}\,(3.156\times10^{15}\,\mathrm s)}{2\times10^{-3}\,\mathrm{kg}\,\mathrm m^{-2}}\\&=1.578\times10^5\,\mathrm m\,\mathrm s^{-1}\\&\simeq158\,\mathrm{km}\,\mathrm s^{-1}\end{aligned}
本模型采用1 Myr=3.156×1013 s1\,\mathrm{Myr}=3.156\times10^{13}\,\mathrm s。压强是每单位面积的力;乘上时间,再除以每单位面积的质量,最后得到速度单位。

对应阅读:Köppen, Jáchym, Taylor & Palouš, Ram Pressure Stripping Made Easy, arXiv:1806.05887v1 (2018) · Original fixed-slice and static-restoring-force materials

让气体的响应时间决定怎样比较

在盘面附近,光滑引力势常可近似成弹簧:高度增加一点,向下的加速度也按比例增加一点。气体于是有自己的上下振荡周期。风压若比这个周期变化得慢,气体有时间逐步寻找新的平衡;变化得很快时,先得到的主要是一份速度增量。

实际穿越常落在两者之间。风还在增强,气体已经向上走;片层速度改变阻力,造星消耗改变气体柱,引力也持续做功。要描述这段中间过程,就沿着运动方程继续积分。

先用静态曲线预测哪些气体更容易抬起,再用完整轨迹跟着它们走。局部最大恢复力、尾部最高处、最后仍受束缚的质量,分别描述同一次事件的不同部分。把它们的定义说清,才能连成一幅有用的物理图像。

对应阅读:Köppen, Jáchym, Taylor & Palouš, Ram Pressure Stripping Made Easy, arXiv:1806.05887v1 (2018) · Original fixed-slice and static-restoring-force materials

练习与解释性答案

  1. 两阵风的峰值压强相同,其中一阵持续时间长十倍,一定移走相同气体吗?

    展开解释

    不一定。短脉冲传来的冲量可能小得多;若两者都比竖直响应慢许多,静态压力比较才更适用。

  2. 气体柱密度减半,而面积和引力保持不变,静态恢复压力怎样变?

    展开解释

    恢复压力也减半,因为它等于柱密度乘最大恢复加速度。同样的入流压力也更容易加速较薄的气体柱。

动手检验预测

跟随迎风星系的气体

  1. 在同一原生时刻比较两条过团路径,先用风压经历预测高度,再看真实气体位置。
  2. 先猜离盘比例是否应等于比机械能为正的比例,再同时检查两项读数。

来源与计算范围

  1. Köppen, Jáchym, Taylor & Palouš, Ram Pressure Stripping Made Easy, arXiv:1806.05887v1 (2018)

    第2.1至2.2节及5.4节用于缓慢压力、冲量极限和可回落气体。自有三分量势、等质量片层与耗尽气体另有明确设定;静态参照没有重现原文全部模拟和观测样本。

  2. Original fixed-slice and static-restoring-force materials

    自有五份精确六十四乘十三字段合成切片、共同初态数组、两种全局光历史/三尘幕/四年龄波段子集,以及独立Python静态峰值和位置、能量、颜色诊断。原轨迹和原生恒星谱卷积是此前的独立计算;公开离线包复算诊断与波段比。没有浏览器引擎、原FITS/谱、论文、作者代码或实测星系数据。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策