计算流体力学--典型算法与算例(附光盘21世纪高等教育建筑环境与能源应用工程系列规划

高歌等编著的《计算流体力学--典型算法与算例 (附光盘)》是按照理工科院校少学时计算流体力学课 程要求编写的一部本科生通用教材。全书分10章，主 要内容包括：计算流体力学的发展简况及基本方程组 、不可压缩流动的典型数值方法、可压缩流动对流项 与时间项处理的典型方法、复杂流动问题的典型处理 方法及Gao-Yong湍流方程组、计算流体力学中网格生 成的典型方法。 本书内容重基础、重实用、少而精，反映了当前 计算流体力学的最新科研成果，深入浅出，还融入了 Gao-Yong理性湍流理论与方程组的相关内容。书中正 文给出了大量的典型算例，并且在第2、3、4篇末还 提供了习题(习题解答见本书所附光盘)，这有利于广 大学生与读者加深对基本格式的理解与认识以及全面 掌握CFD方法。 本书可作为高等院校理工科类计算流体力学课程 的本科生教材，也可作为建筑环境与能源应用工程、 航空航天、动力、机械、水利工程、生物、化工、环 境工程、兵器科学与技术、工程热物理及工程力学等 专业研究生的教材或参考书。

李凤超，1983年出生于黑龙江省哈尔滨市，2006年本科毕业于清华大学热能工程系，2011年获清华大学动力工程及工程热物理专业博士学位。现就职于中国航空研究院新技术研究所，从事航空叶轮机械气动热力性能模拟与优化设计方向研究工作，高级工程师。 吴俊宏，1983年9月生于浙江义乌，2001年考入中国农业大学工学院热能与动力工程专业学习，2010年毕业于北京理工大学宇航学院，获得流体力学博士学位，并进入中国航空研究院工作至今。曾获得过北京市第十九届非物理类专业大学生物理竞赛甲组三等奖、优秀毕业生等荣誉。在攻读博士期间主要从事高精度、高分辨率数值方法的研究，参与过多项国家自然科学基金等项目，发表论文10余篇。目前从事高精度、高分辨率数值格式、叶轮机械气动流场计算以及Gao-Yong湍流方程等研究工作，高级工程师。 王保国，生于山东省临清市，2007年获“北京市教学名师”荣誉称号，现任中国人类工效学学会副理事长、人机工程专业委员会主任。曾在中国科学院力学所和工程热物理所工作16年，后曾分别在清华大学工程力学系和北京理工大学机械与车辆学院、机电学院、宇航学院任教授、博士生导师。另外，还曾在两个国家级航空研究单位从事工程设计与实验工作10年。发表学术论文230余篇，出版专著与国家级规划教材12部。曾两次获国际奖，多次获国家级与省部级奖励。 闫文辉，1979年9月出生于吉林省辉南县，北京航空航天大学博士，中国航空研究院高级工程师。2003年师从北航高歌教授，从事湍流理论及数值模拟研究，2007—2008年，在美国佛罗里达大西洋大学机械系访问交流。现从事计算流体力学及新型动力装置等研究工作。 高歌，1945年1月出生，祖籍山东沂南，北京航空航天大学教授，工程热物理学家，Gao-Yong理性湍流理论及广义宇宙演化论创始人。长期从事流体力学和工程热物理领域的教学和科研工作，曾获国家发明一等奖、首届航空金奖、首届光华特等奖、北京市劳动模范等奖项。

前言 第一篇 计算流体力学的发展简况及基本方程组 第1章 计算流体力学的发展历程及流动控制基本方程组 1.1 计算流体力学发展的主要历程 1.1.1 不可压缩黏性流动计算方法的概述 1.1.2 可压缩流动计算方法的主要发展历程 1.2 流动控制的基本方程组 1.2.1 微分形式的NavierStokes方程组 1.2.2 积分形式的NavierStokes方程组 1.2.3 直角笛卡儿坐标系下的NavierStokes方程组 1.2.4 任意曲线坐标系下的NavierStokes方程组 第二篇 不可压缩流动的典型数值方法 第2章 不可压缩黏性流动的涡量流函数高效算法 2.1 涡量流函数方法的控制方程组 2.2 控制方程的有限差分离散 2.2.1 涡量方程的离散 2.2.2 流函数方程的离散 2.2.3 压力方程的离散 2.3 边界条件的处理 2.3.1 进口边界 2.3.2 出口边界 2.3.3 固体壁面边界 2.4 涡量流函数方法的求解步骤 2.5 涡量流函数方法计算实例 2.5.1 二维方腔顶盖驱动流动模拟 2.5.2 二维外掠后台阶流动模拟 2.5.3 二维方柱绕流模拟 2.5.4 二维层流边界层流动模拟 第3章 一种求解不可压缩黏性流动的SIMPLE通用算法 3.1 交错网格上控制方程的离散 3.1.1 压力波动问题与交错网格 3.1.2 控制方程的有限体积离散 3.2 压力修正技术 3.2.1 压力修正的基本思路和步骤 3.2.2 速度修正方程和压力修正方程 3.2.3 边界条件的处理 3.3 SIMPLE算法的求解步骤 3.4 改进的SIMPLE算法 3.4.1 SIMPLER算法 3.4.2 SIMPLEC算法 3.4.3 PISO算法 3.5 SIMPLE系列算法计算实例 3.5.1 平面Poiseuille流动模拟 3.5.2 二维非定常方柱绕流模拟 3.5.3 SIMPLE系列算法的比较 习题 第三篇 可压缩流动对流项与时间项处理的典型方法 第4章 可压缩流动对流项数值格式的几种处理方法 4.1 模型方程及其数学性质 4.1.1 线性单行波方程 4.1.2 热传导方程 4.1.3 无黏Burgers方程 4.1.4 Burgers方程 4.1.5 KdVBurgers方程 4.1.6 一维Euler方程组 4.2 几个经典的差分格式 4.2.1 基本差分格式 4.2.2 LaxFriedrichs格式 4.2.3 蛙跳格式 4.2.4 LaxWendroff格式 4.2.5 人工黏性方法 4.3 矢通量分裂方法 4.3.1 StegerWarming矢通量分裂方法 4.3.2 Van Leer矢通量分裂方法 4.4 Roe格式 4.5 Godunov间断方法 4.5.1 Godunov格式 4.5.2 MUSCL格式 4.5.3 PPM格式 4.6 TVD格式 4.6.1 Harten的二阶TVD格式 4.6.2 带通量限制因子的二阶TVD格式 4.6.3 YeeRoeDavis对称性型二阶TVD格式 4.7 ENO/WENO格式 4.8 间断Galerkin有限元方法 4.9 数值算例 4.9.1 精度分析 4.9.2 一维线性标量方程数值实验 4.9.3 一维非线性标量方程数值实验 4.9.4 一维气体动力学方程组数值实验 4.9.5 二维气体动力学方程组数值实验 第5章 可压缩流黏性项处理的高阶格式 5.1 物理坐标系下的NS方程 5.2 NS方程组的无量纲化 5.3 计算坐标系下的NS方程 5.4 薄层黏性假设 5.5 黏性项的离散格式 5.6 典型算例分析 5.6.1 激波边界层相互作用算例 5.6.2 后台阶算例 5.6.3 喷管射流算例 第6章 可压缩流时间项处理的几种常用方法 6.1 近似因子分解方法 6.2 显式RungeKutta多步格式 6.3 双时间步方法 6.4 RKDG有限元方法 6.5 典型算例 6.5.1 圆柱绕流算例 6.5.2 NACA0012变攻角绕流算例 6.5.3 前台阶算例 6.5.4 双马赫反射算例 习题 第四篇 复杂流动问题的典型处理方法及GaoYong湍流方程组 第7章 三维黏性流动的有限体积解法 7.1 有限体积法概述 7.2 NavierStokes方程组的几种通用形式 7.2.1 直角坐标系下的守恒型基本方程组 7.2.2 任意曲线坐标系下的守恒型基本方程组 7.3 方程组的无量纲化 7.4 空间离散格式 7.4.1 Jameson中心差分格式 7.4.2 通量差分分裂Roe格式 7.5 时间积分格式 7.5.1 RungeKutta显式时间推进法 7.5.2 双时间步的隐式时间离散方法 第8章 常用的几种湍流模式 8.1 湍流模式的发展 8.2 涡黏模式 8.2.1 零方程湍流模式 8.2.2 一方程湍流模式 8.2.3 两方程湍流模式 8.3 雷诺应力输运模式 8.4 非线性湍流模式 第9章 GaoYong湍流方程组及其典型算例 9.1 时间平均、系综平均及Favre平均 9.2 侧偏平均、加权侧偏平均及其重要性质 9.3 GaoYong可压缩湍流方程组 9.3.1 平均流动基本方程组 9.3.2 加权侧偏平均的漂移流方程 9.4 GaoYong不可压缩湍流方程组 9.5 GaoYong湍流方程组计算的典型算例 9.5.1 圆射流/平面射流异常问题 9.5.2 平板固壁边界层及转捩 9.5.3 二维后向台阶分离流动 9.5.4 翼型大攻角分离流动 9.5.5 翼体角偶流动 9.5.6 激波/湍流边界层干扰流动 9.5.7 拟序结构流动 9.6 小结 习题 第五篇 计算流体力学中网格生成的典型方法 第10章 网格生成的几类方法 10.1 结构化网格生成方法 10.1.1 贴体坐标概念 10.1.2 贴体坐标系与直角坐标系间的变换 10.1.3 代数网格生成法 10.1.4 基于偏微分方程的网格生成法 10.1.5 网格质量要求 10.2 非结构化网格生成方法 10.2.1 Delaunay三角化法 10.2.2 阵面推进法 10.3 自适应Descartes直角网格生成技术 附录 光盘中习题解答的详细图目 参考文献