顶盖驱动流的降阶模型及POD分析.pdf

中国工程热物理学会 传热传质学． 第十二届学术会议 ’． 编号：063197 顶盖驱动流的降阶模型及POD分析 吴学红 丁鹏 陶文铨 矿q安交通大学能源与动力：I：程学院，西安，710049 Tel-029 E—哪ail：wqtao@mail．xjtu．edu．ca 摘要：本文主要采用POD方法对顶盖驱动流的模拟结果进行分析；弗用Oalerkin方法对N-S方程进 行重构建，得到关-3二N．S方程的降阶模型。研究表明：这种方法能很好的捕捉到流动的特征，并且 降阶模型的计算结果与直接数值模拟的结果符合的很好，但节省了更多的计算时间。 关键词：POD方法。Galerkin投影，降阶模型，顶盖驱动流，数值模拟 l、引言 自从上世纪六十年代以来，POD(ProperOrthogonalDecomposition)方法作为～种 分析复杂物理问题的工具，越来越多的得到了重视，并在传热与流动中已经取得了广泛 的应用。在理论上，POD方法足基于统计方法对物理场进行分析。最初在1967年Lumley et 【1】提出POD方法是用来分析紊流场的相干结构。W．Cazemieral‘21用它对Re=22000 的顶盖驱动流紊流模拟结果进行分析，，研究结果显示此方法能很好的提取流动中的相干 结构；并构建出一种降阶的动力学模型，并用它对Re=8000～12000进行求解，计算的 结果与直接数值模拟的结果符合的很好。Ly和Tran”1用它对cVD化学反应器的模拟结 果进行分析，并用Galerkin投影方法对传质方程进行重构建，得到传质方程的降阶模 型。Ly和Tran[4]用此方法对稳态方腔自然对流进行分析，POD方法得出的结果与直接数 值模拟的结果误差在3％以内。 本文主要用POD方法对二维层流稳态顶盖驱动流的数值模拟结果进行分析：并用 Galerkin投影方法对N．S方程重构建进行研究，寻找一种降阶模型，达到能够减少计算 时间的目的。 2数学物理模型及计算方法 计算的物理模型如图1所示，K宽均为L，左右边界及下边界为无滑移边界条件， 上边界给定速度u，计算中不考虑重力的影响。计算的无量纲控制方程如下l，】： V．旷：0 f (1) H ) (细)旷一即+去V2旷 式中的无量纲参数如下： ty’v l Lou ∥=v／‰，X=x／L， a ● 尸=ptpu2,Re=吃LIv 图1顶盏驱动流 基金项目：国家自然科学基金资助项月 1097 控制方程采用有限容积法进行离散，速度与压力采为交错网格的SIMPLIC算法进 行耦合，计算的网格数为52×52。 3构建POD基函数 关于POD的数学描述见文献【3】，从文巾可以得出POD方法主要分为以下四步： Re数每增加100记录一个速度场数据U，把计算所得的数据进行求和、F均，并用每一 个速度场减去这个jr均值，可以得到一组相对丁-jF均速度的脉动值V，也即得到了30 个Snapshots。 j-l∥2—30 (2) D=U=去∑u产(x)·v=U，一D， JV ltI 第二步