赞
踩
根据三大定律建立流体运动的质量、动量和能量方程,包括微分形式和积分形式。
在适当的定解条件下可以获得详细的流场参数分布。
根据方程的积分形式,包括积分形式的质量、动量和能量方程和欧拉方程的伯努利积分,比较方便地解决流动问题。
根据质量守恒定律,系统的质量不随时间变化,对于流体系统则是质量的随体导数为零,即Dm/Dt=0。
以牛顿第二定律表达的动量守恒定律就是F=ma=mDV/Dt。
采用热力学第一定律表述的能量守恒定律为:单位时间内外界传给系统的热量,等于系统的总能量的增加率加上系统对外界输出的功率,即Q·=DE系统/Dt+W·系统输出。或写为:DE系统/Dt=Q·+W·。此处W·为外界对系统的做功率,包括彻体力的做功率和表面力的做功率。
研究流体微团有两种流体模型。
两种方式:描述流场中每一点的流动细节或针对一个有限区域,通过研究某物理量流入和流出的平衡关系来确定总的作用效果。前一种方法称微分方法,后一种称为积分方法。
基本方程可以用两种形式表示。微分形式针对无限小的流体微团,或无限小体积的微元体,描述流场每点细节;积分形式的方程针对有限大体积的系统,或有限大体积的控制体、区域。
要了解流场各点细节是,需要采用微分形式的方程;如果流动过程在某处发生不连续变化,微分形式的方程失效,必须采用积分形式的方程。
当只要求了解流体动力学问题的总体性能关系,而不要求了解流动过程的详细情况时,可用积分形式的基本方程求解,这种方法简单方便。积分形式的方程只要求方程只要求积分本身存在,适用于有间断的流动。
解决将描述“系统”的微分基本方程转化为描述“控制体”的方程的问题。
需要将“系统”物理量的随体导数用“控制体”中物理量的当地时间导数和单位时间从“控制体”界面流出流出的物理量表达。
N=∭η(x,y,z,t)dΩ
D(∭ηdΩ)/Dt=DN/Dt=N(AB^·^CD^·^,t+dt)-N(ABCD,t)/dt
N(AB^·^CD,t+dt)-N(AB^·^CD,t)=dt∭δη(x,y,z,t)/δtdΩ
N(ADCD^·^,t+dt)-N(ABCB^·^,t+dt)=dt∯η(x,y,z,t)V(x,y,z,t)·nδS
DN=dt(∭δηδΩ/δt)+dt∯η(x,y,z,t)V(x,y,z,t)·nδS
D(∭η(x,y,z,t)dΩ)/Dt=∭δη(x,y,z,t)δΩ/δt+∯η(x,y,z,t)V(x,y,z,t)·nδS
用随体观点表述为:封闭“系统”内流体的质量既不能产生,也不能消失,即“系统”的质量随时间的变化率为零。
用当地观点表述为:“控制体”内的流体质量经过某时间段后的变化(增加)应等于该时间段内通过控制面净流入“控制体”的质量。
Dρ/Dt+ρ▽·V=0
Dρ/Dt=0
,即δρ/δt+ρ▽·V=0
,这就是不可压流动的连续方程。▽·V=0
,即δu/δx+δv/δy+δw/δz=0
。δρ/δt+δ(ρu)/δx+δ(ρv)/δy+δ(ρw)/δz=0
,采用散度符号表示为:δρ/δt+▽·(ρV)=0
D(∭ρ(x,y,z,t)dΩ)=0
,可以认为是采用随体观点的积分形式的质量方程,系统的体积和界面积在运动过程中都可能发生变化。∭δρ(x,y,z,t)δΩ/δt+∯ρ(x,y,z,t)V(x,y,z,t)·nδS=0
,式中第一项是对控制体的体积Ω进行积分,第二项是对控制体的外界面进行面积分,控制体的体积和界面积都固定不变。∭[δρ/δt+▽·(ρV)]δΩ=0
。该式对任何大小的控制体均成立,因此被积函数必须为零。∯V(x,y,z,t)·ndS=0
-∬SinV·nδS=∬SoutV·nδS
。式中,n依然是沿封闭曲面S的外法向。上式表明,对于不可压均值流体流动,当不存在内部源时,通过控制面流进和流出控制体的流体体积流量相等。对于定常流动和非定常流均适用。∯ρ(x,y,z,t)V(x,y,z,t)·ndS=0
,可以改写为-∬SinρV·nδS=∬SoutρV·nδS
。以“流体微团”为分析对象建立。
可以表述为:流体微团的质量与加速度的乘积等于该微团所受外力的合力。
FP合=(δPx/δx+δPy/δy+δPz/δz)dxdydz
。FP合=dxdydz[(δσx/δx+δτyx/δy+δτzx/δz)i+(δτxy/δx+δσy/δy+δτzy/δz)j+(δτxz/δx+δτyz/δy+δσz/δz)k]
。ρDV/Dt=ρf+δPx/δx+δPy/δy+δPz/δz
ρDu/Dt=ρfx+δσx/δx+δτyx/δy+δτzx/δz
,ρDv/Dt=ρfy+δτxy/δx+δσy/δy+δτzy/δz
,ρDw/Dt=ρfz+δτxz/δx+δτyz/δy+δσz/δz
。τyx=μδu/δy
。τxy的下标y代表切应力的法向沿y轴,下标x代表哦切应力的方向沿x轴。该流动τxy>0,据此判断应力方向:在固定壁处,固壁的外法向向上(沿正y),固壁受到流体施加的切应力沿正x方向;固定壁处流体微团表面的外法向向下(沿负y),流体受到固壁施加的切应力指向负x方向。τyx=2μSyx
,应力与应变率成比例。上式通常称为牛顿黏性应力公式。P^-^=-(τx+τy+τz)/3=p
,导出b=-2μ▽·V/3-p
,从而有∏==2μS=-(2μ▽·V·E=)/3-pE=
黏性流体运动微分方程(微分形式的黏性流体动量方程),简称N-S方程。
微分形式的动量方程适用于流场中每一点。若所要求的只是总合的作用力,这时就可以转为研究有限大体积的流体系统或控制体,建立积分形式的动量方程来解决问题。
D(∭ρ(x,y,z,t)V(x,y,z,t)dΩ)=∑F
。∭δ(ρV)δΩ/δt+∯(ρV)V·nδS=∑F
。F体=∭ρfδΩ
。FS面=-∯pnδS+FS面,黏性力
,则作用在控制体上的总和力为:∑F=∭ρfδΩ-∯pnδS+F物+FS面,黏性力
。∭δ(ρV)dΩ/δt+∯(ρV)V·nδS=∭ρfδΩ-∯pnδS+F物+FS面,黏性力
。积分形式的动量方程左端的体积分项表示控制体内流体动量的当地时间变化率,一般情况下很难确定,所以方程通常用于定常流动,此时体积分项等于零。
积分形式的动量方程是流体力学中最常用的基本方程之一,其优点在于:流动定常的情况下,只要知道控制体进出口(控制面)的流动情况,就可以求得总合的作用力,无须知道控制体内部的流动细节。
方程有以下几种应用方法。
第一种是在控制体内包有一个物体(如机翼),需求的是作用在这个物体上的合力(既包括压力,也包括黏性力)。由于控制体外表面S上的黏性力可以忽略,这时只要知道外表面S上的流速及压强分布情况,就能理算出物体的受力。
第二种是气流流经一段管道,管内可以有其他物体,也可以是空的,这时管壁和管内机件表面合起来作为控制面S1,进口和出口两个截面合起来为式中的控制面S,只要知道了进出口两个截面上的流动参数,就可以求得管壁和管内机件受到的合力。
第三种是在超声速气流的激波前后。
采用积分形式的方程解决问题时,控制体的选取是关键。一般原则是:控制面必须是封闭的,必须包围研究对象;控制面形状尽可能简单,尽可能利用全部已知量;尽可能便于计算,例如,采用流线或流面做控制面、用对称边界面、取平面边界等。
积分形式的动量方程适用于流体的直线运动或一般曲线运动。
当流体绕定轴转动时,采用动量矩方程更方便,当不仅想知道作用力,还想知道合力作用点时,也需要用到动量矩方程。
动量矩守恒定理表述为:系统对某点的动量矩随时间的变化率,等于此时作用在系统上的所有外力对该点的力矩之和。
对于理想无黏流体,黏性力为零,此时N-S方程简化为欧拉方程:
∂u,v,w/∂t+u,v,w∂u/∂x+v∂u,v,w/∂y+w∂u,v,w/∂z=-∂p/ρ∂x,y,z+fx,y,z
写成矢量形式为:∂V/∂t+(V·▽)V=-▽p/ρ+f
。
这就是理想无黏流体的运动微分方程,对于定常和非定常流动、可压和不可压流动均适用。
欧拉运动微分方程给出了压强的变化、速度的变化以及彻体力之间的关系。需要指出,在空气动力学中,重力作为一种彻体力一般可忽略不计。
兰姆方程显式地包含速度旋度的理想流体运动微分方程,将欧拉方程中速度的迁移导数进行分解,将其中与旋度有关部分分离出来后导出的。
V×(▽×V)=▽(V·V/2)-(▽·V)V
,可以将矢量形式的欧拉方程写为:∂V/∂t+▽(V^2/2)-V×Ω=-▽p/ρ+f
,上式称兰姆型运动方程。Ω为速度的旋度,也是涡量,Ω=2w=rotV。f=-▽∏
,式中,∏是彻体力f的势能函数。∂V/∂t-V×Ω=-▽(P+∏+V^2/2)=-▽H
,其中H=P+∏+V^2/2,H可视为单位质量流体的总机械能。V×Ω=▽(P+∏+V^2/2)
。V·▽(P+∏+V^2/2)=0
,Ω·▽(P+∏+V^2/2)=0
,所以有P+∏+V^2/2=C(沿流线)
,P+∏+V^2/2=C(沿涡线)
。式中,C为常数。上式就是对理想、定常、有旋流动导出的伯努利积分,适用于某条流线或涡线。上式中的积分常数C称为伯努利常数,沿同一流线,积分常数才相等,对不同流线或涡线,积分常数互不相同。p/ρ+∏+V^2/2=C(沿流线)
。z+p/ρg+V^2/2g=C(沿流线)
。上式是水力学应用最广的伯努利方程。它的意义可以从几何和能量两个角度说明。▽(p+∏+V^2/2)=0
,则P+∏+V^2/2=C(全流场)
。对不可压流,有p/ρ+∏+V^2/2=C(全流场)
。上式是适用于全流场的,全流场的C为一个常数;而有旋的伯努利方程积分常数C则是适用于一条流线的,流线改变后,常数C要变。欧拉积分通常也称为伯努利积分,作为其特例。p+V^2/2ρ=C
,上式是在低速气流中经常使用的伯努利公式。C常写为p0,称为总压,则p0=p+V^2/2ρ
。相应地,称p为静压,V^2/2ρ为动压(也称为速度头)。总压可看成在速度滞止为零时全部动能转化为压力势能后的总机械能。Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。