15-2014中文第15章有限元方法在流体力学中的应用.pptVIP

Chapter 15 有限元方法在流体力学中的应用 Applications of FEM in Fluid Mechanics (Fundamentals of finite element analysis by David V. Hutton) 1 引言 不可压缩流体(incompressible flow): 密度不变 可压缩流体(compressible flow) Newton’s law of viscosity μ为绝对粘度(absolute viscosity). 是流体的基本材料参数,与其抗剪应力性能直接相关。 如果流体的粘度很小，则可忽略其剪应力，理想化为无粘性流体。 2 不可压缩流体的控制方程 质量守恒 连续性方程 u, v, w, 是速度在x,y,z方向的分量 质量守恒要求一个体积里面的质量变化率等于该体积的质量净流入速度。 体积dV里的总质量为ρdV, 注意到 dV 为常数, 所以有： 从x,y,z三个方向流入质量造成的控制体积中的质量变化率分别表示为： 所以质量的变化率为： 注意到 dV = dx dy dz, 于是得到连续性方程为： 对于稳态的不可压缩流体(steady flow of an incompressible fluid),密度与时间和空间坐标无关，于是 有旋流和无旋流(Rotational and Irrotational Flow) 把流体流动分为： 有旋流动(rotational) ---平动和转动混合 无旋流动(irrotational) ---仅有平动。流体微元体积没有纯转动。 如果下列三式不能同时满足，则为有旋流动: 我们目前仅考虑无旋流动。 二维流的流函数(steam function) 对于2D,稳态,不可压缩,无旋流动 连续性方程为 无旋条件退化为 如果引入 ψ(x , y) (流函数)，则连续性条件自动满足: 无旋条件变为 (Laplace’s equation) 流函数的物理意义 流线(streamlines) : x-y平面内的曲线，其上的流函数为常数。 流线上的任一点的切向量可以表示为 nt = dx i + dyj， 该点处的流体速度为 V = ui + vj. 则 V × nt = (?v dx + u dy)k=0 因为：两个非零向量的叉积为零意味着这两个向量平行，所以： 在流线上任一点，流体速度正切于流线。 有限元列式 具有M个节点的有限单元内的流函数可以表示为: 利用Galerkin方法, 单元的残差方程为： or 应用Green-Gauss定理，上式变为 其中 S 为单元边界 (nx , ny ) 为边界单位外法线向量。代入流函数表达式，有： 写为矩阵形式 二维流动的流速势函数(Velocity Potential Function) 假定存在流速势函数φ(x , y) 使得 则无旋条件能自动满足: 连续性条件变为: (again, we obtain Laplace’s equation) 在流速势函数为常数的曲线上，有 所以可见，速度向量垂直于流速势为常数的曲线。 于是，流线和流速势等值线形成了正交网格，称为流网 (flow net). 有限元列式与流函数的情形类似： 单元刚度阵与流函数法的相同，但是节点力完全不同。 不可压缩粘性流 假定: 可以看作二维问题 不涉及热 3. 密度和粘度为常数 4. 稳态(不随时间变化) 表示动量守恒的 Navier-Stokes方程为: u,v = x, y方向速度分量 ρ= 密度 p = 压力 μ= 绝对粘度 FBx , FBy = x, y方向单位体积的体力 斯托克斯流(Stokes Flow) Stokes flow (or creeping flow)：流体流动的速度很小，惯性项比粘性项小很多，可以忽略。 高粘度流体，如融化的高分子材料。 动量方程变为 上式和连续性方程一起构成了三个未知数u(x, y), v(x, y), and p(x, y)的三个方程。 这三个节点变量可以表示为： 应用Galerkin方法，残差方程为 最后得到一个 3M 个线性方程系统，求解个3M 未知的节点量 {u}, {v}, {p} : *