基于升华法实验研究后掠翼三维边界层的转捩

在西北工业大学的低湍流度风洞,采用升华法研究不同雷诺数下后掠翼上表面的转捩现象。实验发现雷诺数较低时,后掠翼上的转捩由流向不稳定触发,转捩位置在最小压力点之后,转捩分界为一条直线;当雷诺数足够大时,转捩由横流驻波不稳定触发,转捩提前到最小压力点之前,转捩分界呈现尖楔形状。该结果表明升华法不但能够较准确地分辨出后掠翼上的转捩位置,还能够区分不同的转捩机理,判断转捩是由流向不稳定还是横流驻波不稳定触发。此外,实验中还发现在横流驻波不稳定增长较大时,升华法能够提供转捩上游区域边界层内的横流不稳定信息;当横流驻波不稳定进一步增强时,该不稳定受模型表面粗糙度的影响较大,萘的喷涂有可能会影响到升华法的结果。     

弹性整流罩分离的流固耦合仿真方法

发展了一种流固耦合的动力学数值仿真方法,用以模拟弹性结构在其运动诱导的非定常气动载荷作用下的动态响应.针对给定的运动规律,对不可压流场采用侵入式边界方法求解了非定常气动载荷,其中未计及弹性结构变形对流场的影响.对非定常气动载荷作用下的薄壁弹性结构进行了瞬态动力学数值模拟.以地面试验状态的整流罩分离为例,对刚性体及是否考虑表面气体压力作用的柔性体等情形下的整流罩的动力学行为进行了对比分析,验证了发展的流固耦合动力学仿真方法的可行性.     

预处理特征型边界条件及其回流处理技术

通过特征型预处理Euler方程推导出预处理特征型的边界条件处理方法,并发展了预处理的回流处理技术,增强计算的鲁棒性和精度。通过特征型预处理Euler方程的特征关系得到边界条件处理方法;在亚声速速入口和出口边界,根据其特征方程确定回流边界处理方法。通过数值模拟发现:带回流的特征型边界处理方法比Turkel方法收敛快1.5倍以上,而计算精度与SIMPLE方法相当;在出口带有回流的计算中,其压力误差比Fluent降低68.3%。由此可知带有同流处理的特征型边界处理方法具有较高的鲁棒性和精度,在实际的工程计算中有一定的应用价值。     

折流燃烧室间接点火过程研究

为了研究某折流燃烧室的间接点火过程,对三种不同初始火核位置的燃烧室进行大涡模拟计算,并结合燃烧室点火试验和性能试验进行验证。研究结果表明：折流燃烧室的主燃区可分为三个区域,分别是主回流区、横跨气流区以及次回流区;燃烧室点火成功的关键在于初始火炬进入次回流区;通过大涡模拟计算得到的贫油点火边界油气比与点火试验结果的绝对误差在0.004以内。     

跨声速压气机转子叶尖泄漏涡非定常特性数值研究

为了研究跨声速单转子压气机系统叶尖泄漏涡的非定常特性,选取Rotor 67孤立转子为研究对象,针对不同背压工况与不同转速工况进行了非定常数值模拟。结果表明:在每一转速状态都存在一非定常边界,其将特性线分为定常部分和非定常部分。当转子运行在特性线非定常部分时,随着背压提高,叶尖泄漏涡脉动频率逐渐减小。这是由于背压提高使叶尖前缘负荷变小,叶尖泄漏涡的驱动力也变小。叶尖泄漏涡的频率特征与转子转速息息相关。高转速状态时叶尖泄漏涡主要表现出低频特征,低转速状态时叶尖泄漏涡主要表现出高频特征。这是由于转速不同,叶尖激波的脱体程度不同,激波对于叶尖泄漏涡的激励位置也不同。由非定常叶尖泄漏涡引起的压力波的周向传播速度在各转速下表现出较强的规律性。随着流量系数的减小,压力波波速呈线性减小趋势,且各转速下减小的速率大致相同。且在波速-流量系数曲线中,各转速的非定常状态起始点基本位于同一条直线。      

功能梯度材料热传导的多尺度边界元分析

