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

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

三维可压缩混合层中扰动演化的研究

为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.     

计算气动光学研究进展

计算气动光学是采用数值计算的方法研究高速空气动力流场对光波传输和光学成像影响的一门交叉学科,在红外成像制导设计等工程应用中发挥着重要的作用。总结了计算气动光学研究的三种基本计算方法,包括基于CFD简化方法的光学统计估算、基于RANS流场的计算和基于LES/DNS瞬态流场的计算,对这三种计算方法的基本理论、计算过程进行了介绍,综述了新近的一些计算气动光学研究进展,最后指出了未来的研究重点和方向。     

一个新的可压缩性修正的k-ε模型

考虑结构可压缩性修正的影响,发展了一个同时考虑结构可压缩性修正和膨胀可压缩性修正的k-ε湍流模型,新模型包括Chang可实现性、Heinz湍流动能产生项以及Sarkar可压缩性三部分修正.新模型扩宽了以往发展的可压缩性修正模型的适用范围,适用于高超声速(M>5)复杂湍流流动中.通过对多个复杂超声速横侧射流工况的计算,验证了新模型的预测效果.与实验结果相比表明,几个工况下新模型的预测精度都显著高于标准k-ε模型.流体分离强度越大,新模型的修正效果越显著.与标准k-ε模型相比,新模型计算结果与实验更加接近.     

超声速混合层转捩过程实验研究的进展

超声速湍流机理的实验研究是一件十分困难的工作.在2000年以来,本研究小组在低噪声超声速混合层风洞研究、超声速流动精细结构测量技术研究方面取得了重要进展,这给超声速混合层湍流精细结构的研究奠定了基础.为了研究超声速混合层及其气动光学问题,在研制的超声速混合层风洞中,主要以基于纳米技术的平面激光散射技术(Nano-trace Planar Laser Scattering,简称NPLS)为基础,研究了几种对流马赫数的超声速混合层从层流到湍流转捩过程K-H不稳定涡的空间结构,以及K-H不稳定涡的空间结构随着时间的发展过程.实验结果清晰地反映了湍流混合的不稳定性与转捩的精细结构,以及转捩过程的展向精细结构.     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

超声速混合层在强迫振动下流场结构的实验研究

通过可控振动系统激励双喷管后的薄平板,采用基于纳米粒子的平面激光散射技术(NPLS)观测超声速混合层流场结构。控制薄平板的频率和振幅,结合边缘检测技术研究了超声速混合层大尺度结构的发展和演化规律。实验结果表明:强迫振动对超声速混合层流场结构以及混合效率有重要的影响。与无振动相比,强迫振动增强了混合层的三维特性,混合层的失稳位置提前;同时,强迫振动促进了上下两层流动的掺混,高频振动下混合层增长速度提高约50%,有效增强了混合。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。     

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

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

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。