第一章 计算流体力学动力学基础知识

---

style: candy

第一章 计算流体力学动力学基础知识

CFD基本思想：把原来在时间域及空间域上连续的物理量场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后通过解代数方程组获得场变量的近似解。
计算流体力学：通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。

CFD可以看做是在流动基本方程（质量守恒方程、动量守恒方程、能量守但方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（如速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等

计算流体动力学的工作步骤

采用CFD的方法对流体流动进行数值楼拟，通常包括如下步骤：
(1) 建立反映工程问题或物理问题本质的数学模型，，具体地说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程，动量守恒方程、能量守恒方程以及这些方程相应的定解条件，

(2) 寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法．如有限差分法、有限元法、有限体积法等, 这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是CFD 的核心。

(3) 编制程序和进行计算这部分工作包括计算网格划分、初始条件和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Navier-Stokes方程就是一个十分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲，数值模拟又叫数值试验。应该指出，这部分工作不是轻而易举就可以完成的。

(4) 显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。

有限差分法

将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出方程组的解，就是微分方程定解问题的数值近似解。

粘性

粘性(viscocity)是流体内部发生相对运动而引起的内部柜互作用

流体在静止时虽不能承受切应力，但在运动时，对相邻两层流体间的相对运动，即相对滑动速度却是有抵抗的，这种抵抗力称为粘性应力，流体所具有的这种抵抗两层流体间相对滑动速度，或普遍说来抵抗变形的性质，称为粘性.

注：计算流体力学的复杂性就在于粘性的复杂

流体热传导及扩散

除了粘性外，流体还有热传导及扩散(diffusion)等性质，当流体中存在着温度差时，温度高的地方将向温度低的地方传送热量．这种现象称为热传导。同样地，当流体混合物中存在着组元的浓度差时，浓度离的地方将向浓度低的地方输送该组元的物质，这种现象称为扩散。(温度不是产生热传导的唯一因素)

热传导定义：
热传导（thermal conduction）是介质内无宏观运动时依靠分子、原子及自由电子等微观粒子热运动而进行的传热传递现象，其在固体、液体和气体中均可发生，但严格而言，只有在固体中才是纯粹的热传导，而流体即使处于静止状态，其中也会由于温度梯度所造成的密度差而产生自然对流，因此，在流体中热对流与热传导同时发生。

流体的宏观性质，如扩散、粘性和热传导等，是分子输运性质的统计平均。由于分子的不规则运动，在各层流体间交换着质量、动量和能量，使不同流体层内的平均物理量均匀化。这种性质称为分子运动的输运性质。质量输运在宏观上表现为扩散现象，动量输运表现为粘性现象，能量输运则表现为热传导现象。

层流与湍流

自然界中的流体流动状态主要有两种形式，即层流(laminar)和湍流(turbulence) 。在许多中文文献中，湍流也被译为紊流。层流是指流体在流动过程中两层之间没有相互混掺，而湍流是指流体不是处于分层流动状态。—般说来，湍流是普遍的，而层流则属千个别情况。

流体动力学控制方程

流体流动要受物埋守恒定律的支配，基本的守恒定律包括：质最守恒定律、动量守恒定律、能量守恒定律。如果流动包含有不同成分（组元）的混合或相互作用．系统还要遵守组分守恒定律。如果流动处于湍流状态，系统还耍遵守附加的湍流输运方程。
控制方程(governing equations)是这些守恒定律的数学描述。本节先介绍这些基本的守恒定律所对应的控制方程．有关湍流的附加控制方程将在第4 章中介绍。

质量守恒方程

该定律可表述为：单位时间内流体微元体中质量的增加，等于同一时间间隔内流入该微元体的净质量。质量守恒方程也叫连续方程
$$\frac{\partial \rho }{\partial t}+\frac{\partial (\rho u)}{\partial x}+\frac{\partial (\rho v)}{\partial y}+\frac{\partial (\rho w)}{\partial z}=0$$
引入$div(\mathbf{a})=\partial a_x/\partial x+\partial a_y/\partial y+\partial a_z/\partial z$，用$\nabla$表示散度
KaTeX parse error: No such environment: gather at position 8: \begin{̲g̲a̲t̲h̲e̲r̲}̲ \frac{\partial…
$u$为速度矢量，$u$,$v$,$w$是速度矢量$\mathbf{u}$在$x、y、z$方向的分量

对于不可压流体：
$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}=0$$

动量守恒方程

$$\begin{gathered} \frac{\partial (\rho u)}{\partial t}+div(\rho u\mathbf{u})=-\frac{\partial p}{\partial x}+\frac{\partial {\tau }_{xx}}{\partial x}+\frac{\partial {\tau }_{yx}}{\partial y}+\frac{\partial {\tau }_{zx}}{\partial z}+F_x \\ \frac{\partial (\rho v)}{\partial t}+div(\rho v\mathbf{u})=-\frac{\partial p}{\partial y}+\frac{\partial {\tau }_{xy}}{\partial x}+ \end{gathered}$$