计算流体力学课后题作业

1、实用标准课后习题第一章1 .计算流体动力学的基本任务是什么计算流体动力学是通过计算机数值计算和图像显示， 对包含有流体流 动和热传导等相关物理现象的系统所做的分析。2 .什么叫控制方程？常用的控制方程有哪几个？各用在什么场合？流体流动要受物理守恒定律的支配， 基本的守恒定律包括：质量守恒 定律、动量守恒定律、能量守恒定律。如果流动包含有不同组分的混 合或相互作用，系统还要遵守组分守恒定律。如果流动处于湍流状态， 系统还要遵守附加的湍流输运方程。控制方程是这些守恒定律的数学 描述。常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组 分质量守恒方程。质量守恒方程和动量守恒方程任何流动问题

2、都必须 满足，能量守恒定律是包含有热交换的流动系统必须满足的基本定律。 组分质量守恒方程，在一个特定的系统中，可能存在质的交换，或者 存在多种化学组分，每种组分都需要遵守组分质量守恒定律。4.研究控制方程通用形式的意义何在？请分析控制方程通用形式中 各项的意义。建立控制方程通用形式是为了便于对各控制方程进行分析，并用同一程序对各控制方程进行求解。（夕0）+ 历t （夕日。）= d/t （r grad） + S Cx J各项依次为瞬态项、对流项、扩散项、源项。6.CFD商用软件与用户自行设计的 CFD程序相比，各有何优势？常 用的商用CFD软件有哪些？特点如何？由于CFD的复杂性及计算机软硬件条

3、件的多样性，用户各自的应用程序往往缺乏通用性。CFD商用软件的特点是功能比较全面、适用性强。具有比较易用的前后处理系统和其他 CAD及CFD软件的接口能力, 便于用户快速完成造型、网格划分等工作。具有比较完备的容错机制和操作界面，稳定性高。可在多种计算机、多种操作系统，包括并行环境下运行。常用的商用 CFD 软件有 PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。PHOENICS除了通用CFD软件应该拥有的功能外，PHOENICS软件有自己独特的功能：开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。CFX除了可以使用有限体积法外，还采用基于有限元的

4、有限体积法。用于模拟流体流动、传热、多相流、化学反应、燃烧问题。其优势在于处理流动物理现象简单而几何形状复杂的问题。SRAR-CD基于有限体积法，适用于不可压流体和可压流的计算、热 力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器 三大模块，以良好的可视化用户界面把建模、 求解及后处理与全部的 物理模型和算法结合在一个软件包中。FIDAP完全基于有限元方法。可用于求解聚合物、薄膜涂镀、生物医 学、半导体晶体生长其他领域中出现的各种层流和湍流的问题。对涉及流体流动、传热、传质、离散相流动、自由表面液固相面等问题都 提供了精确而有效的解决方案。FLUENT是一个用于模拟和分析在复杂几何

5、区域内的流体流动与热 交换问题的专用CFD软件。网格划分灵活，FLUENT使用C语言开 发完成，可实现动态内存分配及高效数据结构， 具有很大的灵活性与 很强的处理能力。FLUENT中，解的计算与显示可以通过交互式的用 户界面来完成。第二章1.什么叫离散化？意义是什么？离散化，即对空间上连续的计算区域进行划分，把它划分成多个子区 域，并确定每个区域中的节点，从而生成网格。然后，将控制方程在 网格上离散，即将偏微分格式的控制方程转化为各个节点上的代数方 程组。对于在求解域内所建立的偏微分方程， 理论上是有真解（或称精确解 或解析解）的。但由于所处理的问题自身的复杂性，一般很难获得方 程的真解。因此

6、，就需要通过数值方法把计算域内有限数量位置（网 格节点或网格中心点）上的因变量值当作基本未知量来处理，从而建 立一组关于这些未知量的代数方程组，然后通过求解代数方程组来得 到这些节点值，而计算域内其他位置上的值则根据节点位置上的值来 确定。这样，用变量的离散分布近似解代替了定解问题精确解的连续 数据，当网格节点很密时，离散方程的解将趋近于相应微分方程的精 确解。3 .简述有限体积法的基本思想，说明其使用的网格有何特点？有限体积法的基本思想，将计算区域划分为网格，并使每个网格点周 围有一个互不重复的控制体积；将待解微分方程对每个控制体积积分, 从而得出一组离散方程。节点排列有序，即当给出了一个节

