信息与计算科学丛书（63）：多介质流体动力学计算方法 pdf epub mobi txt 电子书 下载 2025

简体网页||繁体网页

☆☆☆☆☆

贾祖朋，张树道，蔚喜军 著

图书标签:

• 流体动力学
• 计算方法
• 多介质流
• 数值模拟
• 科学计算
• 信息科学
• 工程应用
• 物理学
• 传热学
• 计算流体力学

出版社： 科学出版社

ISBN：9787030409775

版次：1

商品编码：11497175

包装：精装

丛书名： 信息与计算科学丛书63

开本：16开

出版时间：2014-06-01

用纸：胶版纸

页数：344

字数：433000

正文语种：中文

编辑推荐

适读人群 ：《多介质流体动力学计算方法》可供计算数学、计算流体力学以及工程计算等专业的高年级本科生、研究生及相关研究者阅读参考.
　　《信息与计算科学丛书（63）：多介质流体动力学计算方法》适合从事内爆动力学、界面不稳定性、惯性约束聚变、高速冲击、先进常规兵器、天体物理等领域数值模拟研究与应用的科技工作者阅读使用，也可作为高等院校计算数学和计算流体力学专业研究生和高年级本科生的教材与

内容简介

《多介质流体动力学计算方法》详细地介绍作者及其合作者近年来在流体力学多介质大变形问题数值模拟方法研究方面所取得的较新成果, 其中包括一系列具有理论意义及重要应用价值的数值方法. 例如二维和三维显式有限元相容拉氏方法、基于近似RiemAnn 解的有限体积ALE 方法、一种基于level set的欧拉–拉格朗日耦合方法、一种健壮有效的多介质欧拉方法、二维拉氏双曲守恒律方程组的中心型高精度间断GAlerkin 谱有限元方法、基于MOF 界面重构的二维MMALE 方法、三维非结构网格基于MOF 界面重构的MMALE 方法、三维非结构网格基于浓度法的中心型MMALE 方法、几种新的二维滑移线和三维滑移面算法、一种新的非结构网格并行计算方法. 对于流体力学多介质大变形问题数值模拟, 《多介质流体动力学计算方法》是国内第1次比
较系统的讨论, 提出了一系列高精度、健壮、有效、实用的数值方法, 并进行了大量的数值试验.

内页插图

目录

《信息与计算科学丛书》序前言第 1章导论 . 1
1.1流体力学多介质大变形问题的研究背景 1

1.2拉氏方法、欧拉方法和任意拉氏欧拉方法 2

1.3运动界面追踪方法 5

1.4 MMALE方法 . 5

1.5非交错网格 ALE方法8

1.6其他的处理多介质大变形问题的数值方法 . 10

1.7拉氏方法中的接触算法 12

1.8本书第 2章至第 13章的主要内容 13参考文献 . 17第 2章计算流体力学基础 . 26
2.1描述流体运动的基本方法 26

2.2积分形式的流体力学方程 27

2.3微分形式的流体力学方程 29

2.3.1欧拉坐标系中微分形式的流体力学方程 29

2.3.2拉格朗日坐标系中微分形式的流体力学方程组 32
2.4双曲型方程的间断解 34
2.5双曲型方程的唯一解及熵条件 . 36

2.6 RiemAnn问题及其求解方法 .38
2.6.1 RiemAnn问题及其求解的一般概念 .38
2.6.2 HLLC近似 RiemAnn解 . 39

2.6.3人工粘性 42参考文献 . 44第 3章二维流体力学的显式有限元相容拉氏方法 . 46
3.1引言 .46
3.2角质量和子网格质量 46
3.3动量方程半离散格式 50
3.4内能方程半离散格式 51
3.5全离散计算格式及边界条件的处理 52
3.6算例 .53参考文献 . 56第 4章 SALE方法的一种滑移线计算方法 58
4.1引言 .58
4.2 DYNA2D程序的 LAgrAnge滑移线计算方法 59

4.3一种等效的 LAgrAnge滑移线处理方法 63

4.4 SALE方法的滑移线算法 .71
4.4.1对主点作修正 71
4.4.2对从点作修正 77参考文献 . 77第 5章基于近似 RiemAnn解的有限体积 ALE方法 79
5.1引言 .79
5.2控制方程 80
5.3数值方法 81
5.3.1有限体积离散 81
5.3.2移动网格上的 HLLC格式 81
5.3.3 WENO重构 82
5.3.4时间离散 85

5.4算例 .85参考文献 . 90第 6章基于 level set的 Euler-LAgrAnge耦合方法 92
6.1引言 .92
6.2 Euler-LAgrAnge耦合方法的基本框架 . 93
6.3显式有限元相容拉氏方法简介 . 94

