Efficient calculation of the radiation heat flux surrounding a jet fire

Abstract The major objective of the paper is to address the main weakness of the reverse Monte-Carlo method applied to jet fire simulation; its slow convergence rate. This is done by replacing the use of a pseudo-random number generator for calculating ray orientations with a ray direction specification based on Sobol sequences. Sobol sequences produce sequences of points on the unit hemisphere that are not independent of each other in that new points in the sequence avoid previous points generated. This has the property that a more uniform pattern of rays on the unit hemisphere is produced, giving a ray convergence rate for the incident heat flux that is asymptotically equivalent to O ( N Ray −1 ). The use of Sobol sequences to accelerate convergence of the Monte-Carlo method has been applied in mathematical physics and financial modelling but the results presented here are the first study of quasi-Monte-Carlo methods applied to the incident heat flux integral. The use of Sobol sequences to generate ray directions means that the Central Limit Theorem no longer holds. In its place, it has been demonstrated that it is possible to construct a Gaussian variable from the incident intensity distributions calculated using Sobol sequences. This means it is possible to calculate confidence limits for a prediction of incident heat flux and the confidence limits contract with ray number at a rate of O ( N Ray −1  ln( N Ray )).

[1]  Bernard Lapeyre,et al.  Introduction to Monte-Carlo Methods for Transport and Diffusion Equations , 2003 .

[2]  W. Tao,et al.  Ray effect in ray tracing method for radiative heat transfer , 1997 .

[3]  Michael Fairweather,et al.  Evaluation of flame emission models combined with the discrete transfer method for combustion system simulation , 2005 .

[4]  Weeratunge Malalasekera,et al.  An Adaptive Angular Quadrature for the Discrete Transfer Method Based on Error Estimation , 2003 .

[5]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[6]  Michael F. Modest,et al.  Backward Monte Carlo Simulations in Radiative Heat Transfer , 2003 .

[7]  W. A. Fiveland,et al.  Three-dimensional radiative heat-transfer solutions by the discrete-ordinates method , 1988 .

[8]  P. Cumber,et al.  Modelling lifted hydrogen jet fires using the boundary layer equations , 2009 .

[9]  B. Keister.  Multidimensional quadrature algorithms , 1996 .

[10]  M. Fairweather,et al.  Predictions of radiative transfer from a turbulent reacting jet in a cross-wind , 1992 .

[11]  G. Cox,et al.  On Radiant Heat Transfer from Turbulent Flames , 1977 .

[12]  Stephen B. Pope,et al.  An explanation of the turbulent round-jet/plane-jet anomaly , 1978 .

[13]  T. Gautam,et al.  Lift-Off Heights and Visible Lengths of Vertical Turbulent Jet Diffusion Flames in Still Air , 1984 .

[14]  Modelling free and impacting underexpanded jet fires , 1997 .

[15]  S. Jeng,et al.  The structure of buoyant methane and propane diffusion flames , 1982 .

[16]  Michael Spearpoint,et al.  Modeling Lifted Methane Jet Fires Using the Boundary-Layer Equations , 2006 .

[17]  K. D. Lathrop.  Remedies for Ray Effects , 1971 .

[18]  Michael Fairweather,et al.  Predictions of free jet fires from high pressure, sonic releases , 2003 .

[19]  Pierre L'Ecuyer,et al.  Quasi-Monte Carlo methods in finance , 2004, Proceedings of the 2004 Winter Simulation Conference, 2004..

[20]  P. S. Cumber,et al.  Improvements to the discrete transfer method of calculating radiative heat transfer , 1995 .

[21]  P. S. Cumber.  Ray effect mitigation in jet fire radiation modelling , 2000 .

[22]  I. A. Antonov,et al.  An economic method of computing LPτ-sequences , 1979 .

[23]  H. Rushmeier,et al.  Simultaneous optical measurement of soot volume fraction, temperature, and CO2 in heptane pool fire , 1994 .

[24]  K. D. Lathrop.  RAY EFFECTS IN DISCRETE ORDINATES EQUATIONS. , 1968 .

[25]  Michael Spearpoint,et al.  A computational flame length methodology for propane jet fires , 2006 .

[26]  S. Baillie,et al.  A Phenomenological Model for Predicting the Thermal Loading to a Cylindrical Vessel Impacted by High Pressure Natural Gas Jet Fires , 1998 .

[27]  Tao Pang,et al.  An Introduction to Computational Physics , 1997 .

[28]  P. S. Cumber.  Application of Adaptive Quadrature to Fire Radiation Modeling , 1999 .

[29]  M. Fairweather,et al.  Predictions of soot formation in turbulent, non-premixed propane flames , 1992 .

[30]  P. S. Cumber,et al.  Application of wide band radiation models to non-homogeneous combustion systems , 1998 .

[31]  D. Drysdale.  An Introduction to Fire Dynamics , 2011 .

[32]  F. Lockwood,et al.  A new radiation solution method for incorporation in general combustion prediction procedures , 1981 .

[33]  M. Pinar Mengüç,et al.  Thermal Radiation Heat Transfer , 2020 .

[34]  I. Sobol.  Uniformly distributed sequences with an additional uniform property , 1976 .

[35]  Suhas V. Patankar,et al.  RAY EFFECT AND FALSE SCATTERING IN THE DISCRETE ORDINATES METHOD , 1993 .

[36]  Gerard M. Faeth,et al.  Nonluminous Radiation in Turbulent Buoyant Axisymmetric Flames , 1984 .

[37]  William L. Grosshandler,et al.  Radiative heat transfer in nonhomogeneous gases: A simplified approach , 1980 .

[38]  D. Spalding,et al.  GENMIX: A general computer program for two-dimensional parabolic phenomena , 1977 .