高等流体力学 第二章 流体力学基础

第二章 流体力学基础

2.1 连续介质假设

流体力学基本建立在最基本的假设之上——连续介质假设

假设内容：

1. 真实流体由流体质点连续无间隙充满着
2. 每个空间点上、每个时刻都有确定的宏观物理量
3. 这些物理量均为空间，时间的连续函数，除了个别的点线面外

*流体质点是微观无穷大，宏观无穷小的流体单元

使用无量纲量判断：

可使用条件：

$$K_n=\frac{\lambda }{L}<0.001$$

$\lambda$是平均自由程，$L$是所研究对象的特征尺寸。

2.2 流体的性质和分类

流体的性质：

1. 易流动性
2. 粘性
3. 可压缩性

流体的分类：

1. 牛顿流体、非牛顿流体
2. 理想流体、粘性流体
3. 可压缩流体、不可压缩流体

2.3 流体运动的两种描述

2.3.1 拉格朗日法

$$\vec{r}(a,b,c,t)=0$$

2.3.2 欧拉法

定常/非定常是针对欧拉表述下的

$$\vec{r}(x,y,z,t)=0$$

2.4 流体运动的可视化

2.4.1 迹线

$$\frac{\mathrm{d}x}{u}=\frac{\mathrm{d}y}{v}=\frac{\mathrm{d}z}{w}=\mathrm{d}t$$

$t$是自变量。

2.4.2 流线

$$\frac{\mathrm{d}x}{u}=\frac{\mathrm{d}y}{v}=\frac{\mathrm{d}z}{w}$$

$t$是参数，积分时当作常数。

2.4.3 脉线

先求迹线方程，再使用变量替换$（x_0,y_0,z_0)替换(x,y,z),s替换t$

2.4.4 时间线

每一个时刻拍下流场中烟线的变化情况称为时间线。

2.5 应力张量，速度梯度，应变率张量，旋转张量

速度梯度的概念，矢量的梯度是一个二阶张量

$$D_{ij}=\left[\begin{array}{lll}\frac{\partial u}{\partial x}& \frac{\partial u}{\partial y}& \frac{\partial u}{\partial z}\\ \frac{\partial v}{\partial x}& \frac{\partial v}{\partial y}& \frac{\partial v}{\partial z}\\ \frac{\partial w}{\partial x}& \frac{\partial w}{\partial y}& \frac{\partial w}{\partial z}\end{array}\right]={\left[\begin{array}{l}\nabla \vec{v}\end{array}\right]}^T$$

这里写成其转置的形式。

任何一个量都可以分解为对称形式和反对称形似:

$$D_{ij}=\frac{\partial u_i}{\partial x_j}=\frac{1}{2}\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right)+\frac{1}{2}\left(\frac{\partial u_i}{\partial x_j}-\frac{\partial u_j}{\partial x_i}\right)$$

前一项记为$S_{ij}$，为应变率张量，表示剪切变形，第二项记为${\Omega }_{ij}$，为旋转张量，表示旋转变形。

2.6 亥姆霍兹速度分解定理

任一点的速度可以展开求得其他邻域内的点的速度，称为亥姆霍兹速度分解定理。

$$\begin{array}{rl}\vec{v}q& =\vec{v}p+\delta \vec{v}=\vec{v}p+\delta x\frac{\partial \vec{v}}{\partial x}+\delta y\frac{\partial \vec{v}}{\partial y}+\delta z\frac{\partial \vec{v}}{\partial z}\\ & =\vec{v}p+\frac{1}{2}(\nabla \times \vec{v})\times \delta \vec{r}+\mathbf{S}\cdot \delta \vec{r}\end{array}$$

2.7 漩涡的基本概念

2.7.1 漩涡运动的基本概念

定义流体的速度的旋度为流场的涡量$\omega =\nabla \times \mathbf{v}$

$\omega$也是一个矢量场，称为涡量场，涡量场可以给出类似流场的表述形式，如涡线，类比流线：

$$\frac{dx}{{\omega }_x}=\frac{dy}{{\omega }_y}=\frac{dz}{{\omega }_z}$$

涡面，涡管都可以类比流面，流管的概念。

涡通量 在流场中，涡量在某一曲面上的积分称为涡通量。

$$J=\int \int \mathbf{\omega }\cdot \mathrm{d}\mathbf{A}$$

涡管强度 取涡管的一个横截面上的涡通量。

速度环量 对流场中某时刻的封闭曲线做线积分

$$\Gamma ={\oint }_L\mathbf{v}\cdot \mathrm{d}\mathbf{l}$$

由斯托克斯定理可得

$$\Gamma ={\oint }_L\mathbf{v}\cdot \mathrm{d}\mathbf{l}=\int \int \mathbf{\omega }\cdot \mathrm{d}\mathbf{A}$$

微分形式

$$\frac{\mathrm{d}\Gamma }{\mathrm{d}\mathrm{A}}={\omega }_n$$

2.7.2 漩涡运动的基本性质

对欧拉方程取旋得到

$$\frac{\mathrm{d}\mathbf{\omega }}{\mathrm{d}t}-(\mathbf{\omega }\cdot \nabla )\mathbf{v}+\mathbf{\omega }(\nabla \cdot \mathbf{v})=\nabla \times \mathbf{f}+\frac{1}{{\rho }^2}\nabla \rho \times \nabla p$$

这就是无粘流体涡量必须满足的方程

