气相爆轰基元反应模型数值模拟

采用基元反应模型和高精度的ENO格式对气相爆轰进行数值研究，对H2／O2／Ar混合气体起爆和爆轰波传播过程的数值模拟结果表明，计算的爆轰波阵面参数和实验相当符合。数值研究的结果还表明，爆轰波反应区中参与反应的不同组分具有不同类型的变化特征。计算结果的精度随着网格尺寸增加而增加，并能保持较好的收敛性。     

高超声速流动分区对接网格算法研究

本文发展了一种能适应强间断流场的满足通量守恒的分区对接网格算法，并给出了保证通量守恒的通量插值系数应满足的条件。采用有限差分NND格式，求解层流NS方程，数值模拟了马赫数6—20绕圆柱的高超声速粘性流场。数值试验表明：方法具有良好的稳定性和收敛性，强激波能自由地穿越内边界。     

振动尖塔对风洞模拟大气湍流边界层的作用

针对大气边界层风洞模拟中常常出现的湍流度随高度衰减太快的问题，本文尝试了一种新的被动模拟方法——振动尖塔法。该方法主要采用了具有弹性底座的尖塔型旋涡发生器。尖塔群受风洞自由来流驱动，发生随机振动，向下游流场注入了强烈的低频扰动，致使试验区的湍流度、积分尺度和风谱惯性子区的宽度都明显增加，改善了风洞的模拟性能。     

由泄洪水流引起的水舌风现象的数值模拟分析

本文采用细网格非静力能量闭合边界层模式，对水电站大坝泄洪过程中水舌风的生成演变发展过程进行数值模拟。水电站开闸泄洪时，高速水流对气流瞬时的强大冲量作用在坝后形成一个风速较大区，随着高速水流趋向稳定，采用两种模拟方案对水舌风的发展演变过程进行进一步模拟。方案一，通过水气动量交换作用，水流对气流的拖曳力仅来自于大坝溢流面；方案二，拖曳力来自于溢流面及溢流面与消力池接口面。对模拟结果进行分析，发现不考虑水舌风效应时，坝体下风方向出现上下两个方向相反的微弱涡旋；当考虑水舌风作用时，高速水流的瞬时强大冲量作用在坝后形成一个风速较大的逆时针涡旋；水流趋向稳定后，两种方案进行模拟的结果与不考虑水舌风时相比，两个涡旋均得到加强，尤其是处于下面的溢流面处的逆时针涡旋。     

三元收缩段优化设计研究

本文以某低速风洞收缩段为例，对收缩段设计方案的数值计算结果进行了综合分析。文中给出了数值计算结果、分析过程以及典型算例，并给出了该风洞收缩段最佳的气动设计方案，分析了几种常用收缩曲线的特点，同时给出了低速风洞收缩段气动设计的优化过程和可行的方法。结果表明，本文提出的优化设计方法是风洞三元收缩段气动设计较好的有效的气动设计手段。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

复合材料层板低速冲击剩余强度的研究

 针对复合材料层板受低速冲击后的剩余压缩强度问题进行分析计算,把冲击破坏区看作一个含有随机分布裂纹的圆形不均匀体,采用有限元建模分析,结合冲击后层板的试验所得的载荷/ 位移关系,计算得到冲击破坏区的剩余模量。再采用有限元建模分析含圆形冲击损伤区的矩形复合材料层板,求解应力及最大位移,并依据最大应力破坏准则,预测复合材料层板的冲击后压缩强度,计算结果与试验数据的比较表明分析结果可靠。     

轻型高精度卫星的变结构姿态控制器

针对某些小卫星高指向精度和高稳定精度的姿态控制要求,设计了能克服反作用轮转速过零扰动的变结构姿态控制器,并对反作用轮转速过零时低速摩擦对卫星姿态产生扰动的机理进行了分析,建立了仿真用反作用轮低速摩擦动力学模型。同 PID控制器相比,该变结构姿态控制器能有效抑制反作用轮的低速摩擦影响,并具有良好的鲁棒性。数学仿真进一步证明了该变结构姿态控制器的有效性,其指向精度和稳定精度分别可达 0.3°和 0.0 0 1°/