Development of Liquid Capacity Measuring Algorithm Based on Data Integration from Multiple Sensors

This research proposes an algorithm using a process of integrating data from multiple sensors to measure the liquid capacity in real time regardless of the position of the liquid tank.The algorithm for measuring the capacity was created with a complementary filter using a Kalman filter in order to revise the level sensor data and IMU sensor data. The measuring precision of the proposed algorithm was assessed through repetitive experiments by varying the liquid capacity and the rotation angle of the liquid tank. The measurements of the capacity within the liquid tank were precise, even when the liquid tank was rotated. Using the proposed algorithm, one can obtain highly precise measurements, and it is affordable since an existing level sensor is used.


Introduction
Continuous or discontinuous measurements are used to assess the capacity of liquid contained in vehicles, tanker trucks, and industrial tanks.Continuous measurements are divided into the following types: buoyancy, weight measurement, pressure, capacitance, ultrasound, and radiation.Discontinuous measurements are divided into the following types: conductivity, heat transfer, capacitance, optical, ultrasound, and microwave [1][2][3].Some sensors used for measuring liquid level are summarized in Table 1.
The pressure level sensor, a double-continuous measuring type, forms the core of the measurement market owing to its affordability, simplicity of its gauge, and its ability to work at high temperature and high pressure.The ultrasound-and radar-type sensors are more efficient than others in measuring the methods but these sensors are not economical [4,5].
Most demand for current level sensors comes from water treatment plants (sewage, industrial water, wastewater, and black water), a national key industry, and the sector that measures the flow rate for vehicles and vessels.The level sensor is especially adequate for measuring the flow rate of vessels and the fuel of vehicles.
The existing method for measuring the vehicle fuel involves measuring poles and floating box [6].This is structurally limited because it is immeasurable when the fuel tank of the vehicle is below 10%.In addition, measuring errors occur because the fuel gets heeled over when the vehicle turns or accelerates as shown in Figure 1.
This paper proposes a new fuel-measuring algorithm using a pressure level sensor that can be applied for fuel control and measurements of vehicles.The system can accurately measure the fuel capacity below 10% by improving mechanical limits of the existing method in measuring the capacity.Furthermore, the measurement errors are minimized by continuously measuring and revising the slopping or heeling of fuel when the vehicle turns or accelerates.This revised algorithm guarantees precision for laser-and ultrasoundtype sensors.In addition, this algorithm is highly economical because an existing pressure level sensor [7] has been used.

Sensor Design
The proposed sensor should have a characteristic of a measuring algorithm of density and capacity for fuel measurement and a revised method in accordance with the fuel tank's condition, acceleration, and acceleration of gravity.
Figure 2 demonstrates the system which this research has proposed.A pressure level sensor is used to measure capacity, density, and tilting angle, and an inertial measurement unit (IMU) is used to revise the measuring value in accordance  with the rotation, acceleration, and acceleration of gravity [8].
The unique approach of this paper is that the relative value in accordance with the location of the level sensor and the rotation angle data of the IMU are both used to measure the position of a fuel tank.The reasons are as follows.The IMU outputs a precise measurable value for the rotation angle and acceleration measurements, but its disadvantage is that such values are accumulated with time by an integration constant.In the case of a level sensor, it can measure the tilting angle and acceleration in accordance with the measuring value difference located structurally, but the measurements may not be accurate.
The proposed algorithm of this paper guarantees precision and credibility because it uses the Kalman filter to realize a complementary filter with precision despite the error accumulation of the IMU and the imprecision and nonaccumulation of errors of the level sensor.

Sensor Measurement Algorithm for Precise
Measurements of Sensors

