1、理想流体的动量传递一欧拉运动方程冶金机械 流体的运动方程是牛顿第二运动定律在流体流动冶金机械现象中的应用。牛顿第d(mu)描述,即: Fdt定律用动量来,因为作用力与动量率相当,故流体的运动方程也称动量传递方程。对密度为“的理想流体,在直角坐标系下从流场中取一微元六面体,其边长分别为 ax,:和de, 顶角 A 的坐标为(x,:,z) ,如图 1-2-22 所示。根据牛顿第二运动定律,对该微元体有作用在微元体_流出微元体流入微元体_微元体内动_-上的合外力. = 的动量率 的动量率 + 量率的积累(1-2-69)率 L量率的积累 下面首先分析流体沿 x 轴方向流动的动量传递现象。(1)微元体
2、的受力分析作用于该微元体上的力有表面力和质量力两种。对理想流体而言,表面力只有压力,粘性力为零;质量力有重力、电磁力等。设微元体顶角 A 处的压力为?则单位时间内在 x 方向作用于微元体上的净压力为ABCD 面和 EFGH 面所受压力之代数和,即p-(p+#dx)d:dz=-dxayz(1-2-70)2图 1-2-22 理想流体的动量守恒方程推导示意图设单位质量流体的质量力沿 x 轴方向的分量为 X,则微元体在 x 轴方向所受的质量力为pXdxdyz (1-2-71)x 轴方向微元体所受外力之合力为(1-2-72)(2)微元体的动量分析设微元体顶角 A 处的流体质点在 x,y,z 轴方向的对流
3、动量通量分别为 !ux ux, Puxuy和 pUxUz。在 x 轴方向,从 ABCD 面流入和从 EFGH 面流出微元体的动量率分别为: Adxdydz在 x 轴方向流出微元体的净动量率为二者之差,即(PUxUx11(PUxUx)11pUxUx+ -“x-dxjdydz-pUxUxdydz= “x-dxdydz同样,在 y 轴方向,从 ADHE 面流入和从 BCGF 面流出微元体的沿 x 方向的动量率分别为:pUxUydxdz 和pUxUy4(pU xUydxdxdz。因此,在 y 轴方向流出的沿 x 轴方向的净动量率为:(pUxUy)“ydxdydz同理,在 z 轴方向流出的沿 x 轴方向
4、的净动量率为:“(pUxU2)nndxdydz“z三者之和即为 x 轴方向总的净动量率:(pUxUx)(pUxUy)! “V + “y(pUxx)“(U xUy)“(“x又微元体内动量率的积累为:PUxUz“zdxdydz(1-2-73)“tdxdydz(Ul)将式-二-))、!-#-)和-二-)) 代入式(1-2-69),则有(V 1 f“(PUxUx)“(PUxUv)“(PUxUz)l!1 丄“(PUx)11l!入一“x/dxdydz=L “-+ -“+ “ Jdxdydz4dxdydz即“p“(pu xux)“(puxuY)“(puxu2)“(pux)!X_“x=“x“y “x“z“t(
5、1-2-75)将上式右端的各项展开 pX-“p“x=UxL“x“(pUx) “(pUy)“(pu.I“v “y “z“t“Ux “Ux “Ux“Ux,“x +Uy“(+Uz“将连续性方程式(1.57)代入上式,则有“&x1“p_“Uxp“x“t“UxX“r=可 +Ux“x+Uy“yUz“Ux“z同理,可导出沿 y,z 轴方向的分动量守恒方程分别为:“U(“U(“U(丄“#_“Uyp“y“t3“丄“p_“ Uzp“z“t考虑到速度 U 的随体导数为DU_“UDt_“t因此,式(1“2-76) 可改写为XUx“xy“yz“z“Uz“UzUx“x+Uy“(+UZ“Uz(1-2-76a)(1-2-76
6、b)(1-2-76c)“U“UY-U 化 =DtM=Dt 蒂Uax 选办远巧 I“IIPljpIIpUUz“U“z(1-2-77a)1-2-77b)(1-2-77c)上述三式可统一写成向量形式,即r1。DuF-#“p_pFDt(178)式(1-2-76)(1-2-78)即为理想流体的运动微分方程,该方程由欧拉于 1755 年首先提出,因此,也称为欧拉运动方程。它是研究理想流体运动的基础,对稳定流动与不稳定流动、可压缩流体与不可压缩流体都适用。若流体在重力场中作稳定流动,质量力只有重力,则 XzO.YzO.Z:.76)可简化为g,!&=0,式(1“丄 “p=p“x“Uxx“x“UxUy“y+“UxUz“z (1-: -79a)“丄 “p_p“y“Uy&x“x+“Uy&y“y+“Uy&z“z(1-: 1-79b):丄“p “&z “Uz “&z+UzTT- (1-: Z-g-p“z _&x“x+&y“yz“z 79c)上述方程组(1- 79)存在四个未知数,ux,Uy,z 和p,已有三个方程,加上连续性方程,则可以求解。若流体为可压缩性气体,则 p将随 p 而变,需联入气体状态方程才能求解。原文地址:http:/