三维粘性流体内流问题的边界元法的研究中求解奇异积分的一种新的数值解法.pdfVIP

三维粘性流体内流问题的边界元法研究中 求解奇异积分的一种新的数值解法 余流王铁成高永祥刘兴业 (天津大学土木工程系，无律300072) 提要：本文探讨了三维粘蛀流体内流问题的边界元法研究中奇异积分的一种有效的解法． 对边界积分项，采用三角形极坐标来降低奇异积分的维数．从而将整体坐标系下的三角形 单元转换成局部坐标系中的单位正方形单元：肘域嘏蔓卜项．采用四面体樱坐标，不但降髂 了奇异积分的维数，而且将整体坐标系下的四面体单元转换成局部坐标系中的单位立方体 单元，从而使积分域简单化．最后利用高斯积分进行数值求解． 美键词：粘性流体奇异积分三角形极坐标四面体极坐标高斯积分 1前言 不可压缩的粘性流体运动遵循质量守恒定律和牛顿定律，其速度场矿和压力场卢分别 为连续性方程和N--S方程。J．C．Wu“…””“”将流体流动问题划分成运动学和动力学两个 方面，同时以涡量场来代替压力场，利用N-s方程和连续性方程推导出三维不压缩粘性流 体速度场和涡量场积分表达式 _=一击f坠≠‰+击{盟盐等峨㈤ 氍=击I业半丝搬+击{竖宁‰+ m 一1 c[坠虹兰J二譬!!!生坠兰!招 Ir=|F一昂I-【(x一‰)2+(J，一yo)2】+(：一：。)1】2 JF=x7+万+石 l五=x,i+“』+z．k 【s=nji+■，+，‘t m一1j，k三向求和 1，jk按如下轮换可得另外二方向的速度、涡量 i——_J——_t——“ 案漉男74年6月生天津大学建筑工程学院博士生 ·608· ‘工程力学)增刊2001年 2积分表达式离散 2．1离散单元 体积分域划分成四面体单元，边界域(面积分域)划分成三角形单元。 2．2线性分布假定 假定流体内任一点的速度、涡量、总水头为该点体积坐标的线性函数 a．四面体单元 在四面体单元中，引入体积坐标有： 矿 矿 茧；号，…，六=罟，点+{：+￡+￡=1 (3) ．． b．三角形单元 (4) 三角形单元内，引入面积坐标：量=l，--邑=1．点+f：+参=1 j j 最后，积分方程(1)(2)离散为： 吒=∑∑∑吲。舻+∑∑∑州。卯’ (5) Jd r=lp4～1 ，d口BI 厣：：芝主主暇。爿∥+E…Zh：。Ⅳ，。+兰宝主辟訾，且P (6) eel‘2-1，一1 ·，l 0-J o‘Iv…。I ned、neb一体单元、面单元数 3特殊极坐标m 三角形极坐标 将平面三角形变换成单位正方形。在边界元法中，边界单元采用整体坐标系中的空间三角 形单元，单位正方形采用平面局部坐标系，则变换公式为： z，=(1一户；)Jj”+P．(1一p：)xj2’+p，p，xj”，，=f，，，女 (7) 此处，t，，：j，J，^为正交的直角坐标，P．、岛为三角形极坐标，外p：与三角形面积 s毒、乒、￡≤l，点+邑+轰=1)的关系为： 坐标舌、彘、苣(0 (8) f。=1-P．，f：=p，(1一P：)’蠡=p。P： 四面体极坐标。 7，坐标变换公式为： 四面体单元F经过坐标变换．转换成单位边长的立方体F