用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

非合作体附着下的组合体航天器自适应控制

空间中非合作体附着航天器本体后,新组合体惯量参数未知、系统引入的动量变化未知,容易造成航天器姿态控制失稳.针对此问题,提出了一种基于非合作组合体整体惯量估计的自适应控制方案,实现不同情形附着下新组合体姿态的快速和高精度恢复.通过构造李雅普诺夫函数,对自适应控制方案稳定性进行了证明.搭建的某卫星被非合作组合体附着后的姿控仿真平台表明,在非合作体惯量参数未知、引入动量未知以及附着前系统状态不定等恶劣情况下,采用该算法的组合体姿控误差可收敛到一个极小的领域内,可以实现对非合作组合体系统快速而有效的控制.     

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

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

特定涡旋流畸变对跨声速压气机性能的影响

为研究S弯进气道出口和翼身融合体(HWB)飞机发动机进口两种特定涡旋流对压气机气动性能和稳定性的影响,将设计的新型叶片式旋流畸变发生器与跨声速压气机Stage 67进行联合数值仿真,获得进口旋流条件下的压气机特性线和流场分布,并与转子进口为均匀来流时进行对比分析。结果表明:在100%换算转速,S弯进气道产生的旋流使得压气机压比和效率在近峰值效率点分别下降0.06%、0.85%,稳定工作流量范围减小5.97%,压气机稳定裕度降低1.13%,压比和效率特性线均向左下方移动。HWB飞机产生的旋流使得压气机在近峰值效率点的压比增大6.61%,效率下降6.25%,稳定工作流量范围减小26.88%,压气机稳定裕度降低2.84%,压比特性线向右上方移动,效率特性线向右下方移动。     

高速自旋飞行体的气动加热特性

一种典型的弹箭类高速自旋飞行体为例,选取目标在超声速、跨声速、亚声速飞行状态下的弹道计算数据作为来流条件,结合滑移网格和多坐标系法,采用基于密度的耦合隐式算法、Roe-FDS(flux difference splitting)通量格式和SST(shear stress transfer) k -ω湍流模型,对飞行体的外流场进行了热流模拟研究。重点分析了目标在高速自旋飞行条件下的壁面压力、温度、气流密度、热流率、湍流动能等变化规律,并与不考虑飞行体高速自旋的流场进行了对比研究。研究结果表明:在高速自旋飞行条件下,飞行体表面的流线相互干扰,贴壁气流对飞行体产生的扰动更加剧烈,飞行体尾部的气流集聚效应明显,湍流发展与演化过程更加复杂,飞行体壁面的气流速度、压力、温度、热流率均高于无转速飞行情况,尤其在超声速飞行条件下的差异显著。     

白车身轻量化设计分析

针对某轿车减重的项目要求,利用有限元灵敏度优化分析技术,基于经试验验证后的有限元模态分析模型建立灵敏度优化分析模型,用MSC．Nastran求解得到灵敏度矩阵和设计变量的迭代历史过程。结合灵敏度分析对基础模型减重,并对减重后的模型进行刚度、碰撞等相关性能的验证,完成轻量化分析研究。     

基于遥感的汶川震区水体快速提取

＂5.12＂汶川8.0级地震后发生大量滑坡、泥石流,河流部分阻断,形成一系列堰塞湖,急需快速提取水体信息,为救援提供重要依据。本文以震后北川县城为例,利用LBV变换与NDVI（归一化植被指数）图像结合模型对研究区水体进行提取,对＂福卫二号＂影像进行LBV变换,确定B阈值,并提取NDVI指数,两者求交集得到水体影像。该模型理论依据充分可靠,计算机容易实现,具有较高的准确度,可以将水体与低密度覆盖的水植混合体区分开,提取结果符合实际情况,提取速度快,完全能满足灾害背景参数8h内提取的要求。     

小攻角高超声速钝锥边界层失稳特性

利用高阶紧致格式,采用直接数值模拟(DNS)和线性稳定性(LST)分析方法,对高超声速边界层的失稳机制和转捩特点进行了研究.通过对马赫数为6的2?攻角高超声速钝锥边界层的稳定性分析发现:小攻角高超声速钝锥边界层存在多枝不稳定模态;周向速度使钝锥的稳定性特征与不考虑周向速度时有本质的差别;转捩线在接近背风面处出现拐折现象是由失稳模态发生转换引起的.      

基于大幅值Lyapunov轨道的地月转移轨道设计研究

基于平面圆形限制性三体问题模型,利用与绕月轨道相切的大幅值Lyapunov周期轨道,提出了一种新的地月转移轨道设计方法。根据Poincaré截面与限制性三体问题动力学系统对称性计算得到的大幅值Lyapunov轨道,通过与绕月轨道拼接,将地月转移问题转化为地球到大幅值Lyapunov轨道的转移问题。为保证探测器能够从近地轨道（LEO）切向逃逸到达大幅值Lyapunov轨道,通过计算其稳定流形,采用最近点作为Poincaré截面的终止条件求解探测器的初始状态,并根据初始状态完成地月轨道的设计。仿真结果表明,该地月转移策略相比于Hohmann转移,在同样只需要两次速度增量的前提下,约节约100 m/s的速度增量,该研究为地月转移轨道的设计提供了一种新思路。     

从地球到日-火系平动点轨道的转移

平动点轨道特殊的空间位置及动力学特征,使其在深空探测中具有重要的应用。以日-火系平动点轨道(Lissajous与Halo轨道)任务为目标,结合平动点轨道的不变流形理论,研究了小推力转移问题。首先给出了圆型限制性三体动力学模型下平动点附近不变流形(稳定和不稳定流形)高阶分析解以及相应的计算实例。接着以流形分析解为基础,建立了初始小推力轨道优化模型,并利用改进的协作进化算法求解初始小推力轨道。最后将初始轨道离散,采用多点打靶法将最优控制问题转化为参数优化问题,并用序列二次规划方法(SQP)求解。仿真结果证明轨道设计方法的有效性。