6.4 level set函数 . 95
6.5界面表示﹑ Ghost网格及时间和空间耦合格式 96

6.5.1界面表示 96

6.5.2 Ghost网格 . 97

6.5.3时间耦合方式 97
6.5.4空间耦合格式 98
6.6数值算例 .100
6.7结论 107参考文献 107
第 7章一种健壮有效的多介质欧拉方法 109

7.1前言 109

7.2控制方程 .110
7.3 level set方法 110
7.4数值格式 .111
7.4.1界面边界条件的定义 .111
7.4.2物理量控制方程的求解方法 113

7.5算例 113参考文献 116第 8章基于 MOF界面重构的二维多物质 ALE方法 117
8.1引言 117

8.2二维笛卡儿坐标系中的显式有限元相容拉氏方法 118

8.2.1角质量和子网格质量 .118
8.2.2动量方程半离散格式 .121
8.2.3内能方程半离散格式 .123
8.3二维子网格力学模型 . 124
8.3.1混合网格的描述 124
8.3.2几何工具箱 125
n+
8.3.3计算 t12时刻的量 . 126
8.3.4计算 tn+1时刻的量 . 127

8.4 MOF界面重构方法 129
8.5精确积分守恒重映 132

8.6多物质 ALE方法 . 132

8.6.1单元量重构 133
8.6.2单元量重映 133
8.6.3动量重映 134
8.7算例 134参考文献 139第 9章拉氏双曲守恒律的高精度间断 GAlerkin谱有限元方法 . 141
9.1前言 141

9.2控制方程 .143
9.3一些记号、 JAcobi多项式及二维正交谱基函数 144
9.3.1一些记号 144
9.3.2 JAcobi多项式 .145
9.3.3二维正交谱基函数 147
9.4几何变量的定义以及几何守恒律的离散 147
9.5物理守恒律的离散 150

9.6时间离散和时间步长计算 154
9.7节点速度的计算 155
9.7.1守恒关系 156
9.7.2熵不等式 158
9.7.3节点速度和边压力的计算 160

9.7.4向高阶格式 (sp > 0)的推广 162
9.8 HWENO重构 162
9.8.1 sp =1的重构 .163
9.8.2 sp =2的重构 .164
9.9算例 165

9.9.1精度测试 165
9.9.2其他算例 170参考文献 173第 10章三维流体力学的显式有限元相容拉氏方法 . 175
10.1引言 . 175

10.2角质量和子网格质量 175
10.3动量方程半离散格式 178
10.4内能方程半离散格式 180
10.5算例 . 181参考文献 185第 11章三维非结构网格基于 MOF界面重构的 MMALE方法 . 186
11.1前言 . 186

11.2 MMALE方法的计算流程 . 188
11.3拉氏步计算方法 . 189

11.3.1纯网格的算法 189
11.3.2混合网格的算法 . 190
11.3.3更新物质中心点的坐标 190

11.4混合网格界面重构 193
11.5网格重分 193
11.6三维非结构网格高精度精确积分守恒重映 194

11.6.1两个四面体的相交多面体的计算方法 194
11.6.2重映算法 . 198
11.7高精度精确积分守恒重映算法的并行计算 203

11.7.1数据分割 . 203
11.7.2并行环境生成 206
11.7.3高精度精确积分守恒重映算法的并行计算方法 .209
11.8数值算例 216
11.8.1重映方法精度测试 216
11.8.2三维周期漩涡问题 220
11.8.3两种材料的一维 Sod激波管问题 223

11.8.4三维 Sedov问题 224

11.8.5三维 Noh问题 228
11.8.6三维 SAltzmAn问题 233

11.8.7二维 “triple point”问题 234

11.8.8激波与氦气泡相互作用问题 235

11.8.9水中强激波与圆柱形空气泡相互作用问题 237
11.8.10二维 RAyleigh-TAylor不稳定性问题 239
11.8.11三维 RAyleigh-TAylor不稳定性问题 241
11.8.12圆柱形内爆中的 RAyleigh-TAylor不稳定性 246参考文献 249第 12章三维非结构网格基于浓度法的非交错网格 MMALE方法 254
12.1前言 . 254

12.2非交错网格 MMALE方法的计算流程 255

12.3拉氏步计算方法 . 255

12.3.1纯网格的算法 256
12.3.2多介质混合网格的算法 260

12.4网格重分和物理量重映 261

12.5算例 . 261

12.5.1三维周期漩涡问题 261
12.5.2三维 Sedov问题 263

12.5.3三维 RAyleigh-TAylor不稳定性问题 .266参考文献 268第 13章交错型拉氏方法的两种新的三维滑移面算法 270
13.1引言 . 270

13.2三维有限体积相容拉氏方法的简要回顾 274
13.3计算两个三角形的相交多边形的面积的一个简单方法 276

