System and method for determining the amplitude of oscillatory instabilities in fluid mechanical devices

Abstract

Embodiments herein provide a system (100) to estimate the amplitude of oscillations in a turbulent flow system (102) that exhibits oscillatory instabilities. The system (100) comprises of a sensor (102A) mounted on the turbulent flow system (102) to detect an oscillatory variable in the system obtaining a signal, a signal conditioner (104) that conditions the signal from the sensor, an amplitude estimator (110) that estimates the amplitude of the limit cycle oscillations, and also predict the proximity of the system to the oscillatory instability, a processor (108) connected to the amplitude estimator (110) to compare the predicted oscillation amplitude with a threshold value. The amplitude is estimated by estimating the spectral measure of the time series signal obtained from the system.

Claims

1. A system (100) to determine amplitude of oscillatory instability in a turbulent flow system (102), the system (100) comprising: a processor (108); a sensor (102A) provided in communication with the processor (108), configured to obtain at least one signal corresponding to the impending oscillatory instabilities of the turbulent flow system (102); an amplitude estimator (110) provided in communication with the sensor (102A) and the processor (108), wherein the amplitude estimator (110) includes a spectral measure estimator provided in communication with the processor (108), wherein the spectral measure estimator is configured to convert the at least one signal into a frequency spectrum, the amplitude estimator (110) configured to: estimate an amplitude of the impending oscillatory instabilities by estimating the value of spectral measure of the at least one signal, wherein the spectral measures are calculated based on following equation: $M_{f_{0,\Delta f,f_s}}^{m,n,x,y}=\frac{{\left(\underset{i=-\Delta f/f_s}{\overset{\Delta f/f_s}{.Math.}}A\left(f_0-\left(i\star f_s\right)\right)\star {\left(\left(i\star f_s\right)\right)}^m\right)}^x}{{\left(\underset{i=-\Delta f/f_s}{\overset{\Delta f/f_s}{.Math.}}A\left(f_0-\left(i\star f_s\right)\right)\star {\left(\left(i\star f_s\right)\right)}^n\right)}^y}$ wherein A(f) is the amplitude for frequency f in the amplitude spectrum, f.sub.o is the frequency of interest where the oscillatory instability is expected, f.sub.s is a resolution of the amplitude spectrum and Δf is defined so that [f.sub.0−Δf, f.sub.0+Δf] is the range of interest in the amplitude spectrum where the spectral measure is evaluated and m, n, x, y are integers representing the index of the spectral measure; and detecting oscillatory instability upon the estimated amplitude exceeding a threshold amplitude; wherein the processor (108) generates a control signal based on the estimated amplitude that enables an actuator (102B) attached to the turbulent flow system to take control measures.

2. The system as claimed in claim 1, wherein the turbulent flow system (102) is at least one of gas turbine engine, wind turbine, structural elements facing turbulent flows like tall buildings, airplane wings and suspended bridges.

3. The system as claimed in claim 1, wherein a digital to analog convertor (112) is integrated with the system (100) to convert at least one digital signal obtained as an output from the amplitude estimator (110) to at least one analog signal that could be processed by a controller (114).

4. The system as claimed in 1, wherein the sensor (102A) includes a plurality of sensors that are configured to generate the at least one signal corresponding to the dynamics of the turbulent flow system (102).

5. The system as claimed in 1, wherein an analog to digital converter (106) is integrated with the system (100) to convert at least one analog signal generated by the sensor (102A) to at least one digital signal.

6. The system (100) as claimed in claim 1, wherein the system (100) further includes a signal conditioner (104) provided in communication with the sensor (102A) and configured to condition the at least one signal from the sensor (102A) in a form suitable to be transmitted to the analog to digital converter (106).

7. The system (100) as claimed in claim 1, wherein the controller (114) is configured to regulate the functioning of at least one of operating parameters of the turbulent flow system (102) by enabling the actuator (102B) to take appropriate control measures.

