Measurement of mass of fluid in a container

Abstract: The mass of fluid remaining in a container 10 is determined by switching on a series of heating elements 14 and sensing the temperature of the container and contents over a period of time using sensor 16, thereby to deduce the thermal time constant. The thermal time constant is used to determine the mass of the remaining fluid contents from a look-up table or other means. In an alternative arrangement, the thermal response on cooling may be used in a similar manner.

Description: This invention relates to methods and apparatus for determining the mass of fluid in a container and in particular, but not exclusively, to determining the mass of liquid propellant remaining in a storage vessel on board a spacecraft. It is well known that a spacecraft includes thrusters for station-keeping in orbit and that these thrusters are supplied from a finite reservoir. In order to maximise the life of a spacecraft it is important to be able to determine very accurately the mass of propellant remaining at any given instant in time and particularly at the end of life when a small quantity of propellant must be reserved for the final graveyard burn. Traditionally the measurement of the remaining propellant has been accomplished by integrating the predicted propellant consumption for all the liquid propellant gauging and thruster firings over the spacecraft life with check calculations being performed from the tank pressure and temperature telemetry. This method suffers from the fact that the greatest accuracy is found at beginning of life, whilst at end of life, when a high accuracy is most needed, relatively low accuracy is achieved. Also, of course, what appears as a low percentage error when the vessel is full becomes a much larger percentage error when the tank is nearly empty. In practice this means that the spacecraft often will contain more propellant than is required for its service life and it will still have this usable excess even after the final graveyard burn. Thus an improvement in the accuracy of measurement of the mass of the propellant fluid may provide an increase in the service life and/or a decrease in the mass of fluid loaded at launch. In an earlier proposal, the temperature rise over a fixed period immediately after a tank heater is switched on is detected and related in a linear fashion to the thermal capacity of the system. JP-A-61122524 discloses such a system in which the temperature rise is detected and used in conjunction with the specific heat to determine the mass of fluid. We have designed a mass measurement system believed to provide a significant improvement in accuracy of measurement of the mass of the propellant fluid. Accordingly, in one aspect this invention provides a method of determining the mass of fluid remaining in a container which method comprises monitoring the heating or cooling thermal response of said container and/or said fluid to deduce the time constant of the variation of temperature with time for said thermal response, and correlating said time constant with a predetermined relationship between time constant and fluid mass to determine the mass of fluid remaining. In another aspect, this invention provides apparatus for determining the mass of fluid remaining in a container, said apparatus comprising temperature sensor means for sensing the temperature of the container and/or the fluid contained therein, and processor means responsive to the output of said temperature sensor means for monitoring the heating or cooling thermal response of the container and/or the fluid to deduce the time constant of the variation of temperature with time for said thermal response, and correlating said time constant with a predetermined relationship between time constant and fluid mass to determine the mass of fluid remaining. The invention also extends to spacecraft making use of the method or incorporating apparatus defined above. A non-limiting example of the invention will now be described reference being made to the accompanying drawings in which: Figure 1 is a schematic view of a spacecraft propellant storage tank fitted with booster heaters, and Figure 2 shows a typical relationship between the time constant and the liquid mass remaining. Figure 1 illustrates one of several propellant storage tanks provided on a spacecraft. The tanks are loaded with a given mass of fuel before launch and supply propellant fluid to thrusters (not shown) which maintain the spacecraft in the required attitude and orientation during the service life of the spacecraft and are also used at the end of the service life of the spacecraft to send it on a graveyard burn. Each storage tank 10 is insulated from the spacecraft structure and has on its surface six high power electric heating elements 12 and two low power electric heating elements 14. In a typical example, the heaters together dissipate not less than 7 watts. A thermistor 16 is attached to the mid portion of the storage tank 10 and provides a signal representing the temperature of the storage tank and its fluid contents. The heater elements 12 and 14 and thermistor 16 are already provided on most conventional spacecraft to adjust the temperature of the storage tank 10 prior to firing of the thruster. In this example, the thermal response of the tank 20 following switch-on of the heaters 12 and 14 is monitored by a processor 18 which may be on board the spacecraft or supplied with telemetry output at a ground station. The processor 18 determines the time constant (τ) of the temperature versus time characteristic of the storage tank 10 and its fluid content. We have found that the time constant decreases substantially linearly with the mass of the fluid content (see Figure 2), so that the mass of fluid contents can be determined once the time constant is known or calculated. We have developed three different mathematical models from each of which the thermal time constant (τ) may be derived, as follows:- Equation 1 τ = T₂ + T₃ 2 - T₁ + T₂ 2 T₁ - T₂ t₁ - t₂ - T₂ - T₃ t₂ - t₃ where T 1, T 2, T₃ are the temperatures of the container and/or fluid at times t 1, t 2, t₃ respectively. Equation 2 τ = t(a) - t(b) ℓn T∞ - T(b) T∞ - T(a) where T∞ is the steady state temperature of the container and/or fluid attained and T(a),T(b) are the temperatures of the container and/or fluid at times t(a),t(b). Equation 3 τ = -t ℓn 1 - (T(t)-To) (T ∞ -To) where To,T∞ are the initial and steady state temperatures respectively of the container and/or fluid and T(t) is the temperature of the container and/or fluid at time t after start of heating. It will of course be understood that the time constant (τ) or its equivalent may be calculated or determined by monitoring the temperature of the storage tank in ways other than from one of Equations 1 to 3. It is possible to determine the time constant at any time after switch-on of the heaters but the accuracy of measurement of the time constant increases with time and it is preferred to wait for at least 8 or 12 hours before the temperature data from thermistor 16 is processed to determine the time constant. A significant advantage of this system is that the percentage accuracy of the measurement of the time constant - and thus the mass remaining - increases as the mass remaining decreases, so that the measurement technique is most accurate towards the end of the life of the spacecraft when accuracy is most critical. The spacecraft will experience a diurnal temperature variation which is approximately sinusoidal and which is preferably stripped from the raw temperature data from the thermistor before the time constant is calculated. The variation may be detected and processed to form an adaptive sinusoidal model which is used to strip the sinusoidal diurnal component from the sensed temperature data. The accuracy of the value determined for time constant depends on the amount of telemetry output and thus the resolution of the digital temperature signal. The data rate can be increased by switching on the heater elements at such a time so that, during the measurement period, the diurnal variation acts to maximise the increase with time of the sensed temperature. In practice, assuming that the measurement period starts 12 hours after switch-on, the heaters should be switched on when the diurnal temperature is at a maximum, so that there is maximum data output in the second twelve hour period after switch-on. The only first order parameter upon which the time constant is dependent is the full fraction. Other parameters enter the equations since they are functions of the absolute temperature. The variations of these parameters with temperature are accounted for in the post test data processing software and include:- i) radiative coupling from system to environment (αT⁴) ii) specific heat of:- tank + insulation :- liquid within the tank :- vapour within the tank :- cover gas within the tank iii) vapour pressure of the contained liquid iv) latent heat of vapourisation of the contained liquid. The processor 18 will have a look-up table or other means relating calculated time constant to the mass of the remaining fluid contents or the fill fraction and the relationship may be determined empirically before launch taking account of the above parameters. The relation between the time constant τ and the propellant mass is a standard text book relationship and varies from container to container as shown below. The relationship is: τ = Total thermal capacity of container + contents Total thermal coupling from container to environment The remaining liquid mass is then related to the total thermal capacity by:- Liquid mass = Total thermal capacity - tank thermal capacity Liquid specific heat From the above two equations it can be seen that the arithmetic relationship between the time constant τ and the liquid mass is dependent upon a large number of variables e.g. tank mass, tank thermal control, liquid type etc, and Figure 2 shows a typical relationship between the time constant τ and the liquid mass remaining. The processor 18 will also, as mentioned above, be able to determine the diurnal temperature variation and provide appropriate compensating data. Our studies show that the exemplary technique disclosed above should be capable of predicting the fuel/oxidant mass to 6% accuracy which is a significant improvement on existing techniques, so that the excess fuel/oxidant required may be significantly reduced with a resultant increase in life or decrease in fuel/oxidant payload. Another advantage, where there are several fuel and oxidant tanks, is that the technique allows the operator to know the distribution of fuel and oxidant across the tanks. The life of the spacecraft may be further improved by selective tank use to use up all the dynamic residuals (i.e. that usable portion of the fluid as opposed to the unusable static residual). Previously these have not been used as their location is usually unknown. It is also possible to adjust the pressure and temperature in the tanks to preferentially use fuel and oxidant if an excess of either exists, again increasing on-station lifetime.