13.4离散精确匹配法 . 278

13.5离散拉格朗日乘子法 282
13.5.1 计算含有一个碰撞点和一个目标点的一维接触对的接触力的拉格朗日乘子法 282
13.5.2 离散拉格朗日乘子法的接触处理算法 284
13.5.3 与防御节点算法的区别 285

13.6三维滑移面算法的并行计算方法 286
13.7算例 . 287参考文献 318索引 322《信息与计算科学丛书》已出版书目 . 326

精彩书摘

第 1章导论
1.1流体力学多介质大变形问题的研究背景
在内爆动力学、惯性约束聚变 (inertiAl con.ned fusion,ICF)、界面不稳定性、高速冲击、先进常规兵器等国防和高新技术领域 ,大规模、高置信度数值模拟是必不可少的工具 ,正在发挥越来越重要的作用 .数值模拟作为除了理论研究和实验研究之外的第三种研究方法 ,它既是连接理论研究和实验研究的桥梁 ,也是一种新的实验手段 (即数值试验 ).它对研究和发现物理过程中的规律 ,认识其中的机理具有重要的作用 .
上述这些领域所遇到的物理问题常常非常复杂 .从数学上来看 ,主要表现在需要求解可压缩流体力学方程组 ,其计算区域往往是三维复杂区域 ,问题具有动边界、多介质、大变形、强间断、强非线性以及多物理过程强耦合等特点 ,并且要求准确、清晰地刻画物质界面 .其中多介质大变形问题是流体力学计算的难点 ,迄今没有得到很好的解决 .
例如 ,为了探索实现氘氚核聚变 ,提供解决人类清洁能源问题的最终有效途径 ,中国和世界科技强国都开展了 ICF研究 .在激光间接驱动 ICF过程中 ,首先将多束高能激光注入黑腔 ,腔壁在吸收了入射的激光后将其转换成软 X射线辐射 (如图 1.1左图所示 ),之后再利用辐射驱动猛烈压缩置于腔体中间的靶丸 (氘氚燃料小球,外覆碳氢等多层材料 ),使其达到极高能量密度之聚变点火燃烧条件 .在靶丸压缩过程中 ,由于存在复杂波系的相互作用 ,其烧蚀层、氘氚冰层、氘氚气的各个界面上会出现流体力学界面不稳定性 (如图 1.1右图所示 ),包括 RAyleigh-TAylor不稳定性和 Richtmyer-Meshkov不稳定性 .当界面上的压强梯度和密度梯度方向相反时会出现 RAyleigh-TAylor不稳定性 ,当激波经过两种流体的界面时会出现 Richtmyer-Meshkov不稳定性 .界面不稳定性导致物质界面变形、翻转、破碎以及后期的湍流混合现象 ,直接影响靶丸的最终压缩状态 ,对聚变成功与否至为关键 .因此在 ICF中界面不稳定性以及与之相关的多介质大变形是最具挑战性的问题之一 .


1.2拉氏方法、欧拉方法和任意拉氏欧拉方法
基于网格类的流体力学计算方法主要分为采用物质网格的拉氏方法 (LAgrAng-iAn method)、采用固定空间网格的欧拉方法 (EuleriAn method)和采用运动网格的任意拉氏欧拉方法 (ArbitrAry LAgrAngiAn EuleriAn method,简称 ALE方法).
拉氏方法采用物质网格 ,即网格的运动速度取为物质的运动速度 ,从而避免了输运计算及由此造成的误差 .同时用网格边界清晰地刻画、描述物质界面 ,界面处理的精度较高 .其缺点是当流场中发生大变形时网格极度扭曲 ,往往导致计算终止.图 1.2是拉氏方法示意图 ,其中 tn和 tn+1分别表示第 n时间层和第 n +1时间层 .拉氏方法最早是 von NeumAnn和 Richtmyer在文献 [1]中针对一维情况提出来的 . Wilkins在文献 [2]中将该方法推广到二维情况 .这种方法采用人工粘性(Arti.ciAl viscosity)处理激波间断 .采用这种方法的早期的计算格式求解时 ,由于该方法求解内能方程 ,总能量是不守恒的 ,并且会产生虚假的数值振荡 .尽管有这些缺点 ,在过去的几十年中 ,该方法在多介质流动的数值模拟中被广泛应用 .为了解决上述这些问题 ,最近十几年来 ,针对该方法的缺点作了很多改进 . CArAmAnA和 ShAshkov在文献 [3]中提出采用子网格压力消除沙漏运动和虚假漩涡 .在文献 [4,5]中, ShAshkov等基于支撑算子的思想建立了保持总能量守恒性的相容拉氏动力学计算方法 .其主要特点是保证控制方程的主要物理性质和方程中各微分算子之间严格的数学关联在离散后得到保持 ,将控制方程的基本物理性质融入到数值计算中去 ,从而改善了数值模拟结果的可靠性和预测能力 .此外 ,人工粘性的计算方法也取得了明显的进步 [6,7].

