化工管路设计与模拟.pptVIP

化工计算机辅助设计 第四讲 管路设计与模拟 管路设计与模拟 流体流动方程的建立 阻力和管路计算 流体流动方程组的数值求解 管路设计 管路模拟 一. 流体流动方程的建立 流体的流动形式可分为稳定流动与不稳定流动。 稳定流动是指流体流动时的速度、压力等参量只是流体在流动体系中位置的函数。 不稳定流动是指流体流动时的参量值与时间有关。 对于流体在密闭管道内的稳定流动，管道内任两截面处流体的质量流量G相同。只考虑流动中流体的机械能部分，易得伯努利方程的微分形式为： 1. 对于不可压缩流体，比容V为常数，此时将微分式从1截面积分至2截面，有： 2. 对于可压缩流体，比容V可变，但质量流速G仍不变。阻力系数λ只与质量流速相关，此时将微分式两端同时除以V2： 将积分式中各项分别记为： 可压缩流体的多变过程可采用下面方程式同一描述： 将多变过程方程代入压强变化项方程： 对于等温过程，k＝1，则有： 二. 阻力和管路计算 流体在管路中流动时的阻力分为直管阻力和局部阻力。 直管阻力是指流体与直管管壁摩擦所产生的压力降。 局部阻力是指流体流过管路中的管件、阀门等器件时产生的压力降。 流体流动受到的总阻力是上面两种情况之和。 1. 直管阻力计算 圆直管的阻力一般表示为正比于流体动能和相对管长，其比例系数λ便是管路摩擦阻力系数。 λ是雷诺数Re和管壁相对粗糙度ε/D的函数，即： 管路摩擦阻力系数方程形式一般按雷诺数及管壁粗糙度不同分为不同形式，以分段函数形式表示。 Churchill为了方便工程运用，将管路摩擦阻力系数用一个连续函数形式表示。此方程十分有利于不同流速、粗糙度情况下的直管阻力系数计算。 Churchill方程如下： 2. 局部阻力计算 局部阻力的计算常采用两种方法，一种是当量长度法，另一种是阻力系数法。 当量长度法：将由于管件局部阻力引起的压降转化为引起相等压降的同一管径直管的长度值。 阻力系数法：将局部压降表示为管路中流体动能的k倍。 一般当量长度法应用更广泛。 三. 流体流动方程组的数值求解 He= ?PZ + Σhf Σhf= λ(L/D)u2/2 λ=f(Re,ε/D) Re=ρD u/μ u=V/S S=0.785D 2 He =f(V) 变量分析 物性参数: 密度ρ ，粘度μ; 设备参数: 管径D，直管和管件当量长度之和L， 管壁粗糙度ε，流动截面积S; 操作参数: 每小时流体体积流量V，流速u，总势 差?PZ，摩擦阻力损失Σhf，摩擦系数 λ ， 雷诺准数Re，泵对每kg流体提供 的能量He; 1. 简单管径设计 设定:外加机械功He=0，造成流体流动的推动力是总势能差(?PZ0)，即上游的势能大于下游的势能。 管径设计 分析 He=0， ?PZ=－Σhf (?PZ0) 流体种类确定，管路布置确定，管子ε确定 设计D=?使得Vs=要求 优先排序法原理 条件 1．对于某个未知数是可解的，而和其它方程无关，则应首先解这样的方程； 2．含有一个不出现在其它任何方程中的未知数，则应在解出所有其它方程后再来解这个方程求此未知数。 原理 1 首先解只含一个未知数的方程； 2．如果一个未知数只出现在一个方程中，则应在确定这个方程的所有其它未知数后，再来解这个方程求此未知数。这样的方程放在求解决序的最后； 3．将剩下的只含两个未知数的方程（如果有的话）放在一起。选择其中一个方程，经过重新整理后，能将一个未知数（称为解出变量）用另一个未知数（剥离变量）表达出来。 剥离变量应为容易取初值的变量 优先排序法 f1(He, ?PZ, Σhf)=0 f2(Σhf, λ ,L,D,u)=0 f3(λ, ε ,Re,D)=0 f4(Re,D,u,ρ, μ )=0 , f5 (u,V,S)=0 , f6(S,D)=0 优先排序法框图 联立求解 x1=u x2=D x3=Re x4=λ 2. 管路操作模拟 设定:外加机械功He=0，造成流体流动的推动力是总势能差(?PZ0)，即上游的势能大于下游的势能。 f1(He, ?PZ, Σhf)=0 f2(Σhf, λ ,L,D,u)=0 f3(λ, ε ,Re,D)=0 f4(Re,D,u,ρ, μ )=0 , f5 (u,V,S)=0 , f6(S,D)=0 优先排序法框图 模拟泵与管路工作点（介绍） He=A+B ? V2 管路网络模拟（简介） 管路网络模拟 已知 水平管路网络有n个联结节点 第i点和j个节点的管路直径Dij，长度Lij