首页 发现 课程 培训 文章 案例 问答 需求 服务 行家 赛事
消息
首页/文章/ 详情

人物 | D.B.Spalding：在湍流、传热、传质和燃烧等方面的贡献

积鼎科技
积鼎科技
1年前浏览4413


本文概述了D.B.Spalding教授迄今在科学和工程方面的一些贡献。从Spalding早期关于燃烧的工作和他在传质理论方面的独特工作出发，简要介绍了Spalding未发表的“统一理论”。在此之后，帝国理工学院小组在算法方面的发展导致了现代计算流体动力学的诞生，包括著名的SIMPLE算法。还介绍了燃烧、多相流和湍流模型的发展。最后，提到了计算流体动力学的一些学术和工业应用以及随后几年考虑的传热应用。

SIMPOP
翻译自《Internatinal Journal of Heat and Mass Transfer》作者见文末                
               
介绍                  
2008年5月，为纪念BrianSpalding教授85岁生日，由国际热质传递中心主办的国际研讨会(CHT ' 08: Computational Heat Transfer Advances)在摩洛哥马拉喀什举行。会议包括Spalding教授本人的主题演讲和一些介绍他的研究贡献的演讲。会议上决定以Spalding教授的名义出版该杂志的一期特刊。虽然这一期将邀请各种主题的论文，但似乎也应该包括对Spalding研究贡献的广泛回顾。                
               
带着这样的想法，这篇论文的作者(他们中的大多数都与Spalding教授有密切的合作)在会议上开会，规划了这篇论文的性质和结构。那次会议上的讨论和此后几个月的工作最终产生了这一颂词。
在这里，作者不仅给出了Spalding出版物的事实总结。他们试图带领读者踏上一段历史之旅，展示他的主要思想的发展及其随后的影响。作者展示了概括或统一的主题如何贯穿他所有的作品。正是这种宏大的愿景使Spalding做出了如此多的突破性贡献。在本文中，读者将深入了解这位伟人的思维方式以及他对科学和工程的非凡影响。作者们都非常享受撰写这篇论文的过程，并在此过程中重温了与Spalding教授一起工作的经历。                                  
命名法                  

在这种性质的论文中，不需要使用通常的命名法部分。Spalding对科学和诗歌的热情都是广为人知的。

也许正是这些看似不相干的追求之间的相互作用，赋予了他发明奇妙的描述性缩写词的才能，这些缩写词描述了他的许多科学贡献。这个清单很长;这里只提供了几个。其他的则散布在整个论文中。

               

a有时也被称为相间滑移分析仪

早期工作                  
                 
Spalding一生对传热传质和燃烧的贡献可以追溯到他在剑桥大学的博士论文[1]。这篇论文介绍了后来影响他职业生涯的因素。他的目标是了解液体燃料的燃烧。Spalding得出的结论是，气动流动、传热和传质都以协同的方式发挥关键作用。为了理解这些现象，Spalding随后使用了von Kármán[2]的关键流体力学概念、Kruzhilin[3]的传热概念以及Eckert和Lieblein[4]的传质概念。当然，化学动力学和流动过程的相互作用在燃烧中起着关键作用。Spalding还受益于Zeldovich和FrankKamenetskii[5,6]以及Semenov[7]的工作。                
随后，他做出了当时无法预料的预测:化学反应速率常数对燃烧没有影响，直到达到临界传质速率。他从自己的理论中推导出了这些临界速率，并得到了实验的证实。他还根据这一理论开发了一种预测火焰熄灭的方法。这项工作后来被Spalding扩展并以书的形式出版[8]。                
直到20世纪60年代初，Spalding继续研究燃烧和液滴燃烧，并做出了许多创新贡献。其中包括传质驱动力B的概念，定义如下，也称为Spalding数。这些努力在他关于对流传热和传质的杰出著作中达到高 潮[9]。此外，在此期间，他还扩大了他的努力，包括一般工程流程，其中墙体剪力起着重要作用。这将导致他杰出的洞察力和后来被称为“统一理论”的提议。在有限体积方案发展之前，对流热和质量传递问题的解决方案涉及到用速率方程取代扩散项，这是基于热/质量传递系数的概念，Spalding用g表示。因此，获得由此产生的常微分方程的数学解的能力是热和质量传递科学的核心。Spalding对传质的主要贡献是统一(概括)了这个问题，因此本质上相同的系统方法可以应用于从冷却塔到烧蚀再到燃烧的一系列学科。                
以前，在热、质量或动量传递的不同学科中，人们使用不同的方法来制定控制方程。结果，他证明了广义输运方程基本上可以应用于流体力学和热/质传递中的所有状态变量，甚至可以应用于混合分数等“合成”变量。运动方程的普遍化是工程力学问题的普遍化解方法的后续发展的先决条件。Spalding的传质公式是建立在质量分数的基础上的，而不是基于摩尔的平衡，因此在化学反应中质量(而不是摩尔)是守恒的，而且质量加权混合速度通常是在Navier-Stokes方程的解中得到的。事实上，他推导的输运方程的一般形式变成了一个自然方程，所有其他的方程都被强制纳入其中，而不适合实质性导数的额外项，扩散、产生或耗散项被集中为源(或汇)。

1960年，Spalding在本刊第1卷发表了一篇开创性的文章[10]，在这篇文章中，他假设了质量传递的“欧姆定律”，即

=gB                                                   (1)

其中                 是质量流量(“当前”)，g是“电导”或传质系数，它只取决于气动因素而不取决于热力学因素，B是传质驱动力(势)，由                 给出。和分别为流体体和壁面的某个守恒量P的值，PT表示“转移物质状态”或“t状态”。在后来的几年里，符号“P”被“                 ”所取代。对于简单的传质，是                 ，而如果发生化学反应，是                 。式(1)可以应用于焓和动量，也可以应用于化学物质。Spalding在3年后出版的关于对流传质的著作[9]中考虑了各种层流和湍流传质问题，包括传热(包括辐射)、相变和化学反应。该方法随后被许多其他人采用，至今仍是工程课程的标准方法;例如，参见[11,12]。重要的一点是，该方法适用于所有对流扩散系统，即具有守恒性质的一般系统，以及这些变量的组合。                
1960年前后，Spalding提出了一个简单但影响深远的想法[13]。从表面上看，这似乎是一个相对不重要的贡献，但它实际上强调了他将一种新的思维方向应用于他所接触的一切事物的能力。从20世纪30年代开始，研究人员一直试图提供一个连续公式来表示光滑表面附近的平均速度变化，对于简单的剪切流，u+ = f(y+)。紧挨着壁面，在粘性亚层中，变化肯定是线性的，在y+ > 30的完全湍流区域，可以接受对数变化。但是中间的区域，有时被称为缓冲层呢?当时通常的路线是对该地区的不同部分采用分段匹配。对于整个通用速度剖面，提出一个单一的解析表达式是不容易的。Spalding的贡献是将无量纲距离y+表示为量纲速度u+的函数(而不是通常将u+写成y+的形式)。对数部分可以很容易地倒置为y+ =(1/E)exp(ku +)，其中E和k是常数。为了完成这个表达式，他在右边加上层流项u+，然后从指数函数的幂级数展开中减去一些项。结果表达式在非常小和非常大的u+值下确实具有正确的渐近行为(线性和对数)，并且也很好地表示了过渡层。                                  
上述解析表达式的一个有用的结果是，我们可以得到dy+ /du+的解析表达式，然后可以很容易地将其与无量纲湍流(或有效)粘度联系起来。因此,                                  
                                               

右边的第一项代表层流的贡献，其余的是湍流的贡献。一旦用这种方式表示有效粘度，就有可能扩展这一思想来得到有效导热系数的表达式

                                               
式中                 和                 分别为层流普朗特数和湍流普朗特数，                 为有效普朗特数。                
利用这一框架，Spalding提出可以在壁面温度不连续的下游计算穿过已经建立的速度边界层(在湍流中)的热边界层。这篇文章发表在Ref.[14]，他最初将湍流普朗特数设定为统一。Kestin和Persen[15]、Gardner和Kestin[16]以及Smith和Shah[17]利用本文的关键思想，对一系列普朗特数以及壁面温度和热通量边界条件给出了广泛的数值解。
在后来的一篇论文[18]中，Spalding检查了这些解，发现了它们的规律性，为它们提供了近似的解析公式，并引入了非单位湍流普朗特数的思想。他认为= 0.9更为合适。这篇论文还有一个重要贡献。由于层流普朗特数的影响仅限于壁面附近区域，当热边界层比层流亚层更厚时，不同层流普朗特数值的解表现相同。因此，Spalding提出，层流普朗特数的唯一作用是提供“层流亚层的额外阻力”。这种阻力是由于层流普朗特数的值与湍流普朗特数的值不同。
因此对于相关的无量纲变量，可以写成                 。因此,                
               
这个表达式被认为对层流子层之外的u+值有效，在层流子层中，层流普朗特数的影响几乎消失了。这种考虑允许将积分的上限设为                 ，因为当u+的值较大时，被积函数会降为零。                
Spalding将这个积分定义为层流子层的阻力，并称之为P函数，它只取决于比值                 。Jayatilleke[19]与Spalding合作，开发了紊流中光滑壁面和粗糙壁面P函数的综合公式。这些表达中最常用的是                
               
               
图1 . 所考虑的过程类型的说明。来自参考文献。[22]与许可。                                    

对于中等普朗特数，这个表达式有时简化为。该函数至今仍被广泛用于紊流计算中提供壁面边界条件。

统一理论                
Spalding在其著作《对流传质》[9]中试图将传热和传质统一起来。他还考虑了在单一框架内同时进行的传热传质(包括蒸发、烧蚀和燃烧)。大约在同一时间，他被库塔拉泽和列昂捷夫的俄文著作所吸引，并翻译了这本书[20]。在翻译的序言中，他写道:
“以前的教科书作者都没有足够大胆地推荐在高马赫数、压力梯度和同时传质条件下计算阻力和传热的方法……作者设定了一个全面的标准，这将是对后世工作者的启发和挑战。”
Spalding本人似乎接受了这一挑战，并着手为存在压力梯度和传质的湍流剪切层(包括边界层、射流、尾迹、壁面射流和管道流)创建一种单一的计算方法。这就是后来被称为统一理论的理论。在他的整个职业生涯中，流动、热量和质量输运的统一一直是他研究的主要主题，并对CFD社区产生了重大影响。

1964年，Spalding发表了一份航空研究委员会报告[21]，标题为“湍流边界层和壁面射流中摩擦、传热和传质的统一理论”。1964年12月，这份报告被重新打印并做了注释[22]:有87页打印的傻瓜页和29张图。这是Spalding对紊流剪切流的预测作出非凡贡献的真正开端。似乎标题的范围不够宏大，论文的图1(此处转载)揭示了他的视野的范围，尽管正如他在论文的第1页所明确指出的那样，并非所有所代表的特征通常同时出现。值得注意的是，这篇论文从未在期刊上发表过，即使是以缩写形式发表。一种可能的解释是，无论是在他自己的小组内部还是在其他地方，发展都非常迅速，以至于Spalding认为统一理论很快就会被取代，不值得进一步研究。                  

① Spalding自己最近评论说，他觉得这篇论文“对一本书来说太短了，对一篇期刊论文来说也太长了”。

在p4.我们阅读新理论建立在两个主要假设之上。首先，速度、温度和浓度分布可以用两个主要组成部分的公式来描述，一个考虑动量、热量和质量传递到壁面的影响，另一个考虑与主流的相互作用;因此，一般剖面既有“边界层”成分，也有“射流”成分。其次，流体被带入壁层的方式与进入湍流射流的方式相同，并且遵循类似的定量定律。”                
有趣的是，在统一理论(UT)之前，Spalding的出版物对湍流几乎没有兴趣，尽管1961年Ricou和Spalding的论文[23]是一个值得注意的例外，该论文涉及湍流射流的携射。Spalding接着说，一旦“新”理论被构想出来，很快就会发现它的元素根本就不是新的。他建议，将统一理论视为Coles[24]的速度剖面公式Head[25]的夹带假设的结合(尽管这不是该理论的起源)可能会有所帮助。                
该理论在报告的第二节中得到了详细的阐述。给出了边界层动量方程和守恒性质方程的积分形式，并用该关系式“足够精确地”表示速度剖面。                
               

其中为无量纲速度，为自由流速度，壁面函数u+ (y+)的适当形式包含传质效应。这个方程的关键特征，与Coles提出的著名关系非常接近，是参数的出现，据说，“将在后面发挥重要作用”。引入了守恒性质的一个类似方程，它可以表示混合物中化学惰性成分的组成、停滞焓等。在平面自由湍流混合层信息的指导下，引入了初步的夹带规律。Head曾试图将夹带率与形状因子(定义为联系起来，yG是速度边界层厚度，而Spalding则展示了他通常的物理洞察力，意识到夹带率更可能与(1-)成正比，即速度剖面的射流强度或尾流分量。事实上，在最初的统一理论[21]中，他发展了一对更复杂的方程

其中为无因次夹带速率。这两个方程后来都被Spalding简化了[22]，因为他们被认为过于复杂，并且发现高估了夹带速率，并在不久之后被他进一步修改了[26]。统一理论被证明与湍流边界层理论的许多其他方面是一致的，事实上，它表明了如何发展一个统一的方法来涵盖诸如传热、传质、粗糙度甚至性质变化等影响。

发展恒定性质边界层统一理论的纯流体力学方面的责任交给了两个研究生，Escudier和Nicoll。湍流边界层（＜1）的实验数据由前者进行分析，而壁面射流（＞1）的实验数据较少因此需要大量的实验工作。夹带率是用公式推导出来的

其中，                  。虽然                  是通过将速度剖面方程拟合到测量数据来系统地估计的，但该过程中固有的不确定性是明显的。根据统一理论，编制了预测边界层和壁面射流发展的计算机程序。这些代码最初是在五孔打孔磁带上编写的，后来是在IBM卡片上编写的，并在伦敦大学的Atlas3计算机上运行。Escudier-Nicoll合作的结果在[27]中提出，其中夹带速率方程的最终版本为                
               

本文进行了大量的计算和讨论，最后得出结论:“边界层的预测与实验基本一致，除了接近分离条件的地方”和“壁面射流的预测是好的，无论是对零压力梯度还是对逆压力梯度”。另一方面，Escudier和Nicoll对逆压力梯度边界层性能差给出了通常的解释，即将实验数据不满足二维条件归咎于实验数据。这里有些道理;参见Fellouah and Pollard[28]。前面已经提到，在《统一理论》发表的第二年，Spalding修改了无量纲带射与尾迹参数的关系方程[26]。实际上，正是在后一篇论文中，Spalding首次考虑了一种基于求解湍流边界层动量积分方程和动能亏缺方程的边界层预测方法。他表明，在Rotta[29]和Truckenbrodt[30]的成熟理论中，存在带射函数和耗散积分之间的形式关系。他继续推荐改进的，但暂时的表达式，对＜1有效，包括近分离的边界层，以及对＞1有效。Spalding说:“作者目前的观点是，最终使用‘动能’法比使用‘质量守恒’法更可取。”他的理由是“更直接地将与密度变化、曲率半径和其他效应联系起来”。很明显，在Spalding的眼中，到1965年，夹带理论实际上已经成为历史:他说“不应该忘记，基于夹带理论的作用，尽管可能很短暂，但很有价值”，但没有迹象表明，消散积分方法也注定是短命的。

从统一理论到计算流体动力学                    
到1965年，UT得到了很好的建立，这是Spalding的主要关注点和他的期望，UT将普遍适用于所谓的“边界层”流动之外的非常广泛的流动。UT是建立在Taylor[31]的见解基础上的一个“大设计”，它有一个单一的理论，涵盖了边界层、尾迹和射流的典型流动。其本质是使用经典的剖面方法，这种方法已经被波尔豪森和其他人完善了。Spalding意识到，当时流行的侧写法有严重的局限性。必须设计一个能够很好地代表流体的分析剖面，并具有正确的参数数量满足初始条件和边界条件。这是一项艰巨而复杂的任务。因此，Spalding提出用分段表示通用剖面，然后从初始条件和边界条件以及其他数学约束中推导出匹配常数。虽然原则上可以使用任何剖面，但Spalding决定依赖于简单的线性或多项式段。这种分段剖面可以近似于任何描述流的理想剖面的给定精度。他和一群学生一起努力，取得了相当大的成功。他以前的学生已经确定了描述大范围流动所需的最佳夹带函数、对数定律常数和传热传质阻力。一个研究生，S.V. Patankar，在理论方面取得了积极的成果，建立了一个基于分段线性线段的通用“积分剖面”计算机代码。                    

基于动量方程和动能方程解的边界层流动流体力学方面的预测，后来通过使用温度剖面的两个参数:热边界层厚度和尾迹参数，扩展到温度方程。这些是由能量方程和“T-平方”方程得到的，这是动能方程的对应项。该研究发表在1966年芝加哥国际传热会议上[32]。这时，Spalding开始有效地使用加权函数的思想。动能方程和“T-平方”方程不需要物理上可理解的基础;它们可以被认为是动量和能量方程的积分，乘以一个加权函数(可以是因变量本身)。为了寻求更大的灵活性和更广泛的适用性，考虑了许多加权函数族。例如，利用u(或T)的幂得到了动能和相关方程。y的幂(作为加权函数)导致了各种“动量矩”方程。另一个选择是使用子域加权函数。有了这么多的积分方程，轮廓现在可以有大量的参数;它们可以是多项式，或者更好的是，分段线性。在这段时间的探索中，Spalding还提出了边界层中横流坐标的各种自变量。距离y或无量纲距离是明显的候选;但他提出了流函数(及其无量纲对应物)，它简化了横流速度或流量的表示。

Spalding当时很有信心，他的统一理论和分段剖面图可以表示许多工程兴趣流。尽管该理论取得了许多显著的进展，但具有强大的有利或不利压力梯度的流动仍然超出了UT的范围。这些流动包括重要的工业流动类别，如壁面射流和分离流动，包括管道台阶后面或腔体上方的流动。统一理论的主要问题是需要使用一组代数方程来计算轮廓参数。这些方程可以产生奇点，例如，在壁面射流到边界层的过渡中。

剪切在分离流动中起着关键作用，包括那些边界层被不利压力梯度破坏的流动，或者导致分离的几何形状。由于UT的核心依赖于Navier-Stokes方程的抛物线版本，因此这种流动不在UT的范围之内。Spalding的想法是，可以找到能够进化并代表这些流动的概况，否则这些流动就没有“自相似”行为。因此，他指派了两名新学生来解决将这种流动纳入统一理论范围的问题:沃尔夫施因已经在1964年8月加入了这个小组，Spalding让他解决平板上撞击射流的问题;在一个通道中向后台阶后面的流动和在一个驱动的方形腔内的流动是Runchal关注的焦点。这些扩展将牢固地建立统一理论，不仅适用于边界层等抛物流，而且适用于具有强压力梯度、再循环和碰撞的椭圆流。到1965年底，从Runchal和Wolfshtein的初步工作来看，UT在分离或高压梯度流动方面明显存在两个严重缺陷。UT的基本结构是抛物型Navier-Stokes方程，并包含“自相似性”存在的固有假设。Spalding意识到UT和剖面法有一个致命的缺陷。轴向扩散破坏了自相似性，并在分离和反向流动中起着关键作用，因此没有简单的方法来表示轴向扩散的作用。UT不能扩大到这种流动。在尝试将UT扩展到强压力梯度流动的过程中，很明显，当时唯一可用的解决Navier-Stokes方程的完整(椭圆)形式的一般方法是有限差分法。在某种程度上，有限差分法(FDM)与统一理论的中心思想相似:两者都使用分段多项式轮廓。但也有一个本质的区别。UT依赖于求解完成配置文件的匹配常数，而FDM使用简单的通用分段配置文件直接求解控制变量的值。本质上，在计算域中，除了感兴趣点(或节点)的近邻之外，没有明确的轮廓。这可能被认为是一个局部发展的概要，而不是一个通用的概要。

Navier-Stokes方程的有限差分方法已经存在很长时间了:Thom[33]在电子计算机出现之前就已经很好地使用了它们，Burggraf[34]最近将它们应用于方形腔体中的分离流动。很明显，使用FDM方法比修改UT方法更简单、更通用。1966年1月，Runchal和Wolfshtein结合了他们各自的研究目标，开始编写一个通用的FDM计算机程序。由于他们的直接兴趣是解决二维问题，纳维-斯托克斯方程的流函数涡度（φ-ω）形式是他们关注的焦点。φ-ω方法消除了方程中的压力，也产生了一组两个而不是三个耦合的偏微分方程。在早期，很明显存在标准的局限性，然后首选使用FDM的中心差异选项。对于高雷诺数流，解不能收敛。稳定性分析表明，这是因为在高雷诺数下，矩阵系数不再是正定的。Spalding接着做了一个类比:猪圈里吹来的风总是臭烘烘的——或者(北半球的)北风通常带来寒冷。这些讨论导致了“逆风”概念的发展和使用。后来发现，Courant等[35]等人之前也使用了单侧差异。Spalding做出的重要贡献是采用了有限体积方法(FVM)的概念作为数学发展的核心，而不是使用传统的fdm。③Spalding将计算网格比作一系列的“罐”(容纳节点)和“管”(沿着网格输送流体)。这些改进迅速带来了成功。到1966年6月，Runchal和Wolfshtein报告了壁面射流和后台阶后分离流动以及驱动盖方形腔内低和高雷诺数流动的结果。这在帝国理工学院的一系列论文中有报道[36-39]。到1966年末，这种方法已经建立起来，并开始被世界各地的其他研究人员广泛使用。它随后由 Gosman等人以书的形式出版。[40]

③ 直到一段时间后，“有限体积”一词才被采用。

φ-ω方法的一个缺点是，它基本上局限于二维流动，不能很容易地以任何明显的方式扩展到三维。使用矢量涡度和矢量速度势(与流函数对应)的可能扩展在当时被认为过于复杂;人们担心二维领域令人愉快的简单性会消失，因此进一步的进展(如下所述)是基于所谓的原始变量，即速度分量和压力。然而，涡度-矢量势公式实际上后来被Mallinson和de Vahl Davis[41]以及其他许多人成功地使用。

随着FVM的发展和成功，以及无法将UT扩展到具有强压力梯度的流动，UT的进展(基于积分方程)在1966年夏天突然停止，当时Spalding与帕坦卡尔一起前往芝加哥参加1966年国际传热会议，发表了一篇关于UT工作的论文[32]，并在暑假期间进行了一个咨询项目。在此期间，随着椭圆流的FDM已经证明的成功，他们决定基本上放弃UT的“分段剖面”方法，转而使用FDM来处理二维边界层的抛物形式的Navier-Stokes方程。他们引入了一项创新，用一个基于变量的无维流函数(本质上是von-Mises变换)代替与墙的法向距离，作为横流坐标、速度(和温度)的分段线性剖面和子域加权函数。在自由边界处使用合适的夹带公式来确定计算域向周围自由流的增长。该方法具有二维边界层、自由剪切流、壁面射流和管道流的一般程序所需的速度、灵活性、准确性和简洁性。

因此，到1966年底，椭圆型和抛物线型FDMs都牢固地建立起来，整个集团的重点转移到该技术的应用和改进上。随着论文[42,43]的发表，以及后来Patankar和Spalding[44,45]和Gosman等人[40]的著作的出版，椭圆和抛物线方法都获得了相当多的追随者，并被帝国理工学院、斯坦福大学和其他地方的研究人员广泛使用。许多博士论文将该方法应用于各种抛物线和椭圆流。这也是对湍流数学模型的研究如火如荼的时候。CFD代码为湍流和燃烧模型的测试提供了方便的工具。辐射的二通量和六通量模型[46]以及不同粒径范围的粒子浓度公式也在这些方法的框架内实现。

到1968年，越来越多的证据表明，在许多应用中，具有逆风差异的FDM和FMV会导致过度的假(数值)扩散。因此，努力转向寻找更好的方法来表示对流项。研究了许多高阶格式，包括沿流线的矢量差分概念，但没有一个证明是完全令人满意的。Spalding[47]随后提出了一种指数方法来代替上风向，但最终帝国理工学院的研究小组最终确定了一种“混合”差分方案[48]，该方案基于局部Peclet数混合了中心和上风向差分方法。在相当长的一段时间内，这[40,43]成为求解高雷诺数流动的标准方法，但后来被更高阶上旋、QUICK[49]和矢量上旋方案所扩展，并通过使用通量而不是变量的剖面。

人们意识到，“罐管”方法可以被高斯定理所取代。这避免了泰勒级数的使用，并导致在任意形状的控制体积中导出离散方程组的优雅形式。

这种广义的FVM通过产生计算流体力学(CFD)的新领域而彻底改变了工程实践。通用和灵活的计算机代码可以用来解决任何复杂几何形状和任何雷诺数下的流体流动问题。在一个重要的方面，FVM也改变了计算流体动力学背后的基本思想。与FDM和有限元法(FEM)使用的数学方法不同，对FVM的理解是通过物理方法来增强的。感兴趣的焦点不是变量及其与其他节点的相互关系参数;现在的重点是控制体积和跨越它们边界的通量，以及其中包含的源或汇。

进一步的研究和Spalding的持续创新，导致了抛物边界层方法的一个新的改进版本，由Patankar和Spalding在该书的1970年版中发表[45]。其中包含一个名为GENMIX的新计算机程序。1977年，Spalding出版了一本书[50]，本质上是一本GENMIX的用户手册。该方法现在包括一种改进的计算夹带的方法和一种确定受限流动中未知压力梯度的程序。新方法现在包含了先前为高雷诺数流动椭圆模型开发的混合格式。

Patankar-Spalding抛物边界层过程基于一种纯“行进”过程的非迭代格式。Spalding通过“汽车转向类比”解释了在受限流动中压力梯度的非迭代确定。在该方法的许多应用中，自由边界夹带的计算仍然是振荡和不稳定的来源。从形式上讲，这些方程是非线性的，迭代解将以一种直接的方式解决这个问题。(稍后在通用CFD代码PHOENICS中，Spalding[51]解决了这些问题，并在前进的步骤中使用了迭代)夹带速率(以及在受限流动中未知的压力梯度)应该被认为是方程中的特征值，由所得解的某些总体性质确定。根据我们现在所知道的，现在可以构建一个更健壮和有效的二维边界层程序;但在当今时代，几乎没有动力去开发这样的方法。

1967年，Spalding做出了一项看似微不足道的贡献，但对后来CFD工具的思考和发展产生了深远的影响。他引入了“一般输运方程”的概念来表达所有二阶对流-扩散输运方程。他建议把它们写成一个方程，用一个广义变量u来表示速度分量、焓或化学物质。这个广义输运方程有四个组成部分:(1)瞬态或累积，(2)对流，(3)扩散，(4)源项。因此，Navier - Stokes、能量和质量(物种)输运方程都可以用这个方程来表示:

这种表述的早期例子出现在Gosman等人1969年的书中[40]。这是他对数学表达“统一”的努力。这种方法导致了CFD方法的快速发展，因为人们可以集中精力解决单个方程，而不是查看每个变量的细节。

Spalding在1966年左右引入的另一项创新是用“罗盘”符号来表示FD或FV模板。这简化了通量和的可视化和表达，与他对FV的“坦克和管道”方法相吻合。围绕一个中心点(或CV)的六个点(或控制体)被指定为“北”、“南”、“东”、“西”、“高”和“低”。因此，线性代数方程的“标准”形式，即所谓的零残差形式为

到1969年，帝国理工学院团队的兴趣开始转向三维流体。流函数涡度变量在二维流动中具有明显的优势，但如前所述，在三维流动中使用它们是有问题的。在对这种流动进行了一些初步的(φ-ω)扩展之后，Spalding开始寻找直接基于Navier-Stokes方程原始形式的方法。他们尝试了许多备选办法，但没有一种完全令人满意。Caretto等人[52]的耦合SIVA方案带来了第一个成功的迹象。然而，耦合方案被认为是计算昂贵的，重点转移到寻找一种有效的顺序或分离的方法来求解控制方程。此时，Harlow及其同事[53]已经证明了交错网格算法的优势，Chorin[54]发表了关于压力分裂(或压力投影)算法的研究成果。Patankar和Spalding[55]将这些思想与早先开发的FVM技术和混合差分相结合，得出了一种鲁棒且通用的三维Navier-Stokes求解器，称为SIMPLE(压力链接方程的半隐式方法)。SIMPLE算法取得了经典的地位，并引领了新型CFD技术的商业化。                      

SIMPLE算法                      

1970年左右，有人努力开发三维边界层(或三维抛物线流)的程序。人们很快意识到，这不是二维边界层过程的简单或直接的扩展。三维边界层在顺流方向上表现为抛物型，在横流平面上表现为椭圆型。对于这种“椭圆”问题，不能使用流函数-涡度过程，因为流动的三维性排除了使用流函数作为变量。显而易见的选择是直接使用速度分量。

Spalding的突破性贡献在于发现了压力梯度的作用。为了保持该抛物型问题求解过程的行进性，Spalding要求在给定的x点上，由横跨整个边界层的单一压力梯度dp/dx驱动流向动量方程。该压力梯度可以从外部自由流的压力梯度(对于无约束流动)或总体质量守恒(对于受限流动)中获得。然后，横流速度分量由其动量方程和横流平面上满足局部连续性方程所需的压力变化来控制。虽然计算了横流平面上的压力变化，但它没有用于修改用于流向动量方程的流向压力梯度dp/dx。因此，在处理压力时故意接受了一种不一致。Spalding坚持认为，这种不一致是使用推进过程的抛物线流的好处所付出的一个小代价。在实际应用中，横流压力变化很小，不会产生明显的误差。

如果这是一个重要的问题，为什么我们在研究二维边界层时不知道呢?二维边界层的解是否也包含关于压力的不一致？Spalding指出，这个问题一直存在于二维边界层中，但仍然是不可见的;正是三维边界层把这个问题带到了风口浪尖。

在二维边界层中，我们通过求解流向动量方程获得流向速度，其中我们使用单一的dp/dx值。然后由连续性方程求出横流速度。横流动量方程根本没用!如果我们使用它，我们将能够计算横流方向的压力变化，从而暴露与使用的dp/dx值的不一致。

在三维边界层中，有两个横流速度。连续性方程本身不能决定这两个速度分量。水流在两个横流方向上的分布是由横流平面上的压力变化决定的。我们必须求解横流动量方程，并计算横流平面上的压力变化。压力的不一致是藏不住的!Patankar和Spalding[55]在文中描述了压力链接方程的半隐式方法(SIMPLE)。该方法的主要成分是:

  • 通过求解沿横流平面dp/dx值不变的沿流动量方程来计算沿流速度分量。

  • 通过求解横流动量方程和连续性计算横流速度分量和横流压力变化。

后一步需要开发一种新的算法，这是SIMPLE的核心。人们认识到我们总能解出任何特定压力场的动量方程;但是得到的速度场通常不满足连续性方程。然后，我们可以对指定的压力场提出“修正”，使得到的(修正的)速度满足连续性。这导致了所谓的“压力修正方程”的诞生。它实际上是连续性方程的一种重写形式，它提供了一种直接计算压力修正(以及压力)的方法。

尽管在1972年的论文中，Patankar和Spalding在三维抛物流的背景下提出了SIMPLE算法，但很明显，该方法中包含了二维椭圆流的程序，因为这就是处理横向流速度的方法。此外，该方法可以很容易地扩展到三维椭圆流，因为它不依赖于任何特定于二维情况的概念(如流函数)。因此，SIMPLE为求解任何维度的耦合动量和连续性方程提供了一个通用的计算程序。再一次，这种方法在世界范围内被广泛使用。该方法的改进变体，如SIMPLEC[56]、simple[57]和simple[58]后来由不同的研究人员开发，部分是为了加强速度和压力之间的耦合，并提供更快收敛的更鲁棒的算法。                      

多相流                      

在20世纪70年代中后期，Spalding开始将类似simple的方法扩展到多相流。InterPhase Slip Algorithm (IPSA)是一种两相欧拉-欧拉算法[59,60]。它是单相SIMPLE算法的自然而优雅的扩展。虽然在其最初的公式和下面的描述中，只确定了两个阶段，但所有的主要思想都很容易扩展到多相(即两个以上)流程。它首次出现于1977年在杜布罗夫尼克的ICHMT研讨会上，见Spalding[60]。

Spalding认为两相流是由两个相互穿透和混合的连续体组成的。每个位置(因此，每个有限体积电池)都可能存在两种相(例如液体和气体)。各相的量以其体积分数表征。这种两相流模型没有描述两相在局部的分布情况。体积分数表征不能区分大块液体和细小液滴。正是这种特征导致了这个公式的简单性，也导致了它的主要弱点。Spalding后来提出了消除这一弱点的扩展[61]。

对于两个互穿连续体，他提出了两组相的守恒方程。除了多了一对相间输运项外，这些方程与单相方程相似。因此，两相在某一点处的速度差会导致由于一相对另一相施加剪切力而产生的力项，两相之间的温差会导致相间传热。相间传输在很大程度上取决于可用于剪切或热传递的表面积。这反过来又要求颗粒相液滴的大小。体积分数表征不能提供这些信息。事实上，多相流问题是封闭问题的一个经典例子，其中湍流可能是工科学生最熟悉的例子。

用分离的方法求解两相的守恒方程似乎是可行的;也就是说，先求出一个相的速度和温度，然后求出另一个相的速度和温度。因此，在求解一个相位时，假设另一个相位的速度和温度值是已知的。当相间输运占主导地位时，这一过程将非常缓慢地收敛，因为另一相的假设值将主导当前相的行为。Spalding已经在另一种情况下看到了这个问题。在本文其他地方讨论的热交换器分析中，壳侧流体和管侧流体的温度以同样的方式相互影响。因此，Spalding发明了部分消除算法(PEA)来精确地处理这种情况。诀窍是用代数变换的形式来解温度方程，消除了“其他”温度的直接影响。(仍然有一些项间接包含“其他”温度的相邻值，这就是为什么转换是部分而不是全部。)减少了“其他”温度(或速度)的影响，迭代过程得到了令人印象深刻的加速。Spalding立即将两相流问题确定为PEA的候选问题。结果表明，PEA的使用对于大多数两相流计算的成功是必不可少的。还施加了“求和统一”要求，以确保在溶液过程中始终保持相体积分数之间的耦合。

如果将整个连续性方程写成质量守恒方程，则密度较大的流体(液体)的项过多地支配了密度较小的流体(气体)的项。为了克服这一困难，Spalding写了一个“体积连续性方程”[62]，该方程赋予了两相相当的权重。气液分析仪(GALA)采用体积连续性而不是质量连续性，有利于计算含有密度不连续的流体的运动，特别是那些大量级的流体，如气体和液体之间的运动。这种整体的连续性方程是建立两相流压力或压力修正方程的基础。虽然单一“共享”压力的概念是该程序的基础，但IPSA允许存在两种不同的流体压力，只要已知两种单独压力之间的适当代数关系。因此，相位连续性方程用于获得单个相的体积分数，并与SIMPLE类似，将整体连续性转换为压力校正方程。

与SIMPLE一样，IPSA迅速被世界各地的研究人员采用，IPSA的变体构成了当今几乎所有多相有限体积CFD代码中两相欧拉-欧拉公式的基础。两相流的基本框架后来使用“代数-滑移法”(有时称为混合模型或漂移-通量模型)，从概念上扩展到各种多相流。本质上，滑移代数模型由混合物的连续性和动量方程、分散相的连续性方程和滑移速度的代数关系组成。我们可以考虑不同大小范围的液滴，不同浓度的气体包，等等。随后在PHOENICS代码中实现IPSA时，相位连续性方程中也加入了相位扩散项。Spalding用梯度扩散假设对这一项进行了建模，以表示与波动速度和体积分数之间的相关性相关的湍流通量。后来一些研究者主张在动量方程中排除相扩散项，而采用湍流色散力，如Lahey等[63]。在Spalding的双流体湍流模型Malin和Spalding[64]的一些变体中，相扩散项被证明是一个重要的元素。后来(即1980年后)，Spalding利用两相概念设计了先进的湍流模型[64-67](考虑湍流是层流和湍流块的两相混合物)和先进的燃烧模型[68-70](考虑反应气体和未反应气体的混合物)。

湍流建模                        

随着二维薄剪切流的Patankar和Spalding[42,44,45]有限体积求解器的发展,需要(就湍流建模而言)转移到提供在整个计算域中确定有效湍流应力和热通量的途径，Spalding将这一问题简化为获取局部湍流粘度的问题。Patankar-Spalding方法中的初始模型是Prandtl的混合长度假设，该方案在早期的几篇论文中使用，包括他的小组提交给斯坦福会议的Kline等[71]。

然而，Spalding的设想已经在为循环流动提供相应的计算程序，在这里他认识到混合长度方案是相当不充分的。在文献[72]中，他着手解释了为什么在分离流中，壁面热流密度通常在再附着点达到最大值。正如他所说的那样[72]

“回想起来，在再附着点，(时间平均)剪应力为零，这一事实就显得特别有趣;大多数计算传热率的公式都是基于边界层中流动的物理理论，它们会预测，当剪切应力为零时，传热也必然消失。”                        

他解释说，分离流动的明显矛盾行为是由于湍流能量k向壁面的扩散输运。他通过对湍流能量方程的近似解析解，用经验系数证明了这一点，为这一论断提供了证据，从完全不同的流动类型中选择的扩散和耗散项，正确地再现了观察到的斯坦顿数和雷诺数之间的幂律依赖关系。在过去的二十年里，湍流能量方程已经被几个湍流建模者所采用，如Prandtl[73]，但Spalding的论文是第一个展示了将该方程应用于分离流的理想化，以及它如何解决了为什么分离流中的传热遵循不同于边界层中的行为规则的问题。大概正是从这项研究中，人们决定将Prandtl-Glushko湍流能量模型嵌入到稍晚的再循环流数值方案中，Gosman等[40])。

第二个主要贡献，与上述相关，涉及改进近壁对数速度剖面(外部区域的数值解固定在其上，以避免在当时不可能实现的数值要求-更不用说在椭圆求解器中求解粘性子层的建模问题)。如上所述，Spalding认识到通常的公式

以及相应的温度分布

                     
                     

在驻点附近有几面的不足，因为在那里，摩擦速度是一个不合适的速度标在证明了分离流计算中包含湍流能的必要性后，他提出将上述壁面定律公式中出现的无量纲参数改为：

                     

星号量引入了湍流能量k(在平衡状态下相同，相同，但并非相反)。

方程的改变：(12)、(13)至式。(14)和(15)只是简单地用代替。事实上，系数Cμ被选为0.09，这样，在壁面附近的平衡流动中(壁面剪切应力与“对数律”区域的湍流能量成正比)，公式给出的结果与更传统的形式相同，方程(12)、(13)。

Eq.(14)的另一个重要特征是壁面剪切应力与壁面平行速度成正比，这在流动接近分离时很重要。Spalding版本的墙律从20世纪70年代初开始在帝国理工学院的软件中使用，并最终在Launder和Spalding的评论论文中发表[74]，下文进一步参考。在过去的30年里，这种形式一直是大多数商业CFD规范的标准墙体处理方式。与原来的壁面定律不同，它总是给出至少在定性上正确的循环流动中的传热模式，尽管目前在研究级CFD中采用了更可靠的壁面函数方案。

尽管引入湍流能量来计算再循环或其他远离平衡的流动具有非常实际的好处，但Spalding意识到，由于prandtl- Glushko模型(像混合长度假设一样)需要湍流长度尺度的经验规范，因此它的适用范围极其有限。他觉得这更普遍，更完整的模型除了需要一个k的方程外，还需要一个进一步的输运方程来计算流场上湍流长度尺度的分布。受Kolmogorov[75]和Rotta[76]的工作启发，他们分别提出了波动频率ω和k与长度尺度乘积L的长度尺度方程(即kl方程)，Spalding推动了使用该方程的模型的发展，即双方程模型的发展。在早期的报告[77]中，他使用Rotta的积分形式的长度尺度方程表明，可以用同一组常数计算各种自由湍流剪切流(特别是混合层、平面射流和风扇射流)，而混合长度模型则需要不同的混合长度与流宽之比来给出正确的传播速率。然后，为了通过数值计算来支持这些发现，并为他的CFD团队创建的计算机代码开发一个通用模型，他让两名博士生在Patankar-Spalding边界层代码的帮助下，开发和测试一个使用Rotta长度尺度方程、k的输运方程和prandtel - kolmogorov涡流粘度假设的模型。Rodi负责自由流动模型的开发和测试，Ng负责近壁流动模型的开发和测试。这一研究导致了k-kL模型的建立。在Rodi和Spalding[78]中，对于自由流动，证实了该双方程模型优于混合长度模型。这项工作还提出了所谓的平面/圆射流异常，即在停滞环境中平面和圆射流不能用相同常数的模型来计算。Ng和Spalding[79]的近壁计算表明，如果在kL长度尺度方程中加入额外的一项，以产生正确的近壁长度尺度衰减，则该模型在这些流动中也表现得相当好。

在这项工作的同时，Spalding发展了他的k - W模型，在这个模型中，他使用了一个方程W ~k/L2(这是Kolmogorov频率ω的平方)作为长度尺度决定方程，其中量W被解释为涡度波动的时间平均平方。这项工作首先在参考文献[80]中报道，然后在一篇论文中，他也计算了圆形射流中的浓度波动，Spalding[81]，这篇论文成为许多其他论文的基础，特别是在燃烧界。一份内部报告Spalding[82]提供了k - W模型工作的扩展说明，该报告也出现在一本纪念册Spalding[83]中。虽然Spalding采用与学生不同的因变量似乎很奇怪，但他意识到这是一个有一系列选择的领域，需要探索不同路线的优缺点。事实上，他支持由其中一位讲师B.E. Launder领导的一个小组，Hanjalic和Launder [84]， Jones和Launder[85]探索了将耗散率方程ε ~k3/2/L作为长度尺度方程的方法。                      
Spalding在1971年关于k-W模型的报告[82]中对各种长度尺度方程进行了一些比较，并解决了选择问题。他在这份报告中表示，“还不能说任何一种模式明显优于其他模式”，而这种选择可能是基于审美考虑。然而，早在1970年，Spalding就在帝国理工学院写给同事的一份内部报告中分析了这个变量的一般长度尺度方程                      

并考虑哪个长度尺度变量的选择(即哪个指数m和n)最能满足自由流动和近壁流动的要求。他得出结论，汉贾利克和琼斯当时使用的变量ε (m = 3/2, n = -1)比当时帝国理工学院研究的其他长尺度变量更能满足要求。一个类似的分析后来发表在Launder和Spalding[74]上，为优先选择ε方程奠定了基础。然而，1970年的笔记指出，Kolmogorov的[13]变量ω，即湍流波动的频率，也可能是一个合适的变量，这一早期发现得到了最近k-ω模型(Wilcox[86])及其“剪切应力输运”(Shear-Stress Transport, SST)变体Menter[87]的支持。                      
在自由湍流剪切流中，使用不同长度尺度方程的模型产生了几乎相同的结果，正如Spalding在各种论文中所观察到的那样;但是，当时在帝国理工学院研究的模型中，只有ε方程可以免除壁效应修正(附加项或常数的变化)。这一事实导致Spalding在接下来的几年里转而偏爱ε方程。他早期参与的k-ε模型包括1972年NASA-Langley自由剪切流会议的计算[88]，帝国理工学院小组为该会议贡献了23种不同的自由剪切流的计算，使用了从混合长度假设到秒矩闭合的许多湍流模型，这些模型涉及一个尺度决定方程，全部使用ε方程，Launder等人[88]。这些计算坚定地向大部分美国观众展示了双方程模型比使用经验长度尺度处方的自由剪切流模型的优越性和更大的普遍性。在随后的几年里，Spalding也参与了用k-ε模型计算椭圆流的工作，Launder and Spalding的论文[74]对这一早期工作进行了总结，该论文长期以来一直是k-ε模型的主要参考文献之一，在文献中被大量引用。                      
在20世纪80年代，Spalding回到了他的k-W模型，并通过在w方程中加入一个涉及长度尺度梯度的源项，使相同的版本适用于近壁和自由流动。该模型成功地测试了各种流动，Malin和Spalding [64,89]， Ilegbusi和Spalding [90];但到那时，k - ε模型被广泛接受，它仍然是使用的主要模型。大约在同一时间，Spalding还引入了涉及间歇性计算的双流体(湍流/非湍流)模型。一个版本是k-ε模型的修正形式，Malin和Spalding[64]，而另一个版本，Ilegbusi和Spalding[65 - 67]，没有使用涡流粘度，而是使用基于直接计算两种流体的横向流速的扩散通量，而对于长度尺度的确定，则恢复到经验处方。虽然这些想法是新颖的，但与最初的k-ε模型的进步相比，这些模型并没有产生任何重大影响。                      
事实上，k-ε模型迅速成为使用最广泛的双方程模型，虽然现在在研究层面上被更一般的方法所取代，以确定应力-应变相互联系，但它仍然是工业CFD计算的主要工具。虽然Spalding既不是该模型的开发者，也不是该模型的主要倡导者，事实上，他后来对选择e方程的智慧表示怀疑[91]，但他在发展、推广和传播双方程涡流-黏度型模型方面发挥了重要而决定性的作用，他很早就认识到这种模型是在通用CFD计算机代码中使用的最简单的有用级别，因为这种模型不需要在计算域内规定湍流量。这使得湍流建模方面的非专业人员能够在非常广泛的领域内开始计算实际配置的湍流。
湍流燃烧和化学反应                        

Brian Spalding在紊流燃烧模拟领域做出了许多非常重要的贡献。的确，许多当前流行的关键思想在燃烧模型中已经出现在他的开创性发展中，包括混合在分子水平上的经常控制作用，以及火焰前沿表面现象所起的关键作用的问题的多尺度性质。用Spalding自己的话来说， [92]：

“很少有画家发明新的色彩，它们在画布上的排列构成了新奇。”

Spalding在紊流燃烧方面被引用最多的贡献可能是他1971年在第13届燃烧研讨会上发表的论文[93]。根据当时可用的湍流、预混、密闭、崖体稳定火焰的实验观察，他反思了火焰蔓延角度与实验条件(如燃料类型、等效比、接近流的平均速度或温度和压力水平)的相对独立性。根据这一证据，他得出结论，在很大范围的火焰条件下，化学动力学确实是辅助于空气动力学过程的。在1967年发表的一篇更早的前cfd论文中[94]，Spalding分析了假设新鲜混合物在被夹带进入剪切层时立即燃烧的含义(当时称为“夹带-燃烧”模型)。当时提供的假设导致了正确的方向;但是，从多年后第一个CFD计算提供的有利条件来看，它被认为是粗糙的。因此，在1971年的论文中，夹带燃烧模型被“混合燃烧”模型所取代，该模型后来被称为“涡流分解”(EBU)模型。EBU的基本特征是，湍流燃烧的速度是大块新鲜气体分解成小块气体的速度，从而为分子过程的作用创造了界面区域;由于燃烧主要发生在这个界面区域，在足够大的Damköhler数下，假设燃烧速率可以通过求助于涡流破裂速率来模拟。

1971年的开创性论文[92]讨论了湍流剪切流中的火焰，其中，使用平衡假设，可以估计动能衰减率等于产量。因此，在混合控制下，单位体积的反应速率设为:


式中τ为混合物的局部反应度(或反应进度)。为了适应化学动力学是速率控制过程的情况，我们将阿伦尼乌斯速率与这个混合速率结合在一个谐波平均值中。

该模型在1975年发表的一篇论文中得到进一步完善[95]，这是对燃烧室进行全面、多维建模的首次尝试之一。计算完全是椭圆的(与之前实验室火焰的主要抛物线模拟相反);紊流采用双方程k-ε模型;并建立了辐射的四通量模型。利用k和ε预测的优势，将涡流分解模型重新塑造为现在广泛使用的一种形式，即：

其中g是标量(混合分数)方差[81]。

如上所述，许多研究人员使用了几十年的涡流破裂模型(有一些次要的、非基本的变化)，Spalding在很大程度上认为它是第一个粗糙的尝试[93,95]，用于解释火焰空气动力学对燃烧速度的经常控制的影响，以及界面-表面形成和分子混合微观尺度现象的相关性。原始EBU模型在随后的研究中被Spalding和他的同事完善。因此，在[96,97]中，他概述了紊流燃烧的“滑动-拉伸”模型，其中“连贯的气体体在穿过火焰时被挤压和拉伸”。这一思想在[98]中得到了进一步阐述，并随后被赋予了ESCIMO模型的特点。

ESCIMO代表吞噬、拉伸、相干、互扩散和运动观察者，即模型的主要组成部分。Spalding借鉴了许多先前作者的思想，制定了一个模型[100]，该模型可被视为适合在CFD代码中体现的湍流燃烧的第一个多尺度解释。对于EBU，假定化学反应发生在分离反应物(扩散火焰)或反应物与生成物(预混火焰)的界面区域。界面区域被排列成“褶皱”、“包裹”、“三明治”或“层”(图2);它们形成于流体的主要剪切区;它们在对流时伸展;然后当它们被其他漩涡重新吞没时“死亡”。因此，亚太经社会的任务是双重的:第一，计算褶皱的统计数字(即褶皱大小或宽度的局部分布)，这被称为“人口部分”[100];其次，计算每个褶皱内的(一维)属性分布(这被称为“传记部分”)。第一种是欧拉的观点;第二种是拉格朗日，因为瞬态的一维轮廓依赖于(排除简化)流过流时的褶皱历史(这是模型名称的移动观察者部分)。另外，前者描述的是宏观尺度上的流动，对流和拉伸发生在宏观尺度上，而后者关注的是微观尺度上的流动，分子扩散和化学反应发生在微观尺度上。Spalding及其同事将ESCIMO模型应用于几种典型构型，如搅拌反应器、挡板火焰和自由射流[100,101]，可以通过适当的简化假设计算出折叠种群的解析解。后来，它成功地应用于湍流，非预混氢-空气火焰[100]。

尽管ESCIMO在一个框架中体现了后来在广泛接受的模型中发现的关键概念，但整个科学界对ESCIMO的接受程度并不高。这个框架使得用20世纪70年代末的计算机可以解决这个问题。因此，非预混燃烧的小火焰模型使用了薄反应层分布的相同概念，其参数化使用的不是长度尺度(折叠宽度)，而是时间尺度的倒数(耗散率)。褶皱的传记部分的时间依赖性(例如，层内的物种一维轮廓，图2)在所谓的稳态小片中被免除。这是局部计算的，与包裹(或小火焰)历史无关。Kerstein[102]最近独立导出的线性涡模型(LEM)显著地显示了ESCIMO的一些共同建模原则。值得注意的是，LEM还采用了嵌入三维计算的一维计算，以解决三维网格无法解决的小尺度混合和化学反应之间的耦合问题。与Spalding的ESCIMO不同，小尺度是在欧拉框架中计算的;通过一维场的随机搅拌来模拟湍流混合的效果，而不是采用局部的“涡流分布”。

欧拉多相方法的发展，如上文讨论的IPSA[60,103]，为湍流燃烧建模增加物理真实感开辟了新的机会。Spalding很快认识到，在火焰前缘，燃烧(燃烧)和未燃烧的物质混合在一起，就像“陆地和泻湖景观”一样[68]，可以建模为两相流，其中一相代表燃烧流体，另一相代表未燃烧流体。因此，在这个框架中，燃烧主要被视为发生在包裹边界上的界面现象，并由新鲜流体夹带到燃烧的流体中引发。燃烧速率主要由夹带速率决定，这需要建模。夹带模型通常涉及每单位体积的碎片表面积，这是长度尺度的倒数，可以解释为特征碎片尺寸;对于这种碎片大小，已经提出了几种代数和传输模型。燃烧流体内部也提供了燃烧的动力学控制。有趣的是，Magnussen[104]的单流体涡流耗散概念也有类似的概念，其中的“流体”(使用Spalding的术语)是“周围环境”或平均混合物和“精细结构”； 

                     
图2 . ESCIMO模型：  (a)的“折叠” ，来自[99]；  (b)的“滑动和拉伸”。来自参考文献。[92]，许可。

                      
在后者中，反应物在分子水平上混合。Spalding相间传递速率是涡流耗散模型中精细结构产生的速率。
由于其真正的两相特性，Spalding的双流体模型提供了允许两种不同的局部密度(分别为新鲜流体和燃烧流体)的好处。不同密度的流体碎片对共同压力梯度的反应不同，从而产生“非梯度”(有时是“反梯度”)扩散，传统模型无法再现这种扩散，但双流体模型可以[69,70]。
“未燃-燃”(或冷-热)湍流燃烧模型只是Spalding所追求的新型双流体湍流模型的一个例子。虽然"重-轻"模型，如上面概述的燃烧模型，允许表示"筛分"现象(即具有给定密度的流体碎片穿过具有另一密度的碎片)，但在其他模型中，选择密度以外的特性作为区分流体的特性，如上面的湍流建模部分所示。
建立在类似原理之上的越来越复杂的模型谱系中的最新成员是燃烧和化学反应的多流体模型(MFMs)[105,106]。MFMs通常是由Spalding作为涡断模型的演化而引入的[107]。虽然在EBU模型中只允许两种热化学状态(燃烧和未燃烧)，但在MFMs中，通过定义种群定义属性(PDAs)，例如混合物分数和/或反应性，允许任意数量的中间状态。这些pda被离散化，并且允许几种流体在域中共存，每种流体都有一组PDAs的值。多流体模型求解了每种流体局部质量分数的守恒方程。在湍流中，这些不同成分的流体碎片相互碰撞(碎片之间的碰撞最初是雷诺兹和普朗特尔湍流理论背后的一个想法);并在接触一段时间后分离，从而在界面处产生一定比例的中间流体，Spalding称之为“后代”。人口分布就这样改变了。该理论还允许流体中的“连续变化属性”，即不离散的属性。
直观开发的MFMs与用蒙特卡罗方法求解的数学推导概率密度函数(PDF)方程有许多相似之处[108,109]。流体碎片类似于流体颗粒;然而，当前者的定义属性被离散为种群定义属性时，蒙特卡罗粒子的属性被允许连续变化。在这两种方法中，当两个流体碎片碰撞时，微观混合都是建模的:在蒙特卡罗方法中，所谓的分子混合模型，通常将混合视为粒子-粒子相互作用;在MFMs中，混合是通过碰撞的流体碎片相互接触时的传质速率来表示的，从而产生“后代”。在MFMs中，PDF可以从流体种群的质量分数重建。
自20世纪60年代以来，CFD已经被应用于燃烧过程性能的预测，Spalding和帝国理工学院的团队率先将该技术用于燃烧计算。首先，Patankar和Spalding[42]设计了用于计算抛物流的具有里程碑意义的GENMIX CFD代码，1969年，Spalding将燃烧引入GENMIX[45,50]，可以模拟层流和湍流火焰射流，后来被Spalding和Stephenson用于层流火焰传播的预测[110]。Spalding在1969年提出了柴油机的CFD模型[111,112]，但早在1967年，Spalding就已经将简单化学反应系统(SCRS)概念引入到CFD代码中，用于模拟炉内的湍流扩散火焰，Pun和Spalding[113]。
该模型的一个基本特征是，燃烧过程由一个全局的、一步的、无限快速的化学反应来表示，在这个化学反应中，燃料和氧化剂以化学比例反应形成产物。中间物质被忽略，反应物、生成物和惰性物质的质量分数被表示为固定的代数关系，用一个被称为混合物分数的守恒标量变量表示。这种方法在一段时间内成为所有炉和燃烧器CFD建模的基础，很快被Spalding和同事扩展到允许燃料和氧化剂共存于同一位置，尽管在不同的时间，通过使用假定的概率密度函数(pdf)和混合分数方差的传输方程的解。因此，假设pdf模型可能始于CFD中Spalding的“双尖峰”双delta函数假设[114]，随后Lockwood和Naguib[115]提出了“剪切高斯”假设;许多人也紧随其后，包括Kent和Bilger[116]。由Spalding的帝国理工学院团队开创的第三种燃烧模型是使用具有有限反应速率的SCRS，由Arrhenius和Spalding的涡流-破裂模型的最小值表示燃料质量分数的输输方程。Khalil等人[95]在一篇论文中记录了这三种燃烧模型，并将其应用于许多炉子的布置，该论文报告了根据原子能研究机构(AERE)， Harwell (UK)的合同进行的建模工作。本文处理的是二维流动，但实际上，Patankar和Spalding已经模拟了三维炉流[117,118]，Serag-Eldin很快就发表了将其应用于三维燃气轮机燃烧室[119]。
部分抛物线流和复杂流                        

三维抛物线过程的早期应用之一是有曲率的管道。1972年，在抛物计算程序发表后不久，Pratap5[120]在Spalding的监督下，尝试计算螺旋盘管中的流动和传热。这项工作包括对发展中地区和完全发达地区的研究，完全发达地区是空间渐近状态。在热交换器、河流弯道、弯曲动脉、管道弯头等各种应用中，盘管中的流动具有重要意义。当流体流过弯曲的管道时，管道的曲率诱导二次流，在中心区域向外流动，在近壁区域向内流动。这些二次流被称为迪安涡[121,122]。在发展中区域，流向由抛物线型演变为外侧有峰的畸变型。这是动量从弯曲的内部区域转移到外部区域的结果。

在1972年以前，没有对螺旋线圈的发展区域进行过计算。Pratap等[123]发表了第一篇关于弯曲管内层流展开流动的论文。本文采用抛物线计算方法对发展中区域进行求解，并成功地与Austin的数据进行了对比[124]。随后，我们还与Mori和Nakayama[125]以及Dravid等[126]的实验数据进行了比较，分别为Pr = 4.0和0.71。本文揭示了抛物管流动程序预测工业相关复杂流动的能力。这一过程随后由层流扩展到紊流。本文用k-ε湍流模型对曲线坐标系下的动量方程和连续性方程进行了扩充[127]。

④ V.S.Pratap目前被称为S.P万卡。 

以上两个工作考虑了大的半径比。然而，当管道曲率较大时，抛物线假设就会失效。Spalding提出了一种新的过程，称为“半椭圆”或“部分抛物线”过程，它解释了加速流向流动的横向压力变化。Pratap和Spalding[128]跟进了这一想法，他们证明了捕捉这种效应的可能性。将该方法应用于强曲率导管内的流动，并采用Pratap[129]的实验数据进行比较。

绕垂直于管道轴的轴旋转的管道中的流动呈现出另一种情况，物体对流体施加力。Majumdar等[130]应用三维抛物线计算程序求解控制三维流动方程的抛物线形式。动量方程包含由于离心力和科里奥利力而产生的额外项。在随后的研究中，Majumdar和Spalding[131]将部分抛物线过程应用于旋转管道中的流动。与数据的一致性总体良好;然而，观察到某些差异。该算法认识到三维标量压力场确实是椭圆形的，并且在流动的横平面上存在流体再循环，这就需要将稀缺的计算机内存分配给压力和速度分量，但由于在主流方向上没有从下游到上游的流动，因此只需要将压力存储在内存中，因为流动基本上是抛物线流动。Pollard和Spalding[132]采用了这一思想，发展了部分椭圆方法。

Majumdar等[133]研究了圆形风管和围绕平行于风管轴线的水平轴旋转的方形风管完全发达区域内湍流流动的局部和平均换热特性。结果表明，旋转增加了换热。结果与实验数据吻合较好。Majumdar先前的一项研究[134]也考虑了绕垂直轴旋转的径向扩散器。

部分抛物线算法的发展证明了Spalding对流体力学和一般输运现象的掌握。上述关于具有体力的流动——受旋转和曲率影响的管道流动——的工作是开创性的。新开发的部分抛物线法(以及抛物线法)已成功应用。这导致了一套新的计算程序来预测复杂的三维流动，到1975年，帝国理工学院的传热部门正在考虑这种流动(“复杂流动”当时被定义为具有外部施加的应变速率，或者两个或多个基本湍流相互作用[135]，而不是“基本流动”:射流、尾迹、通道、边界层)。在此期间的研究重点之一是建立使用湍流模型预测三维湍流的能力，这在本文的其他地方进行了讨论。在许多情况下，还需要实验数据，如果设备在世界其他地方存在，就把学生送到那里。回顾1970-1980年期间该小组发表的约75篇论文，可以发现大量涉及Spalding的主题，例如:自由和强制湍流对流，Abdelmeguid和Spalding[136];通过和绕过各种形状和障碍物的流动，Caretto等[137]燃烧，Serageldin和Spalding [138];低普朗特数传热，ElHadidy等[139]，t形结流动，Pollard[140];Majumdar等[141]。许多这类研究的一个显著特点是将计算方法和模型的发展结合起来，这些方法和模型能更好地捕捉所研究现象的主要物理特征。

高如今，高性能计算机广泛流行;然而，在那些日子里，计算设备在内存和速度方面都相当有限。因此，今天适合小型手持设备的东西，当时需要庞大的基础设施和大量的支持人员。内存大小和数据存储的限制并不总是允许求解大规模问题的Navier-Stokes方程的愿望得到满足。因此，创新算法开发一直是CFD社区活动的重要组成部分。在20世纪70年代早期到中期，二维流被计算出来，但在今天看来，控制体积的数量非常少。随着具有更多内存和存储空间的新硬件的增加，以及浮点运算的发展，在三维空间甚至时间中求解方程的想法可以被考虑;然而，有了这些，就需要发明更优雅的算法。
热量交换                          
Patankar和Spalding[142]首先提出通过多孔介质或“分布阻力”类比对管壳式换热器中的流体流动和传热进行三维计算。这本质上是多尺度分析的一个非常早期的例子，其中对热交换器挡板进行了详细的计算，但阻力和传热是根据规定的体积摩擦和传热因素“建模”的。这些因素可以通过解析或详细的精细数值计算，LeFeuvre [143]， Beale和Spalding[144]，或从实验数据的相关性中获得。这种方法与传统的基于“假定流量”方法(如Delaware方法[145])的热交换器设计方法有很大不同。如前所述，在这种情况下，部分消除算法(PEA)被证明是最有利的，因为在后来的应用中(见下文)，壳侧和管侧流体的体积平均方程与两相流和传热的体积平均方程非常相似。帕坦卡-Spalding方法随后被许多其他研究人员采用和改进。关于Spalding对换热器分析的哲学方法的完整讨论可以在他的文章《换热器设计手册》[146]中找到。这个主题一直是人们感兴趣的主题之一，Spalding至今仍在继续讲授这个主题[147]。
应力分析                          

工程师通常使用什么数值方法来计算固体物体的应力和应变?通常的答案是“有限元法”(FEM)，它通常应用于结构分析问题的模拟，并形成许多计算机代码的基础。然而，有许多多物理场问题，在这些问题中，固体应力必须同时计算，例如，湍流流体的传热、传质、化学反应等等。这种流固相互作用(FSI)问题具有很大的科学意义和实际必要性。例如，旋转涡轮机中的热应力和离心应力，人体动脉中的血液流动等弹性管中流体流动中的压力波传播，或包含压电膜的打印机腔中的流动，等等。在这些和其他情况下，方法的选择不太清楚，因为与有限体积方法(FVM)相比，FEM算法对于流体动力方程耦合系统的模拟可能没有那么复杂，有限体积方法(FVM)发展得很好，目前几乎所有主要的CFD代码都采用了有限体积方法。此外，此类CFD代码已经包含了许多物理现象(湍流、辐射、化学反应、多相流等)的子组件模型，这些模型是模拟复杂FSI问题的必要组成部分。

Spalding是最早提出FVM可以很容易地用于求解弹性问题的人之一，并且，建议在基于FVM的代码中开发方程可能比将现有的大量物理流体力学和传热模型从已建立的基于FVM的CFD代码移植到基于FEM的代码中更简单。Beale和Elias[148,149]展示了如何将FVM代码PHOENICS应用于应力分析问题，利用蠕变流体流动和固体力学中的流函数和应力函数之间的类比。该方法在理论上类似于Runchal等人[43]的φ-ω方法，因此仅限于二维平面应变/应力问题。Spalding随后考虑了位移，而不是力的类比，后者在三维空间中具有相当普遍的主要优势。

线性热弹性的控制方程以平衡方程为基础，结合广义Hooke定律和应变-位移关系，得到Navier位移方程。Spalding最初将位移方程考虑为如下形式:

                         
那里                          
                         
                       
u是位移，f是体积力。因此，“压力”p与立方膨胀成正比，D=divu,其中λ和G和是lame常数，H与热膨胀系数成正比。式(19)与蠕变流的(Stokes)动量方程非常相似。因此，Spalding[150]在1993年首次提出了FSI问题的FVM解，基于位移(固体)和速度(流体)方程之间的类比，并对通常的FVM进行了以下修改:在流体区域u = (u, v, w)表示速度，而在固体区域u = (u, v, w)表示位移。在后者中，对流项被设为零。压力p、膨胀D和T之间的线性关系，如Eq.(20)，是通过包含压力和温度相关的“质量源”项引入的。应变和应力可以很容易地由位移计算出来。后来，Spalding[151]将式(19)中的压力/膨胀项替换为固相扩散项的各向异性 交换系数Γ。他还指出，其他形式的控制方程，例如基于旋转(对应于涡度的应力方程)是可能的，并且最终可能证明优于位移。对于解算法，可以使用众所周知的SIMPLE算法，但在[151]中，Spalding推测，在这种情况下，诸如SIVA算法[52]之类的耦合方案可能更有效。例如，在多孔介质中，结合固体-热流体的问题在将来很容易适用于这类分析。将应力和应变统一到流体流动和热传递与化学反应和湍流的一般框架中，是建立在物理推理基础上的通用工程方法发展的又一个合乎逻辑的步骤。
集中、热和动量有限公司                          

毫无疑问，Spalding和他的同事们在1969年创立了一家名为燃烧热能和动量有限公司(CHAM)的公司，开创了商业CFD行业。该公司于1974年更名为ConcentrationHeat and Momentum Ltd.，提供商业CFD服务，该服务是基于Spalding在帝国理工学院的研究小组提出的开创性技术。虽然CHAM的活动和CFD代码PHOENICS的起源(见下文)不可避免地交织在一起，但早在PHOENICS出现之前，CHAM基本上是世界上唯一的商业CFD公司，直到20世纪80年代初，Creare公司开始销售FLUENT代码，其早期起源可以追溯到Swithenbank等人[152];参见Runchal[153]。在1974年至1980年间，CHAM为广泛的工业领域的工业客户开发了许多特定应用的CFD计算机代码，包括航空航天，汽车，国防，化学，环境，消防和安全，船舶，制造和工艺，核电和化石燃料发电行业。这些CFD规范(其中一些进入了公共领域)也为中国商会承接CFD咨询合同提供了手段，旨在为行业解决实际问题。CHAM的客户大多是英国、欧洲和美国的大型工程公司。在此期间，在Spalding的指导和监督下，中国商会取得了许多“第一”。例如，它为核工业的蒸汽发生器创建了第一个CFD代码:为此，CHAM必须设计能够模拟同一空间内蒸汽和水的分离流动的CFD程序。另一个例子是，CHAM为铝生产行业编写了第一个电解冶炼厂(ESTER) CFD代码，用于模拟霍尔电解槽型的多阳极电解冶炼厂，其中熔融金属/电解质界面上的电磁效应和重力波可以相互作用，从而限制设备的性能。后来，ESTER被基于PHOENICS的产品Rosten所取代[154]，该产品也被行业广泛使用[155-157]。煤和燃气炉、热交换器、柴油和汽油(汽油)发动机、动力冷凝器、蒸汽发生器和自然通风冷却塔都是CHAM生产的第一批CFD软件包的设备项目。

在20世纪70年代，前面描述的部分抛物线方法为一些3D工业应用提供了比全椭圆代码更实用的路线。以部分抛物型CFD代码FLASH为例，采用曲线坐标系计算船舶和潜艇壳体的三维绕流。这项工作主要由Ishikawa-HarimaHeavy Industries (IHI)(日本)、NationalMaritime Institute (NMI) (Teddington)和AdmiraltyExperimental Works (AEW) (Haslar)资助，Spalding及其同事发表了大量出版物，例如Abdelmeguid等人[158]和Markatos等人[159]。随着计算机技术的日益强大，FLASH被内部的椭圆型CFD代码所取代，然后最终被基于PHOENICS的CFD模型所取代，该模型采用一般的非正交体拟合坐标系，Malin等[160]。部分抛物线概念在涡轮机械应用中也被证明是有用的，1976年，Singhal和Spalding[161]开发了SEMFT部分抛物线CFD代码，用于预测涡轮机械叶栅中的二维稳定湍流。这个代码是1977年提供给苏尔泽兄弟公司的，它能够处理不可压缩、亚音速、超音速或跨音速的气流。CATHY3代码是另一种部分抛物线式涡轮机械代码，提供给许多公司，包括1978年的西屋电气公司。采用曲线坐标系对旋转涡轮机道内的三维定常湍流进行了预测。该代码来源于Pratap和Spalding[128]以及Majumdar和Spalding[162]的早期工作。1979年，它被同名的全椭圆代码Malin等人[163]所取代。该椭圆型CFD代码包含一个空化模型，该软件被提供给NASA Marshall，用于离心叶轮内稳定的三维湍流预测。

如前所述，Patankar和Spalding[142]率先将分布阻力概念应用于壳管式换热器三维流动和传热的CFD建模中。自20世纪70年代初以来，该方法已被CHAM广泛用于进行热交换器，蒸汽发生器和冷凝器在国防，船舶，核能和电力工业中的应用的CFD模拟。这些应用是由许多公司资助的，包括燃烧工程 (美国) 劳斯莱斯(英国)、国家核能公司(英国)、电力研究所(美国)、欧洲原子能共同体(意大利)、北方工程工业帕森斯(英国)、通用电气公司(英国)和大瀑布(瑞典)。Spalding及其同事通过这些研究和其他研究发表了大量论文，例如Singhal等人[164,165]、Al-sanea等人[166]、Jureindini等人[167]和Malin[168]。一些应用包括通过IPSA进行均匀、代数滑移或双流体建模的两相流模拟。最终，PHOENICS成为了CHAM所有此类计算的载体，但著名的前PHOENICS CFD代码是由Spalding和其他CHAM人员开发的ATHOS,CALIPSOS, FOCS, LOBSTER, STELLA和URSULA。

许多工业多相计算的起源可以追溯到GALA[61,62]和IPSA[60,61]，两者都是由Spalding在70年代中期设计的。这些在之前的手稿中有描述，但这里应该提到的是，Spalding的杰出论文[61]也描述了不同的界面跟踪技术，作为使用GALA模拟自由表面流动的替代方案。一种方法涉及跟踪粒子，而另一种方法涉及求解液体高度的守恒方程。在20世纪70年代末和80年代初，CHAM和帝国理工学院的团队都成功地采用了界面跟踪技术，例如，参见Maxwell[169]、Awn[170]和Castrejon和Andrews[171]。IPSA对工业多相流的数值计算产生了巨大的影响。它构成了CHAM及其客户使用PHOENICS和早期CHAM代码(Moult等人[172]、Baghdadi等人[173]、Markatos[174,175]， Marchand等[176]，Al- sanea等[177]，Boisson和Malin[178]等人)进行的大部分两相流模拟的基础。不可避免地，这种方法后来在商业CFD产品中找到了竞争的途径。

Spalding满足于在许多合适的多相流应用中使用代数滑移模型。该模型在20世纪70年代后期在CHAM用于蒸汽发生器应用，Singhal等[179,180]，但在20世纪80年代，它的使用变得更加广泛，用于化学和材料加工行业。在此期间，使用了几种模型的变体，主要限于考虑重力和/或离心力场，例如:Pericleous和Drake[181](旋风中的颗粒分离);Tacke和Ludwig[182](包体分离);Pericleous和Patel[183](搅拌槽式反应器中的固体悬浮物)。最后，Spalding提出了一个更一般的公式，即使用GALA来满足体积连续性。该方法于1992年在PHOENICS实施。近年来，使用代数滑移建模的多相能力开始出现在其他通用CFD软件中。

在20世纪70年代中期到80年代初，另一个新的应用领域是建筑物中火灾和烟雾运动的预测。当时，英国Borehamwood消防研究站(FRS)资助了CHAM在2D和3D椭圆MOSIE CFD代码开发和应用方面的开创性工作，Markatos等[184]。这些特殊用途的火灾和烟雾模拟代码被基于PHOENICS的模型所取代，该模型随后被FRS和建筑研究机构(BRE)使用，以创建著名的CFD火灾模拟代码JASMINE, Markatos & Cox [185]， Kumar等[186]，该代码使用PHOENICS代码的早期版本作为其CFD引擎。通过这种方式，Spalding的贡献可以看作是为FRS和BRE提供了基本的CFD能力，使他们能够开发和验证JASMINE中包含的CFD技术，使CFD现在成为消防安全工程中公认的组成部分。

在CHAM的早期，燃烧建模是业务的重要组成部分，到1973年，该公司为Rocket Propulsion Establishment(Westcott, UK)创建了一个非常先进的椭圆CFD燃烧代码，称为BAFL，参见Tatchell和Spalding[187]和Jensen等人[188]，用于预测以超音速或亚音速移动的火箭下游区域的流动、燃烧和传热。BAFL走在时代的前面，它结合了许多复杂的数学模型，使其能够处理可压缩的、湍流的、气体之间双向化学反应的再循环流动，温度和浓度的波动，特定波长波段的辐射传热(Lockwood和Spalding[46])，以及与不同尺寸范围的颗粒物质的反应(Spalding[189])。事实上，CHAM出售给燃气轮机制造商的第一批CFD代码中，有一些是专门用于燃烧室模拟的。这些燃烧器代码中有几个已经使用了十多年，当然也得到了用户的进一步发展;并且其中至少有一个通过美国陆军进入公共领域，使竞争的CFDcode供应商能够开展业务;他们很快就照做了。1976年，CHAM向NASA Lewis研究中心提供了一个燃烧CFD代码，该代码使用化学平衡模型来表征碳氢化合物反应，并使用化学动力学模型来预测氮的氧化物，Elghobashi等[190]。这可能是预测污染物形成的第一个CFD模型。20世纪70年代另一个值得注意的CHAM燃烧规范是CORA3，它被交付给了许多工业客户。除其他外，该CFD燃烧代码利用了Spalding的“eddybreak”[93191]和“假设PDF”方法[81191]，这些方法考虑了湍流波动。Fiveland提到了Spalding的工作对CFD在发电行业的影响，参见Hanna[192]，他于20世纪70年代在美国的Babcock和Wilcox开始了CFD在锅炉制造业的研究工作:

“当时我知道Spalding教授在英国开创性的CFD工作，我们看到了使用这种‘新’CFD技术改进我们的锅炉和燃烧器设计的价值。”                             

20世纪70年代的计算机内存小，执行速度慢，因此需要巧妙的编程来解决工业问题。此外，在一些圈子里，对于计算机模拟复杂的湍流和多相流是否会有任何实际用途，人们持怀疑态度。1978年的某个时候，Spalding设想了一个能够处理所有流体流动过程的单一CFD代码的想法。因此，CHAM放弃了开发单个应用特定CFD代码的政策，并在1978年底开始创建世界上第一个通用CFD代码:PHOENICS。

PHOENICS代码                            

许多科学家和工程师在他们职业生涯的某个阶段，通过撰写他们专业领域的专著而脱颖而出。然而，Spalding将他的时间和精力投入到一个实用的CFD软件上，这反过来又帮助催生了一个完整的行业领域，并创造了前所未有的就业机会和财富。因此，CFD作为计算机辅助工程的一个子集可以直接追溯到抛物线、双曲或椭圆型数值积分代码系列PHOENICS的发展[51,193-195]。

PHOENICS于1981年作为FORTRAN-input程序商业化推出，后来在1984年作为交互式输入CFD代码推出。这是第一次对所有热流体问题使用单一的代码。这代表了一个重大的哲学变化。实际上，在那个时候，对于通用的CFD代码是可能存在的想法存在着很大的阻力。PHOENICS的成功也要归功于其他一些人，尤其是H.I. Rosten，如果不提及J.C. Ludwig自代码开始以来所做的重要和持续的贡献，那就太失职了。然而，PHOENICS本质上是Spalding的心血结晶，直到今天仍然是他关注的主题。PHOENICS代表了CFD发展的范式转变;工程师们现在能够将之前在帝国理工学院和其他地方开发的算法和封闭假设应用于大量现实世界的工业问题。一旦达到了这一步，CFD就不可逆转地走出了学术界的象牙塔，进入了工程产品和工艺的现实世界。在软件设计、内存管理/优化、可移植性、文档和培训方面，向客户提供工作软件产品的需求将专业水平提升到了一个新的水平，而这在大学环境中是不常见的。由于代码的通用性，CFD超越了热流体问题，发展到了一系列多物理场问题，例如，化学气相沉积。

在20世纪80年代早期，PHOENICS是唯一可用的通用CFD代码，它代表了Spalding小组先前发表的许多科学原理的体现，其中许多原理被精心地复 制到一个案例库中，供用户执行。为层流或湍流编写单相和多相代码，具有复杂的化学动力学，传热和传质(包括热辐射)和可变性质，这是非常独特的概念。该代码被设计为一个被称为“地球”的核心。这包含了所有已建立的硬编码和用户无法访问的实用程序(为了避免没有经验的人损坏，从而便于维护和支持)。还提供了开源组件，其中对物理模型(如湍流闭包、壁函数和求解算法)进行了编码，并且可以由用户检查、更改或替换。“钩子”(线程)在软件的不同位置提供，供用户附加自己的代码，从而允许扩展代码功能，并使其能够用于许多不是专门为其设计的问题。这种“开放源代码”的概念最早出现于1981年，后来被其他一些cfd代码供应商所采用。PHOENICS是基于有限体积法的。代码的核心是单相最简单算法[58]，SIMPLE[55,137]的一种变体。使用IPSA[60](欧拉)和GENTRA(拉格朗日)来调节多相流。用户可以从许多线性方程“求解器”中进行选择，例如逐点，平板和整个场(或提供自己的)，并在每个变量的基础上使用惯性或线性方法来控制松弛。最初采用的是交错速度方案[53]，尽管后来在直线坐标和贴体坐标中添加了并列方案[196]作为选项。边界条件和源项本身被编码为线性源项，在式(11)中，其中G是几何因子，C是源项“系数”，V是源项“值”，用户可以对其进行实质性控制;因此它们是可编程的。代码基于隔离方案，对每个字段变量执行多次“迭代”，并对整个计算域进行多次“扫描”，直到达到收敛，求解器停止。点值和残差最初是打印出来的，在后来的代码版本中，由图形用户界面(GUI)动态绘制，为用户提供对结果的信心度量。

在PHOENICS的早期版本中，用户可以根据需要使用FORTRAN在预处理器' SATELLITE '中或在运行时' GROUND '中定制代码。通过GROUND，用户实际上可以通过执行各种操作的函数和子例程访问存储在单个“F数组”中的所有EARTH变量和数组。代码的性质允许广泛的可编程性，因此如果用户愿意，可以完全绕过所有默认的求解器等等，而支持他们自己的实用程序。需要强调的是，大量的学术和工业CFD计算都是基于PHOENICS代码进行的。1985年第一届国际PHOENICS用户会议记录[197]显示，从一开始就有大量的PHOENICS应用于非常广泛的行业。此后每两年举行一次的用户会议，连同PHOENICS计算流体动力学杂志和公开文献的出版物，将工业应用的数量扩大到真正巨大的比例。Spalding的学生在帝国理工学院(Imperial College)的计算流体动力学部门(Computational Fluid Dynamics Unit)[1988 - 215]和其他地方撰写了许多基于PHOENICS的博士论文。虽然现在研究生在他们的研究工作中使用商业代码的概念几乎没有争议，但在当时这是相当新颖的。

一些代码的发展包括:贴体网格(1983年);生物流动——膀胱、肺等(1987年);可编程输入语言和共轭传热(1989);以专家为基础的系统以改善趋同(1991);固体应力/应变(1992年);并行执行(1993年);细网格嵌入(1997)。后来，引入了新的编程功能，允许用户创建自己的物理模型，而不需要诉诸FORTRAN。随着用户基础的增长，案例库也在增长。多年来，工业界和学术界用户的无数发展都被采用到代码中。此外，特殊版本的PHOENICS被开发用于广泛的应用:电解冶炼厂，化学气相沉积，电子冷却，高炉，冷却塔，外部空气动力学。中国商会成立了PHOENICS专用用户组。

PHOENICS是世界上第一个通用CFD代码;到1981年出版时，基本的科学原理已经被广泛地发表了。因此，一旦CHAM证明了通用CFD代码可以被创建并在市场上盈利，竞争对手很快就会紧随其后。PHOENICS的应用和原理被大量模仿;现在市场上有许多这样的CFD软件包，由竞争厂商提供。大多数这些代码的血统可以追溯到中国商会本身或帝国理工学院团队，中国商会和PHOENICS的影响有时受到员工流动的影响。当然，可以说，中国商会在创建和营销PHOENICS方面的创新开创了一个基于CFD软件和服务销售的全球行业，尽管没有确切的数据，但目前该行业的年营业额已超过六位数。

在第一版PHOENICS发布的几年内，Spalding在CHAM发起了基于PHOENICS的特殊用途产品(SPPs)的创建和营销工作，这些产品采用了特定于应用程序的GUIs。其中第一种是ESTER，它是专门为铝工业开发的。其他spp在20世纪80年代末和90年代初迅速跟进，例如，HOTBOX(用于电子冷却)，FLAIR(用于建筑物中的加热，通风和火灾和烟雾运动)，TACT(用于冷却塔)，PHOENICS-VWT(用于赛车空气动力学)和PHOENICS-CVD(用于化学气相沉积)。Spalding早在1979年就主张将CFD应用于CVD反应器[216]，其中包括使用嵌入式精细网格来更精确地解析晶圆附近流体流动和沉积过程的精细结构的建议。

PHOENICS-CVD, Heritage[217]仍然是SPPs中技术最先进的，能够通过模拟多组分气体中的流体流动和传热来模拟各种CVD反应器的行为，包括气相(均相)和表面(非均相)化学反应、热辐射和等离子体效应。在20世纪80年代和90年代，PHOENICS被CHAM和其他组织大量用于冶金和材料加工应用。此类应用通常具有挑战性，可能涉及流体流动、传热传质、相变、电磁学、自由表面流动、多相和化学的组合。已经提到了在铝工业中成功使用PHOENICS-ESTER。此外，Alcan International与CHAM合作，在PHOENICS中使用IPSA模拟浸入式液态金属射流的混合和凝固，参见Enright等人[218]。CHAM还使用IPSA为NASA Lewis研究中心创建了基于PHOENICS的CFD模型来模拟二元材料的固液相变化，参见Prakash[219,220]。在20世纪80年代中期，CHAM还开始开发基于PHOENICS的HIsmelt过程CFD模型，该模型涉及通过将铁和煤粉直接注入熔融浴来直接熔炼铁矿石，以及随后在顶部空间释放和燃烧熔炼气体。20世纪90年代初，中国商会的员工流动促进了Rio Tinto(澳大利亚)、Davis等人[221]对该模型的持续改进。许多其他公司也开发了基于PHOENICS的材料加工CFD模型，包括Arbed Recherche(卢森堡)、Centro Sviluppo Materiali(意大利)和Hoogevens(荷兰)。在20世纪80年代和90年代初，PHOENICS在麻省理工学院也被Szekely及其同事广泛用于材料加工操作的CFD模型的开发，在此期间发表了大量论文，例如:Ilegbusi和Szekely[222,223](中间包和连铸)，Cartwright等[223](Czochralski系统)，Choo和szkley[224](弧焊)。Saluja等[225](搅拌熔体)。直到今天，PHOENICS还在瑞典冶金研究所(Metallurgical Research Institute AB)和墨西哥国立理工学院(National Polytechnic Institute)高级研究和研究中心(Centerfor Advanced Study and Research)使用了多年，用于进行金属加工研究，例如，参见Solhed和Jonsson[226]和Rodriguez等人[227]。

在20世纪80年代早期，用于表示复杂几何形状的笛卡尔切割细胞方法在CHAM用于特定应用，但直到十年后，Spalding才开始启动一项工作计划，为PHOENICS配备一种普遍适用的设备，称为PARSOL(部分固体)。其动机是隐式地处理导入的CAD几何图形，并大大简化与使用体拟合非正交网格相关的网格生成。尽管在Spalding[228]中进行了一些讨论，但PARSOL从未正式发表。该方法对流域的大部分区域采用了背景笛卡尔网格，并对被实体切割的单元进行了特殊处理，从而保留了符合边界的网格。PARSOL也是部分非结构化的，因为它允许通过嵌入的细网格区域对局部网格进行细化。该方法计算分数面积和体积，并采用了一系列特殊算法来计算界面面积、计算壁面应力、计算固体边界附近的平流和扩散等。笛卡尔切割单元法现在已经成为结构化和非结构化体一致性网格的公认替代方案，并且现在至少在其他三种商业CFD代码中使用。

在1993年，CHAM在PHOENICS中实现了一些采用正态壁面距离阻尼函数的低雷诺数湍流模型。人们很快认识到，使用传统的基于搜索的几何程序来确定壁距分布，对于处理PHOENICS工业应用中使用的任意复杂几何图形来说，计算代价太大。就在那时，Spalding提出了一种巧妙而又极其经济的方法，即通过求解标量变量的泊松型输运方程来计算壁面距离的近代值。Spalding的建议主要是直观的，基于对两块板之间热辐射的类比。这种基于微分方程的壁距算法立即在PHOENICS中用于所有相关的湍流闭包。该方法仅在第10届布莱顿国际传热会议的海报会议上由Spalding本人正式提出，但此后它被Tucker[229,230]等人利用，改进并发表，最终该方法或其某些变体进入了大多数商业CFD代码。

壁面距离法的发明与Spalding提出的零方程低雷诺数湍流模型LVEL同时出现，该模型用于确定流体流过具有许多固体物体的空间时的湍流运动粘度[231]。在这种情况下，附近固体之间的网格密度通常太粗，任何更先进的湍流模型，如双方程k-ε模型，都无法有意义地使用。紊流运动粘度是通过Spalding壁面定律计算的[232]，该定律涵盖了整个层流和紊流状态。该模型的一个关键要素是前面提到的基于泊松的壁面距离微分方程。LVEL模型在PHOENICS中实现，并且在电子冷却应用中具有很大的实际优势，因此几乎不可避免地，该模型开始出现在其他商业CFD代码中，特别是那些致力于电子元件和系统的热设计的代码中。Spalding在他的一篇论文中写道[231]

“LVEL模型应该被视为为传热工程师的问题提供了一个实用的解决方案，而不是提供了一个新的科学见解。” 

当然，该模型的实际用途已被许多工作者证明，包括Rodgers等人[233]、Dhinsa等人[234]和Choi等人[235]。后者报告说，LVEL模型与更先进的双方程模型一样有效，用于模拟机架式服务器，其最大的优势是“可以节省大量的计算时间(三倍或更高)，特别是在进行瞬态CFD模拟或在不同机架设置下进行稳态运行时”。

Spalding认识到工业工程师需要在合理的计算时间内获得适当的解决方案，从而开发了IMMERSOL辐射模型。该模型由Spalding于1994年提出，主要用于经济地表示拥挤空间中的辐射，如电子冷却应用，以及在流体中分布许多导电固体的任何其他问题。IMMERSOL是浸入式固体(immersion SOLids)的缩写，CHAM将其应用于PHOENICS中，作为一种经济的选择，可以在光学厚度和厚度介质中表示复杂几何形状的辐射。Spalding尚未在公开文献中发表该模型，但Rasmussen[236]和Yang等人对其进行了描述。[237]。该模型涉及辐射系数的扩散方程的解，其中扩散系数与介质的吸收和散射系数以及固体壁之间的距离成正比。壁间距离是由泊松方程自动计算的，正如前面在湍流模型的壁距计算中所解释的那样。辐射方程是通过简化Hamakar和Schuster的六通量模型而设计的，参见Spalding[238]，忽略了辐射的方向效应，但将其扩展到透明介质是基于物理直觉。IMMERSOL提供了一个经济可行的近似辐射传递方程，给出了两个无限长的平行板之间的热辐射的精确解决方案;在更复杂的情况下，它的预测总是可信的，而且是正确的数量级。Yang等人[237,239]和Rasmussen[236]都报道了该模型在燃烧应用中的成功应用，Hien和Istiadji[240]也报道了该模型在太阳辐射下的建筑通风研究中的成功应用。到目前为止，除了那些使用PHOENICS的组织外，IMMERSOL似乎对工业CFD影响不大。在Rasmussen[236]的工作中发现了一个例外，他通过用户编程将IMMERSOL实现到另一个商业CFD代码STAR-CD中，然后扩展模型以处理波长依赖性。

PHOENICS站在实用工程CFD的先锋;多年来，它设定了竞争对手试图复 制的标准。不可避免地，它被模仿了，最终另一个CFD代码在销量上超过了它。今天，CFD从业者可以从许多代码套件中进行选择，其中许多代码套件是由帝国理工学院Spalding的前学生编写的，或者是由在CHAM开始其最初职业生涯并离开后组建新公司的企业家编写的。然而，PHOENICS仍然在世界范围内广泛使用，特别是在研究领域。它仍然是唯一适用于“P”(抛物线)和“H”(双曲线)选项的通用代码。实际上，Malin和Spalding[241]扩展了PHOENICS中的抛物过程，以处理释放到亚音速和超音速自由流中的未膨胀自由射流，根据S&Cthermofluid UK的合同，用于火箭排气羽流软件。这种方法产生的解比椭圆解法经济得多。到目前为止，可能没有其他商业代码具有如此广泛的可编程性。因此，PHOENICS是Spalding个人的、不断发展的、关于CFD的活的专著。
在1995年东京举行的PHOENICSCFD用户大会上，Spalding以CFD t恤标语的形式，就CFD技术从帝国理工学院的研究起源向工业分析和设计工具的转变发表了以下轻松的评论:                          
                         
F.C. LockWood，帝国理工学院团队的 原始成员，总结了Spalding在工业CFD中的先驱作用，如下：  [242]：

“Spalding能够将数学天赋与丰富的物理直觉相结合，这导致了一种科学方法，这种方法得到了在实际工业问题的“前沿”工程师的普遍赞赏，”和“商业CFD软件开发人员并不总是以他们的慷慨认可而闻名，因此有必要指出，实际上所有今天销售的商业代码都要归功于Spalding和同事的开创性工作。”

结论                            
这篇文章的准备工作对作者来说是一种热爱的劳动。这是对Brian Spalding在机械和化学工程领域做出的巨大贡献的一个小小的敬意。Spalding的物理洞察力、数学能力和非凡的职业道德使他在描述湍流、传热、传质和燃烧的数学模型的创建和发展方面取得了如此多的进步。他的作品具有高超的创造力和非凡的远见。他的思想经常在我们的理解和预测能力上产生巨大的飞跃，并为新的科学机会打开了大门。他有对给定概念、公式或程序进行概括的愿望和能力。这就是他的研究贡献产生深远影响的原因。知识的非凡深度与他的工作应用的巨大广度相匹配。除了他的技术工作，他能够通过他的清晰和清晰的演讲和写作有效地与他人交流他的想法。他的沟通技巧，加上流利的德语和俄语等语言，使他在世界各地产生了很大的影响。他的领导能力、远见、激 情和执着也为业界提供了CFD实践的多用途工具，并开发了一套用于商业应用的工具和服务。这些贡献在工程软件领域产生了深远的影响。他一直大力倡导国际合作与交流，并对他的许多学生和同事的生活和事业产生了重大影响，我们很幸运地将自己算在其中。在这里，我们试图对这位杰出的人物和他非凡的一生所取得的一些成就表示敬意。毫无疑问，未来还会有更多。


致谢                            
作者感谢Colleen King在准备和收集用于完成本文的材料方面的帮助，并感谢Micha Wolfshtein对早期版本的评论。                            
                           
                           


参考文献
[1] D.B. Spalding, The Combustion of LiquidFuels, Ph.D. thesis, Cambridge University, 1951.                            
[2] T. von Kármán, ÜberLaminare Und Turbulente Reibung, Z. Agnew. Math. Mech. 1921 (1) (1921)233–252.                            
[3] G.N. Kruzhilin, Issledovanie Teplovogo Pogranichnogo Sloia,Zh. Tekh. Fiz. 2 (1936) 183.                             
[4] E.R.G. Eckert, V. Lieblein, BerechnungDes Stoffüberganges an Einer Ebenen, Längs Angeströmten Oberflfläche Bei GroßemTeildruckgefälle, Forschung auf dem Gebiete des Ingenieurwesens 16 (2) (1949)33–42.                             
[5] Y.B. Zeldovich, V.A. Frank-Kamenetskii, Teoriia TeplovogoRasprostraneniia Plameni, Zh. Khim. Fiz. (J. Phys. Chem. (USSR)) 12 (1) (1938)100–105.                             
[6] Y.B. Zeldovich, Selected Works of Yakov Borisovich Zeldovich(A. Granik, E. Jackson, Trans.), Princeton University Press,1992–1993.                             
[7] N.N. Semenov, Teplovaia Teoriia Goreniia I Vzryvov: Iii. TeoriiaNormalnogo Rasprostraneniia Plameni, Usp. Fiz. Nauk (Progress Phys. Sci.(USSR)) 24 (4) (1940) 433–486.                             
[8] D.B. Spalding, Some Fundamentals ofCombustion, Butterworth, London, 1955.                             
[9] D.B. Spalding, Convective MassTransfer; an Introduction, Edward, Arnold, London, 1963.                             
[10] D.B.Spalding, A standard formulation of the steady convective mass transferproblem, Int. J. Heat Mass Transfer 1 (1960) 192–207.                             
[11] W.M. Kays, M.E.Crawford, B. Weigand, Convective Heat and Mass Transfer, fourth ed., McGraw-Hill,New York, 2005.                             
[12] A.F. Mills, Mass Transfer, Prentice Hall, Upper SaddleRiver, NJ, 2001.                             
[13] D.B. Spalding, A single formula for the law of thewall, ASME J. Appl. Math. 28(3) (1961) 455–458.                             
[14] D.B. Spalding, Heattransfer to a turbulent stream from a surface with a stepwise discontinuity inwall temperature, in: Proceedings of the International Developments in HeatTransfer, Part II, 1961, pp. 439–446.                             
[15] J. Kestin, L.N. Persen,Application of Schmidt’s method to the calculation of Spalding’s function andthe skin-friction coeffificient in turbulent flflow, Int. J. Heat MassTransfer 5 (1962) 143–152.                            
[16] G.O. Gardner, J. Kestin, Calculation of theSpalding function over a range of Prandtl numbers, Int. J. Heat Mass Transfer 6(1963) 289–299.                            
[17] A.G. Smith, V.L. Shah, The calculation of wall and flfluidtemperatures for the incompressible turbulent boundary layer, with arbitrarydistribution of wall heat flflux, Int. J. Heat Mass Transfer 5 (1962)1170–1189.                            
[18] D.B. Spalding, Contribution to the theory of heat transferacross a turbulent boundary layer, Int. J. Heat Mass Transfer 7 (1964)743–761.                            
[19] C.L.V. Jayatillake, The inflfluence of Prandtl number and surfaceroughness on the resistance of the laminar sub-layer to momentum and heattransfer, Prog. Heat Mass Transfer 1 (1969) 193–329.                            
[20] S.S. Kutateladze, A.I.Leontev, Turbulent Boundary Layers in Compressible Gases (D.B. Spalding,Trans.), Arnold Press, London, 1964.                            
[21] D.B. Spalding, A Unifified Theory ofFriction, Heat Transfer, and Mass Transfer in the Turbulent Boundary Layer andWall Jet, ARC No. 25 925, Her Majesty’s Stationery Offifice, London, 1964.                            
[22]D.B. Spalding, A Unifified Theory of Friction, Heat Transfer and Mass Transferin the Turbulent Boundary Layer and Wall Jet ARC CP No. 829, Her Majesty’sStationery Offifice, London, 1965.                            
[23] F.P. Ricou, D.B. Spalding, Measurementsof entrainment by axisymmetrical turbulent jets, J. Fluid Mech. 11 (1) (1961)21–32.                            
[24] D. Coles, The law of the wake in the turbulent boundary layer, J.Fluid Mech. 1 (1956) 191–226.                            
[25] M.R. Head, Entrainment in the TurbulentBoundary Layer, ARC R&M No. 3152, Her Majesty’s Stationery Offifice,London, 1960.                            
[26] D.B. Spalding, The Kinetic-Energy Defificit Equation of theTurbulent Boundary Layer, AGARDograph No. 97, Part 1, 1965, pp. 191–244.                            
[27]M.P. Escudier, W.B. Nicoll, The entrainment function in turbulent-boundarylayerand wall-jet calculations, J. Fluid Mech. 25 (2) (1966) 337–366.                            
[28] H.Fellouah, A. Pollard, Reynolds number effects within the development region ofa turbulent round free jet, Int. J. Heat Mass Transfer, in press, doi:10.1016/j.ijheatmasstransfer.2009.03.029.                            
[29]J. Rotta, Uber die Theorie der Turbulenten Grenzschichten Ström. Forsch. Nr. 1.Transl As: On the Theory of the Turbulent Boundary Layer, NACA TM 1344, 1953,1950.[30] E. Truckenbrodt, Ein Qudraturverfahren Zur Berechnung der LaminarenUnd Turbulenten Reibungschicht Bei Ebener Und Rotationssymmetrischer Strömung,Ing. Arch. 20 (1952) 211.                            
[31] B.R. Morton, G.I. Taylor, J.S. Turner, Turbulentgravitational convection from maintained and instantaneous sources, Proc. R. Soc.Lond. A 234 (1956) 1–23.                            
[32] S.V. Patankar, D.B. Spalding, A calculationprocedure for heat transfer by forced convection through two-dimensionaluniform-property turbulent boundary layers on smooth impermeable walls, in:Proceedings of the 3rd International Heat Transfer Conference, Chicago, vol.II, 1966, pp. 50–63.                            
[33] A. Thom, An Investigation of Fluid Flow in TwoDimensions, ARC R&M No. 1194, Her Majesty’s Stationery Offifice, 1928.                            
[34]O.R. Burggraf, Analytical and numerical studies of steady separated flflows, J.Fluid Mech. 24 (1966) 113–151.                            
[35] R. Courant, E. Isaacson, M. Rees, On thesolution of non-linear hyperbolic differential equations byfifinite-differences, Commun. Pure Appl. Math. 5 (1952) 243–255.                            
[36] A.K.Runchal, Three fifinite difference methods for Navier–Stokes equations, in:Proceedings of the 2nd Applied Mechanics Conference, University of Strathclyde,Glasgow, 1967.                            
[37] A.K. Runchal, D.B. Spalding, M. Wolfshtein, The NumericalSolution of the Elliptic Equations for the Transport of Vorticity, Heat andMatter in Two Dimensional Flows, SF/TN/2, Imperial College, Mech. Eng. Dept.,1967.                            
[38] A.K. Runchal, M. Wolfshtein, A Finite-Difference Procedure for theIntegration of the Navier–Stokes Equations, SF/TN/1, Imperial College, Mech.Eng. Dept., 1966.                            
[39] A.K. Runchal, M. Wolfshtein, A Fortran IV ComputerProgram for the Solution of Steady-State, Two-Dimensional Equations of Motion,Energy and Concentrations, SF/TN/10, Imperial College, Mech. Eng. Dept.,1967.                            
[40] A.D. Gosman, W.M. Pun, A.K. Runchal, D.B. Spalding, M. Wolfshtein,Heat and Mass Transfer in Recirculating Flows, Academic Press, 1969.                            
[41] G.D.Mallinson, G. de Vahl Davis, Three-Dimensional Natural Convection in a Box – ANumerical Study, J. Fluid Mech. 83 (1) (1977) 1–31.                            
[42] S.V. Patankar, D.B.Spalding, A fifinite-difference procedure for solving the equations of thetwo-dimensional boundary layer, Int. J. Heat Mass Transfer 10 (1967)1389–1412.[43] A.K. Runchal, D.B. Spalding, M. Wolfshtein, Numerical solutionof the elliptic equations for transport of vorticity, heat, and matter intwo-dimensional flflow, Phys. Fluids 12 (1969) 21–28.                            
[44] S.V. Patankar, D.B.Spalding, Heat and Mass Transfer in Boundary Layers, Morgan-Grampian, London,1967.                            
[45] S.V. Patankar, D.B. Spalding, Heat and Mass Transfer in BoundaryLayers: A General Calculation Procedure, second ed., Intertext Books, London,1970.                            
[46] F.C. Lockwood, D.B. Spalding, Prediction of Turbulent Reacting DuctFlow with Signifificant Radiation, Report No. CCK/TN/A/6, Mech. Eng. Dept.,Imperial College, London, 1971.                            
[47] D.B. Spalding, Convergence and accuracy ofthree fifinite difference schemes for a two dimensional conduction andconvection problem, Int. J. Numer. Methods Eng. 4 (1972) 540–550.                            
[48] A.K.Runchal, Convergence and accuracy of three fifinite difference schemes for atwo dimensional conduction and convection problem, Int. J. Numer. Methods Eng.4 (4) (1972) 541–550.                            
[49] B.P. Leonard, A stable and accurate convectivemodeling procedure based on quadratic upstream interpolation, Comput. MethodsAppl. Mech. Eng. 19 (1) (1979) 59–98.                            
[50] D.B. Spalding, GENMIX: A GeneralComputer Program for Two-Dimensional Parabolic Phenomena, 1st ed., PergamonPress, Oxford, New York, 1977.                            
[51] D.B. Spalding, PHOENICS a general-purposecomputer program for multidimensional one- and two-phase flflow, Math. Comput.Simul. XXIII (1981) 267–276.[52] L.S. Caretto, R.M. Curr, D.B. Spalding, Twonumerical methods for threedimensional boundary layers, Comput. Methods Appl.Mech. Eng. 1 (1) (1972) 39–57.                            
[53] F.H. Harlow, J.E. Welch, Numericalcalculation of time-dependent viscous incompressible flflow of flfluid withfree surface, Phys. Fluids 8 (12) (1965) 2182–2189.[54] A.J. Chorin, Thenumerical solution of the Navier–Stokes equations for an incompressibleflfluid, Bull. Am. Math. Soc. 73 (1967) 928–931.                            
[55] S.V. Patankar, D.B.Spalding, A calculation procedure for heat, mass and momentum transfer inthree-dimensional parabolic flflows, Int. J. Heat Mass Transfer 15 (10) (1972)1787–1806.                            
[56] J.P. Van Doormaal, J.P. Raithby, Enhancements of the SIMPLEmethod for predicting incompressible flfluid flflows, Numer. Heat Transfer BFundam. 7 (2) (1984) 147–163.                            
[57] S.V. Patankar, A calculation fortwo-dimensional elliptic problems, Numer. Heat Transfer 4 (1981) 409–426.                            
[58]D.B. Spalding, Mathematical Modelling of Fluid-Mechanics, Heat-Transfer andChemical-Reaction Processes: A Lecture Course, HTS/80/1, ComputationalFluid Dynamics Unit, Imperial College, University of London, London, 1980.                            
[59]D.B. Spalding, IPSA 1981, New Developments and Computed Results, HTS/81/ 2,Imperial College, London, 1981.                            
[60] D.B. Spalding, Numerical computation ofmulti-phase flflow and heat transfer, in: C. Taylor (Ed.), Recent Advances inNumerical Methods in Fluids, Pineridge Press, 1980, pp. 139–168.                            
[61] D.B.Spalding, The calculation of free-convection phenomena in gas–liquid mixtures,in: D.B. Spalding, N. Afgan (Eds.), Heat Transfer & Turbulent BuoyantConvection, vol. 2, Hemisphere, Washington, DC, 1977.                            
[62] D.B. Spalding, AMethod for Computing Steady and Unsteady Flows Possessing Discontinuities ofDensity, Report 910/2, CHAM Ltd., Wimbledon, London, 1974.                            
[63] R.T. Lahey, M.Lopez de Bertodano, O.C. Jones, Phase distribution in complex geometryconduits, Nucl. Eng. Des. 141 (1993) 177–201.                            
[64] M.R. Malin, D.B. Spalding, Atwo-flfluid model of turbulence and its application to heated plane jets andwakes, Physicochem. Hydrodyn. 5 (5-6) (1984) 339– 361.                            
[65] O.J. Ilegbusi, D.B.Spalding, A two-flfluid model of turbulence and its application to near-wallflflows, Physicochem. Hydrodyn. 9 (1-2) (1987) 127–160.                            
[66] O.J. Ilegbusi, D.B.Spalding, Application of a two-flfluid model of turbulence to turbulent flflowin conduits and free shear layers, Physicochem. Hydrodyn. 9 (1–2) (1987)161–181.                            
[67] O.J. Ilegbusi, D.B. Spalding, Prediction of flfluid flflow andheat transfer characteristics of turbulent shear flflows with a two-flfluidmodel of turbulence, Int. J. Heat Mass Transfer 32 (4) (1989) 767–774.                            
[68] D.B.Spalding, Towards a Two-Fluid Model of Turbulent Combustion in Gases withSpecial Reference to the Spark Ignition Engine, CFD/82/18, Computational FluidDynamics Unit, Imperial College, London, 1982.                            
[69] D.B. Spalding, Thetwo-flfluid model of turbulence applied to combustion phenomena, AIAA J. 24 (6)(1986) 876–884.                            
[70] D.B. Spalding, A turbulence model for buoyant andcombusting flflows, Int. J. Numer. Methods Eng. 24 (1) (1987) 1–23.                            
[71] S.J.Kline, M.V. Morkovin, G. Sovran, D.J. Cockrell, Computation of turbulentboundary layers, in: AFOSR-IFP-Stanford Conference, 1968.                            
[72] D.B. Spalding,Heat transfer from turbulent separated flflows, J. Fluid Mech. 27 (1967)97–109.                            
[73] L. Prandtl, Uber Ein Neues Formelsystem Für die AusgebildeteTurbulenz, Nachr. Akad. Wiss. Göttingen (1945) 6–19.                            
[74] B.E. Launder, D.B.Spalding, The numerical computation of turbulent flflows, Comput. Methods Appl.Mech. Eng. 3 (2) (1974) 269–289.                            
[75] A.N. Kolmogorov, Equations of turbulentmotion of an incompressible flfluid, Izv. Akad. Nauk. SSSR Ser. Phys. V (1–2)(1942) 56–58.                            
[76] J. Rotta, Statistische Theorie Nichthomogener Turbulenz I andII, Z. Phys. 129, 131 (1951) 547–572, 51–77 (Translated into English by W. Rodias Imperial College, Mech. Eng. Dept. Technical Notes Nos. TWF/TN/38 andTWF/TN/39).                            
[77] D.B. Spalding, The calculation of the length scale ofturbulence in some turbulent boundary layers remote from walls, in: T. Irvine(Ed.), Progress in Heat & Mass Transfer, Pergamon Press, Oxford, 1967, pp.255–266.                            
[78] W. Rodi, D.B. Spalding, A two-parameter model of turbulence, andits application to free jets, Wärme Stoffübertrag. 3 (2) (1970) 85–95.                            
[79] K.H.Ng, D.B. Spalding, Some applications of turbulence models for boundary layersnear walls, Phys. Fluids 15 (1) (1972) 20–30.                            
[80] D.B. Spalding, The Predictionof Two-Dimensional, Steady, Turbulent Elliptic Flows, Heat and Mass Transfer inTurbulent Flows with Separated Regions and Measurement Techniques, PaperDelivered at International Center for Heat and Mass Transfer Seminar, Herceg-Novi,Yugoslavia, Imperial College, Mech. Eng. Dept. Report No. EF/TN/A/16, London,1969.                            
[81] D.B. Spalding, Concentration flfluctuations in a round turbulent freejet, Chem. Eng. Sci. 26 (1) (1971) 95–107.                            
[82] D.B. Spalding, The kW Model ofTurbulence, Technical Note No. TM/TN/A16, Mech. Eng. Dept., Imperial College,London, 1971.                            
[83] D.B. Spalding, A two-equation model of turbulence,commemorative lecture for Prof. F. Bosnajakovic, VDI Forschungsheft 549 (1972)5–16.                            
[84] K. Hanjalic, B.E. Launder, A Reynolds-stress model of turbulence andits application to asymmetric shear flflows, J. Fluid Mech. 52 (1972)609–638.                            
[85] W.P. Jones, B.E. Launder, The prediction of laminarization with atwoequation model of turbulence, Int. J. Heat Mass Transfer 15 (1972)301–314.                            
[86] D.C. Wilcox, Reassessment of the scale-determining equation foradvanced turbulence models, AIAA J. 26 (1988) 1299–1310.                            
[87] F.R. Menter,Two-equation eddy-viscosity turbulence models for engineering applications,AIAA J. 32 (1994) 1598–1605.                            
[88] B.E. Launder, A. Morse, W. Rodi, D.B.Spalding, The prediction of free shear flflows – a comparison of theperformance of six turbulence models, in: Proceedings of the NASA Conference onFree Shear Flows, Langley VA, 1972.                            
[89] M.R. Malin, D.B. Spalding, Theprediction of turbulent jets and plumes by the use of the k-w model ofturbulence, Physicochem. Hydrodyn. 5 (2) (1984) 153–198.                            
[90] O.J. Ilegbusi,D.B. Spalding, An improved version of the k-w model of turbulence, ASME J. HeatTransfer 107 (1) (1985) 63–69.                            
[91] D.B. Spalding, The Vorticity-Fluctuations(kW) Model of Turbulence; Early Papers, Report CFD/82/17, Computational FluidDynamics Unit, Imperial College, London, 1982.[92] D.B. Spalding, The ESCIMOTheory of Turbulent Combustion, HTS/76/13, Department of MechanicalEngineering, Imperial College, London, 1976.                            
[93] D.B. Spalding, Mixing andchemical reaction in steady confifined turbulent flflames, in: Proceedings ofthe 13th Symposium (International) on Combustion, Pittsburgh, PA, 1971, pp. 649–657.                            
[94]D.B. Spalding, The spread of turbulent flflames in confifined ducts, in:Proceedings of the 11th Symposium (International) on Combustion, 1967, pp.807–815.                            
[95] E.E. Khalil, D.B. Spalding, J.H. Whitelaw, The calculation oflocal flflow properties in two-dimensional furnaces, Int. J. Heat Mass Transfer18 (6) (1975) 775–791.                            
[96] D.B. Spalding, Mathematical-models of turbulentflflames, Combust. Sci. Technol. 13 (1–6) (1976) 3–25.                            
[97] D.B. Spalding, Somegeneral-comments on papers concerned with theory of reactive turbulence,Combust. Sci. Technol. 13 (1–6) (1976) 245–255.                            
[98] D.B. Spalding, Developmentof the eddy-breakup model of turbulent combustion, in: Proceedings of the 16thSymposium (International) on Combustion, Pittsburgh, PA, 1977, pp. 1657–1663.                            
[99]D.B. Spalding, Computer Modelling Techniques for Laminar and TurbulentCombustion, HTS/77/12, Department of Mechanical Engineering, Imperial College,London, 1977.                            
[100] D.B. Spalding, A general theory of turbulent combustion, J.Energy 2 (1) (1978) 16–23.                            
[101] D.B. Spalding, The inflfluences of laminartransport and chemical kinetics on the time-mean reaction rate in a turbulentflflame, in: Proceedings of the 17th Symposium (International) on Combustion,1979, pp. 431–440.                            
[102] A.R. Kerstein, Linear-eddy modeling of turbulenttransport. II: Application to shear layer mixing, Combust. Flame 75 (1989)397–413.                            
[103] D.B. Spalding, Developments in the IPSA procedure for numericalcomputation of multiphase-flflow phenomena with interphase slip, unequaltemperatures, etc, in: T.M. Shih (Ed.), Proceedings of the 2nd InternationalConference on Numerical Properties and Methodologies in Heat Transfer, MarylandUniversity, 1981, pp. 421–436.                            
[104] B.F. Magnussen, On the structure ofturbulence and a generalized eddy dissipation concept for chemical reaction inturbulent flflow, in: 19th AIAA Aerospace Science Meeting, St. Louis, Missouri,1981.                            
[105] D.B. Spalding, Multi-flfluid models of turbulence progress andprospects, in: Proceedings of the CFD96 4th Annual Conference of the CFDSociety of Canada, Ottawa, 1996, pp. 27–56.                            
[106] D.B. Spalding, Turbulentmixing and chemical reaction; the multi-flfluid approach, in: InternationalSymposium on the Physics of Heat Transfer in Boiling and Condensation, Moscow,1997.                            
[107] D.B. Spalding. Connexions between the multi-flfluid and flflameletmodels of turbulent combustion, A lecture delivered at the University of Leeds,England.<http://www.cham.co.uk/phoenics/d_polis/d_lecs/mfm/flflamelet.htm>, 1998(Accessed 3 March 2008).                            
[108] C. Dopazo, Probability density function approachfor a turbulent axisymmetric heated jet centerline evolution, Phys. Fluids 22(1) (1979) 20–30.                            
[109] S.B. Pope, Monte-Carlo method for the PDF equations ofturbulent reactive flflow, Combust. Sci. Technol. 25 (1981) 159–174.                            
[110] D.B.Spalding, P.L. Stephenson, Laminar flflame propagation in hydrogen and brominemixtures, Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 324 (1558) (1971)315–337.[111] D.B. Spalding, Predicting the performance of diesel enginecombustion chambers, Proc. IMechE 184 (1969) 241–251.                            
[112] S.V. Patankar, D.B.Spalding, Mathematical models of flfluid flflow and heat transfer in furnaces,in: Proceedings of the 4th Symposium of Flames & Industry; PredictiveMethods for Industrial Flames, 1972, pp. 1–5.                            
                           
[113] W.M. Pun, D.B. Spalding, Aprocedure for predicting the velocity and temperature distribution in aconfifined, steady, turbulent, gaseous diffusion flflame, in: Proceedings ofthe International Astronautics Federation Meeting, Belgrade, 1967, pp.3–21.                            
[114] D.B. Spalding, Turbulent, Physically-Controlled CombustionProcesses, Report RF/TN/A/4, Mech. Eng. Dept. Imperial College, University ofLondon, London, 1971.                            
[115] F.C. Lockwood, A.S. Naguib, The prediction of theflfluctuations in the properties of free, round-jet, turbulent, diffusionflflames, Combust. Flame 24 (1975) 109–124.                            
[116] J.H. Kent, R.W. Bilger, Theprediction of turbulent diffusion flflame fifields and nitric oxide formation,in: Proceedings of the 16th Symposium (International) on Combustion Pittsburgh,PA, 1976, pp. 1643–1656.                            
[117] S.V. Patankar, D.B. Spalding, A computer modelfor three-dimensional flflow in furnaces, in: Proceedings of the 14thInternational Symposium on Combustion, 1972, pp. 605–614.                            
[118] S.V. Patankar,D.B. Spalding, Simultaneous predictions of flflow pattern and radiation forthree-dimensional flflames, in: N. Afgan, J. Beer (Eds.), Heat Transfer inFlames, Hemisphere, Washington, 1974, pp. 73–94.                            
[119] M.A. Serag-Eldin, TheNumerical Prediction of the Flow and Combustion Processes in aThree-Dimensional Combustion Chamber, Ph.D. thesis, Imperial College, London,1977.                            
[120] V.S. Pratap, Flow and Heat Transfer in Curved Ducts, Ph.D. thesis,Imperial College, University of London, London, 1975.                            
[121] W.R. Dean, Note onthe motion of flfluid in a curved pipe, Phil. Mag. 14 (1927) 208–223.                            
[122] W.R.Dean, Stream-line motion of a flfluid in a curved pipe, Phil. Mag. 5 (1928)673–695.                            
[123] S.V. Patankar, V.S. Pratap, D.B. Spalding, Prediction of laminarflflow and heat transfer in helically coiled pipes, J. Fluid Mech. 62 (1974)539–551.                            
[124] L.R. Austin, The Development of Viscous Flow within HelicalCoils, Ph.D. thesis, University of Utah, 1971.                            
[125] Y. Mori, W. Nakayama, Studyon forced convective heat transfer in curved pipes (1st report, laminarregion), Int. J. Heat Mass Transfer 8 (1965) 67–82.                            
[126] A.N. Dravid, K.A.Smith, E.W. Merrill, P.L.T. Brian, Effect of secondary flfluid motion onlaminar flflows and heat transfer in helically coiled tubes, AIChE J. 17 (1971)1017–1122.                            
[127] S.V. Patankar, V.S. Pratap, D.B. Spalding, Prediction ofturbulent flflow in curved pipes, J. Fluid Mech. 67 (1975) 583–595.                            
[128] V.S.Pratap, D.B. Spalding, Fluid flflow and heat transfer in three-dimensional ductflflows, Int. J. Heat Mass Transfer 19 (10) (1976) 1183–1188.                            
[129] V.S. Pratap,D.B. Spalding, Numerical computations of flflow in curved ducts, Aeronaut. Q.XXVI (1975) 219–228.                            
[130] A.K. Majumdar, V.S. Pratap, D.B. Spalding, Numericalcomputation of flflow in rotating ducts, ASME J. Fluids Eng. 99 (1977)148153.                            
[131] A.K. Majumdar, D.B. Spalding, A numerical investigation ofthree-dimensional flflows in a rotating duct by a partially-parabolicprocedure, in: Proceedings of the ASME Paper No. 77-WA/FE-7, 1977.                            
[132] A.Pollard, D.B. Spalding, The Prediction of the three-dimensional turbulentflflow fifield in a flflow splitting tee-junction, Comput. Math. Appl. Mech.Eng. 13 (1978) 293–306.                            
[133] A.K. Majumdar, W.D. Morris, D. Skiadaressis, D.B.Spalding, Heat transfer in rotating ducts, Mech. Eng. Bull. (Durgapur) 8 (4)(1978) 87–95.                            
[134] A.K. Majumdar, D.B. Spalding, Numerical investigation offlflow in rotating radial diffusers, Mech. Eng. Bull. (Durgapur) 8 (1–2) (1977)19–28.                            
[135] P. Bradshaw. Complex turbulent flflows, in: W.T. Koiter (Eds.),Proceedings of the 14th Iutam Congress, Theoretical and Applied Mechanics,North Holland, Delft, 1976, pp. 103–113.                            
[136] A.M. Abdelmeguid, D.B. Spalding,Turbulent flflow and heat transfer in pipes with buoyancy effects, J. FluidMech. 94 (2) (1979) 383–400.                            
[137] L.S. Caretto, A.D. Gosman, S.V. Patankar,D.B. Spalding, Two Calculation procedures for steady, three-dimensional flflowswith recirculation, in:Proceedings of the 3rd International Conference onNumerical Methods inFluid Mechanics, vol. 2, 1973, pp. 60–68.                            
[138] M.A.Serageldin, D.B. Spalding, Computations of three-dimensional gasturbinecombustion chamber flflows, Trans. ASME J. Eng. Power 101 (3) (1979)326–336.                            
[139] F. El-Hadidy, F. Gori, D.B. Spalding, Further results on the heattransfer to lowPrandtl number flfluids in pipes, Numer. Heat Transfer A Appl. 5(1) (1982) 107–117.[140] A. Pollard, The numerical calculation of partiallyelliptic flflows, Numer. Heat Transfer A Appl. 2 (3) (1979) 267–278.                            
[141] A.K.Majumdar, A.K. Singhal, D.B. Spalding, Numerical modeling of wet coolingtowers. I. Mathematical and physical models, ASME J. Heat Transfer 105 (4)(1983) 728–735.                            
[142] S.V. Patankar, D.B. Spalding, A calculation procedure forthe transient and steady-state behavior of shell-and-tube heat exchangers, in:N. Afgan, E.U. Schlünder (Eds.), Heat Exchangers: Design and Theory Sourcebook,Scripta Book Company, Washington, DC, 1974, pp. 155–176.                            
[143] R.F. Le Feuvre,Laminar and Turbulent Forced Convection Processes through in-Line Tube Banks,Ph.D. thesis, Imperial College, University of London, 1973.                            
[144] S.B. Beale, D.B. Spalding, Numericalstudy of flfluid flflow and heat transfer in tube banks with stream-wiseperiodic boundary conditions, Trans. CSME 22 (4A) (1998) 394–416.                            
[145] R.W.Serth, Process Heat Transfer: Principles and Applications, Academic Press,2007.                            
[146] G.F. Hewitt (Ed.), Heat Exchanger Design Handbook, vol. 1, BegellHouse, New York, 2008.                            
[147] D.B. Spalding, Numbers count: SFT analysis ofheat-transfer equipment. The James P. Hartnett Lecture (Paper KN-15), in: G. deVahl Davis, E. Leonardi (Eds.), Proceedings of the IHTC-13: International HeatTransfer Conference, Sydney, 2006.                            
[148] S.B. Beale, S.R. Elias, Stressdistribution in a plate subject to uniaxial loading, PHOENICS J. Comput. FluidDyn. Appl. 3 (3) (1990) 255–287.                            
[149] S.B. Beale, S.R. Elias, NumericalSolution of Two-Dimensional Elasticity Problems by Means of a SIMPLE-BasedFinite-Difference Scheme, TR-LT-020, National Research Council, Ottawa,1991.                            
[150] D.B. Spalding, Simulation of flfluid flflow, heat transfer and soliddeformation simultaneously, in: NAFEMS Conference No. 4, Brighton, 1993.                            
[151]D.B. Spalding, Extending the frontiers of computational heat transfer, PaperCht-08-001, in: Proceedings of the CHT-08 ICHMT International Symposium onAdvances in Computational Heat Transfer, Marrakech, 2008.                            
[152] J. Swithenbank,A. Turan, P.G. Felton, D.B. Spalding, Fundamental Modelling of Mixing,Evaporation and Kinetics in Gas Turbine Combustors, Report No. HIC 320,Department of Chemical Engineering & Fuel Technology, University ofSheffifield, 1979.[153] A.K. Runchal, Brian spalding: CFD & reality, Int.J. Heat Mass Transfer, in press,doi:10.1016/j.ijheatmasstransfer.2009.03.058.                            
[154] H.I. Rosten, TheMathematical Foundation of the ESTER Computer Code, CHAM TR/84, Wimbledon,1982.                            
[155] W.E. Wahnsiedler, Hydrodynamic modeling of commercial Hall-Heroultcells, Light Metals (1987) 269–287.                            
[156] V. Potocnik, Modeling of metal–bathinterface waves in Hall-Heroult cells using ESTER/PHOENICS, in: Proceedings ofthe Light Metals, 1989, pp. 227– 235.                            
[157] M. Segatz, D. Vogelsang, C. Droste,P. Baekler, Modeling of transient magnetohydrodynamic phenomena in Hall-Heroultcells, Light Metals (1993) 361– 368.[158] A.M. Abdelmeguid, N.C. Markatos, K.A.Muraoka, D.B. Spalding, A method of predicting three-dimensional turbulentflflows around ships’ Hulls, in: Proceedings of the International Symposium onShip Viscous Resistance, SSPA, Gothenborg, Sweden, vol. 3, 1978, pp. 1–24.                            
[159]N.C. Markatos, M.R. Malin, D.G. Tatchell, Computer analysis of threedimensionalturbulent flflows around ships’ Hulls, Proc. IMechE 194 (24) (1980)239–248.                            
[160] M.R. Malin, H.I. Rosten, D.B. Spalding, D.G. Tatchell,Application of PHOENICS to flflow around ship’s Hulls, in: L. Larsson, C. S.Ohlsson (Eds.), Proceedings of the 2nd International Symposium on Ship ViscousResistance, Gothenburg, Sweden, 1985, pp. 16:1–16:23.                            
[161] A.K. Singhal, D.B.Spalding, A Two-Dimensional Partially-Parabolic Procedure for Axial-FlowTurbomachinery Cascades, ARC R&M No. 3807, Her Majesty’s StationeryOffifice, London, 1978.                            
[162] A.K. Majumdar, D.B. Spalding, A ParametricComputation of Flows and Performance of Rotating Radial Impellers, ARC R&M3810, Her Majesty’s Stationery Offifice, London, 1977.                            
[163] M.R. Malin, H.I.Rosten, D.G. Tatchell, Three-dimensional computations of flflows in centrifugalpumps and compressors; in performance prediction of centrifugal pumps andcompressors, in: Proceedings of the 25th Annual International Gas TurbineConference and Exhibit and 22nd Annual Fluids Engineering Conference, NewOrleans, LA, 1979, pp. 33–45.                            
[164] A.K. Singhal, D.B. Spalding, Mathematicalmodelling of multi-phase flflow and heat transfer in steam generators, in:Proceedings of the 2nd Multi-phase Flow Transport: Fundamentals, Reactor-SafetyApplications, 1980, pp. 373– 406.                            
[165] A.K. Singhal, L.W. Keeton, A.J. Prekwas,J.S. Weems, D.B. Spalding, ATHOS – A Computer Program for Thermal HydraulicAnalysis of Steam Generators, 1982.                            
[166] S.A. Al-Sanea, N. Rhodes, D.G.Tatchell, T.S. Wilkinson, A computer model for detailed calculation of theflflow in power-station condensers, in: Proceedings I. Chem. E. SymposiumSeries No. 75, Condensers: Theory & Practice, 1983, pp. 70–88.                            
[167] R.H.Jureidini, M.R. Malin, M.J. Lord, K.K. Yau, A three-dimensional model for powercondenser design, in: Proceedings of the 4th International PHOENICS UserConference, Miami, FL, 1991.                            
[168] M.R. Malin, Modelling flflow in anexperimental marine condenser, Int. Commun. Heat Mass Transfer 24 (5) (1997)597–608.                            
[169] T.T. Maxwell, Numerical Modelling of Free-Surface Flows, Ph.D.thesis, Imperial College, University of London, 1977.                            
[170] A.G. Awn, NumericalModelling of Flows with Moving Interfaces, Ph.D. thesis, Imperial College,University of London, London, 1979.                            
[171] A. Castrejon, M.A. Andrews. Aprocedure for calculating moving-interface flflows with PHOENICS, in: C.A.Brebbia, S.A. Orszag (Eds.), Proceedings of the 1st International PHOENICSUsers Conference, Dartford, Published in Lecture Notes in Engineering, 18(Springer-Verlag, London, 1985), pp. 433–443.                            
[172] A. Moult, V.S. Pratap, D.B.Spalding, Finite-difference predictions of waves in stratifified gas–liquidflflows, in: Proceedings of the 1st International Conference on PhysicochemicalHydrodynamics, Washington, DC, 1978.                            
[173] A.H.A. Baghdadi, H.I. Rosten, A.K.Singhal, D.B. Spalding, D.G. Tatchell, Numerical modelling of two-phase flflowwith slip, in: Proceedings of the ICHMT seminar on two-phase flflows,Dubrovnik, Yugoslavia, 1978.                            
[174] N.C. Markatos, A. Moult, P.J. Phelps, D.B.Spalding, The calculation of steady, three-dimensional, two-phase flflow andheat transfer in steam generators, in: N. Afgan, F. Durst, G. Tsiklauri (Eds.),Proceedings of the ICHMT Seminar on Two-Phase Flows, Dubrovnik, Yugoslavia,vol. 1, 1978, pp. 485–502.                            
[175] N.C. Markatos, S.M. Rawnsley, D.B. Spalding,Heat transfer during a smallbreak loss-of-coolant accident in a pressurizedwater reactor-a parametric study for a 4-in lower-plenum break, Int. J. HeatMass Transfer 27 (8) (1984) 1379–1394.                            
[176] E.O. Marchand, A.K. Singhal, D.B.Spalding, Numerical modelling of two-phase flflow in steam generators, in: 2ndMulti-phase Flow & Heat Transfer Symposium, Miami, FL, 1978.                            
[177] S.A.Al-Sanea, N. Rhodes, T.S. Wilkinson, Mathematical modelling of twophasecondenser flflows, in: Proceedings of the BHRA 2nd International Conference onMulti-Phase Flow, London, 1985, pp. 169–182.                            
[178] N. Boisson, M.R. Malin,Numerical prediction of two-phase flflow in bubble columns, Int. J. Numer.Methods Fluids 23 (1996) 1289–1310.                            
[179] A.K. Singhal, D.B. Spalding,Mathematical Formulation and Solution Procedure of URSULA2 Computer Program,H3201/1, CHAM, 1979, 1978.                            
[180] A.K. Singhal, L.W. Keeton, D.B. Spalding,Predictions of thermal hydraulics of a PWR steam generator by using thehomogeneous, the two-flfluid, and the algebraic-slip model, in: Proceedings ofthe AIChE Symposium Series No.199, vol. 76, 1980, pp. 45–55.[181] K.A.Pericleous, S.N. Drake, An algebraic-slip model for multi-phase applications,numerical simulation of flfluid flflow and heat/mass transfer processes, in:C.A. Brebbia,S.A. Orszag (Eds.), Lecture Notes in Engineering, Springer-Verlag, London, 1986.[182] K.H. Tacke, J.C. Ludwig, Steel flflow and inclusionseparation in continuous casting tundishes, Steel Res. 58 (1987) 262–270.                            
[183]K.A. Pericleous, M.K. Patel, The modelling of tangential and axial agitators inchemical reactors, Physicochem. Hydodyn. 8 (2) (1987) 105–123.                            
[184] N.C.Markatos, M.R. Malin, G. Cox, Mathematical modelling of buoyancyinduced smokeflflow in enclosures, Int. J. Heat Mass Transfer 25 (1982) 63–75.                            
[185] N.C.Markatos, G. Cox, Hydrodynamics & heat transfer in enclosures containing afifire source, Physicochem. Hydrodyn. 5 (1) (1984) 53–66.                            
[186] S. Kumar, N.Hoffmann, G. Cox, Some validation of JASMINE for fifires in hospital wards, in:C.A. Brebbia, S.A. Orszag (Eds.), Lecture Notes in Engineering,Springer-Verlag, 1986, pp. 159–169.                            
[187] D.G. Tatchell, D.B. Spalding, TheRocket Based-Flow Program – BAFL, CHAM Technical Report 640/1, CHAM Ltd.,Wimbledon, 1973.                            
[188] D.E. Jensen, D.B. Spalding, D.G. Tatchell, A.S. Wilson,Computation of structures of flflames with recirculating flflow and radialpressure gradients, Combust. Flame 34 (1979) 309–326.                            
[189] D.B. Spalding, Noteson the PARMIX Program, Report BL/TN/A/35, Imperial College, Mech. Eng. Dept.,London, 1970.                            
[190] S. Elghobashi, D.B. Spalding, S.K. Srivatsa, Prediction ofHydrodynamics and Chemistry of Confifined Turbulent Methane–Air Mixtures withAttention to the Formation of Oxides of Nitrogen, NASA CR-135170, CHAM TR951/2, 1977.                            
[191] D.B. Spalding, Turbulence Models & Their ExperimentalVerifification: – 12. Turbulence Effects on Combustion, Report HTS/73/27,Imperial College, London, 1973.[192] K. Hanna, Visions of the Future: CFD inthe Power Generation Industry, FLUENT News (Fall Edition), 2004.                            
[193] D.B.Spalding, PHOENICS a General-Purpose Computer Program for MultiDimensional One-and Two-Phase Flow, HTS/81/11, Computational Fluid Dynamics Unit, ImperialCollege, University of London, 1981.                            
[194] D.B. Spalding, Four Lectures on thePHOENICS Code, CFD/82/5, Computational Fluid Dynamics Unit, Imperial College,University of London, 1982.                            
[195] D.B. Spalding, PHOENICS 1984. AMulti-Dimensional Multi-Phase GeneralPurpose Computer Simulator for Fluid Flow,Heat Transfer and Combustion, CFD/84/18, Computational Fluid Dynamics Unit,Imperial College, University of London, 1984.                            
[196] C.M. Rhie, W.L. Chow,Numerical study of the turbulent flflow past an airfoil with trailing edgeseparation, AIAA J. 21 (11) (1983) 1525–1532.                            
[197] N.C. Markatos, D.G.Tatchell, M. Cross, N. Rhodes, Numerical simulation of flfluid flflow andheat/mass transfer processes, in: C.A. Brebbia, S.A. Orszag (Eds.), Proceedingsof the 1st International PHOENICS Users Conference, Dartford, London (1985),Published in Lecture Notes in Engineering, vol. 18, Springer-Verlag, 1986.                            
[198]G.W. Carroll, Experimental and Theoretical Investigation of ThreeDimensionalTurbulent Mixing in Jets & Ducts, Ph.D. thesis, Imperial College,University of London, London, 1980.                            
[199] P.A. Tennant, A Numerical Analysis ofthe Inflfluence of Internal Nozzle Geometry on Diesel Fuel Injection, Ph.D.thesis, Imperial College, University of London, 1981.                            
[200] J.O. Ilegebusi, ARevised Two-Equation Model of Turbulence, Ph.D. thesis, Imperial College,University of London, 1983.                            
[201] A. Castrejon, Unsteady Free-Convection HeatTransfer inside and Annular Cavity, Ph.D. thesis, Imperial College, Universityof London, 1984.                            
[202] H.Q. Qin, Computer and Water-Model Simulation of Flowthrough Poppet Valves, Ph.D. thesis, Imperial College, University of London,1984.                            
[203] S.M. Lo, Development of a New Interphase Friction Law for Air–WaterFlows in Inclined Pipes, Ph.D. thesis, Imperial College, University of London,1985.                            
[204] G.T. Robert, Shock Tube of Convective Heat Transfer, Ph.D. thesis,Imperial College, University of London, 1985.                            
[205] M.J. Andrews, TurbulentMixing by Rayleigh–Taylor Instability, Ph.D. thesis, Imperial College,University of London, 1986.                            
[206] L. Jun, Computer Modelling of Flow with FreeSurface, Ph.D. thesis, Imperial College, University of London, 1986.                            
[207] M.R.Malin, Turbulence Modelling for Flow and Heat Transfer in Jets, Wakes andPlumes, Ph.D. thesis, Imperial College, University of London, 1986.                            
[208] J.Z.Y.Wu, Studies in Unsteady Flame Propagation, Ph.D. thesis, Imperial College,University of London, 1986.                            
[209] S.T. Xi, Transient Turbulent Jets of Miscibleand Immiscible Fluids, Ph.D. thesis, Imperial College, University of London,1986.                            
[210] K.A. Fenech, Analysing Aspects of the Performance of an Iron BlastFurnace, Ph.D. thesis, Imperial College, University of London, 1987.                            
[211] C.A.Miller, An Evaluation of the Performance of the PHOENICS Code, Ph.D. thesis,Imperial College, University of London, 1988.                            
[212] D. Radosavljevik, TheNumerical Simulation of Direct-Contact NaturalDraught Cooling Tower Performanceunder the Inflfluence of Cross-Wind, Ph.D. thesis, Imperial College, Universityof London, 1988.                            
[213] F. Villasenor, The Kelvin–Helholtz Instability in aStratifified Shear Flow, Ph.D. thesis, Imperial College, University of London,1988.                            
[214] N. Fueyo, Two-Fluid Models of Turbulence for Axisymmetrical Jets andSprays, Ph.D. thesis, Imperial College, University of London, 1989.                            
[215] S.B.Beale, Fluid Flow and Heat Transfer in Tube Banks, Ph.D. thesis, ImperialCollege, University of London, 1993.                            
[216] D.B. Spalding, Numerical Modelling ofCVD Reactors, HTS/79/6, Mech. Eng. Dept., Imperial College, University ofLondon, London, 1979.                            
[217] J.R. Heritage, PHOENICS-CVD: a code for the designand development of chemical vapour deposition equipment and processes, PHOENICSJ. Comput. Fluid Dyn. Appl. 7 (4) (1994) 165–180.                            
[218] P.G. Enright, L.Katgerman, J.C. Ludwig, S.C. Rogers, Mixing and solidifification of a turbulentliquid jet in a co-flflowing stream, Int. J. Numer. Methods Eng.24 (1987)231–249.                            
[219] C. Prakash, Two-phase model for binary solid-liquid phase change,part I, Numer. Heat Transfer B Fundam. 18 (1990) 131–145.                            
[220] C. Prakash,Two-phase model for binary solid-liquid phase change, part II, Numer. HeatTransfer B Fundam. 18 (1990) 147–167.                            
[221] M.P. Davis, R.J. Dry, M.P. Schwarz,Flow simulation of the HIsmelt process, in: Proceedings of the 3rdInternational Conference on CFD in the Minerals & Process Industries,Melbourne, 2003, pp. 306–311.                            
[222] O.J. Ilegbusi, J. Szekely, Fluid Flow andtracer dispersion in shallow tundishes, Steel Res. 59 (9) (1988) 399–405.                            
[223]R. Cartwright, J.O. Ilegbusi, J. Szekely, A comparison of order-of-magnitudeand numerical analyses of flflow phenomena in Czochralski and magneticCzochralski systems, J. Crystal Growth 94 (1989) 321–333.                            
[224] R.T.C. Choo, J.Szekely, Vaporization kinetics and surface temperature in a mutually coupledspot gas tungsten arc weld and weld pool, Weld. J. 71 (3) (1992) 77–93.                            
[225] N.Saluja, J.O. Ilegbusi, J. Szekely, On the calculation of the electromagneticforce fifield in the circular stirring of metallic melts, J. Appl. Phys. 68(11) (1990) 5845–5850.[226] H. Solhed, L. Jonsson, An investigation of slagflflotation and entrapment in a continuous-casting tundish using flfluid-flflowsimulations, sampling and physical metallurgy, Scand. J. Metall. 32 (2003)15–32.                            
[227] E. Rodriguez, A.H. Castillejos, F.A. Acosta, Experimental andnumerical investigation of flfluid flflow and mixing in pachuca tanks, Metall.Mater. Trans. B 38B (2007) 641–656.[228] D.B. Spalding, CAD to SFT, withAeronautical applications, Plenary Lecture, in: Proceedings of the 38th IsraelAnnual Conference on Aerospace Sciences, 1998.                            
[229] P.G. Tucker, Assessment ofgeometric multilevel convergence robustness and a wall distance method forflflows with multiple internal boundaries, Appl. Math. Modell. 22 (1998)293–311.                            
[230] P.G. Tucker, Differential equation-based wall distancecomputation for DES and RANS, J. Comput. Phys. 190 (1) (2003) 229–248.                            
[231] D.Aganofer, L. Gan-Li, D.B. Spalding, The LVEL turbulence model for conjugateheat transfer at low Reynolds numbers, in: Proceedings of the EEP Applicationof CAE/CAD to Electronic Systems, ASME International Mechanical EngineeringCongress and Exposition, Atlanta, GA, 1996.                            
[232] D.B. Spalding, A newanalytical expression for the drag of a flflat plate valid for both theturbulent and laminar regimes, Int. J. Heat Mass Transfer 5 (12) (1962)1133–1138.                            
[233] P. Rodgers, J. Lohan, V. Eveloy, C.M. Fager, J. Rantala,Validating numerical predictions of component thermal interaction on electronicprinted circuit boards in forced convection air flflows by experimentalanalysis, Adv. Electron. Packag. 1 (1999) 999–1008.                            
[234] K.K. Dhinsa, C.J.Bailey, K.A. Pericleous, Turbulence Modelling and Its Impact on CFD Predictionsfor Cooling of Electronic Components, in: Proceedings of the Ninth IntersocietyConference Thermal and Thermomechanical Phenomena in Electronic Systems,2004.                            
[235] J. Choi et al., A CFD-based tool for studying temperatures inrack-mounted servers, IEEE Trans. Comput. 57 (8) (2008) 1129–1142.                            
[236] N.B.K.Rasmussen, The composite radiosity and gap (CRG) model of thermal radiation,in: Proceedings of the 6th European Conference on Industrial Furnaces andBoilers (INFUB-6) 2002 Conference, Estoril, Lisbon, Portugal, 2002. (Alsopublished as Danish Gas Technology Centre Report No. CO201, Hørsholm, Demark,2003.)                            
[237] Y. Yang, R.A. de Jong, M. Reuter, Use of CFD to Predict thePerformance of a heat treatment furnace, in: Proceedings of the 4thInternational Conference on CFD in the Oil and Gas, metallurgical and processIndustries, Trondheim, Norway, 2005, pp. 1–9.                            
[238] D.B. Spalding, Idealisationsof Radiation, HTS/80/1, Mech. Eng., Dept., Imperial College, University ofLondon, 1980.                            
[239] Y. Yang, R.A. de Jong, M. Reuter, CFD prediction for theperformance of a heat treatment furnace, Prog. Comput. Fluid Dyn. Int. J. 7(2–4) (2007) 209–218.                            
[240] W.N. Hien, A.D. Istiadji, Effects of externalshading devices on daylighting and natural ventilation, in: Proceedings of the8th International IBPSA Conference, Eindhoven, The Netherlands, 2003, pp.475–482.                            
[241] M.R. Malin, D.B. Spalding, Extensions to the PHOENICS ParabolicSolver for under-Expanded Jets, Report C/4366/1 & 2, CHAM Ltd., London,1998.                            
[242] F.C. Lockwood, T. Abbas, N.H. Kandamby, V. Sakthitharan, CFDexperience on industrial combustors, Prog. Comput. Fluid Dyn. Int. J. 1(1/2/3) (2001) 1–13.                            

                            

本文作者：

V. Artemov , S.B. Beale , G. de Vahl Davis, M.P. Escudier, N. Fueyo, B.E. Launder, E. Leonardi＊ , M.R. Malin, W.J. Minkowycz, S.V. Patankar, A. Pollard, W. Rodi, A. Runchal, S.P. Vanka（注：Eddie Leonardi 于2008年12月14日去世。）  

来源：多相流在线

ACTMechanicalMaxwellSystemMarc碰撞非线性多相流燃烧化学通用航空航天船舶核能理论材料
著作权归作者所有，欢迎分享，未经许可，不得转载
首次发布时间：2023-09-29
最近编辑：1年前
积鼎科技
积鼎科技
联系我们13162025768
获赞 108粉丝 109文章 300课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
哈尔滨工业大学王苹：焊接残余应力场模拟，调控及影响（回放）基于ANSYS Twin Builder连杆结构数字孪生体建模关键技术及应用连接器产品结构/电流/热多场耦合仿真解决方案(回仿)第三届中国仿真技术应用大会主论坛——线上会议特邀主旨报告
相关推荐
CoSim 系列软件 材料非线性问题matlab有限元编程求解 ABAQUS-复合材料工程应用案例四-芳纶纤维增强树脂基复合材料切削损伤失效模拟 抗凹插件开发(压强法)ANSA二次开发工程实例应用12讲 ABAQUS-复合材料工程应用案例二-玻璃纤维复合材料泡沫夹层板弯曲变形损伤失效模拟
最新文章
信测标准华东再落一子，宁波检验检测中心正式落成投产（带一3米法电波暗室）
突发！某锂电厂区发生火情！
全固态电池量产线投产，或将点燃新能源领域，推动相关板块上涨！
2024前三季度中国锂电十大排行榜！
基于领域适应性的数字孪生辅助人工智能框架，用于离心泵轴承缺陷诊断（下）
热门文章
盘点·近十年来国外各公司推出的碳纤维产品
ABAQUS中Cohesive粘聚力模型的2种定义方式（附案例操作步骤）
Abaqus分析常见问题及解决方法（2）：零主元和过约束
几种常见的热仿真软件
仿真工作者必须知道的15款开源软件！
其他人都在看
电磁场仿真 | ChatGPT请回答，我想和你聊聊
Abaqus分析常见问题及解决方法（3）：负特征值（Negative Eigenvalue）
CFX进行共轭传热分析实例
STAR CCM 案例｜电池包散热
春招进行时：仿真人才库电磁仿真工程师内推专场（含社招）
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