固体推进剂能量计算方法.docx

固体推进剂能量计算方法一 固体推进剂能量计算原理1，基本假设在火箭发动机工作时，固体推进剂的化学潜能转换为燃气的动能，经历了推进剂燃 烧和燃烧产物膨胀两个过程发动机的实际工作过程是非常复杂的其复杂性在于：由 于存在热损失，难以保证燃烧过程是等压绝热的；燃烧产物在燃烧室内分布是不均匀的； 对于含铝、含镁、含硼推进剂或含有某些金属化合物的性能添加剂的推进剂，存在凝聚 相产物，这些凝相产物在喷管膨胀过程中导致两相流损失；喷管流动难以保证等熵条件 等等为了反映固体推进剂能量转换过程的本质，抓住主要矛盾，在进行其理论性能预 估时，进行了一些基本假设1) 在燃烧室中，推进剂的燃烧反应达到化学平衡，且燃烧过程为等压绝热过程， 即热力学中的等焓过程；而且燃烧产物的分布是均匀的2) 燃气为理想气体，凝相产物的体积忽略不计3) 喷管中燃气的流动过程为绝热可逆过程，即为等熵过程；燃气在喷管中的流动 为一维定常流，即在喷管的任一截面上，燃气的组成及各性能参数的分布是均匀的4) 不考虑凝聚相燃烧产物的两相流损失2， 基本方程(1) 质量守恒方程常见的固体推进剂是由 C、H、O、N、Cl、Al 等元素构成的某些化学物质的混合 物。

对于这样一个复杂的系统，假设固体推进剂的燃烧产物共有 n 种，而固体推进剂所含有的元素共l种对j元素的质量守恒方程可表达成：£ ax = b (j = 1,2,…,/) (1)ij i ji=1式中，a为混合物系中第i种产物含j种元素的原子摩尔数，它由i燃烧产物的分ij子式得到；x为单位质量燃烧产物中第i种产物的摩尔数；ib为单位质量推进剂中含第j种元素的原子摩尔数，它由推进剂的假想化学j式得到2) 能量守恒方程根据假设(1)，燃烧室内燃烧为等焓过程，则有H = H (2)pc式中，H为单位质量推进剂在初温T时的总焓(通常取T = 298K)；p 0 0H为单位质量推进剂燃烧产物在平衡火焰温度T下的总焓cc实际上，在燃烧室热力计算中，等焓方程式(2)是作为判据使用的即根据式(2)来确定平衡火焰温度，进而计算出在定温(T )和定压(燃烧室工作压强P )条件下单位质量 cc推进剂燃烧产物的平衡组成分布｛x .｝ C = 1,2,…,n)然后可求出平衡燃烧产物在T下ci c具有的总熵Sc根据假设(3)，喷管中燃气的流动为等熵过程，则有S 二 S (3)ce式中，S为单位质量推进剂燃烧产物在T下的总熵;ccS为单位质量推进剂燃烧产物在喷管出口处温度T下的总熵。

ee同样，在喷管热力计算中，等熵方程也是作为判据使用的，即根据式 (3)来确定燃烧产物在喷管出口处的平衡温度T，然后求出喷管出口处平衡燃烧产物的总焓H eesp最后求得固体推进剂的理论比冲以及其它能量参数(特征速度、定压暴热)：(4)(5)式中， R 为燃烧室中燃气的平均气体常数;c比热比函数r=k+4k-1(6)Q = H -£ x AH 298Kp c ci fii=1式中下标 c 表示燃烧室， f 表示生成焓3， 固体推进剂性能计算的一般过程Stepl求单位质量(1kg)推进剂的假定化学式b和推进剂总焓H ； jpStep2燃烧室热力计算，使用温度尝试法，根据等焓方程求解T*、｛x* ｝以及S等; c ci cStep3喷管热力计算，使用温度尝试法，根据等熵方程求解T*、｛x* ｝以及H等;e ei eStep4 推进剂能量参数计算，计算 I 、 C *、 Q 等sp p温度尝试法基本原理：由热力学知识可知，随着温度的增加，燃烧产物总焓的绝对值和总熵是单调上升的以燃烧室热力计算为例，首先假设一个燃烧温度 T ，进而求c出该温度和燃烧室工作压强下的燃烧产物平衡组分(x }(i = 1,2,…,n)，然后求出燃烧ci产物在燃烧室的总焓 H ，进c步用等焓方程式(2)进行比较，H |>Hcp，则可知T

