CN105404774B

CN105404774B - 带太阳能电池的高空气球平飞过程分布温度计算方法

Info

Publication number: CN105404774B
Application number: CN201510771548.0A
Authority: CN
Inventors: 刘强; 祝榕辰; 才晶晶; 张衍垒; 杨燕初
Original assignee: Academy of Opto Electronics of CAS
Current assignee: Academy of Opto Electronics of CAS
Priority date: 2015-11-12
Filing date: 2015-11-12
Publication date: 2018-10-12
Anticipated expiration: 2035-11-12
Also published as: CN105404774A

Abstract

本发明提供了一种带太阳能电池的高空气球平飞过程分布温度计算方法，首先计算大气环境参数及高空气球辐射热环境参数，并基于高空气球几何特征及传热模式，建立高空气球分布温度计算域，然后利用结构化网格离散计算域，建立各微元的质量、动量和能量微分方程，最后根据高空气球球体材料和太阳能电池特性参数，联立求解计算域内所有微元的方程组，计算高空气球平飞过程分布温度。本发明在带太阳能电池的高空气球结构设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义，可以提高带太阳能电池的高空气球设计一次成功率，缩短带太阳能电池的高空气球设计周期，降低带太阳能电池的高空气球设计成本。

Description

带太阳能电池的高空气球平飞过程分布温度计算方法

技术领域

本发明属于高空气球热控制技术领域，尤其涉及一种带太阳能电池的高空气球平飞过程分布温度计算方法。

背景技术

高空气球具有可定点飞行、滞空时间长和分辨率高等优点，在空中预警、监视监测、民用通信等领域具有广泛应用前景，受到世界各主要强国的高度重视。

高空气球在平飞过程中，环境温度、密度、压力、风速、太阳辐射、大气辐射和地面辐射等因素会对高空气球温度特性产生影响。温度过高将提高高空气球内部氦气压力，对高空气球产生重要影响：1、温度过高将改变高空气球球体材料承力特性、增大高空气球球体热应力、增大高空气球球体张力，对高空气球球体的安全构成严重威胁；2、改变高空气球受力状况，导致高空气球飞行高度波动，干扰高空气球执行任务。因此，准确获知高空气球平飞过程中的温度特性，对高空气球结构设计、材料选择、飞行试验规划、规避潜在危险等方面具有重要意义，而目前还没有一个系统性地计算带太阳能电池的高空气球平飞过程分布温度的计算方法。

发明内容

(一)要解决的技术问题

本发明的目的在于，提供一种带太阳能电池的高空气球平飞过程分布温度计算方法，可快速而准确地获得带太阳能电池的高空气球平飞过程分布温度数据。

(二)技术方案

本发明提供一种带太阳能电池的高空气球平飞过程分布温度计算方法，包括：

S1，根据高空气球飞行任务需求，计算高空气球飞行参数及高空气球设计参数；

S2，测量球体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数；

S3，计算高空气球大气环境参数及高空气球辐射热环境参数；

S4，基于高空气球几何特征及传热模式，建立高空气球分布温度计算域，利用结构化网格离散计算域，建立各微元的质量、动量和能量微分方程；

S5，根据高空气球球体材料和太阳能电池特性参数，联立求解计算域内所有微元的方程组，计算高空气球平飞过程分布温度。

(三)有益效果

本发明可以快速和准确地获知带太阳能电池的高空气球平飞过程中的分布温度特性，在带太阳能电池的高空气球结构设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义，可以提高带太阳能电池的高空气球设计一次成功率，缩短带太阳能电池的高空气球设计周期，降低带太阳能电池的高空气球设计成本。

附图说明

图1是本发明实施例提供的带太阳能电池的高空气球结构示意图。

图2是本发明实施例提供的带太阳能电池的高空气球平飞过程分布温度计算方法流程图。

具体实施方式

本发明提供一种带太阳能电池的高空气球平飞过程分布温度计算方法，其根据高空气球飞行参数、高空气球设计参数、球体材料特性参数、太阳能电池特性参数及电池隔热材料特性参数，计算大气环境参数及高空气球辐射热环境参数，并基于高空气球几何特征及传热模式，建立高空气球分布温度计算域，利用结构化网格离散计算域，建立各微元的质量、动量和能量微分方程，根据高空气球球体材料和太阳能电池特性参数，联立求解计算域内所有微元的方程组，计算高空气球平飞过程分布温度。

