有限元法解圆柱绕流问题资料.docxVIP

有限元法求解无限流体中的圆柱绕流问题2016年01月12号一．问题描述 考虑位于两块无限长平板间的圆柱体的平面绕流问题，几何尺寸如下图所示，来流为vx=1,vy=0。由于流场具有上下左右的对称性，只考虑左上角四分之一的计算区域abcde，把它作为有限元的求解区域Ω。要求求解出整个区域中的流函数、vx、vy以及压强值。图1：圆柱绕流模型二．数学建模 在足够远前方选与来流垂直的控制面ae,cd是沿y轴，亦即一流动对称轴，bc是物面，ab亦是流动对称轴，所要考虑的流动区域即由线abcdea所围成的区域，在这一区域中有：1.边界ab为流线，取ψ=0，∂φ∂n=0;2.边界bc也为流线，同样取ψ=0，∂φ∂n=0;3.边界cd，切向速度vτ=0，∂ψ∂n=0，取φ=0；4.边界de为流线,满足ψe=ψa+aevxdy=02dy=2于是在ed上，ψ=2，∂φ∂n=0;5.进口边界ae上，ψ=ψa+ayvxdy=y（本文中采取此条件）也可以提自然边界条件∂ψ∂n=0，∂φ∂n=vx=1我们以流函数ψ作为未知函数来解此问题，流函数所满足的微分方程如下：∇2ψ=0ψ|Γ1=ψ（本质边界条件）∂ψ∂n|Γ2=-vs(自然边界条件)（1）此处Γ1指ab,bc,de和ae四段边界，而Γ2就是就是cd段边界，且切向速度vs=0，Γ1和Γ2合起来是整个边界，并且此二者不重合。下面，按有限元方法的一般步骤来计算