一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

变循环发动机建模技术研究

在对变循环发动机总体结构进行分析的基础上,参考双轴涡扇发动机部件模型的建立方法,建立了变循环发动机部件级数学模型.建立了区分叶根特性和叶尖特性的风扇部件模型,单外涵和双外涵模式的核心驱动风扇级数学模型.根据变循环发动机的特点,建立了反映变几何部件变化的稳态和动态共同工作方程.数字仿真结果表明:所建的变循环发动机模型能够实现工作模式之间的相互转换,并且在低空低马赫数双外涵模式下表现出了涡扇发动机特性,即高推力与低耗油率,而在高空高马赫数单外涵模式下相比涡扇模式提供的推力更大、耗油率更低,符合变循环发动机特点,验证了建模方法的可行性.      

固体推进剂/衬层界面脱粘裂纹的三区域界面层模型

将固体推进剂/衬层粘结界面层按照剪切模量的变化分成三部分,由此建立了其界面脱粘裂纹的三区域界面层模型。应用Fourier积分变换和传递矩阵方法,推导并数值求解Cauchy奇异积分方程,得到了反平面剪切状态下裂纹问题的半解析解,并讨论了界面层模量和厚度对断裂参量的影响。结果表明,界面层模量随贮存时间的延长而降低时,应力强度因子逐渐减小,但裂纹张开位移变化不大;界面层厚度对应力强度因子的影响则不明显。     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程,控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型,空间离散采用中心有限体积格式,时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷.方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混...     

离心泵双密封环平衡系统泄漏特性预测方法研究

为了准确预测离心泵双密封环平衡系统的泄漏特性,通过分析该平衡系统的几何特征和液体流动特点,对其主要过流部件进行简化,将泵腔液体流动简化为由圆周剪切流和径向压差流叠加而成的轴对称二维粘性层流运动,同时将密封环间隙液体流动简化为两平行平板间的定常不可压缩粘性流体层流运动.然后,利用N-S方程分别建立了泵腔和密封环间隙的液体...     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

高速流非平行边界层稳定性研究

导出了可压缩流的抛物化稳定性方程（PSE）。针对高速流动,特别是超声速和高超声速流动的非平行边界层稳定性问题进行了研究。引入高效的边界层变换、全流场高精度的差分格式及预估校正迭代和推进求解法来求解PSE方程,使得PSE方法中至关重要的正规化条件得到了满足,确保了数值计算的稳定。采用高马赫数下对稳定性起支配作用的第二模式,研究了高速流边界层稳定性的演变和特征,分析了流动的非平行性、压缩性,以及壁面冷却等因素对流动稳定性的影响．所得结果与相关实验数据吻合较好。     

热完全气体的热力学特性及其N-S方程的求解

首先用五次多项式拟合给出了温度在50～3000K范围内的热完全空气的焓值与温度之间的函数关系式，导出了其它热力参数e，cp，cv和γ的表达式。接着提出了热完全空气总温、总压的计算方法，并将其计算结果与量热完全空气的结果进行了比较。最后，将该热完全空气模型用于N-S方程求解，对NAPA软件进行了改进，并用该软件时高马赫数钝体绕流流场和乘波体流场进行了计算，分析了气体真实效应时流场结构及参数的影响。结果表明，本文提出的总温、总压计算方法及时NAPA软件的改进是成功的，可较准确地模拟高超声速流动的主要特征。