Abstract:
The present invention relates to a sensor for measuring the instantaneous rate of mass flow and the cumulative mass flow in steady or unsteady flows of single-phase liquids or gases in a duct. By measuring the shear stress or the streamwise pressure gradient at the duct wall, and relating it to mass flow through solutions to the Navier Stokes equations of fluid mechanics, information on mass flow through the entire duct cross sectional area is deduced.

Description:
This application claims benefit of Provisional Application No. 60/099108 Sep. 3, 1998. 
    
    
     BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     The present invention relates to a sensor for measurement of the instantaneous rate of mass flow and the cumulative mass flow, in steady or unsteady flows, of single-phase liquids or gases along ducts. 
     2. Discussion of the Prior Art 
     The principal criterion for a successful design of a mass flow sensor is that it provides accurate measurements (e.g. ±0.5%) of the momentary mass-flow rate in steady and unsteady flows. In certain applications, the criterion that it measure accurately the cumulative mass flow (e.g. ±1.0%) over a specified time interval is more appropriate. 
     For example, in today&#39;s automobile engine, the device used most widely to measure the mass flow rate {dot over (m)} of air into the engine is the mass airflow sensor. The known mass airflow sensor generally comprises a “hot-wire” sensor housed in the center of a section of a duct as illustrated in FIG. 1 wherein the hot-wire sensor indicates the fluid speed |V| at the duct center. 
     When the momentary fluid speed is multiplied by the duct cross sectional area A and by the average fluid density ρ (inferred from other measurements of air temperature and pressure), a measure of the mass flow rate m=ρA|V| can be made. In this approach, the implicit assumption that the density across the flow is constant is quite reasonable in engine intake flows. A second implicit assumption, i.e., that the velocity at the duct center is also the velocity over the entire cross sectional area, is not completely valid and so A is replaced with an “effective”area A eff  determined from a steady flow calibration experiment. More significantly, the hot wire sensor infers fluid speed according to the convective heat transfer from the wire and cannot distinguish the direction of the flow. 
     Since many ducts flows, and particularly those in intakes to internal combustion engines, are not unidirectional but can momentarily reverse, this existing method may result in incorrect readings of mass flow rate. Consequently, control of the proportions of air and fuel supplied to the engine can be compromised. 
     Variations on the above described method involve mounting the hot wire sensor within a bypass tube positioned in the duct to damp out unsteadiness in part of the flow, in an effort to remove the directional ambiguity. The presence of the tube obstructs the flow and changes the flow field, necessitating calibration experiments to determine A eff , though in reality A eff  then becomes a function of time and of the unsteadiness of the particular flow. A schematic of a prior art mass airflow sensor with a bypass duct is shown in FIG.  1 . 
     In addition to the foregoing, venturi and orifice meters are known which measure mass flow rates based on an independent measurement of density and assumptions of steady flow. However, these types of meters are considered to be inadequate in non-steady flows. Turbine meters and moving vane devices also require an independent measure of density and have frequency responses which are too low for automotive intake applications, as well as being intrusive to the flow. Coriolis meters, which have high frequency responses and measure mass flow rate in direct proportion to the Coriolis force, impose large pressure drops on the flow and are most suitable for liquid flow applications. 
     As such, there is a need in the art for sensors capable of measuring instantaneous rates of mass flow and cumulative mass flow which are non-intrusive. 
     SUMMARY OF THE INVENTION 
     According to the present invention, the flow rate is determined by measuring a surface force such as the fluid shear stress τ w  or the streamwise pressure gradient ∂p/∂x with a sensor of the present invention. While the sensor is preferably enclosed within a duct, the duct itself may be of various cross-sectional shapes, i.e., round, square, eliptical, convoluted, etc. By deducing the flow rate {dot over (m)} or cumulative mass flow m cum  from a fluid-mechanics relationship between the mass-flow rate and the measured shear stress, the fluid stress and/or streamwise pressure gradient can be obtained. For fully developed, unsteady, laminar, axisymmetric flow of a Newtonian fluid, these relationships are of the form:            m   .            ∫   0   t              (       -       ∂   p       ∂   x              (     t   ′     )       )     ·     f        (     t   -     t   ′       )                   t   ′                                  m   cum     =         R   2     v            ∫   0     t   cum                (       -       ∂   p       ∂   x              (     t   ′     )       )     ·     g        (     t   -     t   ′       )                   t   ′                   or           m   .     =     R          ∫   0   t                τ   w          (     t   ′     )       ·     h        (     t   -     t   ′       )                   t   ′                       m   cum     =         R   3     v            ∫   0     t   cum                  τ   w          (     t   ′     )       ·     j        (     t   -     t   ′       )                   t   ′                                    
     where t is time, R is the duct radius, ν is the fluid kinematic viscosity, ∂p/∂x is the streamwise pressure gradient, and f g, h, and j are smoothly varying mathematical functions. Processing of shear-stress or pressure gradient measurements with either of these relationships, which involve convolution integrals, are carried out in real-time using an integrated analog-electronics circuit. In specific applications of this kind, where functions f, g, h, and j are known, analog circuits can have higher processing speeds than digital ones. While either shear stress or the pressure gradient could be used as the sensed variable in flows which are fully developed at the sensor location, it is necessary to sense shear stress when the unsteadiness in the flow past the sensor is wave-like. 
     When compared to the hot-wire mass-airflow sensor used in today&#39;s automobile engine, the present has the advantages that: i) it is wall mounted and therefore non-intrusive to the flow; ii) it can be combined easily with other sensors; iii) it can measure mass-flow rate without requiring an independent measurement of density; and iv) it measures the direction of the flow. Further, when contrasted with orifice meters, the present invention has the advantages of performing accurately in unsteady flow and being able to measure flow direction. When compared with Coriolis meters, the present invention has the advantage that it does not impose a large pressure drop on the flow. 
     One application for the sensor and methods of the present invention is that automotive engine manufacturers could use the sensor in place of existing mass-airflow sensors in engines as well as for measuring flow rates of mixtures of air, exhaust gases and hydrocarbon vapors in engine intake systems. Still another perceived application for the sensor of the present invention is the measurement of pulmonary flows and other unsteady duct flows in which a surface mounted sensor is preferred. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a side elevation view of a central section through a prior art “hot wire” sensor; 
     FIG. 1A is an end view of the prior art hot wire sensor of FIG. 1; 
     FIG. 2 is a side elevation view illustrating the sensor of the present invention; 
     FIG. 2A is an end view of the sensor of FIG. 2; 
     FIG. 3 is an elevational view illustrating the shear stress sensor of the present invention; 
     FIG. 3A is a cross sectional view of the sensor of FIG. 2 taken at line A—A; 
     FIG. 4 is an elevational view of the sensor of FIG. 2 illustrating the position of the pressure sensors in relation to the velocity and shear force vectors; 
     FIG. 5 is a schematic diagram illustrating a convolution integral circuit for single voltage inputs. 
     FIG. 6 is a schematic diagram illustrating a convolution integral circuit for differential voltage inputs; and 
     FIG. 7 is a schematic diagram of a circuit which converts the difference in currents in reverse-bias junctions into voltage. 
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     Referring to FIGS. 2 and 2A, the sensor  10  of the present invention, which preferably is formed from a silicon or polysilicon, comprises a support member  12  and a floating pad  14 . The floating pad  14  is etched out of the support member  12  utilizing known micromachining techniques to ensure that an exposed surface  16  of the sensor along both the support member  12  and the floating pad  14  are disposed along the same plane. As a result of the etching, the sensor  10  includes a lateral gap  18  occurring between the support member  12  and floating pad  14  which is on the order of 0.1 mm and a vertical gap  20 . Thus, the floating pad  14  is suspended by tethers or coupling members  22  occurring in relative proximity to the perimeter or edge  24  of the pad. 
     To implement the sensor  10  as a shear stress sensor, i.e., for flows in which the streamwise velocity varies with streamwise distance at the plane of the sensor such as, for example, along an interior surface  28  of a duct  26 , the fluid shear stress acting on the floating pad  14  is measured as a function of the difference between the extensional and compressional strain of the tethers  22 , using area ratios and the Young&#39;s modulus for the tether material to relate strain to stress. To measure strains, the tethers  22  are equipped with a semiconductor piezo-resistive strain gauge as shown in FIG. 3 or optionally with a piezo-tunneling (reverse-bias n ++ /p +  or p ++ /n +  junction) strain sensor (not shown) adhered to the top surface of or fabricated within each tether at position  24  as shown in FIG. 3. A suitable commercially available piezoelectric pressure sensor is Data Instruments Corp., Model DUXL 20D 7C21-7. 
     When there is flow past the sensor  10 , a shear force acts on the exposed surface  16  of the floating pad. This shear force is equal to the shear stress imparted by the streamwise fluid velocity when integrated over the pad&#39;s upper surface area. Thus, larger pads produce greater shear forces but reduce the spatial resolution of the sensor. In equilibrium, the shear force acting on the sensor exposed surface is balanced by equal and opposite forces of tension/compression, shared almost equally between each tether arm  22 . The forces in each tether arm  22  divided by the tether cross-sectional area determine the tensile/compressive stresses in each tether arm, with smaller cross-sectional areas resulting in larger stresses. The strains measured by strain gauges attached to the tether arms  22  are the stresses divided by the Young&#39;s modulus of the tether-arm material. The cross-sectional area of the tether arms is chosen to be small enough to produce a measurable strain, yet large enough to resist failure through overloading, vibration and other external forces. Thus, the overall ratio of measured strain to applied shear stress is proportional to the ratio of pad surface area to tether-arm cross-sectional area. Other design considerations favor using as small a pad  14  as possible, to provide the fastest response to changing shear forces. 
     Preferably, in a strain gauge implementation, two strain gauges would be employed as the active arms of a Wheatstone bridge (temperature compensated). The remainder of the bridge circuit could be located on the support member  12  or external to the entire sensor assembly. In the event that a reverse bias junction sensor is employed, the difference in the currents through each junction would be measured by circuitry (described below) positioned predominately along the support member  12 . This circuitry is described below (see attached pages). 
     The sensor  10 , when utilized as a stress-strain instrument, can be employed singly to analyze axisymmetric flows or in an array to analyze circular asymmetric flows in a duct wall by way of non-limiting example. 
     Referring to FIG. 4, the sensor  10  is shown as being employed as part of a pressure-gradient sensor system for the measurement of mass-flow rates. According to this embodiment, two sensors  11  and  11 A are positioned a pre-determined distance apart in order to measure the difference in pressure between the sensors  11  and  11 A. As should be understood by those skilled in the art, pressure-gradient sensors are useful for flows in which streamwise velocity past the sensors does not vary with streamwise distance. 
     Referring to FIG. 2, the sensor  10  is shown being employed as part of a shear-stress sensor system for measuring mass flow rates. As should be understood by the foregoing description, the fluid stress acting on the pad is measured as a function of the difference between the extensional and compressional strains on the tethers using area ratios and Young&#39;s modulus for the sensor material, i.e., silicon or polysilicon to relate stress to strain. 
     To process the sensor signal, analog circuits which perform real time convolution integrals of sensor outputs for τ ω  and ∂P/∂κ with functions f g, h and j are utilized. An example is given of the circuit which would perform the convolution integral for functions f and g, whose shapes approximate exponential decays or piecewise-exponential decays. The convolution integral circuit is designed for possible integration within the sensor microstructure. With minor variations, it can accept as its input signal either: i) a single absolute voltage between ground and the supply voltage (FIG.  5 ), as provided by reverse-bias piezo-tunneling strain sensors whose current difference is converted to a voltage; or ii) a differential voltage (FIG.  6 ), as supplied by a shear stress sensor with piezoresistive strain gages in two arms of a Wheatstone bridge, or supplied as the outputs from two pressure transducers. 
     The circuit shown in FIG. 5 carries out the convolution integral of input voltage V in  with an exponentially decaying function, providing the result of the integral as voltage V out . The circuit comprises two transconductance amplifiers and a capacitive element  30 . The capacitive element  30  may be integrated within the sensor chip if design considerations allow it to be less than a few picofarads. Otherwise, a larger capacitance would require external mounting and connections. The standard transconductance amplifier, in CMOS or bipolar circuits, comprises a differential pair and a current mirror block. Typically, this can be achieved by 3 transistors for the differential pair and 2 to 6 transistors for the mirror(s). 
     The bias voltage of transconductance amplifier  1 , applied at node  7  in FIG. 5, serves as a scaling factor for the input signal V in . The bias voltage applied at node  8  of transconductance amplifier  2  controls the convolution circuit&#39;s impulse response. 
     The model equations governing this convolution circuit are          C               V   out            t         =         I   1     +     I   2       =         G   1          (       V   ε     -     V   ref       )       -       G   2          (       V   out     -     V   ref       )                                  
     where C is the capacitance. The “nonlinear” transconductance amplifier functions G 1  ( ) and G 2  ( ) approximate hyperbolic tangents with amplitudes that depends on the bias voltages at nodes  7  and  8  respectively, when operated in the subthreshold regions of CMOS transistors. By keeping the amplitudes below about 80 millivolts, functions G 1  and G 2  are restricted to the linear regions of the hyperbolic tangent functions and the bias voltage controlling G 2  serves to change the slope of G 2 , effectively changing the exponent of the exponentially decaying impulse response. When the sensor output is a differential voltage, convolutions of the voltage difference V in 1 −V in 2  with decaying exponential functions can be achieved with a similar circuit shown in FIG.  4 B. These circuits are very efficient and consume only a fraction of a micro- or milli-watt when they operate in the subthreshold region. 
     If a wall shear stress is measured in terms of the difference between a compressive and a tensile strain, using reverse-bias piezo-tunneling strain sensors implanted on tether arms, a circuit is needed to convert the difference in small sensor currents ΔI=I 1 −I 2  to the voltage input required for the convolution integral circuit. This current conversion circuit is shown in FIG.  7  and is also designed so it may be integrated as part of the sensor structure on a single chip or module. In this design, the piezo-tunneling junctions, which act as strain sensors, are connected via a current mirror which feeds the current difference ΔI to the sense-amp block. Two cascade transistors control the performance of the circuits via the bias voltages V cas 1  and V cas 2 . In practice, one of the bias voltages (V cas 1 ) may be fixed while the second one would be adjusted to optimize circuit performance. The measured current difference ΔI is fed to the sense amp block, which includes an operational amplifier and a large resistance R to scale the current difference. The governing equation of the overall circuit is: 
     
       
           V   sense   =V   ref   −RΔI . 
       
     
     Thus, the wall shear stress is measured in terms of the sensor currents I 1 , and I 2 , whose difference ΔI is processed as an amplified and offset voltage V sense . This voltage may then be fed to the convolution integral circuit of FIG. 4A, the output voltage of which will be proportional to the momentary mass flow rate. 
     The sensors described herein can be used in conjunction with other types of sensors such as oxygen and hydrocarbon sensors for medical applications and semi-conductors or solid state sensors for the measurement of chemical composition, by way of non-limiting example. 
     As an example of the calculation of mass flow rate {dot over (m)}(t) from measurements of wall shear stress τ ω (t) in an unsteady flow, we consider a pipe flow in which the measured shear stress oscillates back and forth as a sine wave at angular velocity ω, or τ ω =sin (ωt). The {dot over (m)}˜τ ω  relation is written as: 
     
       
         
           
             
               
                 m 
                 . 
               
                
               
                 ( 
                 t 
                 ) 
               
             
             = 
             
               R 
                
               
                 
                   ∫ 
                   0 
                   t 
                 
                  
                 
                   
                     
                       τ 
                       ω 
                     
                      
                     
                       ( 
                       
                         t 
                         ′ 
                       
                       ) 
                     
                   
                    
                   
                     h 
                      
                     
                       ( 
                       
                         t 
                         - 
                         
                           t 
                           ′ 
                         
                       
                       ) 
                     
                   
                    
                   
                      
                     
                       t 
                       ′ 
                     
                   
                 
               
             
           
         
                 
         
             
         
      
     
     but can be equivalently written as            m   .          (   t   )       -     R          ∫   o   t              ∂     τ   w         ∂   t            (     t   ′     )            h   1          (     t   -     t   ′       )                 t   ′                                    
     where the Laplace transform of h 1  is equal to the Laplace transform of h divided by s. If h 1 (t−t′) is set equal to πR 2 /(4ν)+πH 1 (t−t′) then            m   .          (   t   )       =           π                   R   3         4      v              τ   ω          (   t   )         +     π                 R          ∫   0   t              ∂     τ   w         ∂   t            (     t   ′     )            H   1          (     t   -     t   ′       )                 t   ′                                      
     Finally, if G is chosen so that its Laplace transform is equal to the Laplace transform of H 1  multiplied by s            m   .          (   t   )       =           π                   R   3         4      v              τ   ω          (   t   )         +     π                 R          ∫   0   t              τ   ω          (     t   ′     )            G        (     t   -     t   ′       )                 t   ′                                      
     The first term in this expression is the steady flow result while the second term represents the extra effect of unsteadiness. For laminar fully-developed pipe flow, the function G will not have a simple analytical form, so we approximate it as a decaying exponential for the purposes of this example calculation. Thus we set G(t−t 1 )=e −(t−t     1     )ν/R     3   , where the normalization by R 2 /ν assures a proper dimensionless exponent. For this shear-stress history and function G, the {dot over (m)}˜τ ω  relation becomes:            m   .          (   t   )       =           π                   R   3         4      v              τ   ω          (   t   )         +     π                 R          ∫   0   t            sin        (     ω                   t   ′       )                   -     (     t   -     τ   ′       )            v   /     R   2                     t   ′                                      
     The integration can be performed analytically to yield the result:                  m   .          (   t   )       =                      π                   R   3         4      v              τ   ω          (   t   )         +                                  π                   R   3       v            1     1   +       (     ω                     R   2     /   v       )     2              [       sin        (     ω                 t     )       -         ω                   R   2       v          cos        (     ω                 t     )         +            -   v                     t   /     R   2             ]                                      
     This result, which includes unsteady effects, has an additional group of terms which makes it significantly different from the result which would be obtained using only steady-flow principles, which is:            m   .          (   t   )       =         π                   R   3         4      v              τ   ω          (   t   )                                
     If the exact function G, derived from solution of the Navier Stokes equations, were used, a similar result and a similar distinction between unsteady- and steady-flow analysis would be deduced. 
     While it will be apparent that the preferred embodiments of the invention disclosed are well calculated to fulfill the objects stated, it will be appreciated that the invention is susceptible to modification, variation and change without departing from the spirit thereof.