三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

计算带有层流分离的粘性/无粘流相互作用的一个方法

基于粘性/无粘流相互作用模型和边界层积分方程,推导出计算二维和三维不可压层流粘性/无粘流相互作用的方程组。该方程组可以用于带有中等程度的分离和再附情况下。文中给出了计算方法和算例。所得的二维计算结果与其他计算和实验结果相符。本法具有简便迅速的优点。     

边界层横流驻定扰动稳定性问题

采用直接对稳定性方程进行数值积分的方法,计算和研究三维可压（可带边界层控制）边界层横流驻定扰动（ｓｔａｔｉｏｎａｒｙｄｉｓｔｕｒｂａｎｃｅ）稳定性问题。导出积分初始边界处用于数值计算的矩阵表达式。通过改进的正则正交化并与误差判断相结合的方法,有效地控制正交化进程,克服了稳定性方程是刚性方程在积分求解中的困难。算例分析研究了跨音速后掠翼前部横流驻定扰动稳定性变化的关键因子和主要影响因素,结果是满意的     

基于非Favre过滤的动力学非线性亚格子模型

采用直接过滤的Navier-Stocks（NS）方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,并仅保留级数的一阶导数项。因为截断级数产生的误差相当于丢失了模型的部分耗散作用,采用一种动力学非线性亚格子模型来补偿丢失的耗散影响。用模型模拟了超音速湍流平板边界层及绕双椭球高超声速湍流,并将所得结果与经验公式和计算结果进行了对比。对比表明模型可以有效的模拟可压缩湍流流场。     

气冷叶片表面局部对流换热系数的数值计算

本文提出了联立求解二维可压缩边界层方程组和稳定导热方程的迭代法,从而修正了文献[1]假定物体表面温度等于来流总温的局限;同时提供了用这种方法计算气冷叶片表面局部对流换热系数的完整的FORTRAN程序。文中用近似积分法求解边界层方程组,用有限元素法计算叶片的稳定温度场。算例的结果与有关的理论及实验结果相接近。     

干扰剪切流动稳定性理论及其对高雷诺数流动数值模拟方法的改进

在干扰剪切流(Interacting Shear Flow,ISF)理论的基础上,提出ISF稳定性理论并把它用于改进高雷诺(Re)数流动计算方法。(1)高Re数内外绕流的RANS计算及工业标准PNS计算中,流动转捩的预测均基于经典边界层理论;然而转捩并非总是最早发生在边界层中,例如发生在壁面小突起、小凹坑、小窄缝等局部粘性/无粘强干扰区,这些强干扰区可能位于边界层内,但边界层理论并不适用于它们,又如转捩发生在分离点邻域强干扰区等。(2)ISF理论表明:高Re数内外绕流为一复杂ISF,转捩总是最早发生在该ISF的层流区中。(3)ISF稳定性理论表明:作者提出的干扰剪切扰动流(Interacting Shear Perturbed Flow,ISPF)方程组可以计算ISF层流中非湍流扰动运动演化并预测转捩;ISF方程组和ISPF方程组分别与PNS和抛物化稳定性方程(PSE)为同类方程组,PSE分析计算边界层稳定性的众多成功实践,说明用ISPF(即PSE)方程组计算ISF层流扰动流并预测转捩完全可行。(4)RANS和PNS方法经ISF稳定性理论改进后,在转捩前用ISF方程组(即PNS)计算ISF层流基本流,用ISPF方程组(即PSE)计算ISF层流扰动流并预测转捩位置;转捩后RANS方法计算RANS或RANS/LES,PNS方法计算干扰剪切湍流(ISTF)方程组即抛物化RANS(PRANS)方程组。改进后的两方法,理论合理正确,方程体系完备、自洽,ISF方程组只能用ISPF方程组相配对,因此是高Re数内外绕流计算的理想且可持续发展的两种方法。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

激波与壁面粉尘相互作用的数值研究

