基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式，具有十分优异的分辨率特性。但在求解多维方程组时，HCW-R格式需要求解块状三对角方程组，因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P，构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式，其可沿迎风方向、由边界处向内推进求解，避免了处理三对角或块状三对角方程组，从而其计算代价与显式格式无异。虽然就分辨率而言，HCW-E稍逊于HCW-R，但前者的计算效率要显著高于后者。因此，当花费相同的计算代价，HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

SAS与PANS模型在圆柱绕流中的应用比较

分别运用SST(Shear Stress Transport)、SST-SAS(Scale-Adaptive Simulation)、两种变fk(模化湍动能的比例)函数的SST-PANS(Partially Averaged Navier-Stokes)湍流模型对Re=3900的圆柱绕流进行了数值研究,重点从湍流结构捕捉、气动力计算、涡黏性控制等方面,比较了SAS与PANS两类RANS/LES混合模型的计算能力,并通过不同网格计算分析了模型的网格敏感性。数值结果表明:SAS及两种变fk方法的PANS模型均具有求解小尺度涡运动的能力,并能较好地反映出绕流尾迹的三维非定常特性,同时PANS模型能捕捉到更多的非定常结构;SAS模型中自适应尺度Lvk立足于当地流动,对网格依赖较小,计算的湍动粘度分布更合理,能够更好地计算剪切层及回流区;两种PANS模型网格独立性较差,出现了雷诺应力不足的现象;类DES可变fk函数构造相对简单,所得fk分布更准确,使用tanh函数计算的尾迹区fk值偏低,对流场调控能力稍差。     

多梁复合结构的动力学建模及全局解析振型分析

多梁复合结构广泛用于工程领域，其动力学特性对相关装备的工作状况有着重要影响，建立准确适用的动力学模型对后续的动力学响应参数影响规律研究有决定性作用。针对工程中常见的多梁复合结构，给出了位移场的一般表达式，利用哈密顿原理得出多梁复合结构的运动控制方程和边界条件，对其全局解析振型进行理论计算，给出了计算多梁复合结构数值振型的一般过程，基于数值分析和实验分析以L型梁和Z型梁为例利用频率相对误差、振型定性比较和MAC定量分析说明了全局解析振型的求解验证方法，利用伽辽金方法和全局解析振型给出了截断原系统得到离散系统运动方程的一般形式。提出的多梁复合结构的建模方法和全局解析振型求解验证方法，对工程中此类结构动力学响应的参数影响规律研究有重要意义。     

圆弧激波缓聚点火及后续燃烧传播

采用激波管实验和准一维数值模拟的方法，对预混可燃气体中圆弧汇聚激波的自点火现象及后续燃烧波的传播特性进行研究。其中圆弧汇聚激波由平面运动激波通过精确设计的弧形过渡管段转变得到。研究表明:收缩段中圆弧汇聚激波波后的非均匀梯度环境由激波在平直段、弧形过渡段和扇形收缩段中传播所分别诱导的3个梯度区共同构成。随着圆弧汇聚激波的不断增强，圆弧激波后某处首先形成一个无激波的温和反应区。该反应区逆流锋面的初期运动速度远超Chapman-Jouguet（CJ）爆轰波速，而反应产物区流动则呈现出一定弱爆轰波特征。进一步分析发现，该反应锋面本质上是一种"自发反应波"（spontaneous reaction wave），而非常规意义上的动力学波，其速度与汇聚激波波后气流点火时间梯度的倒数吻合。而后，反应区的扩张速度很快降至CJ爆轰波速以下，伴随反应锋面附近激波的产生以及激波-火焰复合结构的形成。激波-火焰结构最终加速演变为反向传播的爆轰波。在一定的条件下，由于入射激波转变过程和汇聚所构造的特定点火环境，自发反应波可再次赶超爆轰波，成为新的燃烧波前；而当自发反应波速度再次低于CJ爆轰波速时，它将再次转变为爆轰波；在此过程中，原先的爆轰波阵面蜕变为反应产物中传播的激波。     

类IXV飞行器初期再入制导与姿态控制方法研究

摘要： 针对类IXV飞行器无翼式升力体构型及采用RCS/尾襟翼组合控制特点，研究其初期再入段的高精度制导与姿态控制方法.设计带过程约束的数值预测 校正制导律以提高制导系统的鲁棒性及精确性；根据其执行机构配置特点，设计基于RCS的航向控制律及基于RCS/气动舵的纵、横向复合控制律，并采用鲁棒伺服LQR技术进行控制参数快速设计.通过蒙特卡洛打靶仿真来验证算法的精度及鲁棒性.     

多种工质下螺旋波等离子体源波场结构数值模拟

针对螺旋波等离子体放电机理，开展了多种工质条件下的螺旋波放电等离子体内波场结构数值模拟研究。计算发现:氦气等离子体的Er分量在径向边界处的峰值更为突出，有利于等离子体在径向的输运，波电场径向分量决定了电流密度径向分量在内部的表现。在0.266 Pa和1.064 Pa两种气体压强条件下，通过波场结构验证了气压对于波阻尼影响的结论。波场结构是螺旋波在等离子体内传播以及能量沉积的微观体现，研究螺旋波波场结构是揭示其高电离效率的重要途径。初步探索了功率耦合机制，为实验系统优化及实验方案设定奠定理论基础。     

运载火箭推进剂复杂流动传热问题数值模拟中的模型简化方法

在运载火箭改进优化或新型运载火箭设计过程中,会遇到各种复杂的流动与传热问题。按照真实产品和物理过程直接建立的数值计算模型过于复杂,不同尺度的结构、流动与传热、长时间的飞行过程相互耦合,出现计算时间太长(数年)、收敛困难、程序调试困难等问题。针对初始条件复杂、飞行过程复杂、边界条件复杂3类具体问题,总结了3种相应的模型简化方法,分别是工程经验法、极限参数法、低维近似法,为推进剂复杂流动与传热的数值模拟提供参考。     

背景波系下的隔离段激波串运动特性及其流动机理研究进展

对高超声速进气道-隔离段激波串在复杂背景波系下的突跳运动特性及其流动机理的最新研究进展进行了综述，涵盖了背景波系作用下的激波串运动特性、突跳机理和突跳运动特性的数学描述方法，以期对高超声速进气道相关研究工作提供一定的参考。首先，对固定背景波系和变化背景波系下的激波串运动特性和突跳机制进行阐述，指出隔离段壁面压力顺压力梯度和逆压力梯度的交替变化是激波串突跳特性产生的内在物理机制。其次，对背景流场下隔离段激波串突跳运动的触发机理和触发条件进行了讨论。最后，基于对运动特性和突跳机制的认识，尝试给出了背景波系作用下的隔离段激波串运动特性的数学模型，为激波串前缘位置控制提供参考。     

预测悬停旋翼跨音速流动的一种组合自由尾流/CFD方法

一种新的组合自由尾流 / CFD方法用于悬停旋翼流动的 CFD解中以考虑实际尾流的作用 .用文中描述的尾流分析方法研究了螺旋尖涡的运动 .首先从广义尾流模型开始 ,用半经验公式模化了涡核对旋翼尾流的作用 ;然后在环量收敛和尾迹收敛的条件下完成了自由尾迹计算 ;最后应用 Jameson有限体积龙格 -库塔推进格式求解了欧拉方程 .所得结果与相关文献和实验数据进行了比较