椭圆截面喷管排气辐射热流数值分析

本文对直升机红外抑制器椭圆截面喷管排气的辐射热流进行数值分析,并预估排气温度分布速度和燃气成分分布对辐射热流的影响。为了正确计算辐射热流,本文采用六热流辐射模型来考虑尾焰对周围环境传热的影响。紊流模型采用修正的双方程k—ε模型。由于喷管出口截面为椭圆形,因此排气为三元自由射流,控制方程为抛物型偏微分方程,可采用前进积分求解。计算结果与试验结果基本一致,说明本计算方法是可行的。     

用预处理方法及可压Navier-Stokes方程计算低马赫数流动

用可压方法及预处理方法求解了低马赫数紊流流动,其中空间离散格式采用了Roe格式,紊流模型采用了DES模型.以绕NACA0012的紊流流动为例,分析了该方法对提高流动速度和精度的特性,并将其用到多段翼型及绕圆柱的非定常流动中.数值模拟结果表明,该方法能有效求解低速流动及混合流场问题,且具有较高的数值精度.     

SOLUTION OF 3-D TURBULENCE NAVIER-STOKES EQUATIONS USING HYBRID GRIDS

INTRODUCTIONWith the development of computer,consid-erable process has been made in the developmentof numerical methods for the solution of Navier-Stokes equations.However,most of these me-thods are impractical for complicated flows in adesign environment.The primary reason is thatwe can′t generate a suitable mesh easily for acomplicated configuration.Although there aremany techniques,such as structured grid or un-structured grid,which have been successfullyused for numerical simulation…     

旋翼与紊流场干扰噪声计算

本文给出了旋翼与紊流场相互作用产生噪声的理论计算方法。该方法从翼型段在紊流场中作直线运动时的噪声计算分析出发，由理论分析推广到旋翼旋转运动的情况，其中考虑了桨叶上非定常载荷的弦向及展向的非紧致性，叶－叶载荷之间的相关，特别是考虑了由旋翼运动造成的素流场畸变所引起的非各向同性和非均匀性。本文计算了一架模型直升机在＂准悬停＂状态下紊流场的畸变以及相应的干扰噪声谱，并对考虑畸变与非畸变、均匀畸变与非均匀畸变的结果作了比较。结果表明，紊流场的收缩畸变对干扰噪声有显著影响，并对此给出了合理的解释。     

盒式机翼布局气动特性研究

采用数值模拟的方法研究了盒式机翼布局的气动特性。在数值模拟过程中,控制方程为雷诺平均N-S方程,采用基于非结构网格的有限体积方法进行求解。空间离散采用了Osher逆风格式,时间离散则采用五步Runge-Kutta格式。紊流模型采用了S-A方程模型。对盒式机翼布局外形在不同状态下的流场进行了数值模拟,给出了每个状态下的气动参数,并进行了分析,结果显示盒式机翼布局具有较好的气动特性。     

槽道湍流相干结构的研究

采用高阶差分方法求解非定常不可压Navier-Stokes方程，进行槽道流相干结构的数值研究，与常用的谱方法相比，不需限制边界条件，可用于更一般的粘性流体的流动问题，文中所用时间分裂法和三维耦合高阶差分格式适用于包括邻近边界点在内的整个计算区域，克服了三维各自用高阶中心差分格式不适用邻近边界点的困难。通过算例模拟了槽道流动中单个相干结构的演化发展规律，结果与实验一致。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。     

从壁面脉冲到湍流相干结构演化的数值模拟

发展了高效求解基于Navier-Stokes方程的扰动方程高阶紧致差分方法，研究了壁面强脉冲到湍流相干结构的复杂演化过程。该方法采用高精度、高分辨率的三维耦合中心差分格式，以及和高阶时间分裂法相结合因而更为有效，比一般谱方法更适用于较为广泛的流动区域和边界条件的湍流相干结构研究。中模拟了槽道流动中壁面强脉冲诱导的湍流相干结构的产生、发展和演化过程，显示了湍流相干结构的基本特征和发展规律．计算结果表明，所用壁面强脉冲模型比共振三波模型更有效．     

矩形翼飞机对大气紊流的动态响应

本文采用时间历程法和谱分析技术,计算了一架典型矩形翼飞机对大气紊流的纵向动态响应,讨论了展弦比对动态响应的影响,分析了带有理想驾驶仪方程时对响应的改善。由于计算是从基本飞行运动方程出发,故对紊流的离散型和连续型有统一的求解方法。     

基于格子Boltzmann方法的可压缩翼型绕流模拟

