FLASHCAT for tray failures
FLASH-CAT Method
The initial fire source may ignite electrical raceways and cables in the compartment, which are termed as secondary combustibles in this report. Despite the term “secondary,” these objects may not combust simultaneously but rather at different timings. The secondary fires can add to the heat released by the initial fires, which eventually changes the fire and temperature calculations. In order to estimate the total heat from the initial and secondary fires, the FLASH-CAT empirical methodology developed by the National Institute of Standards and Technology (NIST) [3] is adopted in this work. Circled area 1 in Figure 3 shows the application of FLASH-CAT in the scenario generation flow.
Calculation Process The FLASH-CAT methodology estimates the HRR from the burning cables, the timing and duration of fire, and the length of cable burned during the initial combustion. These values are approximated from controlled fire experiments conducted by NIST. For that reason, this methodology requires certain assumptions that fit the experiment conditions to be met. The FLASH-CAT assumptions are: The cable trays are horizontal and stacked vertically with a spacing of less than 0.45 m (18 in) The cables burn in the open; that is, they are away from walls and well below the ceiling The cables are not exposed to elevated temperature sources except for the ignition source below There are no barriers separating the trays, and the tray tops and bottoms are open The cables are not protected with coatings, armor shielding, or thermal blankets of any kind There is a fire beneath the lowest tray The initial extent of the fire in the lowest tray is equal to the width of the source fire Each tray has at least a single row of cables or roughly 25 % of the NRC limit. Under these assumptions, the fire is assumed to propagate upwards through the array of cable trays according to an empirically determined timing sequence. The length of cables within a given tray that ignite initially increases as the fire spreads upwards. Lateral spread of the fire begins as soon as the cables within the tray ignite. This produces a solid V-shaped burning pattern that expands laterally with time. As the mass of combustible material within the center of the V is consumed, the V-shape becomes an expanding, open wedge of burning cable. The fires in each tray continue to spread until the end of the tray is reached. This process is illustrated in Figure 4.
Figure 4. V-Shaped burning pattern in FLASH-CAT model. The initial length of fire on the raceways is determined by the raceway’s vertical distance from the fire, and the order of the raceway among other raceways. For the first raceway (i.e., the raceway closest to the fire), the initial length is the same with the length of the fire itself. Once this first raceway combusts, the other raceways on top of it will eventually burn. The initial length of fire on those other raceways are calculated as:
Where L(i) denotes the fire length of the current raceway, L(i+1) is the fire length for the next raceway in the vertical order, and the hi is the vertical distance between the two raceways. The fire spreads laterally with a rate of 1.1 m/hour for thermoset cables and 3.2 m/hour for thermoplast cables. If there are multiple cable types in the same raceway, the fire spread rate is assumed to be the fastest among the cables (i.e., thermoplast). Cables consist of metal and plastic elements. The amount of combustible mass in a raceway depends on the plastic mass of cables in it. A variable called the combustible mass per unit area mc” is calculated as: m_c^"=(∑_j▒(Y_p (1-v)m^' )_j )/W (2)
Where j is the index of a cable in the raceway, Yp is the plastic mass fraction of the j-th cable, v is the char yield of the cable, m’ is the mass per unit length of the cable, and W is the tray’s width. The amount of combustible mass in a raceway determines how long the fire burns, which is denoted as Δt: Δt=(m_c^" ΔH)/(5q ̇_avg^"/6) (3)
Where ΔH is the heat of combustion, and q ̇_avg^" is the HRR per unit area which is 150 kW and 250 kW for thermoset and thermoplast cables respectively. To calculate the time-dependent HRR profile, the raceway is divided into spatial grids, and the distance from the center of the fire is measured as x. The time when a spatial grid begins to combust, tign,i(x), is calculated as: t_(ign,i) (x)=t_(ign,i,0)+max(0,(|x|-L_i/2)/S ) (4)
Where S is the fire spread rate in m/s, and t_(ign,i,0) is the time when the tray first ignites. If the grid is located within the original section of raceway that first ignites, the two ignition times will be the same. After calculating the combustible mass m_c^", the fire burn period Δt, and the ignition time t_(ign,i) (x), each of the raceway grid will have a local HRR per unit area as shown in the Figure 5. HRR starts to increase from t_(ign,i) (x) until reaching the peak Heat Release Rate Per Unit Area (HRRPUA) at one-sixth of Δt later, remains stable until 5Δt/6, and then decreases linearly to zero.
Figure 5. Time history of the local HRRPUA for each raceway grid. Once all the HRRPUA profiles for each raceway grid have been calculated, the total HRR of the secondary combustible is calculated by summing all the local HRRPUA profiles from all raceways as follows: Q ̇(t)=Q ̇_fire+W∑(i=1)^(N_trays)▒[∫(x_0)^(x_j)▒〖q ̇"(t-t_(ign,i) (x))dx 〗] (5)
Where Q ̇(t) is the total secondary combustible’s HRR, Q ̇_fire is the HRR of the initial fire, W is the raceway’s width, q ̇" is the local HRRPUA, and t is the discrete observation time. The process to initiate FLASH-CAT in our work is as follows: Model the compartment, fire source(s), vent(s), cables, and raceways in FRI3D. Run the initial CFAST simulation.
Get CFAST output of actual HRR, fire area, heat flux, and temperature at each raceway. Iterate through each raceway. Check raceway geometrical coordinate if it is located above the fire source. (If a raceway is above it is conservatively included vs. the 0.45 m spacing requirement). If raceway is on top of a fire source, calculate the initial length of cable that combusts. Iterate on each cable in the raceway. Check if the maximum temperature or heat flux at the raceway exceeds the cable’s combustion threshold. If it does, check the cable material properties. If the properties are not defined by the modeler in FRI3D, the default material properties are used for the cable, typically thermoplastic unless changed by the user (see Appendix B). Create a one-dimensional spatial mesh along the cable’s length (dx) and calculate which mesh segments are burning at certain times (dt) by taking into account the horizontal fire spread rate. Integrate the spatio-temporal heat release rate over the cable’s length to obtain the time-dependent HRR for that cable. Sum the time-dependent HRR for all the cables within the raceway. Repeat Step 5 until Step 7 for other raceways in the vertical stack. Add the cable HRR to the initial fire’s HRR, update the CFAST input file with the new HRR, and re-run CFAST simulation. Use the new CFAST results to get the temperature and heat flux at the next raceway in the iteration. Repeat from Step 4. After all raceways have been evaluated, get the temperature and heat flux from the last CFAST simulation to estimate component failures. Data from the 3D model When modeling raceways, a single raceway is typically broken up into sections wherever a cable enters or exits, which enables simple linking to the FRANX style of logic modeling. However, to calculate FLASH-CAT, the entire tray data and all cables need to be used. To enable this, each raceway has links to other raceway pieces. When the simulation indicates a particular raceway section will have secondary combustibles, the data from all linked raceways is used in the FLASH-CAT calculation. FRI3D automatically adds the links to raceways if they are touching and distinguishes between horizontal and up or down connections. Currently only horizontal connections are calculated for secondary combustion. The linking of local raceways described above is the solution chosen to ensure FLASH-CAT fire can propagate along the whole physical raceway. From the initial fire source, data for each route or path away from the fire is gathered, branches in the raceway cause additional fire propagation paths. The raceways sections that are linked programmatically typically each contain different cables than adjacent raceways. For more realism, these differences should be accounted for to calculate the FLASH-CAT HRR. However, in the current version, a conservative simplified approach is used. Instead of calculating the HRR contribution for each raceway section, the conservative method uses the section with the largest combustible mass and HRR per unit area for the entire raceway. This could be updated in a future revision but was done currently to reduce code complexity and computation time. The combustible mass is calculated for each local raceway following Equation (2). The maximum value among all local raceways is taken as the combustible mass for all pieces of raceway in the FLASH-CAT calculation. This value is then used to calculate the burning period Δt following Equation (3). Therefore, the burning period for each raceway grid is conservatively longer than what is expected. These local HRRPUA profiles are then summed following Equation (5) to form the overall HRR profile. As a result, the HRR profile has conservative higher gradient and the fire burns longer.