[发明专利]一种氧化亚氮储箱自增压过程预测方法

权利要求书

1.一种氧化亚氮储箱自增压过程预测方法，其特征在于：包括以下几个步骤：

步骤一、构建储箱物理模型：

将储箱内从上到下依次划分为气相区、饱和液相层区和过热液相区，各区内部温度均匀，饱和液相层区温度为恒处于氧化亚氮气体压力下的饱和温度，饱和液相层区的物理厚度视为零，进出饱和液相层区的物质量及能量均是守恒的，过热液相区与饱和液相区间换热为沸腾换热，沸腾换热率远大于“液-液”导热和“气-液”对流的换热率，气相区与饱和液相层区间换热为自然对流换热，沸腾换热率与自然对流换热率成正比，气相区、饱和液相层区和过热液相区的氧化亚氮液相压力与氧化亚氮气相压力相等，氧化亚氮液相及气相总焓仅与温度相关；氧化亚氮气相符合真实气体状态方程；储箱的内壁满足绝热边界条件；

步骤二、建立微分代数方程组：

按照整个储箱的气相区、饱和液相层区、过热液相区三部分控制容积，根据开口系统能量方程、质量守恒方程及步骤一中建立的储箱物理模型中的条件，获得包含九方程及九未知数的微分代数方程组：

质量守恒方程：mG·=mLV·-mprop·---(1)]]>

mL·=-mLV·---(2)]]>

其中为气相区的质量变化率，为液体蒸发速率，为氧化亚氮质量流率，为过热液相区的质量变化率，单位为kg/s；

对于气相区和过热液相区适用的能量方程Qin·=dEdt+Σme·he-Σmi·hi+W·---(3)]]>

其中为进入储箱中气相区或过热液相区的净热量，E为系统总能量，为流出储箱中气相区或过热液相区的氧化亚氮质量流率，为进入储箱中气相区或过热液相区的氧化亚氮质量流率，h_e为流出储箱中气相区或过热液相区的氧化亚氮比焓，h_i为流入储箱中气相区或过热液相区的氧化亚氮比焓，为储箱中气相区或过热液相区与外界交换净功率；

根据储箱物理模型获得的公式：TS=A2A1-lg(PG100000)-A3---(4)]]>

mL·=a′QLShLV---(5)]]>

P_L＝P_G                        (6)

其中T_S为饱和液相区的温度，P_G为气相区压力，P_L过热液相区压力，a′为经验因子，A₁、A₂、A₃和a′为常数，取值分别为A₁＝4.80716087，A₂＝967.819748，A₃＝19.6368887，a′＝2.1×10⁴；h_LV为过热液相区氧化亚氮汽化潜热，Q_LS表示过热液相区和饱和液相层区的热交换率；

气体状态方程：PG=(RTGρG44.02-B2ρG-B1(ρG44.02)2)×1.01325×105---(7)]]>

其中，T_G为气相区温度，ρ_G为气相区密度，R、B₁、B₂为常数，取值分别为R＝0.08206，B₁＝3.782，B₂＝0.04415；

补充方程：ρL=1000×C1C2-(1-TL/TC)n---(8)]]>

V_G+V_L＝V_total         (9)

其中ρ_L为过热液相区密度，T_L为过热液相区温度，T_C为氧化亚氮临界温度，V_G气相区的体积，V_L过热液相区的体积，V_total储箱的总体积；T_C、n、C₁、C₂为常数，取值分别为T_C＝309.57，n＝0.2882，C₁＝0.44927，C₂＝0.27244；

步骤三、求解上述的微分代数方程组，获得储箱内部各区域压力、温度、密度、质量等参数随时间变化情况。