三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

基于泊松方程的二元翼型网格生成技术

利用求解Posson 方程的方法来划分计算网格，使用Hilgenstock 方法来 控制网格间距和网格正交度。对于安装在风洞中的二元翼型，在进行Navier－Stokes 方程 计算时，需要在物面附近的网格间距很小，这导致网格划分时计算收敛速度很慢。 本文对 Hilgenstock 方法进行了改进，大大提高了网格划分的计算收敛速度。     

三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

一种适合全流速的流场数值方法

发展一种适合于所有流速的压力修正方法,在考虑了密度的可压缩性后,压力修正方程夸为对流扩散模型的,从而与般标量一样使用二阶中心差分离散。整个算法具有二阶精度,进一步使用多步压力修正来增强网格非正交时的收敛性态,在对流项离散中,设计了独特的人工耗散,也同样适用于所有流速,利用本算法求解全三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,计算采用了非交错网格方案,使用Ｒｈｉｅ－Ｃｈｏｗ方法消除了速度－压力的不关联,     

求解N-S方程SIMPLE算法的半交错网格研究

运用半交错网格设置对不可压缩流动的Navier-Stokes方程进行数值离散求解，文中对不合理压力场产生的原因及抑制方法进行了探讨，导出了压力控制体界面上速度的修正表达式。作为算例．本文采用区域湍流模化方法，在核心区采用标准k-ε模型，在近壁区采用单方程低雷诺数模型．对封闭转-静盘腔、具有径向出流的转-静盘腔和非规则环形通道内的流动进行了数值模拟．并与实验数据进行了对比分析，结果表明本文的计算方法是合理的。     

直升机旋翼翼型及桨叶气动外形反设计分析

建立了一个基于余量修正思想的直升机旋翼反设计计算方法,用于直升机桨叶的气动外形设计研究.使用Poisson方程为控制方程生成围绕桨叶的贴体网格;在悬停状态下建立了以Euler方程为主控方程的旋翼流场求解方法,并采用了嵌套网格方法进行数值计算;在流场计算及网格生成基础上,采用MGM方程作为翼型反设计方程,建立了一套直升机旋翼翼型及桨叶气动外形的反设计方法.应用该方法,分别对二维翼型以及悬停状态下的旋翼桨叶进行了反设计分析.反设计结果表明,在给定的目标压力分布条件下,使用本文方法分别获得了满足要求的二维翼型及直升机桨叶外形,并与目标压力吻合良好.     

侧壁约束效应对三维方腔自持振荡和噪声辐射影响的实验研究

实验研究了侧壁约束效应对三维方腔流动结构和噪声辐射特性的影响,固定方腔长深比为2∶1,使用麦克风阵列测量了方腔宽长比从0.1变化至0.5过程中流致噪声在不同指向性下的强度变化规律,并使用脉动压力传感器测量了不同宽长比方腔内部壁面压力分布,同时采用TR-PIV（Time-Resolved Particle Image Velocimetry）测量了方腔内流动结构的发展。实验结果表明:对于宽长比为0.5的方腔,当来流马赫数大于0.03时,方腔流动开始出现振荡并向上游辐射噪声;当来流马赫数增大至0.20时,方腔流动发展为对应Rossiter三阶模态的自持振荡,并辐射出尖频噪声。减小方腔宽度,当宽长比小于0.3时,方腔流动的自持振荡和尖频噪声被大幅度抑制甚至消除,来流马赫数为0.20和0.25时,方腔上游总声压级能够降低3 dB以上。通过对比壁面压力分布和PIV流场测量结果,发现减小方腔宽长比时,方腔内主回流涡向上游移动,涡强度降低,使得方腔的流动反馈不足以形成自持振荡,从而降低了辐射噪声。     

处理动边界问题的无网格/直角网格混合算法

