基于本征正交分解的并列双圆柱绕流尾流场分析.pdf
《基于本征正交分解的并列双圆柱绕流尾流场分析.pdf》由会员分享,可在线阅读,更多相关《基于本征正交分解的并列双圆柱绕流尾流场分析.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 卷 第 期 年 月气 体 物 理 .:./.基于本征正交分解的并列双圆柱绕流尾流场分析王军辉 田书玲(南京航空航天大学航空学院空气动力学系 江苏南京)()摘 要:并列双圆柱绕流尾流场流动复杂 在不同的间距下会呈现多种复杂尾迹模式 基于有限体积方法 针对低 数 下不同间距并列二维双圆柱绕流流动开展了数值模拟研究 并采用本征正交分解()法对仿真得到的尾流场进行模态对比分析 分析结果表明:间距.时 上下两圆柱间干扰较弱 尾迹呈现周期性同步同向模式 间距.时 上下两圆柱间干扰较强 流动现象较为复杂 尾迹呈现交替翻转模式 流场不再具有明显周期性规律 同时 方法能得到尾流流动在空间上的主要分布特征 较
2、少的模态包含流场的绝大部分能量 且间距.的尾流模态结构较间距.的模态结构更加紧凑有序这和不同间距 下的尾流流动形态密切相关关键词:圆柱绕流 尾迹 数值模拟 本征正交分解 模态 中图分类号:.文献标志码:.().:收稿日期:修回日期:第一作者简介:王军辉()男 硕士 主要研究方向为复杂尾流场降阶分析:.引用格式:王军辉 田书玲.基于本征正交分解的并列双圆柱绕流尾流场分析.气体物理 ():.:.():.引 言具有复杂流动分离现象的单个或多个钝体绕流现象广泛存在于航空航天等工程应用领域 因此准确描述和理解其流动机理对工程研究来说至关重要 早在 世纪中下期 等、等等通过实验手段对并列双圆柱绕流流动现象
3、开展了一系列研究 研究结果表明 流动分离现象随着两圆柱距 的改变 尾涡脱落也会随之改变研究了 下并列双圆柱绕流的尾迹脱落模式 并讨论了 种不同尾迹现象及流动机理 年 聂德明等采用不可压格子 模型对不同 数下并列双圆柱绕流进行了更为细致的数值模拟研究 分析了间距 以及第 期王军辉 等:基于本征正交分解的并列双圆柱绕流尾流场分析 数对圆柱尾流及升力、阻力的影响 并在此基础上得到了 种尾迹模式 同时研究了流场的初始扰动对流动分岔现象的影响 发现在适当的扰动下可以很快得到同步同相的尾涡脱落模式随着计算机性能的提高以及高效算法的不断提出 复杂流动现象的模拟仿真研究成为了可能然而要准确提取复杂流场流动结构
4、 研究流动内在机理仍然具有一定的难度 近年来 本征正交分解()方法作为一种流场结构提取方法被国内外一众学者广泛使用 方法是目前常用的一种高效数据降维方法 其核心是在空间域中寻找一组最优标准正交基 样本数据在该正交基上的投影迅速递减 保留投影中含能较高的前几阶模态 使用较少的模态基就能对高维原始样本数据进行低维近似描述 方法的本质与数据处理中的主成分分析()算法大致趋同 其最早是由 于 年提出来的一种高效数据降维模态分解方法 随后 国内外学者将其应用于流场结构的分析中 并在原有算法的基础上进行优化 年 首次将 方法应用于湍流相干结构的提取中 其基于时间序列的流场空间样本数据 得到了前几阶含能较高
5、的流动模态以此为基础研究了流场的流动机理 随流场数据更加复杂化与高维化 基础算法对于高维流场空间矩阵难以处理 随后 提出了快照 算法 是对原始 算法的优化与改良其核心在于将原始 算法的空间维度矩阵通过矩阵变换转换到时间维度上 这样大大节省了计算时间 解决了原始 算法难以求解超大型维度矩阵的问题 等采用快照 方法针对机翼翼型设计优化问题建立了降阶模型 大大降低了不同翼型的数值仿真计算成本并列双圆柱绕流是流体力学中常见的流动问题之一 本文首先对两种不同间距下并列双圆柱绕流尾流的形成与发展过程进行了非定常数值模拟 分别研究了低 数下两圆柱之间间距 对尾涡脱落模式的影响 然后使用快照 方法对仿真得到的
6、非定常涡量场数据进行模态分析 得到了各阶模态所占流场能量大小以及各阶模态所对应的流场结构 最后分别使用前几阶高能占比模态对两种不同间距 双圆柱绕流涡量场进行了重构 验证了快照 方法在分析复杂流场内在流动机理方面的作用 数值模拟.数值求解方法二维非定常 方程的积分守恒形式为()()()式中 为控制体 为控制体单元边界 为守恒变量()为无黏通量()为黏性通量 采用有限体积法对式()进行离散 如下所示 式中 为流体密度 为单元体积 为单元面的个数 为面的面积 为面的质量流率 为面 的梯度 空间离散格式采用 阶迎风格式 时间离散项采用 阶隐式格式进行离散 即()式中 为当前时间步长 为前一个时间步长
7、为下一个时间步长、分别为当前时间步长的标量值()为残差项 由于本文计算低 数圆柱绕流问题 故采用层流模型进行数值模拟.双圆柱绕流数值模拟计算 区 域 如 图 所 示 两 圆 柱 直 径 均为 为两圆柱间间距 计算模型区域大小为.计算域网格采用多块结构化网格对两中心圆柱采用 形网格进行加密 两圆柱间距.和.的网格单元数均为.计算边界条件如下所示:入口:出口:/上下边界均为对称边界:/两圆柱面:图 计算区域及边界条件.气 体 物 理 年 第 卷采用有限体积法进行数值模拟 为保证计算精度 非定常计算时间步长为.来流 数 选用层流模型 计算总步数设置为 步 即分别模拟 内两圆柱间距 .以及.的双圆柱非
8、定常绕流流动过程 方法 方法的核心思想是在最小二乘意义下从时间序列的空间数据样本中寻找最能代表样本数据的一组最优正交基函数来 使用较少的正交基函数来对原始样本数据进行低维近似 具体原理推导如下:令 为流场变量 如速度、温度、涡量等 由指定的时间点获得的快照集合()()构成如下的 样本矩阵 其中 为网格节点数 为快照数快照集合的时间均值可表示为 快照脉动矩阵 其中 方法的核心是寻找一组最优正交基函数()来最大化下式 满足 其中、以上最优问题则转化为以下的特征值问题其中 ()()本文采用 提出的快照 方法求解 模态 该方法主要是将所要求解的矩阵 的空间维度转换到时间维度上 从而利用时间维度矩阵特征
9、值线性组合表征特征模态 即 令 即求 式中 为对应特征值的特征函数 由于所求特征值维度只有 阶 故与原始 算法所要计算的 阶矩阵特征值问题相比极大减小了计算量 计算结果及 分析.尾迹模式及升阻力系数图 为双圆柱间距.、.对应的尾流涡结构图 图 为.、.对应尾流流线图图、图 分别为双圆柱间距 .、.对应的升力系数和阻力系数随时间的变化曲线图 图 中两圆柱升阻力系数说明 当两圆柱间距.时 尾涡脱落无明显周期性规律 两圆柱的阻力系数随时间不断交替变化 两圆柱间流场互相干扰较强 此种间距下的圆柱间流场在某一个时间段偏向某个圆柱则随时间推移慢慢又偏向另外一个圆柱不断交替循环 可以看出偏向的圆柱受尾流影响
10、其阻力系数明显大于另一圆柱 再结合图 的涡结构以及流线图可得 偏流的存在使得流场偏向的圆柱后方形成的旋涡较为窄小 此时尾迹脱落模式是以交替翻转模式进行 图 中两圆柱升阻力系数说明 当两圆柱间距.时 尾涡脱落具有明显周期性规律 两圆柱的升力系数随时间变化具有相反的趋势 二者升力系数曲线存在 的相位差再结合图 的涡结构以及流线图可得 旋涡在两圆柱表面以反向模式进行 圆柱后上下尾流涡关于中轴线对称 在圆柱后方很长的一段区域内 尾迹流场结构仍然保持着不变 此时尾迹脱落模式是以同步反向模式进行对比两种不同间距.、.尾迹涡结构以及升阻力系数不难得出 随两圆柱间距 发生改变 尾迹模式也有着相应变化 当间距
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 正交 分解 并列 圆柱 绕流尾流场 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。