用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

三维超声速混合层内标量混合的大涡模拟

对三维对流马赫数0.62的超声速混合层流动的标量混合进行大涡模拟, 控制方程对流项采用五阶精度的WENO格式求解, 小尺度涡的作用采用一方程LDKM亚格子模型处理, 过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了混合层流场大尺度拟序结构以及标量场的演化过程, 研究表明标量混合过程受混合层内涡系演化所控制, 标量场具有明显的三维特征.模拟得到的速度、组分及其脉动的统计时均结果和实验结果相符较好.      

跨音速翼型试验侧壁干扰数值模拟研究

采用分区网格技术，数值模拟NACA0012翼型在二维风洞中的侧壁干扰特性。控制方程采用全场三维可压、雷诺平均的Navier—Stokes方程；计算格式采用Jameson中心有限体积、多步Runge—Kutta时间步长格式进行求解；湍流模型采用Baldwin—Lomax模型。通过分析不同迎角、马赫数和展弦比情况下的流动状态，反映了由于侧壁干扰而产生的二维风洞流场的三维性。     

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

AUSM+方法在 粘性内流模拟中的应用

本文阐述了三维RANS方程求解程序中使用的AUSM+离散方法以及Spalart-Allmaras湍流模型，然后对两个算例—Sajben收缩扩散通道内跨声速流动和NASA TP-3252文献中给出的透平静叶排内部流动，进行了数值模拟，并把数值结果与实验结果进行了比较，说明三维RANS方程求解程序可以准确模拟复杂的内部流场。     

突扩管中受限射流的计算方法研究

突扩型燃烧室的结构简单可靠,近年来已引起广泛注意。为了计算时均流场和湍动流场,必须考虑迴流区及附面层全方位的压力梯度,给计算方法带来困难。 本文由时均的N—S方程和—ε湍流模型出发,将Patankar解耦技术推广到可压缩流,采用牛顿迭代法和稀疏矩阵技术,由连续和动量方程直接求解压力和速度场,继而获得温     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

高超声速再入飞行器的多孔流动控制

针对高超声速再入飞行器流动控制问题,提出一种头部进气、尾部排气的多孔流动控制概念。采用Fluent软件针对高超声速再入飞行器无孔模型、排气孔模型在马赫数为6,高度为25km状态下开展了三维雷诺平均N-S方程(RANS)数值模拟研究。结果表明:排气孔模型阻力系数略有增加,下排气孔(SH-DB)模型产生正的升力系数增量和低头力矩增量,上排气孔(SH-UB)模型产生负的升力系数增量和抬头力矩增量,实现了无活动气动面产生控制力矩的目的,为高超声速再入飞行器复合控制提供了参考。      

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。     

Application of Coolant Diffusion Equation in Numerical Simulation of Flow in Air-cooled Turbines

This article deals with time-averaged coolant diffusion equation in arbitrary, body-fitted coordinates. Schmidt number and Lewis number are used to obtain the mass diffusion coefficient. Other than the existing N-S solver, a new solver for the diffusion equation and the influences of cooling air on the thermal parameters, such as gas constant, constant-pressure specific heat, is established. Capable of simulating the flow of multi component fluid, this solver is used to simulate the flow in air-cooled two-stage turbines. The comparison of test results from the modified solver to those from the traditional N-S solver has evidenced that a marked difference present is in thermal parameters and flow field structure.     

叶片安装角槽式轮毂处理的三维数值模拟

采用三维定常不可压Ｎ－Ｓ方程对叶片安装角槽式的轮毂处理的叶栅通道进行了数值模拟,分析了叶片安装角槽式的轮毂处理的扩稳机理,在计算中采用了非交错网格的ＳＩＭＰＬＥ算法,并解决了因使用非交错网格离散方程而带来的压力场振荡的抑制问题,紊流模型采用ｋ－ε两方程模型。     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

跨声压气机内确定应力耦合问题研究

基于三维定常Denton程序发展了二维定常、二维非定常、三维非定常计算程序,以二维跨声压气机叶栅、三维跨声单级压气机近失速点为算例,进行确定应力耦合问题的研究.结果表明:在确定应力应用到时均方程的耦合处理中存在一些因素导致时均方程的解和非定常解的时均值有一定偏差,它们分别是空间差分离散引入的误差、黏性相关项引入的误差、掺混面处理引入的误差.在多排耦合问题中掺混面的处理方式至关重要,对耦合计算结果有较大的影响,简单的掺混面处理并不能获得理想的效果.      

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

壁面振荡对边界层分离的影响

流动分离的控制是当今流体力学领域的一个十分重要的研究课题,本文研究外激发小扰动对基本定常分离流的二阶时均整流效应,结果表明:在合适的St/Re比值下,一个足够大振幅的激发扰动可以推迟甚至抑制分离的发生。     

湍流通用物理方程的推导及讨论

湍流的通用方程已经获得。利用变形率张量与摄动变形位移的点积分解湍流脉动速度,并使用密度加权的侧偏系综平均方法取代Reynolds平均方法,对Navier-Stokes方程进行平均。在二阶精度上没有出现未知相关量,从而获得了具有二阶精度的封闭的湍流通用物理方程。方程中无相关量,不引入任何经验系数,无需构造任何湍流模型,即可描述任意复杂湍流场并解释湍流物理本质。方程在长时标上描述湍流平均行为,在短时标上描述相应尺度层次的湍流结构现象,从而统一了湍流的统计及结构二大学派。方程简化后得到湍流的规范方程,是混合的Burgers-Korteweg-de Vries方程,对该方程的分析给出了湍流转捩准则。文中对湍流的间歇性、拟序结构、各向异性、湍流能量逆转等现象的物理本质进行了初步讨论。      

不稳定飞机气动参数辨识的一种实用方法

地是飞行器气动参数辨识最广泛采用的一种方法。但对于气动高度不稳定的飞机，辨识算法往往因状态方程组和灵敏度方程组积分发散而无法进行。为克服这一数值困难，将方程解耦技术引入输出误差法，发展了高度不稳定飞机气动参数辨识的一种实用方法－－基于方程解耦的误差法。并对某鸭式布局飞机的纵向劝仿真数据进行了辨识仿真计算，计算结果证实了所述方法的有效性。     

稀薄流到连续流的气体运动论统一数值算法初步研究

从非线性模型Ｂｏｌｔｚｍａｎｎ方程出发,引入简化速度分布函数、使用离散速度坐标法对速度空间进行离散、降维,去掉分布函数对速度分量的连续依赖性;采用时间分裂法,将简化速度分布函数松驰变化方程分解为源项碰撞变化方程、对流运动方程,进行耦合计算,应用ＮＮＤ耗散差分方法直接模拟气体分子速度分布函数;发展离散速度数值积分法,通过宏观取矩获取物理空间各点的流动参数,从而建立一套能有效模拟各流域气动问题的简化的     

基于DMD方法的超声速进气道喘振特性分析

采用非定常数值仿真方法对典型超声速进气道的喘振现象进行了研究,并引入动力模态分解(DMD)方法对小喘和大喘流动特性进行了分析,获取了小喘及大喘的流场振荡特征,其中DMD得到的1阶模态反映了时均流场特征、2阶模态反映了主频振荡流场特性。在此基础上,针对小喘与大喘的相互关系进行了研究,结果表明:进气道内小喘流动包含大喘的流场振荡特性,小喘状态是进气道由不喘到大喘状态的中间状态,由小喘向大喘演化过程中,进气道内一些流动特征逐渐减弱并趋于稳定收敛,大喘的流场结构整体上比小喘状态更为稳定。