高超声速进气道内激波－边界层干扰研究

为了研究高超声速进气道内激波干扰和流场结构，对一个二维进气道模型在中国空气动力研究与发展中心超高速所风洞流场中进行了实验（Ｍ＿∞＝５．０，Ｒｅ＿∞＝７．５×１０￣６ｍ￣（－１））。采用纹影方法，结合压力测量对进气道内流结构进行诊断，同时用Ｙｅｅ的ＴＶＤ有限差分格式通过求解二维Ｎ─Ｓ方程对相应实验状态下的进气道内流场进行了数值模拟。风洞实验与数值模拟取得了较好的一致。进一步的分析还表明高超声速进气道内的激波─边界层干扰对进气道性能有很大的影响。     

翼身—平面角区湍流场研究

实验在低湍流风洞中完成．对翼身－平面角区流场里诸参数，如翼型表面和平板上的压力分布、绕翼型及尾迹区内的平均速度、脉动速度、湍动能及二阶关联量ｕ′ｖ′及ｕ′ｗ′进行了测量，然后应用拟压缩性方法，在曲线坐标系下求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，选取Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模式，对角区湍流场进行数值模拟。最后把计算结果与实验结果进行了对比，二者符合较好。     

跨声速三维非线性洞壁干扰的数值计算

分别以固壁条件和洞壁附近的压力分布模拟各类实壁和透气壁试验段的洞壁边界条件，利用Ｅｕｌｅｒ方程和Ｎ－Ｓ方程数值求解模型在风洞中的绕流场，得出洞壁干扰对跨声速模型绕流和气动力的影响，初步的研究结果表明，该方法能较有效地模拟模型在跨声速风洞中的绕流场，经洞壁干扰修正后的ＧＢＭ－０４Ａ模型在０．６ｍ风洞中的试验结果与无干扰参考结果吻合较好     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

真空羽流场的 N-S 和 DSMC 耦合数值模拟

根据调姿发动机真空羽流场的特点,采用差分求解Ｎ－Ｓ方程数值模拟喷管内流场,用ＤＳＭＣ方法数值模拟喷管外羽流场,得到了符合流动规律的数值模拟结果。为姿控发动机喷管设计提供了参考,并为真空羽流污染研究打下了良好的基础。     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件,采用数值求解Ｅｕｌｅｒ方程的方法模拟模型在风洞中的绕流场,然后将洞壁边界条件用远场边界条件替换,用同样方法算出模型在自由流中的绕流场,从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性,采用重叠网格技术,对飞机大迎角标模和ＳＢ－０３模型进行了亚声速大迎角的洞壁干扰计算与修正。     

透平进气阀箱的数值模拟

用数值方法模拟了透平进气阀箱的内部流动。为了改善计算的收敛性和提高计算精度,程序采用了ＳＩＭＰＬＥＳＴ迭代求解方法,并引入ｋ－ε湍流模型、壁面函数法和空度的概念进行了数值计算。计算结果表明,本文的方法能够较好的模拟进气阀箱内部复杂的流动,并与已有的一些实验结果吻合良好。     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

基于次流喷射控制推力矢量喷管的实验及数值研究

应用实验和数值模拟的方法,对一种新型的推力矢量喷管—基于次流喷射控制的二维推力矢量喷管的推力矢量性能和流场进行了研究。实验是在西北工业大学小型超高速吹气式风洞中进行,测量了在不同的二次喷流情况下,推力矢量和流场的变化规律;采用时间推进求解N-S方程的方法数值模拟了二维推力矢量喷管内流场和性能。研究结果表明,应用次流喷射控制主流流动可以实现较大的推力矢量转折,但是,二次喷流必须具有足够的压力值;如何从推力矢量工作方式恢复到轴向流动工作方式则是需要进一步研究的问题。      

跨声速开槽壁超临界翼型洞壁干扰的N－S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一步，其结果就符合得更好，这为下一步用Ｎ－Ｓ方程进行洞壁干扰修正打下了基础。     

压气机叶片排内三维紊流流动的隐式高分辨率数值分析

对跨音速压气机叶片排内具有复杂三维激波的紊流流场所进行的数值分析，采用了空间三阶精度无波动格式以保证对紊流特征的正确模拟，利用ＬＵ－ＳＧＳ算法避免了一般隐式算法的繁琐运算，使之与显式方法推进一步的运算工作量基本相当，从而构成了一种既准确又高效的跨音速压气机紊流流动数值求解系统。计算与实验进行了对比，结果是令人满意的。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

喷管二维跨声速两相湍流流场的数值模拟