本文采用厚边界层理论建立了考虑两相耦合效应的气固两相湍流边界层方程用以研究激波与壁面粉尘的相互作用。方程中考虑了压力梯度，Ｓａｆｆｍａｎ升力和湍流扩散等因素对颗粒飞行轨迹的影响，并用Ｍｉｒｅｌｓ方程作为边界条件使方程组封闭可解。数值计算时，经适当变换的气相方程组用ＢＯＸ格式，拉格朗日坐标中的颗粒相方程用Ｒｕｎｇｅ－Ｋｕｔｔａ法，而源项则用颗粒源法和求质量加权的平均值获得。计算结果与某些实验结果较为     

高精度差分格式在机翼分离流计算中的应用

用一个模型方程分析了边界层方程的数值稳定性 ,指出稳定性问题会随计算雷诺数减小而变得严重 ,因而对较低的Re数 ,计算分离边界层流动采用高阶精度差分格式十分必要。本文同时给出用四阶精度差分格式求解在流动分离情况下边界层方程的技术方法 ;对一个大展弦比后掠翼在攻角等于 1 4°时计算了机翼上表面的分离线 ,并与文献上公布的用别的测算方法测出结果做了比较 ,两者总体上符合良好     

来流马赫数对附面层的影响

本文从附面层积分方程出发, 定性地分析了当表面压力系数 C_p分布相同时, 在保证Re数相等的条件下进口 Ma数对附面层发展的影响。并且用数值方法进行了模拟。与此同时用数值方法对由于表面曲率和自由流紊流度所引起的影响作出估计。结果表明; 在表面压力系数分布和 Re相同的情形下, 进口Ma数, 直到临界值附近时, 所引起的附面层的变化量是很小的, 一般在工程允许的误差范围之内。表面曲率所引起的变化同马赫数引起的变化相近, 而进口自由流紊流度引起的变化较之前两者较为明显。      

基于Newton-Raphson算法的径向轴承瞬态热流体动力润滑的研究

建立了径向轴承在载荷和速度突然变化时的三维数学模型，模型中考虑了轴瓦的热变形，在油膜和轴瓦交界面采用热流连续的理想边界条件，数值模拟轴承的瞬态温度场，并对轴承的瞬态性能进行分析。在每一瞬时，用Newton—Raphson算法同时求解Reynolds方程、膜厚方程和轴颈运动方程获得轴承油膜的压力分布和轴颈中心的运动速度，然后数值积分压力分布得到轴承的油膜力，差分运动速度得到轴颈中心位置和运动加速度。用一有效的有限差分法同时求解油膜和轴瓦的温度控制方程。最后将Reynolds方程和能量方程通过节点压力和温度相耦合获得轴承的瞬态三维温度场。结果表明本所介绍的方法收敛快，大大节约计算时间。     

Effects of adverse pressure gradient on Reynolds stresses inturbulent boundary layers

The effects of adverse pressure gradient (APG) on Reynolds stresses in turbulent boundary layers (TBLs) with APG were analyzed. The difficulty of this work was attributable to the Reynolds stresses in TBLs with APG under two combined effects, i.e.: effect of upstream flow and effect of APG. The effect of upstream flow is an inherent effect no matter pressure gradient exists or not. The individual effect was analyzed from absolute developments of Reynolds stresses in TBLs with zero pressure gradient (ZPG) firstly. Effect of APG was then analyzed from absolute developments of Reynolds stresses in TBLs with APG. Result showed that, for absolute development of mean streamwise Reynolds stresses, APG accelerated its development in TBL with ZPG; for absolute development of mean normal or shear Reynolds stresses, APG increased their magnitude in the outer part, and decreased their extent of large value region.      

非零压差边界下间隙流动雷诺方程近似解析解

为了获得间隙密封、滑动轴承等的液膜压力分布解析表达式,将间隙流动雷诺方程的解析解表示为特解和通解之和,其中特解表示为周向坐标函数和轴向坐标函数之和,通解表示为周向坐标函数和轴向坐标函数之积,结合分离变量法和非零压差边界条件,推导获得间隙流动雷诺方程的近似解析解。基于近似解析解计算间隙密封的液膜承载力和偏位角等参数,并与采用有限差分法的计算结果对比,其中液膜承载力偏差小于5%,偏位角偏差小于3°,而计算效率提高了数百倍。      

