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

基于小尺度方程亚网格模式的大涡模拟方法’ 昊锤结1，2关晖2赵红亮2 (1北京大学湍流研究国家重点实验室，北京100871； l 2中国人民解放军理工大学流体力学研究中心，南京2ll01) 摘要本文提出大、小尺度相互辆合计算的大涡模拟方法，其中小尺度计算采用基于最冼基的低维动力系 统．通过采用基于小尺度方程亚网格模式的大涡模拟数值计算了三维槽道湍流，并将结果与陆n鲫r吣模 式和动力学模式相比较．这一新方法的主要优点是．不包含任何可调参数；可适甩于各类复杂湍流流场计 算；可得到正确的亚网格流场信息，包括大、小尺度阃的逆向能量传输等；可获得其他方法无法得到的近似 亚同格流场，并可分析其动力学特性． 关奠词太涡模拟；PoD方法；动力系统 O引 言 湍流问题是经典物理学中最后一个未被攻克的堡垒，同时又对国民经济、人类活动各方面影响极 大．湍流数值模拟是研究、揭示湍流基本规律的重要工具，对湍流大涡模拟(LEs)的研究既可以克服 因计算机能力不足所造成的困难，又可以加深对湍流本身的认识，揭示湍流的基本规律．大涡模拟的 关键是给出正确的亚网格(sGs)模式． 通常，湍流大涡模拟方法包含如下步骤：对Na“er．stok*方程和连续性方程进行滤波 r，(j)=，(，)+，’(z) 17(J)：fG(¨，)“一)“ (” 得到大涡模拟方程 aM。 au。虫 l a卫 护u。 强 瓦+百一i茬”硒a即 (2) 挚：o dZJ 其中，G(J，一)为给定核函数．对上式中的湍流亚网格应力项o：司孑一哥可必须进行模拟，建立 scs模式．同时在计算过程中还必须对边界条件等进行特殊处理．上述Navi—s10ks方程的大涡模拟 方法基于如下假设：亚网格小尺度结构比大尺度结构更为一般，对边界条件等的影响不敏感，因此， 更便于建立模式． 目前，sGs模式主要有：涡粘模式(如蛳dmky模式…和动力学模式【2—3])，尺度相似模式【4，“， 速度估计模式№，基于涡方法的sGs模式【…，重正化群模式㈨o·“】，及其他模式，如Klaic}㈣的 * 中国科学院力学研究所非线性力学国家重点实验室资助项目 ——一——————虿■孺蕊—— simulati∞)方法L16-等．亦可参见’17]． 及最近由spal越提出的D岱(DeIached．Eddv 我们认为，Ⅱs研究中的主要问题在于对湍流本质认识不足，所采用的sGs模式无法适用于各类 复杂湍流．其具体表现为，在各种sGs模式中大都存在着可调参数，并且对边界附近流动依赖于壁面 率的使用．因此，通过研究湍流的物理和数学本质，提出在数学上严格、物理上合理、适用于各种复杂 湍流的、不存在任何可调参数的sGS模式理论，对湍流数值模拟及湍流本质的认识都将是重要的贡 献．为此，我们提出太、小尺度相互耦合计算的、基于sn印sIlots．PoD基sGs模式的大涡模拟新方法． 1小尺度方程模式的基本思想与解法 I．I基本思想 小尺度方程模式的基本思想是，首先通过滤波将大、小尺度分离，得到大涡模拟方程和小尺度亚 网格流场的控制方程；小尺度亚网格模式为基于等网格的小尺度方程，利用snapsho*PoD方法构造 满足亚网格流场边界条件的最优基，将小尺度方程在最优基上投影而得到小尺度动力系统；将其与大 涡模拟方程耦合求解，得到近似小尺度速度场，进而得到大涡模拟方程中的RevI“ds应力项，这样即 可构成新的大涡模拟方法．构造能反映小尺度场动力学特性的基函数——小尺度场的拟序结构—— 是该方法的关键所在． 1．2解 法 首先，通过滤波将流动变量，(z)分解为大尺度量7(r)与小尺度量，7(j，．这里，滤波仅在zy 平面内进行．对不计质