低温推进剂液体火箭发动机超导电磁泵压循环系统初步设想

为解决大型液体火箭发动机现有循环系统性能和技术难度及成本之间的矛盾，本研究基于低温推进剂液体火箭发动机工作环境特点，并结合磁流体发电机和直流超导电机的原理和工作特性，提出了一种不同于涡轮泵循环的超导电磁泵压循环系统。文中阐述了系统中超导电磁泵和磁流体发电装置的基本结构与工作原理，分析了系统的优缺点，并对超导电磁泵压循环系统进行了初步的可行性论证，最后对超导电磁泵压循环系统的应用前景进行了分析和论述。结论认为．作为一种全新的火箭发动机循环系统，超导电磁泵压循环系统具有非常广阔的发展及应用前景。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

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

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

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

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

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

挤压油膜阻尼器的均匀两相流模型

在挤压油膜阻尼器内气穴现象实验研究的基础上，对挤压油膜因气穴现象而发生物理特性的变化进行了讨论。气液均相流模型作为挤压油膜的近似描述，导出了上述模型的密度和粘度与油膜压力之间的关系式。定性地分析了当挤压油膜采用这一可压缩模型时，其压力分布的描述更能与实际接近，对阻尼器的力学特性的准确预测有重要的参考价值。     

超声速自由剪切层流动的数值模拟和理论分析

本文从二维完全Ｎ－Ｓ方程出发，采用ＮＮＤ２Ｍ差分格式，计算了超声速自由剪切层流场，并且研究了上游扰动在剪切层内的传播规律。在数值模拟的基础上对流动涡结构进行了分析。     

轴对称射流矢量喷管的试验和数值模拟

对基于激波来实现推力矢量的轴对称射流矢量喷管的缩比模型进行了测力和测压试验，用推广到可计算可压缩流的SIMPLE方法其内外流场进行了数值模拟，根据试验和数值模拟结果分析了喷管主流与次流相互作用产生的复杂流场结构、二次流流量和落压比对气动矢量角的影响，在落压比3～6范围内，二次流流量和喷管主流比值增大，气动矢量角增大，两者比值相同时，落压比增大，气动矢量角减小。     

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

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

表面磁流体气动激励控制楔面激波规律数值研究

针对表面磁流体(MHD)气动激励对高超声速进气道在非设计状态下激波控制问题,从唯象学的角度出发,基于低磁雷诺数假设,将电磁作用简化为Navier-Stokes方程组中的源项处理,同时考虑到低气压、低磁场环境下电子回旋效应引起的Hall效应,并联立Ohm定律,建立磁流体动力学模型,通过与实验纹影对比验证了模型的合理性,并利用该模型研究了表面MHD加/减速激励作用位置与宽度、磁场强度、电导率和能量转化率等参数对楔面激波的影响规律。结果表明:表面MHD气动激励包括焦耳热与洛伦兹力作用,当放电功率密度为3.8×10~9k W/m~3,磁场强度为0.34T时,MHD加/减速激励分别使激波位置前移6mm与10mm;而磁场强度较低时,由于焦耳热的主导作用,将会出现激波前后压力比增大,激波强度增加等负面效应;根据激励参数的影响规律,激励器电极应靠近尖端布置,增大磁场强度,并改善等离子体源以提高气体电导率,同时适当增大激励区域宽度。     

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

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