Method and device for filling or withdrawing from a pressurized gas tank

Patent metadata

Field	Value
Publication number	US-10704737-B2
Application number	US-201615573246-A
Country	US
Kind code	B2
Filing date	May 4, 2016
Priority date	May 12, 2015
Publication date	Jul 7, 2020
Grant date	Jul 7, 2020

Abstract

A method for filling or withdrawing from a pressurized gas tank. The tank having a wall having a cylindrical overall shape with dimensions and thermophysical properties that are given and known. The method including the regulation of the flow rate of the introduced or withdrawn gas, and/or of the temperature of the introduced gas, to avoid a situation in which the tank reaches a given high temperature threshold or a given low temperature threshold. The method including a step of estimating, by calculating in real time, at least one tank temperature from: the average temperature of the tank wall, the maximum temperature reached by the tank wall, the minimum temperature reached by the tank wall, and in that the flowrate of gas or the temperature of the gas is regulated depending on the calculated tank temperature.

First claim

What is claimed is:

1. A method for filling with, or bleeding off, a pressurized gas tank containing a gas, the tank being delimited by a wall of cylindrical general shape having a longitudinal axis and determined and known dimensions and thermo-physical properties and comprising a liner, comprising a gas, wherein a control valve is connected to the tank,
the method comprising;
calculating, in real time at least one temperature of the tank selected from the group consisting of:
an average temperature of the wall of the tank Twall,average(r,t) as a function of time (t) and r being a radius coordinate starting from a longitudinal axis of the tank,
a maximum temperature reached by the wall of the tank Twall,max(t) as a function of time, and
a minimum temperature reached by the wall of the tank Twall, min(t) as a function of time; and

using the control valve to regulate an inlet flow rate of the gas into the tank, a bleed off flow rate of the gas out of the tank, and/or an temperature of the gas used to fill the tank, thereby preventing the pressurized tank from reaching a determined high temperature threshold or a determined low temperature threshold, wherein the regulation of the inlet flow rate, the bleed off flow rate and/or the temperature of the gas used to fill the tank is carried out as a function of said at least one calculated temperature of the tank.

2. The method of claim 1 , wherein, during filling, when the calculated temperature of the tank reaches a determined high threshold (HT), the inlet flow rate is decreased and/or the inlet temperature is decreased by thermal exchange with a source of cold.

3. The method of claim 1 , wherein, during filling, when the calculated temperature of the tank is less than a high threshold (HT) by a determined value, the inlet flow rate is increased and/or the inlet temperature and/or of the temperature of the tank is increased by thermal exchange with a source of heat.

4. The method of claim 1 , further comprising;
calculating a Richardson Number (Ri) for the gas in the tank as a function of time,
comparing the Richardson Number (Ri) calculated with a determined reference value (Vr) lying between 0.05 and 1.5 and, when the Richardson Number (Ri) calculated is less than determined reference value (Vr) the temperature of the gas in the tank is considered to be homogeneous, wherein the maximum temperature reached by the wall of the tank Twall, max(t) as a function of time is equal to the average temperature of the wall of the tank Twall,average(r,t) in contact with the gas as a function of time (t): Twall,max(t)=Twall,average(r=r_liner,t).

5. The method of claim 4 , wherein, during filling, when the Richardson Number (Ri) is greater than the reference value (Vr), the inlet flow rate is increased.

6. The method of claim 4 , further comprising regulating the inlet flow rate, wherein when the Richardson Number number (Ri) calculated is greater than determined reference value (Vr), the temperature of the gas in the tank is considered to be heterogeneous wherein the maximum temperature reached by the wall of the tank Twall,max(t) as a function of time is not equal to the average temperature of the wall of the tank in contact with the gas as a function of time (t) wherein at the liner, Twall,average(r=radius of the liner of the tank, t), the method comprises a step of increasing the inlet flow rate thereby decreasing the value of the Richardson Number (Ri) calculated below the determined reference value (Vr) and thus render the gas homogeneous in temperature.

