(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 201710146487.8
(22)申请日 2017.03.13
(71)申请人 北京环境特性研究所
地址 100854 北京市海淀区永定路50号
(72)发明人 段然 朱希娟 马静 吴杰
(74)专利代理机构 北京君恒知识产权代理事务
所(普通合伙) 11466
代理人 黄启行 张璐
(51)Int.Cl.
G06F 17/50(2006.01)
(54)发明名称一种固体火箭喷焰氧化铝固液相变的数值模拟方法(57)摘要公开了一种固体火箭喷焰氧化铝固液相变的数值模拟方法，包括：获取氧化铝在当前喷焰高度和喷焰温度下的固相比热值和液相比热值；确定糊状相氧化铝中固相氧化铝和液相氧化铝的含量；依据固相氧化铝和液相氧化铝的含量、以及氧化铝的固相比热值和液相比热值，确定喷焰中氧化铝的有效比热值；利用氧化铝的有效比热值模拟固体火箭喷焰中氧化铝的固液相变过程。本发明能够有效的模拟固体火箭喷焰中氧化铝颗粒固液相变过程所引起的热力学变化，提高喷焰流场的数值模拟精度和以喷焰为目标的红 外特性模拟的精度。权利要求书1页 说明书4页 附图3页CN 107092715 A 2017.08.25
C N 107092715
A
1.一种固体火箭喷焰氧化铝固液相变的数值模拟方法，其特征在于包括：
获取氧化铝在当前喷焰高度和喷焰温度下的固相比热值和液相比热值；
确定糊状相氧化铝中固相氧化铝和液相氧化铝的含量；
依据固相氧化铝和液相氧化铝的含量、以及氧化铝的固相比热值和液相比热值，确定喷焰中氧化铝的有效比热值；
利用氧化铝的有效比热值模拟固体火箭喷焰中氧化铝的固液相变过程。
2.如权利要求1所述的数值模拟方法，其特征在于，固体火箭喷焰中每个氧化铝颗粒的固液相变过程都时时处于平衡状态。
3.如权利要求2所述的数值模拟方法，其特征在于，固体火箭喷焰中，单个氧化铝颗粒表面上的温度梯度为零。
4.如权利要求3所述的数值模拟方法，其特征在于，固体火箭喷焰中，单个氧化铝颗粒通过热传导达到温度平衡的速度高于氧化铝颗粒与空气对流传热达到平衡的速度。
5.如权利要求4所述的数值模拟方法，其特征在于，糊状相氧化铝中液相氧化铝和固相
氧化铝的质量组成为：
式中，f l 为糊状相氧化铝中液相氧化铝的质量含量；T m 为糊状相氧化铝的温度中值，T为喷焰温度，ε
为糊状相氧化铝的温度中值与最低温度或者最高温度的差值的绝对值，三者的单位相同，(T m -ε,T m +ε)表示糊状相氧化铝的温度区间。
6.如权利要求5所述的数值模拟方法，其特征在于，固体火箭喷焰中，氧化铝颗粒在喷焰中流动和相变时不发生破碎或者重聚。
7.如权利要求6所述的数值模拟方法，其特征在于，当固体火箭喷焰中的氧化铝为糊状相时，
其有效比热值满足如下关系：
式中，Cp为喷焰中氧化铝的有效比热值；Cp l 为氧化铝在当前喷焰高度和喷焰温度下的固相比热值；Cp s 为氧化铝在当前喷焰高度和喷焰温度下的液相比热值；f l 为糊状相氧化铝中液相氧化铝的质量含量；L m 为单位质量氧化铝的潜热值；T为喷焰温度。
8.如权利要求7所述的数值模拟方法，其特征在于，单位质量氧化铝的潜热值为1.138×106J/kg。
9.如权利要求6所述的数值模拟方法，其特征在于，当固体火箭喷焰中的氧化铝全部为固相氧化铝时，喷
焰中氧化铝的有效比热值与氧化铝在当前喷焰高度和喷焰温度下的固相比热值相等；当固体火箭喷焰中的氧化铝全部为液相氧化铝时，喷焰中氧化铝的有效比热值与氧化铝在当前喷焰高度和喷焰温度下的液相比热值相等。
10.如权利要求1-9任一所述的数值模拟方法，其特征在于，所述比热为定压比热。
权 利 要 求 书1/1页CN 107092715 A
一种固体火箭喷焰氧化铝固液相变的数值模拟方法
技术领域
[0001]本发明涉及固体火箭尾喷焰流场及目标红外特性建模，尤其涉及一种固体火箭喷焰氧化铝固液相变的数值模拟方法。
背景技术
[0002]以下对本发明的相关技术背景进行说明，但这些说明并不一定构成本发明的现有技术。
[0003]固体火箭尾喷焰中含有三氧化二铝(Al2O3)颗粒，是由推进剂中的金属铝(Al)粉末氧化形成的，Al2
O3颗粒的形状类似球体，其粒径与发动机的尺寸和工作状态有关。由于尾喷焰中的Al2O3颗粒是重要的辐射源且有很强的散射作用，而且其存在会对气体流场产生较大的影响，所以在模拟固体火箭尾喷焰时必须考虑含有Al2O3颗粒的气固两相流动。
[0004]推进剂中的金属铝质量含量约为10％至20％，经过燃烧室内的氧化过程，金属铝转化为Al2O3颗粒的比例超过90％，因此可以认为尾喷焰中的颗粒成分为单一的Al2O3。在数值计算中，通常假设喷管和尾喷焰中的Al2O3完全为固态且不发生相变。然而实际上，Al2O3的熔点约为2327K，由于燃烧室出口的温度很高，Al2O3颗粒一般以液态的形式存在，因此其形状多为球形，其颗粒直径受到燃烧室的内弹道状态和喷管内加速过程以及发动机尺寸的影响。液态的Al2O3颗粒会在空气中发生凝固，由于火箭尾喷焰会在空气中发生二次燃烧，使得喷焰剪切层的温度大幅上升，当温度高于Al2O3的熔点，部分已经转化为固态的颗粒会再次融化形成液态，因此尾喷焰中事实上会同时存在颗粒液态转化为固态、固态转化为液态两个过程。固液相变的过程伴随着相变潜热的释放和吸收，且固、液两相的颗粒液具有不同的物理性质，因此相变会影响气体和颗粒的的温度，进而影响目标的红外特性。目前Al2O3颗粒相变的模拟方法集中于燃烧室部分，涉及了金属铝的氧化过程，然而在尾喷流段，对于存在的相变过程则多数予以了忽略，有关尾喷焰中Al2O3固液相变的数值模拟方法也是一个空白。
发明内容
[0005]本发明的目的在于提出一种固体火箭喷焰氧化铝固液相变的数值模拟方法，能够有效的模拟固体火箭喷焰中氧化铝颗粒固液相变过程所引起的热力学变化，提高喷焰流场的数值模拟精度和以喷焰为目标的红外特性模拟的精度。
[0006]本发明的固体火箭喷焰氧化铝固液相变的数值模拟方法，包括：
[0007]获取氧化铝在当前喷焰高度和喷焰温度下的固相比热值和液相比热值；[0008]确定糊状相氧化铝中固相氧化铝和液相氧化铝的含量；
[0009]依据固相氧化铝和液相氧化铝的含量、以及氧化铝的固相比热值和液相比热值，确定喷焰中氧化铝的有效比热值；
[0010]利用氧化铝的有效比热值模拟固体火箭喷焰中氧化铝的固液相变过程。[0011]优选地，固体火箭喷焰中每个氧化铝颗粒的固液相变过程都时时处于平衡状态。
[0012]优选地，固体火箭喷焰中，单个氧化铝颗粒表面上的温度梯度为零。
[0013]优选地，固体火箭喷焰中，单个氧化铝颗粒通过热传导达到温度平衡的速度高于氧化铝颗粒与空气对流传热达到平衡的速度。
[0014]优选地，糊状相氧化铝中液相氧化铝和固相氧化铝的质量组成为：
[0015]
[0016]式中，f l为糊状相氧化铝中液相氧化铝的质量含量；T m为糊状相氧化铝的温度中值，T为喷焰温度，ε为糊状相氧化铝的温度中值与最低温度或者最高温度的差值的绝对值，三者的单位相同，(T m-ε,T m+ε)表示糊状相氧化铝的温度区间。
[0017]优选地，固体火箭喷焰中，氧化铝颗粒在喷焰中流动和相变时不发生破碎或者重聚。
[0018]优选地，当固体火箭喷焰中的氧化铝为糊状相时，其有效比热值满足如下关系：
[0019]
[0020]式中，Cp为喷焰中氧化铝的有效比热值；Cp l为氧化铝在当前喷焰高度和喷焰温度下的固相比热
值；Cp s为氧化铝在当前喷焰高度和喷焰温度下的液相比热值；f l为糊状相氧化铝中液相氧化铝的质量含量；L m为单位质量氧化铝的潜热值；T为喷焰温度。
[0021]优选地，单位质量氧化铝的潜热值为1.138×106J/kg。
[0022]优选地，当固体火箭喷焰中的氧化铝全部为固相氧化铝时，喷焰中氧化铝的有效比热值与氧化铝在当前喷焰高度和喷焰温度下的固相比热值相等；当固体火箭喷焰中的氧化铝全部为液相氧化铝时，喷焰中氧化铝的有效比热值与氧化铝在当前喷焰高度和喷焰温度下的液相比热值相等。
[0023]优选地，所述比热为定压比热。
附图说明
[0024]通过以下参照附图而提供的具体实施方式部分，本发明的特征和优点将变得更加容易理解，在附图中：
[0025]图1是飞行高度15km的火箭在有无考虑Al2O3颗粒相变时喷焰中心线温度分布对比图；
[0026]图2是飞行高度15km的火箭在有无考虑Al2O3颗粒相变时喷焰高温区域的分布对比图；图2a是图2中a处的局部放大图；图2b是图2中b处的局部放大图；图2c是图2中c处的局部放大图；图2d是图2中d处的局部放大图。
具体实施方式
[0027]下面参照附图对本发明的示例性实施方式进行详细描述。对示例性实施方式的描述仅仅是出于示范目的，而绝不是对本发明及其应用或用法的限制。
[0028]本发明固体火箭喷焰氧化铝固液相变的数值模拟方法，包括：
[0029]获取氧化铝在当前喷焰高度和喷焰温度下的固相比热值和液相比热值；
[0030]确定糊状相氧化铝中固相氧化铝和液相氧化铝的含量；
[0031]依据固相氧化铝和液相氧化铝的含量、以及氧化铝的固相比热值和液相比热值，确定喷焰中氧化铝的有效比热值；
[0032]利用氧化铝的有效比热值模拟固体火箭喷焰中氧化铝的固液相变过程。[0033]固体火箭喷焰中的氧化铝同时存在颗粒液态转化为固态、固态转化为液态两个过程，固液相变的过程伴随着相变潜热的释放和吸收，因此相变会影响气体和颗粒的的温度，且固、液两相的颗粒液具有不同的物理性质。这些都会对目标红外特性模拟的准确性产生影响。本发明依据固相氧化铝和液相氧化铝的含量、以及氧化铝的固相比热值和液相比热值确定喷焰中氧化铝的有效比热值，并以氧化铝的有效比热值模拟固体火箭喷焰中氧化铝的固液相变过程，能够有效的模拟固体火箭喷焰中氧化铝颗粒固液相变过程所引起的热力学变
化，提高喷焰流场的数值模拟精度和以喷焰为目标的红外特性模拟的精度。
[0034]本发明中提及的比热，可以是定压比热、即单位质量的物质在压力不变的条件下温度升高或下降1℃或1K所吸收或放出的能量，也可以是定容比热、即单位质量的物质在容积(体积)不变的条件下温度升高或下降1℃或1K吸收或放出的内能，也可以是饱和状态比热、即单位质量的物质在某饱和状态时温度升高或下降1℃或1K所吸收或放出的热量，本发明对此不做具体限定。为了便于描述，以下以定压比热为例进行详细说明。
[0035]本发明中，可以假设固体火箭喷焰中每个氧化铝颗粒的固液相变过程都时时处于平衡状态，氧化铝颗粒的固相与液相的质量分配只与颗粒温度有关。如此可以在保证模拟准确性和精确性的基础上，大大简化数值模拟的数据处理。
[0036]由于喷焰中氧化铝颗粒的直径很小，而且具有很高的速度，所以即使喷焰中存在很高的温度梯度，仍然可以认为氧化铝颗粒表面的所有部位处于相同的气体环境。为了简化数值模拟的数据处理，在一些实施例中，固体火箭喷焰中，单个氧化铝颗粒表面上的温度梯度为零。进一步优选地，可以假定固体火箭喷焰中单个氧化铝颗粒通过热传导达到温度平衡的速度高于氧化铝颗粒与空气对流传热达到平衡的速度，因此可以忽略处于固液相变的颗粒的动力学过程。“高于”的程度可以根据数值模拟的精度要求进行设置，例如当数值模拟的精度要求较高时，可以要求单个氧化铝颗粒通过热传导达到温度平衡的速度是氧化铝颗粒与空气对流传热达到平衡的速度的10倍以上。
[0037]本发明中所述的糊状相，是指物质处于固液共存状态。优选地，糊状相氧化铝中液相氧化铝和固相氧化铝的质量组成可以表示为：
[0038]
[0039]式中，f l为糊状相氧化铝中液相氧化铝的质量含量；T m为糊状相氧化铝的温度中值，T为喷焰温度，ε为糊状相氧化铝的温度中值与最低温度或者最高温度的差值的绝对值，(T m-ε,T m+ε)表示糊状相氧化铝的温度区间。从上述公式可以看出，只要T m、T和ε的单位相同，三者的单位在计算过程中便可以抵消，因此T m、T和ε采用何种形式的单位对本发明数值模拟方法没有影响，本领域技术人员可以根据实际需要进行设置。
[0040]为了简化数值模拟的数据处理，可以假定：固体火箭喷焰中，氧化铝颗粒在喷焰中流动和相变时不发生破碎或者重聚，因此氧化铝颗粒的粒径、密度和比热值就只与温度有
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议