渐扩通道内湍流流动与混合的数值模拟

应用标准形式的ｋ－ε湍流模型,对具有逆向压力梯度的渐扩通道内两股等温或不等温射流的湍流流动与混合进行了数值模拟。结果表明,时均速度场与温度场的计算结果与实验测量数据相符合。     

环形水射流流场的实验研究与统计分析

在11、12和15MPa射流压力下对环形自由水射流流场进行实验研究,以获得该射流场的能量特征与液滴尺寸分布规律。运用相位多普勒粒子测速（PDPA）技术对射流流场中的速度分布和液滴粒径分布进行测量,并对通过不同位置控制体的单个液滴行为进行了统计分析。研究表明：环形射流流束中心存在着较宽的高速区域,且射流能量沿着射流方向衰减缓慢;距离喷嘴的轴向距离越大,射流横断面上的速度与液滴粒径分布越平坦;射流压力对流束中心的轴向速度变化影响较大,对液滴粒径分布的影响不明显;射流流束中心区的湍流脉动较弱,但通过位于射流中心位置的控制体的液滴粒径谱较宽。     

等离子射流在圆柱充液室中扩展、掺混的实验及数值模拟

为了研究等离子射流与液体工质相互作用特性,设计了圆柱型充液室,运用数字高速录像系统观察了等离子射流的扩展过程,重点研究了放电电压对等离子射流扩展特性的影响。实验结果表明,等离子射流在液体中扩展时,两相流体存在较大的速度差,Taylor-helmholtz不稳定效应导致强烈的两相湍流掺混,Taylor空腔的亮度在时空上呈非单调分布;放电电压越大,Taylor-helmholtz不稳定效应越强。在实验基础上,建立了等离子射流扩展的二维轴对称非稳态可压缩流的数学模型,并进行了数值模拟,获得了射流场中Taylor空腔界面的形态变化,以及流场中压力、速度、温度、湍动能和湍流耗散率的分布特性。由Taylor空腔界面扩展云图计算得到的Taylor空腔轴向扩展位移与实测值吻合较好。     

三维亚声速冲击射流流场的数值模拟

利用有限体积法对三维可压缩的N－S方程进行离散，对喷嘴的亚声速垂直冲击射流和斜冲击射流进行了数值模拟。网络的划分采用非结构性网络，湍流模型为RNGκ-ε湍流模型，在近壁处采用壁面函数进行修正，流场计算结果与实验值吻合一致。计算可得到垂直冲击流场结构的三个区域：自由射流、滞止区和壁面射流区域，并显示出垂直冲击点附近区域的卷吸回流、速度分布。冲击射流近壁处存在着回流现象，由于边界层流动和近壁回流引起距离固壁2D内的速度梯度变化较大。斜冲击在近壁处平行壁面速度的最大值，以及壁面上的总压最大值都偏离冲击点，并且相对应。     

室内气流分布的大涡模拟

准确地模拟室内气流分布对优化室内通风设计及预测舒适和健康的室内环境有着重要的意义，传统的空内气流分布的数值模拟大都采用к-ε两方程湍流模型，该方法求解出的是流动变量的平均值，无法给出流场结构的详细信息，体现不出湍流流动的瞬时性的特点。本文采用大涡模拟（LES）法求解室内气流分布，使HSMAS法保证计算的稳定性和收敛性。通过对室内等温贴附射流和非等温受限射流两种类型的射流作用下的室内气流分布计算，以及导入к-ε两方程湍流模型对计算结果进行比较与分析表明，大涡模拟求解室内气流分布是一种确定可行和有效的方法。     

液膜射流和圆射流雾化过程的湍流演化特性

