三级压缩楔前缘半径对壁面静压及热流分布影响研究

为了研究三级压缩楔前缘半径对壁面静压和热流分布的影响,在0.6m激波风洞上开展了试验测量,模型长约0.6m前缘半径为0mm3mm试验名义马赫数为5.98。研究结果表明,试验得到的壁面静压和热流数据重复性很好。采用Fluent软件进行了二维和三维流场参数辅助模拟分析,三种不同湍流模型获得的壁面静压分布差别不大,均与试验结果吻合较好;不同湍流模型获得的壁面热流分布差异较大,采用标准k-ε模型得到的结果与试验吻合较好。试验和数值模拟结果均表明,在第二和第三压缩面上,经过激波后,壁面静压逐渐上升到一个压强平台;壁面热流逐渐上升到一个局部极大值,然后在同一压缩面持续下降。随前缘半径增加,壁面静压和壁面热流整体减小,压强平台值和热流局部极大值也减小,而达到压强平台和热流局部极大值需要的长度增加,显示激波边界层干扰影响区域增大。     

旋翼桨涡干扰噪声开环桨距主动控制研究

直升机在斜下降飞行时旋翼产生的桨涡干扰(BVI)噪声十分严重,桨距主动控制是降低旋翼BVI噪声的有效手段之一。为摸索其对旋翼BVI噪声的影响规律并阐释其机理,开展了开环桨距主动控制对旋翼BVI噪声的影响研究。建立能够计入开环桨距主动控制的旋翼自由尾迹模型,并结合翼型气动力模型及基于FW-H(Ffowcs Williams-Hawkings)方程的旋翼载荷噪声计算模型,建立旋翼BVI噪声开环主动控制模型。以40%缩比的4桨叶BO-105直升机模型旋翼为算例,在风洞配平状态下开展开环桨距主动控制对旋翼BVI噪声的影响研究。通过分析算例旋翼在不同相位、幅值的桨距主动控制下的BVI噪声声压级、桨盘气动载荷及桨盘迎角分布,总结出开环桨距主动控制影响旋翼BVI噪声的规律,并初步阐释了其机理:适当的桨距主动控制可改善桨盘迎角分布,降低桨涡干扰位置附近的桨叶气动载荷,从而降低BVI噪声。     

考虑螺旋桨滑流影响的机翼气动优化设计

涡桨飞机的机翼、短舱等部件在滑流作用下周围的流场特性与无滑流状态下截然不同。所以,应该在涡桨飞机的机翼气动设计过程中考虑螺旋桨滑流的影响,从而使得机翼在真实飞行时滑流作用下表现出更好的气动特性。采用基于雷诺平均Navier-Stokes方程的多重参考坐标系(MRF)方法对螺旋桨滑流进行高精度准定常数值模拟,通过自由变形(FFD)技术实现螺旋桨飞机机翼的参数化构建,应用径向基函数(RBF)插值的动网格技术进行网格自动生成,获得样本机翼在滑流影响下的气动数据后,建立目标函数和状态函数的Kriging代理模型,结合随机权重粒子群优化(PSO)算法、Kriging代理模型和对应的EI(Expected Improvement)函数加点准则进行加样本点以及代理模型重建,从而建立滑流影响下机翼气动优化设计系统。使用该系统对某型螺旋桨飞机进行了滑流影响下的优化设计,结果表明,优化后的构型机翼和短舱在巡航状态下减阻达3.98counts,升阻比提高了3.325%。因此,建立的考虑滑流影响下的机翼优化设计方法是可行的。     

高速子弹穿透充液油箱的数值模拟

高速子弹击穿充液油箱产生水锤效应,造成油箱全面且灾难性的破坏,因而对子弹穿透油箱的过程进行模拟分析具有重要价值。使用ANSYS/LS-DYNA有限元软件,模拟900 m/s的子弹穿透装有60%水的4340钢制油箱以及空油箱的过程,模拟再现了水锤的产生过程,得到了不同时刻的物理图像和各种瞬变现象,确定了水锤效应的成因以及对油箱的破坏方式。研究表明,油箱中水压峰值出现在子弹入射板内侧,油箱底板和子弹出射板变形更严重。     

Simulation of 3D parachute fluid–structure interaction based on nonlinear finite element method and preconditioning finite volume method

A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual(GMRES) method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor(HHT) time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC(Harten–Lax–van Leer contact) scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel(LUSGS) approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.     

非线性谐波法在双级对转压气机中的进一步校检

