# Patent

## Title
Electroactive polymer-based flow sensor and methods related thereto

## Abstract

An apparatus includes an electroactive polymer member and a first electric terminal and a second electric terminal. The electroactive polymer member has a free portion configured to change position relative to a fixed portion in response to an external stimulus corresponding to at least one of a flow parameter or a fluid parameter. The first electric terminal and the second electric terminal are coupled to the electroactive polymer member. The first electric terminal and the second electric terminal provide an electrical signal in response to the change.

## Background
This application is a U.S. National Stage Filing under 35 U.S.C. 371 from International Application No. PCT/US2011/034374, filed Apr. 28, 2011, and published in English as WO 2011/137259A2 on Nov. 3, 2011, which claims the benefit under 35 U.S.C. 119 (e) of U.S. Provisional Application Ser. No. 61/330,149 filed on Apr. 30, 2010, which applications and publications are hereby incorporated by reference in their entireties.


STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT
This invention was made with government support under Grant No. 0547131 from the National Science Foundation and Grant No. N000140810640 from the Office of Naval Research. The Government has certain rights in this invention.


BACKGROUND
Accurate measurements of air and fuel flows of automobile engines are critical because they govern the fuel-air mixing processes in the combustion chamber that could directly impact the combustion efficiency and exhaust emissions.
A traditional flow sensor may be unsatisfactory for certain applications. For example, a calibration error can lead to inaccurate data.
SUMMARY
An apparatus includes at least one electroactive polymer (EAP) member (e.g., an ionic polymer-metal composite (IPMC) beam), with each of EPA members having a first electric terminal and a second electric terminal Each of the at least one EAP members further has a free portion configured to change position relative to a fixed portion in response to an external stimulus corresponding to at least one of a flow parameter or a fluid parameter. The first electric terminal and the second electric terminal are coupled to the EAP member. The first electric terminal and the second electric terminal provide an electrical signal in response to the change.



BRIEF DESCRIPTION OF THE DRAWINGS
In the drawings, which are not necessarily drawn to scale, like numerals may describe similar components in different views. Like numerals having different letter suffixes may represent different instances of similar components. The drawings illustrate generally, by way of example, but not by way of limitation, various embodiments discussed in the present document.
FIG. 1 illustrates a schematic of a beam sensor according to an example.
FIG. 2 illustrates a beam sensor in a flow channel according to an example.
FIG. 3 illustrates a beam sensor according to an example.
FIG. 4 illustrates a beam sensor according to an example.
FIG. 5 illustrates a beam sensor system according to an example.
FIG. 6 illustrates a model of a sensor system according to an example.
FIGS. 7A and 7B illustrate a sensing principle according to an example.
FIG. 8 illustrates a flow chart of a method according to an example.
FIG. 9A illustrates a beam according to an example.
FIG. 9B illustrates the beam of FIG. 9A modeled by five rigid elements according to an example.
FIG. 10 illustrates a free body diagram of the ith beam element of the beam of FIGS. 9A and 9B according to an example.
FIG. 11 is a process flow diagram showing an algorithm for converting a polymer signal into flow information according to an example.
FIG. 12 shows beam sensor responses after a flow pulse start according to an example.
FIG. 13 shows beam sensor responses after a flow pulse end according to an example.
FIG. 14 shows data selection for water at injection start according to an example.
FIG. 15 shows a tracking error of the identified result of water at injection start according to an example.
FIG. 16A shows current signal of water at injection start according to an example.
FIG. 16B shows current signal of n-heptane at injection start according to an example.
FIG. 16C shows current signal of n-heptane at injection end according to an example.



DETAILED DESCRIPTION
An example of the present subject matter can be configured for use as a sensor of flow characteristics and/or fluid characteristics. As used herein, the term “fluid” is intended to refer to a substance which can continually deform or flow under an applied shear stress. A fluid can therefore include a material in a liquid or gaseous state.
Flow characteristics (or properties) can include flow rate, pressure, duty cycle, or duration. Fluid characteristics (or properties) can include density, viscosity (liquid state), or temperature. As a flow sensor, the present subject matter can be configured to monitor flow of a fluid medium within a flow channel, such as a lumen, or configured to monitor flow relative to another structure, such as air flow over a surface of an aircraft or a rotating body (turbine).
The sensor can include a member fabricated of an electroactive polymer (EAP). For example, a conjugated polymer, dielectric elastomer, or an ionic polymer-metal composite (IPMC) can be used as part of the sensing element. In one example, the sensor includes a piezoelectric polymer or other material, such as a polyvinylidene fluoride (PVDF).
Two or more electrical terminals (electrodes) can be coupled to the sensor to provide an electrical signal corresponding to a detected physical change in the member. The electrodes can be affixed to a common side of the member or to different (or opposing) sides of the member. For example, two electrodes can be attached on the two sides of an IPMC member, respectively. In one example, each surface of the IPMC member is conductive and an electrode layer of the IPMC member can be patterned to provide multiple outputs from the member.
The member can be configured as a beam or other form. For example, movement of a free portion of the member can be measured relative to a fixed portion of the member. The fixed portion of the member is attached to a substrate or other relatively stationary structure.
In various examples, the member includes a single cantilevered beam having one portion fixed (or securely clamped) or a double cantilevered beam having opposing ends fixed and a center portion allowed to move freely. For instance, the member can include a flag-shaped portion or a round center portion. The perimeter of a round member can be fixed and the central portion can be free to move. In one example, the central portion is fixed and the perimeter is free to move. More complex-shaped members can also be used including plates or membranes.
In one example, an algorithm can be executed to “invert” the sensor dynamics to infer information on deformation. The algorithm can include computation elements to allow modeling of the flow sensor in order to compensate for temperature-dependent behavior of the sensor.
In one embodiment, an electroactive polymer (EAP) member in the form of a cantilever beam, such as an ionic polymer-metal composite (IPMC) cantilever beam (hereinafter “IPMC beam”) is used for sensing characteristics of the surrounding fluid medium, such as temperature. The temperature-dependent dynamics of an IPMC beam can be experimentally characterized and modeled. The IPMC beam can be submersed in a surrounding fluid medium comprising a water bath, the temperature of which can be regulated. The empirical frequency response of the sensor, with the tip displacement as input, and the short-circuit sensing current as output, shows dependence on temperature of the surrounding fluid medium. The sensing dynamics can be modeled with a transfer function with temperature-dependent coefficients. By fitting the values of the coefficients to a set of test temperatures, the temperature-dependence of coefficients of a polynomial function can be determined The model can be used in predicting the sensing behavior at various temperatures.
As described herein, in one embodiment, an IPMC beam can be used for sensing with applied mechanical stimulus causing redistribution of the ions inside an IPMC beam, resulting in a detectable electrical signal and thus the sensing effect (See FIGS. 7A and 7B). An example of the present subject matter considers the dependence of the IPMC beam sensing dynamics on the temperature. Accounting for such dependence, however, can be used to improve the accuracy of interpretation of the IPMC beam sensor output in various applications.
Testing for characterizing IPMC beam sensor dynamics at different temperatures was completed in which a (cantilevered) IPMC beam, with length 18 mm, width 6.5 mm, and thickness 0.21 mm, was soaked in a water bath sitting on a digital hotplate (Thermo Scientific, HP131225), which was feedback-controlled using a T-type thermal couple (Omega, HTTC36-T-116G-6) with an accuracy of 0.5° C. The IPMC sample was fabricated in-house, using Nafion®-117 (DuPont) as the raw material. For mechanical excitation, the tip of the IPMC beam was inserted into a slit cut in a wooden piece about 10 cm by 1 cm in size, which was mechanically coupled to the center of the membrane of a subwoofer speaker (Polk/MOMO, MM2154). The speaker was driven by a power amplifier (Biamp Systems, MPA-250). A voltage applied to the speaker induced vibration of the membrane and thus the prescribed tip displacement of the IPMC beam at the order of a few millimeters. A laser sensor (Baumer Electric, OADM 2016441/S14F) was used to measure the actual tip displacement. The short-circuit current measured with a current-amplifier circuit consisting of two cascaded operational amplifiers between the two electrodes of the IPMC beam was at the order of 20-100 micro amps. This measurement was then used as the sensor output.
Completed testing revealed various measured time responses at different bath temperatures, such as about 25° C., about 45° C., and about 65° C., although other temperatures, including all ranges there between can be utilized in the various embodiments. Both the gain and the phase shift can vary significantly with temperature. In one example, one sensor exhibited an increased sensor gain with the temperature within a frequency range of about 40 to about 60 Hz, with the trend reversing outside this range. The phase lag was shown to decrease with an increasing temperature. A physical explanation of the observed trend cannot be readily provided since the sensing dynamics depend on a number of physical factors (such as diffusivity, resistivity, and dielectric constant). Additionally, and while not wishing to be bound by this proposed theory, it is possible that temperature of the water baths may have influenced each of these factors in a complex way.
In one embodiment an empirical model comprising a fourth-order transfer function G(s) is adopted to approximate the measured time response for IPMC beam sensors, according to the following mathematical formula:)




  
    G
    ⁡
    
      (
      s
      )
    
  
  =
  
    
      
        
          b
          4
        
        ⁢
        
          s
          4
        
      
      +
      
        
          b
          3
        
        ⁢
        
          s
          3
        
      
      +
      
        
          b
          2
        
        ⁢
        
          s
          2
        
      
      +
      
        
          b
          1
        
        ⁢
        s
      
      +
      
        b
        0
      
    
    
      
        s
        4
      
      +
      
        
          a
          3
        
        ⁢
        
          s
          3
        
      
      +
      
        
          a
          2
        
        ⁢
        
          s
          2
        
      
      +
      
        
          a
          1
        
        ⁢
        s
      
      +
      
        a
        0
      
    
  




