具有迎角／过载限制器的俯仰增稳系统控制律的研究

本文以JJ-6飞机为例,研究了采用具有迎角/过载限制器的俯仰增稳系统对该机在大迎角下稳定性和操纵性的影响.并根据该机在大迎角下的偏离/失速运动机理,探讨系统中限制器边界值确定的原则。研究结果表明,该系统具有抗偏离/失速的能力,进而可以防止尾旋发生。但在大迎角下过大的横侧操纵,飞机仍能进入尾旋运动。     

二维间断有限元混合网格混合算法研究

针对二维流场进行结构／非结构混合网格划分,在靠近物面处生成结构化网格,在结构化网格以外和远场之间以非结构化网格填充。使用间断有限元方法（DGM）离散Euler方程。在混合网格混合算法中,结构化网格区域使用结构DGM,非结构化网格区域使用非结构DGM,保证结构／非结构对接面处的通量守恒,并加入moment限制器消除激波前后伪震荡。对NACA0012和RAE2822翼型跨音速无粘流动的数值模拟结果说明了混合算法的有效性以及在捕捉间断方面的能力。     

二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.      

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。     

轴对称短化喷管型面设计及流动分析

对最短长度喷管(MLN)设计方法进行了改进和优化,在喷管喉部扩张段实现了光滑过渡,采用特征线方法设计轴对称短化喷管型面。采用高阶、高分辨率WENO格式对设计的喷管进行流场数值模拟,结果表明:短长度喷管和喉部加圆弧过渡的短化喷管出口流场均匀,能够达到设计要求,经过圆弧过渡的喷管型面出口流场品质更好。     

高精度混合型格式的构造及数值模拟

在守恒变量形式的无波动无自由参数耗散(NND)格式的基础上,引入样条插值的思想,并采用具有一致二阶精度格式的限制器函数,构造了一个在光滑区具有四阶精度,在激波附近降为二阶精度的混合型格式。在计算中采用具有总变差减少(TVD)性质的四阶R unge-K u tta方法进行时间离散,通过对多个典型算例的计算,并与NND格式作了比较,分析结果表明:混合四阶格式引起的耗散和波动较小,并且能够高精度地分辨间断。     

一种基于间断Galerkin格式的新型限制器设计

使用高阶间断Galerkin格式求解守恒律方程组时,激波附近的Gibbs效应容易导致非物理解的产生。为抑制这一现象,必须构造合理的限制器对数值解进行处理。目前间断Galerkin格式中的限制器多源于有限体积法,在非结构网格上只对低阶导数项进行限制,对高阶导数项则很难给出普适判据。文章对间断Galerkin解进行广义Fourier展开,实现不同频域范围内的谱分解;在新的模板坐标系下描述各阶方向导数的变化规律;结合当前单元和相邻单元的信息,分层限制各阶方向导数,实现对非物理解的抑制。通过求解Euler方程,对二维Riemann问题、翼型周围的亚、跨声速流动问题、前台阶问题以及超燃冲压发动机内流场激波反射问题进行数值模拟,检验了新型限制器的可靠性以及向高阶格式推广的可行性。     

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。     

应用维数分裂方法推广MUSCL和WENO格式的若干问题

首先对有限差分法和有限体积法的差异进行了讨论,在已有文献观点的基础上补充了二者在边界条件处理及网格需求等方面存在差异的新论据。介绍了逐维推导的MUSCL和WENO格式计算控制体界面通量的过程,认为此类格式计算界面通量的方法直接应用于高斯积分型有限体积法不够严谨,从而得到了应用维数分裂方法构造的MUSCL格式和WENO格式不属于高斯积分型有限体积法的观点,“积分格式”这一定义更能准确反映这类格式的特点。此外,讨论了MUSCL格式和WENO格式在曲线坐标系下不能保证守恒的原因,并简单介绍了消除方法。