机械能方程的耗散积分解析表示式

发展了二维可压湍流边界层的耗散积分方法。采用两层涡粘性模型和速度剖面的指数律导出机械能方程的耗散积分解析表示式。它是局部摩阻系数、形状因子,Ma数和基于动量厚度的Re数的函数。该方法运算简捷,使用方便,大大简化了湍流边界层计算。和其他理论结果相比,本文结果更接近于实验数据。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

TRIP2．0软件湍流模块的开发与确认

基于TRIP2.0_SOLVER数值模拟软件,开发了湍流模型方法库,目前该CFD软件中集成的湍流模型主要包括B-L代数湍流模型、SA一方程模型和SST两方程模型.为了考核不同的湍流模型工程适用性,本文采用对接网格技术,通过有限体积法数值离散三维任意坐标系下的RANS方程组,应用LU分解、MUSCL差分格式和低雷诺数SA和SST两种湍流模型,数值模拟了二维NLR-7301两段翼型和三维DLR-F6翼身组合体的绕流流场,计算与试验比较的内容包括了表面的压力系数分布、典型站位的速度型和气动特性曲线等内容.通过计算与试验的比较,在本文的计算范围内,采用两种湍流模型均可以得到与试验结果相吻合的压力分布和升力曲线,但在边界层内的速度型、粘性阻力和力矩特性等方面,不同的湍流模型具有明显的差异.     

CFD/CSD方法分析动力效应对民机气动特性影响

基于Reynolds average Navier-Stockes(RANS)的三维Navier-Stokes流场控制方程耦合结构静力学方程时域分析方法，研究了带有发动机的民用飞机其动力效应对全机气动性能的影响。首先采用数值方法对发动机进排气边界条件进行了模拟，分析了带动力的涡扇发动机模型的流场，并将计算结果与实验进行比较，验证边界条件处理的准确性；以此为基础，考虑结构弹性变形，采用计算流体动力学/计算结构动力学(CFD/CSD)耦合的方法，分别对通气和带动力的翼吊发动机全机的气动性能进行了研究。结果表明:基于通气构型预测的升阻力系数，气动载荷和压心位置与考虑动力效应后的计算结果存在明显不同。弹性变形又会加剧这一差异，使得全机的升阻比下降约12.6%，升力系数下降约8.9%,压心位置后移。数值算例显示，在靠近发动机区域气动载荷受动力效应影响显著，远离该区域，弹性变形效应占主要影响因素，因此在进行带动力效应的民机气动性能分析时，考虑弹性变形的影响是十分必要的。      

边界层特性对雷诺数变化的敏感性分析

通过对有压力梯度下的边界层进行研究，模拟低压涡轮叶片表面边界层的发展过程，并分析了叶片表面不同负荷分布形式情况下吸力面边界层特性对雷诺数的敏感程度。计算中，分别采用了三种典型的涡轮气动负荷分布形式：后加载、前加载、均匀加载等，并分别模拟了三种不同的雷诺数下的边界层发展情况。结果表明：在同一雷诺数情况下，不同负荷分布形式吸力面边界层发展呈现出明显差异；而不同负荷分布形式下叶片表面边界层特性对雷诺数的敏感程度也不一样，相对而言，均匀加载形式对雷诺数的敏感程度最小，更适应低雷诺数条件下工作。     

再入飞行器和模型转换雷诺数计算

本文采用近似的激波关系、球锥体压力分布、平衡气体状态方程、边界层厚度计算公式及转捩准则等,按照流管法计算了4种再入飞行器的外流转捩雷诺数和转捩距离,并给出了高硅氧、碳石英、碳碳球锥模型的转捩雷诺数范围,对再入飞行器和模型的计算结果作了分析和比较。