欧拉方法采用固定空间网格进行计算 ,能够处理大变形问题 .其缺点是由于输运计算带来的误差 ,欧拉方法很难给出精确的物质界面 .图 1.3是欧拉方法示意图 .近几十年来欧拉方法取得了长足的进步 ,发展了基于 RiemAnn间断解的 Godunov方法 [8]和 Roe方法 [9]、积分平均型间断解方法 MUSCL[10]和 PPM[11],高分辨率方法 TVD[12]、ENO[13]、WENO[14]、紧致格式 [15],以及间断有限元方法 [16-19]等.其中 ENO方法采用逐次扩展的节点模板来构造高阶的逼近多项式 ,利用比较差商绝对值大小进行取舍的方法判定、选择扩展的节点模板 ,尽量避免在所选择的模板中包含间断 ,以此提高插值方法精度并实现高分辨率和无振荡的效果 .在 ENO方法的实施过程中 ,有许多中间计算结果被可惜的丢弃不用 .文献 [14]提出了 WENO方法 ,对这个不足作了弥补 ,因此得到了更为广泛的应用 .间断有限元方法是提高数值逼近精度的一个有效方法 .间断有限元方法的出现 ,最早可以追溯到 1973年 Reed和 Hill关于中子输运方程问题的论文 [20]. 20世纪 90年代以来 ,以 Cockburn和 Chi-WAng Shu为代表提出的 Runge-KuttA间断 GAlerkin有限元方法[16-19]在解决含有间断现象的问题中发挥着越来越大的作用 ,卓有成效地应用到了水动力学、气动力学和波传播等问题 .间断有限元方法的基本思想是用检验函数乘以原方程并在网格单元上积分 ,然后通过分步积分获得原问题的弱形式 .选取有限元逼近函数空间 (允许其在单元边界不连续 ),选择合适的数值流通量构造单元边界连接条件 ,形成可解的封闭线性代数方程组 .间断有限元方法具有很多优点：易于处理复杂边界和边值问题 ;具有灵活处理间断的能力 ,克服了一般有限元不适合间断问题的缺点 ;可以通过提高单元插值函数的次数来提高精度 ,不必像有限体积方法那样首先扩大节点模板 ,再做复杂的多项式重构 ;容易实施自适应策略 ;可以适当选取基函数 ,使得质量矩阵是分块对角的 ,容易求逆 ;为了求解给定单元的自由度,只需要相邻单元的自由度 ,易于实现并行计算 .用间断有限元等高阶格式计算含有强间断的问题会出现非物理振荡甚至非线性不稳定 (Gibbs现象 ). Qiu JiAnxiAn和 Shu ChiwAng提出了一个小模板的 HWENO(Hermite weighted essentiAlly non-oscillAtory)重构方法 [21,22]来抑制非物理振荡 ,并用作间断 GAlerkin有限元方法的限制器 .此外 , H. Luo等也提出了一个类似的方法 [23].

为了克服拉氏方法和欧拉方法的缺点并结合它们的优点 ,人们提出并发展了 ALE(ArbitrAry LAgrAngiAn-EuleriAn)方法 [24-29],其主要特点是采用网格速度可人为控制的运动网格 ,物质内部网格可以通过优化保持较好的几何品质 ,能够适应大变形问题的数值模拟 ,同时物质界面仍用网格边界显式刻画 ,界面处理清晰、准确 .图 1.4是 ALE方法示意图 .由于上述优点 ,近几十年来 ALE方法一直是流体力学计算方法研究的热点 ,具有重要的学术价值和工程应用价值 ,在相关基础研究领域和应用基础研究领域 ,特别是内爆动力学和 ICF等领域得到广泛应用 .
根据物理量定义位置的不同 , ALE方法大致可以分为两类：一类为运动量 (如位移、速度和加速度 )定义在网格节点、状态量 (如密度、压力、内能等 )定义在网格中心的交错网格 ALE方法 (stAggered grid ALE method),一类为所有物理量 (包括运动量和状态量 )均定义在网格中心的非交错网格 ALE方法 (cell-centered ALE method).
另一方面 ,根据物理量更新方式的不同 , ALE方法大致可以分为两种：比较常见的一种 ALE方法通常由三个步骤组成 ,即拉氏步计算、网格重分和物理量重映 (如文献 [24-27]).当网格变形不大、形状较好时 ,仅进行拉氏步计算 ,把 tn时刻的拉氏网格上的物理量更新为 tn+1时刻的拉氏网格上的物理量 ,高精度地追踪物质界面 ;当 tn+1时刻的拉氏网格变形较大、网格品质较差时 ,首先对该网格进行优化 ,生成品质较好的新网格 ,然后进行重映计算 ,把 tn+1时刻的拉氏网格上的物理量映射到 tn+1时刻的新网格上 .另一种 ALE方法首先确定 tn时刻的网格速度 ,得到 tn+1时刻的新网格 ,然后直接求解含有网格速度的 ALE形式的流体力学方程组 (如 [28,29]),把 tn时刻的拉氏网格上的物理量直接更新为 tn+1时刻的新网格上的物理量 .文献 [30, 31]中指出 ,这两种 ALE方法的精度是相当的 .