8. A method to determine amplitude of oscillatory instabilities in a turbulent flow system (102), the method comprising: obtaining, by an amplitude estimator (110), at least one signal corresponding to the oscillatory instabilities of the turbulent flow system (102) using a sensor (102A) connected to the amplitude estimator (110) and a processor (108) connected to the sensor (102A) and the amplitude estimator (110), wherein the sensor (102A) is mounted on the turbulent flow system (102); estimating, by the amplitude estimator (110), an amplitude of the impending oscillatory instabilities by estimating a spectral measure value of the at least one signal, wherein the spectral measures are calculated based on following equation: $M_{f_{0,\Delta f,f_s}}^{m,n,x,y}=\frac{{\left(\underset{i=-\Delta f/f_s}{\overset{\Delta f/f_s}{.Math.}}A\left(f_0-\left(i\star f_s\right)\right)\star {\left(\left(i\star f_s\right)\right)}^m\right)}^x}{{\left(\underset{i=-\Delta f/f_s}{\overset{\Delta f/f_s}{.Math.}}A\left(f_0-\left(i\star f_s\right)\right)\star {\left(\left(i\star f_s\right)\right)}^n\right)}^y}$ wherein A(f) is the amplitude for frequency f in the amplitude spectrum, f.sub.o is a frequency of interest where the oscillatory instability is expected, f.sub.s is the resolution of the amplitude spectrum and Δf is defined so that [f.sub.0−Δf, f.sub.0+Δf] is a range of interest in the amplitude spectrum where the spectral measure is evaluated and m, n, x, y are integers representing the index of the spectral measure; converting the at least one signal into a amplitude spectrum, and detecting, by the amplitude estimator (110), oscillatory instability upon the estimated amplitude exceeding a threshold amplitude; wherein the processor (108) generates a control signal based on the estimated amplitude that enables an actuator (102B) attached to the turbulent flow system (102) to take control measures.

9. The method as claimed in claim 8, wherein a digital to analog convertor (112) communicably coupled to the amplitude estimator (110) is configured to convert at least one digital signal obtained as an output from the amplitude estimator (110) to at least one analog signal that could be processed by a controller (114).

10. The method as claimed in claim 8, wherein an analog to digital converter (106) is communicably coupled to the amplitude estimator (110) to convert at least one analog signal generated by the sensor (102A) to at least one digital signal.

11. The method as claimed in claim 8, wherein a controller (114), communicably coupled to the processor (108) and the amplitude estimator (110), is configured to regulate the functioning of at least one of operating parameters of the turbulent flow system (102) by enabling the actuator (102B) to take appropriate control measures.

Description

BRIEF DESCRIPTION OF FIGURES

(1) This invention is illustrated in the accompanying drawings, throughout which like reference letters indicate corresponding parts in the various figures. The embodiments herein will be better understood from the following description with reference to the drawings, in which:

(2) FIG. 1 is a block diagram illustrating a system used to estimate the amplitude of limit cycle oscillations, according to embodiments as disclosed herein;

(3) FIG. 2 is a graphical diagram illustrating the variation of spectral measure with amplitude of oscillations in a log-log plot, according to embodiments as disclosed herein; and

(4) FIG. 3 is a graphical diagram illustrating the estimation of amplitude of oscillations for a system approaching oscillatory instability using spectral measure, according to embodiments as disclosed herein.

(5) While embodiments of the present disclosure are described herein by way of example using several illustrative drawings, those skilled in the art will recognize the present disclosure is not limited to the embodiments or drawings described. It should be understood the drawings and the detailed description thereto are not intended to limit the present disclosure to the form disclosed, but to the contrary, the present disclosure is to cover all modification, equivalents and alternatives falling within the spirit and scope of embodiments of the present disclosure as defined by the appended claims.

DETAILED DESCRIPTION OF INVENTION

(6) The embodiments herein and the various features and advantageous details thereof are explained more fully with reference to the non-limiting embodiments that are illustrated in the accompanying drawings and detailed in the following description. Descriptions of well-known components and processing techniques are omitted so as to not unnecessarily obscure the embodiments herein. The examples used herein are intended merely to facilitate an understanding of the ways in which the embodiments may be practiced and to further enable those of skill in the art to practice the embodiments. For example, although, certain embodiments herein are related to the system and method for early detection of instabilities in devices such as combustors in gas turbines, and industrial processing devices such as furnaces and burners for the ease of understanding the invention, it should be noted that the system according to the present invention may also be used for any other devices in which the transition to oscillatory instability from chaotic behavior happens though intermittent bursts.

(7) Further, although, certain embodiments herein are related to the system and method for detecting and preventing oscillatory instabilities in combustion systems, it should be noted that the system and method according to the present invention could also be used for preventing oscillatory instability in any noisy or chaotic or turbulent systems; for example, structural instabilities that may arise due to aeroelastic flutter or flow induced vibration, aerodynamic and aeromechanical instabilities such as surge and, or instabilities arising in magnetohydrodynamics, or aeroacoustic instabilities in gas transport systems exhibiting pipe tone instabilities. Accordingly, the examples should not be construed as limiting the scope of the embodiments herein.