with coefficients {ai}i=03 and {bi}i=04 dependent on the surrounding fluid medium temperature. It is desirable to capture such dependence with simple functions, so that the model at any temperature (within a selected range) can be obtained. To achieve this result, the coefficients corresponding to a set of temperatures are identified. Thereafter, through curve-fitting, low-order polynomial functions of the temperature “T” for individual coefficients can be obtained. The effectiveness of the modeling approach can be established by comparing the measured and predicted frequency responses at temperatures (e.g., 35° C. and 45° C.) that were not used in the data-fitting process. The model can also be validated in free, damped oscillations. Additional details of this multi-link dynamic modeling can be found in FIGS. 9A and 9B discussed below and in D. Hung, X. Yang, G. Zhu, T. Ganely, X. and Tan, “Experimental Studies on Using Ionic Polymer-Metal Composite Materials for Automotive Flow Sensing Applications,” published in the proceedings of the 5th IFAC Symposium on Mechatronic Systems, Cambridge, Mass., September, 2010 (hereinafter “Hung Paper”), which is hereby incorporated by reference in its entirety herein.

It is also possible to model the drag force (Fdi) produced by the surrounding fluid medium as a concentrated load applied at the gravity center of each element, with its direction determined by the flow direction of the surrounding fluid medium. (See FIG. 10 and the Example).
A model of the motion can be established using high speed photographic techniques for the IPMC beam vibrating in air and water (stagnant). The sensor output can include an open-circuit voltage or a short-circuit current and can be expressed as a function of tip displacement for the beam.
In one example, the IPMC beam has a length of between about 10 mm and about 30 mm In one example, the IPMC beam can be modeled in three dimensions in which the IPMC beam width can be monitored for torsional bending.
Opportunities to improve vehicle fuel economy and reduce emissions depend on improved sensors. A beam sensor as described herein can be used to control and monitor power train (ignition, fuelling and mass air flow), and thus improve fuel economy and reduce the emissions. For instance, precise metering of the air and fuel (including any type of biofuel, such as ethanol, E85, other alcohol blends, biodiesel, and the like) entering into the combustion chamber can affect the combustion efficiency and exhaust emissions (including hydrocarbon) during engine cold-start and warm-up conditions.
Traditional sensors are inadequate because they exhibit slow response time and typically provide only low-pass filtered data compiled over many engine cycles. In an automotive application, the engine operating conditions change rapidly.
The nature of cycle-by-cycle engine combustion can lead to cycle-to-cycle pulsating fuel flow. According to one example of the present subject matter, a beam sensor provides good calibration to ensure a proper fuel-air ratio and further provides real-time monitoring of fuel metering for each individual cylinder. An example of the present subject matter provides direct mechano-electrical transduction, flexibility (able to measure low-speed flow), lightweight, low cost and compact size.
FIG. 1 illustrates a schematic of a sensor having an electroactive polymer (EAP) member 10A according to one example. The EAP member 10A includes a free portion 15A and a fixed portion 20A. In this example, the EAP member 10A is configured as a cantilever (beam) with the free portion 15A capable of moving relative to the fixed portion 20A. Electrodes 25 and 27 are coupled to surfaces of the EAP member 1 OA and provide an electrical signal based on movement or deflection of the EAP member 10A.
The EAP member 10A is responsive to a fluid medium as well as to the flow properties of the fluid medium, represented by arrow 35. Both the flow properties and the fluid medium can act on the EAP member 10A and cause deflection, as indicated by arrow 40A.
An aspect ratio for a beam can be defined as a ratio of the beam length to beam width. In FIG. 1, the EAP member 10A has a length that is greater than the width and in one embodiment, the aspect ratio is greater than two. A high aspect ratio may exhibit greater sensitivity. Also as shown in FIG. 1, the beam length has an axis that is substantially perpendicular to the flow direction. The thickness of the beam (measured in a dimension substantially perpendicular to the width) can be uniform or non-uniform. For example, the thickness can be tailored to correspond with non-uniform flow within the cross-section of a flow channel.
The EAP member 10A can be configured to provide an output signal corresponding to static deflection or a dynamic frequency response.
FIG. 2 illustrates a flow sensor having an EAP member 10B disposed within a flow channel 29A according to one example. The flow channel 29A can carry a fluid (such as a liquid or a gas) in direction denoted by arrow 35. The EAP member 10B is electrically coupled by electrodes terminating at nodes 28 and 26.
The flow channel 29A can be configured for installation in a hydraulic system, a fuel system, an airway, a pneumatic line, a biological channel (such as an arterial system) or other fluid carrying conduit. The EAP member 10B, as with EAP member 10A, illustrates a cantilevered element.
FIG. 3 illustrates a sensor having an EAP member 10C suspended between two endpoints. In particular, the EAP member 10C includes a free portion 15B suspended between a fixed portion 20B and a fixed portion 20C. Fixed portions 20B and 20C, respectively, are rigidly coupled to structure 30B and 30C, respectively, leaving the free portion 15B free to oscillate.
FIG. 4 illustrates a sensor having a central disk-shaped portion with an EAP member 10D, having a free portion 15C. The free portion 15C is circular and, in this example, is free to oscillate in an upwardly or downwardly direction and/or twist in a torsional manner, as depicted by opposing arrows 40B and 40C. The EAP member 10D can also move in a complex manner, depending on flow detected relative to the free portion 15C.
FIG. 5 illustrates a sensor system 50 according to one example. The sensor system 50, in the example shown, includes an array of four sensors and a corresponding number of electrical connections for those sensors. In particular, the sensor system 50 can be configured for use in a four-cylinder engine. In this embodiment, sensors 29C, 29D, 29E, and 29F are coupled to electrical nodes 23C, 23D, 23E, and 23F, respectively.
FIG. 6 illustrates a model of a sensor system 60, according to one example. The sensor system 60 includes a single EAP member 10E, although more than one EAP member is also possible. In this embodiment, the EAP member 10E is electrically coupled to a processor 62. The processor 62 executes instructions stored in memory 64 and generates data for storage in the memory 64. In addition, the processor 62 is coupled to an output 66, which can include, but is not limited to, a printer, a display, a server, a transceiver, or other output device.
An example of two different states of an IPMC member, 700A and 700B, is illustrated in FIGS. 7A (without deformation) and 7B (under deformation), respectively. As shown, the IPMC without and with deformation (700A and 700B respectively), has a polymer layer 1 sandwiched by outer metal surface layers (electrodes) 2. In one embodiment, the polymer layer 1 is an ion-exchange polymer membrane (e.g., Nafion®) containing hydrated mobile cations 3 having mobile cations 4 therein, fixed anions 5 and mobile (neutral) water molecules 6. Fixed anions 5 within the IPMC member without deformation700A in FIG. 7A and covalently fixed to polymer chains in the polymer layer 1 are substantially balanced by the mobile cations 4 within the hydrated mobile cations 3. In contrast, the IPMC member under deformation 700B, such as a deformation under a mechanical stimulus redistributes the hydrated mobile cations 3 (and thus the mobile cations 4) (as well as the waater molecules 6 to a certain extent), producing a detectable electrical signal 7 correlated with the deformation.
FIG. 8 illustrates a flow chart of a method 80 according to one example. Method 80 represents a process for configuring a sensor. At 82, the method includes providing a flow channel. The flow channel can include a lumen or a pathway by which a fluid medium passes a sensor. At 84, the method includes positioning an electroactive polymer (EAP) member in the flow channel. This can include placing a cantilever in a manner in which the member axis is aligned perpendicular to an axis of the flow channel. At 86, the method includes electrically coupling the EAP member to a processor. At 88, the method includes configuring the processor to generate data. Configuring can include calibrating or programming the processor.
An EAP member can provide an electrical signal in response to relatively small motion of the EAP member. The electrical signal from the EAP member will change based on a detected change in oscillation frequency or amplitude. In particular, a large deflection of the member will produce a large electrical signal.
The EAP member can be used to detect a change in the fluid medium. For example, air, water, and fuel each produce a different signal pattern. This difference in signal pattern can be used to identify the fluid medium. In one example, a sensor in a fuel line of a flexible fuel vehicle can provide information to adjust engine performance based on a detected change in the fluid medium. In particular, when the sensor detects that a fuel known as E-85 is passing the member, a first engine parameter can be operative and when the sensor detect that a different fuel is passing the EAP member, then a second engine parameter can be operative. The output signal is distinct for different types of fluid media or flow conditions.
The EAP member can detect a change in fluid medium based on a change in viscosity, density, or other characteristic.
In addition, the EAP member can detect a pressure or temperature of the fluid medium and adjust calibration accordingly. In one example, the EAP member provides a signal for use by the engine control unit (ECU) that is used to adjust engine performance. Other sensors can also be integrated with the EAP member of the present subject matter in a system for operating an engine.
The present subject matter can be used to provide near real time data as to the fluid medium or flow characteristics. In particular, the sensor output is nearly instantaneous with a change detected by the EAP member.
The present subject matter can be configured to provide a signal based on a property of the surrounding fluid medium. This can include, for example, density or viscosity. In addition, the present subject matter can be configured to provide a signal based on a flow property. This can include, for example, a pulsation (duty cycle), turbulence, direction of flow, or other operational characteristic.
The EAP member can be configured to respond to a flow of a surrounding fluid medium confined to a lumen in a flow channel or flow of a fluid over a surface. For example, air flow over an airfoil or surface of an aircraft structure can be monitored using the present subject matter. Various examples can be configured as an angle of attack sensor, an airspeed indicator, a turbulence (laminar flow) detector, a stall detector, or other such sensor in an aircraft system.
As noted above, the dynamic response of an EAP member (or beam) can be determined In one embodiment, the dynamic responses of an IPMC beam subject to an evenly distributed external force due to the flow of the surrounding fluid medium can be captured. In this embodiment a full cantilever IPMC beam is discretized into a finite number of rigid-body elements as shown in FIG. 9A. FIG. 9B illustrates an example of discretizing the IPMC beam into five rigid elements. These elements have equal length, with each element linked with its neighbor through the connecting joints. The flexible nature of the material is modeled by using a spring (Ki) and a damper (Hi) at each connecting joint. The drag force (FDi) produced by the surrounding medium can be modeled as a concentrated load applied at the gravity center of each element, and its direction is determined by the flow direction of the fluid.
After discretization, the original cantilever IPMC beam can be converted into a multi-link dynamical system as discussed in the Hung Paper. In one embodiment, a five-element approximation can be used. See, for example, FIG. 10, which shows the free body diagram of one element. For each element, the governing equation can be written in the form as follow,




  
    
      
        
          
            
              J
              i
            
            ⁢
            
              
                θ
                 ̈
              
              i
            
          
          +
          
            
              H
              i
            
            ⁡
            
              (
              
                
                  
                    θ
                    .
                  
                  i
                
                -
                
                  
                    θ
                    .
                  
                  
                    i
                    -
                    1
                  
                
              
              )
            
          
          +
          
            
              K
              i
            
            ⁡
            
              (
              
                
                  θ
                  i
                
                -
                
                  θ
                  
                    i
                    -
                    1
                  
                
              
              )
            
          
          -
          
            
              H
              
                i
                +
                1
              
            
            ⁡
            
              (
              
                
                  
                    θ
                    .
                  
                  
                    i
                    +
                    1
                  
                
                -
                
                  
                    θ
                    .
                  
                  i
                
              
              )
            
          
          -
          
            
              K
              
                i
                +
                1
              
            
            ⁡
            
              (
              
                
                  θ
                  
                    i
                    +
                    1
                  
                
                -
                
                  θ
                  i
                
              
              )
            
          
        
        =
        
          
            
              F
              
                
                  (
                  
                    i
                    ,
                    
                      i
                      -
                      1
                    
                  
                  )
                
                ⁢
                x
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              θ
              i
            
            ⁢
            
              l
              2
            
          
          +
          
            
              F
              
                
                  (
                  
                    i
                    ,
                    
                      i
                      -
                      1
                    
                  
                  )
                
                ⁢
                y
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              θ
              i
            
            ⁢
            
              l
              2
            
          
          +
          
            
              F
              
                
                  (
                  
                    
                      i
                      +
                      1
                    
                    ,
                    i
                  
                  )
                
                ⁢
                x
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              θ
              i
            
            ⁢
            
              l
              2
            
          
          +
          
            
              F
              
                
                  (
                  
                    
                      i
                      +
                      1
                    
                    ,
                    i
                  
                  )
                
                ⁢
                y
              
            
            ⁢
            cos
            ⁢
            
              
            
            ⁢
            
              θ
              i
            
            ⁢
            
              l
              2
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




