高等流体力学第三章不可压粘性流体运动

夏天，我来了

已于 2022-03-01 22:06:47 修改

阅读量1.4k

点赞数 1

分类专栏：高等流体力学文章标签：流体力学学习经验分享

于 2022-03-01 22:05:18 首次发布

本文链接：https://blog.csdn.net/weixin_41474955/article/details/115142565

版权

高等流体力学专栏收录该内容

5 篇文章 4 订阅

订阅专栏

第三章不可压粘性流体运动

从我的笔记搬运过来的，有些格式图片没办法直接搬过来，因此后续内容可以直接链接到我的Notion笔记中
第三章不可压粘性流体运动

上一章中给出了流体力学的控制方程，可以看到方程组是非常复杂的，因此在实际应用中我们往往进行一些简化，不可压缩是其中最常用的一种方法，通过不可压缩的条件，可以得到最经典的N-S方程的形式。

3.1 控制方程

不可压缩，连续方程变为：

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

动量方程转变为：

$\frac{{\partial \vec v}}{{\partial t}} + (\vec v \cdot \nabla )\vec v = \vec f - \frac{1}{\rho }\nabla p + \frac{\mu }{\rho }\Delta \vec v$

左端第一项是非定常加速度，第二项是对流加速度，右端第一项是体积力，第二项是压强梯度，第三项是粘性项。

此时能量方程结合傅里叶定律，不可压，均匀流，方程可以转变为

$\frac{1}{2} \rho \frac{\partial}{\partial t} \int_{V} v^{2} d V=\rho \int_{V} v_{i} \frac{\partial v_{i}}{\partial t} d V$

状态方程

$f(p,\rho ;T) = 0$

3.2 确定N-S方程解的条件

微分方程得到唯一定解还需要结合边界条件和初始条件，因此本节主要给出常见的几种边界条件。

固体边界

固体边界要满足不穿透、无滑移的条件

$无穿透\vec v \cdot \vec n = \vec U \cdot \vec n$

$无滑移\begin{array}{l}{{\vec v}_f} = {{\vec U}_s}\\{T_f} = {T_s}\end{array}$

*对于多孔界面，存在质量交换，则不满足无穿透条件。

2. 自由表面

运动边界条件：

$t)=\frac{\partial \zeta}{\partial t}+u \frac{\partial \zeta}{\partial x}+v \frac{\partial \zeta}{\partial y}$

动力边界条件：

无表面张力： $p_1=p_2$

有表面张力： $p_1-p_2=\sigma (\frac{1}{R_1}+\frac{1}{R_2})$

能量条件： ${\left. {k\frac{{\partial T}}{{\partial n}}} \right|_1} = {\left. {k\frac{{\partial T}}{{\partial n}}} \right|_2}$

*不考虑交界面不连续的情况

对于无对流项、均匀流、体积力有势的情况下

N-S方程如下：

$\left\{\begin{array}{l} \nabla \cdot \vec{v}=0 \\ \rho \frac{\partial \vec{v}}{\partial t}=\rho \vec{f}-\nabla p+\mu \nabla^{2} \vec{v} \\ \rho c_{v} \frac{\partial T}{\partial t}=k \nabla^{2} T+\Phi \end{array}\right.$

对二式取散度，得到

${\nabla ^2}p = 0$

这是一个椭圆方程。

若取散度，得到

$\frac{{\partial \vec \omega }}{{\partial t}} = \nu \Delta \vec \omega$

这是一个抛物线方程。

对于定常欧拉方程，马赫数不同，方程类型不同，所需要的边界条件和初值条件也是不同的。

在这里插入图片描述

3.3 无量纲数（相似参数）

无量纲化控制方程：

方程中本含有 $\rho, T_{0} ,T_{0}=T_{\infty} \text { or } T_{w}-T_{\infty}$ ，对这些量进行无量纲化，需要使用到特征长度 $L$ ，特征速度 $U$ ，特征密度 $\rho$

$\begin{aligned} x^{*}=\frac{x}{L}, y^{*} &=\frac{y}{L}, z^{*}=\frac{z}{L}, t^{*}= \frac{t}{L / U} \\ &\vec{v}^{*}=\frac{\vec{v}}{U}, p^{*}=\frac{p+\rho g z}{\rho U^{2}}, \\ &T^{*}=\frac{T}{T_{0}} \text { or } \frac{T-T_{\infty}}{T_{0}}\end{aligned}$

由此变形代入回原方程中，便可以得到无量纲的方程形式：

连续方程

${\nabla ^*} \cdot {{\vec v}^*} = 0$

动量方程

$\frac{{\partial {{\vec v}^*}}}{{\partial {t^*}}} + ({{\vec v}^*} \cdot {\nabla ^*}){{\vec v}^*} = - {\nabla ^*}{p^*} + \frac{\mu }{{\rho UL}}{\Delta ^*}{{\vec v}^*}$

能量方程

$\frac{{\partial {T^*}}}{{\partial {t^*}}} + ({{\vec v}^*} \cdot {\nabla ^*}){T^*} = \frac{k}{{\rho UL{c_v}}}{\Delta ^*}{{\vec v}^*} + \frac{{\mu U}}{{\rho {c_v}L{T_0}}}(2\overrightarrow {\vec S} :\overrightarrow {\vec S} )$

由此得到了一些常用的无量纲数：

Reynolds number： ${\mathop{\rm Re}\nolimits} = \frac{{\rho UL}}{\mu }$
Prandtl number： $\Pr = \frac{\mu }{{{c_p}k}}$
Eckert number： $\frac{{{U^2}}}{{{c_p}{T_0}}}$
比热： $\gamma = {c_p}/{c_v}$

方程可以进一步简化：

这里少了重力一项，还有Fr弗劳德数
在这里插入图片描述

无量纲化边界条件：

在这里插入图片描述

其他无量纲数：

在这里插入图片描述

在所有无量纲数中，最重要的就是雷诺数

雷诺数

雷诺数从小到大，圆柱绕流会出现不同的尾部结构。
雷诺数和阻力的关系，在小雷诺数下成线性关系，大雷诺数下成二次关系，与阻力系数的关系先是下降再是水平。不同结构对阻力系数影响也很大，影响流场何时与物理表面发生脱落。
雷诺数和弗劳德数的关系，当存在自由表面时，需要关注弗劳德数，因为会产生兴波阻力，船鼻首就是为了减小兴波阻力。
雷诺数和马赫数关系，马赫数低于0.4时，雷诺数占据主导，马赫数大于0.4时，马赫数占据主导。

相似定理

几何相似

运动相似

动力相似

3.4 N-S方程的精确解

精确解基本上都是省略掉对流项，因此，通常都是层流下的平行流动，低雷诺数流动等。本章主要内容是层流运动。

对流体问题的解答方法论

选择合适的坐标系
写下控制方程
分析流动状态
简化方程
边界条件、初始条件
求解
检验结果

常见的几种能够得到解析解的情况

1. 线性，定常

定常，可以将时间项直接去掉，线性可以将非线性项去掉，从而将动量方程转化为线性偏微分方程

（1）平行平板间的不可压缩运动

控制方程转换为：

$\begin{array}{l} 0 = G + \mu \frac{{{\partial ^2}u}}{{\partial {y^2}}}\\ 0 = k\frac{{{\partial ^2}T}}{{\partial {y^2}}} + \mu {(\frac{{\partial u}}{{\partial y}})^2} \end{array}$