高等流体力学复习02

流体变形包括哪几种?变形(线变形/角变形)、平移、旋转

一 连续性方程

1.取控制体(任意)
2.质量守恒定律
3.高斯公式-将面积分转化为体积分(一个函数的面积分等于其散度的体积分)
4.积分定理-连续函数,任意积分体,其体积分为零,则函数必为零
V [ ρ t + ρ V n ] d V = 0 \int_V[\frac{ \partial\rho}{\partial t} +\bigtriangledown ·(\rho V_n) ]{\rm d}V=0
5.定常流
ρ t = 0 \because \frac{ \partial\rho}{\partial t}=0
ρ V n = 0 \Rightarrow \bigtriangledown ·(\rho V_n)=0
6.定常不可压缩流
ρ t = 0 ρ = c o n s t \because \frac{ \partial\rho}{\partial t}=0,\rho=const
u x + v y + w z = 0 \Rightarrow\frac{ \partial u}{\partial x}+\frac{ \partial v}{\partial y}+\frac{ \partial w}{\partial z}=0
思考题:
1.试判别下列流动是否存在 u = x 2 y 2 , v = 2 x y , w = 0 u=x^2-y^2,v=-2xy,w=0
2.试给出圆管中的连续性方程
3.试给出喷管中的连续性方程

二 运动方程(动量方程)

1.动量定理:
(1)冲量–力在时间过程中的累积效应,等于力与其作用时间的乘积,为矢量,方向与力一致。
(2)动量:是表征物体运动状态的量,等于质量与速度的乘积,为矢量,方向与速度一致。
(3)动量定理
2.取控制体(任意)
3.受力分析
(1)体积力 V ρ F d V \int_V \rho·F {\rm d}V
(2)表面力 A τ n d A \int_A \tau_n{\rm d}A
(3)动量变化率 d d τ V ρ V n d V \frac{d}{d\tau}\int_V \rho·V_n{\rm d}V
4.由动量定理
d d τ V ρ V n d V = V ρ F d V + A τ n d A \frac{d}{d\tau}\int_V \rho·V_n{\rm d}V=\int_V \rho·F {\rm d}V+\int_A \tau_n{\rm d}A
5.由高斯公式得:
d d τ V ρ V n d V = V [ ρ F + ( τ n ) ] d V \frac{d}{d\tau}\int_V \rho·V_n{\rm d}V=\int_V [ \rho·F+\bigtriangledown·(\tau_n) ] {\rm d}V
6.函数积分的随体导数
d d τ V ρ V n d V = V d ( ρ V n ) d τ d V \frac{d}{d\tau}\int_V \rho·V_n{\rm d}V=\int_V\frac{{\rm d}(\rho·V_n)}{{\rm d}\tau}{\rm d}V
7.积分定理-控制体任意且连续
d ( ρ V n ) d τ ρ F τ n = 0 \frac{{\rm d}(\rho·V_n)}{{\rm d}\tau}-\rho·F-\bigtriangledown·(\tau_n)=0
8.方程的封闭性
9.纳维-斯托克斯方程(N-S方程)
10.N-S方程讨论:
(1)方程的封闭性
(2)理想流体
(3)静止流体
练习:已知u=yzt,v=xzt,w=0,求质点在(1,2,1)处,t=1s时的加速度。

另一种推导动量方程的方法 —微分形式

牛顿第二定律 (受力平衡)
合力=体积力+表面力=ma
F = F b + F s = m a = ρ d x d y d z d V d t \vec F=\vec F_b+\vec F_s=m\vec a=\rho{\rm d}x{\rm d}y {\rm d}z\frac{{\rm d }\vec V}{{\rm d }t}
其中,体积力为: F b = f b ρ d x d y d z \vec F_b=\vec f_b·\rho {\rm d}x{\rm d }y{\rm d }z
表面力: F s , r i g h t F s , l e f t = ( τ x + τ x x d x ) d y d z τ x d y d z = τ x x d x d y d z \vec F_{s,right}-\vec F_{s,left}=(\vec\tau_x+\frac{\partial\vec\tau_x}{\partial x} {\rm d}x) {\rm d}y {\rm d}z-\vec\tau_x {\rm d}y {\rm d}z=\frac{\partial\vec\tau_x}{\partial x} {\rm d}x{\rm d}y {\rm d}z

F s = ( τ x x + τ y y + τ z z ) d x d y d z \Rightarrow \vec F_s=(\frac{\partial\vec\tau_x}{\partial x} +\frac{\partial\vec\tau_y}{\partial y} +\frac{\partial\vec\tau_z}{\partial z} ) {\rm d}x{\rm d}y {\rm d}z

