传递过程原理讲课提纲04粘性流体运动的微分方程及其应用2百精.docx
- 文档编号:11530438
- 上传时间:2023-03-18
- 格式:DOCX
- 页数:28
- 大小:62.11KB
传递过程原理讲课提纲04粘性流体运动的微分方程及其应用2百精.docx
《传递过程原理讲课提纲04粘性流体运动的微分方程及其应用2百精.docx》由会员分享,可在线阅读,更多相关《传递过程原理讲课提纲04粘性流体运动的微分方程及其应用2百精.docx(28页珍藏版)》请在冰豆网上搜索。
传递过程原理讲课提纲传递过程原理讲课提纲04粘性流体运动的微分方程及其粘性流体运动的微分方程及其应用应用2百精百精24NS方程在柱坐标及球坐标中的表示1.柱坐标中的表示x=rcosy=rsinz=z在r分量方向zuuruururuuurzrrrr+-+2=(+-+-2222222111zuururrurrrrpXrrrr在分量方向zuuruuururuuuzrr+=(+-2222222111zurururrurrrprXr在z分量方向zuuururuuuzzzzrz+=1(1122222zuurrurrrzpXzzzz+-2.球坐标中的表示x=(rsincosy=(rsinsinz=rcosr分量:
ruuuruururuuurrrrr22sin+-+sin1(2sin1sinsin11122222222+-+-uctguuurururrurrrrpXrrrrr=yy图19分量:
rctguruuuruururuuurr2sin-+sincos2sin2sin1sinsin1111222222222222-+-urruurururrurrrprXr=分量:
rctguuruuuruururuuurrr+sinsincos2sinsin2sin1sinsin11sin112222222222222+-+-urruurururrurrrprXr=3流体运动方程的应用3-1平壁间的稳定层流设平板无限大,相互平行,作层流运动一维稳定流动,不可压缩于是uy=uz=0(1由连续性方程有0=xux(2又对稳定流动0=xu(3故NS方程简化为:
(12222zuyuXxpxx+=对无限大平板可认为0=zux故022=zux,又在x方向X=0x图20于是222211dyuddxdpyuxpxx=边界条件(BoundaryConditiony=y0,ux=0初始条件(Initialconditiony=0,0=dydux于是(2221yydxdpu-=及-=20max1yyuu式中:
dxdpyu-=20max21又设平板的宽度为w,则流体流过二平行平板间的体积流量为-=AyvydxdpwudywudAQ0300322另一方面,若设平板间的主体流速(即截面平均流速为ub则有QV=ubB(2y0可得:
2031ydxdpub-=故32=xbuu及203yudxdpb-=3-2圆形直管内的稳定层流化工原理中已得出了相应的结论。
若采用NS方程可知:
在沿轴的一维流动时,若处于重力场中:
000=zuuuuzzxy及00=xpypz,故NS方程可简化为:
+=22221yuxuzpyxA图21应用柱坐标系则为:
dzdpdrdurdrudzz=+1122或dzdpdrdurdrdrz11=边界条件r=Ruz=0初始条件r=0uz=umax故:
-=2214RrdzdpRuz,dzdpRu-=42max及-=2max1Rruuz,max21uub=,28Rudzdpb-=上式均为哈根泊稷叶(HagenPoiseuille方程。
根据以上关系,不难证明,在稳定的层流时,管壁处的流体剪应力drduzw-=R-=dzdpZR及(Re64Re1641=f3-3环形套管中的稳态层流圆形直管中的柱坐标系NS方程仍可适应于环形套管中:
dzdpdrdurdrudzz=+1122或dzdpdrdurdrdrz11=边界及初始条件:
r=R1,uz=0;(2r=R2,uz=0;(3r=Rmax,uz=umax;及0=drduz对前式积分得:
xz图22z图23-=12max212ln221RrRRrdzdpurz或-=22max222ln221RrRRrdzdpurz根据以上二式有:
21221222maxln(RRRRR-=此即为对数平均“半径平方”,于是:
-=2max2max222maxmaxln221RRRRRdzdpu及:
-=1max2max212maxmaxln221RRRRRdzdpu环形面积中流体的主体流速(平均流速ub(2max2122281222121RRRdzdprdrdrrudAdAuuRRRRrrb-+-=或:
-+-=2max2122218RRRudzdpb4爬流(蠕动流1、概念爬流非常低Re数下的流动。
具体来说指Re1时的流动,在此情况下,流体流动过程的粘性力超过其惯性力。
2、根据NS方程可知,爬流时,惯性力的影响比粘性力小得多而可忽略。
故有:
对不可压缩流体0=+zuyuxuzyx(不可压缩流体连续性方程X,Y及Z0(惯性力DDux,DDuy及DDuz0(惯性力故NS方程可简化为:
+=222222zuyuxuxpxxx+=222222zuyuxuypyyy及+=222222zuyuxuzpzzz当流体以极慢的速度绕过半径为r0的球体时,采用球坐标方程并结合球坐标系中的连续性方程及初始条件:
r=r0时,0000=uuur可解得+-=300021231cosrrrruur-=300041431sinrrrruu2000023cos-=rrrupp3、球体在流体中所受阻力包括:
a.由于压力分布在球体表面所产生的形体阻力Fdfb.球体表面处的剪应力所产生的表面阻力Fds对形体阻力,其方向在“-Z”方向,其大小为(=-=20000202sincosuddrpFrdf0:
20:
对表面阻力,其方向在+Z方向(即与球面法向应力相反的方向其大小为:
(=20020sinsinddrFrdfr-作用于球面的剪应力,因pr,故-pcos恰与rsin在同一方向P0r图24xy图25对于球坐标系,牛顿粘性定律dydu-=可表示为+-=rrurrur1由上述已知的解,代入方程得sin2300rur=故:
Fds=4r0u0总阻力:
Fd=Fdf+Fds=6r0u0此即为斯托克斯方程(注意与化工原理中的斯托克斯方程相联系通常对粒子的自由沉降,定义220uCAFDd=,式中CD为阻力系数(根据以上结论可得:
CD=Re24220=AuFd式中A球形粒子在运动方向的投影面积A=r025流线及流函数1、概念流线欧拉法考察流体运动的结果。
即指任一瞬间,在流场(流动空间的一条曲线,处于该曲线上的各质点流动方向与该点处曲线的切线相重合(流线既可以为直线,也可为曲线。
2、流线的几个属性一个流动空间是由许多条流线组成的,这许多流线通常称为流线束(曲线族;不同时间,流动空间中的流线束有可能不同(如不稳定流动过程;同一时间,同一空间位置,由于流体质点速度大小及方向的唯一性,故通过该空间点的流线唯一。
亦即,同一时间、空间中的不同流线不可能相交。
3、流线方程对质点A,在时间d内,有:
dx=uxddy=uydxdz=uzd故zyxudzudyudx=u质点A此即为三维流动系统的流线微分方程。
特殊地,对二维流动,则有:
uxdy-uydx=0根据这一关系,可以求出流体在空间流动的流线方程。
即y=f(x关系(见例题4、流函数概念:
相对于某基准流线而言的流体体积流量,设有一函数满足:
xuyuyx=-=则d=dxudyudyydxxyx-=+即定义为流函数。
对二维不可压缩流体,则:
0(022=+=-yuxuxyyxyx即(不可压缩流体连续性方程用流函数表示的结果亦可表述为:
穿过由基准流线和任一流线及垂直于纸面方向上的单位厚度流道所构成的体积流量。
即dQ=A-B=d不难看出:
当=常数时,其所表示的即为流线xyudxudy=;流速越大的地方,流线越密集;流体流过曲线c的单位厚度流量等于曲线c上A、B两点的流函数差;5柱坐标中流函数的定义式定义:
rdur=ru-=6势流及势函数1概念势流函数:
即速度势(VelocityPotential函数。
根据势能的概念可知:
在重力场中,单位质量流体(固体势能的变化d等于将流体升举一个微分高度所做的功,即:
d=-gdz或-g=d/dz(注意:
流体势能的意义p+gz线图27将此概念引入速度问题,规定一个函数,并定义:
xux=可见,流体沿x方向的运动速度是其在x方向的速度势梯度。
2用势函数表示的不可压流体的二维流动连续性方程直角座标中,不可压流体的二维流动连续性方程为:
0=+yuxuyx引入速度势概念则有02222=+yx(不可压流体连续性方程用势函数表示此即为Laplace方程。
通过引入Laplace变换,并已知适当的边界条件,即可求解。
3根据势函数xux=及yuy=故有:
xuyuyx=(xyyx=22即由右图28分析可知:
a、平面上微元流体在经过d时间后,由于x方向速度梯度使上层流体运动所引起的位移为dydyux(顺时针方向;b、同理,在y方向引起的位移为dxdxuy(逆时针方向。
这二者使流体微元绕z轴方向发生旋转。
故总的旋转速度为二者之和的平均值。
即:
yuxuxuyunxyyx-=+-=20(以逆时针方向为正计当=0时,称此流体为无旋运动,或说为有势运动。
即此时:
yuxuxy=4势函数的应用举例泊努利方程的导出x由欧拉方程(二维xpXuzuuyuuxuuxxzxyxx-=+1(1ypYuzuuyuuxuuyyzyyyx-=+1(2ypZuzuuyuuxuuzzzzyzx-=+1(3在稳态、二维运动及不可压缩流体条件下有:
xpXyuuxuuxyxx-=+1(1ypYyuuxuuyyyx-=+1(2又由有势(无旋运动xuyuyx=有:
xpXxuuxuuyyxx-=+1(1”ypXyuuyuuyyxx-=+1(2”即:
puuyx+-+222121=常数或:
-+pu22=常数此即为著名的柏努利(Bernoulli方程。
5讨论由有势运动(无旋运动可知:
只有对于理想流体流动才会有势函数(无旋运动存在。
对实际流体,由于其粘性,将会使流体微元发生有旋运动。
尽管实际流体不存在势函数,但其流函数仍然存在。
只要该流体满足连续性方程,而并不需要满足流体不可压缩条件。
6流网、势线与流线的关系流网:
由一组等势线和一组流线所交织成的图称为平面有势(无旋运动流网。
流网的特性:
等势线与流线必然正交。
证明:
对等势线,其方程为(x,y=常数或dx+dy=0uxdx+uydy=0(法线向量)xy(切线向量)对流线:
uxdyuydx=0即:
=ux(流线)y而对势线:
=ux(势线)x又在直角座标中:
xy故:
=即与正交。
xy7应用举例例A.对于稳态二维流动的流体,已知某流场的速度向量形式为u(x,y=3x2yi+4xy3j试求出过点(1,3)的流线方程;判别其无旋性(有势性)。
解:
流线方程的一般形式为由已知向量形式可知:
对稳定二维流动d=0即:
3x2ydy4xy3dx=0故d=uxdyuydxux=3x2yuy=4xy3ydy2=4dx3x即:
114=lnx+c或y=4y3lnx+c3314lnx又该流线过点(2,3)故c=1/3y=此即为过点(1,3)的流线方程。
判别有旋性(有势性)若流体有势,则应满足即:
uxdx+uydy=0udy=xdxuy或uyx=uxy34对方程求导:
4dy3(14lnx43x=y=22dx(14lnx(14lnx3x(14lnx2又:
22ux3x2y3x3(14lnxx(14lnx=2=x=uy49124xy34y可见在一般情况下,该流体均为有旋运动,即为非理想流体。
例B.已知二维流场中,稳态流动下的速度向量为u(x,y=3xy2i+3x2yj,且流线过点(1,。
2)试问其是否作无旋运动。
若无旋,试求其势函数及流函数,并证明与正交。
解:
由题可知:
ux=3xy2,uy=3x2y若作无旋运动,则应满足:
uyx=uyux=6xyyxxux=6xyy322xy+c(y2由此可判断其作无旋运动(有势运动)又:
ux=3xy2=故=3x2y+c(y=uy=3x2yy故:
于是:
c(y=0即c(y=常数=c=322xy+c势函数2对于稳定的二维运动,其流函数应满足:
uxdyuydx=0即:
3xy2dy3x2ydx=0即:
ydy=xdxy2=x2+c当其过点(1,2)时有:
22=12+cc=3即:
y2=x2+3或y=x2+3正交性若流线的切线与势线的法线二者斜率相等,则与正交对流线:
y2=x2+3故2yy=2x即:
k1=y=x/y流线切线斜率对势线:
=3/2x2y2+C在等势线(=C1)上,有:
C1=3/2x2y2+C即:
x2y2=C2故2yy=2xcc1=422xx35即:
y=c2又由于x=1时y=2故C2=4yx3即:
y=44y4yy=22=3c2xxyxyxx1x=yyx则势线的法线斜率k2=k1*k2=-1,即与正交。
36
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 传递 过程 原理 讲课 提纲 04 粘性 流体 运动 微分方程 及其 应用 百精