基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。     

一类非线性微分方程组奇点稳定的条件

本文讨论了一类非线性电路微分方程组奇点稳定的条件,用线性近似的方法,讨论了方程组(1)中参数P*对奇点稳定性的影响,得到了奇点(P*/Ф,Ф)在P*相似文献   

改进的小波均匀化方法求解小周期系数抛物型方程

对传统的小波数值均匀化方法通过加入逐层修正的步骤进行了改进,并将改进后的方法用于求解具有小周期系数的抛物型方程。数值算例结果表明：在保持较高计算效率的同时,改进后的方法比改进前方法在精度方面有明显的提高。因此,改进后方法是小周期系数抛物型方程的一种高效率、高精度的数值求解方法。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

基于N-S方程的旋翼翼型优化设计方法

将N—S方程流场分析程序和序列二次规划结合起来，提出了一种工程实用的直升机旋翼翼型设计方法，以提高旋翼翼型在悬停、前飞和机动等多个设计点及在多种约束条件下的气动性能。由N—S方程计算得到的升力、阻力等气动参数构成目标函数，应用数值优化程序对其进行最优化。计算实例表明：本方法设计质量高，所需机时少，易于实施，有较大的工程应用价值。     

轴对称流中音速线附近的特征线方程

 从Ｏｓｗａｔｉｔｓｃｈ的无粘、跨音速、轴对称无旋流的基本方程出发,在细长体条件下,通过新变量将物理面上的运动方程变成形式简单的速度面方程。利用音速线附近条件对速度面方程进行积分,可得解析解；由此可按Ｃｏｕｒａｎｔ理论推导出音速线附近的特征线方程,结果表明,特征线有两族,均为２／３次幂曲线,其数学形式与平面流相同。     

Navier—Stokes（NS）方程组差分计算中的物理和网格尺度效应及NS方程组的简化

对于高Re数流动计算，在通常二阶精度NS差分格式和网格数条件下，存在某些粘性项落入修正微分方程截断误差项的问题。这类NS方程组计算实际是计算某种简化NS方程组，而且重复计算误差物理粘性项既浪费机时和内存，误差积累又会对数值解产生不可预测的影响，避免外述缺陷的办法一个提高NS差分格式的精度，另一个是丢掉可能落入截断误差项的物理粘性项，把NS方程组简化为广义NS方程组，广义NS计算避免了误差物理粘性项误差积累对数值解的不可知影响，又可节省内存和机时，对高Re数流体工程计算很有好处。利用广义NS方程组计算超声速绕前向和后向台阶流动的结果表明：广义NS方程组与NS方程组的数值结果很好相符符。     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

基于CFD的电磁散射数值模拟

计算流体力学(CFD)技术中的时域有限体积法(FVTD)被推广到计算电磁学中,直接数值求解时变麦克斯韦方程组.FVTD使用电磁场矢量共置于网格单元中心、散射场积分形式的麦克斯韦方程组守恒形,空间离散使用基于特征值的近似黎曼解构建网格单元边界通量,时间推进采用四阶龙格-库塔法.TM和TE波极化下几种二维完全导电体的表面诱导电流密度和雷达散射截面(RCS)计算验证表明,时域有限体积法是一种高精度有效的时域方法.     

带座舱飞船高超声速粘性绕流的分离形态及流场结构分析

通过求解薄层近似的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带座舱飞船高超声速粘性绕流，首次系统地给出了攻角为５°、１０°、２０°、３０°和４０°情况下的物面分离形态、横戴面流线形态、纵剖面流线形态以及旋涡沿涡轴的演变形态。分析和比较表明，本文采用数值可视化方法给出的流场拓扑结构与相应的拓扑规律完全吻合。在飞船后体流场横截面流线中，本文发现了稳定的极限环流态，该极限环涡心处沿涡轴方向的速度是超声速