此时，使T增加一个步长h，c c c c即T = T + h，然后重新计算该温度和P下的平衡组成和总焓H，再将H与推进剂 c c c c c总焓H比较，此计算过程可以构成一个迭代计算过程步长h可以变化，如取100K、 p10K、IK、0.1K等，根据所需的精度要求去逼近推进剂总焓H二 最小自由能法求解定温定压下产物的平衡组成经过上面的阐述，剩下的问题是如何求解定温定压下燃烧产物的平衡组成目前工 程上经常使用的是吉布斯最小自由能法最小自由能原理：对于一个混合物系，求该物系的产物平衡组成，就是求该系统的 自由能函数达到最小时的产物分布1，自由能函数为了使问题得到一定简化,只考虑无凝相的混合物系的化学平衡组成，即只考虑C、 H、 O、 N、 Cl 五种元素组成的体系下面简要推导混合物系自由能函数的表达式设一个混合物系共含有l中元素，由n种气态产物组成若该物系在给定温度和压强下，单位质量的混合物系中有x摩尔的第i种产物，则单位质量该物系的总自由能函i数为(7)g o + R T xln P + ln x - ln X-i 0 i gi=1 i=1式中，X =Hx，将(7)式除以RT进行无量纲化，又由吉布斯自由能的定义式 g i 0i=1可得-S0(8)(9)g 0 i—RT0式中，H0为第i种产物的摩尔焓，可由式(9)求得;iS0为第i种产物的摩尔熵，可由式(10)求得。

iH = R C T + d T2/2 + d T3 3 + d T4 4 + d T5 5 + d )i 0 i1 i 2 i 3 i 4 i 5 i 6S = R (d lnT + d T + d T2.. 2 + d T3 3 + d T4 4 + d )i 0 i1 i 2 i 3 i 4 i 5 i 7式中，R为普适气体常数；T为温度;d、d、…、d为第i种物质的温度系0 i1 i 2 i 7数将(8)式代入(7)式，又令(11)C 三 £ + In Pi RT0该物系的总自由能函数可写成G (x )=工 x(Cii i=1+ ln x - ln Xig(12)2，Lagrange 乘数法求平衡组分(1) 自由能函数的近似表示首先假设满足质量守恒方程式(1)的解为｛y ｝ (i = 1,2,…,n)，此时该混合物系的总i自由能函数为(13)G(y)=Hy(Ci i i gi=1式中Y =工ygi i=1令A = x - y，即A为体系平衡时第i种产物的摩尔数(真解)与假设值只差又令 i i i iA = X -Y，则A=LAg g ii=1假设A很小，取G(x)在x = y处Taylor级数展开的前三项Q(x)来近似表示 iG(x),利用多元函数的微分原理，整理可得(14)G(x)u Q(x)= G(y)+工(C + ln y -ln Y)A +1i i g i 2i =1(2) 线性方程组的建立现在面临的问题是求出混合物系的自由能函数式 (12)在满足质量守恒条件系达到 最小时的产物组成｛xj,这是一个多元函数的条件极值问题。

写出Lagrange函数，即L (x ) = G (x )+ 工九jj=1b -工 a x u Q (x)+工九 I j i=i j ‘丿jj=1(15)当反应达到平衡时，混合物系的自由能函数达到最小值，即该点处 Lagrange 函数的一阶导数应为零将Lagrange函数对x求导，有 iC + ln y — ln Y iix X Ll+ ―^ — —g — /y九a二0j ijj=1i对上式进行整理，可改写成Xx = #yi Y ig(17)L b 九=L g (y )j j ij =1 i=1另一方面，将式(18)代入元素质量守恒方程式(1)，有厶a a yik ij i\ 丿(18)整理上式，a y + ik i i=1jj=1 i=1—La g (y)=b (k =1,2,---,l) (19)ik ii=1同时令Ykj jki=1a - a - y ,ik ij iXW =戸—1，贝I」式(20)可写成g+ y工九ai j ijj=1将上式两边对i求和，进行整理可得(y) (k = 1,2,…,l)(20)b -W + LY 入=La -gk kj j ik ij =1 i=1联立式(19)和(21)构成一个(l +1)元一次线性方程组，解该方程组可以求出l个待定系数九,九，…,九和未知数W，然后利用式(18)可求得n个平衡组分｛x ｝。

1 2 l i(3) 负值修正推进剂的燃烧产物大致可分为含量较多的主要产物和含量甚微的次要产物一般情况下，假设的初值｛y ｝ (i = 1,2,…,n)往往与真值偏差较大，其结果是利用不合理初值i计算出来的结果｛x ｝部分出现负值，显然在下一次迭代前，需对负值进行修正i设本次计算用假设值为｛y｝，本次计算结果为｛x ｝，修正后的值为｛x '｝负值修i i i正的公式为x' = y +0 (x - y ) (21)i i i i式中0为修正系数负值修正的具体作法是：丄 § — yStep 1令° = »，式中§取很小(如§ = 10-15)0i x —yiiStep 2在心｝(i = 1,2,…,n)中选最小值° ，取° = ° ，利用式(22)对｛x ｝进行i min min i修正，这样可使全部 x 均修正为正值iStep 3检验GC')-G(y)< 0是否成立,若不成立，令©二0帥，重新对｛x ｝进i行修正，直至G (x') < G (y )为止Step 4令y = x'(i = 1,2,…,n)，进行下一步迭代计算iiStep 5精度检验要求|x - y.| <£ i = 1,2,…,n。

满足上式时的｛x ｝即为给定T、1 1 max 1P 下的平衡组成。