CFD控制方程离散方法：有限容积法
cfd控制方程的离散方法
cfd控制方程的离散方法
CFD（Computational Fluid Dynamics，计算流体力学）是一种利用数值方法解决流体力学问题的技术。
在CFD中，控制方程是描述流体运动的基本方程，包括质量守恒方程、动量守恒方程和能量守恒方程。
离散方法是将连续的物理方程转化为离散的代数方程，以便通过计算机进行求解。
离散方法常用的有有限差分法（Finite Difference Method）、有限体积法（Finite Volume Method）和有限元法（Finite Element Method）。
对于CFD中的控制方程，离散方法的选择取决于问题的性质和所需的精度。
以下是几种常用的离散方法：
1. 有限差分法：将微分算子近似为差分形式，通过在网格上进行逐点近似来离散化方程。
有限差分法简单易用，适用于规则网格和简单几何形状的问题。
2. 有限体积法：将控制方程应用到一个控制体积（Control Volume）上，使用积分形式得到离散化的方程。
有限体积法适用于复杂几何形状和非结构网格，能够保持物理量的守恒性。
3. 有限元法：将求解域划分为离散的有限元，使用基函数对方程进行近似。
有限元法适用于复杂几何形状和非结构网格，能够处理不规则网格以及局部自适应网格细化。
这些离散方法各有优缺点，需要根据具体问题的性质和要求选择合适的方法。
同时，为了保证数值解的准确性和稳定性，还
需要考虑网格的划分方式、边界条件的处理以及迭代求解算法等因素。
CFD控制方程离散方法_有限容积法
有限容积法
离散与数值解
方程非线性
各个变量通过系数等相互影响：内耦合
把原来在空间与时间坐标中连续的物理量的场，用一系列有限个离散点上 的值的集合来代替，建立起离散方程，求解所建立起来的代数方程以获得所求 解变量的近似值。 分离式求解：在一个方程中认为系数（包括速度）为常数，采用上一次迭代 算出来的值。 迭代求解线性方程组，从一个初始的（猜测的）变量分布开始，不断更新直至 收敛。
数值传热学基本思想
有限差分法
一维稳态有源项的对流-扩散方程 泰勒展开
得到线性方程组
有限容积法
在一个控制体内做积分
积分之后对方程进行离散
有限容积法
扩散项线性化处理（中心差分） 源项认为在整个控制体内不变，局部线性化处理 两种型线假设
有限容积法
求解线性方程组
有限容积法
有限容积法
有限容积法
迎风算法 扰动不向上游传递 一阶迎风的数值误差较大，在调试阶段可采用。
知识回顾 Knowledge Review
放映结束 感谢各位的批评指导！
谢 谢！
让我们共同进步
3、输运性Transportiveness
邻点W和E有两个恒定源，画出等值线
纯扩散
对流扩散
纯扩散使源的影响向各个方向同等地传播；纯对流时，P点只受上游影响不受下游影响。
无源时场随时间变化
3、输运性Transportiveness
n时刻 n+1时刻
扰动被均匀向两侧传递
中心差分
对流项的中心差分不合理，因为aE为负，使得下游增大会使上游 减小
在均分网格情况下与Taylor展开法的结果一致。 在FVM中所谓不同的格式就是指不同的型线。
分离式求解过程
有限元素法有限体积法有限差分法有限容积法的区别
1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。
该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。
有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。
该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。
1.2 差分格式（1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。
（2）从差分的空间形式来考虑，可分为中心格式和逆风格式。
（3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。
目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。
差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
1.3 构造差分的方法构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。
其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。
通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。
2. FEM2.1 概述有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
采用不同的权函数和插值函数形式，便构成不同的有限元方法。
2.2 原理有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。
在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。
cfd离散的四项法则
CFD离散的四项法则1.离散化方法离散化是计算流体动力学（CFD）中的核心步骤，它涉及到将连续的物理空间和时间转化为离散的数值网格。
离散化的目的是将偏微分方程转换为数值求解的差分方程，以便在计算机上进行数值模拟和分析。
常见的离散化方法包括有限差分法、有限元法和有限体积法等。
这些方法各有优缺点，适用于不同的流动和几何形状。
2.离散格式在离散化过程中，需要对偏微分方程中的各个导数项进行离散化。
不同的离散格式会导致不同的数值精度和稳定性。
常见的离散格式包括中心差分格式、前向差分格式、后向差分格式和混合差分格式等。
选择合适的离散格式对于保证数值模拟的精度和稳定性至关重要。
3.时间积分方案时间积分方案决定了如何推进求解的进程，即在离散的时间步长上逐步求解离散的差分方程。
常见的时间积分方案包括隐式方案、显式方案和半隐式方案等。
隐式方案具有较高的稳定性和精度，但计算量较大；显式方案稳定性和精度较低，但计算量较小；半隐式方案则结合了隐式和显式的优点，具有较好的稳定性和精度，同时计算量也相对较小。
4.离散方程的求解方法在CFD中，离散方程的求解方法通常包括迭代法和直接法。
迭代法是通过不断迭代来逼近方程的解，常见的迭代法包括Jacobi迭代法和Gauss-Seidel迭代法等。
直接法则是通过一定的算法直接求解方程的解，常见的直接法包括高斯消去法和LU分解法等。
选择合适的求解方法可以提高计算效率，并保证数值模拟的准确性。
以上是CFD离散的四项法则中各重要元素的简单概述。
在实际应用中，需要根据具体问题选择合适的离散化方法、离散格式、时间积分方案和离散方程的求解方法。
在保证数值模拟的精度和稳定性的同时，提高计算效率是CFD模拟的关键。
随着计算机技术的不断发展，CFD的应用范围越来越广泛，CFD技术也面临着新的挑战和机遇。
未来，CFD技术将不断发展和完善，为流体动力学、气象学、环境科学等领域提供更加精确和可靠的数值模拟和分析工具。
cfd控制方程理解
CFD控制方程是指流体力学中的基本方程，用于描述流体运动的规律。
这些方程包括质量守恒方程、动量守恒方程和能量守恒方程。
质量守恒方程，也称为连续性方程，表示单位时间内流体质量的增加等于该时间间隔内流入该控制体的净质量。
动量守恒方程，也称为动量方程或N-S方程，描述了流体动量的变化规律。
能量守恒方程则描述了流体能量的变化规律，包括机械能、内能和热能等。
这些控制方程可以用来描述各种复杂的流体流动现象，例如湍流、分离流、波动等。
通过求解这些控制方程，可以得到流场内各点的速度、压力、温度等参数的分布情况，从而对流体运动进行预测和分析。
在实际应用中，CFD控制方程通常需要使用数值方法进行求解。
常用的数值方法包括有限差分法、有限元法和有限体积法等。
通过这些数值方法，可以将控制方程离散化，将其转化为一系列代数方程组，然后使用计算机进行求解。
总之，CFD控制方程是流体力学中的基本方程，用于描述流体运动的规律。
通过求解这些方程，可以对流体运动进行预测和分析，为工程设计和优化提供重要的依据。
03_控制方程的离散化方法
7xxxxΔΔΔCartesian 坐标系下输运方程的表达式•有限差分的概念是从导数的定义中得到的：•几种常用的差分格式812有限元法的插值函数•第一类：只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；•第二类：不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。
有限元法的单元坐标•单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。
•常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状。
13有限元法的求解步骤•建立积分方程•区域单元剖分•确定单元基函数•单元分析•总体合成•边界条件的处理•解有限元方程14有限元法在很多地方与FV 非常类似，有限元法也将空间分成连续的控制体，但是通常采用非结构网格。
FE 法的最明显特点是方程在整个域内积分之前被乘上了一个权函数，在最简单的FE 法中，每个控制体内的函数被假设成线性的，并且保证解在边界上的连续性。
有限元法的求解速度比有限差分法和有限体积法，因此，商用CFD 软件中应用并不普遍。
而有限元法目前在固体力学分析中占绝对比例，几乎所有软件都采用有限元法。
15有限体积法(Finite Volume Method)——控制体内的平均近似有限体积法是近年发展非常迅速的一种离散化方法，在CFD 领域广泛应用。
大多数商用CFD 软件都采用这种方法。
特点是计算效率高，不紧表现在对控制方程的离散结果上，还表现在所使用的网格上。
163.3 有限体积法基本思想有限体积法又称控制体积法(Control Volume Method)，基本思路是：将计算区域划分一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解微分方程(控制方程)对每一个控制体积积分，得到一组离散方程。
从积分区域的选取方法来看，有限体积法属于加权余量法中的子域法。
从未知解的近似方法来看，有限体积法属于采用局部近似的离散方法。
有限差分,有限元,有限体积等离散方法的区别介绍
有限差分，有限元，有限体积等等离散方法的区别介绍一、区域离散化所谓区域离散化，实质上就是用一组有限个离散的点来代替原来连续的空间。
实施过程是；把所计算的区域划分成许多互不重迭的子区域，确定每个子区域的节点位置及该节点所代表的控制容积。
节点：需要求解的未知物理量的几何位置；控制容积：应用控制方程或守恒定律的最小几何单位。
一般把节点看成是控制容积的代表。
控制容积和子区域并不总是重合的。
在区域离散化过程开始时，由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。
网格是离散的基础，网格节点是离散化物理量的存储位置。
大家都知道，常用的离散化方法有：有限差分法，有限元法，有限体积法。
1. 有限差分法是数值解法中最经典的方法。
它是将求解区域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程（控制方程）的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。
这种方法发展比较早，比较成熟，较多用于求解双曲线和抛物线型问题。
用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。
2. 有限元法是将一个连续的求解域任意分成适当形状的许多微小单元，并于各小单元分片构造插值函数，然后根据极值原理（变分或加权余量法），将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为各单元极值之和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组就得到各节点上待求的函数值。
对椭圆型问题有更好的适应性。
有限元法求解的速度较有限差分法和有限体积法慢，在商用CFD软件中应用并不广泛。
目前的商用CFD软件中，FIDAP采用的是有限元法。
3. 有限体积法又称为控制体积法，是将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积，将待解的微分方程对每个控制体积积分，从而得到一组离散方程。
其中的未知数十网格节点上的因变量。
子域法加离散，就是有限体积法的基本方法。
就离散方法而言，有限体积法可视作有限元法和有限差分法的中间产物。
计算流体力学常用的五大类数值方法简介
计算流体力学常用的五大类数值方法简介流体力学数值方法有很多种，其数学原理各不相同，但有二点是所有方法都具备的，即离散化和代数化。
总的来说其基本思想是：将原来连续的求解区域划分成网格或单元子区域，在其中设置有限个离散点(称为节点)，将求解区域中的连续函数离散为这些节点上的函数值；通过某种数学原理，将作为控制方程的偏微分方程转化为联系节点上待求函数值之间关系的代数方程(离散方程)，求解所建立起来的代表方程以获得求解函数的节点值。
不同的数值方法，其主要区别在于求解区域的离散方式和控制方程的离散方式上。
在流体力学数值方法中，应用比较广泛的是有限差分法、有限元法、边界元法、有限体积法和有限分析法，现简述如下。
一、有限差分法这是最早采用的数值方法，它是将求解区域划分为矩形或正交曲线网格，在网格线交点(即节点)上，将控制方程中的每一个微商用差商来代替，从而将连续函数的微分方程离散为网格节点上定义的差分方程，每个方程中包含了本节点及其附近一些节点上的待求函数值，通过求解这些代数方程就可获得所需的数值解。
有限差分法的优点是它建立在经典的数学逼近理论的基础上，容易为人们理解和接受；有限差分法的主要缺点是对于复杂流体区域的边界形状处理不方便，处理得不好将影响计算精度。
二、有限元法有限元法的基本原理是把适定的微分问题的解域进行离散化，将其剖分成相连结又互不重叠的具有一定规则几何形状的有限个子区域(如：在二维问题中可以划分为三角形或四边形；在三维问题中可以划分为四面体或六面体等)，这些子区域称之为单元，单元之间以节点相联结。
函数值被定义在节点上，在单元中选择基函数(又称插值函数)，以节点函数值与基函数的乘积的线性组合成单元的近似解来逼近单元中的真解。
利用古典变分方法(里兹法或伽辽金法)由单元分析建立单元的有限元方程，然后组合成总体有限元方程，考虑边界条件后进而求解。
由于单元的几何形状是规则的，因此在单元上构造基函数可以遵循相同的法则，每个单元的有限元方程都具有相同的形式，可以用标准化的格式表示，其求解步骤也就变得很规范，即使是求解域剖分各单元的尺寸大小不一样，其求解步骤也不用改变，这就为利用计算机编制通用程序进行求解带来了方便。
试采用有限容积法推导出关于内节点p的显式离散化方程 题目
试采用有限容积法推导出关于内节点p的显式离散化方程题目采用有限容积法推导内节点p的显式离散化方程，首先需要明确几个关键概念和步骤。
1. 内节点定义：在计算区域内，位于边界之内的节点称为内节点。
2. 控制容积：对于内节点p，其控制容积是指围绕该节点的封闭几何区域。
3. 质量守恒方程：对于控制容积，质量守恒方程为：流入控制容积的质量 - 流出控制容积的质量 = 控制容积内原有质量的变化。
4. 离散化：将连续的质量守恒方程转化为离散的形式，以便于数值计算。
基于以上概念，推导过程如下：假设内节点p的控制容积为V，其周围有m个网格单元。
对于每个网格单元，假设其体积为ΔV，质量为Δm。
根据质量守恒方程，我们有：流入控制容积的总质量 = 流出控制容积的总质量 + 控制容积内原有质量的变化即：m1ΔV + m2ΔV + ... + mkΔV = mk+1ΔV + mk+2ΔV + ... + mnΔV + Δm其中，mk+1到mn表示从其他网格单元流入控制容积p的网格单元，m1到mk表示从控制容积p流出的网格单元。
由于我们只关心内节点p，所以可以简化为：m1ΔV + m2ΔV +... + mkΔV = Δm进一步简化，得到：(m1 + m2 + ... + mk)ΔV = Δm由于Δm = ΔV (ρ1Δx,t + ρ2Δx,t + ... + ρkΔx,t)，其中ρi表示第i个网格单元的密度，Δx表示网格大小。
代入上式，得到：(m1 + m2 + ... + mk) ΔV = ΔV (ρ1Δx,t + ρ2Δx,t + ... + ρkΔx,t)进一步整理，得到：(m1 + m2 + ... + mk) = (ρ1Δx,t + ρ2Δx,t + ...+ ρkΔx,t)这就是关于内节点p的显式离散化方程。
有限容积法和有限体积法
有限容积法和有限体积法有限容积法和有限体积法是计算流体力学中常用的两种数值方法，它们在流体动力学的数值计算中占有非常重要的地位。
本文将从概念、原理、特点、应用等方面，对这两种方法进行详细介绍。
一、有限容积法1.概念有限容积法（Finite Volume Method，FVM）是一种离散化的数值方法，它将连续的物理量离散化为有限个体积元，在每个体积元内计算其平均值，进而求解整个流体系统的物理量。
FVM方法的核心是质量守恒原理，即物质的进出必须平衡，这种保证了物理量在每个体积元内的守恒关系，从而保证了数值计算的准确性。
2.原理FVM方法的数值计算是基于网格的，它将流体动力学问题离散化为一个由有限体积元组成的系统，将原问题转化为流量守恒方程的求解，即$$\frac{\Delta m}{\Delta t}=\Sigma_{faces}\rho uA$$其中，$\Delta m$是在$\Delta t$时间内通过一个表面的质量变化量，$\rho$是介质的密度，$u$是速度，$A$是面积。
对于每个有限体积元，上式可以写为其中，$F_{ij}^p$和$F_{ij}^n$分别是流向有限体积元内部和外部的通量，$i,j$是有限体积元的编号。
3.特点（1）FVM方法基于质量守恒原理，具有非常强的数值稳定性和保真性；（2）FVM方法的计算结果具有局部守恒性，能够准确反映流场内部的物理现象；（3）FVM方法可以处理非结构化网格，适用范围广泛；（4）FVM方法求解的是面积分，所需的时间和空间存储相对较少。
4.应用（1）流体力学领域，如空气动力学、水力学、燃烧问题等；（2）材料科学领域，如薄膜生长、材料变形等。
有限体积法（Finite Element Method，FEM）是一种离散化的数值方法，它将求解的物理场离散化为有限个单元，然后在每个单元内进行近似计算。
相比于FVM方法，FEM方法更加精确，适用于需要高精度计算的问题。
CFD仿真模拟【转】
一．介绍1.什么是CFD?简单地说，CFD就是利用计算机求解流体流动的各种守恒控制偏微分方程组的技术，这其中将涉及流体力学(尤其是湍流力学)、计算方法乃至计算机图形处理等技术。
因问题的不同，CFD技术也会有所差别，如可压缩气体的亚音速流动、不可压缩气体的低速流动等。
对于暖通空调领域内的流动问题，多为低速流动，流速在10m/s以下；流体温度或密度变化不大，故可将其看作不可压缩流动，不必考虑可压缩流体高速流动下的激波等复杂现象。
从此角度而言，此应用范围内的CFD和数值传热学NHT(Numerical Heat Transfer)等同。
另外，暖通空调领域内的流体流动多为湍流流动，这又给解决实际问题带来很大的困难。
由于湍流现象至今没有完全得到解决，目前HVAC内的一些湍流现象主要依靠湍流半经验理论来解决。
总体而言，CFD通常包含如下几个主要环节：建立数学物理模型、数值算法求解、结果可视化。
2．建立模型建立数学物理模型是对所研究的流动问题进行数学描述，对于暖通空调工程领域的流动问题而言，通常是不可压流体的粘性流体流动的控制微分方程。
另外，由于暖通空调领域的流体流动基本为湍流流动，所以要结合湍流模型才能构成对所关心问题的完整描述，便于数值求解。
如下式为粘性流体流动的通用控制微分方程，随着其中的变量f的不同，如f代表速度、焓以及湍流参数等物理量时，上式代表流体流动的动量守恒方程、能量守恒方程以及湍流动能和湍流动能耗散率方程。
基于该方程，即可求解工程中关心的流场速度、温度、浓度等物理量分布。
3．数值算法上述的各微分方程相互耦合，具有很强的非线性特征，目前只能利用数值方法进行求解。
这就需要对实际问题的求解区域进行离散。
数值方法中常用的离散形式有：有限容积，有限差分，有限元。
目前这三种方法在暖通空调工程领域的CFD技术中均有应用。
总体而言，对于暖通空调领域中的低速，不可压流动和传热问题，采用有限容积法进行离散的情形较多。
有限容积法
有限容积法有限容积法（Finite Volume Method）又称为控制体积法。
其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。
其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。
从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。
简言之，子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解，并能得出直接的物理解释。
离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。
限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。
这是有限体积法吸引人的优点。
有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。
就离散方法而言，有限体积法可视作有限单元法和有限差分法的中间物。
有限单元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。
有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。
在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。
有限容积法（FVM）是计算流体力学（CFD）和计算传热学（NHT）中应用最广泛的数值离散方法。
它通常包括如下五个部分：1.? ? ? ? 网格生成2.? ? ? ? 对流项的离散化3.? ? ? ? 边界条件的离散化4.? ? ? ? 压力速度耦合5.? ? ? ? 离散方程的求解对以上五个部分的处理将直接影响到最准结果的SIMPLE算法自1972年问世以来在世界各国计算流体力学及计算传热学界得到了广泛的应用，这种算法提出不久很快就成为计算不可压流场的主要方法，随后这一算法以及其后的各种改进方案成功的推广到可压缩流场计算中，已成为一种可以计算任何流速的流动的数值方法。
CFD的最新发展
计算流体力学(Computational Fluid Dynamics，简称CFD)是当代迅速发展的一门学科，是利用高速计算机求解流体流动的偏微分方程组，目的是为了更好的从定性上和定量上了解流体流动的物理现象，改进设计的一门学科。
目前在航空航天、造船、气象、海洋、水利、液压和石油化工等工程领域都有广泛的应用。
作为流体力学和发动机设计的新手段，CFD是一种令人鼓舞的模拟流体流动的方法，它大大缩短了设计的时间，节省了设计费用。
它相对于理论方法来说，具有假设限制少、应用范围广的特点，其方法也容易应用。
相对于实验来说，计算流体力学很少有马赫数和物体尺寸的限制，并且具有较高的经济价值。
数值仿真优于实验的地方还在于：计算机仿真的诊断“探测”并不干扰流动且不使所研究的现象变得不可捉摸。
CFD已经代替了许多环境发动机的试验项目，而试验的目的也逐渐从验证设计参数的合理性，改变为对CFD数值模拟的正确性及最终设计的校核。
CFD不仅可以为固体环境发动机提供快速而经济的设计依据，并且可以观测到一些试验中无法观测到的物理现象，还可以为新型发动机的设计提供理论依据。
CFD技术艰深的理论背景与流体力学问题的复杂多变阻碍了它向工业界推广。
一般工程技术人员很难较深入地了解这门学科，由专家编制的程序用起来也不容易，因为总有不同条件、参数要根据具体问题以及运算过程随时做出修改调整，若不熟悉广法和程序，往往会束手无策，此外，前、后处理也显得十分棘手。
CFD研究成果与实际应用的结合成为极大难题，这一切曾使人们对CFD的工程应用前景产生疑虑。
在此情况下，通用软件包应运而生，使CFD计算变得方便、简单。
CFD软件一般包括三个主要部分：前处理器（建模，网格生成等），解算器（具体的数值运算）和后处理器（运算结果的具体演示）。
常见的CFD软件有：FLUENT，PHOENICS，CFX，STAR-CD，FIDAP等。
以FLUENT公司开发的大型CFD软件FLUENT为例，它可计算从不可压缩（低亚音速）到轻度可压缩（跨音速）直达高度可压缩（超音速）流体的复杂流动问题。
第4讲 控制方程离散方法简介
1 V
dAx
A
1 V
N
i Aix
i1
在 x 方向
y
1 V
V y
dV
1 V
dAy
A
1 V
N
i Aiy
i1
在 y 方向
13
离散化例子I
在结构化均匀网格中，利用有限体积法对二维连续性方程离散化。
连续性方程
u v 0 x y
14
离散化例子II
在x方向
w W
N
1 u dV 1
V V x
V
控制方程的离散方法
7、数值积分变换法(Integral transformation method，ITM) 对不具备分析解的非线性偏微分方程，把它的解表示成一
个特征值问题的解及一个降维的定解问题解的组合 优点是计算时精度可以较高，但不容易形成通用程序，特
征值问题的选取有任意性，对强非线性问题，计算量较大。
第四讲 控制方程离散方法简介
屠基元 教授 清华大学
墨尔本皇家理工大学
三个关键环节
• 建立恰当的数学模型 Proper Mathematical Modelling
• 对求解区域进行离散化处理 Discretization of Computational Domain
• 对数学模型进行离散化处理 Discretization of Mathematical Model
(x)
i + 1
i i – 1
中心差分 A
C 精确解
向前差分 B
向后差分
i–1
i
i+1
x
有限差分的图像表示
11
有限体积法I
控制体体积
计算流体动力学(CFD)简介
图3-2 各软件之间的关系图
2.3.1 Fluent软件包的安装
Fluent的安装顺序如下： (1) 安装Exceed。推荐安装Exceed6.2版本。 (2) 安装Gambit。单击Gambit的安装，按照提示就可以完成安装，推 荐安装Gambit2.2.30。 (3) 安装Fluent。单击Fluent安装文件，按照提示就可以完成安装, 推 荐安装Fluent6.2..016。 一般来说，Fluent和Gambit的安装推荐使用默认安装设置。当按照以 上的安装步骤安装完毕以后，还要对Fluent和Gambit的环境变量进行设置 。
图3-6 矩形截面管道示意图
图3-7 流体计算区域示意图
2.4.2 实例分析
当利用Fluent解决某一工程问题时，要详细考虑以下几个问题： (1) 确定计算目标； (2) 选择计算模型； (3) 确定物理模型； (4) 确定解的程序。
2.4.3 实例操作步骤
1. 利用Gambit建立计算区域和指定边界条件类型 步骤1：文件的创建及其求解器的选择 (1) 启动Gambit软件 Gambit设置好环境变量以后，可以选择“开始”→“运行”打开如 图3-8所示的对话框，单击“确定”按钮。接着又会弹出如图3-9所示的 对话框，单击Run按钮可以启动Gambit软件，它的窗口布局如图3-10所 示。
根据控制方程离散方式，分为 有限差分法（FDM） 有限元法（FEM） 有限分析法（FAM） 有限体积法或者控制体积法（FVM或CVM）。 有限体积法导出的离散方程可以保证守恒特性， 而且离散方程的系数物理意义明确，是目前计算 流体力学中应用最广的一种方法。
优势 1.可得流动问题满足工程需要的数值解 2.可利用计算机进行各种数值试验 局限性 1.是一种离散近似算法 2.需充分了解所求解问题 3.程序编制、正确使用等要求较高
容积法
有限容积法有限容积法（Finite Volume Method）又称为控制体积法。
其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。
其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。
从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。
简言之，子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解，并能得出直接的物理解释。
离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。
限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。
这是有限体积法吸引人的优点。
有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。
就离散方法而言，有限体积法可视作有限单元法和有限差分法的中间物。
有限单元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。
有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。
在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。
有限容积法（FVM）是计算流体力学（CFD）和计算传热学（NHT）中应用最广泛的数值离散方法。
它通常包括如下五个部分：1.? ? ? ? 网格生成2.? ? ? ? 对流项的离散化3.? ? ? ? 边界条件的离散化4.? ? ? ? 压力速度耦合5.? ? ? ? 离散方程的求解对以上五个部分的处理将直接影响到最准结果的SIMPLE算法自1972年问世以来在世界各国计算流体力学及计算传热学界得到了广泛的应用，这种算法提出不久很快就成为计算不可压流场的主要方法，随后这一算法以及其后的各种改进方案成功的推广到可压缩流场计算中，已成为一种可以计算任何流速的流动的数值方法。
有限容积法
有限容积法有限容积法是一种能够精确地求解复杂热力学系统的分子模拟方法。
此法的基本思想是采用一定的容积，在一定的温度和压力下将分子进行势能最小化，以达到复杂物质以最小体积放到容积中的效果，并讨论及分析其结构及性质，从而实现对复杂系统的模拟。
有限容积法利用模拟技术建立热力学系统，引入热力学模型及计算方法，将历史上实验测量众多的温度及压力值进行分析，给出热力学模型，并将模型量化，建立模拟程序。
在建立模拟程序的过程中进行能量最小化，将原有的容积空间按照一定的程序处理，即改变容积，使其朝向更低的能量，即使容积很小，但是能量也不会增加，当整个容积的能量达到最小时，模型程序就可以终止，这时候就得到了系统最终的结构。
有限容积法可以让实验者将系统放在一定体积中，然后在其中寻找最有效的结构。
有限容积法得到的系统结构极为准确，因此成为复杂热力学系统的重要分析方法，如化学传输流动泵、有序多孔材料、有序分子分子装配，包括复杂结构的单分子层等。
此外，有限容积法还可以应用于药物分子的设计与优化，预测新材料的性能等领域。
有限容积法不仅可以提高计算效率，而且可以避免模拟结果的偶然性。
其应用前景十分广阔。
尽管有限容积法能够实现准确的模拟系统，但是也存在一些问题。
首先，有限容积法需要较高的计算能力。
其次，有限容积法的分子模拟容积有限，对于过度的模拟会受到一定的限制；此外，每个步长的时间要求仍然比较长，会影响计算的运行效率。
最后，实际的计算模型会受到外部的环境的影响，当外部环境发生变化时，这些模型可能会出现误差。
由于有限容积法具有上述缺点，在计算复杂热力学系统时仍需注重加以考虑。
计算流体力学中常用的控制方程离散化方法概述
计算流体力学中常用的控制方程离散化方法概述计算流体力学是现代流体力学的一种数值计算方法，最早出现是在20世纪50年代。
它主要应用于流体的流动、传热、化学反应、物质转移等方面的数值计算，成为了工程和科学界不可或缺的工具。
计算流体力学中的控制方程离散化方法则是其中重要的一部分，本文将就此进行概述。
一、控制方程离散化在计算流体力学中，控制方程是解决问题的基础，主要包括连续性方程、动量方程和能量方程等。
这些方程通过离散化方法进行处理，变成可以计算机进行处理的数学模型。
离散化的基本思想是将时间和空间分成有限个点来处理，利用简单的数值运算方法计算每个时间步长中的各个物理量。
常用的离散化方法包括有限差分方法、有限体积方法、有限元方法等。
二、有限差分方法有限差分方法是计算流体力学中常用的一种离散化方法，它是一种基于差分的数值方法，利用有限差分近似代替微分方程，求解微分方程数值解的方法。
它的主要思想是将一个连续的空间域区间划分为一些点，对连续波动函数的任意一阶导数代替为该点处差分的近似，从而把原问题转化为一个差分方程组，通过解这个方程组来求解微分方程的近似解。
三、有限体积方法有限体积方法是一种对控制方程离散化方法，它是一种基于控制方程积分形式的方法。
该方法基于微积分的思想，通过对空间区域划分成有限的体积单元来进行数值计算。
在有限体积方法中，我们通常选择一个体积单元V，然后从该体积单元周围的表面积进行积分，得到控制方程的离散形式。
四、有限元方法有限元方法是计算流体力学中另一种常用的离散化方法，它能够适应各种复杂流动情况。
该方法可以将连续问题变为离散问题，进而离散化求解成一些小片断组成的离散问题，并且可以在不同的片段上使用不同阶次的多项式进行近似，从而得到更为准确的结果。
在有限元方法中，我们通常需要先对区域进行剖分，然后利用插值法来构造近似解。
五、总结综合来说，计算流体力学中常用的控制方程离散化方法有有限差分方法、有限体积方法和有限元方法三种。
有限容积法
有限容积法有限容积法（Finite Volume Method）又称为控制体积法。
其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。
其中的未知数是网格点上的因变量的数值。
为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。
从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。
简言之，子区域法属于有限体积发的基本方法。
有限体积法的基本思路易于理解，并能得出直接的物理解释。
离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。
限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。
这是有限体积法吸引人的优点。
有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。
就离散方法而言，有限体积法可视作有限单元法和有限差分法的中间物。
有限单元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。
有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。
有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。
在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。
有限容积法（FVM）是计算流体力学（CFD）和计算传热学（NHT）中应用最广泛的数值离散方法。
它通常包括如下五个部分：1. 网格生成2. 对流项的离散化3. 边界条件的离散化4. 压力速度耦合5. 离散方程的求解对以上五个部分的处理将直接影响到最准结果的SIMPLE算法自1972年问世以来在世界各国计算流体力学及计算传热学界得到了广泛的应用，这种算法提出不久很快就成为计算不可压流场的主要方法，随后这一算法以及其后的各种改进方案成功的推广到可压缩流场计算中，已成为一种可以计算任何流速的流动的数值方法。
3.3 微分方程的有限容积法离散
常用型线选择 – 空间
阶梯分布
分段线性分布
常用型线选择 – 时间
阶梯分布
分段线性分布
例：一维有源对流扩散方程的离散
( u ) S t x x x
对 流 项：需对时间分布取假设
a) 阶梯显式：
b) 阶梯隐式：
c) 取随 t 分段线性分布
算术平均格式 Crank-Nicolson 格式
扩散项：需对时间分布取假设 Nhomakorabeaa) 阶梯显式：
b) 阶梯全隐式：
c) 取随 t 分段线性分布
算术平均格式 Crank-Nicolson 格式
源项：需对时间、空间分布取假设
1) 积分：
t t e
t
dxdt t w
t t e
t
u dxdt x w
t t e
t
dxdt x x w
t t e
Sdxdt
t w
2) 型线假设
• 非稳态项：
需对空间做假设，取阶梯分布
空间取阶梯分布，而时间可取： a) 阶梯显式： b) 阶梯全隐式：
c) 取随 t 分段线性分布
算术平均格式 Crank-Nicolson 格式
界面上的函数值和导数值
阶梯分布
分段线性分布
界面上的函数值和导数值
离散格式：各种情况的组合
• 显式
• 全隐式
Crank-Nicolson 格式
• 均匀网格，不可压
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数值传热学基本思想
有限差分法
一维稳态有源项的对流-扩散方程
泰勒展开
得到线性方程组
有限容积法
在一个控制体内做积分
积分之后对方程进行离散
有限容积法
扩散项线性化处理（中心差分）
源项认为在整个控制体内不变，局部线性化处理
两种型线假设
有限容积法
求解线性方程组
有限容积法
有限容积法
有限容积法
在均分网格情况下与Taylor展开法的结果一致。 在FVM中所谓不同的格式就是指不同的型线。
求 扩 散 通 量
二次插值：有问题
二次插值格式QUICK 守恒
2、有界性Boundedness
没有源项时内部节点的参数值应该位于边界节点的范围限制内
迭代收敛的充分条件：
系数矩阵对角占优
1、源项的线性化系数应该为负 若源项的线性化系数为负：T增大→S增大→T增大→S增大：不稳定 2、离散方程里的所有系数应该有相同的符号（通常为正）。一个节点参数的增加 应该导致相邻节点参数的增加。 中心节点的参数值为相邻节点的加权平均。 如果离散格式不满足局限性，那么解可能不收敛，或者有“摆动” 3、没有源项时，中心点系数为相邻点系数之和。 当控制方程只包含微分项时，T和T+C都满足微分方程。故参 数加上常数时离散方程仍然成立。
CFD控制方程离散方法：有限容积法
控制方程
连续性方程
质量流动
动量方程
动量流动 表面力
能量方程
热流通量
RANS动量方程
湍流动能k 方程
是k的扩散系数, G 是湍流动能的产生速率， 是耗散速率 , 是方程的净源项。
一般形式
守恒形式与非守恒形式
守恒形式
随体导数
连续性方程
非守恒形式 CFD有限容积法计算中用守恒形式，只有守恒型的方程才能 保证有限大小体积内守恒定律的成立
举例
节点1： 节点5：
小于1在边界处达到，因为ap包含所有边界点（包括值 已知）的边界点，已知值的边界点作为源项出现，使ap 更大。
3、输运性Transportiveness
邻点W和E有两个恒定源，画出等值线
纯扩散
对流扩散
纯扩散使源的影响向各个方向同等地传播；纯对流时，P点只受上游影响不受下游影响。
有限容积法
离散与数值解
方程非线性
各个变量通过系数等相互影响：内耦合 把原来在空间与时间坐标中连续的物理量的场，用一系列有限个离散点上 的值的集合来代替，建立起离散方程，求解所建立起来的代数方程以获得所求 解变量的近似值。 分离式求解：在一个方程中认为系数（包括速度）为常数，采用上一次迭代 算出来的值。 迭代求解线性方程组，从一个初始的（猜测的）变量分布开始，不断更新直至 收敛。
无源时场随时间变化
3、输运性Transportiveness
n时刻 n+1时刻 扰动被均匀向两侧传递 中心差分
对流项的中心差分不合理，因为aE为负，使得下游增大会使上游 减小
迎风算法
扰动不向上游传递
一阶迎风的数值误差较大，在调试阶段可采用。
分离式求解过程
初始
系数
解方程
系数
解方程 ……
总结
离散格式的要求
不合理的解
控制方程离散后求解能够收敛但是得到物理上不合理的解。
1、守恒性Conservativeness
跨过同一个界面离开一个控制体和进入相邻控制体的通量必须相等。相邻控制 体计算跨过界面的通量的式子必须相同（为同一个式子）。
中心差分：没问题