US20150308621A1 | 2015-10-29 | |||
US20180335181A1 | 2018-11-22 | |||
US20140311622A1 | 2014-10-23 | |||
US20200346554A1 | 2020-11-05 | |||
US20190271439A1 | 2019-09-05 | |||
US20190184847A1 | 2019-06-20 |
WHAT IS CLAIMED IS: 1. A method of filling a tank with gaseous fuel, comprising: delivering a gas from a filling station to the tank; while delivering the gas: communicating a gas temperature measurement and a gas pressure measurement from the tank to the filling station; based on the gas temperature measurement, computing a time remaining (tremain_Tgas) for a temperature of the gas in the tank to reach a predetermined target temperature; based on the gas temperature and gas pressure measurements, computing a time remaining (tremain_SOC) for a state of charge of the gas in the tank to reach a predetermined target state of charge; determining a difference between (tremain_Tgas) and (tremain_SOC); and modulating a pressure ramp rate of the gas being delivered so as to reduce the difference between (tremain_Tgas) and (tremain_SOC). 2. The method according to claim 1, wherein, the tank is a hydrogen vehicle gas tank. 3. A method of filling a tank with gaseous fuel, comprising: delivering a gas from a filling station to the tank, at a ramp pressure; communicating a gas temperature measurement and a gas pressure measurement from the tank to the filling station; increasing the ramp pressure at a pressure ramp rate (PRR_calculated) which is calculated using a preselected fueling protocol (PRR_calculated), until the ramp pressure exceeds a threshold value (P_threshold); once the ramp pressure exceeds P_threshold: based on the gas temperature and gas pressure measurements, computing a state of charge (SOC) based on density using an equation of state; computing a rate of change of SOC (SOCRR) over a lookback period (t_lookback), and based on SOCRR, computing a time remaining (tremain_SOC) for a state of charge of the gas in the tank to reach a predetermined target state of charge; computing the rate of change of the ramp pressure (PRR_lookback) over a lookback period (t_lookback), and based on PRR_lookback, computing a time remaining (tremain_PRR) for the ramp pressure to reach a predetermined value (P_ramp_target); determining a difference between (tremain_PRR) and (tremain_SOC); in response to this difference being a negative value, computing a new pressure ramp rate PRR_SOC based on tremain_SOC, so as to reduce an absolute value of the difference between tremain_PRR and tremain_SOC; comparing PRR_SOC to (PRR_calculated); and in response to PRR_SOC being lower than PRR_calcuated, setting the pressure ramp rate to PRR_SOC. 4. A method of filling a tank with gaseous fuel, comprising: delivering a gas from a filling station to the tank, at a ramp pressure: communicating a gas temperature measurement and a gas pressure measurement from the tank to the filling station; increasing the ramp pressure at a pressure ramp rate (PRR_calculated) calculated using a predetermined fueling protocol); once the ramp pressure exceeds P_threshold: based on the gas temperature and gas pressure measurements, computing a state of charge (SOC) based on density using an equation of state; computing a rate of change of SOC (SOCRR) over a lookback period (t_lookback), and based on SOCRR, computing a time remaining (tremain_SOC) for a state of charge of the gas in the tank to reach a predetermined target state of charge; computing a new pressure ramp rate PRR_SOC based on tremain_SOC; comparing PRR_SOC to PRR_calculated; and in response to PRR_SOC being lower than PRR_calcuated, setting the pressure ramp rate to PRR_SOC. 5. The method according to claim 3 or 4, wherein P_threshold is a preset value. 6. The method according to claim 3 or 4, wherein P_threshold is a value based on a percentage of a nominal working pressure (NWP) of the tank; 7. The method according to claim 3 or 4, wherein P_threshold is a function of P_ramp_target and a pressure drop measured between the ramp pressure and the pressure in the tank. 8. The method according to claim 3 or 4, wherein P_threshold is a function of P_ramp_target and a pressure drop measured between a station pressure and the pressure in the tank. 9. The method according to claim 3 or 4, wherein P_threshold is a function of P_ramp_target and a maximum pressure drop measured between the ramp pressure and the pressure in the tank. 10. The method according to claim 3 or 4, wherein P_threshold is a function of P_ramp_target and a maximum pressure drop measured between the station pressure (P_station) and the pressure in the tank. 11. The method according to claim 3 or 4, wherein P_threshold is a function of P_ramp_target, the current ramp pressure as a percentage of the P_ramp_target, and the SOC_target. 12. The method according to claim 3 or 4, further comprising computing SOCRR based on a linear regression fit of the SOC data over the lookback period t_lookback. 13. The method according to any of the previous claims, wherein the gaseous fuel storage tank is a hydrogen vehicle gas tank. 14. A method of filling a tank with gaseous fuel, comprising: delivering a gas from a filling station to the tank, at a ramp pressure: communicating a gas temperature measurement (Tgas_high) and a gas pressure measurement from the tank to the filling station; setting values for a set of parameters (a, b, Tgas_max and Tgas_target); increasing the ramp pressure at a pressure ramp rate (PRR_calculated) calculated using a predetermined fueling protocol; computing a pressure drop (DeltaP) as the difference between the ramp pressure (P_ramp) and the tank pressure; computing a maximum value of the pressure drop (DeltaP_max); computing a temperature threshold value T_threshold as (DeltaP_max) subtracted from Tgas_target; in response to the gas temperature being greater than T_threshold: computing a pressure ramp rate (PRR_threshold); computing an adaptable denominator (AD) as the maximum value of either "b" or the product of "a" and DeltaP; computing a new pressure ramp rate (PRR_throttle) as the maximum value of either zero or the product: PRR_threshold times (Tgas_target minus Tgas_high) divided by AD; comparing PRR_throttle to PRR_calculated; and in response to PRR_throttle being lower than PRR_calculated, setting the pressure ramp rate to PRR_throttle. 15. The method of claim 14, wherein: The tank is part of a compressed gas storage system having two or more tanks; and the gas temperature measurement (Tgas_high) represents the highest gas temperature of all of the two or more tanks in the compressed gas storage system. 16. The method of claim 14, wherein Tgas_target is an offset (DeltaT) lower than Tgas_max, or is a percentage of Tgas_max. 17. The method of claim 14, wherein DeltaP is calculated as the station pressure minus tank pressure; 18. The method of claim 14, further comprising computing T_threshold as DeltaP subtracted from the target temperature (Tgas_target); 19. The method of claim 14, further comprising computing T_threshold as a percentage of Tgas_target. 20. A method of filling a tank with gaseous fuel, comprising: delivering a gas from a filling station to the tank, at a ramp pressure; setting values for a set of parameters (a, b_min, TMAL, Tgas_max, Tgas_target, Tgas_smooth_threshold, Tgas_diff_factor, and Tgas_diff_multiplier); communicating a gas temperature measurement (Tgas_high) and a gas pressure measurement from the tank to the filling station; increasing the ramp pressure at a pressure ramp rate (PRR_calculated) calculated using a predetermined fueling protocol; computing a pressure drop (DeltaP) as the difference between the ramp pressure and the tank pressure; computing a maximum value of the pressure drop (DeltaP_max); computing a triple moving average of Tgas_high with each moving average having a length of TMAL, referred to as Tgas_smooth; computing a temperature threshold value T_threshold as DeltaP_max subtracted from Tgas_target; in response to Tgas_smooth being greater than T_threshold: computing a pressure ramp rate (PRR_threshold); in response to Tgas_smooth being greater than Tgas_smooth_threshold: computing Tgas_diff as the difference between Tgas_high and Tgas_smooth; computing a maximum value of Tgas_diff, defined as Tgas_diff_max; computing Tgas_offset as the product of Tgas_diff_factor and Tgas_diff_max; computing Tgas_target as Tgas_max minus Tgas_offset; computing parameter b as the maximum value of either b_min or the product of Tgas_offset_multiplier and Tgas_offset; computing an adaptable denominator AD as the maximum value of either parameter b or the product of: parameter a and DeltaP; computing a new pressure ramp rate (PRR_throttle) as the maximum value of either zero or the product of PRR_threshold multiplied by (Tgas_target minus Tgas_smooth) divided by AD; comparing PRR_throttle to PRR_calculated; and in response to PRR_throttle being less than PRR_calculated, setting the pressure ramp rate to PRR_throttle. 21. The method of claim 20, wherein: The tank is part of a compressed gas storage system having two or more tanks; and the gas temperature measurement (Tgas_high) represents the highest gas temperature of all of the two or more tanks in the compressed gas storage system. 22. The method of claim 20, wherein DeltaP is calculated as the station pressure minus the tank pressure. 23. The method of claim 20, wherein T_threshold is computed as DeltaP subtracted from Tgas_target. 24. The method of claim 20, wherein T_threshold is computed as a percentage of Tgas_target. 25. The method of claim 20, wherein Tgas_smooth is computed using a triple moving average, wherein each moving average length is an independent value. 26. The method of claim 20, wherein Tgas_smooth is computed using an alternative smoothing function of Tgas_high. 27. The method of claim 14 or claim 20, wherein PRR_threshold is computed PRR_calculated at a time that Tgas_high or Tgas_smooth initially exceeds T_threshold. 28. The method of claim 14 or claim 20, wherein PRR_threshold is computed as (P_final - P_min)/t_final. |
i = i +1 and Loop Back to Eq 1.4 ELSE (Eq 1.10): Calculate PRR according to the fueling protocol being utilized, e.g. PRR(i) = PRRcalculated(i) [0060] Flow Chart: FIG.2 is a flow chart showing the implementation of Eq.1.1 through 1.10 set forth above. It is noted that the algorithm in the flow chart may have a calculation frequency suitable for the particular application. One example calculation frequency is 1 Hz. [0061] FIG.3 illustrates the operation of the Lookback T gas Throttle Method. The illustration shows the current point of the fill as time i and fill history up to this point is represented by solid lines. The dashed portion of the line labeled (3), is a projection of T gas_high based on the rate of change, represented by TGASRR lookback . The dashed portion of the line labeled (1), is a projection of the SOC based on the rate of change, represented by SOCRRlookback. As can be seen, based on these projections, Tgas_high reaches T gas_target at a time sooner than when the SOC reaches SOC target . This time difference is represented by comparing tremain_Tgas to tremain_SOC. Because tremain_Tgas is shorter than tremain_SOC, the pressure ramp rate needs to be reduced. As long as tremain_Tgas is shorter than tremain_SOC, the pressure ramp rate is reduced by ΔPRR each time step until it reaches PRR min . When t remain_Tgas is longer than t remain_SOC , the pressure ramp rate is increased by ΔPRR each time step until it reaches PRRcalculate. The PRR throttle can never exceed PRR calculated. This process essentially regulates the pressure ramp rate so that Tgas_high reaches Tgas_target at the same time as SOC reaches SOCtarget. [0062] SOC Throttle: [0063] SOC throttle is a fueling method intended to reduce the fueling time for fills that are constrained by the pressure drop between the dispenser pressure and the CHSS pressure. This is especially important for fueling of heavy-duty vehicles where the mass flow rates can be substantially higher than with light duty vehicle fueling. The SOC Throttle algorithm is intended to be used in conjunction with fueling protocols that use pressure ramp rate control strategies, for example the table-based protocol and the MC Formula protocol in SAE J2601, or any other pressure ramp rate based control. The control strategies regulate the rate of change of the ramp pressure, which is the pressure the dispenser targets for its measured pressure at each point in time (or time step) during the fill. This control strategy can increase or decrease the pressure ramp rate based on current conditions, at least for approaches like MC Formula that dynamically calculate the pressure ramp rate throughout the fill. Typically, the pressure ramp rate used in these control strategies is derived by ensuring the following four conditions are satisfied: a) the dispenser pressure does not exceed the maximum allowable pressure; b) the gas temperature within the CHSS does not exceed the maximum allowable temperature; c) the mass flow rate does not exceed the maximum value; and d) the SOC within the CHSS reaches an acceptable level, typically 95 to 100%. The problem with this approach is that often times, the derived pressure ramp rate is constrained by the requirement to meet the SOC target without the dispenser pressure exceeding the maximum dispenser pressure. The pressure ramp rate must be constrained such that pressure drop between the dispenser pressure and CHSS pressure is low enough to achieve the SOC target in the CHSS while at the same time not allowing the dispenser pressure to exceed the maximum value. This typically occurs under fueling conditions where the gas temperature in the CHSS does not approach the maximum allowable temperature, for example, under colder ambient temperatures or colder fuel delivery temperatures. If the pressure drop between the dispenser pressure and the CHSS pressure could be reduced, it would allow for faster fueling. However, the pressure drop is typically a function of the flow coefficients (or flow resistance) of the components utilized (i.e., the breakaway fitting, hose, nozzle, receptacle, and the tubing / manifolds / valves that make up the fuel delivery system of the CHSS). One way to reduce the pressure drop at the end of the fill is to fuel at a faster pressure ramp rate earlier in the fill so that more hydrogen is dispensed, allowing for the pressure ramp rate to be reduced towards the end of the fill. The SOC Throttle fueling methods implements this approach. [0064] The SOC Throttle operates by utilizing both the CHSS gas temperature and CHSS gas pressure communicated from the vehicle to the dispenser to calculate an SOC within the CHSS for each timestep throughout the fill, and then to throttle or reduce the pressure ramp rate (only if needed) so that the desired target SOC can be achieved at the end of the fill without exceeding the operating process limit for the dispenser pressure. The SOC is directly a function of the gas density. Using an equation of state, the gas density can be calculated from the gas pressure and gas temperature within the CHSS. The SOC is expressed as a percentage of the measured gas density in the CHSS to the gas density at reference conditions, typically at the nominal working pressure of the storage vessel and a temperature of 15 °C. As an example, the density of hydrogen at 100% SOC for a 70 MPa storage vessel is 40.2 g/l. [0065] The SOC Throttle fueling method allows the derivation of the average pressure ramp rate or the parameters used in calculating a dynamic pressure ramp rate (for example, the tfinal equation and associated a, b, c, d parameters used in SAE J2601 MC Formula) to be much faster than otherwise would be the case. This is made possible by eliminating one of the constraints in the derivation of the average pressure ramp rate or the parameters for a dynamic pressure ramp rate. The constraint that is eliminated is the requirement that the dispenser pressure not exceed the maximum allowable pressure. The remaining three constraints are that the gas temperature within the CHSS does not exceed the maximum allowable temperature, the flow rate does not exceed the maximum value, and the SOC within the CHSS reaches an acceptable level (typically 95 to 100%). For conditions where the CHSS gas temperature does not exceed the maximum allowable temperature, the peak flow rate becomes the constraining factor. Without the SOC Throttle fueling method applied, however, this pressure ramp rate would be much too fast at the end of the fill, causing a very large pressure drop between the dispenser pressure and the CHSS pressure, resulting in poor SOC. The SOC Throttle approach throttles back or reduces the pressure ramp rate once a pressure threshold value has been exceeded, only when required, and only sufficiently that the SOC target can be achieved without the dispenser pressure exceeding its maximum allowable pressure. [0066] There are two implementations of the SOC Throttle approach. These two implementations have the same objective, but differ slightly in their implementation. [0067] Definitions (Common to both example implementations): [0068] Tgas_low: the lowest bulk average gas temperature in the CHSS (communicated from the vehicle to the station) [0069] Pgas: the gas pressure in the CHSS (communicated from the vehicle to the station) [0070] P ramp : The ramp pressure for each time step. This is the pressure the dispenser is targeting and is a function of the ramp pressure one time step ago and the pressure ramp rate. [0071] Pramp_target: The end of fill maximum ramp pressure, typically 1.25 X NWP, although it can be a lower value than this to provide necessary margin. [0072] ΔP: The difference between the ramp pressure Pramp(i) and the pressure in the OR The difference between the station pressure Pstation(i) and the pressure in the CHSS PCHSS(i). [0073] ΔPmax: The maximum value of ΔP measured during the fill. [0074] P threshold : A threshold pressure above which the SOC Throttle approach is applied. P threshold can be a fixed value or a function. Examples of functions are shown below: SOC [0075] PRRcalculated: A calculated pressure ramp rate based on the fueling protocol utilized. This could be an average pressure ramp rate from a lookup table or a dynamically calculated pressure ramp rate based on a control parameter such as t final in the MC Formula protocol (i.e. PRR MC ) [0076] tlookback: A lookback time period in seconds. It can be a fixed value, a function of the pressure ramp rate, or be based on some other criteria such as a percentage of tfinal [0077] i: Represents the current time step. i – tlookback is the time step of the lookback period (i.e. tlookback seconds ago). i + 1 represents the next time step. A timestep is typically one second. [0078] SOC: The state of charge in the CHSS based on 100 times the calculated density divided by the density at nominal conditions (i.e. nominal working pressure and 15 °C temperature). The density is calculated based on an equation of state and is a function of Pgas and Tgas_low (these parameters are communicated from the vehicle to the station) [0079] SOCt arget: The desired end of fill SOC (eg 97%) [0080] SOCRR lookback : The rate of change of SOC over the lookback period [0081] tremain_SOC: The time remaining for SOC to reach SOCtarget based on SOCRRlookback [0082] PRR SOC : A pressure ramp rate calculated such that P ramp reaches the P ramp_target at the same time that the SOC reaches SOC target [0083] PRRlookback: The rate of change of the ramp pressure Pramp over the lookback period [0084] t remain_PRR : The time remaining for P ramp to reach P ramp_target based on PRR lookback [0085] Equations for Example Implementation 1: (Eq 2.1) : Calculate PRR calculated(i) : PRR calculated(i) = PRR calculated for the next time step according to the PRR control strategy being utilized (e.g. APRR for table-based method, PRRMC based on tfinal for MC Formula method) (Eq 2.2) IF Pramp(i) ≥ Pthreshold THEN (Eq 2.3) Calculate Lookback Period (rounded to nearest integer value) (Eq 2.4) (Eq 2.5) END IF (Eq 2.12) Calculate the pressure ramp for the next time step i + END IF [0086] Flow Chart: FIGS.4A and 4B contain a flow chart showing the implementation of Eq.2.1 through 2.14 set forth above. It is noted that the algorithm in the flow chart may have a calculation frequency suitable for the particular application. One example calculation frequency is 1 Hz. [0087] FIG.5 illustrates the operation of the SOC Throttle fueling method example implementation 1. The illustration shows the current point of the fill as time i and fill history up to this point is represented by solid lines. The dashed portion of the line labeled (2), is a projection of P ramp based on the pressure ramp rate calculated over the lookback period, i.e. PRRlookback. The dashed portion of the line labeled (1), is a projection of SOC based on the rate of change, represented by SOCRRlookback. As can be seen, based on these projections, P ramp reaches P ramp_target at a time sooner than when SOC reaches SOCtarget. If the fill continues at this pressure ramp rate, the ending SOC will be below the target SOC which is not desirable. This requires the pressure ramp rate to be reduced such that P ramp reaches P ramp_target at the same time that SOC reaches SOC target This is accomplished by calculating a new pressure ramp rate PRRSOC based on tremain_SOC. Finally, PRRSOC is compared to PRRcalculated and the slower of the two values is used as the pressure ramp rate PRR for the next time step. This approach allows the target SOC to always be reached in the fastest time possible, while also facilitating a slower pressure ramp rate if the calculated pressure ramp rate dictates that, for example, due to the fuel delivery temperature of the hydrogen warming (causing the enthalpy to increase and requiring an extension of the fill time). [0088] Equations for Example Implementation 2: (Eq 3.1): Calculate PRRcalculated(i): PRRcalculated(i) = PRR calculated for the next time step according to the PRR control strategy being utilized (e.g. APRR for table-based method, PRRMC based on tfinal for MC Formula method) (Eq 3.2): IF P ramp(i) ≥ P threshold THEN (Eq 3.3): Calculate Lookback Period (rounded to nearest integer value): (Eq 3.5): Calculate time remaining based on lookback SOC ramp rate: END IF [0089] Flow Chart: FIG.6 is a flow chart showing the implementation of Eq.3.1 through 3.10 set forth above. It is noted that the algorithm in the flow chart may have a calculation frequency suitable for the particular application. One example calculation frequency is 1 Hz. [0090] FIG.7 illustrates the operation of the SOC Throttle fueling method example implementation 2. The illustration shows the current point of the fill as time i and fill history up to this point is represented by solid lines. At this point, P ramp has just reached P threshold .. The rate of change in SOC is calculated based on a lookback period tlookback. This value is called SOCRRlookback. From this the time remaining for the SOC to reach its target value can be calculated. This value is t remain_SOC . A pressure ramp rate PRR SOC is calculated based on t remain_SOC , which is the pressure ramp rate at which Pramp will reach Pramp_target at the same time that SOC reaches SOCtarget Finally, PRRSOC is compared to PRRcalculated and the slower of the two values is used as the pressure ramp rate PRR for the next time step. This approach allows the target SOC to always be reached in the fastest time possible, while also facilitating a slower pressure ramp rate if the calculated pressure ramp rate dictates that, for example, due to the fuel delivery temperature of the hydrogen warming (causing the enthalpy to increase and requiring an extension of the fill time). [0091] Another fueling method may be referred to as an "Adaptable Tgas Throttle Method". This method provides a means to reduce and modulate the pressure ramp rate of a hydrogen fueling protocol to ensure that the maximum gas temperature limit of the CHSS is not exceeded, while also doing so in a manner that maximizes fueling performance (i.e. minimizes the fueling time). The method does this by calculating a threshold CHSS gas temperature above which the method is applied. Prior to this threshold being exceeded the fill utilizes the normal fueling protocol control (which can be a variety of approaches, including a fixed average pressure ramp rate, a dynamic pressure ramp rate such as used in the MC Formula protocol, etc.). This threshold temperature is dynamically calculated as a function of the pressure drop between the dispenser ramp press or the dispenser pressure measurement (station pressure) and the in-tank or CHSS gas pressure measurement. The reason the pressure drop is utilized is that when this pressure drop is high, changes in the dispenser pressure ramp rate have less effect on the mass flow rate or pressure ramp rate within the CHSS, which in turn affect the gas temperature rise. Therefore, when the pressure drop is high, it is necessary to utilize a lower threshold temperature to begin the pressure ramp rate reduction and modulation earlier in the fill. When the pressure drop is low, the threshold temperature is higher and the pressure ramp rate reduction / modulation occurs later in the fill. [0092] The Adaptable T gas Throttle method uses an equation which reduces and modulates the pressure ramp rate after the threshold temperature has been exceeded. This equation also uses an adaptable parameter in the denominator. The adaptable parameter is again a function of the pressure drop between the dispenser ramp pressure or station pressure and the CHSS gas pressure. The pressure ramp rate equation causes the pressure ramp rate to be reduced as the gas temperature rises. When the adaptable parameter in the denominator is large due to a high pressure drop, the reduction in the pressure ramp rate in relation to the rise in gas temperature is larger. As the pressure drop decreases later in the fill, the adaptable parameter in the denominator gets smaller, which makes the reduction in the pressure ramp rate in relation to the rise in gas temperature smaller. In essence, this equation is designed so that the pressure ramp rate is reduced sufficiently to avoid overshooting the maximum gas temperature in the CHSS, which can occur after the threshold temperature is breached and the pressure drop is still high (typically in the early to middle part of the fill), while at the same time keeping the pressure ramp rate high enough in the latter part of the fill so that the gas temperature approaches the target gas temperature. In summary, this approach uses an adaptable parameter to determine when the pressure ramp rate throttling begins and uses an adaptable parameter to regulate the pressure ramp rate once throttling begins. [0093] There are two implementations of the Adaptable Tgas Throttle methodology described. The first implementation is referred to simply as "Adaptable Tgas Throttle" and the second implementation is referred to as "Adaptable Tgas Throttle with Self- Adjusting Parameters." This second implementation calculates adjustments to some of the parameters based on the measured "noise" or "fluctuations" in the CHSS gas temperature Tgas_high. Tgas_high is supposed to represent the highest "bulk-average" gas temperature in each of the tanks in the CHSS. However, in practice the measurement of T gas_high is not a perfect representation of the bulk-average gas temperature and is influenced by many factors, including the placement of the temperature sensor in the tank and the amount of mixing of the gas in the tank, which is in turned influenced by the gas injector geometry and the aspect ratio of the tank. The problem these fluctuations in T gas_high present to the Adaptable Tgas Throttle control algorithm is that as Tgas_high approaches Tgas_target, Tgas_high can momentarily spike above Tgas_target. If T gas_target is set at or near T gas_max (the maximum allowed gas temperature in the CHSS), then the fill must stop. To counter this, T gas_target must be set lower, which causes the fueling time to be longer. There are two approaches to dealing with fluctuations in Tgas_high. In the first implementation of the Adaptable Tgas Throttle, the user defines sufficiently conservative parameter values (a, b, and T gas_target ) to deal with the inherent fluctuations in Tgas_high. The second approach is to implement a method whereby the fluctuations in Tgas_high are measured, and adjustments are automatically made to the parameters based on the amplitude of these fluctuations. This second approach is best suited to a hydrogen fueling station with many different makes and models of vehicles fueling, since the fluctuations in Tgas_high for each vehicle can be different, and fueling performance is automatically optimized for each vehicle. [0094] To reduce the fluctuations in the pressure ramp rate (PRR) caused by fluctuations in Tgas_high, a noise filter is implemented to smooth these fluctuations. A number of different noise filters were investigated and can be used, but a triple moving average (TMA) of the CHSS gas temperature demonstrated the best combination of effectiveness and simplicity. FIGS.8A and 8B illustrate the TMA (solid line), which is a reasonable approximation of the bulk average CHSS gas temperature (dashed line). The TMA does introduce a time lag, but this doesn’t have a material effect on the control since the precision of the PRR throttling is most important after the CHSS gas temperature has reached an asymptote during the latter part of the fill. The TMA is simply a moving average of a moving average of a moving average of the CHSS gas temperature. The length or periodicity of each moving average can be different, e.g.15, 10, 5, or it can be the same, e.g.10, 10, 10. Extensive computer fueling simulations were conducted and it was determined that a TMA of equal periodicities of 10 worked well. [0095] The effectiveness of the TMA in smoothing the CHSS gas temperature and consequently the PRR is illustrated in FIG.9. Note in FIG.9 that the Tgas_target value had to be set lower than 85 °C to avoid the CHSS gas temperature from momentarily spiking above the temperature limit. In this case, T gas_target was set to 83.25 °C, which resulted in a peak Tgas_high value of 84.5 °C. Although this methodology is effective, a shortcoming of it is that each vehicle will have different levels of fluctuations in T gas_high . So how should T gas_target be set if the level required depends on the magnitude of these fluctuations? One way to deal with this is to set Tgas_target sufficiently low so that it is effective against all expected fluctuation levels in Tgas_high (e.g. a noise level of +/- 5 °C). However, this approach increases the fueling time substantially, because the lower Tgas_target is set, the more the PRR is reduced, lengthening the fueling time. One option is for the vehicle OEM to measure the fluctuations in Tgas_high under a variety of fueling conditions and determine an appropriate T gas_target value, which can then be communicated from the vehicle to the dispenser. An alternative approach is to utilize an approach whereby the Tgas_target value automatically adjusts to the fluctuations inherent in T gas_high . [0096] There are two control parameters in the Adaptable Tgas Throttle algorithm which need to be adjusted due to fluctuations in Tgas_high: Tgas_target and the adaptable denominator value AD. Tgas_target is explained above. AD is the denominator in the PRR throttle equation. The smaller the value of AD, the more sensitive PRR is to changes in Tgas_high, or Tgas_smooth (which is the TMA applied to Tgas_high). Therefore, with higher amplitude in the fluctuations of Tgas_high, Tgas_target needs to be reduced and the minimum value of AD needs to be increased. The minimum value of AD is determined by the parameter "b". [0097] To determine the inherent fluctuations in Tgas_high, the difference between T gas_high and T gas_smooth is measured. This parameter is named T gas_diff . T gas_diff is measured after T gas_smooth crosses above a threshold temperature named Tgas_smooth_threshold. Once Tgas_diff begins to be measured, the maximum value is recorded as T diff This process is illustrated in FIG 10. [0098] The reason that T gas_diff is measured only after T gas_smooth rises above Tgas_smooth_threshold is because early in the fill, the CHSS gas temperature is rising rapidly. As noted previously, Tgas_smooth lags due to the TMA smoothing function. If T gas_diff is measured from the beginning of the fill, T gas_diff_max will be artificially high due to this lag. Therefore, the objective is to set the T gas_smooth_threshold value at a value where the CHSS gas temperature is naturally beginning to asymptote. Fueling simulations show that a value between 75 °C and 80 °C works well (although other values lower and higher could be used). Fueling performance is relatively insensitive to the Tgas_smooth_threshold value utilized within this range. This asymptote behavior in the CHSS gas temperature is illustrated in FIG.11. This region where Tgas_high and T gas_smooth are relatively flat is referred to as the throttling region and it is where the pressure ramp rate throttling is most critical to avoid exceeding the maximum gas temperature. [0099] T gas_diff_max is a measurement of the magnitude of fluctuation inherent in Tgas_high. Its purpose is to determine an appropriate setting for Tgas_target and the parameter "b", which determines the minimum value of AD. To utilize Tgas_diff_max in this manner, a derivative parameter T gas_offset is calculated. T gas_offset is calculated by multiplying T gas_diff_max by a parameter named T gas_offset_factor , i.e. T gas_offset = T gas_diff_max x Tgas_offset_factor. Tgas_target is then calculated as follows: Tgas_target = Tgas_max - Tgas_offset. Tgas_max is the maximum CHSS gas temperature allowed (typically 85 °C). When T gas_smooth_threshold has not yet been reached, T gas_target = T gas_max . "b" is calculated as follows: b = MAXIMUM [bmin, (Tgas_offset_multiplier x Tgas_offset)]. bmin is a user defined parameter, but simulations show a value of 4 works well. When Tgas_smooth_threshold has not yet been reached, b= b min . T gas_offset_factor and T gas_offset_multiplier are both user defined tuning parameters. Multiple fueling simulations were conducted to determine appropriate settings for these two parameters. These simulations demonstrated that T gas_offset_factor = 0.6 and T gas_offset_multiplier = 5 work well, but other values can be used. To illustrate the self-adjusting parameters in an example fill, see FIG.12. This graph shows how Tgas_target changes based on Tgas_diff_max and thus Tgas_offset. [0100] To confirm the robustness of the Adaptable Tgas Throttle with Self-Adjusting Parameters algorithm, approximately 150 computer fueling simulations were conducted at different ambient temperatures, fuel delivery temperatures, noise amplitudes in Tgas_high, initial CHSS pressures, CHSS Cv values, CHSS type 4 liner thermal conductivity values, CHSS type 3 liner properties, and CHSS surface to volume ratios. In other words, to test the robustness of this approach, all relevant parameters were varied over a wide range. In every simulation conducted, the peak CHSS gas temperature T gas_high was kept below 85 °C (which was used as the T gas_max setting). The highest peak gas temperature observed was 84.7 °C. [0101] Definitions applicable to both Adaptable Tgas Throttle and Adaptable Tgas Throttle with Self-Adjusting Parameters: [0102] T gas_high : the highest bulk average gas temperature in the CHSS (communicated from the vehicle to the station) [0103] Tgas_smooth: a triple moving average of Tgas_high using a moving average length of TMAL. This parameter is only calculated in the self-adjusting parameters approach. [0104] TMAL: The moving average length used in the triple moving average of T gas_high to calculate T gas_smooth. This parameter is only used in the self-adjusting parameters approach. [0105] Tgas_max: the maximum gas temperature allowed for the CHSS (can be a user setting or can be communicated from the vehicle to the station) [0106] T gas_target : the practical allowable gas temperature in the CHSS (accounting for any needed tolerances – can be a user setting or can be a function of Tgas_max) [0107] Tthreshold: Threshold temperature –when exceeded by Tgas_high or Tgas_smooth the PRR throttling algorithm is applied [0108] T gas_smooth_threshold : Another threshold temperature – when exceeded by Tgas_smooth Tgas_diff and Tgas_diff_max are calculated. This parameter is only used in the self-adjusting parameters approach. [0109] T gas_diff : A measurement of the difference between the current T gas_high and the current Tgas_smooth. This measurement is only conducted in the self-adjusting parameters approach. [0110] T gas_diff_max : The maximum T gas_diff measured during the fill up to and including the current timestep. This calculation is only conducted in the self-adjusting parameters approach. [0111] T gas_diff_factor : A user defined parameter utilized only in the self-adjusting parameters approach. [0112] Tgas_offset: A parameter which is calculated as a function of Tgas_diff_max and utilized only in the self-adjusting parameters approach. [0113] T gas_offset_multiplier : A user defined parameter utilized only in the self-adjusting parameters approach. [0114] PCHSS: The gas pressure in the CHSS. In a multi-tank CHSS, PCHSS is the lowest pressure in all of the tanks. (communicated from the vehicle to the station) [0115] Pstation: The dispenser pressure [0116] Pramp: The dispenser ramp pressure – the pressure the dispenser is targeting for each point in time during the fill [0117] ^P: A value representing the calculation of the current ramp pressure Pramp or current station pressure P station minus the current CHSS pressure P CHSS [0118] ^Pmax: The maximum ^P measured during the fill. [0119] Pfinal: The final pressure used in the derivation of the tfinal parameter used in MC Formula protocol, typically set at 1.25 x NWP [0120] P min : The initial pressure used in the derivation of the t final parameter used in MC Formula protocol [0121] t final : A parameter used in the MC Formula protocol. It is a calculation of the minimum time required to fill from Pmin to Pfinal without exceeding the CHSS maximum gas temperature limit and the maximum flow rate limit. It is typically derived using computer fueling simulations. [0122] PRR threshold : The pressure ramp rate at the time that T gas_high = T threshold or a pressure ramp rate based on a control parameter such as tfinal, e.g. PRRthreshold = (Pfinal – P min )/t final [0123] PRR: The pressure ramp rate for each time step throughout the fill. This is the control pressure ramp rate used to calculate the ramp pressure Pramp for each time step. [0124] PRRcalculated: The pressure ramp rate calculated by the fueling protocol, for example, PRR MC using the MC Formula fueling protocol algorithm. PRR calculated is calculated throughout the fill and is used as PRR prior to the point in the fill when Tgas_high or Tgas_smooth exceeds Tthreshold [0125] PRR throttle : The pressure ramp rate calculated by the Adaptable Tgas Throttle algorithm. [0126] a: A dimensionless user defined parameter, which is multiplied by ΔP. Parameter "a" is a user defined input which can be tuned for the best performance in a particular application. It can also be a function of another parameter such as the initial CHSS pressure. [0127] b: A dimensionless user parameter used to calculated a minimum value for AD. In the non-self adjusting parameter implementation "b" is a user defined input which can be tuned for the best performance in a particular application. [0128] bmin: A dimensionless user parameter used to calculate a minimum value of "b" and utilized only in the self-adjusting parameters approach. [0129] An example implementation of the Adaptable T gas Throttle fueling method is as follows: [0130] Equations: (Eq 4.1) Set a (user input – one embodiment IF P 0 < 10, a = 3, ELSE a = 4) Set b (user input – one embodiment b = 4) Set Tgas_max (user input or received via communications from vehicle – one embodiment Tgas_max = 85 ºC) Set T gas_target (user input or a function of T gas_max – one embodiment T gas_target = Tgas_max, another embodiment Tgas_target = Tgas_max – ^T where ^T is a user input) Begin Fueling (Eq 4.2) Calculate PRR calculated(i) according to the fueling protocol being utilized (Eq 4.3) ΔP(i) = Pramp(i) - PCHSS(i) (one embodiment); or ΔP(i) = Pstation(i) - PCHSS(i) (another embodiment) IF ΔP (i) > ΔP max , THEN ΔP max = ΔP (i) END IF (Eq 4.4) T threshold = T gas_target – aΔP max (one embodiment)
END IF Calculate the ramp pressure for the next time step i + 1: (Eq 4.12) i = i +1 and Loop Back to Eq 4.2 [0131] Flow Chart: FIGS.13A & 13B contain a flow chart showing the implementation of Eq.4.1 through 4.12 set forth above. It is noted that the algorithm in the flow chart may have a calculation frequency suitable for the particular application. One example calculation frequency is 1 Hz. [0132] FIG.14 illustrates the operation of the Adaptable T gas Throttle method. The illustration shows station pressure (which in this case also represents the ramp pressure), tank pressure, and tank temperature versus time. The line labeled (1), represents the pressure ramp rate. The line labeled (2), represents the station or ramp pressure. The line labeled (3), represents the tank pressure (PCHSS). The line, labeled (4), represents the threshold temperature. The line, labeled (5), represents the tank gas temperature (T gas_high ). The line, labeled (6), represents ΔP. The line, labeled (7), represents the product of a and ΔP. The line, labeled (8), represents AD. [0133] An example implementation of the Adaptable Tgas Throttle with Self- Adjusting Parameters fueling method is as follows: [0134] Equations: (Eq 5.1) Set T gas_smooth = T gas_high Set a (user input – one embodiment IF P0 < 10, a = 3, ELSE a = 4) Set bmin (user input – one embodiment bmin = 4) Set TMAL (user input – one embodiment TMAL = 10) Set T gas_max (user input or received via communications from vehicle – one embodiment Tgas_max = 85 ºC) Set T gas_target (user input or a function of T gas_max – one embodiment T gas_target = T gas_max ) Set Tgas_smooth_threshold (user input – one embodiment Tgas_smooth_threshold = 77 ºC) Set Tgas_diff_factor (user input – one embodiment Tgas_diff_factor = 0.6) Set T gas_offset_multiplier (user input – one embodiment T gas_offset_multiplier = 5) Begin Fueling (Eq 5.2) Calculate PRR calculated(i) according to the fueling protocol being utilized (Eq 5.3) embodiment) IF ΔP(i) > ΔPmax, THEN ΔPmax = ΔP(i) END IF (Eq 5.5) T threshold = T gas_target – aΔP max (one embodiment) Tthreshold = Tgas_target – aΔP (another embodiment) (Eq 5.6) IF T gas_smooth(i) ≥ T threshold THEN (Eq 5.7) PRR threshold(i) = PRR at the timestep when T gas_high ≥ T threshold for the first time (one embodiment); OR ^^ ^^ ^^ ( ^^ ^ ^ℎ ^^ ^^ ^^ℎ ^^ ^^ ^^( ^^) = ^^ ^^ ^^ ^^ ^^− ^^ ^^ ^^ ^^) ^ ^ ^^ ^^ ^^ ^^ ^^( ^^) (another embodiment) (Eq 5.8) IF Tgas_smooth(i) ≥ Tgas_smooth_threshold (Eq 5.9) THEN Tgas_diff(i) – Tgas_high(i) – Tgas_smooth(i) (Eq 5.10) ELSE
(Eq 5.19) END IF Calculate the ramp pressure for the next time step i + 1: (Eq 5.20) i = i +1 and Loop Back to Eq 5.2 [0135] Flow Chart: FIGS.15A & 15B contain a flow chart showing the implementation of Eq.5.1 through 5.20 set forth above. It is noted that the algorithm in the flow chart may have a calculation frequency suitable for the particular application. One example calculation frequency is 1 Hz. [0136] Implementation of The Fueling methods [0137] The fueling methods described above include algorithms that can be implemented by a hydrogen station dispenser. The elements that are required to practically implement these fueling methods are illustrated in FIG.16. [0138] FIG.16 illustrates a representative example filling station 10 in conjunction with a vehicle 12. The vehicle 12 includes a gaseous fuel storage tank 14 having a known volume "V" and equipped with a fill receptacle 15. The vehicle 12 includes a pressure sensor 16 operable to sense a gas pressure in the tank 14 and produce a signal representative thereof. The vehicle 12 includes a temperature sensor 18 operable to sense a gas temperature in the tank 14 and produce a signal representative thereof. The vehicle 12 includes an electronic controller 20 operable to receive the signals from the sensors and transmit that information along with other pertinent information such as the tank volume V using a wired or wireless data transmitter 22 such as the illustrated infrared transmitter. [0139] While FIG.16 illustrates a single tank, it will be understood that there can be multiple storage vessels that comprise the compressed hydrogen storage system. In the rendition with multiple storage vessels, the gas temperature measurement "T" would be used to determine the lowest gas temperature Tgas_low and the highest gas temperature T gas_high . [0140] Furthermore, it is noted that the use of the vehicle is merely an illustrative example, and the method described herein are suitable for application to a filling process for any gas storage container. [0141] The filling station 10 includes a fuel supply 24 (e.g. hydrogen). The fuel may be stored as liquid, low-pressure gas, or high-pressure gas, and is dispensed in gaseous form through a nozzle 26 which is disposed at a distal end of a fill hose 28 and which is configured to be coupled to the fill receptacle 15. It will be understood that a filling station 10 of this type may include conventional ancillary equipment for handling the fuel such as heat exchangers, pumps, compressors, and/or valves. The filling station 10 includes an electronic controller 30. [0142] The controller 30 includes one or more processors capable of executing ladder logic, programmed instructions, or some combination thereof. For example, it may be a general-purpose microcomputer of a known type, such as a PC-based computer, or may be a custom processor, or may incorporate one or more programmable logic controllers (PLC). [0143] The filling station 10 is equipped with a wired or wireless data receiver 32 configured to receive data from the data transmitter 22 of the vehicle 12, such as the illustrated infrared receiver mounted on the nozzle 26. The data receiver 32 is operably connected to the controller 30. [0144] The filling station 10 includes a gas property sensor 34 which may be disposed in a proximate end of the fill hose 28. The gas property sensor 34 is operable to sense one or more properties of the gaseous fuel flowing through the fill hose 28 generate a signal representative thereof. One example of a sensed property is pressure. Another example is flow rate (volume or mass). And another example is the gas temperature. [0145] The filling station 10 includes at least one throttling device operable to affect some aspect of the flow of gaseous fuel. One example of a throttling device is a controllable valve, shown schematically at 36. This could be, for example a pressure regulating valve operably connected to the controller 30 in such a way that the controller 30 can change the pressure set point of the pressure regulating valve. Another example of a controllable valve 36 would be a flow metering valve operated by an actuator operably connected to the controller 30 such that the controller 30 can change the flow rate through the flow metering valve. [0146] Another example of a throttling device is a variable-speed pump 38 operating control by the controller 30. [0147] Operation of the throttling device, e.g. modulation of a pressure, flow rate, or pump speed is effective to control a property of the gaseous fuel flow 34 measured by the gas property sensor 34. The controller 30 may be programmed to execute one or more of the throttle control algorithms described above and to operate the throttling device in a feedback loop to produce the desired pressure ramp rate (PRR) computed within the throttle control algorithm. [0148] The foregoing has described a gaseous fueling method. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and/or all of the steps of any method or process so disclosed, may be combined in any combination, except combinations where at least some of such features and/or steps are mutually exclusive. [0149] Each feature disclosed in this specification (including any accompanying claims, abstract and drawings) may be replaced by alternative features serving the same, equivalent or similar purpose, unless expressly stated otherwise. Thus, unless expressly stated otherwise, each feature disclosed is one example only of a generic series of equivalent or similar features. [0150] The invention is not restricted to the details of the foregoing embodiment(s). The invention extends to any novel one, or any novel combination, of the features disclosed in this specification (including any accompanying claims, abstract and drawings), or to any novel one, or any novel combination, of the steps of any method or process so disclosed.