振荡边界层大涡模拟.pdfVIP

振荡边界层的大涡模拟掌 张强 (北京理工大学宇航学院力学系， 100081) 摘要 振荡流研究在理论和实际应用上都有很强的重要性和研究价值．这类问题在自然界 中普遍存在，如：海洋中表面重力渡与海洋底床之间的相互作用，海底床面的形成及变化， 浅水中的泥沙输运问题等等(李家春，周济福等．1998；李家春和吴承康1998；李家春 和谢正桐1999；周济福和李家春2000；赵子丹，蒋昌波等200l；蒋昌渡，白玉川等 2003)．因为此类问题具有很强的非定常性，在简单边界务件下的情况对于研究湍演膜型能 否用于非定常流动很有帮助(SpalartBaldwin1989)． 海洋中的波向着岸边传播的过程中，由于海洋深度变化，导致表征边界层特性的参数，雷 诺款：忍=Un参，l，，发生从小到大的变化过程．其中：Ua为外部波动的速度振幅， 覆=√2y／国为Stokes厚度，缈为渡动的角频率，y为流体的运动粘性系数．振荡流的边界 层也会出现由层流向湍流发展的过程．在湍流的状态下，流动将表现与单向流动很大的差 异．因此，振荡流的研究，特别是层流向湍流的转捩，湍流状态下速度以及剪应力分布， 流场相干结构的研究引起了人们的很大兴趣． 关键词：湍流．振荡边界层．大涡模拟． 考虑简谐振荡的平板诱导的不可压流动，并引入以下无量纲变量： t=f‘缈，薯=i／4，％=强‘／砜，P=—等?万 (1) pUq 其中带‘为有量纲量，而最：,／T万g为Stokes厚度，uo为平板振荡的速度幅值。对 不可压缩流动的Navier-Stokes方程和连续性方程进行无量纲化后可以写成： 塑+堕v．(uu)：一生vp+!v2u、 ‘ (2) 研 2 2 2 一 V．u=0 (3) 计算区域如图l所示，底部为周期振荡的平板。垂寅方向(z)上尺度根据计算流动状 态不同选取为足够大的Stokes层厚度的倍数以满足顶部为自由流动的边界条件。流向(x) 和展向(y)的尺度选择满足在这两个方向上采用周期边界条件要求。 这样的物理模型经过简单的坐标变换，将坐标固结于振荡平板，模型就成为外部自 由流动为一U。COs(f)周期振荡诱导振荡流边界层。此时得到的整个流场和(Vittori Verzicco 1998)的直接数值模拟情形类似，因此非常便于与其结果进行比较。 ’国家自然科学基金资助项目 7 图1振荡流计算区域模型 在流向和展向采用均匀网格，而在垂向采用非均匀网格，壁面附近的网格越加密。 流向和展向速度以及压力的网格点则是间于垂向速度网格点之间，采用交错网格。 在大涡模拟中-，被分解成大尺度部分，和小尺度部分．，，即有：，2，+，。。大尺 度部分，通过计算直接得到，而小尺度部分的贡献采用Smagorinsky模型处理： 6 ．r—=。=一 zH一专f磕=-2(C,A／,)z心2s^S口 t4) ) 在振荡流中，对于计算中的一个物理量八x，J，，：，，J的统计平均用如下方法获得： 1—坚1墨丝 (5) f(x,y,z，t)=百j矗F∑∑∑厂(x，Y，z，f+2肼) f2l