槽道剪切流动中湍流斑问题研究.pdfVIP

槽道剪切流动中湍流斑问题研究 杨颖朝 (南京航空航天大学空气动力学系， 南京 210016) 摘 要 通过对 Navier－Stokes 方程的直接数值模拟，研究槽道剪切流动中湍流斑的形 成和发展。采用高效的高阶紧致差分方法，模拟从壁面脉冲到湍流斑的演化过程,展示 了湍流斑内物理量脉动和涡系的形成等湍流的基本特征及其湍流斑之间的相互作用。 关键词 湍流斑，Navier－Stokes 方程，直接数值模拟 在近壁剪切流动的转捩过程的后期，会出现具有一定特征的局部湍流区域，尽管它 的周围流动依然为层流，但在这个区域内部却具有随机脉动等湍流特征，该区域被称为 湍斑。这些湍斑随机的出现在近壁区域的不同位置，并随着它向下游的传播而逐渐增长 和抬升，然后相互混合形成完全湍流。Emmons 在水流实验中首先发现了这种湍斑，后 来发现这是转捩过程中的重要特征[1] ，从而引起了人们的广泛关注，采用了实验[ 2, 3]和数 值[4]等不同方法对其进行研究，并取得了一定的成果。但是，由于问题的高度复杂性， 关于湍斑的形成和发展机理等许多问题还不完全清楚，有待进一步深入研究。 壁面粗糙（凹凸）形成的脉冲是一种常见的扰动形式，是诱发湍斑的重要成因之一。本 文以壁面脉冲作为初始扰动场，采用直接数值模拟的方法研究槽道流动中湍流斑的发展和演 化过程，展示湍流斑演化及其湍流斑间的相互作用。 1．控制方程 这里的控制方程采用无量纲形式的非定常不可压 Navier－Stokes 方程和连续方程: ∂U 1 2 +( •∇) −∇ + ∇ （1） U U p U ∂t Re ∇•U 0 （2） T 式中：U (u, v, w) ，u u / U ，v v / U ，w w / U ，为无量纲速度矢量，这里的是流 2 U 向最大速度;P p / ρU ,为无量纲压力；雷诺数 Re Uh /ν ； t tU / h ，为无量纲 2 时间；▽是梯度算子；▽ 是 Laplace 算子。其中 h 为槽道半高度；u ，v ，w ，分别是流 向、法向和展向有量纲速度；ρ 为流体密度，p 为静压，t 为实际时间，ν 是运动粘性系数。 2．数值方法 发展四阶时间分裂法，离散方程 （1）和（2 ），得半离散方程组。为了保持精度的一 致性，空间离散均采用四阶精度：对于化得的 Poisson 方程形式，用中心差分方法离散并 保留四阶精度；化得的 Hemlholtz 方程形式，用四阶精度三维耦合紧致有限差分格式。 最后构成了不可压 Navier-Stokes 方程的三维耦合高阶紧致差分解法[5] ，其初值和边界条 件根据具体情况给出 3． 初始模型和边界条件 1 3．1 初始模型 湍流斑的初始模型是决定湍流斑形成和发展的重要因素之一. 根据实验观察, 湍流斑的 出现无论从时间上还是空间上都具有随机性, 但是在满足一定的条件后, 湍流斑就可能出现. 湍流斑的形成原因复杂多样, 而壁面扰动是一重要的触发因素. 这里合理地选取因壁面粗糙 而诱发的壁面速度脉冲作为湍流斑的初始模型, 它可表示为 u′ δu δΩ, t 0 其中δu 为扰动强度。 3．2 边界条件 入口边界：采用无扰动边界条件 ∂U   速度： 