始于实验，陷于仿真，忠于流体

实验和仿真就好像推动产品研发的两副车轮，他们时而相互扶持，时而相互鄙视，像极了一对相爱相杀的伙伴。不过有一位流体力学的大神却用自己的实际行动完美诠释了实验和仿真是可以手牵手，一起做朋友，那便是提出经典的k-epsilon湍流模型的Brian Edward Launder教授。

01从实验流体力学开始

1939年，生于伦敦的Launder很早便展示出自己的学霸属性。1961年，年仅22岁的Launder在伦敦帝国理工学院获得机械工程学士学位，并以全班第一名的身份获得Bramwell奖章。随后，他进入麻省理工学院攻读硕士和博士学位，并以实验流体力学作为自己研究的起点。

Launder在MIT的燃气轮机实验室主要从事边界层实验工作，其主要的贡献是通过实验研究湍流边界层的加速效应。Launder在1964年发表了经典的论文《Laminarization of the turbulent boundary layer in a severe acceleration》阐述了他的研究内容。简而言之，如果加速足够强，湍流边界层可以变成层流，从而大大减少传热速率，这种现象被称为层流化，对于涡轮的流动至关重要。

02又回到最初的起点

在MIT完成博士学位后，Launder于1964年回到帝国理工学院担任讲师，后来成为流体机械专业的教授，直到1976年去了加州大学戴维斯分校。Launder在帝国理工学院供职12年，从25岁到37岁，Launder把最才华横溢的12年给了自己的母校。

在帝国理工学院，Launder加入了当时已经功成名就的Brian Spalding教授的团队，而就在同年，CFD行业的另外一位大佬Patankar进入帝国理工学院，师从Spalding教授攻读博士学位。

也许当时的人们很难想象，Spalding的团队将在接下来的二十年收获一系列引领CFD行业的成就。作为现代CFD行业的创始人之一，Spalding一直致力于CFD的工程化，其中很重要的一个研究方向就是为工业领域的流体仿真应用搭建合适的湍流模型，而Spalding将这一重要的研究方向交给Lauder来领导，于是便有了后来大名鼎鼎的k-epsilon模型。

03k-epsilon模型的诞生

还记得我们在前面的文章中提到过，早在1924年，普朗特基于“混合长度理论”提出了CFD领域第一个实用的涡粘模型，封闭了雷诺平均的N-S方程的求解。作为湍流模型的起点，普朗特将混合长度定义为与到壁面的距离成正比，在今天看来当然不够严谨，他不仅忽略了边界层内的变化，也忽略了湍流本身的演化。

或许是得益于早年在MIT做实验的经历，Launder亲眼目睹过湍流的生成和耗散，尤其是流体在边界层附近的演化。另外，Launder比普朗特幸运的是他站在更多的巨人肩膀上，尤其是1940年代，以柯尔莫果洛夫为首的统计学派为流体力学带来了全新的视野，让人们在湍流的无序中看到有序，为无规则的湍流制定了规则。于是，Launder很敏锐的发现，如果用湍动能和耗散率来表示混合长度能够更准确的从物理层面上反映湍流的生成和耗散。

04k-epsilon模型的推导

正如前面所说的，k-epsilon模型属于涡粘性模型，不过它变换了另外一种更加湍流的思想来求解混合长度，进而求解涡粘性系数：混合长度定义为流体微团耗散前所经过距离，因此，也可以定义为耗散时间和流体微团速度的乘积，当然流体微团的速度和脉动速度是相当的，通过经验系数来表示。于是我们可以用湍动能和湍流耗散率来表示湍流粘性系数。

而面对多出来的两个未知量k和epsilon，则需要补充两个微分方程来让方程的求解封闭。这便是我们熟悉的k方程和epsilon方程。不过需要特别指出的是补充的epsilon方程中包含了经验系数，这些系数一般由实验得出，因此限制了k-epsilon模型的适用范围。实际工程应用中，也有很多工程师会根据特定的应用场景修改相应的经验参数来匹配实验的结果。

k-epsilon模型本质上是在湍流各向同性假设基础上，对混合长度理论的进一步发展，因此相对于普朗特的模型，k-epsilon模型适用于更复杂的流动。但是熟悉k-epsilon模型的朋友都知道，它不适用于逆压梯度较高的流动，不适用于喷流和二次流。于是，Launder在后续提出了一种更彻底的湍流模型。

05最彻底的湍流模型

文章提出了一种湍流模型，不再使用涡粘性假设，其中雷诺应力由这些变量的输运方程和湍流能量耗散率epsilon的解来确定。模型特别注意压力和应变关系的近似；所采用的形式对壁面附近和自由剪切流中的应力进行了合理的分配。对于应变均匀剪切流和二维非均匀剪切流，包括射流、尾迹和尾流，新的湍流模型相对于k-epsilon模型在求解精度方面有了明显的改善。此外，研究还表明，对于曲面上的流动，新的模型确实预测了二次应变项的强烈影响。这种方法称为雷诺应力输运模型或者更一般地称为二阶矩封闭。

二阶矩封闭的一个主要优点是，表示雷诺应力产生的项在它们各自的输运方程中精确地表示出来。当然，这些方程中的其他项必须是近似的，与两方程模型相比，现在有七个输运方程需要求解——需要更多的计算工作。当然，这种方法能够捕捉到比k-epsilon模型更复杂的流动和现象。雷诺应力输运模型由Launder和他的同事G. J. Reece和W. Rodi共同提出，因此也称为LRR模型，很多人也将它称为最彻底的湍流模型。

06辗转于高校与科研机构之间

虽然Launder在近二十年的时间里并没有停下科研的脚步，不过相对于在帝国理工大学期间取得影响CFD发展的成就，Launder后来的工作更多的是对之前工作的一些延续和总结，尤其是针对雷诺应力输运模型的延续。为了表彰Launder对湍流建模的贡献，英国皇家学会和皇家工程院都将他选为Fellow。当然，他还获得了许多国际荣誉，包括很多大学的荣誉学位。

结语

作为CFD行业的先驱， Brian Launder教授一直致力于湍流模型的开发，而由他提出来的k-epsilon模型及其变种更是CFD走向工业化的基石。时至今日，各种商用CFD软件中，下拉湍流模型选项，k-epsilon都是当之无愧的王者。

然而成名于CFD的Launder却一直记得当年在MIT做实验的时光，也一直坚定的认为实验和仿真是相辅相成的，他无疑是“知行合一”的杰出践行者。最后就用Launder的一段话作为总结吧：“我总是看到实验和建模并驾齐驱。虽然我关于建模的论文被引用最多，但实验无疑给了我洞察力，也为测试建模思想提供了数据基础。”

-END-