计算流体力学5.pptVIP

 系数可由三个边界节点定出，如 单元内解拉普拉斯方程 以中心节点代入 取h=k=1 则边中点与角点 若边中点0.25，角点0.0，则为5点格式。 而本例为9点格式，更为精确。 其中 2.矩形单元 线性插值（4节点单元），设 则 二次插值（9节点单元），略 曲线四边形单元，可采用双线性变换 3.三维情形 (1)四面体: 4节点单元 10节点单元 (1)正六面体: 8节点单元 20节点单元（略） 例（同前）： 步骤（3） 积分整理得 (4)解出 (5)近似解 4.矩量法 取权函数 则 例（同前）： 步骤（3）取  积分整理得 (4)解出 (5)近似解 5.伽辽金法（Galerkin） 取基函数本身作权函数，即取 则 例（同前）： 步骤（3）取 构造二级近似解 余量 代入伽辽金 加权积分式 积分 (4)解出 (5)近似解 此例题精确解为 由各种解比较结果可见，伽辽金法精度最高。 因此，该方法获得广泛应用。 3 切削叶轮外径以调节泵的性能 切削后的叶轮与切削前不再符合几何相似. 叶轮切削前后的速度比 叶轮切削前后性能参数的关系如下 （近似认为切削前后容积效率相等，叶片排挤系数相等， 水力效率相等，涡流系数相等）： 1.低比转数的泵（ns≤150） 叶片出口宽度b2小，近似认为叶轮切削前后 流量比 扬程比 轴功率比 上述三个关系式称为第一切削定律。 2.中、高比转数的泵（ns＞150） 近似认为切削前后叶轮出口面积不变 上述三个关系式称为第二切削定律。 切削叶轮调节工况分析 根据第一切削定律，可以推导得出其切削曲线是一条过原点的直线； 根据第二切削定律，可以推导得出其切削曲线是一条过原点的 二次抛物线。B点和C点是满足第一切削定律的工况点， B点和D点是满足第二切削定律的工况点。 切削叶轮的调节方法，其切削量不能太大， 否则效率明显下降。泵的最大切削量与比转数有关。 泵的比转数 60，允许最大切削率为20％ 泵的比转数120，允许最大切削率为15％ 泵的比转数200，允许最大切削率为11％ 泵的比转数300，允许最大切削率为9％ 泵的比转数350，允许最大切削率为7％ 泵的比转数350以上，允许最大切削率为0 所谓切削率，即 切削叶轮的调节方法，不增加额外的能量损失，设备效率下降少， 是一种节能的调节方法，但需要停机换轮，因此常用于 水泵的季节性调节。 * 第五章 边界元法、有限体积法和有限分析法 基本思想：将控制微分方程转化为边界积分方程，再用有限元法来处理边界积分方程。 特点： 1.区域内满足微分方程，边界上近似满足边界条件。 2.维数减少一个，可以简化计算。 3.精度一般高于有限元法。 4.系数矩阵不对称并为满阵，需要解析函数的基本解， 目前主要适用于线性问题。 一、边界元法（Boundary Element Method或BEM） 1.边界积分方程的建立 例：Laplace方程 其伽辽金方程为 （在Ω内） （在Γ1上） （在Γ2上） 当近似解不要求满足边界条件时，由格林公式可得: 引入权函数w=δu代入上式 ? 而 代入上式得 --------（a） 称为边界元出发方程。 若选权函数w满足 而不要求满足边界条件时，则 选择另一权函数w，使其对区域内任一点i满足 （在i点） （不在i点处） 则 由（c） 则（a）式变为 说明:内点的函数值可用边界上的函数值及其 法向导数值沿区域的边界积分来表示。 满足(c)式的解称为基本解。对于各向同性介质 r---区域内任意点i到边界积分点的距离。 将i点移到边界上，则 （d）式称为边界积分方程。 （d）式也可写为为 3.数值离散 (1)边界上剖分和插值 a.常数单元（1节点） 取单元中点为节点,则 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) j i j i j j i j i j j i j i j j j i j i j y y x x q n u q u u j j j j = = = ÷ ? ? ? è ? ? ? = = b.线性单元（2节点） 取单元两端点为节点,则 c.二次单元（3节点） 取中点及两端点为节点,则 （2）数值化 边界积分方程（e）的离散形式 若采用常数单元： 引入记号 则 ? 将已知值 或写成矩阵形式 等代入，整理成代数方程组形式 式中 解方程组后，则全部边界上u,q为已知， 再求出区域内任意点u值为 二、有限体积法（Finite Volume Method或FVM） 1.基本思想： 将计算区域划分为一系列不重复的控制体积，使每个网格点周围有一个控制体积，