圆外Stokes问题的边界积分公式研究.pdfVIP

固体力学学报 S．Issue 第27卷专辑 V01．27 ACTAMECHANICASOLlDASINICA 2006 2006年12月 December 圆外Stokes问题的边界积分公式。 彭雏红 董正筑 赵慧明 曹国华 (中国矿业大学理学院，徐州，221008) 擒要 针对圆外区域的Stokes问题．利用Fourier展开法，通过自然边界归化得到了一个只与边界速度有 关的Stokes问题的边界积分公式．根据此公式及边界速度值，求得区域内速度及压力分布的解析表达式，并通过数 值积分的方法进行求解，计算量小，所得速度及压力分布图曲线光滑．最后借助流体软件进行数值计算，结果验证 了边界积分公式的正确性，可行性． 关■硼 圆外区域，Stokes同题．Fourier晨开法，边界积分公式 0 引育 力粘滞系数． 迄今为止，还没有一般的解析方法可以用来积 Stokes方程组的边界条件通常有以下两类： 分Navier-Stokes方程【1矗】，而只是在某些特殊情形1)第一类边界条件为约束边界条件，规定了边 下，得出了对所有粘性系数都成立的解．求解 界上的速度，即口=uo，它表示流体在边界上的速度 Navier-Stokes方程的精确解的问题，一般都会遇到必须与边界本身的运动速度相等． 难以克服的数学困难．这主要是因为方程是非线性 2)第二类边界条件为自然边界条件，流体在自 的，而不能应用叠加原理的缘故．因而许多学者转向 2 由边界上必须满足岛=∑O071，=o，或更一般些，t 对线性Stokes方程的研究．祝家麟用单层位势表示 =g．其中，露=(以。．n2)为界面外法向单位向量，它 定常Stokes方程的Dirichlet问题的解[3】．在二维情 况下，对于几何形状特殊的问题，余德浩采用与区域 指向所考察流体的外部．譬一(g。，＆)为流体在边界 上所受到的力． 形状相关的Green函数形成的自然边界元法求 解H～‘]．张耀明等以此为基础，归化出并建立了相应 2 Stokes问怒的Poisson积分公式 的数值求解体系[7]．半解析法半数值方法也是当前 设有如下Stokes方程组边值问题 n内 研究Stokes流的主流，其中四种为近年发展起来的 r一可△雎+gradp=0， 新方法【|qoj(多极子配点法，边界积分方程法，体内 ldivu=0， n内 (1) 奇点分布法和多极子矩法)．本文考虑了在定常运动 【H—l|o 11上 情况下，Reynolds数较小时不可压缩粘性流体的运单位圆外区域n上的Stokes方程组(1)的任意一 动，即略去Navier-Stokes方程的非线性项及局部惯组解必可表示成如下复变函数的实部或虚部形式嘲 性力项，讨论了平面Stokes问题．根据圆外区域 ． Stokes问题的自然边界元法，利用相应的边界积分 (2) 公