回流区中对涡与扩散火焰相互作用的数值模拟

对涡与扩散火焰的相互作用是涡轮喷气发动机燃烧室中十分重要的问题.本研究通过数值计算揭示了在对涡形成的流场中扩散火焰的结构以及扩散火焰对流场的影响.当不考虑释热引起密度变化时,数值解与相应的解析解符合较好.     

对流-扩散型方程的一种简单、有效的欧拉-拉格朗日分裂格式

本文提出了一种简单而有效的数值求解一维对流-扩散型方程的欧拉-拉格朗H(E.-L.)分裂格式。无论对于线性对流-扩散方程和非线性Burgers方程,或者是对流占优或非对流占优,它都能给出较好的结果。文中给出了四个典型而又较困难的计算例子,与精确解的比较表明,数值结果是相当令人满意的。     

逆风因子法在非正交曲线坐标系中解不可压N-S方程中的应用

所谓的逆风因子法就是分别用二阶精度公式计算出对流项和扩散项的贡献,然后按此两项贡献的比值来确定一个逆风因子,来保证对流扩散方程数值解的精度和稳定性。这种逆风因子法是作者在用有限元素法解不可压N-S方程中提出来的。本文将逆风因子法应用于在非正交曲线坐标系中用有限差分法解不可压N-S方程。      

关于病态线性方程组的迭代解法

本文基于简单迭代法的基础上，提供解病态线性方程组的一种迭代方法，并给出一些数值例子及应用于高雷诺数下计算对流扩散方程。     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

非定常动静叶干涉作用的二维数值模拟

采用涡方法研究叶轮机内动、静叶间的相互作用,它利用全三维的、通道时间平均的定常解为基础解,并由此给出非定常扰动场的初始解。为计算叶片对扰动场的响应过程,采用拉格朗日方法追踪扰动涡团的对流流动过程,用确定性涡方法来描述流体的粘性扩散过程。用该方法对某型压气机的一级（动、静叶栅距比为１．７４）的动、静叶相互作用,进行了二维有粘问题的数值模拟。数值结果表明,该方法具有较好的收敛性和收敛速度;转子出口的尾迹和各种旋涡等非定常因素在静叶通道中的切割、迁移过程,非定常性对静叶攻角的影响以及脉动压力的振幅衰减特性也在文中作了分析。     

数字全息法测量CO_2在RP-3燃油中的扩散系数

首先,介绍了数字全息激光干涉测量实验原理,设计并搭建扩散系数测量实验平台,运用MATLAB软件编写数字图像处理程序,利用傅里叶相位分析方法进行数字全息图像分析,得到扩散系数解。然后,分别测量温度为298.15K时0.33mol/L的KCl溶液和CO_2气体在水中扩散系数值并与文献研究数据进行对比,进行实验正确性验证,结果显示实验与文献标准值误差分别为1.63%和2.29%,满足实验要求。最后,运用此装置测试了278.15~323.15K时CO_2在RP-3航空燃油中的扩散系数随着温度升高扩散系数增大,并且不同温度下扩散系数值符合阿伦尼乌斯方程。     

一类新型自适应反扩散近似Riemann求解器及其应用

对于包含激波、剪切层等复杂结构的流动问题，为了精确模拟剪切层等精细结构，且保证激波计算的稳定性，必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大，Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解，容易导致不稳定。针对这一问题，本文在Riemann求解器中通过合理设计反扩散矩阵，发展了一类具有自适应反扩散的新型Riemann求解器，并将其应用到高阶加权紧致格式，实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性，结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷，既能准确识别剪切层等精细结构，又能保证激波解的稳定性。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

求解输流管道动态响应问题的一种有效方法

研究了两端铰支输流管道在脉动内流作用下的参数共振问题。根据数值分析中发现的参数共振响应特性,给出了第一和第二振型运动相互解耦的控制方程,并用谐波平衡法导出了参数共振响应曲线的近似解析解。算例表明,第一振型近似解与数值解吻合得很好,而第二振型解在共振点附近与数值解略有差距,但二者的整体走向完全相同。     

A low-diffusion robust flux splitting scheme towards wide-ranging Mach number flows

This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows. Based on Toro-Vázquez splitting approach, the new scheme splits inviscid flux into convective system and pressure system. This method introduces Mach number splitting function and numerical sound speed to evaluate advection system. Meanwhile, pressure diffusion term, pressure momentum flux, interface pressure and interface velocity are constructed to measure pressure system. Then, typical test problems are utilized to systematically assess the robustness and accuracy of the resulting scheme. Matrix stability analysis and a series of numerical cases, such as double shear-layer problem and hypersonic viscous flow over blunt cone, demonstrate that TVAP scheme achieves excellent low diffusion, shock stability, contact discontinuity and low-speed resolution, and is potentially a good candidate for wide-ranging Mach number flows.     

控制体有限元方法中的一种高阶插值格式

针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟采用各种迎风格式都容易产生假扩散的问题,本文给出了一种由插值单元及其迎风相邻单元所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟,表明该格式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩散。      

二维粘性跨音叶栅流场计算中逆风因子法的应用

本文是二维亚音叶栅粘性流场计算的推广。在非正交贴体坐标中用有限差分法解非定常运动方程;用逆风因子法确保其数值稳定性。采用总温处处相等这个能量方程。时间步长由运动方程的离散方程的系数非负性确定。采用k-ε紊流模型。对DFVLR跨音叶型叶栅分别在进口马赫数为1.05和0.82的条件下作了数值试验。计算和实验的对比是满意的。     

二维对流—扩散方程的控制体有限元法数值模拟

推导了二维对流—扩散方程的控制体有限元公式，其主要特点是所得离散化方程适合于逐点迭代法及共轭梯度法进行求解，使得该方法更适合于求解复杂的几何域及节省更多的计算机内存。通过几个算例表明，该方法对于特别大的Ｐｅｃｌｅｔ数其精度也是很高的     

平面叶栅不可压粘性流场数值解

本文旨在发展一种二维叶栅不可压粘性流场计算的工程数值法。采用了主变量控制方程。在贴体非正交曲线坐标中求解运动方程。应用了一种新的叠加网格。用逆风因子法来确保数值解的稳定。应用压力修正代数方程来修正压力场,以满足连续方程。应用等效粘性系数紊流模型。最后对四种不同的叶栅在攻角-13°～+11°的范围内进行了数值试验,并将落后角以及分离点的位置与平面叶栅吹风数据作了比较,相互吻合是满意的。      

NS方程激波计算的摄动有限差分方法

摄动有限差分(PFD)方法从一阶迎风差分格式出发,将差分系数展开为网格步长的幂级数,通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式,因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导,将摄动有限差分格式引入时间相关法的计算,并构造了守恒形式的摄动有限差分格式,然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明:本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率,而且保持了一阶迎风格式的无振荡性质。     

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式.该格式在激波过渡点降为一阶迎风格式.并从单个线形方程推广到非线性方程及方程组情形.通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的.