wherein l is the length of one beam element. F(i,i−1)x, F(i+1,i)x, F(i,i−1)y, F(i+1,i)y are the reaction forces each on the joint of (i, i−1) and (i+1,i), and they have the following relationship,

mi·aix=F(i,i−1)x−F(i+1,i)x  (2)

mi·aiy=FDi+F(i,i−1)y−F(i+1,i)y  (3)

wherein mi is the mass of beam i, aix and aiy are defined by equations (4) and (5),





  
    
      
        
          
            a
            ^
          
          i
        
        =
        
          
            
              
                a
                ix
              
              ⁢
              i
            
            +
            
              
                a
                iy
              
              ⁢
              j
            
          
          =
          
            
              
                a
                ^
              
              
                (
                
                  
                    i
                    -
                    1
                  
                  ,
                  i
                
                )
              
            
            +
            
              
                
                  θ
                  
                     ̈
                    ^
                  
                
                i
              
              ×
              
                
                  l
                  ^
                
                2
              
            
            +
            
              
                
                  θ
                  
                    .
                    ^
                  
                
                i
              
              ×
              
                
                  θ
                  
                    .
                    ^
                  
                
                i
              
              ×
              
                
                  l
                  ^
                
                2
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            a
            ^
          
          
            (
            
              
                i
                -
                1
              
              ,
              i
            
            )
          
        
        =
        
          
            
              a
              ^
            
            
              (
              
                
                  i
                  -
                  2
                
                ,
                
                  i
                  -
                  1
                
              
              )
            
          
          +
          
            
              
                θ
                
                   ̈
                  ^
                
              
              
                i
                -
                1
              
            
            ×
            
              l
              ^
            
          
          +
          
            
              
                θ
                
                  .
                  ^
                
              
              
                i
                -
                1
              
            
            ×
            
              
                θ
                
                  .
                  ^
                
              
              
                i
                -
                1
              
            
            ×
            
              l
              ^
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




