非惯性非正交曲线坐标系中的边界层方程

应用张量分析方法将 N-S方程写于作非惯性运动的非正交曲线坐标系中 ,然后用匹配渐近展开法导出此坐标系中的边界层方程。消去压强项后 ,与惯性坐标系中的边界层方程相比 ,非惯性系中的边界层方程在 x1和 x2 方向的动量方程中分别只多了一哥氏力项 ,它们是由物体转动角速度的物面法向分量产生的 ,其它惯性力项均不进入边界层方程。方程中所含曲线坐标系的拉梅系数及其导数只需在物面上取值 ,物面曲率影响不进入边界层方程     

高性能歼击机进气道进口流场的数值模拟

应用数值模拟技术，采用“模拟前机身”的方法，把模拟对象按某种规律在长度方向上缩短1/2后，可将模拟流场与基准流场间的参数差控制在10%左右。吸除边界层后，可直接使用三维Euler法来判断机身缩短后的影响。而不必求解雷诺平均薄层N-S方程组，理论研究业已证明，用此方法，可在较小尺寸的自由射流式超音速风洞中开展带有飞机前部操纵面的进气道/前机身组合大模型的实验研究，从而解决了进/发匹配课题中的一大技术关键。     

绕扁长椭球体三维分离的积分法

本文利用三维积分边界层法并结合流线法决定涡层型流动分离。由三维不可压湍流边界层方程求得边界层参数后,按流线法计算粘性流线。作为算例,研究了一带迎角的扁长椭球体在不可压缩湍流里的绕流,经与水洞里所观察到的α=15°流型的对比,二者基本一致。     

高超声速边界层的转捩及预测

首先描述了边界层转捩的基本过程及研究内容。在此基础上,指出了高超声速边界层不同于不可压缩边界层的流动不稳定性特性,并介绍了边界层的转捩机理与感受性特征;给出了高超声速三维边界层中预测转捩的常用方法,并着重介绍了多用于工程实际的e N方法以及对e N方法的理性改进,同时列举了在高超声速三维边界层中应用e N方法实现转捩预测的多个实例。最后,分析并总结了高超声速边界层转捩预测所存在的困难及需要解决的问题。     

三维边界层的湍流特性及其二阶矩模拟

根据现有实验数据,本文进一步分析了三维边界层内纵向涡对雷诺应力分布的影响,从而建立了一个足以模拟三维边界层湍流特性的二阶矩封闭模式。在新模式中,着重考虑了三维平均涡场引起的雷诺应力输运的非线性作用。以文[3]给出的模拟无限后掠翼实验为例,本文分别应用基准的二阶矩模式和新建立的非线性二阶矩模式对三维边界层内的雷诺应力分布作了数值预计。结果表明:应用基准二阶矩模式所得雷诺应力诸分量的预计值均比相应的测量值大得多,且完全不能预计切应力方向相对于速度梯度方向的滞后;而非线性二阶矩模式对雷诺应力分布的预计则基本上与测量数据符合,且能定性地模拟切应力方向相对于速度梯度方向的滞后。     

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

再入飞行器迎风面边界层转捩的数值预估

基于求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,首先数值预估了高超声速三维流动的边界层特性;采用已有的高超声速三维边界层转捩准则,即在转捩位置Ｒｅθ／Ｍｅ＝常数,以及通过修正由ＯＮＥＲＡ提出的代数转捩模型,建立了预估三维边界层转捩的半经验数值技术并应用于三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的求解。利用ＨＡＬＩＳ模型及相关的实验数据,首次进行了在Ｍ∞＝６时不同雷诺数及迎角的验证试验。预估结果表明：在不考虑     

带有前缘气膜冷却燃气轮机叶片的气热耦合研究(英文)

使用气热耦合数值计算方法对带有"淋浴喷头"式前缘气膜冷却的C3X燃气轮机叶片进行数值模拟。通过对使用不同湍流模型的计算结果与实验测量值进行对比,使用转捩模型能够比较准确的模拟带有气膜冷却射流的边界层内的流动和传热过程。气膜冷却射流在叶片吸力面上游层流边界层区域内表现出比较明显的三维湍流特性,但在吸力面下游边界层转捩为湍流以后,湍流使得气膜冷却射流与主流燃气掺混过程加速而使其三维流动特性迅速消失;在吸力面由于边界层流动始终为层流状态,冷却射流与主流掺混较缓慢而使其仍能保持较强的三维湍流特性。边界层内由高温燃气与被冷却叶片之间温度差形成的温度梯度有益于增强边界层流动的稳定性。     

