二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

一种绕曲面物体的可压缩湍流混合格式(英文)

提出了一种中心-迎风型混合格式。在该混合格式中,中心差分格式和Roe通量差分裂格式进行混合,它们之间的切换通过一个二进制开关函数实现。为了验证该混合格式在计算绕曲面物体可压缩湍流问题时的可靠性,尤其是带激波的流动问题,采用分离涡模拟方法计算了3个典型的问题。研究结果表明,当前数值结果与已有的实验数据相符较好,这说明该混合格式可以用来研究带激波和湍流的曲面物体可压缩绕流问题。     

Stiffened状态方程下γ-model方法跟踪二维多介质流动

在Stiffened状态方程下，运用Level—set方程跟踪界面运动变化，把界面捕捉的等效方程、Level-set函数和欧拉方程组耦合，求解包含比热比和材料参数的耦合形式的流体力学方程组。计算方法采用二阶精度Wave-propagation算法，通过对气体-液体两种流体的一维Riemann问题、激波和气泡相互作用以及气液两种流体Richtmyer—Meshkov界面不稳定性问题进行数值模拟，抑制了不同介质界面两侧的非物理振荡。     

激波冲击下Air/SF_6斜界面不稳定性实验研究

激波在不同密度介质上的交互作用在可压缩湍流上具有重要的基础价值。激波在界面上的作用会引起Richtmyer-Meshkov不稳定性。激波不正规折射时,流场存在更多复杂的涡。研究马赫数为1.23、1.41的激波在初始倾角β=60°的Air/SF6界面上非正规折射的情况。入射激波的切向冲击和法向冲击的相互作用,在界面处产生涡,折射波在壁面发生马赫反射。利用阴影显示技术,给出了界面演化和混合的过程。     

基于矩阵误差的流动自适应求解

为了构造一个网格沿流场特征方向进行分布,以保证整个流场的误差尽量均匀化的网格系统,本文对可压缩流动的自适应求解进行了研究.研究内容包括基于误差矩阵的流场特征探测方法和网格自适应方法.网格自适应方法包含网格移动、网格加密、网格稀疏以及这3种方法的混合.基于误差矩阵的流场探测方法能找到流场特征的主方向, 通过网格自适应使得网格按主方向进行分布,可以减少该区域的误差而得到高精度的解.以边为单位能使得该流场的误差尽可能均匀化.本文以捕捉激波为例,对该方法的过程及应用进行了描述,结果验证了该方法的可靠性.     

空腔可压缩流致噪声问题研究进展

空腔是一种常见的飞行器结构,高速来流条件下,空腔流致噪声成为重要噪声源,制约飞行器性能的提高。针对空腔可压缩流致噪声问题的研究进展进行综述,分析空腔流激振荡及其诱导噪声的产生机制、参数影响规律与噪声控制技术发展趋势。基于空腔非定常流动及其诱导噪声特性,总结空腔可压缩流致噪声产生、传播规律以及相关研究方法,从来流参数、空腔形态、结构振动等关键参数出发,阐述空腔可压缩流致噪声的参数影响规律,总结空腔噪声控制技术的研究现状以及存在不足。最后,对空腔可压缩流致噪声问题研究面临的挑战和发展趋势进行展望。     

一种求解多介质流问题的正保护数值方法

在计算爆炸或高速射流等多介质流问题时,会遇到低密度低内能的情况.由于数值误差等原因,计算过程中会出现负的密度、压力或内能,从而导致计算无法继续推进.针对该问题,给出了一种基于rGFM界面处理方法和WENO高精度格式的正保护算法.算法首先使用具有正保护功能的求解器在界面处求解双激波近似黎曼问题,并用近似解更新界面附近真实...     

用热流探针测量激波速度

