抛物化稳定性方程在曲面边界层中的应用

为了研究抛物化稳定性方程在曲面边界层中的应用,选取了等曲率圆环管道、等曲率板和NACA0012翼型3种典型的曲面边界层.通过计算小幅值扰动波的演化,抛物化稳定性方程的计算结果和线性稳定性理论、数值模拟扰动方程结果相符,说明可以使用抛物化稳定性方程研究曲面边界层的小幅值扰动波的演化和稳定性分析;通过计算有限幅值扰动波的演化,线性阶段抛物化稳定性方程计算结果和扰动方程相符,在非线性很强的阶段,抛物化稳定性方程计算发散,发散的位置可作为转捩位置.      

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

用椭圆型、抛物型两种方法计算平面叶栅湍流尾迹

本文用连续方程、稳定不可压缩椭圆型和抛物型时均N-S方程及K-ε双方程湍流模型计算了平面叶栅湍流尾迹。计算表明,上述两种方法都能得到较好的数值结果,但抛物方法可节省大量计算时间和内存空间。     

三维不可压边界层抛物化稳定性方程的椭圆特性研究

针对三维平行不可压边界层流动的抛物化稳定性方程,研究其剩余的椭圆特性,提出了消除该特性的具体方法,并对抛物化稳定性方程解与全椭圆问题解的一致性进行了证明。     

抛物化Navier-Stokes方程(PNS)和高氏PNS理论及其应用的评述（英文）

抛物化NS方程得到广泛应用,已经成为工业标准气动计算的基础。现有的八种抛物化NS方程有不同的名称,方程中粘性项的形式略有不同,其中的PNS和薄层(TL)NS方程应用最多。但是这些方程都具有类似的数学性质,例如,当流向方向上马赫数大于1时,他们都是抛物型方程,可以采用空间推进算法(SMA)进行求解。与采用时间推进算法求解的NS方程或雷诺平均(RA)NS方程相比,PNS-SMA计算降低了空间的维数,节省了大量的存储空间和CPU计算时间。PNS-SMA算法也获得了巨大的进展。但是,早期PNS研究在理论上是相当模糊的,高智在1990年提出的粘性/无粘干扰剪切流理论(ISF)弥补了这一不足。ISF理论概括了PNS方程所能描述的基本流动,提出了其流动的运动规律及数学定义式,所导出的ISF方程组也属于PNS方程的一种。为了不增加新的名称,我们将ISF方程组也称为高氏PNS理论和方程组。这一理论在NS方程和RANS方程的计算中均有重要的应用。例如,计算最优坐标系的选择以减少伪扩散,网格尺度选择及局部网格加密设计以捕捉高超声速流动中物体表面热流等的急剧变化,壁面压力边界条件的选择以及由高PNS导出的壁面判据来进行NS和RANS近壁数值解可信度评估。本文评述了一些初步的应用,进一步的应用和综合PNS-SMA,RANS-SMA以及PSE-SMA计算值得深入研究,这里PSE指抛物化稳定性方程。     

超音进气道的亚音扩压器流场的数值计算和实验研究

本文对二维、亚音扩压管道中的粘性流动进行了数值计算,并对同样的扩压器模型进行了风洞吹风实验。 根据要计算的二维、亚音、壁面扩张角不太大的管道中的粘流,我们对N-S方程进行了抛物化假设,化成可沿流程推进求解的抛物型方程,再用有限差分进行计算。 经与实验结果对比,证明本文采用的结合零方程紊流模型,在几何流线坐标系下求解抛物型方程的方法对解决二维、亚音扩压管中分离区以前的流动问题,可得到满意的结果。     

扩压器内二维湍流流动的时间推进数值解

ＡＦ有限差分理论。指出薄层Ｎ－Ｓ方程抛物化后右端粘性项是主元素变量矢的一个简单线性函数,这对实际计算很有用。应用隐式ＡＦ方法计算了为发动机进气道设计的典型扩压器内的复杂湍流流场。研究表明,计算结果与实验结果基本一致,这为三维非定常湍流流动进一步研究提供了基础。     

缓变过渡段几何外形对实验段内气动力特性的影响

本文提出了一种计算高速风洞三维过渡段的气动特性的方法。本方法是以一维管流解为基础,并通过坐标变换,将变系数的线化速度势方程化为波动方程(对超音速流)或拉普拉斯方程(对亚音速流),然后用有限基本解法求解。本文着重讨论了超音速流求解的过程。在超音速流时,对矩形实验段区域,则采用分离变数法求解。作为例子,计算了由圆变方的超音速流过渡段,得出了三维变化对实验段气流参数的影响。     

改进的小波均匀化方法求解小周期系数抛物型方程