基于线性稳定性理论的转捩控制试验

 使用线化稳定性理论研究了某后掠机翼试验模型在风洞试验条件下可能的试验结果。在当前试验条件下,三维边界层中横流驻波的不同波长扰动发展趋势不同,根据理论分析和数值计算给出了最大放大率波长为4.0 mm。通过升华法试验验证了理论分析与数值计算结果。在试验中使用丝网印刷技术在后掠机翼前缘添加粗糙带,成功地引入了理论计算的最大放大率波长。根据预期影响了三维边界层的转捩。     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

后掠平板横流不稳定性

本文研究了无限后掠平板三维边界层的横流不稳定性，证明行进波和定常波几平同时失稳，且行进波比定常波更不稳定；     

方管内超声速湍流流动及旋涡的测定技术研究

本文在Ｍ∞＝１．５８及湍流边界层厚度δ／Ｄ＝０．１条件下，用实验方法研究了管长Ｌ／Ｄ＝６方管内的三维超声波湍流场及其沿程的发展和反压的影响，讨论了壁面边界层沿流向和展向的主要变化特征；指出管内横截面的三维静压场特点，文中最后探讨了超声速湍流流动的方管内旋涡的测量技术。     

高超声速三维边界层转捩数值研究进展及预测软件

高超声速边界层转捩问题已经成为飞行器设计需要考虑的一个关键因素，也是高超声速流动研究的热点和难点。三维边界层中存在横流、流向涡、接触线三种典型的流动结构和失稳特征，他们诱导转捩的机理不同。本文总结了低雷诺数湍流、γ-Re_θt、k-ω-γ、层流动能4种常用转捩模型的构造、适用性及应用情况，介绍了适用于三维边界层的BiGlobal、 PSE3D、 TriGlobal三种全局稳定性分析方法，分析了基于转捩预测e^N方法的e Malik、 LASTRAC、 LILO、HyTEN等转捩软件的优缺点，展望了高超声速三维边界层转捩研究及软件开发的未来。     

在三维湍流边界层中摩擦应力的测量研究

本文给出了V形壁面热膜测量物面摩擦应力的方法研究。实验表明该技术对三维湍流边界层中的摩阻测量是一类有效的方法。从本实验的数据可以看到,壁面热膜在二维湍流边界层中得到的校测函数能够直接用于三维湍流边界层的摩擦应力测量。     

超声速边界层三维扰动引起小激波的数值研究

通过直接数值模拟的方法,研究了M=4.5的超声速边界层中三维扰动的演化。以某一剖面为入口,加入一个及一对三维T-S波,发现随展向波数的增加,扰动幅值的增长率逐渐减小,证实了M=4.5的超声速边界层中,当三维扰动达到一定幅值时会有小激波出现,为建立可压缩流动稳定性理论提供了依据。     

高速三维边界层稳定性问题的数值方法

本文把四阶精度紧致差分格式用于求解声速、高超声速三维边界层稳定性问题。提出了适用地求解高速边界层特征值问题的Ｃｈｅｂｙｓｈｅｖ多项式配置点方法,并将各种数值方法进行了比较。数值结果表明上述两种方法是求解高速三维边界层稳定性问题的有效数值方法。     

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

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

利用κ—ε湍流能量闭合方法对城市街渠内气流结构的模拟

利用非静力、能量闭合的大气边界层模式分别绎二维情况下不同高宽比的街渠内部的气流结构和三维情况下城市街渠在来流方向不同的条件下的内部气流结构进行了数值模拟。并且对三维情况下街渠的上风向建筑物顶部的风向与街渠中心近地面的风向之间的关系和街渠内部上层的风向与街渠中心近地面的风向之间的关系进行了模拟研究，将这些结果与G T Johnson等人的现场观测和模式模拟结果进行了比较，表明两者有较好的一致性。     

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

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

用非线性涡粘性模式计算三维湍流边界层

针对非线性涡粘性模式在求解三维湍流边界层流动时的不足，本文从压力－变形率关联项中的快速项出发，考虑流动非均匀影响，在显式代数应力模式的雷诺应力表达式中引入了反映雷诺应力“松驰”效应的速度二阶导数项，构造了一个新的非线性涡粘性模式。通过对典型算例－翼体角隔流动的计算结果表明，新模式能较好地再现出三维湍流边界层内雷诺切应力方向的发展滞后于速度梯度方向发展这一流动特性。