Research Papers

Benjamin Franklin and Computational Fluid Dynamics: The Population Approach to Turbulence

[+] Author and Article Information
D. Brian Spalding

Concentration Heat & Momentum (CHAM) Limited,
Bakery House, 40 High Street,
Wimbledon Village,
London, SW19 5AU, England
e-mail: brianspalding@cham.co.uk

Manuscript received November 4, 2010; final manuscript received September 6, 2012; published online December 6, 2012. Assoc. Editor: Gerard F. Jones.

J. Heat Transfer 135(1), 011005 (Dec 06, 2012) (12 pages) Paper No: HT-10-1516; doi: 10.1115/1.4007652 History: Received November 04, 2010; Revised September 06, 2012

Turbulent fluids can be usefully considered as populations, the members of which are distinguished by their positions in non-geometric spaces of which the dimensions may be arbitrarily chosen. These dimensions can be discretised in the same manner as is customary for distance and time. In a distance-time “cell,” variations in mass per unit volume are computed from balances of mass, energy, etc., in which the influences of neighbouring cells are represented in terms of convection and diffusion. Just so can distributions in population space be deduced from balances of mass migration effected by other physical or chemical processes. One such process is movement in population space, as when exothermic chemical reaction moves fluid from a lower to a higher temperature interval; it is akin to convection in geometric space. Another is coalescence, by way of which population members having differing attributes merge and are replaced by elements of intermediate attribute. A third is differential convection, whereby the differing body forces experienced by population elements of differing densities, by reason of the pressure gradients which they share, cause “sifting” and “filtering” to change the population composition. Just as turbulent diffusion is expressed by man-created hypotheses, so hypotheses must be invented for movement in population space, coalescence, and differential convection. The above ideas are explained in the present paper and represented mainly by reference to the age distribution in human populations, to swirling flows and to the distributions of mixture ratio and temperature within turbulent combusting gases.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Daskopoulos, P., and Harris, C. K., 1996, “Three-Dimensional CFD Simulations of Turbulent Flow in Baffled Stirred Tanks: An Assessment of the Current Position,” ICHEME Series No. 140.
Spalding, D. B., 1996, “Progress Report on the Development of a Multi-Fluid Model of Turbulence and Its Application to the Paddle-Stirred Mixer/Reactor,” Invited Lecture at 3rd Colloquium on Process Simulation, Espoo, Finland, June 12–14.
Chesters, J. H., Howes, R. S., Halliday, I. M. D., and Philip, A. R., 1949, “Flow Patterns in Open-Hearth Furnaces,” J. Iron Steel Inst., 163, p. 385.
Longwell, J. P., and Weiss, M. A., 1955, “High-Temperature Reaction Rates in Hydrocarbon Combustion,” Ind. Eng. Chem.47, pp. 1634–1643. [CrossRef]
Patankar, S. V., and Spalding, D. B., 1974, “Simultaneous Prediction of Flow Patterns and Radiation for Three-Dimensional Flames,” Heat Transfer in Flames, N.Afgan and J. M.Beer, eds., John Wiley, New York.
Spalding, D. B., 1971, “Mixing and Chemical Reaction in Confined Turbulent Flames,” 13th International Symposium on Combustion, The Combustion Institute, pp. 649–657.
Spalding, D. B., 1971, “Concentration Fluctuations in a Round Turbulent Free Jet,” J. Chem. Eng. Sci., 26, pp. 95–107. [CrossRef]
Spalding, D. B., 1987, “A Turbulence Model for Buoyant and Combusting Flows,” Int. J. Numer. Methods Eng., 24, pp. 1–23. [CrossRef]
Spalding, D. B., 1979, “The Influences of Laminar Transport and Chemical Kinetics on the Time-Mean Rate in a Turbulent Flame,” Proceedings 17th International Combustion Symposium, Combustion Institute, pp. 431–439.
Peters, N., 1986, “Laminar Flamelet Concepts in Turbulent Combustion,” Proceedings 21st International Combustion Symposium, The Combustion Institute, Pittsburgh, pp. 1231–1250.
Dopazo, C., and O'Brien, E. E., 1974, “An Approach to the Auto Ignition of a Turbulent Mixture,” Acta Astronaut., 1, pp. 1239–1266. [CrossRef]
Pope, S. B., 1982, “An Improved Turbulent Mixing Model,” Combust. Sci. Technol., 28, pp. 131–145. [CrossRef]
Spalding, D. B., 1995, “Models of Turbulent Combustion,” Proceedings of 2nd Colloquium on Process Simulation, Helsinki University of Technology, Espoo, Finland, pp. 1–15.
Poinsot, T., Staffelbach, G., Boileau, M., Wolf, P., 2010, “Using LES to Study Reacting Flows and Instabilities in Annular Combustion Chambers,” Proceedings of 8th International Symposium on Engineering Turbulence Modelling and Measurements, Vol. 2. pp. 301–308.
Freeman, D. J., and Spalding, D. B., 1995, “The Multi-Fluid Turbulent Combustion Model and Its Application to the Simulation of Gas Explosions,” PHOENICS J., 8, pp. 255–277.
Lockwood, F. C., and Naguib, A. S., 1975, “The Prediction of the Fluctuations in the Properties of Free, Round-Jet, Turbulent, Diffusion Flames,” Combust. Flame, 24, pp. 109–124. [CrossRef]
Magnussen, B. F., and Hjertager, B. H., 1976, “On Mathematical Modelling of Turbulent Combustion With Special Emphasis on Soot Formation and Combustion,” 16th International Symposium on Combustion, The Combustion Institute, pp. 719–729.
Patankar, S. V., and Spalding, D. B., 1972, “A Calculation Procedure for Heat, Mass and Momentum Transfer in Three-Dimensional Boundary Layers,” Int. J. Heat Mass Transfer, 15, pp. 1787–1806. [CrossRef]
Spalding, D. B., 1998, “The Simulation of Smoke Generation in a 3-D Combustor, by Means of the Multi-Fluid Model of Turbulent Chemical Reaction,” Paper presented at the “Leading-Edge-Technologies Seminar” on “Turbulent Combustion of Gases and Liquids,” Energy-Transfer and Thermofluid-Mechanics Groups, Institution of Mechanical Engineers, Lincoln, England, Dec. 15–16.


Grahic Jump Location
Fig. 1

The TriMix diagram fraction

Grahic Jump Location
Fig. 6

Contour lines illustrating TriMix reaction rates

Grahic Jump Location
Fig. 7

Visual representation of five turbulent-combustion models

Grahic Jump Location
Fig. 9

Interrelation of geometric and population grids

Grahic Jump Location
Fig. 10

Numerical diffusion

Grahic Jump Location
Fig. 11

Determining Pi (π) by discretisation

Grahic Jump Location
Fig. 12

Contour diagrams of product concentration (left) and PDFs (right) for an acid + base reaction




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In