📚 课程目录第3章空间外热流建模与计算

3. 空间外热流建模与计算

空间外热流是遥感载荷热控设计的输入边界条件，它决定了载荷外表面在轨接收到的能量水平，也直接影响散热面面积、加热回路功率和温度波动范围。工程上通常把外热流拆成四个独立分量：太阳辐射、地球反照、地球红外和深空冷背景。前三项是“热源”，最后一项是“冷汇”。本节围绕这四项的物理来源、工程计算方法以及典型轨道下的外热流曲线特征展开。

一、四个分量的物理来源与量级

太阳辐射是外热流中最主要的能量来源。在日地平均距离处，太阳常数约为 1367 W/m²，工程计算中常取 1367 W/m² 或 1353 W/m²（视任务约定）。太阳辐射到达载荷表面时，实际吸收的能量取决于表面太阳吸收率 α_s 和太阳入射角 θ，即 q_sun = S·α_s·cosθ。当 θ 超过 90° 时，表面处于阴影区，该项为零。

地球反照是太阳辐射经地球大气和地表反射后到达载荷的能量。它的大小与地表反射率、云层分布、太阳—地球—载荷几何关系密切相关。工程上常用反照率 ρ 乘以太阳常数再乘以一个与轨道高度和视角相关的几何因子来估算，典型反照率取 0.3 左右，但沙漠、云顶等区域可显著偏高。

地球红外是地球自身向外辐射的长波能量，主要来自地表和大气的热辐射。在低轨高度上，可近似把地球看作一个等效黑体，其红外辐射强度约为 237 W/m² 左右。载荷表面吸收的地球红外能量为 q_IR = E_IR·α_IR·F，其中 F 为地球对表面的角系数，α_IR 为红外吸收率。

深空冷背景是外热流中的“冷端”。深空等效温度约 4 K，辐射能量极低，工程上通常将其视为理想冷汇，即载荷向深空辐射散热时不受反向热流影响。这一项在散热计算中体现为辐射散热项 σ·ε·T⁴。

二、工程计算方法

外热流计算的核心是“几何关系 + 表面属性 + 时间积分”。对于任意一个外表面，其在某一时刻接收的总外热流可写为：

q_total(t) = S·α_s·cosθ(t)·δ_sun(t)
           + S·ρ·α_s·F_albedo(t)
           + E_IR·α_IR·F_IR(t)
           - σ·ε·T_surface⁴

其中 δ_sun(t) 为太阳可见性函数，进入阴影时取 0。F_albedo 和 F_IR 分别是地球反照和地球红外对目标面的角系数，它们随轨道位置变化。工程上常用两种方法求解：一是解析法，把地球简化为圆盘，用球面三角计算角系数；二是数值法，用射线追踪或有限元热分析软件（如 Thermal Desktop、NX TMG）逐面计算。

对于低轨遥感载荷，轨道高度通常在 500～800 km，轨道周期约 90～100 分钟，其中约 60% 时间处于日照区，40% 处于阴影区。太阳入射角在一个轨道周期内可从 0° 变化到 180°，导致太阳辐射项剧烈波动。地球反照和地球红外则相对平稳，但会随载荷姿态和地球视角变化。

工程上常把外热流曲线做成“轨道周期内的热流—时间”图，用于确定散热面朝向、加热回路启停阈值和热控涂层选型。曲线峰值通常出现在太阳直射且地球反照叠加的时刻。

三、典型轨道外热流曲线分析

以一颗 700 km 太阳同步轨道遥感载荷为例，假设其一个外表面法线指向天顶，太阳吸收率 α_s=0.2，红外发射率 ε=0.8，反照率取 0.3，地球红外强度取 237 W/m²。在一个轨道周期内，该表面的外热流曲线会呈现如下特征：

日照区初期：太阳入射角较大，cosθ 较小，太阳辐射项约 100～300 W/m²；地球反照和红外项合计约 80～120 W/m²。
正午附近：太阳接近直射，cosθ 接近 1，太阳辐射项可达 270 W/m² 以上；若同时处于地球反照高值区，总外热流可超过 400 W/m²。
阴影区：太阳辐射项归零，仅剩地球红外约 50～80 W/m²，表面温度开始下降。
出阴影后：太阳辐射迅速回升，表面温度开始爬升，形成周期性的温度波动。

这种周期性波动是热控设计的核心挑战。如果散热面朝向阳照区，其外热流峰值高，散热能力被削弱；如果朝向阴影区，则散热能力稳定但可能过冷。因此，工程上常通过调整散热面法线方向、选用低太阳吸收率涂层（如 α_s/ε 比值低的二次表面镜）来平衡。

四、计算中的关键约定与注意事项

外热流计算涉及方向、正负号和坐标系，必须先明确约定。通常定义表面法线方向为外法线，太阳入射角 θ 为太阳方向与表面法线的夹角，cosθ 为正表示太阳照射到表面，为负表示表面背向太阳。地球反照和地球红外的角系数 F 定义为地球辐射到表面的能量占地球总辐射的比例，取值范围 0～1。

工程计算中还需注意以下几点：

太阳常数取值：不同任务可能采用 1367 W/m² 或 1353 W/m²，需在热控设计文件中统一约定，避免计算偏差。
反照率时变性：反照率随季节、纬度、云量变化，工程上常取保守高值（如 0.35～0.4）用于热设计，取低值用于冷设计。
地球红外均匀性：实际地球红外分布并不均匀，极区、沙漠、云顶差异明显，高精度分析需采用地球红外地图数据。
表面属性退化：热控涂层在轨会因紫外辐照、原子氧侵蚀导致 α_s 增大、ε 变化，计算时应考虑寿命末期参数。

外热流计算结果是热控设计的输入，不是最终结论。实际在轨温度还受内部功耗、热容、接触热阻和加热回路控制策略影响。建议在初步设计阶段用外热流曲线确定散热面面积和加热功率，在详细设计阶段用热分析软件进行全轨道瞬态仿真验证。

总结来说，空间外热流建模的核心是：明确四个分量的物理来源，建立几何关系与表面属性的数学表达，结合轨道参数生成时间曲线，并在工程约定下完成计算。典型轨道外热流曲线呈现“日照区高、阴影区低”的周期性波动，峰值和谷值决定了热控设计的边界条件。掌握这一方法，才能为后续散热面布局、涂层选型和加热回路设计提供可靠依据。