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

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

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

再入飞行器非平衡气动加热工程计算方法研究

综合比较了现有的非平衡热环境工程计算方法,发现采用不同方法给出的计算结果相互之间差别较大。本文基于边界层中的原子复合反应主要发生在靠近壁面薄层内的特点,将边界层中的气相反应等价到表面上,建立了同时考虑边界层内非平衡反应和表面催化特性的非平衡边界层气动加热工程计算新方法。计算结果表明,非平衡效应主要集中在飞行器头部区域,下游边界层热流逐渐趋于平衡值。本文工程方法计算结果与有关飞行测量结果吻合较好。     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

气膜冷却涡轮叶栅气动损失数值计算

在二维气膜冷却边界层流动和传热计算方法 [1] 基础上 ,引入边界层积分参数法 [2 ] ,计算叶栅的气动损失。结果表明 ,这种方法较适用于低吹风比和边界层较厚处有气膜冷却的情况。作者将边界层积分参数法与叶栅损失计算的另一种方法进行数值计算的比较 ,表明两种方法各有长处与适用范围 ,实际应用中应考虑两者的结合使用。     

多级轴流压气机端壁边界层中叶片力亏损的研究

本文基于对有经向间隙和无任向间隙两种情况下叶栅二次流造成边界层内叶片力的偏转的分析建立了计算叶片力亏损的模型。文中给出了由瑞壁边界层叶片力亏损计算压气机效率的方法。用上述方法对两台轮流压气机和一个进口导向叶片排计算了端壁边界层厚度,得到了环面堵塞系数,计算了一台压气机的效率。对计算结果的分析表明,用这种方法可见得到较合理的计算结果。      

钝尾物体底部压力计算

 提出了一种简便可靠的计算钝尾物体表面压力分布及底部压力的方法。将钝尾物体绕流流场分为无粘区、边界层区和分离区。无粘区和边界层区分别用位流理论和边界层理论计算。分离区由根据尾流函数等建立的替代模型模拟。利用奇点法解出整个模型的压力分布。计算釆用迭代法进行。对Vanwagenen模型和鲁尔大学模型进行了计算,结果表明该方法是成功的。     

三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。     

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

自适应多重网格法在跨声速势流计算中的应用

用全近似格式求解了轴向大扰动，径向小扰动的跨声速势流方程。并采用自适应多重网格技术控制粗细网格层间的迭代转换，采用嵌套迭代技术较好地选取了初值。数值计算表明：自适应多重网格法可使计算效率提高８倍以上。     

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。     

利用雷诺方程数值模拟波形壁湍流扰动场

为了研究波形壁湍流及由波形壁引起的扰动场,本文将N-S方程的联立解法应用于湍流。扰动场是从时间平均的湍流场中直接分解出来的。由此计算的类Stokes层仅局限于壁面附近的薄层内,它的特征可以通过对扰动场的分析加以解释。     

可压缩非平行流边界层稳定性研究

 采用线性稳定性理论和多重尺度方法,研究三维可压缩的非平行流边界层稳定性问题。分析了可解条件的特征,导出精确计算所需的伴随问题方程渐近外边界条件的矩阵表达式,给出有控制的重正化方法,以有效地克服刚性方程在积分求解中的困难。探讨与非平行性作用相关的方程和影响因素,特别是新的特征函数畸变对扰动增加率的作用。通过算例,清楚地显示了流动的非平行性对边界层稳定性的影响。     

超声速尾迹-剪切流的混合增强

超声速条件下燃料和空气之间的高效混合是超然冲压发动机技术上的主要挑战。基于大涡模拟和流动稳定性分析,针对超声速尾迹-剪切流动开展了混合增强方法研究。尾迹的存在改变了混合层流动的速度剖面,对流动稳定性产生了重要影响,使混合层由三维最不稳定变为二维最不稳定,最不稳定扰动波频率和增长率增大。基于流动稳定性结果引入扰动的混合增强方式依然有效,根据稳定性结果设计了波纹隔板。数值结果表明：二维波纹壁引入的扰动未能增长,不具备混合强化效果,而三维波纹壁引入的扰动能够快速增长,具有混合强化效果,且波纹壁参数越接近最不稳定扰动波参数,混合强化效果越明显。     