1.3运动界面追踪方法
前面已经提到 ,由于输运计算带来的误差 ,欧拉方法很难给出精确的物质界面 .例如 , TVD、Roe等方法在光滑区域可以达到三阶、四阶精度 ,而在关键的间断区域却只能是一阶 (物质界面也是一种间断 ).为此 ,近几十年来人们发展了一类界面追踪方法 ,用于模拟、描述和追踪自由面和运动界面的轨迹、特征和发展 .目前比较流行的是 VOF(volume of .uid)方法 [32-33]、等值面 (level set)函数方法 [34,35]等.其中 VOF方法是在整个流场中定义一个函数 (称之为流体体积函数 ),在每个网格中 ,这个函数定义为一种流体 (称之为目标流体 )的体积与网格体积的比值 .在任意时刻 ,知道了这个流体体积函数在每个网格上的值 ,就可以通过某种途径显式地构造出运动界面 .然后在求解物理方程时可以在界面附近作特殊的精细处理 ,以提高分辨率和精度 .另一种较理想的做法是利用所谓的等值面函数 (level set function). (x,t)代替 VOF方法中的流体体积函数 (x是空间变量 , t表示时间 ).让 .以适当的速度移动 ,使其零等值面就是物质界面 .在任意时刻 ,只要知道 .,然后求出其零等值面,就知道了此时的运动界面 .等值面函数法不需要显式地追踪运动界面 ,从而可以较容易地处理复杂的物质界面及其拓扑结构发生变化的情形 ;而且界面的一些特征 (如法向、曲率等 )直接隐含在 level set函数中 ,便于精细地描述界面 ;此外 ,它还易于向高维推广 .由于这些优点 ,等值面函数法已经被用来处理几何、流体力学、工艺过程等许多方面的问题 .

1.4 MMALE方法
传统的 ALE方法在一个网格中只允许含有一种介质 ,在计算多介质问题时 ,一般将物质界面取为计算子区域的边界 ,与边界上的网格边重合 ,界面在法向保持拉氏运动 .这种情况下界面至少应该保持连续 ,最好不要发生较大的变形 ,只有这样才能保证在进行网格重分时能够在子区域中生成质量较好的新网格 .但是在内爆压缩等一些实际问题的数值模拟中 ,不仅网格随流体运动发生大变形 ,而且界面也会发生较大的变形 ,甚至出现破裂、合并等拓扑结构改变的情形 ,这时传统的拉氏方法和 ALE方法将无法继续应用 .美国 SAndiA国家实验室的 Peery等在 2000年提出的 MMALE方法 (multi-mAteriAl ArbitrAry LAgrAngiAn EuleriAn method)[36]是解决这个问题的一个比较有效的方法 ,其基本思想是把传统的 ALE方法和界面追踪方法相结合 ,处理包含网格扭曲以及界面破裂、合并等复杂现象的多介质大变形问题.其实现途径是在传统的 ALE方法的基础上 ,通过引进混合网格 (即允许一个网格内含有多种介质 ,物质界面可以穿过网格 ),用界面追踪方法在混合网格中刻画、重构运动界面 ,来进一步提高 ALE方法处理多介质大变形问题的能力 .他们将该方法应用于 ALEGRA程序中 ,计算 ICF研究中遇到的含有强剪切变形的多介质大变形问题 .图 1.5是 MMALE方法示意图 .

图 1.5 MMALE方法示意图
粗实线表示用网格边刻画的物质界面 ,虚线表示在混合网格中用界面追踪方法刻画、重构物质界面

