索罗·苏奇(Sauro Succi)是一位物理学家,1987年博士毕业于洛桑联邦理工学院等离子体物理研究中心,现任意大利理工学院首席研究员。他在物理学领域取得了诸多成就,特别是在流体物理和计算物理方面。他的研究广泛涉及非平衡统计力学,从核聚变等离子体到工业设计的湍流,以及多孔介质、微纳米流体等。他还发表了多篇高引用率论文,并获得了包括Alexander von Humboldt奖、Gordon Bell超级计算竞赛入围奖在内的多项国际奖项。他不仅在学术界享有盛誉,还积极参与科学传播,为高中生讲授复杂系统科学。他现任多个国际科研机构的科学顾问,并在欧洲科学界发挥着重要作用。
◎序言
本书所介绍的格子玻尔兹曼(Boltzmann)方程(the Lattice Boltzmann Equation,简称 LBE)是玻尔兹曼方程极为独特的一种形式,专门用于解决流体动力学及其他学科中的问题。
LBE 最初由格子气元胞自动机(Lattice Gas Cellular Automata,简称 LGCA)衍生而来,迅速发展成为统计力学、特别是(离散)动理学理论框架下的独立研究对象。
大约十年前,主要在格子气元胞自动机的影响下,人们意识到,利用相空间可以对动量空间进行非常大胆的、在物理上又很合理的离散化处理,并且只要适当选择为数不多的几个离散动量,就能在连续介质极限下得到流体动力学方程。
这个发现不仅天然散发着智慧之光,还一举扫清了“真正的” 玻尔兹曼方程所涉及的计算障碍,为基于离散动理学理论的一系列相当高效的流体动力学(及其他学科中的)数值模拟方案铺平了道路。
无论是从分析角度还是从计算角度来说,玻尔兹曼方程的数学结构(“移动–碰撞”)是 LBE 成功的关键。
计算无小事。
LBE 计算中的大部分运算来自碰撞算子。因为粒子之间的碰撞在本质上是点状的,所以碰撞算子在位形空间中是完全局部的,这意味着可以同时处理物理空间的不同区域,而这正是并行计算的理想场景,也是 LBE 的关键优势。
这种内在的计算特性使 LBE 在某种程度上有别于更早的离散动理学理论模型,即 Broadwell 离散速度模型。
LBE 的主要关注点在于深刻理解复杂流动(通常通过数值模拟),而不是简化模型方程的精确解析解。因此,它的目标是流体动力学(以及更广泛学科)中最紧迫的一些问题,例如湍流、无序介质中的多相流以及含有悬浮物的流动。
这与近年来动理学理论的复兴完美契合,见 M. Ernst、E. van Beijeren 和 F. Rochelle 的优美表述(Dynamics: models and kinetic methods for non-equilibrium many-body systems, STATPHYS 20 Satellite Meeting, Leiden, July, 1998):“现代动理学理论提供了一个统一的理论框架,使我们能够一并研究各种看似毫无关联,并且表现出复杂动力学行为的物理系统。例如,这些方法被应用于胶体悬浮液动力学、颗粒流、介观系统中的电子输运等各种领域,还可用来计算经典多体系统的李雅普诺夫(Lyapunov)指数以及其他一些混沌特征,可谓包罗万象,异彩纷呈。”
没有词汇能更好地概括本书的“其他学科”部分。考虑到所有这些情况后自然会问,与现有的数值技术相比,LBE 应该带来哪些真正的优势呢?这个问题本身就带来了一小章——谁需要 LBE?我在这一章中(非常主观地)划分了LBE 的四类应用:DU、CU、SU、MU,分别代表“勿用(Don’t Use)”、“可用(Can Use)”、“该用(Should Use)”和“必用(Must Use)”。DU 类应用包含那些(在本书作者看来)更适合用其他现有技术处理的问题。(每一种值得尊敬的数值格式都必然有自己的 DU 类应用!)相反,MU 类应用是突破的摇篮,其中所包含的问题无法用任何其他方法解决。就目前而言,我相信 SU 类应用已经初具规模,而备受期待的 MU 类应用虽然尚显“稀薄”,但可能也并非一无所有。这个主题已经发展得足够成熟了,值得汇总成书,但是其研究前沿仍然不断涌现出新的进展。因此,希望读者在阅读本书后能充满自信地迎接挑战,在 MU类应用中留下自己的身影。
本书更多依赖于物理直觉和启发性论证,而不是数学上的详细分析,其目的是为更广泛的读者群体提供一份轻松的“入门”,而不是为希望在该领域深入钻研的专业人士提供全面的论述[1]。后一种类型的读者可以持续关注当前的文献,而本书没有取代这些文献的打算。
本书可供物理学、数学、工程技术和计算机科学等学科的研究生和高年级本科生参考。更广泛地说,我们希望格子气领域以外的研究者也能在本书中找到足够多的材料,从而自行判断是否值得在他们自己的研究项目中尝试这种方法。为了让他们对这种方法的态度偏向积极的一面,本书还提供一个正在开发的入门性质的初级计算机程序。
对统计力学和动理学理论的基础知识有所了解的读者不必阅读第 1 章。类似地,熟悉格子气元胞自动机的读者不必阅读第 2 章。撰写这两章的唯一目的是为了保证本书内容完备。
致谢
本书汇集了我十年以来关于格子玻尔兹曼方程这一主题的愉快研究。这十年让我有幸与许多杰出的同行建立了个人联系,他们给我诸多教益,多到远远超出我能够用语言表达的范围。对于以下致谢名单中的无意疏漏,我深表歉意:F. Abraham、B. Alder、G. Amati、M. Anile、G. Bella、R. Benzi、M. Bernaschi、J. Bernsdorf、M. Bertsch、B. Boghosian、J. P. Boon、A. Brandt、G. Brenner、H. Cabannes、N. Cabibbo、P. Carnevali、F. Castiglione、C. Cercignani、H. Chen、S. Chen、Y. Chen、M. L. Chiofalo、B. Chopard、G. Ciccotti、S. Ciliberto、E. Cohen、P. Coveney、J. Dongarra、G. Doolen、M. Droz、D. d’Humi`eres、T. Dupont、M. Ernst、B. Favini、O. Filippova、P. Franchi、D. Frenkel、U. Frisch、R. Gatignol、I. Ginzbourg、D. Gosman、P. Grosfils、J. Jimenez、F. Hayot、X. He、M. Henon、F. Higuera、A. Hoekstra、W. Hoover、L. Kadano、R. Kapral、I. Karlin、E. Kaxiras、W. Kohn、M. Krafczyk、A. Ladd、P. Lallemand、J. Lebowitz、L. S. Luo、M. Mareschal、N. Margolus、E. Marinari、F. Massaioli、W. Miller、P. Moin、K. Molvig、B. Nadiga、R. Natalini、H. C. Oettinger、I. Ohashi、S. Orszag、G. Parisi、L. Pietronero、R. Piva、Y. Pomeau、E. Presutti、M. Pulvirenti、Y. H. Qian、A. Quarteroni、C. Rebbi、S. Remondi、D. Rothman、D. Ruelle、P. Santangelo、G. Smith、J. Somers、D. Stauffer、C. Teixeira、T. Toffoli、F. Toschi、M. Tosi、R. Tripiccione、C. Tsallis、M. Vergassola、A. Vulpiani、V. Yakhot、J. Yeomans、J. Yepez、S. Zaleski、P. Zanella和 G. Zanetti。
特别感谢 Piero Sguazzero 博士,他建议我使用 HPP 自动机作为学习 IBM- 370 汇编语言的试验平台,从而开启了我的 LBE 探索之旅,这段旅程比预期更加趣味盎然。特别感谢 Roberto Benzi 教授,他分享了 LBE 早期研究中充满乐趣和兴奋的“欢乐时光”。
Hudong Chen(陈沪东)博士与我进行了无数次令人振奋的讨论,他与 EXA公司、波士顿大学和哈佛大学的其他多位朋友让波士顿成为了我的第二故乡,我向他们致以深深的谢意。我还要特别感谢 E. Kaxiras 教授,他为我安排的哈佛大学访问学者奖学金令我受益匪浅。
Steven Orszag 教授和 R. Piva 教授多年以来一直提供友好的支持,Carlo Cercignani 教授在本书得以面世的过程中发挥了至关重要的作用,我向他们致以最真挚的谢意。我也非常感谢 B. Boghosian 教授、J. P. Boon 教授、C. Cercignani 教授、P. Coveney 教授、U. Frisch 教授和 L. S. Luo 教授对本书初稿的批判性审阅。
在这份致谢名单中还必须补充我在动理学理论方面最早的几位老师,即 V. Boffi 教授、V. Molinari 教授和 G. Spiga 教授,多年以前从他们那里学到的基本知识至今仍然很有用。我同样要感谢 K. Appert 博士和 J. Vaclavik 博士,他们教会了我怎样把优美的理论转化为同样优美的数值方法,以便解决相关问题。
我要感谢 M. Adamo 博士负责本项目的图形处理,感谢我的同事 G. Amati、Y. Chen、O. Filippova、O. Inamuro、I. Karlin、A. Ladd、L. S. Luo、A. Masselot 和 J. Yeomans 慷慨提供了其图形的原始 postscript 文件,也非常感谢牛津大学出版社的 Soenke Adlung 博士和 Anja Tschoertner、R. Lawrence 和 J. Harris 在整个项目中的鼎力协助。
我还要向我的妻子 Claudia 和女儿 Caterina 致以深深的谢意,这或许是一种标准态度,但我是真心实意的。
即使是薄薄的一本书也难免出错,而这本书不算……
预先感谢指出书中各种漏洞、错误和其他任何不妥之处的热心读者,请使用以下电子邮件地址:succi@iac.rm.cnr.it。
S.
罗马
2000 年 9 月
[1] 在最后校对本书清样的过程中,我注意到由 D. A. Wolf-Gladrow 撰写的一本非常好的书:Lattice Gas Cellular Automata and Lattice Boltzmann Models. Lecture Notes in Mathematics. Berlin: Springer, 2000. 这本书对更侧重数学的读者来说是对本书的有益补充。