BEST MODE FOR CARRYING OUT THE INVENTION

[0030] Referring to FIG. 1 , an apparatus, generally shown as 10 , is provided to measure gas volume fraction in liquids and mixtures (e.g. paper and pulp slurries or other solid liquid mixtures) having entrained gas therein (including air). The apparatus 10 in accordance with the present invention determines the speed at which sound propagates through the fluid 12 within a pipe 14 to measure entrained gas in liquids and/or mixtures 12 . To simplify the explanation of the present invention, the flow 12 propagating through the pipe will be referred to as a mixture or slurry with the understanding that the flow may be a liquid or any other mixture having entrained gas therein.

[0031] The following approach may be used with any technique that measures the sound speed of a flow or speed at which sound propagates through the flow 12 . However, it is particularly synergistic with flow meters using sonar-based array processing, such as described in U.S. patent application Ser. No. ______ (Cidra's Docket No. CC-0122A) and U.S. patent application Ser. No. 09/729,994 (Cidra's Docket No. CC-0297), filed December 4, 200, now U.S. Pat. No. 6,609,069, which are incorporated herein by reference. While the sonar-based flow meter using an array of sensors to measure the speed of sound of an acoustic wave propagating through the mixture is shown and described, one will appreciate that any means for measuring the speed of sound of the acoustic wave may used to determine the entrained gas volume fraction of the mixture/fluid.

[0032] FIG. 2 is a block diagram 1 of the apparatus 10 of FIG. 1 that includes a device 2 for measuring the speed of sound (SOS) propagating through the flow 12 within a pipe 14 . A pressure sensor and/or temperature sensor 3 , 4 measures the pressure and/or temperature, respective, of the mixture 12 flowing through the pipe. In response to the speed of sound signal 5 and the characteristics 6 of the flow (e.g., pressure and temperature), an entrained gas processing unit 25 determines the gas volume fraction (GVF) of the flow 12 . The pressure and temperature sensors enables the apparatus 10 to compensate or determine the gas volume fraction for dynamic changes in the pressure and temperature of the flow 12 . Alternatively, the pressure and/or temperature may be estimated rather than actually measured.

[0033] A flow chart 13 shown in FIG. 3 illustrates the function of the entrained gas processing unit 25 . As shown in FIG. 2 , the inputs to the processing unit 25 includes the speed of sound (SOS) 5 within the mixture 12 in the pipe 14 , and the pressure and/or temperature of the mixture. The fluid properties of the mixture (e.g., SOS and density) are determined knowing the pressure and temperature of the mixture. The gas volume fraction of the mixture (GVF) is determined using the SOS measurement and fluid properties, which will be described in greater detail hereinafter.

[0034] Other information relating to the gas volume fraction in a fluid and the speed of sound (or sonic velocity) in the fluid, is described in “Fluid Mechanics and Measurements in two-phase flow Systems”, Institution of mechanical engineers, proceedings 1969-1970 Vol. 184 part 3C, Sept. 24-25, 1969, Birdcage Walk, Westminster, London S.W. 1, England, which is incorporated herein by reference.

[0035] FIG. 1 illustrates a schematic drawing of one embodiment of the present invention. The apparatus 10 includes a sensing device 16 comprising an array of pressure sensors (or transducers) 18 - 21 spaced axially along the outer surface 22 of a pipe 14 , having a process flow propagating therein. The pressure sensors measure the unsteady pressures produced by acoustical disturbances within the pipe, which are indicative of the SOS propagating through the mixture 12 . The output signals (P ₁ -P _N ) of the pressure sensors 18 - 21 are provided to the processor 24 , which processes the pressure measurement data and determines the speed of sound and gas volume fraction (GVF).

[0036] In an embodiment of the present invention shown in FIG. 1 , the apparatus 10 has at least pressure sensors 18 - 21 disposed axially along the pipe 14 for measuring the unsteady pressure P ₁ -P _N of the mixture 12 flowing therethrough. The speed of sound propagating through the flow 12 is derived by interpreting the unsteady pressure field within the process piping 14 using multiple transducers displaced axially over 2 diameters in length. The flow measurements can be performed using ported pressure transducers or clamp-on, strain-based sensors.

