基于自由尾迹和升力面方法的双旋翼悬停气动干扰计算

建立了一个同时适合于共轴式、横列式及纵列式直升机双旋翼气动干扰分析的计算方法。在该方法中,为更好地模拟气动干扰特性及桨尖三维效应,桨叶模型采用了升力面/涡格法,尾迹则使用畸变的自由尾迹模型。通过旋翼下洗速度的计算值与可得到的实验值对比,验证了计算方法的有效性。应用该方法,以横列、共轴、纵列式双旋翼为例,分别计算了悬停状态双旋翼的尾迹特性及诱导速度变化,给出了部分尾迹边界,并对计算结果进行了分析。最后,给出了几点结论。     

动力增升飞机气动力计算方法研究

分析了国外动力增升飞机气动力计算方法,提出了一套新的动力增升飞机快速计算方法。新方法解决了传统方法在小速度下气动力发散,多速度点求解气动力带来计算效率低等问题。以某动力增升运输机为算例,验证了新方法的正确性。     

一种新的适用于叶栅绕流计算的网格生成法及其应用

本文对平面叶栅气动计算网格生成方法作了简要评述。在此基础上,根据叶栅槽道流动的特点,提出一种以等比级数为基础,用几何方法构作H网格的数值方法。用这种网格进行流场计算的结果表明,所生成的网格在叶片前后缘区有很高的分辨能力,能够算出速度分布的细微结构。本文的计算结果与亚音速及跨音速叶栅速度图解符合得很好。     

带有双层壁扩压器的波瓣喷管混合流场的数值计算

采用三维贴体曲线坐标网格,边界网格加密且正交,在整个计算区域进行全场计算。为避免因波瓣造成的网格强烈的非正交引起解的发散,采用了Chen和Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法,对带有双层壁扩压器的波瓣喷管混合流动进行了数值计算和分析。对波瓣及双层壁,采用大粘性的方法解决流固耦合。计算结果表明:在双层壁间有外界冷气流被引射进入,形成了壁面的冷却气流;自波瓣出口截面沿流向产生的环流速度场,强化了主次流的掺混,速度分布渐趋均匀。计算结果与实验数据符合良好,二者在离开波瓣52 mm的混合管内相对主流速度的最大误差为21.65 %。     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

双旋翼自由尾迹计算分析中的加速算法

对快速多极展开算法(FMM)在共轴式和纵列式双旋翼自由尾迹互诱导速度计算中的应用进行了研究。介绍了FMM算法的基本原理,建立了FMM算法的数学计算模型,并对其加速特性和准确性进行了分析;将方法引入到基于直涡元的自由尾迹互诱导速度的计算分析中,推导了诱导速度的多极展开表达式。在此基础上,以某旋翼模型为例,通过对自由尾迹中涡元诱导速度的计算,验证了加速算法的加速度特性和准确性。使用方法对共轴式和纵列式双旋翼悬停和前飞状态下的自由尾迹几何形状进行了计算,分析了双旋翼之间的相互干扰。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

波瓣喷管双层壁扩压器流场的数值分析

采用三维贴体曲线坐标网格，波瓣边界网格加密且正交，在整个计算区域进行全场计算。文中采用了Chen-Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法，对带有双层壁扩压器的波瓣喷管流动进行了冷热态数值计算和分析。对波瓣及双层壁，采用大粘性的方法解决流固耦合，计算结果表明：在双层壁间有外界冷气流被引射进入，形成了壁面的冷却气流，相对单层壁扩压器，双层壁扩压器的壁面温度明显降低；自波瓣出口截面沿流向产生的环流速度场，强化了主次流的掺混，速度分布渐趋均匀。通过试验测试，计算结果与实验数据符合良好，二者在离开波瓣420mm的混合管内相对主流速度的最大误差为15．35％。     

基于递推傅里叶变换的飞行器参数在线辨识方法

在线辨识在现代飞行控制系统设计中扮演越来越重要的角色,飞行器模型的在线更新使得人们可以采用更智能的控制方法。基于计算精度和速度的考虑,在线辨识方法通常以递推方式进行,主要分为时域和频域两大类方法。在建立飞行器系统模型结构的基础上,利用频域递推傅里叶变换及最小二乘方法,实现对气动及控制偏导数的在线辨识。针对某飞机纵向通道的在线辨识仿真验证该方法有效,且计算速度和收敛速度快,辨识结果与参数真实值之间的一致性好,方法对传感器噪声有较强的适应性。最后,分析比较了各种典型激励信号对辨识结果的影响,为进行实际在线辨识试验提供了参考依据。     

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

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

弹道式地地导弹精度计算公式及其拟合方法

归纳和总结厂弹道式地地导弹的精度变化规律,根据这些规律和谐波分析理论建立了导弹精度计算公式。文中根据函数逼近理沦,提出了确定谐波函数阶次的方法和具体的拟合步骤。应用这种方法导出了实际导弹的精度计算拟合公式.采用这些公式和采用弹道计算方法分别计算任意相同条件下的导弹精度值.其相对误差小于1%。本文提出的精度计算拟合公式是简便实用的.     

量子计算及其在空气动力学中的应用前景综述

量子计算是最重要的后摩尔计算技术之一，拥有经典计算机无可比拟的超强计算能力，未来能够对各行业应用产生颠覆性的影响。针对量子计算给空气动力学带来的机遇和挑战，详细综述了量子计算机、量子算法、量子底层软件栈等方面的研究进展。结合空气动力学领域常用的基础方法，在综述量子计算线性方程组求解、插值操作、数值积分、优化搜索等最新进展的基础上，结合典型量子算法深入分析了量子计算在空气动力学领域的应用前景，并指出了需要重点关注的研究方向。     

加油机尾流场建模与仿真分析

提出了等效气动效应法,建立了加油机尾流场的等效扰动模型。等效气动效应法采用逐点积分方法和加权平均方法,将尾流的影响等效为受油机运动方程直接可用的平均风速度和风梯度。对不同状态下受油机所受尾流的扰动作用进行了仿真。结果表明,扰动模型与实际情况比较符合。     

颤振板的极限环振动

本文研究了在匀速气流中悬臂板的极限环振动。采用冯·卡门大挠度板理论和准定常空气动力,应用瑞雷-里兹法导出悬臂颤振板的非线性振动的计算方程。并提出一组模态函数展开式用于拉格朗日方程,从而确定随时间历程的极限环振动和颤振边界。计算表明,模态展开式具有收敛性,板的长宽比对极限环振幅有显著影响。     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

方向交替隐含法并行计算的一种模型

介绍了流体力学计算中运用的方向交替隐含法和求解三对角方程组的追赶法，提出了利用方向交替法并行求解流体力学问题的一种模型。该模型将计算中的通信复杂度降到了很低的情况。给出了该模型的节点机互连结构及其数据的划分及存储方式，确定了共享数据和其主、从节点，讨论了共享数据的交换方法，并给出了并行计算的实现细节。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

Navigation Computation in Terrestrial Strapdown Inertial Navigation Systems

The problem of computing the translational velocity and position relative to earth, which has to be solved by the processor of a strapdown inertial navigation system is discussed. Several approaches are briefly examined with consideration given to the form in which the sensor data are generated and to the computational burden involved in each approach. A computational scheme is finally selected in which the computation is divided into three rate levels. The differential equations of this scheme are developed and the assumptions on which the development is founded are stated. Three variants of the basic scheme are presented, each based on a different level of simplifying assumptions. The main purpose of this work is to develop the differential equations to be solved at each stage of the computation, rather than the numerical implementation of the solution. This work supplies the theoretical background for some of the numerical methods which are now being used.     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。