激波速度测量是激波管和激波风洞运行状态的一个重要参数,压电传感器或光学方法测速系统成本高,而传统电离探针在激波马赫数较低、波后温度达不到空气电离程度的情况下无法满足实验要求。提出了一种使用同轴热电偶作为测速探针来测激波速度的方法,弥补了电离探针在激波马赫数较低时的不足。介绍了同轴热电偶探针测速原理,并设计了测量激波速度的系统电路。通过信号放大电路锁定激波冲激信号,触发脉冲信号发生电路,实现了一种单通道、多测点的激波风洞测速系统。分别开展以温度与热流为触发信号的风洞实验,结果表明只有使用热流信号才能满足激波测速的时间要求。     

圆锥激波与平板边界层干扰的数值模拟研究

为进一步研究参数变化以及三维效应对圆锥激波/平板边界层之间相互干扰的影响,使用两方程Menter-SST模型,针对来流马赫数为2时的三维圆锥激波与平板边界层的相互干扰现象进行了数值模拟与定性、定量分析。分别研究了圆锥激波发生器半锥角和来流单位雷诺数变化对干扰区流动的影响,总结了参数变化引起的流动分离变化规律;此外,还计算了与三维计算的中心对称面上的入射激波等效的二维情形,并将三维结果与二维情形进行对比,对比结果显示中心对称面上的壁面压力系数、分离涡尺寸、涡量分布等与相应的二维情形存在明显差异。     

超声速喷流激波噪声基础问题数值模拟研究进展

超声速喷流问题是一个包含激波、旋涡、湍流和声波的多尺度复杂流动问题,其数值模拟方法及激波噪声产生机制是相关研究的长期热点和难点。本文简要回顾了超声速喷流激波噪声研究进展,重点介绍了针对激波噪声计算方法的非物理噪声消除技术和光滑因子设计准则,针对超声速喷流激波噪声研究设计的模型问题（包括旋涡–旋涡相互作用、激波–旋涡相互作用和激波–剪切层相互作用等）,以及轴对称和三维超声速喷流的研究进展。本文还介绍了作者最近针对超声速喷流开展的三维直接数值模拟、实验验证工作和初步分析结果（包括轴对称模态定位、束缚波演化和摆动模态发展等）。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

高精度混合型格式的构造及数值模拟

在守恒变量形式的无波动无自由参数耗散(NND)格式的基础上,引入样条插值的思想,并采用具有一致二阶精度格式的限制器函数,构造了一个在光滑区具有四阶精度,在激波附近降为二阶精度的混合型格式。在计算中采用具有总变差减少(TVD)性质的四阶R unge-K u tta方法进行时间离散,通过对多个典型算例的计算,并与NND格式作了比较,分析结果表明:混合四阶格式引起的耗散和波动较小,并且能够高精度地分辨间断。     

斜激波入射V形钝前缘溢流口激波干扰研究

针对内转式进气道溢流口这一关键部位所面临的三维复杂激波干扰问题,将溢流口提炼简化为V形钝前缘平板,采用激波风洞实验观测结合数值模拟的方法,研究了前体斜激波与V形钝前缘溢流口相对位置变化引起的激波干扰的演化规律。结果表明:由于V形钝前缘自身的激波干扰,其驻点前弓形激波的脱体距离较大,波后存在大范围的亚声速区。当斜激波入射在该弓形激波接近正激波的部分时,发生Edney第Ⅳa类激波干扰,该流动结构与V形钝前缘自身带来的三维激波干扰相互耦合,形成多处超声速射流区域;当斜激波入射在该弓形激波亚声速区的声速点附近时,呈现出不同于Edney第Ⅲ类激波干扰的波系结构;当斜激波入射在该弓形激波的超声速部分时,形成的波系结构与Edney第Ⅱ、Ⅵ类激波干扰类似。     

激波风洞气动力试验不确定度影响因素分析

参考北大西洋公约组织和AIAA推荐的风洞试验数据不确定度计算方法,结合激波风洞运行特点,确定激波风洞气动力试验的主要误差源,计算激波风洞B-2标模气动力测量结果的不确定度.采用改变单一变量的方法计算主要误差源对测量结果不确定度的影响程度,辨析对不确定度起主要作用的基本参数.计算结果表明:皮托压力和总压的测量结果对流场参数影响显著,皮托压力的测量结果比总压测量结果对流场参数与气动力测量结果影响更大;降低皮托压力和总压的偏离极限,有利于提高激波风洞气动力试验数据的质量.     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

