圆外区域sopes问题的边界积分法.docxVIP

圆外区域sopes问题的边界积分法 0 基于半解析法的stkies方程 到目前为止，一般分析方法无法用于积分navier变量，只是在某些特殊情况下获得所有粘性系数的解。解决navier-stokes方程的精确解通常是困难的。这主要是因为方程是非线性的，不能应用叠加原理。因此，许多科学家正在研究线性stokes方程。朱立军用单调势表示定常stokes方程的线性问题的解。在二维情况下，在几何形状特殊的问题上，余德浩采用了基于区域形状的green函数形成的自然边界元法。基于此，张耀明等人对方程进行了归一分析，并建立了相应的数值解算方法。半分析法和半值法也是当前stokes流研究的主流。其中四种是近年来发展起来的新方法（多极点法、边界积分法、体内点分布法和多极子矩法）。在这项工作中，我们考虑到了在定常运动中，为了不压缩低粘度流的运动，即省略了davier-stokes方程的非线性项和局部惯性力项。本文根据圆外区域stokes问题的自然边界元法，使用相应的边界积分公式和边界速度值，对stikes方程的边值问题的解提供了解。通过计算数值方法进行了解算方法，并显示了该方法的准确性和有效性。 1 es方程组的边界条件 Stokes方程组为 {-ηΔu+gradp=0divu=0 其中,Δ为拉普拉斯算子,u为流速,p为压力,η为动力粘滞系数. Stokes方程组的边界条件通常有以下两类: 1) 第一类边界条件为约束边界条件,规定了边界上的速度,即u=u0,它表示流体在边界上的速度必须与边界本身的运动速度相等. 2) 第二类边界条件为自然边界条件,流体在自由边界上必须满足ti=2∑i=1σijnj=0,或更一般些,t=g.其中,n=(n1,n2)为界面外法向单位向量,它指向所考察流体的外部.g=(g1,g2)为流体在边界上所受到的力. 2 复系数的测定 设有如下Stokes方程组边值问题 {-ηΔu+gradp=0?Ω内divu=0?Ω内u=u0Γ上 (1) 单位圆外区域Ω上的Stokes方程组(1)的任意一组解必可表示成如下复变函数的实部或虚部形式 {u1(x,y)=Re[-φ′(z)ˉz+φ(z)-ψ(z)]u2(x,y)=Ιm[φ′(z)ˉz+φ(z)+ψ(z)]p(x,y)=-4ηReφ′(z) (2) 其中φ(z),ψ(z)为Ω上的两个解析函数,z=x+iy,ˉz=x-iy.当单位圆内边界的合速度为零时,设 {φ(z)=∞∑0α-nz-nψ(z)=∞∑0β-nz-n (3) 其中,α-n=ˉα-n?β-n=ˉβ-n?n=1?2??,且α-n及β-n为复系数。 将(3)式代入(2)式即可得区域内以及单位圆边界上速度与压力的表达式.此外,设单位圆边界上速度为 {u1(1,θ)=∞∑-∞aneinθ?a-n=ˉa-nu2(1,θ)=∞∑-∞bneinθ?b-n=ˉb-n (4) 其中a0,b0为实数,ai,bi为复数,i≠0,将(4)式与所得单位圆边界上速度与压力的表达式进行比较,得出未知系数并代入区域内速度及压力的表达式,可得 u1(r,θ)=∞∑-∞anr-|n|einθ+(1-1r2)[12cos2θ?∞∑-∞(|n|an-inbn)r-|n|einθ+12sin2θ?∞∑-∞(inan+|n|bn)r-|n|einθ]u2(r,θ)=∞∑-∞bnr-|n|einθ+(1-1r2)[12sin2θ?∞∑-∞(|n|an-inbn)r-|n|einθ-12cos2θ?∞∑-∞(inan+|n|bn)r-|n|einθ] p(r,θ)=2ηr[cosθ∞∑-∞∞∑-∞(|n|an-inbn)r-|n|einθ+sinθ∞∑-∞(inan+|n|bn)r-|n|einθ] 引入公式 12π∞∑-∞r-|n|einθ=r2-12π(1+r2-2rcosθ)?Ρ(r,θ),r＞1(5) 利用公式(5)及坐标变换公式,容易求得极坐标分解下的圆外区域Stokes问题的边界积分公式 ur(r,θ)={cosθΡ(r,θ)+r2-R22r2[cosθ(-r???rΡ(r,θ))+sinθ??θΡ(r,θ)]}*ur(R,θ)+{sinθΡ(r,θ)+r2-R22r2[sinθ(-r??rΡ(r,θ))-cosθ??θΡ(r,θ)]}*uθ(R,θ)uθ(r,θ)={-sinθΡ(r,θ)+r2-R22r2[sinθ(-r???rΡ(r,θ))-cosθ??θΡ(r,θ)]}*ur(R,θ)+{cosθΡ(r,θ)-r2-R22r2[cosθ(-r??rΡ(r,θ))+sinθ??θΡ(r,θ)]}*uθ(R,θ)p(r,θ)=2ηr{[cosθ(-r??rΡ(r,θ))+sinθ??θΡ(r,θ)]* uθ(R,θ)} 其中*表示卷积. 3 边界条件及