根据本发明的一种实施方式，温度计算方法包括：

根据本发明的一种实施方式，高空气球飞行参数包括高空气球飞行时间、高空气球飞行地点经度Lon、高空气球飞行地点纬度Lat、高空气球飞行海拔高度h和高空气球飞行空速v；

高空气球设计参数包括高空气球体积V、高空气球长度L、高空气球最大直径D、高空气球表面积A和太阳能电池面积A_S。

根据本发明的一种实施方式，球体材料特性参数包括球体材料表面吸收率α、球体材料表面发射率ε、球体材料面密度ρ和球体材料比热容c；

太阳能电池特性参数包括太阳能电池效率η、太阳能电池表面吸收率α_S、太阳能电池表面发射率ε_S、太阳能电池面密度ρ_S和太阳能电池比热容c_S；

电池隔热材料特性参数隔热材料特性参数包括隔热材料厚度δ_{S_I}和隔热材料导热系数λ_{S_I}。

根据本发明的一种实施方式，高空气球大气环境参数包括高空气球飞行海拔高度h处的大气温度T_Atm、大气压力P_Atm和大气密度ρ_Atm，

其中，大气温度T_Atm的数学表达式为：

大气压力P_Atm的数学表达式为：

大气密度ρ_Atm的数学表达式为：

高空气球热环境参数包括高空气球辐射热环境参数和对流换热环境参数，所述高空气球辐射热环境参数包括太阳直接辐射热流q_{D_S}、大气散射太阳辐射热流q_{A_S}、地面反射太阳辐射热流q_{G_S}、大气长波辐射热流q_{A_IR}和地面长波辐射热流q_{G_IR}，

太阳直接辐射热流q_{D_S}的数学表达式为：

q_{D_S}＝I₀·τ_Atm，

其中，I₀为大气层上界太阳辐射强度，τ_Atm为太阳直接辐射衰减系数；

所述大气散射太阳辐射热流q_{A_S}的数学表达式为：

q_{A_S}＝k·q_{D_S}，

其中，k为大气散射系数；

地面反射太阳辐射热流q_{G_S}的数学表达式为：

q_{G_S}＝I_Ground·r_Ground·τ_{IR_G，}

其中，I_Ground为抵达地球表面太阳直接辐射强度，r_Ground为地球表面反射系数，τ_{IR_G}为地球表面辐射衰减系数；

所述大气长波辐射热流q_{A_IR}的数学表达式为：

其中，σ为辐射常数，T_Atm为大气温度；

地面长波辐射热流q_{G_IR}的数学表达式为：

其中，T_Ground为地面温度，ε_Ground为地面发射率；

根据本发明的一种实施方式，步骤S4包括：

建立高空气球及其外流场区域，利用结构化网格将计算域划分为多个微元，分析高空气球球体、太阳能电池、太阳能电池隔热材料、内部氦气微元传热过程，建立所有微元的质量、动量和能量微分方程；

其中，计算域内质量、动量和能量微分方程为：

质量微分方程：

动量微分方程：

能量微分方程：

其中，T是温度，ρ是密度，c_p是定压比热容，t代表时间，u代表流体速度矢量，k是导热系数，S_u代表动量广义源项，S_T代表能量广义源项，μ是流体的粘度系数，P是流体压力，X指代坐标向量；

其中，太阳能电池微元i的能量广义源项表达式：

S_{T_S,i}＝Q_{S,i_D}+Q_{S,i_Atm}+Q_{S,i_IR_Atm}+Q_{S,i_IR}+Q_{S,i_Cond}，

Q_{S,i_D}是吸收太阳直接辐射热量，Q_{S,i_Atm}是吸收大气散射辐射热量，Q_{S,i_IR_Atm}是吸收大气长波辐射热量，Q_{S,i_IR}是对外界环境长波辐射热量，Q_{S,i_Cond}是通过隔热层与球体的传导换热热量。

太阳能电池微元i的能量广义源项表达式中各项热量计算式列述如下：

