平流层飞艇热敏感因素分析
Thermal sensitivity factors analysis of stratospheric airships
通讯作者:
收稿日期: 2018-11-20
Received: 2018-11-20
作者简介 About authors
程晨(1994—),女,硕士生,从事浮空器热特性研究.orcid.org/0000-0003-4657-116X.E-mail:
以平流层飞艇为研究对象,建立飞艇热特性模型并分析主要的热影响因素. 通过有限拆分法,建立完整的平流层飞艇热特性模型. 为了使建立的热模型能够适用于不同类型的飞艇,将其外形建模和表面离散化处理,分别采用分布参数法和集总参数法,对蒙皮单元和内部填充气体进行热特性数值计算. 通过相关试验数据对建立的模型及其求解方法的可靠有效性进行验证. 分别计算8种不同热影响因素作用下飞艇表面蒙皮温度及内部填充气体温度的变化,对每种因素的影响规律进行定性和定量的分析,得到飞艇热特性的主要敏感因素及不敏感因素,归纳出飞艇蒙皮出现极值温度的条件. 比较不同热源对飞艇热特性的影响和规律,给出飞艇热特性模型中的可简化部分.
关键词:
A thermal characteristics model of the airship was established and the main thermal influence factors were analyzed by taking stratosphere airship as the research object. A complete thermal model of the stratospheric airship was established with the finite split method. The shape was modeled and the surface was discretized in order to make the established thermal model applicable to different types of airships. The distribution parameter method and the lumped parameter method were used to calculate the thermal characteristics of the skin cell and the internal filling gas. The reliability and effectiveness of the established model and its solution method were verified by relevant experimental data. The changes in the surface temperature of the airship and the temperature of the internal filling gas under the influence of 8 different thermal influence factors were calculated, and the qualitative and quantitative effects of each factor were analyzed. The conditions for extreme temperature of the airship skin were summarized. The effects and laws of different heat sources on the thermal characteristics of airships were compared, and a simplified part of the thermal characteristics model of airships was given.
Keywords:
本文引用格式
程晨, 王晓亮.
CHENG Chen, WANG Xiao-liang.
平流层太阳辐射强而对流换热弱,因此飞艇蒙皮和内部气体的温度场受太阳辐射的影响极大. 驻空期间,飞艇蒙皮与内部填充气体之间的热平衡研究能够促进飞艇设计和控制的发展. 近年来,关于飞艇热力学模型的大多数研究都是以计算蒙皮及内部填充气体昼夜温差为主的结果导向型研究,所建立的模型相似,但是最终计算结果往往相差较大,例如内部填充气体的昼夜最大温差从35 K到60 K[1-5],蒙皮昼夜最大温差从65.8 K到98 K[1-5]. 计算结果的偏差是由飞艇热力学模型中的相关经验公式不一致、环境参数及所用材料不同以及分析方法不同等因素造成的. 在建立热模型的同时,要分析不同因素对飞艇热特性的影响规律,模型才具备工程通用价值和意义.
目前,飞艇热特性的主要研究方法可以分成两大类. 第1类为集总参数法,研究者将飞艇看成一个或几个温度均匀分布的部分进行研究,飞艇蒙皮及内部填充气体的热特性仅仅随着时刻的变化而变化,与具体的坐标位置无关,集总参数法计算结果较粗略,但是计算速度较快,适用于计算上升、下降时的飞艇状态,例如Carlson等[6]建立高空气球的轨迹和热特性模型;Farley[7]采用集总参数法,研究零压和超压浮空气球的热特性;Yao等[5]研究带有太阳能电池板的平流层飞艇在上升和下降过程中的热特性. 第2类研究方法为分布参数法,将研究对象网格化,分别计算每一个单元的热特性,这种分析方法相较于集总参数法的计算量更大,但是计算结果更准确,更适用于计算飞艇驻空状态时的热特性,例如Dai等[8]计算了高空气球的热特性,Zheng等[9]在计算时省略了飞艇内表面的辐射换热.
由于飞艇内部填充气体的温度分布相对均匀统一,表面蒙皮不同处的温差较大,本文根据不同对象不同的热特性差异,结合已有的2种热模型分析方法,提出有限拆分法,即采用集总参数法分析填充气体的热特性,采用分布参数法分析表面蒙皮的热特性. 有限拆分法因地制宜,优势诸多,相较于传统的集总参数法,提高了计算精度,相较于分布参数法,提高了计算效率. 该方法首先由计算流体动力学(CFD)软件建立所设计飞艇的三维几何模型,将飞艇蒙皮划分成若干三角形单元,采用建立的热模型计算飞艇驻空状态时的热特性.
1. 研究对象和数学模型
如图1所示,研究对象为平流层飞艇的蒙皮以及内部填充气体,影响飞艇热特性的主要外部环境因素可以分为3部分,分别为太阳短波辐射、长波辐射及对流换热.
图 1
1.1. 太阳短波辐射
飞艇表面蒙皮单元吸收的太阳短波辐射由3部分构成,分别是太阳直射辐射
式中:下标i表示蒙皮的第i个单元;
1.2. 长波辐射
长波辐射可以分成地球长波辐射
飞艇蒙皮内表面的辐射换热为
其中,
1.3. 对流换热
飞艇对流换热模型由内表面的自然对流换热
式中:
1.3.1. 内表面自然对流换热
飞艇内表面自然对流换热系数可以通过下式[16]计算:
式中:
其中
1.3.2. 外表面混合对流换热
飞艇外表面与大气既存在一定的相对运动,又有温度上的差别,因此飞艇外表面与空气之间的换热为混合对流换热,要综合考虑自然对流和强迫对流2种换热. 飞艇外表面对流换热系数可以通过下式[17]计算:
式中:
1.3.2.1. 外部自然对流换热
在计算飞艇外部自然对流换热时,可以将飞艇按截面分成许多个微元柱进行分析,根据流动特性的不同,将换热分为层流和湍流. 在自然对流初期,边界层厚度远远小于飞艇截面直径,换热近似为层流传热. 随着流动距离的增加,边界层厚度越来越大,与此相对应,表面传热系数越来越小. 当边界层厚度达到临界值的时候,层流传热将转变为湍流,边界层厚度保持恒定,不再与流动距离相关.
1)层流与湍流分界点 层流湍流的分界点对于外部自然对流换热的计算至关重要. 流体介质为空气,普朗特常数约为0.733 6.
图 2
图 2 圆柱体层流湍流转化分界点与瑞利数关系[21]
Fig.2 Relationship between demarcation point of laminar-turbulent flow and Rayleigh number
2)层流换热与湍流换热 对于直径为D的圆柱,壁面温度和来流温度均为常数,x为沿流线的距离. 如图3所示为外部自然对流截面参数图示.
图 3
外部自然层流换热与湍流换热的对流换热系数及努赛尔数的计算式如表1所示.
表 1 层流及湍流对流换热系数及努赛尔数计算式
Tab.1
换热方式 | 外部自然对流换热系数 | 努赛尔数 |
层流换热 | | |
湍流换热 | | |
1.3.2.2. 外部强迫对流换热
式中:
1.4. 飞艇热力学方程
根据能量守恒微分方程,可以推导出飞艇蒙皮单元的热力学方程:
内部填充气体的热力学方程为
式中:
式(15)、(16)建立了完整的飞艇热特性模型,可以用于迭代计算得到飞艇表面各蒙皮单元及内部填充气体随时间的变化关系.
2. 计算方法与模型验证
2.1. 飞艇热力学分析实现方法
通过CFD软件建立所设计飞艇的3D几何模型,将飞艇蒙皮划分成若干三角形单元;采用Tecplot软件对单元数据进行转换,形成Matlab仿真程序需要的飞艇模型数据. 整个热模型的仿真通过Matlab编程实现. 基于建立的飞艇热力学特性的数学模型,使用Matlab软件编程,实现整个热模型的仿真. 采用Tecplot软件及Origin软件处理计算结果,对结果进行一系列分析. 飞艇的热分析流程如图4所示.
图 4
2.2. 网格无关性及时间步长无关性分析
在数值仿真的计算程序中,有蒙皮单元网格数以及时间步长2个离散量. 为了在保证计算结果的准确性的前提下提高计算的效率,分析蒙皮网格数与时间步长对飞艇热特性计算的影响.
图 5
图 5 不同网格数下飞艇热特性48 h内的变化
Fig.5 Thermal characteristics of airship in 48 h under different grid numbers
图 6
图 6 不同时间步长下飞艇热特性48 h内的变化
Fig.6 Thermal characteristics of airship in 48 h under different calculation time steps
由图6可知,网格数对计算结果的影响极小;当时间步长在45 s以内时,时间步长对飞艇表面蒙皮及内部填充气体热特性的影响极小,但当时间步长达到60 s时,蒙皮最高超热温度在白天开始出现了振荡,可以判断60 s的时间步长对于计算模型来说过大,会导致计算结果的不准确.
在之后的计算中采用45 s作为计算的时间步长,4 666作为网格数,既能够保证计算结果的准确性,也能够提高计算效率.
2.3. 模型验证及结果分析
表 2 文献[15]实验的基本几何参数
Tab.2
基本参数 | 参数值 | 基本参数 | 参数值 | |
飞艇总长 | 1.42 m | 太阳辐射热流 | 972 W/m2 | |
中间圆柱直径 | 0.47 m | 蒙皮发射率 | 0.81 | |
蒙皮厚度 | 0.1 mm | 太阳吸收比 | 0.45 |
图7给出不同外部对流条件下,飞艇内部空气平均温度计算值Ts与李德富实验中飞艇内部上、下部分空气温度测量值T18、T17(第18和17号传感器)的对比.
图 7
图 7 内部填充气体平均温度的计算值与实验值对比
Fig.7 Comparison of calculated value and experimental value of temperature of internal filling gas
由图7可以得出,通过数值计算得出的飞艇内部空气温度的平均值在17、18号传感器测量的空气温度之间,与测量值吻合较好且变化趋势基本相同.
当来流速度为0 m/s时,300 s时上部空气温度(18号传感器)实验值为303.7 K,下部空气温度(17号传感器)实验值为290.0 K,因此实验中空气平均温度为296.9 K,空气平均温度计算值为298.7 K,比实验值高1.8 K,相差7.5%;600 s时上部空气温度实验值为311.7 K,下部空气温度实验值为292.6 K,因此实验中空气平均温度为302.2 K,空气平均温度计算值为304.6 K,比实验值高2.4 K,相差8.2%.
当来流速度为7 m/s时,600 s时上部空气温度实验值为294.3 K,下部空气温度实验值为289.0 K,因此实验中的空气平均温度为291.7 K,空气平均温度计算值为291.6 K,比实验值低0.1 K,相差0.5%;1 200 s时上部空气温度实验值为294.3 K,下部空气温度实验值为288.9 K,因此实验中的空气平均温度为291.6 K,空气平均温度计算值为292.3 K,比实验值高0.7 K,相差3.8%.
当来流速度为0 m/s时,计算值与实验值的最大偏差为8.2%;当来流速度为7 m/s时,计算值与实验值的最大偏差为3.8%. 考虑到李德富模拟实验中,太阳模拟器的误差约为±5%,风速仪的误差约为±5%,热电偶的测量误差约为±1%,综合考虑实验中存在的其他环境与人员操作误差,可以认为数值计算结果与试验结果吻合较好.
3. 飞艇热敏感因素分析
主要研究平流层飞艇的热敏感因素,以典型的椭球外形飞艇为研究对象,分析季节、云层覆盖率、飞艇尺度、填充气体种类、蒙皮厚度、空速及飞艇所处经纬度对于飞艇蒙皮温度及内部浮升气体热特性的影响. 本文数值计算的相关参数取值如表3所示.
表 3 飞艇的热敏感因素算例基本参数取值
Tab.3
基本参数 | 参数值 | 基本参数 | 参数值 | |
飞行高度 | 20 000 m | 表面单元个数 | 4 666 | |
飞艇长度 | 50 m | 飞艇最大直径 | 16 m | |
飞艇体积 | 6 704.5 | 蒙皮太阳辐射吸收率 | 0.33 | |
蒙皮比热容 | 3 600 J/(kg·K) | 蒙皮长波辐射吸收率 | 0.33 | |
蒙皮发射率 | 0.8 | 初始时刻蒙皮温度 | 216.5 K | |
地球发射率 | 0.92 | 初始时刻浮升气体温度 | 216.5 K |
3.1. 季节对于飞艇热特性的影响
太阳辐射极大地影响着平流层飞艇热特性,由于一年四季地球相对太阳的位置不断在变化,许多与太阳热源相关的角度在不同季节的变化特性各不相同,影响热源及飞艇热特性. 分析季节对于飞艇蒙皮及内部填充气体热特性的影响程度. 图8中,α为太阳方位角,β为太阳高度角.
图 8
图 8 不同太阳角度在不同季节的48 h内的变化
Fig.8 Change of solar angles in 48 h in different seasons
表 4 不同季节飞艇蒙皮的部分热特性指标
Tab.4
季节 | 最大温度 | 最大温度时段 | 最小温度 | 最小温度时段 |
春 | 289.86 K | 9 h-15 h | 193.3 K | 19 h-5 h |
夏 | 285.41 K | 9 h-15 h | 193.3 K | 19 h-5 h |
秋 | 289.12 K | 9 h-15 h | 193.3 K | 19 h-5 h |
冬 | 289.86 K | 9 h-15 h | 193.3 K | 19 h-5 h |
图 9
图 9 不同季节飞艇内部填充气体48 h内的温度变化
Fig.9
3.2. 云层对于飞艇热特性的影响
云层覆盖越多,会使得飞艇受到的地面及云层反射辐射越强,ρ为地面和云层反射系数,一般认为,晴天时ρ可取0.2,云层全覆盖时ρ可取0.7. 由图10可以看出,内部气体温度受云层的影响较大,云层全覆盖时的内部气体昼夜温差比晴天时高20.36 K,接近40%.
图 10
图 10 不同云层覆盖情况下浮升气体48 h内的温度变化
Fig.10
云层不仅影响内部氦气的热特性,而且影响蒙皮的热特性. 蒙皮的最低超冷温度不随着云层覆盖程度的变化改变(均为193.3 K),但是蒙皮的最高超热温度随着云层覆盖率的增加而增加.
由图11可知,云层全覆盖时的蒙皮最大超热温度比晴天时高13 K,因此设计飞艇及选取蒙皮材料时,应考虑到云层全覆盖的极端情形.
图 11
图 11 蒙皮最大超热随云层覆盖程度的变化
Fig.11
Relationship of
3.3. 飞艇尺寸对于飞艇热特性的影响
图 12
图 12 不同尺寸飞艇内部浮升气体48 h内的温度变化
Fig.12
图 13
图 13 内部气体达到最高温度时刻与飞艇尺寸的关系
Fig.13 Relationship between time when inner gas (helium) reaches maximum temperature and size of airship
由图13可知,随着飞艇尺寸的增大,氦气温度到达最大值的时刻从11.99时变成了12.83时,可知尺寸越大,热平衡实现得越慢,到达氦气超热最大值的时间越长. 这是由于飞艇尺寸变大时,到达换热平衡状态所需的时间越长.
除了影响内部气体的热特性,飞艇尺寸对蒙皮温度也有一定的影响.
图14表明,随着尺寸的变大,蒙皮最低超冷温度
图 14
图 14 最低超冷温度与飞艇尺寸的关系
Fig.14
Relationship between
3.4. 空速对于飞艇热特性的影响
图 15
图 15 不同空速时飞艇内部气体48 h内的温度变化
Fig.15
图 16
图 16 内部气体昼夜温差随空速的变化
Fig.16
Relationship between airspeeds and
从图17可以看出,随着空速的增加,飞艇蒙皮单元最高超热温度单调递减,最低超冷温度单调递增. 这是由于空速越大,外部对流换热越大,致使飞艇的温度越接近环境温度.
图 17
可以看出,空速是影响飞艇热特性的重要因素,在选择蒙皮材料的时候,需要综合考虑最小空速时蒙皮单元的超热极限及超冷极限温度.
3.5. 纬度、经度、蒙皮厚度及内部填充气体种类对于飞艇热特性的影响
纬度越低,太阳辐射越强,对飞艇热特性的影响越大. 计算结果表明,纬度变化对蒙皮热特性的影响不大,但对填充气体的热特性有一定的影响.
分析图18可得,随着纬度的升高,填充气体的最高温度和昼夜温差逐渐下降,纬度为0°时内部填充气体温差为64.79 K,纬度为60°时内部填充气体温差为53.14 K.
图 18
图 18 不同纬度时飞艇内部气体48 h内的温度变化
Fig.18
如图19所示,经度对于飞艇内部填充气体的热特性及蒙皮热特性基本没有影响. 算例中的飞艇所处经度为42E,分析发现飞艇在不同经度时具有完全重合的热源曲线及内部填充气体热特性曲线,在同一时间有最高蒙皮温度(290.56 K)及最低蒙皮温度(193.3 K),昼夜内部填充气体温差均为58.72 K. 经度不是飞艇的热敏感因素.
图 19
图 19 不同经度时飞艇内部气体48 h内的温度变化
Fig.19
如图20所示,蒙皮厚度t对于飞艇内部气体以及蒙皮热特性的影响极小. 分析发现不同蒙皮厚度的飞艇有着完全重合的热源曲线及内部填充气体热特性曲线,在同一时间有最高蒙皮温度(289.86 K)及最低蒙皮温度(193.3 K),昼夜内部填充气体温差均为61.5 K.
图 20
图 20 不同蒙皮厚度时飞艇内部气体48 h内的温度变化
Fig.20
除上述影响因素以外,分析不同种类的内部填充气体对飞艇热特性的影响. 由图21可见,当内部填充气体分别为氢气和氦气时,填充气体的温度基本相同,只在18 h左右气体温度曲线的斜率有微小的偏差. 可以得出,填充气体的种类对昼夜温差几乎没有影响,但是对温度变化的速率有细微影响. 当内部填充氦气和氢气时,飞艇蒙皮和内部填充气体的所有热特性指标的差值均小于1 K,可以初步判断内部填充气体的种类对于飞艇来说不是重要的热敏感因素.
图 21
图 21 飞艇不同内部填充气体48 h内的温度变化
Fig.21
4. 环境热源特性分析及模型简化
在影响飞艇热特性的全部8个热源中,太阳直射、天空散射以及地面及云层反射3个热源直接受到太阳的影响,其余热源不与太阳辐射直接相关,其中地面长波、大气长波、蒙皮内表面单元间辐射及外部对流4个热源主要由蒙皮表面温度决定,内部对流这个热源不仅与蒙皮温度相关,还与内部填充气体温度相关.
图 22
图 22 热源对飞艇蒙皮热特性影响占比
Fig.22 Percentage of contribution of heat source of airship
简化模型能够提高计算效率,同时基本不影响计算精度,适用于对计算速度要求很高的情况. 简化后的飞艇蒙皮单元的热力学方程为
省略天空散射辐射及蒙皮内表面热辐射后的简化模型的计算结果如图23所示. 内部填充气体温度较完整计算模型稍微偏低,但是相差始终保持在1 K之内,证明了简化模型的合理性.
图 23
图 23 简化模型与完整模型内部填充气体温度变化图
Fig.23
图 24
图 24 春季12 h主要热源及蒙皮温度云图
Fig.24 Main heat sources and skin temperature at 12 h in spring
5. 结 论
(1)飞艇最敏感的热影响因素为空速及云层,对于飞艇填充气体的昼夜温差,在空速为30 m/s 时比1 m/s时低33.46 K,云层全覆盖时比晴天时高20.36 K,接近40%. 季节、尺寸及纬度对飞艇热特性均有一定的影响,飞艇填充气体的昼夜温差,夏季比冬季高6 K,纬度为60度时比纬度为0度时低12 K;飞艇蒙皮最低超冷温度,6倍尺寸时比0.5倍尺寸时低6 K. 经度、蒙皮厚度以及填充气体种类对飞艇热特性基本无影响.
(2)当选取蒙皮材料时,应取云层全覆盖且空速最小时蒙皮的最高温度以及空速最小时蒙皮的最低温度作为材料的极值温度.
(3)在影响飞艇热特性的热源中,天空散射辐射及蒙皮内表面热辐射的热量贡献极小,仿真中可以忽略,形成新的飞艇热特性简化模型.
参考文献
平流层飞艇热仿真初步探讨
[J].DOI:10.3969/j.issn.1009-8518.2007.02.002 [本文引用: 1]
Preliminary discussion on thermal simulation of stratospheric airship
[J].DOI:10.3969/j.issn.1009-8518.2007.02.002 [本文引用: 1]
平流层飞艇的热数值分析
[J].
Thermal numerical analysis of stratospheric aerostats
[J].
Transient thermal behavior of stratospheric balloons at float conditions
[J].
A heat transient model for the thermal behavior prediction of stratospheric airships
[J].DOI:10.1016/j.applthermaleng.2014.05.050 [本文引用: 3]
New thermal and trajectory model for high-altitude balloons
[J].
Performance simulation of high altitude scientific balloons
[J].DOI:10.1016/j.asr.2011.12.026 [本文引用: 2]
A simplified thermal model and comparison analysis for a stratospheric lighter-than-air vehicle
[J].DOI:10.1115/1.4037194 [本文引用: 1]
蒙皮热辐射特性对平流层飞艇氦气温度影响
[J].
Effect of skin heat radiation characteristics on helium temperature of stratospheric aerostats
[J].
太阳电池对平流层飞艇热特性的影响分析
[J].
Analysis of the influence of solar cells on the thermal characteristics of stratospheric airship
[J].
A general method of obtaining approximate solutions to laminar and turbulent free convection problems
[J].
Extension of Emmons spot theory to flow on blunt bodies
[J].
/
〈 |
|
〉 |
