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

计算物理 ›› 2015, Vol. 32 ›› Issue (1): 13-19.

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

孙晨扬, 李启良, 杨志刚

同济大学上海地面交通工具风洞中心, 上海 201804

收稿日期:2013-12-25 修回日期:2014-04-01 出版日期:2015-01-25 发布日期:2015-01-25
作者简介:孙晨扬(1991-),男,江苏兴化人,硕士生,主要研究车身与空气动力学,E-mail:sun-chenyang@163.com
基金资助:
国家自然科学基金青年科学基金(11302153);同济大学中央高校基本科研业务费专项资金(20123315)资助项目

Least-squares Finite Element Method for Unsteady Stress Formulation of Navier-Stokes Equations

SUN Chenyang, LI Qiliang, YANG Zhigang

Shanghai Automobile Wind Tunnel Center, Tongji University, 201804 Shanghai, China

Received:2013-12-25 Revised:2014-04-01 Online:2015-01-25 Published:2015-01-25

摘要/Abstract

摘要： 为精确求解非定常层流问题,发展一种非定常速度-应力-压力的方法.采用牛顿法对非线性对流项进行线性化处理和预处理共轭梯度法,实现了非定常应力形式Navier-Stokes方程的求解.方腔层流流动比较发现,非定常应力形式比涡量形式与试验结果更加吻合,精度更高.该方法有效地解决亚格子应力项的问题,实现基于最小二乘有限元法的湍流求解.比较方腔湍流流动的试验与仿真结果,证明本文的方法具有可行性,为湍流大涡模拟计算打下基础.

关键词: 最小二乘有限元法, 非定常速度-应力-压力, 大涡模拟, 方腔流

Abstract: To solve unsteady laminar flow problems,a method of velocity-stress-pressure formulation instead of velocity-vorticitypressure formulation is developed.With Newton's linearized method to linearize convective terms and preconditioned conjugate gradient method to solve equations,unsteady stress formulation of Navier-Stokes equations is solved.Comparison between numerical and experimental results of cavity laminar flow shows that result of stress formulation fits experiment better and has higher accuracy than vorticity formulation.The stress formulation can deal with subgrid stress with least squares finite element method.Comparison with experimental results of cavity turbulent flow reveals feasibility of the method.It lays a firm foundation for large eddy simulation computation.

Key words: least squares finite element method, unsteady velocity-stress-pressure, large eddy simulation, cavity flow

中图分类号:

TP301.6

孙晨扬, 李启良, 杨志刚. 最小二乘有限元法求解非定常应力的Navier-Stokes方程[J]. 计算物理, 2015, 32(1): 13-19.

SUN Chenyang, LI Qiliang, YANG Zhigang. Least-squares Finite Element Method for Unsteady Stress Formulation of Navier-Stokes Equations[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2015, 32(1): 13-19.

[1]	刘梅, 李曙光, 吴正人, 王松岭, 邓宇涵. 带沟槽结构的平板流动中的熵产分析[J]. 计算物理, 2020, 37(2): 182-188.
[2]	方乐, 王楚涵. 理性亚格子建模的实践与反思[J]. 计算物理, 2018, 35(3): 253-261.
[3]	潘宏禄, 李俊红, 程晓丽, 马汉东. 气动光学效应湍流脉动模型系数修正[J]. 计算物理, 2018, 35(2): 194-204.
[4]	李为君, 张云伟, 顾兆林, 段翠娥, 张力元, LU Weizhen Jane. 时变来流条件下大气边界层流动大涡模拟[J]. 计算物理, 2016, 33(6): 691-697.
[5]	上官燕琴, 王娴, 李跃明. 基于格子Boltzmann方法的平板射流大涡模拟[J]. 计算物理, 2015, 32(6): 669-676.
[6]	潘宏禄, 关发明, 袁湘江, 卜俊辉. 高超声速平板边界层/圆柱粗糙元非定常干扰[J]. 计算物理, 2015, 32(5): 537-544.
[7]	刘同新, 马宝峰. 低阶格式在大涡模拟计算中的适用性[J]. 计算物理, 2014, 31(3): 307-313.
[8]	洪宁, 李体云, 李秋香, 施保昌. 两相方腔流的格子Boltzmann模拟[J]. 计算物理, 2013, 30(4): 520-526.
[9]	陈瑜, 夏振华, 蔡庆东. 基于树网格的格子Boltzmann方法以及曲线边界的处理[J]. 计算物理, 2010, 27(1): 23-30.
[10]	邱剑, 顾兆林, 王赞社. 一种改进后的亚格子特征尺度表达式[J]. 计算物理, 2007, 24(6): 711-716.
[11]	谢志刚, 许春晓, 崔桂香, 王志石. 方柱绕流大涡模拟[J]. 计算物理, 2007, 24(2): 171-180.
[12]	许春晓, 吴超, 崔桂香. 壁面在展向作周期运动的槽道湍流的大涡模拟[J]. 计算物理, 2006, 23(5): 537-544.
[13]	王小华, 何钟怡. 肋片对后台阶湍流流动影响的数值模拟[J]. 计算物理, 2005, 22(4): 344-350.
[14]	夏健, 刘锋. 各向同性紊流能量衰减的大涡模拟[J]. 计算物理, 2005, 22(1): 61-64.
[15]	崔桂香, 周海兵, XU Chun-xiao, 张兆顺, L. Shao. 新型大涡数值模拟亚格子模型及应用[J]. 计算物理, 2004, 21(3): 289-293.