采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

涵道尾桨的CFD模拟与验证

采用CFD(计算流体力学)方法分析在悬停和侧飞时直升机涵道尾桨的流场与性能。在轴对称的圆柱坐标系中,用有限容积法和SIMPLE(压力耦合方程的半隐式法)算法求解定常不可压的湍流N-S方程。在分析中,旋转的螺桨被描绘成沿螺桨桨叶展向分布的、与本地流动参数相关的、时间平均的动量源项,通过在N-S方程中加入此动量源项来替代螺桨对流体的作用。计算方法还包括涵道壁的阶梯型近似,原始变量的交错网格,k-ε湍流模型和涵道壁上的壁面函数法等措施。涵道尾桨的流场分析和性能预测同实验测量数据的良好的一致性表明,这种CFD方法可以有效的分析涵道尾桨的具体设计问题。      

湍流模型对发动机排气系统流场和红外辐射特征计算的影响

计算分析了采用k-ω SST(shear-stress transport)模型、RNG(renormalization group)k-ε模型配合标准壁面函数(RNG k-ε+标准壁面函数)、RNG k-ε模型配合增强型壁面函数法(RNG k-ε+增强壁面函数)3种湍流模拟方法,对典型涡扇发动机轴对称收敛喷管排气系统流场、温度场、组分质量浓度场及红外辐射特征计算结果的影响.研究表明:对于这3种方法,排气系统腔体内壁面压力分布基本一致;但在温度场、组分质量浓度分布以及红外辐射特征方面的差别较大.RNG k-ε+增强壁面函数计算得到的内涵壁面与燃气流的温度、CO₂与H₂O的质量浓度和红外辐射强度值最大,RNG k-ε+标准壁面函数的计算结果其次,而k-ω SST 模型的计算结果最小.通过与试验结果对比,发现k-ω SST模型对排气系统红外辐射特征的模拟更加合理一些.      

高超声速流动壁面催化复合气动加热特性

针对高超声速流动壁面催化特性,计算了不同壁面催化复合系数条件下的球锥驻点热环境。引入了经验证的数值求解Navier Stokes方程的方法,在不同壁温500K～2500K的条件下分别分析了O 2和N 2气体在壁面处的催化复合气动加热特性,得到如下结论：（1） 原子复合放热将提高近壁面温度梯度,改变近壁面组分分布;原子复合放热一部分加热飞行器形成组分扩散热流,一部分加热近壁气体提高近壁温度梯度。（2） 在壁面催化复合系数较小时,原子复合放热主要转化为组分扩散加热,对于不同壁面温度,壁面催化复合系数α<0.1时,单一气体反应组分扩散热流小于总热流的20%。     

磁屏蔽霍尔推力器技术的发展与展望

磁屏蔽霍尔推力器技术是近年来霍尔推进领域最具影响的创新突破,对于拓展霍尔推力器的应用范围,提高推力器的寿命具有重要意义。介绍了磁屏蔽霍尔推力器的原理及优缺点,从磁屏蔽的提出与验证、不同功率量级霍尔推力器的磁屏蔽技术以及磁屏蔽霍尔推力器热设计、背景压力敏感性、振荡模式转换等方面介绍了磁屏蔽的研究现状,并对未来磁屏蔽霍尔推力器技术的发展进行了展望。     

竖直平板Blasius流层流边界层流动与传热耦合解

对考虑辐射传热、流体黏度随温度变化、有和无滑移边界条件下的可渗透竖直平板Blasius流层流边界层的无量纲速度场与温度场进行了深入研究.经相似变换将描述速度场与温度场耦合的偏微分方程组转换成非线性常微分方程组,用Runge-Kutta法对常微分方程组进行了数值求解.研究了无量纲参数对无量纲速度场及温度场的影响,着重分析了滑移边界条件下速度和温度随无量纲参数的变化规律.结果表明:吸入时边界层变薄,喷注时边界层加厚;对比于无滑移边界条件,滑移边界条件下速度、温度边界层变薄;随着变黏度参数a或喷注与吸入参数的增大,壁面摩擦因数、局部努塞尔数Nu增大,速度和温度边界层变薄;随着普朗特数Pr、热辐射参数R的增加,或毕渥数Bi,布林克曼数Br的减小, 温度边界层变薄.      

渗透壁面竖直平板Blasius流熵产特性

对渗透壁面竖直平板Blasius流层流边界层流动熵产及传热熵产与总熵产的比值(Bejan数)进行了深入研究,将耦合的动量与能量偏微分方程组通过相似变换转换为非线性常微分方程组,用Runge-Kutta法求解获得了无量纲流函数与无量纲温度的相似解,进而研究了无量纲参数对熵产及Bejan数的影响.结果表明:吸入时最大熵产位于壁面处,吸入速度越大最大熵产越大,熵产随相似变量的增加单调下降;喷注速度越大壁面处熵产越小;近壁面处熵产随变黏度参数增大而增大,远离壁面处随喷注/吸入参数或变黏度参数的增大而减小;远离壁面处Bejan数随喷注/吸入参数或变黏度参数的增大而增大,近壁面处与此相反;熵产和Bejan数随着辐射参数和Reynolds数增大而减小,随着Biot数增大而增大;Biot数在0~2范围内对传热熵产的影响很剧烈,Biot数较小时总熵产由流动熵产控制,Biot数较大时总熵产由传热熵产控制;熵产随滑移参数的减小而增大,Bejan数随滑移参数的增大而增大;Grashof数对熵产的影响较大,在较大的Grashof数条件下熵产随相似变量的变化较剧烈;滑移参数越大,流动熵产越小,Bejan数越大.      

MEMS 壁面剪应力传感器阵列水下标定实验研究

壁面剪应力的精确测量对于研究水下物体边界层流动、寻求有效的减阻增效措施至关重要。MEMS 壁面剪应力传感器的标定首先是最基本的静态标定,其决定了其测量的精度和数据的可信度。本文在分析已有标定方法的基础上研发1种新型水下壁面剪应力给定装置,并采用数值方法计算分析不同流速下的壁面剪应力给定条件,进而设计壁面剪应力传感器静态标定方案,开展了一种 MEMS 热膜式壁面剪应力传感器阵列的水下静态标定实验,获得了各传感单元的标定系数。     

浮动环涡动对浮环轴承油膜压力分布影响的研究

研究了考虑浮动环涡动时,内、外层挤压油膜力的分布情况。从广义不可压Reynolds方程出发,推导了浮环轴承内、外层油膜的瞬态控制方程,结合Hirs整体流动理论和Moody壁面摩擦系数方程,获取浮动环的运动方程。通过算例,较为深入地研究了浮动环涡动情况以及对油膜压力分布的影响：涡动半径随轴颈偏心率增大而略呈减小趋势,内层油膜其正压力区将沿周向平移。     

重叠网格洞面优化技术的改进与应用

针对传统割补法在壁面狭缝处网格重叠失败的问题,提出了以单元顶点状态进行洞面优化的vertex方法及以单元大小进行洞面优化的volume方法的两种改进方法。这两种方法将传统割补法中切割和填补两步搜索方法改为一步搜索方法进行重叠网格的洞面优化。通过数值试验证明:两种洞面优化方法都能很好地解决壁面狭缝网格重叠问题,网格重叠区域流场数值传递正确,计算结果与试验值吻合,优化效率高,收敛速度快。volume优化方法与vertex优化方法相比可以减少由网格大小不对等引起的数值误差,计算结果和计算精度均优于vertex优化方法。