把无网格/直角网格混合算法发展用于涉及动边界的非定常流动计算问题.采用无网格/直角网格混合算法,计算区域整体采用直角网格,只在物面附近嵌入局部无网格.对于动边界问题,本文让主体直角网格保持固定,用原混合算法中的局部无网格区跟随运动边界,使得动边界物面附近有不变的点云结构,避免了点云重构;无网格区与直角网格交界面随边界运动跟踪确定,由于主体基于直角网格,方便了交界面的跟踪.用发展的混合算法,结合非定常Euler方程双时间推进方法,先对翼型作俯仰振荡的非定常跨声速流动进行了数值模拟,计算结果与现有文献和实验数据进行了比较与验证;结合刚体运动方程,成功模拟了二维外挂物投放过程,算例展示出发展的混合算法适合处理大位移运动边界问题.     

基于旋转坐标系的高阶间断有限元方法非定常湍流数值模拟

基于高阶间断有限元方法(Discontinuous Galerkin method,DGM),对旋转非惯性系下耦合了修正的一方程S-A模型的RANS方程进行了离散求解。为了在稀疏网格上获得更贴近真实的物面形状,使用了多层高阶弯曲网格方法对物面进行拟合。非定常时间推进采用了隐式双时间步方法,每个时间步产生的线性系统采用预处理的方法,即广义最小残差方法(Generalized minimal residual method,GMRES)来求解。计算了旋转圆柱绕流以及经典翼型振荡算例的升力和力矩迟滞曲线,与实验结果以及前人的计算结果对比验证了本文方法的正确性和有效性。     

模型燃烧室湍流亚网格尺度模型

分别采用Smagorinsky-Lilly和k方程亚网格尺度模型,在任意曲线坐标系下对模型燃烧室湍流流场进行大涡模拟.分别用QUICK和二阶隐式差分格式离散控制方程及其时间项,并用SIMPLE算法进行求解.计算结果与PIV测量数据比较表明:大涡模拟方法能较好模拟旋流器后面回流区的瞬态变化,两种亚网格尺度模型都能较好模拟湍流流动,但k方程亚网格尺度模型稍优于Smagorinsky-Lilly亚网格尺度模型,因此采用k方程亚网格模型更适用于大涡模拟实际燃烧室复杂湍流流场.     

平面扩压叶栅流场犘犐犞与三孔尾迹探针对比测试研究

针对PIV技术在暂冲式高亚声速平面叶栅流场中遇到的示踪粒子投放问题,通过采用高压雾化式粒子发生器以及安装在稳定段前的撒播器,有效地使示踪粒子均匀地与主流混合,并成功开展了某扩压平面叶栅叶片槽道及出口尾迹可视化测量,获得了零迎角、进口马赫数0.2~0.8状态下的二维速度矢量场。为了验证PIV在叶栅流场测试结果的可靠性,在相同工况下,将PIV测量结果分别与数值计算结果和三孔尾迹探针测量结果进行了对比分析。结果表明:采用PIV技术测得的叶栅中截面二维速度矢量场合理地反映了叶片槽道及尾迹的流动结构,与数值模拟结果较为接近;PIV与三孔楔形尾迹探针在叶栅出口尾迹的测量所获得的气流速度和主流区的出口气流角重合性较好;尾迹分离区的出口气流角重合性略差,主要原因是尾迹区气流角超出了探针校准范围,这也说明了PIV测试技术优势。本文提出的PIV测量技术也可用于连续式叶栅风洞中。     

矩形隔离段内激波串结构及其影响因素分析

隔离段是超燃冲压发动机的重要组成部分,主要作用是隔绝燃烧与进气道的相互干扰。隔离段中存在的复杂流动现象一直是人们研究和关注的重点。利用三维数值模拟方法对矩形隔离段激波串特性影响因素进行了研究,主要分析了不同来流马赫数、单侧和对称扩张角以及壁面凹腔等因素影响下的激波串特性。结果表明:在高来流马赫数条件下,隔离段内激波串长度变短,隔离段抗反压能力增强,总压损失增大;在单侧和对称扩张隔离段内的激波串结构存在差异,且隔离段后的流场总压损失与扩张形式无关;隔离段添加壁面凹腔后,在不同反压下会出现2种模态（亚临界凹腔模态和超临界凹腔模态）,2种模态下隔离段内激波串结构及流场参数特性有所不同,超临界凹腔模态下隔离段抗反压能力下降,总压损失增大。本文的研究结果可为隔离段和燃烧室设计及试验提供参考。     