[0037] The apparatus 10 has the ability to measure the gas volume fraction by determining the speed of sound of acoustical disturbances or sound waves propagating through the flow 12 using the array of pressure sensors 18 - 21 . While the apparatus of FIG. 1 shows at least four pressure sensors 18 - 21 , the present invention contemplates an apparatus having an array of two or more pressure sensors and having as many as sixteen (16) pressure sensors.

[0038] Generally, the apparatus 10 measures unsteady pressures created by acoustical disturbances propagating through the flow 12 to determine the speed of sound (SOS) propagating through the flow. Knowing the pressure and/or temperature of the flow and the speed of sound of the acoustical disturbances, the processing unit 24 can determine the gas volume fraction of the mixture, as described and shown in FIG. 3 .

[0039] The apparatus in FIG. 1 also contemplates providing one or more acoustic sources 27 to enable the measurement of the speed of sound propagating through the flow for instances of acoustically quiet flow. The acoustic source may be a device the taps or vibrates on the wall of the pipe, for example. The acoustic sources may be disposed at the input end of output end of the array of sensors 18 - 21 , or at both ends as shown. One should appreciate that in most instances the acoustics sources are not necessary and the apparatus passively detects the acoustic ridge provided in the flow 12 . The passive noise includes noise generated by pumps, valves, motors, and the turbulent mixture itself.

[0040] The apparatus 10 of the present invention may be configured and programmed to measure and process the detected unsteady pressures P ₁ (t)-P _N (t) created by acoustic waves propagating through the mixture to determine the SOS through the flow 12 in the pipe 14 . One such apparatus 110 is shown in FIG. 4 that measures the speed of sound (SOS) of one-dimensional sound waves propagating through the mixture to determine the gas volume fraction of the mixture. It is known that sound propagates through various mediums at various speeds in such fields as SONAR and RADAR fields. The speed of sound propagating through the pipe and mixture 12 may be determined using a number of known techniques, such as those set forth in U.S. patent application Ser. No. 09/344,094, entitled “Fluid Parameter Measurement in Pipes Using Acoustic Pressures”, filed Jun. 25, 1999, now U.S. Pat. No. 6,354,147; U.S. patent application Ser. No. 09/729,994, filed Dec. 4, 2002, now U.S. Pat. No. 6,609,069; U.S. patent application Ser. No. 09/997,221, filed Nov. 28, 2001, now U.S. Pat. No. 6,587,798; and U.S. patent application Ser. No. 10/007,749, entitled “Fluid Parameter Measurement in Pipes Using Acoustic Pressures”, filed Nov. 7, 2001, each of which are incorporated herein by reference.

[0041] In accordance with one embodiment of the present invention, the speed of sound propagating through the mixture 12 is measured by passively listening to the flow with an array of unsteady pressure sensors to determine the speed at which one-dimensional compression waves propagate through the mixture 12 contained within the pipe 14 .

[0042] As shown in FIG. 4 , an apparatus 110 embodying the present invention has an array of at least three acoustic pressure sensors 115 , 116 , 117 , located at three locations x ₁ ,x ₂ ,x ₃ axially along the pipe 14 . One will appreciate that the sensor array may include more than three pressure sensors as depicted by pressure sensor 118 at location x _N . The pressure generated by the acoustic waves may be measured through pressure sensors 115 - 118 . The pressure sensors 15 - 18 provide pressure time-varying signals P ₁ (t),P ₂ (t),P ₃ (t),P _N (t) on lines 120 , 121 , 122 , 123 to a signal processing unit 130 to known Fast Fourier Transform (FFT) logics 126 , 127 , 128 , 129 , respectively. The FFT logics 126 - 129 calculate the Fourier transform of the time-based input signals P ₁ (t)-P _N (t) and provide complex frequency domain (or frequency based) signals P ₁ (ω),P ₂ (ω),P ₃ (ω),P _N (ω) on lines 132 , 133 , 134 , 135 indicative of the frequency content of the input signals. Instead of FFT's, any other technique for obtaining the frequency domain characteristics of the signals P ₁ (t)-P _N (t), may be used. For example, the cross-spectral density and the power spectral density may be used to form a frequency domain transfer functions (or frequency response or ratios) discussed hereinafter.