(8) According to embodiments disclosed herein, FIG. 1 is a block diagram of an apparatus 100. The apparatus 100 comprises of a sensor 102A attached to a turbulent flow system 102. The turbulent flow system 102 can be, but not limited to, gas turbine engines, wind turbines, structural elements facing turbulent flows like tall buildings, airplane wings, suspended bridges and the like. The sensor 102A obtains a signal corresponding to an oscillatory quantity of the turbulent flow system 102. The signal from the sensor 102A is passed on to a signal conditioner 104 which in turn conditions the signal into a suitable form which can be used for the analysis. An analog to digital (A/D) converter 106 converts the analog signals to digital signals which can be processed in a processor 108. In the processor 108, the signal is sent to an amplitude estimator 110 which estimates the amplitude of the oscillatory instability in the turbulent flow system 102. The estimated amplitude is passed back into the processor 108 which in response produces a control signal. The control signal is send to a digital to analog (D/A) converter 112 which in turn acts upon a controller 114 to enable an actuator 102B attached to the turbulent flow system 102 to take appropriate control measures.

(9) Estimation of Spectral Measures:

(10) The amplitude estimator 110 works as follows. The time series, t(k) of the oscillatory quantity, obtained from the turbulent flow system 102 using the sensor 102A, is Fourier transformed into a corresponding frequency domain and the Fourier series, A(f) is obtained. The spectral measures are then calculated based on the following formula:

(11) $M_{f_{0,\Delta f,f_s}}^{m,n,x,y}=\frac{{\left(\underset{i=-\mathrm{\Delta f}/f_s}{\overset{\mathrm{\Delta f}/f_s}{.Math.}}A\left(f_0-\left(i*f_s\right)\right)*{\left(\left(i*f_s\right)\right)}^m\right)}^x}{{\left(\underset{i=-\mathrm{\Delta f}/f_s}{\overset{\mathrm{\Delta f}/f_s}{.Math.}}A\left(f_0-\left(i*f_s\right)\right)*{\left(\left(i*f_s\right)\right)}^n\right)}^y}$
Where, A(f) is the amplitude of the frequency spectrum for frequency f, f.sub.0 is the frequency of interest where the oscillatory instability is expected, f.sub.s is the resolution of the frequency spectrum and Δf is defined so that [f.sub.0−Δf, f.sub.0+Δf] is the range of interest in the frequency spectrum where the spectral measure is evaluated and m, n, x, y are integers representing the index of the spectral measure. In another embodiment, the can also be evaluated using an integral form of equation,

(12) $M_{f_{0,}\Delta f}^{m,n,x,y}=\frac{{\left({\int }_{f_0-\Delta f}^{f_0+\Delta f}A\left(f\right){\left(f_0-f\right)}^{m}df\right)}^x}{{\left({\int }_{f_0-\Delta f}^{f_0+\Delta f}A\left(f\right){\left(f_0-f\right)}^{n}df\right)}^y}$
In an embodiment, is evaluated for various fluid dynamic systems such as an aero acoustic system, and two configuration of gas turbine combustors as they approach oscillatory instabilities. versus peak amplitude of power spectra at different states of the system is illustrated in FIG. 2. It can be observed that in the log-log plot, the variation of versus peak amplitude is linear. Further, the slope of the line is same for all the cases. As we approach oscillatory instability, the amplitude of oscillation increases, and the value of the spectral measure reduces. The measure can thus be used as a precursor to an impending oscillatory instability. At the onset of oscillatory instability, the spectral measure crosses the threshold value. The threshold value can be appropriately fixed using conventional methods available in the existing state of the art.
Estimation of Amplitude of Instability:

(13) From FIG. 2, it can be seen that for all transitions, in a log-log plot, the spectral measure varies linearly with the peak amplitude of oscillations. It is shown in FIG. 2 that all systems have the same slope for the linear variation. In FIG. 3, the peak amplitude of oscillations during oscillatory instability is estimated corresponding to the aforementioned threshold value of the spectral measure by extrapolating the graph obtained in FIG. 2 towards the threshold spectral measure.

(14) Once the amplitude of the oscillations during oscillatory instability is estimated, the estimated amplitude is compared with a threshold amplitude. If the estimated amplitude exceeds the threshold amplitude, oscillatory instability is detected and the necessary control actions are taken to drive the system away from an onset of oscillatory instability.

(15) Further, it is to be understood that the phraseology or terminology employed herein is for the purpose of description and not of limitation. Therefore, while the embodiments herein have been described in terms of preferred embodiments, those skilled in the art will recognize that the embodiments herein can be practiced with modification within the spirit and scope of the embodiments as described herein.