最小二乘有限元法求解非定常应力的Navier-Stokes方程.pdfVIP

第32卷 第 1期 计 算 物 理 Vo1．32，No．1 2015年 1月 CHINESE JOURNAL OFCOMPUTATIONAL PHYSICS Jan．，2015 文章编号：1001—246X(2015l01-0013-07 最小二乘有 限元法求解非定常应力的 Navier-Stokes方程 孙晨扬 ， 李启 良， 杨志刚 (同济大学 上海地面交通工具风洞 中心 ，上海 201804) 摘 要 ：为精确求解非定常层流 问题 ，发展一种非定常速度一应力一压力的方法．采用牛顿法对非线性对流项进行 线性化处理和预处理共轭梯度法，实现了非定常应力形式 Navier．Stokes方程的求解．方腔层流流动 比较发现，非定 常应力形式 比涡量形式与试验结果更加吻合 ，精度更高．该方法有效地解决亚格子应力项的问题 ，实现基于最小二 乘有限元法的湍流求解．比较方腔湍流流动的试验与仿真结果，证明本文的方法具有可行性，为湍流大涡模拟计算 打下基础． 关键词：最小二乘有限元法；非定常速度一应力一压力；大涡模拟；方腔流 中图分类号 ：TP301．6 文献标志码 ：A 0 引言 流体流动均为非定常流动，包括层流和湍流．为有效模拟流体流动，常用的有限差分法和有限体积法均 发展了相应的非定常形式 ¨ ．最／b-乘有限元法作为与之相对应的方法 ，已发展了非定常的速度 一涡量 一 压力形式 ，并在层流流动中有所应用 ． 近年来基于有 限体积法 的大涡模拟在工程和研究中得到广泛应用．然而 ，这些计算结果均与试验结果差 距较大 J．主要原 因是有限体积法的精度 限制在二阶．而最tJx--乘有限元法具有通用性强、计算效率高、 收敛性好和计算准确并且能够使用高阶精度等特点 ，成为获得精确湍流结果的新方法． 由于传统的速度 一涡量 一压力形式无法处理大涡模拟 出现的亚格子应力项 ，本文首先发展非定常的速 度 一应力 一压力形式 ，在此基础上 ，引入亚格子应力 ，并 由此发展可以求解湍流的方法．通过方腔流动的数值 计算与试验对比，评估本文方法的准确性 ，为最小二乘有限元方法在湍流计算 中的应用提供支撑． 1 数值方法 1．1 非定常应力形式的N．S方程 非定常不可压无量纲应力形式 的N-S方程 ， Ou j — — ： 0． a 塑 + 十op at ax 眠a％；鲁Ox；～ (1) 式中，ll为速度 ，p为压力 个未知量 ，三维问题有 10个未知量． 收稿 日期 ：2013—12—25；修 回 日期 ：2014—04—0l 基金项 目：国家 自然科学基金青年科学基金；同济大学中央高校基本科研业务费专项资金资助项 目 作者简介 ：孙晨扬(1991一)，男，江苏兴化人，硕士生，主要研究车身与空气动力学 ，Em·ail：sun-chenyang@163．com 14 计 算 物 理 第 32卷 1．2 大涡模拟控制方程及亚格子应力模型 按照大涡模拟理论 ，经滤波函数处理后的不可压缩大涡模拟控