超声速非平行边界层的非线性稳定性研究(英文)

用非线性抛物化稳定性方程(NPSE),研究超声速非平行边界层的扰动波的非线性演化过程。发展了高精度的数值方法,对NPSE进行有效求解和稳定推进。计算得到了各谐波的振幅增长曲线和不同的扰动剖面函数,特别是展示了谐波之间的非线性作用的物理特征,包括流向涡、展向涡和Λ结构的形成和演化等,以进行转捩过程的机理研究和分析。NPSE的计算结果显示了这个新方法是研究超声速边界层非线性稳定性的有效方法。     

超声速平面剪切层多模态混合扰动数值分析

对Ma_c=1.2二维空间发展超声速平面自由剪切层进行扰动模态及流动结构数值分析.使用线性稳定理论和扰动直接数值模拟方法对超声速剪切层中不同流向站位的流动剖面进行扰动参数分析和局部扰动直接数值模拟,获得了剪切层不同站位处的扰动特性.通过特定单频扰动下的直接数值模拟,在剪切层的不同区域内,无相差地促发了慢模态、混合模态和快模态等多种声涡模态,得到了特定单频扰动促发的多模态无相差不稳定波.      

合成射流微扰动对后台阶湍流分离流动控制的实验研究

后台阶流动是流体力学中一个经典的研究课题,代表着工程中一类横截面突扩的钝体绕流问题。后台阶流动分离会导致一些不利的影响,如高速旋涡的形成、流动损失、压力脉动以及气动噪声等。基于阵列式合成射流激励器对二维矩形后台阶湍流分离再附流动控制进行了研究,综合应用表面测压、七孔探针、粒子图像测速仪(PIV)和热线等多种实验手段,获取了后台阶的表面压力分布和非定常流场结构。结果表明:利用在台阶前缘形成的合成射流微扰动可使无量纲再附点长度降低25%,合成射流控制使得沿台阶下游的湍动能和雷诺应力增强,提高了台阶下游流场的混合效率。热线结果表明,频率是后台阶分离流动控制的重要参数,当频率为260 Hz,扰动频率与剪切层涡脱落频率之比为1.32时,合成射流控制可使位于1/2倍频的剪切层能量增强,仅需消耗较小的能量即可实现流动控制的目的。     

凹型粗糙元对边界层稳定性的影响

为了能够利用eN方法对带有凹型粗糙元的平板边界层进行转捩预测.通过数值模拟和流动稳定性的方法研究了马赫数为4.5的超声速边界层中凹型粗糙元对扰动演化的影响.对两种尺度的凹型粗糙元下扰动沿流向的幅值和增长率分布进行了计算.结果表明:凹型粗糙元对基本流的影响只是局限在凹型粗糙元及其附近很小的范围内.凹型粗糙元对扰动幅值的演化有抑制作用,尺度越大抑制作用越明显.在凹型粗糙元的后面,凹型粗糙元对扰动幅值的增长率的影响很小,相对变化量在2%以内.验证了对壁面有凹型粗糙元的情况,可以通过流动稳定性分析加N值修正的方法进行转捩预测.      

空中交通相依网络的脆弱性研究

机场、航路和管制扇区构成空中交通相依网络,节点在受到扰动时,网络运输性能下降,引发网络脆弱性。提出一种空中交通相依网络脆弱性的度量方法,首先构建机场网络、航路网络和管制扇区网络3个层网络,基于空中交通管理规则与层网络间的逻辑连接关系,建立空中交通相依网络模型。在随机扰动和蓄意扰动节点两种不同失效模式下,采用最大连通度和结构熵两个指标,并给出了相应算法,分析相依网络的结构脆弱性;创建流量熵和交通流损失比指标,设计了相应算法,研究相依网络的功能脆弱性。研究结果表明,随机扰动对空中交通相依网络影响更大;网络的结构脆弱源与功能脆弱源是机场网络;网络的脆弱性与层网络间的连接方式和层网络的交通量不匹配有关。