对传统的小波数值均匀化方法通过加入逐层修正的步骤进行了改进,并将改进后的方法用于求解具有小周期系数的抛物型方程。数值算例结果表明：在保持较高计算效率的同时,改进后的方法比改进前方法在精度方面有明显的提高。因此,改进后方法是小周期系数抛物型方程的一种高效率、高精度的数值求解方法。     

小波变换数值均匀化方法求解振荡系数的方程

将小波数值均匀化方法用于求解二维快速振荡系数的抛物型方程.小波数值均匀化方法基于小波的多分辨分析,在细尺度空间上得到原方程的半离散方程,然后利用小波映射得到粗尺度空间的数值均匀化方程,即可在粗尺度上求解原方程,有效地减小了计算量.计算结果表明:小波数值均匀化方法与精细剖分的有限差分法相比较,既大大地节省了计算时间又获得了较好的精度.     

90°圆截面弯管内三维紊流场实验研究

用CTA热线风速仪和五孔探针, 对一曲率直径比R_c/D=0.87, 直径D=280mm的90°圆截面弯管内部三维紊流流场进行了详细测量和分析。给出了不同雷诺数下弯管内不同截面上时均速度场、紊流动能、总压及静压沿曲率半径方向的分布。      

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

圆截面U形管流有限元模拟

弯曲圆管内流动具有很强的工业应用背景,管内二次涡的大小、位置及发生的频率会对流动系统的阻力及热质交换产生重大作用,研究弯曲圆管内的流动具有重要的理论意义.本文应用湍流大涡模拟方法对弯道流动进行了数值模拟.用时间推进法,给出了隐含流线迎风耗散作用的二阶完全展开Taylor-Galerkin有限元离散格式,数值计算的结果与实验结果符合较好.表明大涡模拟方法适用于边界形状复杂,存在各向异性的大尺度涡的流动仿真.     

埋入式进气道流场数值计算与实验研究

应用有限体积法对埋入式弯曲进气道内外流场进行了数值计算。用ＴＴＭ法生成贴体网格;用ｋ－ε模型模拟复杂湍流流动。计算结果表明,前方气流一部分被吸入管道,另一部分沿机体表面流向下游。管道入口背风侧存在一个分离区,管道内部气流畸变严重。计算结果和实验数据吻合较好。     

瞬时起动圆柱的初期流动

 本文将ADI方法和Poisson方程直接法结合起来求解涡量-流函数形式的N-S方程,研究瞬时起动圆柱的初期流动(Re=3000,5000,9500)。成功地算出实验观察到的α流动结构、β流动结构,以及这些结构在流动发展过程中交替出现的情况。本文方法效率高,计算结果与实验结果及他人的数值计算结果相符较好。可期望用该方法计算其它复杂流动。     

绕任意轴旋转的矩形截面弯管流动特性研究(续)

利用数值计算的方法研究了绕任意轴旋转的矩形截面弯管内的粘性流动情况，分析了在不同入口距离处离心力和科氏力共同作用下的二次流动、轴向主流的发展变化情况。计算结果表明：F=-4．0时，流动呈现出4一涡结构；F=-4．0的轴向速度矢量图在垂直面内与F=-4．0的轴向速度矢量图正好相反；在大曲率情况下管道内会出现三个轴向速度极值点。     

包含平均旋流的圆环管道声传播的模态分析

通过对三维线化欧拉方程的特征模态分析,研究包含平均旋流情况下声波在敷设声衬的圆环管道中的传播.特征模态分析采用谱配置方法计算轴向波数谱及扰动量的分布.当有平均旋流存在时,特征模态谱主要包括近/纯声模态和近/纯对流模态,旋流的存在甚至会改变传播模态个数,从而改变声波传播截止特性.在验证程序正确性的基础上,研究流动参数及声衬对声模态及对流模态的影响.      

不等熵一维定常管流数值解法

不等熵一维定常管流控制方程的数值解法,已有许多讨论,但大多局限于某些特定情况,特别是控制方程在奇点(M=1)附近的计算,还没有一个简便的方法.本文引进新变量,把方程变为非奇异的,提出了对各种不等熵流都适用的数值解法.由于引进了喉部条件,保证了数值解的唯一性.经过对模型方程的计算精度分析,证明该数值解除了在奇点附近稍有误差外,其它地方与精确解一致.本文最后给出了该方法在无喷管发动机内弹道计算中的应用.     

某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

非均匀来流对新型高超弯曲激波二维进气道的影响

针对一种采用新型压升规律的曲面压缩面所设计的高超弯曲激波二维进气道,利用数值模拟方法研究了6种不同长度平板发展而来的非均匀流对其性能参数的影响,并和两种常规二维进气道吞非均匀来流的能力进行了比较,着重对比了设计Ma数下三种不同压缩形式进气道喉道截面的流场畸变程度。数值模拟结果表明,采用新型压升规律的高超弯曲激波二维进气道性能受来流气流不均匀性的影响较小,且喉道截面的流场畸变指数小,对非均匀来流具有一定的校正作用。