S. V. Patankar 和 D. B. Spalding 提出了分布阻力的概念，来表征壳程固体表面对流体流动的影响，即把壳程内的换热管束、支撑板以及折流挡板等固体结构看成是多孔介质的骨架，用体积多孔度来表示流体流动所占有的空间占到整个换热器名义空间的百分比，形成了简化的换热器流体流动和传热的数学模型。这种方法既保留了用微分方程来描述壳程流体流动与传热的特点，又简化了换热器的复杂的几何结构，从而实现了管壳式换热器壳程流场的数值模拟。尽管 S. V. Patankar 和 D. B. Spalding 所研究的换热器较简单，网格划分较粗，计算精度也较低，未准确地得到换热管束间的流体流动情况。但是其模拟过程阐述了流动的数值模拟技术，为常用换热器的三维流场数值模拟提供了一定基础。

正是基于 S. V. Patankar 和 D. B. Spalding 的分布阻力概念，后来的研究者们提出了更加精细的数值模拟方法，以更准确地模拟管壳式换热器壳程流场。例如，基于计算流体力学（CFD）的方法可以通过网格划分、数值方法和求解算法等手段，对流场进行模拟和计算，得到更加准确的结果。这些方法为换热器设计和优化提供了重要的工具和手段。同时，这些研究也有助于深入理解流体力学的基本原理和数值模拟的技术方法。

相关问题

Orkiszewski J.气液两相流方法用matlab编程实现

实现气液两相流方法的编程需要掌握以下几个方面：

1. 了解气液两相流模型及其基本方程式：气相方程、液相方程、质量守恒方程、动量守恒方程、能量守恒方程等。

2. 选择数值方法：通常情况下，气液两相流模型采用有限体积法（FVM）或有限元法（FEM）进行求解。选择合适的数值方法可以保证计算结果的准确性和稳定性。

3. 编写程序：选择合适的编程语言，如MATLAB，并编写程序实现气液两相流模型的求解。具体实现过程包括将方程离散化，建立网格，计算通量、速度和压力等，最后求解方程组得到结果。

参考文献：

[1] Orkiszewski J. Two-phase flow calculations using an Eulerian-Lagrangian method[J]. International Journal of Multiphase Flow, 1985, 11(3): 345-360.

[2] Patankar S V. Numerical heat transfer and fluid flow[M]. CRC press, 1980.

[3] Versteeg H K, Malalasekera W. An introduction to computational fluid dynamics: the finite volume method[M]. Pearson Education Limited, 2007.