亚、超音喷管中声传播的数值模拟

采用计算气动声学方法研究了管道近音速区的声传播和声波与激波的相互作用问题。空间离散采用高阶DRP格式,时间离散采用低耗散低频散的龙格库塔方法,边界条件采用Giles的特征变量一维非定常无反射边界条件。计算结果和解析解的精确吻合表明,计算中应用的这种高阶的DRP格式有很高的分辨率和很好的稳定性,适合声传播的计算。      

基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

阵风与叶栅干涉噪声的数值模拟

采用计算气动声学方法求解二维线化欧拉方程对阵风与叶栅的干涉噪声进行了数值模拟。为减少频散和耗散误差,采用七点四阶频散相关保持有限差分格式进行空间离散;时间推进采用优化的4/6低频散低耗散龙格库塔法;在各边界处均选择或建立了适当的边界条件。首先通过模拟阵风在自由空间中的传播验证了无反射边界条件的可靠性;然后对阵风与平板叶栅的干涉问题进行了模拟,数值结果与精确解符合得很好;最后给出了阵风与NACA-0012翼型叶栅干涉的数值模拟结果,显示了数值求解方法的可行性。      

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

声波与剪切流作用的数值模拟

基于完整的可压缩流体力学方程 ,利用特征分析的方法给出了气动声学中数值计算扰动的边界处理方法。采用所提出的边界处理方法 ,模拟了静止流场、不同来流速度的均匀流场及均匀剪切流场中扰动的演化 ,结果表明 :给出的边界条件的处理方法是合适的。文中还模拟了混合层中扰动的演化 ,给出了扰动在混合层中传播的特征      

高阶格式及其在内外流场计算中的应用

高阶格式这里包含两种不同的含义,一种是人们常说的高精度格式;另一种是高分辨率格式,因为它常与高精度有密切联系,但两者毕竟有严格的区别.在格式的构建上,相比之下高分辨率格式较高精度格式难得多.本文从4个方面讨论了高精度格式与高分辨率格式在相关领域中的应用并给出了大量算例,这对澄清人们对高阶格式的认识,弄清它们各自的适用范围是十分必要的.大量的数值计算与研究表明:对于复杂流场,尤其是既考虑激波与边界层相互作用,又要考虑壁面传热的流场,用低阶数值格式计算是无法得到与实验较接近的数值结果的,因此发展高阶格式对处理高温涡轮的气膜冷却问题以及高速飞行器的气动热计算是至关重要的.      

非均匀横流作用下冲击射流冷却的数值研究

运用数值计算的方法对非均匀横流作用下冲击射流冷却的流动和传热过程进行了三维数值研究,揭示出射流冲击孔和横流孔相对排列方式、射流雷诺数和横流雷诺数相对大小、射流孔中心至横流孔出口截面距离等因素对冲击射流冷却流动传热特性的影响。结果表明:若保持总的冷却流量不变,当横流雷诺数与射流雷诺数之比等于或大于0.5后,有横流时的冲击射流冷却局部努塞尔数比无横流时有较为显著的下降;在非均匀的横流作用下,对于叉排和顺排方式,由于射流和横流的相互作用,壁面温度分布呈现出复杂的特征;射流孔和横流孔呈顺排布置时,峰值努塞尔数所对应的射流驻点区向下游偏移。      

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

自适应直角坐标网格方法在对流换热数值模拟中的应用

复杂几何边界是对流换热计算面临的主要难点之一。自适应直角坐标网格方法可以比较方便的处理复杂几何边界,但其应用目前主要集中于可压无粘流的Euler方程、可压粘性流的Navier-Stokes方程的外部绕流计算。本文应用自适应直角坐标网格方法计算了两个不可压内流对流换热算例,并把计算结果同商用CFD软件STAR-CD的计算结果进行了对比。自适应直角坐标网格方法在本文的两个算例中得到了同STAR-CD一致的结果。这说明采用该方法数值模拟对流换热问题是适用的。      

跨声压气机周向槽处理机匣非设计工况研究

设计了某单级跨声压气机周向槽处理机匣,利用全三维定常数值模拟方法研究了三种典型转速下的机匣扩稳机理;对比分析得出处理机匣扩稳性能与压气机工况、处理槽位置的关系:转速越高,处理槽扩稳效果越好,而且位于叶尖弦长中部的处理槽更有效.当转速降低时,扩稳效果减弱,靠近尾沿的处理槽更为明显.研究表明,该压气机在降转速过程中存在一个扩稳机理的转变:在高转速下扩稳以抑制间隙涡破裂为主,而低转速下则以抑制叶尖吸力面附面层分离为主.      

特征无反射边界条件应用于二维圆柱绕流模拟

采用特征无反射边界条件和吸收层技术相结合的边界处理方法,研究了二维圆柱绕流模拟的边界反射问题。近年来在计算气动声学中提出了将大涡模拟和声学传播方程结合的混合方法求解声场,基于计算量的考虑,须在较小的计算区域内进行大涡模拟,这就需要在计算区域较小时,边界的反射尽可能小。作为大涡模拟研究的技术储备,将该方法应用于二维圆柱绕流的模拟中。数值结果表明这种方法在计算区域较小时,能起到很好的无反射效果,计算结果与前人在计算区域取得很大时的结果比较吻合。     

二维平行剪切层声波产生和辐射的数值模拟

 在二维线化Euler 方程的基础上解决了第3 届计算气动声学专题讨论会的第5 类标准问题: 二维平行剪切层声波的产生和辐射。频散相关保持有限差分格式用于空间离散, 低频散低耗散的龙格库塔法用于时间积分。重点使用了远场无反射边界条件。计算结果与NASA GLENN 研究中心的解析解符合得很好。     

平板附面层等离子体流动控制的数值模拟

通过求解电场中的拉普拉斯方程和赫姆霍兹方程的变形形式,成功地将等离子体激励对平板附面层流动的影响,以体积力向量的形式引入到NS方程之中.借助求解的电荷密度,通过将计算结果与实验数据的对比,可确定德拜长度、最大电荷密度、形状因子等可调参数的取值原则,建立起平板附面层等离子体流动控制的数值模拟方法,为将等离子体流动控制方法应用于外流及内流场中的强剪切流动控制,奠定关键性的技术支撑.      

抽吸对高超声速进气道抗反压能力的影响

对不同抽吸开孔率下某典型高超声速进气道进行了二维数值模拟,简要分析了压比变化对隔离段内部流场的影响,在该分析的基础上讨论了附面层抽吸对进气道最大抗反压能力的影响,给出了最大抗反压能力随抽吸流量之间的变化规律,分析了采用附面层抽吸技术提高进气道抗反压能力的原因.分析结果表明:采用附面层抽吸技术能够提高进气道的抗反压能力,相当于增加了隔离段的长度.抽吸孔面积越大,相对抽吸流量越大,进气道的抗反压能力越大.