[0043] The frequency signals P ₁ (ω)-P _N (ω) are fed to an array processing unit 138 which provides a signal to line 40 indicative of the speed of sound of the mixture a _mix , discussed more hereinafter. The a _mix signal is provided to an entrained gas processing unit 142 , similar to the processing unit 25 , which converts a _mix to a percent composition of a mixture and provides a gas volume fraction or % Comp signal to line 44 indicative thereof (as discussed hereinafter).

[0044] The data from the array of sensors 115 - 118 may be processed in any domain, including the frequency/spatial domain, the temporal/spatial domain, the temporal/wave-number domain or the wave-number/frequency (k-ω) domain. As such, any known array processing technique in any of these or other related domains may be used if desired, similar to the techniques used in the fields of SONAR and RADAR.

[0045] One such technique of determining the speed of sound propagating through the flow 12 is using array processing techniques to define an acoustic ridge in the k-ω plane as shown in FIG. 9 . The slope of the acoustic ridge is indicative of the speed of sound propagating through the flow 12 . This technique is similar to that described in U.S. Pat. No. 6,587,798 filed Nov. 28, 2001, titled “Method and System for Determining The Speed of Sound in a Fluid Within a Conduit”, which is incorporated herein by reference. The speed of sound (SOS) is determined by applying sonar arraying processing techniques to determine the speed at which the one dimensional acoustic waves propagate past the axial array of unsteady pressure measurements distributed along the pipe 14 .

[0046] The signal processor 24 performs a Fast Fourier Transform (FFT) of the time-based pressure signals P ₁ (t)-P _N (t) to convert the pressure signal into the frequency domain. The power of the frequency-domain pressure signals are then determined and defined in the k-ω plane by using array processing algorithms (such as Capon and Music algorithms). The acoustic ridge in the k-ω plane, as shown in the k-ω plot of FIG. 9 , is then determined. The speed of sound (SOS) is determined by measuring slope of the acoustic ridge. The gas volume fraction is then calculated or otherwise determined, as described hereinafter.

[0047] The flow meter of the present invention uses known array processing techniques, in particular the Minimum Variance, Distortionless Response (MVDR, or Capon technique), to identify pressure fluctuations, which convect with the materials flowing in a conduit and accurately ascertain the velocity, and thus the flow rate, of said material. These processing techniques utilize the covariance between multiple sensors 18 - 21 at a plurality of frequencies to identify signals that behave according to a given assumed model; in the case of the apparatus 10 , a model, which represents pressure variations 20 convecting at a constant speed across the pressure sensors comprising the flow meter monitoring head 12 .

[0048] To calculate the power in the k-ω plane, as represent by a k-ω plot (see FIG. 9 ) of either the pressure signals, the processor 58 determines the wavelength and so the (spatial) wavenumber k, and also the (temporal) frequency and so the angular frequency ω, of various spectral components of the acoustic waves created passively or actively within the pipe. There are numerous algorithms available in the public domain to perform the spatial/temporal decomposition of arrays of sensor units 18 - 21 .

[0049] In the case of suitable acoustic pressures being present, the power in the k-ω plane shown in a k-ω plot of FIG. 9 so determined will exhibit a structure that is called an acoustic ridge 61 associated with sound propagating with the flow and one associated with sound propagating against the flow. The acoustic ridge represents the concentration of the disturbances that propagate with and against the flow and is a mathematical manifestation of the relationship between the spatial variations and temporal variations described above. Such a plot will indicate a tendency for k-ω pairs to appear more or less along a line with some slope, the slope indicating the speed of sound traveling in both directions, as is described in more detail below. The power in the k-ω plane so determined is then provided to a acoustic ridge identifier, which uses one or another feature extraction method to determine the location and orientation (slope) of any acoustic ridge present in the k-ω plane. Finally, information including the acoustic ridge orientation (slope) is used by an analyzer to determine the speed of sound.

[0050] The array processing unit 23 uses standard so-called beam forming, array processing, or adaptive array-processing algorithms, i.e. algorithms for processing the sensor signals using various delays and weighting to create suitable phase relationships between the signals provided by the different sensors, thereby creating phased antenna array functionality. In other words, the beam forming or array processing algorithms transform the time domain signals from the sensor array into their spatial and temporal frequency components, i.e. into a set of wave numbers given by k=2πr/λ where λ is the wavelength of a spectral component, and corresponding angular frequencies given by ω=2πν.