文献 [36]中的 MMALE方法包含以下三个关键技术 .一是对混合网格提出一个热力学封闭模型 ,这个模型用来计算经过一个拉氏步以后混合网格中每种物质的体积分数和物理量的变化 .文献 [36]采用平均应变率混合模型 [37].在这个模型中 ,混合网格中各种介质的应变率就取为网格的平均应变率 ,这实际上就是假设经过一个拉氏步以后各种介质的体积份额保持不变 .当混合网格中界面两侧密度比或压力比很大时 ,应用这个模型会得到非物理的结果 [38].二是用重分重映处理网格的大变形.文献 [36]采用对流重映方法 [39-41],这种方法通过计算网格边上的通量来得到新网格上的积分量 ,因而要求新网格充分接近老网格 .但是这种方法不需要计算老网格和新网格的交点 ,计算效率较高 .三是用 VOF方法处理界面的大变形 .
近年来 ,由于内爆动力学和 ICF等领域中多介质大变形问题研究的需求牵引 , MMALE方法受到学术界的广泛关注 ,成为计算流体力学的一个研究热点 .
在混合网格的封闭模型方面 ,为了克服平均应变率混合模型的缺点 ,提出了几个新的封闭模型 ,例如： 1压力平衡模型 [42].这种模型假设混合网格中各种介质的压力在瞬间达到平衡 .压力平衡是非线性问题 ,需要设计复杂的迭代格式进行计算 .
2压力松弛模型 [43-46].这种模型引进一种类似于粘性的松弛机制使得混合网格中的压力趋于平衡 . 3一维子网格力学模型 [45,47-49].这个模型用声波 RiemAnn解估计界面速度 ,并由此计算每种物质的体积分数和物理量的改变量 . 4质量分数模型[50,51].这个模型也被称为浓度法 .就计算效率和有效性而言 ,压力松弛模型由于

前言/序言

<div>　　在内爆动力学、惯性约束聚变、界面不稳定性、高速冲击、先进常规兵器等国防和高新技术领域，高性能、大规模数值模拟是必不可少的工具，正在发挥越来越重要的作用。在上述这些领域遇到的物理问题常常非常复杂，从数学上来看，主要表现在需要求解可压缩流体力学方程组，其计算区域往往是三维复杂区域，问题具有动边界、多介质、大变形、强问断、强非线性，以及多物理过程强耦合等特点，并且要求准确、清晰地刻画物质界面。其中多介质大变形问题是流体力学计算的难点，迄今没有得到很好的解决。为了解决这些困难，国内外的计算流体力学工作者进行了不懈的努力，各种新思想、新方法应运而生。</div><div>　　对于流体力学多介质大变形问题数值模拟，本书在国内第一次比较系统地进行了讨论，提出了一系列高精度、健壮、有效、实用的数值方法，并进行了大量的数值试验。本书是作者及其课题组成员多年来在流体力学多介质大变形问题数值模拟方法研究方面所取得的成果总结。课题组曹雄研究员、刘军助理研究员等为本书的内容作出了重要贡献。</div><div>　　作者撰写本书的目的有两个：一是和同行进行交流，希望得到同行专家的指正；二是期望为各行各业同行的教学、科研工作起到一点助益，为国家的人才培养、科技进步有所奉献。由于本书涉及的内容比较广泛，且各种计算方法都在不断发展中，加上作者的学术水平有限，书中难免会有缺点和错误之处，恳切希望同行专家及读者批评指正。</div><div>　　作者感谢国家出版基金和国家自然科学基金面上项目(11072040)以及中国工程物理研究院科学技术基金重点项目(2012A0202010)的资助。</div><div>　　本书适合从事内爆动力学、惯性约束聚变、界面不稳定性、高速冲击、先进常规兵器、天体物理等领域数值模拟研究与应用的科技工作者阅读使用，也可作为高等院校计算数学和计算流体力学专业研究生和高年级本科生的教材与参考书。</div><div>　　贾祖朋 张树道 蔚喜军</div><div>　　2014年3月于北京</div>