吸收太阳直接辐射热量Q_{S,i_D}：

Q_{S,i_D}＝α_S·q_{D_S}·A_S,i·F_S-S，

其中，F_S-S是太阳能电池微元i外表面与太阳直接辐射的辐射角系数，A_S,i是太阳能电池微元i外表面面积。

吸收大气散射辐射热量Q_{S,i_Atm}：

Q_{S,i_Atm}＝αS·q_{IR_Atm}·A_S,i，

吸收大气长波辐射热量Q_{S,i_IR_Atm}：

Q_{S,i_IR_Atm}＝ε_S·q_{IR_Atm}·A_S,i，

对外界环境长波辐射热量Q_{S,i_IR}：

其中，T_S,i是太阳能电池微元i的温度。

通过隔热层与球体的传导换热热量Q_{S,i_Cond}：

其中，T_{Enup_S,j}是球体微元j的温度，球体微元j被太阳能电池微元i遮盖；

其中，球体上半部分被太阳能电池遮盖部分微元j的能量广义源项表达式：

S_{T_Enup_S,j}＝Q_{Enup_S,j_IR}+Q_{Enup_S,j_Cond}，

其中，Q_{Enup_S,j_IR}是吸收球体内部辐射换热热量，Q_{Enup_S,j_Cond}是通过隔热层与太阳能电池的传导换热热量。

球体上半部分被太阳能电池遮盖部分微元j的能量广义源项表达式中各项热量计算式列述如下：

吸收球体内部辐射换热热量Q_{Enup_S,j_IR}：

Q_{Enup_S,j_IR}＝A_{Enup_S,j}·(G_{Enup_S,j}-J_{Enup_S,j})，

其中，G_{Enup_S,j}是投射到球体上半部分被太阳能电池遮盖部分微元j的辐射热流，J_{Enup_S,j}是离开微元j的辐射热流。

其中，J_{Enup_S,j}可以表达为微元辐射热流和反射热流之和，其表达式：

其中，X_k,j是球体内表面微元k到球体上半部分被太阳能电池遮盖部分微元j的辐射角系数。

通过隔热层与太阳能电池的传导换热热量Q_{Enup_S,j_Cond}：

其中，T_{Enup_S,j}是球体上半部分被太阳能电池遮盖部分微元j的温度，A_{Enup_S,j}是球体上半部分被太阳能电池遮盖部分微元j的面积。

球体上半部分未被太阳能电池遮盖部分微元l的能量广义源项表达式：

S_{T_Enup_R,l}＝Q_{Enup_R,l_D}+Q_{Enup_R,l_Atm}+Q_{Enup_R,l_IR_Atm}+Q_{Enup_R,l_IR_E}+Q_{Enup_R,l_IR_I}，

其中，Q_{Enup_R,l_D}是吸收太阳直接辐射热量，Q_{Enup_R,l_Atm}是吸收大气散射辐射热量，Q_{Enup_R,l_IR_Atm}是吸收大气长波辐射热量，Q_{Enup_R,l_IR_E}是对外界环境长波辐射热量，Q_{Enup_R,l_IR_I}是与球体内部长波辐射换热热量。

球体上半部分未被太阳能电池遮盖部分微元l的能量广义源项表达式中各项热量计算式列述如下：

吸收太阳直接辐射热量Q_{Enup_R,l_D}：

Q_{Enup_R,l_D}＝α·q_{D_S}·A_{Enup_R,l}·F_{Enup_R,l-S}，

其中，A_{Enup_R,l}是微元l的面积，F_{Enup_R,l-S}是微元l与太阳直接辐射的辐射角系数。

Q_{Enup_R,l_IR_I}是与球体下半部分长波辐射换热热量。

吸收大气散射辐射热量Q_{Enup_R,l_Atm}：

Q_{Enup_R,l_Atm}＝α·q_{A_S}·A_{Enup_R,l}，

吸收大气长波辐射热量Q_{Enup_R,l_IR_Atm}：

Q_{Enup_R,l_IR_Atm}＝ε·q_{A_IR}·A_{Enup_R,l}，

其中，ε是球体材料发射率；

对外界环境长波辐射热量Q_{Enup_R,l_IR_E}：

