高超音速粘性流的快速算法

由 Lombard 等人提出的 CSCM-S 算法是求解 Euler 及 Navier-Stokes 方程的高效率方法。但是,由于该法在整个计算区域上采用总体迭代,实用上是不经济的。本文提出一种修正方案,即对超音速流为主的区域用一次推进解法与复杂流动区域用多次扫描,可大大节省 CPV 时间,更适于工程应用。     

张量方程在三维粘性绕流计算中的应用

本文采用张量形式的粘性激波层方程,用空间推进的数值方法计算了球锥、椭球锥有攻角高超音速绕流问题,并计算了组合椭球锥的粘性绕流,从而说明了本文的方法可推广应用于一般外型飞行器的小攻角绕流计算问题。文中考虑了在高超音速流动条件下空气的非平衡化学反应,认为化学反应的速率是有限的     

粘性金属中的正激波稳定性问题

G.H.Miller等把高压金属中的粘性激波作为强间断面处理,解析推论出:在大粘性系数条件下小扰动激波是不稳定的,物质粘性是导致失稳的因素。本文中针对平面正激波,认为高压金属中的粘性激波的物理量是连续变化的,利用线性稳定性理论,用数值解推论出:在有粘性条件下小扰动激波都是稳定的,物质粘性是致稳的因素。指出G.H.Miller等获得错误结论的原因在于:从无粘流动解推出的小扰动边界条件导致粘性激波小扰动增长。给出实验确定的小扰动速度梯度的边界条件,这样既可以把粘性正激波作为强间断面处理,也能够保证粘性正激波的稳定性。     

激波管氢氧爆轰驱动技术的发展进程

回顾了爆轰驱动的出现与认识过程及对其性能的观察结果. 采用卸爆管消除爆轰波反射高压和利用变截面缓解Taylor波有害影响, 使得反向和前向爆轰驱动模式能用来产生高焓高压试验气源. 双爆轰驱动段彻底消除了爆轰波后的Taylor稀疏波, 不仅提高了前向爆轰的驱动品质, 而且为进一步提高驱动能力开辟了新途径.      

二维激波与剪切层相互作用的直接数值模拟研究

采用五阶weighed esseritially non-oscillatory (WENO) 格式和三阶total variation diminishing (TVD) Runge-Kutta 格式, 通过求解二维非定常Navier-Stokes 方程, 直接数值模拟了激波与剪切层相互作用, 目的在于揭示激波与剪切层相互作用过程中噪声产生的机理. 研究发现:(1) 当入射激波穿过剪切层时, 剪切层中心位置向下层区域偏移;(2) 入射激波穿过剪切层产生小激波, 在小激波与剪切层接触点处产生声波并向外辐射;(3) 反射激波穿过剪切层后形成了分段弧状激波;(4) 当反射激波穿过剪切层时, 激波在鞍点处泄漏并向外辐射声波, 这是一种激波泄漏机制.     

机动弹头高超声速粘性绕流数值模拟

采用空间二阶精度的交替方向隐式分解的NND格式求解完全气体假定下的非定常薄层近似NaVierStokes方程,并采用抛物化的椭圆型方程生成复杂带翼弹头的空间网格.最后给出了机动弹头在M∞=7.3,α=20°下的流场计算结果.计算结果表明本文所采用的计算网格及计算方法可以适用于高超声速复杂流场的数值模拟.     

两相正激波松弛流的数值研究

对气、固两相正激波松弛流动进行了数值研究,方程中考虑了压力梯度及非匀速运动引起的附加质量力。给出了完整的计算方法,对三种不同固相容积比、不同固相颗粒直径和不同马赫数分别进行了数值计算,所得计算结果给出了气、固两相的速度、温度、压力、密度、固相容积比及熵值沿流向分布情况。文章对气、固两相参数在松弛区中的变化规律进行了讨论。     

关于超音速流中运动激波若干问题的研究

研究了在无粘完全气体流中的运动激波 ,讨论了激波运动速度D和来流速度U对激波后气流参数的影响 ,包括对激波后的总焓比值和总压比值以及对流转角的影响。计算结果表明它们不同于通常静止激波下所得到的结果。该内容涉及到超音速射流与障碍物或空腔体相互作用时出现的失稳状态下激波的振动和空腔体底部的反常加热问题。     

粘性指进显示及其分维测量

考虑到液体驱油时粘性指进可能造成原油回采严重减少，利用图像采集系统对多孔介质中由粘性指进不稳定性导致的指进花样进行显示，测算其分维指数并研究聚合物浓度与分维的关系。     

两相正激波松弛流的数值研究

对气、固两相正激波松弛流动进行了数值研究,方程中考虑了压力梯度及非匀速运动引起的附加质量力。给出了完整的计算方法,对三种不同固相容积比、不同固相颗粒直径和不同马赫数分别进行了数值计算,所得计算结果给出了气、固两相的速度、温度、压力、密度、固相容积比及熵值沿流向分布情况。文章对气、固两相参数在松弛区中的变化规律进行了讨论。     

分层流体中gKdV型孤立波的迎撞

