不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

用三维N—S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松驰迭代来实现方程的求解，湍流的计算采用了标准的双方程模型，本文给出了钝体和类小汽车外形的计算结果，钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。     

用N-S方程计算三维叶栅流

发展了基于微元中心有限体积方法,采用雷诺平均N-S方程求解三维叶栅流动的程序,特别注意了格式中人工粘性的加法;程序采用Baldwin-Lomax模型。格式经过平板层流及湍流流动的校验,层流与Blasius解、湍流与实验经验公式及Wilcox的计算值符合很好,最后以Goldman涡轮导向器为例进行了叶栅流动计算校核,结果符合于实验以及Chima等人的计算。      

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

μ理论在电液负载模拟器中的应用

 研究μ理论在电液负载模拟器中的应用，提出电液负载模拟器的鲁棒控制策略。电液负载模拟器是一个复杂的机-电-液复合系统，且是一个强耦合、时变受控对象。综合考虑参数变化、模型变动和干扰等不确定性的影响，利用μ综合控制理论设计电液负载模拟器的鲁棒力控制器，并通过μ分析使用鲁棒力控制器时系统的鲁棒稳定性和鲁棒性能。给出使用鲁棒力控制器和经典力控制器时的实验结果，实验结果表明所设计鲁棒力控制器的有效性和优越性。     

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

基于马赫数分布可控曲面外/内锥形基准流场的前体/进气道一体化设计

提出了一种高超声速飞行器乘波前体的外锥形基准流场设计方法,在锥面马赫数分布规律给定的条件下,通过有旋特征线法实现反设计,提高了基准流场设计的灵活性。该基准流场通过锥形"下凹"弯曲激波和波后等熵压缩波系压缩气流,可以在较短的长度内完成高效压缩。基于反正切马赫数分布外锥形基准流场设计的乘波前体具有较高的容积率,乘波特性良好且出口均匀,设计点时有黏升阻比为1.89。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。     

Enhancing the efectiveness of film cooling

Advanced gas turbine stages are designed to operate at increasingly higher inlet temperatures to increase thermal efficiency and specific power output.To maintain durability and reasonable life,film cooling is needed in addition to internal cooling,especially for the first stage.Film cooling lowers material temperature by forced convection inside film-cooling holes and by forming a layer of coolant about component surfaces to insulate them from the hot gases.Unfortunately,each cooling jet forms a pair of counter-rotating vortices that entrains hot gas and causes the film-cooling jet to lift off from the surface that it is intended to protect.This paper gives an overview of efforts to enhance the effectiveness of film-cooling.This paper also describes two new design concepts.One design concept seeks to minimize the entrainment of hot gases underneath of film-cooling jets by using flow-aligned blockers.The other design concept shifts the interaction between the approaching hot gas and the cooling jet to occur further above the surface by using an upstream ramp.For both design concepts,computational fluid dynamics results are presented to examine their usefulness in enhancing film-cooling effectiveness.      

Potentials of manufacture and repair of nickel base turbine components used in aero engines and power plants by laser metal deposition and laser drilling

IntroductionExpensive turbine parts like HPT(HighPressure Turine)blades or vanes are replaced bynew parts in case of damage.For example theburn through of the inner side of a blade or vane(Figure 1)is a frequently appearing damage,which cannot be repaired…