应用维数分裂方法推广MUSCL和WENO格式的若干问题

首先对有限差分法和有限体积法的差异进行了讨论,在已有文献观点的基础上补充了二者在边界条件处理及网格需求等方面存在差异的新论据。介绍了逐维推导的MUSCL和WENO格式计算控制体界面通量的过程,认为此类格式计算界面通量的方法直接应用于高斯积分型有限体积法不够严谨,从而得到了应用维数分裂方法构造的MUSCL格式和WENO格式不属于高斯积分型有限体积法的观点,“积分格式”这一定义更能准确反映这类格式的特点。此外,讨论了MUSCL格式和WENO格式在曲线坐标系下不能保证守恒的原因,并简单介绍了消除方法。     

粘性超声速绕体流动的有限谱法计算

首先推导出粘性可压缩流体动力学基本方程组在广义坐标下矢通量分裂的一般形式，并将有限谱ENO格式用于求解超声速流动问题。通过典型算例，验证了这一格式的精度和稳定性。然后进行了圆柱和机翼超声速绕流的数值模拟，得到了与实验和其他计算相符的结果。这表明有限谱法不单具有精度高，应用灵活等多方面的优点，而且突破了传统谱方法全域性的局限，可以成功地应用到一般形状流场可压缩流动的数值计算中。     

可压缩后台阶流的数值模拟

基于可压缩Navier-Stokes方程,采用Beam-Warming矢通量分裂格式以及vanLeer的矢通量分裂方法,分别通过Johnson-King及Baldwin-Lomax湍流模型对亚音速和超音速的后台阶流进行了计算,计算结果与实验符合较好。     

一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

用于二维无粘跨音速流计算的隐式TVD有限体积法(英文)

本文结合Harten的TVD(Total variation diminishing)概念和有限体公式建立一种求解守恒型Euler方程组的有效方法,采用通量分裂的迎风格式对方程组的隐式算子部分离散化,得到系数矩阵是主对角线占优的线性方程组,采用线Gauss-Seidel迭代推进法求解,提高了数值收敛速度,给出了二维跨音速内流的计算结果。     

用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

一种混合格式的构建及其在多尺度流动中的应用

面向非定常湍流数值模拟方法精度提升问题,构建了实用的矢通量分裂(FVS)形式高阶精度混合格式。该格式平滑区具有五阶精度线性紧致格式的高分辨率特证,而强间断和强梯度区呈二阶耗散型格式特征,进而保证全流场计算的稳定性。基于所建立的混合格式开展有限体积大涡模拟(LES)计算,将其应用于三维高Reynolds数非定常湍流预测中。结果显示,所建立的数值格式能够可靠预测湍流多尺度涡结构特征,且计算量小,具有较高分辨率和鲁棒性。     

跨、超音速流动的区域分解方法与并行算法

研究二维跨、超音速无粘性流动的 Euler方程区域分解方法、并行算法及其应用。通过内边界耦合条件实现相邻子区域解的光滑过渡,以得到总体流场的数值解。发展了一种多块区域之间守恒型的有效内边界耦合方法,对二维翼型跨音速流动和钝头体超音速流动等进行了分区数值求解,分区计算结果与其他单区计算结果作了比较,并讨论了多种区域分解数目的分区计算效率。并行计算采用纯结点并行编程方式和“先进先出”的同步控制等待机制,利用 PVM并行环境对二维绕翼型跨音速流动做了二区和四区分区并行计算。     

流体注入的轴对称矢量喷管三维流场计算

采用Roe通量差分分裂格式对基于流体注入控制的轴对称矢量喷管内流场进行了数值模拟。流体注入的位置分别为前孔和后孔,注气压强比为0.75～2.0,注气流量比为2.5％～10.0％,矢量角变化范围为2.8°～7.8°。计算结果表明:随着注气流量和注气压强增加,流体注入所产生的喷管矢量角相应增加;注气位置对喷管矢量角影响较大,注气位置靠近喷管尾沿(后孔注气)比注气位置靠近喷管喉部(前孔注气)所产生的矢量角明显增大。     

直射式喷咀垂直跨流喷油与沙丘驻涡稳定器匹配的实验研究

试验研究了直射式喷咀垂直跨流喷射时, 进口气流速度、喷咀压降与稳定器径向位置等参数对沙丘驻涡(BD)稳定器截面燃油浓度分布的影响, 并与沿z向喷射情况作了比较。结果表明, 垂直跨流喷油与BD稳定器的匹配比沿z向(即水平方向)喷油好。为使良好气动特性与燃烧稳定性的BD稳定器获得推广与应用, 必须考虑喷油方式、气流速度、喷咀压降与稳定器径向位置等参数对燃油浓度分布的影响。      

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

基于CATIA V5系统的非金属曲面零件展开方法研究

提出了一种应用于飞机非金属曲面零件简便实用的展开方法。基于CATIA V5系统环境下,应用可展直纹面分块拟合零件定义复杂曲面,将零件三维轮廓按拟合直纹面进行分块展开,所得的分块展开轮廓在二维工程图纸环境中依据给定的规则拼接成零件完整展开图形,实现此类非金属曲面零件的展开。分析了展开方法产生的误差,给出了实用的展开精度控制方法,实现了在未进行二次开发情况下,获得满足此类零件工程制造精度要求的二维展开图形数据。     

Box法用于三维边界层计算的数值稳定性分析

本文对 Box法用于三维边界层计算的数值稳定性做了 Neumman分析,并借助于数值计算,证实了目前公认的一些稳定性问题和三维边界层计算中所发现的现象,得到了一些有参考价值的结论。      

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

用热流探针测量激波速度

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

二元喷管流体矢量控制方案数值研究

采用基于雷诺平均的三维N-S方程和RNGk-ε湍流模型对两种流体注入方式下的某型二元收敛-扩散喷管全流场进行了数值模拟.计算结果表明:在激波诱导和喉道倾斜矢量控制方案下, 扩散段射流位置对喷管矢量效率的影响趋势是一致的, 射流位置越靠近喷口, 喷管获得的矢量角越大.相同射流流量在同一位置注入时, 由于喉道倾斜方案下的喷管主流可以实现亚声速偏转, 所以其总压损失较激波诱导方案要小.