与球体内部长波辐射换热热量Q_{Enup_R,l_IR_I}：

Q_{Enup_R,l_IR_I}＝A_{Enup_R,l}·(G_{Enup_R,l}-J_{Enup_R,l})，

其中，G_{Enup_R,l}是投射到微元l的辐射热流，J_{Enup_R,l}是离开微元l的辐射热流；

其中，球体下半部分微元m的能量广义源项表达式：

S_{T_End,m}＝Q_{End,m_D}+Q_{End,m_Atm}+Q_{End,m_G}+Q_{End,m_IR_Atm}+Q_{End,m_IR_G}+Q_{End,m_IR_E}+Q_{End,m_IR_I}，

其中，Q_{End,m_D}是吸收太阳直接辐射热量，Q_{End,m_Atm}是吸收大气散射辐射热量，Q_{End,m_G}是吸收地面反射辐射热量，Q_{End,m_IR_Atm}是吸收大气长波辐射热量，Q_{End,m_IR_G}是吸收地面长波辐射热量，Q_{End,m_IR_E}是对外界环境长波辐射热量，Q_{End,m_IR_I}是与球体内部长波辐射换热热量。

球体下半部分微元m的能量广义源项表达式中各项热量计算式列述如下：

吸收太阳直接辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{D_S}·A_End,m·F_End,m-S，

其中，A_End,m是微元m的面积，F_End,m-S是微元m与太阳直接辐射的辐射角系数；

吸收大气散射辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{A_S}·A_End,m，

吸收地面反射辐射热量Q_{End,m_G}：

Q_{End,m_G}＝α·q_{G_S}·A_End,m，

吸收大气长波辐射热量Q_{End,m_IR_Atm}：

Q_{End,m_IR_Atm}＝ε·q_{A_IR}·A_End,m，

吸收地面长波辐射热量Q_{End,m_IR_G}：

Q_{End,m_IR_G}＝ε·q_{G_IR}·A_End,m，

对外界环境长波辐射热量Q_{End,m_IR_E}：

与球体内部长波辐射换热热量

Q_{End,m_IR_I}＝A_End,m·(G_End,m-J_End,m)，

其中，G_End,m是投射到微元m的辐射热流，J_End,m是离开微元m的辐射热流。

根据本发明的一种实施方式，步骤S5包括，加载微元的热边界条件，通过微元之间能量数据传递，联立求解微元能量方程组，计算高空气球平飞过程分布温度分布数据。

综上所述，本发明可以快速和准确地获知带太阳能电池的高空气球平飞过程中的分布温度特性，在带太阳能电池的高空气球结构设计、材料选择、飞行试验规划、规避潜在危险等方面具有指导意义，可以提高带太阳能电池的高空气球设计一次成功率，缩短带太阳能电池的高空气球设计周期，降低带太阳能电池的高空气球设计成本。

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。

如图1所示，本发明实施例提供的带太阳能电池的高空气球包括高空气球由球体上半部分1、球体下半部分2、太阳能电池3、太阳能电池隔热层4、尾翼5和推进装置6构成。

其中，高空气球主体由球体上半部分1和球体下半部分2构成，球体上半部分顶部铺设有太阳能电池3，太阳能电池与球体上半部分之间安装有隔热层4，尾翼5呈倒Y型安装于高空气球尾部，推进装置6左右对称安装于高空气球两侧。

如图2所示，带太阳能电池的高空气球平飞过程分布温度计算方法，包括：

根据高空气球飞行任务需求，计算出本实施例中的高空气球主要飞行参数如表1所示，主要设计参数如表2所示。

表1高空气球主要飞行参数

表2高空气球主要设计参数

测量拟采用的高空气球球体材料特性参数如表3所示；测量太阳能电池特性与太阳能电池隔热材料特性参数如表4所示。

表3球体材料特性参数

表4太阳能电池与太阳能电池隔热材料特性参数

计算高空气球热环境：大气压力、温度、密度。其中，高空气球在海拔高度h处的大气温度T_Atm(K)、大气压力P_Atm(Pa)、大气密度ρ_Atm(kg/m³)可由公式计算：

大气温度随海拔高度h变化的数学表达式为：