本文采用约化摄动法和PLK方法并通过双参数摄动展开,讨论了分层流体中以推广的Korteweg-de vries方程(gKdV方程)描述的孤立波的迎撞问题,求得了二阶近似解。分析结果表明,gKdV型孤立波碰撞后保持原来的形状不变,在碰撞时最大波幅为两个来碰孤立波的最大波幅的线性叠加。     

泰氟龙烧蚀对黏性激波层电子密度的影响

采用黏性激波层基本方程组对有无泰氟龙烧蚀两种情况下的钝锥体化学非平衡绕流作了数值求解,以研究泰氟龙烧蚀对流场电子密度的影响规律.算例结果表明:泰氟龙烧蚀确有降低流场电子密度的效应.平衡催化壁工况下这种效应的强度,显著大于非催化壁工况下的强度;远下游截面处的这种效应的强度,显著大于驻点区的强度.此外,通过计算结果分析,对形成上述影响规律的原因作了初步讨论.     

冲击波阵面上粘度的一种唯象描述

本文从粘度的基本定义出发,导出了冲击波阵面上粘度系数、应变率和粘应力的一种唯象表达式。并以45~#钢和L Y-12铝材料实验测量的冲击波剖面为例,计算出在这些材料中冲击波阵面上的上述参量及平均粘度系数。所得结果与W.C.Moss和张若棋等给出的结果进行了比较和讨论。     

LOCAL BENDING OF THIN FILM ON VISCOUS LAYER

Effects of deposition layer position film are systematically investigated. Because the and number/density on local bending of a thin deposition layer interacts with the thin film at the interface and there is an offset between the thin film neutral surface and the interface, the deposition layer generates not only axial stress but also bending moment. The bending moment induces an instant out-of-plane deflection of the thin film, which may or may not cause the socalled local bending. The deposition layer is modeled as a local stressor, whose location and density are demonstrated to be vital to the occurrence of local bending. The thin film rests on a viscous layer, which is governed by the Navier-Stokes equation and behaves like an elastic foundation to exert transverse forces on the thin film. The unknown feature of the axial constraint force makes the governing equation highly nonlinear even for the small deflection chse. The constraint force and film transverse deflection are solved iteratively through the governing equation and the displacement constraint equation of immovable edges. This research shows that in some special cases, the deposition density increase does not necessarily reduce the local bending. By comparing the thin film deflections of different deposition numbers and positions, we also present the guideline of strengthening or suppressing the local bending.     

在超声速流动中激波与边界层的干扰特性

激波与物面边界层的干扰涉及可压缩流动的稳定性、转捩、分离等问题,直接影响到飞行器的阻力、表面热防护和飞行性能等工程技术问题。首先总结了前人对于激波与边界层的干扰所做的工作,之后重点研究和对比分析了超声速与跨声速流动中,正激波、斜激波以及头部激波对于飞行器层流和湍流边界层的干扰影响。激波强度的不同对边界层干扰作用不同,在强干扰情况下将会引起边界层分离和翼型失速。     

INTERFERENCE OF SHOCK AND BOUNDARY LAYER IN SUPERSONIC FLOW

激波与物面边界层的干扰涉及可压缩流动的稳定性、转捩、分离等问题,直接影响到飞行器的阻力、表面热防护和飞行性能等工程技术问题。首先总结了前人对于激波与边界层的干扰所做的工作,之后重点研究和对比分析了超声速与跨声速流动中,正激波、斜激波以及头部激波对于飞行器层流和湍流边界层的干扰影响。激波强度的不同对边界层干扰作用不同,在强干扰情况下将会引起边界层分离和翼型失速。     

非惯性坐标系中的边界层方程

本文用匹配渐近展开法导出二维翼型在非惯性坐标系中的一阶和二阶边界层方程,消去压强项后,一阶边界层方程与经典边界层方稃相同;惯性力的作用在二阶边界层中才出现。     

爆轰驱动高焓激波膨胀管性能研究

介绍了中国科学院高温气体动力学实验室JF-16爆轰驱动高焓激波膨胀管的研制进展及其性能测试结果. 性能研究主要测量了入射激波速度和压力曲线. 根据入射激波速度,应用Mirels的黏性激波管理论计算了超高速试验气流的运动速度,然后应用Gaseq软件计算了试验气流的热力学参数. 研究结果表明:在16.35\,m长的JF-16激波膨胀管中获得了流速7\,000\,m/s以上的稳定试验气流,试验时间为50$\sim$100$\mu$s, 试验气流总焓30\,MJ/kg左右.      

激波卷扬的气体粉尘边界层的实验与数值研究

采用瞬态阴影技术和激光消光技术对激波卷扬粉尘的现象进行了实验研究，理论研究采用ＰＳＩＣ（ＰａｒｔｉｃｌｅＳｏｕｒｃｅｉｎＣｅｌｌ）模型，模型同时考虑了湍流效应及气固两相间的输运效应。结果表明，激波卷扬的粉尘颗粒在壁面附近湍流区内先被加速，一定高度后减速，从而在边界层内形成一个高颗粒浓度的区域。数值结果与实验结果基本一致。