7、点的编号后，立即可以得出其相邻 节点的编号。称为结构网格，是一种传统的网格形式，网格自身利用 了几何体的规则形状。非结构网格的节点以一种不规则的方式布置在 流场中，网格生成复杂，但却有着极大地适应性，尤其对具有复杂便 捷的流场计算问题特别有效。4 .简述在空间域上离散控制方程的基本做法， 说明对流项及扩散项在 离散处理时的异同，给出所生成的二维稳态对流 -扩散问题的离散方 程的形式。生成计算网格，包括节点、控制体积；将守恒型控制方程在每个控制 体积上作积分，得到关于节点未知量的代数方程组；求解代数方程组, 得到个计算节点的未知量值。对流项与扩散项而言，则通过采用控制节点上的值选取合适的插值函数

8、来表示，其中扩散项通常采用中心差分格式来表示， 而对流项的形式则是多种多样的。二维稳态对流-扩散问题的离散方程6.分析比较中心差分格式、一阶迎风格式、混合格式、指数格式、二 阶迎风格式、QUICK格式各自的特点及适用场合。中心差分格式(central differencing scheme )：就是界面上的物理 量采用线性插值公式来计算，即取上游和下游节点的算术平均值。 它 是条件稳定的，在网格Pe数小于等于2时，中心差分格式的计算结 果与精确解基本吻合，在不发生振荡的参数范围内，可以获得较准确 的结果。如没有特殊声明，扩散项总是采用中心差分格式来进行离散。 但中心差分格式因为有限制而不能作为

9、对于一般流动问题的离散格 式，必须创建其他更合适的离散格式。一阶迎风格式(first order upwind scheme):即界面上的未知量恒取上游节点(即迎风侧节点)的值。这种迎风格式具有一阶截差，因 此叫一阶迎风格式。无论在任何计算条件下都不会引起解的振荡，是绝对稳定的。但是当网格Pe数较大时，假扩散严重，为避免此问题， 常需要加密网格。研究表明，在对流项中心差分的数值解不出现振荡 的参数范围内，在相同的网格节点数条件下，采用中心差分的计算结 果要比采用一阶迎风格式的结果误差小。因此，随着计算机处理能力 的提高，在正式计算时，一阶迎风格式目前常被后续要讨论的二阶迎 风格式或其他高阶格式

10、所代替。混合格式(hybrid scheme )：综合了中心差分和迎风作用两方面的 因素，当|Pe|<2时，使用具有二阶精度的中心差分格式；当|Pe|>2 时，采用具有一阶精度但考虑流动方向的一阶迎风格式。该格式综合了中心差分格式和一阶迎风格式的共同的优点，其离散系数总是正的，是无条件稳定的。计算效率高，总能产生物理上比较真实的解，且是 高度稳定的。但缺点是只具有一阶精度。指数格式(exponential scheme ):将扩散与对流的作用合在一起来 考虑，绝对稳定。在应对于一维的稳态问题时，指数格式保证对任何 的Pelclet数以及任意数量的网络点均可以得到精确解。缺点是指数

11、运算较为费时，对于多维问题以及源项不为零的情况此方案不准确。 乘方格式(power-law scheme ):绝对稳定，与指数格式的精度较 接近，但比指数格式省时。主要适用于无源项的对流-扩散问题。对 有非常数源项的场合，当Pe数较高时有较大误差。二阶迎风格式(second order upwind scheme ): 二阶迎风格式与 一阶迎风格式的相同点在于，二者都通过上游单元节点的物理量来确 定控制体积界面的物理量。但二阶格式不仅要用到上游最近一个节点 的值，还有用到另一个上游节点的值。它可以看作是在一阶迎风格式 的基础上，考虑了物理量在节点间分布曲线的曲率影响。在二阶迎风格式中，实际上只