非线性气体振荡整流效应在矩形机翼上的应用研究

在对圆管内气体的非线性振荡进行理论和实验研究的基础上,利用开口圆管非线性气体振荡的整流效应,通过流动显示及压力测量,对翼面分离流控制进行实验研究。研究结果表明,翼面开缝引入气体振荡有利于改善大迎角下翼面的压力分布,提高升力,增大机翼的稳定性,同时,不论采取何种开缝模式,这种输入能量的方式对提高其气动性能都有一定作用。     

跨声速翼型风洞开孔壁的简化模型数值模拟

以开孔壁翼型风洞为研究对象,构建简化的仿真模型,模拟开风洞孔壁附近及小孔内的流动,研究开孔壁对风洞试验的影响。研究了开孔壁流动的主要特征参数并建立多孔板模型,为进一步建立数值风洞模型及研究洞壁干扰提供参考。通过简化的孔壁模型,研究了开闭比等特征参数对风洞流场和翼型绕流的影响。构建了二维简化孔壁模型和多孔介质孔壁模型,并验证了以多孔板模型模拟孔壁风洞流场的可行性。本文建立了一种研究跨声速孔壁风洞的孔壁效应的数值方法,为跨声速孔壁风洞流场的模拟研究提供参考,为进一步构建可靠的风洞孔壁数值模拟数学模型提供一种研究思路。     

带辅助进气的二元超声速进气道低速特性

通过飞行试验,研究了带辅助进气门的二元超声速进气道在亚声速低马赫数飞行时的内部流场分布及进气畸变等特点。结果表明,辅助进气流与进气道内主流掺混会在进气道内辅助进气门下游区域产生低总压区,引起畸变增加、总压恢复下降,并且进气量及引起的畸变随发动机转速的增大而增大,随马赫数的增加而减小,随飞行高度的变化则无显著差异;同时,侧滑角向左及向右增大时辅助进气产生的低压区范围以及进气畸变也会增大,而迎角变化时进气畸变及总压恢复变化不明显。另外,通过分析各试验点进气道出口低总压区的变化与流量系数的关联,确定了该型进气道辅助进气门打开及关闭状态对应的工作范围。     

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究.采用“双时间步”方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化.结果表明:盲腔构型在高超声速飞行中会出现周期性的激波振荡现象.数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好.     

某组合发动机进气道抽吸作用分析

对某组合发动机进气道设计马赫数下的三维流场进行了数值模拟。研究了不同位置抽吸对进气道流场及性能的影响,研究结果表明,外压部分抽吸能在较小的抽吸流量下有效改善进气道的起动性能和气动性能,而又不改变进气道的内通道特性;内通道抽吸均化了喉道出口流场,提高了总压恢复系数,但减弱了压缩程度,造成压比、温升比较大幅度地下降,喉道出口马赫数增加,且抽吸流量越大,压比、温升比的降幅和喉道出口马赫数的升幅越大。     

含化学反应流动的TVD格式及其在燃气射流中的应用

以Harten标准TVD格式为基础，结合固体火箭燃气射流的特点，以数学方法系统地推导了适合于高温、高压和高速流体流动的数值格式，给出了压力偏导数的合理计算公式；利用化学动力学知识，对火箭燃气射流流场中存在的多组分、含有限速率的化学反应问题进行了论述，阐明了化学反应质量源项的求解方法。以12组分9反应方程模型为例，利用编制的计算程序，对某火箭燃气自由射流流场进行了模拟。通过对结果的分析，肯定了数值格式的正确性。     

运输飞机非常态俯仰振荡分析

对某运输飞机在大襟翼偏度下的非常态俯仰振荡原因进行了分析。研究结果表明,该运输飞机在大襟翼偏度下飞行时因平尾受机翼的强下洗作用导致平尾下翼面发生局部分离,该分离未达到失速分离,从而引起飞机非常态的俯仰振荡。分离的主要原因是机身与平尾融合处抗流场逆压梯度能力下降,导致在大襟翼偏度飞行时出现局部分离,使得飞机出现纵向非常态俯仰振荡。该结论在风洞试验及该运输飞机试飞中得到验证。     

非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。