有限元法解圆柱绕流报告.docx

有限元法解圆柱绕流计算流体力学期末大作业工学院熊思 00986080 2012/6/3 【摘要】有限元法是求解计算流体力学问题的重要方法。用有限元法求解具体问题时, 首先需要将求解区域进行离散化, 即将求解区域划分为许多几何形状简单规则的单元, 在二维一般是三角形或四边形, 在三维是四面体或六面体。然后, 在每个单元内, 用一个比较简单的解析函数来逼近微分方程的解, 此函数在单元内用一组选定的单元基函数的线性组合表示, 而其中的系数通常是节点参数,它是待定的。这样,每个单元只要有适当数量的节点参数值,就可以满足对插值函数的光滑性和精度的要求。第三, 在满足微分方程和相应的初边值条件下, 对全部子域进行积分。对每个单元分别进行积分, 形成“单元方程”, 通过总体合成, 得到总体有限元方程组。最后, 用适当的方法解方程组, 可得节点参数值, 进而可求得各单元内的近似解。本文将以圆柱绕流问题为例,展示有限元法求解的一般步骤。一. 有限元法概述有限元法( finite element method )是一种高效能、常用的计算方法。有限元法在早期是以变分原理为基础发展起来的, 所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中, 这类场与泛函的极值问题有着紧密的联系。有限单元法最早可上溯到 20 世纪 40 年代。 Courant 第一次应用定义在三角区域上的分片连续函数和最小位能原理来求解 扭转问题。现代有限单元法的第一个成功的尝试是在 1956 年, Turner 、 Clough 等人在分析飞机结构时, 将钢架位移法推广应用于弹性力学平面问题,给出了用三角形单元求得平面应力问题的正确答案。 1960 年, Clough 进一步处理了平面弹性问题,并第一次提出了" 有限单元法" ,使人们认识到它的功效。自从 1969 年以来,某些学者在流体力学中应用加权余量法中的迦辽金法(Galerkin) 或最小二乘法等同样获得了有限元方程。具体而言, 参考差分法中网格化的做法, 把求解区域划分为有限多子区域, 称这些子区域为单元, 在每个单元上构造解的近似分布,将 Ritz 法或加权余量法应用到分块的逼近函数上。因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系。实质上,有限元法就是 Ritz 法或加权余量法。二. 问题描述考虑位于两块无限长平板间的圆柱体的平面绕流问题,几何尺寸如下图所示,来流为。由于流场具有上下左右的对称性,只考虑左上角四分之一的计算区域 abcde , 把它作为有限元的求解区域Ω。要求求解出整个区域中的流函数、以及压强值。 1. 边界 ab 为流线,取ψ=0, 2. 边界 bc 也为流线,同样取ψ=0, 3. 边界 cd ,切向速度=0 ,取 4. 边界 de 为流线, 满足于是在 ed 上, ψ=2, 5. 进口边界 ae 上, ψ= (本文中采取此条件) 也可以提自然边界条件我们以流函数ψ作为未知函数来解此问题,流函数所满足的微分方程如下: (1) 此处就是就是 cd 段边界,且切向速度=0,Γ 1和Γ 2 合起来是整个边界,并且此二者不重合。下面,按有限元方法的一般步骤来计算此问题。 1 .建立有限元积分表达式根据求解问题的基本控制方程, 应用变分法或加权余量法将求解的微分方程定解问题化为等价的积分表达式,作为有限元法求解问题的出发方程式。对于方程(1) ,它是一椭圆型方程,具有正定性,可以用变分法,这里直接给出泛函令其变分δ J=0 ,可以得到自然边界条件已经包含在变分表达式中( 其名称的由来) ,而本质边界条件必须强制ψ满足( 因此称其为本质边界条件,也称为强制边界条件)。如果根据原微分方程中无法给出泛函 J ,则可以用 Galerkin 加权余量方法得到积分方程,这相当于将原来的微分方程写为如下变分形式: 这里的δψ是函数ψ的改变量, 是一种“虚位移”, 在本质边界条件。因此, 上式做分部积分后,边界积分仅剩下。具体为即( 3 )式。可见, 如果ψ满足原来的微分方程和边界条件, 那么, 必然有ψ满足(4) 式, 进而满足(5) 式。注意,在(5) 式中,包含的边界Γ 2 上的边界条件信息,对边界Γ 1 的部分,仅知道它是给定了函数值的边界, 却不知道边界上的值是多少, 为了确定这些值, 还需要额外的处理方法。正是因为Γ 2 上的边界信息可以包含在积分表达式中, 这种边界条件也称为自然边界条件。 2 .区域剖分根据物理问题的特点以及区域的形状, 把计算区域分成许多几何形状规则但大小可以不同的单元, 确定单元节点的数目和位置, 建立表示网格的数据结构。采用的单元形状和