Measuring Algorithm of Capacity.
As demonstrated in Figure 3, the capacity can be computed in accordance with the container volume and the input capacity of the liquid.The mass can be computed using the pressure applied to the bottom of a container with liquid.First, the calculation equation of the capacity is as follows: Pressure applied to the bottom of the container ∝ liquid level × liquid density.
(1)   The equation to compute the force applied to the bottom in accordance with the liquid is as follows [17]: Therefore the pressure  applied to the unit area of the bottom side can be defined as in the equation below [2]: (If there is no internal pressure in the tank) ,  −  0 =  ×  × ℎ (If there is an internal pressure in the tank) .
(3) In this control system, the following method has been used to measure the capacity.
The sensor attached toward the -axis and -axis from the center of the target container is located on the opposite side of the container.The obtained measurement value is used to calculate the average value of each axis, and the = The average value of the output voltage → The measuring comparison by the input data. ( Here the output value of each sensor does not accumulate, but the precision is low.This sensor output value changes with turning and accelerating and is used in the revised algorithm to measure the capacity.

Measuring Tilting Angle with Level Sensors.
In Figure 4, the sensor installed on the opposite side of each axis receives In the equation above, Sensor 1-Sensor 6 are digital values obtained by an ADC board used for implementation, and ADC/ ref is the ADC value (constant) per tilting angle obtained from experiments.Sensors 5 and 6 are only used as a means to increase the accuracy because they duplicate axis of rotation.Therefore Sensors 5 and 6 are not used in the data analysis.

Measuring Algorithm of Density.
where  is liquid density,  is gravitational acceleration, ℎ is height, and  is pressure.The equation above is a basic equation to compute liquid density, and  1 and  2 from the above equation are Sensor 7 and Sensor 4 in Figure 5, respectively [18,19].

Multiple Sensors Algorithm to Improve the Precision of the Result Value of the Sensor
Figure 6 is a block diagram of the algorithm to accurately measure the capacity when the liquid tank has rotated by the force applied to it.An attitude Kalman filter computes accurate values, while accumulating almost no errors, by receiving the rotation angle data from the IMU and level sensor.Furthermore, a level Kalman filter precisely measures the capacity by using the obtained rotation angle data and level increment.The holistic drift reduction (HDR) filter used is excellent in avoiding drift in the obtained data from the IMU.The acceleration of gravity is distributed to the gravity value of each axis in accordance with the fuel tank conditions.

Revision of IMU Sensor.
When the output of the IMU sensor has an angular speed of ω ,   with a revised drift and scale is computed using the HDR algorithm: Here the HDR algorithm is used to compute the bias value   .

Attitude Kalman
Covariance matrix of Eulerian angles is as follows: State variables and covariance matrix are initialized as follows:

Attitude Prediction
(i) System Model.The angular speed measured by sensors has differentials of the rotation matrix and the following relation: Here, Ω  is the skew-symmetric matrix of   and is defined as follows: In addition, the differentials of Eulerian angles have angular speed and the following relation: Here,  −1 1 is a matrix that converts the angular speed measured by the gyro sensor into Eulerian angles: If  ≈ ±/2 in the matrix above, it is cos  ≈ 0, which generates the denominator of matrix element to be close to 0 at times.However, such a case refers to the fact that when the sensor stands vertically, this does not occur if one considers the location and direction when attaching the sensor.

Attitude Updates by Gravity (i) Attitude Measurement Model of Sensors of Gravity.
The matrix that indicates the attitude of sensors utilizing 0 and θ is revised as follows: Furthermore, Eulerian angles are revised as follows: Here,  −1 0 is as follows: The attitude measurement model function ℎ(⋅) of the sensor utilizes Eulerian angles without changes as follows: The attitude measurement value   is computed as follows, utilizing 0 and θ: (ii) Attitude Update.One has Here,  0 is column 1 row 1 element of   , and   is column 2 row 2 elements of   : Here,  = diag( 2 0 ,  2  ,  2  ) is computed as below because the benefit must increase when acceleration is near 1 g, and it should decrease when the size of acceleration exceeds 1 g: (26)

Attitude Update by Level Sensors (DIFF).
The fuel tank of a vehicle has its direction of motion set to the axis (heading), so the attitude may be updated with the rotation angle data obtained by level sensors.First, velocity and rotation angle DIFF obtained from the IMU are brought in.Then, from the rotation matrix    , the yaw angle of the vehicle can be computed as follows: Here,   is the matrix "" of the rotation matrix    and the element of column "." The matrix indicating the attitude of the fuel tank from the angle and direction of vehicle Ψ  measured by level sensors is revised as follows: Furthermore, Eulerian angles are revised as follows: The attitude measurement model function ℎ( * ) of a fuel tank still utilizes the Eulerian angles as follows: The attitude measurement value   is computed as follows, utilizing Ψ: Here,  Ψ is a 3 × 3 element of   : Here,  = diag( 2 0 ,  2  ,  2  ) is computed as follows: ,  = 10 −6 ,  level = 10. (34)

Attitude Measurement of Vehicles by Gravity.
Since the acceleration of gravity is always toward the center of the earth,  = (0, 0, −9.81 m/s 2 ) is measured if another force is not applied to the acceleration sensor.The sensor attitude can be revised by comparing the acceleration measured by the acceleration sensor and the acceleration of gravity.However, this is possible if gravity is the only force applied to the acceleration sensor.Since it is not possible to separate the acceleration of gravity from other forces and measure it, the simplest method is to assess if the condition of the applied forces other than the acceleration of gravity is ‖‖ ≈ ‖‖. = [      ], the acceleration measured by the acceleration sensor, includes various types of acceleration in accordance with the acceleration of gravity and the acceleration of the sensor.As an equation, this would be as follows:

Capacity Measuring Algorithm through Revised Rotation Angle Data and Measurement Value of Level Pressure
Capacity is computed as a digital value through the ADC board.When the acquired ADC and capacity are = Sensor adc /(ℓ/adc), the equation to compute the capacity is as follows.

System Model.
The difference value of the level sensor measured in the coordinate system of the fuel tank is converted into a navigation system: The formula above is made into an equation of state as follows: State variables are    (39)

Experiment
Table 2 shows the specifications of the pressure sensor applied to the system.The display device of the algorithm built on this research has a level sensor (Figure 7) to measure absolute pressure and the tilting angle; an IMU (Figure 8) for measuring the rotation angle, acceleration, and so forth; a programmable AMP (Figure 9) that can consistently program the increase rate of pressure between sensors by amplifying the values of pressure sensors in the desired offset; and a test bed (Figure 10).

Capacity Measurement of Distilled Water Using Level Sensors.
In order to measure capacity using level sensors, the sensor output must have a linear output in accordance with the capacity.To assess it, the research team sequentially inserted 1 L-11 L of water (distilled water, density ≈ 1) into the test bed and assessed each sensor value.In order to avoid continuity   and obtain objective data, the experiment was conducted once a day over three days.After the experiment, the team removed all distilled water from the test bed and conducted the experiment again.
Figure 11 demonstrates a 1 L/75 ADC value and 0.25 V/ 75 ADC value.Figure 11 indicates the difference of values in the first, second, and third measurements and each sensor value with a margin of error of ±1%.The output graph of Sensors numbers 5 and 6 differs from the output graph of Sensors numbers 1, 2, 3, and 4 owing to the installed location.Sensor number 7 reacts when it is over 11 L because it was located on the upper part of the test bed to measure density.
The level sensors used in this research are linear overall, demonstrating that capacity and density computation is possible.However, the margin of error increases with rotation.Figure 12 shows the output value measured three times at a 0 ∘ , 15 ∘ , 30 ∘ , and 45 ∘ rotation, and it indicates that the value increased or decreased significantly in accordance with the rotation angle.When tilting toward Sensor 3, the value of Sensor 3 and Sensor 5 on the Sensor 3-axis increases, whereas the value of Sensor 4 and Sensor 6 on the opposite side decreases.When tilting toward Sensor 4, the value of Sensor 4 and Sensor 6 on the Sensor 4-axis increases, whereas the value of Sensor 0 and Sensor 2 on the opposite side decreases.Sensor 1 and Sensor 2, which are irrelevant to the tilting direction, demonstrated no change, and so the team has not included it here.

Revision of rotation angle and acceleration of gravity
The experiment results indicate that measuring capacity and density only with a level sensor generates a large margin of error when there is tilting involved.Therefore, it is inadequate to measure the capacity and density of a liquid using only level sensor.
When revising the rotation angle and acceleration of gravity with the algorithm this paper proposes, the measurement results are obtained in Figure 13.
The research team inserted 11 L of distilled water into the test bed, tilted it by 5 ∘ , and measured the sensor value.Figure 13 demonstrates a stable result compared with the sensor graph before the application of the algorithm.This The number of repetitions is because the accurate estimation of the rotation angle and its revision value were designed with a Kalman filter.Thus estimating the rotation angle is important in this algorithm.
The comparative results of the proposed algorithm that estimates the rotation angle are as shown in Figure 14.
Figure 14 demonstrates the measurements of rotation angles using only the IMU; the measurement is carried out once every minute.As indicated in Figure 14, the margin of error continuously accumulated, demonstrating a large difference in value after 17 min from the beginning.Furthermore, when measuring the rotation angle with only level sensors (Figure 15), the margin of error does not accumulate and is drifting as with the IMU data, but the margin of error is significantly large.
Figure 16 shows the rotation angle using the algorithm proposed by this paper.The margin of error is within ±5%, which is relatively accurate compared with existing methods.The results of measuring capacity and density using the algorithm this paper has proposed are as shown in Figure 17.
Figure 17 shows the output values of the experiment, where 1 L was inserted as a sensor value without tilting.However, the measurements were not conducted sequentially because each sensor has a slightly different minimum pressure value that computes the sensor output and because the sensor locations differ as well.For precise measurements, the research team inserted 3 L of distilled water, and when all sensors were reacting sufficiently, the team initialized the value and measured again.As demonstrated in Figure 17, accurate results from the measuring capacity could be obtained because the sensor value was close to the primary equation and each sensor had less difference in value.
Figure 18 is the computation of density and capacity, which is very precise because it has a margin of error of only ±0.05 L compared with the actual capacity.The density was also very precise, with a margin of error of ±0.02 kg/m 3 .
Sensor number 4 and Sensor number 7 are used to measure the density.Sensor number 7 is placed on the upper part so it reacts only when the solution is more than 10 liters.Therefore the density can be measured from 10 liters.But if the sensor can be placed lower, we can measure the density by using less fluid.
In addition, Figure 19 demonstrates the results of the measurements with the same method as above, but tilting by 5 ∘ .
As indicated in Figure 19, the margin of error of measurement increases in accordance with the tilting, and the maximum margin of error was only 0.5 L. Since it can be presumed that the larger the target capacity for measurement, the lesser the error ratio, this can be said to be a precise algorithm.
As demonstrated in Figure 20, compared with the results when there was no tilting, there has rarely been any difference in density, and the result was ±0.05 kg/m 3 .
For Figure 21, the program was implemented with a graphical user interface (GUI), and its resulting values were continuously obtained and observed.

Conclusion
This research proposed an algorithm using a level sensor and IMU sensor to improve the limitations of the existing methods of measuring the liquid capacity.For this, the research    team installed level sensors and IMU sensors in a liquid tank and conducted sensing in real time.The proposed algorithm utilized the strengths of each sensor and supplemented their weaknesses, which were assessed through experiments.
A measurement evaluation of the liquid level sensor was conducted by increasing and decreasing the liquid in the tank.The team increased the water level by units of 1 L and analyzed the output characteristics of up to 12 L.As a result, the system demonstrated a highly linear precision that was within ±0.5%.The offset voltage of the liquid level sensor was 93.5 mV, with a sensitivity of 363.5 mV/L.The output characteristic in inclined environments resulted in a value of 1% FS.In addition, the sensitivity to density measurements was 1085.28 [kg/m 3 ], which is consistent with the computed value of the sensor output.The reproducibility of the sensor showed minute hysteresis, but since there was repeatability, a high-precision pressure level sensor was implemented as it can ignore this repeatability.
The research team believes that, with the proposed algorithm, precise liquid capacity measurements in real time would be possible even during movement or tilting of the liquid tank.

Figure 1 :
Figure 1: Disadvantages of current fuel-measuring methods for vehicles.(a) Structural limits result in inability to measure the remaining 10% of fuel.(b) Fuel-measuring error on slopes.(c) Fuel-measuring error when accelerating.

Figure 2 :
Figure 2: Algorithm applied in this assignment.

Figure 4 :
Figure 4: Sensor location to measure the sensor.

Figure 5 :
Figure 5: Liquid tank and sensor mimetic diagram to measure density.

Figure 5
is the schematic diagram of the attached sensor in the container in order to compute the density.The density is computed with the signals outputted from the sensor attached to the container: Sensitivity of the sensor × Pressure = Output voltage → Output Voltage Sensitivity of the sensor = Pressure (pa) ,

Figure 6 :
Figure 6: Measuring methods of the rotation angle of the liquid tank and acting forces. Filter

4. 2 . 1 .
Initialize.The state variables to be used in the Kalman filter are indicated as a rotation matrix and as Eulerian angles, and the covariance matrix is indicated as the covariance of Eulerian angles.State variables are as follows:  =    (State variables for the rotation matrix)  = Λ  (State variables for the Eulerian angles) .

Figure 13 :
Figure 13: Measurement result graph using algorithm that has revised rotation angle and acceleration of gravity.

Figure 14 :
Figure 14: Result graph of measuring the rotation angle using IMU data.

Figure 15 :
Figure 15: Result graph of measuring the rotation angle using level sensor data.

Figure 16 :
Figure 16: Result graph of measuring the rotation angle with the proposed algorithm.

Figure 17 :Figure 18 :
Figure 17: Measurement graph of capacity and density using the proposed algorithm.

Figure 19 :
Figure 19: Comparative graph 2 to assess the precision of measurements of capacity and density.

Figure 20 :
Figure 20: Comparative graph 3 to assess the precision of measurements of capacity and density.

Figure 21 :
Figure 21: GUI to observe the measurement values in real time.

Table 1 :
Comparison of sensors used for measuring liquid level.