大气压力随海拔高度h变化的数学表达式为：

大气密度随海拔高度h变化的数学表达式为：

计算太阳直接辐射热流q_{D_S}，大气散射太阳辐射热流q_{A_S}，地面反射太阳辐射热流q_{G_S}，大气长波辐射热流q_{A_IR}，地面长波辐射热流q_{G_IR}；对流换热环境参数包括高空气球与外部环境的对流换热系数h_Ex，高空气球与内部氦气的对流换热系数h_In。

太阳直接辐射热流q_{D_S}是大气层上界太阳辐射强度I₀与太阳直接辐射衰减系数τ_Atm的乘积，计算式如下：

q_{D_S}＝I₀·τ_Atm (4)

大气散射太阳辐射热流q_{A_S}是太阳直接辐射热流q_{D_S}与大气散射系数k的乘积，计算式如下：

q_{A_S}＝k·q_{D_S} (5)

地面反射太阳辐射热流q_{G_S}是抵达地球表面太阳直接辐射强度I_Ground、地球表面反射系数r_Ground与地球表面辐射衰减系数τ_{IR_G}的乘积，计算式如下：

q_{G_S}＝I_Ground·r_Ground·τ_{IR_G} (6)

大气长波辐射热流q_{A_IR}计算式如下：

其中，σ是辐射常数，T_Atm是大气温度。

地面长波辐射热流q_{G_IR}计算式如下：

其中，T_Ground是地面温度。

计算域内质量、动量和能量微分方程为：

质量微分方程：

动量微分方程：

能量微分方程：

其中，T是温度；ρ是密度；c_p是定压比热容；t代表时间；u代表流体速度矢量；k是导热系数；S_u代表动量广义源项；S_T代表能量广义源项；μ是流体的粘度系数；P是流体压力；X指代坐标向量。

建立各微元质量、动量和能量微分方程。其中，针对质量和动量微分方程，固体微元域内无流动，质量和动量微分方程退化；流体微元质量和动量微分通过联立能量微分方程一起求解。针对能量微分方程，固体微元的辐射热量、导热热量、内热源是其广义能量源项，添加广义能量源项作为边界条件即可建立完整的能量微分方程；流体微元与固体微元边界的对流换热通过质量微分方程、动量微分方程和能量微分方程联立求解。

太阳能电池微元i的能量广义源项表达式：

S_{T_S,i}＝Q_{S,i_D}+Q_{S,i_Atm}+Q_{S,i_IR_Atm}+Q_{S,i_IR}+Q_{S,i_Cond} (12)

吸收太阳直接辐射热量Q_{S,i_D}：

Q_{S,i_D}＝α_S·q_{D_S}·A_S,i·F_S-S (13)

吸收大气散射辐射热量Q_{S,i_Atm}：

Q_{S,i_Atm}＝α_S·q_{IR_Atm}·A_S,i (14)

吸收大气长波辐射热量Q_{S,i_IR_Atm}：

Q_{S,i_IR_Atm}＝ε_S·q_{IR_Atm}·A_S,i (15)

对外界环境长波辐射热量Q_{S,i_IR}：

其中，T_S,i是太阳能电池微元i的温度。

通过隔热层与球体的传导换热热量Q_{S,i_Cond}：

其中，T_{Enup_S,j}是球体微元j的温度，球体微元j被太阳能电池微元i遮盖。

球体上半部分被太阳能电池遮盖部分微元j的能量广义源项表达式：

S_{T_Enup_S,j}＝Q_{Enup_S,j_IR}+Q_{Enup_S,j_Cond} (18)

吸收球体内部辐射换热热量Q_{Enup_S,j_IR}：

Q_{Enup_S,j_IR}＝A_{Enup_S,j}·(G_{Enup_S,j}-J_{Enup_S,j}) (19)

通过隔热层与太阳能电池的传导换热热量Q_{Enup_S,j_Cond}：

S_{T_Enup_R,l}＝Q_{Enup_R,l_D}+Q_{Enup_R,l_Atm}+Q_{Enup_R,l_IR_Atm}+Q_{Enup_R,l_IR_E}+Q_{Enup_R,l_IR_I} (23)