[0051] The prior art teaches many algorithms of use in spatially and temporally decomposing a signal from a phased array of sensors, and the present invention is not restricted to any particular algorithm. One particular adaptive array processing algorithm is the Capon method/algorithm. While the Capon method is described as one method, the present invention contemplates the use of other adaptive array processing algorithms, such as MUSIC algorithm. The present invention recognizes that such techniques can be used to determine speed of sound propagating through the fluid 12 .

[0052] Also, some or all of the functions within the processor 130 may be implemented in software (using a microprocessor or computer) and/or firmware, or may be implemented using analog and/or digital hardware, having sufficient memory, interfaces, and capacity to perform the functions described herein.

[0053] It is within the scope of the present invention that the pressure sensor spacing may be known or arbitrary and that as few as two sensors are required if certain information is known about the acoustic properties of the process flow 12 . The pressure sensors are spaced sufficiently such that the entire length of the array (aperture) is at least a significant fraction of the measured wavelength of the acoustic waves being measured. The acoustic wavelength is a function of the type or characteristics of flow 12 .

[0054] Based on the above discussion, one may use a short length scale aperture to measure the sound speed. For example, the characteristic acoustic length scale is: λ=c/f; where c is the speed of sound in a mixture, f is frequency and λ is wavelength.

[0055] If Aperture=L and if L/λ is approx. constant.

[0056] Then Lwater/λwater=Lwater*f/C _water ≈L _GVF *f/c _GVF

[0057] Therefore: L _GVF =Lwater (C _GVF /C _water ); where GVF is gas volume fraction.

[0058] Thus for SOS of water (Cwater=5,000 ft/sec), and SOS of the Gas volume fraction (C GVF=500 ft/sec) and a length aperture of L water=5 ft (which we have shown is sufficient to accurately measure the SOS of water), the length aperture for a gas volume fraction L _GVF would be about 0.5 feet.

[0059] The entrained gas processing unit 25 assumes a nearly isothermal condition for the flow 12 . As such the gas volume fraction or the void fraction is related to the speed of sound by the following quadratic equation:

Ax ² +Bx+C= 0

[0060] wherein x is the speed of sound, A=1+rg/rl*(K _eff /P−1)−K _eff /P, B=K _eff /P−2+rg/rl; C=1−K _eff /rl*a _meas {circumflex over ( )}2); Rg=gas density, rl=liquid density, K _eff =effective K (modulus of the liquid and pipewall), P=pressure, and a _meas =measured speed of sound.

[0061] Effectively,

[0062] Gas Voulume Fraction (GVF)=(−B+sqrt(B{circumflex over ( )}2−4*A*C))/(2*A)

[0063] Alternatively, the sound speed of a mixture can be related to volumetric phase fraction (φ _i ) of the components and the sound speed (a) and densities (ρ) of the component through the Wood equation. 1 $\frac{1}{{\rho }_{\mathrm{mix}}a_{{\mathrm{mix}}_{\infty }}^2}=\sum _{i=1}^N\frac{{\phi }_i}{{\rho }_ia_i^2}\mathrm{where}{\rho }_{\mathrm{mix}}=\sum _{i=1}^N{\rho }_i{\phi }_i$

[0064] One dimensional compression waves propagating within a mixture 12 contained within a pipe 14 exert an unsteady internal pressure loading on the pipe. The degree to which the pipe displaces as a result of the unsteady pressure loading influences the speed of propagation of the compression wave. The relationship among the infinite domain speed of sound and density of a mixture; the elastic modulus (E), thickness (t), and radius (R) of a vacuum-backed cylindrical conduit; and the effective propagation velocity (a _eff ) for one dimensional compression is given by the following expression: 2 $\begin{array}{cc}{a}_{\mathrm{eff}}=\frac{1}{\sqrt{\frac{1}{a_{{\mathrm{mix}}_{\infty }}^2}+{\rho }_{\mathrm{mix}}\frac{2R}{\mathrm{Et}}}}& \left(\mathrm{eq}1\right)\end{array}$

