广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况,开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质,应用合适的高阶精度格式对各项分别进行离散求解,并采用经典几何声学射线法计算声爆传播路径,实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟,验证了本文高阶离散方法的可靠性。进一步的结果分析表明,采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性,在同等网格数目上可获得更高分辨率的计算结果,且其网格收敛性远优于传统二阶精度格式,同时在计算效率上也有一定优势;在Burgers方程的各效应项中,非线性项的影响更为明显,因而采用高阶离散格式时其优势更为突出;热黏性吸收项对数值结果影响较小,采用高阶离散格式对计算精度的提升效果也不明显,实际计算中仍可以采用传统离散格式,甚至忽略该项的贡献。     

高效率的特征型紧致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格式的效率优势在求解高维问题时更为明显。      

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

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

基于离散多项式相位变换的SAR高阶相位误差估计

离散多项式相位变换（DPT）是分析恒定幅度多项式相位信号的有力工具，其主要用途是估计相位多项式的系数，具有估计精度高、算法简单等特点。在SAR的回波数据中相位误差往往可以用幅度恒定的多项式相位来表示。补偿这种类型的相位误差，常用的方法是多孔径地图漂移法、多孔径相位差分法和相位梯度自聚焦算法。在SAR自聚焦中，运用DPT算法估计二阶和高阶相位误差系数是可行的。本文详细介绍了DPT算法的原理，通过实测数据验证了DPT算法的性能，并与多孔径地图漂移法、多孔径相位差分法和相位梯度自聚焦算法进行了比较，结果表明在信噪比较高的条件下，DPT算法的估计精度远远高于多孔径地图漂移法、多孔径相位差分法和相位梯度自聚焦算法。     

自适应数字波束形成的格式滤波算法

从矢量空间分析的角度，分析了格式滤波算法在自适应波束形成中的应用，以空间正交化原理为基础，得到了一种新的递推型格式滤波算法，使格式滤波器的阶次降为原来的一半，从而较少地减少了算法的运算量以及简化了硬件结构。     

非结构网格高超声速MHD流场逆风格式数值模拟

应用逆风通量分裂格式在非结构混合网格上对二维高超声速理想磁流体绕钝头体的流场进行了数值模拟。控制方程为Euler方程耦合Maxwell方程的理想MHD方程,空间离散采用AUSM格式,时间推进采用显式5步Runge-Kutta格式。引用双曲型散度清除技术加强·Ｂ＝0的条件。计算模型为二维钝头体,在高超声速来流条件下,对不同磁感应强度的均匀分布磁场干扰下的流场进行了计算,得到了较满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法可用于高超声速MHD方程的求解。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

基于PIV速度场测量的压强梯度计算

介绍了基于PIV速度场测量结果计算压强梯度的基本原理和3种不同的计算方法,并通过对流高斯涡模型系统研究了离散格式选取、PIV参数设置和流场特征等对压强梯度计算的影响.结果表明:PIV速度测量的精度是压强梯度计算的关键性影响因素,1%噪声对结果的影响比其他过程高约2个数量级;时间分辨率和空间分辨率的设置并非越高越好,而是有个合适的中间值;拉格朗日方法在大多数情况下优于欧拉方法,但当流动结构较复杂时性能较差,应根据流动类型和流场特征合理选取压强梯度计算方法.      

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

守恒双曲型方程组的一类高效差分格式

本文将文［１］所构造的差分格式推广到方程组情形。数值实验结果表明,格式具有高效分辨激波的能力。