The invention relates to a method according to the preamble of claim 1 for determining the flame shape of a swirling flame in a closed combustion chamber, using third-octave band acoustic analysis.
Pollutant emission regulations for heat engines, boilers, and combustion systems are continuously stringent. Nowadays, these requirements for steady-operating combustion applications are often fulfilled by using a lean mixture, using burners developed explicitly for such operation. It is known from Glassman, Yetter: “Combustion” (Burlington: Academic Press, 2008) that the lean flammability limit in terms of air-to-fuel equivalence ratio is typically the double of the stoichiometric mixture, showing marginal variations for hydrocarbon fuels. Nowadays, the most important pollutants are nitrogen oxides (NOx), which concentration in the flue gas is mitigated by using a leaner mixture. It is known from Dunn-Rankin “Lean Combustion: Technology and Control” (Academic Press, 2007.) that if there is excess oxygen in the combustion zone compared to the required for complete combustion, the flame temperature will be lower, ultimately leading to a reduced NOx emission. The leaner mixture can be characterized with a reduced flame propagation speed, however, to provide a suitably wide operating range, the flame is usually stabilized by swirl vanes. Swirling flames can be classified into strongly and weakly swirling flames, depending on the swirl intensity, discussed by Beer, Chigier “Combustion aerodynamics” (London, Robert E. Krieger Publishing Company, Inc. 1972.). While the flame shape is helical at low swirl, it forms a V shape when the swirl intensity is high. There is a transitory regime between the mentioned two, in which both straight and V-shaped flames can be observed, see Jozsa, Kun-Balog “Stability and emission analysis of crude rapeseed oil combustion” “Fuel Processing Technology”, (156, pp. 204-210. 2017.). The mentioned flame shapes are qualitatively shown in
The preferred flame shape is the V since it is characterized by a notably reduced NOx emission. Consequently, there are numerous journal papers available in the public literature dealing with such flame shape, such as Stair, Boxx, Carter, Meier “Experimental study of vortex-flame interaction in a gas turbine model combustor” (Combustion and Flame, 159(8), pp. 2636-2649. 2012.); Durox, Moeck, Bourgouin, Morenton, Viallon, Schuller, Candel: “Flame dynamics of a variable swirl number system and instability control” (Combustion and Flame, 160(9), pp. 1729-1742. 2013.); Noiray, Denisov: “A method to identify thermoacoustic growth rates in combustion chambers from dynamic pressure time series” (Proceedings of the Combustion Institute, 36(3), pp. 3843-3850. 2017.); Idahosa, Basu, Miglani: “System Level Analysis of Acoustically Forced Nonpremixed Swirling Flames” (Journal of Thermal Science and Engineering Applications, 6(3), pp. 1-15. 2014.)
The flame shape is critical from the point of view of the operation of an engine; inappropriate design may lead to severe damage. Hence, the requirement for time-to-time or online diagnostics of the flame shape is becoming increasingly important. Such methods are used explicitly for low reaction time flame diagnostics or control.
It is known that various flame shapes feature distinct acoustic characteristics. Based on the requirements above, acoustical and mathematical investigation of flame shapes and analysis of combustion, finding the governing and modifying phenomena by control has become a highlighted research area in the past decades. This is supported by the following comprehensive publications: Huang, Y., Yang, V. “Dynamics and stability of lean-premixed swirl-stabilized combustion” (Progress in Energy and Combustion Science, 35(4), pp. 293-364. 2009.), Lieuwen, T. “Unsteady Combustor Physics” (New York, N.Y.: Cambridge University Press. 2012.), and Reed, R. J. “North American Combustion Handbook: A Basic Reference on the Art and Science of Industrial Heating with Gaseous and Liquid Fuels” (Claveland, Ohio: North American Mfg. Co. 1986.). The known flame acoustic analysis methods principally evaluate the Fourier-transform of the pressure signal, however, the spectral data is not processed further. Hence, they do not realize the goal of the present method. Nevertheless, time series analysis is another known method in the literature; this method similarly fails to provide the goals discussed earlier.
The characteristic frequencies of the combustion chamber can be well-localized; the mentioned book of Reed tells that the phenomenon of turbulent combustion can be described by broadband noise in the spectrum. A typical acoustic spectrum originated from combustion is presented in
Our goal with this invention is to make simpler, more robust, and accurate determination of the identification of notable flame shapes during combustion possible.
We have recognized that in the case of a steady-operating combustion chamber, in which all three flame shapes may occur, there are such characteristic frequencies, which helps with apparently distinguishing the three flame shapes by calculating the ratio of the sound pressure levels of the nearby band center frequencies. The frequency analysis can be performed online by the currently available microphones and software.
Our goal, namely, determining the ratio of the sound pressure levels of the band center frequencies has been reached by a method according to claim 1. Several main advantageous realizations of the method are disclosed in dependent claims.
A main advantage of the method according to the invention is that a third-octave analysis can be performed online with common devices, hence microphones and widespread software. Since these frequency bands are sufficiently wide, robust diagnostics can be performed as the method is insensitive to the temporal evolution of frequency peaks. The result can be directly determined by the method according to the invention, there is no need for further technology or method.
Additional features and advantages are described herein and will be apparent from the following detailed description and the figures wherein
It can be observed that the peaks are weakly localized in the spectral regimes, and the changing of the acoustic impedance also means temporal variation, described by, e.g., Kabiraj, L., Sujith, R. I. “Nonlinear self-excited thermoacoustic oscillations: intermittency and flame blowout” (Journal of Fluid Mechanics, August 2015, pp. 1-22. 2012) and Sampath, R., Chakravarthy, S. R. “Investigation of intermittent oscillations in a premixed dump combustor using time-resolved particle image velocimetry” (Combustion and Flame, 172, pp. 309-325. 2016).
These public sources are also supporting that the researchers try to make conclusions on flame shapes based on the detailed frequency regime, which makes the above-detailed excessive effort essential, and confines the available information for further processing.
Taking apart from this well-established, common practice a new method has been developed in which sound pressure levels of various frequency intervals are investigated at various spectral bands instead of using a detailed frequency resolution. Third-octave analysis is a possible method, which, however, contains frequency resolution detailed in the ISO 18405:2017 standard, its intentional use in practical combustion was not published according to our best knowledge, and also, its use was not discussed anywhere as a possibility.
In the following example, there were five third-octave bands identified, originated from measurement results for maximum efficiency. However, if the spectrum of a similar burner to the subject one is known, then a prediction can also be performed, which might work without correction based on the experimental results, but it is strongly advised.
A lean premixed prevaporized burner is analyzed as an example to the application of the present invention, in which the combustion air inlet had a tangential component, making the flow swirling. For the analysis, there must be at least one sensor put into the combustion chamber, which is capable of sensing combustion noise, hence, the acoustic fluctuations created by combustion. Such a sensor can be, e.g., a pressure sensor or a microphone, which is sensitive in the above-detailed spectral range; hence, its output signal can be used for further processing.
The sampling frequency of the microphone, used as a pressure sensor, should be at least double of the largest frequency component, like in the present case, according to the Nyquist-Shannon sampling theorem; however, a factor of three to five is recommended in practice. Consequently, using at least 10 kHz sampling frequency is advised. It should be noted that most of the commercially available noise analyzer systems operate at 20 kHz by default; hence, the application of this method can be performed by an expert in a familiar environment; it does not require special, expensive technology. Combustion noise at lower frequencies is not compact; consequently, it cannot be assumed as a point source. However, the location of the microphone is limited by practical reasons as the device may defect at high temperatures. To enhance its thermal resistance capability, a cooled sensor socket can be used if the operating temperature given by the manufacturer of the microphone cannot be met by proper placement. In the presented example, the microphone was placed at the height of the burner lip, 1 m sideways. Nevertheless, this is not possible in the case of an industrial application. Hence, by knowing the propagation of the noise, one should select a sensor with the right sensitivity. The noise inside the flame maybe 150-170 dB in the case of industrial burners, while this is 50-70 dB in the case of domestic appliances. The decrease of noise intensity is quadratic as the function of the distance measured from the source in the case of free noise propagation. Fundamentally, combustion noise shows low directional dependence; hence, the placement of the sensor is nearly arbitrary. The data are evaluated by the aforementioned third-octave method, provided by the ISO 18405:2017 standard.
The measurement can be performed from time-to-time, however, for maximum efficiency, continuous data evaluation is recommended.
To identify the characteristic frequency bands, a few different operating conditions are set to analyze the flame spectrum within the boundaries of physically possible parameter ranges given by the burner. The system incorporating the burner determine a minimum and maximum airflow rate, the fuel system sets an upper limit of thermal power, and the air delivery system with the swirl vanes limits the range of possible swirl numbers. During our investigation, the notable regimes are selected, and the analysis is performed on these.
The various frequency bands are respective to the burner design. For instance, in the case of the burner taken as an example, the third-octave band center frequencies were 200, 250, 400, 500, and 3150 Hz. However, these frequencies depend on burner design, similar values were expected, based on, e.g., Singh, A. V. et al. “Investigation of noise radiation from a swirl stabilized diffusion flame with an array of microphones” and Candel, S., Durox, D., Schuller, T., Bourgouin, J.-F., Moeck, J. P. “Dynamics of Swirling Flames” (Annual Review of Fluid Mechanics, 46(1), pp. 147-173. 2014).
The sound pressure level ratio of straight and V-shaped flames notably differs. It was shown in
It can be seen that the sound pressure level ratios are nearly identical. A similar trend characterized V-shaped flames, while there is a continuous transition in the transitory state. Based on this ratio, we can properly determine the flame shape and use this information for creating control algorithms if the combustion system characteristics are now.
The method can be used for measurement by the third-octave method, and the results can be stored within the physically possible operation range of a burner, where combustion can be sustained. Typically, we can go beyond the factory limitations, which might mean less favorable operation. Based on the diagram of
It should be noted that even the name of the third-octave method refers to the fact that the spectrum is analyzed in bands, more precisely, in logarithmic order. Hence, the spectral resolution is not perfect; the sound pressure level is corresponding to a broader band. However, combustion is a broadband phenomenon. Therefore this simplification does not lead to the loss of information. The first band center frequency that meets the above-detailed condition is 200 Hz, 500 Hz, and 6300 Hz for V-shaped flames, as it is known from the ISO 18405:2017 standard. Based on our combustion-related experiences, the last identified peak is an outlier since there is no such characteristic frequency as all the significant frequency components were below 6 kHz. Consequently, the third result was omitted. Nevertheless, precision would require the use of 6.3 kHz; the 300 Hz difference does not mean deviation in practice due to the logarithmic nature of the third-octave frequencies. The 6 kHz limitation was determined by the Fourier-transform technique since only background noise can be found from this frequency up to 20 kHz. In conclusion, there is no need to check higher frequency ranges.
The first local peak for straight flames is located at 500 Hz, then at 3150 Hz, which is also the global maximum. Since the spectra of both the straight and the V-shaped flame appears in the spectra of transitory flames hence, the following frequency peaks were found: 500 Hz, which was already present for both straight and V-shaped flames, and 3150 Hz, which was the global peak of the straight flame. Also, there was a local maximum located at 16 kHz, however, it was omitted due to the above reasons.
The method presented in the above example can be used for all spectra. It is possible that both the number and frequency of the detected notable peaks vary within a single flame shape as the operating conditions change. If there are technically relevant criteria for the burner to meet, such as pollutant emission, these data should be evaluated along with the spectra. The flame shape should be generally understood as a parameter range in which there is no sudden jump in the characterizing properties, meaning even favorable or unfavorable conditions. The sudden jump can be defined here as a 10% change in a single parameter on a relative scale. Following this, it is possible to find multiple characteristic regimes in the case of a single combustion chamber for, e.g., uniformly V-shaped flames. It should be noted that this method also works for combustion chambers that do not feature V-shaped flame.
It should be highlighted that a flame can be of any physical shape; the distinction between flame shapes should be made based on their behavior since their visible appearance might be similar. Therefore, it might be beneficial to provide optical access to the flame during the measurement procedure to see the flame structure; having such access is not mandatory but greatly helps the processing of measurement data. Finally, the number of the monitored frequencies and the relevant frequency ratios should be determined in a way to have a correlation with the notable characteristic parameter or parameters of the flame. For instance, these can be pollutant emissions or system efficiency.
A person skilled in the art may gather the same information by using a technique other than the above-detailed third-octave method. The information content of the evolution of the frequency ratios may be derived by using another spectral resolution technique, consequently, the third-octave approach is only a possible solution that can be used. In other cases, an octave-based, coarser approach and finer spectral resolution than that offered by the third-octave method can also lead to success.
A further possible approach, different from the detailed above, is the adaptive parameter analysis by using Artificial Intelligence. During this process, the system response to the change of the substantial parameters is investigated. Based on the results, a ‘learning database’ is created, then the correlation model is tested on data outside of the learning database. By using this approach, the same result can be achieved as detailed above, not necessarily relying on acoustical data. This technique is usually less favorable in industrial applications since its operation is not transparent, and the generated database by machine learning is not well-defined and not searchable for the operator. The statistical nature of this technique usually means an excessive risk for critical areas; hence, they are seldom used.
Number | Date | Country | Kind |
---|---|---|---|
P2000063 | Feb 2020 | HU | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/HU2021/050010 | 2/5/2021 | WO |