12、有对流项采用了二阶迎风格式，而扩散项仍采用中心差分格式。二阶迎风格式具有二阶精度的截差，但仍有假扩散的问题。QUICK格式：是“对流项的二次迎风插值”，是一种改进离散方程截 差的方法，通过提高界面上插值函数的阶数来提高格式截断误差的。 对流项的QUICK格式具有三阶精度的截差，但扩散项仍采用二阶截 差的中心差分格式，QUICK格式具有守恒特性。对于与流动方向对 齐的结构网格而言，QUICK格式将可产生比二阶迎风格式等更精确 的计算结果。QUICK格式常用于六面体（二维中四边形）网格。对 于其它类型的网格，一般使用二阶迎风格式。第三章2.可压流动与不可压流动，在数值解法上各有何特点，为何不可压流

13、 动在求解时反而比可压流动有更多的困难？如果流动是可压的，可以把密度视作连续方程中的独立变量进行求解， 即以连续方程作为一个普通的关于密度的输运方程生成相对简单的 离散方程组，压力根据气体状态方程得到。对于不可压流动，密度是 常数，就不可能把密度与压力相联系，因此将密度作为基本未知量的 方法不可行。7. SIMPLE算法的基本思想是什么？动量方程和连续性方程在其中 是如何得到满足的？在交错网格上如何实施 SIMPLE算法？ SIMPLE算法的基本思想是对于给定的压力场（它可以是假定的值， 或是上一次迭代计算所得到的结果），求解离散格式的动量方程，得 出速度场。因为压力场是假定的或不精确的，这样

14、得到的速度场一般不满足连续 方程，因此，必须对给定的压力场加以修正。修正的原则是：与修正 后的压力场相对应的速度场能满足这一迭代层次上的连续方程离散 形式。据此原则，我们把由动量方程的离散形式所规定的压力与速度 的关系带入连续方程的离散格式，从而得到压力修正方程，由压力修 正方程得出压力修正值。接着，根据修正后的压力场，求的新的速度 场。然后检查速度场是否收敛。若不收敛，用修正后的压力值作为给 定的压力场，开始下一层次的计算。直至获得收敛的解。SIMPLE算法计算步骤1 .假定一个速度分布，记为u0, v0,用于计算离散方程系数2 .假定一个压力场p*3 .依次求解两个动量方程，得到u* ,

15、v*4 .求解压力修正方程，得到p'5 .据p '得速度修正方程改进速度值6 .利用改进后的速度场求解那些通过与动量方程耦合的其它变量方 程，如果其它变量方程不影响动量方程， 那么应该在速度场收敛后再 求解7 .利用改进后的速度场和改进后的压力场作为初值进行下一轮的迭 代，返回3.8 . SIMPLEC 算法与SIMPLE算法相比有何改进？效果在哪里 ？请给 出SIMPLEC算法的具体实施条件和步骤。 一、一一 ' 一 ' SIMPLE算法中速度修正值万程略去I：?E?%?使速度修正完全归结为压差项直接作用SIMPLEC算法没有像SIMPLE算法那样将汇？?荻忽

16、略，因此得到 的压力修正值p' 一般是比较合适的，因此，在SIMPLEC算法中可不 再对p'进行欠松弛处理。与SIMPLE算法基本相同SIMPLEC算法计算步骤1 .假定一个速度分布，记为u0, v0,用于计算离散方程系数2 .假定一个压力场p*3 .依次求解两个动量方程，得到u* , v*4 .求解压力修正方程，得到p'5 .据p '得速度修正方程改进速度值6 .利用改进后的速度场求解那些通过与动量方程耦合的其它变量方程，如果其它变量方程不影响动量方程， 那么应该在速度场收敛后再 求解7 .利用改进后的速度场和改进后的压力场作为初值进行下一轮的迭代，返回3.10. SIMPLER和PISO的主要特点是什么？与 SIMPLE的区别在哪里？SIMPLE算法的各种改进算法，主要是提高了计算的收敛性，从而缩 短计算时间。改进后的SIMPLER算法只用压力修正值p'来修正速度，另外构建 一个更加有效的压力方程来产生正确的压力场。由于在推导 SIMPLER算法的离散化压力方程时，没有任何项被省略，因此所得到的压力场和速度场相适应。在 SIMPLER算法中，正确的速度场将 导致正确的压力场，而在SIMPLE算法中则不是这样。所以SIMPLER 算法是在很高的效率下正确计算压力场的。SIMPLER算法的计算量比SIMPLE算法高出3