功能梯度材料是物理性能连续变化的非均匀复合材料,其结构特点消除了材料内部严重的热应力界面,在航空、航天及核反应堆等工程领域具有广阔的应用前景。针对功能梯度材料的稳态热传导问题,给出一种基于边界元模型的统计多尺度分析方法,并针对典型问题进行数值模拟。首先,在功能梯度材料微结构随机分布且体积分数随着位置变化的表征基础上,建立等效热传导系数的多尺度边界元分析模型;然后,给出统计意义下的数值计算方法;最后,研究颗粒随机分布功能梯度材料的热传导性能,并通过与实验结果对比对算法进行验证。结果表明:统计多尺度边界元方法对于预测功能梯度材料的热传导性能是有效的。     

壁面粗糙度湍流扩展模型及流动数值模拟

采用有限体积法数值求解雷诺平均Navier-Stokes流动控制方程,数值模拟表面粗糙度对气动特性和热力学特性的影响。通过修正Knopp的基于切应力输运k-ω湍流模型的粗糙度湍流扩展模型中边界条件ω的值,提高了壁面附近的涡黏度,增加了总壁面切应力,从而引入壁面粗糙度的影响。通过对不同粗糙度的平板流动和NACA 65_2215翼型绕流进行数值模拟,分析粗糙度对表面摩阻系数、斯坦顿数和升力系数的影响。计算结果表明:粗糙度使得平板表面摩阻系数和斯坦顿数增大,翼型流动分离提前,升力系数减小。流经粗糙壁面的流动处于过渡状态时,改进的粗糙度湍流扩展模型改善了原模型速度偏移量与经验公式相比偏小的缺点,同时,表面摩阻系数、斯坦顿数和升力系数与实验结果吻合得更好。该改进模型不需要对壁面附近网格进行加密,能够快速、准确地模拟流经粗糙壁面的流动。     

基于孔边距约束和Shepard插值的孔位修正方法

在飞机装配中制孔的位置精度直接影响飞机的最终装配质量,自动化制孔以飞机数模为基准,受装配误差影响,在机翼翼盒等封闭结构中,难以确保制孔位置与骨架沟槽关键特征的边距满足工艺要求。提出一种基于孔边距约束和Shepard插值的自动化制孔位置评价及修正方法,通过激光跟踪仪及扫描仪获取实测骨架点云数据并进行去噪和配准,将骨架理论孔位映射到点云模型获得实际制孔位置,提取骨架沟槽关键特征的边界点并构造Ferguson曲线作为边界曲线,基于此计算制孔约束区域中实际孔位的边距及孔边距偏差。针对孔边距可能超差的情况,提出一种先基于孔边距偏差修正约束孔位再利用Shepard插值法修正其余孔位的孔位修正策略,在满足孔边距要求的同时保留了原孔位的行列分布特性。运用所提出的方法对某型飞机机翼骨架进行测量和数据处理,验证了该方法的有效性。     

Suppressing unsteady motion of shock wave by high-frequency plasma synthetic jet

Three Plasma Synthetic Jet Actuators (PSJA) under the high-frequency actuation are used to control the Shock Wave Boundary Layer Interaction (SWBLI), a high-speed schlieren image processing method based on spatial Fourier transform as well as snapshot proper orthogonal decomposition were used to study the control effect of high-frequency plasma synthetic jet on low-frequency unsteadiness of SWBLI. The analysis of the base flow shows that the separated shock wave actually has both large- and small-amplitude vibrations at low frequency. And the results revealed that the PSJA with an operating frequency of 2 kHz has the ability to reduce the energy of low-frequency component of shock wave motion, indicating that the 2 kHz actuation can effectively suppress low-frequency unsteadiness of the separated wave. Compared with the actuation frequency of 2 kHz, the energy of low-frequency component of the shock wave is enhanced under the 8 kHz actuation, which aggravates the low-frequency unsteady motion of the shock wave. It is likely that the actuation frequency is too high, thus the intensity of the precursor shock wave induced by PSJA becomes weaker. Additionally, as the 4 kHz actuation is applied, the pulsation of the separation region was enhanced, it is speculated that the actuation frequency is coupled with the oscillation frequency of the separation region.