[0065] Note: “vacuum backed” as used herein refers to a situation in which the fluid surrounding the pipe externally has negligible acoustic impedance compared to that of the mixture internal to the pipe 14 . For example, meter containing a typical water and pulp slurry immersed in air at standard atmospheric conditions satisfies this condition and can be considered “vacuum-backed”.

[0066] The mixing rule essentially states that the compressibility of a mixture (1/(ρa ² )) is the volumetrically-weighted average of the compressibilities of the components. For gas/liquid mixtures 12 at pressure and temperatures typical of paper and pulp industry, the compressibility of gas phase is orders of magnitudes greater than that of the liquid. Thus, the compressibility of the gas phase and the density of the liquid phase primarily determine mixture sound speed, and as such, it is necessary to have a good estimate of process pressure to interpret mixture sound speed in terms of volumetric fraction of entrained gas. The effect of process pressure on the relationship between sound speed and entrained air volume fraction is shown in FIG. 5 .

[0067] Conversely, however, detailed knowledge of the liquid/slurry is not required for entrained air measurement. Variations in liquid density and compressibility with changes in consistency have a negligible effect on mixture sound speed compared to the presence of entrained air. FIG. 6 shows the mixture sound speed as a function of entrained air volume fraction for two slurries, one with 0% wood fiber and the other with 5% wood fiber by volume. As shown, the relationship between mixture sound speed and gas volume fraction is essentially indistinguishable for the two slurries. Furthermore, mixture sound speed is shown to an excellent indicator of gas volume fraction, especially for the trace to moderate amounts of entrained air, from 0 to 5% by volume, typically encountered in the paper and pulp industry.

[0068] For paper and pulp slurries, the conditions are such that for slurries with non-negligible amounts of entrained gas, say <0.01%, the compliance of standard industrial piping (Schedule 10 or 40 steel pipe) is typically negligible compared to that of the entrained air.

[0069] FIGS. 7 and 8 above show the relationship between sound speed and entrained air for slurries 12 with pulp contents representative of the range used in the paper and pulp industry. Referring to FIG. 7 , two slurry consistencies are shown; representing the lower limit, a pure water mixture is considered, and representing the higher end of consistencies, a 5% pulp/95% water slurry is considered. Since the effect of entrained air on the sound speed of the mixture is highly sensitive to the compressibility of the entrained air, the effect of the entrained air is examined at two pressures, one at ambient representing the lower limit of pressure, and one at 4 atmospheres representing a typical line pressure in a paper process. As shown, the consistency of the liquid slurry 12 , i.e., the pulp content, has little effect on the relationship between entrained air volume fraction and mixture sound speed. This indicates that an entrained air measurement could be accurately performed, within 0.01% or so, with little or no knowledge of the consistency of the slurry. The chart does show a strong dependence on line pressure. Physically, this effect is linked to the compressibility of the air, and thus, this indicates that reasonable estimates of line pressure and temperature would be required to accurately interpret mixture sound speed in terms of entrained air gas volume fraction.

[0070] FIG. 7 also shows that for the region of interest, from roughly 1% entrained air to roughly 5% entrained air, mixture sound speeds (a _mix ) are quite low compare to the liquid-only sound speeds. In this example, the sound speed of the pure water and the 5% pulp slurry were calculated, based on reasonable estimates of the constituent densities and compressibilities, to be 1524 m/s and 1541 m/s, respectively. The sound speed of these mixtures with 1% to 5% entrained air at typical operating pressure (1 atm to 4 atms) are on the order of 100 m/sec. The implication of these low sound speed is that the mixture sound speed could be accurately determined with a array of sensors, ie using the methodology described in aforementioned U.S. patent application Ser. No. 10/007,749 (Cidra's Docket No. CC-0066B), with an aperture that is similar, or identical, to an array of sensors that would be suitable to determine the convection velocity, using the methodology described in aforementioned U.S. patent application Ser. No. ______ (Cidra's Docket No. CC-0122A), which is incorporated herein by reference.

[0071] For the sound speed measurement, the apparatus 1 - 10 utilizes array processing algorithms. The temporal and spatial frequency content of sound propagating within the process piping is related through a dispersion relationship. 3 $k=\frac{\omega }{a_{\mathrm{mix}}}$

