文档介绍：天津大学
博士学位论文
基于空间模式的平板边界层层流到湍流转捩的研究
姓名:李宁
申请学位级别:博士
专业:流体力学
指导教师:罗纪生
20070501
中文摘要本文基于空间模式的角度,采用直接数值模拟姆椒ǎτ肕并行技术,模拟了不可压缩平板边界层层流到湍流转捩的过程,继续研究了转捩过程中“惫痰哪谠诨恚氐憧疾炝似骄髌拭嫘拚诓懔髯^娴“惫讨兴鸬淖饔谩流动的雷诺数选定为玖魑2懔鰾相似性解,计算域入口扰动形式为一个二维ê鸵欢远猿频娜琓波。控制方程为直角坐标系下无量纲化的扰动形式匠蹋奔淅肷⒉捎萌拙ǘ然旌舷砸嘟岷的分裂格式,空间离散则采用高精度紧致有限差分逼近与逼追椒ㄏ嘟合的方法。法向采用非等距网格坐标变换,出口边界条件采用嵌边函数法。应用具体算例,验证了本文并行程序不仅符合线性理论,而且与串行计算结果一致。按照扰动幅值的大小,本文计算了亚谐转捩的三种情况,随后从壁面摩擦系数、平均脉动动能、谱空间扰动演化、矫嫠俣仁噶糠植肌⑵骄俣绕拭娴演化、流动稳定性分析、边界层的各种厚度及形状因子的演化等方面对计算结果进行了分析,得出以下结论:煌某跏既哦担^娣⑸谙掠尾煌奈恢谩H绻哦姆值较小,只有在“前的大约一、两个二维流向波长的范围内有明显三维扰动,且展向尺度基本上与入口处的三维波的展向尺度一致。转捩后更高的高次谐波被快速激发,扰动的展向尺度变得更小,同时平均流修正也增大的很快。凇癰”的前期,用当地层流边界层排移厚度无量纲化的壁面摩擦系数是相同的。在“暮笃冢诿婺Σ料凳旧显大到同样的大小。“钡墓涛个计算域入口处层流边界层排移厚度。骄龆芩婵瘴恃莼墓媛捎氡诿婺Σ料凳嗨疲扔诒诿婺Σ系数开始快速增长。总的脉动能量可以达到笥遥渲辛飨蚵龆量占五分之三,法向和展向脉动能量各占五分之一。另外,初始扰动越小,转捩越靠下游,各个方向的平均脉动能量都越大,但平均脉动能量快速增长区间的大小是基本一致的。煌跏既哦鹱^妫敝痢癰”前,二维波的幅值都是不一样的,但平均流修正开始快速增长的位置是相同的,平均流修正快速增长的过程也是大致相同的。平均流修正开始快速增长的位置与二维波
关键词:不可压缩平板边界层开始不规则的位置是一致的。在这个位置有规则的三维波达到最大,约为笥遥缓笕ǹJ急涞梦薰嬖颉U庖黄诩淙ㄔ龀ず芸欤也许是流动稳定性的二次稳定性在起作用。在转捩前,二维扰动幅值的演化与流动稳定性线性理论的预测基本符合,波长也与入口处二维波的波长近似相等,从某种意义上讲,这也许为转捩预测中的方法提供了理论依据。匝亓飨虿煌恢玫钠骄魉俣冉械奈榷ㄐ苑治霰砻鳎刈帕鞫向,不稳定区域先逐渐减小,缩小到很小以后,再逐渐的增大到最大范围,然后再一次逐渐减小,趋于很小。不稳定区域开始逐渐增大的位置对应于壁面摩擦系数开始快速增长的位置,不稳定区域达到最大的位置大致位于壁面摩擦系数曲线快速增长段的处,不稳定区域再一次逐渐减小,趋于零的位置对应于壁面摩擦系数达到最大的位置。当层流中的扰动幅值增大后,由于非线性作用将修正平均流剖面。当平均流剖面被修正到一定程度,不稳定性区域会扩得很大,这将导致更多的谐波被快速激发并增长,开始了“惫獭T凇癰”过程中,由于大量谐波快速增长,导致扰动能量快速增长。扰动能量的快速增长又进一步加速了平均流剖面的快速修正。这样相互的快速增长使得流动快速地转变为湍流。这里的结论与其他流动中的研究结果是一致的。谧^婀讨校乓坪穸鹊脑龀け认嘤Φ牟懔髋乓坪穸鹊脑龀ひB转捩直接数值模拟空间模式稳定性分析发卡涡
.、....甌..,,.、,·
,甒瑃甀甀,甌琤琽玛甀,瓺,、,.,畉猟..琲瑆瑃甌.
图目录情勋中个整周期时扰动速度婵占涞难莼肌情础中稣芷谑比哦俣萓的模随空间的演化图⋯⋯⋯情黝稳定性分析矫嫔系闹行郧摺情黝中稣芷谑比哦俣萓随空间的演化图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.·烟线演化示意图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯计算域示意图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯..狻:;嵌边函数的形状示意图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.局部超松弛迭代示意图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.压力迭代并行示意图⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.三维ㄈ肟谌哦俣确植肌计算结果与理论值对比⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯..串行、并行程序计算结果的比较⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.计算域入口处的特征函数⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯..⋯⋯⋯⋯.情况一个整周期时扰动速度材在驴矫娴牡戎迪咄肌情况稣芷谑比哦俣取暗哪K婵占涞难莼肌三种情