[发明专利]一种RP‑3航空煤油空化的数值模拟方法

权利要求书

1.一种RP-3航空煤油空化的数值模拟方法，其特征在于：具体步骤包括：

步骤一：提供RP-3航空煤油物理替代模型的配方，各成分的摩尔百分比为：

正十二烷：52～58％

癸烷：28～32％

辛烷：1～3％

甲基环己烷：1～5％

甲苯：8～12％

步骤二：采用美国国家标准与技术研究院(NIST)开发的Supertrapp物性分析软件获取RP-3航空煤油物理替代模型配方的各成分的各类物质属性；

步骤三：采用GERG-2004方程将各成分混合为一种新的虚拟物质，实现所述虚拟物质的属性尽可能逼近真实RP-3航空煤油的属性，用来代表RP-3航空煤油混合物，所述虚拟物质具有完整的物质属性数据，所述虚拟物质即为RP-3航空煤油物理替代模型，其基本结构形式可以表示为：

α(δ,τ,x)＝α⁰(ρ,T,x)+α^r(δ,τ,x) (1)

α0(ρ,T,x)=Σa=1Nxa(α0a0(ρ,T)+lnxa)---(2)]]>

αr(δ,τ,x)=Σa=1Nxaα0ar(δ,τ)+Σa=1N-1Σb=a+1NxaxbFabαabr(δ,τ)---(3)]]>

δ=ρρr(x)---(4)]]>

τ=Tr(x)T---(5)]]>

其中下标a,b分别表示两组分，多组分时则依次类推，下标0和r分别表示初始量和相对量，上标0和r分别表示理想项和余项，α⁰为Helmholtz自由能方程的理想项，α^r为余项，δ为对比密度，τ为对比温度，F_ab为调节因子，ρ和T分别为各组分所对应的密度和温度，N为组分的种类个数，ρ_r(x)和T_r(x)分别为混合物的密度和温度函数，x为各组分所占混合物的摩尔分数；利用上述方程组可以通过由无量纲Helmholtz自由能方程偏微分求解得到混合物的物质属性；如压强p的求解可由公式(6)和公式(7)得到：

p(δ,τ,x)ρRT=1+δαδr---(6)]]>

αδr=(∂αr∂δ)τ,x---(7)]]>

上式中，R为各组分所对应的气体常数；

步骤四：将RP-3航空煤油物理替代模型的物质属性导入商用CFD软件材料库或自编CFD程序，定义为一种新的流体介质，分别设置液相和汽相数据；

步骤五：采用三维造型软件进行流场区域的建模，采用网格划分软件进行网格划分，将网格文件导入商用CFD软件或自编CFD程序；

步骤六：对计算流域的边界条件和初始化条件进行设置；

步骤七：采用基于有限元的有限体积法对方程组进行离散，其中对流项采用高精度差分格式，其他项采用中心差分格式，对方程组的求解采用全隐式耦合技术；求解方程组包括：

连续性方程：

∂ρm∂t+∂(ρmuj)∂xj=0---(8)]]>

动量方程：

∂(ρmui)∂t+∂(ρmuiuj)∂xj=-∂p∂xi+∂∂xj[(μ+μt)(∂ui∂xj+∂uj∂xi-23∂ui∂xjδij)]---(9)]]>

能量方程：

∂∂t(ρmCpT)+∂∂xj(ρmujCpT)=∂∂xj[(μPrL+μtPrt)∂h∂xj]-{∂∂t[ρm(fvL)]+∂∂xj[ρmuj(fvL)]}---(10)]]>

质量守恒方程：

∂ρ1α1∂t+∂(ρ1α1uj)∂xj=m·++m·----(11)]]>

蒸发源项：

m·+=Fe3αnuc(1-αv)ρvRB23|Pv-P|ρl---(12)]]>

凝结源项：

m·-=Fc3αvρvRB23|Pv-P|ρl---(13)]]>

湍流动能k方程：

d(ρk)dt=Pt-ρϵ+∂∂xj[(μ+μtσk)∂k∂xj]---(14)]]>

湍流动能耗散率ε方程：

dϵdt=Cϵ1ϵkPt-Cϵ2ϵ2k+∂∂xj[(μ+μtσϵ)∂ϵ∂xj]---(15)]]>

湍流粘性：

μt=Cμρmk2ϵ---(16)]]>

其中，ρ_m＝ρ_lα_l+ρ_v(1-α_l)，u和p分别为混合相的密度、速度和压强，μ和μ_t分别为混合介质的层流和湍流动力粘度,f_v为氮蒸汽的质量分数，L为汽化潜热，Pr_L和Pr_t分别为层流和湍流的普朗特数，h为焓，α_l液相体积分数，能量方程中最后一项为能量源项；分别为凝结和蒸发源项，下标i和j分别代表坐标方向，下标m,l和v分别代表混合相，液相和汽相；采用标准k-ε两方程湍流模型实现了方程组的封闭；

求解器计算达到给定的精度要求后求解停止，即可获得RP-3航空煤油在给定的几何和边界条件下的空化流动情况；以给定温度下流体介质液相和汽相的物质属性作为数据支撑的条件下实现RP-3航空煤油空化的数值模拟。