拔刀斋作者
15年1个月前IP:未同步
342910
9楼
公式推导过程与Mathematica计算结果
按《“气态氧-汽油”火箭发动机设计与制作》
“以往的小型水冷式火箭发动机设计经验已经表明,我们可以预计铜燃烧室和喷嘴平均传热率Q约为3 Btu/(in2*sec)。”
1 btu = 1 055.05585 joules
1 in2 = 6.4516 cm2
Q=490.6J/(cm2*sec)=490.6W/cm2=4.906*10^6W/m2=4.906MW/m2

燃烧室壁厚1cm,材料:40Cr

用偏微分方程计算:热传导方程
x处微元(x-0.5dx,x+dx)温度为T(x)
流入x处厚度为dx的微元的热流率=(dT/dx(x-0.5dx)-dT/dx(x+0.5dx))*导热系数=d2T/dx2*dx/导热系数
温升dT/dt=流入微元x的热量/单位长度比热/dx=(d2T/dx2*dx*导热系数)/单位长度比热/dx
dT/dt=(d2T/dx2)*(导热系数/单位长度比热)=(d2T/dx2)*(导热系数/比热/密度)

40Cr导热系数约为30W/m*k、比热为450J/kg*k、密度为7800kg/m^3
导热系数/比热/密度=8.547*10^-6(W/m/k/J*kg*k/kg*m^3)=8.547*10^-6(m^2/s)

微分方程:dT/dt=(d2T/dx2)*8.547*10^-6(m^2/s)

边界条件:给定温度梯度dT/dx=热流率/导热系数
内壁面x=0,dT/dx=-(4.906*10^6W/m2)/(30W/(m*k))=-1.633*10^5(k/m),外壁面x=0.01m,dT/dx=0

初始条件:室温25度

Mathematica计算:
为了保证边界条件与初始条件连续,内壁面dT/dx数值乘上(1 - Exp[-1000 t])作起点平滑过渡。

结果:
燃烧6秒后内壁温度达到钢材熔点1400度;外壁温度达到600度红热,40Cr或不锈钢接近失去强度。

参考资料:
《“气态氧-汽油”火箭发动机设计与制作》https://www.kechuang.org/t/37433
钢的比热为450J/(kg*k):XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
40Cr的导热系数约为30W/(m*k):XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
高温强度特性:XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
钢材熔点:XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX

heat1.png
+1
学术分
虎哥
2011/12/09
这个计算帮了大忙。
评论
文号 / 342910

万流景仰
名片发私信
学术分 15
总主题 213 帖总回复 6946 楼拥有证书:学者 机友 笔友
注册于 2006/08/19 15:42:00最后登录 2022/11/01 09:18:28
主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:0

个人简介

暂未填写
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的