火箭垂直回收中发动机布局与喷流壁面效应影响研究

针对重复使用火箭垂直着陆过程的喷流流场问题开展研究,利用计算流体力学(Computational Fluid Dynamics, CFD)方法研究了壁面效应和发动机布局对超声速喷流的影响。研究表明,着陆距离(L)在2.24D～11.2D(D为喷管出口的直径)的范围内,地面效应对喷管出口中心处的温度分布影响较小;在当前计算条件下,当L<2.24D时,超声速喷流撞击地面会形成强烈的激波,随着离地高度的降低,该激波位置往喉部方向移动,由于壁面效应,喷管内部形成斜激波,导致中心喷管壁面处的温度升高;中心喷管相对外侧喷管往外突出增大了壁面流动速度,导致外侧喷管出口的温度降低;研究还表明子级火箭底部端面的喷管数量增加后,会导致喷管的温度升高。研究结果将为火箭发射及回收方案选取提供参考。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

平行平板内连续布置肋块流场与传热的数值模拟

为了研究平行平板内连续布置8个肋块的流场结构和传热现象, 利用计算流体力学软件FLU-ENT6.2内置的标准k-ε湍流模型、RNGk-ε湍流模型和Realizablek-ε湍流模型对此进行了数值模拟并对结果作了比较.流场结果显示在相同来流条件下, 由Realizablek-ε模型预测的肋块间逆时针方向回流区最大, 标准k-ε模型最小.数值结果表明从第2个肋块开始加热段表面压力系数Cp和努塞尔数Nu开始呈周期性变化.同时将第7个肋表面的努塞尔数数值解与实验值进行了比较, 结果为三种湍流模型都能准确模拟出努塞尔数的变化趋势.计算结果可为工程设计和使用提供参考.      

湍流模型对叶片表面换热计算的影响

为了比较湍流模型对涡轮叶片外换热计算结果的影响,采用五种湍流模型及两种壁面函数计算了NASA-MarkII导向叶片表面换热并与实验数据进行了对比,结果表明:在相同边界条件下不同湍流模型的计算结果有很明显的差别,即使是同一种湍流模型,如果采用不同的壁面处理函数其计算结果也是有很大差别的;某些湍流模型的计算值只是在某个区域较为理想,还不能找到在整个叶片表面计算结果与实验数据较为接近的湍流模型.在尚无普适性较好的湍流模型的情况下,研究在不同计算域采用不同湍流模型和壁面处理函数的计算技术,是一种较好的可行的方案.      

温度边界层对壁湍流多尺度相干结构的影响

用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中壁面加热平板湍流边界层不同法向位置流向速度分量的时间序列信号, 与常温情况相比, 热对流加速了壁湍流中流体的动量、质量和能量的交换, 加快了湍流边界层的发展;用子波分析方法研究了壁面加热对壁湍流多尺度相干结构喷射和扫掠条件相位平均波形的影响, 显示热对流明显增大了缓冲层相干结构扫掠过程和喷射过程的强度, 是增强壁湍流中流体的动量、质量和能量交换的根本原因.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

Flight Flutter Modal Parameters Identification with Atmospheric Turbulence Excitation Based on Wavelet Transformation

针对基于大气紊流激励的飞机颤振试验数据具有信噪比低以及测量数据以加速度响应形式给出等特点,将小波变换与随机减量法相结合对飞机颤振模态参数进行识别,首先利用随机减量法获得结构在非零初始加速度条件下的自由衰减响应,然后对该自由衰减响应进行Morlet连续小波变换,为使变换简单,本文采用了Parseval定理和留数定理,利用Morlet的带通滤波特性,通过搜索小波变换系数在不同的尺度区域内的极大值,再根据小波变换系数的极大值、相角与固有振动频率及阻尼系数间的关系,即可识别出各模态参数。并对利用小波变换识别颤振模态参数时,实现模态解耦应满足的条件进行了分析研究。通过仿真和飞机颤振试验数据分析,验证了该识别方法不仅简单、有效和可行,而且具有较强的抗噪性。     

平板湍流摩阻系数的显式表达式

 湍流摩阻系数的计算是工程技术中经常遇到的问题。计算时所依据的公式多为半经验半理论的结果,如Karman公式,Van Driest公式和Prandtl公式等。这些公式是以隐函数给出的,给计算带来不便。但是不难发现,这类公式经过适当的变换,可以化为统一形式的超越代数方程,这种方程没有有限形式的解析解,本文给出了一种近似解析解。根据这种方法除了导出上述公式的显式表达式外,还对White公式作了修正。当然解超越代数方程最根本的方法是迭代法,本文对此也作了简单讨论。     

对称翼型近场尾流及后缘边界层流动特性的试验研究

本文简述了用热线风速仪测量“NACA 63-012”对称翼型近场尾流及后缘边界层流动特性,讨论了流场的时均特性及雷诺应力等湍流量的分布特性,提出了无因次速度分布的半经验公式、相关参数及相关曲线解析式,并得出了几点重要结论。     

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。