可压不可压缩流体计算详解.pptVIP

为使方程组封闭，还需要引入气体状态方程和内能与焓的关系式。对于完全气体有： 控制方程组及求解方法 或 对于不可压缩流体： 得 展开为 可压流体连续性方程 状态方程不再成立 压力P不再作为热力学变量 连续方程退化为速度的散度为零的方程,在可压缩流动的计算中可用于求解密度和压力的连续方程在不可压缩流动求解中仅是动量方程的一个约束条件，由此求解不可压缩流动的压力成为一个困难。 动量方程和能量方程可以分开求解 不可压和可压方程的差别 不同类型的偏微分方程在特性上主要区别是它们的倚赖域与影响域不同 双曲型：初值、边值问题，有两个实的特征线，依赖域是两个特征线之间的部分，由于倚赖域具有方向性，所以是一个初边值问题； 抛物型：初值、边值问题，因为有一个实的特征线，依赖域取决于特征线前面的区域，从而是初边值问题； 椭圆型：边值问题，需要注意的是倚赖域和影响域是整个流场，设置边条的时候需要给定完整的边条，所以认为是边值问题； 控制方程的数学分类及其对数值解的影响 不含时间变量的不可压缩流方程NS方程组是椭圆型方程，在椭圆型方程中信号传播速度无穷大，即声速无穷大，在边值问题求解中需要全部封闭的边解条件。在实际计算中可看到下游对上游的强烈影响。 含时间变量的不可压缩粘性流NS方程是抛物型方程，但它含有椭圆型方程的强烈次特征。当人们用时间相关法来求解含时间的NS方程组时，在达到稳态解时完全是椭圆型方程问题，方程的边界条件也应按椭圆型方程来提。 不可压流方程的数学性质 原始变量方法 将NS方程写成速度和压力的形式，进行直接求解 非原始变量方法 不直接求解速度、压力的方法 求解不可压缩流动的各种方法主要在于求解不同的压力过程 不可压缩流动求解主要有两类方法 非原始变量法 Fasel提出的流函数——涡函数方法 Aziz和Hellums提出的势函数——涡函数方法 对于二维问题，具有计算上的优势；对于三维题需要求解三个Poisson方程，非常耗时，反而不及原始变量法。 不可压缩流动求解主要有两类方法 原始变量法 Harlow和Welch提出的压力Poisson方程方法 Patanker和Spalding的SIMPLE（Semi-Implicit Method for Pressure-Linked Equation）法 Chorin于1967年提出的虚拟压缩方法 鉴于求解时间，后两种方法使用较多。 不可压缩流动求解主要有两类方法 DNS（直接数值模拟 ） 准确，计算耗费最大，难以应用于工程实际 LES（大涡模拟 ） 较准确，计算耗费较大，工程上用的较少 RANS（雷诺平均的湍流模型法 ） 相对准确，计算耗费小，应用广泛，但引入了湍流模型问题 NS方程的求解方法 雷诺平均 ， 将瞬时压力和密度按时间平均法（雷诺平均）分解为平均值与脉动值之和，即： 其它参数q采用质量加权的时间平均，即： RANS方程形式 对NS方程组进行雷诺平均，得到RANS方程组 连续方程无变化 动量方程粘性应力项多出了一项，类比于粘性应力，称为雷诺应力，引入了封闭问题 雷诺应力 RANS方程形式 工程上计算雷诺应力主要的三种方法：雷诺应力输运模型（RSTM） 、涡粘性模型（EVM） 和代数雷诺应力方程模型（ARSM） 中，应用最多的要算涡粘性模型。根据实际情况选取——综合考虑算例流动种类以及计算资源。 涡粘性模型： 零方程模型——混和长度、C-S 、B-L模型 一方程模型——Splart-Allmaras 模型 两方程模型——k -ε模型、k -ω模型 四方程模型——V2-f模型 没有哪一个湍流模型具有普适性！ 封闭问题以及湍流模型 结构化网格——数值计算方法简单，计算效率高，但复杂结构网格生成困难 非结构化网格——网格生存容易，但数值方法实现复杂，计算效率低，数值精度也容易受到影响 网格生成 拓朴结构主要有H型、C型、O型三类 生成技术一般可分为代数方法、微分方程方法（椭圆型、抛物型、双曲型） 网格质量：贴体性、正交性、光滑性、加密程度（与使用湍流模型相关） 流场参数梯度大的地方网格一定要加密，比如边界层、激波附近 结构化网格 有限差分法——将 求解域划分为差分网格，用有限个网格节点代替连续的求解域。该方法是一种直接将微分问题变为代数 问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法 。主要缺点是对复杂区域的适应性较差及数值解的守恒性难以保证。 有限体积法——计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。有限体积法的基本思路易于理解，并能得出直接的物理解释。具有守恒