以双级对转压气机为研究对象,在与试验结果及常规非定常方法对比的基础上,对新型非定常模拟方法——非线性谐波法(NHM)进行了时域和频域中的校检分析,探索其相比常规非定常方法的适用性.研究结果表明:非线性谐波法在设计工况下能够较为可靠地模拟确定应力的大小及其分布规律,具有较高的建模精度.对于近失速状态,由于流场中的固有周期性被破坏,导致非线性谐波法在该工况下不再适用,其与常规非定常模拟结果存在明显差异.      

安溢活门动态流场数值模拟研究

为了研究安溢活门主阀开启过程中的动态流场特性,采用Fluent动网格法和用户自定义函数方法对其进行流固耦合数值模拟。通过Fluent 3D非定常显式耦合求解器以及虚拟挡板技术求解某时刻主阀腔内的气压力,并通过有限气容充放气过程方程确定背压腔充放气过程对膜片的气压力,得到膜片上的总和力。进而通过求解质量弹簧阻尼系统动力学方程确定主、副活阀的运动速度和位移。结果表明:活门开启过程中开度越大,阀道内流场越复杂,流动越不稳定;活门开度增大,阀座处流量线性增加,入口腔流量变化平稳,而出口腔流量波动较大;随着开度增加,活门不同位置和膜片上的气压力变化规律不同。     

一种局部动态数据交换方法在流固耦合分析中的应用

 准确有效的数据交换方法是流固耦合(FSI)计算的关键,为此,提出了一种局部动态的耦合界面数据交换方法,一方面避免了曲面拟合方法由于病态矩阵引起的数值不稳定性和整体插值计算量大的问题,另一方面克服了表面跟踪方法引起的表面不光顺、出现锯齿状的问题。为了验证该方法,对三维柱面和某临近空间螺旋桨及某高空气球进行了分析,并与常体积四面体(CVT)方法做了比较。结果表明该方法在数值稳定性和表面光顺性方面有很大优势,是一种适合流固耦合分析的界面数据交换方法。     

基于ALE流固耦合方法的刷式密封泄漏特性理论与实验研究

现有刷式密封流固耦合求解模型因刷丝变形接触导致网格畸变而难以计算。本文提出了基于ALE (Arbitrary Lagrange-Euler)流固耦合方法考虑刷丝变形接触的三维瞬态双向流固耦合求解模型,设计搭建了刷式密封泄漏特性和可观测刷丝运动轨迹的实验装置。实验研究了刷式密封刷丝运动状态对泄漏特性的影响,不同压比条件下的泄漏特性变化规律,在验证刷式密封流固耦合求解模型准确性基础上,数值分析了不同时刻刷式密封速度分布与刷丝运动状态,并研究了不同结构参数对泄漏特性的影响规律。研究结果表明:基于ALE流固耦合方法解决了传统刷式密封求解模型因刷丝变形接触导致网格畸变而难以计算的问题,可准确计算刷式密封流场特性和力学特性;刷式密封刷丝吹下效应会降低泄漏量,刷丝颤振现象主要发生在气流入口处的前排刷丝域,会引起泄漏量增大,在气流力作用下,刷丝束轴向厚度变小,有利于减小泄漏量;在本文研究条件下,刷丝固定端与自由端的中间部位最先发生变形,随着时间的变化,气流速度逐渐趋于稳定,最大速度出现在末排刷丝与后挡板接触位置,刷丝自由端变形最大;刷式密封泄漏量随径向密封间隙的增大而增大,与径向密封间隙0mm时相比,径向密封间隙0.2mm时泄漏量增加了51.6%～62.8%。随着刷丝排数的增大,泄漏量逐渐减小,当刷丝排数大于25排时,其对泄漏量影响较小。随着刷丝间隙、刷丝与后挡板之间的轴向间隙、后挡板保护高度的增大,泄漏量先逐渐增大后趋于平缓。     

Maxwell模型中keV能量的Xe粒子对铜表面的适应系数

使用热流传感器对氙离子推力器轴向距离为500、700 mm和900 mm,径向角度为0°~15°(推力器出口平面中心为圆心,推力器出口轴线为0°)范围内羽流热流密度的分布进行了实验研究,获得了热流随角度和半径变化的实验数据。采用PIC-DSMC (particle in cell direct simulation of Monte Carlo)算法在不同适应系数下对实验条件进行仿真分析,对比仿真结果和实验结果得到适应系数。结果表明:keV能量的Xe粒子对热流传感器表面(铜)的适应系数接近1。