[0072] As before, k is the wave number, defined as k=2π/λ, ω is the temporal frequency in rad/sec, and a _mix is the speed at which sound propagates within the process piping. Unlike disturbances, which convect with the flow, however, sound generally propagates in both directions, with and against the mean flow. For these cases, the acoustic power is located along two acoustic ridges, one for the sound traveling with the flow at a speed of a _mix +V _mix and one for the sound traveling against the flow at a speed of a _mix −V _mix . FIG. 9 shows a k-ω plot generated for acoustic sound field recorded from water flowing at a rate of 240 gpm containing ˜2% entrained air by volume in a 3 in, schedule 10 , stainless steel pipe. The k-ω plot was constructed using data from an array of strain-based sensors attached to the outside of the pipe. Two acoustic ridges are clearly evident. Based on the slopes of the acoustic ridges, the sound speed for this for this mixture was 330 ft/sec (100 m/s), consistent with that predicted by the Wood equation. Note that adding 2% air by volume reduces the sound speed of the bubbly mixture to less than 10% of the sound speed of single phase water.

[0073] In one embodiment of the present invention as shown in FIG. 1 , each of the pressure sensors 18 - 21 may include a piezoelectric film sensor to measure the unsteady pressures of the mixture 12 using either technique described hereinbefore.

[0074] The piezoelectric film sensors include a piezoelectric material or film to generate an electrical signal proportional to the degree that the material is mechanically deformed or stressed. The piezoelectric sensing element is typically conformed to allow complete or nearly complete circumferential measurement of induced strain to provide a circumferential-averaged pressure signal. The sensors can be formed from PVDF films, co-polymer films, or flexible PZT sensors, similar to that described in “Piezo Film Sensors Technical Manual” provided by Measurement Specialties, Inc., which is incorporated herein by reference. A piezoelectric film sensor that may be used for the present invention is part number 1-1002405-0, LDT4-028K, manufactured by Measurement Specialties, Inc.

[0075] Piezoelectric film (“piezofilm”), like piezoelectric material, is a dynamic material that develops an electrical charge proportional to a change in mechanical stress. Consequently, the piezoelectric material measures the strain induced within the pipe 14 due to unsteady pressure variations (e.g., vortical and/or acoustical) within the process mixture 12 . Strain within the pipe is transduced to an output voltage or current by the attached piezoelectric sensor. The piezoelectrical material or film may be formed of a polymer, such as polarized fluoropolymer, polyvinylidene fluoride (PVDF). The piezoelectric film sensors are similar to that described in U.S. patent application Ser. No. ______ (CiDRA Docket No. CC-0676), which is incorporated herein by reference.

[0076] Another embodiment of the present invention include a pressure sensor such as pipe strain sensors, accelerometers, velocity sensors or displacement sensors, discussed hereinafter, that are mounted onto a strap to enable the pressure sensor to be clamped onto the pipe. The sensors may be removable or permanently attached via known mechanical techniques such as mechanical fastener, spring loaded, clamped, clam shell arrangement, strapping or other equivalents. These certain types of pressure sensors, it may be desirable for the pipe 12 to exhibit a certain amount of pipe compliance.

[0077] Instead of single point pressure sensors 18 - 21 , at the axial locations along the pipe 12 , two or more pressure sensors may be used around the circumference of the pipe 12 at each of the axial locations. The signals from the pressure sensors around the circumference at a given axial location may be averaged to provide a cross-sectional (or circumference) averaged unsteady acoustic pressure measurement. Other numbers of acoustic pressure sensors and annular spacing may be used. Averaging multiple annular pressure sensors reduces noises from disturbances and pipe vibrations and other sources of noise not related to the one-dimensional acoustic pressure waves in the pipe 12 , thereby creating a spatial array of pressure sensors to help characterize the one-dimensional sound field within the pipe 12 .