若体力有势，且流体正压，则可以写成

$$\frac{\mathrm{d}\mathbf{\omega }}{\mathrm{d}t}-(\mathbf{\omega }\cdot \nabla )\mathbf{v}+\mathbf{\omega }(\nabla \cdot \mathbf{v})=0$$

由此可以得到开尔文定律——正压无粘流体在体力有势时，沿着封闭流体线的速度环量在流动过程中保持不变

$$\frac{\mathrm{d}\Gamma }{\mathrm{d}t}={\oint }_L\frac{\mathbf{v}}{\mathrm{d}t}\cdot \mathrm{d}\mathbf{l}=0$$

由此条性质可以得到Lagrange定理：

在体力有势、流体正压的条件下，无粘流体若在某一时刻无旋，则该流体运动就永远无旋。漩涡不会产生，也不会消失。

和亥姆霍兹涡线、涡面、涡管保持定理（亥姆霍兹第一定理）：

在某一时刻组成涡面、涡管的流体质点在其他时刻也仍在组成涡面和涡管。

$$\frac{\mathrm{d}\Gamma }{\mathrm{d}t}=\frac{\mathrm{d}J}{\mathrm{d}t}$$

涡管强度保持不变时亥姆霍兹第二定理，表面涡管只能在流场的边界开始或结束。

2.8 势函数和流函数

1 势函数

当流体无旋时，满足 $\nabla \times \mathbf{v}=0$

由张量分析得知 梯度的散度等于0，即 $\nabla \times (\nabla \phi )=0$

因此，流体无旋时，速度有势。

$$u=\frac{\partial \phi }{\partial x}，v=\frac{\partial \phi }{\partial y},w=\frac{\partial \phi }{\partial z}$$

且势函数是调和函数，流动为不可压缩时，满足Laplace方程

$$\Delta \phi =0$$

2 流函数

流函数只适用于二维，不可压缩。

天然满足不可压缩的条件

$$\nabla \cdot \mathrm{v}=0$$

$$u=\frac{\partial \psi }{\partial y}、v=-\frac{\partial \psi }{\partial x}$$

同时当流体无旋时，满足Laplace方程

$$\Delta \psi =0$$

当流体有旋时，不满足Laplace方程

$$\Delta \psi =-\omega$$

2.9 应力张量、主应力、主平面、主轴

某一个面的应力张量的求法

$${\vec{T}}_i={\sigma }_{ij}{n}_j$$

求特征值就是求主应力

$$\begin{array}{l}\vec{T}(\vec{n})=\sigma \vec{n}\\ {\sigma }_{ij}{n}_j=\sigma {\vec{n}}_i\\ ({\sigma }_{ij}-\sigma {\delta }_{ij})n_j=0\\ \left|{\sigma }_{ij}-\sigma {\delta }_{ij}\right|=0\end{array}$$

2.10 （雷诺）输运定理

$$\begin{array}{rl}\frac{\mathrm{d}}{\mathrm{d}t}{\int }_{V}L(\vec{x},t)\mathrm{d}V(\vec{x})& =\frac{\mathrm{d}}{\mathrm{d}t}{\int }_{V_0}L(\vec{x}(\vec{a}),t)J\mathrm{d}V(\vec{a})\\ & ={\int }_{V_0}\frac{\mathrm{d}}{\mathrm{d}t}[L(\vec{x}(\vec{a}),t)J]\mathrm{d}V(\vec{a})\\ & ={\int }_{V_0}\left(J\frac{\mathrm{d}}{\mathrm{d}t}L+L\frac{\mathrm{d}J}{\mathrm{d}t}\right)\mathrm{d}V(\vec{a})\\ & ={\int }_{V_0}\left(J\frac{\mathrm{d}}{\mathrm{d}t}L+LJ\nabla \cdot \vec{v}\right)\mathrm{d}V(\vec{a})\\ & ={\int }_V\left(\frac{\mathrm{d}L}{\mathrm{d}t}+L\nabla \cdot \vec{v}\right)\mathrm{d}V(\vec{x})\\ & ={\int }_V\left[\frac{\partial L}{\partial t}+\nabla \cdot (L\vec{v})\right]\mathrm{d}V(\vec{x})\end{array}$$

连续方程、动量方程、动能方程可以统一左边的格式。

2.11 流体运动的控制方程

1. 连续方程 $L(\vec{x},t)=\rho$
2. 动量方程 $L(\vec{x},t)=\rho v$
3. 能量方程 $L(\vec{x},t)=(\frac{1}{2}\rho v^2+\rho \epsilon )$

加上状态方程、本构关系、热力学方程即可组成完备的流体力学控制方程

$$\begin{array}{l}{\partial }_t\rho +\nabla \cdot (\rho \vec{v})=0\\ {\partial }_t{v}_i+v_j{\partial }_j{v}_i=f_i+{\partial }_j{\sigma }_{ij}/\rho \\ {\partial }_t\epsilon +v_j{\partial }_j\epsilon =\left({\sigma }_{ij}{\partial }_j{u}_i-{\partial }_j{q}_j\right)/\rho +Q\\ f(p,\rho ,T)=0\\ {\sigma }_{ij}=-p{\delta }_{ij}+2\mu S_{ij}+\lambda S_{kk}{\delta }_{ij}\\ \epsilon \sim T\\ \vec{q}\sim \nabla T\end{array}$$