高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

二阶滑模控制在十字梁控制实验系统中的应用

为了削弱传统滑模控制在十字梁实验系统中应用所产生的抖振问题,运用了高阶滑模的设计思想设计了多输入多输出十字梁试验系统的输出跟踪二阶滑模控制器,从而实现了高阶滑模在多输入多输出系统中的应用,并对该系统进行了仿真,仿真结果表明,二阶滑模控制器能够有效地削弱系统产生的抖振.      

重力波波包在时变背景场中传播的数值研究

采用二阶算子分裂格式对高斯型重力波波包在时变背景温场和时变背景风场中的非线性传播过程分别进行了数值模拟和分析.研究表明,时变背景场对波包空间结构演变影响不明显,各物理量的偏振关系与线性理论偏振关系较一致;在时变背景温场中传播,水平群速度和垂直群速度的变化趋势与线性理论结果相反,波频率的变化也与线性理论色散方程计算结果有一定的差别,这表明了线性重力波理论的局限性;在时变背景风场中传播,波频率的变化比时变背景温场中更显著.传播路径与在稳恒背景场中比较,更接近线性理论计算出的路径.      

一般状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型和高精度多介质Piecewise Parabolic Method(PPM)算法模拟一般形式Mie-Grüneisen状态方程的复杂多流体可压缩流动.根据多流体接触界面无振荡原则设计高精度计算格式,采用双波近似和两层迭代算法求解一般状态方程的Riemann问题,利用Strang分裂方法将一维算法直接推广到高维情况,不需要动能校正.二维多流体接触间断平移问题的模拟结果表明本算法在接触间断附近压力和速度没有振荡,界面数值耗散很小,被控制在两到三个网格之内.本文还模拟了二维Riemann问题和激波同多流体交界面的相互作用,算例结果表明本文方法可以有效地处理接触间断、激波等物理问题,且具有耗散小精度高的特点.     

一类单控制方向已知的随机不确定二阶系统反演控制研究

针对于一类单控制方向已知的给定区间随机不确定二阶线性系统进行了反演控制研究,从反演控制原理说明反演控制的有效性,并且通过与PID控制方法进行仿真对比,说明了针对标称模型难以确定的随机分布被控系统,反演控制方法对模型参数在给定区间完全随机不确定的情况下有更强的鲁棒性。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

基于五阶容积-二阶平滑变结构滤波算法的捷联惯导系统初始对准技术

为了改进捷联惯导系统大方位失准角初始对准情况下的方位失准角对准精度和稳定性问题,提出了一种组合对准方法。该方法将Kalman滤波对准模型与罗经对准模型相结合构建了一种组合对准机制,并在此模型基础上提出了基于非线性滤波和鲁棒滤波的五阶容积-二阶平滑变结构滤波算法,将五阶容积-二阶平滑变结构滤波算法与组合对准机制相结合,从而实现大方位失准角的初始对准。仿真结果表明,五阶容积-二阶平滑变结构滤波对准方法得到的方位失准角对准精度明显优于常见的几种对准方法,且其对准结果的重现性更高。因此,提出的五阶容积-二阶平滑变结构滤波算法能够很好地适用于大方位失准角的初始对准。     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

沿流向微结构沟槽流场直接数值模拟

采用高阶格式对覆有V型对称沟槽表面的槽道湍流流动进行了直接数值模拟,数值方法采用有限差分法。为精确求解沟槽壁面的湍流流动,对流项的离散采用7阶WENO（Weighted Essentially Non-Oscillatory）格式;时间推进采用分数步时间推进与低耗散、低色散Runge-Kutta方法（LDDRK方法）结合的格式;黏性项的离散采用6阶中心格式。模拟的雷诺数为5 000（基于槽道高度的1/2）,计算的沟槽宽度范围为13~44,沟槽斜壁与水平面夹角为60°。数值模拟结果表明,与平板相比,沿流向沟槽表面的阻力最高降低了9%。数据分析发现出现减阻效果时,沟槽减少了近壁面处顺流向涡的数目,并且减阻机理与微沟槽阻碍大尺度流向涡与沟槽壁面的直接碰撞,使沟槽表面湍流脉动得到削弱有关。