第八章-瞬变流1分析课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,输油管道设计与管理,输油管道设计与管理,1,第八章 输油管道瞬变流,第一节 管道的瞬变流动过程,第二节 管道瞬变流动过程的数学描述与特征线解法,第三节 输油管道中的水力瞬变流动,第四节 输油管道瞬变流动过程的控制,第八章 输油管道瞬变流第一节 管道的瞬变流动过程,2,稳定流动：,压力管道中任一点的流速和压力仅与该点的位置有关，而与时间无关的流动称为稳定流动。,不稳定流动,瞬变流动（Transient flow）：,流体从一种稳定流动状态过渡到另一种新的稳定流动状态的过程。,水击（Water Hammer）：,又称“水锤”。有压管道中因流速、压力急剧变化而引起压力波在水中沿管道传播的现象。,第一节 管道的瞬变流动过程,稳定流动：压力管道中任一点的流速和压力仅与该点的位置有关，而,3,水击现象（Water-hammer Phenomena）：,在有压管道系统中，由于某一管路中的工作状态的突然改变，引起管内液体流速的急剧变化，同时液体压强大幅度波动，这种现象称为水击现象。,正水击（Positive Water-hammer）：,当管道阀门迅速关闭时，管中流速迅速减小，压强显著增大，这种水击称为正水击。,负水击（Suction Water-hammer）：,当管道阀门迅速开启时，管中流速迅速增大，压强显著减小，这种水击称为负水击。,第一节 管道的瞬变流动过程,水击现象（Water-hammer Phenomena）：,4,一、水击的研究历史,1897年，茹可夫斯基（俄）发表了他的经典报告，全面而正确的阐述了水击的机理，给出了惯性水击压力和水击波传播速度的计算公式。,1902年，阿列维（意）建立了不稳定流动的微分方程，奠定了水击分析的理论基础。,20世纪60年代以后，借助于电子计算机的发展，数值法求解阿列维方程成为可能。斯特里特（美）1963年创建的特征线法，开创了水击分析的新局面。,一、水击的研究历史1897年，茹可夫斯基（俄）发表了他的经典,5,二、水击产生的原因流量突变,计划地调整输量或改变输送流程,开关阀门：阀门最后关闭的0.10.3秒，水击最强。,油品的注入与分输,泵站节流、调速、启停泵,改变输送压力、流量,油品切换,混油通过中间泵站,二、水击产生的原因流量突变计划地调整输量或改变输送流程,6,二、水击产生的原因流量突变,事故引起的流量变化,停电,阀门失灵,管道泄漏,二、水击产生的原因流量突变事故引起的流量变化,7,三、水击形成过程,水击过程四个阶段：,增压逆波,减压顺波,减压逆波,增压顺波,三、水击形成过程 水击过程四个阶段：,8,四、管道瞬变流的基本公式,直接瞬变压力公式,根据儒可夫斯基的水击理论，由流速瞬间变化直接产生的压力脉动值，可用下式计算,如果阀门阶段性关闭，则,四、管道瞬变流的基本公式直接瞬变压力公式,9,四、管道瞬变流的基本公式,压力波传播速度,管内压力波的传播速度取决于液体的可压缩性和管子的弹性。,含气液体中的波速计算（略）,四、管道瞬变流的基本公式压力波传播速度,10,五、水击危害,1、正水击超压，管道破裂,2、负水击管道失稳、液柱分离、气体逸出,3、振动,五、水击危害,11,六、水击分析的意义,在优化设计中的应用,在编制运行方案中的应用,预测混油界面及检漏,在操作培训中的应用,六、水击分析的意义在优化设计中的应用,12,七、水击的控制,对于管道瞬变流动过程，控制的目的是：,避免管道超压（包括超高压和超低压）；,减轻管道运行参数的脉动，维持管道的平稳运行。用于控制管内瞬变流动过程的装置和措施很多，根据其作用原理可分为两类：,改变流速变化过程（幅度和时间）,使用各种压力保护措施,七、水击的控制对于管道瞬变流动过程，控制的目的是：,13,压力保护装置和措施,调节阀,压力保护,泄压阀,回流保护,泵机组自动停运,超前保护,七、水击的控制,压力保护装置和措施七、水击的控制,14,第二节 管道瞬变流的数学描述与特征线法,分析管道的瞬变流动过程，使用弹性理论比较符合客观实际。描述管道瞬变流动过程的弹性理论以液体可以压缩和管壁产生变形为基础，分析受压力波作用的微元体的受力与质量守恒，建立,运动方程,和,连续方程,的微分形式。将此两方程联立，对应相应的边界条件，是用数值方法即可分析管道的瞬变流动过程。,第二节 管道瞬变流的数学描述与特征线法分析管道的瞬变流动过,15,一、管道瞬变流动过程的数学描述,运动方程,连续方程,一、管道瞬变流动过程的数学描述运动方程,16,（一）,运动方程,运动方程是从处于瞬变流动的管内流体中选取隔离体，应用牛顿第二定律建立的。推导方程是有两点假设：,（1）管内流体为均质一维流动，这就意味着管道横截面上的压力、流速和密度是均匀分布的；,（2）管子和流体变形均在线性弹性变形范围内。,（一）运动方程运动方程是从处于瞬变流动的管内流体中选取隔离体,17,第八章-瞬变流1分析课件,18,19,20,把压力P换算成液柱高度：,由于：,隔离体内视为均匀密度。,由于隔离体足够小，,（8-9）,把压力P换算成液柱高度：由于：隔离体内视为均匀密度。由于隔离,21,（二）连续方程,连续方程描述了时刻t流入及流出隔离体的流体净质量等于隔离体内的质量随时间的变化率。,对隔离体，考虑时间的质量守恒，应存在,展开，忽略高阶小量，整理后可得,（二）连续方程连续方程描述了时刻t流入及流出隔离体的流体净质,22,由于,（8-10）,由于（8-10）,23,二、,特征线解法,分析液体管道不稳定流的两个基本方程,一阶拟线性双曲型偏微分方程组,由于二次摩擦项的存在，目前用分析方法还无法得到它们的解析解，只能用数值方法求解。,求解液体管内不稳定流问题的最常用的一种数值方法特征线解法,二、特征线解法,24,（一）特征方程,特征方程是用一对特征值，将两个偏微分方程线性组合，即,（一）特征方程 特征方程是用一对特征值，,25,（一）特征方程,（一）特征方程,26,（二）有限差分方程,特征方程虽然是常微分方程的形式，但由于摩阻项是非线性的，仍不能用微分方法得到解析解。下面介绍用有限差分方法求数值解。,使用数值方法求解特征方程，必须把管道沿线长度方向离散成若干管段，把瞬变过程离散成若干个时间步长，在小范围内进行近似计算。,（二）有限差分方程特征方程虽然是常微分方程的形式，但由于摩阻,27,沿管道（简单管道）长度等间距分成N段，每段长 ，称为空间步长，管道的结点数为N+1；压力波传播距离所需时间 ，称为时间步长。这样，把x-t平面画成网格。,沿管道（简单管道）长度等间距分成N段，每段长,28,第八章-瞬变流1分析课件,29,在一个网格范围内，沿特征线，在i-1和i之间，对方程分离变量后积分,在一个网格范围内，沿特征线，在i-1和i之间，对方程分离变量,30,由于无法知道瞬变过程中V与x之间的函数关系，后面一项积分只能作数值计算。对于摩阻较小的管道，后面一项积分可采用一阶近似方法计算。差分方程为,由于无法知道瞬变过程中V与x之间的函数关系，后面一项积分只能,31,在进行管道瞬变流动计算时，往往用流量代替流速更方便。重新整理成下面的形式：,（8-18）,（8-19）,在进行管道瞬变流动计算时，往往用流量代替流速更方便。重新整理,32,带有p下标的表示当前时刻的计算值，不带p下标的表示已知参数。显然，是分别由i-1和i+1点已知的参数计算的常系数。联立求解（8-18）和（8-19）式可得,对于i结点，时间终了时的参数和可由 时间开始时i-1和i+1结点的参数计算求得。,（8-20）,（8-21）,带有p下标的表示当前时刻的计算值，不带p下标的表示已知参数。,33,【例】站间管道如图所示，管材 ，长25.5km，稳态下流量Q=0.025m,3,/s，,（V1.47m/s），m=0.125，f=14.13，a=1214m/s，稳态时末端进站压头为434米。末端阀门在0s时瞬间关闭，上游罐液位恒定。,【例】站间管道如图所示，管材,34,0s,1s,2s,3s,4s,Q=0.025,502,485,468,451,434,i=0.014,18,19,20,21,22,0s1s2s3s4sQ=0.025502485468451,35,（三）扰动的影响区,沿管道不同的位置，受到扰动的影响的时间是不同的。,发生扰动的位置不同，管道沿线受到扰动影响的过程也不同。,以P点为起点，两特征线上方包络的范围内的状态受到影响。因此，该区间也被称为P点扰动的影响区。而两特征线下方区间管内没有受到P点扰动的影响，只取绝于初始状态，所以该区间也称为初始条件支配区。,（三）扰动的影响区沿管道不同的位置，受到扰动的影响的时间是不,36,（四）初始条件和边界条件,初始条件,稳定工况管内的压力分布和流量分配，作为瞬变过程计算的初值,边界条件,内部边界条件,外部边界条件,（四）初始条件和边界条件初始条件,37,三、内部边界条件,管道变径点,内部局部损失,中间恒速泵站,减速离心泵,中间泵站失电,三、内部边界条件,38,四、外部边界,上游边界,起点恒液位罐,首站恒速泵,下游边界,四、外部边界,39,五、,基本的计算程序,五、基本的计算程序,40,