wherein âi is the acceleration at the mass center of beam i; â(i−1,i) is the acceleration at the joint of beam i−1 and beam i. Note the acceleration â(0,1) at joint (0, 1) is zero (fixed). The force applied on each element is exerted by the surrounding medium consists of two terms, namely





  
    
      
        
          F
          Di
        
        =
        
          
            b
            ·
            l
            ·
            
              C
              D
            
            ·
            ρ
            ·
            
              
                V
                2
              
              2
            
          
          +
          
            M
            ·
            
              
                d
                ⁢
                
                  
                
                ⁢
                V
              
              
                d
                ⁢
                
                  
                
                ⁢
                t
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




wherein the first term is the drag force and the second term is the hydrodynamic force caused by the added-mass effect due to the beam vibration under a fluid medium. Parameters b and c are the beam element width and thickness. Note that V=V0+Vi, wherein V0 is fluid velocity, Vi=Vyi cos θi+Vxi sin θi and Vxi and Vyi are the tangential and normal components of the element linear velocity at its center of mass, respectively. Assume that angle θi is relatively small, and Vi can be approximated by Vyi. Also note that CD in equation (6) is the drag coefficient and ρ is the density of the surrounding medium. M represents the added mass that is being pushed by the beam element, and it can be calculated by M=πρb2l/4.

The beam rotational inertial Ji, rotational stiffness Ki, and the damping coefficient Hi can be calculated from the beam dimensions and the properties of the material, as defined below




  
    
      
        
          J
          i
        
        =
        
          
            
              m
              i
            
            ·
            
              l
              2
            
          
          12
        
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          K
          i
        
        =
        
          
            γ
            ·
            b
            ·
            
              c
              3
            
          
          
            12
            ⁢
            l
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          H
          i
        
        =
        
          
            ξ
            ·
            2
          
          ⁢
          
            
              
                J
                i
              
              ·
              
                K
                i
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




wherein γ is the Young's modulus of the IPMC beam material and ξ is the critical damping ratio of the material. The value of critical damping ratio ξ can be calibrated using the experimental results. Note that, each element can have the same rotational inertial, stiffness, and the damping coefficient.

In order to obtain a linearized model, a small angle assumption can be made such that

sin θi≈0, cos θi≈1  (10)

Additional assumptions which can be made include 1) that the added mass is very small and 2) that the fluid velocity is much larger than that of the beam at the center of the mass. Such assumptions imply that dV/dt is also small, therefore the added mass term can be ignored during the calculation. Based upon the above assumptions, equation (1) can be written in a matrix form as below,

