单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,计算水力学,第一章 基本方程,李光炽,计算水力学,计算水力学 李光炽计算水力学,1,1.守恒定律,1.1,微分方程的意义,连续方程的矢量式,:,对于不可压缩流体，上式简化为,divU=0,物理意义可表述为：流体密度的增加等于位体积内流入的净质量。,李光炽,计算水力学,1.守恒定律1.1 微分方程的意义李光炽计算水力学,2,设 表示水流和输运现象中某种物理量在单位质量中的含量。当为,i,方向的速度分量,Ui,时，其物理意义为单位质量流体所含的,i,方向的动量。当,=1,时，表示单位质量流体所含的质量,李光炽,计算水力学,设 表示水流和输运现象中某种物理量在单位质量中的含量。当,3,设,J,表示对因变量发生影响的通量。J在,x,方向的分量即通量 ，表示进入面积为,dydz,表面的通量，而离开对立面的通量可表示为,，因而,x,方向的净通量为,李光炽,计算水力学,设J表示对因变量发生影响的通量。J在x方向的分量即通量,4,单位体积的净通量,(1),微分方程中的因变量可看作单位质量流体所含的物理量，(2),微分方程中的各项均表示某种物理机制，对单位体积中因变量的含量发生影响，(3),整个微分方程，是上述各项的总和，表示各种影响之间的平衡或守恒。,李光炽,计算水力学,单位体积的净通量 李光炽计算水力学,5,1.2,守恒定律,一、输运物质守恒和热量守恒,设水流的速度场为U，水流挟带着泥沙悬移质或某种污染物质，其质量浓度(即单位体积中输运物质质量与混合物质量之比)为C，则C的守恒可表示为,李光炽,计算水力学,1.2守恒定律 一、输运物质守恒和热量守恒李光炽计算水力学,6,物质守恒,第一项可理解为单位体积内所含输运物质随时间的变化率。2-1是该物质的对流通量，即由流速场引起，随水流运动的通量。2-2表示扩散通量，通常由,C,的梯度所引起对流通量和扩散通量的散度，即单位体积的净通量，构成微分方程的第二顷。方程右端的是输运物质在单位体积中的产生率，称为源项。,李光炽,计算水力学,物质守恒第一项可理解为单位体积内所含输运物质随时间的变化率。,7,物质守恒,费克,(Fick),扩散定律表示扩散通量,J,C,李光炽,计算水力学,为,扩散系数,物质守恒费克(Fick)扩散定律表示扩散通量 JC李光炽计算,8,张量形式为,李光炽,计算水力学,描述水流中热量守恒,张量形式为 李光炽计算水力学描述水流中热量守恒,9,二、动量守恒,(NavierStokes,方程,),李光炽,计算水力学,写出水流动量守恒的微分方程，即著名的纳维埃,-,斯托克斯方程,式中,P,为压力，,v,为运动粘性系数，,B,为单位体积的体积力,二、动量守恒(NavierStokes方程)李光炽计算水,10,三、紊流的时均,(Reynolds),方程,对纳维埃,斯托克斯方程作时均演算,并采用鲍辛捏斯克,(Boussinesq),关于紊动粘性系数的假设，便可得出时均流的动量方程：,李光炽,计算水力学,式中 为紊动粘性系数，字母上方的横线表示时均值,三、紊流的时均(Reynolds)方程 对纳维埃斯托克斯方,11,2.通用微分方程,如果用 表示通用变量，则有各方程的通用形式：,李光炽,计算水力学,称为水流和输运现象的通用微分方程。通用微分方程包含变化率项、对流项、扩散项和源项。令因变量代表不同的物理量，并对扩散系数 和源项,S,作相应的调整,2.通用微分方程 如果用 表示通用变量，则有各方程的通,12,连续方程、,浓度输运方程,李光炽,计算水力学,得出连续方程,得出浓度输运方程,连续方程、浓度输运方程 李光炽计算水力学得出连续方程 得出浓,13,温度输运方程、动量方程,李光炽,计算水力学,得出温度输运方程,得出动量方程,温度输运方程、动量方程 李光炽计算水力学得出温度输运方程 得,14,通用微分方程的张量形式为,李光炽,计算水力学,爱因斯坦求和约定（,Einstein summation convention,）,通用微分方程的张量形式为 李光炽计算水力学爱因斯坦求和约定（,15,3 圣维南方程,一、连续方程,质量守恒定律：,通过控制面进到控制体的质量控制体内质量增量,李光炽,计算水力学,3 圣维南方程 一、连续方程李光炽计算水力学,16,李光炽,计算水力学,1,1,2,2,t,t+,t,Q,X,质量守恒原理示意图,李光炽计算水力学1122tt+t QX质量守恒原理示意,17,由质量守恒定律可得：,李光炽,计算水力学,化简后得：,简化过程中用到了为常数的假设。式中：为过水断面面积，为断面流量,由质量守恒定律可得：李光炽计算水力学化简后得：,18,二、动量方程,动量方程的推导主要依据是动量守恒定律。由于动量是一矢量，推导是建立水流流动方向的动量方程。,通过控制面流进到控制体内的动量,作用于控制体外力的冲量,控制体积内动量的增量,李光炽,计算水力学,二、动量方程动量方程的推导主要依据是动量守恒定律。由于动量是,19,李光炽,计算水力学,1,1,2,2,t,t+,t,X,Qu,动量守恒原理示意图,李光炽计算水力学1122tt+t XQu动量守恒原理,20,经过时间后,控制体的动量发生了变化,其增量为:,李光炽,计算水力学,在时间内通过控制面流入到控制体的动量为：,经过时间后,控制体的动量发生了变化,其增量为:李光炽,21,作用于断面之间的控制体上的力有：重力、摩阻力及水压力,重力在水流方向上的分量为：,李光炽,计算水力学,在非恒定流情况下摩阻损失与恒定流情况下差别不大，仍可用曼宁公式、谢才公式或流量模数公式，,由此可得出作用于控制体上的摩阻力为,作用于断面之间的控制体上的力有：重力、摩阻力及水压力,22,压力：压力对断面之间水体的作用，分为两部分：断面压力差和,侧壁上的压力,李光炽,计算水力学,P,Fa,Fa,控制体上压力分布图,压力：压力对断面之间水体的作用，分为两部分：断面压力,23,李光炽,计算水力学,h,B,断面压力积分示意图,李光炽计算水力学hB断面压力积分示意图,24,断面上的压力为,李光炽,计算水力学,断面与断面的压力差为：,断面上的压力为 李光炽计算水力学断面与断面的压力差为：,25,李光炽,计算水力学,1,1,2,2,侧壁压力分力计算示意图,李光炽计算水力学1122侧壁压力分力计算示意图,26,作用于断面1-2之间的侧壁压力在水流方向上的分量等于断面与断面的差值部分面的压力，其大小为,李光炽,计算水力学,作用于控制体上的总压力在水流方向上的分量为,作用于断面1-2之间的侧壁压力在水流方向上的分量等于断面与,27,可得作用于控制体上的总外力为,李光炽,计算水力学,由动量守恒定律可得：,可得作用于控制体上的总外力为 李光炽计算水力学由动量守恒定律,28,化简可得动量方程,李光炽,计算水力学,可写成等价形式,化简可得动量方程 李光炽计算水力学可写成等价形式,29,明渠非恒定流的基本方程圣维南方程组,李光炽,计算水力学,明渠非恒定流的基本方程圣维南方程组 李光炽计算水力学,30,4 定解条件和定解问题,非恒定问题必须给出初始条件和边界条件。初始条件与边界条件是确定微分方程解的必不可少的条件,合称为定解条件。,定解条件少给了,就无法定解，即所谓欠定;多给了可能产生矛盾,即所谓过定;,只有给出适当个数的定解条件才能求解,即所谓适定。,适定,性：,解的存在性、唯一性和稳定性,李光炽,计算水力学,4 定解条件和定解问题 非恒定问题必须给出初始条件和边界条件,31,对流方程,李光炽,计算水力学,称为特征线,即的值沿特征线保持不变,对流方程 李光炽计算水力学称为特征线 即的值沿特征线保持不,32,李光炽,计算水力学,x,t,C,D,a,b,A,B,特征线示意图,李光炽计算水力学xtCDabAB特征线示意图,33,将特征线视为有向线段,其方向与增加时的方向一致。定解域边界(包括初始时刻)上的某点是否要给出定解条件,可从该点处作特征线,若其方向指向定解域内,则需给出定解条件,否则不必给。,在计算域的边界,(请注意，这里没有必要区别什么是初值条件边界，什么是边值条件边界),处，边界条件的个数就等于穿越边界而进入定解域的特征线的条数,。,李光炽,计算水力学,将特征线视为有向线段,其方向与增加时的方向一致。定解域边界,34,圣维南方程组,沿顺特征线,李光炽,计算水力学,沿逆特征线,圣维南方程组 沿顺特征线 李光炽计算水力学沿逆特征线,35,缓流,和急流,由于是缓流，因而水流的佛汝德数 ，故域内一点处的两条特征线分别指向上、下游。,由于是急流，因而水流的佛汝德数 ，故域内一点处的两条特征线均指向下游。,李光炽,计算水力学,缓流 和急流由于是缓流，因而水流的佛汝德数 ，,36,李光炽,计算水力学,t,x,C,D,L,B,A,T,缓流情况下边界处特征线方向示意图,李光炽计算水力学txCDLBAT缓流情况下边界处特征线方向示,37,李光炽,计算水力学,x,t,C,D,L,B,A,T,急流情况下边界处特征线方向示意图,李光炽计算水力学xtCDLBAT急流情况下边界处特征线方向示,38,