5 计算流体力学基础一
计算流体力学基础及其应用
计算流体力学基础及其应用计算流体力学（CFD）是计算机运用精确的数学模型和算法来研究流体力学物理过程的一种技术。
它利用计算机模拟方法处理流体流动和相互作用的过程，以更准确、更快捷的方式研究热流体流动、传热、传质和湍流等物理过程的问题。
CFD的基础是数学方面的流体力学，应用计算机模拟的基本方法是数值方法，用于分析各种流体流动问题以及相关热传导、传质等热力学现象。
此外，计算流体力学还集成有计算机动力学，流体动力学，热力学，结构力学，能量方法，计算工程和多物理场的数值模拟技术，可以更加精准地研究流体动力学，热传递，流体机械，复杂流动等问题。
CFD在工程实践中具有重要作用，其应用领域非常广泛，包括空气、液体、气体和粘性流动等各种固体表面及流体体系的运动和相互作用。
例如，可以用来分析大气环境中污染物的扩散，水力学中河流水流的流动性能和可能形成的机械，风能资源的开发利用，以及气体控制元件的设计等。
CFD技术的研究和应用对改善工业和生活的质量起着重要作用，具有重大的经济效益。
它可以帮助工程师进行快速和准确的表征及设计，从而大大缩短研发和评估的周期，并节省大量的研发费用，从而提高产品的质量和可靠性。
例如，可以用CFD模拟来分析火力发电厂泄漏物介质的运动和湍流，从而确定阀门及其参数，进行管道设计，抑制烟气污染，提高系统效率，实现节能减排等。
此外，CFD还可以用于水工工程，海洋工程，气候变化，大气和海洋环境监测，飞机设计，汽车行业和其他工程方面的问题，有助于数字信息的可视化，预测及避免工程问题，提高效率。
因此，CFD既可以用于重要的实际问题的研究，也可以用于开发新产品，从而为工程实践提供可靠的计算技术，有效地改善系统质量和可靠性，提高经济效益。
综上所述，CFD的研究和应用具有重要的实际意义，可以显着提高工程的质量和可靠性，并带来可观的经济收益。
未来，CFD技术将逐步发展壮大，有效地改善人们的生活和工作环境。
流体力学与计算流体力学基础
第1章流体力学与计算流体力学基础流体力学是力学的一个重要分支，它主要研究流体本身的静止状态和运动状态，以及流体和固体界壁间有相对运动时的相互作用和流动的规律，在生活、环保、科学技术及工程中具有重要的应用价值。
计算流体力学或计算流体动力学（Computational Fluid Dynamics，CFD），是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。
本章先介绍流体力学中支配流体流动的基本物理定律，然后在此基础上介绍用数值方法求解流体力学问题的基本思想，进而阐述计算流体力学的相关基础知识，最后简要介绍常用的计算流体力学商业软件。
学习目标：•学习流体力学的基础知识，包括基本概念和重要理论；•学习计算流体力学的相关理论和方法；•了解CFD软件的构成；•了解常用的商业CFD软件。
1.1 流体力学基础流体力学是连续介质力学的一个分支，是研究流体（包含气体及液体）现象以及相关力学行为的科学。
1.1.1 流体力学概述1738年，伯努利在他的专著中首次采用了水动力学这个名词并作为书名；1880年前后出现了空气动力学这个名词；1935年以后，人们概括了这两方面的知识，建立了统一的体系，统称为流体力学。
在人们的生活和生产活动中随时随地都可遇到流体，因此流体力学是与人类日常生活和生产事业密切相关的。
大气和水是最常见的两种流体，大气包围着整个地球，地球表面的70%是水面。
大气运动、海水运动（包括波浪、潮汐、中尺度涡旋、环流等）乃至地球深处熔浆的流动都是流体力学的研究内容。
20世纪初，世界上第一架飞机出现以后，飞机和其他各种飞行器得到迅速发展。
20世纪50年代开始的航天飞行，使人类的活动范围扩展到其他星球和银河系。
航空航天事业的蓬勃发展是同流体力学的分支学科——空气动力学和气体动力学的发展紧密相连的。
这些学科是流体力学中最活跃、最富有成果的领域。
石油和天然气的开采、地下水的开发利用，要求人们了解流体在多孔或缝隙介质中的运动，这是流体力学分支之一——渗流力学研究的主要对象。
计算流体基础
X 方向净流出量为：
Y 方向净流出量为：
Z 方向净流出量为：
控制体质量流出量
该控制体的质量为密度 X 体积，即：
则其质量随时间的变化率为：
对于该控制体来说， 质量守恒就是流出该控制体的质量流量等于控制体内的质量 减少量.将质量减少定义为负，则质量守恒可以表达为：
即：
或者写为散度形式：
综上，控制体在 x 方向上的所受表面力为：
同理，可得所有表面力在 x 方向所做功率为：
所有的体积力和表面力所做功功率为：
流入控制体净热流量主要来自于体积加热，比如吸收或释放的辐射热和热传 导等.定义 为单位质量的体积加热率，则控制体的体积加热为：
假设 为热传导在单位时间内通过单位面积在 x 方向上输运的热量，那么对于控 制体在 x 方向上热传导的热量为：
同理可得整个控制体的热流量为：
根据傅里叶热传导定律，热传导产生的热量与当地温度梯度成正比：
所以总热流量可表示为：
Байду номын сангаас
物体的总能量为动能与内能之和，为 e v 2 / 2 则控制体能量变化率为：
则能量方程可表示为：
这五个方程统称为 NS 方程，然后再引入两个补充状态方程：
这样，整个方程组就封闭了，CFD 的任务就是求解这一组方程.
当 x 足够小时，我们可以近似的认为
用线性代数式表示非线性微分方程，这就是求解微分方程的基本思想. 根据泰勒展开公式有：
可以看到最低阶项是 x 的一次方项，我们可以将上式写为：
舍弃 (x) ，则可得到有限差分表达式：
其中舍弃的部分 (x) 称为截断误差.可以看出该表达式的截断误差最低项是 x 的一次方，则称该有限差分格式具有一阶精度，由于该表达式仅仅用到了(i,j)和 (i+1,j)的信息，该格式称为一阶向前差分. 再看(i‐1,j)这个点的泰勒展开式为：
计算流体力学及其并行算法
计算流体力学及其并行算法一、引言计算流体力学（Computational Fluid Dynamics, CFD）是研究流体运动和相互作用的一门学科，广泛应用于工程、天文、地球科学等领域。
随着计算机技术的发展，CFD的数值模拟方法也得到了极大的发展，其中并行算法在加速CFD计算过程中起到了重要的作用。
二、计算流体力学基础1. 流体力学基本方程计算流体力学的基础是流体力学的基本方程，包括质量守恒方程、动量守恒方程和能量守恒方程。
这些方程描述了流体的运动、力学性质和能量转换。
2. 数值离散化方法为了将流体力学方程转化为计算模型，需要对连续域进行离散化。
常用的数值离散化方法包括有限差分法、有限体积法和有限元法等。
这些方法将连续的流体域离散为网格，通过在网格上的节点上进行数值计算，得到流体的各个物理量。
三、并行算法在计算流体力学中的应用1. 并行计算的需求计算流体力学涉及大规模的计算，需要处理大量的数据和复杂的计算操作。
传统的串行计算方式往往难以满足计算需求，因此并行算法成为加速CFD计算的重要手段。
2. 并行算法分类并行算法根据不同的并行计算方式，可以分为共享内存并行和分布式内存并行两大类。
共享内存并行算法使用多个处理器共享同一块内存，通过线程间的数据共享和同步来实现并行计算；分布式内存并行算法则将计算任务分配到不同的处理器上，通过消息传递来实现并行计算。
3. 并行算法的优势并行算法在加速CFD计算中具有显著的优势。
首先，通过并行计算，可以将计算任务分配到多个处理器上，实现计算资源的充分利用。
其次，并行算法可以处理大规模的计算问题，提高计算效率和精度。
此外，并行算法还可以实现实时计算和交互式计算，提供更好的用户体验。
四、并行算法的挑战和发展方向1. 数据通信和负载均衡在并行计算过程中，处理器之间需要进行数据通信，这涉及到数据传输和同步操作。
数据通信的效率和负载均衡是并行算法面临的挑战之一，需要合理设计算法和优化通信过程。
第二章--计算流体力学的基本知识
第二章--计算流体力学的基本知识第二章计算流体力学的基本知识流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。
这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。
2.1计算流体力学简介2.1.1计算流体力学的发展流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。
20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。
数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。
从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。
数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。
数值计算方法最近发展很快，其重要性与日俱增。
自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。
最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。
航空技术的发展强烈推动了流体力学的迅速发展。
*流体运动的规律由一组控制方程描述。
计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。
但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。
计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力学这门交叉学科。
计算流体力学是一门用数值计算方法直接求解流动主控方程(Euler 或Navier-Stokes方程)以发现各种流动现象规律的学科。
计算流体力学入门
并 且 数 值 上 等 于 格 心 处 的 流 场 参 数 值 ， Fi 是 每 个 控 制 面 上 F 的 平 均 值 ， 即 记
u
1 V
1 ， F d u V i C.V Si
u V F 。那相当于求解 F dS i Si 0 。这个方程就 c.si t i
第 i-1 点:
(1.3)
ui 1 (u u ) ai 1 i 1/ 2 i 3/ 2 0 t x
以上三式相加，得：
(1.4)
ui 1 ui ui 1 ai 3/ 2 ai 3/ 2 (ui 3/ 2 ui 3/ 2 ) ( ) t x 3 6 a u ) a (u (u ) a a u i 1/ 2 i 3/ 2 i 3/ 2 i 1/ 2 i 3/ 2 i 1/ 2 i 1/ 2 i 3 / 2 6 6 x x
这里 a 不仅限于一个数，对于最复杂的情形——多维、多个方程，a 是一个 Jacobi 矩阵。 这个变形可以描述为方程可以化为拟线性形式。这里只分析最简单的情况，即 a 是一个数。 对这种问题的分析可以很类似地推广到方程组和多个空间坐标的情形， 比如欧拉方程， 甚至 是完整的 NS 方程，RANS 方程。 对于方程
u u a(u, x) 0 t x
以中心差分方法为例来说明。 对于第 i 点:
(1.1)
ui (u u ) ai i 1/ 2 i 1/ 2 0 t x
其中 ai 第 i+1 点:
(1.2)
1 (ai 1/ 2 ai 1/ 2 ) 2
ui 1 (u u ) ai 1 i 3/ 2 i 1/ 2 0 t x
计算流体力学 有限体积法基础及其应用
一、计算流体力学简介1.1 计算流体力学的定义1.2 计算流体力学的研究对象1.3 计算流体力学的发展历史二、有限体积法基础2.1 有限体积法的理论基础2.1.1 有限体积法的基本原理2.1.2 有限体积法的数学模型2.2 有限体积法的数值求解2.2.1 离散化2.2.2 迭代求解三、有限体积法在计算流体力学中的应用3.1 有限体积法在流体流动模拟中的应用 3.1.1 管道流动模拟3.1.2 自由表面流动模拟3.2 有限体积法在传热问题中的应用3.2.1 对流传热3.2.2 辐射传热四、有限体积法在工程领域中的应用4.1 有限体积法在航空航天领域中的应用 4.2 有限体积法在汽车工程中的应用4.3 有限体积法在建筑工程中的应用五、有限体积法的发展趋势5.1 高性能计算技术对有限体积法的影响5.2 多物理场耦合对有限体积法的挑战5.3 人工智能在有限体积法中的应用六、结论一、计算流体力学简介1.1 计算流体力学的定义计算流体力学（Computational Fluid Dynamics, CFD）是利用计算机模拟流体力学问题的一门学科。
它通过对流动流体的数值解，来研究流体在各种情况下的运动规律和性质。
1.2 计算流体力学的研究对象计算流体力学的研究对象包括流体的流动、传热、传质、振动等现象，以及与流体相关的各种工程问题，如飞机、汽车、建筑等的气动特性分析与设计。
1.3 计算流体力学的发展历史计算流体力学的发展可以追溯到20世纪50年代，当时计算机技术的进步为流体力学问题的数值模拟提供了可能。
随着计算机硬件和软件的不断发展，CFD的应用领域不断扩大，成为现代工程领域不可或缺的工具之一。
二、有限体积法基础2.1 有限体积法的理论基础2.1.1 有限体积法的基本原理有限体积法是求解流体动力学问题的数值方法之一，它基于质量、动量和能量守恒的控制方程，将求解域离散化为有限数量的体积单元，通过对控制方程进行积分，将方程转化为代数方程组。
计算流体力学课件-part1
2024/2/28
19
❖Computational Fluid Dynamics
计算流体流体力学
第二讲 典型模型方程的数学性质
模型方程的概念
➢完整方程
连续方程
动量方程
能量方程
2024/2/28
20
❖Computational Fluid Dynamics
沿特征线，扰动波的幅值不变，传播速度为c
则在t>0时，传播过程如下图：
2024/2/28
27
❖Computational Fluid Dynamics
计算流体流体力学
第二讲 典型模型方程的数学性质
模型方程的特征
➢单波方程
➢c>0时，传播沿x正向 ➢C<0时，传播沿x负向 ❖扰动波以有限速度传播是双曲型方程的重要 特征（波形和波幅可能会变化，此处为什么不 变？）
如何表达初始形状三角形
如何存储数据 如何积分
数值积分，HOW？
如何显示结果
TECPLOT
尝试改变几个常数，看看结果有何变化，常数反映了什么？
2024/2/28
22Biblioteka ❖Computational Fluid Dynamics
回顾
控制方程
模型方程
➢NS ➢EULER ➢Impressible NS ➢RANS
➢单波方程可以模拟EULER方程的一些特征
2024/2/28
28
❖Computational Fluid Dynamics
计算流体流体力学
第二讲 典型模型方程的数学性质
模型方程的特征
计算流体力学基础
物理模型与数学模型在概念上的区别
数学模型：对物理模型的数学描写。
比如N-S方程就是对粘性流体动力学的一种数学描写，值得注意的 是，数学模型对物理模型的描写也要通过抽象，简化的过程。
物理模型是指把实际的问题，通过相关的物理定律概括和抽象出来并满足 实际情况的物理表征。
比如，我们研究管道内的流体流动，抽象出来一个直管，和粘性流体模型， 或者我们认为管道内的液体是没有粘性的，使用一个直管和无粘流体模型. 还有，我们根据热传导定律，认为固体的热流率是温度梯度的线形函数， 相应的傅立叶定律就是导热问题的物理模型。因此，不难理解物理模型是 对实际问题的抽象概念，对实际问题的一种描述方式，这种抽象包括了实 际问题的几何模型，时间尺度，以及相应的物理规律。
确定边界条件与初始条件 初始条件与边界条件是控制方程有确定解的前提，控制方程与 相应的初始条件、边界条件的组合构成对一个物理过程完整的数学 描述。 初始条件是所研究对象在过程开始时刻各个求解变量的空间分 布情况。对于瞬态问题，必须给定初始条件。对于稳态问题，不需 要初始条件。 边界条件是在求解区域的边界上所求解的变量或其导数随地点 和时间的变化规律。对于任何问题，都需要给定边界条件。例如， 在锥管内的流动，在锥管进口断面上，我们可给定速度、压力沿半 径方向的分布，而在管壁上，对速度取无滑移边界条件。 对于初始条件和边界条件的处理，直接影响计算结果的精度。
划分计算网
采用数值方法求解控制方程时，都是想办法将控制方程在空 间区域上进行离散，然后求解得到的离散方程组。要想在空间域 上离散控制方程，必须使用网格。现已发展出多种对各种区域进 行离散以生成网格的方法，统称为网格生成技术。 不同的问题采用不同数值解法时，所需要的网格形式是有一 定区别的，但生成网格的方法基本是一致的。目前，网格分结构 网格和非结构网格两大类。简单地讲，结构网格在空间上比较规 范，如对一个四边形区域，网格往往是成行成列分布的，行线和 列线比较明显。而对非结构网格在空间分布上没有明显的行线和 列线。
流体计算理论基础讲解
流体计算理论基础讲解(总32页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--流体计算理论基础1 三大基本方程连续性方程连续性方程也称质量守恒方程，任何流动问题都必须满足质量守恒定律，该定律可表示为：单位时间内流体微元中质量的增加等于同一时间间隔内流入该微元体的净质量，其形式如下：()()()0u v w t x y zρρρρ∂∂∂∂+++=∂∂∂∂ 可以写成：()0div u tρρ∂+=∂ 其中ρ密度，t 为时间，u 为速度矢量，u ，v 和w 为速度矢量在x ，y 和z 方向上的分量。
若流体不可压缩，密度为常数，于是：0u v w x y z∂∂∂++=∂∂∂ 若流体处于稳态，则密度不随时间变化，可得出：()()()0u v w x y zρρρ∂∂∂++=∂∂∂ 动量守恒定律该定律可以表述为：微元体中流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和，该定律实际是牛顿第二定律，按照这一定律，可导出x ，y 和z 三个方向上的动量守恒方程：()()()()()()yx xx zx x xy yy zy y yz xz zz z u p div uu F t x x y z u p div uv F t y x y z u p div uw F tz x y z τττρρτττρρτττρρ∂⎧∂∂∂∂+=-++++⎪∂∂∂∂∂⎪⎪∂∂∂∂∂⎪+=-++++⎨∂∂∂∂∂⎪⎪∂∂∂∂∂+=-++++⎪∂∂∂∂∂⎪⎩式中，p 为微元体上的压力，xx τ，xy τ和xz τ等是因分子粘性作用而产生的作用在微元体表面上的粘性应力τ的分量。
x F ，y F 和z F 是微元体上的体力，若体力只有重力，且z 轴竖直向上，则：0,0x y F F ==，z F g ρ=-。
对于牛顿流体，粘性应力τ与流体的变形率成比率，有：x yy x 2();==()2();==()2();==()xx xy y xz z zz yz zy u u v div u x y x v u w div u x z x w v w div u x z y τμλττμτμλττμτμλττμ∂∂∂⎧=++⎪∂∂∂⎪∂∂∂⎪=++⎨∂∂∂⎪∂∂∂⎪=++⎪∂∂∂⎩其中，μ为动力粘度，λ为第二粘度，一般可取23λ=-，将上式代入前式中为：()()()()()()()()()u v w u p div uu div gradu S t x v p div uv div gradv S ty w p div uw div gradw S tz ρρμρρμρρμ⎧∂∂+=-+⎪∂∂⎪∂∂⎪+=-+⎨∂∂⎪⎪∂∂+=-+⎪∂∂⎩ 其中：()()/()/()/grad x y z =∂∂+∂∂+∂∂μ为动力粘度(dynamic viscosity)，λ为第二粘度(second viscosity),一般可取：23λ=-(参考文献：,Boundary Layer Theory,8th ed,McGraw Hill, New York,1979)。
计算流体力学课件完整版
●实验要受测量技术限制，实验周期长、费用高。
☆ 理论研究 ●在研究流体流动规律的基础上，建立了流体流动基 本方程。 ●对于一些简单流动，通过简化求出研究问题的解析 解。
计算流体力学
●对于实际流动问题，通常需运用流体力学基本方程， 借助于计算机求数值解（计算机数值模拟）— 计算流体力学CFD。
Z
skirt.plt X Y
75 50 25
0 -25 -50 -75
-2
Y(M) 0
2
0 2 4 6 10 8 X(M) 12 14
D) 16 Feb 2003 Velocity Vectors
4.5
4 velocity.plt
3.5
3
2.5
2
1.5
Z
Z
(3D) 16 Feb 2003 IJK-Ordered DZ ata
ijkcyl.plt X Y
Z
-0.4 -0.2 Y0 0.2 0.4
1
0.8
0.6
0.4
0.2
0 -0.4 -0.2 0 X 0.2 0.4
Z
jetflow.plXt Y
0.6 0.5 0.4 0.3 0.2 0.1
0 0 Y0.1 0.2
-0.6 -0.4 -0.2 0 X 0.2 0.4 0.6
轴流叶轮计算与实验叶片表面极限流线
计算流体力学
轴流叶轮计算与实验性能比较
计算流体力学
轴流叶轮计算与实验流场结构比较
计算流体力学
第二章 流体力学数值计算数学模型及定解条件
☆本章所涉及的基本方程有两类： ●流体力学基本方程，基本出发点：质量守恒、动量守恒和能
计算流体力学(中科院力学所)_第讲-基本方程ppt课件
YF23
7
● 90年代， CFD 在飞机设计中发挥了主力作用 波音777， CFD占主角
● 2000 之后， CFD 取代了大部分风洞实验 波音787：全机风洞实验仅3次
● 航天领域，CFD发挥着实验无法取代的作用 实验难点：复现高空高速流动条件
波音777
Copyright by Li Xinliang
s
s
控制体内的动量增加=流入的动量+表面力的冲量+体积力的冲量
t V d [ (V V ) F P ]d
V (V V )F P
t
Copyright by Li Xinliang
12
基本概念： 应力 （张量）
pn Pn
pn
根据本构方程（广义牛顿粘性定律）
Pijpijij ：静止部分+运动部分
✓基本概念: 随体导数 dV
dt t
11
2） 动量守恒律
单位时刻内，流出面元ds的动量为：
d V d m V V n dS
总流出动量为:
d （ V V ） n d S ( V V ) d
S
s
外力的合力：
质量力：Fd 表面力：
根据动量守恒：
p nd SP n d S P d
控制体
单位时刻表面微元ds的流出质量为： dm V n dS
V
总质量流出为 d m V n d S (V )d
n
s
s
根据质量守恒： 控制体内质量的增加=流入控制体的质量
dS
控制体的任意性
td d m (V )d
s
(V)0
t
（1） Copyright by Li Xinliang
计算流体力学基础ppt课件
它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性， 能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、 易燃等真实条件和实验中只能接近而无法达到的理想条件。
8
数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适 用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果 不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并 有一定的计算误差。
对于初始条件和边界条件的处理，直接影响计算结果的精度。
16
划分计算网 采用数值方法求解控制方程时，都是想办法将控制方程在空间区
域上进行离散，然后求解得到的离散方程组。要想在空间域上离 散控制方程，必须使用网格。现已发展出多种对各种区域进行离 散以生成网格的方法，统称为网格生成技术。
不同的问题采用不同数值解法时，所需要的网格形式是有一定区 别的，但生成网格的方法基本是一致的。目前，网格分结构网格 和非结构网格两大类。简单地讲，结构网格在空间上比较规范， 如对一个四边形区域，网格往往是成行成列分布的，行线和列线 比较明显。而对非结构网格在空间分布上没有明显的行线和列线。
数学模型就好理解了，就是对物理模型的数学描写。 比如N-S方程就是对粘性流体动力学的一种数学描写，值得注意的是，数学 模型对物理模型的描写也要通过抽象，简化的过程。
14
建立控制方程 确立初始条件及边界条件 划分计算网格，生成计算节点
建立离散方程
离散初始条件和边界条件
给定求解控制参数
解收敛否
否
显示和输出计算结果
21
给定求解控制参数 在离散空间上建立了离散化的代数方程组，并施加离散化的
计算流体力学part基础知识PPT课件
①在直角坐标系中：A P(x, y, z)i Q(x, y, z) j R(x, y, z)k
i rotA
x P
jk y z Q R 第18页/共56页
（21）
一、向量分析初步
5、向量场的环量及旋度
rot A 0 有旋运动， rot A 0 无旋运动。应当指出，流体微团 是否作有旋运动，需视微团是否围绕着通过流体微团的瞬时 轴旋转，而并非决定于流体微团轨迹的几何形状。
a(t) ax (t)i ay (t) j az (t)k （10） 结论：
向量导数在坐标轴上的投影等于相应的向量投 影的导数。
向量的导数在几何上为一切向矢量。
da(t) a(t) dt
第10页/共56页
一、向量分析初步
2、向量函数对于数变量的导数
一个流体微团在空间的位置可用坐标 x, y, z 确定，也可用向径确定：
一、向量分析初步
2、向量函数对于数变量的导数
da(t) lim a(t) dt t0 t
lim
t0
ax (t t
)
i
ay (t) t
j
az (t) t
k
dax (t) i day (t) j daz (t) k
dt
dt
dt
第9页/共56页
一、向量分析初步
2、向量函数对于数变量的导数
dx i dt
dy dt
j
dz dt
k
vxi vy j vzk
第11页/共56页
（11）
一、向量分析初步
3、数量场的梯度
若在数量场 x, y, z 中的一点 p
处，存在着矢量 G ，其方向为函数
流体力学与计算流体力学基础
f = f xi + f y j + f z k
(1-14)
式中： f x ， f y ， f z 为单位质量力在 x ， y ， z 轴上的投影，或简称为单位质量分力。 表面力： 大小与表面面积有关而且分布作用在流体表面上的力称为表面力。 表面力按其作 用方向可以分为两种：一种是沿表面内法线方向的压力，称为正压力；另一种是沿表面切向的 摩擦力，称为切应力。 作用在静止流体上的表面力只有沿表面内法线方向的正压力。 单位面积上所受到的表面力 称为这一点处的静压强。静压强有两个特征： 静压强的方向垂直指向作用面。 流场内一点处静压强的大小与方向无关。
α =−
dρ ρdT
(1-12)
α 的单位为 1/ K 。 这里的负号是考虑到随着温度的增高， 体积必然增大， 而密度必然减小； 体积膨胀系数的物理意义为当压强不变时，每增加单位温度所产生的流体体积的相对变化率。 气体的体积膨胀系数可由气体状态方程求得：
α = 1/ T
(1-13)
在研究流体流动过程时，若考虑到流体的压缩性，则称为可压缩流动，相应地称流体为可 压缩流体，如相对速度较高的气体流动。 若不考虑流体的压缩性，则称为不可压缩流动，相应地称流体为不可压缩流体，如水、油 等液体的流动。 （6）液体的表面张力 液体表面相邻两部分之间的拉应力是分子作用力的一种表现。 液面上的分子受液体内部分
6
Fluent 17.0 流体仿真从入门到精通
图 1-1
边界层示意图
大雷诺数边界层流动的性质： 边界层的厚重较物体的特征长度小得多，即 δ / L （边界层相对厚度）是一个小量。边界 层内粘性力和惯性力同阶。 对于二维平板或楔边界层方程，通过量阶分析得到：
∂u ∂v + =0 ∂x ∂y ∂U ∂ 2u ∂u ∂u ∂u ∂U +v 2 +u +v = +U ∂x ∂y ∂t ∂x ∂y ∂t
流体力学 第14章 计算流体力学基础
行求解。近似公式应用在空间和时间的小域上,从而通过求解微分方程的数值解,
得到离散空间各个小域上具体物理量的数值,给出数值结果,这就是计算流体力学
的基本数学指导思想。
计算流体力学的基本原理
利用计算流体力学对流动问题进行数值模拟时,通常包括如下四个步骤:
• 有限体积法（FVM）——控制体内的平均近似
出发点是守恒型方程的积分形式,求解域被分成若干连续的控制体。在每
一个控制体上满足守恒方程。在每一个控制体的中心作为计算节点,计算该点
上的物理量。控制体边界上的函数值用节点函数值的插值获得。体积分和面
积分用适当的求积公式近似。结果在每个控制体上都有一个代数方程,未知数
网格线的任意一条有且只有一个交点。
(2)非结构化网格
指网格区域内的内部点不具有相同的毗邻单元,即与
网格剖分区域内的不同内点相连的网格数目不同。
计算流体力学的基本要素
• 有限近似
在选定数值网格以后,还必须确定数值离散过程中的近似方法。
近似程度决定了数值求解的精度以及求解的难度和费用。高精度格式的
方程中包含了更多网格节点数,因此求解的工作量和难度也相应地增加。
• 坐标和矢量系统
流体力学的基本方程与坐标无关,但在不同的坐标系下有不同的表达形式,
因此在数值计算时必须选择合适的坐标系,此外,矢量在该坐标系下的表达形
式也必须事先予以确定。
• 数值网络
数值网格定义了所求物理量在空间的位置。数值网格的种类大致可分为:
(1)结构化网格
由多族网格线构成,同族的网格线互不相交,且和其他族
返回目录
3.计算流体力学的控制方程组
计算流体力学基础
For personal use only in study and research; not for commercial use一、计算流体力学的基本介绍一、什么是计算流体力学(CFD)?计算流体力学(Computational Fluid Dynamics)是流体力学的一个新兴的分支，是一个采用数值方法利用计算机来求解流体流动的控制偏微分方程组，并通过得到的流场和其它物理场来研究流体流动现象以及相关的物理或化学过程的学科。
事实上，研究流动现象就是研究流动参数如速度、压力、温度等的空间分布和时间变化，而流动现象是由一些基本的守恒方程（质量、动量、能量等）控制的，因此，通过求解这些流动控制方程，我们就可以得到流动参数在流场中的分布以及随时间的变化，这听起来似乎十分简单。
但遗憾的是，常见的流动控制方程如纳维一斯托克斯(Navier-Stokes)方程或欧拉(Euler)方程都是复杂的非线性的偏微分方程组，以解析方法求解在大多数情况下是不可能的。
实际上，对于绝大多数有实际意义的流动，其控制方程的求解通常都只能采用数值方法的求解。
因此，采用CFD方法在计算机上模拟流体流动现象本质上是流动控制方程（多数情况下是纳维一斯托克斯方程或欧拉方程）的数值求解，而CFD软件本质上就是一些求解流动控制方程的计算机程序。
二、计算流体力学的控制方程计算流体力学的控剖方程就是流体流动的质量、动量和能量守恒方程。
守恒方程的常见的推导方法是基于流体微元的质量、动量和能量衡算。
通过质量衡算可以得到连续性方程，通过动量守恒可以得到动量方程，通过能量衡算可以得到能量方程。
式(1)一(3)是未经任何简化的流动守恒微分方程，即纳维一斯托克斯方程( N-S方程)。
N-S方程可以表示成许多不同形式，上面的N-S方程是所谓的守恒形式，之所以称为守恒形式，是因为这种形式的N-S方程求解的变量p、pu、pv、pw、pE是守恒型的，是质量、动量和能量的守恒变量。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
经典差分格式及其性质
?u ? c ?u ? 0 ?t ?x
u ?x,0 ?? φ?x?
精确解
u ?x, t ?? φ ?x ? ct ?
FTBS Lax 格式
un?1 j
?
u
n j
?
c
u
n j
?
un j?1
?
0
?t
?x
稳定条件 CFL ? 1
? ? u
n j
?
1
?
1 2
??u
n j?1
?
u
n j?1
CDS (二阶)
一阶格式
达到与均匀网格同样的收敛速度：O (D x 2)
离散格式精度（阶数）的含义
? 当Δx? 0(即对网格不断加密)时，截断误差 减小的程度(收敛速度)。
? 高阶格式的优势是在网格不断加密的过程 中体现出来的，所以：
? 只有在网格足够密时，高阶格式才具有与 之相称的高精度。
均匀网格二分加密：
有限体积法基本思想
? 多项式插值
三点二次多项式插值：
? 格式精度等于插值次数减一 ? 偶数次插值、均匀网格：精度加一
五点四次多项式插值：
二阶导数差分格式
扩散项处理：
高精度格式
? 高精度格式需要更多节点 ? 方程更难求解 ? 边界条件更难处理 ? 二阶格式常用
耗散与色散
耗散
数值扩散系数ν数
? 耗散只存在于对流项的一阶离散格式 ? 正比于网格步长
离散方法
? 有限差分 ? 有限元 ? 有限体积
有限差分
有限差分基本思想
结构化网格：
?f ruj ?x
微分-> 差分：
有限差分：Taylor 级数展开
FDS： BDS： CDS：
截断误差：
有限差分：多项式插值
二阶CDS：
均匀网格高阶格式：
迎风(Upwind)格式
非均匀网格
? 截断误差 ~ 网格尺寸*导数
董f , D x ? 动f , D x ?
误差均匀分布
? 节点数给定，非均匀网格可使误差降到最小 ：网格自适应
非均匀网格：加密与截断误差
CDS：
一阶
CDS均匀网格：
Dx2
6
二阶
对非均匀网格进行系统加密可以
达到与均匀网格同样的收敛速度：O (D x 2)
非均匀网格：加密与截断误差
加密方法
? 简单二分法 ? 固定增长因子
计算流体力学
CFD基本要素
? 数学模型 ? 离散方法 ? 坐标系 ? 计算网格 ? 有限近似 ? 求解方法 ? 收敛准则
CFD方法性质
? 方程的相容性(Consistency) ? 数值算法的稳定性 ? 解的收敛性
? Lax等价定理
? 守恒性 ? 有界性 ? 模型的可实现性 ? 精度
? 模型误差、离散误差、迭代误差
?
σ
u
n j?
1
?
u
n j?1
??
σ ? c? t / ? x
Leap-Frog
un?1 j
?
u n?1 j
?
f
n j?1
?
fn j?1
?
0,
(f ? cu)
2? t
2? x
Lax-Wandrof
u
n? j
1
?
u
n j
?
c
?t ?x
u
n j?1
?
u
n j?1
2
?
1
?
?
t
2
?
2 ?? ? x ??
δx? cδx? cu j
非均匀网格：加密与截断误差
( )e
r o t 2h
( ) t
et h
非均匀网格CDS (一阶)收敛速度：
ì??í??>=
4, 4,
re = 1 re ? 1
对非均匀网格进行系统加密可以
均匀网格：
r t
=
ì???????í???????? ((((DDDDxxxx))))2222hhhh
= =
4, 2,
u ?u ?
α ?u ? g ?x?
?n ?
α1
?x ?? u
?n
?
? α2 ?x?u ?
?
h ?x ?
? 边界附近的差分可能要用到边界之外的节点 ? 如果不降低精度，需采用不同的差分格式
四次多项式插值：
边界导数差分格式
FDS：
内部节点三次插值：
离散后的代数方程
有限差分
微分方程
代数方程
?f ru
un j?1
?
u
n j
? ? ? ? ? f
? ? a j?1 2
?
?? ?
u j?1 ? f u j u j?1 ? u j
, u j?1 ? u j ? 0
? ?
a
?u
j
?,
? ? u j?1 ? u j ? 0
经典差分格式性质比较
，精确解
1 0.5
u u0, 0
-0.5
u0 FTBS lf lw lpfrog bw roe bwsplit
j ?x
计算网格单元与节点存储
代数方程组
合适的节点排序
稀疏矩阵 带矩阵
离散误差
微分方程真实解 离散方程真实解
线性问题：
截断误差 离散误差
离散误差估计与Richardson外插
? 网格足够细 ? 单调收敛
Richardson
外插:
fR h
=
f h
+
ed h
: O (hr ), r > p
有限体积法
? ? f ? ? c? u, c? ? c ? c / 2
f ? f ? ? f ? , c ? c? ? c?
? ? un?1 j
?
u
n j
?
?t ?x
f?jn?1 2 ? f?jn?1 2
? ? ? ? ? ? f?jn?1 2
?
1 2
??
f
u
n j
?f
un j?1
?? ?
1 2
an j?1 2
网格Peclet数与网格Reynolds数
数值扩散系数ν数
数值耗散不能大于物理耗散： 或
色散
色散关系： CDS approximation of a wave:
差分格式只能正确模拟低频分量（相对于kmax）
色散与振荡
CDS
UDS
? 网格较粗时产生较大的虚假低频分量 ? 一阶格式的强耗散抹平了振荡
r t
=
ì???????í???????? ((((DDDDxxxx))))222hh2hh
= =
2, 4,
一阶格式
CDS (二阶)
r =9 t r =9 t
r = 16 t
二阶导数差分格式构造方法
? 对一阶导数差分
二阶导数差分格式构造方法
? Taylor 级数展开：量 ，
一阶
二阶导数差分格式构造方法
-1
-1
-0.5
0
0.5
1
x
经典差分格式性质比较
，精确解
1
0.8
0.6
u
0.4
Initial ftbs lf lw lpfrog bwsplit roe bw
0.2
0-1
-0.5
0
0.5
1
x
边界条件差分格式
偏微分方程边界条件：
1、Dirichlet问题 2、Neumann问题 3、Robin问题
经典差分格式及其性质
Beam-Warming 流通矢量分裂
Roe格式
? ? un? 1 ? ? j
?
u
n j
?
c? t
3u
n j
?
4u
n j?
12? x?来自un j? 2?
c2
?t
2
u
n j
?
2u
n j?1
?
un j?2
2
?x 2
uin ? 1
?
u
n j
?
?
? x
f
? j
?
?
? x
f
? j
?
0
?t
?x