吸收太阳直接辐射热量Q_{Enup_R,l_D}：

Q_{Enup_R,l_D}＝α·q_{D_S}·A_{Enup_R,l}·F_{Enup_R,l-S} (24)

Q_{Enup_R,l_IR_I}是与球体下半部分长波辐射换热热量。

吸收大气散射辐射热量Q_{Enup_R,l_Atm}：

Q_{Enup_R,l_Atm}＝α·q_{A_S}·A_{Enup_R,l} (25)

吸收大气长波辐射热量Q_{Enup_R,l_IR_Atm}：

Q_{Enup_R,l_IR_Atm}＝ε·q_{A_IR}·A_{Enup_R,l} (26)

其中，ε是球体材料发射率。

对外界环境长波辐射热量Q_{Enup_R,l_IR_E}：

与球体内部长波辐射换热热量Q_{Enup_R,l_IR_I}：

Q_{Enup_R,l_IR_I}＝A_{Enup_R,l}·(G_{Enup_R,l}-J_{Enup_R,l}) (28)

其中，G_{Enup_R,l}是投射到微元l的辐射热流，J_{Enup_R,l}是离开微元l的辐射热流。

球体下半部分微元m的能量广义源项表达式：

S_{T_End,m}＝Q_{End,m_D}+Q_{End,m_Atm}+Q_{End,m_G}+Q_{End,m_IR_Atm}+Q_{End,m_IR_G}+Q_{End,m_IR_E}+Q_{End,m_IR_I} (29)

吸收太阳直接辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{D_S}·A_End,m·F_End,m-S (30)

其中，A_End,m是微元m的面积，F_End,m-S是微元m与太阳直接辐射的辐射角系数。

吸收大气散射辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{A_S}·A_End,m (31)

吸收地面反射辐射热量Q_{End,m_G}：

Q_{End,m_G}＝α·q_{G_S}·A_End,m (32)

吸收大气长波辐射热量Q_{End,m_IR_Atm}：

Q_{End,m_IR_Atm}＝ε·q_{A_IR}·A_End,m (33)

吸收地面长波辐射热量Q_{End,m_IR_G}：

Q_{End,m_IR_G}＝ε·q_{G_IR}·A_End,m (34)

对外界环境长波辐射热量Q_{End,m_IR_E}：

与球体内部长波辐射换热热量

Q_{End,m_IR_I}＝A_End,m·(G_End,m-J_End,m) (36)

氦气压力控制范围为：

0≤ΔP_He＝P_He-P_Atm≤300Pa (37)

其中，ΔP_He是氦气超压量，P_He是氦气绝对压力，P_Atm是大气环境压力。

氦气质量控制：当高空气球内部氦气超压超过300Pa时候，氦气阀门打开，排出部分氦气，至超压量等于300Pa时阀门关闭。

氦气质量流量计算式为：

其中，ρ_He是氦气密度，A_{v_He}是氦气阀门面积，k_{v_He}是氦气阀门流量系数。

内部氦气温度和速度通过求解球体内部流体微元内质量、动量以及能量微分方程获得。

输入高空气球设计参数、飞行任务参数，加载微元的热边界条件，通过微元之间能量数据传递，联立求解微元能量方程组，计算高空气球平飞过程分布温度分布数据。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims

1.一种带太阳能电池的高空气球平飞过程分布温度计算方法，其特征在于，包括：

S5，根据高空气球球体材料和太阳能电池特性参数，联立求解计算域内所有微元的方程组，计算高空气球平飞过程分布温度,

所述步骤S4包括：

其中，计算域内质量、动量和能量微分方程为：

质量微分方程：

动量微分方程：

能量微分方程：

其中，太阳能电池微元i的能量广义源项表达式：

S_{T_S,i}＝Q_{S,i_D}+Q_{S,i_Atm}+Q_{S,i_IR_Atm}+Q_{S,i_IR}+Q_{S,i_Cond}，

Q_{S,i_D}是吸收太阳直接辐射热量，Q_{S,i_Atm}是吸收大气散射辐射热量，Q_{S,i_IR_Atm}是吸收大气长波辐射热量，Q_{S,i_IR}是对外界环境长波辐射热量，Q_{S,i_Cond}是通过隔热层与球体的传导换热热量；