[0078] The pressure sensors 18 - 21 of FIG. 1 described herein may be any type of pressure sensor, capable of measuring the unsteady (or ac or dynamic) pressures within a pipe 14 , such as piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers (or geophones), velocity measuring devices, displacement measuring devices, etc. If optical pressure sensors are used, the sensors 18 - 21 may be Bragg grating based pressure sensors, such as that described in U.S. patent application Ser. No. 08/925,598, entitled “High Sensitivity Fiber Optic Pressure Sensor For Use In Harsh Environments”, filed Sep. 8, 1997, now U.S. Pat. No. 6,016,702, and in U.S. patent application Ser. No. 10/224,821, entitled “Non-Intrusive Fiber Optic Pressure Sensor for Measuring Unsteady Pressures within a Pipe”, which are incorporated herein by reference. In an embodiment of the present invention that utilizes fiber optics as the pressure sensors 14 they may be connected individually or may be multiplexed along one or more optical fibers using wavelength division multiplexing (WDM), time division multiplexing (TDM), or any other optical multiplexing techniques.

[0079] In certain embodiments of the present invention, a piezo-electronic pressure transducer may be used as one or more of the pressure sensors 15 - 18 and it may measure the unsteady (or dynamic or ac) pressure variations inside the tube 14 by measuring the pressure levels inside of the tube. These sensors may be ported within the pipe to make direct contact with the mixture 12 . In an embodiment of the present invention, the sensors 14 comprise pressure sensors manufactured by PCB Piezotronics. In one pressure sensor there are integrated circuit piezoelectric voltage mode-type sensors that feature built-in microelectronic amplifiers, and convert the high-impedance charge into a low-impedance voltage output. Specifically, a Model 106B manufactured by PCB Piezotronics is used which is a high sensitivity, acceleration compensated integrated circuit piezoelectric quartz pressure sensor suitable for measuring low pressure acoustic phenomena in hydraulic and pneumatic systems. It has the unique capability to measure small pressure changes of less than 0.001 psi under high static conditions. The 106B has a 300 mV/psi sensitivity and a resolution of 91 dB (0.0001 psi).

[0080] The pressure sensors incorporate a built-in MOSFET microelectronic amplifier to convert the high-impedance charge output into a low-impedance voltage signal. The sensor is powered from a constant-current source and can operate over long coaxial or ribbon cable without signal degradation. The low-impedance voltage signal is not affected by triboelectric cable noise or insulation resistance-degrading contaminants. Power to operate integrated circuit piezoelectric sensors generally takes the form of a low-cost, 24 to 27 VDC, 2 to 20 mA constant-current supply. A data acquisition system of the present invention may incorporate constant-current power for directly powering integrated circuit piezoelectric sensors.

[0081] Most piezoelectric pressure sensors are constructed with either compression mode quartz crystals preloaded in a rigid housing, or unconstrained tourmaline crystals. These designs give the sensors microsecond response times and resonant frequencies in the hundreds of kHz, with minimal overshoot or ringing. Small diaphragm diameters ensure spatial resolution of narrow shock waves.

[0082] The output characteristic of piezoelectric pressure sensor systems is that of an AC-coupled system, where repetitive signals decay until there is an equal area above and below the original base line. As magnitude levels of the monitored event fluctuate, the output remains stabilized around the base line with the positive and negative areas of the curve remaining equal.

[0083] It is also within the scope of the present invention that any strain sensing technique may be used to measure the variations in strain in the pipe, such as highly sensitive piezoelectric, electronic or electric, strain gages and piezo-resistive strain gages attached to the pipe 12 . Other strain gages include resistive foil type gages having a race track configuration similar to that disclosed U.S. patent application Ser. No. 09/344,094, filed Jun. 25, 1999, now U.S. Pat. No. 6,354,147, which is incorporated herein by reference. The invention also contemplates strain gages being disposed about a predetermined portion of the circumference of pipe 12 . The axial placement of and separation distance ΔX ₁ , ΔX ₂ between the strain sensors are determined as described herein above.

[0084] It is also within the scope of the present invention that any other strain sensing technique may be used to measure the variations in strain in the tube, such as highly sensitive piezoelectric, electronic or electric, strain gages attached to or embedded in the tube 14 .

[0085] While a number of sensor have been described, one will appreciate that any sensor the measures the speed of sound propagating through the fluid may be used with the present invention, including ultrasonic sensors.

[0086] It should be understood that any of the features, characteristics, alternatives or modifications described regarding a particular embodiment herein may also be applied, used, or incorporated with any other embodiment described herein.

[0087] Although the invention has been described and illustrated with respect to exemplary embodiments thereof, the foregoing and various other additions and omissions may be made therein and thereto without departing from the spirit and scope of the present invention.