利用PIV技术,在相同喷射压力下对分属液膜射流和圆形射流的两种理论喷射类型的典型喷油器———涡旋喷孔和双喷孔喷油器的射流雾化的湍流演化过程进行了测量和分析,研究表明:在整个射流过程中,两种射流的空间相对湍流强度分布比速度分布发散,且集中分布在射流的上下边缘和根部区域;随着喷射时间的延长,射流流场的湍流强度值逐渐增大,且主要分布在涡旋射流的非前锋区域以及双孔射流的非前锋和非中心线区域;喷射开始后1.5ms和4ms时,涡旋射流中心区域的湍流积分尺度分量值较大,分别达到约6mm和7mm的量级,8ms时,涡旋射流流场的湍流积分尺度分布的发散程度比双孔射流明显,且雾化粒子速度方向与积分方向不同时得到的两个湍流积分尺度分量大小以及连续性等分布情况恰好相反;整个喷雾过程,雾化粒子y方向速度分量沿x方向积分得到的湍流积分尺度分量值均小于其他3个湍流积分尺度分量值且分布的不连续性明显。     

亚声速旋拧射流噪声中的温度效应

采用大涡模拟（LES）方法模拟亚声速旋拧射流,着重考察温度效应对旋拧射流近场流动演化过程、湍流脉动空间发展和远场噪声的影响。线性稳定性分析表明,旋拧射流中提高射流中心温度会增加剪切层的扰动增长率;数值结果显示,加热会促进剪切层中大尺度结构的产生及相互作用,促使流动更快进入湍流状态,并缩短射流势核区的长度。在初始层流发展阶段,加热会提高中心线上的流向速度脉动峰值,但是对剪切层中的流向速度脉动峰值几乎没有影响;在湍流发展阶段,提高射流中心温度会提高流向速度脉动衰减率,并降低脉动幅值。此外,在非等温射流中,密度脉动幅值要远高于等温射流。在30°方位角附近,等温射流的总声压级幅值最高,冷射流的噪声幅值最低。方位角大于50°时,加热使总声压级降低,且随着方位角幅值的增大,降低越明显;而冷却则会提高总的声压级幅值。     

受限射流的数值模拟

对柱形燃烧室内的轴对称等温射流流场进行了数值模拟。回流是这类流场的普遍特征。它对于燃烧室火焰的稳定起着重要作用。本文研究了影响回流特征的多种因素。以压力-速度为主要变量,采用K-ε二方程湍流模型进行计算。计算方法的特征之一是膨胀流动的阶梯形边界表示法。结果包括回流区特征及预示的平均流线分布图。计算结果与实验资料相当吻合。这说明了此方法可用于燃烧室性能的分析。     

大速差射流燃烧室内三维湍流气固两相流的实验研究

本文用二维激光多普勒测速仪和探针取样法对大速差射流燃烧室中冷态三维湍流气固两相流动进行了实验研究,测量了气相速度,颗粒速度、颗粒质量流分布及两相各自的湍流特性,探讨了高速射流作用下两相流场的主要物理特征及其对火焰的稳定和强化燃烧的作用。     

超声速流场中侧向射流的数值研究

通过数值方法研究拱柱体上侧向射流与超声速来流的相互干扰,为侧向射流喷管出口横截面几何形状的设计提供依据,实现高速飞行器机动飞行的姿态控制。采用多块嵌套网格并行算法,数值求解Navier-Stokes方程以及k-ω湍流模型,分别对四种不同横截面形状(相同面积的圆形、椭圆形、楔形以及倒置楔形)的射流喷管出口形成的流场进行了研究。圆形喷管出口的流场数值计算结果与实验结果相一致,表明本文采用的数值方法可靠以及网格分布合理。对比物面压力分布和射流干扰放大因子发现,对于长轴与来流垂直的椭圆形喷管出口,在射流上游一侧的物面附近形成较大的分离区以及较强的激波,由此产生的高压区有利于侧向射流对飞行器姿态的控制。     

基于能量法和特征值法的颤振预测数值方法研究

为研究适用于工程设计的颤振数值预测方法,采用能量法和特征值法对出现颤振现象的真实叶片进行预测。通过将2种数值方法的预测结果与试验结果进行对比,验证数值预测方法的精度和工程适用性。结果表明:2种方法都较准确的判断了某典型风扇转子叶片的气弹稳定性,与试验结果基本吻合,具有较高的工程实用价值。     

水滴撞击特性的粒子统计法研究