高超声速进气道流场的数值模拟

运用自主开发的数值模拟软件对高超声速进气道复杂内流场计算进行了数值模拟,以验证其对高超声速进气道流动模拟的适用程度。通过某二元高超进气道模型内流场,德国RWTH Aachen所公布的高超进气道模型内流场计算,对该软件描述高超声速进气道内波系反射相交、激波/边界层干扰、高反压下隔离段内激波串等复杂流动现象的能力进行了考核与分析。计算结果表明该软件能够描述高超声速进气道内复杂的流动现象,即使在进气道承受高的燃烧室反压时,仍具有较高精度。但该软件中的紊流模型不能较好地预报边界层分离,特别是对分离区的大小及其诱导激波的强度等的预测存在一定的误差,因此需要进一步的改进和完善。     

液体压力激波加工技术研究

机械产品的零部件，通常需要通过成形与改性才具有所需形状及实用功能，本文提出了一种新型材料成形与改性加工的方法——液体压力激波加工。采用脉冲超声波聚焦技术，通过压电陶瓷激波发生器在液体介质中产生高能瞬时压力激波，对材料进行冲击，产生变形或表面改性。经过理论研究和试验验证，初步掌握了液体压力激波的产生、激波能量的控制和激波加工的动力学机制。这种加工方法具有能量可控性好、成本低廉和安全高效等优点，非常适合具有复杂结构和型面的中小零件加工，是一种颇具竞争力的先进制造技术手段。     

高速流动PIV示踪粒子跟随响应特性实验研究

示踪粒子的跟随响应能力是影响高速流动PIV测量精度的重要因素。针对法向马赫数大于1.4的高速流动所提出粒子松弛特性分析模型,结合理论分析与数值模拟方法,发展了高速流动下的示踪粒子布撒技术,提高了PIV技术定量化测量能力。基于上海交通大学多马赫数风洞,以不同粒径的氧化钛颗粒作为示踪粒子,利用PIV技术观测Ma4的高速流动诱导的一道22°激波,结果显示30nm粒径的示踪粒子有更优秀的跟随响应能力;并以该粒子进行了不同条件下（包括斜激波与脱体激波）的跟随性实验验证,为高速流动PIV示踪粒子选择提供了实验支撑。     

犖犘犔犛技术及其在高速飞行器气动研究中的应用

近年来,与高速飞行器相关的超声速/高超声速流动受到了极大关注。这类流动所具有的非定常性、强梯度和可压缩性对试验研究提出了挑战。纳米示踪的平面激光散射技术(NPLS)是2005年由作者所在的研究团队研发的非接触光学测试技术。它能够获得超声速三维流场的某个剖面的瞬态流动结构,并且具有较高的时空分辨率。目前,许多研究结果表明NPLS是研究超声速湍流的一项非常有效的技术。近年来,作者应用 NPLS 技术在超声速湍流研究中取得了较大的进展,并且基于NPLS开发了其它几种技术,比如基于 NPLS 的密度场测量技术(NPLS-DT),能够获得超声速流动的密度场信息并还能进一步得到雷诺应力分布。本文介绍了NPLS技术并回顾了其在超声速边界层、激波/边界层相互作用等流动中的应用。由于能够获得雷诺压力和湍动能等统计量, NPLS技术有望在发展可压缩湍流模型的研究中发挥作用。     

多孔介质传热传湿过程多层不连续问题的数值分析

提出了一种基于有限体积法预测非线性边界条件下多层多孔介质内的传热传湿过程的数值分析方法。求解过程中考虑了瞬态边界条件，从而避免了通常处理中边界条件设定为常数而给计算带来的误差，对于多层多孔介质每一层物性参数的非连续性，采用了有效的有限差分逼近处理。利用该处理对典型的三层墙体层与层界面处相对湿度的瞬态值进行了预测，计算结果与Liesen R J等的传递函数求解方法符合很好。