低雷诺数牛顿流体的边界积分法.docxVIP

低雷诺数牛顿流体的边界积分法 stokes问题广泛应用于采矿、沉淀、化工、土木工程等领域。今年是快速发展的流量分支。然而，到目前为止，很少有人能够准确地分析地壳运动。这两种方法主要用于研究相互依存的颗粒之间的弱流量动态干扰。为了考虑到颗粒之间的薄流量动态干扰，文献成功解决了一系列复杂的强干扰问题，即多极子配点法。然而，如果水流场中的一些边界状况是不规则的，无法找到合适的坐标系，并且是坐标的表面，则不适合使用多极子配点法进行计算。因此，有必要发展一种新的方法来灵活解决各种不规则形状下的stokes流问题。在文献中，边界积分方程法被广泛应用于各种不规则情况下的stokes流问题。该方法在70年代中期用于各种三维问题，主要因为二维stoks问题的解是奇怪的。后来，在文献的基础上，我们研究了一种更普遍的复合函数法，并使用了基于区域形状的自然边界规律。在任何边界条件下，都可以使用边界问题的边界积分公式。在文献的基础上，这项方法提出了满足圆外区域stokes方程解的更一般复变方法，并给出了具有一定距离间隔的stokes方程的速度解的边界积分公式。一旦边界问题是无限的，则只能计算边界速度为零的stokes运动。根据文献，本文提出了满足圆外区域stokes方程解的更一般复变函数，并给出了任意速度边界条件下stokes方程的速度解的边界积分公式。 1 应力分量范围的复变函数分析 在连续介质假设下,当流体的迁移惯性力远远小于黏性力时,不可压缩流体的运动在Ω内服从Stokes方程 {-ηΔ→u+gradp=0,div→u=0,(1) 式中:→u=(ux,uy)为速度矢量;p为压力;η为动力黏滞系数;区域Ω?R2. 平面区域Ω上的Stokes方程组的任意一组解可表示成复变函数的实部或虚部形式 {ux(x,y)=Re[-φ′(z)ˉz+φ(z)-ψ(z)],uy(x,y)=Ιm[φ′(z)ˉz+φ(z)+ψ(z)].(2) 应力分量也可表示为 {σx(x,y)=2ηRe[2φ′(z)-φ″(z)ˉz-ψ′(z)],σy(x,y)=2ηRe[2φ′(z)+φ″(z)ˉz+ψ′(z)],τxy(x,y)=2ηΙm[φ″(z)ˉz+ψ′(z)],(3) 式中φ(z),ψ(z)为Ω上的两个解析函数,z=x+iy,ˉz=x-iy. 圆外区域作为一多连体区域,其解析函数φ(z)和ψ(z)可能表现为多值的.在一般情况下,多连体可能具有m个内边界,本文考虑仅有一个内边界和一个外边界的情形.根据应力分量及速度分量的单值性条件分析可知,对于内边界外的任一点z,复变函数φ(z)及ψ(z)可表示为 {φ(z)=(α+iβ)lnz+φ**(z),ψ(z)=(-α+iβ)lnz+ψ**(z),(4) 式中φ**(z) ,ψ**(z) 为在内边界之外解析的复变函数,按罗朗级数展开, φ**(z)=∞∑n=0(anzn+a-nz-n),ψ**(z)=∞∑n=0(bnzn+b-nz-n),(5) 式中an,bn为复系数,n∈Z. 根据应力分量式(3)可知 σx+σy=8ηReφ′(z),(6)σy-σx+2τxy=4η[φ″(z)ˉz+ψ′(z)].(7) 将式(4)及式(5)中的第一式代入式(6),有 σx+σy=4η[α+iβz+α+iβzˉˉ+∑n=-∞∞n(anzn-1+anzn-1ˉ)].(8) 对于无界区域,考察式(8)的右边,随着|z|的增大,z的一次方及其以上的项会迅速增大.为了使应力为有限值,有下列条件 an=0,bn=0(n≥2).(9) 于是,复变函数φ(z)及ψ(z)可以表示为 {φ(z)=(α+iβ)lnz+(α′+iβ′)z+φ0(z),ψ(z)=(-α+iβ)lnz+(α″+iβ″)z+ψ0(z),(10) 式中:α′,α″,β′,β″为任意实常数;φ0(z),ψ0(z) 为区域内解析函数,且 {φ0(z)=∑n=0∞a-nz-n,ψ0(z)=∑n=0∞b-nz-n,(11) 式中a-n,b-n为复系数,n=0,1,2,……. 由式(2)可知 ux+iuy=φ(z)-zφ′(z)ˉ-ψ(z)ˉ=(α+iβ)lnz+(α′+iβ′)z+φ0(z)-z[α+iβz+(α′+iβ′)+φ′0(z)ˉ]-[(-α+iβ)lnz+(α″+iβ″)z+ψ0(z)ˉ].(12) 为了使速度保持为有限大,要求式(12)中α,β,β′=0以及α″+iβ″=0.因此,满足Stokes方程组解的复变函数表达式(2)的圆外无界区域解析函数φ(z)及ψ(z)可表示为 {φ(z)=α′z+φ0(z),ψ(z)=ψ0(z).(13) 要求速度ux和uy满足在无穷远处为零的条件来求解Stokes方程组,内边界速度之和就必须为零.因此,为了求解边界合速度不为零时扰动场在有限远处的Stok