一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

自身场磁等离子体推力器数值仿真流场特性分析

磁等离子体推力器以其推力大、比冲高等特点,成为未来深空探测、星际航行任务首选的电推力器类型,研究推力器内部等离子体流场特性,有利于解释推力器出现的物理现象。针对特定自身场磁等离子体推力器,建立磁流体模型,使用TVD Lax-Friedrich格式以及ADI方法对推力器内部及羽流进行数值求解,得到等离子体流场及等离子体参数分布情况,仿真结果显示在高电流工况下,等离子体羽流更加集中,轴向加速效果更加显著,但高电流模式下阴极温度较高,不利于阴极的使用寿命。     

高超声速变外形飞行器建模与有限时间控制

针对高超声速变外形飞行器变形带来的参数摄动大、变形过程建模难、外界干扰复杂等大不确定问题,研究了一类可变后掠飞行器建模与姿态控制问题,设计了一种有限时间控制方案。针对变外形飞行器建立了带有变形量的面向姿态控制的三自由度模型,该模型能够反映出变外形飞行器的内在影响。分析了变外形飞行器在典型状态下的气动特征,并给出了连续变形关键气动数据可行处理方案。针对可连续变形的飞行器设计了一套有限时间控制方案,并证明了系统稳定性。进一步考虑控制律中用到的指令微分项,设计了有限时间指令收敛滤波器。利用扩张状态观测器,估计不易测量状态和“综合扰动”。以考虑复杂干扰下的高超声速变外形飞行器为对象进行仿真,结果表明：所设计的控制方案可解决不同变形速率下、存在复合干扰的飞行器姿态控制问题。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。     

基于截尾概率-非概率混合模型的可靠性优化算法

针对工程中截尾概率变量与非概率变量同时存在的情况,给出一种新的截尾概率与非概率混合可靠性模型。在该混合可靠性模型基础上,按照可靠性指标(RIA)法给出双层嵌套可靠性优化模型,并采用改进搜索策略后的ST-Powell优化算法在外层搜索设计变量的最优值,内层采用能保证收敛的改进的有限步长迭代法求解混合可靠性指标。数值算例表明,改进搜索策略后的ST-Powell优化算法的全局寻优性得到显著提升;改进搜索策略后的ST-Powell优化算法与改进的有限步长迭代法相结合求解双层嵌套混合可靠性优化模型的正确性得到验证,且对于非线性程度较高的极限状态函数同样能够得到满足截尾概率与非概率混合可靠性模型指标要求的最优解,并对工程结构算例具有很好的适应性。     

铣槽结构液体火箭发动机推力室壳体热应力分析

介绍了液体火箭发动机推力室铣槽结构热应力的数值分析方法,通过建立液体火箭发动机推力室的流场燃烧和导热理论模型,运用有限体积法考虑液膜冷却计算出发动机工作时的燃气、燃烧室壳体和冷却工质的温度场,将得出的结果作为壳体热应力计算模型的边界条件进行热应力场有限元分析。内、外壁温度的计算数据与实验结果基本相符。     

一类高阶有理差分动力系统

讨论了高阶有理差分方程xn 1=(αxn-k 1)/(-1 xnxn-1…xn-k 1),for n=0,1,…,其中α<0,k∈{1,2,…},x-k 1,…,x0都是实数,并且使得x-k 1…x0≠1。给出了这类方程解的明确表达式,并对其解的稳定性和收敛性也做了进一步讨论。     

小波变换数值均匀化方法求解振荡系数的方程

将小波数值均匀化方法用于求解二维快速振荡系数的抛物型方程.小波数值均匀化方法基于小波的多分辨分析,在细尺度空间上得到原方程的半离散方程,然后利用小波映射得到粗尺度空间的数值均匀化方程,即可在粗尺度上求解原方程,有效地减小了计算量.计算结果表明:小波数值均匀化方法与精细剖分的有限差分法相比较,既大大地节省了计算时间又获得了较好的精度.     

利用偏振图像加权融合及CLAHE算法的水下成像方法

针对传统基于偏振差分原理的水下光学成像方法中目标退偏振特性差异引起的图像中局部反射光损失的问题,本文提出了基于偏振图像加权融合与限制对比度自适应直方图均衡（Contrast Limited Adaptive Histogram Equalization, CLAHE）算法的水下成像方法。一方面,将原始水下图像分解为偏振光强图像和非偏振光强图像,根据不同退偏振特性目标在两幅图像中的灰度值分布特点,设计相应的权重因子,对两幅图像进行加权融合。从而实现在压缩原始水下图像中散射光的同时,保留更多目标反射光,提升整体目标反射光在融合图像中所占的比例。另一方面,为了进一步提升融合图像的对比度,利用限制对比度自适应直方图均衡算法对融合图像进行处理。该算法能够在提升融合图像对比度的同时,有效避免图像噪声的放大。实验结果表明,相比于传统的偏振差分方法以及独立的直方图均衡化算法,本文提出的算法能够有效提升水下图像的清晰度和对比度。