《信息与计算科学丛书（63）：多介质流体动力学计算方法》 —— 深入探索复杂流体系统的数值模拟前沿 丛书定位与背景 《信息与计算科学丛书》一直致力于汇集信息科学、数学、物理学交叉领域的前沿研究成果与经典理论。本分册，作为丛书中的第63部力作，聚焦于当前计算科学领域极具挑战性与应用价值的方向之一：多介质流体动力学（Multiphase Fluid Dynamics）的数值模拟方法。 在现代工程、环境科学和基础物理研究中，流体系统很少是单一成分的。从自然界的气-液、液-固两相现象（如波浪、雾化、沉积）到工业生产中的复杂混合过程（如石油天然气开采、化工反应器、先进制造中的金属凝固），都涉及多相流体的相互作用。精确地捕捉和预测这些系统的行为，对提高效率、保障安全以及推动新材料研发至关重要。然而，由于界面追踪的复杂性、不同相间物质交换的非线性耦合以及极端条件下的物态变化，多介质流体动力学问题的解析解几乎不可得，这使得高效、稳定且高精度的计算方法成为核心瓶颈。 本书正是为填补这一领域内理论与实用方法之间鸿沟而设计，旨在系统梳理和深入剖析当前主流及新兴的多介质流体动力学计算框架。 --- 核心内容与技术深度解析 本书的结构经过精心设计，从基础理论出发，逐步深入到前沿的先进算法与大规模并行实现。全书的核心内容涵盖以下几个关键板块： 第一部分：多相流基础理论与数学模型重述 本部分首先回顾了连续介质力学的基础，重点阐述了如何将单相流体方程（Navier-Stokes方程）推广至多相环境。内容包括： 1. 相变与界面动力学基础： 详细分析了相界面（如气液界面、固液界面）的物理描述，包括表面张力、Marangoni效应的引入，以及热力学驱动的相变过程（蒸发、凝结、熔化）。 2. 宏观描述模型比较： 深入对比了描述多相流的两种主要范式： 欧拉-欧拉（Eulerian-Eulerian）框架： 侧重于平均化处理，介绍相际动量、质量和能量交换的建模（如MFIX模型、子网格模型）。 欧拉-拉格朗日（Eulerian-Lagrangian）框架： 适用于稀相分散流，重点阐述粒子追踪、曳力模型和湍流对离散相的影响。 3. 界面捕捉方法（Interface Capturing）： 详细讲解了用于描述清晰界面的先进技术，如水平集（Level Set Method, LSM），包括其演化方程、曲率计算以及处理界面拓扑变化的挑战。 第二部分：先进的界面解析与追踪技术 在涉及显著界面结构（如自由表面流、液滴破碎）的问题中，直接追踪界面至关重要。本部分着重介绍了高精度界面解析方法： 1. 体积平均法与守恒性： 深入讨论了体积平均方法（Volume of Fluid, VOF），强调其在质量守恒方面的优势。详细介绍了VOF的重构技术（如PLIC方案），以确保界面在网格尺度上的锐利性。 2. 守恒型界面方法（Conservative Interface Methods）： 对体积平均拉格朗日（VAL）等混合方法的原理进行剖析，特别是在处理界面移动和拓扑变化时的数值稳定性。 3. 几何与拓扑处理： 探讨了如何结合LSM和VOF的优点，实现既能精确追踪界面形状，又能保证质量守恒的混合方法，例如相场模型（Phase-Field Model）在介观尺度上的应用及其向宏观流体动力学方程的过渡。 第三部分：求解器的设计与数值实现 计算方法的效率和精度高度依赖于底层离散化方案和求解器的构建。本部分是本书的工程核心： 1. 空间离散化技术： 系统对比了有限差分（FDM）、有限体积（FVM）和有限元（FEM）在处理多相流问题时的适用性与局限性。特别关注了高分辨率格式（如WENO、TVD）在抑制相界面处的数值振荡方面的应用。 2. 时间推进策略： 讨论了隐式、显式和半隐式时间步进方案。针对强耦合和刚性方程组（如涉及高速流动或强热传导时），详细分析了压力-速度耦合算法（如SIMPLE族算法的扩展）在多相流求解中的改进，尤其是如何处理不同相间的压力跳跃。 3. 网格技术与自适应： 介绍了处理剧烈界面运动所需的动网格技术（Adaptive Mesh Refinement, AMR），如何根据界面曲率或速度梯度自动细化计算区域，从而在保证精度的同时节省计算资源。 第四部分：特定物理过程的计算挑战与前沿应用 最后，本书将理论与方法应用于几个典型的复杂物理场景，展示计算方法的威力： 1. 气液/液液界面湍流耦合： 讨论了如何将雷诺平均纳维-斯托克斯（RANS）或大涡模拟（LES）与界面追踪方法结合，模拟界面处的湍流生成、传输和耗散。 2. 固液相变与冶金模拟： 重点解析了用于模拟凝固过程的计算模型，包括枝晶生长、孔隙率演化，以及界面移动与潜热释放的耦合求解。 3. 多相流的并行计算与高性能实现： 针对大规模仿真需求，探讨了域分解、负载均衡以及如何高效地在GPU或CPU集群上实现VOF/LSM等复杂算法的并行化（如使用MPI和OpenMP）。 --- 本书的目标读者 本书适合于流体力学、计算数学、化学工程、航空航天、材料科学等领域的高年级本科生、研究生、科研人员以及相关领域的工程技术人员。它不仅是理解多介质流体动力学数值模拟最新进展的理论参考，更是指导实际工程计算工具开发与应用的高阶教材。阅读本书需要具备流体力学基础知识和熟悉数值方法的基本概念。 通过对这些复杂计算方法的系统梳理，本书旨在帮助读者建立起一套全面的、可操作的工具箱，以应对从微观尺度到工程尺度的多介质流体动力学难题。

评分☆☆☆☆☆