吸收太阳直接辐射热量Q_{S,i_D}：

Q_{S,i_D}＝α_S·q_{D_S}·A_S,i·F_S-S，

其中，α_S是太阳能电池的太阳辐射吸收率，q_{D_S}是太阳直接辐射热流，F_S-S是太阳能电池微元i外表面与太阳直接辐射的辐射角系数，A_S,i是太阳能电池微元i外表面面积；

吸收大气散射辐射热量Q_{S,i_Atm}：

Q_{S,i_Atm}＝α_S·q_{IR_Atm}·A_S,i，

其中，q_{A_S}是大气散射太阳辐射热流；

吸收大气长波辐射热量Q_{S,i_IR_Atm}：

Q_{S,i_IR_Atm}＝ε_S·q_{IR_Atm}·A_S,i，

其中，ε_S是太阳能电池表面发射率；q_{A_IR}是大气长波辐射热流；

对外界环境长波辐射热量Q_{S,i_IR}：

其中，T_S,i是太阳能电池微元i的温度；

通过隔热层与球体的传导换热热量Q_{S,i_Cond}：

其中，λ_{S_I}是隔热层导热系数，δ_{S_I}是隔热层厚度，T_{Enup_S,j}是球体微元j的温度，球体微元j被太阳能电池微元i遮盖；

S_{T_Enup_S,j}＝Q_{Enup_S,j_IR}+Q_{Enup_S,j_Cond}，

其中，Q_{Enup_S,j_IR}是吸收球体内部辐射换热热量，Q_{Enup_S,j_Cond}是通过隔热层与太阳能电池的传导换热热量；

吸收球体内部辐射换热热量Q_{Enup_S,j_IR}：

Q_{Enup_S,j_IR}＝A_{Enup_S,j}·(G_{Enup_S,j}-J_{Enup_S,j})，

其中，G_{Enup_S,j}是投射到球体上半部分被太阳能电池遮盖部分微元j的辐射热流，J_{Enup_S,j}是离开微元j的辐射热流；

其中，X_k,j是球体内表面微元k到球体上半部分被太阳能电池遮盖部分微元j的辐射角系数；

通过隔热层与太阳能电池的传导换热热量Q_{Enup_S,j_Cond}：

其中，λ_{S_I}是隔热层导热系数，δ_{S_I}是隔热层厚度，T_{Enup_S,j}是球体上半部分被太阳能电池遮盖部分微元j的温度，A_{Enup_S,j}是球体上半部分被太阳能电池遮盖部分微元j的面积；

其中，Q_{Enup_R,l_D}是吸收太阳直接辐射热量，Q_{Enup_R,l_Atm}是吸收大气散射辐射热量，Q_{Enup_R,l_IR_Atm}是吸收大气长波辐射热量，Q_{Enup_R,l_IR_E}是对外界环境长波辐射热量，Q_{Enup_R,l_IR_I}是与球体内部长波辐射换热热量；

吸收太阳直接辐射热量Q_{Enup_R,l_D}：

Q_{Enup_R,l_D}＝α·q_{D_S}·A_{Enup_R,l}·F_{Enup_R,l-S}，

其中，α是球体材料太阳辐射吸收率，A_{Enup_R,l}是微元l的面积，F_{Enup_R,l-S}是微元l与太阳直接辐射的辐射角系数；

Q_{Enup_R,l_IR_I}是与球体下半部分长波辐射换热热量；

吸收大气散射辐射热量Q_{Enup_R,l_Atm}：

Q_{Enup_R,l_Atm}＝α·q_{A_S}·A_{Enup_R,l}，

吸收大气长波辐射热量Q_{Enup_R,l_IR_Atm}：

Q_{Enup_R,l_IR_Atm}＝ε·q_{A_IR}·A_{Enup_R,l}，

其中，ε是球体材料发射率；

对外界环境长波辐射热量Q_{Enup_R,l_IR_E}：

其中，T_{Enup_R,l}是微元l的温度；

与球体内部长波辐射换热热量Q_{Enup_R,l_IR_I}：