采用能够恢复可压缩Navier-Stokes方程的耦合双分布(DDF)格子Boltzmann方法(LBM)模拟了贴体网格下的NACA0012翼型绕流.首先在Ma＝0.5,α＝0.0,Re＝5 000状态下通过模拟低雷诺数流动检验了该方法;然后分别在Ma＝0.5,0.85和1.2,α＝-0.05,1.0和0.0状态下进行了模拟,并在Ma＝0.85时通过加密网格得到了更好的结果.为了避免转捩问题采用了无粘边界条件进行模拟,但得到的压力系数与实验结果吻合良好,证明了该方法的可靠性并显示了其模拟可压缩的潜力.该方法在翼型绕流上的应用为湍流流动的数值模拟提供了基础.     

求解织物热湿耦合方程的控制体-时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方程的算法。首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化连续性的非线性微分方程组为一系列的递归形式的线性代数方程组。文中给出了问题的求解步骤和算例。结论表明该方法的预测结果具有不依赖时间步长的特点。     

Highly-Efficient Aerodynamic Optimal Design of Rotor Airfoil Using Viscous Adjoint Method

In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型，采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式，对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟，探讨了内外喷流及喷管内的流动特征，并与实验数据进行了比较，计算结果与实验结果吻合较好，验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性，从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

飞翼无人机低速着陆状态抖振响应研究

基于RANS方程和SST湍流模型,采用飞翼无人机升力系数判据、俯仰力矩判据和表面极限流判据等对抖振始发迎角度进行了预测估计。基于较为详细的结构模型和气动模型,构造了结构与气动的耦合求解技术,同时采用弹簧近似光滑和局部重构组合方法对气动的动网格技术进行更新;然后分别在时频域内分析了飞翼无人机的抖振载荷响应。研究结果表明:基于CFD的RANS方程、SST湍流模型及各类判据预计刚性飞翼无人机在低速大迎角状态下的抖振始发迎角可以得到较为合理的结果;采用升力系数判据、俯仰力矩系数判据预测的抖振始发迎角与表面极限流判据预测的始发迎角相比要保守一些;而表面极限流判据方法能给出翼面气流流动的细节,采用此判据可以分析对比不同迎角下的翼面气流流动变化情况;与气动弹性及嗡鸣响应相比,抖振是一种强迫振动,其响应频率并不单一。     

软管-锥套式空中加油动态建模与性能分析

提出了一种基于凯恩方程的自主空中加油软管-锥套动态模型及性能分析方法。该方法将软管视为由有限段以铰链形式连接的刚性 杆组成，锥套为软管末端的一个质点。本文定义了描述系统状态的广义坐标及广义速率，导出了软管段位置多级递推公式和系统动态方程，估算了软管在加油机尾流、定常流和大气扰动下的气动载荷。通过数值仿真分析了平稳大气中锥套在加油机不同飞行条件下的稳态阻力和软管拖拽轨迹，验证了模型的正确性及系统稳定性。最后研究了大气扰动对锥套运动的影响及不同软管段的受扰运动。     

基于空间相关函数的二维紊流场数值生成法

大气紊流是复杂的扰动形式，影响飞行品质、飞行安全。对于复杂的飞行任务，有必要生成二维紊流场模型，从而更准确地反映扰动情况。传统方法是以时间频谱为基础，利用谱分解定理得到所需的成形滤波器结构，进而生成紊流递归模型。这种方法需要经过对模型进行简化，以便利用谱分解定理，使仿真模型的统计特性受到影响；另外，从理论上讲，传统方法只能生成一维大气紊流递归模型。本文直接从空间相关函数入手，提出建立离散自递归模型的新方法，给出了二维紊流场的递归模型     

一种基于EMD的低湍流度信号处理分析方法

采用热线风速仪在3座典型低速风洞中进行了流场湍流度测量，这3座低速风洞包括1个低湍流风洞、1个常规闭口风洞和1个开口射流风洞。针对湍流信号通常受噪声干扰的问题，在湍流度值处理中引入了经验模式分解（EMD)自适应滤波和HHT时频谱分析方法。将EMD方法与其他几种湍流度值处理方法进行了比较，包括带通滤波方法（BP)、电磁噪声解耦方法（ENC)和高通惯性衰滤波方法（HPIA)。采用EMD方法测得低湍流风洞的湍流度值，在流场速度30~100 m/s的范围内小于0.05%。采用HHT方法完成了脉动速度信号的时频分析，分析发现开口风洞试验段的脉动速度HHT时频谱有突出的低频信号。所构建的EMD自适应滤波器可以有效控制噪声对热线输出信号的影响，是一种有效的低湍流度信号处理方法。     

多变量自回归模型的快速辨识方法

本文提出了多变量AR(Autoregressive)模型建模的快速算法,该算法采用解超定矩阵的最小二乘方法,并利用正交变换技术,从而避免了最小二乘估计中矩阵求逆的病态问题,保证了数值计算的稳定性。还介绍了一种可转化为单变量建模的多变量AR模型,它可以直接利用单变量建模模块。一般情况下,可将序列进行典则分解。最后简要列出了有关结论。