基于小尺度方程网络模式的大涡模拟方法研究.pdfVIP

基于小尺度方程亚网格模式的大涡模拟方法 ‘ 吴锤结’,2 关 晖“赵红亮2 (’北京大学湍流研究国家重点实验室，北京100871; 2中国人民解放军理工大学流体力学研究中心，南京211101) 摘要 本文提出大、小尺度相互辆合计算的大涡模拟方法，其中小尺度计算采用荃于最优基的低维动力系 统.通过采用基于小尺度方程亚网格模式的大涡棋拟数值计算了三维借道湍滚，并将结果与Smagai麟》模 式和动力学模式相比较.这一新方法的主要优点是 不包含任何可调参数;可适用于各类复杂湍流流场计 算;可得到正确的亚网格流场信息，包括大、小尺度间的逆向能量传输等;可获得其他方法无法得到的近似 亚网格流场，并可分析其动力学特性. 关扭词 大涡模拟-POl〕方法:动力系统 0 引 言 湍流问题是经典物理学中最后一个未被攻克的堡垒，同时又对国民经济、人类活动各方面影响极 大.湍流数值摸拟是研究、揭示湍流基本规律的重要工具，对湍流大涡模拟 L〔ES)的研究既可以克服 因计算机能力不足所造成的困难，又可以加深对揣流本身的认识，揭示湍流的基本规律.大涡模拟的 关键是给出正确的亚网格(SGS)模式. 通常，湍流大涡模拟方法包含如下步骤:对Navier-Stokes方程和连续性方程进行滤波 r . f ( X . f(x)+f(x) . J I J f ( (1) X =}(二‘，二，)f(二，)dx 得到大涡模拟方程 丛 两 ,32又 +”dxjdx,ad=x}j (2) 其中，G(x,x)为给定核函数.对上式中的揣流亚网格应力项rd二u;可一可可必须进行模拟，建立 SGS模式.同时在计算过程中还必须对边界条件等进行特殊处理，上述Navier-Stok。方程的大涡模拟 方法基于如下假设:亚网格小尺度结构比大尺度结构更为一般，对边界条件等的影响不敏感，因此， 更便于建立模式. 目前，SGS模式主要有:涡粘模式(如Smagorinsk)模式川和动力学模式[[2,31)，尺度相似模式a[[,5] 速度估计模式[6.，基于涡方法的SGS模式[81，重正化群模式}9,10,111，及其他模式，如Kraie11rian的 EDQNM方法{12]，M81arst,P5lenr的结构函数模式[131,Jim6二的惯性律模式!141，Reynolds应力模式[15] 中国科学院力学研究所非线性力学国家重点实验室资助项目 t_If抓日 及最近由Spalart提出的DES(Detached-EddySimulation)方法16〔一等.亦可参见1〔7] 我们认为，LES研究中的主要问题在于对湍流本质认识不足，所采用的SGS模式无法适用于各类 复杂揣流.其具体表现为，在各种SGS模式中大都存在着可调参数，并且对边界附近流动依赖于壁面 率的使用.因此，通过研究湍流的物理和数学本质，提出在数学上严格、物理上合理、适用于各种复杂 湍流的、不存在任何可调参数的SGS模式理论，对揣流数值模拟及揣流本质的认识都将是重要的贡 献.为此，我们提出大、小尺度相互藕合计算的、基于Snapshots-POD基SGS模式的大涡模拟新方法. 1 小尺度方程模式的基本思想与解法 1.1基本思想 小尺度方程模式的基本思想是，首先通过滤波将大、小尺度分离，得到大涡模拟方程和小