基于线化稳定性理论的后掠翼边界层内横流稳定性研究

采用Mack方法求解忽略曲率和压缩性的三维线化稳定性方程（Orr-Sommerfeld方程）,针对无限展长后掠翼,计算分析边界层内不同波长横流驻波沿弦向的放大率,确定最不稳定波的波长。计算结果与参考文献中的实验结论和计算结果符合较好,表明该方法能够准确的预测出最不稳定波的波长,可以用来指导横流驻波主导下的后掠翼流动稳定性问题研究。     

一种三维弹性机翼气动优化新方法

基于控制理论方法和有限元分析软件,通过将网格生成、流场计算、有限元分析、粘性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,研究与发展了计及静气动弹性影响的三维机翼气动优化设计方法,其中,雷诺平均Navier-Stokes方程为主控方程,气动载荷和结构静弹性变形量由气动/结构方程的耦合迭代求解得到,目标函数梯度信息由共轭方程数值求解得到。典型大展弦比机翼气动减阻设计结果表明：研究及发展的计及静气动弹性影响的三维机翼气动优化设计方法是有效的,能够有效考虑静气动弹性的影响。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

战损飞机动力学特性分析及杀伤研究

现代战争中,防空武器性能不断提升,对飞机造成的威胁越来越大,确定飞机受到打击时的杀伤状态及杀伤等级成为飞机设计时需要考虑的重要问题。飞机被防空武器击中时,可能发生机翼脱离损伤,其质量、面积产生明显变化,将导致飞机的重心、气动特性的变化,严重危害飞机的稳定性。针对上述情况,以某型作战飞机为例,建立损伤情况下的飞机模型,综合分析机翼脱离损伤对飞机重量、重心和机翼面积带来的影响;建立损伤飞机的气动模型,计算损伤飞机气动特性的变化情况。基于这些变化,分析损伤对飞机动力学特性的影响,建立损伤状态下的动力学方程,计算损伤飞机的失控时间,依此确定飞机的杀伤等级。本文所做研究可用于在飞机设计过程中分析其杀伤状态及杀伤等级,为飞机设计时飞机杀伤的分析评估提供一种有效的解法。     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

Modelling and observing Jovian electron propagation times in the inner heliosphere

The propagation of Jovian electrons in interplanetary space was modelled by solving the relevant transport equation numerically through the use of stochastic differential equations. This approach allows us to calculate, for the first time, the propagation time of Jovian electrons from the Jovian magnetosphere to Earth. Using observed quiet-time increases of electron intensities at Earth, we also derive values for this quantity. Comparing the modelled and observed propagation times we can gauge the magnitude of the transport parameters sufficiently to place a limit on the 6 MeV Jovian electron flux reaching Earth. We also investigate how the modelled propagation time, and corresponding Jovian electron flux, varies with the well-known ∼13 month periodicity in the magnetic connectivity of Earth and Jupiter. The results show that the Jovian electron intensity varies by a factor of ∼10 during this cycle of magnetic connectivity.     

基于迹向控制量的编队星群相位控制

为了得到控制量与控制时机对相位改变的直观公式,基于Hill方程以及迹向控制的约束条件,定性分析共面绕飞编队中控制量、控制时机与编队卫星相位之间的关系,分析结果表明:当需要改变的相位角小于90°时,可实现最省燃料的控制;当需要改变的相位角大于90°时,可进行趋于最小控制量的控制;当在相对运动椭圆上下点实施控制时,相位无变化.     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。