二氧化硅气凝胶隔热复合材料的性能及其瞬态传热模拟
国 防 科 技 大 学 学 报
第31卷第2期 J OUR NAL OF NA TIONAL UNIVER SITY OF DEFENSE TECHNOL OGY
文章编号:1001-2486(2009) 02-0036-05Vol . 31No . 22009
二氧化硅气凝胶隔热复合材料的性能及其瞬态传热模拟
周祥发1, 2 , 冯 坚, 肖汉宁, 张长瑞, 姜勇刚1211
(1. 国防科技大学航天与材料工程学院, 湖南长沙 410073; 2. 湖南大学材料科学与工程学院, 湖南长沙 410082)
摘 要:高马赫数、新型航天飞行器承载严重的气动加热环境, 其隔热材料的性能与传热计算对飞行器热
防护结构的优化和预测具有一定的指导意义。研究了二氧化硅气凝胶隔热复合材料的结构与性能, 并进行了
瞬态传热模拟。结果表明复合材料热导率仅为0. 018W m ·K , 瞬态传热模拟结果与考核测试值吻合, 为预测和
优化隔热材料提供了一定的依据。
关键词:隔热; 二氧化硅气凝胶; 热导率; 有限差分法
中图分类号:TQ343 文献标识码:A
Performance and Heat Transfer Simulation of Silica
Aerogel Composites
ZHOU Xiang -fa 1, 2, FE NG Jian 1, XIAO Han -nin g 2, ZH ANG Chang -rui 1, JIANG Yong -gang 1
(1. College of Aerospace and M aterial Engineering , National Univ . of Defense Technology , Changsha 410073, China ;
2. College of Materials and Engineering , Hunan Univers ity , Changsha 410082, China )
A bstract :With the development of high -speed spacecrafts , the properties and heat transfer simulation of thermal protection materials used by the spacecrafts are of great significance . This study aims to stud y the properties of silica aerogel composites and heat transfer si mulation . The results indicate that thermal conductivit y of the silica aerogel composites is only 0. 018W m ·K and transient heat transfer si mulation shows good agreement with the experimental data . Thus the simulation method can forecast the experimental conclusions and optimize the geometric structure of the composites .
Key words :insulation ; silica aerogel ; thermal conductivity ; finite difference methods
新型航天飞行器向高马赫数、长航时方向发展, 飞行器表面累积气动加热非常严重。因此, 对高效隔热材料提出了更高的要求。对于被动式热防护系统, 二氧化硅气凝胶隔热复合材料展现了优异的隔
2热性能, 因为气凝胶具有三维纳米颗粒骨架(骨架颗粒2~5nm ) , 高比表面积(500~1200m g ) , 纳米级孔
洞(约20nm ) 等独特的结构[1-4]
[5], 能够有效抑制固态热传导和气体对流传热, 是目前热导率最低的固态材料(常温为0. 015W m ·K ) 。
为使隔热材料获得更为优异的隔热性能, 需要对材料结构与性能进行研究, 同时, 隔热材料的传热模拟计算使材料更具设计和预测性。Dar yabeigi 应用二热流近似方法分析了辐射换热, 建立了多层隔
[7]热结构的数值分析模型。Kamran 等采用能量方程和二热流密度方法对纤维隔热材料建立了传热模
型, Spinnler 等[8, 9][6]采用Scaling 模型对辐射传热进行求解并建立了模型, 程远贵等
[11-17][10]利用分形理论建立了纤维隔热材料的传热模型, 但是这些模型没有与实验相结合, 模型的解还需进一步确证。近年来, 美国军方为了解决飞行器承载与防热问题, 开发了具有承载和隔热双重功能的金属面板TPS 防热结
构, 并对金属面板TPS 防热结构做了一系列的重要测试和传热建模计算, NASA 的研究者编制了简化的一维程序用来设计金属及与之竞争的其它材料的热防护性能, 这个程序包含了每种热防护结构的一维非线性有限元热流近似方法, 分析辐射换热, 建立了多层隔热结构数值分析模型。
本研究测试和分析了二氧化硅隔热复合材料热物理性能; 采用有限差分隐式方法进行瞬态传热模 收稿日期:2008-11-07
基金项目:国防科技重点实验室基金资助项目1975, ,
周祥发, 等:二氧化硅气凝胶隔热复合材料的性能及其瞬态传热模拟37
拟计算并进行实验验证。
1 实验
1. 1 二氧化硅气凝胶隔热复合材料的制备
采用酸 碱两步法制备二氧化硅溶胶, 即将正硅酸乙酯(tetraethyl orthosilicate , TE OS ) 与一定量的乙醇混合后搅拌一段时间, 加入盐酸、水和乙醇进行酸催化, 静置一段时间后加入氨水、水和乙醇(ethanol , EtOH ) 进行碱催化, 制得二氧化硅溶胶。然后, 将一定体积分数的陶瓷纤维(主要成分:SiO 2, Na 2O ) 热处理后与氧化硅溶胶充分混合, 待其凝胶后老化一段时间, 使凝胶缩聚反应继续进行, 网络结构趋于完整, 再以乙醇为干燥介质(乙醇临界点243℃,6. 3MPa ) 经超临界干燥制得二氧化硅气凝胶隔热复合材料。
1. 2 测试
采用JSM -6360LV 扫描电镜观察复合材料的微观形貌; 用Hot Disk 导热仪测试二氧化硅气凝胶隔热复合材料热物性参数, 仪器型号:TPS 2500S , 样品尺寸:(Υ180mm ×20mm ) , 测量材料在不同温度下的热导率和比热容; 实验测试:热面温度800℃,隔热层冷面温度随时间的变化。
2 二氧化硅气凝胶隔热复合材料传热计算
2. 1 问题的描述及基本假设
飞行器在飞行过程中, 已知舱体内外壁温度随时间变化的曲线, 给定初始环境及条件, 要求计算不同厚度的隔热层内温度, 使得在工作时间内舱内空气温度在规定厚度的隔热层不高于给定温度。假定采用单隔热材料的隔热层结构在工作时间为600s 过程中最高温度不超过200℃。在单隔热材料结构中, 假设只存在简单的导热过程, 同时材料的厚度与另外两个方向的尺寸相比小得多。所以, 隔热层内的传热可以简化为沿厚度方向的一维瞬态导热问题。
在建立隔热层内的传热模型时, 提出以下假设:仅考虑厚度方向的热传导; 隔热层在加热之前各处温度相同; 加热过程中舱壁整体温度一致, 其升温方式如图1
。
Fig . 1 Temperature rising modes of the cabin wall Fig . 2 Mathematics -ph y s ical model 图1 舱壁升温方式图2 数学物理模型
2. 2 数学模型的建立
数学物理模型如图2所示, 隔热层与舱壁的接触处为原点O , 隔热层的厚度为δ, 用x 表示厚度方向坐标, t 表示加热时间, T (x , t ) 表示隔热层内温度。
隔热层内数学模型的导热微分方程式为[18]
T ρc T k x 其中T =T (x , t ) ; ρ为隔热材料的密度, c 为隔热材料的比热容, k 为隔热材料的热导率。
·初始条件
在t =0时, 隔热层处于均匀温度, 故可表示为, 02(1)
38国防科技大学学报 2009年第2期式中, t 0为隔热层内的初始温度。
·边界条件
前面假设隔热层表面与舱壁接触处温度相等, 舱壁温度曲线已知, 故隔热层外侧温度满足第一类边界条件, t ≥0时, 关系式可表示为
t wf =f (t )
式中, t w f 为隔热层表面温度。
隔热层内侧与舱内空气之间存在对流热交换, 满足第三类边界条件, 可表示为
-k [19]x =δ(3) =h (t w -t f ) (4) 式中, h 为对流换热表面传热系数, t w 为隔热层内侧温度; t f 为舱内空气温度。2. 3 导热问题数值解法
工程应用上要获得解析解是一个复杂的过程, 通常的处理方法是在工程允许的误差范围内采用数值计算得到它的近似解。本研究将用到有限差分法隐式差分格式来求解导热微分方程, 隐式差分格式不仅可提高计算的精度, 而且它是无条件稳定的。首先对导热的时间、空间区域进行离散化, x 为空间坐标, 可将计算区域划分为i 等份, 得到i +1个空间节点; t 为时间坐标, 可将时间坐标上的计算区域划分为j 等份, 得到j +1个节点。空间网格线和时间网格线的交点(i , j ) 代表了时间-空间区域中的
t i +1, j -t i , j , 其导数的有限差分表一个节点。用有限差商代替微商, 将微分方程转化为差分方程达式为t i +1, j =t i , j
t i -1, j =t i , j t t x
Δx i , j x t Δx i , j x i , j 22Δx i , j 2! x 2 t Δx x i , j 2! 323Δx +…i , j 3! Δx +…i , j 3! 33(5) t x t i +1, j -t i , j +O (Δx ) Δx (6)
(7)
(8) t x x
隐式差分格式向后差分式为2i , j i , j t i , j -t i -1, j +O (Δx ) Δx t i +1, j -2t i , j +t i -1, j 2+O (x ) Δx
t i -t i t i , k Δτk k -1
则
t i -t i t i -1-2t i +t i +1=a ΔτΔx
等价于:
t i k +1k k -1k k k (9) -t i t i -1-2t i +t i +1k +1a Δk +1k +1k =a 1+2t i t i -1+t i +1+t i Δx Δx Δx k k +1k +1k +1(10)
3 结果与讨论
3. 1 二氧化硅气凝胶隔热复合材料的微观结构
纤维经600℃热处理(去除表面有机浸润剂及吸附水) 后, 再将二氧化硅溶胶浸渍到定型的纤维模型中, 实验发现, 复合体的结构保持了较好的完整性, 二氧化硅溶胶渗入纤维空隙中, 填满了纤维与纤维,
周祥发, 等:二氧化硅气凝胶隔热复合材料的性能及其瞬态传热模拟39
结构如图3所示, 由于热处理改善了纤维与硅溶胶的浸润
性, 纤维中大部分空隙被气凝胶所填充, 同时陶瓷纤维表面
被气凝胶包裹, 形成了一个有机的整体。
3. 2 二氧化硅气凝胶隔热复合材料的热物性参数
经测试, 二氧化硅气凝胶隔热复合材料的密度仅为
0. 23g cm 。表1是hot disk 测试的不同温度下复合材料的
热导率和比热值, 从这些数据可知, 复合材料的常温热导率
低, 仅为0. 018W m ·K , 复合材料随温度的升高, 热导率变化
缓慢, 400℃热导率较常温时上升仅为0. 007W m ·K 。因
此, 制备的二氧化硅气凝胶复合材料不仅轻质而且具有优图3 二氧化硅气凝胶隔热复合材料的微观形貌异的隔热性能。Fig . 3 Schematic views and SEM micrographs
of aerogel composite
表1 二氧化硅气凝胶隔热复合材料的物性参数
Tab . 1 Thermal physical parameter of the silica aerogel composites
温度(℃)
热导率(W m ·K )
比热(J g ·K ) 100℃0. 018(常温) 0. 549200℃0. 0190. 528300℃0. 0210. 500400℃0. 0250. 4533
3. 3 气凝胶隔热复合材料的瞬态传热模拟
3. 3. 1 有限差分计算
为定量计算飞行器单隔热层复合材料瞬态热传导问题, 首先分析隔热层外壁温度的变化方式, 由舱内壁升温方式曲线(见图1) :升温阶段, 在120s 舱壁升温至800℃,即隔热层外壁温度在120s 升温至800℃;定温阶段, 隔热层外壁温度保持在800℃。首先, 利用有限差分隐式方法模拟分析升温阶段隔热层内不同节点(厚度) 温度随时间的变化。
2经计算, 热流密度q 0″为3160W m , 表面温度从初始温度
(28℃) T (0, 0) 经过120s 升至T (0, 120s ) , 即800℃,与实际的
隔热层外壁升温条件相符。所以, 对隔热层内进行节点划分,
如图4, 每个节点的厚度Δx 为2mm , Δt 为1. 2s 。初始温度
T (x , 0) 为28℃。编写计算程序, 用迭代法数值求解一维瞬态
导热问题, 软件计算采用复合材料的热物性数据, 见表1。
软件计算各节点温度随时间的变化曲线, 如图5。各节点
代表了隔热层内不同厚度处温度随时间的变化。图5中节点图4 隔热层内节点划分Fig . 4 Nodes for the composites [20]
T 0是隔热层外壁(0mm ) 温度随时间的变化, 可以发现在时间120s 壁面温度升至800℃,符合实际情况, T 1-T 9是隔热层内2mm -18mm 厚度处温度随时间的变化曲线。当隔热材料的厚度为18mm (节点
图5 节点T 0-T 9温度随时间的变化曲线Fig . 5 Temperature change curves with time of nodes T 0-T 9
T 9) , 升温阶段对其厚度的影响很小, 温度变化小于22℃,所以, 对于隔热材料厚度为小于20mm 要考虑升温阶段的影响, 但是当隔热材料厚度为大于或等于20mm , 为便于简化计算隔热层冷面温度与时间的关系, 可以忽略升温阶段对其厚度的影响。
3. 3. 2 定温阶段瞬态传热
经升温阶段的软件模拟分析, 当隔热材料的厚度在大于或等于20mm 时, 将条件假定为外壁定温阶段, 即隔热层外壁温度T s 保持在800℃。模型如图6, 编写计算程序, 用迭代法数值求解一维瞬态导热问题, 针对不同温度选取表1中的隔热材料热物性参数, 计算各时刻的冷面温度值并与实验测试值比较
。
图6 定温表面计算模型图7 软件计算值与测试值比较 Fig . 6 Modelin g of constant surface temperature Fig . 7 Comparison of calculated and tested
隔热层冷面各时刻的温度计算值与测试值比较如图7, 这种单隔热层结构的隔热材料在传热时间为600s , 800s , 1000s 时, 实验测试值分别为87℃,145℃,200℃;计算值分别为115℃,161. 4℃,201. 2℃,所以计算值与测试值吻合较好, 并且在给定工作时间(600s ) 内舱内空气温度在规定厚度(20mm ) 的隔热层不高于给定温度(200℃) , 该计算方法对隔热材料实际应用预测及结构优化奠定了一定的基础。4 二氧化硅气凝胶复合材料隔热性能的理论分析及计算模型的评价
气凝胶低密度、高孔隙率有效地降低了固体传导, 同时, 低于分子平均自由程的平均孔径有效地抑制了空气对流传热, 因此, 在常温下其热导率较低。但是, 氧化硅气凝胶对于波长为3~5μm 高温近红外热辐射具有较强的透过性, 而红外辐射产生的热导与温度的3次方成正比, 因此, 温度升高, 气凝胶的热导率上升较快。而当氧化硅气凝胶与纤维复合后, 一方面, 低热导的气凝胶降低了固体传热, 气凝胶纳米孔径抑制空气对流传热, 气凝胶与纤维的复合减少了纤维-纤维接触传热, 导致大部分固体传热通过低热导率的气凝胶传递, 同时气凝胶充分填充纤维间的空洞, 有效抑制了自由气体分子的传热; 另一方面, 纤维本身的挡红外能力较强, 提高了复合材料抗红外辐射能力。因此, 气凝胶复合材料在常温、高温都具有优异的隔热性能。
根据飞行器内气凝胶隔热复合材料的传热特点, 计算模型首先采用有限差分法分析计算了隔热层内温度与时间的变化规律(飞行器内单隔热材料传热模型忽略了辐射换热) , 然后, 模拟计算了冷面温度随时间的变化, 对飞行器内隔热部件在服役环境下传热行为的预测具有一定的参考和指导意义。5 结论
(1) 制备了二氧化硅气凝胶隔热复合材料, 其密度和常温热导率低, 分别为0. 23g cm , 0. 018W m ·K ; 复合材料随温度的升高, 热导率变化缓慢, 400℃热导率为0. 025W m ·K ;
(2) 计算模型采用有限差分法分析计算了隔热层内温度与时间的变化规律, 单层隔热材料传热模拟计算值与实验测试值吻合得较好, 对隔热材料性能的预测和设计提供了一定的依据。
(下转第69页) 3
参考文献:
[1] Currie A , Bro wn M A . Wide -s wath SAR [C ] IEE Proc . Inst . Elect . Eng . F . , 1992, 139(2) :122-135.
[2] Callaghan G D , Longs taff I D . Wide -s wath Space -borne SAR Us ing a Quad -element Array [C ] IEE Proc . -R adar Sonar Navig . , 1999, 146(3) :
159-165.
[3] Suess M . Graf mueller B , Az hn R . A Novel High R esolution Wide Swath SAR Sys tem [C ] Proc . IGAR SS ' 01, 13-19Jul y 2001, Sidney ,
Australia , 2001.
[4] Heer C , Soualle F , Zahn R , et al . Investigation on a Ne w High R es olution Wide Swath SAR Concept [C ] Proc . IGAR SS ' 03, Toul ouse , France ,
2003.
[5] 李世强, 杨汝良. 天线相位中心偏移方位多波束合成孔径雷达的误差分析[J ]. 电子学报, 2004, 32(9) :1436-1440.
[6] 马晓岩, 吴顺华, 向家彬. 速度与PR F 失配对MPCSAR 成像的影响及补偿方法研究[J ]. 电子学报, 2005, 33(12) :2130-2134.
[7] Aguttes J P . The SAR Train Concept :R equired Antenna Area Distributed over N Smaller Satellites , Increase of Perfor mance by N [C ] IGAR SS
2003, Toulous e , France , 2003:542-544.
[8] 杨凤凤, 王敏, 梁甸农. 基于非均匀采样的小卫星分布式多通道SAR 无模糊成像[J ]. 电子学报, 2007, 35(9) .
[9] 马仑, 李真芳, 廖桂生. 一种稳健的利用分布式小卫星获取宽域、高分辨SAR 图像的方法[J ]. 航空学报, 2007, 28(5) .
[10] Li Z F , Wang H Y , Su T , et al . Generation of Wide -s wath and High -resoluti on SAR Images fro m Multichannel Small Spaceborne SAR Systems [J ].
IEEE Geoscience and R emote Sensing Letters , 2005, 2(1) :82-86.
[11] Krieger G , Gebert N , Moreira A . Unambiguous SAR Si gnal Recons truction from Nonuniform Displaced Phas e Center Sampling [J ]. IEEE
Geoscience and Re mote Sens ing Letters , 2004, 1(4) :260-264.
[12] Goodman N A , Lin S C , et al . Proces sing of Multiple -receiver Spaceborne Arrays for Wide -area SAR [J ]. IEEE Trans . on GR S , 2002, 40:841
-852.
[13] 邢孟道, 等. 分布式小卫星雷达空时频成像方法研究[J ]. 宇航学报, 2005, 26:70-77.
[14] Lin Y P , et al . Periodicall y Nonuniform Sampling of Bandpass Signals [J ]. IEEE Transactions on Circuits and Sys tems —Ⅱ:Analog and Digital
Signal Process ing , 1998, 45(3) :340-351.
(上接第40页)
参考文献:
[1] Hrubes h L W . Aerogel applications [J ]J . Non -cryst . Solids , 1998, 225:335-342.
[2] Fes mire J E . Aerogel Insulation Systems for Space Launch Application [J ]. Cryogenics , 2006, 46:111-117.
[3] Reim M , Korner W , Manara J , et al . Sil ica AerogelGranulate Material for Thermal Ins ulation and Daylighting [J ]. Solar Energy , 2005, 79:131-
139.
[4] 沈军, 周斌, 吴广明, 等. 纳米孔超级绝热材料气凝胶的制备与热学特性[J ]. 过程工程学报, 2002, 2(4) :341-345.
[5] Carta D , Corrias A , Mountjoy G , et al . Structural Study of Highl y Porous Nanocompos ite Aerogels [J ]. J . Non -crys t . Solids , 2007, 353:1785-
1788.
[6] Kamran D . Thermal Analysis and Des ign of M ulti -layer Insul ation for Reentry Aerodynamic Heating [R ]. AIAA2001-2834, 2001.
[7] Kamran D . Desi gn of High Temperature Multilayer Insul ation for Reus able Launc h Vehicles [D ]. U niversit y of Virginia , 2000.
[8] Spinnler M . Studies on High -temperature Multilayer Thermal Ins ulations [J ]. Internati onal Journal of Heat and Mas s Transfer , 2004, 47:1305-
1312.
[9] Spinnler M . Theoretical Studies of High -temperature Multilayer Thermal Insul ations Us ing R adiation Scaling [J ]. J ournal of Quantitative Spectroscopy
and Radicalize Transfer , 2004, 84:477-491.
[10] 程远贵, 周勇, 朱家骅. 耐火纤维高温热导率的分形[J ]. 化工学报, 2002, 53(11) :1193-1192.
[11] 姚草根, 吕宏军, 贾新潮, 等. 金属热防护系统材料与结构研究进展[J ]. 宇航材料工艺, 2005, 35(2) :10-13.
[12] 夏德顺. 重复运载器金属热防护系统的述评[J ]. 导弹与航天运载技术, 2002, 256(2) :21-26.
[13] 关春龙, 赫晓东. 可重复使用热防护系统防热结构及材料的研究现状[J ]. 宇航材料工艺, 2003, 33(6) :7-11.
[14] 曹义, 程海峰, 肖加余, 等. 美国金属热防护系统研究进展[J ]. 宇航材料工艺, 2003, 33(3) :9-12.
[15] 韩杰才, 陈贵清, 孟松鹤, 等. 新型AR MOR 热防护系统[J ]. 宇航学报, 2004, 25(3) :350-353.
[16] 赵颖. 2000年世界运载器发展综述[J ]. 导弹与航天运载技术, 2001, 249(1) :16-22.
[17] Brewer W D , Wallace T A , Bird R K , et al . All oys and Coating Development for MetallicTPS for Reus able Launch Vehicles [C ] 2000National
Space Miss ile Materials Symposium , SanDiego , California , February 28-March 2, 2000.
[18] Cengel Y A . Heat Transfer , a Practical Approach (2nd edition ) [M ]. Boston :MCB McGra w -hill , 2003.
[19] Lienhard J H . A Heat Trans fer Textbook (3rd edition ) [M ]. Cambridge :Phlogis ton Pres s , 2003.
[20] Incropera F P , DeWitt D P . Fundamentals of Heat and Mass Trans fer (6th Edition ) [M ]. New York :J ohn Wiley &Sons , 2007.