计算机辅助氧弹量热计探索火箭燃料的推进性能
Propulsion Performance of Rocket Fuel through Bomb Calorimeter and Computer Simulation
通讯作者:
收稿日期: 2019-11-18 接受日期: 2019-11-20
Received: 2019-11-18 Accepted: 2019-11-20
本实验基于现行“综合化学实验”进行创新性改进。将一定压力的氮气充入氧弹,在以硝酸钾为氧化剂、三氧化二铁为催化剂的条件下,混合一定比例的蔗糖得到硝糖燃料,进行燃烧,采用环境恒温量热计记录蔗糖燃烧过程中的温度变化,再通过雷诺作图法校正产生的ΔT偏差,最终计算得到硝糖燃烧过程的恒容热效应,即为硝糖燃料的爆热值。当蔗糖:硝酸钾:氧化铁的质量比为39:59:2时,其爆热值最大,比冲最大。在上述实验的基础上,选择最佳配比的硝糖燃料,用于计算机仿真模拟火箭发射系统。通过改变燃烧室压力和燃料流速,计算得到该最佳配比硝糖燃料的比冲和火箭飞行高度,进而得到采用单级火箭将东方红一号卫星送入预定轨道的发动机参数。改变每一级火箭的燃料类型与比例,设计得到适于推进不同卫星进入预定轨道的二级火箭和三级火箭。通过设计氧弹量热计的拓展应用与计算机仿真模拟相结合的实验,达到已有实验的创新设计、启发学生创造性和引起科研兴趣的目的。
关键词:
This work has been designed as an example in current comprehensive chemistry laboratory course. Firstly, nitrogen was filled into a bomb calorimeter, and sucrose was combusted in this calorimeter using KNO3 as the oxidant and ferric oxide as the catalyst. During combustion of sucrose, change of the temperature was measured through using isothermal calorimeter. Next, the deviation of ΔT was corrected using Reynolds method, and the heat of combustion of sucrose at constant volume was calculated, that is the heat of explosion of nitrosaccharose. When the mass ratio of sucrose-KNO3-ferric oxide was 39:59:2, both the heat of explosion and specific impulse were the largest. Based on the above results, the best proportion of nitrosaccharose fuel was selected for computer simulation of rocket launching system. A series of the specific impulse of nitrosaccharose fuel and the flying height of the rocket were obtained through changing the pressure of combustion chamber and flow rate of nitrosaccharose fuel. In addition, the engine parameters of single-stage rocket for transporting Dongfanghong-1 into the predetermined orbit were calculated through using this simulation system. The engine parameters of multiple-state rockets were also obtained by changing the types and ratio of rocket fuel. As a result, combination of bomb calorimeter experiments and computer simulation system greatly helped innovative design of the current experiment, enlightening creative thought of students.
Keywords:
本文引用格式
王玉琪, 钟子婷, 盘盈滢, 刘帅, 林昕元, 林耿忠, 庞浩然, 章伟光, 何广平.
Yuqi Wang.
1 引言
利用氧弹量热计测定燃烧热的实验,在国内外各高校的本科基础物理化学实验中普遍采用,如国外牛津大学[1]、华盛顿大学[2]等的物理化学实验教材,国内北京大学[3]、南京大学等[4]主编的物理化学实验教材均通过氧弹量热计测定萘的燃烧热。国内外高校也有将氧弹量热计应用于一些为高年级本科生开设的综合性实验,比如加州大学伯克利分校[5]将Lab view软件与氧弹量热计结合使用,成功校准并测定了糖样的燃烧热;朗沃德大学[6]将有机实验与物理化学燃烧实验相结合,指导学生合成生物柴油,测定生物柴油燃烧热,并使用UV-Vis分光光度计测量其浊点;南京大学出版的《大学化学实验综合与探索性实验》[7]中,通过开设呋喃甲醇、呋喃甲酸的合成及燃烧热测定实验,将有机合成与物理化学实验相结合。本参赛单位主编和使用的《综合化学实验》中,应用氧弹量热计观察添加二茂铁对柴油的助燃和消烟作用,并融入燃烧尾气的处理方法、计算机仿真模拟[8]等实验以创新实验内容。
航天飞行器的关键推进设备是火箭发动机,其中燃料又是关键之中的关键,是导弹、空间飞行器等各类发动机的动力源,其性能优劣直接影响发动机的工作性能。从黑火药到现在的高能燃料,燃料技术已经历了数百年的发展。一些新型燃料如端羟基聚丁二烯燃料(HTPB)等[9]的成功研制,在火箭燃料的发展史上具有划时代意义。
本创新实验基于本参赛单位主编和使用的《综合化学实验》中实验三“二茂铁及其衍生物对柴油的助燃和消烟作用”进行创新性改进,先采用氧弹量热计测定硝糖燃料的爆热值,找出最佳配比的硝糖燃料。再结合计算机仿真模拟,探索研究作为火箭发动机燃料的推进性能。改变燃烧室压力和燃料流速,计算偏二甲肼-四氧化二氮和煤油-液氧等不同种类和配比的燃料燃烧后推动火箭进入不同轨道所需要的比冲。最后,通过仿真模拟系统进行人工智能模拟计算与操作,判断装有该类型燃料的多级火箭能否把相关航天器送入预定轨道。
2 实验部分
2.1 实验原理
2.1.1 氧弹量热计
氧弹量热计按照工作方式可分为绝热、恒温和环境恒温三种,本实验所使用的氧弹量热计为环境恒温式热量计。图1是氧弹量热计测量装置,被测样品置于氧弹内,体系与外界隔以空气层绝缘,下方设有绝热缘垫片,上方设有绝热缘的胶板,以减少对流和蒸发。为了减少热辐射,稳定环境温度,体系与环境间包有温度与体系相近的空气隔热套,减少与环境的热交换。
图1
根据能量守恒定律,样品燃烧放出的热量促使量热计本身及周围的介质(本实验用水)温度升高,测量介质燃烧前后温度的变化,可求得该样品的恒容燃烧热,其关系式为:
式(1)中,QV为恒容燃烧热(J·g-1),W为量热计和周围介质的热当量(J·K-1),W一般用经恒重的标准物(苯甲酸)来标定,m1为点火丝的重量(g),m2为样品的重量(g),Q1为点火丝在恒容条件下燃烧放出的热量(J·g-1),ΔT为燃烧过程温度的变化。
利用热量计多功能控制箱记录燃烧过程中的温度变化,经过雷诺作图法校正由于系统热漏等原因产生的ΔT的测量偏差,通过式(1)可计算得出样品的恒容热效应。
2.1.2 固体燃料的爆热值与测定方法
燃料爆热值是表征固体燃料能量特性的重要物理参量之一,是燃料进行燃烧反应的热效应。
根据定容爆热的定义,1 kg的燃料在0 K、惰性介质(如N2)中绝热定容下燃烧,然后将生成的气体冷却到0 K,并假定在冷却过程中没有二次反应和燃烧产物的凝结过程(即成分不变),则所发出的热量定义为定容爆热(QV0)。通过弹式量热计可测定固体燃料的爆热值,大量实验结果证明,采用一般实验方法测定燃料爆热时,在常温下(20 ℃或室温)所测定得到的QV(水)可作为定容爆热(QV0) [10]。
2.1.3 爆热值与比冲的关系
燃料比冲的定义[11]:单位质量的燃料所能带来的冲量(动量的改变),单位为米/秒(m·s-1)或牛·秒/千克(N·s·kg-1),工程上习惯使用秒(s)。
燃料的爆热值与比冲之间具有重要关系,它反映了燃料能量的水平,关系到火箭发动机可以把多大的荷载推送到多高的轨道。爆热值和比冲一般呈正比关系,爆热值越高,比冲越大。装有这种燃料的火箭发动机的性能越好,即可以把更重的荷载推送到更高的地球轨道上。
2.1.4 计算机仿真模拟火箭发射系统
本模拟系统采用Labview (一种用图标代替脚本的框图程序语言)仿真软件。通过CpropepShell_cn软件(发动机热力学计算的基本工具)计算出部分相同燃料种类和配比,以及不同燃烧室压强的比冲数据。再将数据导入origin中进行数据拟合出新的函数方程。
利用不同比冲数据可以计算出火箭与荷载卫星的两段飞行高度,一是主动段高度,即从火箭发射到燃料燃尽,火箭脱落时的高度;二是被动段高度,即火箭脱落后,卫星抛射段高度。这两段高度可依据齐奥尔科夫斯基公式[12]:VmtIgln(m0/mF),进行推导计算。对多级火箭,按照一级火箭初速度V0为零,二级火箭初速度为一级火箭的末速度,三级的初速度为二级的末速度。只有当燃料推力大于火箭的总重,火箭才能起飞。在此条件下,依据比冲的定义,可以计算出火箭燃料燃烧的最小质量流速。具体公式的推导与计算过程如下。
由齐奥尔科夫斯基公式得:Vmt=Igln(m0/mF)
随着时间的增加,燃料不断被消耗,因此设为一定时间t内,消耗燃料的质量为mt × t,因此列出下式:
${V_{{\rm{mt}}}} = Ig\ln [{m_0}/({m_0} - {m_{\rm{t}}}t)]$
主动段升空高度H1计算如下:
$\begin{aligned} H_{1} &=\int_{0}^{t} V_{\mathrm{mt}} \mathrm{d} t \\ &=I g \int_{0}^{t_{0}} \ln m_{0}-\ln \left(m_{0}-m_{\mathrm{t}} t\right) \\ &=\operatorname{Igm}_{0}\left\{\left(1-m_{\mathrm{燃}} / m_{0}\right) \ln \left[\left(m_{0}-m_{\mathrm{燃}}\right) / m_{0}\right]+m_{\mathrm{燃}} / m_{0}\right\} \end{aligned}$
燃料燃尽后抛射段高度推导如下。
该级火箭燃料燃尽后,末速度为:
${V_末} = Ig\ln ({m_0}/{m_{\rm{F}}}) = Ig\ln [{m_0}/({m_{\rm{火箭}}} + {m_{\rm{荷载}}} - {m_{\rm{燃}}})]$
之后火箭做抛射运动,对应高度为H2,计算如下:
$\begin{array}{l}0.5mV_{\rm{末}}^2 = mg{H_2}\\{H_2} = 0.5V_{\rm{末}}^2/g = 0.5g{\{ I\ln [{m_0}/({m_{\rm{火箭}}} + {m_{\rm{荷载}}} - {m_{\rm{燃}}})]\} ^2}\end{array}$
故火箭的升空总高度H:
$H = {H_1} + {H_2}$
式中:I为燃料比冲(s);mt为单位时间内燃料消耗量(kg·s-1);g为重力加速度(N·kg-1);m0为初始质量(m0= m火箭+ m荷载) (kg);mF为末质量(kg);m火箭为火箭总质量(kg);m燃料为火箭燃烧初始质量(kg);Vmt为火箭瞬时速度(m·s-1)。
通过选择不同燃料种类、配比和质量流速,系统可自动计算比冲,根据这些参数,通过仿真模拟系统进行人工智能模拟计算与操作,判断装有该类型燃料的多级火箭能否把相关航天器送入预定轨道。
2.2 仪器与试剂
2.2.1 试剂
硝酸钾(AR,广州化学试剂厂),蔗糖(AR,广州化学试剂厂),三氧化二铁(AR,上海麦克林生化科技有限公司)。
2.2.2 仪器
电子天平(AL型,精确至0.1 mg,梅特勒-托利多仪器上海有限公司),Pro’s Kit MT-1820万用电表(中国台湾宝工工具有限公司),氮气瓶(纯度99.99%,广州市华诺工业气体有限公司),HR-15A型数显氧弹量热计(湖南长沙长兴高教仪器设备公司),HR-15B氧弹量热计多功能控制箱(长沙长兴高教仪器设备公司)。
2.3 实验步骤
2.3.1 实验预习
图2
①样品的配制:实验组学生分为A、B、C三小组,每小组按照下述比例分别称量总质量约为20 g的燃烧样品,其中蔗糖:硝酸钾:氧化铁质量比分别为:A组29 : 69 : 2,B组39 : 59 : 2,C组49 : 49 : 2。各小组将称量好的样品倒入研钵中混合,研磨8 min,再装入称量瓶中备用。
②装样与称重:剪取一段长约30 cm的铁丝,将铁丝中间用玻璃棒绕成螺旋数圈,准确称量其质量记为m1。将点火丝置于坩埚内,用药勺将样品装入坩埚,盖住铁丝,注意铁丝不得触碰坩埚壁,样品称取5.0000 ± 0.0500 g,精确至0.0002 g,记质量为m2。
①充氮气:将装有样品的坩埚放入氧弹中,铁丝的两端分别接在氧弹的两个接线柱上,旋紧盖子,氧弹内充入1.4 MPa氮气,保持10–15 s,排出氧弹内的气体,再向氧弹内充入1.4 MPa氮气,达到所需压力后保持10–15 s。
②回路电阻检测:充好氮气后,用万用表检查两接线柱间的电阻,检查是否通路,若通路则进行下一步操作,否则放气打开氧弹,检查断路原因,重新进行此步操作。
③氧弹量热计安装调试:将氧弹置于内筒的固定座上,将内筒放入外筒的绝缘支架上,将热量计多功能控制箱的电极引线接入氧弹的两个电极上。准确量取2800 mL自来水倒入内筒中,浸没氧弹排气阀门螺母高度的三分之二处,盖上筒盖,将温度计的热敏探头插入内筒水中。
温度稳定后开始读点火前最初阶段的温度,每隔0.5 min记录一次温度示数,读取6个数据后,在记录第6个读数的同时,迅速通电点火,指示灯亮后瞬时熄灭,温度迅速持续上升,表明点火成功。自点火后,每隔0.5 min读一次温度读数,至两次读数之差小于0.005 ℃后,仍为每0.5 min读一次数,继续记录约3 min,作为燃烧反应的后期温度。取出燃烧后剩余的铁丝,准确称量,记为m3。将实验数据记录到表1中。每个实验配比进行3次平行实验。
表1 计算机辅助氧弹量热计探索火箭固体燃料推进性能实验记录表
小组成员:________________实验日期:________________ | ||||
实验室温度(℃)________________大气压(kPa):________________ | ||||
硝糖固体燃料样品配比(蔗糖:硝酸钾:氧化铁)________________ | ||||
实验-1 | 实验-2 | 实验-3 | ||
点火丝质量m1/g | ||||
样品质量m2/g | ||||
剩余点火丝质量m3/g | ||||
记录(次) | 时间(分) | 温度/℃ | 温度/℃ | 温度/℃ |
1 | 0.5 | |||
2 | 1.0 | |||
… | … | |||
n | t |
根据2.3.2小节不同配比硝糖爆热值的测定,筛选出硝糖燃料最佳配比。在计算机模拟系统中,通过调节火箭燃烧室压力以及燃料质量流速(由火箭发动机喷嘴大小决定),探索火箭运行情况,能否起飞,起飞后可把东方红一号送入多高的轨道。将数据记录到表2中。
表2 计算机模拟参数记录表
记录:_________________实验日期:_________________ | |||||
硝糖固体燃料样品最佳配比(蔗糖:硝酸钾:氧化铁) _________________ | |||||
序号 | 燃烧室压力/MPa | 质量流速/(kg·s-1) | 比冲/s | 飞行高度/km | 运行情况 |
1 | |||||
...... | |||||
n |
通过选择一、二级火箭的燃料,调节相应燃料的含量、相关燃料质量流速,探索火箭运行状况,能否起飞,起飞后可把亚太五号通讯卫星送入多高的轨道。将计算机模拟数据记录到表3中。
通过选择一、二、三级火箭的燃料,调节相应燃料的含量、相关燃料质量流速,探索火箭运行状况,能否起飞,起飞后可把嫦娥一号送入多高的轨道。将计算机模拟数据记录到表4中。
表4 三级火箭模拟参数记录表
序号 | 一级火箭 | 二级火箭 | 三级火箭 | 质量流速/ (kg·s-1) | 飞行高度/km | |||||
燃烧剂种类 (含量/%) | 比冲/s | 燃烧剂种类 (含量/%) | 比冲/s | 燃烧剂种类 (含量/%) | 比冲/s | |||||
1 | ||||||||||
... | ||||||||||
n |
将表1记录的数据导入Origin中,绘制T–t图并进行雷诺校正,求出由于固体燃料燃烧引起的卡计温度变化值ΔT,代入公式(1)中即可求得该固体燃料的爆热值QV。
重复2.4.1步骤,求出不同配比的硝糖燃料的爆热值,并以硝糖燃料中蔗糖含量为横坐标,爆热值为纵坐标,绘制不同配比的硝糖燃料与爆热值关系曲线图,得出硝糖燃料最佳配比。
3 结果与讨论
3.1 硝糖爆热值的计算
将测定数据进行雷诺校正后代入公式(1),计算求得不同比例的硝糖燃料爆热值,并进行误差分析,如表5所示。实验中每组样品平行测定3次,组内相对误差均在2%以下,具有较好的平行性,实验结果可靠。
表5 不同比例的硝糖燃料爆热值表
蔗糖含量ω/% | 次数 | 燃料质量m/g | 爆热值QV/(J·g-1) | 绝对误差/J | 相对误差/% |
29 | 1 | 5.0084 | 2292 | 28 | 1.3 |
29 | 2 | 5.0013 | 2220 | -42 | -1.9 |
29 | 3 | 5.0055 | 2277 | 14 | 0.6 |
平均值 | 2264 | ||||
39 | 1 | 5.0020 | 2580 | 22 | 0.9 |
39 | 2 | 5.0162 | 2529 | -28 | -1.1 |
39 | 3 | 5.0036 | 2563 | 5 | 0.2 |
平均值 | 2558 | ||||
49 | 1 | 5.0009 | 2257 | -1 | -0.04 |
49 | 2 | 5.0005 | 2260 | 2 | 0.1 |
49 | 3 | 5.0000 | 2257 | -1 | -0.04 |
平均值 | 2258 |
3.2 硝糖配比对其爆热值的影响
图3
3.3 不同燃烧室压力对硝糖燃料比冲的影响
表6 单级火箭模拟参数记录表
序号 | 燃烧室压力/MPa | 比冲I/s | 质量流速/ (kg·s-1) | 飞行高度/km | 运行情况 |
1 | 1.529 | 120.00 | 25 | 0 | 不能起飞 |
2 | 1.529 | 120.00 | 40 | 179 | 无法到达近地轨道 |
3 | 2.753 | 130.00 | 25 | 232 | 达到近地轨道 |
4 | 2.753 | 130.00 | 40 | 207 | 达到近地轨道 |
5 | 5.203 | 140.00 | 25 | 263 | 达到近地轨道 |
6 | 5.203 | 140.00 | 40 | 236 | 达到近地轨道 |
7 | 7.346 | 145.00 | 25 | 280 | 达到近地轨道 |
8 | 7.346 | 145.00 | 40 | 252 | 达到近地轨道 |
9 | 10.636 | 150.00 | 35 | 274 | 达到近地轨道 |
10 | 10.636 | 150.00 | 40 | 267 | 达到近地轨道 |
硝糖燃料样品最佳配比(蔗糖:硝酸钾:氧化铁) = 39:59:2
图4
3.4 多级火箭燃料性能的计算机模拟
表7 二级火箭模拟参数记录表
序号 | 一级火箭燃料 | 二级火箭燃料 | 质量燃速/ (kg·s-1) | 飞行高度/km | 运行情况 | |||
燃烧剂种类 (含量/%) | 比冲/s | 燃烧剂种类 (含量/%) | 比冲/s | |||||
1 | UDMH (9.3) | 200 | 煤油(6.5) | 200 | 3500 | 2683 | 达到中地球轨道 | |
2 | UDMH (9.3) | 200 | 煤油(6.5) | 200 | 5000 | 2637 | 达到中地球轨道 | |
3 | UDMH (25.3) | 285 | 煤油(24.8) | 325 | 3500 | 4989 | 达到中地球轨道 | |
4 | UDMH (9.3) | 200 | UMDH (9.3) | 200 | 3500 | 2680 | 达到中地球轨道 | |
5 | UDMH (9.3) | 200 | UMDH (9.3) | 200 | 5000 | 2634 | 达到中地球轨道 | |
6 | UDMH (25.3) | 285 | UMDH (25.3) | 285 | 3500 | 4375 | 达到中地球轨道 | |
7 | 煤油(6.5) | 200 | UMDH (9.3) | 200 | 3500 | 2680 | 达到中地球轨道 | |
8 | 煤油(6.5) | 200 | UMDH (9.3) | 200 | 5000 | 2634 | 达到中地球轨道 | |
9 | 煤油(23.0) | 320 | UMDH (25.3) | 285 | 3500 | 4625 | 达到中地球轨道 | |
10 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 3500 | 2674 | 达到中地球轨道 | |
11 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 5000 | 2634 | 达到中地球轨道 | |
12 | 煤油(23.0) | 320 | 煤油(23.0) | 320 | 3500 | 5172 | 达到中地球轨道 |
UDMH表示偏二甲肼
表8 三级火箭模拟参数记录表
序号 | 一级火箭 | 二级火箭 | 三级火箭 | 质量燃速/ ((kg·s-1) × 104) | 飞行高度/km | 运行情况 | |||||
燃烧剂种类 (含量/%) | 比冲/s | 燃烧剂种类 (含量/%) | 比冲/s | 燃烧剂种类 (含量/%) | 比冲/s | ||||||
1 | UDMH (9.0) | 200 | UDMH (9.0) | 200 | UDMH (9.0) | 200 | 10 | 0 | 无法起飞 | ||
2 | UDMH (9.0) | 200 | UDMH (9.0) | 200 | UDMH (9.0) | 200 | 12 | 15439 | 到达中地球轨道 | ||
3 | UMDH (25.3) | 285 | UDM (25.3) | 285 | UDMH (25.3) | 285 | 10 | 30981 | 到达中地球轨道 | ||
4 | UMDH (25.3) | 285 | UDM (25.3) | 285 | UDMH (25.3) | 285 | 12 | 30921 | 到达中地球轨道 | ||
5 | 煤油(25.0) | 325 | UDMH (9.0) | 200 | UDMH (9.0) | 200 | 12 | 20684 | 到达中地球轨道 | ||
6 | 煤油(25.0) | 325 | 煤油(25.0) | 325 | UDMH (9.0) | 200 | 12 | 27113 | 到达中地球轨道 | ||
7 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 8 | 0 | 无法起飞 | ||
8 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 煤油(6.5) | 200 | 12 | 15382 | 到达中地球轨道 | ||
9 | 煤油(25.0) | 325 | 煤油(25.0) | 325 | 煤油(25.0) | 325 | 8 | 40194 | 到达登月轨道 | ||
10 | 煤油(25.0) | 325 | 煤油(25.0) | 325 | 煤油(25.0) | 325 | 12 | 40099 | 到达登月轨道 |
UDMH表示偏二甲肼
图5
4 结语
在充有一定氮气的氧弹中,采用环境恒温量热计,对蔗糖在硝酸钾为氧化剂、三氧化二铁为催化剂条件下的燃烧过程进行研究,记录蔗糖燃烧过程的温度变化,通过雷诺作图法校正ΔT的测量偏差,最终计算得到蔗糖燃烧过程的恒容热效应。研究了硝糖燃料配比对其爆热值和比冲的影响,结果表明,在蔗糖:硝酸钾:三氧化二铁的质量比为39 : 59 : 2时,其爆热值最大。
进而在火箭燃料的计算机模拟系统,通过改变燃烧室压力和燃料燃速,得到硝糖燃料的比冲、火箭飞行高度、单级火箭将东方红一号卫星送入预定轨道的发动机参数。还在该计算机模拟系统,对二级火箭和三级火箭中燃料与火箭飞行状态的关系进行研究,探究装有对应火箭燃料的飞行器能否把相关的航天飞行器送入预定轨道。
5 创新性/特点/特色声明
1)氧弹量热计充氮测定硝糖燃料的爆热值,筛选出硝糖燃料最佳配比;
2)将实验理论数据与人工智能结合,实现计算机对火箭发射仿真模拟;
3)融合航天技术的创新性实验改进,启发学生创造性和对科研的兴趣。
参考文献
DOI:10.1021/ed083p260 [本文引用: 1]
/
〈 |
|
〉 |