基于水滴撞击特性研究的拉格朗日方法,提出了适合于考察非定常情况的粒子统计法.通过统计非定常流场中不同时刻撞击在防冰表面上不同区域的粒子数,获得防冰表面上粒子浓度分布随时间的变化情况,与基于染色法水滴撞击试验获得的色度分布具有直接对应关系;采用基于染色法的动态图像测量方法对某型发动机帽罩缩比模型开展了水滴撞击特性试验,对所提粒子统计法进行了验证.基于拉格朗日方法的粒子统计法获得的计算结果与试验结果符合较好,与采用拉格朗日轨迹法获得的结果具有较好的一致性.粒子统计法所模拟的物理过程与试验过程一致,在处理湍流随机性和来流液态水含量分布不均等方面具有优势.     

基于变域变分有限元的翼型反设计

反问题作为一种自由(未知)边界问题近年来得到了广泛的研究。变域变分通过把自由边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而自由边界求解和流场分析可以完全耦合进行。本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型，计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果。文中还采用了两种计算网格(H型和O型)来设计翼型，比较了它们的精度，计算表明本方法完全能用于翼型的反设计。     

连续变迎角测力试验技术在大型暂冲式跨声速风洞中的应用

由于暂冲式高速风洞运行时间短暂,普遍采用阶梯变迎角方式进行静态测力试验,其试验信息量难以满足先进飞行器研制的试验需求。为在暂冲式高速风洞中获得更为详尽的气动力信息,在2．4m跨声速风洞中进行了连续变迎角测力试验技术应用研究。主要介绍了该项试验技术的基本特点,给出了J7标模的主要试验结果。结果表明,该项试验技术获得的气动力数据与常规阶梯方式具有很好的一致性,可以满足工程实用的要求。     

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

仿真结果检验和状态分析方法

提出了一种仿真结果检验和状态分析方法,不仅可以有效地判断仿真结果正确与否,而且当仿真结果不正确时,可以进一步对仿真结果进行状态分析,确定哪些状态的仿真结果正确,哪些状态的仿真结果错误,从而指导仿真软件系统的编制和修改.此外,还给出了仿真结果检验表和状态分析表,便于工程应用.      

某粉末冶金高温合金涡轮盘破裂转速分析

介绍了平均应力法、平均径向应力修正法、残余变形法、局部塑性应变法等破裂转速预测方法,并采用这些方法对某粉末冶金高温合金涡轮盘破裂转速进行了分析.通过理论预测与试验结果的对比得知:采用平均径向应力修正法、残余变形法、局部塑性应变法预测的破裂转速与试验结果吻合较好,而采用基于径向应力分布均匀假设的平均应力法推算的破裂转速与试验结果相差较大.      

确定性仿真结果检验和因子分析方法

提出一种确定性仿真结果(不考虑随机因素影响的仿真结果)检验方法, 能够对仿真结果与试验结果之间的误差进行分析, 确定该误差是由于仿真不正确引起的系统误差, 还是由于试验结果的分散性引起的偶然误差, 进而判断仿真结果正确与否.同时建立一种确定性仿真结果多因子分析方法, 能够对仿真结果的各个影响因素(因子)进行分析, 确定哪些因素已被正确仿真, 哪些因素尚未被正确模拟, 从而指导仿真软件的编制.      

航空管路补偿器耐久性计算方法对比分析

本文研究了管路补偿器补偿位移疲劳寿命和耐久振动寿命的有限元分析计算方法。通过对某型管路补偿器有限元仿真计算结果和传统工程计算法计算的结果进行对比分析,研究管路补偿器补偿位移疲劳寿命的有限元计算方法,探索耐久振动寿命的有限元计算流程,并将结果同试验数据对比,以验证其设计方法的计算精度。结果表明,补偿位移的疲劳寿命有限元法的计算精度优于工程设计法;耐久振动寿命有限元法能够给出确定的振动应力分布情况,并预测出耐久振动寿命时间,工程设计法仅能计算出管路补偿器的自振频率,不能明确振动应力和寿命。     

确定性仿真结果检验与状态分析方法

为判断多状态的仿真结果正确与否,建立了变异系数相等情况下的确定性仿真结果检验方法.同时,在变异系数相等或方差相等的情况下,分别建立了一种确定性仿真结果状态分析方法,可对仿真结果的各个状态逐一进行分析,以确定哪些状态已被正确仿真,哪些状态尚未被正确模拟,从而指导仿真软件的编制.