一直以来，我对流体动力学的研究都充满了热情，尤其是在湍流、多相流等领域，总觉得还有许多未解之谜等待探索。我最近在进行一个关于海底油气开采模拟的项目，其中涉及到复杂的多相流体（油、气、水）以及与固体（井筒、储层）的相互作用，这让我深切体会到数值模拟在解决这类复杂工程问题中的关键作用。我非常希望能找到一本能够系统性地介绍如何处理多介质流体动力学问题的书籍，它最好能够涵盖从理论模型建立到具体计算方法实现的整个流程。我期望这本书能够提供一些关于如何选择合适的数值离散格式、如何处理流体界面、如何耦合不同物理过程（如传热、传质）以及如何进行大规模并行计算的详细指导。如果书中能够包含一些实际案例的分析，或者能够提供一些开源代码的参考，那将是更棒的事情，这对于我实际应用中所遇到的挑战将会有极大的帮助。

评分☆☆☆☆☆

作为一个对计算科学领域充满好奇的学习者，我一直以来都很喜欢《信息与计算科学丛书》系列。我印象深刻的是其中一本关于现代数值方法在偏微分方程求解中的应用，它系统地介绍了特征线法、有限元法以及谱方法等，并分析了它们的收敛性、稳定性和适用性。同时，我也对书中关于数据结构和算法效率的讨论颇有心得。我非常希望这本关于多介质流体动力学计算方法的新书，能够在延续该系列严谨的学术风格的同时，为我开启一个全新的研究方向。我期待它能解释如何将那些基础的数值技术巧妙地应用于处理具有复杂物理行为的多介质流体系统，例如，如何将这些数值方法与多相流模型（如两相流、三相流）以及非牛顿流体等概念结合起来，并给出具体可行的计算流程和算法实现。

评分☆☆☆☆☆

我对流体力学，尤其是那些涉及到真实世界复杂场景的模拟，一直抱有浓厚的兴趣。我曾经尝试阅读过一些介绍计算流体力学基础的书籍，它们详细地讲解了Navier-Stokes方程的推导、数值离散化方法（如有限差分、有限体积）以及一些基本的流动模拟案例。然而，当我开始研究一些涉及多相流、界面现象甚至相变的课题时，我发现现有的知识体系显得有些单薄，迫切需要更专业、更深入的指导。我特别希望能在这本新书中找到关于如何精确捕捉和处理流体界面（例如，自由表面、液滴、气泡），如何模拟不同介质之间的相互作用（如表面张力、相界面的动量和能量传递），以及如何建立和求解描述多介质耦合动力学行为的数学模型。如果书中能提供一些关于如何选择和设计适应复杂物理过程的数值算法的建议，那对我将是巨大的帮助。

评分☆☆☆☆☆

这套《信息与计算科学丛书》一直是我非常钟爱的系列，尤其是当我知道第六十三卷将聚焦“多介质流体动力学计算方法”时，我内心充满了期待。我之前已经拜读过该丛书中关于数值分析、微分方程求解以及高性能计算等几本，可以说，它们都为我打下了坚实的理论基础，并提供了许多实用的算法思路。我记得有一本详细讲解了有限元方法在边值问题中的应用，让我对如何离散化连续域有了深刻的理解；另一本则深入探讨了并行计算的策略，为我解决大规模科学计算问题提供了重要指导。因此，我对于这本新书所能带来的知识增量，以及它将如何将这些前沿的计算科学技术与我一直感兴趣的多介质流体动力学领域相结合，抱有极大的好奇心。我猜想，它应该会进一步拓展我在数值模拟方面的视野，或许会介绍一些更复杂、更贴近实际工程应用的算法，比如如何处理界面追踪、相变以及复杂边界条件下的流体行为。这本书的到来，无疑是我学习道路上的一大福音，我迫不及待地想深入其中，汲取新的知识养分。

评分☆☆☆☆☆

我近期一直在关注计算流体动力学（CFD）领域的最新进展，特别是那些能够处理复杂边界和多相流的先进技术。我之前阅读过一本关于计算流体力学基础的书籍，它详细介绍了有限体积法、有限差分法等传统数值方法，也初步涉及了拉格朗日-欧拉方法（ALE）的概念。但是，在面对多介质流体时，我感觉现有的知识体系还有些不足，尤其是在处理不同相之间的界面捕捉和相变模拟方面。我希望这本新书能够在这方面提供更深入的讲解，例如，它是否会介绍一些更先进的界面追踪技术，如Level Set方法、VOF（Volume of Fluid）方法，或者更通用的相场方法？同时，我也很关心在多介质耦合模型中，如何有效地处理不同介质之间的物理属性差异以及相界面上的守恒律问题。如果书中能对此有详尽的阐述，并结合一些实际的工程应用场景进行分析，那我将受益匪浅。