Q_{Enup_R,l_IR_I}＝A_{Enup_R,l}·(G_{Enup_R,l}-J_{Enup_R,l})，

其中，球体下半部分微元m的能量广义源项表达式：

其中，Q_{End,m_D}是吸收太阳直接辐射热量，Q_{End,m_Atm}是吸收大气散射辐射热量，Q_{End,m_G}是吸收地面反射辐射热量，Q_{End,m_IR_Atm}是吸收大气长波辐射热量，Q_{End,m_IR_G}是吸收地面长波辐射热量，Q_{End,m_IR_E}是对外界环境长波辐射热量，Q_{End,m_IR_I}是与球体内部长波辐射换热热量；

吸收太阳直接辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{D_S}·A_End,m·F_End,m-S，

吸收大气散射辐射热量Q_{End,m_Atm}：

Q_{End,m_Atm}＝α·q_{A_S}·A_End,m，

吸收地面反射辐射热量Q_{End,m_G}：

Q_{End,m_G}＝α·q_{G_S}·A_End,m，

吸收大气长波辐射热量Q_{End,m_IR_Atm}：

Q_{End,m_IR_Atm}＝ε·q_{A_IR}·A_End,m，

吸收地面长波辐射热量Q_{End,m_IR_G}：

Q_{End,m_IR_G}＝ε·q_{G_IR}·A_End,m，

其中，q_{G_IR}是地面长波辐射热流，

对外界环境长波辐射热量Q_{End,m_IR_E}：

其中，T_End,m是微元m的温度；

与球体内部长波辐射换热热量

Q_{End,m_IR_I}＝A_End,m·(G_End,m-J_End,m)，

2.根据权利要求1所述的温度计算方法，其特征在于，所述高空气球飞行参数包括高空气球飞行时间、高空气球飞行地点经度Lon、高空气球飞行地点纬度Lat、高空气球飞行海拔高度h和高空气球飞行空速v；

所述高空气球设计参数包括高空气球体积V、高空气球长度L、高空气球最大直径D、高空气球表面积A和太阳能电池面积A_S。

3.根据权利要求2所述的温度计算方法，其特征在于，所述球体材料特性参数包括球体材料表面吸收率α、球体材料表面发射率ε、球体材料面密度ρ和球体材料比热容c；

所述太阳能电池特性参数包括太阳能电池效率η、太阳能电池表面吸收率α_S、太阳能电池表面发射率ε_S、太阳能电池面密度ρ_S和太阳能电池比热容c_S；

所述电池隔热材料特性参数隔热材料特性参数包括隔热材料厚度δ_{S_I}和隔热材料导热系数λ_{S_I}。

4.根据权利要求3所述的温度计算方法，其特征在于，所述高空气球大气环境参数包括高空气球飞行海拔高度h处的大气温度T_Atm、大气压力P_Atm和大气密度ρ_Atm，

其中，大气温度T_Atm的数学表达式为：

大气压力P_Atm的数学表达式为：

大气密度ρ_Atm的数学表达式为：

所述高空气球辐射热环境参数包括太阳直接辐射热流q_{D_S}、大气散射太阳辐射热流q_{A_S}、地面反射太阳辐射热流q_{G_S}、大气长波辐射热流q_{A_IR}和地面长波辐射热流q_{G_IR}，

所述太阳直接辐射热流q_{D_S}的数学表达式为：

q_{D_S}＝I₀·τ_Atm，

所述大气散射太阳辐射热流q_{A_S}的数学表达式为：

q_{A_S}＝k·q_{D_S}，

其中，k为大气散射系数；

所述地面反射太阳辐射热流q_{G_S}的数学表达式为：

q_{G_S}＝I_Ground·r_Ground·τ_{IR_G}，

所述大气长波辐射热流q_{A_IR}的数学表达式为：

其中，σ为辐射常数，T_Atm为大气温度；

所述地面长波辐射热流q_{G_IR}的数学表达式为：

其中，T_Ground为地面温度，ε_Ground为地面发射率。

5.根据权利要求1所述的温度计算方法，其特征在于，所述步骤S5包括，加载微元的热边界条件，通过微元之间能量数据传递，联立求解微元能量方程组，计算高空气球平飞过程分布温度分布数据。