涡旋演化的小波自适应模拟-应用数学和力学.pdf

摇 , 32 1 Applied Mathematics and Mechanics摇 摇 2011年1月15 日出版 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 Vol.32,No.1,Jan.15,2011 文章编号:1000鄄0887(2011)01鄄0033鄄11 訫 应用数学和力学编委会,ISSN 1000鄄0887 涡旋演化的小波自适应模拟* 1 1 2 赵摇 勇 宗摇 智 邹文楠 ,摇 ,摇 (1.大连理工大学 运载工程与力学学部 船舶工程学院 工业装备结构分析国家重点实验室,辽宁 大连 116024; 2.南昌大学 工程力学研究所,南昌330029) (林建忠推荐) 摘要:摇 该文考察了小波自适应方法用于涡旋运动的演化过程.首先,通过两个初边值问题,说明 小波方法具有可精度可控和局部结构自动捕捉的能力.然后,计算了涡旋的合并过程,结果表明, 小波方法可以准确高效的应用于流动涡旋的演化预测,进而,讨论了小波方法在湍流数值模拟中 的应用. 关摇 键摇 词:摇 小波自适应;摇 涡旋演化;摇 相干结构;摇 涡量方程;摇 湍流 中图分类号:摇 O357摇 摇 摇 文献标志码:摇 A DOI:10.3879/ j.issn.1000鄄0887.2011.01.004 引摇 摇 言 小波基函数具有许多优秀性质,比如有限支撑基函数、对称性、正交或双正交性质等,使得 [1鄄4] 小波方法可以用于微分方程数值模拟 .相比于Fourier方法,小波方法的最大的优点是基函 数有限支撑性质.在Fourier方法中,基函数是三角函数,具有无穷域支撑,它的展开系数是整 个域分布的平均效果,局部结构信息很难集中展现.而在小波方法中,小波基函数是有限域支 撑的,通过伸缩和平移变换,小波分析可以细致反映局部结构信息,据此小波方法是适合求解 [5鄄7] 局部结构问题的 . 能准确地模拟涡旋演化在流动数值模拟特别是湍流数值模拟中非常紧要,因为涡旋普遍 [8] 分布在流动中,随流动迁移并发生伸缩扭转变形 .正是由于涡旋形态的多样性和涡旋的多 尺度分布特征,极大地增加了湍流数值模拟的计算量和难度.湍流本质上是三维的,其数值模 拟应当通过三维流动控制方程来计算.本文旨在探讨小波用于湍流数值模拟的可行性,将通过 一个二维涡旋演化的模拟算例,考察基于小波数值方法的精度及其自动捕捉强脉动的能力,为 后续实现三维湍流数值模拟积累一些经验. 小波系数临界值确定的方案选择对于湍流的小波自适应模拟的效果和效率非常关键. [5鄄6] Farge和Schneider 基于去噪理论获得临界值公式,提出相干涡模拟(coherent vertex simula鄄 * 收稿日期:摇 2010鄄07鄄02;修订日期:摇 2010鄄10鄄29 基金项目:摇 创新研究群体基金资助项目;973资助项目(2010CB832700) 作者简介:摇 赵勇(1981—),男,博士生(联系人.E鄄mail:fluid@mail.dlut.edu.cn); 宗智(1964—)男,教授,博士(E鄄mail:zongzhi@dlut.edu.cn); 邹文楠(1968—),男,教授,博士(E鄄mail:zouwn@ncu.edu.cn). 33 34 赵摇 摇 勇