采用ＳＩＭＰＬＥ算法,结合分散颗粒轨道（ＰＳＩＣ）模型,在同位网格上实现了喷管中跨声速两相湍流流场的数值模拟。根据Ｒｈｉｅ,Ｃｈｏｗ提出的动量插值概念,以协变物理速度分量作为求解变量,完成跨声速气相流场的数值计算。湍流模式采用ｋ－ε双方程模型,计算结果与提供的实验数据相符合。     

运输机翼型结冰的计算和实验

采用数值计算和结冰风洞实验相结合的手段,对某运输机机翼剖面缩比翼型的结冰特性进行了研究.结冰实验在气动中心0.3m×0.2m结冰风洞中进行,实验模型为弦长0.18m的层流翼型.翼型的绕流流场通过求解低速黏流的时均Navier-Stokes(N-S)方程得到,采用拉格朗日法计算水滴撞击特性,在此基础上,求解结冰热力学模型...     

大飞机布局模型跨声速风洞实验尾支撑干扰研究

对某大飞机布局风洞实验尾支撑干扰开展了数值模拟和实验研究,发展的数值方法计算结果与风洞实验结果有很好的一致性。对于类似构型的飞机,在迎角-2°～6°范围,可认为尾支撑干扰量随迎角呈线性变化,采用前位叶片支撑作为辅助支撑带来的二次干扰量可以忽略,新型双天平辅助支撑系统试验进一步验证了这一结果;尾支撑对机身、尾翼、机翼等部件的绕流都有影响,干扰量随构型而变,对阻力、力矩影响较大,且随Ma数变化,因此不同构型实验数据需要单独修正。所发展的带风洞支撑系统的数值模拟软件能够满足工程应用要求,可用于支撑干扰修正研究以及风洞实验支撑系统优化设计。     

带进气道的隔离段流场实验研究与数值模拟

通过风洞实验与数值模拟研究了超燃冲压发动机带进气道的隔离段流动。结合RNGk ε模型隐式求解了三维N S方程,并将计算结果与实验结果进行了比较。研究发现,RNGk ε模型能较好地模拟出超声速流动的激波与分离现象;进气道喉道流动的非均匀性使隔离段内激波/附面层干扰流场与均匀来流条件下的流场有显著差别。风洞实验表明,同一个模型,风洞马赫数为3 85试验的隔离段入口压力畸变大于马赫5 3;但前者隔离段出口截面压力分布比后者更均匀;隔离段入口畸变度大,隔离段实际能达到的压升就低。研究表明隔离段内的总压损失在整个进气道 隔离段组合体总压损失中占了相当大的比重。     

高温燃气流风洞非稳态运行过程实验和数值模拟研究

针对热考核用高温燃气流风洞运行过程中的非稳态过程开展研究,采用风洞实验与数值计算相结合的方法研究了风洞整体起动、扩压器背压抬升、关机三种过程中喷管、实验舱和扩压器内瞬态流动特性。实验结果揭示了燃气风洞诸多有趣的瞬态现象,且披露了一手实验数据。而通过数值计算能较好地复现上述瞬态现象;借助数值计算能合理地解释试验现象产生的原因;还能捕捉实验中无法观测到的现象;此外还评估了该风洞扩压器的抗反压裕度在10kPa以上。因此,数值方法是研究大型燃气风洞瞬态流动特性最重要的辅助手段。该研究可为类似风洞运行调试提供借鉴,并为风洞实际运行提供数据支撑。     

跨声速开槽壁超临界翼型洞壁干扰的N—S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一     

数值模拟模型尺度与来流条件对实验数据的影响

采用商业计算软件CFD-Fastran,数值模拟了高超声速一体化飞行器在真实飞行以及风洞实验来流条件下的流场,通过对壁面压力与摩擦力进行积分,得到全机的阻力系数、升力系数以及俯仰力矩系数。对不同来流条件以及模型缩比尺度的气动数据相关性进行研究,其中风洞来流条件的马赫数与静压与真实飞行条件相同,但总焓以及模型尺度不同。数值模拟结果揭示了风洞来流条件以及模型缩比尺度对飞行器整体气动性能的影响,表明发动机内流道摩阻系数计算结果的不同是发动机通流状态下全尺寸模型与缩比模型的阻力系数不同的主要原因。该结果为风洞实验数据的修正提供了参考。     

再入飞行器迎风面边界层转捩的数值预估

基于求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,首先数值预估了高超声速三维流动的边界层特性;采用已有的高超声速三维边界层转捩准则,即在转捩位置Ｒｅθ／Ｍｅ＝常数,以及通过修正由ＯＮＥＲＡ提出的代数转捩模型,建立了预估三维边界层转捩的半经验数值技术并应用于三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的求解。利用ＨＡＬＩＳ模型及相关的实验数据,首次进行了在Ｍ∞＝６时不同雷诺数及迎角的验证试验。预估结果表明：在不考虑