建筑结构周围流场分析的一种新型数值模拟方法研究.pdfVIP

第十四届全国结构风工程学术会议论文集 建筑结构周围流场分析的一种新型数值模拟方法 周岱 何涛 包艳 (上海交通大学船舶海洋与建筑工程学院，上海200240) 摘要针对建筑结构风场以及水体流动等不可压缩粘性牛顿流体流场的数值模拟，提出一种基于双线性和线 隐式化并进行空间离散，运用三步有限元法实施时间离散．通过对不同几何形状腔内的流体流场的数值模拟显示， 本方法的计算精度和稳定性好，同时与其他数值算法比较，计算效率更高． 关键词半隐式三步法、SUPG有限元法、非定常不可压缩流、腔内拖拽流 1．引言 风场流体域的高效、高精度数值计算是建筑结构风场模拟的关键问题之一。计算流体动力学 是，工程中常面临求解高雷诺数流场问题，而N．S方程中存在非线性对流项，致使采用传统Galerkin 流体有限元法难以获得稳定的数值解。为克服这一困难，许多学者提出各种求解流体流场的稳定化 通过引入SUPG稳定化技术，突破迭代求解时对时间步长的限制，从而提高计算效率，同时保持满 意的计算精度和稳定性。 2．不可压缩Navier-Stokes方程 程和连续方程为： 一‘一dui．{-ubut，，=呀，．，+Z，j20。，j+ji U】 tubu-{tO)1( ui,f=0 (2) 式中，“，”表示求导，嘶是f方向的速度分量，P是流体密度，石为外力，t是时间；a0为应力张量， 由下式给出： 吩=一PSij／p+vlui，，+叶，fJ (3) 式中，P和D分别为压力和动力粘度，磊，为Kronecker符号。边界条件为： 坼=gt， 在0上 (4) 在n上 (5) q，nj=tl,， 式中，厂g Fh为不相交二维边界厂的非重叠子段，gf和hi为给定位移和给定表面力，nj为边界 ．，方向的法向单位向量。 3．半隐式三步法 运用半隐式三步法对控制方程进行时间积分，需显式化对流项，而隐式化粘性项。因此动量方 程在每个时间子步，离散为如下形式： 《≯=母扩∥屈+y伽护(1嘞Ⅶn+l，3)+∥(6a) (6b) 910 第十四届全国结构风工程学术会议论文集 (6c) 竺：=!云j筹竺=一Ⅳ；剧3甜嚣2，3一刀+1屈+yL岛嵋箩73+(1一历甜搿)+∥捌3 式中，△f是时间增量，上标刀为时间步数。文中对粘性项采用二阶隐式梯形法积分，且取∥=o．5， 这样时间步长取值可不受粘性条件限制。 4．控制方程的有限元离散 4．I Navier-Stokes方程的SUPG有限元格式 将动量方程(1)乘以速度检验函数，在整个计算域内积分，并加上SUPG稳定项，则其弱式为1)： L毛u粥越，钆√鲁+％％一％一z)搠=。(7) L【跳(鲁+巧％一z)+％跳√】扣一￡。^她矿+兰e=l 式中，net为单元数，瓯是速度检验函数。龟uPG为稳定化参数，定义为13)’14)’15)： h8 ％UPG ～：譬 (8) 2研死弘 其中，h。、t／e和Re8分别为单元特征尺寸、单元