振动斜板上下落液体薄膜的数值模拟.pdfVIP

维普资讯 计 算 物 理 第23卷 第 1期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS V0I．23．No．1 2006年 1月 Jan．．2006 [文章编号] 1001．246X(2006)O1．0057．04 振动斜板上下落液体薄膜的数值模拟 胡国辉 ， 胡 军 ， 尹协远 (1．上海大学应用数学和力学研究所，上海 200072；2．中国科技大学力学和机械工程系，安徽 合肥 230027) [摘 要] 用 Galerkin有限元法数值模拟了振动斜板上的下落液体薄膜流动 ，研究了流场中不同强迫频率和初 始扰动频率的影响．结果表明由于它们之间的相互作用 ，流场中出现低频信号，其在流场 中的非线性效应使得 表面波发生合并，加速了表面波的发展．对较高的强迫频率 ，这种效应更明显，且可观察到孤立包的产生 ． [关键词] 薄膜流动；Galerkin有限元；自由表面 [中图分类号] 0359 [文献标识码】 A 0 引言 斜板上下落薄膜的流动是流体力学中一个经典的理论问题 ，在工程实际中也有广泛的应用背景 ．该领域 的理论研究始于 Benjamin 和 Yih 的线性稳定性研究．Chang对流动中的线性和非线性稳定性问题和表面 波的稳定性演化研究做 了详细的评述 ．目前大量 的工作是基于长波近似 ，将 Navier-Stokes方程简化为关于 薄膜厚度的Benny方程 ，并进行求解和分析 ．尽管这样 的处理对许多物理现象给 出了合理 的解释，最近 的研 究表明即使在流动的线性阶段，长波理论的描述也不是足够准确的H．因此，近年来陆续有一些基于 Navier- Stokes方程的稳定性和数值模拟工作发表 ． 在许多工程问题中，出于流动控制的考虑或其他因素，流动 中可能牵涉到斜板振动的问题 ．美国Clarkson 大学un的课题组用 Chebyshev谱方法结合 Floquet定理对这种流动的稳定性进行了研究 ，分析了重力波、 剪切波和参数共振 Faraday波的相互作用情况，并认为在某些参数下，大振幅的振动可能抑制流动不稳定的 产生 ]．本文旨在从Navier-Stokes方程出发，用数值计算研究流场波动的非线性演化过程． 1 物理模型和数值方法 考虑粘性不可压二维液体薄膜在重力作用下沿水平倾角为 0的斜板流动(图1)．引入直角坐标系，设沿 斜板方向为 轴 ，Y轴为垂直斜板方向．无量纲后雷诺数定义为 Re=gH3sinO／2v，佛罗德数 Fr=u／,／gH，其 中g是重力加速度 g的大小，u是运动粘性系数，日是一维均匀 Nusseh流动的薄膜厚度 H=(3vRe／ sin)”，U是平均Nusselt速度 ，U=gsinOH213v．自由表面高度随 和时间t的变化是h(，t)．流动可以用 有周期外力作用的Navier．Stokes方程和连续性方程描述为 刮 + = 吉p+1 +专g一 )， (1) ． ： 0． (2) 式中 =(u，)是速度矢量，P是压力，』D是密度，A(t)=(dsin(2rcfet)，0)是单位质量流体受到的外部激励． 在平板上使用无滑移和无穿透条件 ， u： ：0． (3) 在自由面上使用运动学和动力学边界条件 ， 【收稿 日期】2004—10—21；修【 回日期 】2005—02—06 [基金项目]国家 自然科学基金 10472062)和上海市高等学校青年基金资助项 目 [作者简介]胡181~(1969一)，男，江西南 昌，副教授 ，博