7. The method of claim 4 , further comprising, before filling, a step of determining or detecting, by sensor(s), an initial temperature T(0) of the gas in the tank, an initial pressure P(0) of the gas in the tank, an initial average temperature of the wall of the tank Tw, average(0) and a step of determining an initial mass of gas in the tank m(0) and then, during filling
when the Richardson Number (Ri) calculated is less than determined reference value (Vr) the temperature of the gas in the tank ( 1 ) is considered to be homogeneous, and under these conditions, the method comprises, in the course of filling, a step of calculating an average temperature Tgas,average(t) of the gas in the tank in real time as a function of time and the average temperature of the wall Twall,average(r,t) in real time as a function of time (t) on the basis of a mass and enthalpy balance applied to the gas in the tank and on the basis also of an energy balance in the wall of the tank, of the equation of state of the gas, and of a balance of the thermal exchanges between the gas and the wall, and between the wall of the tank is the exterior, and/or
when the Richardson Number (Ri) calculated is greater than determined reference value (Vr) the temperature of the gas in the tank is considered to be heterogeneous, and under these conditions, the method comprises, in the course of filling, a step of calculating the average temperature Tgas,average(t) of the gas in the tank in real time as a function of time (t) and the average temperature of the wall Twall,average(r,t) in real time as a function of time on the basis of a mass and enthalpy balance applied to the gas in the tank and of an energy balance in the wall of the tank, of the equation of state of the gas, and of a balance of the thermal exchanges between the wall of the tank is the exterior,

method comprising a step of calculating the maximum temperature reached by the wall of the tank Twall,max(t) as a function of time, this maximum temperature reached by the wall of the tank Twall,max(t) as a function of time being obtained by correlation on the basis of the average temperature Tgas,average(t) of the gas in the tank calculated in real time as a function of time and as a function of the average temperature of the wall Twall,average(r,t) in real time as a function of time (t).

8. The method of claim 7 , further comprising
calculating the enthalpy hin(t) of the gas entering or exiting the tank as a function of time,
measuring or calculating the mass of gas m(t) introduced or bled off from the tank as a function of time or, respectively,
determining the pressure P(t) in the tank as a function of time,
determining the average temperature of the gas Tgas,average(t) at the instant t in the tank in degrees K, this average temperature Tgas,average(t) being expressed as a first-degree function of the average temperature of the gas T(t−1) at the previous instant (t−1) and of a coefficient of convective heat exchange between the gas and the internal wall of the tank at the instant (t−1) in W·m−2·K−1, in which the heat exchange coefficient kg(t−1) is given by the relation kg=(λg/Dint)·Nuint in which λg is the thermal conductivity of the gas in the tank in W·m−1·K−1, Dint is the internal diameter of the tank in meters and NuDint the Nusselt number of the gas in the tank (dimensionless), and in which the Nusselt number of the gas is expressed as a function of the Reynolds number (Redin) (dimensionless) relating to the forced convection in the tank and of the Rayleigh number (RaDint) (dimensionless) relating to the internal natural convection in the tank according to a formula NuDint=a·RaDintb+c·Redind in which a and c are dimensionless coefficients dependent on the ratio (Lint/Dint) between the internal length of the tank Lint in meters and the internal diameter of the tank Dint in meters and on the ratio (Dint/di) between the internal diameter of the tank Dint in meters and the diameter of the injector di in meters, a, b, c and d being dimensionless positive real numbers, a lying between 0 and 1, b lying between 0.2 and 0.5, c lying between 0 and 1 and d lying between 0.5 and 0.9.

9. The method of claim 7 , further compris

Assignees

• Air Liquide

Inventors

• Bourgeois Thomas
• Ammouri Fouad
• Weber Mathilde
• Vyazmina Elena

Classifications

• Y02E60/32

  Hydrogen storage · CPC title

• F17C5/06 Primary

  for filling with compressed gases · CPC title

• F17C2250/0495

  the indicated parameter is a converted measured parameter · CPC title

• F17C2250/0443

  Flow or movement of content · CPC title

• F17C2203/0663

  in form of fibers or filaments · CPC title