广义HWD-POD方法及湍流耦合低维动力系统.pdfVIP

广义HWD．PoD方法及湍流耦合低维动力系统 吴锤结1t2赵红亮2 (‘北京大学湍流研究国家重点实验室，北京lo嘴7I； 2中国人民解放军理工大学流体力学研究中心，南京2lll0¨ 摘要吴锤结(昊锤结与赵红亮，19卯)提出基于最优广义螺旋波基低维动力系统建模方法，试图通过分析 其特性，以最大限度地揭示湍流涡结构的动力学演化规律．奉文根据涡矢量国的手征性将流场分为一对螺 旋量符号相反的流动数据库，分别对两组数据库进行POD分析，得到广义Hwnp0D基；在其上建立了耦合 低维动力系统，并对其三维复杂流动的空同特性和动力学特性进行分析 关键词螺旋波基；动力系统；湍流 0引言：螺旋波基——湍流的基元 我们认为湍流基元至少必须具有如下性质：在动力学方面，它必须深刻反映湍流的动力学特性； 上，它必须具有集聚性，且在远场衰减得足够快；在时间演化方面，必须具有缓变性，即：基元的时间 演化比流场脉动过程慢得多．通过缓变基元间的非线性相互作用导致流态的瞬息万变． 经过多年研究，人们认识到湍流拟序结构中的涡结构——涡元——对湍流的产生和发展均起着 至关重要的作用．但至今为止，对涡本身还没有公认的统一定义．鉴于认识的不统一，加上涡元的时间 演化与流场脉动往往是同一量级的(不具有缓变性)，我们认为涡元不是数学上精确、物理意义明确、 实际可操作的湍流的基元． 间中定义螺旋波基为 『y}(置，J)=(cos七·工)J士(sin七一工)f iwft，．r)：f。jn量．JH干(。。。七．J)f ‘1’ 它们分别满足 fVx=±^ 1vx聪：±^y} ‘2’ 对应的螺旋量分别为±^．它们均为满足Euler方程的Beltmmi流． 现：涡量场的快速变化是由一对极性相反的、缓变的HwD涡相向运动，通过线性迭加而成的；在最大 涡与最小涡之问存在相互作用；并且观察到小尺度极性相同的HwD涡相互配对的反级串过程，它是 大尺度拟序结构内在间歇性的基本原因之一． 可见}fwD具有如下优点：它深刻反映了涡动力学特性，是揭示其物理本质的最自然、最基本的 正交分解；mm基具有空问集聚性、时间缓变性，可作为湍流的基元．但是，经典的}1wD基具有如下 ·7D- 缺点：F叭i”空间中的}IWD无法用于复杂几何、物理边界情形；对湍流低维动力系统建模而言，它不 是晟优的． 1广义HwD．P(m数据库分析方法 结合p。o数据库分析方法，求出满足mm基部分基本性质的广义m巾．PoD基．具体做法如下： 首先，根据A的符号将流场分为满足手征性的Ⅳ+，Ⅳ一两组数据库．设一已知速度场Ⅱ’(z，I)，其 stokes．Hd出0ltz分解为 H：Ⅳ7一甲乒=Ⅱ++“一 相应的涡量场∞(z，￡)为 ∞=V×口’=可ד=V×口++V×Ⅳ 根据mm分解，即： V×{}=馗{ (3) 对上式两边同时点乘{f，得 {l·V×{}=鸠￡·{f (4) 因为{}·{}≥0，所以A的符号与流场中各点处螺旋量的符号一致． HWD．P(ID基． 2广义HWmp0D基的数值分析 以标准三维槽道流的湍流大涡模拟数据库作为分析对象．数据库的网格数为32×32×32，基于 壁面摩擦速度蜥和半槽宽△的ReyT·ol出数为3∞，流动数据的时间范围为无量纲计算时间O．5～1．5+ 共计200卟计算步长．因在三维槽道流计算时已扣除势流部分，所以它可直接分解为Ⅱ+与Ⅳ一阿部 分．应用广义HwD．POD方法即可求得两组广义m∞．POD基．下面考察其所具有的独特性质． 2．1 广义HwD．P()D基俘获能量的