Fire Modeling
Fire Modeling
Fire Modeling
EPRI 1023259
Final Report
NUREG-1934
EPRI 1023259
Final Report
November 2012
NOTE
For further information about EPRI, call the EPRI Customer Assistance Center at 800.313.3774 or e-mail
[email protected].
Electric Power Research Institute, EPRI, and TOGETHERSHAPING THE FUTURE OF ELECTRICITY
are registered service marks of the Electric Power Research Institute, Inc.
Copyright 2012 Electric Power Research Institute. All rights reserved.
ABSTRACT
There is a movement to introduce risk-informed and performance-based (RI/PB) analyses into
fire protection engineering practice, both domestically and worldwide. This movement exists in
both the general fire protection and the nuclear power plant (NPP) fire protection communities.
The U.S. Nuclear Regulatory Commission (NRC) has used risk-informed insights as a part of its
regulatory decision making since the 1990s.
In 2001, the National Fire Protection Association (NFPA) issued NFPA 805, Performance-Based
Standard for Fire Protection for Light-Water Reactor Electric Generating Plants, 2001 Edition.
In July 2004, the NRC amended its fire protection requirements in Title 10, Section 50.48 of the
Code of Federal Regulations (10 CFR 50.48) to permit existing reactor licensees to voluntarily
adopt fire protection requirements contained in NFPA 805 as an alternative to the existing
deterministic fire protection requirements. In addition, the NPP fire protection community has
been using RI/PB approaches and insights to support fire hazard analysis in general.
One key element in RI/PB fire protection is the availability of verified and validated (V&V) fire
models that can reliably estimate the effects of fires. The U.S. NRC, together with the Electric
Power Research Institute (EPRI) and the National Institute of Standards and Technology
(NIST), conducted a research project to verify and validate five fire models that have been used
for NPP applications. The results of this effort are documented in a seven-volume NUREG
report, NUREG-1824 (EPRI 1011999), Verification & Validation of Selected Fire Models for
Nuclear Power Plant Applications.
This report describes the implications of the V&V results for fire model users. The features and
limitations of the five fire models documented in NUREG-1824 (EPRI 1011999) are discussed
relative to NPP fire hazard analysis (FHA). Finally, the report provides information on the use of
fire models in support of various commercial NPP fire hazard analysis applications.
iii
CONTENTS
ABSTRACT ................................................................................................................................. iii
CONTENTS .................................................................................................................................. v
FIGURES ..................................................................................................................................... xi
TABLES ...................................................................................................................................... xv
REPORT SUMMARY ............................................................................................................... xvii
PREFACE .................................................................................................................................. xix
CITATIONS................................................................................................................................ xxi
ACKNOWLEDGEMENTS ....................................................................................................... xxiii
LIST OF ACRONYMS .............................................................................................................. xxv
NOMENCLATURE ................................................................................................................. xxvii
1 INTRODUCTION .................................................................................................................... 1-1
1.1 Background ..................................................................................................................... 1-1
1.2 Objective ......................................................................................................................... 1-2
1.3 Scope .............................................................................................................................. 1-2
1.3.1 User Capabilities ..................................................................................................... 1-2
1.3.2 Training Resources ................................................................................................. 1-3
1.4 Fire Modeling Theory ...................................................................................................... 1-4
1.5 Fire Modeling Tools ......................................................................................................... 1-7
1.5.1 Algebraic Models ..................................................................................................... 1-7
1.5.2 Zone Models ............................................................................................................ 1-7
1.5.3 CFD Models............................................................................................................. 1-9
1.5.4 Fire Model Verification and Validation (V&V) ........................................................ 1-11
1.6 Fire Modeling Applications ............................................................................................ 1-11
1.6.1 License Amendments and Exemptions ................................................................. 1-11
1.6.2 Fire Induced Circuit Failures .................................................................................. 1-11
1.6.2.1 Deterministic Application ............................................................................. 1-12
1.6.2.2 Risk-Informed, Performance-Based Application.......................................... 1-12
1.6.3 NFPA 805 Performance-Based Applications ......................................................... 1-12
1.6.3.1 Fire Modeling ............................................................................................... 1-13
1.6.3.2 Fire Risk Evaluations ................................................................................... 1-14
1.6.4 Fire Modeling in Support of Fire PRA .................................................................... 1-14
1.7 Organization of the Guide ............................................................................................. 1-15
2 THE FIRE MODELING PROCESS ......................................................................................... 2-1
2.1 Step 1: Define Fire Modeling Goals and Objectives........................................................ 2-2
2.2 Step 2: Characterize Fire Scenarios ............................................................................... 2-4
2.2.1 General Considerations ........................................................................................... 2-4
2.2.2 Enclosure Details..................................................................................................... 2-5
2.2.3 Fire Location ............................................................................................................ 2-6
2.2.4 Credited Fire Protection........................................................................................... 2-7
2.2.5 Ventilation Conditions .............................................................................................. 2-8
2.2.6 Target Locations ...................................................................................................... 2-8
2.2.7 Secondary Combustibles ......................................................................................... 2-9
2.2.8 Source Fire .............................................................................................................. 2-9
2.3 Step 3: Select Fire Models ............................................................................................ 2-10
2.3.1 Fire Dynamics Tools (FDTs) .................................................................................. 2-14
2.3.2 Fire Induced Vulnerability Evaluation (FIVE-Rev1) ............................................... 2-16
2.3.3 Consolidated Fire Growth and Smoke Transport (CFAST) Model ........................ 2-18
2.3.4 MAGIC ................................................................................................................... 2-19
2.3.5 Fire Dynamics Simulator (FDS) ............................................................................. 2-19
2.3.6 Verification and Validation (V&V) .......................................................................... 2-22
2.3.7 Fire Modeling Parameters Outside the Validation Range ..................................... 2-28
2.3.7.1 Sensitivity Analysis ...................................................................................... 2-29
2.3.7.2 Additional Validation Studies ....................................................................... 2-31
2.4 Step 4: Calculate Fire-Generated Conditions................................................................ 2-32
2.5 Step 5: Conduct Sensitivity and Uncertainty Analyses.................................................. 2-33
2.6 Step 6: Document the Analysis ..................................................................................... 2-34
2.7 Summary ....................................................................................................................... 2-35
3 GUIDANCE ON FIRE MODEL SELECTION AND IMPLEMENTATION................................ 3-1
3.1 Model Implementation of Fire Scenario Elements........................................................... 3-1
3.1.1 Heat Release Rate .................................................................................................. 3-1
3.1.2 Plant Area Configuration ......................................................................................... 3-2
3.1.3 Ventilation Effects .................................................................................................... 3-3
3.1.4 Targets .................................................................................................................... 3-4
3.1.5 Intervening Combustibles ........................................................................................ 3-5
3.2 Guidance on Model Selection and Analysis .................................................................... 3-5
3.2.1 Targets in the Flames or Plume .............................................................................. 3-8
3.2.1.1 General Objective .......................................................................................... 3-8
3.2.1.2 Modeling Strategy.......................................................................................... 3-8
3.2.1.3 Recommended Models.................................................................................. 3-9
3.2.1.4 Detailed Examples....................................................................................... 3-10
3.2.2 Scenario 2: Targets Inside or Outside the Hot Gas Layer ..................................... 3-11
3.2.2.1 General Objective ........................................................................................ 3-11
3.2.2.2 Modeling Strategy........................................................................................ 3-11
3.2.2.3 Recommended Modeling Tools ................................................................... 3-12
3.2.2.4 Detailed Examples....................................................................................... 3-13
3.2.3 Scenario 3: Targets Located in Adjacent Rooms .................................................. 3-14
3.2.3.1 General Objective ........................................................................................ 3-14
3.2.3.2 Modeling Strategy........................................................................................ 3-14
3.2.3.3 Recommended Modeling Tools ................................................................... 3-15
3.2.3.4 Detailed Examples....................................................................................... 3-15
3.2.4 Scenario 4: Targets in Rooms with Complex Geometries ..................................... 3-16
3.2.4.1 General Objective ........................................................................................ 3-16
3.2.4.2 Modeling Strategy........................................................................................ 3-16
3.2.4.3 Recommended Modeling Tools ................................................................... 3-17
3.2.4.4 Detailed Examples....................................................................................... 3-17
3.2.5 Scenario 5: Main Control Room Abandonment ..................................................... 3-18
3.2.5.1 General Objective ........................................................................................ 3-18
3.2.5.2 Modeling Strategy........................................................................................ 3-18
3.2.5.3 Recommended Modeling Tools ................................................................... 3-18
3.2.5.4 Detailed Examples....................................................................................... 3-19
3.2.6 Scenario 6: Smoke Detection and Sprinkler Activation ......................................... 3-20
3.2.6.1 General Objective ........................................................................................ 3-20
3.2.6.2 Modeling Strategy........................................................................................ 3-20
3.2.6.3 Recommended Modeling Tools ................................................................... 3-21
3.2.6.4 Detailed Examples....................................................................................... 3-21
3.2.7 Scenario 7: Fire Impacting Structural Elements .................................................... 3-22
3.2.7.1 General Objective ........................................................................................ 3-22
vi
viii
ix
FIGURES
Figure 1-1.
Figure 1-2.
Figure 1-3.
Figure 2-1.
Figure 3-1.
Figure 3-2.
Figure 3-3.
Figure 3-4.
Figure 3-5.
Figure 3-6.
Figure 3-7.
Figure 3-8.
Figure 3-9.
Figure 4-1.
Figure 4-2.
Figure 4-3.
Figure 4-4.
Figure 4-5.
Figure 4-6.
Figure A-1.
Figure A-2.
Figure A-3.
Figure A-4.
Figure A-5.
Figure A-6.
Figure A-7.
Figure A-8.
Figure A-9.
Figure A-10.
Figure A-11.
Figure A-12.
Figure A-13.
Figure A-14.
Figure B-1.
Figure B-2.
Figure B-3.
xi
Figure B-4.
Figure B-5.
Figure B-6.
Figure B-7.
Figure B-8.
Figure B-9.
Figure B-10.
Figure B-11.
Figure B-12.
Figure B-13.
Figure B-14.
Figure B-15.
Figure B-16.
Figure C-1.
Figure C-2.
Figure C-3.
Figure C-4.
Figure C-5.
Figure C-6.
Figure C-7.
Figure C-8.
Figure C-9.
Figure C-10.
Figure C-11.
Figure C-12.
Figure C-13.
Figure C-14.
Figure C-15.
Figure C-16.
Figure D-1.
Figure D-2.
Figure D-3.
Figure D-4.
Figure D-5.
Figure D-6.
Figure D-7.
Figure D-8.
Figure D-9.
Figure D-10.
Average HGL temperature (red line) from FPA correlation and HRR (blue
line) for SWGR cabinet fire scenario. ................................................................ B-12
Average CFAST/Smokeview rendering of SWGR. ........................................... B-12
CFAST inputs for compartment geometry in SWGR scenario. ......................... B-13
CFAST fire specification inputs for the SWGR scenario. .................................. B-14
CFAST mechanical ventilation inputs for the SWGR scenario. ........................ B-15
CFAST inputs for cabinet and cable targets for the SWGR scenario. .............. B-16
FDS/Smokeview rendering of the SWGR. ........................................................ B-18
FDS/Smokeview rendering of the SWGR fire showing localized ignition of
extinction of secondary cable fires resulting from initial cabinet fire. ................ B-18
Heat release rate inputs to CFAST and FDS for a SWGR cabinet fire
scenario............................................................................................................. B-20
Estimated temperatures for Cable Tray A directly above the fire source for
a SWGR cabinet fire scenario. .......................................................................... B-21
Estimated temperature and heat flux to a cabinet adjacent to the fire
source in a SWGR cabinet fire scenario. .......................................................... B-22
Distribution of HRR for an electrical cabinet fire. .............................................. B-23
Distribution of flame heights for the entire range of cabinet fires. ..................... B-23
Geometry of the pump room. .............................................................................. C-2
MAGIC view of the pump room. .......................................................................... C-9
Heat release rate curve for lubricating oil fire.................................................... C-10
MAGIC fire input screen for pump fire............................................................... C-11
Modeling multi-conductor cables in MAGIC. Source: NUREG-1824
(EPRI 1011999), Volume 6, Figure 3-3. ............................................................ C-12
Representation of the ERFBS protected cable raceway for MAGIC. ................ C-13
MAGIC input screen for ERFBS layers. ............................................................ C-13
MAGIC input screen for ventilation through door. ............................................. C-14
FDS/Smokeview rendering of the pump room scenario at the early stage of
the fire, before the compartment becomes under-ventilated............................. C-15
HRR predicted by algebraic methods, MAGIC, and FDS for the pump room
fire scenario....................................................................................................... C-17
HGL Temperature Predicted by MAGIC and FDS for the pump room fire
scenario............................................................................................................. C-17
Cable surface temperature predicted by MAGIC and FDS for the pump
room fire scenario. ............................................................................................ C-19
Cable surface temperature predicted by MAGIC for changes to insulation
wrapping............................................................................................................ C-20
Temperature predicted by MAGIC for increased door size. .............................. C-21
HRR for base case and HRR sensitivity case. .................................................. C-22
Temperature for base case and HRR sensitivity case. ..................................... C-22
Geometry of the MCC/SWGR in a BWR. ............................................................ D-3
Typical electrical cabinet in the lower part of the SWGR. ................................... D-4
View of the high ceiling space. ............................................................................ D-4
Smokeview rendering of the geometry of two-height ceiling SWGR, as
modeled in CFAST. ............................................................................................. D-9
CFAST inputs for compartment geometry for SWGR. ...................................... D-10
CFAST fire inputs for two-height ceiling SWGR scenario. ................................ D-10
CFAST target inputs for two-height ceiling SWGR scenario. ............................ D-11
CFAST inputs for vents connections in two-height ceiling SWGR scenario. .... D-12
FDS/Smokeview representation of the MCC/SWGR scenario. ......................... D-13
Heat flux and temperature predictions for the adjacent cabinet. ....................... D-15
xii
Figure D-11. Summary of the cable temperature and heat flux predictions for the
MCC/SWGR. ..................................................................................................... D-16
Figure E-1. Geometry of the CSR. ......................................................................................... E-4
Figure E-2. Photograph of the CSR used for analysis. Note that the cables are located
in the trays in the overhead. ................................................................................ E-5
Figure E-3. Geometric detail of the CSR. .............................................................................. E-5
Figure E-4. Schematic diagram of transient trash fire in CSR (not to scale). ........................ E-9
Figure E-5. Plume temperatures at cable trays located above a transient trash fire. .......... E-10
Figure E-6. Average HGL temperature (red line) from FPA correlation and heat release
rate (blue line) for the CSR trash fire scenario. ................................................. E-11
Figure E-7. CFAST/Smokeview rendering of the CSR scenario.......................................... E-12
Figure E-8. CFAST inputs for compartment geometry for the CSR scenario. ..................... E-13
Figure E-9. CFAST inputs for the fire in the CSR scenario.................................................. E-14
Figure E-10. CFAST inputs for cable targets above the fire source in the CSR scenario. .... E-15
Figure E-11. CFAST inputs for natural ventilation for the CSR scenario. Note the
opening of the doorway at 170 s, the time estimated for smoke detectors to
activate. ............................................................................................................. E-16
Figure E-12. CFAST inputs for smoke detectors in the CSR scenario. ................................. E-17
Figure E-13. FDS/Smokeview rendering of the CSR scenario. ............................................. E-17
Figure E-14. HRR and estimated HGL temperature for CSR scenario. ................................. E-19
Figure E-15. Estimated cable conditions for the CSR............................................................ E-22
Figure E-16. Distribution of HRR for a trash fire. ................................................................... E-23
Figure E-17. Distribution of plume temperatures at Trays 3 and 6, respectively. .................. E-24
Figure F-1. Geometry of the turbine building. .........................................................................F-2
Figure F-2. Main turbine lubricating oil tanks in the turbine building.......................................F-3
Figure F-3. Typical steel column in the turbine building. ........................................................F-3
Figure F-4. Schematic diagram of the fire impinging on the ceiling. .....................................F-10
Figure F-5. Detail from Figure F-1 with estimated flame extension beneath ceiling
superimposed.....................................................................................................F-11
Figure F-6. FDS geometry for the turbine building fire scenario. ..........................................F-13
Figure F-7. Temperatures of columns for Curb Location 1. ..................................................F-15
Figure F-8. Temperatures of columns for Curb Location 2. ..................................................F-16
Figure G-1. Geometry of the multi-compartment corridor. ..................................................... G-2
Figure G-2. Geometry details of the multi-compartment corridor. .......................................... G-3
Figure G-3. HRR for the corridor fire scenario. ...................................................................... G-4
Figure G-4. Effective corridor layout for implementation in zone models (not to scale). ........ G-8
Figure G-5. MAGIC rendering of the corridor scenario. ......................................................... G-9
Figure G-6. MAGIC screen capture of the fire in corridor scenario. ..................................... G-10
Figure G-7. MAGIC screen capture for supply vent specification. ....................................... G-11
Figure G-8. MAGIC screen capture of the smoke detector specification. ............................ G-12
Figure G-9. HGL temperature predictions by MAGIC for the corridor scenario. .................. G-14
Figure G-10. HGL temperature for reduced ceiling height by MAGIC. ................................... G-15
Figure G-11. Detector temperature prediction by MAGIC for fire corridor.............................. G-15
Figure G-12. Detector temperature for two radial distances predicted by MAGIC. ................ G-16
Figure H-1. Geometry of the annulus. .................................................................................... H-3
Figure H-2. Geometry details of redundant cable trays located in the annulus. .................... H-4
Figure H-3. HRR for a cable fire in the annulus. .................................................................... H-5
Figure H-4. FDS/Smokeview rendering of the annulus fire scenario. .................................. H-10
Figure H-5. Summary of simulation results for the annulus. ................................................ H-12
Figure H-6. Predicted sprinkler link temperature for the annulus fire scenario. ................... H-12
xiii
TABLES
Table 2-1.
Table 2-2.
Table 2-3.
Table 2-4.
Table 2-5.
Table 3-1.
Table 3-2.
Table 4-1.
Table 4-2.
Table 4-3.
Table 6-1.
Table A-1.
Table A-2.
Table A-3.
Table A-4.
Table B-1.
Table B-2.
Table B-3.
Table B-4.
Table B-5.
Table C-1.
Table C-2.
Table C-3.
Table C-4.
Table C-5.
Table C-6.
Table D-1.
Table D-2.
Table D-3.
Table E-1.
Table E-2.
Table E-3.
Table E-4.
Table E-5.
Table E-6.
Table F-1.
Table F-2.
Table F-3.
xv
Table F-4.
Table G-1.
Table G-2.
Table G-3.
Table G-4.
Table H-1.
Table H-2.
Table H-3.
xvi
REPORT SUMMARY
Background
Beginning in the 1990s, when the U.S. Nuclear Regulatory Commission (NRC) adopted the
policy of using risk-informed methods to make regulatory decisions whenever possible, the
nuclear power industry has been moving from prescriptive rules and practices toward the use of
risk information to supplement decision making. Several initiatives have furthered this transition
within the fire protection field, including risk-informed, performance-based fire protection
programs (FPPs) compliant with Title 10, Section 50.48(c) of the Code of Federal Regulations
(10 CFR 50.48(c)) and FPP change evaluation under the existing Title 10 Section 50.48 and
Regulatory Guide 1.189, Fire Protection for Nuclear Power Plants, October 2009. RI/PB fire
protection often relies on fire modeling to estimate the effects of fires.
Objectives
To provide guidance on the proper application of fire models to NPP fire scenarios
To fulfill the need as a teaching tool and support the NRC/EPRI Fire PRA training
Approach
There were five different fire models evaluated in NUREG-1824 (EPRI 1011999), Verification
and Validation of Selected Fire Models for Nuclear Power Plant Applications: (1) the NRCs Fire
Dynamics Tools (FDTs), (2) the Electric Power Research Institutes (EPRI) Fire-Induced
Vulnerability Evaluation Revision 1 (FIVE-Rev1), (3) the National Institute of Standards and
Technologys (NIST) Consolidated Model of Fire Growth and Smoke Transport (CFAST),
(4) Electricit de Frances (EdF) MAGIC, and (5) NISTs Fire Dynamics Simulator (FDS). To
support consistency, the same fire models were used in the development of this report. The
project team developed guidance on the selection and application of each model and treatment
of uncertainty and/or sensitivity as part of the fire modeling analysis. Based on this guidance,
the project team selected appropriate models and conducted fire modeling analyses of eight
different fire scenarios of interest in nuclear power plants (NPPs). The results of each analysis
were to be documented in a format appropriate for fire model calculation submittals.
Results
This report presents a step-by-step process for using fire modeling in NPP applications. The
recommended methodology consists of a six step process: (1) define fire modeling goals, (2)
characterize fire scenarios, (3) select fire models, (4) calculate fire-generated conditions, (5)
conduct sensitivity and uncertainty analyses, and (6) document the results.
This report is designed to assist fire model users in applying this technology in the NPP
environment. There are a number of unique construction and fire hazard attributes associated
with NPPs. It was the authors goal to explore and demonstrate the use of different models for
this application. In addition, a fifth module, Advanced Fire Modeling, has been added to the
annual NRC/EPRI Fire PRA training workshop in 2011. This report expands upon the
information provided in NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA
Methodology for Nuclear Power Facilities, and will serve as the training material for the
Advanced Fire Modeling module.
xvii
The report contains a CD to assist the reader in reproducing the examples discussed in the
appendices. The CD contains the installation files for CFAST and FDS. It also contains the
latest versions (November 2012) of the NUREG-1805 spreadsheets. The CD also contains the
input files for each model used to analyze the eight fire scenarios.
EPRI Perspective
The use of fire models requires a good understanding of their limitations and predictive
capabilities, and also presents challenges that should be addressed if the fire protection
community is to realize the full benefit of fire modeling and performance-based fire protection.
EPRI, in partnership with NRC under a Memorandum of Understanding (MOU), will continue to
provide training to the fire protection community, using this document to promote fire modeling
and gain feedback on how the results of this work may affect known applications of fire
modeling. In the long term, model improvement and additional experiments should be
considered.
This report supersedes EPRI 10002981, Fire Modeling Guide for Nuclear Power Plant
Applications, August 2002, as guidance for fire modeling practitioners in NPPs. The report has
benefited from the insights gained since 2002 on the predictive capability of selected fire models
in improving confidence in the use of fire modeling in NPP applications.
Note: Due to EPRI limitations on distribution of software, the FDTs spreadsheets, CFAST, and
FDS are not included on the EPRI version of the CD. The FDTs can be obtained from the NRC
web site (www.nrc.gov), and CFAST and FDS can be obtained from the NIST web site
(www.nist.gov).
Keywords
Performance-Based
Fire Hazard Analysis (FHA)
Fire Protection
Probabilistic Risk Assessment (PRA)
Fire
Verification and Validation (V&V)
Risk-Informed Regulation
Fire Safety
Nuclear Power Plant (NPP)
Fire Modeling
xviii
PREFACE
This report describes research sponsored jointly by EPRI and U.S. Nuclear Regulatory
Commission, Office of Nuclear Regulatory Research (RES). The main purpose of this report is
to provide guidance on the application of fire models to nuclear power plant (NPP) fire scenarios
and to serve as a teaching tool and support the Advanced Fire Modeling module of the NRCRES/EPRI Fire PRA Course. The fire modeling analyses presented in this report represent the
combined efforts of individuals from RES and EPRI. Both organizations provided specialists in
the use of fire models/fire hazard analysis tools to support this work. These results are intended
to provide technical analysis of the predictive capabilities of five fire modeling calculation tools.
This report is the fifth in a series designed to assist those responsible for performing fire
modeling in NPP applications.
In August 2002, EPRI published EPRI 1002981, Fire Modeling Guide for Nuclear Power Plant
Applications. This report offered step-by-step guidance that analysts could follow when using
fire modeling to support nuclear power plant fire protection applications. It also included FIVERev1, an Excel-based library of fire models previously documented by EPRI, and additional
models from fire protection literature.
In December 2004, the NRC published NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative
Fire Hazard Analysis Methods for the U.S. Nuclear Regulatory Commission Fire Protection
Inspection Program. This report provided an introduction to the principles of fire dynamics, and
included an Excel-based library of fire models comparable to EPRI FIVE-Rev1.
In a follow-up effort as a part of the NRC/RES-EPRI Memorandum of Understanding (MOU),
NRC/RES and EPRI jointly conducted a verification and validation of selected fire models for
use in nuclear power plant fire modeling to gain insight into the predictive capabilities of these
models. The results of this work were published in NUREG-1824 (EPRI 1011999), Verification
and Validation of Selected Fire Models for Nuclear Power Plant Applications, May 2007. Using,
in part, the findings of this work, the NRC conducted a Phenomena Identification and Ranking
Table (PIRT) study to evaluate the current state of knowledge for fire modeling for NPP
applications. The results of this work were published in NUREG/CR-6978, A Phenomena
Identification and Ranking Table (PIRT) Exercise for Nuclear Power Plant Fire Modeling
Applications, November 2008.
This document does not constitute regulatory requirements. RES participation in this
study does not constitute or imply regulatory approval of applications based upon this
methodology.
xix
CITATIONS
This report was prepared by:
U.S. Nuclear Regulatory Commission,
Office of Nuclear Regulatory Research (RES)
Washington, D.C. 20555-0001
Principal Investigator:
D. Stroup
University of Maryland
Principal Investigator:
J. Milke
This report describes research sponsored jointly by the U.S. Nuclear Regulatory Commissions
(NRC) Office of Nuclear Regulatory Research (RES) and the Electric Power Research Institute
(EPRI).
The report is a corporate document that should be cited in the literature in the following manner:
Nuclear Power Plant Fire Modeling Analysis Guidelines (NPP FIRE MAG), U.S. Nuclear
Regulatory Commission, Office of Nuclear Regulatory Research (RES), Washington, D.C.,
2012, and Electric Power Research Institute (EPRI), Palo Alto, CA, NUREG-1934 and EPRI
1023259.
xxi
ACKNOWLEDGEMENTS
This report is the result of a multi-year effort by the authors, during which it has been
significantly rewritten in response to comments from various reviewers. The authors wish to
thank those reviewers for taking the time to provide their thoughts on this document. The
authors express appreciation to the members of the peer review panel who provided comments
on the original draft of this document in 2009: Professor Jose L. Torero and his students at BRE
Centre for Fire Safety Engineering, The University of Edinburgh, Scotland; Professor Frederick
Mowrer, formerly of the Department of Fire Protection Engineering, University of Maryland, and
currently with the Department of Fire Protection Engineering, California Polytechnic State
University; Mr. Patrick Finney, NRC Resident Inspector at Susquehanna Nuclear Plant;
Mr. Naeem Iqbal, Fire Protection Engineer, Office of Nuclear Reactor Regulation (NRR), NRC;
and Mr. Thomas Gorman, Project Manager, Pennsylvania Power and Light, Susquehanna
Nuclear Plant.
Drafts of this report were Noticed twice in the Federal Register for public comment: first on
December 29, 2009 (74 FR 68873), and second on August 2, 2011 (76 FR 46331). The authors
thank those members of the public who provided comments during the two public comment
periods, specifically Robert M. Brady, Schirmer Engineering Corporation; Patricia L. Campbell,
GE Hitachi Nuclear Energy; Jason E. Floyd, Hughes Associates, Inc.; Pablo Guardado, Entergy;
Michael D. Jesse, Excelon Nuclear; Daeil Kang, Korea Atomic Energy Research Institute;
Nancy McNabb, National Fire Protection Association; Daniel Orr, NRC, Region 1; Mark
Schairer, Engineering Planning and Management (EPM); and Robert Webster, AREVA, all of
whom provided insightful comments on the two drafts of the document.
This report was used as a textbook for pilot offerings of an advanced fire modeling course held
as part of the annual NRC/EPRI Fire PRA training in 2011 and 2012. The authors would like to
thank the students in those two classes for their many constructive comments that significantly
improved the final report.
Meetings were held before the Reliability and PRA Subcommittee of the Advisory Committee on
Reactor Safeguards (ACRS) and before the full ACRS committee on March 21, 2012 and
July 11, 2012, respectively. We would like to thank the members of the ACRS for their time and
efforts in supporting the meetings and reviewing and providing comments and suggestions for
improvements to the final report.
Finally, the authors express their thanks to Laurent Gay and Eric Wizenne, Electricit de France
(EdF), for reviewing the MAGIC calculations. The authors would also like to acknowledge Mr.
Bryan Klein of Thunderhead Engineering and Drs. Nathan Siu and Raymond Gallucci of the
NRC for their valuable contributions to this report. The authors also express appreciation to Ms.
Aixa Belen, Mr. Ken Canavan, Mr. Stuart Lewis, Mr. Nicholas Melly, Ms. Carolyn Siu, and Mr.
Robert Vettori for their reviews of and comments on various drafts of this document. Publication
of this document would not have been possible without assistance from the NRCs Office of
Administration. Specifically, the authors thank Tojuana Fortune-Grasty and Guy Beltz for their
efforts in publishing this report.
xxiii
LIST OF ACRONYMS
ACH
Air Changes per Hour
ACRS
Advisory Committee on Reactor Safeguards
AGA
American Gas Association
AHJ
Authority Having Jurisdiction
ANS
American Nuclear Society
ASET
Advanced Science and Engineering Technologies
ASME
American Society of Mechanical Engineers
ASTM
American Society for Testing and Materials
ATF
Bureau of Alcohol, Tobacco, and Firearms
BE
Benchmark Exercise
BFRL
Building and Fire Research Laboratory
BRE
Building Research Establishment
BWR
Boiling Water Reactor
CAROLFIRE Cable Response to Live Fire
CDF
Core Damage Frequency
CFAST
Consolidated Fire Growth and Smoke Transport Model
CFD
Computational Fluid Dynamics
CFR
Code of Federal Regulations
CHRISTIFIRE Cable Heat Release, Ignition, and Spread in Tray Installations during Fire
COR
Code of Record
CSR
Cable Spreading Room
DBD
Design Basis Document
ECCS
Emergency Core Cooling Systems
EdF
Electricit de France
EPRI
Electric Power Research Institute
ERFBS
Electrical Raceway Fire Barrier System
FDS
Fire Dynamics Simulator
FDTs
Fire Dynamics Tools (NUREG-1805)
FHA
Fire Hazard Analysis
FIVE-Rev1
Fire-Induced Vulnerability Evaluation, Revision 1
FFT
Fast Fourier Transform
FLASH-CAT Flame Spread over Horizontal Cable Trays
FM/SNL
Factory Mutual & Sandia National Laboratories
FPA
Foote, Pagni, and Alvares
FPRA
Fire Probabilistic Risk Assessment
FRA
Fire Risk Analysis
GRS
Gesellschaft fr Anlagen- und Reaktorsicherheit (Germany)
HGL
Hot Gas Layer
HRR
Heat Release Rate
HRRPUA
Heat Release Rate Per Unit Area
HVAC
Heating, Ventilation, and Air Conditioning
IAFSS
International Association of Fire Safety Science
iBMB
Institut fr Baustoffe, Massivbau und Brandschutz
ICFMP
International Collaborative Fire Model Project
ID
Identification
IEEE
Institute of Electrical and Electronics Engineers
xxv
IPEEE
LERF
LES
LFS
LLNL
LOL
MCC
MCR
MEFS
MOVs
MQH
MOU
NBS
NEI
NFPA
NIST
NPP
NRC
NRR
PE
PMMA
PRA
PVC
PWR
RCP
RES
RG
RI/PB
RIS
RTE
RTI
SBDG
SDP
SFPE
SNL
SWGR
THIEF
TP
TS
UL
V&V
WTC
XPE
ZOI
xxvi
NOMENCLATURE
area (m2)
floor area (m2)
opening area (m2)
surface area of enclosure boundary (m2)
specific heat, solid material (kJ/kg/K)
specific heat, gas, constant pressure (kJ/kg/K)
fire diameter (m) or optical density (m-1)
characteristic diameter
experimental measurement
acceleration of gravity (m/s2)
heat transfer coefficient (kW/m2/K)
ceiling height (m)
height of fire base above floor (m)
opening height (m)
thermal conductivity (kW/m/K)
constant
volumetric entrainment coefficient
light extinction coefficient (m-1)
mass-specific extinction coefficient (m2/kg)
compartment length (m)
flame height (m)
total mass per unit length of a single cable (kg/m)
combustible mass per unit area of tray (kg/m2)
mass loss or flow rate (kg/s)
mass loss rate per unit area (kg/s/m2)
model prediction
number of cables per tray
probability, wall perimeter (m)
heat flux (kW/m2)
critical heat flux (kW/m2)
heat release rate (kW)
fire Froude number
radial distance (m)
ceiling jet distance (m)
specific area (m2/g)
burnout time (s)
time (s)
temperature (C)
volume (m3)
volume flow rate (m3/s)
compartment width (m), width of tray (m)
cell dimension (m)
critical value
product yield (kg/kg)
mass fraction (kg/kg)
xxvii
Subscripts:
,
xxviii
INTRODUCTION
1.1 Background
In 2001, the National Fire Protection Association (NFPA) issued the first edition of NFPA 805,
Performance-Based Standard for Fire Protection for Light-Water Reactor Electric Generating
Plants, 2001 Edition1. Effective July 16, 2004, the U.S. Nuclear Regulatory Commission (NRC)
amended its fire protection requirements in Title 10, Section 50.48(c) of the Code of Federal
Regulations (10 CFR 50.48(c)) to permit existing reactor licensees to voluntarily adopt fire
protection requirements contained in NFPA 805 following a performance-based approach as an
alternative to the existing deterministic2 fire protection requirements. One important element in
a performance-based approach is the estimation of fire hazard using mathematical fire models.
Fire modeling is often used in constructing Fire PRAs to determine the effects of fire hazard so
that the associated risk can be quantified.
As part of its fire modeling requirements, NFPA 805 states that fire models shall be verified and
validated (section 2.4.1.2.3) and that only fire models that are acceptable to the authority
having jurisdiction (AHJ) shall be used in fire modeling calculations (section 2.4.1.2.1). This is
an important requirement because the verification and validation (V&V) of fire models is
intended to ensure the correctness, suitability, and overall quality of the method. Specifically,
verification is the process used to determine whether a model correctly represents the
developers conceptual description (i.e., whether it was built correctly), while validation is used
to determine whether a model is a suitable representation of the real world and is capable of
reproducing phenomena of interest (i.e., whether the correct model was built).
In 2007, the NRCs Office of Nuclear Regulatory Research (RES) and the Electric Power
Research Institute (EPRI) completed a collaborative project for the V&V of five select fire
modeling tools. The results of this study, which was performed under the NRC/RES-EPRI
Memorandum of Understanding (MOU), are documented in NUREG-1824 (EPRI 1011999),
Verification and Validation of Selected Fire Models for Nuclear Power Plant Applications. The
National Institute of Standards and Technology (NIST) was also an important partner in
developing this publication, providing extensive fire modeling and experimentation expertise.
The V&V effort is intended to support the use of fire modeling for various NPP fire hazard
analysis applications.
This report builds on the V&V research described earlier by incorporating the results into a set
of guidelines and recommendations for conducting fire modeling studies in support of
1
All references in this report to NFPA 805 are specific to the 2001 edition of the standard, which is the code of record (COR)
required by 10 CFR 50.48(c).
2
In nuclear fire protection, the term deterministic is typically used to refer to prescriptive requirements while deterministic is often
used as an adjective with fire model in the general fire protection field. Within this report, the usage should be clear from the
context in which the word is used.
1-1
INTRODUCTION
commercial nuclear industry applications. When the NRCs Advisory Committee on Reactor
Safeguards (ACRS) issued a letter to Luis Reyes, Executive Director for Operations,
recommending publication of NUREG-1824 (EPRI 1011999), they identified two major items to
be included in the users guide (Wallis, 2006). Specifically, the ACRS recommended that the
users guide include:
The ACRS indicated that quantitative estimates of the intrinsic model uncertainty would be a
valuable input in risk-informed as well as non-risk-informed applications. Chapters 2 and 3
address the first ACRS recommendation. Chapter 4 specifically addresses the second ACRS
recommendation, that is, the development of V&V results into quantitative estimates of model
uncertainty. Finally, the appendices contain examples that illustrate the entire process for
several nuclear power plant (NPP) scenarios.
1.2 Objective
The objective of this guide is to describe the process of conducting a fire modeling analysis,
principally for commercial NPP applications. The process described in this guide addresses
most of the technical elements relevant to fire modeling analysis, such as the selection and
definition of fire scenarios and the determination and implementation of input values, sensitivity
analysis, uncertainty quantification, and documentation. In addition, requirements associated
with fire modeling analyses and analytical fire modeling tools are addressed through generic
guidance, recommended best practices, and example applications.
1.3 Scope
1.3.1 User Capabilities
This guide should be used as a complement to, not a substitute for, users manuals for specific
fire modeling tools, fire dynamics textbooks, technical references, education, and training. This
guide only compiles information and organizes it procedurally for NPP applications. Analysts
are encouraged to review the references identified throughout the guide for in-depth coverage of
the advantages and the range of applicability of specific models. Once a fire scenario has been
selected, this guide will help the fire model user define the necessary modeling parameters,
select an appropriate model, and properly interpret the fire modeling results. Since all models
are merely approximations of reality, this guide also provides useful insights for translating real
configurations into modeling scenarios. Due to the technical nature of this guide, users with the
following areas of expertise will benefit the most from it:
INTRODUCTION
This guide focuses on the capabilities of the models selected for V&V. However, some generic
guidance is also provided, and most of the discussion is applicable to any fire model of the
respective type (algebraic model, zone model, or computational fluid dynamics (CFD) model).
Five specific models are discussed in this guide:
(1) The NRCs Fire Dynamics Tools (FDTS), NUREG-1805 and Supplements
(2) EPRIs Fire-Induced Vulnerability Evaluation, Revision 1 (FIVE-Rev1)
(3) NISTs Consolidated Model of Fire Growth and Smoke Transport (CFAST) Version (6)
(4) Electricit de Frances (EdF) MAGIC code Version (4.1.1)
(5) NISTs Fire Dynamics Simulator (FDS) Version (5)
Finally, the user of this document would benefit from a familiarity with the use of fire modeling in
a fire protection performance-based approach and NUREG/CR-6850 (EPRI 1011989),
EPRI/NRC-RES Fire PRA Methodology for Nuclear Power Facilities.
1-3
INTRODUCTION
ASTM E135505a (2005), Standard Guide for Evaluating the Predictive Capability of
Deterministic Fire Models, American Society for Testing and Materials, West
Conshohocken, PA, 2005.
Buchanan, A. H. (2001), Structural Design for Fire Safety, John Wiley and Sons, LTD,
Chichester, England, 2001.
Babrauskas, V., Ignition Handbook, Fire Science Publishers/Society of Fire Protection
Engineers, Issaquah WA (2003).
Drysdale, D., An Introduction to Fire Dynamics, 3rd Ed., John Wiley, 2011.
Karlsson, B. and Quintiere, J.G., Enclosure Fire Dynamics, CRC Press, 2000
Quintiere, J.G., Principles of Fire Behavior, Delmar Publishers, 1998.
Quintiere, J.G., Fundamentals of Fire Phenomena, John Wiley, 2006.
NFPA, Fire Protection Handbook, National Fire Protection Association, 20th Ed., A.E.
Cote, (Editor) 2008.
NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, Nuclear Regulatory Commission, 2007.
SFPE, SFPE Engineering Guide to Assessing Flame Radiation to External Targets from
Pool Fires, SFPE Engineering Guide, Society of Fire Protection Engineers,
Bethesda, MD, March, 1999.
SFPE, SFPE Engineering Guide to Fire Exposures to Structural Elements, SFPE
Engineering Guide, Society of Fire Protection Engineers, Bethesda, MD,
November, 2005.
SFPE, SFPE Engineering Guide to Piloted Ignition of Solid Materials Under Radiant
Exposure, SFPE Engineering Guide, Society of Fire Protection Engineers,
Bethesda, MD, January, 2002.
SFPE, SFPE Engineering Guide to Predicting Room of Origin Fire Hazards, SFPE
Engineering Guide, Society of Fire Protection Engineers, Bethesda, MD,
November, 2007.
SFPE, SFPE Handbook of Fire Protection Engineering, 4th Ed., P. DiNenno (Editor),
National Fire Protection Association, 2008.
1-4
INTRODUCTION
In an ideal situation, once the ceiling jet reaches the enclosing walls, a Hot Gas Layer4 (HGL)
develops. As a result of the continuing supply of smoke mass and heat via the plume, the HGL
becomes deeper, and its temperature increases. Other properties of the smoke in the HGL also
increase (including concentration of gas species and solid particulates).
Radiant heat from the HGL to other combustibles not involved in the fire increases their
temperature. Similarly, the temperature of non-burning combustibles will also increase as a
result of receiving thermal radiation from the burning item(s). As the other combustibles reach
their respective ignition temperatures, they will also ignite. In some cases, the ignition of many
other combustibles in the space caused by heating from the HGL occurs within a very short time
span. This is commonly referred to as flashover.
Several aspects of fire behavior may be of interest when applying fire models, depending on the
purpose of the modeling application. Analysts may seek to determine the effects associated
with heating of targets submerged in smoke or receiving radiant heat from the flames, the
response of ceiling-mounted detectors or sprinklers to the fire environment, or other
phenomena.
Hot Gas Layer or HGL is also called smoke layer or hot upper layer in other publications in fire protection engineering.
1-5
INTRODUCTION
The most common aspects of fire behavior that typically are of interest in such analyses include,
but are not limited to:
A detailed review of each of these aspects is provided in texts on fire dynamics. A brief review
of each is provided here.
Rate of smoke production
Smoke is defined as a combination of the gaseous and solid particulates resulting from the
combustion process, plus the air that is entrained into the flame and/or smoke plume.
Consequently, the rate of smoke production at a particular height in the plume is the
combination of the generation rate of combustion products and air entrainment rate into the
flame and/or smoke plume between the top of the fuel and the height of interest. In most cases,
the air entrainment rate greatly exceeds the generation rate of fuel volatiles. Thus, the
correlations used to estimate the rate of smoke production are usually taken from experimental
research on entrained air.
Rate of smoke filling
The rate of smoke filling is dependent on the rate of smoke production, the heat release rate
(HRR), floor area, height and configuration of the space, and time from ignition. For a fire with a
steady HRR, the rate of smoke filling in a compartment will decrease with time due to a
decrease in the smoke production rate, which decreases as the height available to entrain air
decreases when the HGL deepens.
Properties of the ceiling jet
The ceiling jet transports smoke and heat horizontally away from the plume after it impacts with
the ceiling. The response of ceiling-mounted fire detectors or sprinklers is governed primarily by
their interaction with a ceiling jet. The temperature and concentration of smoke in a ceiling jet is
principally dependent on the height and configuration of the space, distance to the ceiling
impact point of the smoke plume, and the HRR of the fire.
Properties of the HGL
As smoke and heat are transported to the HGL via the smoke and fire plumes, the properties of
the HGL will change. The principal properties of interest include the depth, temperature, and
gas concentrations in the HGL. The magnitude of the properties depends on the HRR of the
fire, geometry of the space, ventilation openings (permitting material from the HGL to leave the
space, providing air to the fire, and/or causing a stirring action), yields of combustion products,
1-6
INTRODUCTION
and the elapsed time after ignition. These changes can be tracked by considering the
conservation of energy, mass, and species relative to the HGL.
Target response to incident heat flux via either thermal radiation or convection
The targets temperature will increase as a result of receiving heat via either thermal radiation or
convection. Radiation heat transfer is dependent on the intensity of thermal radiation emitted by
a source, the size of the source, and the proximity of the target to the source. For this
application, the flame height, the portion of heat released from the fire as radiation, and the
distance separating the target from the flame are the dominant parameters. Convective heating
occurs whenever the target is submerged in the smoke plume or HGL.
1-7
INTRODUCTION
HGL
Plume mass
flow
Mass
outflow
Mass inflow
Figure 1-2. A two-zone enclosure fire with an HGL above and a cool lower layer below.
Zone models are most applicable in situations involving simple geometries or where spatial
resolution within a compartment is not important. The preparation of input for a zone model, the
computation time, and the amount of output data generated are slightly more extensive than a
simple algebraic model; however, the overall computational time cost is still low.
Zone models can easily analyze conditions resulting from fires involving single compartments or
compartments with adjacent spaces, and are often used to compute the HGL temperature, HGL
composition, and target heat fluxes. They are also capable of modeling some effects of natural
and mechanical ventilation in both horizontal and vertical directions. Some zone models allow
the user to select a thermal plume model, which may assist in better characterization of a known
fire scenario, while others use an axisymmetric smoke plume. Other features of a zone model
may include a user-specified single zone or multiple fire plumes.
Simulations of spaces with complex ceilings or numerous compartments can be challenging
with a zone model. Because zone models specify uniform conditions in the HGL and lower
layer, results cannot be distinguished between locations at different distances from the fire. Due
to the zone approach, smoke transport time lags are not considered in the simulation, which is
an acceptable approximation in relatively small spaces but may lead to significant error in largevolume spaces or spaces with large aspect ratios.
Smoke production, fire plume dynamics, ceiling jet characteristics, heat transfer, and ventilation
flows are all algebraic models embedded within zone models. Other parameters that can be
calculated with a zone model include thermal behavior, detection response, and suppression
response. The output of a zone model is typically simple to understand and is generally
presented through an automatic user interface.
1-8
INTRODUCTION
Most zone models have default values that must be recognized and adjusted as necessary to
obtain an accurate solution. The model user must understand and justify the relevance of the
default values used in any application. Fire model users are expected to assess the
appropriateness of default values provided in the fire models and make changes or adjust
values as necessary. User manuals and technical references for each zone model outline such
values and may provide recommended ranges for the parameters.
Spatial resolution is important, relative to either the locations of fuel packages or targets.
Large compartments relative to the fire size are involved.
Compartments have complex geometries, flow connections, or numerous obstructions in
the upper part of the compartment.
Large numbers of compartments are within the area of interest and the presence of each
compartment is expected to affect the fire environment in the area of interest.
An example of a CFD simulation of a fire experiment is shown in Figure 1-3. The purpose of the
calculation was to simulate an experiment that was part of the validation study described in
NUREG-1824 (EPRI 1011999). In the experiment, a pan fire was placed in a relatively small
1-9
INTRODUCTION
compartment, and temperatures and heat fluxes were measured at various locations. The CFD
simulation is able to describe the changing behavior of the fire as it interacts with its
surroundings.
While CFD models provide a detailed analysis of a space, they are costly to create, simulate,
and maintain. The input files created in the preprocessing stage require a significant effort to
create. The user must understand the code syntax and the implications and approximations
embedded in the model. A firm understanding of fire dynamics is important in providing input
data that is relevant to the application. Most CFD models have default values that must be
recognized and adjusted as necessary to obtain an accurate solution. The model user must
understand and justify the relevance of the default values used in any application. Fire model
users are expected to assess the appropriateness of default values provided in the fire models
and make changes or adjust values as necessary. User manuals and technical references for
each CFD model outline such values and may provide recommended ranges for the
parameters.
Depending on the complexity of the scenario and the computers computational power, the
solver within the model can take anywhere from a few hours to weeks to complete all the
calculations. This time cost depends on the measured parameters, the size of the geometry,
and the mesh size of the calculations. Outputs of CFD models are visualized through a post-
1-10
INTRODUCTION
processing program. The CFD model developed at NIST, Fire Dynamics Simulator (FDS),
employs the program Smokeview to represent distributions of temperature, mass, heat flux,
burning rate, etc. throughout the geometry. These parameters can be described through point
locations, isocontours, or vector diagrams. Output data may also be stored in a commaseparated value file format that can be read by a standard spreadsheet program.
1-11
INTRODUCTION
The term fire-induced circuit failure is intended to address loss of function due to fire damage to cables and equipment as well as
spurious operation(s) of that equipment either singly or in multiples.
1-12
INTRODUCTION
The use of fire models acceptable to the authority having jurisdiction (AHJ), i.e., the U.S.
Nuclear Regulatory Commission
In the context of this application, the specific analytical capabilities within the fire model need to
be verified and validated. Model capabilities not invoked in a specific calculation are outside the
scope of this requirement. NUREG-1824 (EPRI 1011999) is an example of a verification and
validation study for fire models specifically developed for NPP applications.
Section 4.2.4.1 is subdivided by process element as follows:
Identify Targets (NFPA 805 4.2.4.1.1): This subsection requires the description of the
targets (e.g., equipment or cables) and target locations (specific locations of raceways/conduits
containing the cables, electrical cabinets, or equipment) associated with them needed to
achieve the nuclear safety performance criteria.
Establish Damage Thresholds (NFPA 805 4.2.4.1.2): This subsection requires the
description of damage thresholds for the equipment and cables needed to achieve the nuclear
safety performance criteria. The damage threshold (i.e., target vulnerability) for cables exposed
to fire is expressed in most cases in the form of an incident heat flux on the cables or the cables
surface temperature shall be established in accordance with Section 2.5 of NFPA 805.
Determine Limiting Conditions (NFPA 805 4.2.4.1.3): This subsection requires the
description of the combination of equipment or required cables with the highest susceptibility to
any fire environment.
Establish Fire Scenarios (NFPA 805 4.2.4.1.4): This subsection requires the description of
a given areas fire conditions resulting from the identified and analyzed fire scenarios. It should
be noted that the scenario definition is consistent with the requirements listed under 2.4.1.3 of
NFPA 805. Appendix C of NFPA 805 provides two categories of fire scenarios used in the
standard as follows:
Maximum Expected Fire Scenario: The maximum expected fire scenario (MEFS) is
defined in NFPA 805 Section C.3.2, as the scenario that is used to determine by fire
modeling whether the performance criteria are met in the fire area being analyzed. The
input data for the fire modeling of the MEFS should be based on the following:
1-13
INTRODUCTION
Heat release and fire growth rates for actual in-situ and transient combustibles
that are realistic and conservative based on available test data and applicable
fire experience
Ventilation within normal operating parameters with doors in the open or closed
position
Limiting Fire Scenario: The limiting fire scenario (LFS) is defined in NFPA 805 Section
C.3.3 as fire scenario(s) in which one or more of the inputs to the fire modeling
calculation (e.g., heat release rate, initiation location, or ventilation rate) are varied to the
point that the performance criterion is not met. The intent of this scenario(s) is to
determine that there is a reasonable margin between the expected fire scenario
conditions and the point of failure [T]he LFS can be based on a maximum possible,
though unlikely, value for one input variable or an unlikely combination of input variables.
The values used for LFS input should remain within the range of possibility but can
exceed that determined or judged to be likely or even possible.
For each fire scenario, the environmental conditions resulting from each MEFS are compared to
the damage thresholds for the targets in the fire area. If damage thresholds are not exceeded,
the targets in the fire area can be considered free of fire damage under the conditions of the
postulated MEFS.
By definition, the effects of the LFS include damage to the targets in the fire area under
consideration. Fire modeling parameters that have been varied to establish the LFS conditions
are identified and described.
INTRODUCTION
predetermined damage criteria. Damage states between the first and final states capture target
sets compromised as the fire propagates through intervening combustibles.
Each scenario progression postulated in a Fire PRA is quantified to determine its contribution to
fire risk. The most common fire risk metrics are the reactor Core Damage Frequency (CDF)
and Large Early Release Frequency (LERF) of radioactive material. NUREG/CR-6850 provides
guidance on the quantification of fire risk.
The Fire PRA standard (ASME/ANS RA-Sa-2008 and Addenda RA-Sa-2009) lists requirements
for all the technical elements associated with a Fire PRA, includes specific requirements for the
use of fire models. The standard addresses: (1) the selection of appropriate fire modeling tools
for estimating fire growth and damage behavior, considering the physical behaviors relevant to
the selected fire scenarios, and (2) the application of fire models that are sufficiently capable of
modeling the conditions of interest within known limits of applicability. In the case of analytical
fire models, the standard requires the use of appropriate fire modeling tools with the ability to
model the conditions of interest within known limits of applicability.
Chapter 2 presents a qualitative overview of the process for conducting fire modeling,
including the basic principles of fire simulation, advantages and limitations of the
technology, and brief descriptions of the five models
Chapter 3 provides specific guidance on selecting models to address typical scenarios in
commercial NPPs
Chapter 4 contains information on determining the sensitivity and uncertainty associated
with fire modeling calculations
Chapter 5 contains the list of references identified throughout this document
Chapter 6 contains general information on the eight scenarios documented in the
appendices
Appendices A through H provide detailed examples of fire modeling analyses of typical
NPP scenarios:
o Appendix A Cabinet Fire in Main Control Room
o Appendix B Cabinet Fire in Switchgear Room
o Appendix C Lubricating Oil Fire in Pump Compartment
o Appendix D Motor Control Center Fire in Switchgear Room
o Appendix E Trash Fire in Cable Spreading Room
o Appendix F Lubricating Oil Fire in Turbine Room
o Appendix G Transient Fire in Multi-Compartment Corridor
o Appendix H Cable Tray Fire in Annulus
CD containing the following files:
o Appendix A Files
Algebraic models
CFAST and FDS input files
o Appendix B Files
1-15
INTRODUCTION
o
o
o
o
o
Algebraic models
CFAST and FDS input files
Appendix C Files
FDS and MAGIC input files
Appendix D Files
CFAST and FDS input files
Appendix E Files
Algebraic models
CFAST and FDS input files
Appendix F Files
Algebraic model files
FDS files
Appendix G Files
Algebraic model files
MAGIC input files
Appendix H Files
FDS input files
Fire Modeling Software
CFAST Software Guides
FDS Software Guides
NUREG-1805 Spreadsheets
Related Publications
1-16
Definefiremodelinggoals
andobjectives(Section2.1)
Characterizethefire
scenarios(Section2.2)
Selectthefiremodel(s)
(Section2.3)
Yes
Canthe
problem
be
redefined
?
Yes
No
Is
another
model
available
?
No
No
Isthe
model
suitable?
Yes
Analternateanalysismaybe
required
Calculatethefiregenerated
conditions(Section2.4)
End
Conductasensitivityand
uncertaintyanalysis(Section
2.5)
Arethe
goalsand
objective
satisfied?
Yes
Documenttheanalysis
(Section2.6)
End
2-1
No
Evaluating whether or when a fire could damage a single electrical cable or component
Evaluating the potential for fire propagation through or across a fire barrier
Evaluating the potential for fire propagation between fire zones or fire areas, or to
secondary combustibles
The performance criteria are an important aspect for developing the fire modeling objectives
and will be specific to the fire modeling application. The performance criteria will often include,
but are not limited to one or more of the following:
Maximum acceptable incident heat flux for a cable, component, structural element, or
secondary combustible
2-2
The performance criteria may also involve sequences of events, such as detection or sprinkler
actuation before cable damage, which occurs when the surface temperature exceeds a
specified threshold value. NUREG/CR-6850 (EPRI 1011989) provides some performance
thresholds for common NPP targets (see Appendix H) as well as for habitability (see Section
11.5.2.11).
A few simple examples will illustrate the various ways in which a fire modeling goal may be
stated. In many NPP fire modeling applications, the motivation for a fire modeling analysis is a
need to know whether or not an electrical cable or a component remains free of damage from a
fire. This could be very specific (i.e., a particular exposure fire exposing a particular cable) or
general (i.e., the maximum distance from which a particular type of fire could damage cables).
In addition, it may only be necessary for a single fire to damage a single cable, or it may be
necessary to simultaneously damage two particular cables with a known separation. The
following are examples of goals that would lead to a fire modeling analysis:
Ensure that Panel X and Cable Y are not subject to failure given a fire in Room Z.
Tabulate the distance through which an electrical cabinet ignition source could
adversely affect electrical targets for use as a screening tool in support of a PRA.
Ensure that Cable X and Cable Y both do not fail in Room Z when exposed to a
single transient fuel package fire.
Quantify the potential benefit of the sprinkler system in Room X for preventing damage
to Cable Y.
Each of the goals explicitly states the purpose of the analysis and the means by which success
is determined in terms that can be achieved by a fire modeling analysis. The objectives are
more specific and may include specific steps that will be followed to satisfy the goal. Some
examples of objectives that support the previously cited example goals are as follows:
Evaluate the maximum distance from any surface of an electrical cabinet at which a 98th
percentile heat release rate (HRR) fire in Fire Zone X could cause the surface
temperature of a cable to exceed a specified threshold value.
Determine if any ignition sources in Fire Zone X could damage cables in Raceway Y
by causing the surface temperature of any cables in the raceway to exceed a specified
threshold value before the sprinkler system actuates.
2-3
A number of the fire scenario elements may also be viewed as fire model input. This section
provides a broad perspective on the considerations that apply when formulating the appropriate
fire scenario, given a fire modeling goal. Chapter 3 provides additional guidance on specific fire
scenario elements as they apply to various fire modeling goals and objectives evaluated with a
particular fire model.
Note that when characterizing the fire scenarios, preliminary consideration should also be given
to how many scenarios are needed to address a particular goal and which specific fire event
characteristics each scenario should capture (i.e., which scenarios are needed). In general, at
least one fire scenario would be necessary to assess the effects for a single ignition sourcetarget set pair. The analyst should become familiar with the information necessary to develop
input files for the fire modeling tools. In practice, this information should be collected during the
process of selecting and describing fire scenarios to minimize the number of walkdowns and
document/drawing reviews.
2-4
Targets in the fire plume or ceiling jet. Locating a source on top of a cabinet ignition
source usually results in the most severe fire conditions, since the cabinet walls will not
affect fire-generated conditions. Furthermore, since the fire is located in the highest
possible position, flames are expected to be higher, and temperatures in the plume and
ceiling jet will also be high. The user should judge whether this is conservative based on
2-6
Targets affected by flame radiation. Combustible materials that are not fixed, such as
transient fuel packages and unconfined liquid spills, should be located so that there is an
unobstructed (when no passive fire protection system is credited) view between the
source and the target. A horizontal path between flame and target provides the highest
heat flux to the target.
Targets immersed in the HGL. The fires elevation may influence how far down the HGL
will develop as predicted by some fire models, although other important scenario
characteristics will also be influential.
Fire detection systems. These include smoke detectors, heat detectors, or high
sensitivity detection systems
Passive fire protection systems. These include structural fire barriers, fire doors,
Electrical Raceway Fire Barrier Systems (ERFBSs), radiant shields, and fire stops
When assessing the performance of a system against the postulated fire hazard, it is necessary
to consider the conditions under which the system is designed (fire size, fuel load, exposure
temperature, plant operation mode, etc.). For example, an ordinary hazard sprinkler system
may not have a sufficient water spray delivery to protect against a large hydrocarbon pool fire.
Another example would be passive fire protection systems that are rated against an ASTM
2-7
2-9
When fire modeling is used to support a Fire PRA, the HRR for a source fire may be
represented as a conditional probability distribution. In this case, depending on the type of
analysis, a screening value may be selected (e.g., a 98th percentile peak HRR as recommended
in NUREG/CR-6850), or the effects may be represented using multiple points on the conditional
probability distribution.
As was the case with secondary combustibles, combustible materials in sealed or rated
containers may be excluded from consideration if the container is capable of resisting the
effects of the fire. In addition, self-ignited cable fires are generally postulated only for non-IEEE383 qualified power cables (NUREG/CR-6850/EPRI 1011989).
2-10
2-12
Examples
Typical Applications
Advantages
FDTS
FIVERev1
Screening calculations;
ZOI; target damage by
thermal radiation, HGL,
or thermal plume acting
in isolation.
Zone Model
CFAST
MAGIC
Computation
Fluid
Dynamics
Model
FDS
Disadvantages
Limited
application range;
treats phenomena
in isolation;
typically
applicable only to
steady state or
simply defined
transient fires
(e.g., proportional
to the square of
time or t fires).
Error increases
with increasing
deviation from a
rectangular
enclosure; large
horizontal flow
paths not well
treated.
Significant effort
to create input
files and postprocess the
results; long
simulation times;
difficult to model
curved geometry,
smoke detector
performance, and
conditions after
sprinkler
actuation.
An important consideration in the fire model selection process is the type of analysis performed.
Because of the large number of potential scenarios in a typical NPP, it is usually not practical to
default to the most sophisticated tool available. The analysis typically begins with a series of
screening analyses (NUREG/CR-6850/EPRI 1011989) performed to identify a subset of fire
scenarios and targets that require further analysis with greater resolution. The screening
process will typically use fairly simple fire modeling tools, such as algebraic models or generic
solutions. When such screening is conducted, it is important to remain within the model
limitations and the verification and validation (V&V) basis for the screening model. Section 2.3.6
and Chapter 4 provide additional guidance on the significance of the fire model V&V basis and
steps that the user should take to ensure that the fire model is used within acceptable limits.
2-13
02.1_Temperature_NV.xls
02.2_Temperature_FV.xls
02.3_Temperature_CC.xls
03_HRR_Flame_Height_
Burning_Duration_Calculation.
xls
2-14
NUREG-1824
Verification
and Validation
Status
V&V provided
V&V provided
V&V provided
NUREG-1824
Verification
and Validation
Status
04_Flame_Height_Calculations
.xls
05.1_Heat_Flux_Calculations_
Wind_Free.xls
05.2_Heat_Flux_Calculations_
Wind.xls
05.3_Thermal_Radiation_From
_Hydrocarbon_Fireballs.xls
06_Ignition_Time_Calculations.
xls
07_Cable_HRR_
Calculations.xls
08_Burning_Duration_Soild.xls
09_Plume_Temperature_
Calculations.xls
10_Detector_Activation_
Time.xls
13_Compartment_Flashover_
Calculations.xls
14_Compartment_Over_
Pressure_Calculations.xls
15_Explosion_Claculations.xls
2-15
16_Battery_Compartment_
Flammable_Gas_Conc.xls
17.1_FR_Beams_Columns_
Substitution_Correlation.xls
17.2_FR_Beams_Columns_
Quasi_Steady_State_Spray_
Insulated.xls
17.3_FR_Beams_Columns_
Quasi_Steady_State_Board_
Insulated.xls
17.4_FR_Beams_Columns_
Quasi_Steady_State_
Uninsulated.xls
18_Visibility_Through_Smoke.
xls
NUREG-1824
Verification
and Validation
Status
2-16
CThrr
Function Description
Heat release rate profile considering t growth and four stages.
Estimates flame irradiation at distance r from the fire source.
Point source approximation for REMOTE targets.
Flame height based on Heskestad flame height correlation.
Plume temperature at a specific height based on Alpert plume
temperature correlation.
Plume temperature at a specific height based on McCaffrey plume
temperature correlation.
Plume temperature at a specific height based on Heskestad
plume temperature correlation.
Estimates convective heat flux in the fire plume.
Plume velocity at a specific height based on Alpert's plume
temperature correlation.
Plume velocity at a specific height based on McCaffrey plume
temperature correlation.
Plume velocity at a specific height based on Heskestad plume
temperature correlation.
Air entrainment into plume based on Zukoski plume entrainment
correlation.
Air entrainment into plume based on Thomas plume entrainment
correlation.
Air entrainment into plume based on Heskestad plume
entrainment correlation.
Estimates plume diameter based on Heskestad plume correlation.
Unconfined ceiling jet temperature based on Alpert ceiling jet
correlation.
Confined ceiling jet temperature based on Delichatsios ceiling jet
correlation.
Estimates convective heat flux in the ceiling jet.
Unconfined ceiling jet velocity based on Alpert ceiling jet
correlation.
Compartment temperature after a specified time, given a steady
HRR based on MQH approach.
Heat release rate required for flashover after a specified time
based on MQH approach.
Estimates compartment temperature based on FIVE-Rev1.
Activation time of heat detection devices based on heat release
rate profiles.
Time required by Hot Gas Layer to reach a specific height based
on heat release rate profiles and openings at the bottom of the
enclosure.
Estimates heat release rate from cable trays. The correlation is
based on 14 experiments with a stack of 12 horizontal cable trays
and 2 experiments with a combination of 12 horizontal cable trays
and 3 vertical trays.
2-17
NUREG-1824
Verification
and Validation
Status
V&V not provided
V&V provided
V&V provided
V&V not provided
V&V not provided
V&V not provided
V&V not provided
V&V not provided
V&V provided
V&V provided
V&V not provided
V&V not provided
V&V not provided
V&V not provided
V&V provided
V&V not provided
V&V not provided
V&V not provided
V&V provided
V&V not provided
V&V provided
V&V not provided
V&V not provided
Function Description
Estimates the length of a visible path in a smoke environment.
The correlation applies to light-reflecting signs.
Estimates target temperature under constant heat flux.
Time to target damage under constant heat flux.
NUREG-1824
Verification
and Validation
Status
V&V not provided
V&V not provided
V&V not provided
2.3.4 MAGIC7
MAGIC is another two-zone computer fire model, developed and maintained by Electricit de
France (EdF) specifically for use in NPP analysis. MAGIC is supported by three EdF
publications, including: (1) the technical manual, which provides a mathematical description of
the model (Gay et al., 2005b); (2) the users manual, which details how to use the graphical
interface (Gay et al., 2005a); and (3) the validation studies, which compare MAGICs results to
experimental measurements (Gay et al., 2005c). These three proprietary publications and the
MAGIC software are available through EPRI to EPRI members. Additional V&V results for
MAGIC are documented in Volume 6 of NUREG-1824 (EPRI 1011999).
MAGIC is fundamentally the same type of model as CFAST and thus solves the same basic set
of differential equations. The combustion model and vent flow models are similar as well.
Despite this, MAGIC differs from CFAST in that it does not have the corridor or shaft submodels, and the ceiling jet and wall jet treatments are different. The user should consult the
technical manual for a complete description of the MAGIC sub-models (Gay et al., 2005b).
Once a given simulation is completed, MAGIC generates an output file with all of the solution
variables. Through a post-processor interface, the user selects the relevant output variables
for the analysis. Typical outputs include the temperatures of hot and cold zones, concentrations
of oxygen and unburned gases, smoke migration into each compartment, the mass flow rates of
air and smoke through the openings and vents, the pressures at the floor level of each
compartment, the temperatures at the surfaces of the walls, and the thermal fluxes (radiative
and total) exchanged by the targets placed by the user.
2-19
2-20
(2-1)
where is the fire HRR (kW), is the ambient density of air (kg/m3),
is the specific heat of
air (kJ/kg/K), is the ambient air temperature (K), and is the acceleration of gravity (m/s2).
In many instances,
is comparable to the physical diameter of the fire. FDS employs a
numerical technique known as large eddy simulation (LES) to model the unresolvable or subgrid motion of the hot gases. The effectiveness of the technique is largely a function of the ratio of
the fires characteristic diameter, , to the size of a grid cell, . In short, the greater the ratio
/ , the more the fire dynamics are resolved directly, and the more accurate the simulation.
Past experience has shown that a ratio of 5 to 10 usually produces favorable results at a
moderate computational cost for problems where gross smoke movement is of interest.
As an example, suppose the HRR of the fire were 700 kW. The characteristic diameter may
then be calculated as follows:
/
700 kW
1.2 kg/m
1.012 kJ/kg/K
0.83 m
(2-2)
To perform a grid sensitivity analysis, a good place to start might be a cell resolution of 15 cm
5.5. Then choose a grid of 10 cm (4 in), and then 5 cm (2 in).
(6 in), which means that /
At this point, the calculation time will have increased by a factor of roughly three hundred,
making it potentially impractical to compute; however, if it can be shown that there is little
difference between the 5 cm (2 in) and 10 cm (4 in) grids, then the objective has been achieved.
The meaning of little difference can be interpreted in several ways. Given that NUREG-1824
(EPRI 1011999), the fire model V&V study, lists the relative error expected of the various
models for the various quantities, it is reasonable to interpret the difference in results on
different grids in light of what is expected of the model accuracy.
Although the fire size and dimensions often determine the optimum grid resolution, there are
other factors that can influence the selection of the grid resolution. These include the number of
cells used to resolve a flow path dimension, the number of cells used to describe the fire
dimension, and the number of cells used to resolve the conditions in a partially isolated volume.
These considerations are related in that it is generally advisable to include at least three cells
across any flow path, such as a door or a window, and fire dimension, regardless of the
minimum number of cells computed using the fire characteristic diameter. In some cases,
partially isolated volumes are created by various obstructions; if the temperature and flow
conditions are of interest in these areas, a minimum of three cells across any dimension should
be provided. Another consideration that could influence the grid resolution is the dimension of
2-21
2-22
Fire Model
FDTs
FIVE-Rev1
CFAST
MAGIC
FDS
YES
YES
YES
YES
YES
NO
NO
YES
YES
YES
NO
YES
YES
YES
YES
Plume Temperature
YES
YES
NO
YES
YES
Flame Height
YES
YES
YES
YES
YES
YES
YES
YES
YES
YES
NO
NO
YES
YES
YES
NO
NO
YES
YES
YES
Wall Temperature
NO
NO
YES
YES
YES
Target Temperature
NO
NO
YES
YES
YES
Smoke Concentration
NO
NO
YES
YES
YES
Oxygen Concentration
NO
NO
YES
YES
YES
Room Pressure
NO
NO
YES
YES
YES
Table 2-5 identifies normalized parameters that may be used to compare NPP fire scenarios
with validation experiments. A full description of each parameter follows Table 2-5. The
validation range for the normalized parameters shown in Table 2-5 were derived from NUREG1824 (EPRI 1011999), Table 2-4, and are intended to provide guidance on which groups of
validation experiments to consider when evaluating a certain attribute based on the validation
results. These parameters may not be the only ones appropriate for evaluating the applicability
of a specific experiment; Table 2-5 of NUREG-1824 (EPRI 1011999) Volume 1 lists the ranges
of values for different physical characteristics and normalized parameters based on the
experiments considered in the validation study.
It is seen in Table 2-5 that the fire diameter plays a role in three of the normalized parameters:
the fire Froude Number, the flame length ratio, and the radial distance ratio. The fire diameter
may be estimated from the total plan area of the burning fuel using the following equation
(SFPE, 2011):
4
(2-3)
where is the effective fire diameter (m) and is the plan area of the burning fuel (m) equal to
the fuel package length dimension multiplied by the width dimension. Equation 2-3 produces a
diameter of a circle having an equivalent area to that of the burning fuel and is therefore
2-23
Quantity
Normalized Parameter
Fire Froude
Number
Flame Length
Ratio
3.7
1.02
cj
Ceiling Jet
Distance Ratio
Equivalence
Ratio
O
Compartment
Aspect Ratio
General Guidance
0.23
1
2
0.23
(Natural)
Mechanical
or
Radial Distance
Ratio
2-24
Ratio of characteristic
velocities. A typical accidental
fire has a Froude number of
order 1. Momentum-driven fire
plumes, like jet flares, have
relatively high values.
Buoyancy-driven fire plumes
have relatively low values.
A convenient parameter for
expressing the size of the fire
relative to the height of the
compartment. A value of 1
means that the flames reach
the ceiling.
Ceiling jet temperature and
velocity correlations use this
ratio to express the horizontal
distance from target to plume.
The equivalence ratio relates
the energy release rate of the
fire to the energy release that
can be supported by the mass
flow rate of oxygen into the
compartment, O . The fire is
considered over- or underventilated based on whether
is less than or greater than 1,
respectively.
This parameter indicates the
general shape of the
compartment.
This ratio is the relative
distance from a target to the
fire. It is important when
calculating the radiative heat
flux.
NUREG1824
Validation
Range
0.4 2.4
0.2 1.0
1.2 1.7
0.04 0.6
0.6 5.7
2.2 5.7
is the Fire Froude number (non-dimensional), is the fire heat release rate (kW),
where
is the heat capacity at constant pressure for ambient air
is the density of ambient air (kg/m),
(kJ/kg-C), is the ambient temperature (K), is the fire diameter or the effective fire diameter
(m), and is the acceleration of gravity (9.81 m/s). The heat release rate that is used to
generate the validation range is the peak heat release rate; thus, when comparing the model
application range to these values, the peak heat release rate should be used. The effective fire
diameter may be computed using Equation 2-3. The ambient density and temperature are
typically coupled through the following equation:
352
(2-5)
where all terms have been defined. The heat capacity for air varies slightly with temperature,
and, at 298 K, it is equal to about 1.012 kJ/kg-C.
Flame Length Ratio
The flame length ratio is a measure of the flame height relative to the upper horizontal boundary
(ceiling). The NUREG-1824 validation tests all involved fires with a flame height that was at or
below the ceiling. In situations where the flame height is greater than the ceiling, the flames will
extend radially outward from the impingement point. This configuration is not well characterized
by simple algebraic models, or even by zone models that use plume and flame height
correlations.
Flame Length Ratio
(2-6)
where
is the base height of the fire (m),
is the flame height (m), and
is the enclosure
height (m). Note that the flame length ratio does not apply to fires modeled outside an
enclosure. The flame height is computed using the following equation:
3.7
1.02
(2-7)
where is the fire diameter or the effective fire diameter (m) and
is the fire Froude Number
calculated using Equation 2-4 (non-dimensional). The effective fire diameter may be computed
using Equation 2-3.
2-25
(2-8)
is the
where cj is the horizontal distance within the ceiling jet from the fire centerline (m),
is the base height of the fire (m). The ceiling jet ratio is applicable
enclosure height (m), and
primarily when sprinkler or heat detector actuation is calculated. It is not applicable when the
fire is modeled outside an enclosure.
Equivalence Ratio
This quantity is the ratio of the generation rate of fuel to the supply rate of oxygen. When the
equivalence ratio is equal to one, the exact amount of oxygen required for complete combustion
is available. When the ratio is greater than one, the environment is fuel rich and the fire is
considered to be under-ventilated. The reverse is true when the ratio is less than one. The
upper limit for the equivalence ratio in enclosure fires is about three (Gottuk et al., 2008).
Modeling under-ventilated fires is challenging in part because the fuel mass loss rate is typically
specified by the user and is not adjusted by the fire models to the conditions actually present.
The equivalence ratio may be estimated for natural and forced ventilation by the expression:
(2-10)
O
where is the equivalence ratio (non-dimensional), is the heat release rate of the fire (kW),
is the heat of combustion for oxygen (kJ/kg), and O is the mass flow rate of oxygen
into the enclosure (kg/s). The heat of combustion for oxygen is the energy released per unit
mass of oxygen and varies by less than five percent among nearly all carbon-based fuels
(Janssens, SFPE Handbook, 4th edition). The value for
is typically taken to be 13,100
kJ/kg, an average value over a wide range of common fuels. A notable exception to this
convention is hydrogen, whose heat of combustion based on oxygen consumption is
approximately 18,000 kJ/kg. For other pure fuels, it is possible to calculate a more accurate
value than 13,100 kJ/kg, in which case Eq. (2-10) can be modified accordingly.
The mass flow rate of oxygen into the enclosure is given by the following equation:
0.23
0.23
1
2
(Natural)
Mechanical
2-26
(2-11)
(2-12)
In many compartments, both natural and forced ventilation is present. The recommended
procedure to follow in this case is to select the dominant ventilation mode (i.e., the mode that
produces the highest mass flow of oxygen). If the two modes are comparable, the oxygen
masses may be added. As was the case for the ceiling jet ratio and the flame length ratio, the
equivalence ratio applies only to fires that are modeled within an enclosure.
Compartment Aspect Ratio
The compartment aspect ratio is a measure of the deviation of the enclosure dimensions from a
cube. When at least one of the compartment aspect ratios is large, the enclosure takes on the
characteristics of a corridor. In such cases, the transport time of the combustion products and a
non-uniform layer can both become significant parameters that require consideration. When at
least one of the compartment aspect ratios is low, the enclosure takes on the characteristics of
a shaft. In these cases, stratification of the combustion products, the interaction of the fire
plume and the enclosure boundaries, or choked flow could become parameters that influence
the results. These situations may lie outside the development basis for algebraic and zone fire
models. The compartment aspect ratio is computed using the following equation:
2-27
(2-14)
(2-15)
is
For a given set of experiments and NPP fire scenarios, the user can calculate the relevant
normalized parameters. If the fire scenario parameters fall within the ranges evaluated in the
study, then the results of the study offer appropriate validation for the scenario. If they fall
outside the range, then a validation determination cannot be made based on the results from
the study. For any given fire scenario, more than one normalized parameter may be necessary
for determining the applicability of the validation results.
The V&V study provides valuable insight into the predictive capability of the five fire models.
This insight is ultimately characterized in terms of a bias and a standard deviation for a number
of output parameters. The closer the bias is to unity, the more accurate the fire model tends to
be in predicting the given parameter; the smaller the standard deviation, the smaller the
expected scatter around the mean bias. Chapter 4 of this guide describes how the V&V
uncertainty information can be used to assign a probability function to the output data.
NUREG-1824 (EPRI 1011999) provides V&V documentation for specific versions of fire models.
Because the fire models considered are under active development, new releases are expected
to and do occur. The user has the option of using the model version that has been verified and
validated in NUREG-1824 (EPRI 1011999) or re-evaluating cases in NUREG-1824 (EPRI
1011999) to demonstrate that the predictive capability of the model has not decreased for the
application at hand. It is expected that NUREG-1824 (EPRI 1011999) will be updated from time
to time, as the need arises.
Froude number: The two parameters that can be practically varied are the fire diameter
and the HRR. For fire sizes (i.e., HRR) that are small for the postulated diameter, the
resulting Froude number can fall under the low end of the applicability range. Similarly,
for fires that are relatively large for the postulated diameter, the Froude number can fall
above the applicability range. In the former situation, the analysts may consider
reducing the fire diameter and keeping the HRR profile unchanged. In most fire
modeling tools, the fire diameter is simply used to determine HRRs or to calculate the
fire plume conditions, such as the flame height or plume temperature. Considering that
the HRR is fixed in this sensitivity study, the fire diameter may not be a relevant
parameter in the analysis, with the important exception of scenarios where the fire plume
conditions are relevant. A similar approach could be used for the latter situation.
Increasing the fire diameter can force the dimensionless term into range. It should be
stressed that fire diameter is often a parameter that influences predicted flame height
and fire plume conditions, and that the effects of diameter variations should be explicitly
addressed in the analysis. This includes other dimensionless terms where the fire
diameter is a key input (e.g., target distance to diameter (r/D), etc.).
Flame length relative to ceiling height: This is a convenient parameter for expressing the
size of the fire relative to the height of the compartment. A value of 1 means that the
flames reach the ceiling. The validation range extends up to a value of 1.0, which
should cover most of the scenarios of interest in commercial NPPs. Scenarios that are
expected to fall out of the range are:
o
Those associated with relatively short flames. Typical ceiling heights in NPP
scenarios range from about 3 to 6.1 m (10 to 20 ft), excluding the containment and
turbine buildings, which have relatively large openings between elevations.
Consequently, flame lengths shorter than 0.6 to 1.2 m (2 to 4 ft) will be considered
outside of validation range. A sensitivity analysis increasing the HRR values should
provide a conservative estimate of fire conditions within the validation range. In
cases where the conclusion of the analysis does not change given the increased fire
intensity (e.g., no damage within the flame length of fire plume), the suggested
sensitivity analysis can be used as the justification for the evaluation of a
compartment that falls outside the validation range.
2-29
Flame extensions under ceilings. In this particular case, not only are such flame
lengths out of the range of validation, but also the models for predicting this
phenomenon have not been verified or validated with a process similar to the one
documented in NUREG-1824 (EPRI 1011999).
Ceiling jet radial distance relative to the ceiling height: Ceiling jet temperature and
velocity correlations use this ratio to express the horizontal distance from target to
plume. Ceiling jet applications in commercial NPPs should be carefully evaluated due to
the numerous obstructions near the ceiling (e.g., cable trays, HVAC ducts, piping, etc.).
Most of its applications include determination of time to detection and sprinkler
activation, in which the ceiling jet velocity is a sub-model in the analysis. An alternative
option is a sensitivity analysis consisting of moving the fire location to distances that
would fall within the validation range; it is recognized, however, that in many situations
the fire location cannot be altered, particularly in the case of fixed ignition sources or
transient fires postulated near areas where redundant targets are in close proximity
(pinch-points). In general, longer horizontal distances will result in longer activation time
results; by contrast, shorter horizontal distances would result in conservative time-todamage results. In situations where the ceiling jet geometry deviates significantly from
the idealized flat horizontal surface, as may be the case when there are large numbers
of obstructions or bays, a CFD model may be the better choice for calculating detector
response times.
Equivalence ratio as an indicator of the ventilation rate: The validation available is for
well-ventilated fires: that is, no model validation information is available for underventilated compartment fires, including fire extinction due to lack of oxygen. In general,
fires that are considered well ventilated in the enclosure should result in bounding
conditions as long as the HRR profile is appropriate. Conditions in the enclosure are not
expected to be worse in a fire where the combustion process is affected by lack of
oxygen than they would be under fire conditions where the combustion process is
unaffected. However, under-ventilated fire conditions should be considered carefully as
sudden air inflows into compartments with under-ventilated fire conditions could produce
relatively severe fire conditions.
2-30
Radial distance relative to the fire diameter: This ratio is the relative distance from a
target to the fire, and is important when calculating the radiative heat flux. Note that the
validation range starts at a distance approximately twice the fire diameter. In practice,
targets at a very close distance to the fire (approximately two fire diameters or less)
should be expected to fail, given the relatively low damage threshold levels for cables.
An alternative option is a sensitivity analysis, which consists of moving the fire location to
distances that would fall within the validation range; it is recognized, however, that in
many situations the fire location cannot be altered, particularly in the case of fixed
ignition sources or transients fires postulated near pinch-points. In general, shorter
horizontal distances will result in higher heat flux levels.
Scenarios involving targets within the fire plumes: A useful discussion of fire plumes is
contained in Gunnar Heskestads chapter in the SFPE Handbook of Fire Protection
Engineering, 4th ed., Fire Plumes, Flame Height, and Air Entrainment. The plume
correlations used in the empirical and zone models are described, as well as their range
of applicability. NUREG-1824 (EPRI 1011999) contains experimental measurements of
fire plumes, but the range is somewhat limited. The plume correlations used by the
models have a much wider range of applicability than that exercised in NUREG-1824
(EPRI 1011999).
Scenarios involving targets within the ceiling jet: Similarly, Ronald Alperts chapter
Ceiling Jet Flows in the SFPE Handbook contains a description of the various
correlations used to estimate the temperature and gas velocity of ceiling jets. There are
extensive references to the original experimental test reports from which the correlations
were derived.
Scenarios involving electrical failure of cables: The Cable Response to Live FIRE
(CAROLFIRE) program led to the development and validation of the Thermally-Induced
2-31
Scenarios involving cable burning: The Cable Heat Release, Ignition, and Spread in
Tray Installations in Fire (CHRISTIFIRE) program led to the development and validation
of the Flame Spread in Horizontal Cable Trays (FLASH-CAT) model (NUREG/CR-7010,
Volume 1). This model addresses the growth and spread of fire within vertical stacks of
horizontal, open-back cable trays.
In addition to NUREG-1824 (EPRI 1011999) and the various documents cited above, the
individual model developers typically maintain a collection of validation cases that are included
as part of the model documentation. The algebraic spreadsheet models, FDTs and FIVE-Rev1,
are based directly on experimental correlations. Validation of these models is typically not part
of the model documentation; rather, there are references to source material like the SFPE
Handbook or the original test reports. Validation studies by the CFAST and FDS developers are
contained within:
NIST Special Publication 1086, CFAST Consolidated Model of Fire Growth and Smoke
Transport, Software Development and Model Evaluation Guide, 2008.
NIST Special Publication 1018, Fire Dynamics Simulator, Technical Reference Guide,
Volume 3, Validation, 2007.
In summary, the purpose of the sensitivity analysis is to re-shape the scenario with parameters
that fall within the V&V range and result in more severe fire generated conditions (e.g., higher
HGL or plume temperature, higher incident heat flux, etc.). Depending on the application, one
or more parameters may need to be varied affecting multiple dimensionless parameters. It is
recommended that the results from the sensitivity calculations always be compared to those
resulting from the base case to ensure that the input parameter manipulation produces more
severe fire generated conditions.
2-32
2-33
Models are developed based on idealizations of the physical phenomena and simplifying
approximations, which unavoidably introduces the concept of model uncertainty (i.e.,
model error) into the analysis.
The concepts of model and parameter uncertainty have traditionally been addressed in fire
modeling using uncertainty and/or sensitivity analysis. The uncertainty in a variable represents
the lack of knowledge about the variable, and is often represented with probability distributions.
Its objective is to assess the variability in the model output, that is, how uncertain the output is
given the uncertainties related to the inputs and structure of the model. By contrast, the
sensitivity of a variable in a model is defined as the rate of change in the model output with
respect to changes in the variable. A model may be insensitive to an uncertain variable.
Conversely, a parameter to which a model is very sensitive may not be uncertain.
Details of the uncertainty and sensitivity analysis are included in Chapter 4.
Marked-up plant drawings. Plant layout, detection, suppression, cable tray, HVAC, and
conduit drawings are often marked to highlight the location of the compartment, the
ignition sources, the targets, the ventilation flow paths, and the fire protection features.
The drawings also serve as sources of fire model input values, such as compartment
2-34
Design basis documents (DBDs). Design basis documents (DBDs) provide in-depth
assessments of plant features in various operation modes, such as the HVAC system.
Sketches. Sketches are perhaps one of the most useful ways of documenting a fire
scenario. A sketch typically consists of a drawing illustrating the ignition source,
intervening combustibles, targets, and fire protection features. A first draft of the sketch
is usually prepared during walkdowns. The analyst should take the opportunity to
include details such as raceways and conduit identifications (IDs), and other information
relevant to the fire modeling analysis. Pictures often supplement sketches.
Write-ups and input tables. Write-ups and input tables are used to compile the
information collected from drawings and walkdowns in an organized way. The write-up
should include a brief scenario description and detailed documentation supporting
quantitative inputs to the fire modeling analysis, as well as any relevant sketches or
pictures associated with each scenario.
Software versions, descriptions, and input files. The documentation package should
include the version numbers of any software, brief descriptions of the software, and
copies of the input files.
The examples presented in Appendices A through H of this guide illustrate techniques for the
proper documentation of fire modeling calculations using the format described. In conclusion, a
properly documented analysis should enable someone else to reproduce the results from the
information contained within the documentation.
2.7 Summary
This chapter described a recommended process for conducting and documenting a fire
modeling analysis. Chapter 3 provides guidance on selecting the appropriate fire modeling tool
and input parameters for typical commercial NPP applications. Fire model uncertainty is
addressed in Chapter 4 of this document. Specific fire modeling examples evaluated using the
process described in this Chapter are provided in Appendices A through H.
2-35
2-36
Ventilation Parameters
Targets
Intervening Combustibles
3-1
Fire elevation: The elevation of the base of the fire, measured from the floor. It is important
in scenarios involving targets in the fire plume where the relative distance between the fire
and the target strongly influences the exposing temperature. It is also important because
the height of the fire relative to the hot gas layer (HGL) influences the air entrainment into
the plume, the position of the HGL, and, potentially, the actual HRR (since air entrained from
the HGL is oxygen-depleted).
Fire location: In scenarios where the fire is located near a wall or corner, the plume is
expected to entrain less air, resulting in higher plume temperatures (Karlsson and Quintiere,
2000, p. 72).
Fuel mass: This parameter is an important factor in determining the burning duration
Soot and product yields: The yield of a combustion product is the mass of the product
generated per unit mass of fuel consumed. In particular, the soot yield is an important factor
in radiative heat transfer (e.g., targets immersed in the HGL), visibility calculations, and
smoke detector response estimates. The yields of toxic gases can also be important in
habitability calculations
Radiative fraction: The fraction of energy emitted in the form of thermal radiation. For most
materials, the radiative fraction is approximately one third. That is, one third of the total
HRR radiates in all directions and two thirds is convected upwards into the smoke plume
3-2
Thermal
Conductivity
(W/m/K)
0.8
1.6
386
0.17
0.12
0.192
54
0.235
Density
(kg/m3)
Specific Heat
(kJ/kg/K)
Source
2600
2400
8954
960
540
1380
7850
1375
0.8
0.75
0.38
1.1
2.5
1.289
0.465
1.390
,
,
(3-1)
Here Ao,i and Ho,i are the individual door areas and heights. The effective width of multiple
vertical openings can be estimated by the ratio Ao/Ho.
Regarding doors (and other operable openings), consideration should be given to the doors
being opened (or closed) during a fire. For example, when the fire brigade arrives, they will
open the doors to the fire area to gain access, which will affect the ventilation and possibly
result in smoke spread.
3-3
3.1.4 Targets
A target is an object of interest that can be affected by the fire-generated conditions and
typically consists of cables in conduits, cables in raceways, or plant equipment. Targets are
characterized by their location, damage criteria, and thermophysical properties.
A targets location simply refers to its location relative to the fire. The location is represented by
three-dimensional coordinates within the volume of the room in which the fire conditions are
simulated. Where the target faces in a particular direction, an orientation vector to indicate that
direction needs to be entered.
The damage criteria refer primarily to a damage/response threshold. In general, the damage
criteria for scenarios involving cable damage is expressed in terms of damage temperature or
incident heat flux.
The models within the scope of this Guide require specification of the targets thermophysical
properties, primarily the density, specific heat, and thermal conductivity, for the analysis. These
parameters are used to estimate heat conducted into the targets. The predicted time for the gas
temperature surrounding a target to reach a specific limit is usually less than the time it takes
the target to reach the same limit because the heat conduction inside the target will delay the
temperature rise at the surface during the heating process.
3-4
Chapter Section
3.2.1
3.2.2
3.2.3
3.2.4
3.2.5
6
7
3.2.6
3.2.7
Scenario Description
Scenarios consisting of determining time to damage of cables above the ignition
source located inside the flames or the fire plume.
Scenarios consisting of determining time to damage of cables located inside or
outside the HGL. This scenario also includes a secondary fuel source (i.e.,
propagation to cable trays).
Scenarios consisting of determining time to damage of cables located in a room
adjacent to the room of fire origin.
Scenarios consisting of determining time to damage of cables located inside or
outside the HGL in rooms with complex geometries.
Scenarios consisting of determining time to loss of habitability of the main control
room.
Scenarios consisting of determining time to smoke or heat detector activation.
Scenarios consisting of determining temperature of structural elements.
3-5
A sketch capturing most of the technical elements relevant to the analysis. A legend
summarizing the different elements presented in the sketches is provided in Figure 3-2.
A scenario objective stating the purpose of the modeling exercise in engineering terms.
A modeling strategy section summarizing the recommended steps for performing the
calculation.
A section referencing relevant detailed fire modeling examples documented in the Appendix
section of this guide.
3
7
4
1
6
Figure 3-1. Pictorial representation of the fire scenario and corresponding technical
elements described in this section.
3-6
Cable
Tray
1. Exposed
2. Target
1.
Open
Door
2.
Junction
Box and
Conduit
Closed
Door
Electrical
Cabinet
Pump
Fire
Plant
Employee
Sprinkler
Control
Board
HVAC
Smoke
Detector
1. Supply
2. Exhaust
Structural
Element
1.
2.
Obstruction
Figure 3-2. Legend for fire modeling sketches presented in this chapter.
3-7
3-9
3-10
The more detailed models, such as FDS, can handle the ignition and contribution of multiple
fires, provided that the ignition criteria and source HRR characteristics are provided as input.
Other models, especially the algebraic models, only accept the total HRR as a function of time,
which is found by summing up the individual HRRs.
In the present example, consider a cable tray directly above the fire. The time to ignition of the
cable tray can be determined via algebraic models that estimate the flame height and plume
temperature as a function of time for the initial cabinet fire (see scenario previously discussed).
Once the flames from the cabinet reach the cable tray, the cable ignites. The same is true when
the plume temperature at the elevation of the cable tray reaches the ignition temperature of the
cables. Both calculations should be completed, and the shorter time used as the ignition time.
The HRR from the cable tray can be added to the HRR of the cabinet to determine a combined
HRR as a function of time. The resulting HRR profile can take into consideration both the fuel
consumption and propagation to additional intervening combustibles as a function of time. In
cases where fuel consumption is not considered, the resulting HRR profile is expected to
overestimate the fire intensity. This total rate can then be used in the various models as an
approximation of the HRR profile.
Appendix R of NUREG/CR-6850 (EPRI 1011989) addresses cable fires, including methods for
calculating the HRR for a variety of cable configurations.
It should be noted that the simple summation of the two HRRs is a simplification of a complex
phenomenon and only provides an approximation of the conditions created by the two separate
fires.
3-12
3-13
3-14
3-15
3-16
3-17
3-19
3-20
3-21
3-22
3-23
3-24
4-1
FDTs
1.44
0.25
HGL Depth*
N/A
N/A
FIVE-Rev1
1.56
0.32
N/A
CFAST
MAGIC
FDS
Exp
1.06
0.12
1.01
0.07
1.03
0.07
0.07
1.04
0.14
1.12
0.21
0.99
0.07
0.07
1.84
0.29
1.15
0.24
1.01
0.08
1.04
0.08
0.08
0.73
0.24
0.94
0.49
1.25
0.28
1.01
0.07
1.15
0.11
0.07
Flame Height**
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
I.D.
Oxygen Concentration
N/A
N/A
0.91
0.15
0.90
0.18
1.08
0.14
0.05
Smoke Concentration
N/A
N/A
2.65
0.63
2.06
0.53
2.70
0.55
0.17
N/A
N/A
1.13
0.37
0.94
0.39
0.95
0.51
0.20
N/A
N/A
1.00
0.27
1.19
0.27
1.02
0.13
0.07
1.32
0.54
1.07
0.36
1.10
0.17
0.10
2.02
0.59
1.42
0.55
N/A
N/A
0.81
0.47
1.18
0.35
0.85
0.22
0.10
N/A
N/A
1.25
0.48
1.38
0.45
1.13
0.20
0.07
N/A
N/A
1.05
0.43
1.09
0.34
1.04
0.21
0.10
Models: Five fire models were selected for the study, based on the fact that they are commonly
used in fire analyses of NPPs in the U.S. Two of the models consist of simplified engineering
4-2
700
FDS Target Temperature
600
500
400
300
200
100
0
0
100
200
300
400
500
600
700
Consider again Figure 4-1. To make better use of results such as these, two statistical
parameters8 have been calculated for each model and each predicted quantity. The first
parameter, , is the bias factor. It indicates the extent to which the model, on average, underor over-predicts the measurements of a given quantity. For example, the bias factor for the data
shown in Figure 4-1 is 1.02. This means that the model has been shown to slightly
8
The statistical parameters listed in Table 4-1 are based on the versions of the fire models used in the V&V study, circa 2006. As
the models are improved and new validation data introduced, these values may change.
4-4
(4-1)
(4-2)
For some models/quantities, there was an insufficient amount of data to calculate the relative standard deviation of the distribution,
in which case I.D. is shown in the Table.
10
Excel 2007 does not evaluate erfc
for negative values of , even though the function is defined for all real . In such cases,
use the identity erfc
2 erfc .
4-5
ln
(4-3)
ln
ln
(4-4)
For a given model prediction, , the true value of the quantity of interest, is a normally
/ .
distributed random variable with a mean of / and a standard deviation of
4-6
(4-5)
0.14
0.13
0.16
0.14
0.09
0.33
0.40
0.80
0.20
0.14
11
The Kolmogorov-Smirnov test for normality has been applied using the software package SigmaPlot 10, Systat Software, Inc.
The default P value of 0.05 was used.
4-7
3
Plume Temperature, FDTs
Normality Test: Fail
Wall Temperature
Normality Test: Pass
40
30
20
10
0
1
10
10
Interval Number
Interval Number
Figure 4-2. Two examples demonstrating how the validation data is tested for normality.
20 C
280
(4-6)
Step 2: Refer to Table 4-1, which indicates that, on average, FDS over-predicts target
temperatures with a bias factor, , of 1.02. Calculate the adjusted model prediction:
280 C
275 C
1.02
Referring again to Table 4-1, calculate the standard deviation of the distribution:
0.13
280
1.02
36 C
(4-7)
(4-8)
Step 3: Calculate the probability that the actual cable temperature would exceed 330 C12:
12
In the result of Equation 4-9, the precise value for and were used rather than the rounded values shown.
4-8
1
erfc
2
330
1
erfc
2
20
36
275
0.16
(4-9)
The process is shown graphically in Figure 4-3. The area under the bell curve for
temperatures higher than 330 C (625 F) represents the probability that the actual cable
temperature would exceed that value. Note that this estimate is based only on the model
uncertainty.
0.014
0.012
Model Prediction
0.010
0.008
0.006
0.004
P(T>330 C)
0.002
0.000
100
150
200
250
300
350
400
450
500
Temperature (C)
Figure 4-3. Normal distribution of the true value of the cable temperature in a
hypothetical fire.
(4-10)
Step 2: Refer to Table 4-1, which indicates that, on average, MAGIC over-predicts radiant heat
flux with a bias factor, , of 1.07. Calculate the adjusted model prediction:
6
1.07
5.6 kW/m
4-9
(4-11)
6
1.07
1
erfc
2
(4-12)
2.0 kW/m
1
6 kW/m
5.6 kW/m
erfc
2
2.0 kW/m 2
0.42
(4-13)
Even though the model predicts a peak radiant heat flux equal to the critical value, there is only
a 42% chance that the actual heat flux would exceed this value. This is mainly due to the fact
that MAGIC has been shown to overestimate the heat flux by about 7%.
It is important to note that this calculation of model uncertainty does not take into account the
input parameters, such as the HRR of the fire. It is only an assessment of how well the model
MAGIC can estimate the radiant heat flux to a target. Another key point is that the damage
threshold, temperature in the first example and heat flux in the second, is known exactly, i.e.,
without uncertainty.
4-10
13
NUREG/CR-6850 specifies gamma distributions for the various types of combustibles found within an NPP. Microsoft Excel
provides a built-in function (GAMMA.DIST) that calculates the probability density function given the parameters and .
4-11
0.005
HRR Distribution
0.004
0.003
0.002
0.001
0.000
0
200
400
600
For convenience, a spreadsheet can be used to take each value of , from 1 kW to 600 kW,
and compute a corresponding flame height, f (in meters), using Heskestads correlation:
f
0.235
1.02
(4-14)
The diameter of the fire, , is fixed at 0.48 m (1.6 ft), based on the equivalent diameter of the
vent. Whereas Appendix E of NUREG/CR-6850 recommends dividing the range of into 15
bins, it is just as easy for this example to compute the flame height for 600 values of (each
by the
bin has a width of 1 kW). The pdf for the flame height,
f , is related to the pdf for
following expression:
; ,
f
; ,
0.094
(4-15)
This distribution is shown in Figure 4-5. Note that when the derivative in Equation (4-15) is not
easily written in closed form, it is sufficient to calculate the bin width of the model output divided
by that of the model input. In this example, the bin width of the model input parameter, , is
1 kW and the bin width of the model output parameter is the difference in flame heights for two
successive values of .
4-12
1.0
Flame Height Distribution
0.8
0.6
0.4
0.2
0.0
0.0
0.5
1.0
1.5
2.0
2.5
3.0
Figure 4-5. Distribution of flame heights for the entire range of cabinet fires.
Once the pdf for the flame height is calculated, it can be used to determine the probability of the
flames reaching a certain height. In this case, the cable tray is 1.5 m (4.9 ft) above the top of
the cabinet. The probability that the flames from a randomly chosen fire will reach the cables is
given by the area beneath the curve in Figure 4-5 for flame heights greater than 1.5 m (4.9 ft).
In this example, it is approximately 0.31.
Power
(4-16)
The exact value of the Constant is not important; rather, it is the Power that matters. The larger
its absolute value, the more important the Input Parameter. According to the McCaffrey,
Quintiere, and Harkleroad (MQH) correlation, for example, the hot gas layer (HGL) temperature
rise in a compartment fire is proportional to the HRR raised to the two-thirds power:
/
(4-17)
What is important here is the amount that the HGL temperature, , changes due to a shift in
the HRR, . It is the two-thirds power dependence, as found in Table 4-3, that matters. To
4-13
with respect to
2
3
(4-18)
This is a simple formula with which one can readily estimate the relative change in the model
output quantity, /
, due to the relative change in the model input parameter, / .
The uncertainty in a measured quantity is often expressed in relative terms14. Suppose that the
uncertainty in the HRR of the fire, / , is 0.15, or 15%. The expression above indicates that a
15% increase in the HRR should lead to a 2/3 x 15 = 10% increase in the prediction of the HGL
temperature. The result is equally valid for a reduction; if the HRR is reduced by 15%, the HGL
temperature is reduced by 10%.
Table 4-3. Sensitivity of model outputs from Volume 2 of NUREG-1824 (EPRI 1011999).
Important Input
Parameters
Power Dependence
HRR
Surface Area
Wall Conductivity
Ventilation Rate
Door Height
2/3
-1/3
-1/3
-1/3
-1/6
Door Height
HRR
Production Rate
1/2
1
HRR
Soot Yield
1
1
Pressure
HRR
Leakage Rate
Ventilation Rate
2
2
2
Heat Flux
HRR
4/3
Surface/Target
Temperature
HRR
2/3
Output Quantity
HGL Temperature
HGL Depth
Gas Concentration
Smoke Concentration
This relationship is based on experimental data, and has nothing to do with any particular
model; however, an effective way to check a fire model is to take a simple compartment fire
simulation, vary the HRR, and ensure that the change in the HGL temperature agrees with the
correlation. Consider the two curves shown in Figure 4-6. For Benchmark Exercise #3 of the
International Collaborative Fire Model Project (ICFMP), Test 3 was simulated with FDS, using
HRR values of 1000 kW and 1150 kW. An examination of the peak values confirms that the
relative change in the HGL temperature (10%) is two-thirds the relative change in the HRR
14
Note that a differential relationship is only approximate. This method of relating input parameters to output quantities is valid for
relative differences that are less than approximately 30% in absolute value.
4-14
Temperature (C)
200
100
0
0
10
15
20
25
30
Time (min)
Figure 4-6. FDS predictions of HGL Temperature as a function of time due to a 1,000 kW
fire (solid line) and a 1,150 kW fire (dashed).
This section illustrates the usefulness of sensitivity analysis. As an example, consider that
NFPA 805 uses the term Maximum Expected Fire Scenario (MEFS) to describe a severe fire
that could be reasonably anticipated to occur within a compartment and the term Limiting Fire
Scenario (LFS) to describe a severe fire that exceeds one or more performance criteria. The
analyst is often asked to determine the model inputs for both of these scenarios. For MEFS,
input parameters can be chosen from distributions for a particular percentile value. Gallucci
(2011) discusses the issues important in determining the appropriate choice of HRR distribution.
The development of the LFS is essentially a sensitivity analysis performed to identify which
combinations of input parameters or variables are critical to the analysis. The particular
variables to be evaluated depend entirely on the problem being analyzed. At a minimum, the
following parameters should be varied until failure conditions result: HRR, the fire growth rate or
the flame spread rate, the flame radiative fraction or the radiative power, and the location of the
fuel package relative to the target (if variable).
Suppose, for example, that as part of an analysis the problem is to determine the minimum HRR
needed to cause damage in a particular compartment whose HGL temperature is not to exceed
500 C (930 F). The geometrical complexity of the compartment rules out the use of the
algebraic and zone models, and that FDS has been selected for the simulation.
Step 1: Determine a reasonable, but conservative, estimate of what might be the maximum fire
that could occur in the compartment. Using data from NUREG/CR-6850, for this example,
suppose that a 98th percentile HRR for the electrical cabinet fire, 702 kW, has been chosen for
this representative estimate. Choose a model and calculate the peak HGL temperature.
Step 2: The model chosen is FDS and it predicts 450 C (840 F) for the selected fire scenario.
Adjust the prediction to account for the model bias, (see Table 4-1):
4-15
450
20
adj
3
2
adj
500
437
702 kW
437
20
159 kW
(4-20)
This calculation suggests that adding an additional 159 kW to the original 702 kW will produce
an HGL temperature in the vicinity of 500 C (930 F). This result can be double-checked by rerunning the model with the modified input parameters.
Table 4-3 lists several other parameters besides the HRR that can affect the HGL temperature.
Following the example just discussed, similar calculations can be performed in which these
other parameters are varied to determine how else the minimum damage threshold might be
reached. For example, suppose that the surface area, , of the compartment is 400 m2 (4300
ft2). How much would the surface area have to increase or decrease to raise the HGL
temperature to 500 C (930 F)? If the thermal conductivity of the walls, , is 0.1 W/m/K, how
much would it have to change? If the ventilation rate is 1 m3/s (2100 cfm), how much would it
have to change? If the door height, , is 2 m (6.6 ft), how much would it have to change?
Following the example for the HRR, the required changes in these parameters can be
calculated as follows:
3 400 m
adj
adj
(4-21)
500 437
417
0.045 W/m/K
(4-22)
3 1 m s
500 437
417
0.45 m /s
(4-23)
1.8 m
(4-24)
adj
181 m
3 0.1 W/m/K
adj
500 437
417
6 2m
500 437
417
For this example, to increase the HGL temperature by 63 C (145 F), one could increase the
HRR by 159 kW, decrease the surface area of the compartment by 181 m2 (1948 ft2), decrease
the thermal conductivity of the walls by 0.045 W/m/K, decrease the ventilation rate by 0.45 m3/s
(950 cfm), or decrease the door height by 1.8 m (5.9 ft). Of course, some of these options are
not physically possible. Room dimensions and thermal properties are not subject to significant
change, but the HRR and ventilation rates can vary significantly. Also note that if the relative
change in the parameter values exceeds 30 %, it is recommended that further calculations be
performed to confirm the estimated quantity changes.
4-16
4-17
4-18
REFERENCES
10 CFR 50, Voluntary Fire Protection Requirements for Light-Water Reactors, 10 CFR Part 50,
Section 50.48(c), RIN 3150-AG48, Federal Register, Volume 69, Number 115, U.S. Nuclear
Regulatory Commission, Washington, D.C., June 16, 2004.
AISC, Steel Construction Manual, 13th edition, American Institute of Steel Construction, New
York, New York, 2006.
ASME/ANS RA-Sa-2009, Standard for Level 1/Large Early Release Frequency Probabilistic
Risk Assessment for Nuclear Power Plant Applications, Addendum A, ASME/ANS RASa-2009, American Society of Mechanical Engineers, 2009.
ASTM E119-10a, Standard Test Methods for Fire Tests of Building Construction Materials,
American Society for Testing and Materials, West Conshohocken, PA, 2010.
ASTM E135505a (2005), ASTM Standard Guide for Evaluating the Predictive Capability of
Deterministic Fire Models, American Society for Testing and Materials, West
Conshohocken, PA, 2005.
Babrauskas, V., Ignition Handbook, Fire Science Publishers/Society of Fire Protection
Engineers, Issaquah, WA, 2003.
Buchanan, A. H., Structural Design for Fire Safety, John Wiley and Sons, LTD, Chichester,
England, 2001.
Drysdale, D., An Introduction to Fire Dynamics, John Wiley and Sons, 3rd Edition, Chichester, pp.
14, 283-284, 2011.
EPRI 1002981, Fire Modeling Guide for Nuclear Power Plant Applications, Electric Power Research
Institute, Palo Alto, CA, August 2002.
Forney, G. P., Users Guide for Smokeview Version 5 - A Tool for Visualizing Fire Dynamics
Simulation Data, Volume II: Technical Reference Guide, NIST Special Publication 10172, National Institute of Standards and Technology, Gaithersburg, MD, October 2010.
Gallucci, R.H.V., Examination of the Efficacy of the NFPA-805 Fire Modeling Approach
(Comparison between Maximum Expected and Limiting Fire Scenarios), ANS PSA
2011 International Topical Meeting on Probabilistic Safety Assessment and Analysis,
Wilmington, NC, March 13-17, 2011, American Nuclear Society, LaGrange Park, IL,
2011.
Gay, L., User Guide of MAGIC Software V4.1.1, EdF HI82/04/022/B, Electricit de France, France,
April 2005.
Gay, L., C. Epiard, and B. Gautier, MAGIC Software Version 4.1.1: Mathematical Model, EdF
HI82/04/024/B, Electricit de France, France, November 2005.
Gay, L., J. Frezabeu, and B. Gautier, Qualification File of Fire Code MAGIC Software Version 4.1.1:
Mathematical Model, EdF HI82/04/024/B, Electricit de France, France, November 2005.
5-1
REFERENCES
Grove, B. S., and J.G. Quintiere, Calculating Entrainment and Flame Height in Fire Plumes of
Axisymmetric and Infinite Line Geometries, Journal of Fire Protection Engineering, Volume
12, August, 2002, pp. 117 - 137.
Gottuk, D. T., and Lattimer, B. Y., Section 2, Chapter 2-5, Effect of Combustion Conditions on
Species Production, SFPE Handbook of Fire Protection Engineering, 4th Edition (P. J.
DiNenno, Editor-in-Chief), National Fire Protection Association and The Society of Fire
Protection Engineers, Quincy, MA, 2008.
Heskestad, G., Section 2, Chapter 2-1, Fire Plumes, Flame Height, and Air Entrainment, SFPE
Handbook of Fire Protection Engineering, 4th Edition (P. J. DiNenno, Editor-in-Chief),
National Fire Protection Association and The Society of Fire Protection Engineers, Quincy,
MA, 2008.
Holman, J. P., Heat Transfer, 7th edition, McGraw-Hill, New York, 1990.
Janssens, M., Section 3, Chapter 3-2, Calorimetry, SFPE Handbook of Fire Protection
Engineering, 4th Edition (P. J. DiNenno, Editor-in-Chief), National Fire Protection Association
and The Society of Fire Protection Engineers, Quincy, MA, 2008.
Jones, W., R. Peacock, G. Forney, and P. Reneke, CFAST: An Engineering Tool for Estimating Fire
Growth and Smoke Transport, Version 5 - Technical Reference Guide, NIST Special
Publication 1030, National Institute of Standards and Technology, Gaithersburg, MD, 2004.
Karlsson, B., and J. Quintiere, Enclosure Fire Dynamics, CRC Press, Boca Raton, Florida, 2000.
Kylu, U. O., and G.M. Faeth, Carbon Monoxide and Soot Emissions from Liquid-Fueled Buoyant
Turbulent Diffusion Flames, Combustion and Flame, 87:61-76, 1991.
McCaffrey, B. J., J. G. Quintiere, and M. F. Harkleroad, Estimating Compartment Temperature and
Likelihood of Flashover Using Fire Test Data Correlation, Fire Technology, Volume 17, No.
2, pp. 98-119, Quincy, MA, 1981.
McGrattan, K. et al., Fire Dynamics Simulator (Version 5) Technical Reference, Volume 3:
Validation, NIST Special Publication 1018-5, National Institute of Standards and
Technology, Gaithersburg, MD, 2010.
McGrattan, K., B. Klein, S. Hostikka, and J. Floyd, Fire Dynamics Simulator (Version 5) Users
Guide, NIST Special Publication 1019-5, National Institute of Standards and Technology,
Gaithersburg, MD, 2009.
McGrattan, K., and B. Toman, Quantifying the predictive uncertainty of complex numerical models,
Metrologia, 48:173-180, 2011.
Mulholland, G. W., and C. Croarkin, Specific Extinction Coefficient of Flame-Generated Smoke,
Fire and Materials, 24:227-230, 2000.
NFPA, Fire Protection Handbook, National Fire Protection Association, 20th Ed., A. E. Cote
(Editor), 2008.
NFPA 70 (NEC 2008), National Electric Code, National Fire Protection Association, Quincy, MA,
2008.
NFPA 805, Performance-Based Standard for Fire Protection for Light Water Reactor Electric
Generating Plants, National Fire Protection Association, Quincy, MA, 2001.
NEI 00-01 (2009), Guidance for Post Fire Safe Shutdown Circuit Analysis, Revision 2, Nuclear
Energy Institute, Washington, D.C., May, 2009.
5-2
REFERENCES
NEI 00-01 (2010), Guidance for Post Fire Safe Shutdown Circuit Analysis, Draft Revision 3,
Nuclear Energy Institute, Washington, D.C., 2010.
NEI 04-02 (2009), Guidance for Implementing a Risk-Informed Performance-Based Fire
Protection Program Under 10 CFR 50.48(c), Rev. 1, Nuclear Energy Institute,
Washington, D.C., September, 2009.
NIST NCSTAR 1-5F, Federal Building and Fire Safety Investigation of the World Trade Center
Disaster: Computer Simulation of the Fires in the World Trade Center Towers, National
Institute of Standards and Technology, Gaithersburg, MD, 2005.
NRC, Regulatory Guide 1.189, Fire Protection for Nuclear Power Plants, U.S. Nuclear
Regulatory Commission, Bethesda, MD, April, 2009.
NUREG-1805, Fire Dynamics Tools (FDTs): Quantitative Fire Hazard Analysis Methods for the U.S.
Nuclear Regulatory Commission Fire Protection Inspection Program, U.S. Nuclear
Regulatory Commission, Washington, D.C., December 2004.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 1: Main Report, U.S. Nuclear Regulatory Commission, Office of
Nuclear Regulatory Research (RES), Washington, D.C., 2007, and EPRI 1011999, Electric
Power Research Institute (EPRI), Palo Alto, CA.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 3: Fire Dynamics Tools (FDTs), U.S. Nuclear Regulatory Commission,
Office of Nuclear Regulatory Research (RES), Washington, D.C., 2007, and EPRI 1011999,
Electric Power Research Institute (EPRI), Palo Alto, CA.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 4: Fire-Induced Vulnerability Evaluation (FIVE-Rev1), U.S. Nuclear
Regulatory Commission, Office of Nuclear Regulatory Research (RES), Washington, D.C.,
2007, and EPRI 1011999, Electric Power Research Institute (EPRI), Palo Alto, CA.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 5: Consolidated Fire Growth and Smoke Transport Model (CFAST),
U.S. Nuclear Regulatory Commission, Office of Nuclear Regulatory Research (RES),
Washington, D.C., 2007, and EPRI 1011999, Electric Power Research Institute (EPRI), Palo
Alto, CA.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 6: MAGIC, U.S. Nuclear Regulatory Commission, Office of Nuclear
Regulatory Research (RES), Washington, D.C., 2007, and EPRI 1011999, Electric Power
Research Institute (EPRI), Palo Alto, CA.
NUREG-1824, Verification and Validation of Selected Fire Models for Nuclear Power Plant
Applications, Volume 7: Fire Dynamics Simulator (FDS), U.S. Nuclear Regulatory
Commission, Office of Nuclear Regulatory Research (RES), Washington, D.C., 2007, and
EPRI 1011999, Electric Power Research Institute (EPRI), Palo Alto, CA.
NUREG-1855 Volume 1, Guidance on the Treatment of Uncertainties Associated with PRAs in
Risk-Informed Decision Making, U.S. Nuclear Regulatory Commission, Office of Nuclear
Regulatory Research (RES), Washington, D.C., 2009.
NUREG/CR-4680, Heat and Mass Release Rate for Some Transient Fuel Source Fires: A Test
Report, U.S. Nuclear Regulatory Commission, Washington, D.C., December 2004.
5-3
REFERENCES
NUREG/CR-6738, Risk Methods Insights Gained from Fire Incidents, U.S. Nuclear Regulatory
Commission, Office of Nuclear Regulatory Research (RES), Washington, D.C., 2001.
NUREG/CR-6850, EPRI/NRC-RES Fire PRA Methodology for Nuclear Power Facilities: Volume 1:
Summary and Overview, U.S. Nuclear Regulatory Commission, Office of Nuclear
Regulatory Research (RES), Washington, D.C., 2005, and EPRI 1011989, Electric Power
Research Institute (EPRI), Palo Alto, CA.
NUREG/CR-6850, EPRI/NRC-RES Fire PRA Methodology for Nuclear Power Facilities: Volume 2:
Detailed Methodology, U.S. Nuclear Regulatory Commission, Office of Nuclear Regulatory
Research (RES), Washington, D.C., 2005, and EPRI 1011989, Electric Power Research
Institute (EPRI), Palo Alto, CA.
NUREG/CR-6850 Supplement 1, Fire Probabilistic Risk Assessment Methods Enhancements, U.S.
Nuclear Regulatory Commission, Office of Nuclear Regulatory Research (RES),
Washington, D.C., 2010, and EPRI 1019259, Electric Power Research Institute (EPRI), Palo
Alto, CA.
NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 1: Test Descriptions and
Analysis of Circuit Response Data, U.S. Nuclear Regulatory Commission, Washington,
D.C., 2007.
NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 2: Cable Fire Response
Data for Fire Model Improvement, U.S. Nuclear Regulatory Commission, Washington, D.C.,
2007.
NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 3: Thermally-Induced
Electrical Failure (THIEF) Model, U.S. Nuclear Regulatory Commission, Washington, D.C.,
2007.
NUREG/CR-7010, Cable Heat Release, Ignition, and Spread in Tray Installations during Fire
(CHRISTIFIRE), Phase 1: Horizontal Trays, National Institute of Standards and
Technology, Gaithersburg, MD, 2012.
Peacock, R., K. McGrattan, B. Klein, W. Jones, and P. Reneke, CFAST Consolidated Model
of Fire Growth and Smoke Transport (Version 6) Software Development and Model
Evaluation Guide, Special Publication 1086, National Institute of Standards and
Technology, Gaithersburg, MD, 2008.
Peacock, R., W. Jones, P. Reneke, and G. Forney, CFAST: An Engineering Tool for Estimating
Fire Growth and Smoke Transport, Version 6 Users Guide, Special Publication 1041,
National Institute of Standards and Technology, Gaithersburg, MD, 2008.
Quintiere, J. G., Principles of Fire Behavior, Delmar Publishers, 1998.
Quintiere, J. G., Fundamentals of Fire Phenomena, John Wiley, 2006.
Salley, M.H., J. Dreisbach, K. Hill, R. Kassawara, B. Najafi, F. Joglar, A. Hamins, K. McGrattan,
R. Peacock; and B. Gautier, Verification and Validation--How to Determine the
Accuracy of Fire Models, Fire Protection Engineering, vol. 34, pp. 34 - 44, Spring 2007.
Schiffiliti, R., and W. Pucci, Fire Modeling, State of the Art, Fire Detection Institute, 1996.
SFPE, SFPE Engineering Guide to Assessing Flame Radiation to External Targets from Pool
Fires, SFPE Engineering Guide, Society of Fire Protection Engineers, Bethesda, MD,
March 1999.
5-4
REFERENCES
SFPE, SFPE Engineering Guide to Fire Exposures to Structural Elements, SFPE Engineering
Guide, Society of Fire Protection Engineers, Bethesda, MD, November, 2005.
SFPE, SFPE Engineering Guide to Piloted Ignition of Solid Materials Under Radiant Exposure,
SFPE Engineering Guide, Society of Fire Protection Engineers, Bethesda, MD, January,
2002.
SFPE, SFPE Engineering Guide to Predicting Room of Origin Fire Hazards, SFPE Engineering
Guide, Society of Fire Protection Engineers, Bethesda, MD, November, 2007.
SFPE, SFPE Engineering Guide to Substantiating a Fire Model for a Given Application, SFPE
Engineering Guide, Society of Fire Protection Engineers, Bethesda, MD, June, 2010.
SFPE, SFPE Handbook of Fire Protection Engineering, 4th Edition (P. J. DiNenno, Editor-in-Chief),
National Fire Protection Association and The Society of Fire Protection Engineers, Quincy,
MA, 2008.
SFPE, SFPE Engineering Standard on Calculating Fire Exposures to Structures, SFPE S.01.2001,
Society of Fire Protection Engineers, Bethesda, MD, 2011.
Underwriters Laboratories, Single Station Fire Alarm Device, UL 217, Underwriters Laboratories,
Northbrook, Illinois.
Wallis, G. B., Draft Final NUREG-1824, Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, Memorandum from Advisory Committee on Reactor
Safeguards to Luis Reyes, ML062980154, U.S. Nuclear Regulatory Commission, October,
2006.
Yuan, L., and G. Cox, An Experimental Study of Some Line Fires, Fire Safety Journal, Volume 27,
1996, pp. 123 139.
5-5
REFERENCES
5-6
Chapter 3
Generic Fire
Scenario
3.2.5
3.2.1, 3.2.2
3.2.2
3.2.4
3.2.1
3.2.7
3.2.3
3.2.1, 3.2.6
The structure of each appendix reflects the analysis sequence as well as the recommended
means by which to document each type of fire modeling analysis. In particular, the appendices
use the following format
Fire modeling objective. This section clearly presents the issue that is being evaluated
and the reason the analysis is being performed.
Description of the fire scenario. This section provides a detailed description of the fire
scenario being evaluated, including important parameters such as the room geometry,
6-1
Selection of the fire models. This section presents the basis for selecting the different
types of fire models used to address the fire modeling objective. This includes
considerations of the model capability as well as the available verification and validation
(V&V) basis. Generally, two or more models are selected in order to highlight the
strengths and weaknesses of the models for the given application. The selected models
are not necessarily the only models that can be used, but the basis provided is broadly
applicable to a given class of fire models.
Estimation of the fire-generated conditions. This section presents model specific inputs
and supporting calculations and identifies applicable model nuances for the scenario
considered.
Evaluation of results. The results of the fire models are presented in terms of the data
necessary to support the fire modeling objective in this section. Model sensitivity and
uncertainty are also described for each output parameter and model considered.
Conclusions. A summary of the fire modeling analysis, including the outcome relative to
the fire modeling goal is presented in this section.
The authors strove to create fire scenarios based on actual commercial NPP configurations in
order to provide realistic fire scenario development and results. Be cautioned that the fire
modeling results are not generic, and different applications will require careful consideration of
the various input parameters, the fire model applicability, and the dominant exposure
mechanism(s). However, the structure presented is generic, and when followed it is expected to
lead to a fire modeling analysis that could be used to support safety related applications.
This report also contains a CD. On the CD you will find the complete NUREG-1934 final report
in pdf format and the input files used for the fire model runs discussed in each appendix. In
addition, the current versions of the FDTs, CFAST and FDS that were used for each appendix
are on the CD. The complete layout and contents of the CD are described in the README file
included on the CD.
6-2
APPENDIX A
A-1
A-2
A-3
A-4
600
500
400
300
200
100
0
0
600
1200
1800 2400
Time (s)
3000
3600
Figure A-5. Time history of the HRR used by all models in the MCR scenario.
Based on a physical assessment of the cabinet, it is determined that the exterior panels of the
burning cabinet do not open before or during the fire. The smoke, heat, and flames are
exhausted from an air vent in the side of the cabinet. The top of the air vent is 0.3 m (1 ft) below
the top of the cabinet. The air vent is 0.6 m (2 ft) wide and 0.2 m (8 in) high. The cabinet is
2.4 m (8 ft) tall.
The heat of combustion and product yields for XPE/neoprene cable are taken from Table 3-4.16
of the SFPE Handbook, 4th edition and are listed in Table A-1. When estimating the
composition of the fires combustion products, the jacket and insulation material of the cable are
taken as an equal-parts mixture of polyethylene (C2H4) and neoprene (C4H5Cl), with the
effective chemical formula C3H4.5Cl0.5.
Table A-1. Data for MCR fire based on XPE/neoprene electrical cable.
Parameter
Effective Fuel Formula
Peak HRR
Time to reach peak HRR
Heat of Combustion
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Mass Extinction Coefficient
Value
C3H4.5Cl0.5
702 kW
720 s
10,300 kJ/kg
0.63 kg/kg
0.175 kg/kg
0.082 kg/kg
0.53
8700 m2/kg
Source
Combination of polyethylene and neoprene
NUREG/CR-6850 (EPRI 1011989), App. G
NUREG/CR-6850 (EPRI 1011989), App. G
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
Mulholland and Croarkin (2000)
A-5
The original edition of NUREG/CR-6850 (EPRI 1011989) contains an error in the specification of the optical density (NRC ADAMS
-1
Accession Number ML061630360). The value of 3 m is correct.
A-6
A.3.4 Validation
The principal source of validation data justifying the use of the fire models discussed above for
this scenario is the U.S. Nuclear Regulatory Commission/Electric Power Research Institute
(NRC/EPRI) verification and validation (V&V) study documented in NUREG-1824 (EPRI
1011999). The National Institute of Standards and Technology (NIST) has expanded the
NRC/EPRI V&V to include the latest versions of CFAST (6.1.1) (Peacock, 2008) and FDS
(5.5.3) (McGrattan, 2010). In particular, a Factory Mutual/Sandia National Laboratories
(FM/SNL) test series was designed specifically as a mock-up of a control room in a nuclear
power plant (NPP). One of these experiments (Test 21) involves a fire within a hollow steel
cabinet.
Table A-2 lists the important non-dimensional parameters that characterize the fire scenario and
the ranges for which the NRC/EPRI validation study is applicable. A few parameters fall outside
the validation parameter space and are addressed individually:
The Fire Froude Number falls outside the range. This parameter is essentially a measure of
the fires heat output relative to its base area. In this example, the fire is specified as
emanating from the side of the cabinet, with the vent opening serving as its base. This
leads to a higher value of
than would be calculated if the fire were burning completely
outside of the cabinet. Thus, the high value of
is the represents more severe fire
conditions than would be expected if the fire were specified as burning partially within the
cabinet.
The relatively low Equivalence Ratio for the compartment is a result of the relatively large
amount of air forced into the room during the smoke purge mode. Twenty-five air changes
per hour is a considerable flow rate, and no validation experiment in NUREG-1824 (EPRI
1011999) involved such a high ventilation rate. However, the results of all the model
simulations indicate that the scenario in which the ventilation is turned off is most likely to
compromise human habitability, and the presence of any level of ventilation reduces room
temperature and heat flux and increases visibility.
For the scenario with no ventilation, the classic definition of the Equivalence Ratio does not
apply because there is no supply of oxygen in the room. However, it can be shown that
there is sufficient oxygen in the room to sustain the specified fire. The total mass of oxygen
in the room is the product of the density of air, , the volume of the room, , and the mass
fraction of oxygen in the air, O :
O ,tot
1.2 kg/m
1945 m
A-7
0.23
537 kg
(A-1)
12
8
3
13,100 kJ/kg
60 s/min
19
min
2
69 kg
(A-2)
These calculations show that the quantity of oxygen in the room would be able to sustain the
specified cabinet fire.
The ratio of the Target Distance relative to the Fire Diameter, / , exceeds the range of the
validation study. However, this parameter is only relevant to the point source radiation heat
flux calculation, which is by definition more accurate, as the target moves further from the
source. Thus, although the parameter is outside the validation range, it is not outside of the
methodologys range of validity.
A-8
Quantity
Fire Froude
Number
702 kW
1.2 kg/m
Flame
Height,
,
relative to the
Ceiling
Height,
Ceiling Jet
Radial
Distance, cj ,
relative to the
Ceiling
Height,
Equivalence
Ratio, , of
the room,
based on
Forced
Ventilation of
Purge Mode
Compartment
Aspect Ratio
Target
Distance, ,
relative to the
Fire
Diameter,
3.7
1.02
0.9
0.4 m 3.7
6.2
1.02
702 kW
13,100 kJ/kg 3.7 kg/s
0.23
24.6 m
5.2 m
4.7
1.2 kg/m
13.4 m /s
22
No
0.2 1.0
Yes
1.2 1.7
N/A
0.04 0.6
No
0.6 5.7
Yes
2.2 5.7
No
0.014
3.7 kg/s
16.2 m
5.2 m
8.8 m
0.4 m
0.4 2.4
2.7 m
In
Range?
6.2
9.8 m/s
2.1 m 2.7 m
5.2 m
NUREG1824
Validation
Range
3.1
Notes:
(1)
(2)
A-9
(A-3)
The parameters for the FIVE-Rev1 implementation of the FPA are listed in Table A-3. Note that
the specified time-dependent HRR is used in the calculation, but that the calculation does not
include either the fires elevation above the floor or any other information about the fire. The
walls, ceiling, and floor are all specified as gypsum board rather than concrete because the FPA
correlation only accounts for one type of lining material. Gypsum board was chosen because it
has a lower thermal conductivity than concrete, which results in a slightly higher HGL
temperature.
Table A-3. Summary of input parameters for the FPA calculation of the MCR.
Parameter
Room height (H)
Room effective length (Le)
Room effective width (We)
Room boundary material
Mech. ventilation rate ( )
Ambient temperature (Ta)
Fire parameters
Value
5.2 m
27.1 m
13.8 m
Gypsum board
13.4 m3/s
20 C
A-10
Source
Figure A-1
Equation (A-3)
Equation (A-3)
Table 3-1
Specified (25 ACH)
Specified
Table A-1
Uniform
compartment
temperature
Pointsourcefirewithtime
dependentheatreleaserate
butnospecifiedlocationor
height
Closeddoor
Figure A-6. Schematic diagram of the FPA calculation for the MCR smoke purge scenario.
Heat Flux
The point source model is used to estimate the heat flux from the flames to the operator when
the fire is at its peak HRR. The peak HRR, , is 702 kW, the radiative fraction, , is 0.53, and
the distance from the cabinet vent to the operator is approximately 8.8 m (29 ft). The heat flux
is calculated:
0.53
4
702 kW
8.8 m
0.38 kW/m
(A-4)
While this heat flux prediction is well below the critical value of 1 kW/m2, it does not account for
the thermal radiation from the HGL. Thus, the point source method can be used as a screening
tool, and further analysis can be performed by CFAST and FDS.
Smoke concentration and visibility
Neither the FDTs nor FIVE-Rev1 include methods to calculate smoke concentrations or visibility
in mechanically ventilated enclosure fires, but calculation methods provided in Section 3,
Chapter 9, of the SFPE Handbook, 4th ed. are relatively simple to apply and are based on the
same principles and concepts embodied in zone models. These hand calculations provide an
estimate of the fire-generated smoke concentrations and visibility conditions for this scenario
and will indicate whether more detailed modeling is warranted.
The soot mass generation rate,
, is the product of the soot yield, , and the mass burning
rate of fuel,
. The latter quantity is obtained by dividing the HRR, , by the heat of
combustion, :
0.175
702 kW
10,300 kJ/kg
A-11
0.012 kg/s
(A-5)
tot
, is then calculated:
0.012 kg/s
1.2 kg/m
13.4 m s
0.00075 kg/kg
(A-6)
, is calculated:
1.2 kgm
0.00075 kg/kg
7.8 m
(A-7)
Here
is the mass-specific extinction coefficient listed in Table A-1. By definition, the optical
density of the smoke is related to the extinction coefficient via the expression:
ln 10
7.8 m
2.3
3.4 m
(A-8)
This calculated optical density is then compared with the tenability limit, 3 m-1. This calculation
indicates that the visibility criterion for tenability would be exceeded for the MCR smoke purge
scenario based on the specified parameters and the MCR uniformly filling with smoke.
However, because this analysis is based only on the peak rather than the time-dependent HRR,
further modeling with zone and CFD models is warranted.
A-12
Figure A-7. Compartment geometry and surface material selection in CFAST for the MCR
Fire.
While there are numerous cabinets and tables in the compartment, they do not significantly
change the overall volume of the room and may be neglected as obstructions. There are no
mechanisms within CFAST to account for the open-grate ceiling. This ceiling provides a
negligible flow resistance to the heat and air that go through it, so it can be ignored. It is
expected that neglecting it will lead to slightly higher HGL temperatures because there is less
resistance for the rising smoke and hot gases.
Fire: In CFAST, a fire is described as a source of heat placed at a specific point within a
compartment that generates combustion products according to user-specified combustion
chemistry. Consistent with typical practice for the use of zone fire models for electrical cabinet
fires, the fire is positioned at the top of the air vent, 0.3 m (1 ft) below the top of the cabinet, at
A-13
Note: Values for Total Mass, Heat of Gasification, and Volatilization Temperature are set at default values.
Figure A-8. Specification of the fire in CFAST for the MCR Fire.
Materials: CFAST does not include the ability to model individual walls of different materials,
but it can model different materials for walls, floors and ceilings. For this example, the
compartment walls are specified as being entirely made of gypsum wallboard which has a
thermal conductivity lower than the other wall materials. The floor and ceiling are modeled as
0.5 m (1.6 ft) thick concrete.
Ventilation: For the smoke-purge calculation, air is supplied to the MCR via the six supply
vents and exhausted through the two returns. The total ventilation rate is 25 air changes per
hour, or 13.4 m3/s. Mechanical ventilation inputs for CFAST are shown in Figure A-9. A
snapshot of the CFAST simulation is shown in Figure A-10.
A-14
Note: Values for Begin Dropoff At and Zero Flow At are set at default values.
Figure A-9. Mechanical ventilation inputs in CFAST for the MCR fire.
Figure A-10. Smokeview rendering of the CFAST simulation of the MCR fire with
mechanical ventilation.
A-15
200 kg/m
2.0 kJ/kg/K
0.68 kW/m /K
The
used to estimate the individual thermal parameters is derived from an empirical ignition
model, and it may differ somewhat significantly from the true value. However, in this case, the
heat losses to the carpet are not expected to play a significant role in the simulation and the
choice of its thermal properties is not expected to affect the results. If the boundary is expected
to be an important heat loss surface, then a more accurate determination of the thermal
parameters is necessary.
Ventilation: Air is supplied to the MCR via the six supply vents and exhausted through the two
returns. The supply rate is divided equally among the six supply vents, and the return rate is
divided equally among the two returns. The leakage from the compartment is modeled by
specifying a small vent located at the base of the door through which air escapes at a rate
determined by the pressure difference between the MCR and ambient. Note that the door crack
itself is not modeled explicitly, as the numerical grid is not fine enough. Rather, the leak is
spread over a slightly larger area. The volume flow through the leakage area, , is estimated
via the equation:
A-16
A-17
Model
Bias
Factor,
Standard
Deviation,
Ventilation
Critical
Value
Probability
of
Exceeding
70
95
0.000
61
48
82
70
95
95
95
95
0.000
0.000
0.009
0.000
0.4
0.1
0.2
0.6
0.4
1
1
1
1
1
0.000
0.000
0.000
0.228
0.000
7.6
0.5
54
31
3
3
3
3
0.471
0.000
0.912
0.909
Predicted
Value
1.56
0.32
1.06
1.03
1.06
1.03
0.12
0.07
0.12
0.07
Purge
No Vent.
Heat Flux (kW/m2)
FIVE-Rev1
CFAST
FDS
CFAST
FDS
1.42
0.81
0.85
0.81
0.85
0.55
0.47
0.22
0.47
0.22
Purge
No Vent.
Optical Density (m-1)
CFAST
FDS
CFAST
FDS
2.65
2.7
2.65
2.7
0.63
0.55
0.63
0.55
Purge
No Vent.
A-18
HGL Height
6
FPA (Purge)
120
100
CFAST (Purge)
FDS (No Vent.)
80
Height (m)
Temperature (C)
140
FDS (Purge)
60
40
4
3
2
1
20
0
0
0
600
1200
1800
2400
3000
3600
600
1200
1800
2400
3000
3600
Time (s)
Time (s)
Figure A-12. Hot Gas Layer Temperature and Height for the MCR scenario.
The bias factor for the CFAST predictions of HGL temperature is 1.06. Taking this into account, the CFAST prediction of HGL
temperature for the 35% radiative fraction case would be 90 C (194 F) rather than 95 C (203 F). This is close enough to the
critical value to warrant further modeling.
A-19
1.0
FIVE (Purge)
CFAST (No Vent.)
0.8
CFAST (Purge)
FDS (No Vent.)
0.6
FDS (Purge)
0.4
0.2
0.0
0
600
1200
1800
2400
3000
3600
Time (s)
A-20
60.0
SFPE (Purge)
CFAST (No Vent.)
50.0
CFAST (Purge)
FDS (No Vent.)
40.0
FDS (Purge)
30.0
20.0
10.0
0.0
0
600
1200
1800
2400
3000
3600
Tim e (s)
A.6 Conclusion
A fire modeling analysis has been performed to assess the habitability of the MCR in the event
of a fire within an isolated electrical cabinet. The fire is not expected to spread to other
cabinets. Of the three MCR abandonment criteria, it is most likely that the operators would be
forced to abandon the MCR because the optical density would surpass 3 m-1 approximately 12
minutes after the fire ignites if the smoke purge system is not activated before this time,
according to the FDS analysis. A simple analytical method and the zone model CFAST indicate
that the optical density would exceed the critical value with the smoke purge system on and with
the ventilation system turned off. However, these analyses are based on the use of several
important conservative parameters. For the smoke purge case, the analytical method predicts
that the smoke fills the entire compartment uniformly, even though the FDS analysis shows that
the supply vents maintain visibility in the vicinity of the operator location. CFAST reports the
optical density of the upper layer, but does not predict that the upper layer would descend to the
level of the operator in either the purge or no-ventilation scenario based on the conservative
specifications, at least for a fire having a base height of 2 m (6.6 ft).
A-21
A.7 References
1. Mulholland, G.W., and C. Croarkin, Specific Extinction Coefficient of Flame Generated
Smoke, Fire and Materials, 24:227230, 2000.
2. NIST SP 1018-5, Fire Dynamics Simulator (Version 5), Technical Reference Guide,
Volume 3, Experimental Validation, 2010.
3. NIST SP 1030. CFAST: An Engineering Tool for Estimating Fire Growth and Smoke
Transport, Version 5 - Technical Reference Guide, National Institute of Standards and
Technology, Gaithersburg, Maryland, 2004.
4. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
5. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for Nuclear
Power Plant Applications, 2007.
6. NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
7. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
A-22
APPENDIX B
The cable interior temperature is taken to be that of the insulation material surrounding the conductors.
B-1
In this scenario, the fire propagates upward through the array of cable trays according to an
empirically determined timing sequence. First, ignition of the cables in the lowest tray occurs
when the internal temperature of a target cable within that tray reaches the failure temperature
of 205 C (400 F). This analysis is based on guidance provided in NUREG/CR-6850 (EPRI
1011989), Appendix R. The calculation of the cables internal temperature is based on the
Thermally-Induced Electrical Failure (THIEF) methodology (NUREG/CR-6931, Volume 3).
Following ignition, the cables in the lowest tray burn at a rate of 250 kW/m2, a value appropriate
for TP cables (NUREG/CR-7010, Volume 1). The width of the burning cable is the same as the
B-2
tan 35
(B-1)
where is the length of tray and is the distance (bottom to bottom) between tray
1 and
tray . The 35 upward spread angle is described in NUREG/CR-6850 (EPRI 1011989),
is equal to the length of the vent in the cabinet (0.6 m (2 ft))
Appendix R. In this example,
and is the distance between trays (bottom to bottom).
Following ignition, the fire in the first tray spreads laterally at a rate of 3.2 m/h (NUREG/CR6850 (EPRI 1011989), Appendix R). The fire in the second tray ignites 4 min after the first, and
the lateral extent of the initial fire in the second tray is widened based on Equation (B-1) to
1.3 m. The burning and spread rates of the fire in the second tray are the same as the first.
The fire in the third tray ignites 3 min after the fire in the second, and the initial lateral extent of
the fire is widened yet again following the 35 spread angle to 2.0 m.
Local burnout of the fire occurs when the cable plastic is consumed. The time to burnout is
calculated as follows. First, determine the combustible mass per unit area of tray,
:
1
50
0.32
1 0
0.8 m
0.38 kg/m
7.6 kg/m
(B-2)
is the mass fraction of combustible (i.e., nonwhere is the number of cables per tray,
is the total mass per unit
metallic or plastic) material in the cable, is the residue yield,
length of a single cable, and
is the tray width. Next, calculate the burnout time, :
5 avg /6
7.6 kg/m
20,900 kJ/kg
5/6 250 kW/m
762 s
(B-3)
where is the heat of combustion18, and avg is the average HRR per unit area of tray. The
FLASH-CAT model asserts that the HRR per unit area ramps linearly to its average value over a
time period of /6, remains steady for a time period of 2 /3, and then decreases linearly to
zero over a time period of /6. The linear ramp-up and ramp-down are typical ways of
approximating the time history of an items HRR. Further details of the FLASH-CAT model are
provided in NUREG/CR-7010, Volume 1.
The heat of combustion and product yields for PE/PVC cables are taken from Table 3-4.16 of
the SFPE Handbook, 4th edition. Note that five different types of PE/PVC cables are listed in
the chapter. The values listed in Table B-1 are for a cable with relatively high soot and CO
yields, typical of an under-ventilated fire burning within a closed cabinet. Note also that the nonmetallic components of the cables are a mixture of PE (C2H4) and PVC (C2H3Cl). Because the
mixture consists of approximately the same mass of each, the cable materials with an effective
chemical formula of C2H3.5Cl0.5 have been selected. Table B-1 summarizes the fuel and
combustion parameters for this scenario.
18
By default, the FLASH-CAT model uses a heat of combustion of 16,000 kJ/kg unless there is experimental data that is more
appropriate. In this case, the chosen heat of combustion is based on a measurement of PE/PVC cable.
B-3
Parameter
Effective Fuel Formula
Peak HRR
Heat of Combustion
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Value
C2H3.5Cl0.5
464 kW
20,900 kJ/kg
1.29 kg/kg
0.136 kg/kg
0.147 kg/kg
0.49
Source
Combination of polyethylene and PVC
NUREG/CR-6850 (EPRI 1011989), App. G
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
B-4
B-5
B.3.3 Validation
NUREG-1824 (EPRI 1011999) contains experimental validation results for CFAST and FDS that
are appropriate for this scenario. The National Institute of Standards and Technology (NIST)
has expanded the U.S. Nuclear Regulatory Commission/Electric Power Research Institute
B-6
B-7
Quantity
NUREG1824
Validation
Range
In
Range?
0.4 2.4
No
0.2 1.0
Yes
1.2 1.7
N/A
0.04 0.6
Yes
0.6 5.7
Yes
2.2 5.7
Yes
464 kW
1.2 kg/m
Fire Froude
Number
2.6
9.8 m/s
1600 kW
Flame Length, Lf ,
relative to the
Ceiling Height, Hc
Ceiling Jet Radial
Distance,rcj ,
relative to the
Ceiling Height, Hc
Equivalence Ratio,
, as an indicator of
the Ventilation Rate
Compartment
Aspect Ratio
Target Distance, r,
relative to the Fire
Diameter, D
3.7
1.4
1.2 kg/m
1.02
0.48 m 3.7
2.6
1.02
2.1 m
0.23
kW
kJ/kg
. kg/s
0.23
26.5 m
6.1 m
1.2 kg/m
1.4 m /s
18.5 m
6.1 m
4.3
1.5 m
0.48 m
3.1
0.4 kg/s
3.0
Notes:
(1)
(2)
B-8
Cable
Tray C
Cable
Tray B
Cable
Tray A
Bank A
Bank B
Bank C
The first step in using algebraic models is to determine whether the cables in the cable trays
located directly above the cabinet fire are likely to be damaged and potentially ignited by the
cabinet fire. As shown in Figure B-1, the top of the cabinet is located at an elevation of 2.4 m
(7.9 ft), and the lowest overhead cable tray is located at an elevation of 3.9 m (12.8 ft), which is
1.5 m (4.9 ft) above the top of the cabinet. As shown in Table B-2, the flame length of the
cabinet fire is calculated to be 2.1 m (6.9 ft) at the cabinet peak heat release rate of 464 kW, so
this empirical correlation for flame length can be used to confirm that the overhead cables would
likely be damaged in this scenario. This calculation supports the scenario that the overhead
cable trays would be ignited by the cabinet fire, as discussed in Section B.2.
The next step is to calculate the fire plume temperatures that develop from the cabinet fire to
determine which of the three cable trays located above the cabinet would be damaged by the
cabinet fire. The Heskestad plume temperature correlation included in the FDTs and the FireInduced Vulnerability Evaluation (FIVE-Rev1) was used to calculate the plume centerline
temperature above the cabinet fire. The results of this calculation are shown in Figure B-3.
These results show that the plume temperature at all three cable trays would exceed the cable
damage threshold temperature of 205 C (400 F). However, the Heskestad plume temperature
B-9
Plumetemperatures Heskestadcorrelation
1000
900
Temperature(C)
800
700
600
Bottom
Middle
Top
500
400
300
200
100
0
0
600
1200
1800
Time(s)
2400
3000
3600
Figure B-3. Plume temperatures at cable trays located above a cabinet fire.
The FPA forced ventilation correlation is used to estimate the average HGL temperature of the
SWGR resulting only from the cabinet fire only, based on the parameters described in the
following subsections and summarized in Table B-3.
Geometry: The FPA correlation requires room dimensions to be specified in terms of length,
width, and height. For this example, the selected compartment is a rectangular parallelepiped,
so its length, width, and height are specified directly from dimensions shown in Figure B-1.
Fire: As applied to this scenario, the FPA correlation is used with the time-dependent HRR
specified for the cabinet fire only. This HRR history is shown in Figure B-4.
Materials: The walls, ceiling, and floor are all specified as concrete, with the thermal properties
specified in Table 3-1.
Ventilation: The ventilation rate of the smoke purge mode is 1.42 m3/s (3,000 cfm). This value
is used as a direct input parameter in the FPA correlation.
B-10
Table B-3. Summary of input parameters for FPA analysis of switchgear room scenario.
Parameter
Room height (H)
Room length (L)
Room effective width (We)
Room boundary material
Mech. Ventilation rate ( )
Fire elevation (Hf)
Value
5.2 m
26.5 m
18.5 m
Concrete
1.42 m3/s
2.4 m
Source
20C
See Table B-1
Figure B-1
Figure B-1
Calculation
Figure B-1. See Table 3-1 for properties.
From scenario description
From scenario description of cabinet height
and vent location.
Specified
Temperature: The FPA HGL temperature correlation for mechanically ventilated spaces is
expressed in non-dimensional terms as:
(B-4)
0.63
B-11
500
450
60
350
300
40
250
30
200
HRR(kW)
Temperature(C)
400
50
FPA[C]
HRR[kW]
150
20
100
10
50
0
0
600
1200
1800 2400
Time(s)
3000
0
3600
Figure B-4. Average HGL temperature (red line) from FPA correlation and HRR (blue line)
for SWGR cabinet fire scenario.
B-12
Fire: CFAST requires a user-specified, time-dependent HRR and stoichiometry for the
combustion of fuel and oxygen. The HRR is the combined cabinet/cable fire described above.
Figure B-7 shows the CFAST inputs for the fire taken directly from Table B-1.
B-13
Note: Values for Lower Oxygen Limit and Gaseous Ignition Temperature are set at default values.
Figure B-7. CFAST fire specification inputs for the SWGR scenario.
Materials: The walls, floor, and ceiling are specified as concrete in CFAST with properties as
previously described. The target properties are provided directly as input.
Ventilation: Mechanical ventilation and leakage are specified as input to CFAST directly from
the scenario description. CFAST uses three inlet and three outlet vents for the mechanical
ventilation at the heights specified in the scenario description. Horizontal placement of the
mechanical ventilation within the compartment does not affect the zone model calculation and is
not part of the input. Figure B-8 shows the CFAST inputs for the mechanical ventilation. Since
pressure for this scenario peaks approximately at the lower fan cutoff pressure threshold,
default values for these pressures are used as they should not impact the calculations.
B-14
Note: Values for Begin Dropoff At and Zero Flow At are set at default values.
Figure B-8. CFAST mechanical ventilation inputs for the SWGR scenario.
B-15
Figure B-9. CFAST inputs for cabinet and cable targets for the SWGR scenario.
B-16
B-17
Figure B-11. FDS/Smokeview rendering of the SWGR fire showing localized ignition of
extinction of secondary cable fires resulting from initial cabinet fire.
B-18
Model
Bias
Factor,
CFAST
FDS
CFAST
FDS
1.00
1.02
1.00
1.02
CFAST
FDS
0.81
0.85
Standard
Deviation,
Location
Predicted
Value
Critical
Value
Probability
of
Exceeding
205
205
205
205
0.937
1.000
0.177
0.000
6
6
0.576
0.159
B-19
1800
500
450
400
350
300
250
200
150
100
50
0
1600
CFAST
CFAST
1400
HRR (kW)
HRR (kW)
FDS
1200
FDS
1000
800
600
400
200
0
600
1200
1800
2400
3000
3600
600
Time (s)
1200
1800
2400
3000
3600
Time (s)
Figure B-12. Heat release rate inputs to CFAST and FDS for a SWGR cabinet fire scenario.
There are several ways to assess the initial fire sources potential to ignite the lowest cable tray
(Cable Tray A). Cables are considered damaged when the temperature reaches 205 C
(400 F) or the exposure heat flux reaches 6 kW/m2 (NUREG-1805, Appendix A). These criteria
are intended to indicate electrical failure, but are routinely also applied as ignition criteria. In
newer studies in NUREG/CR-7010, cable ignition was not observed at fluxes below 25 kW/m2,
and most often only with direct flame impingement. Handbook values for minimum ignition flux
for power and communication cables are reported in the range of 15 kW/m2 to 35 kW/m2 (SFPE
Handbook, Table 3-4.2). For this scenario, CFAST predicts that the flame height reaches the
cable tray in approximately 490 s, quite similar to the temperature-based prediction. Table B-5
shows the lowest cable trays estimated time to ignition for a variety of ignition criteria. For this
simulation, 490 s was chosen.
Table B-5. Estimated time to ignition of lowest cable tray by CFAST for the SWGR cabinet fire.
Ignition Criterion
Gas temperature 205 C
Cable temperature 205 C
Heat flux 6 kW/m2
Heat flux 15 kW/m2
Flame impingement
Time
270 s
860 s
490 s
740 s
490 s
The CFAST and FDS temperature predictions for the Cable Tray A cables are shown in Figure
B-13. FDS predicts cable failure in Tray A at about 495 s, CFAST in about 600 s. Peak
temperatures from both models are well above the failure criteria for the cables, so it can be
expected that the cables will ignite and provide an additional source of fire.
B-20
Temperature (C)
700
CFAST
600
500
FDS
400
300
200
100
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure B-13. Estimated temperatures for Cable Tray A directly above the fire source for a
SWGR cabinet fire scenario.
Qualitatively, the results of the CFAST and FDS predictions are quite different. The radiation
from the fire source in CFAST is calculated based on a point source fire positioned at the base
of the fire. Thus, once the fire grows and the flame height approaches the target cable tray,
CFAST can be expected to underestimate the local cable temperature and heat flux, since the
cable would actually be immersed within the flames. CFAST does include an estimate of the
flame height, which can also be used as an indicator of damage to the cable. For this scenario,
CFAST predicts that the flame height will reach the cable tray in approximately 490 s, quite
similar to the temperature-based prediction. Past this point, CFAST estimates of the local target
temperature are expected to be under-predictions. FDS predictions include the impact of direct
flame impingement and immersion of the target in flames. Thus, the higher temperatures
predicted by FDS are expected.
Upon ignition of the bottom cable tray (Cable Tray A), the higher cable trays are ignited,
consistent with the FLASH-CAT model.
B-21
Cabinet A Temperature
6
180
CFAST
140
Temperature (C)
160
120
FDS
100
80
60
40
20
0
0
600
1200
1800
2400
3000
CFAST
4
FDS
3
2
1
0
3600
Time (s)
600
1200
1800
2400
3000
3600
Time (s)
Figure B-14. Estimated temperature and heat flux to a cabinet adjacent to the fire source
in a SWGR cabinet fire scenario.
19
NUREG/CR-6850 specifies gamma distributions for the various types of combustibles found within an NPP. Microsoft Excel
provides a built-in function (GAMMA.DIST) that calculates the probability density function given the parameters and . In this
case, these parameters are 2.6 and 67.8, respectively.
B-22
HRR Distribution
0.005
0.004
0.003
0.002
0.001
0.000
0
100
200
300
400
500
600
Applying Heskestads flame height correlation to the entire range of HRR, now taken as a
random variable, leads to a distribution of flame height shown in Figure B-16.
Flame Height Distribution
1.0
0.8
0.6
0.4
0.2
0.0
0
0.5
1.5
2.5
Figure B-16. Distribution of flame heights for the entire range of cabinet fires.
The cable tray is 1.5 m (4.9 ft) above the top of the cabinet. The probability that the flames from
a randomly chosen fire will reach the cables is equal to the area beneath the curve in Figure
B-16 for flame heights greater than 1.5 m (4.9 ft), or approximately 0.31. Consistent with the
guidance in NUREG/CR-6850, this resulting probability can be used as the severity factor for
the quantification of corresponding fire ignition frequencies. It is interesting to note that an
analysis that uses the 50th percentile fire would lead to the conclusion that the flame from the
B-23
B.6 Conclusion
This analysis has considered the potential that a fire in an electrical cabinet in a 4160 V SWGR
will damage overhead cables and adjacent electrical cabinets. Algebraic equations from the
FDTs and FIVE-Rev1, including the Heskestad flame height correlation and the Heskestad
plume temperature correlation, were used for screening purposes, to evaluate the potential for
damage as well as to determine whether more detailed analysis with CFAST and FDS was
warranted. The algebraic equations demonstrate that the calculated flame height from the
cabinet fire would be high enough to potentially ignite the lowest of the three horizontal cable
trays located directly above the cabinet fire. They also demonstrate that the calculated fire
plume temperatures are high enough at all three horizontal cable trays located directly above
the cabinet fire to potentially damage cables in all three trays. As applied in this scenario, the
algebraic equations demonstrate that a more detailed analysis with CFAST and FDS is
warranted.
The more detailed analyses with CFAST and FDS demonstrate that the cabinet fire is likely to
ignite the electrical cables in the lowest cable tray directly above the cabinet fire in
approximately 10 min. The additional cable trays directly above the lowest tray would then
ignite in turn. In addition, the models indicate that it is possible that the combined cabinet and
cable fire could damage adjacent cabinets. These conclusions are based on an analysis of a
severe (i.e., 98th percentile) fire. A subsequent analysis involving the propagation of the entire
HRR distribution through a simple flame height model indicates that 31 % of all possible fires
within the cabinet would ignite the overhead trays and the resulting fire could damage adjacent
cabinets.
B.7 References
1. Jones, W., R.D. Peacock, G.P. Forney, and P.A. Reneke, CFAST - Consolidated Model of
Fire Growth and Smoke Transport (Version 6), Technical Reference Guide, SP 1026,
National Institute of Standards and Technology, Gaithersburg, MD, 2009.
2. NIST SP 1018-5, Fire Dynamics Simulator (Version 5), Technical Reference Guide,
Volume 3, Experimental Validation, 2010.
3. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
4. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
5. NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
6. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
B-24
3.
B-25
B-26
APPENDIX C
C-1
C-2
Material
Ceramic Fiber
Insulation
Cable
Parameter
Thickness (2 layers)
Thermal conductivity
Density
Specific heat
Emissivity
Diameter
Jacket thickness
Insulation/jacket conductivity
Insulation/jacket density
Insulation/jacket specific heat
Mass per unit length
Conductor mass fractions
Value*
5 cm
0.06 W/m/K
128 kg/m3
1.07 kJ/kg/K
0.9
15 mm
2 mm
0.192 W/m/K
1380 kg/m3
1.289 kJ/kg/K
0.4 kg/m
33% PE/PVC, 67% copper
*Source: Product literature (ERFBS) and NUREG/CR-6850 (EPRI 1011989), Volume 2, Appendix R
(PVC cable insulation).
Fire Protection Systems: As shown in Figure C-1, a smoke detector and a sprinkler are
located in the pump room. However, to determine whether the barrier system alone will provide
the required protection, the fire detection and suppression systems are not credited in the fire
scenario under consideration.
Ventilation: There is one supply and one return air vent, each with an area of 0.25 m2 (2.7 ft2),
providing a volume flow rate of 0.25 m3/s (530 cfm). The locations are shown in Figure C-1.
The ventilation system continues to operate during the fire, with no changes brought about by
fire-related pressure effects. This does not imply that the fire does not impact the ventilation
system, but rather that there is typically limited information about the ventilation network that
connects to a given compartment. The pump compartment has one door; it is 1.1 m (3.6 ft)
wide and 2.1 m (6.9 ft) tall. The door is normally closed, but it is opened 10 min after ignition by
C-3
Parameter
Value
Source
Effective Fuel Formula
CnH2n+2
Specified as C14H30
Mass burning rate
0.039 kg/s.m2 NUREG-1805 Table 3-4
Fuel volume
190 L
Specified
3
Fuel density
760 kg/m
NUREG-1805 Table 3-4
Heat of Combustion
46,000 kJ/kg NUREG-1805 Table 3-4
Heat of Combustion per unit
13,100 kJ/kg Huggett 1980, Average value
mass of oxygen consumed
CO2 Yield
2.64 kg/kg
SFPE Handbook, 4th ed., Table 3-4.16*
Soot Yield
0.059 kg/kg
SFPE Handbook, 4th ed., Table 3-4.16*
CO Yield
0.019 kg/kg
SFPE Handbook, 4th ed., Table 3-4.16*
Radiative Fraction
0.34
SFPE Handbook, 4th ed., Table 3-4.16*
2
Mass Extinction Coefficient
8700 m /kg
Mulholland and Croarkin (2000)
*Material identified as Hydrocarbon in SFPE Handbook was used to derive the properties.
C-4
C.3.3 Validation
Table C-3 lists various important non-dimensional parameters and the ranges for which the
validation study NUREG-1824 (EPRI 1011999) is applicable. The only parameters that fall
outside the validation parameter space are the equivalence ratios for the mechanically
ventilated portion of the scenario (first 10 minutes, while the door is closed) and for the natural
ventilation portion of the scenario (after 10 minutes when the door is opened). In both cases,
the high equivalence ratios for the compartment are a result of the relatively large fire and low
airflows.
For MAGIC, a sensitivity case is run, in which double doors, rather than a single door, are
opened after 10 minutes to determine whether an increase in airflow would cause higher
temperatures in the room. With the enlarged opening, the equivalence ratio is 0.5, putting it
within the verification and validation (V&V) range. The results of the sensitivity case, presented
in Section C.5.2, show that the HGL temperatures predicted in the room are not sensitive to the
size of the door opening.
C-5
Quantity
Fire Froude
Number
4934 kW
1.2 kg/m
Flame Length,
L , relative to the
Ceiling Height,
Hc
3.7
1.02
9.8 m/s
0. 93
1.02
Equivalence
Ratio, , based
on natural
ventilation
0.23
0.23 0.5
Target Distance,
r, relative to the
Fire Diameter, D
1.2 kg/m
0.5
1.2 1.7
N/A
0.04 0.6
No
0.04 0.6
No
0.6 5.7
Yes
2.2 5.7
N/A
0.38 kg/s
0.6
N/A
C-6
Yes
0.99
Notes:
(1)
0.2 1.0
0.07 kg/s
2.31 m 2.1 m
2.8 m
4.9 m
1.9
Yes
5.5
0.25 m /s
4934 kW
13,100 kJ/kg 0.38 kg/s
0.23
9.4 m
4.9 m
Compartment
Aspect Ratios
4934 kW
13,100 kJ/kg 0.07 kg/s
0.23
0.4 2.4
4.8 m
N/A
Equivalence
Ratio, , based
on the
mechanical
ventilation rate
In
Range?
0.9
0.99
1.9 m 3.7
Ceiling Jet
Radial
Distance,rcj ,
relative to the
Ceiling Height,
Hc
NUREG1824
Validation
Range
(2)
(3)
0.039 kg/m s
46,000 kJ/kg
2.75 m
4,934 kW
(C-1)
The oxygen needed to sustain the fire is calculated from the following equation:
4934 kW
13,100 kJ/kg
0.377 kg/s
(C-2)
where
is the heat of combustion per unit mass of oxygen consumed. The quantity of
oxygen provided by the ventilation system is calculated by multiplying the oxygen content (0.23)
by the density and the ventilation rate of the air:
0.23
0.23
1.2 kg/m
0.25 m /s
0.069 kg/s
(C-3)
The oxygen provided by the ventilation system is much lower than the amount needed to
sustain the fire. The oxygen initially in the room can provide the additional oxygen needed for
combustion for a short time. The available oxygen in the room, calculated from the room
dimensions (Table C-4), is:
0.23
0.23
1.2 kg/m
2.81
9.39
4.9 m
35.7 kg
(C-4)
The oxygen initially in the room can sustain the fire for an amount of time equal to the oxygen
quantity in the room divided by the consumption rate minus the ventilation supply rate, as shown
below:
35.7 kg
0.377 kg/s 0.069 kg/s
C-7
116 s
(C-5)
13,100 kJ/kg
904 kW
(C-6)
These results show that the oxygen supply available to the room will only allow a fire of reduced
size to burn until the door is opened (under-ventilated condition).
Parameter
Effective Length
Effective Width
Fire Diameter
Peak Heat Release Rate
Fire Duration
Mass of Fuel (kg)
Stoichiometric Mass-Oxygen
Specific Area
C-8
Value
2.81 m
9.39 m
1.87 m
4,934 kW
1345 s
144.4
3.5 g/g
513 m2/g
Fire: For this scenario, the lubricating oil is preheated prior to the spill, such that the HRR
reaches the peak immediately upon fire initiation, as shown in the HRR curve plotted in Figure
C-3. The lower oxygen level is specified to be 10%. Using the specified spill area and volume,
the spill depth is calculated as 0.069 m (0.23 ft).
The fire is modeled as a single circular area of equivalent diameter. The actual entrainment for
the pool fire is proportional to the perimeter of the fire, which is significantly greater than the
perimeter of the circular area. However, the enclosure is small and the smoke filling rates are
expected to be short regardless of the fire shape.
The fire duration, , is determined from the pool depth, , density, , and burning rate,
C-9
1345 s
22.4 min
:
(C-7)
HRR (kW)
Input HRR
4000
3000
2000
1000
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure C-3. Heat release rate curve for lubricating oil fire.
The location of the fire is placed at the edge of the dike, closest to the target cable. The total
mass of fuel is 144.4 kg, calculated from the volume multiplied by the density from Table C-2.
The stoichiometric mass-oxygen-to-fuel ratio, , is calculated using Equation 22 from Chapter
3-4 of SFPE Handbook, 4th Edition and the values from Table C-1, as follows:
46,000 kJ/kg
13,100 kJ/kg
3.5 kg/kg
(C-8)
where is heat of combustion per unit mass of fuel consumed and O is the heat of
combustion per unit mass of oxygen consumed. One of the inputs required by MAGIC is the
specific area, s, a measure of the smoke generation, which is calculated as described in
NUREG-1824 (EPRI 1011999), Volume 6, Section 3.2.7:
0.059
8,700 m /kg
513 m /kg
(C-9)
where
is the mass extinction coefficient and is the soot yield, as listed in Table C-2. The
pyrolysis rate (g/s) is calculated for input to MAGIC by dividing the HRR values (4934 kW) at
each time step by the heat of combustion (46,000 kJ/kg). Other inputs needed for MAGIC are
listed in Tables 3-1, C-1, C-2, and C-4. Figure C-4 is a screenshot of the source fire in the
MAGIC input file.
C-10
C-11
The radius of the conductor in an equivalent single-conductor cable is needed for input to
MAGIC. The mass of the conductor per unit length is calculated from the mass fraction of the
conductor multiplied by the mass per unit length of the multi-conductor cable: 0.67 x 0.4 kg/m =
0.27 kg/m (values from Table C-1). To determine the conductor radius, , the mass per unit
length (0.27 kg/m) is set equal to the cross-sectional area times the density of copper, 8954
. Rearranging the equation to solve for conductor radius
kg/m3 (Table 3-1) or 0.27
results in the following:
0.27
8954
3.1 mm
(C-10)
The insulation thickness is calculated by cable radius jacket thickness copper thickness =
7.5 mm 2 mm 3.1 mm = 2.4 mm. Since the jacket and the insulation are both composed of
PE/PVC, the thicknesses are added together for a total thickness of 4.4 mm. As a result, the
ERFBS protected cable raceway is modeled with three layers (Figure C-6): ceramic fiber
blanket (5 cm), PE/PVC (4.4 mm), and copper (3.1 mm). The input screen for the layers of the
ERFBS and the cables is shown in Figure C-7.
C-12
Figure C-6. Representation of the ERFBS protected cable raceway for MAGIC.
Ventilation: Mechanical ventilation is maintained constantly during the simulation, using the
values provided above. MAGIC uses circular ducts, so the rectangular ducts seen in Figure C-1
are modeled as circular areas with equivalent diameter of 0.56 m (1.8 ft). As noted above, the
door is normally closed, but it is opened 10 minutes after ignition by the arriving fire brigade.
Before the door opens, leakage due to the doorway occurs via a 1.3 cm (0.5 in) gap under the
door. The MAGIC input screen for the doorway is shown in Figure C-8.
C-13
C-14
Fire: Due to the limited amount of validation data available for scenarios of this type and the
considerable uncertainties involved, the approach taken is to specify, rather than attempt to
predict, the burning rate of the fuel, even though the FDS model does provide the physical
mechanisms to estimate burning rates. The fire is specified in the diked area surrounding the
pump. Although FDS has a liquid fuel burning model, it is not used here because there is not
enough information about the fuel, and, more importantly, it lacks the exact geometry of the
pump and diked area. Typically, FDS would expect that the oil has formed a relatively deep
pool with relatively little influence by the surrounding solids. This is not the case here. Instead,
the specified burning rate, 0.039 kg/m2/s, is applied directly to the model over an area of
2.75 m2 (29.6 ft2), yielding a burning rate of 0.107 kg/s. The density of the oil is 0.76 kg/L, which
means that the oil burns at a rate of 0.141 L/s. At this rate, 190 L (50 gal) will require 1,348 s to
burn out. Note that this is slightly different from the burning duration of 1345 s computed for the
MAGIC input. The fire duration computation for FDS converts the mass data to volumetric data,
C-15
Model
Bias
Factor,
Standard
Deviation,
Predicted
Value
Critical
Value
Probability of
Exceeding
1.19
0.27
135
205
0.000
FDS
1.02
0.13
145
205
0.000
C-16
HRR (kW)
16000
6000
MAGIC
5000
MAGIC
FDS
4000
FDS
3000
Hand Calculation
HRR (kW)
20000
12000
8000
2000
1000
4000
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure C-10. HRR predicted by algebraic methods, MAGIC, and FDS for the pump room
fire scenario.
Temperature (C)
HGL Temperature
1000
900
800
700
600
MAGIC
FDS
ASTM E-119
500
400
300
200
100
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure C-11. HGL Temperature Predicted by MAGIC and FDS for the pump room fire
scenario.
C-17
(C-11)
is the integrated heat flux received by the ERFBS, is the HGL temperature,
is the
Here,
ambient temperature (20 C), is the convective heat transfer coefficient (about 0.025 kW/m2/K
in fully developed fires), and is the Stefan-Boltzmann constant (5.67 10-11 kW/m2/K4). Note
that Equation C-11 is the total energy transferred from a thermally thick (emissivity of 1) hot gas
to a cold target. It is intended only to compare the different temperature curves. In reality, the
net heat transferred to a target in the compartment decreases as the target heats up.
Applying Equation (C-11) to each of the HGL temperature curves in Figure C-11 yields values of
346 MJ/m2 for the ASTM E119 curve and approximately 40 MJ/m2 for both FDS and MAGIC.
This 40 MJ/m2 exposure corresponds to an approximately 14 min exposure within the standard
test furnace. Table C-6 lists the thermal exposure as a function of time in the standard test
furnace. It is also significant to note that the maximum predicted exposure temperature remains
lower than the maximum exposure temperature that the ERFBS protected raceway was
exposed to during the ASTM E119 fire test.
C-18
Thermal
Time
Exposure
(min)
(MJ/m2)
5
6
10
23
15
47
20
75
25
104
30
135
35
167
40
200
45
235
50
270
55
307
60
346
Direct Calculation of Cable Temperature
MAGIC and FDS have heat conduction algorithms to account for the multiple layers of insulation
and cable materials. The surface temperature predictions of the cables protected by the
ERFBS (ceramic fiber insulation in this case) are shown in Figure C-12. MAGIC predicts a
maximum cable surface temperature of approximately 135 C (275 F). FDS predicts a
maximum cable surface temperature of approximately 145 C (293 F). Note that although the
HGL temperature drops and then increases dramatically when the door opens, as shown in
Figure C-11, the cable temperature slowly rises. This is due to the thermal lag caused by the
ERFBS.
Cable Surface Temperature
Inside ERFBS
Temperature (C)
160
140
MAGIC
120
FDS
100
80
60
40
20
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure C-12. Cable surface temperature predicted by MAGIC and FDS for the pump room
fire scenario.
C-19
Temperature (C)
140
120
100
80
60
40
25% Tighter
Wrapping
20
25% Thinner
Wrapping
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure C-13. Cable surface temperature predicted by MAGIC for changes to insulation
wrapping.
0.23 0.5
0.23
0.5
4.62 m 2.1 m
4934 kW
13,100 kJ/kg 0.77 kg/s
0.77 kg/s
0.5
(C-12)
(C-13)
Figure C-14 shows the temperature comparison for the HGL and the cable surface temperature
(measured inside the ERFBS) for the base case and for the case with double doors. The plots
show that the results for both cases are very similar, indicating that the door size does not
significantly affect the results. Nevertheless, it is consistent with experimental data that the
scenario with the equivalence ratio closest to unity produces the highest enclosure temperature.
C-20
HGL Temperature
700
160
500
140
Double Doors
120
Double Doors
Temperature (C)
Temperature (C)
600
400
300
200
100
100
80
60
40
20
0
0
600
1200
1800
2400
3000
3600
Time (s)
600
1200
1800
2400
3000
3600
Time (s)
C-21
6000
6000
Input HRR
5000
Extended
HRR
5000
4000
3000
2000
MAGIC
4000
HRR (kW)
HRR (kW)
MAGIC
3000
2000
1000
1000
0
0
600
1200
1800
2400
3000
3600
600
1200
Time (s)
1800
2400
3000
3600
Time (s)
Figure C-15. HRR for base case and HRR sensitivity case.
HGL Temperature
250
600
MAGIC Base
Case
500
Extended
HRR
Temperature (C)
Temperature (C)
700
400
300
200
100
0
0
600
1200
1800
2400
3000
200
Extended HRR
150
100
50
0
3600
Time (s)
600
1200
1800
2400
Time (s)
Figure C-16. Temperature for base case and HRR sensitivity case.
C-22
3000
3600
C.6 Conclusion
This analysis has considered the potential for a relatively large lubricating oil spill fire in a
relatively small enclosure to damage a cable tray protected by an ERFBS. Algebraic
calculations, the zone model MAGIC, and the CFD model FDS were all used to evaluate the fire
conditions within the enclosure. MAGIC and FDS were used to calculate the thermal response
of the cables to these calculated fire conditions.
Based on the specified lubricating oil spill area and burning characteristics, a fire of
approximately 5 MW is expected. However, after the rapid consumption of the limited quantity
of air in the room, the mechanical ventilation to the enclosure could only support a HRR of less
than 1 MW before the door to the enclosure opens after 10 min. This analysis suggests that to
avoid rapid fire escalation, doors to such rooms should not be opened until firefighters are
prepared to suppress the fire, and, even then, the potential for rapid fire escalation should be
considered.
Two different strategies were applied to assess the integrity of the ERFBS. Because the
thermal and chemical properties of the insulating material are only partially known, it is
necessary to implement an alternative technical approach of comparing the predicted HGL
temperatures from the models with the standard temperature curve under which the ERFBS
received an hour rating. Since the predicted HGL temperatures do not lie completely within the
standard curve, a simple integrated heat flux calculation was performed to demonstrate that the
ERFBS received approximately 10 times the thermal exposure in the standard fire endurance
test than is predicted by the two models.
A second strategy for assessing the integrity of the ERFBS was to directly calculate the heat
penetration through the insulating blankets using the thermal material properties of the cables
and the ERFBS. Both models predicted cable temperatures below the reported critical values.
Based on the two approaches to determine its performance, the ERFBS is expected to prevent
the cables from reaching temperatures that would limit their functionality in the event of a fire
involving burning spilled lubricating oil. This conclusion is based on certain expected burning
behavior of the lubricating oil during the under-ventilated stages. A sensitivity study on the
burning behavior of the lubricating oil concluded that the results could change if the burning rate
decreases during the under-ventilated stage. The results are also shown to be sensitive to the
thermal properties of the ERFBS material. Further research or testing of the ERFBS thermal
properties may be necessary to confirm the initial conclusion.
C.7 References
1. Gay, L., C. Epiard, and B. Gautier, MAGIC Software Version 4.1.1: Mathematical Model,
EdF HI82/04/024/B, Electricit de France, France, November 2005.
2. Gay, L., and J. Frezabeu, Qualification File of Fire Code MAGIC version 4.1.1, EdF HI82/04/022/A, Electricit de France, France, December 2004.
3. Hugget, C., Estimation of Rate of Heat Release by Means of Oxygen Consumption
Measurements, Fire and Materials, 4:61-65, 1980.
4. Mulholland, G. W., and Croarkin, C., Specific Extinction Coefficient of Flame Generated
Smoke, Fire and Materials, 24:227230, 2000.
C-23
Pump_Room.cas
Pump_Room_tighter_wrapping.cas
Pump_Room_thinner_wrapping.cas
Pump_Room_2Doors.cas
Pump_Room_extendedHRR.cas
C-24
APPENDIX D
20
The cable failure temperature is based on experiments conducted with Cable #14, an XLPE/CSPE, 3 conductor control cable.
D-1
Parameter
Effective Fuel Formula
Peak HRR
Time to reach peak HRR
Heat of Combustion
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Value
C3H4.5Cl0.5
702 kW
720 s
10,300 kJ/kg
0.63 kg/kg
0.175 kg/kg
0.082 kg/kg
0.53
Source
Combination of polyethylene and neoprene
NUREG/CR-6850 (EPRI 1011989), App. G
NUREG/CR-6850 (EPRI 1011989), App. G
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
SFPE Handbook, 4th Ed., Table 3-4.16
Ventilation: The compartment is normally supplied with three air changes per hour (ACH).
The supply and return vents are indicated on the drawing. The two doors are normally closed.
Normal heating, ventilation, and air conditioning (HVAC) operations continue during the fire, and
the doors remain closed. The volume of the compartment is 882 m3 (31,150 ft3); thus, three air
changes per hour is equivalent to a volume flow rate of 0.735 m3/s.
21
4 / , where
D-2
D-3
Figure D-2. Typical electrical cabinet in the lower part of the SWGR.
D-4
D.3.3 Validation
The principal source of validation data justifying the use of the above-listed fire models for this
scenario is the U.S. Nuclear Regulatory Commission/Electric Power Research Institute
(NRC/EPRI) verification and validation (V&V) study documented in NUREG-1824 (EPRI
1011999). The National Institute of Standards and Technology (NIST) has expanded the
NRC/EPRI V&V to include the latest versions of CFAST (6.1.1) (Peacock, 2008) and FDS
(5.5.3) (McGrattan, 2010). Also, CFAST and FDS utilize the cable failure algorithm, THIEF,
which was developed and validated in NUREG/CR-6931 (Vol. 3).
D-5
D-6
Quantity
Fire Froude
Number
702 kW
1.2 kg/m
Flame
Length,
,
relative to the
Ceiling
Height,
Ceiling Jet
Radial
Distance,rcj ,
relative to the
Ceiling
Height, Hc
Equivalence
Ratio, , as
an indicator
of the
Ventilation
Rate
Compartment
Aspect Ratio
(Lower |
Upper)
Target
Distance, r,
relative to the
Fire
Diameter, D
2.4 m 2.5 m
3.0 m
/
3.7
1.02
0.48 m 3.7
9.8 m/s
8.5 m
3.0 m
0.23
2.8 ;
0.23
8.5 m
3.0 m
3.6
1.02
0.735 m /s
8.6 m
9.1 m
2.8
1.1 m
0.5 m
2.2
0.4 2.4
No
0.2 1.0
No
1.2 1.7
N/A
0.04 0.6
Yes
0.6 5.7
Yes
2.2 5.7
Yes
1.6
702 kW
13,100 kJ/kg 0.2 kg/s
1.2 kg/m
In
Range?
3.6
2.5 m
NUREG1824
Validation
Range
0.9 ;
0.3
0.2 kg/s
8.5 m
9.1 m
0.9
Notes:
(1)
(2)
(3)
D-7
0.53
4
702 kW
1.1
24.5 kW/m
(D-1)
This estimate does not include any contributions to the heat flux by the HGL or by the flame
extension beneath the ceiling. However, it does indicate that the heat flux to the adjacent
cabinet could exceed the critical heat flux by a relatively large margin. Consequently, this
scenario would warrant more detailed analysis with either a zone model or a CFD model.
D-8
2264 kg/m
(D-2)
D-9
Figure D-6. CFAST fire inputs for two-height ceiling SWGR scenario.
D-10
Figure D-7. CFAST target inputs for two-height ceiling SWGR scenario.
Ventilation: The two individual compartments used to model the entire room are connected by
a single large vent. Although the size of this vent relative to the compartment size is not typical
of a zone model application, the simple two-compartment geometry of the space and the more
dominant mechanical ventilation flow from one side of the SWGR to the other should minimize
any uncertainty in the calculation resulting from the large connecting vent. Mechanical
ventilation is included at the specified height and with the specified volume flow applied to the
single supply (in the low-ceiling space) and return (in the high-ceiling space). Additionally, since
zone fire models use compartments that are completely sealed unless otherwise specified, a
typical leakage vent, 13 mm (0.5 in) in height, is included at the bottom of each closed doorway
to reflect the fact that the doorways are not totally airtight. Figure D-8 shows the CFAST inputs
for these natural vents.
D-11
Figure D-8. CFAST inputs for vents connections in two-height ceiling SWGR scenario.
D-12
D-13
Model
Bias
Factor,
CFAST
FDS
CFAST
FDS
CFAST
FDS
CFAST
FDS
1
1.02
1
1.02
1
1.02
1
1.02
CFAST
FDS
CFAST
FDS
CFAST
FDS
CFAST
FDS
0.81
0.85
0.81
0.85
0.81
0.85
0.81
0.85
Standard
Deviation,
Target
Predicted
Value
Critical
Value
Probability of
Exceeding
0.460
0.000*
0.950
0.997
0.112
0.000
0.000
0.000
0.911
0.006*
0.974
1.000
0.823
0.997
0.000
0.000
D-14
Cabinet Temperature
500
25
400
20
Temperature (C)
CFAST
FDS
15
FDTs
10
5
CFAST
FDS
300
200
100
0
0
0
600
1200
1800
2400
3000
3600
600
1200
1800
2400
3000
3600
Time (s)
Time (s)
Figure D-10. Heat flux and temperature predictions for the adjacent cabinet.
3
2
3
702 kW
2
100 C
300 C 20 C
376 kW
(D-3)
Both FDS and CFAST predicted cable temperatures for Tray C indicate that the cables are
unlikely to fail.
D-15
800
CFAST
100
FDS
FDS
400
200
80
60
40
20
0
0
600
1200
1800
2400
3000
3600
600
1200
Tim e (s)
2400
3000
3600
300
CFAST
1800
Tim e (s)
400
FDS
200
100
CFAST
20
FDS
15
10
5
0
0
0
600
1200
1800
2400
3000
3600
600
1200
Tim e (s)
1800
2400
3000
3600
Tim e (s)
Cable C Temperature
100
2.5
CFAST
80
CFAST
Temperature (C)
CFAST
600
FDS
60
40
20
0
0
600
1200
1800
2400
3000
3600
2.0
FDS
1.5
1.0
0.5
0.0
0
Time (s)
600
1200
1800
2400
3000
3600
Time (s)
Figure D-11. Summary of the cable temperature and heat flux predictions for the
MCC/SWGR.
D-16
D.6 Conclusion
The purpose of the calculations in this example is to predict if and when various components
within a compartment will become damaged due to a fire in the MCC. The fire model analyses
performed for this scenario indicate that the fire would damage the cables in Tray A because all
the models (FDTs, CFAST, FDS) predict that the flames would directly impinge on the cables
themselves.
CFAST and FDS predict that the cables in Tray B are likely to be damaged based on the
heat flux criterion. However, neither model predicts that the interior cable temperatures are
likely to be high enough to cause failure.
Neither FDS nor CFAST predicts that the cables in Tray C would be damaged.
A point source heat flux analysis indicates that the adjacent cabinet housing would be
exposed to a heat flux that would cause damage. Even though FDS does not predict
damage, its predictions of heat flux to surfaces very near the adjacent cabinet are
sufficiently high to cast doubt on the conclusion that the cabinet would not be damaged.
Small changes in the positions of various obstructions could easily change the predicted
heat flux by an order of magnitude. Even though the point source method tends to overpredict the heat flux to targets close to the fire, there is too much uncertainty in the
geometric configuration to accept the validity of the more detailed calculation.
D.7 References
1. NIST SP 1018-5, Fire Dynamics Simulator (Version 5), Technical Reference Guide, Volume
3, Experimental Validation, 2010.
2. NIST SP 1019, Fire Dynamics Simulator (Version 5) Users Guide, 2007.
3. NIST SP 1086, CFAST Consolidated Model of Fire Growth and Smoke Transport (Version
6), Software Development and Model Evaluation Guide, 2008.
4. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
5. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
6. NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
7. NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 2: Cable Fire
Response Data for Fire Model Improvement, 2007.
8. NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 3: ThermallyInduced Electrical Failure (THIEF) Model, 2007.
D-17
D-18
APPENDIX E
E-1
317 kW
480
480 s
3
480 s
(E-1)
Here,
is the peak HRR of 317 kW, which is reached following a t-squared growth curve in
480 s. Solving for yields a total burning time of 799 s, which is rounded up to 800 s.
While the exact composition of the trash is unknown, typical waste at the plant includes wood
scraps and polyethylene protective suits. For this scenario, the trash is specified as being
comprised of equal parts of each. The radiative fraction22 and product yields, like the heat of
combustion, are taken to be averages of values for red oak and polyethylene (SFPE Handbook,
4th ed., Table 3-4.16) and shown in Table E-1. For the purpose of modeling, the fuel molecule is
specified as C4H7O2.5.
22
E-2
Parameter
Effective Fuel Formula
Peak HRR
Time to reach peak HRR
Heat of Combustion
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Value
C4H7O2.5
317 kW
480 s
30,400 kJ/kg
2.0 kg/kg
0.038 kg/kg
0.014 kg/kg
0.40
Source
Specified
NUREG/CR-6850 (EPRI 1011989), App. G
NUREG/CR-6850 (EPRI 1011989), App. G
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
E-3
E-4
Figure E-2. Photograph of the CSR used for analysis. Note that the cables are located in
the trays in the overhead.
E-5
E.3.3 Validation
NUREG-1824 (EPRI 1011999) contains experimental validation results for CFAST and FDS that
are appropriate for this scenario. These experiments include ventilation effects on, heat fluxes
to, and temperatures of various targets, particularly cables. Fire sizes in these experiments
E-6
E-7
Quantity
Fire Froude
Number
317 kW
1.2 kg/m
0.8 m 1.6 m
4.0 m
Flame Length,
Lf , relative to the
Ceiling Height,
Hc
Ceiling Jet
Radial
Distance,rcj ,
relative to the
Ceiling Height,
Hc
3.7
1.02
0.6 m 3.7
0.23
9.8 m/s
0.23
40 m
4.0 m
Compartment
Aspect Ratio
Target Distance,
r, relative to the
Fire Diameter, D
1.2 kg/m
1.0
1.02
2.3 m
0.6 m
3.8
Yes
0.2 1.0
Yes
1.2 1.7
N/A
0.04 0.6
Yes
0.6 5.7
No
2.2 5.7
Yes
1.6 m
0.06
1.4 m /s
18.5 m
4.0 m
10
0.4 2.4
0.6
317 kW
13,100 kJ/kg 0.4 kg/s
In
Range?
1.0
Equivalence
Ratio, , as an
indicator of the
Ventilation Rate
NUREG1824
Validation
Range
4.6
0.4 kg/s
Notes:
(1)
, is the sum of the height of the fire from the floor and
E-8
Air Supply
(2 Vents)
Smoke
Detector
Tray B
Air Return
(2 Vents)
Tray A
Bottom
Figure E-4. Schematic diagram of transient trash fire in CSR (not to scale).
The first step to estimate the fire-generated conditions using algebraic models is to determine
whether the cables in the cable trays located directly above the trash fire are likely to be
damaged and potentially ignited by the trash fire. As shown in Figure E-3, the top of the trash
receptacle is located at an elevation of 0.8 m (2.6 ft), and the lowest overhead cable tray is
located at an elevation of 1.8 m (5.9 ft), which is 1.0 m (3.3 ft) above the top of the trash
receptacle. As shown in Table E-2, the flame length of the cabinet fire is calculated to be 1.6 m
(5.2 ft) at the peak heat release rate of 317 kW, so this empirical correlation for flame length can
be used to confirm that at least some of the overhead cables would likely be damaged in this
scenario. However, due to the metal barrier on the bottom of the lowest cable tray, algebraic
models cannot be used to predict whether the cables in the lowest cable tray would ignite.
The next step is to calculate the fire plume temperatures that develop from the trash fire to
determine whether cable trays A or B located above the trash fire would be damaged by the
trash fire alone. The Heskestad plume temperature correlation included in the Fire Dynamics
Tools (FDTs) and the Fire-Induced Vulnerability Evaluation (FIVE-Rev1) was used to calculate
the plume centerline temperature above the trash fire. The results, shown in Figure E-5, show
that the plume temperature of all the cable trays would exceed the cable damage threshold
temperature of 205 C (400 F). However, the Heskestad plume temperature correlation is
based on an unobstructed plume. The obstruction caused by the position of the cable trays
within the fire plume would alter the actual fire plume entrainment and temperatures.
Nonetheless, these results demonstrate that the potential for damage and ignition of the cable
trays located above the transient trash fire warrants more detailed analysis.
E-9
Plumetemperatures Heskestadcorrelation
1400
Temperature(C)
1200
1000
800
Bottom
TrayA
TrayB
600
400
200
0
0
600
1200
1800
Time(s)
2400
3000
3600
Figure E-5. Plume temperatures at cable trays located above a transient trash fire.
The FPA forced ventilation correlation is used to estimate the average HGL temperature of the
CSR resulting only from the transient trash fire, based on the parameters described in the
following subsections.
Geometry: The FPA correlation requires room dimensions to be specified in terms of length,
width, and height. For this example, the selected compartment is a rectangular parallelepiped,
as shown in Figure E-1, so its length, width, and height are specified directly.
Fire: As applied to this scenario, the FPA correlation is used with the time-dependent HRR
specified for the transient trash fire. This HRR history is shown in Figure E-6.
Materials: The walls, ceiling, and floor are all specified as concrete, with the thermal properties
specified in Table 3-1.
Ventilation: The ventilation rate in the smoke purge mode is 1.4 m3/s (3,000 cfm). This value
is used as a direct input parameter in the FPA correlation.
Table E-3. Summary of input parameters for FPA analysis of CSR scenario.
Parameter
Room height (H)
Room length (L)
Room width (W)
Room boundary material
Mech. Ventilation rate ( )
Fire elevation (Hf)
Ambient temperature (Ta)
Fire parameters
Value
4.0 m
40.0 m
18.5 m
Concrete
1.4 m3/s
0.8 m
20C
See Table E-1
Source
Figure E-1
Figure E-1
Figure E-1
Figure E-1 See Table 3-1 for properties.
From scenario description
From scenario description
Specified
E-10
60
350
50
300
250
40
200
30
150
20
HRR(kW)
Temperature(C)
HGLtemperaturecalculation FPAcorrelation
FPA[C]
HRR[kW]
100
10
50
0
0
600
1200
1800 2400
Time(s)
3000
0
3600
Figure E-6. Average HGL temperature (red line) from FPA correlation and heat release rate (blue
line) for the CSR trash fire scenario.
E-11
Geometry: Since zone models are concerned with volumes and not physical length and width,
the volume of the computer compartment, as well as the numerous cable trays, ductwork, and
beams, was modeled in CFAST with a cross-sectional area that varies with height. Table E-4
shows the cross-sectional area as a function of the height calculated from the compartment
geometry shown in Figure E-1. Figure E-8 shows the CFAST inputs for compartment geometry.
Table E-4. Cross-sectional area as a function of height used for CFAST calculation.
Floor Level
Bottom of Cable Trays
Bottom of Obstructions
HVAC Ductwork
Top of Obstructions
Top of Cable Trays
Ceiling Level Obstructions
Ceiling Level
E-12
Height
(m)
0
1.8
2.2
2.4
2.8
3.2
3.6
4
Area
(m2)
700.04
635.74
483.74
514.89
634.74
699.04
291.46
291.46
Figure E-8. CFAST inputs for compartment geometry for the CSR scenario.
Fire: The specified fire is input directly. The combustion chemistry in CFAST is described, at a
minimum, by the production rates of CO and soot. Figure E-9 shows the CFAST inputs for the
fire location in the scenario taken directly from Table E-1.
Cables: In CFAST, target temperatures are calculated using a one-dimensional cylindrical heat
transfer calculation based on the material properties and cable diameter, as specified in the
scenario description. Following the Thermally-Induced Electrical Failure (THIEF) methodology
in NUREG/CR-6931, Vol. 3, electrical functionality is lost when the temperature just inside the
jacket of a thermoplastic cable reaches 205 C (400 F). Specific heat and thermal conductivity
for the cables are taken to be 1.5 kJ/kg/K and 0.2 W/m/K, as specified in the THIEF model.
Cable density is determined from the mass per unit length and cross-sectional area of the cable
from the scenario description, 0.4 0.015 2 2264 kg/m3. To account for the shielding
of the cables on the lower surface of the cable tray, the CFAST input for the normal vector from
the cable surface is directed upwards. This effectively shields the cables from the fire below
while exposing them to the surrounding gas temperature for convection and to the hot upper
gas layer for radiation. Figure E-10 shows the CFAST inputs for the three cable targets above
the fire.
2
E-13
Figure E-9. CFAST inputs for the fire in the CSR scenario.
E-14
Figure E-10. CFAST inputs for cable targets above the fire source in the CSR scenario.
Ventilation: The supply and return airflow rates are input directly into CFAST. Upon smoke
detector activation, mechanical ventilation fans stop and dampers are closed. Thus, before a
stop time for the fans can be specified, the time to smoke detector activation is needed. This
requires that CFAST be run with the fans on for the entire time to find the first smoke detector
activation. The model is then re-run using the smoke detector activation time as the fan stop
time. Figure E-11 shows the CFAST inputs for the natural ventilation.
E-15
Figure E-11. CFAST inputs for natural ventilation for the CSR scenario. Note the opening of the
doorway at 170 s, the time estimated for smoke detectors to activate.
Fire Detection: Although there are multiple smoke detectors in the space, the closest detector
is the only one that needs to be modeled to determine time to detection, based on the fact that
the nearest detector will be exposed to the greatest concentration of smoke products and the
highest gas temperatures, thus leading to the earliest response. There are no geometric or
ventilation features that would prevent this from being the case in the example considered. In
CFAST, there is no direct way of calculating smoke density for smoke detector activation.
Instead, the smoke detector is modeled as a sprinkler with a very low activation temperature
and response time index (RTI). An activation temperature of 30 oC (86 oF) and an RTI of
5 (m/s)1/2 were used for this scenario. Figure E-12 shows the CFAST inputs for the detectors in
the scenario.
E-16
Figure E-12. CFAST inputs for smoke detectors in the CSR scenario.
Geometry: The interior of the compartment is modeled, and all obstructions have been
included. To get increased resolution in the area of interest, multiple meshes are used. The
finest mesh has a 10 cm (4 in) resolution and spans a volume surrounding the trash can, which
E-17
E-18
HGL Temperature
400
50
300
40
Temperature (C)
HRR (kW)
CFAST
200
FDS
100
0
CFAST
30
FDS
20
10
600
1200
1800
2400
3000
3600
Time (s)
600
1200
1800
2400
3000
3600
Time (s)
Figure E-14. HRR and estimated HGL temperature for CSR scenario.
Because of the uncertainty in all models smoke detector activation predictions and the
uncertainty associated with the possible ignition of cables in the trays just above the fire, it is
difficult to predict whether or not the CO2 suppression system would be activated in time to
prevent possible cable ignition. No validation results are available that address time to detector
activation. Thus, the analysis does not consider activation of the suppression system.
Table E-5 summarizes the results of all the models for the chosen damage criteria. For each
predicted value, a calculation is performed to determine the probability of exceeding the critical
value. The procedure for calculating this probability is given in Chapter 4, and it accounts for
the model bias and scatter. The purpose of this table is to highlight the criterion that is most
likely to be exceeded so that further analysis can be focused on this criterion and the model or
models that predict it. Each criterion is discussed in greater detail in the following subsections.
E-19
Model
CFAST
FDS
CFAST
FDS
CFAST
FDS
CFAST
FDS
CFAST
FDS
CFAST
FDS
Bias
Factor,
1
1.02
1
1.02
1
1.02
Standard
Deviation,
M
Location
Predicted
Value
Critical
Value
Probability
of
Exceeding
205
205
205
205
205
205
0.893
0.000
0.472
0.000
0.003
0.000
6
6
6
6
6
6
0.367
0.001
0.091
0.000
0.000
0.001
0.81
0.85
0.81
0.85
0.81
0.85
0.47
0.22
0.47
0.22
0.47
0.22
Model
CFAST
FDS
Time (s)
170 s
160 s
E-20
E-21
Flame Height
Bottom CableTemperature
2.0
350
CFAST
FDS
200
CFAST
1.5
250
Temperature (C)
300
150
100
50
1.0
0.5
0.0
0
0
600
1200
1800
2400
3000
3600
600
1800
3000
3600
Cable A Temperature
250
CFAST
150
FDS
200
100
50
0
0
600
1200
1800
2400
3000
CFAST
FDS
2
1
0
3600
600
1200
Time (s)
1800
2400
3000
3600
Time (s)
Cable B Temperature
140
2.0
120
CFAST
CFAST
1.5
100
Temperature (C)
2400
Time (s)
Time (s)
Temperature (C)
1200
FDS
80
60
40
20
0
FDS
1.0
0.5
0.0
600
1200
1800
2400
3000
3600
Time (s)
600
1200
1800
2400
Time (s)
E-22
3000
3600
0.010
0.008
0.006
0.004
0.002
0.000
0
100
200
300
400
500
23
NUREG/CR-6850 specifies gamma distributions for the various types of combustibles found within an NPP. Microsoft Excel
provides a built-in function (GAMMA.DIST) that calculates the probability density function given the parameters and . In this
case, these parameters are 1.8 and 57.4, respectively.
24
A value of 0.35 was used for the radiative fraction in determining the convective heat release rate.
E-23
0.010
0.008
0.006
0.004
0.002
0.000
0.008
0.006
0.004
0.002
0.000
100
200
300
400
500
600
100
200
300
400
500
E.6 Conclusion
This analysis considers the potential impact of a transient trash fire on a stack of horizontal
cable trays located directly above the transient fire source; it also considers the potential for the
cables in these trays to ignite.
Algebraic equations from the FDTs and FIVE-Rev1, including the Heskestad flame height
correlation and the Heskestad fire plume temperature correlation, were used to evaluate the
potential for the transient trash fire to ignite or damage cables in the cable trays located directly
above this fire source. These calculations showed that the unobstructed flame height would
reach multiple cable trays at the bottom of the stack; they also showed that the unobstructed
plume temperatures would exceed the cable damage temperature in all trays located above the
fire source. However, these calculations do not account for the impact of the cable trays on the
actual flame height and plume temperatures, so more detailed analysis with CFAST and FDS is
needed.
The FPA correlation included in the FDTs and FIVE-Rev1 was used to evaluate whether the
transient trash fire alone would cause damage to cable trays not located directly above this fire
source. This calculation showed that the HGL temperature from the transient trash fire alone is
well below the damage temperatures for the cables.
The more detailed analysis with CFAST and FDS shows that a 317 kW waste bin fire beneath a
vertical array of horizontal cable trays is unlikely to damage cables in the trays located three
(Tray A) and six (Tray B) levels above the fire. Both CFAST and FDS estimate peak
temperatures and heat fluxes below the failure criteria for cables in the third tray from the
bottom. Based on FDS calculations, temperatures and heat fluxes for the protected lowest
cable tray are well below critical values. Estimates from CFAST for unprotected cables
demonstrate the importance of the protection afforded by the solid metal lower surface of the
cable trays.
E-24
600
E.7 References
1. NIST SP 1086, Consolidated Model of Fire Growth and Smoke Transport, CFAST (Version
6), Software Development and Model Evaluation Guide, 2008.
2. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
3. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
4. NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
5. NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE) Volume 3: ThermallyInduced Electrical Failure (THIEF) Model, 2008.
6. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
E-25
E-26
APPENDIX F
F-1
F-2
Figure F-2. Main turbine lubricating oil tanks in the turbine building.
F-3
Member
Maximum Cross-Section
Average Temperature
C (F)
Beam
593 (1,099)
Column
538 (1,000)
Parameter
Effective Fuel Formula
Value
CnH2n+2
Source
Developed from fuel chemistry (n in range of
12-15)
NUREG-1805 Table 3-4
Specified
NUREG-1805 Table 3-4
NUREG-1805 Table 3-4
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
Mulholland and Croarkin (2000)
0.039 kg/s.m2
3,000 L
760 kg/m3
46,000 kJ/kg
2.64 kg/kg
0.059 kg/kg
0.019 kg/kg
0.34
8,700 m2/kg
The peak heat release rate (HRR), , is computed from the fuel mass burning rate,
of combustion, , and the specified area of the spill, :
0.039 kg/m s
46,000 kJ/kg
28.1 m
, the heat
50,400 kW
The fire duration, , is determined from the pool depth, , density, , and burning rate,
F-5
2,144 s
35.7 min
(F-1)
:
(F-2)
F.3.2 Validation
Table F-3 provides a summary of the normalized parameter calculations for the turbine building
fire scenario. A number of parameters fall outside of the parameter space of the NRC/EPRI
V&V study (NUREG-1824 (EPRI 1011999)), including the flame length/ceiling height ratio, the
compartment aspect ratios and some of the target distance to fire diameter ratios. In addition,
the equivalence ratio cannot be calculated with the simple equations that have been used for
various other examples.
The calculation of the equivalence ratio is challenging because natural ventilation is provided
through the 18 roof vents located around the perimeter of the turbine deck level. To evaluate
the potential impact of ventilation on the fire for this scenario, the quantity of oxygen available in
the turbine building is compared to the amount of oxygen that would be consumed by the
specified lubricating oil fire. Given a total volume of approximately 209,600 m3, the mass of
oxygen within the turbine building is estimated to be:
O ,tot
1.1 kg/m
209,600 m
0.23
53,030 kg
(F-3)
The mass of oxygen required to burn all the fuel is estimated to be:
O ,req
50,400 kW
2,144 s
13,100 kJ/kg
F-6
8,249 kg
(F-4)
F-7
F-8
8.5 m
6.0 m
28.0 m
6.0 m
4.7
1.4
1.02
1.1 kg/m
N/A
6.0 m 3.7
0. 52
1.02
4.5
1.2
18.8 m
6.0 m
8.8 m
6.0 m
100.3 m
4.6 m
1.5
3.1
21.8 ;
18.3 m
6.0 m
3.9 m
6.0 m
0.7
3.1
99.5 m
4.6 m
7.2
6.1
11 m
2.2 5.7
0.6 5.7
0.04 0.6
1.2 1.7
0.2 1.0
0.4 2.4
NUREG1824
Validation
Range
Yes/No
No
Yes
N/A
No
Yes
In
Range?
13.5
13.1
0.524
78. m
6.0 m
80. m
6.0 m
9.8 m/s
4 / , where
36.5 m
6.0 m
43.3 m
6.0 m
21.6
3.7
7.2 m
6.0 m
26.9 m
6.0 m
50,400 kW
Notes: The effective diameter of the fire is determined from the formula,
Target
distance to
fire diameter
(Columns
A,B,C,D,E,F)
Compartment
aspect ratios
Equivalence
ratio based
on opening
area
Ceiling jet
radius
relative to the
ceiling height,
Hc
Flame length,
Lf, relative to
ceiling height,
Hc
Fire Froude
Number
Quantity
Table F-3. Normalized Parameter Calculations for the Turbine Building Fire Scenario.
F-9
Unobstructed
flame height
Effective flame
length beneath
ceiling
Figure F-5 shows the approximate radial flame extension of 6.4 m (21 ft) superimposed (yellow
circle) on a detail of Figure F-1 near the two proposed curb locations. As shown in Figure F-5,
Column B is at the edge of the estimated flame extension for Curb Location 1, and Column D is
well within the estimated flame extension for Curb Location 2, which means that the upper part
of Column D would be engulfed in flames. Column A is near the edge of the estimated flame
extension for Curb Location 1, and Column C is near the edge of the estimated flame extension
for Curb Location 2.
F-10
Figure F-5. Detail from Figure F-1 with estimated flame extension beneath ceiling
superimposed.
The heat flux from the fire to a nearby column can be estimated using one of the FDTs
(NUREG-1805, 2005), specifically:
For the point source method, the estimated peak HRR is 50,400 kW, the radiative fraction is
0.33, and the horizontal distance from the center of the lubricating oil pool to the nearest column
(Column D) is approximately 4.2 m (13.8 ft):
0.34
4
50,400 kW
4.2 m
75.0 kW/m
(F-5)
The solid flame calculation in the FDTs is based on the flame being approximated as a
cylindrical, blackbody, homogeneous radiator with an average emissive power (Beyler, SFPE
Handbook, 4th edition). However, Beyler points out that the solid flame calculation represents
the average emissive power over the whole flame and is significantly less than the emissive
power that can be attained locally. Thus, the solid flame calculation is inappropriate for
estimating the heat flux to columns that are potentially engulfed in flames. In fact, neither
calculation is appropriate for estimating the heat flux to columns that are potentially in flames.
Neither of these simple radiation calculations account for the flame extension beneath the
ceiling, so it is likely that they underestimate the actual heat flux at the column, particularly for
the case of the solid flame heat flux calculation.
F-11
/
To calculate the time, crit , when the steel failure temperature is reached, this equation is
rearranged, with the critical steel temperature, crit , inserted for the steel temperature.
/
crit
crit
crit
(F-8)
The term / is sometimes called the section factor, and is the effective thickness of the steel
member; it is calculated as the cross-sectional area of a steel member divided by the heated
perimeter of the member. In the U.S., it is more common to use a parameter referred to as the
W/D ratio, which is simply the section factor multiplied by the steel density. For a W14x145
steel column, the W/D ratio has a value of approximately 96.2 kg/m2 (1.64 lb/ft/in). Using this
value for the W/D ratio, the time to reach the critical steel temperature for the column can be
estimated using the radiant heat flux estimated in Equation F-5:
crit
538 C
36 C
300 s
(F-9)
Because this time is much shorter than the estimated burning duration of the lubricating oil fire
(~2,100 s), this analysis suggests that it may be necessary to protect the columns nearest the
proposed curb locations with fire resistant coatings to prevent the columns from reaching the
specified steel failure temperature.
Because a number of conservative approximations were used for this analysis using algebraic
models, FDS is also used to analyze this scenario.
F-12
Hatchesandstairs
betweenlevels
Upperlevel
HVACroom
Lowerlevel
ColumnE
ColumnB
ColumnA
ColumnD
ColumnC
Plant
north
Ventstoexterior
Curb
Location1
Curb
Location2
Stairopening
Figure F-6. FDS geometry for the turbine building fire scenario.
Materials: The material properties are applied directly as specified to the walls, floor, and
ceiling.
F-13
Model
Bias
Factor,
FDS
FDS
FDS
FDS
FDS
FDS
1.02
1.02
1.02
1.02
1.02
1.02
FDS
FDS
FDS
FDS
FDS
FDS
1.02
1.02
1.02
1.02
1.02
1.02
Standard
Deviation,
Target
Predicted
Value
Critical
Value
Probability of
Exceeding
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.000
0.001
0.828
0.000
0.000
F-14
Temperature of Column B
Temperature of Column A
400
400
350
FDS
300
Temperature (C)
Temperature (C)
350
250
200
150
100
50
FDS
300
250
200
150
100
50
0
0
0
600
1200
1800
2400
3000
3600
600
1200
Time (s)
3000
3600
Temperature of Column D
Temperature of Column C
400
350
350
FDS
300
Temperature (C)
Temperature (C)
2400
Time (s)
400
250
200
150
100
50
FDS
300
250
200
150
100
50
0
0
600
1200
1800
2400
3000
3600
600
1200
Time (s)
1800
2400
3000
3600
Time (s)
Temperature of Column E
Temperature of Column F
400
400
350
350
FDS
300
Temperature (C)
Temperature (C)
1800
250
200
150
100
50
FDS
300
250
200
150
100
50
0
0
600
1200
1800
2400
3000
3600
0
0
600
1200
1800
Time (s)
Time (s)
F-15
2400
3000
3600
Temperature of Column A
Temperature of Column B
800
800
700
FDS
600
Temperature (C)
Temperature (C)
700
500
400
300
200
100
FDS
600
500
400
300
200
100
0
0
0
600
1200
1800
2400
3000
3600
600
1200
Time (s)
3000
3600
Temperature of Column D
Temperature of Column C
800
700
700
FDS
600
Temperature (C)
Temperature (C)
2400
Time (s)
800
500
400
300
200
100
FDS
600
500
400
300
200
100
0
0
600
1200
1800
2400
3000
3600
600
1200
Time (s)
1800
2400
3000
3600
Time (s)
Temperature of Column E
Temperature of Column F
800
800
700
700
FDS
600
Temperature (C)
Temperature (C)
1800
500
400
300
200
100
FDS
600
500
400
300
200
100
0
0
600
1200
1800
2400
3000
3600
0
0
600
1200
1800
Time (s)
F-16
2400
3000
3600
F.6 Conclusion
This analysis has addressed the potential for a relatively large lubricating oil fire to damage
exposed structural steel in a turbine building. The analysis is complicated by the significant
flame impingement on the ceiling caused by an oil fire spread over a relatively large area. This
type of fire behavior is beyond the validation range addressed in NUREG-1824 (EPRI 1011999).
Algebraic calculations were performed to estimate the extent of flame extension beneath the
ceiling. These algebraic calculations indicate that at least one of the columns (Column D) would
be engulfed in the flames extending from the fire at Curb Location 2 These calculations also
indicate that other columns would be located near the outer extent of flames from Curb
Locations 1 and 2. Algebraic calculations were also performed to estimate the time to reach the
critical steel temperature of the nearest column. These calculations indicate that damage could
occur within a time frame of approximately five minutes. These calculations indicate that a more
detailed analysis is warranted. The CFD model, FDS, was used to perform this more detailed
analysis because zone models do not have the necessary physical models to simulate the
postulated fire.
Based on the FDS simulation of this scenario, a 50 MW lubricating oil fire in Curb Location 1 is
not predicted to cause the steel columns to exceed a temperature of 538 C (1,000 F). This is
not the case for the proposed Curb Location 2, which is located closer to Column D.
Consequently, the recommendation for the design package is to install the curbed area at Curb
Location 1.
F.7 References
1. American Institute of Steel Construction, Steel Construction Manual, 13th Edition, New York,
2006.
2. ASTM E119-10a. Standard Test Methods for Fire Tests of Building Construction Materials,
American Society for Testing and Materials, West Conshohocken, PA, 2010.
3. Heskestad, G., and Hamada, T., Ceiling jets of strong fire plumes, Fire Safety Journal, Vol.
21, No. 1, pp. 69-82, 1993.
4. Karlsson, B., and J.G. Quintiere, Enclosure Fire Dynamics, CRC Press LLC, Boca Raton,
FL, 2000.
5. Mulholland, G.W., and C. Croarkin, Specific Extinction Coefficient of Flame Generated
Smoke, Fire and Materials, 24:227-230, 2000.
6. NIST NCSTAR 1-5F, Federal Building and Fire Safety Investigation of the World Trade
Center Disaster: Computer Simulation of the Fires in the World Trade Center Towers, 2005.
7. NIST SP 1018-5, Fire Dynamics Simulator (Version 5), Technical Reference Guide, Vol. 3,
Experimental Validation, 2010.
8. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
9. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
10. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
F-17
F-18
APPENDIX G
G-1
G-2
G-3
Parameter
Effective Fuel Formula
Peak HRR
Time to reach peak HRR
Heat of Combustion
Heat of Combustion per unit
mass of oxygen consumed
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Value
C6H10O5
2,500 kW
420 s
17,100 kJ/kg
13,100 kJ/kg
Source
Cellulose
SFPE Handbook, 4th Ed., Figs. 3-1.65, 3-1.100
SFPE Handbook, 4th Ed., Figs. 3-1.64
SFPE Handbook, 4th Ed., Table 3-4.16
Hugget 1980, Average value
1.27 kg/kg
0.015 kg/kg
0.004 kg/kg
0.37
HRR (kW)
2500
2000
1500
1000
500
0
0
600
1200
1800
2400
3000
Time (s)
G-4
3600
G.3.2 Validation
A source of validation data justifying the use of the fire models discussed above for this scenario
is the U.S. Nuclear Regulatory Commission/Electric Power Research Institute (NRC/EPRI)
verification and validation (V&V) study documented in NUREG-1824 (EPRI 1011999). The
National Bureau of Standards (NBS, now the National Institute of Standards and Technology
(NIST)) Multi-Room Test Series consists of 45 fire experiments in two relatively small rooms
connected via a relatively long corridor. The fire source, a gas burner, was located against the
rear wall of one of the small compartments, and fire tests of 100 kW, 300 kW, and 500 kW were
conducted. The present scenario has a larger fire in a larger, longer compartment. The V&V
study assessed the MQH HGL temperature correlation, the Alpert ceiling jet correlation, and the
zone model, MAGIC. The unconfined ceiling jet correlation of Delichatsios is included in this
analysis to assess the sensitivity of various input parameters.
Table G-2 lists various important model parameters and the ranges for which the NRC/EPRI
validation study is applicable. The calculations in Table G-2 are for the fire room
(compartment 8, as shown in Figure G-4). The ceiling jet radial distance relative to the ceiling
height is not within the range of validation. As a sensitivity case, an additional detector was
added at a distance that would fall within the validation range (6.8 m (22.3 ft)). In addition, the
room-width geometry is not within the range of validation. A sensitivity case with the height
reduced to a value within the validation range (5 m (16.4 ft)) was performed. The remaining
model parameters are within the applicable validation ranges.
G-5
Table G-2. Normalized parameter calculations for the multi-compartment corridor fire scenario.
Quantity
Fire Froude
Number
2500 kW
1.2 kg/m
Flame Length,
L , relative to
the Ceiling
Height, Hc
3.7
1.02
Ceiling Jet
Horizontal
Radial
Distance,rcj ,
relative to the
Ceiling Height,
Hc
Equivalence
Ratio, , as an
indicator of the
Ventilation Rate
Compartment
Aspect Ratios
Target
Distance, r,
relative to the
Fire Diameter,
D
0.23
9.8 m/s
0.7
1.29 m 3.7
1.2
1.02
1.2 kg/m
1.67 m /s
3.0 m
6.1 m
2.49 ;
N/A
In
Range?
0.4 2.4
Yes
0.2 1.0
Yes
1.2 1.7
No
0.04 0.6
Yes
0.6 5.7
No
2.2 5.7
N/A
1.2
3.8 m
0.8
2500 kW
13,100 kJ/kg 0.46 kg/s
0.23
15.2 m
6.1 m
0.44 m 3.8 m
6.1 m
4.46 m
6.1 m 0.44 m
NUREG1824
Validation
Range
0.4
0.46 kg/s
0.49
Notes:
(1)
(2)
(3)
G-6
G-7
5
4
7
1
Figure G-4. Effective corridor layout for implementation in zone models (not to scale).
Length
m (ft)
Width
m (ft)
Area
m2 (ft2)
8.1 (26.6)
4.1 (13.5)
33.2 (357)
2.0 (6.6)
23.4 (76.8)
46.8 (504)
45.1 (148)
4.1 (13.5)
184.9 (1990)
8.1 (26.6)
6.0 (19.7)
48.6 (523)
10.3 (33.8)
6.6 (21.7)
68.0 (732)
10.3 (33.8)
6.6 (21.7)
68.0 (732)
12.2 (40)
8.2 (26.9)
100.0 (1076)
3 (9.8)
15.2 (49.9)
45.6 (491)
G-8
Fire: The fire is modeled as a 1.3 m2 (14 ft2) source (equivalent diameter of 1.29 m (4.2 ft)) at
an elevation of 0.44 m (1.4 ft) (see Figure G-2). The stoichiometric oxygen to fuel ratio, , is
calculated using Equation 22 from Chapter 4 in Section 3 of the SFPE Handbook, 4th edition
and the values from Table G-1, as follows:
17.1 kJ/g
13.1 kJ/g
1.3 g/g
(G-1)
where is heat of combustion per unit mass of fuel consumed and O is the heat of
combustion per unit mass of oxygen consumed. One of the inputs required by MAGIC is the
specific area, (NUREG-1824, Volume 6, Section 3.2.7):
0.015
7,600 m
g
114 m /g
(G-2)
is a constant and is the soot yield, as listed in Table G-1. Products of combustion
where
for a wood pallet fire.. The pyrolysis rate (g/s) is calculated for input to MAGIC by dividing the
HRR values (kW) at each time step (as shown in Figure G-3) by the heat of combustion
G-9
G-10
G-11
G-12
Model
Bias
Factor,
Standard
Deviation,
Ventilation
Predicted
Value
Critical
Value
Probability
of
Exceeding
1.56
0.32
Natural
256
330
0.001
1.01
0.07
Mechanical
240
330
0.000
G-13
300
200
MAGIC
160
FIVE (MQH)
Temperature (C)
Temperature (C)
250
200
150
100
50
120
80
40
0
0
600
1200
1800
2400
Time (s)
600
1200
1800
2400
Time (s)
Figure G-9. HGL temperature predictions by MAGIC for the corridor scenario.
Chapter 4, Model Uncertainty, provides guidance on how to express the uncertainty of the
MAGIC predictions. MAGIC predicted a temperature rise of M = 240 20 = 220 C (428 F).
As shown in Table 4-1, the NRC/EPRI V&V study (NUREG-1824 (EPRI 1011999)) found that
MAGIC predictions of HGL temperature rise are, on average, 1% greater than corresponding
measurements, with a relative standard deviation of 7%. The adjusted model prediction is =
220/1.01 = 218, and the standard deviation is = 0.07 (220/1.01) = 15.2 C. Therefore, the
probability that the cable temperature would exceed 330 C (625 F) is:
330
1
330 20 218
erfc
2
15.22
(G-3)
The same uncertainty calculation for the MQH correlation in FIVE-Rev1 results in a probability
of exceeding the damage temperature of 0.001, or 0.1%. In other words, there is a near-zero
probability of exceeding the damage temperature threshold for cables within the compartment of
fire origin based on a surrounding HGL temperature, according to the FIVE-Rev1 and MAGIC
predictions. This demonstrates that detailed analyses of the cables outside the compartment of
origin are not warranted.
Sensitivity of the Room Width to Height Ratio
As shown in Table G-2, the room width to height ratio is not within the validation bounds. A
sensitivity case is run in which the room height is reduced until the ratio is within the bounds. All
inputs are the same, except that the height is reduced from 6.1 m (20 ft) to 5 m (16.4 ft) in the
fire room. The results, plotted in Figure G-10, show that although the temperatures are slightly
higher in the fire room, there is no change in the adjacent room, which contains the target. This
shows that the room width to height ratio is not a significant influence on the results in this
scenario and that the MAGIC results are applicable, even though the ratio falls outside of the
validation bounds.
G-14
250
Base Case
(H=6.1 m)
200
Ceiling
Height=5 m
Temperature (C)
Temperature (C)
300
150
100
160
Base Case
(H=6.1)
120
Height=5 m
Fire Room
50
0
0
600
1200
1800
80
40
0
2400
600
1200
Time (s)
1800
2400
Time (s)
Temperature (C)
80
FIVE (Alpert)
70
FIVE (Delichatsios)
60
MAGIC
50
40
30
20
10
0
0
50
Time (s)
100
150
G-15
70
Radial distance=6.8 m
Temperature (C)
60
50
40
30
20
10
0
0
50
Time (s)
100
150
Figure G-12. Detector temperature for two radial distances predicted by MAGIC.
G.6 Conclusion
This analysis addresses a transient fire in a relatively long, narrow room, identified as a corridor,
which is connected with other rooms in a multi-room arrangement. The targets include cables in
rooms located beyond the corridor where the fire is specified to occur. The MQH correlation
included in FIVE-Rev1 was used to estimate temperature conditions within the corridor, and it
was found that the HGL temperature predicted by the MQH correlation was lower than the cable
damage temperature. Because the cables are located in rooms beyond the corridor where the
fire occurs and it can be reasoned that the temperatures in these rooms will be lower than in the
corridor, cable damage would not be expected.
The zone model MAGIC was used to predict the HGL temperatures in all of the interconnected
compartments, resulting from the stack of burning pallets and trash bags. These MAGIC
calculations also demonstrate that the HGL temperature is not high enough to cause cable
damage in any compartment or corridor, including the corridor of fire origin. These predictions
account for model uncertainty in the temperature predictions and the sensitivity to variations in
the heat release rate. Based on a simplified method for smoke detector activation, smoke
detector operation is expected to occur between 40 s and 50 s after the ignition of the fire.
G-16
G.7 References
1. Gay, L., C. Epiard, and B. Gautier, MAGIC Software Version 4.1.1: Mathematical Model,
EdF HI82/04/024/B, Electricit de France, France, November 2005.
2. Huggett, C., Estimation of Rate of Heat Release by Means of Oxygen Consumption
Measurements, Fire and Materials, 4:61-65, 1980.
3. Karlsson, B., and Quintiere, J. G., Enclosure Fire Dynamics, CRC Press, Boca Raton, 2000.
4. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
5. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
6. NUREG/CR-6850 (EPRI 1011999), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
7. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
8. UL 217, Single and Multiple Station Smoke Alarms, 6th edition, Underwriters Laboratories,
Inc., 2006.
G-17
G-18
APPENDIX H
H-1
Parameter
Heat of Combustion
CO2 Yield
Soot Yield
CO Yield
Radiative Fraction
Value
20,900 kJ/kg
1.29 kg/kg
0.136 kg/kg
0.147 kg/kg
0.49
Source
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
SFPE Handbook, 4th ed., Table 3-4.16
Local burnout of the fire occurs when the cable plastic is consumed. The time to burnout is
:
calculated as follows. First, determine the combustible mass per unit area of tray,
1
120
0.33
1 0
0.6 m
0.4 kg/m
26.4 kg/m
(H-1)
5 avg /6
26.4 kg/m
20,900 kJ/kg
5/6 250 kW/m
2648 s
(H-2)
where is the heat of combustion and avg is the average HRR per unit area of tray. The HRR
per unit area ramps linearly to its average value over a time period of /6, remains steady for a
time period of 2 /3, and then decreases linearly to zero over a time period of /6. The linear
ramp-up and ramp-down are typical ways of approximating the time history of an items HRR.
H-2
H-3
Figure H-2. Geometry details of redundant cable trays located in the annulus.
H-4
HRR (kW)
600
1200
1800
2400
3000
3600
Time (s)
Suppression System: Standard response sprinklers are located on the inner wall, as shown in
Figure H-2. The sprinklers have a response time index (RTI) of 130 (ms)1/2 and activate at a
temperature of 100 C (212 F) (NUREG-1805, Chap. 10). Each sprinkler is topped by heat
collectors designed to trap heat from a fire. The ambient temperature within the annulus region
is typically 35 C (95 F).
H-5
H.3.3 Validation
Only one of the non-dimensional parameters that have been used to characterize the fire
scenarios is applicable here (Table H-2), mainly because the other parameters address
phenomena unique to compartment fires. The only relevant parameter, / , indicates the
relative distance separating the fire from the target cables. However, the effective fire diameter,
, is not well-defined when the fire is expected to spread vertically and horizontally within the
cable trays. The lateral distance between the burning tray and the target tray is 1.7 m (5.6 ft).
The trays are each 0.6 m (2 ft) wide. For International Collaborative Fire Model Project (ICFMP)
Benchmark Exercise #3, which is described in NUREG-1824, the relative distance between the
fire and the heat flux gauges in some of the tests was similar to this scenario. However, the
accuracy of the solid flame and point source algorithms included in the FDTs varied
considerably in the validation study. The point source heat flux model loses accuracy when the
target is relatively close to the fire. The solid flame model is designed to improve the accuracy
in the near-field. However, the implementation within the FDTs is based on the fire being
modeled as a pool fire with a nearly circular base and cylindrical shape. This is not the case in
the scenario under consideration. In short, the simple heat flux models included in the FDTs
have been assessed in NUREG-1824 for experiments of comparable scope to the given
scenario, but the models have not been demonstrated to be particularly accurate. For this
reason, the models will be used to provide screening estimates, but the predictions will be
supplemented by CFD calculations.
The FDS sprinkler activation algorithm was validated using a variety of experimental test series
(NIST SP 1018-5). The plume algorithm was also assessed in NUREG-1824. The cable failure
algorithm, THIEF, was developed and validated in NUREG/CR-6931 (Vol. 3).
H-6
NUREG-1824
Validation
Range
In
Range?
Fire Froude
Number
N/A The fire does not conform to classic fire plume theory.
0.4 2.4
N/A
N/A The fire does not conform to classic fire plume theory.
0.2 1.0
N/A
1.2 1.7
N/A
Equivalence Ratio,
, as an indicator
of the Ventilation
Rate
0.04 0.6
N/A
Compartment
Aspect Ratio
0.6 5.7
N/A
Target Distance, r,
relative to the Fire
Diameter, D
2.2 5.7
Yes
Fire Height,
, relative
to the Ceiling
Height,
For the point source method, the estimated peak HRR, , is 945 kW, the radiative fraction, ,
is 0.49, and the horizontal distance from the center of the burning tray to the edge of the target
tray, , is 2 m (6.6 ft). The calculated heat flux is:
ps
0.49
4
945 kW
2.0 m
9.2 kW/m
(H-3)
The solid flame calculation is based on the fire having a roughly circular base with which it
employs a flame height correlation to estimate the vertical extent of the luminous flame region.
H-7
0.49 255
4
2.9
172.5
2.4
172.5
2.0
172.5
2.2
172.5 kW
2.9
m
6.2 kW/m
(H-4)
Much like the solid flame method, the distributed point source method provides a refined
estimate of the heat flux, based on the fact that the fire is distributed over the two parts of the
tray and is not concentrated at a single point. The end result is a lower estimate of the heat flux
to a given target point. The fire could be distributed over more than five points, but the answer
would not significantly change. Neither the point source nor the distributed point source method
accounts for changes in the orientation of the target. The target cables are modeled as having
a direct view of the fire with no obstructions, such as the sides of the trays, that can reduce the
heat flux to the target cables. The CFD model can provide a more refined estimate.
H-8
H-9
H-10
Model
Bias
Factor,
Standard
Deviation,
Predicted
Value
Critical
Value
Probability
of
Exceeding
1.42
9.2
6.0
0.553
1.42
0.55
6.2
6.0
0.248
6.0
0.000
205.0
0.000
100.0
0.001
Point Source
Distributed Point
Source
FDS
1.10
FDS
1.02
FDS
1.15
Target
Cables
0.17
2.5
Target Temperature (C)
0.13
Cables
120.0
Plume Temperature (C)
0.11
Sprinkler
90.0
H-11
150
120
Temperature (C)
90
FDS
60
7
6
Point Source
Distributed Point Source
FDS
5
4
3
2
30
1
0
0
0
600
1200
1800
2400
3000
3600
600
1200
1800
2400
3000
3600
Time (s)
Time (s)
Temperature (C)
100
80
60
40
20
0
0
600
1200
1800
2400
3000
3600
Time (s)
Figure H-6. Predicted sprinkler link temperature for the annulus fire scenario.
It is possible to determine how large a fire is needed to increase the plume temperature by
10 C. Table 4-3 indicates that the HGL temperature rise is proportional to the HRR to the 2/3
power. The plume temperature behaves similarly. Following the methodology in Section 4.4, in
H-12
3
2
3
945 kW
2
10 C
90 C 35 C
258 kW
(H-5)
In other words, the peak HRR of the fire would have to be approximately 945 + 258 = 1203 kW
to bring the plume temperature into a range to cause sprinkler activation.
H.6 Conclusion
Simple point source heat flux calculations indicate that a fire in one of the cable trays within the
annulus region of the containment building might damage the cables in an adjacent tray.
However, an additional analysis using FDS indicates that cable damage is unlikely due to the
orientation of the target cables and the blockage of thermal radiation by the cable tray itself.
This suggests that the details of the cable tray location, orientation, and configuration can
significantly impact potential for damage.
FDS predicts that sprinkler activation above the fire is unlikely. However, its prediction is
sensitive to the exact location of the sprinkler relative to a fire plume that may be subject to
unpredictable air movements throughout the entire facility. Alternative protection strategies,
such as shielding between trays or other thermal barriers, should be considered to ensure the
protection of the redundant cables.
H.7 References
1. NIST SP 1018-5, Fire Dynamics Simulator (Version 5), Technical Reference Guide, Volume
3, Experimental Validation, 2010.
2. NUREG-1805, Fire Dynamics Tools (FDTs) Quantitative Fire Hazard Analysis Methods for
the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program, 2004.
3. NUREG-1824 (EPRI 1011999), Verification and Validation of Selected Fire Models for
Nuclear Power Plant Applications, 2007.
4. NUREG/CR-6850 (EPRI 1011989), EPRI/NRC-RES Fire PRA Methodology for Nuclear
Power Facilities, 2005.
5. NUREG/CR-6931, Cable Response to Live Fire (CAROLFIRE), Volume 3: ThermallyInduced Electrical Failure (THIEF) Model, 2008.
6. NUREG/CR-7010, Cable Heat Release, Ignition, and Spread In Tray Installations during
Fire (CHRISTIFIRE), Volume 1: Horizontal Trays, 2012.
7. SFPE Handbook of Fire Protection Engineering, 4th edition, 2008.
H-13
H-14
INDEX
Administrative controls ............................................................................................................... 2-7
Algebraic models .................................... 1-7, 1-15, 1-16, 2-10, 2-11, 2-13, 3-18, 3-21, 3-23, 4-17
ASTM E1355...................................................................................................................... 1-4, 5-1
Ceiling jet ................................................................................ 2-24, 2-30, A-9, B-8, E-8, F-8, F-17
CFAST .xvii, xviii, xxv, 1-3, 1-7, 1-15, 1-16, 2-11, 2-13, 2-18, 2-19, 2-23, 2-32, 2-33, 3-4, 3-21, 42, 4-3, 5-2, 5-3, 5-4, 6-2, A-6, A-7, A-11, A-13, A-14, A-15, A-18, A-19, A-20, A-21, A-22, B-6,
B-11, B-12, B-13, B-14, B-15, B-16, B-19, B-20, B-21, B-24, B-25, D-5, D-8, D-9, D-10, D-11,
D-12, D-14, D-15, D-16, D-17, D-18, E-6, E-11, E-12, E-13, E-14, E-15, E-16, E-17, E-19, E20, E-21, E-24, E-25, F-9
CFD model ... 1-9, 1-10, 2-11, 2-19, 2-20, 2-22, 2-30, 2-33, 3-2, 3-4, 3-10, 3-12, 3-13, 3-15, 3-16,
3-17, 3-19, 3-21, 3-23, 4-11, A-7, A-12, B-6, C-5, C-7, C-16, C-23, D-8, E-6, E-20, F-6, F-17,
H-8
Compartment aspect ratio................................................................................................ 2-30, F-8
Compartment geometry ................................................................................................... 3-2, A-13
Detection ..................................................................... 3-20, 5-4, E-2, E-16, E-18, E-20, F-4, G-15
Documentation ...................................................................................................... 2-14, 2-18, 2-34
Equivalence ratio ............................................................................................................. 2-30, F-8
Fire Dynamics Simulator (FDS) . xvii, 1-3, 1-11, 2-11, 2-19, 4-3, 5-3, A-6, B-6, C-6, D-5, E-6, H-6
Fire Froude number ................................................................................................................. 2-25
Fire location ............................................................................................................................... 3-2
Fire PRA ... xvii, xix, xxiii, 1-1, 1-3, 1-14, 1-15, 2-5, 2-10, 4-10, 4-11, 5-4, A-22, B-24, C-24, D-17,
E-25, G-17, H-13
Heat release rate ................................................................................................. 2-17, B-20, C-10
Hot Gas Layer (HGL) ....................................................................................................... 2-6, 2-23
Intervening combustibles ........................................................................................................... 2-9
Limiting Fire Scenario (LFS) .................................................................................................... 4-15
MAGIC .. xvii, xxiii, 1-3, 1-7, 1-16, 2-11, 2-13, 2-19, 2-23, 2-33, 3-21, 4-2, 4-3, 4-9, 4-10, 5-1, 5-3,
A-6, C-3, C-5, C-7, C-8, C-9, C-10, C-11, C-12, C-13, C-14, C-15, C-16, C-17, C-18, C-19, C20, C-21, C-23, C-24, D-5, G-5, G-7, G-9, G-10, G-11, G-12, G-13, G-14, G-15, G-16, G-17
Main control room .................................................................................................. xxvi, 1-15, 3-18
Maximum Expected Fire Scenario (MEFS) .............................................................................. 4-15
Model uncertainty....................................................................................................................... 4-1
NFPA 805 ................................................ iii, 1-1, 1-11, 1-12, 1-13, 1-14, 2-4, 2-22, 4-2, 4-15, 5-2
Normalized parameter .......................................................2-22, A-9, B-8, C-6, D-7, E-8, G-6, H-7
Parameter uncertainty ...................................................................................................... 4-1, 4-10
Plume .......................................................... 2-15, 2-17, 2-23, 3-8, 4-2, 4-7, B-10, E-7, E-10, H-11
Secondary combustibles ............................................................................................................ 2-9
Sensitivity analysis ................................................................................................................... 4-10
Smoke production ...................................................................................................................... 1-8
Smokeview.1-10, 1-11, 2-19, 5-1, A-15, A-16, B-12, B-18, C-14, D-9, D-12, D-13, E-12, E-17, H8
Structural steel .......................................................................................................................... F-4
Suppression ...................................................................................................................... E-2, H-5
Index-1
INDEX
Target. 1-6, 1-7, 2-8, 2-15, 2-23, 3-12, 3-13, 4-2, 4-7, 4-8, 4-14, A-8, A-9, B-8, C-6, D-7, D-14, E8, F-6, F-8, F-14, G-6, H-7, H-11
Uncertainty ................................................................... 2-33, 4-1, 4-6, 4-8, 4-10, B-22, E-23, G-14
Ventilationxxv, 1-14, 2-8, 2-14, 3-1, 3-3, 3-4, 4-14, A-1, A-9, A-14, A-16, A-18, A-22, B-1, B-8, B10, B-11, B-14, B-17, C-3, C-13, C-16, D-2, D-7, D-11, D-13, E-2, E-8, E-10, E-15, E-18, F-4,
F-14, G-1, G-6, G-11, G-13, H-7
Zone model ............................... 1-8, 2-11, 3-10, 3-13, 3-15, 3-17, 3-19, 3-21, 3-23, A-7, D-6, G-5
Index-2
(12-2010)
NRCMD 3.7
1. REPORT NUMBER
NUREG-1934
Nuclear Power Plant Fire Modeling Analysis Guidelines (NPP FIRE MAG)
YEAR
November
2012
5. AUTHOR(S)
6. TYPE OF REPORT
Technical
7. PERIOD COVERED (Inclusive Dates)
8. PERFORMING ORGANIZATION - NAME AND ADDRESS (If NRC, provide Division, Office or Region, U. S. Nuclear Regulatory Commission, and mailing address; if
contractor, provide name and mailing address.)
U.S. Nuclear Regulatory Commission, Office of Nuclear Regulatory Research, Washington, DC 20555-0001
Electric Power Research Institute, 3420 Hillview Avenue, Palo Alto, CA 94303
National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, MD 20899
Kleinsorg Group Risk Services, a Division of Hughes Associates, 3610 Commerce Drive, Baltimore, MD 21227
9. SPONSORING ORGANIZATION - NAME AND ADDRESS (If NRC, type "Same as above", if contractor, provide NRC Division, Office or Region, U. S. Nuclear Regulatory
Commission, and mailing address.)
There is a movement to introduce risk-informed and performance-based (RI/PB) analyses into fire protection engineering
practice, both domestically and worldwide. This movement exists in both the general fire protection and the nuclear
power plant (NPP) fire protection communities. The Nuclear Regulatory Commission (NRC) has used risk-informed
insights as a part of its regulatory decision making since the 1990s. In 2001, the National Fire Protection Association
(NFPA) issued NFPA 805, Performance-Based Standard for Fire Protection for Light-Water Reactor Electric Generating
Plants, 2001 Edition. In July 2004, the NRC amended its fire protection requirements in Title 10, Section 50.48 of the
Code of Federal Regulations to permit existing reactor licensees to voluntarily adopt requirements contained in
NFPA 805 as an alternative to the existing deterministic fire protection requirements. One key element in RI/PB fire
protection is the availability of verified and validated (V&V) fire models that can reliably estimate the effects of fires.
The U.S. NRC, together with the Electric Power Research Institute (EPRI) and the National Institute of Standards and
Technology (NIST), conducted a research project to verify and validate five fire models that have been used for NPP
applications. The results of this effort are documented in the seven-volumes of NUREG-1824 (EPRI 1011999),
Verification & Validation of Selected Fire Models for Nuclear Power Plant Applications. NUREG-1934 (EPRI 1023259)
describes the implications of the V&V results for fire model users. The features, limitations, and uses of the five fire
models, documented in NUREG-1824, are discussed relative to NPP fire hazard analysis applications.
12. KEY WORDS/DESCRIPTORS (List words or phrases that will assist researchers in locating the report.)
Fire, Performance-Based, Verification and Validation (V&V), Fire Hazard Analysis (FHA),
Risk-Informed Regulation, Fire Safety, Fire Protection, Nuclear Power Plant (NPP),
Probabilistic Risk Assessment (PRA), Fire Modeling
unlimited
14. SECURITY CLASSIFICATION
(This Page)
unclassified
(This Report)
unclassified
15. NUMBER OF PAGES
16. PRICE
NRC FORM 335 (12-2010)
UNITED STATES
NUREG-1934
Nuclear Power Plant Fire Modeling Analysis Guidelines (NPP FIRE MAG)
November 2012