u: s .3. 9 23‘th .1 v 1.13 251.25. ! ‘ .7: I, 19!; . . H.371 913...... a .431 .n. I. 3.. 5 n.0, ‘ . X mam. y Xo-n‘lu-Q - . . a .K: . r E h. . ‘ ‘ V . ‘ 11%... If... 97:! , 3 i. . , . 3.. $7 . .13.. . . 5.22qu f... . .11.. I o p ,i. : may! . ‘ 0 I .ll,“. II a, fin:fiuwnn.flotzu: fitrTlMFf : (Hanan: ngfia. i. . . a 1:. . . cc- 16 .- .agt: . p. 60.1!“ 1 .3 «twin». 1...... L c L 1:5..hihxfig a v Pu.» it. l 1.". fifty! t. .52.. .31.... *5». E it. A - 3 i ; 2... «sq... ithfl. 3:2, ‘ I’d-v.5 In: A is: . , 2 ".943. . s 1 >§zfli§§ Eng. J '(wl-.$ r. x. I at ‘ ‘25» hd....-&émfif. . 3%. i 1:}: N :75 03:93.; I03...)- 1...! . an 3...... in. 3§ ‘w -—-> Differential Pressure Flow Manometer 90° Direction Junction Figure 2.3 Blower and elbow meter (including the differential pressure manometer). 20 In order to know what the volume flow rate of air is (for the low volume flow rates examined in this study) the inlet cross-section must be constricted down to a rather small area, thus negating the advantage of downstream induction altogether. Air is injected into the settling chamber using a Dayton 1/2 hp blower (Dayton PSC 4C831A). The flow rate is regulated with a sliding gate at the blower inlet. The blower is attached to 7.6cm diameter PVC pipe that directs flow into the settling chamber. In order to measure the flow rate, an elbow meter was constructed. The pipe leading to the settling chamber is fitted with a 90° junction. Pressure taps were drilled into the inner and outer bends of the 90° junction. Tygon tubing was extended from these pressure taps to a differential pressure manometer (MKS Baritron model #225AD-00010AABS, capable of measuring pressure differences from 0 to 1244 Pa, calibrated to in-HZO). The Baritron pressure transducer is powered by a Power-One linear power supply (Power-One, HB28-1-A , 28 volts at 1 amp). The pressure transducer, 90° junction, pressure taps, etc. are collectively known as the elbow meter. The elbow meter was calibrated using a correlation between the pressure difference displayed by the output of the manometer and the velocity in the combustion chamber. This correlation (Fig. 2.4) was developed over a range of velocities by adjusting the damper plate at the blower inlet. The fin (discussed in Chapter 3) was attached to the trailing edge of the sample holder). 21 Velocities were measured in the center of the combustion chamber using a TSI IFA 300 Constant Temperature Anemometer. A least squares linear regression analysis was performed. The equation for the linear curve fit is y: 0.9755x + 0.3446. The R2 value is 0.99. .5 O) _L .b .4. N O h O N Measured velocity (m/s) . . .0 .0 . . . O l I I l 0.2 0.4 0.6 0.8 1 Pressure difference (in.-l-IZO) Figure 2.4 Calibration curve for the elbow meter. 22 1.2 Honeycomb was placed within the PVC pipe between the blower and the elbow meter. The 2.5cm thick honeycomb consists of a bundle of plastic tubes (similar to drinking straws) bound together. This suppresses any vorticity and inconsistencies in the airflow generated by the blower. As the air passes through the elbow meter it enters the sealed settling chamber where its velocity drops. The settling chamber also serves as the structural base for the MSU Fire Tunnel. Four Unistrut beams extend from the settling chamber’s plywood sides and run up the entire length of the MSU Fire Tunnel. These beams are bolted to a wood frame within the settling chamber. A steel platform beside the MSU Fire Tunnel allows for access to the combustion chamber (Fig. 2.2). 2.3 Converging Duct and Conditioner As the oxidizer flow emerges from the settling chamber it passes through a series of conditioning structures. Honeycomb and screens (wire cloth) are positioned inside the conditioner duct orthogonal to the mean flow direction. This is intended to suppress vorticity and turbulent eddies that may have formed within the settling chamber. The cross-sectional area of the conditioner duct is different than that of the settling chamber, requiring the use of a lip at the inlet of the conditioner. Two- inch PVC pipe was cut in half and secured to the edge of the conditioner inlet, thus providing a smooth transition for the flow from the settling chamber to the conditioner. 23 In order to accelerate, direct, and decrease turbulence in the oxidizer flow a contraction was built. The geometry of this contraction was chosen based on the study of [11], which outlines the design of two-dimensional wind-tunnel contractions. A two dimensional contraction was chosen because it provided the required plane symmetric flow within the combustion chamber and was relatively easy to construct. A one-dimensional contraction was not used in this study because of the desire for a plane symmetric velocity profile within the combustion chamber. A three-dimensional contraction was not chosen due to the complexity of its construction and an axisymmetric velocity profile was not required within the combustion chamber. The size of the two dimensional contraction was determined based on the dimensions of the combustion chamber. The outlet of the contraction has the same dimensions as the inlet of the combustion chamber. The contraction ratio (defined as the area of the contraction inlet divided by the area the contraction outlet) was chosen as five [12]. The mathematical method for determining the geometry of a contraction involves piecing together two cubic polynomials. The location at which these cubic polynomials meet (defined as X",) is determined by the chosen contraction ratio, the virtual origin of the boundary layer, the maximum wall pressure, the separation point of a turbulent boundary layer, as well as other factors. The quantity Xm is a normalized distance from the contraction inlet representing the location where the cubic polynomials meet. Xm is defined as zero at the contraction inlet and unity at the outlet. Once the point at which the two cubic 24 polynomials met was chosen, a computer program was employed to generate the geometry [13]. This computer program generates a series of points that were then plotted in Excel. The points depicted in Figure 2.5 describe the geometry of the converging duct. Converging Duct Geometry (X1, = 0.61) 2’ 0.3 E, 0.25 w A E E 0.2 3 *5 015 3 ‘3 0.1 5 '5', 0.05 I 0 0 0.2 0.4 0.6 Length of converging duct (m) Figure 2.5 Converging duct geometry, X,,n = 0.61. The contraction was constructed using an aluminum frame and flexible plastic sheets. The desired contraction geometry was cut out of a piece of 25 1.6mm thick aluminum sheet metal using a computer-controlled mill that interpolated between the points provided to create a continuous cut. A frame was created from these aluminum shapes, and the plastic sheets were glued into place with an epoxy called Plastic Welder. 2.4 Combustion Chamber The combustion chamber rests atop the contraction and is the location of the flame-spread tests. The aluminum frame of the combustion chamber is mounted to the Unistnrt supports. A non-combustible Marinite board mounted within the combustion chamber acts as a sample holder. Strip heaters bolted to a sliding aluminum frame (collectively known as the radiant panel) can be heated up outside of the combustion chamber and then rapidly slid into place, thus providing the sample with a relatively time invariant impulsive heat flux. The wall behind the sample holder is a 1.6mm thick aluminum plate. The sides of the combustion chamber are 6.4mm thick Pyrex glass. 2.4.1 Sample Holder The sample holder consists of an aluminum tray with a 90° bend at one end. A Marinite sheet with a thickness of 1.9cm was placed within the tray. The Marinite sheet had a 25.4cm x 22.80m x 1.30m square opening cut into its center using a router. The samples that were burned fit into this space. 26 1 .90m Figure 2.6 The sample holder of the MSU Fire Tunnel. The 10.50m section is the airfoil leading edge, out from aluminum using the EDM. The three sections are joined by a threaded rod. Shown also are the sample, 21 .6cm downstream of the leading edge. 27 The sample holder was fitted with an airfoil (bolted to the 90° bend) that was designed to create a stagnation point at its leading edge. The airfoil has a cross-section that is half of a 10:1 ellipse. In terms of the fluid dynamics, this is the most important feature of the combustion chamber because it fixes the location of the boundary layer leading edge at the stagnation point on the ellipse [14]. The equation used (an equation for an ellipse) was’5y=[b2(1-x2/a2)]"2, where a is the minor axis and b is the major axis of the ellipse. This equation was used in a program called Ez-Cam to generate a series of points for a computer controlled Electrical Discharge Machine (EDM). The EDM cut three 6” long, solid aluminum sections of the airfoil from three blocks of aluminum that were secured together using threaded rod. The airfoil was then bolted to the sample holder. Dry wall filler was used to fill in the small gap between the trailing edge of the airfoil and the leading edge of the Marinite board. It was sanded to create a smooth transition for the oxidizer flow. A sample (1.3 cm thick) is secured onto the void in the Marinite sheet and the entire sample holder was placed into the combustion chamber. The solid fuel samples used were 25.4cm x 22.8cm x 1.3cm clear Plexiglas G. This type of Plexiglas is cast, making it ideal for ignition and flame spread experiments because it does not melt or drip much. The fuel sample has mounting holes pre-drilled and is secured to the sample holder with nine steel mounting screws. These screws penetrate halfway into the fuel sample. 28 When the sample is placed into the sample holder, small gaps exist between the Marinite edge and the fuel sample edge. All the gaps except the one at the leading edge were filled with fiberglass insulation. The gap at the leading edge was filled with dry wall filler and sanded until smooth. This created a smooth transition for the flow from the sample holder to the fuel sample. 2.4.2 Ignition Assembly Nickel-Chromium resistive heating wire was used to ignite the sample. The wire (Omega Engineering, Part # Nl80—015-50ft) is 80% Nickel and 20% Chromium with a resistance of 8.435E-2 ohms/cm. The heating wire was wrapped around a hollow ceramic rod with a 3mm diameter. Two steel wires insulated by hollow ceramic rods (6mm outer diameter, not shown in Fig. 2.7) carry current through the combustion chamber walls to the Nickel-Chromium wire and also act as structural support for the ignition assembly. The ignition assembly was placed near the trailing edge of the sample, directly on the surface of the sample. A current of 5 amps raises this wire to a temperature of approximately 800° C. A variable voltage transformer was used to supply this current to the ignition wire. 29 Nickel-Chromium Wire Ceramic Rod Figure 2.7 Ignition wire assembly. Placement of the ignition assembly is given in Figure 2.8. 30 Figure 2.8 The Combustion Chamber of the MSU Ignition Assembly Radiant Panel Control System Fuel Sample Pyrex Sides of the Test Section Airfoil Fire Tunnel, showing the sample holder with the Inlet airfoil leading edge, fuel sample, ignition system, Flow radiant panel, and inlet flow. (U..) 3] 77.5cm 26.7cm +— 10.20m Figure 2.9 The general dimensions of the combustion chamber, showing airfoil with the sample, ignition assembly, and radiant panel. 32 2.4.3 Combustion Chamber Design Considerations It is desirable for the oxidizer flow in the region between the radiant panel and the sample holder to behave like that of free stream flow over an external flat plate. In order to obtain the dimensions of the combustion chamber that would allow this, boundary layer calculations were performed. Within these calculations it was assumed that the flow within the chamber was laminar, steady, incompressible, flow over a flat plate. By using the displacement thicknesses a determination of the desired size of the combustion chamber (H) was made. The distance from the leading edge of the airfoil to the leading edge of the sample (L,) and the volume flow rate necessary to produce a given free stream velocity were also calculated. Minimum chamber dimensions were sought so that the sum of the displacement thicknesses due to the sample holder and the radiant panel walls (5'u and 8'. respectively) was small compared to the distance between the radiant panel to the sample holder (H). C =IH-(6I+6Z)I (1) q H 33 ////////. //////////////////////////////////////V//////// Figure 2.10 The combustion chamber, airfoil, sample holder and relevant 34 parameters. Equation (1) can be used to determine the effective cross-sectional area for free stream flow. Cq was chosen to be 0.88. While the choice of Cq was arbitrary, such a value ensures that the boundary layer near the radiant panel will not significantly affect the boundary layer near the fuel sample. The leading edge of the fuel sample was used as the reference point along the length of the sample holder from which calculations were performed. By referring to figure (5) an approximation of the boundary layer thickness 655 ’ , (2) was established, where v is the kinematic viscosity, U0, is the free stream velocity, and L,( is the distance from the leading edge of the sample holder to the leading edge of the sample. The velocity U... was chosen to be 0.5 m/s (the lowest free stream velocity anticipated within the chamber). Thermal expansion of the gases in the combustion chamber occurs and a correction can be made to compensate for it. The kinematic viscosity (v) was used as the correction parameter. The fluid inside the combustion chamber was 35 approximated as air and the value for the kinematic viscosity of air at 1000 K was used in equation (2). The distance from the leading edge of the airfoil to the leading edge of the fuel sample (L) was then determined using (2). It is desirable for the boundary layer thickness to change as little as possible over the length of the fuel sample. The difference in thickness is minimized as the sample is moved downstream with respect to the leading edge of the airfoil. Size constraints prevent this distance from being too large so a compromise was reached. The percentage differences between the boundary layer thicknesses at the trailing and leading edge of the fuel sample for a U,o of 0.5 m/s and a fuel sample length of 23cm for various sample positions with respect to the airfoil are given in Table 2.1. Position of fuell B.L. B.L. Percentage leading edge thickness thickness difference (Lx) (cm) (leading) (trailing) (cm) (cm) 5 0.61 1.42 57.0 10 0.86 1.54 44.1 15 1.06 1.66 36.3 20 1.22 1.77 31.0 22 1.28 1.81 29.3 25 1.36 1.87 27.1 30 1.49 1.97 24.0 Table 2.1 Calculated percentage differences between the boundary layer thicknesses at the leading and trailing edges of a 230m fuel sample, U..,=0.5 m/s. From this table, a distance of 22cm gives a percentage difference between the boundary layers of less than 30%. This choice of less than a 30% 36 difference though arbitrary, represents a compromise between size constraints and the need to make the flow over the fuel sample as regular as possible. Once L,‘ (Fig. 2.10) is known it fixes the boundary layer thicknesses and the displacement thicknesses for a given U... for both the fuel sample wall and the radiant panel wall. It was assumed that the radiant panel wall has the same boundary layer behavior as the fuel sample wall (inexact, but adequate for this calculation). 5' = -—5 (3) Using (1),(2),and(3), the distance H between the fuel sample wall and the radiant panel wall can be calculated. H= 7cm was obtained. The general design of the sample holder within the combustion chamber (Fig. 2.10) is credited to [12]. Modifications were made from the original design and are discussed in greater detail in section 3.1. 2.4.4 Radiant Panel The radiant panel consists of seven OT-1801 Omegalux resistive strip heaters with an iron oxide finish. The strip heaters are mounted to a rigid aluminum frame and coated with a low gloss, black, high temperature paint. Each strip heater requires approximately 8 amps of current in order to achieve 37 full power. This paint enables the surface of the strip heaters to be approximated a blackbody radiator (ez0.95). This frame fits into tracks cut into the frame of the combustion chamber. This was done to allow the radiant panel to be rapidly slid in and out of the combustion chamber in order to avoid any power up phase of radiant sample heating. A control system was constructed for the MSU Fire Tunnel. The strip heaters are controlled using a CN9120A Temperature Controller. A closed loop control scheme is used (Figure 2.13) with “E” type thermocouples used for temperature monitoring. The thermocouple is mounted to a strip heater with Omegabond 300 to produce good thermal contact and high electrical resistance. Four circuits are controlled with the temperature control unit using four DC input/AC output solid state relays (SSRs). One of the SSRs is an Omega SSR240DC75 with a load current capability of 75 amps. The remaining three SSRs are Crozet Gordos GAS-6D25 with load capabilities of 25 amps each. All four SSRs can be controlled with a pulsed DC current between 3 and 32 volts. One circuit is powered with a 125 Volt 15 amp line. The 15 amp line was used because only one strip heater is connected to it and because the lab only has three 30 amp lines. This 15 amp circuit powers one strip heater and is controlled via the SSR. The remaining six strip heaters draw power from three 125-Volt 30- amp lines and are also controlled with SSRs. High temperature wire connects the SSRs to the strip heaters. High temperature ring terminals were used to secure the high temperature wire to the strip heaters. The 30-amp lines come in from a three way disconnect box with slow blow 25-amp fuses. Each 30-Amp 38 line then has a second “fast blow" fuse in the control box. The 15-Amp circuit has a 10-amp fuse. The Temperature Controller provides a DC pulse to operate the SSRs. Each SSR has a control side impedance of 10009 with a minimal voltage requirement of 3 Volts DC. The Temperature Controller can produce 5 Volts DC (pulsed) at 25mA. By connecting the SSRs in parallel (Fig. 2.11), the resistance “seen” by the temperature controller becomes smaller. The combined impedance is 2500 and the minimum current necessary to switch the SSRs is 20mA, 5mA lower than that provided by the Temperature Controller (5mA margin of safety). SSR #l SSR #2 SSR #3 SSR #4 Temp. Controller Figure 2.11 Control side wiring of the SSRs, which are connected in parallel in order to reduce the overall resistance required by the temperature controller to control the system. 39 During operation the SSRs will heat up. For this reason they are mounted to finned heat sinks. The SSRs and the rest of the control system is mounted within a metal enclosure that has a fan (capable of delivering 30 CFM) to provide cooling. Temperature Over-temperature controller sensor Cut-off switch for three 30- amp lines. .....-..o-... ...,.|-UD"”fl.t- ”‘ Controller/7th j? strip heater - T} power switch 40 The design of this control system has an inherent danger. Seven strip heaters are controlled yet only one is monitored. It is possible (however unlikely) that the strip heater that is being monitored could burn out. This would fool the control system into treating the strip heaters as if they were cooling, thus causing a possible overheat situation. Therefore as a safety precaution, a separate over-temperature sensor has been included in the control system (not included in Figure 2.13). A CN355-KF High Limit Controller monitors the temperature of the radiant panel via a thermocouple (again bonded to a strip heater with Omegabond 300). The CN355-KF has a trip light on it. If the measured temperature exceeds the set temperature the green light will go off. It can be connected to an audible alarm designed to signal the user in the event of an over-temperature potential emergency. Once the control system and radiant panel were constructed and in place an auto-tune feature was used to optimize the PID control parameters of the Temperature Controller. The controller goes through a succession of heat up and cool-down cycles to determine the system characteristics. This algorithm is built into the controller and needs to be performed only once. The parameters obtained are stored in the controller even when the power is turned off. 41 CO Th ac fro an 2.5 Exhaust Hood The exhaust hood is positioned downstream of the outlet of the combustion chamber and hangs from four chains secured to overhead supports. The steel exhaust hood is 91.50m x 91.50m at its inlet and 61cm deep. A 12.7cm diameter round pipe extends from the exhaust hood to an overhead duct that leads to a separate chemical hood. When the fan on the roof of the building is turned on the glass slide for the chemical hood must be closed in order to achieve the desired flow rate in the exhaust hood. The round duct extending from the exhaust hood is fitted with a butterfly valve in order to shut off the flow and increase the flow rate to the chemical hood. 42 Control System IAAAAAAAAA 1AAAAAAAAA SSR Strip + heaters 30A Fuse 125V IAAAAAAAAA ' E — IAAAAAAAAA |—- SSR Strip + heaters 30A Fuse 125V AAAAAAAAA ' r | . \. 1AAAAAAAAA SSR. SHITPI + , . heaters 30A Fuse 125V H + L wwww 15 A? SSR? Stn'p heater 125V Fuse Temperature Controller Figure 2.13 Basic outline of the radiant panel control system. 43 CHAPTER 3 TEST METHODOLOGIES AND DATA ACQUISITION TECHNIQUES Tests on various components of the MSU Fire Tunnel were performed. Because this is a brand new experimental apparatus, considerable time was also spent refining the flame spread test methodologies and the geometrical configuration of the test section. Ultimately, full-scale flame-spread tests were performed and the flame-spread rates were recorded. 3.1 Various Experimental Configurations Tested In the process of refining the experiment a variety of experimental configurations were tested and examined. Changes (from what was initially built) were made to the elbow meter, the inlet to the settling chamber, and the sample holder. A fin attached to the trailing edge of the sample holder was tested and eventually not used for flame spread tests (for practical reasons). Initially, the PVC pipe leading from the blower to the elbow meter contained three pieces of honeycomb and four wire cloth sections placed within it to suppress turbulence. It was found that the blower could not produce the desired free stream velocity within the combustion chamber with these pieces in place, so all but one piece of honeycomb was removed. Another part that was tested was a diffusion pipe within the settling chamber. As air entered the settling chamber it passed through a pipe with a 44 series of holes in it designed to diffuse the initial momentum of the incoming air over a larger area. Again, it was found that this impeded the blower to the degree to which it could not produce the required free stream velocity within the combustion chamber. Several refinements were made to the sample holder and the type of sample used. The sample that was used for the first flame-spread test was a piece of molded therrno-set acrylic 6.4mm thick. When this sample was exposed to heat it did burn; it also melted and dripped into the combustion chamber. It was found that the use of cast acrylic is the standard for flame spread tests and cast acrylic was used in all subsequent flame spread experiments. The first few flame-spread tests were done with 6.4mm thick cast acrylic mounted in the Marinite board. The first Marinite sample holder was had the same dimensions as the one described in Chapter 2 but differed in the space that the fuel sample was mounted. A 1.3cm gap was left between the back of the sample (the side not exposed to the radiant flux) and the aluminum sample holder tray. As flame spread occurred over the surface, the fire would creep around to the back of the sample and not only ruin the structural integrity of the sample but also create irregular surface flame fronts. A new Marinite sample holder board was made for a different sample thickness (1 .Scm). Also the number of screws that mounted the sample was increased to nine. This solved the problems of sample sagging and flame- spread to the back of the sample. The new Marinite board was made such that a 1 .3cm deep recessed region (not all of the way through the 1.7cm thick sheet) 45 just big enough for the fuel sample was cut in its center. The 1.30m thick sample was placed into this void. Flow profile and flame spread tests were also performed with a fin mounted to the sample holder. While Figure 2.10 depicts the primary test configuration used, an alternate test configuration was also utilized. The design for this configuration was the original design discussed in [12]. A fin was placed on the back end of the sample holder (Fig. 3.1). Because the sample holder is offset with respect to the centerline of the combustion chamber the pressure drop along a streamline on the front of the sample holder is different than on the back of the sample holder. It is thought that this moves the hydrodynamic stagnation point from the tip of the airfoil to some location behind it (near the back wall). This may be one of the causes of non-uniforrnity in the free stream velocity profile over the sample holder. Use of the fin at the trailing edge of the sample holder has proved to be problematic during flame spread tests. A first version of the fin was made from 0.8mm thick aluminum sheet and was secured on the trailing edge. Fire entered the gap between the fin and the sample holder burned a hole into it. A second fin was made from 1.6mm thick aluminum but similar problems occurred. While it is not clear that the presence of the fin significantly affects flame-spread, further study of this issue is recommended. 46 77f r////¢/J ........ v ///// ............ ///// ......... ////////////////// Figure 3.1 The sample holder including the fin on the trailing edge. 47 3.2 Component Testing The first of the components tested were the contraction and conditioner. Flow profile testing of the contraction was done to verify that the velocity profile was uniform at the outlet. It is desirable for the leading edge of the airfoil to be situated in the free stream flow. Therefore, the thickness of the boundary layer at the outlet of the converging duct determines the minimum distance that the sample holder can be placed with respect to the combustion chamber wall. For this reason, tests were performed to obtain an estimate of the boundary layer thickness at the contraction outlet. The velocity profile was measured across the short dimension of the converging duct outlet using a pitot tube probe and manometer. The outlet of the converging duct was fixed to a plenum held at a pressure lower than atmospheric. This test configuration included the flow conditioner, which was attached to the converging duct at one end and left open to ambient air at the other end (Fig. 3.2). Another component tested was the ignition wire system. Nickel-Chromium wire was wrapped around short lengths of ceramic rod and placed near small PMMA samples. In order to make ignition at the trailing edge of the fuel sample uniform, several parameters were adjusted until it was found that ignition was insensitive to their effect. The voltage per unit length of ignition wire, the spacing between 48 successive wraps and the proximity to the sample were varied until it was found that there was spanwise uniformity in ignition. Measurements taken here \ N Plenum pressure < atmospheric Figure 3.2 The configuration for the contraction/conditioner component tests. The next component tested was the radiant panel. Spatial variations in the surface temperature of the strip heaters will occur due to manufacturing inconsistencies, surface emmisivity variations, and conduction losses. This, in 49 turn, will cause variations in the surface heat flux. Thermal imaging of the radiant panel was performed in order to characterize this variation. Additionally, any discrepancies in temperature between that measured by the control system and that measured by the thermal imaging camera were sought. An lnframetrics 600L infrared camera capable of measuring infrared radiation wavelengths between 8 and 12 microns was used to image the radiant panel surface. In order to minimize natural convection the radiant panel was oriented such that it was horizontally with respect to the ground. Additionally, 8cm long baffles were secured along the edges to inhibit airflow around the edges of the radiant panel. Images were taken at two different temperatures and recorded on a VCR. Image processing software (lmagePro 3.0) was used to pseudo-color the images. The effective heating area was determined using the image processing software. The software was spatially calibrated using structures recognizable in the thermal image (i.e. lines that exist on the image were recognizable on the thermal image and represent the junction between successive strip heaters). The distance between the outer edges of four of the strip heaters is known to be 15.3cm and this was used as a calibration. Next, the length and width of the region within the thermal image that fell within the measured range of the camera was measured. 50 ‘ 2 ' 3 Position . 4 Number Figure 3.3 Locations of hot wire probe traverses. 51 3.3 Combustion Chamber Flow Measurements Tests of the velocity profile in the combustion chamber were performed under the standard operating configuration (MSU Fire tunnel Fig. 2.1). Point measurements of the velocity within the combustion chamber were taken with a TSI lFA 300 Constant Temperature Anemometer. The data-sampling rate was 2000 Hz and 16 Kpts were collected for each probe location. Thus each data point represents the average of 16,384 velocity measurements. Thermal Pro software was used to enable the calculation of the turbulence intensity to be made from the velocity data. The hot wire probe was traversed across the width of the test section between the sample holder and the radiant panel. The probe was inserted into holes cut in a panel just below the strip heaters. Five such holes were drilled, thus enabling data to be gathered at five different spanwise locations in the oxidizer flow. 3.4 Flame-Spread Tests Flame-spread tests were conducted within the combustion chamber. The purpose of these tests was to determine the relationship between the set radiant panel temperature and the flame-spread rate. Secondary goals were to determine the repeatability of the tests and to gain an understanding of the general nature of the test apparatus. The test apparatus was set up in the configuration described in Chapter 1. No fin was used in these flame-spread tests. A PMMA fuel sample was mounted 52 in the sample holder and the ignition wire was placed directly on it. The blower was turned on thus setting the free stream velocity to 2.2 m/s 10.1. The radiant panel was heated to the desired temperature and slid into the combustion chamber. A sample heat-up time of 2:00 minutes (chosen arbitrarily) was used before the ignition wire was powered up. The entire event was filmed with a Sony CCD-TRV75 NTSC video camera from one side of the combustion chamber through the glass side. As the ignition Ignition Wire \ l J b a C Flame Front b a c Figure 3.4 Solid fuel sample. Flame-spread measurements were made along a-a not b-b or c-c. 53 wire was powered up and subsequent combustion occurred, the aperture of the camera was closed in order to compensate for the brightness of the flame. The experiment continued until the entire sample burned or until 14 minutes elapsed. The flames were then extinguished using an injection of bottled Nitrogen gas. The videotape was analyzed using the image analysis software lmagePro 3.0 to determine the rate of flame-spread over the sample surface. The software was configured to grab an image from the videotape at a constant time interval. The distance from some reference point at the top of the sample to the flame front was measured for each picture acquired. The software enables the user to define a reference length with regard to the image. A “calibration” image was obtained and the length of the sample was chosen as the reference dimension. Using this method, each subsequent measurement on the sample can be made with respect to the defined dimension, allowing absolute measurements to be obtained fairly easily. An origin at the top center of the sample was chosen for each test run and was referenced in terms of its pixel location. This pixel was then used as the starting point for the measurements of the flame front location. In order for this to be a valid technique it must be assumed that the camera does not move significantly during one test. Measurements of the flame front location were made “visually", i.e. by just drawing a line from the origin pixel to the leading edge of the flame. This appears to be an inaccurate method but because the aperture of the camera was 54 adjusted to compensate for the brightness of the flame the difference between the flame and the fuel sample was stark and could easily be visually resolved to two or three pixels. Measurements of the location of the flame front were always made along the centerline of the sample (Fig. 3.4). The flame fronts on the fuel surface were non-uniform in the spanwise direction but the center always spread the furthest. 55 CHAPTER 4 RESULTS AND DISCUSSION 4.1 Introduction Presented are the results of the various component tests and full scale flame-spread tests. Measurements of the velocity profile at the outlet of the converging section (in a testing configuration) for various free stream velocities were performed. Tests of the flow profile in the combustion chamber were conducted with no combustion and no radiant heating. Thermal images of the radiant panel were taken. Finally, flame-spread tests were conducted at several radiant panel set temperatures. 4.2 Converging Duct Testing Flow profile testing of the converging section was done to verify that the velocity profile was uniform at the outlet of the contraction. These tests were also perf0rmed in order to obtain an estimate for the boundary layer thickness. The error bars represent errors accumulated from uncertainties in the reading of the pressure transducer as well as errors associated with its calibration. The Figures 4.1 through 4.4 represent measurements taken in the manner described in section 3.3. They indicate a uniform velocity profile in the bulk of the outlet of the conditioner as well as relatively thin boundary layers. 56 insocnyr~unksntcankacrczztnaaer {iiiiiiiii i i I i i I i I i if “” I :0 PM 0 n i-o-r no ”-1 Velocity (mIs) 7‘ (.71 ql-l GD T T I 1 I -1.5 05 25 4.5 65 8.5 105 firms pos'ilm (on) Figure 4.1 The velocity profile across the outlet of the converging duct at a bulk velocity of approximately 3.7 m/s as measured by a pitot tube. The error bars represent errors accumulated from uncertainties in the reading from the pressure transducer as well as errors associated with its calibration. 57 VelocitdeileRQIletdCa'lrafim C0 ‘~I Iiiiiiiii i i i i l i Hit, CD U'I Velocity (mls) l. I‘D unl CD I 1 1 I -1.5 05 25 4.5 as as 105 Fiche position (err) Figure 4.2 The velocity profile across the outlet of the converging duct at a bulk velocity of approximately 6.7 We as measured by a pitot tube. 58 do VelocityProfileNOuletofConiraction 7 mm l i i I i mm; ”i i 8? g 3 ii i i i 0:. 2‘. .1. .5 .‘5 would"? 10.5 Figure 4.3 The velocity profile across the outlet of the converging duct at a bulk velocity of approximately 6.7 m/s as measured by a pitot tube. 59 VelocityProfileAtQJtIetofContraclion 2.25 2.00 175 {Hi i i i i i i H; “1.50 V1.25 1.00 0.75 0.50: I 0.25 0.00 l i l i i 0 2 4 6 8 10 Probe position (cm) HH mls f-O-I Velocity HH Figure 4.4 The velocity profile across the outlet of the converging duct at a bulk velocity of approximately 1.8 m/s as measured by a pitot tube. 4.3 Radiant Panel Thermal Images Thermal images were taken at two radiant panel set temperatures (232° C and 315° C) using an infrared camera. Figure 4.5 represents the gray-scale image acquired at the set temperature of 232° C along with the dimensions. It is difficult to determine the relative intensity of the thermal radiation using a gray- scale image so a pseudo-colored image was also produced. Figure 4.6 represents the pseudo-colored version of figure 4.5. Each color division represents 3.4° C. The highest temperature level displayed is 298° C and the lowest is 264° C. It is important to note that the camera has a finite range and everything that is blue is 264° C or lower but obviously not necessarily equal to 264° C. The camera (operating in this mode) has an error of approximately 13° C. Figure 4.7 represents the gray-scale image of the radiant panel set to 315° C. The dimensions are also given. Figure 4.8 represents the pseudo-colored version of figure 4.7. Each color division represents 3.4° C. The highest temperature level displayed is 386° C and the lowest is 342° C. 61 u.) 4—————— 30.50111 #1 + 3540: these NUDE sens: +24.3°c + 314°C Figure 4.5 Thermal image of the radiant panel set to 232° C (gray-scale). 62 53 524 , '96 ifiiiFRmE'rmcs II “If. llll Illl. II “ ' IMHGE MODE EC I-C..= +24 8°C + 3114“: Figure 4.6 Pseudo-colored thermal image of the radiant panel set to 232° C. 63 Figure 4.7 Thermal image of the radiant panel set to 315° C (gray-scale). ‘61-. Figure 4.8 Pseudo-colored thermal image of the radiant panel set to 315° C. 65 In order to determine the heat flux variation across the width of the radiant panel a plot was generated for a simplified geometry (Fig. 4.9). D" .- Y V , -;-...wegh.~§‘~;._-u ~ -;'}-_-A‘ . mi. 2.3“. 1.9%, a-ifi.fi&:§4fi'mfig.fi rent-H-502; . g. 4 Ki, . . ‘ . . - . . . . . . . \ . « . . ‘ , ,, ‘ . . , - 5 ,. ..,. _._ .. ‘ . . a . ~. ..' . . . . . , l . . : '“E , a: v.1" -fi'me1£”.\. . 13" r U, ‘—;:~ 3 , ‘ 1' e.‘ "9"\‘A“"‘ ‘V‘A' . (3".43‘ ,-- n. 'g ,.:.A- "rV‘w-w ”v -‘ ~ 1...: L7. Vv-‘~.l-" ‘Avu l‘ “’ >| Figure 4.9 Parallel, isothermal, infinitely long plates for radiant heat flux calculation. The geometrical configuration depicted in Figure 4.9 is solved in [15]. The system considered is a set of two infinitely long parallel plates held at a constant temperature of the same width and placed in a large, cold environment. The emssivities of the two plates are considered to be the same (e=0.95) and in this case h = w = 1. The solution was obtained using a 1m successive approximation technique. As given in [15] the nondimensionalized heat flux loss to the surroundings, ‘P, can be represented by the following equation = _E: W": ___§____ 2 _ 2 4 W) s 2iflm—:)2+1+52I0[8(1£)WI () Where W = w/h = 1 and the nondimensionalized width parameter 5 = x/h. As 5 —> 1 the term o[s"’(1-s)W2] —> 0 so in this case of a “high” a it maybe discarded. Figure 4.10 represents the nondimensionalized heat flux as a function of the location along the plates. From figure 4.10 it can be seen that the heat flux losses are greatest near the edges and minimal near the center. This calculation is obviously not a perfect representation of the true radiation problem but it represents adequately the conditions of the radiant panel. This calculation could be modified to more accurately model the problem. 67 Nondimensional Heat Flux vs. Plate Loaction .9 co .° 00 A v A V .0 O) o (0 Nondimensional heat flux loss, We: q/eo T"4 p c or \I I I I I I I I I I 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Location along plates, §/x P .h 0 Figure 4.10 Non-dimensional heat flux vs. plate location. Note that the losses to the surroundings are the highest near the edges and that the radiant flux is the highest in the middle. 68 4.4 Combustion Chamber Flow Profile Tests Tests of the velocity profile in the combustion chamber were performed under the standard operating configuration (MSU Fire tunnel Fig. 2.1). Point measurements of the velocity within the combustion chamber were taken with a TSI lFA 300 Constant Temperature (hot wire) Anemometer. Also presented are data representing the turbulence intensity across the combustion chamber. The hot wire probe was built and calibrated using another test facility (The Ford MSU Test Facility, FMTF). The calibration was performed in a range from 2.5m/s to around 10m/s. Error can be generally estimated from these sorts of calibrations as approximately 2% when measurements are made within the calibrated range. In this case an extrapolation was made so error bars with 14% were used. A large number of data points (16,384) were used for each velocity data point represented on the figures. The reliability of 7c (the sample mean) as a )°'5 where o is the standard deviation and measure of u (the true mean) is 1- o/(n n is the sample size [16]. This error vanishes with the sample size used. The hot wire probe was traversed across the width of the test section between the sample holder and the radiant panel. The probe was inserted into a hole cut in a panel just below the strip heaters. Five such holes were drilled, thus enabling data to be obtained at five different spanwise locations in the oxidizer flow. The first sets of tests were designed to determine the uniformity of the flow within the test section. These tests were conducted without the fin attached to 69 the trailing edge of the sample holder. Figure 4.11 represents the flow profile through the test section for this first set of tests. The zero position is at the fuel sample wall. The measurements of the flow depicted in Figure 4.11 indicates non- uniformity in the velocity profile when no fin is attached to the trailing edge of the sample holder. This non-uniformity appears to be consistent through the width of the combustion chamber. Next, measurements of the velocity profile in the combustion chamber were performed with the fin in place at the trailing edge of the sample holder. The fin impedes the flow and thus drops the velocity compared to the case with no fin (for the same blower throttle plate position). Figure 4.12 indicates the velocity through the test section, at position #3, for the case with the fin. Around the probe locations of about 7cm the velocities are high. Since the pressure in the combustion chamber is higher with the fin in place compared to the case with no fin (compared to atmospheric), a greater pressure differential exists. This causes a high flow velocity near the hole for the hot wire probe. Because there was no hole in the wall of the combustion chamber during normal operation, the near wall velocity data was ignored when comparing profiles with and without the fin. 70 mumenmmmmommmmzemo Fin) 35 stillntlllfgg“ 3.4° . 3‘WAL u t' O O I A25 A+ a O E , . E I o 2 omits % + Imam“ > AMOS 21.5 +me#2 3 I 8 2 1 05 o I I I I I 1 RT 0 1 2 3 4 5 6 7 8 Prunloedlaum') Figure 4.11 Velocity profile through the combustion chamber at four position locations, with no fin at trailing edge of the sample holder. The 1:4% error bars are not displayed. 71 Use of the sample standard deviation of the velocity data is an effective way of determining the degree of variation across the bulk flow. The sample variances and standard deviations of the velocity through the test section for the cases with and without the fin are compared (Figures 4.11 and 4.12 for position #3). The first three and last three data points were discarded for this calculation (for both cases, with and without the fin). | |Variance Standard deviation | [Without fin [1.405-02 1205-01 | |wrth fin |1.70E-03 4.105—02 | Table 4.1 The variance and standard deviation of velocity measurements through the test section with and without a fin. There is an order of magnitude difference between the two standard deviations. The bulk velocity tested with the fin in place was lower, but this cannot account for such a large change in the standard deviation of the velocity profile. Clearly, the fin is effective for reducing variations in the bulk velocity profile. While knowledge of the bulk flow behavior is useful, the subject of most interest is in the near-fuel region. For this reason, tests of the velocity profile near the fuel sample were conducted. 72 veodiynuiienmgiuemmommmm 4.5 I-"O—l 3.5 OD 8 N I Measured Velocity (mls) a .ngnununinnni I —L D 01 O I I I I 0 1 2 3 4 5 6 7 Rwelocaiaflar) Figure 4.12 Velocity profile through the combustion chamber (position #3) with the fin. 73 For reasons stated earlier, the existing fin design is impractical for flame spread tests. Because of this, the near wall tests were conducted without the fin in place. A micrometer was used to achieve the fine degree of spatial resolution necessary to make measurements in the near wall boundary layer. It must be noted that for Figures 4.13, 4.14, and 4.15 the location of the sample wall was not precisely known at the time the measurements were taken. The wall location was determined by comparing the data to the Blasius profile. Figures 4.13 and 4.14 indicate that the velocity profile is quite similar to the Blasius profile. The deviation that does exist may be due to a pressure gradient in the free-stream. The Reynolds number based on the momentum thickness (Ree) was calculated at the leading edge of the fuel sample. The momentum thickness is defined as m a = 0.664%; (5) for the Blasius solution. For a free stream velocity of 2.25rn/s Ree is 180. 74 mmrveocirypmiire Wflflbfin) 2.5 {III 2 €15 % i 0 m g 1 l —aasius 0.5 0 T i l 0 5 10 15 Z) Riemann!" Figure 4.13 Near sample wall velocity profile, position #3, with no fin. 75 MMIVdodtdeilaPosifimflMfin) 25 WI i i g 2 5' 31.5 E. g o Nbarerrerts 'g 1 —Ba9ts o. 05 0 l r 0 10 15 20 Wodtymifs) Figure 4.14 Near sample wall measured velocity profile and Blasius profile, position #4, with no fin. 76 MMITuulemelrtensitthimffl, nofin) 4.5 4 35 I S I} 5 3 f 525 I .3 I I § 2 IE; WWI—I I .9 I 3 e 15 3 .— 1 0.5 0 I I F l l I T I 0 2 4 6 8 10 12 14 16 18 Pithelocdiaimn') Figure 4.15 Turbulence intensity in the near sample wall boundary layer for position #3 (velocity data from figure 4.13 were used). 77 4.5 Flame-Spread Tests Flame-spread tests were conducted in the manner presented in section 3.5. Presented are the results of these tests. The set radiant panel temperature was varied in increments of 50° F for the first six tests. A jump up of 100° F was taken for the seventh experiment. The remaining experiments (#8-#11) involved repeating the test for a radiant panel set temperature of 260° C in order to determine the consistency of the experiment. A least squares linear regression analysis was conducted for each experiment. The lines generated are displayed on the corresponding figures. In order to make effective comparisons the least squares line was forced to go through the origin. The equation for the linear fit and the R2 value are displayed on each figure. During the beginning of each test, the ignition wire was warming up and began to light the sample. Because the spatial proximity of the ignition wire to the sample surface will strongly affect the heat absorbed by the sample, variations will occur in the time to ignition. In other words, a change in ignition wire placement from one test to another as little as one millimeter will alter the time to ignition (elapsed time from when the ignition wire was turned on). Because the video data was fed into the computer at regular intervals and the first image was “snapped” when the ignition wire was tumed on, an error exists in the location of the vertical axis from test to test. While this error is present in all the data it is not large and doesn’t significantly affect the slope of the linear fit. 78 2201224PH Figure 4.16 The flame front within the combustion chamber (measurement line shown). 79 4.5.1 Error Analysis Error analysis was performed to determine the error in estimating the slope and intercept of the linear regression applied to the flame spread data. The standard error (8,) can be expressed as Sy = {2,-0.- ‘mx; *bI/(M-ZHUZ (6) where m and b are the slope and intercepts of the least-squares linear fit. M is the total number of measurements. The error in estimating m and b can be expressed as Sm = (N /A)"2Sy _ (7) S, = (2x? /A)“23y (8) Where A = (2219.)2 — m3 (9) It was found that the errors in estimating m and b were small. Table 4.2 indicates the errors in the estimation of the slopes and intercepts of the linear 80 least-squares fit to the data shown in Figures 4.17 through 4.24 (excluding 4.23, where no least-squares fit was applied). The slopes that are reported in Figures 4.17 through 4.24 are given only to the significant figure that is reliable. 10 1 11 -05 Table 4.2 Errors in estimating the slope and intercept of the least-squares linear fit for the flame-spread tests. The error in resolving the location of the flame front on the pictures acquired into the image processing software was around two pixels. This translates into an error of one or two percent for each data point shown in the flame spread tests. Error bars indicating this error would be about the same size as the data point itself so they were not included. 81 4.5.2 Flame-Spread Test Results Figure 4.17 (corresponding to a set radiant panel temperature of 177° C, the lowest panel temperature tested) indicates a fairly regular spread rate. It appears as though the spread rate is accelerating (with respect to the least squares fit) over time. Figure 4.18 indicates the same general trend. The spread rate accelerates with respect to the least squares fit and then begins to level off near the end of test. Figure 4.19 indicates the same general trend of acceleration of the flame front near the top of the sample. It can be seen cleariy that this acceleration wanes and the flame front begins to decelerate. This inflection point appears to occur at around the 14.5cm point along the sample length. Figure 4.20 displays this same general trend as the previous figures. The flame-spread rate is higher and the flame front spread further than the previous tests. The inflection point appears to be at the same point (around 14.5 or 15cm) as in the previous tests. This seems to indicate that this is a characteristic of the experimental configuration. I Figure 4.21 (corresponding to experiment #5) indicates a vigorous flame spread rate. It follows the general trend of early acceleration followed by a deceleration phase. This time the inflection point appears to have moved up to approximately 15 to 15.50m. 82 Ramfiutposifimvsfinefipmmpmel setterp=1770 .a co _s O) —L .s N 1h \ .5 o O) 00 \ 0 Flame front posltlon (cm) A \ O 4 .\ TING”) Figure 4.17 Flame front position vs. time for a set radiant panel temperature of 177° C, Experiment #1. 83 Rarefiutposiflmvsfinefipmmpa'elsetterwq .L .h .s N O a It? I56 89 Flame front position (cm) to «h C» 00 \ O \ Figure 4.18 Flame front position vs. time for a set radiant panel temperature of 204° 0, Experiment #2. 84 RarefiortposiflmvsfineExpanPa'elsetMC) .i. a: —L a) —l .5 .L N =.0021x /: R2: =090/ .s C O) 00 \ 4 Flame front posltlon (cm) N -k \ O The (see) Figure 4.19 Flame front position vs. time for a set radiant panel temperature of 232° C, Experiment # 3. 85 RarefiutposffimvsfineEmMWpa'elsettawq l3 e 5 \ 63 63 1% II p 8 .s .h 0 4 —L N —L C \ Flame front position (cm) C N «5 O) m .\ O {3‘ 8 8 6 6 8 a Figure 4.20 Flame front position vs. time for a set radiant panel temperature of 260° C, Experiment #4. 86 Ramfiutposifimisfineawmmsetmq 23 26 y=0CB7x / :4 Ri=096 / , 22:21 / C ./" .5“ . 7 E16 ./ 314 , 512 , / i: / 010 5 a / u. 6 / 4 //’ 0 2 e 0 r r . . . . T Figure 4.21 Flame front position vs. time for a set radiant panel temperature of 288° C, Experiment #5. 87 RamirorlposifimvsfineEipilemsettemhm6Q 28 26 Y=QI3L 24 / 22 #:09 / E2) /o ° 0 ~€18 e ./ %16 j / 3.14 / E12 . / E10 / E 8 ./ 6 / 4 / 2 / 0 . . . . o 8 MG”) Figure 4.22 Flame front position vs. time for a set radiant panel temperature of 316° C, Experiment #6. 88 1(D 150 Z!) 250 Figure 4.22 (corresponding to experiment #6) indicates a much faster flame spread rate. Because the spread rate was much faster than the previous tests the video capture interval was reduced to 20 seconds. The inflection point has again moved up to around 16cm. As the flame spread rate increases the linearity of the data seems to decrease. The acceleration in the beginning of the test (near the top of the sample) may be explained by recognizing the fact that the flame front is moving into a region of higher incident heat flux. Four effects possibly cause this higher heat flux near the middle of the test period. The first is the increase in the radiative view factor near the center of the sample, resulting in a higher heat flux. The second is an increased radiant panel temperature near the center compared to the edges (as indicated by the radiant panel thermal images) due to conduction losses near the top. The third explanation for this phenomenon could be that there is an increase in the total sample area that is burning. As more of the sample burns it produces radiation that is absorbed by the unburned fuel sample both directly and from reflection within the combustion chamber. The fourth factor could be a transient issue. The center of the sample has most likely not reached thermal equilibrium as the flame front begins at the top of the sample. As the flame front progresses the sample absorbs more and more heat causing a variation (specifically an acceleration) in the flame-spread rate. As the flame front moves through this region of high heat flux toward the leading edge of the sample it begins to slow down. This slowdown could be due 89 to several effects. The first reason for this could be that the heat flux impinging on the unburned surface decreases. This occurs as a result of a decrease in both the radiant panel temperature (due again to conduction losses at the edges) and in the view factor. Additionally, a decrease in temperature at the leading edge of the radiant panel is sure to occur due to forced convection losses, also contributing to a decrease in the heat flux. An additional explanation in the slowdown of the flame spread rate is the thinning of both the hydrodynamic and thermal boundary layers near the leading edge of the sample (in comparison to the trailing edge). This causes an increase in the oxidizer velocity in the near fuel region where combustion is occurring. The increased oxidizer velocity carries away both more gasified fuel and more heat from the sample’s surface, which produces a resultant decrease in the flame-spread rate. As the temperature of the radiant panel increases the effects that contribute to a variation in flame spread rate also increase. This causes a greater deviation from the linear fit resulting in a lowered R2 value. Figure 4.23 indicates a much greater increase in the flame-spread rate compared to the previous tests. This test exhibits (while impossible to glean from the figure) the “Oscillatory Flame Advancement (OFA)” behavior discussed in Chapter 5. Because the points are so deviant from any linear curve fit one was not applied and the test was not used in the later comparisons of flame spread rate. RarefiortposifimvsfineExpff/(Raipaelsetterrpfiflq 81892 63 I .5 O) O .5 A O .5 O O Flame front Position (cm) '16 ore-50100 Time (see) Figure 4.23 Flame front position vs. time for a set radiant panel temperature of 371° C, Experiment #7. 91 Figure 4.24 represent repeated versions of experiment #4 (radiant panel set temperature of 260° C). It indicates that the repeatability of the experiment is fairly good. The slope of the linear fit variation is small (although the number of tests is not high enough to calculate a reliable standard deviation). Figure 4.25 represents the linear regression applied to the averaged data at each time level from Figure 4.24. This is the center dotted line shown. The two lines around represent the upper and lower maximum errors associated with estimating the slope and intercept for these data. The error bars are quite close to the central line. This indicates that the variation of the data seen in Figure 4.24 is most likely not due to errors in estimating the slope and intercept of the linear curve fit. 92 Rarefratpos’llmvstirreEmMB—11,(Mmel setterp=2880 81852 .5 m d O) .5 .5 .5 o 0 5004 I 59” A 59” X 50m 0 Bpfl1 Flame front position (cm) 1'6 010me The(sec) Figure 4.24 Flame front position vs. time for a set radiant panel temperature of 260° C, Experiments #4, 8-11. 93 Umfigeesfmfittohmfirp. #5811)vrilhlllhia'mmEirors inSlopea'idlrleroen prererrory= 0(2841+02 / 81833 UnearegessimfittoAredEmM8-11y=0m(dwiedlire) / Lmererrory= 0029(— 0.2 Flame frgnt gosltlonJcm) N Figure 4.25 Slope of least squares fit for various radiant panel temperature settings, Experiments #1 -6, 8-11. 94 Figure 4.26 is a plot of the slope of the linear curve fits (roughly the average flame spread rate) versus the set radiant panel temperature. This indicates that there is marked increase in the flame-spread rate with the increase in radiant panel temperature. Experiments #1 and #2 may be called into question because they don’t necessarily follow the trend. The lack of variation from experiment to experiment is exhibited well by the data taken at the radiant panel set temperature of 260° C. 95 SopedLquaethvsSetFniatPaidTerpa-aue 0.14 0.12 p 8 Slope (cm/s) p 8 004 . 0C2 I Figure 4.26 Slope of least squares fit for various radiant panel temperature settings, Experiments #1 -6, 8-11. 96 CHAPTER 5 OSCILLATORY FLAME ADVANCEMENT During the course of testing an interesting phenomenon was discovered. This so-called Oscillatory Flame Advancement (OFA) occurs only under a unique set of operating conditions and was found to involve a large portion of the sample surface area. A conceptual description of this phenomenon and a preliminary f theory for its behavior are given. The area of interest is given in Figure 5.1. For OFA to occur, a portion of the fuel sample’s surface must begin to degrade into a highly viscous melt layer. In order to form this melt layer the radiant panel must be heated up to a set temperature above 320° C for a oxidizer flow rate of around 2.25m/s. The panel must be slid into the combustion chamber and a long heat-up period must be given to allow the sample’s surface to begin to “boil” or actively degrade. In this state the temperature at which a balance is struck between total heat losses and the heat gained due to the external radiation at the sample surface is higher than the boiling temperature of the fuel. Therefore the sample’s surface begins to boil. After a significant fraction of the surface begins to boil (approximately 2/3 to 3/4 of the surface area) the ignition wire is turned on. A flame front appears near the ignition source and very quickly sweeps over the entire degraded surface. As quickly as it appears it recedes back to the ignition source, completely extinguished over approximately 95% of the sample’s surface. This cycle then repeats itself but for this cycle the time interval between full 97 advancement and regression is smaller. With each cycle the process continues and the time interval between advancement and regression is considerably smaller. Additionally, the line to which the flame fully regresses begins to creep down the sample’s surface with each successive cycle until the advancement front and regression front meet. Like a rubber ball bouncing between two parallel converging walls that finally leave no room for movement, the flame front oscillations reach a steady state and ultimately a continuous flame engulfs the sample surface. What causes this unstable behavior? As with other complex transient phenomena it results from a delicate balance between several interrelated competing mechanisms. The rate of fuel evaporation from the liquid state to the gaseous state is related to the temperature of the liquid state. As the temperature of the liquid state increases so does the rate of fuel evaporation to the gaseous state such that rick,” = f (Tswm) In order for combustion to occur an adequate concentration of fuel (known as the lean limit) must be present in some region Ax (the quenching distance) away from the fuel surface of thickness Ax: (the flame thickness). In order to produce the lean limit at a distance Ax, therefore, a minimum surface temperature must be achieved. Since the flow field tends to carry away both heat and evaporated fuel, this surface temperature must be higher in the presence of a forced oxidizer flow. Both the external heat flux and the heat flux produced by the flame flash increase this surface temperature. In this way the 98 heat flux is directly related, through the fuel surface temperature, to the presence of a lean limit. When the concentration of evaporated fuel reaches the lean limit at a distance Ax and a pilot flame is present at the top edge of the decomposing sample combustion occurs a flame sheet flashes across the surface. The pilot is situated near the trailing edge of the fuel sample where the boundary layer of evaporated fuel is the thickest. The concentration at Ax is above the lean limit and the flame begins to travel down the sample, still in a region where the concentration is above the lean limit. As it travels downward it uses up the fuel in the volume of thickness Mb The temperature at the surface of the fuel is not high enough to produce a fuel evaporation rate that will replace this fuel quick enough. So after the fuel is all burned up within the critical volume the flame extinguishes. The fuel is all burned near the leading edge of the flame front before trailing edge because the thermal, hydrodynamic, and evaporated fuel boundary layers are at their thinnest in this region. As the flame front burns over the fuel surface (and even after it extinguishes) a heat flux from the flame flash is transferred to the fuel. This heat flux incrementally increases the temperature (T -» T + 67') of the liquid layer and therefore changes the evaporation rate, m": n't(T + 67') . The time required for the concentration of the evaporated fuel to reach the lean limit within the region Ax + Ax. now becomes smaller. In other words as the evaporation rate increases the time delay between successive flame front advancements decreases (Fig. 5.2). 99 This effect creates an acceleration in the oscillatory cycle until a steady state is reached. The evaporation rate has increased to such a degree by the time steady state is reached that it is now delivering an ample amount of fuel to the burning region even to the point of causing fuel rich combustion (noted by its yellowish-orange appearance and soot laden products of combustion). The change in surface temperature necessary to produce the large increase in fuel evaporation seen, however, is most likely not caused by the external heat flux. Since the time required for this whole process (from first flash to continuous flame) to occur is between 10 and 15 seconds the heat flux from the flame front seems to be the more likely candidate. Because the heat flux from the flame is most likely the driving mechanism behind the acceleration in the oscillations perhaps these two parameters could be linked. This may enable measurements of the oscillation frequency to be used as a rough predictor of the heat flux from the flame to the fuel surface (a rather elusive quantity). The creation of a rigorous theoretical model of OFA could prove difficult due to the complexity of the system. An easier (and perhaps more fmitful) path would make used of some simplifications in the governing physical system. One such simplification would eliminate the reaction equations by simply replacing the combustion with a transient volumetric heat generation term, inserted into the energy equation, but it is not yet clear how accurate such a model might be. The development of such a model is beyond the scope of this thesis. 100 1. Flame front heatfluxto fu I Thin melt layer —> S‘ ++++++ mfiwlmp /' Partially premixed fuel/air mixture $155533?“ Oxidizer flow +X Triple flame front of flash External heat flux 4 ext {— Figure 5.1 The region under consideration for the OFA phenomenon. 101 Continuous Flash interval z constant flame / r W l - n q to _ surface <——> At; time t;-)Oast—>oo Figure 5.2 The heat flux to the surface as a function of time for the OFA theory. 102 CHAPTER 6 CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK Preliminary testing of the MSU Fire Tunnel and its components indicate that nearly all of the original design goals have been met. The flow field within the chamber is predictable, the radiant panel is effective, and the flame-spread tests were successful and appear highly repeatable. It has been shown that the flow field within the combustion chamber is largely laminar and predictable. The hydrodynamic boundary layer near the fuel sample surface has been shown to be thin and nearly unaffected by the rest of the combustion chamber surfaces. The radiant panel has proven to be a durable radiant energy source that is easy to regulate, control and use. While there seems to be some non- uniforrnities across its surface they are not large and with a little effort could be easily corrected. Flame-spread tests indicate that the MSU Fire Tunnel is capable of producing the predictable, repeatable environmental conditions necessary for precise solid fuel ignition and flame spread studies. While the tests need to be refined the initial results look promising and indicate that this device could contribute a great deal to the knowledge of solid fuel combustion. The MSU Fire Tunnel is a versatile device. It can perform combustion as well as basic heat transfer experiments. For example, temperature profiles along the surface of the sample have been generated for use in other studies using a 103 data acquisition unit and a series of thermocouples. The oxygen concentration of the oxidizer flow, the flow velocity, the level of radiant flux, the type of sample tested, the location of the ignition source, as well as other parameters can all be vaned. While the MSU Fire Tunnel has exhibited promising preliminary results, several improvements could be made to increase its effectiveness. The major areas were improvements need to be made are in the flow system, the radiant panel, and the sample holder. While not really in the scope of this project the method of acquiring data could also be improved quite easily. A few improvements in the flow system would improve the precision and operating range of the device. The blower is a bit undersized and a more powerful blower would allow for a greater range of flow velocities to be studied. It may even allow for the study of a forced turbulent flow field, an interesting problem. Also, examining the flow profiles with radiant heating (and possibly even a burning sample) would prove very useful. Particle seeding techniques would most likely be necessary to obtain this data. These techniques, while expensive and somewhat complicated, enable measurements to be obtained under high temperature conditions. The uniformity of the radiant panel could be improved by using an additional set of radiant strip heaters. In order to compensate for the lack of a uniform heat flux (lower levels near the edges) due to the view factor and conduction losses the temperature of the edge of the panel could be raised using another set of strip heaters. The heaters could be placed on the back of the 104 existing heaters and bound with a highly conductive paste thus creating a temperature bulge near the edges of the radiant panel. This bulge could be fine tuned using the thermal imaging system and theoretical calculations in order to compensate for the view factor difference. This may allow a uniform flux to be produced from the perspective of the fuel sample. Making the heat flux more uniform will most likely correct the u-shaped non-uniformity present in nearly all of the tested samples. Also, isolation of the various effects that may cause the velocity variation in the flame spread (as indicated by the inflection point on the flame-spread rate curve) would allow for a greater knowledge of the relative intensities of these effects. The radiant panel could be moved with respect to the fuel sample to determine the relative effect of the radiant panel view factor. The fuel sample could be moved with respect to the elliptical shape to determine the relative effect of the thinning boundary layers. Further investigation into the near-fuel boundary layer is necessary. Investigation into the validity of the flat plate assumption would be useful. In order to determine if there is no pressure gradient in the free stream (implicit in the flat plate assumption) velocity measurements could be made at various streamwise locations. Additionally, the stability of the boundary layer could be investigated by performing more hot wire tests. Improvements in the sample holder may correct the smaller wave-like non-unifonnities along the flame front. As the flame-spread tests progressed the dry wall filler used to smooth the transition between the Marinite sample holder 105 and the leading edge of the fuel sample began to break dOWn. This may cause some unpredictable fluid flow over the sample surface. Better, high temperature filler could be used to correct this. Better measurement techniques would certainly improve the resolution of the data. Instead of simply measuring the furthest point of flame spread it may be possible to characterize an average flame front location as well as characterize the degree of non-uniformity along the flame front using image processing tools and some simple mathematical models. Measurements of the fuel surface (as well as within the surface) could easily be performed using thermocouples and the existing data acquisition unit. The addition of the fin on the trailing edge of the sample holder would improve the uniformity of the flow profile within the combustion chamber. Several designs were tested and they failed for a number of reasons. Because the flames impinge on the surface of the fin it gets very hot. This causes the fin to expand and warp, causing sealing problems. If the seal between the trailing edge of the sample holder and the leading edge of the fin is not perfect, the fire tends to go through it causing the fin to warp even more. One possible solution would be to avoid using the fin altogether and introduce a separate blower. This blower could be arranged so that it sucked air out of the combustion chamber behind the end of the sample holder. Several hoses could be used so as to avoid point source suction that may compromise spanwise uniformity. Further investigation into the Oscillatory Flame Advancement phenomenon could lend additional insight into nature of flame spread over solid 106 fuels. The frequency of the oscillations in the experimental device could be obtained by simply monitoring the audio signal generated. As the flame advances over the fuel surface it creates a noticeable sound. This sound corresponds with the frequency of the advancement. The intensity of the audio signal could be monitored and the locations of the peaks in the intensity may correlate with the frequency of the oscillations. If the analytical model of OFA could be advanced further, a hybrid analytical/experimental analysis may allow for the calculation of the heat flux from the flame to the sample surface for a given set of operating conditions. Because the heat flux is the driving mechanism in this phenomenon, it may be possible to relate it to the acceleration in the advancement frequency. Once the advancement frequency is known it might be able to be plugged into an analytical model to determine what heat flux is necessary to produce this frequency for a given set of conditions. The analytical model for this phenomenon is in the initial stages of development. 107 APPENDIX 108 Operating Instructions for The MSU Fire Tunnel 1.1 Introduction These instructions are designed to allow a user to operate the MSU Fire Tunnel who has no previous experience with it. It is assumed at the beginning of the instruction that the MSU Fire Tunnel is in a storage configuration. The MSU Fire Tunnel is in the storage configuration when the sample holder is out of the combustion chamber, the radiant panel is mounted onto the hood, and (obviously) all power is off and all of the components are unplugged. 1 .2 Performing Flame-Spread Tests with the MSU Fire Tunnel A fuel sample must be mounted into the sample holder. The back of the sample holder has nine holes for mounting screws. Pre-drill holes (appropriate for the mounting screws) approximately 1/4” into the fuel sample. Once the fuel sample is secured into place the edges should be filled. This can be done with interior spackling paste but it is not ideal. A filling compound that does not crack and expand when heated would be better. Fiberglass insulation can be used to fill the gaps that are not as critical to the flow over the sample holder. The back of the sample holder has four spacers and bolts. The bolts can be inserted into the holes on the combustion chamber and secured with washers and nuts. The radiant panel can now be mounted to the combustion chamber. The (safety) bolts used to mount the radiant panel to the hood are also used to prevent the panel from sliding out of position. Remove the radiant panel from its 109 position on the hood and slid it partially into the slots cut into the combustion chamber. Now remove the safety bolts and slide the panel in. Replace the safety bolts and secure with nuts. Imp rtant! Check each terminal connecting the radiant panel to the high- temperature lead wires. The wires leading to the terminals become loose with time and must be replaced periodically. The ignition wire system can now be installed. Insert the hollow ceramic tubes from the outside of the combustion chamber. Next, insert the igniter assembly into the ceramic tubes. The steel wires of the ignition assembly should stick out of the ceramic tubes about 1/4 inch. The leads from the variable voltage supply can now be attached to the ignition wire. Electrical tape (while not ideal) works well. The components can now begin to be plugged in and turned on. Plug in and tum on the blower (Fig. A.1), the pressure transducer power supply, and the voltmeter. Adjust the damper on the blower to the desired flow-rate (referring to the pressure transducer output). Mount the camera or other measuring equipment. Be sure that the butterfly valve at the top of the hood is open. Turn on the hood and shut the glass slide. Slide out the radiant panel. 110 ‘ Over-temperature sensor ' Power switch I ' for the blower Cut-off switch for - . ‘ ; the three I. l 5' - - 1 ' V' . - lines. __ 3 Power switch for the temperature ' , controller, over temperature sensor, and the 7‘h strip heater. Figure A.1 The radiant panel control system, component switches, and supporting equipment. Plug in the variable voltage supply and the power supply for the temperature controller. Turn on the power to the temperature controller. Note: The top strip heater is powered from the same circuit that powers the temperature controllers. When the controller is on that line is live! The SSRs are off when there is no control voltage. If the set temperature on the temperature controller is lower than the room temperature the SSR should not trip. However, 111 don’t maintenance the radiant panel if any power is on to the control system. Ground potential, thermocouple noise, ect. could cause the SSR to trip even if the set temp is less than the room temperature. Turn on the power to the power to the controller. Set the desired radiant panel temperature (°F) by holding down the * and simultaneously pressing the T or the 1. Flip the big disconnect switch into the on position. The radiant panel is now heating up. The controller will display the measured temperature. For more detail on the temperature controller see the CN9120A manual. The temperature is measured on the bottom of the radiant panel so the response is slow. The controller could be re-tuned (see manual) to improve the response time. Set the desired high-limit temperature. Usually it is set to around 100° F higher than the set temperature. Push the reset button. A green light should appean Important: If this light goes off it means that he temperature has exceeded the set temperature and the system should be shut down immediately. Slide the radiant panel into the combustion chamber once it has reached the desired temperature. Usually a sample heat-up time period is waited before the ignition wire is turned on. Turn on the ignition wire by dialing in a voltage on the variable voltage power supply. A calculation can be made to determine the voltage necessary to run the desired current through the wire. Tables on the relationship between the current through the wire and the temperature of the wire can be found in the Omega Temperature catalog. Usually the wire is heated until 112 glows a bright orange. Care must be taken so that the wire is not heated too much and broken. After the fuel sample is lit the ignition wire must be turned off in order to prevent it from breaking. Turn off the radiant panel when the test is complete. In order to extinguish the fire an injection of nitrogen is used. Open the main valve for the bottled nitrogen cylinder and open the pressure regulator valve. Turn off the blower. The fire should be extinguished in about ten to twenty seconds. Turn off and unplug all equipment. 113 BIBLIOGRAPHY 114 (1) (2) (3) (4) (5) (6) (7) (8) (9) BIBLIOGRAPHY Wichman, IS, “A Review of the Literature of Material Flammability, Combustion, and Toxicity Related to Transportation“, unpublished. ASTM E 84 — 97a, “Standard Test Method for Surface Burning Characteristics of Building Materials”, Annual Book of ASTM Standards, Vol. 4.07, 1997. ASTM E 1354 — 97, “Standard Test Method for Heat and Visible Smoke Release Rates for Materials and Products Using an Oxygen Consumption Calorimeter”, Annual Book of ASTM Standards, Vol. 4.07, 1997. ASTM E 1740 - 97a, “Standard Test Method for Determining the Heat Release Rate and Other Fire-Test-Response Characteristics of Wallcovering Composites Using a Cone Calorimeter", Annual Book of ASTM Standards, Vol. 4.07, 1997. ASTM E 648-97, “Standard Test Method for Critical Radiant Flux of Floor- Covering Systems Using a Radiant Heat Energy Source”, Annual Book of ASTM Standards, Vol. 4.07, 1997. ASTM E 1321 - 97a, “Standard Test Method for Determining Material Ignition and Flame Spread Properties”, Annual Book of ASTM Standards, Vol. 4.07, 1997. J.L. Cordova, J. Ceamanos and AC Femandez-Pello, “Piloted Ignition of a Radiatively Heated Solid Combustible Material in a Forced Oxidizer Flow”, 7th AlAA/ASME Joint Therrnophysics and Heat Transfer Conference, June 15—1 8 1998. J. Ouintiere, M. I-Iarkleroad and D. Walton, “Measurement of Material Flame Spread Properties,” Combustion Science and Technology, 1983, Vol32, pp. 67-89. J.G. Ouintiere, Principals of Fire Behavior, Delmar Publishers, 1998, p. 98 115 (10) (11) (12) (13) (14) (15) (16) J.L. Cordova, J. Ceamanos and AC Femandez-Pello and FLT. Long, J.L. Torero, J.G. Ouintiere, “Flow Effects on the Flammability Diagrams of Solid Fuels", Paper Proceedings of the 1997 4th International Microgravity Combustion Workshop, 1997. T., Morel, “Design of Two-Dimensional Wind Tunnel Contractions“, Journal of Fluids Engineering, Trans. ASME, Vol. 99, 1977, pp. 371 -378. Foss, J., Personal communication. Design Procedure For Optimum 2-D Wind Tunnel Contraction, B. Voss, 1990, Dos based computer program. Schlichting, H., Boundary Layer Theory, Pergamon Press, New York (1955). R. Siegel, J.R. Howell, Thermal Radiation Heat Transfer, Hemisphere Publishing Corp., 2"“' Edition, 1981. I. Miller, J.E. Freund, Probability and Statistics for Engineers, Prentince- Hall, lnc. 1965. 116 "Illll‘lllllllllllilies