用加权余量法建立三维涡流场的有限元方程

本文应用加权余量法推导了三维涡流场的有限元离散化方程,使得三维涡流场的有限元方程表达方便,概念清楚,计算较为简捷。     

利用配重改善平尾颤振特性方法的研究

针对某型飞机平尾颤振速度偏低的实际问题，详细探讨了利用配重提高颤振速度的方案。分析方法及所得结果将有助于改善平尾的颤振特性，从而提高颤振速度。     

防雹火箭整体尾翼结构研究

提出防雹火箭尾翼整体注塑成型新技术。介绍整体尾翼的结构设计和应用效果检验。这种新的整体尾翼成型技术，可以大幅度降低尾翼制造成本，在其它民用火箭及小型野战火箭研制中具有推广价值。     

平直槽道内湍流运动的大涡模拟

为在中型机上实现对湍流运动的大涡模拟,采用Moin的基本模型,在数值计算中采用了诸如高斯滤波函数、速度压力交错网格、谱方法以及一种简单的压力计算的修正等方法和技巧,计算结果与Moin结果比较,表明在我国现有条件下也有可能进行大涡模拟方面的工作。     

涡发生器控制平板边界层分离的大涡模拟

为了研究涡发生器(VGs)间距λ对控制边界层分离效果的影响,选取了4种涡发生器间距,λ/H(H为涡发生器高度)分别为5,7,9,11.采用大涡模拟(LES)方法对带逆压梯度的平板边界层分离流动及VGs控制分离流动进行了数值模拟.分析了有无VGs控制时,湍流场中大尺度相干结构及其演化规律,分别从旋涡间距、边界层内流体动能、压差损失等方面考察了VGs间距对控制流动分离效果的影响.研究结果表明当λ/H为5时,VGs间距过小抑制了旋涡的展向发展,λ/H为9,11时,VGs间距过大边界层内流体动能偏低,当间距λ/H为7时流动控制效果更优,此时计算域压差损失最小,相比较无VGs控制时,压差损失降低了30.95%.      

Transonic wing stall of a blended flying wing common research model based on DDES method

Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation(DDES) and unsteady Reynolds-averaged Navier-Stokes(URANS) equations methods based on the shear stress transport(SST) turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10~6. A joint time step/grid density study is performed based on power spectrum density(PSD) analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x_(cr)= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.     

基于被动控制的D型钝体减阻大涡模拟研究

钝体减阻在航空航天等多个领域内具有重要的应用潜力,也是重要的基础研究热点问题。为了精确捕捉流动控制的流动细节,并发展效果优良的流动控制减阻方案,论文围绕抽象出的D型钝体,采用大涡模拟方法,开展了被动控制减阻高精度数值模拟研究。首先基于前期研究成果,对D型钝体尾迹区剪切层附近放置一个光滑小圆柱的被动减阻方法开展了数值模拟,发现总阻力减小17.7%,与试验结果吻合很好,同时数值预测的速度场分布也与试验结果吻合良好。在此基础上,进一步提出了采用齿槽型表面结构的小圆柱对D型钝体尾迹区进行扰动,并开展了数值验证,发现总阻力减小21.4%,优于前期的减阻方法。最后研究了在三种雷诺数工况下两种小圆柱扰动情况下的减阻效果,均表现出良好的减阻效果,两种小圆柱扰动下总阻力最大降幅分别为19.6%和23.1%,同时基于大涡模拟计算结果对减阻流动机理进行了探讨。上述研究结果表明,通过进一步优化流场结构,可以得到更优的流动减阻方案。     

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

气体状态方程对正十二烷射流燃烧的影响

选取正十二烷作为航空煤油的替代燃料,应用大涡数值模拟方法和详细化学反应动力学相结合的方法对该燃料在超临界环境下的射流燃烧进行计算,研究其射流、燃烧及排放特性。重点考察了气体状态方程对射流形态、贯穿距、温度、混合分数等的影响。结果表明:理想气体和真实气体状态方程对仿真得到的射流贯穿距和宏观形态影响不大,而真实气体状态方程可以更好地反映在超临界环境下的液体燃料射流内部特征,如混合分数、温度;燃烧方面,仿真的点火延迟与实验存在10%的误差,火焰浮起长度、火焰形态、碳烟等均得到很好的预测,NOx与OH分布与温度场具有高度一致性。     

旋翼干扰下的直升机尾桨气动噪声计算研究

针对旋翼干扰下的直升机尾桨气动噪声计算问题,将CFD方法与Farassat 1A(F1A)公式相结合,提出了一种适用于旋翼干扰下的尾桨气动噪声数值计算方法。首先,应用所提出的方法进行了噪声的算例验证,并与参考文献结果进行了对比;然后,着重针对悬停和前飞状态旋翼干扰下的尾桨噪声及气动特性进行了计算分析,并与孤立尾桨状态进行了对比。计算分析结果表明:悬停状态下,旋翼干扰下的尾桨气动噪声水平显著大于孤立尾桨;而在前飞状态,由于旋翼尾迹对尾桨桨盘平面产生更强的气动干扰,这一现象会更加明显;同时,在旋翼干扰作用下,尾桨噪声的主传播方向也会发生明显变化。