能量方程
这里仍然采用随流体运动的无穷小流体微团为研究对象推导能量方程。根据热力学第一定律,对于随流体运动的微团模型而言,其满足能量守恒定律,即:流体微团内能量的变化率=流入微团内的净热流量+体积力和表面力对微团做功的功率
首先来看体积力以及表面力对微团做功的功率。 可知作用在一个运动物体上的力,对物体做功的功率等于这个力乘以速度在此力作用方向上的分量。所以作用在速度为
V
V
V的流体微团上的体积力,做功的功率为:
ρ
f
.
V
(
d
x
d
y
d
z
)
\rho f.V(dxdydz)
ρf.V(dxdydz)
对于表面力(压力加上正应力和切应力)做功的功率,只考虑作用在x方向上的力。在x方向上压力和切应力对流体微团做功的功率,等于x方向上速度的分量u乘以力(比如,在面abcd上的切应力为
τ
x
y
d
x
d
y
\tau_{xy}dxdy
τxydxdy,相应的功率表示为
μ
τ
x
y
d
x
d
y
\mu \tau_{xy}dxdy
μτxydxdy)。在进行推导之前,这里约定作用在x正向上的力做正功,在x负方向上的力做负功。与前面一节推导动量方程的方法类似,这里我们根据上面的推导可得到作用在面adhe以及bcgf上的压力在x方向上做功的功率为:
[
u
p
−
(
u
p
+
∂
(
u
p
)
∂
x
d
x
)
]
d
y
d
z
=
−
∂
(
u
p
)
∂
x
d
x
d
y
d
z
\left[u p-\left(u p+\frac{\partial(u p)}{\partial x} \mathrm{~d} x\right)\right] \mathrm{d} y \mathrm{~d} z=-\frac{\partial(u p)}{\partial x} \mathrm{~d} x \mathrm{~d} y \mathrm{~d} z
[up−(up+∂x∂(up) dx)]dy dz=−∂x∂(up) dx dy dz
类似的,在面abcd以及面efgh上,切应力在x方向上做功的功率为:
[ ( u τ y x + ∂ ( u τ y x ) ∂ y d y ) − u τ y x ] d x d z = ∂ ( u τ y x ) ∂ y d x d y d z \left[\left(u \tau_{y x}+\frac{\partial\left(u \tau_{y x}\right)}{\partial y} \mathrm{~d} y\right)-u \tau_{y x}\right] \mathrm{d} x \mathrm{~d} z=\frac{\partial\left(u \tau_{y x}\right)}{\partial y} \mathrm{~d} x \mathrm{~d} y \mathrm{~d} z [(uτyx+∂y∂(uτyx) dy)−uτyx]dx dz=∂y∂(uτyx) dx dy dz
进一步,可推导出所有的表面力对运动流体微团做功的功率为:
[ − ∂ ( u p ) ∂ x + ∂ ( u τ x x ) ∂ x + ∂ ( u τ y x ) ∂ y + ∂ ( u τ z z ) ∂ z ] d x d y d z \left[-\frac{\partial(u p)}{\partial x}+\frac{\partial\left(u \tau_{x x}\right)}{\partial x}+\frac{\partial\left(u \tau_{y x}\right)}{\partial y}+\frac{\partial\left(u \tau_{z z}\right)}{\partial z}\right] \mathrm{d} x \mathrm{~d} y \mathrm{~d} z [−∂x∂(up)+∂x∂(uτxx)+∂y∂(uτyx)+∂z∂(uτzz)]dx dy dz
上述式子仅仅考虑了x方向上的表面力。接下来,再考虑y以及z方向上的表面力。最终得到体积力以及应力的合力对流体微团做功的功率总和为:
C
=
[
−
(
∂
(
u
p
)
∂
x
+
∂
(
v
p
)
∂
y
+
∂
(
w
p
)
∂
z
)
+
∂
(
u
τ
x
x
)
∂
x
+
∂
(
u
τ
y
x
)
∂
y
+
∂
(
u
τ
z
x
)
∂
z
+
∂
(
v
τ
x
y
)
∂
x
+
∂
(
v
τ
y
y
)
∂
y
+
∂
(
v
τ
z
y
)
∂
z
+
∂
(
w
τ
z
)
∂
x
+
∂
(
w
τ
y
z
)
∂
y
+
∂
(
w
τ
z
)
∂
z
]
d
x
d
y
d
z
+
ρ
f
⋅
V
d
x
d
y
d
z
\begin{aligned} C=&\left [ -\left(\frac{\partial(u p)}{\partial x}+\frac{\partial(v p)}{\partial y}+\frac{\partial(w p)}{\partial z}\right)+\frac{\partial\left(u \tau_{x x}\right)}{\partial x}+\frac{\partial\left(u \tau_{y x}\right)}{\partial y}+\frac{\partial\left(u \tau_{z x}\right)}{\partial z}+\frac{\partial\left(v \tau_{x y}\right)}{\partial x}+\right.\\ &\left.\frac{\partial\left(v \tau_{y y}\right)}{\partial y}+\frac{\partial\left(v \tau_{z y}\right)}{\partial z}+\frac{\partial\left(w \tau_{z}\right)}{\partial x}+\frac{\partial\left(w \tau_{y z}\right)}{\partial y}+\frac{\partial\left(w \tau_{z}\right)}{\partial z}\right] \mathrm{d} x \mathrm{~d} y \mathrm{~d} z+\rho f \cdot V \mathrm{~d} x \mathrm{~d} y \mathrm{~d} z \end{aligned}
C=[−(∂x∂(up)+∂y∂(vp)+∂z∂(wp))+∂x∂(uτxx)+∂y∂(uτyx)+∂z∂(uτzx)+∂x∂(vτxy)+∂y∂(vτyy)+∂z∂(vτzy)+∂x∂(wτz)+∂y∂(wτyz)+∂z∂(wτz)]dx dy dz+ρf⋅V dx dy dz
接下来继续考虑流入微团的总热流量。主要由两部分组成:一)来自于体积加热,如吸收或者释放的辐射热;2)由温度梯度导致的跨过表面的热输送,即热传导。定义
q
˙
\dot{q}
q˙为单位质量的体积加热率。运动微团的质量为
ρ
d
x
d
y
d
z
\rho dxdydz
ρdxdydz,从而可以得到:
微团的体积加热
=
ρ
q
˙
d
x
d
y
d
z
微团的体积加热=\rho \dot{q}dxdydz
微团的体积加热=ρq˙dxdydz
热传导从面adhe输运给微团内的热量为
q
x
˙
d
y
d
z
\dot{q_x}dydz
qx˙dydz。其中
q
x
˙
\dot{q_x}
qx˙是热传导在单位时间内通过单位面积在x方向上输运的热量。因此,经过面bcgf输运到微团外的热量是:
[ q ˙ x − ( q ˙ x + ∂ q ˙ ∂ x d x ) ] d y d z = − ∂ q ˙ x ∂ x d x d y d z \left[\dot{q}_{x}-\left(\dot{q}_{x}+\frac{\partial \dot{q}}{\partial x} \mathrm{~d} x\right)\right] \mathrm{d} y \mathrm{~d} z=-\frac{\partial \dot{q}_{x}}{\partial x} \mathrm{~d} x \mathrm{~d} y \mathrm{~d}z [q˙x−(q˙x+∂x∂q˙ dx)]dy dz=−∂x∂q˙x dx dy dz
进一步,加上其他面在y,z方向上的热输运量,可以得到:
热传导对流体微团的加热 = − ( ∂ q ˙ x ∂ x + ∂ q ˙ y ∂ y + ∂ q ˙ z ∂ z ) d x d y d z \text { 热传导对流体微团的加热 }=-\left(\frac{\partial \dot{q}_{x}}{\partial x}+\frac{\partial \dot{q}_{y}}{\partial y}+\frac{\partial \dot{q}_{z}}{\partial z}\right) \mathrm{d} x \mathrm{~d} y \mathrm{~d} z 热传导对流体微团的加热 =−(∂x∂q˙x+∂y∂q˙y+∂z∂q˙z)dx dy dz
综上所述,总的热量=辐射总热量+热传导热量
[ ρ q ˙ − ( ∂ q ˙ x ∂ x + ∂ q ˙ y ∂ y + ∂ q ˙ z ∂ z ) ] d x d y d z \left[\dot{\rho q}-\left(\frac{\partial \dot{q}_{x}}{\partial x}+\frac{\partial \dot{q}_{y}}{\partial y}+\frac{\partial \dot{q}_{z}}{\partial z}\right)\right] \mathrm{d} x \mathrm{~d} y \mathrm{~d} z [ρq˙−(∂x∂q˙x+∂y∂q˙y+∂z∂q˙z)]dx dy dz
根据傅里叶热传导定律,热传导产生的热流与当地的温度梯度成正比,其中k为热导率。
q ˙ x = − k ∂ T ∂ x q ˙ y = − k ∂ T ∂ y q ˙ z = − k ∂ T ∂ z \dot{q}_{x}=-k \frac{\partial T}{\partial x} \quad \dot{q}_{y}=-k \frac{\partial T}{\partial y} \quad \dot{q}_{z}=-k \frac{\partial T}{\partial z} q˙x=−k∂x∂Tq˙y=−k∂y∂Tq˙z=−k∂z∂T
所以,总的能量公式可进一步推导为:
B
=
[
ρ
q
˙
+
∂
∂
x
(
k
∂
T
∂
x
)
+
∂
∂
y
(
k
∂
T
∂
y
)
+
∂
∂
z
(
k
∂
T
∂
z
)
]
d
x
d
y
d
z
B=\left[\rho \dot{q}+\frac{\partial}{\partial x}\left(k \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial y}\left(k \frac{\partial T}{\partial y}\right)+\frac{\partial}{\partial z}\left(k \frac{\partial T}{\partial z}\right)\right] \mathrm{d} x \mathrm{~d} y \mathrm{~d} z
B=[ρq˙+∂x∂(k∂x∂T)+∂y∂(k∂y∂T)+∂z∂(k∂z∂T)]dx dy dz
最后讨论流体微团能量变化的时间变化率。根据热力学第一定律中内能的物理意义:气体系统的内能就是系统内每个分子和原子能量的总和。因此,运动流体微团的能量主要有两个来源:
1) 由于分子随机运动而产生的单位质量内能 e e e;
2) 流体微团平动时具有的动能。单位质量的动量为 V 2 / 2 V^2/2 V2/2文章来源:https://www.toymoban.com/news/detail-477560.html
因此,运动着的流体微团既有动能又有内能,两者之和就是总能量。由于是跟随着一个运动的流体微团,单位质量的总能量变化的时间变化率由物质导数给出。微团的质量为
ρ
d
x
d
y
d
z
\rho dxdydz
ρdxdydz,所以有:
A
=
ρ
D
D
t
(
e
+
V
2
2
)
d
x
d
y
d
z
A=\rho \frac{\mathrm{D}}{\mathrm{D} t}\left(e+\frac{V^{2}}{2}\right) \mathrm{d} x \mathrm{~d} y \mathrm{~d} z
A=ρDtD(e+2V2)dx dy dz
综上所述,根据能量守恒定律,总的计算公式为:
ρ
D
D
t
(
e
+
V
2
2
)
=
ρ
q
˙
+
∂
∂
x
(
k
∂
T
∂
x
)
+
∂
∂
y
(
k
∂
T
∂
y
)
+
∂
∂
z
(
k
∂
T
∂
z
)
−
∂
(
u
p
)
∂
x
−
∂
(
v
p
)
∂
y
−
∂
(
w
p
)
∂
z
+
∂
(
u
τ
x
x
)
∂
x
+
∂
(
u
τ
y
x
)
∂
y
+
∂
(
u
τ
z
x
)
∂
z
+
∂
(
v
τ
x
y
)
∂
x
+
∂
(
v
τ
y
y
)
∂
y
+
∂
(
v
τ
z
y
)
∂
z
+
∂
(
w
τ
x
z
)
∂
x
+
∂
(
w
τ
y
z
)
∂
y
+
∂
(
w
τ
z
z
)
∂
z
+
ρ
f
⋅
V
\begin{aligned} \rho \frac{\mathrm{D}}{\mathrm{D} t}\left(e+\frac{V^{2}}{2}\right)=& \rho \dot{q}+\frac{\partial}{\partial x}\left(k \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial y}\left(k \frac{\partial T}{\partial y}\right)+\frac{\partial}{\partial z}\left(k \frac{\partial T}{\partial z}\right)-\\ & \frac{\partial(u p)}{\partial x}-\frac{\partial(v p)}{\partial y}-\frac{\partial(w p)}{\partial z}+\frac{\partial\left(u \tau_{x x}\right)}{\partial x}+\frac{\partial\left(u \tau_{y x}\right)}{\partial y}+\frac{\partial\left(u \tau_{z x}\right)}{\partial z}+\\ & \frac{\partial\left(v \tau_{x y}\right)}{\partial x}+\frac{\partial\left(v \tau_{y y}\right)}{\partial y}+\frac{\partial\left(v \tau_{z y}\right)}{\partial z}+\frac{\partial\left(w \tau_{x z}\right)}{\partial x}+\frac{\partial\left(w \tau_{y z}\right)}{\partial y}+\frac{\partial\left(w \tau_{z z}\right)}{\partial z}+\rho \boldsymbol{f} \cdot \boldsymbol{V} \end{aligned}
ρDtD(e+2V2)=ρq˙+∂x∂(k∂x∂T)+∂y∂(k∂y∂T)+∂z∂(k∂z∂T)−∂x∂(up)−∂y∂(vp)−∂z∂(wp)+∂x∂(uτxx)+∂y∂(uτyx)+∂z∂(uτzx)+∂x∂(vτxy)+∂y∂(vτyy)+∂z∂(vτzy)+∂x∂(wτxz)+∂y∂(wτyz)+∂z∂(wτzz)+ρf⋅V
需要注意的是,上述能量方程是非守恒的。文章来源地址https://www.toymoban.com/news/detail-477560.html
到了这里,关于流体力学控制方程——能量方程的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!