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

涡桨飞机的机翼、短舱等部件在滑流作用下周围的流场特性与无滑流状态下截然不同。所以,应该在涡桨飞机的机翼气动设计过程中考虑螺旋桨滑流的影响,从而使得机翼在真实飞行时滑流作用下表现出更好的气动特性。采用基于雷诺平均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。     

小尺寸多级轴流压气机级间干扰现象试验

针对某小尺寸多级轴流压气机,利用Kulite动态压力传感器以及高频压力探针,测量不同转速下转子叶顶动态静压分布以及级间流场,并对测得的动态压力信号进行分析,研究多级压气机转/静干扰作用。结果表明:高转速下,首级转子叶顶测得的流场能够反映典型的超声叶栅流动特征,多级压气机存在的级间干扰作用只影响到首级转子出口流场;中等转速下,级间干扰作用非常显著,影响到首级转子进口,并随压气机节流级间扰动强度增大。对这些级间扰动的对应频率进行模态分析,从而获得了非定常级间扰动的传播速度和模态数。     

马赫数6柱-裙激波/边界层干扰直接模拟

高超声速激波/边界层干扰比超声速工况下具有更强的可压缩效应,再附之后会形成极高的局部力/热载荷,严重影响飞行器飞行安全。而激波/湍流边界层干扰区附近流动的三维特性使得流动更加复杂而难以预测。采用直接数值模拟对高超声速条件下的柱-裙激波/湍流边界层干扰进行了详细研究,特别是对Görtler涡结构对分离泡、物面压力和热流造成的展向差异开展了定性和定量分析。研究发现,干扰区附近的分离泡结构呈现出明显的三维效应,且Görtler涡展向分离位置截面的分离泡要明显小于再附位置,而这两个截面上分离泡的运动基本同步,没有明显的延迟或超前现象。物面压力和热流在展向出现显著的不均匀性,展向再附位置的平均压力和热流要比展向分离位置分别高13%和16.2%,脉动压力和热流比展向分离位置分别高28%和20%。截面流向速度特征正交分解结果显示两个位置上的能量都集中在剪切层附近,并且展向再附位置上低频模态占有更高的能量。在采用模态重构流场分析分离区面积发现,展向分离位置重构误差更小,而展向再附位置上的重构误差收敛更快。