计算流体力学完整
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3
计算流体力学(CFD)：通过数值方法求解流体力学控制 方程，得到流场的离散的定量描述，并以此预测流体运 动规律的学科。
在CFD中， 首先，把控制方程中的积分、微分项近似地表示为离散的代数形 式，把积分、微分形式的控制方程转化为一组代数方程，这个过 程称为控制方程的离散化(discretization)；所采用的离散化方法 称为数值方法或数值格式。
The Elements of Computational Fluid Dynamics
1
第一章 绪论
§1.1 计算流体力学的概念与意义 §1.2 流体力学的基本方程 §1.3 流体力学方程组的类型判别
2
§1.1 计算流体力学的概念与意义
1、流体运动遵循3个基本定律： 1) 质量守恒定律；2) 动量守恒定律；3) 能量守恒定律
6
第六，数值解的显示和评估
计算感兴趣的力、力矩等； 应用流场可视化软件对流场进行显示、分析； 对数值方法和物理模型的误差进行评估等。
7
计算流体力学典型流程
物
数
网
理
学
格
模
模
生
型
型
成
结
流
果
场
分
显
析
示
验 证 与 确 认
离 散 方 法 选 择
时、空离散
解 代 边界条件离散 数 方 程 组
8
举例：自然循环回路内的流动与传热特性
优点：原则上可以研究流体在任何条件下的运动，使得我们研究流体运动的范围和 能力都有本质的扩大和提高。费用低，周期短。
16
§1.2 流体力学基本方程
守恒型积分方程
t
d
Ò V
s
ds
0
t
Vd
乙 VV
s
ds
Fd
s
t
*
然后，通过电子计算机求解这些代数方程组，得到流场在离散的 时间/空间点上的数值解(numerical solution)。
CFD也被称作流场的数值模拟、数值计算、数值仿真等。
4
计算流体力学的研究步骤
第一，问题的界定和流动区域的几何描述。
流场的几何形状：源于对已有流动区域的测量或者新的产品和 工程的设计结果。
第四，程序设计和调试。
在网格划分策略和数值方法的基础上，编制、调试数值求解流体运动方程 的计算机程序或软件。
第五，程序验证和确认。
验证(Verification)：The process of determining that a model implementation accurately
represents the developer’s conceptual description of the model and the solution to the
10
数学模型：
控制方程
r
(V
)
0
t
r
V
rr
(VV )
r
P [( V )]
gr
t
E
r
(VE)
(T
)
T
P
q
t
t
定解条件
初始条件：
r V
0;
T T0
边界条件：固体壁面上无滑移；
恒温热源、恒温热沉，
其余为绝热壁面。
11
网格划分：
12
数值算法：
离散方法： FDM、FVM、FEM……
空间离散： 对流项，粘性项，源项……
时间离散： 显式、隐式
边界离散： 来流、出流、固壁、远场、周期性……
求解代数方程组
13
数值解的验证与确认：
14
流场显示及结果分析：
15
计算流体力学的特点及意义
实验研究
优点：借助各种先进仪器，给出多种复杂流动的准确、可靠的观测结果，这些结果 对于流动机理的研究和与流体运动有关的机械和飞行器的设计具有不可替代的作用。
边界条件通常依赖于控制方程。
固体壁面条件，来流、出流条件，周期性条件，对称条件等
附加的物理模型：湍流模型，化学反应等。
5
第三，确定网格划分策略和数值方法。
网格划分：结构网格、非结构网格、组合网格、重叠网格。 网格可以是静止的，也可以是运动的，还可以根据数值解动态调整
（自 适应网格）。
数值方法：有限差分、有限体积、有限元、谱方法等。 数值方法和网格划分策略是相互关联的。
2、流体的本构模型和状态方程
控制方程（Governing equations）
偏微分方程（方程组）或积分形式的方程（方程组）
流体运动的复杂性主要表现为控制方程的高度非线性和流动 区域几何形状的复杂性等，导致对绝大多数流动问题无法得 到解析解。 高速计算机的发展，使得计算流体动力学（Computational Fluid Dynamics,CFD）逐渐成为一门独立学科。
9
物理模型：
(1) 空间维数：1D、2D、3D (2) 时间特性：定常、非定常 (3) 流动性质：无粘/粘性、可压缩/不可压缩、层流/湍流 (4) 流体物性：常物性、变物性
Geometric parameter: Height H Width W Length of heat sink (source) L Tube diameter d Rayleigh number Ra Heat source temperature Th Heat sink temperature Tc Operation pressure P
缺点：费用高昂，周期很长，有些流动条件难以通过实验手段来模拟。
理论研究
优点：可以给出具有一定适用范围的简洁明了的解析解或近似解析解，这些解析解对 于分析流动的机理和预测流动随参数的变化非常有用。
缺点：只能研究简单流动问题，能够得到解析解的流动问题为数不多，远远不能满足 工程设计的需要。
计算流体力学
发展CFD的主要动机：利用高速电子计算机，克服理论研究和实验研究的缺点，深 化对于流体运动规律的认识并提高解决工程实际问题的能力。
流动条件：雷诺数、马赫数、边界处的速度及压力等 对数值模拟的要求：精度、所花费的时间。
第二，选择控制方程和边界条件。
在牛顿流体范围内，用Navier-Stokes方程描述。 根据问题的特点，可以考虑定常或非定常，可压或不可压的流动模型。
简化的数学模型：势流方程，Euler方程，边界层方程， 薄层近似的Navier-Stokes方程等。
nds
t
Ed
乙 Βιβλιοθήκη sEVds
F
Vd
s
(t
*
V
q)
nds
牛顿流体本构关系 ti*j pij
计算流体力学(CFD)：通过数值方法求解流体力学控制 方程，得到流场的离散的定量描述，并以此预测流体运 动规律的学科。
在CFD中， 首先，把控制方程中的积分、微分项近似地表示为离散的代数形 式，把积分、微分形式的控制方程转化为一组代数方程，这个过 程称为控制方程的离散化(discretization)；所采用的离散化方法 称为数值方法或数值格式。
The Elements of Computational Fluid Dynamics
1
第一章 绪论
§1.1 计算流体力学的概念与意义 §1.2 流体力学的基本方程 §1.3 流体力学方程组的类型判别
2
§1.1 计算流体力学的概念与意义
1、流体运动遵循3个基本定律： 1) 质量守恒定律；2) 动量守恒定律；3) 能量守恒定律
6
第六，数值解的显示和评估
计算感兴趣的力、力矩等； 应用流场可视化软件对流场进行显示、分析； 对数值方法和物理模型的误差进行评估等。
7
计算流体力学典型流程
物
数
网
理
学
格
模
模
生
型
型
成
结
流
果
场
分
显
析
示
验 证 与 确 认
离 散 方 法 选 择
时、空离散
解 代 边界条件离散 数 方 程 组
8
举例：自然循环回路内的流动与传热特性
优点：原则上可以研究流体在任何条件下的运动，使得我们研究流体运动的范围和 能力都有本质的扩大和提高。费用低，周期短。
16
§1.2 流体力学基本方程
守恒型积分方程
t
d
Ò V
s
ds
0
t
Vd
乙 VV
s
ds
Fd
s
t
*
然后，通过电子计算机求解这些代数方程组，得到流场在离散的 时间/空间点上的数值解(numerical solution)。
CFD也被称作流场的数值模拟、数值计算、数值仿真等。
4
计算流体力学的研究步骤
第一，问题的界定和流动区域的几何描述。
流场的几何形状：源于对已有流动区域的测量或者新的产品和 工程的设计结果。
第四，程序设计和调试。
在网格划分策略和数值方法的基础上，编制、调试数值求解流体运动方程 的计算机程序或软件。
第五，程序验证和确认。
验证(Verification)：The process of determining that a model implementation accurately
represents the developer’s conceptual description of the model and the solution to the
10
数学模型：
控制方程
r
(V
)
0
t
r
V
rr
(VV )
r
P [( V )]
gr
t
E
r
(VE)
(T
)
T
P
q
t
t
定解条件
初始条件：
r V
0;
T T0
边界条件：固体壁面上无滑移；
恒温热源、恒温热沉，
其余为绝热壁面。
11
网格划分：
12
数值算法：
离散方法： FDM、FVM、FEM……
空间离散： 对流项，粘性项，源项……
时间离散： 显式、隐式
边界离散： 来流、出流、固壁、远场、周期性……
求解代数方程组
13
数值解的验证与确认：
14
流场显示及结果分析：
15
计算流体力学的特点及意义
实验研究
优点：借助各种先进仪器，给出多种复杂流动的准确、可靠的观测结果，这些结果 对于流动机理的研究和与流体运动有关的机械和飞行器的设计具有不可替代的作用。
边界条件通常依赖于控制方程。
固体壁面条件，来流、出流条件，周期性条件，对称条件等
附加的物理模型：湍流模型，化学反应等。
5
第三，确定网格划分策略和数值方法。
网格划分：结构网格、非结构网格、组合网格、重叠网格。 网格可以是静止的，也可以是运动的，还可以根据数值解动态调整
（自 适应网格）。
数值方法：有限差分、有限体积、有限元、谱方法等。 数值方法和网格划分策略是相互关联的。
2、流体的本构模型和状态方程
控制方程（Governing equations）
偏微分方程（方程组）或积分形式的方程（方程组）
流体运动的复杂性主要表现为控制方程的高度非线性和流动 区域几何形状的复杂性等，导致对绝大多数流动问题无法得 到解析解。 高速计算机的发展，使得计算流体动力学（Computational Fluid Dynamics,CFD）逐渐成为一门独立学科。
9
物理模型：
(1) 空间维数：1D、2D、3D (2) 时间特性：定常、非定常 (3) 流动性质：无粘/粘性、可压缩/不可压缩、层流/湍流 (4) 流体物性：常物性、变物性
Geometric parameter: Height H Width W Length of heat sink (source) L Tube diameter d Rayleigh number Ra Heat source temperature Th Heat sink temperature Tc Operation pressure P
缺点：费用高昂，周期很长，有些流动条件难以通过实验手段来模拟。
理论研究
优点：可以给出具有一定适用范围的简洁明了的解析解或近似解析解，这些解析解对 于分析流动的机理和预测流动随参数的变化非常有用。
缺点：只能研究简单流动问题，能够得到解析解的流动问题为数不多，远远不能满足 工程设计的需要。
计算流体力学
发展CFD的主要动机：利用高速电子计算机，克服理论研究和实验研究的缺点，深 化对于流体运动规律的认识并提高解决工程实际问题的能力。
流动条件：雷诺数、马赫数、边界处的速度及压力等 对数值模拟的要求：精度、所花费的时间。
第二，选择控制方程和边界条件。
在牛顿流体范围内，用Navier-Stokes方程描述。 根据问题的特点，可以考虑定常或非定常，可压或不可压的流动模型。
简化的数学模型：势流方程，Euler方程，边界层方程， 薄层近似的Navier-Stokes方程等。
nds
t
Ed
乙 Βιβλιοθήκη sEVds
F
Vd
s
(t
*
V
q)
nds
牛顿流体本构关系 ti*j pij