Abstract:
A method for determining the flow of an anode gas out of an anode sub-system. The method includes providing pressure measurements at predetermined sample times over a predetermined sample period and using the pressure measurements to calculate a slope of a line defining a change of the pressure from the beginning of the time period to the end of the time period. The slope of the pressure line is then used in a flow equation to determine the amount of gas that flows out of the anode sub-system, which can be through a valve or by system leaks.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    This invention relates generally to a method for determining flow out of a closed system and, more particularly, to a method for determining flow of anode gas out of an anode sub-system in a fuel cell system, where the method determines the slope of a pressure line formed from a plurality of pressure measurements and uses the slope in a flow calculation. 
         [0003]    2. Discussion of the Related Art 
         [0004]    Hydrogen is a very attractive fuel because it is clean and can be used to efficiently produce electricity in a fuel cell. A hydrogen fuel cell is an electro-chemical device that includes an anode and a cathode with an electrolyte therebetween. The anode receives hydrogen gas and the cathode receives oxygen or air. The hydrogen gas is dissociated at the anode catalyst to generate free protons and electrons. The protons pass through the electrolyte to the cathode. The protons react with the oxygen and the electrons at the cathode catalyst to generate water. The electrons from the anode cannot pass through the electrolyte, and thus are directed through a load to perform work before being sent to the cathode. 
         [0005]    Proton exchange membrane fuel cells (PEMFC) are a popular fuel cell for vehicles. The PEMFC generally includes a solid polymer electrolyte proton conducting membrane, such as a perfluorosulfonic acid membrane. The anode and cathode typically include finely divided catalytic particles, usually platinum (Pt), supported on carbon particles and mixed with an ionomer. The catalytic mixture is deposited on opposing sides of the membrane. The combination of the anode catalytic mixture, the cathode catalytic mixture and the membrane define a membrane electrode assembly (MEA). MEAs are relatively expensive to manufacture and require certain conditions for effective operation. 
         [0006]    Several fuel cells are typically combined in a fuel cell stack to generate the desired power. For example, a typical fuel cell stack for a vehicle may have two hundred or more stacked fuel cells. The fuel cell stack receives a cathode input gas, typically a flow of air forced through the stack by a compressor. Not all of the oxygen is consumed by the stack and some of the air is output as a cathode exhaust gas that may include water as a stack by-product. The fuel cell stack also receives an anode hydrogen input gas that flows into the anode side of the stack. 
         [0007]    A fuel cell stack includes a series of bipolar plates positioned between the several MEAs in the stack, where the bipolar plates and the MEAs are positioned between two end plates. The bipolar plates include an anode side and a cathode side for adjacent fuel cells in the stack. Anode gas flow channels are provided on the anode side of the bipolar plates that allow the anode reactant gas to flow to the respective MEA. Cathode gas flow channels are provided on the cathode side of the bipolar plates that allow the cathode reactant gas to flow to the respective MEA. One end plate includes anode gas flow channels, and the other end plate includes cathode gas flow channels. The bipolar plates and end plates are made of a conductive material, such as stainless steel or a conductive composite. The end plates conduct the electricity generated by the fuel cells out of the stack. The bipolar plates also include flow channels through which a cooling fluid flows. 
         [0008]    It is necessary to accurately determine the flow rate through bleed valves, drain valves, and possibly other valves, in the anode sub-system of a fuel cell system to know when to close the particular valve, as is well understood by those skilled in the art. Traditional valve orifice models work fairly well, but are subject to part-to-part variations because the models assume an effective area of the orifice. Further, the orifice model calculation also requires a difference between an inlet and outlet pressure to determine the flow. For certain known systems, this pressure differential is on the same order of magnitude as the error of the pressure sensors, which could lead to large estimation errors. 
         [0009]    U.S. patent application Ser. No. 12/636,276, titled Injector Flow Measurement for Fuel Cell Applications, filed Dec. 11, 2009, assigned to the assignee of this application and herein incorporated by reference, discloses a method for determining flow through a valve in a fuel cell system. An anode sub-system pressure is measured just before an injector pulse and just after the injector pulse, and a difference between the pressures is determined. This pressure difference, the volume of the anode sub-system, the ideal gas constant, the anode sub-system temperature, the fuel consumed from the reaction of the fuel cell stack during the injection event and the fuel cross-over through membranes in the fuel cells of the fuel cell stack are used to determine flow through a valve. The determination of the flow in the &#39;276 application uses two pressure measurements that could be affected by noise in the pressure readings and provides no feedback as to the validity of the pressure measurements. 
       SUMMARY OF THE INVENTION 
       [0010]    In accordance with the teachings of the present invention, a method is disclosed for determining the flow of an anode gas out of an anode sub-system. The method includes providing pressure measurements at predetermined sample times over a predetermined sample period and using the pressure measurements to calculate a slope of a line defining a change of the pressure from the beginning of the time period to the end of the time period. The slope of the pressure line is then used in a flow equation to determine the amount of gas that flows out of the anode sub-system, which can be through a valve or by system leaks. 
         [0011]    Additional features of the present invention will become apparent from the following description and appended claims, taken in conjunction with the accompanying drawings. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0012]      FIG. 1  is a schematic block diagram of a fuel cell system; and 
           [0013]      FIG. 2  is a graph with time on the horizontal axis and pressure on the vertical axis showing a change in pressure over time of an anode sub-system of a fuel cell system. 
       
    
    
     DETAILED DESCRIPTION OF THE EMBODIMENTS 
       [0014]    The following discussion of the embodiments of the invention directed to a method for determining the flow of an anode gas out of an anode sub-system is merely exemplary in nature, and is in no way intended to limit the invention or its applications or uses. Particularly, the method discussed herein determines a flow of an anode gas out of the anode sub-system. However, the method will have application for determining the flow out of other closed systems. 
         [0015]      FIG. 1  is a schematic block diagram of a fuel cell system  10  including a fuel cell stack  12 . Hydrogen gas from a high pressure hydrogen gas source  14 , such as a tank, is provided to the anode side of the fuel cell stack  12  on an anode input line  16 . The hydrogen gas from the source  14  is injected into the stack  12  by an injector  20 , where the injector  20  is intended to represent a single injector or a bank of injectors suitable for the purposes described herein. Anode exhaust gas from the fuel cell stack  12  is output on an anode exhaust line  22 . A bleed valve  24  is provided in the anode exhaust gas line  22 , and is periodically opened to bleed nitrogen from the anode side of the fuel cell stack  12  in a manner that is well understood by those skilled in the art. The system  10  also includes an anode recirculation line  34  that directs the anode exhaust in the anode exhaust line  22  back to the injector  20  to be provided on the anode input line  16  when the bleed valve  24  is closed. A drain valve  36  is provided in the recirculation line  34  and is operable to periodically drain water from the recirculated anode exhaust gas, also in a manner that is well understood by those skilled in the art. A pressure sensor  26  is provided in the recirculation line  34  to measure the pressure of the anode sub-system. The location of the pressure sensor  26  is intended to represent any suitable location for measuring the pressure at the anode inlet, the anode outlet or a recycle line between the anode inlet and outlet. 
         [0016]    The system  10  also includes a high temperature pump  28  that pumps a cooling fluid through a coolant loop  30  external to the stack  12  and through cooling fluid flow channels in the stack  12  in a manner that is well understood by those skilled in the art. A temperature sensor  32  measures the temperature of the cooling fluid flowing through the coolant loop  30 , and can be provided at any suitable location in the coolant loop  30 , such as at an inlet to the stack  12  where the cooling fluid is typically the coolest or at an outlet of the fuel cell stack  12  where the cooling fluid is typically the hottest. A controller  38  receives a pressure signal from the pressure sensor  26  and a temperature signal from the temperature sensor  32 , and controls the duty cycle of the injector  20  and the position of the bleed valve  24  and the drain valve  36 . It is necessary to know how much anode gas flows out of the valves  24  and  36  when they are opened so that it is known when to close the valves  24  and  36  for optimal system performance. The discussion below provides an accurate determination of the flow based on anode sub-system pressure. 
         [0017]      FIG. 2  is a graph with time on the horizontal axis and pressure on the vertical axis showing the anode sub-system pressure P on graph line  42  during and between injection events shown at graph locations  44  and  46 . An injection event at the locations  44  and  46  is identified by a sharp drop in the anode sub-system pressure P when the injector  20  is open and then a sharp rise in the anode sub-system pressure P when the injector  20  is closed. A pressure P 1  is defined just after the injector  20  is closed at one injection event at the location  44  and a pressure P 2  is defined just before the injector  20  is opened at the next injection event at the location  46 . A decay of the pressure P along the line  42  between the pressures P 1  and P 2  typically occurs for any number of reasons including a valve being opened and normal loss of pressure in the anode sub-system. For example, the bleed valve  24  can be opened during this time, the drain valve  36  can be opened during this time, there can be a leak in the closed anode sub-system if the valves  24  or  36  are closed, etc. 
         [0018]    As will be discussed in detail below, the slope of the line  42  between the injection events is used to determine the flow rate out of the anode sub-system for any of these occurrences. The controller  38  will open the bleed valve  24  or the drain valve  36  based on various system controls and operations know to those skilled in the art, and use the flow through these valves to accurately know when to close them. Also, if the anode sub-system is closed where the valves  24  and  36  are closed, the controller  38  can determine whether there is too large of a leak from the anode sub-system. Each point  40  along the line  42  represents a pressure measurement by the pressure sensor  26 . Those measurements can be taken at any suitable interval, such as every three milli-seconds. 
         [0019]    The slope {dot over (P)} of the pressure line  42  can be calculated along any segment of the line  42 , such as a sample time between P 1  and P 2 , as: 
         [0000]    
       
         
           
             
               
                 
                   
                     P 
                     . 
                   
                   = 
                   
                     
                       
                         ∑ 
                         
                           i 
                           = 
                           1 
                         
                         n 
                       
                        
                       
                         
                           ( 
                           
                             
                               t 
                               i 
                             
                             - 
                             
                               t 
                               _ 
                             
                           
                           ) 
                         
                          
                         
                           ( 
                           
                             
                               P 
                               i 
                             
                             - 
                             
                               P 
                               _ 
                             
                           
                           ) 
                         
                       
                     
                     
                       
                         ∑ 
                         
                           i 
                           = 
                           1 
                         
                         n 
                       
                        
                       
                         
                           ( 
                           
                             
                               t 
                               i 
                             
                             - 
                             
                               t 
                               _ 
                             
                           
                           ) 
                         
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Where t i  is the time that a pressure measurement is taken,  t  is an average of the sample times, P i  is a pressure reading at a particular sample time and  P  is the average of all of the measured pressures during the sample period. 
         [0020]    In order for the algorithm to determine whether the pressure measurement data is accurate enough to determine the flow out of the anode sub-system, the algorithm uses a pressure fit value r, as determined by equation (2) below, where the value r determines how close the pressure measurement data is to forming a straight line, and if r is less than a predetermined threshold, the pressure measurement data does not form a straight enough line, and thus is not accurate enough, and is not used to determine flow. 
         [0000]    
       
         
           
             
               
                 
                   
                     r 
                     2 
                   
                   = 
                   
                     
                       ( 
                       
                         
                           
                             ∑ 
                             
                               i 
                               = 
                               1 
                             
                             n 
                           
                            
                           
                             
                               ( 
                               
                                 
                                   t 
                                   i 
                                 
                                 - 
                                 
                                   t 
                                   _ 
                                 
                               
                               ) 
                             
                              
                             
                               ( 
                               
                                 
                                   P 
                                   i 
                                 
                                 - 
                                 
                                   P 
                                   _ 
                                 
                               
                               ) 
                             
                           
                         
                         
                           
                             
                               ∑ 
                               
                                 i 
                                 = 
                                 1 
                               
                               n 
                             
                              
                             
                               
                                 
                                   ( 
                                   
                                     
                                       t 
                                       i 
                                     
                                     - 
                                     
                                       t 
                                       _ 
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                                
                               
                                 
                                   ∑ 
                                   
                                     i 
                                     = 
                                     1 
                                   
                                   n 
                                 
                                  
                                 
                                   
                                     ( 
                                     
                                       
                                         P 
                                         i 
                                       
                                       - 
                                       
                                         P 
                                         _ 
                                       
                                     
                                     ) 
                                   
                                   2 
                                 
                               
                             
                           
                         
                       
                       ) 
                     
                     2 
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
           
         
       
     
         [0021]    Using this method, a number of variables are obtained to determine the validity of the calculated slope {dot over (P)}. The two major variables are the number of pressure measurement points n collected and the r 2  value of the linear fit compared to the data. This is important because the pressure response does not become linear until some finite time after the injector  20  was commanded close. This decay time is variable depending on a number of factors affecting the injector  20 , such as supply pressure, supply voltage, coil temperature, etc. The closing time of the injector  20  is not modeled and therefore if it changes, the algorithm collects an incorrect pressure point so it is important to ignore the resulting value. 
         [0022]    The slope {dot over (P)} of the pressure line  42 , as calculated by the equation (1), can be used to determine the flow through a valve, such as the bleed valve  24  and the drain valve  36 , when those valves are open, or can be used to determine if there is a leak in the anode sub-system where a reduction in anode pressure indicates a flow of the anode gas out of the anode sub-system over time. Thus, if that flow out of the system is greater than some threshold indicating abnormal stack operation, a diagnostic flag can be set that a leak is occurring. Of course, the example herein is for flow out of the anode sub-system of the fuel cell system  10 . However, using the pressure measurement to determine the slope of the pressure line can be used in any closed system where pressure can be measured. 
         [0023]    Using equation (3) below, the flow {dot over (n)} out of the anode sub-system can be determined as: 
         [0000]    
       
         
           
             
               
                 
                   
                     n 
                     . 
                   
                   = 
                   
                     
                       
                         
                           P 
                           . 
                         
                         · 
                         V 
                       
                       
                         R 
                         · 
                         T 
                       
                     
                     - 
                     
                       
                         I 
                         · 
                         
                           N 
                           cell 
                         
                       
                       
                         2 
                         · 
                         F 
                       
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   )