A fully coupled fluid/solid model for open air combustion of horizontally-oriented PMMA samples

Abstract The predictive capability of computational fluid dynamics (CFD) fire models is highly dependent on the accuracy with which the source term due to fuel pyrolysis can be determined. The pyrolysis rate is a key parameter controlling fire behavior, which in turn drives the heat feedback from the flame to the fuel surface. The main objectives of the present study were twofold. First, an in-depth pyrolysis model of a semi-transparent solid fuel (here, clear poly-methyl-methacrylate or PMMA) with in-depth radiation and a moving gas/solid interface was coupled with a CFD code which included turbulence, combustion and radiation for the gas phase. Second, experiments were conducted in order to validate coupled model results. A combined genetic algorithm/pyrolysis model was used with Cone Calorimeter data from a non-flaming pyrolysis experiment to estimate a unique set of kinetic parameters for PMMA pyrolysis. Flaming experiments were conducted on square slabs of PMMA with side dimensions of 10, 20 and 40 cm. From data collected at the center of the slabs, it was found that i) for any sample size, the experimental regression rate remains almost constant with time, with average values of 5.8, 8.6 and 10.9 µm s −1 for the PMMA slabs with side lengths of 10, 20 and 40 cm respectively, and ii) although the radiative and total heat transfers increase significantly with sample size, the radiative contribution to the total heat flux remains almost constant (∼80%). Coupled model results show a fairly good agreement with the literature and with current measurements of the heat fluxes, gas temperature and regressing surface rate at the center of the slabs. Predicted flame heights based on a threshold temperature criterion were found to be close to those deduced from the correlation of Heskestad. However, discrepancies between predicted and measured total pyrolysis rates are observed, which result from the underestimation of the flame heat feedback at the edges of the slab, as confirmed by the comparison between predicted and observed topography of burned samples.

[1]  B. Porterie,et al.  Steady and transient pyrolysis of thick clear PMMA slabs , 2015 .

[2]  Paul E. DesJardin,et al.  An embedded upward flame spread model using 2D direct numerical simulations , 2009 .

[3]  G. D. Raithby,et al.  A Finite-Volume Method for Predicting a Radiant Heat Transfer in Enclosures With Participating Media , 1990 .

[4]  P. Beaulieu,et al.  Flammability characteristics at applied heat flux levels up to 200 kW/m2 , 2008 .

[5]  Hugues Pretrel,et al.  Numerical simulations of fire-induced doorway flows in a small scale enclosure , 2015 .

[6]  S. Patankar,et al.  Finite volume method for radiation heat transfer , 1994 .

[7]  Edwin R. Galea,et al.  The numerical simulation of the noncharring pyrolysis process and fire development within a compartment , 1999 .

[8]  T. Kashiwagi,et al.  Transient Modeling of Thermal Degradation in Non-Charring Solids , 1999 .

[9]  L. Rigollet,et al.  Verification and validation of a CFD model for simulations of large-scale compartment fires , 2011 .

[10]  Francis H. Harlow,et al.  Transport Equations in Turbulence , 1970 .

[11]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[12]  B. Maury.  Characteristics ALE Method for the Unsteady 3D Navier-Stokes Equations with a Free Surface , 1996 .

[13]  B. Hjertager,et al.  On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion , 1977 .

[14]  A. Kulkarni,et al.  Heat Loss to the Interior of a Free Burning Vertical Wall and its Influence on Estimation of Effective Heat of Gasification , 1990 .

[15]  B. Porterie,et al.  Optical and radiative properties of clear PMMA samples exposed to a radiant heat flux , 2014 .

[16]  Jun Yu Li,et al.  A flame spread simulation based on a comprehensive solid pyrolysis model coupled with a detailed empirical flame structure representation , 2015 .

[17]  M. Modest.  Radiative heat transfer , 1993 .

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

[19]  Guillermo Rein,et al.  The application of a genetic algorithm to estimate material properties for fire modeling from bench-scale fire test data , 2006 .

[20]  Vb Novozhilov,et al.  Computational fluid dynamics modeling of compartment fires , 2001 .

[21]  Jose L. Torero,et al.  Average centreline temperatures of a buoyant pool fire obtained by image processing of video recordings , 1995 .

[22]  B. Porterie,et al.  A theoretical and numerical evaluation of the steady-state burning rate of vertically oriented PMMA slabs , 2008 .

[23]  James G. Quintiere,et al.  Fundamentals of Fire Phenomena , 2006 .

[24]  Michel Bellet,et al.  An arbitrary Lagrangian–Eulerian finite element approach to non-steady state turbulent fluid flow with application to mould filling in casting , 2000 .

[25]  Experimental study on radiation attenuation by a water film , 2014 .

[26]  A. Snegirev,et al.  Statistical modeling of thermal radiation transfer in buoyant turbulent diffusion flames , 2004 .

[27]  Jose L. Torero,et al.  Experimental observations on the steady-state burning rate of a vertically oriented PMMA slab , 2008 .

[28]  Nicholas A. Dembsey,et al.  Effect of oxygen on flame heat flux in horizontal and vertical orientations , 2008 .

[29]  Christopher W. Wilson,et al.  Genetic algorithms for optimisation of chemical kinetics reaction mechanisms , 2004 .

[30]  K. McGrattan,et al.  Modeling Solid Sample Buming , 2005 .

[31]  M. Gollner,et al.  Estimation of local mass burning rates for steady laminar boundary layer diffusion flames , 2015 .

[32]  T. F. Smith,et al.  Evaluation of Coefficients for the Weighted Sum of Gray Gases Model , 1982 .

[33]  Jean-Louis Delfau,et al.  Experimental and Numerical Study of the Thermal Degradation of PMMA , 1987 .

[34]  R. Chitty,et al.  A study of the deterministic properties of unbounded fire plumes , 1980 .