J·{umlaut over (θ)}+H·{umlaut over (θ)}+K·θ=Γ  (11)

wherein J, H, and K are the matrices of moment of inertia, damping coefficient, and stiffness, respectively, and for a five-element case, they are given below





  
    
      
        
          J
          =
          
            [
            
              
                
                  
                    52
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    42
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    30
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
              
              
                
                  
                    42
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    40
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    30
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
              
              
                
                  
                    30
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    30
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    28
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
              
              
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    18
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    16
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
              
              
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    6
                    ⁢
                    
                      J
                      i
                    
                  
                
                
                  
                    4
                    ⁢
                    
                      J
                      i
                    
                  
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          H
          =
          
            [
            
              
                
                  
                    2
                    ⁢
                    
                      H
                      i
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      H
                      i
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      H
                      i
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      H
                      i
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    -
                    
                      H
                      i
                    
                  
                
                
                  
                    H
                    i
                  
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          K
          =
          
            [
            
              
                
                  
                    2
                    ⁢
                    
                      K
                      i
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      K
                      i
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      K
                      i
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      K
                      i
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
              
              
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  
                    -
                    
                      K
                      i
                    
                  
                
                
                  
                    2
                    ⁢
                    
                      K
                      i
                    
                  
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          
            θ
            =
            
              [
              
                
                  
                    
                      θ
                      1
                    
                  
                
                
                  
                    
                      θ
                      2
                    
                  
                
                
                  
                    
                      θ
                      3
                    
                  
                
                
                  
                    
                      θ
                      4
                    
                  
                
                
                  
                    
                      θ
                      5
                    
                  
                
              
              ]
            
          
          ,
          
            Γ
            =
            
              
                1
                2
              
              ⁡
              
                [
                
                  
                    
                      
                        
                          F
                          
                            D
                            ⁢
                            
                              
                            
                            ⁢
                            1
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              2
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              3
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              4
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              5
                            
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          F
                          
                            D
                            ⁢
                            
                              
                            
                            ⁢
                            2
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              3
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              4
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              5
                            
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          F
                          
                            D
                            ⁢
                            
                              
                            
                            ⁢
                            3
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              4
                            
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              5
                            
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          F
                          
                            D
                            ⁢
                            
                              
                            
                            ⁢
                            4
                          
                        
                        +
                        
                          2
                          ⁢
                          
                            F
                            
                              D
                              ⁢
                              
                                
                              
                              ⁢
                              5
                            
                          
                        
                      
                    
                  
                  
                    
                      
                        F
                        
                          D
                          ⁢
                          
                            
                          
                          ⁢
                          5
                        
                      
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




In one embodiment, a least squares parameter identification is employed. In this embodiment, all the IPMC beam model parameters and fluid flow velocity are assumed to be known, with the only unknown parameters left in equation (6) are CD and ρ. As such, Equation (11) can be rewritten into the following state space form,

{dot over (x)}=Ax+Bu, u=V2CDρ

y=Cx, y={dot over (θ)}5  (13)

wherein





  
    x
    =
    
      [
      
        
          
            θ
          
        
        
          
            
              θ
              .
            
          
        
      
      ]
    
  
  ,
  
    A
    =
    
      [
      
        
          
            
              0
              
                5
                ×
                5
              
            
          
          
            
              I
              
                5
                ×
                5
              
            
          
        
        
          
            
              
                -
                
                  J
                  
                    -
                    1
                  
                
              
              ×
              
                K
                
                  5
                  ×
                  5
                
              
            
          
          
            
              
                -
                
                  J
                  
                    -
                    1
                  
                
              
              ×
              
                H
                
                  5
                  ×
                  5
                
              
            
          
        
      
      ]
    
  
  ,
  
    B
    =
    
      [
      
        
          
            
              0
              
                5
                ×
                5
              
            
          
        
        
          
            
              
                J
                
                  -
                  1
                
              
              ×
              
                Q
                
                  5
                  ×
                  5
                
              
            
          
        
      
      ]
    
  
  ,




and

C=[01×91];

also I5×5 is an identity matrix and





  
    
      Q
      
        5
        ×
        5
      
    
    =
    
      
        
          bl
          2
        
        2
      
      ⁡
      
        [
        
          
            
              
                1
                /
                2
              
            
            
              1
            
            
              1
            
            
              1
            
            
              1
            
          
          
            
              0
            
            
              
                1
                /
                2
              
            
            
              1
            
            
              1
            
            
              1
            
          
          
            
              0
            
            
              0
            
            
              
                1
                /
                2
              
            
            
              1
            
            
              1
            
          
          
            
              0
            
            
              0
            
            
              0
            
            
              
                1
                /
                2
              
            
            
              1
            
          
          
            
              0
            
            
              0
            
            
              0
            
            
              0
            
            
              
                1
                /
                2
              
            
          
        
        ]
      
    
  
  ,
  
    
      V
      2
    
    =
    
      
        [
        
          
            
              
                
                  (
                  
                    
                      V
                      0
                    
                    +
                    
                      V
                      1
                    
                  
                  )
                
                2
              
            
          
          
            
              
                
                  (
                  
                    
                      V
                      0
                    
                    +
                    
                      V
                      2
                    
                  
                  )
                
                2
              
            
          
          
            
              
                
                  (
                  
                    
                      V
                      0
                    
                    +
                    
                      V
                      3
                    
                  
                  )
                
                2
              
            
          
          
            
              
                
                  (
                  
                    
                      V
                      0
                    
                    +
                    
                      V
                      4
                    
                  
                  )
                
                2
              
            
          
          
            
              
                
                  (
                  
                    
                      V
                      0
                    
                    +
                    
                      V
                      5
                    
                  
                  )
                
                2
              
            
          
        
        ]
      
      .
    
  




Note that Γ=CDρV2Q. In this embodiment, fluid velocity V0 can be obtained based upon the pressure difference in a tube and other measurement approaches. Therefore, in this embodiment, V0 is assumed to be known. Beam element velocities V1 to V5 at the center of mass and the beam tip velocity Vtip are defined below assuming small angle,

V1=l{dot over (θ)}1/2

V2=l{dot over (θ)}1+l{dot over (θ)}2/2

V3=l{dot over (θ)}1+l{dot over (θ)}2+l{dot over (θ)}3/2

V4=l{dot over (θ)}1+l{dot over (θ)}2+l{dot over (θ)}3+l{dot over (θ)}4/2

V5=l{dot over (θ)}1+l{dot over (θ)}2+l{dot over (θ)}3+l{dot over (θ)}4+l{dot over (θ)}5/2

Vtip=l{dot over (θ)}1+l{dot over (θ)}2+l{dot over (θ)}3+l{dot over (θ)}4+l{dot over (θ)}5  (5)

As mentioned above, the beam short current signal, ishort, is proportional to the beam tip velocity, Vtip, namely:

Vtip=η·ishort  (15)

wherein η is the proportional gain, which can be determined according to known methods or obtained using experimental data. In one embodiment, the parameter η is obtained based on the high speed imaging and the beam short current signal (See also

In addition, note that the damping matrix H is proportional to stiffness matrix in equation (11), or H=βK, where β is a scalar. The natural frequencies of each mode in system (11) can be obtained by using Rayleigh damping theory.
Equation (16) shows the system in modal coordinate,

{umlaut over (q)}+H{dot over (q)}+Kq=PTΓ  (16)

wherein q=PTθ is the uncoupled coordinates, H=PTHP and K=PTKP are transformed damping and stiffness matrices, respectively. The coordinate transformation matrix P=[P1 P2 P3 P4 P5] consists of five system mode shapes, Pi(i=1,2, . . . ,5). Since the force applied to each beam element is assumed to be evenly distributed at the same direction, we assumed that the first mode dominates the beam vibration motion. This can also be observed by the high speed imaging in the Hung Paper. Therefore, the first mode vibration is used for analysis. Note that the first mode shape of:

P1=[P11 P21 P31 P41 P51]T  (17)

is obtained from the coordinate transformation from beam physical coordinate to its modal coordinate. Under this assumption, beam element velocities V1 to V5 can be expressed in terms of the function of beam tip velocity Vtip as well as the beam tip angular velocity {dot over (θ)}5. Equation (18) provides the detailed solutions.





  
    
      
        
          
            V
            1
          
          =
          
            
              
                
                  P
                  11
                
                /
                2
              
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  P
                  51
                
              
            
            ⁢
            
              V
              tip
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            V
            2
          
          =
          
            
              
                
                  P
                  11
                
                +
                
                  
                    P
                    21
                  
                  /
                  2
                
              
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  P
                  51
                
              
            
            ⁢
            
              V
              tip
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            V
            3
          
          =
          
            
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  
                    P
                    31
                  
                  /
                  2
                
              
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  P
                  51
                
              
            
            ⁢
            
              V
              tip
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            V
            4
          
          =
          
            
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  
                    P
                    41
                  
                  /
                  2
                
              
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  P
                  51
                
              
            
            ⁢
            
              V
              tip
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            V
            5
          
          =
          
            
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  
                    P
                    51
                  
                  /
                  2
                
              
              
                
                  P
                  11
                
                +
                
                  P
                  21
                
                +
                
                  P
                  31
                
                +
                
                  P
                  41
                
                +
                
                  P
                  51
                
              
            
            ⁢
            
              V
              tip
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              θ
              .
            
            5
          
          =
          
            
              
                P
                51
              
              ·
              
                V
                tip
              
            
            
              
                (
                
                  
                    P
                    11
                  
                  +
                  
                    P
                    21
                  
                  +
                  
                    P
                    31
                  
                  +
                  
                    P
                    41
                  
                  +
                  
                    P
                    51
                  
                
                )
              
              ·
              l
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




The continuous time beam system model (13) was discretized into the following discrete state space model:

x(k+1)=Âx(k)+{circumflex over (B)}u(k)

y(k)=Ĉx(k)  (19)

with a sample step of 0.5 ms. The time response of the discrete model can be written as follows,

x(k)=Âkx(0)+Σj=0k−1Âk−j−1{circumflex over (B)}V2(j)CDρ

y(k)=Ĉx(k)  (20)

Let x0=x(0) and the following equation can be obtained based upon equation (20)





  
    
      
        
          
            
              
                y
                _
              
              =
              
              ⁢
              
                [
                
                  
                    
                      
                        y
                        ⁡
                        
                          (
                          1
                          )
                        
                      
                    
                  
                  
                    
                      
                        y
                        ⁡
                        
                          (
                          2
                          )
                        
                      
                    
                  
                  
                    
                      ⋮
                    
                  
                  
                    
                      
                        y
                        ⁡
                        
                          (
                          n
                          )
                        
                      
                    
                  
                
                ]
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  [
                  
                    
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              A
                              ^
                            
                            1
                          
                        
                      
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              A
                              ^
                            
                            0
                          
                          ⁢
                          
                            B
                            ^
                          
                          ⁢
                          
                            
                              V
                              2
                            
                            ⁡
                            
                              (
                              0
                              )
                            
                          
                        
                      
                    
                    
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              A
                              ^
                            
                            2
                          
                        
                      
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              ∑
                              
                                j
                                =
                                0
                              
                              1
                            
                            ⁢
                            
                              
                                
                                  A
                                  ^
                                
                                
                                  1
                                  -
                                  j
                                
                              
                              ⁢
                              
                                B
                                ^
                              
                              ⁢
                              
                                
                                  V
                                  2
                                
                                ⁡
                                
                                  (
                                  j
                                  )
                                
                              
                            
                          
                        
                      
                    
                    
                      
                        ⋮
                      
                      
                        ⋮
                      
                    
                    
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              A
                              ^
                            
                            n
                          
                        
                      
                      
                        
                          
                            C
                            ^
                          
                          ⁢
                          
                            
                              ∑
                              
                                j
                                =
                                0
                              
                              
                                n
                                -
                                1
                              
                            
                            ⁢
                            
                              
                                
                                  A
                                  ^
                                
                                
                                  n
                                  -
                                  j
                                  -
                                  1
                                
                              
                              ⁢
                              
                                B
                                ^
                              
                              ⁢
                              
                                
                                  V
                                  2
                                
                                ⁡
                                
                                  (
                                  j
                                  )
                                
                              
                            
                          
                        
                      
                    
                  
                  ]
                
                ⁡
                
                  [
                  
                    
                      
                        
                          x
                          0
                        
                      
                    
                    
                      
                        
                          
                            C
                            D
                          
                          ⁢
                          ρ
                        
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  Φ
                  ⁡
                  
                    (
                    n
                    )
                  
                
                ⁡
                
                  [
                  
                    
                      
                        
                          x
                          0
                        
                      
                    
                    
                      
                        
                          
                            C
                            D
                          
                          ⁢
                          ρ
                        
                      
                    
                  
                  ]
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




Since system matrices (Â,{circumflex over (B)},Ĉ) are known and velocity vector, V2, defined after equation (13) can be calculated using equations (15) and (18), assuming that the beam electronic current signal is measurable and the fluid speed V0 is known. Therefore matrix Φ(n) can be obtained based upon the beam current signal. To calculate the fluid property parameter CDρ, the least squares minimal solution of equation (21) is




  
    
      
        
          [
          
            
              
                
                  x
                  0
                
              
            
            
              
                
                  
                    C
                    D
                  
                  ⁢
                  ρ
                
              
            
          
          ]
        
        =
        
          
            
              [
              
                
                  
                    Φ
                    T
                  
                  ⁡
                  
                    (
                    n
                    )
                  
                
                ⁢
                
                  Φ
                  ⁡
                  
                    (
                    n
                    )
                  
                
              
              ]
            
            
              -
              1
            
          
          ⁢
          
            
              Φ
              T
            
            ⁡
            
              (
              n
              )
            
          
          ⁢
          
            y
            _
          
        
      
    
    
      
        (
        22
        )
      
    
  




The solution minimizes the following cost function





  
    
      
        
          
          
            
              y
              _
            
            -
            
              
                Φ
                ⁡
                
                  (
                  n
                  )
                
              
              ⁡
              
                [
                
                  
                    
                      
                        x
                        0
                      
                    
                  
                  
                    
                      
                        
                          C
                          D
                        
                        ⁢
                        ρ
                      
                    
                  
                
                ]
              
            
          
          
        
        2
      
    
    
      
        (
        23
        )
      
    
  




Then, fluid property parameter CDρ can be obtained from equation (22).

In one embodiment, a method comprising converting a polymer signal into flow information is provided. As shown in FIG. 11, the method can comprise sampling 102 a polymer short-circuit signal in fluid and using 104 the polymer (short-circuit) signal of the fluid to process the start and end of pulsating flow information. The next step can comprise selecting 106 a “useful” (short-circuit) signal (i.e., a signal that has a corresponding tip angular movement magnitude of less than about 10 degrees) from start and end of pulsating flow with a given length “n.” The method can continue by calculating 108 a drag coefficient based on selected data and a polymer beam model and then comparing 110 the identified parameter with the parameter obtained with previous “n.” A determination 112 is then made as to whether or not the difference between the two parameters is smaller than a calibration threshold. If not, the process continues with the selecting step 106 described above. If it is, the process is complete 114, and the flow information obtained can now be used to provide information as a sensor output signal to be used as a feedback signal for a control system. If the fuel blend content of E85 and gasoline is detected, it can be used to adjust the fuel injection quantity to regulate the engine air-to-fuel ratio accurately.
Embodiments of the invention will be further described by reference to the following examples, which are offered to further illustrate various embodiments of the present invention. It should be understood, however, that many variations and modifications may be made while remaining within the scope of the present invention.
EXAMPLE 1
In this example, fluid property parameter, CDρ, of the IPMC beam submerged in a fluid tube was calculated based upon the electronic signals of two fluid media, namely, distilled water and n-heptane (a standard test fuel of the straight-chain alkane with the chemical formula H3C(CH2)5CH3 or C7H16), as shown in FIGS. 12 and 13.
Table 1 defines the parameters for the IPMC beam used in this testing. Note that, the beam element length l is 1⁄5 of the total beam length.







TABLE 1







Parameters for Tested IPMC Beam
















Current








ishort and





Vtip


Young's
Material
Density
Proportion


Thick-


Modulus γ
damping
ρbeam
Gain η
Length l
Width b
ness


(N/mm2)
ratio ξ
(g/mm3)
(m/s/μA)
(mm)
(mm)
c (mm)





500
0.14
3
0.22
2
4
0.25









Table 2 provides the properties of the two fluids (water and n-heptane), as well as their measured average velocity. The Reynolds number, Re, of the two fluids is obtained by solving the following equation:




  
    R
    e
  
  =
  
    
      
        V
        0
      
      ⁢
      b
    
    v
  




(See, M. Potter and D. Wiggert, 1997, Mechanics of Fluid, 2nd ed., Prentice Hall, pp. 339-343 (hereinafter “Potter Reference”). Based upon the obtained Reynolds number, the actual drag coefficient CD can be found using the table in the Potter Reference. Note that the obtained drag coefficient assumes that the actual beam shape is a round cylinder. Considering the beam is rectangular, an interval value is given in Table 2. Next, the drag coefficient was calculated based upon the data for different liquids.








TABLE 2







Parameters for Tested Liquid














Kinematic
Average

Drag



Density ρ
viscosity υ
velocity
Reynolds
coefficient


Liquid
(g/mm3)
(mPa · s)
V0 (m/s)
number Re
CD















Water
1
 1.007 * 10−6
0.15
595.8
1.1 to 1.3


n-heptane
0.684
0.5643 * 10−6
0.22
1583.4
0.8 to 1  









The data shown in FIGS. 12 and 13 was used to estimate the fluid property CDρ of water and n-heptane. Note that parameter n, defined in equation (21) (above), affects the accuracy of the least squares estimation. Using the water test data of FIG. 12 and referring now to FIG. 14, the beam current signal under water at the start of pulse flow, and its corresponding beam tip velocity obtained through simulation can be seen. Due to the nonlinearity at the start of pulse flow, the beam current signal cannot match with the simulated beam tip velocity. The linear model also assumes a small displacement angle of the IPMC beam. Therefore, in order to precisely extract the fluid property, the current signal used for parameter estimation is small. Since the sensor transfer function between the current signal and tip displacement angle is known, beam tip angle displacement is configured to be less than 10 degrees. FIG. 14 shows the signal range selected for parameter estimation, where the current signal matches with the simulation beam tip velocity well.
FIG. 15 illustrates the estimation error in terms of the signal length n. The estimation error is respect to the final identified value of CDρ shown in Table 3. The estimation error reduces as the signal length increases because the estimation method is mainly dependent of the beam composite damping (beam and fluid damping), which requires certain signal length to be estimated precisely. FIGS. 16A, 16B and 16C show the selected data regions of the other three cases, with the final estimation results were summarized in Table 3.







TABLE 3







Minimal Least Squares Solution














Predicted Value




Injection
Fluid
CDρ
Identified CDρ
















Start
Water
1100-1300
1103.7



End
Water
1100-1300
1137.0



Start
n-Heptane
547.2-684  
1377.3



End
n-Heptane
547.2-684  
613.2










The predicted value (obtained based upon the beam and fluid medium physical properties) for CDρ was calculated based on the data in Table 2, with the identified CDρ obtained using the data together with the proposed least squares approach. The result shows that the significant deviation between the predicted value shown in Table 3 and estimated CDρ is from n-heptane data at the start of the pulse flow. The remaining results in Table 3 show that the estimated values are close to their corresponding predicted values. The error of the n-heptane suggests that the IPMC beam has different behavior under n-heptane from water, although this preliminary conclusion cannot be confirmed without further testing. These results also suggest using the end of pulse flow data for parameter estimation.
EXAMPLE 2
Prophetic
Future testing will determine possible differences in behavior of beam sensors, such as an IPMC beam sensor, when exposed to fluid mediums, such as flowing fluid mediums having different viscosities. In addition to various solvents, such as water and n-heptane, a variety of fuels such as alcohol and biodiesel blend may be tested. The beam sensors may also be tested in gaseous fluid mediums. Additional testing will also be performed to confirm possible benefits in using the end of pulse flow data for parameter estimation.
Conclusion
The various embodiments described herein allow, for the first time, estimation of fluid properties, including fluid flow properties, based on the current signal of an electroactive polymer (EAP) member surrounded by the fluid medium using a finite element model. In one embodiment, a five element finite model is used.
The above detailed description includes references to the accompanying drawings, which form a part of the detailed description. The drawings show, by way of illustration, specific embodiments in which the invention can be practiced. These embodiments are also referred to herein as “examples.” Such examples can include elements in addition to those shown or described. However, the present inventors also contemplate examples in which only those elements shown or described are provided. Moreover, the present inventors also contemplate examples using any combination or permutation of those elements shown or described (or one or more aspects thereof), either with respect to a particular example (or one or more aspects thereof), or with respect to other examples (or one or more aspects thereof) shown or described herein.
All publications, patents, and patent documents referred to in this document are incorporated by reference herein in their entirety, as though individually incorporated by reference. In the event of inconsistent usages between this document and those documents so incorporated by reference, the usage in the incorporated reference(s) should be considered supplementary to that of this document; for irreconcilable inconsistencies, the usage in this document controls.
In this document, the terms “a” or “an” are used, as is common in patent documents, to include one or more than one, independent of any other instances or usages of “at least one” or “one or more.” In this document, the term “or” is used to refer to a nonexclusive or, such that “A or B” includes “A but not B,” “B but not A,” and “A and B,” unless otherwise indicated. In the appended claims, the terms “including” and “in which” are used as the plain-English equivalents of the respective terms “comprising” and “wherein.” Also, in the following claims, the terms “including” and “comprising” are open-ended, that is, a system, device, article, or process that includes elements in addition to those listed after such a term in a claim are still deemed to fall within the scope of that claim. Moreover, in the following claims, the terms “first,” “second,” and “third,” etc. are used merely as labels, and are not intended to impose numerical requirements on their objects.
Method examples described herein can be machine or computer-implemented at least in part. Some examples can include a computer-readable medium or machine-readable medium encoded with instructions operable to configure an electronic device to perform methods as described in the above examples. An implementation of such methods can include code, such as microcode, assembly language code, a higher-level language code, or the like. Such code can include computer readable instructions for performing various methods. The code may form portions of computer program products. Further, the code can be tangibly stored on one or more volatile or non-volatile tangible computer-readable media, such as during execution or at other times. Examples of these tangible computer-readable media can include, but are not limited to, hard disks, removable magnetic disks, removable optical disks (e.g., compact disks and digital video disks), magnetic cassettes, memory cards or sticks, random access memories (RAMs), read only memories (ROMs), and the like.
The above description is intended to be illustrative, and not restrictive. For example, the above-described examples (or one or more aspects thereof) may be used in combination with each other. Other embodiments can be used, such as by one of ordinary skill in the art upon reviewing the above description. The Abstract is provided to comply with 37 C.F.R. §1.72(b), to allow the reader to quickly ascertain the nature of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. Also, in the above Detailed Description, various features may be grouped together to streamline the disclosure. This should not be interpreted as intending that an unclaimed disclosed feature is essential to any claim. Rather, inventive subject matter may lie in less than all features of a particular disclosed embodiment. Thus, the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment, and it is contemplated that such embodiments can be combined with each other in various combinations or permutations. The scope of the invention should be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled.

## Claims

1. A sensor comprising:
an electroactive polymer (EAP) member deformable in response to at least one of a flow characteristic of a fluid medium or a fluid characteristic of the fluid medium; and
first and second electric terminals coupled to the electroactive polymer member and configured to provide a current electrical signal based on deformation of the member, the current electrical signal corresponding to the flow or fluid characteristic.

2. The sensor of claim 1 wherein the electroactive polymer (EAP) member is selected from a conjugated polymer member, a dielectric elastomer member, an ionic polymer-metal composite (IPMC) member, and combinations thereof.

3. The sensor of claim 1 wherein the fluid medium is confined to a flow channel.

4. The sensor of claim 3 wherein the member comprises a single cantilevered beam having a free portion configured to be disposed within the fluid medium and to deform relative to a fixed portion or a double cantilevered beam having opposing ends fixed.

5. The sensor of claim 4 wherein the free portion has a length to width ratio greater than about 2:1.

6. The sensor of claim 4 wherein the flow channel has a flow axis and the free portion has a length aligned substantially perpendicular to the flow axis.

7. The sensor of claim 4 wherein the free portion is substantially circular.

8. The sensor of claim 1 wherein the fluid medium is a liquid medium.

9. The sensor of claim 1 wherein the electroactive polymer member comprises at least one beam.

10. The apparatus of claim 9 wherein said beam has a non-uniform thickness.

11. The sensor of claim 9 wherein the fluid characteristic is selected from density, viscosity and temperature, and the flow characteristic is selected from flow rate, flow type, pressure, duty cycle and duration.

12. The sensor of claim 11 wherein the current electrical signal is a short-circuit current electrical signal.

13. The sensor of claim 12 wherein the flow type is a pulsating flow having a start and an end and the short-circuit current electrical signal provides information on the pulsating flow by using a linear model.

14. The sensor of claim 12 comprising one beam, wherein the short-circuit current electrical signal has a corresponding tip angular movement magnitude of less than 10 degrees from the start to the end of the pulsating flow for a given beam length.

15. The sensor of claim 1 wherein the first electric terminal is coupled to a first surface of the electroactive polymer member and the second electric terminal is coupled to an opposing second surface of the electroactive polymer member.

16. The sensor of claim 1 wherein the fluid medium is a fuel.

17. The sensor of claim 16 wherein the fuel is a gaseous fuel or a liquid fuel.

18. A method comprising:
allowing a sensor comprising an electroactive polymer (EAP) member to deform in response to at least one of a flow characteristic of a fluid medium or a fluid characteristic of the fluid medium; and
coupling first and second electric terminals to the electroactive polymer member, wherein the first electric terminal and the second electric terminal are configured to provide a current electrical signal based on deformation of the member, the current electrical signal corresponding to the flow or fluid characteristic.

19. The method of claim 18 further comprising providing the fluid medium.

20. The method of claim 18 wherein the fluid medium is confined to a flow channel.

21. The method of claim 20 wherein the electroactive polymer member comprises at least one beam.

22. The method of claim 21 wherein said beam has a non-uniform thickness corresponding with a non-uniform flow within a cross-section of the flow channel.

23. The sensor of claim 1 wherein the deformation is a change in position which is a stationary displacement of the sensor or a dynamic response of the sensor.

24. The sensor of claim 1 wherein the fluid medium flows relative to another structure.

25. The sensor of claim 24 wherein the structure is a surface or a rotating body.