Claims: A method of determining the mass of fluid remaining in a container (10) which method comprises monitoring the heating or cooling thermal response of said container (10) and/or said fluid to deduce the time constant of the variation of temperature with time for said thermal response, and correlating said time constant with a predetermined relationship between time constant and fluid mass to determine the mass of fluid remaining. A method according to Claim 1, wherein monitoring of the thermal response is carried out after a predetermined delay period following the start of said heating or cooling. A method according to Claim 1 or Claim 2, wherein the heating thermal response is monitored and said heating is performed by one or more heating elements (14) associated with said container (10), and the temperature of said container (10) is sensed by a temperature sensor (16). A method according to any preceding Claim, wherein an estimate of said time constant (τ) is deduced substantially in accordance with the following formula: τ = T₂ + T₃ 2 - T₁ + T₂ 2 T₁ - T₂ t₁ - t₂ - T₂ - T₃ t₂ - t₃ where T₁,T₂,T₃ are the temperatures of the container (10) and/or fluid at times t 1. t₂,t₃ respectively. A method according to Claim 3, wherein an estimate of said time constant (τ) is deduced substantially in accordance with the following formula: τ = t(a) - t(b) ℓn T∞ - T(b) T∞ - T(a) where T∞ is the steady state temperature attained by the container (10) and/or fluid and T(a), T(b) are the temperatures of the container and/or fluid at times t(a),t(b). A method according to Claim 3, wherein an estimate of said time constant (τ) is deduced substantially in accordance with the following formula: τ = -t ℓn 1 - (T(t)-To (T ∞ - To where To,T∞ are the initial and steady state temperatures respectively attained by the container (10) and/or fluid and T(t) is the temperature of the container (10) and/or fluid at time t after start of heating. A method according to any preceding Claim, wherein the thermal response is compensated for variations in ambient conditions. A method according to any preceding Claim, wherein the container (10) is mounted on a spacecraft and contains a propellant fluid. Apparatus for determining the mass of fluid remaining in a container (10), said apparatus comprising temperature sensor means (16) for sensing the temperature of the container and/or the fluid contained therein, and processor means (18) responsive to the output of said temperature sensor (16) means for monitoring the heating or cooling thermal response of the container (10) and/or the fluid to deduce the time constant of the variation of temperature with time for said thermal response and correlating said time constant with a predetermined relationship between time constant and fluid mass to determine the mass of fluid remaining. A spacecraft including apparatus according to Claim 9. A spacecraft having a plurality of containers (10) containing propellant gas under pressure and each having associated therewith apparatus according to Claim 9, wherein fluid may be transferred from one container to another on the basis of the values obtained for the masses of fluid remaining in the containers (10).

Inventors: HONSTVET IAN ANDREW BRITISH AE, KENDREW MARK STUART BRITISH AE, HONSTVET, IAN ANDREW BRITISH AEROSPACE PUBLIC, KENDREW, MARK STUART BRITISH AEROSPACE PUBLIC
Applicants: BRITISH AEROSPACE, BRITISH AEROSPACE PUBLIC LIMITED COMPANY
Application year: 1990
Main classification: G01F23
IPC codes: G01F23, G01N25, B64G1