Prediction of Plume Formation and Dispersion from Gas Flares

Results from three different applications are presented to illustrate how plumes from various types of flares have been analyzed with different levels of detail. The first application considers an enclosed flare plume. Results from a thermodynamic estimate of combustion temperature and species are used to calculate plume characteristics using buoyancy flux and Gaussian dispersion. The second application considers an air assisted flare and employees a steady state CFD analysis to investigate the impact flare tip geometry and a cross wind has on plume formation and dispersion. The third application considers a single tip ground flare and uses a transient LES based CFD analysis investigate the impact that a flare plume has on nearby equipment. Results from these applications provide examples of how flare plumes may be analyzed and the various types of information available from each approach.

Prediction of Plume Formation and Dispersion from Gas Flares Joseph D. Smith, Ph.D. Ahti Suo-Antilla, Ph.D., and Robert Jackson Systems Analyses & Solutions, Inc. Tulsa, OK Scot K. Smith, P.E. Zeeco, Inc., Broken Arrow, OK ABSTRACT Results from three different applications are presented to illustrate how plumes from various types of flares have been analyzed with different levels of detail. The first application considers an enclosed flare plume. Results from a thermodynamic estimate of combustion temperature and species are used to calculate plume characteristics using buoyancy flux and Gaussian dispersion. The second application considers an air assisted flare and employees a steady state CFD analysis to investigate the impact flare tip geometry and a cross wind has on plume formation and dispersion. The third application considers a single tip ground flare and uses a transient LES based CFD analysis investigate the impact that a flare plume has on nearby equipment. Results from these applications provide examples of how flare plumes may be analyzed and the various types of information available from each approach. INTRODUCTION AND BACKGROUND Plume formation from large ground flares, elevated flares, off shore flares, and other types of flares is a subject of interest due to potential safety hazards for workers and potential damage to surrounding buildings and equipment. Also, air pollutants are dispersed into the atmosphere downwind of flare stacks from the flare plume. Understanding how flare plumes form and are impacted by cross winds, how these plumes affect personnel and equipment, and how they disperse air pollutants are the main topics of this paper. An introduction to plume formation and how plumes are analyzed will first be discussed. Next, several examples of various flares and the associated plume analysis will be discussed. Lastly, a summary of plume analysis with conclusions on the appropriate analysis techniques will be given. Plume Formation and Dispersion Modeling Plumes from an elevated flare forms at the stack exit while plumes from a multi-tip ground flare form as plumes from individual burner tips merge into a large composite plume (Figure 1). Although plumes can be considered as starting as a point source subject to atmospheric conditions (i.e., stability class, wind conditions), detailed flare plume analysis must also consider the overall flare geometry. Questions of plume shape, rise height, and downwind concentration profiles must be answered. In general, taller flare stacks disperse pollutants better than shorter ones (i.e., ground flares verses elevated flares) because the resulting plume must disperse further before it reaches ground level as it spreads and disperses. For an elevated flare, flare gases mix with surrounding air and combust above the flare tip. The plume evolves from the visible flame and travels downwind of the flare stack. Due to the temperature difference between the flare plume and the surrounding air, buoyance causes the plume to rise and expand as it mixes with cooler surrounding air. As the plume rises and expands, pollutants (and heat) disperse into the surrounding atmosphere so pollutants have further to travel to reach ground level. The final "effective plume height (Hp) is the sum of the physical stack height (hs), the flame height (hf) and the average plume rise (Δh) as shown in Figure 1. Two factors mainly control plume rise: momentum and buoyancy. Plume momentum is determined by flare gas volumetric flow rate through the flare tip geometry. Plume buoyancy is determined by flare gas composition and resulting combustion products and effluent gas temperature. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 2 Figure 1- Flare plume formation and rise Plume momentum carries the flare effluent downwind from the stack to a point where atmospheric conditions control its trajectory. Once emitted, the initial stack exit velocity is quickly reduced by ambient air entrainment and wind effects which cause the plume to bend over as it acquires horizontal momentum. Since wind speed usually increases with distance above ground level, as plume elevation increases its horizontal momentum increases and it bends over until it appears level with the ground. Depending on ambient conditions, the point where the plume becomes horizontal may be a considerable distance from its origin. Buoyance mainly depends on the temperature difference between combustion products temperature (Tg) and the surrounding ambient air temperature (Ta). In an unstable atmosphere, the plume buoyance decreases as it rises while in neutral conditions the plume buoyancy remains constant with increasing elevation. Wind also reduces plume buoyancy by increasing the entrainment rate of cooler surrounding air into the plume. The greater the wind speed the faster the plume mixes with the atmosphere and cools thus losing its buoyancy. Thus, in "windy" regions a flare plume will not rise as high as in regions with little wind. FLARE PLUME DISPERSION ANALYSIS Flare plume dispersion models range in complexity from standard Gaussian dispersion models to transient reacting flow multi-physics CFD models. Each approach provides different amounts of dispersion information but also requires more input data and more effort to complete. The following illustrative examples demonstrate various levels of flare dispersion analysis for different types of gas flares. 1-D Plume Dispersion Analysis Flare plume dispersion may be analyzed using standard formulas for plume rise and distribution. The first example illustrates application of this approach. To set the stage for this analysis, the model is briefly discussed with application illustrated with an analysis of plume dispersion from an enclosed flare. The simplest form of plume dispersion model consists of a set of mathematical equations that predict pollutant concentrations downwind of the plume origin. These models use empirical equations to determine effective plume height and shape and estimate pollutant dispersion based on various meteorological conditions including temperatures, wind speeds, stabilities, and topography. A common model based on the plumeFlare StackΔhhfhsHpwindGround FlareTaTg Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 3 Gaussian dispersion equation estimates dispersion of nonreactive pollutants as a function of complexity of source and topography surrounding the flare. Plume Rise Model To calculate plume rise to estimate effective plume height, Brigg's standard plume rise formulas [1, 2, 3] relating "buoyancy flux" and wind speed to plume rise can be used: uxFhb3231)(6.1=Δ (1) Here, u is the average wind speed, x is the downwind distance from the plume origin (i.e., flare) and Fb is the "buoyancy flux" factor defined as: −=aagsgbTTTdgvF42 (2) where g is the acceleration due to gravity, vg is the gas velocity at the flare exit, Tg is the plume gas temperature and Ta is the ambient air temperature. Eq. 1 estimates plume rise for buoyancy dominated plumes. Some plumes are momentum dominated (i.e., high tip velocity from elevated flares). In this case, plume rise is more appropriately calculated from: 31223=ΔuxFhjmβ (3) where jβ is the jet entrainment coefficient defined as: gjvu+=31β (4) For Eq. 4, the "momentum flux" parameter is defined as: =gasgmTTdvF422 (5) with vg as the gas velocity at the flare exit and ds as the flare diameter (Ta and Tg defined as before). To determine whether a plume is buoyancy or momentum dominated, the critical temperature difference (ΔTc) must be determined. For unstable or neutral atmospheric conditions, the following equations are recommended: ()32310297.0sggcdvTT=Δ for 55<bF (6a) ()313200575.0sggcdvTT=Δfor 55≥bF (6b) When cTTΔ≥Δ the plume is buoyancy dominated and when cTTΔ<Δthe plume is momentum dominated. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 4 Plume Dispersion Model Source complexity (i.e., flare plume formed from a wake stabilized flame) is an important issue in determining the plume shape so the 1-D model must account for this issue. In addition, local topography and buildings play a key role in plume dispersion. Elevated plumes may impact airplane flight paths1 while elevated terrains (buildings) may affect plume formation and lead to higher surface concentrations since the predicted plume centerline is closer to the surface. The general Gaussian distribution equation is: +=  +−  −−  −2222121212zzyHzHzyzyeeeuQσσσσπσχ (7) where: χ = ground level pollutant concentration (g/m3) Q = mass emitted per unit time yσ = standard deviation of pollutant concentration in y (horizontal) direction zσ = standard deviation of pollutant concentration in z (vertical) direction u = wind speed y = distance in horizontal direction z = distance in vertical direction H = effective stack height As discussed elsewhere [4], the Gaussian distribution equation uses relatively simple calculations using two dispersion parameters (i.e. yσfor horizontal dispersion and zσ for vertical dispersion) to identify the variation of plume concentration away from plume center. Plume distribution predicted by Eq. 7 determines ground level plume concentration using time-averaged atmospheric temperature and wind speed (and direction). This approach limits predicted plume concentration downwind of the stack in terms of normally distributed "average" concentration as illustrated in Figure 2: Figure 2 - Normal distribution of concentration for a Gaussian plume dispersion analysis Gaussian based plume analysis predicts plume rise and effective size as a function of atmospheric stability using vertical and horizontal dispersion coefficients, which are functions of wind, cloud cover and 1 For example, in Tulsa, OK the flight path for the Tulsa International airport passes directly over the John Zink flare test facility. When flares are tested, John Zink must inform the airport of their testing to avoid potential hazards to low flying aircraft on their final approach. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 5 solar radiation. This approach assumes: 1) relatively flat surface along plume path, 2) plume reflected from surface, 3) constant plume emission rate and 4) uniform wind speed and direction. In this case, the horizontal and vertical dispersion coefficients are obtained graphically as a function of atmospheric stability2 with each increasing with decreasing atmospheric stability. One-dimensional plume rise formulas discussed previously calculate the imaginary centerline of a flare plume trajectory where the maximum concentration exists and assumes a uniform distribution from the centerline. However, these simple one-dimensional formulas are not able to estimate how pollutant concentrations vary non-uniformly from the centerline path. Gaussian dispersion predicts pollutant concentrations at points of interest away from the source (i.e., downwind ground level, surface of an elevated building platform, etc.). Existing plume models developed by the US Environmental Protection Agency (USEPA) are available to perform this type of analysis. [5] Application: Enclosed Flare Plume Formation and Dispersion To illustrate the Gaussian approach for plume analysis, the formation and dispersion of a plume from the enclosed flare shown in Figure 3 was analyzed. The flare system included up to 3 units similar to the one shown in Figure 3 with the system processing up to 37,000 kg/hr flare gas with a molecular weight of approximately 25 containing saturated and unsaturated hydrocarbons plus halogens with heavy metal contaminates and a Lower Heating Value between 1200 and 6200 KJ/m3. The analysis estimated the initial plume conditions using the NASA-Glenn Chemical Equilibrium with Applications code (NASAcea2) to calculate expected plume composition at various stoichiometries and temperatures [6]. This information was then used to determine the flare tip gas exit velocity. For this plume analysis, the objective was to estimate ground level concentration of certain emissions at a predetermined location to support the air permit application process. Using NASAcea2, the combustion effluent composition with solid particulates and the expected gas temperature as a function of combustion stoichiometry were calculated. Given operating conditions and flare diameter, the flare tip velocity was estimated (see Figure 4). With the tip velocity, effluent composition, and the atmospheric stability, plume dispersion was estimated for the downwind location. Figure 3 - Enclosed flare with burner system 2 Atmospheric stability depends on temperature difference between plume and surrounding air. Stability levels can be stable, conditionally stable, neutral, conditionally unstable, or unstable which for plume modeling are classified by five surface wind speed categories, three types of daytime insolation, and two types of nighttime cloudiness and are referred to as Pasquill-Gifford stability classes. Stabilities A, B, and C refer to daytime hours with unstable conditions while stability D represents overcast days or nights with neutral conditions and stabilities E and F refer to nighttime, stable conditions based on amount of cloud cover. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 6 Figure 4 - Estimated Flare Tip Velocity as a function of gas temperature and composition Results shown in Table 1 include estimated levels of SO2, NO2, CO, Total Solid Particulate (TPS) and concentration of particles with mean diameters less than 2.5 microns (PM2.5). This information was then reported to the client and used as part of their air permit application for the enclosed flare system. Table 1 - Estimated Flare Plume Dispersion 3-D Steady State Plume Dispersion Model To more accurately predict plume formation, rise and dispersion, computational fluid dynamics (CFD) can be used. To analyze the steady state plume dispersion, the CFD code Fluent has been used. This code describes the reactive turbulent flow associated with the combustion of flare gas and the subsequent formation and dispersion of the flare plume. A detailed discussion of the CFD code is not given here but a summary is given for comparison of the complexity of the CFD approach compared to the Gaussian dispersion approach. Steady State RANS Model In general, the general Reynolds Averaged Navier-Stokes (RANS) steady state equations describing the conservation of mass, momentum, and energy written in Cartesian tensor form are [7, 8]: ∂ρ ˜ u j()∂xj=0; mass conservation (8) ∂ρ ˜ u i˜ u j()∂xj=−∂P ∂xi+∂∂xjτ ij−ρ′ ′ u i′ ′ u j()+ρ fi; momentum conservation (9) ∂ρ Φ uj()∂xj=∂∂xjΓ φ∂Φ ∂xj−ρ′ ′ φ ′ ′ u j      +S Φ; general conserved scaler Φ (i.e., energy) (10) where: τij = viscous stress tensor Φ = transported variable (chemical species, enthalpy, etc.) ΓΦ = molecular diffusion coefficient for Φ 4.004.505.005.506.006.507.007.50800100012001400160018002000Temperature (K)Velocity (m/s)Eq.Rat=0.5Eq.Rat=0.75Eq.Rat=1.0 Base Elevation (m)Stack Height (m)Stack Diameter (m)Exit TemperatureExit VelocityStack Capped or Not X (m)Y (m)(H)(H)(K)(m/s)(Yes/Not)SO2 (g/s)NO2 (g/s)COa (g/s)TSP (g/s)PM2.5b (g/s) 007.0016.763.511144.267.10Not0.142140.0406252.44E-050.03950.027 007.0016.763.51939.2610.94Not0.0895570.04830200.0520.036 007.0016.763.51693.717.17Not0.0007740.01229900.0260.018 007.0016.763.51916.4811.01Not2.6246590.0408300259180 007.0016.763.51682.047.16Not0.0198210.010775013090Source CoordinatesPollutants Emission Rates at Normal Operating Conditions Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 7 SΦ = source term for Φ ρ = density of the fluid P = static pressure g = acceleration of gravity Individual species transport is described by: iiiiiYYDYuYt+∇⋅∇=⋅∇+∂∂)()(ρρρ (11) where the species source term, iY accounts for the consumption and production of each specie described by: kKkikiiiirWYWtcΣ===1ν∂∂ (12) together with the total heat released during combustion which is written as the sum of the product of reaction enthalpy "h" and concentration change for the species: qqHdcdthrkkKkikikkKiiNi==−⋅=−===ΣΣΣ111Δνω, (13) Turbulent pre-mixed combustion is simulated using either the Eddy-Breakup model (EBU)[9] or the more general Eddy-Dissipation Concept (EDC). [10] Several turbulence models have been proposed [11], but the ε−k turbulence model, originally proposed by Harlow and Nakayama [12] remains the most widely used to describe practical flow systems [13]. Using the modified Boussinesq hypothesis [14] the transport equations for the standard k-ε model are: ∂∂+∂∂+−∂∂=∂∂jkTjjiijjjxkxxUxkUσμμρετρ: Turbulence Kinetic Energy (14) ∂∂+∂∂+−∂∂=∂∂jkTjjiijjjxxkCxUkCxUεσμμερτεερεε221: Dissipation Rate (15) ερμμ2kCT=: Eddy Viscosity (16) with the empirical closure coefficients included in Eqs. 14-16 taken as: 3.1,0.1,09.0,92.1,44.121=====εμεεσσkCCC (17) Equations 8 through 16 include several steady-state, second-order, non-linear, elliptic partial differential equations (PDE). Solving these continuous PDE's requires that they be transformed into discrete finite difference equations (FDE's). The set of FDE's are then solved using various numerical schemes (i.e., SIMPLE = Semi-Implicit Method for solving Pressure Linked Equations) to obtain a "converged" solution for the flow, temperature, and concentration profiles of an operating flare and its resulting plume. Application: Air-Flare Tip Design Impact on Plume Shape A recent analysis of an air-assisted gas flare illustrates the use of a RANS based CFD code to predict flare plume formation and dispersion. Typically, air-assisted flares employ complex geometries to enhance Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 8 mixing between the flare gas and surrounding air. The present air-assisted flare included alternating air and gas tubes (see Figure 5) through which approximately 170,000 lbs/hr flare gas and 29,000 CFM is flowing at 60ºF. The plume formed with this condition was predicted using the Fluent CFD model for no-wind and a 3-20 mph wind (see Figure 6). As shown, no wind (or low wind) has little effect on this plume while at higher wind velocities (~20 mph) the plume clearly is deflected and begins to form the expect barrel vortices which are commonly observed in plumes as they rise from a stack into the wind above the stack. Recalling the Assist Air Tubes (hydrocarbon flows around air tubes)Flare Tip GeometryFlare Gas (red) at Tip Exit Figure 5 - Air-assisted flare tip geometry Flare Exit3'6'9.5'12'Flare ExitFlare ExitNo Wind3 mph Wind20 mph Wind Figure 6 - Wind effect on air flare plume formation earlier discussion about plume dynamics and the relationship between stack exit velocity and cross wind velocity, this behavior is not unexpected. As shown at no or low wind conditions, plume momentum dominates as the plume rises well above the stack exit and maintains a cylindrical shape. As wind velocity increases, the outer regions of the plume are deformed and the shape changes into a horse shoe shape with vortices forming on the outer edges perpendicular to the wind. Plume rise is reduced as the wind bends the plume over and deforms its shape. This example shows the limitation of a Gaussian plume dispersion analysis which predicts an average plume path with a constant concentration along this path. 3-D Transient Plume Dispersion Model An even more completed modeling approach utilizes an LES based turbulent reacting flow CFD model called C3D. This tool has previously been applied to analyze plumes from pool fires [15, 16, 17] and gas flares. [18, 19] To simulate flare plumes, C3D uses an LES based turbulent mixing model. Turbulent mixing sets the combustion rate and controls plume mixing with surrounding air. The energy equation combines combustion with turbulent mixing and radiation heat transfer. Individual species transport equations predict the distribution and concentration of fuel, oxygen, intermediate species, soot, and combustion products. The combustion model (described below) provides source and sink terms as a function of species concentrations, local gas temperature, and turbulent diffusivity. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 9 The radiation model estimates flame emissivity as a function of molecular gas composition, soot volume fraction, flame size, shape and temperature distribution. Radiation flux between the flame and external surfaces is predicted assuming the internal flame is optically thick so radiation transport is diffuse. Diffuse radiation is modeled employing an effective, temperature-dependent flame thermal conductivity. Outside the flame (defined as local soot volume fraction < specified minimum soot volume fraction) the gaseous medium does not participate in radiation heat transfer. The radiation model uses view factors between all points on the flame surface and all points on surrounding objects in the computational domain3. Solid objects in the domain also radiate to the environment when they are not engulfed in flames. The combustion scheme used is based on a proposed hybrid EBU method. [20] Species included in this hybrid combustion model include fuel vapor (F) from the flare tip, oxygen (O2), products of combustion (PC) including water vapor and carbon dioxide, radiating carbon soot (C), and non-radiating intermediate species (IS). The general combustion reactions involving these species include: 1kg F + (2.87-2.6S1) kg O2 → S1 kg C + (3.87-3.6S1) kg PC + (50-32S1) MJ (18) which describes incomplete combustion of fuel vapor (F) and produces carbon soot (C) and products of combustion (PC) plus energy. The standard combustion soot stoichiometric parameter (S1) is set to 0.05 but can be adjusted based on fuel type. For natural gas a value of 0.005 is used. The endothermic fuel pyrolysis or cracking reaction (soot producing) consumes fuel vapor (F) and energy and produces radiating carbon soot (C) plus the intermediate species (IS): 1kg F + 0.3 MJ → S2 kg C + (1-S2) kg IS (19) This step includes the cracking parameter (S2) which is set as 0.15 but can also be adjusted based on fuel type. Finally, soot combustion is described by: 1kg C + 2.6 kg O2 → 3.6 kg CO2 + 32 MJ (20) which consumes carbon soot (C) and oxygen (O2) and produces carbon dioxide (CO2) plus some energy. Combustion of intermediate species (IS) is described by: MJPCOIS2222222S-1S3250 kgS-1S 6.387.3 kgS1S6.287.2 kg 1−+−→−−+ (21) where the coefficients are selected so that complete combustion of carbon soot (C) and intermediate species (IS) produce the same species and thermal energy as direct combustion of the fuel. The coefficients in the formula are in terms of mass weight and not moles. The advantage of the three-step reaction is that the first reaction has a low activation energy, which allows the partial burning and heat release of the flare gas. This maintains combustion since the partial heat released allows the second reaction, which produces most of the heat and all of the soot, to proceed. The flare gas Arrhenius combustion time scale is combined with the turbulence time scale to yield an overall reaction time scale for the reaction: diffturbxCετ2Δ= (22) 3 The view factor calculation does not include shadowing effects from irregular fire surfaces. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 10 where Δx is the characteristic cell size, C is a user input constant (0.2E-4), εdiff is the eddy diffusivity from the turbulence model, and τturb is the turbulence time scale (i.e. characteristic time required to mix contents in computational cell). The reaction rates are combined by simple addition of the time scales. All reactions are based on Arrhenius kinetics: TTNiRRAefCdtdf−−=Π11 (22) where the "C" coefficients and the Activation Temperatures, TA, are supplied for each reaction. This hybrid combustion approach efficiently simulates turbulent hydrocarbon combustion to estimate the flare plume composition and then uses an LES based mixing model to simulate plume transport. Application: Single Tip Ground Flare Plume Impact on Surrounding Equipment Sometimes flare plumes impact surrounding equipment and may pose safety hazardous to operations staff. The next application considers a single tip ground flare tip located 2 meters above grade and 23 meters below an overhead cable. The problem considered was whether the flare plume temperature would exceed the safe operating temperature (250ºC) of the nearby equipment (see Figure 7). To address this question, a transient LES based combustion analysis of the flare was completed using the C3D code describe earlier. 23 m Figure 7 - Single point ground flare located below cable The ground flare burned ethylene flowing at 200 Nm3/hr at 20ºC. The combustion model used for this analysis included three reaction steps: The computational domain and respective mesh structure are shown in Figure 8: Soot MJ/Kg32CO 3.66O 2.66SootteIntermedia MJ/Kg34.1 Soot0.25OH 0.7CO 2.7O 2.58HCEthylene MJ/Kg0.94OH 0.64HC 0.93O 0.57HC2222222222242+→++++→+++→+ Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 11 Cable 30 m above gradeFlare tip80 m80 m60 m23 mGrid horizontal slice through cable Figure 8 - C3D computational domain and mesh structure As shown above, the mesh is refined near the cable and around the flare tip. This allows improved resolution of the flare plume and better prediction of temperature profile. Performing transient LES analysis of the flare system using the 3-step combustion model on the computational mesh provided estimates of plume temperature up to 50 meters above the flare tip (see Figure 9). The elevation shown in Figure 9 represents total elevation above ground level. Given the flare tip was located approximately 2 meters above grade, the actual cable elevation is 25 meters above grade. Temperature (°C)45 m elevMax Temp: 120 °C40 m elevMax Temp: 165 °C50 m elevMax Temp: 100 °C35 m elevMax Temp: 235 °C30 m elevMax Temp: 260°C25 m elevMax Temp: 340°C Figure 9 - Predicted gas temperature at various elevations above flare tip Since C3D provides a transient analysis of the flare plume, a time averaged calculated temperature profile is possible (see Figure 10). This clearly shows that the plume temperature at the cable elevation is 221ºC but actually ranges between about 100ºC and 350ºC. This type of plume analysis also provides the same level of detail for gas species which allows one to assess O2 and CO levels around nearby buildings. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 12 Mean temperature221ºC171ºC140ºC110ºC89ºC62ºC50 m62ºC25 m221ºCCenterline Temperature at 50 sCenterline Plume Temperature averaged over 50 s Figure 10 - Time averaged plume temperature over a 50s period SUMMARY AND CONCLUSIONS Flare plume analysis is an important topic that must be addressed when designing a flare system located in an operating plant. This paper presented a summary of three different types of dispersion models used to analyze flare plumes. Depending on the amount and type of information required, the selected model may range from a 1-dimensional Gaussian dispersion model to a 3-dimensional steady state RANS based CFD model to a 3-dimensional LES based transient combustion model. To estimate the concentration of flare emissions downwind of a flare stack required for environmental air permits, a simple Gaussian dispersion analysis may be appropriate. This may be appropriate when complex chemistry and solid particulates are involved. The Gaussian dispersion model predicts the average plume path and assumes a constant source location with uniformly distributed dispersion from the plume along its path. As shown, a flare plume can be dramatically affected by windy conditions and flare tip geometry. In this case, a steady state CFD model was used to estimate plume shape. Results of this analysis showed that the flare plume may not have a uniformly distributed dispersion profile along its path. In this case, the plume shape changed as it transitioned from a momentum dominated plume to a buoyancy dominated plume. The plume momentum was also shown to be a strong function of flare tip geometry and assist media flow rate. This is important as the plume dynamics will change as a function of operating condition for a given tip geometry. Thus, for an accurate estimate of plume dispersion downwind of the flare a detailed CFD plume analysis should be conducted. In some cases, plume dynamics are also critical when evaluating potential hazards to nearby equipment and operations personnel. For the case shown, the average predicted temperature was below the critical operating temperature but dynamic fluctuations in plume temperature exceeded the safe operating condition for the equipment. Thus, when considering safety issues related to flare plume dispersion, a full transient LES based CFD plume dispersion analysis will be required. Flare Plume Formation and Dispersion 2012 Annual American Flame Research Committee Meeting Salt Lake City, Utah September 5 - 7, 2012 13 REFERENCES [1] Briggs, G.A., "A Plume Rise Model Compared with Observations", JAPCA, 15:433-438 (1965) [2] Briggs, G.A., "Plume Rise", USAEC Critical Review Series (1969) [3] Briggs, G.A., "Discussion: Chimney Plumes in Neutral and Stable Surroundings", Atmos. Envir., 6:507-510 (1972) [4] Irwin, J.S., "Estimating Plume Dispersion: A Comparison of Several Sigma Schemes", J. Appl. Meteorol., 22:1 (1983) [5] User's guidelines for EPA SCREEN3 and AERSCREEN dispersion model: http://www.epa.gov/ttn/scram/dispersion_screening.htm) [6] http://www.grc.nasa.gov/WWW/CEAWeb/ [7] Bird, R.B., Stewart, W.E., and Lightfoot, E.N. Transport Phenomena. John Wiley & Sons, New York, NY (1960) [8] Anderson, J. D., Computational Fluid Dynamics-The Basics With Applications. McGraw-Hill, NY (1995) [9] Spalding, D.B. "Mixing and chemical reaction in steady confined turbulent flames," 13th Symp. (Int.) on Combustion, 649-657 (1971) [10] Magnussen, B.F.; Hjertager, B.H.; Olsen, J.G.; Bhaduri, D., "Effects of turbulent structure and local concentrations on soot formation and combustion in C2H2 diffusion flames," 17th Symp. (Int.) on Combustion, 1383 - 1393 (1978) [11] Nallasamy, M. "A Critical Evaluation of Various Turbulence Models as Applied to Internal Fluid Flows," NASA TP - 2474, Marshall Space Flight Center, Alabama, May (1985) [12] Harlow, F.H. and Nakayama, P.I. "Transport of Turbulence Energy Decay Rate," Rpt. # LA-3854, Los Alamos Scientific Laboratory, Los Alamos, NM, January (1968) [13] Speziale, C.G. "On Nonlinear k-l and k-ε Models of Turbulence," J. Fluid Mech., 178, 459 (1987) [14] Smoot, L.D. and Smith, P.J. Coal Combustion and Gasification. Plenum Press, New York, NY (1985) [15] Greiner, M, Are, N., Lopez, C., and Suo-Anttila, A., "Effect of Small Long-Duration Fires on a Spent Nuclear Fuel Transport Package," Institute of Nuclear Materials Management 45th Annual Meeting, Orlando, FL, July 18-22 (2004) [16] Greiner, M., and Suo-Anttila, A., "Radiation Heat Transfer and Reaction Chemistry Models for Risk Assessment Compatible Fire Simulations," Journal of Fire Protection Engineering, 16, 79-103 (2006) [17] Fuss S.P., Hamins, A., "An estimate of the correction applied to radiant flame measurements due to attenuation by atmospheric CO2 and H2O", Fire Safety Journal, 37, 181-190 (2002) [18] Suo-Anttila, A., Smith, J.D., "Application of ISIS Computer Code to Gas Flares Under Varying Wind Conditions," 2006 American Flame Research Committee International Symposium, Houston, TX, October 16-18 (2006) [19] Smith, J.D., Suo-Anttila, A., Smith, S.K., and Modi, J., "Evaluation of the Air-Demand, Flame Height, and Radiation Load on the Wind Fence of a Low-Profile Flare Using ISIS-3D," AFRC-JFRC 2007 Joint International Combustion Symposium, Marriott Waikoloa Beach Resort, Hawaii, October 21-24, (2007) [20] Said, R., Garo, A., and Borghi, R., "Soot Formation Modeling for Turbulent Flames," Combustion and Flame, 108, 71-86 (1997)