d V d t = f s + 1 ρ ( τ x x + τ y y + τ z z ) \Rightarrow \frac{{\rm d }\vec V}{{\rm d }t}=\vec f_s+\frac{1}{\rho}(\frac{\partial\vec\tau_x}{\partial x} +\frac{\partial\vec\tau_y}{\partial y} +\frac{\partial\vec\tau_z}{\partial z})
任一表面应力可分解为3个方向的应力张量:一个正应力,二个切应力
τ x = τ x x + τ x y + τ x z = τ x x i + τ x y j + τ x z k \vec \tau_x=\vec \tau_{xx}+\vec \tau_{xy}+\vec \tau_{xz}=\tau_{xx}\vec i+\tau_{xy}\vec j+ \tau_{xz}\vec k
d V d t = f s + 1 ρ ( τ x x x + τ y x y + τ z z ) + 1 ρ ( τ x x + τ y y + τ z z ) + 1 ρ ( τ x x + τ y y + τ z z ) \Rightarrow \frac{{\rm d }\vec V}{{\rm d }t}=\vec f_s+\frac{1}{\rho}(\frac{\partial\vec\tau_{xx}}{\partial x} +\frac{\partial\vec\tau_{yx}}{\partial y} +\frac{\partial\vec\tau_z}{\partial z})+\frac{1}{\rho}(\frac{\partial\vec\tau_x}{\partial x} +\frac{\partial\vec\tau_y}{\partial y} +\frac{\partial\vec\tau_z}{\partial z})+\frac{1}{\rho}(\frac{\partial\vec\tau_x}{\partial x} +\frac{\partial\vec\tau_y}{\partial y} +\frac{\partial\vec\tau_z}{\partial z})
张量表示: d u i d t = f b , i + 1 ρ ( τ j i x j ) \frac{{\rm d}u_i}{{\rm d }t}=f_{b,i}+\frac{1}{\rho}(\frac{\partial \tau_{ji}}{\partial x_j})
又因为: τ x y = τ y x , τ x z = τ z x , τ y z = τ z y \tau_{xy}=\tau_{yx},\tau_{xz}=\tau_{zx},\tau_{yz}=\tau_{zy}
则应力分量应该有6个独立分量:
τ x x , τ y y , τ z z , τ x y , τ x z , τ y z \tau_{xx},\tau_{yy},\tau_{zz},\tau_{xy},\tau_{xz},\tau_{yz}
另外:欧拉认为流体的粘性力相比与压力是非常小的,假设忽略切应力,正应力等于压力,即: τ x x = τ y y = τ z z = p , τ x y = 0 , τ x z = 0 , τ y z = 0 \tau_{xx}=\tau_{yy}=\tau_{zz}=-p,\tau_{xy}=0,\tau_{xz}=0,\tau_{yz}=0
d V d t = f s 1 ρ p \Rightarrow \frac{{\rm d }\vec V}{{\rm d }t}=\vec f_s-\frac{1}{\rho}·\bigtriangledown p----无粘流动欧拉方程

三 能量方程

1.取控制体(任意)
2.能量分析
(1)动能变化率(分子热运动+质点运动)
d E v d τ = d d τ V ρ ( e + V n 2 2 ) d V \frac{{\rm d}E_v}{{\rm d }\tau}=\frac{{\rm d} }{{\rm d }\tau}\int_V \rho(e+\frac{V_n^2}{2}) {\rm d }V
(2)功(体积力功+表面力做功)
W = V ρ F V n d V + A τ n V n d A W=\int_V\rho·F·V_n{\rm d}V+\int_A\tau_nV_n{\rm d}A
(3)传热(导热项+源项) Q = A λ T n d A + V ρ q d V Q=\int_A\lambda\frac{ \partial T}{ \partial n}{\rm d }A+\int_V\rho q{\rm d}V
3.由能量守恒定律的:
ρ C p d T d τ = ϕ + λ T + ρ q \rho C_p\frac{{\rm d}T}{{\rm d}\tau}=\phi+\lambda\bigtriangledown·T+\rho q
其中, ϕ = \phi=
λ = λ μ \lambda^·=\lambda^·(\mu)

思考题:
(1).试给出常见的理想流体能量方程—(伯努利方程)
(2).当沿程存在阻力损失时的能量方程
(3).当存在损失和外功输入的能量方程
练习:
1.水泵
2.
3.
4
虹吸现象
毛细现象

猜你喜欢

转载自blog.csdn.net/weixin_43466418/article/details/85050679