Impact of Sloshing on Fossil Fuel Loss during Transport

: This study attempts to uncover the most common issue of fuel shortage faced by the oil and transportation industry worldwide. In Pakistan, petroleum is transported to the northern areas from the south coast. Currently, this is done using road tankers as a pipeline is still under construction. However, even after the pipeline becomes operative, road tankers would still be used for intra-city transport. Findings from this study can be used to determine the inter-city transport losses faced by oil companies. This study determines the hydrocarbons lost to the environment during inter-city road transport of petroleum. It takes nearly 2–3 days to complete a one-way trip with the fully loaded tank. Much work has been reported worldwide on hydrocarbon emissions, but nearly all of it has been done either for storage tanks / vessels or fuel tanks in rails / cars. The aim of this study was to investigate the actual amount of fuel lost to the environment due to the sloshing of liquid. Also, the results were expected to help in determining the extent of hazardous emissions resulting from road transport of petroleum. Hence, measures could be taken by the concerned authorities to mitigate the emissions if they exceeded the acceptable range. The sloshing was not found to contribute much in terms of vapor loss. Valve location was found to be important as no loss was obtained from the third compartment because the valve is to the right in this chamber. A negligible amount of fuel was lost from the ﬁrst and second compartments per application of the brakes. Over a whole trip of 2–3 days, if the tanker braked 500 times, a total of 9–10 L would be lost to the environment.


Introduction
This study explicitly calculated the amount of fuel lost from a 48,000 L tank on a one-way trip from Karachi to Sheikhupura (Machike). The tank was comprised of three individual compartments of 16,000 L. Each compartment was fitted with a pressure/vacuum breather vent valve at the top of the tank to prevent the tank from creating overpressure/vacuum.
Besides the evaporation losses, the movement of the liquid inside the tank when it is subjected to harsh acceleration or brakes can also contribute to fuel loss. This is because the pressure exerted by the liquid on the tank walls can be substantial enough to reach the opening pressures of the P/V vent valve of the tank. The movement of the fluid after an impact produces large shifts of weight which build momentum and do not settle down easily. This movement of the liquid is termed 'sloshing'.
One horizontal and two vertical baffles Out of these, the most reliable i.e., the one with least slosh was with three vertical baffles [7]. Joshi et al., compared the average and maximum impact pressures of different fill levels using CFD. This was done for a chamfered LNG road tank. The dimensions of the chamfered tank were 1.3 m × 0.9 m with a 45 degree chamfer, 30 cm offset chamfer. The simulation was performed with and without baffles for each fill-level. Five different monitor points were selected each on the vertices of the tank to find the maximum impact locations. Fill levels of 30, 50, and 80 were used. Baffles seemed to mitigate the impact pressures at all points and fill levels. However, for low filling levels, baffles seemed to mitigate the impact pressures nearly by a factor of 10 [8].
Močilan et al. performed a CFD analysis on a fuel tank using ANSYS FLUENT. The proportion of fuel to air was 60:40 and time-varying acceleration was applied. The free surface was modeled. Also, the maximum stress in the tank was found. The research found that reduced amounts of stress levels were found with baffles as opposed to a tank without baffles [9].
In another study, Siddapureddy and Satapathy performed CFD on a rectangular tank with dimensions 5 × 4 × 4 m 3 . The software used was FLUENT. The tank was partially filled with water. The tank started from rest. In the first two seconds, it rose to 30 km/h and remained at 30 km/h for another 4 s. Then in the last two seconds, the velocity decreased from 30 km/h to 0 km/h. In the study, three different fill levels of 30%, 50%, and 70% were analyzed. The exercise of different fill levels was repeated for tanks without baffles, with horizontal baffles, with vertical baffles, and with multiple baffles. Also, the simulation results were validated by experiments. It also introduced a new concept of percentage reduction (PR) which was calculated using mean deviation (MD). MD is the mean of the deviations from the mean position i.e., the position of the free surface at the start of the simulation. Hence, PR can be used to assess the performance of different baffle designs. Results found that vertical baffles were more efficient in contrast to horizontal ones. However, multiple baffles showed better results than horizontal/vertical ones where minimization of the sloshing effect was more crucial than the complexity of the design of the tank [10].

Valve Model
The tank is fitted with a Breather Vent Valve in each compartment to prevent overpressure/vacuum [11,12]. The tank is composed of three compartments each of 16,000 L capacity. Hence, three P/V Vent valves are present in total. The threshold absolute pressures of the valve are in compliance with EN 14595 [13,14] as shown in Figure 1: The valve is also fitted with a spill over protection which prevents the spillage of fuel when a rollover occurs.
The flow rate of the valve was modeled. The cracking pressure of the valve is 111,000 Pa. Hence, the model was divided into three regions:

1.
Valve closed: This is when the pressure is below 111,000 Pa. Valve partially open: This is when pressure is between 111,000 and 131,323 Pa. Two data points were provided by the manufacturer [15] which are shown in Table 1 below: These data points were used to form the equation of this region:

3.
Valve fully open: This is when the pressure is equal to or greater than 131,323 Pa. Hence, Q = 132,000 L/h The flow rates of each region are summarized in Table 2. The below graph ( Figure 2) shows the behavior of the valve in all three regions. A curve fit was done to determine the equation of region 2, however, flow rates of region 1 and region 3 were directly supplied by the manufacturer.
The flow rate of the valve was modeled. The cracking pressure of the valve is 111,000 Pa. Hence, the model was divided into three regions: 1. Valve closed: This is when the pressure is below 111,000 Pa. Q = 0  [15].

Computational Fluid Dynamics
Sloshing is the movement of the free surface when the tank is subjected to acceleration/deceleration. The CFD method solves Navier Stokes equations of momentum and heat and mass transfer. This was performed in ANSYS CFX. CFX uses the finite volume technique to solve these equations. In this technique, the area of interest is divided into smaller regions called control volumes. The equations are solved numerically for each control volume [16].
A multiphase sloshing simulation was set up to determine the hydrodynamic pressures generated. These pressures were then used to estimate vapor loss.

Sensitivity Analysis
A sensitivity analysis was performed on the elliptical tank to determine the spatial and temporal coordinates where the desired convergence was achieved. A general Gaussian waveform of acceleration with a peak of 2.6 m/s 2 i.e., 0.265 g was applied as shown in Figure 3. An initial velocity of 0.2 m/s was applied. The period of the wave was kept at 5 s. Peak acceleration of 2.6 m/s 2 was applied because it was obtained from the tracker report of the truck. The simulation was run for 1 s in each case. The mesh was generated using a curvature size function. The curvature size function was used because of the curvature in the geometry. The function examines curvature on edges and faces and forms elements without exceeding the maximum element size and the curvature normal angle specified. The CFX solver calculates the curvature normal angle, growth rate, and minimum edge

Computational Fluid Dynamics
Sloshing is the movement of the free surface when the tank is subjected to acceleration/deceleration. The CFD method solves Navier Stokes equations of momentum and heat and mass transfer. This was performed in ANSYS CFX. CFX uses the finite volume technique to solve these equations. In this technique, the area of interest is divided into smaller regions called control volumes. The equations are solved numerically for each control volume [16].
A multiphase sloshing simulation was set up to determine the hydrodynamic pressures generated. These pressures were then used to estimate vapor loss.

Sensitivity Analysis
A sensitivity analysis was performed on the elliptical tank to determine the spatial and temporal coordinates where the desired convergence was achieved. A general Gaussian waveform of acceleration with a peak of 2.6 m/s 2 i.e., 0.265 g was applied as shown in Figure 3. An initial velocity of 0.2 m/s was applied. The period of the wave was kept at 5 s. Peak acceleration of 2.6 m/s 2 was applied because it was obtained from the tracker report of the truck. The simulation was run for 1 s in each case.

Computational Fluid Dynamics
Sloshing is the movement of the free surface when the tank is subjected to acceleration/deceleration. The CFD method solves Navier Stokes equations of momentum and heat and mass transfer. This was performed in ANSYS CFX. CFX uses the finite volume technique to solve these equations. In this technique, the area of interest is divided into smaller regions called control volumes. The equations are solved numerically for each control volume [16].
A multiphase sloshing simulation was set up to determine the hydrodynamic pressures generated. These pressures were then used to estimate vapor loss.

Sensitivity Analysis
A sensitivity analysis was performed on the elliptical tank to determine the spatial and temporal coordinates where the desired convergence was achieved. A general Gaussian waveform of acceleration with a peak of 2.6 m/s 2 i.e., 0.265 g was applied as shown in Figure 3. An initial velocity of 0.2 m/s was applied. The period of the wave was kept at 5 s. Peak acceleration of 2.6 m/s 2 was applied because it was obtained from the tracker report of the truck. The simulation was run for 1 s in each case. The mesh was generated using a curvature size function. The curvature size function was used because of the curvature in the geometry. The function examines curvature on edges and faces and forms elements without exceeding the maximum element size and the curvature normal angle specified. The CFX solver calculates the curvature normal angle, growth rate, and minimum edge The mesh was generated using a curvature size function. The curvature size function was used because of the curvature in the geometry. The function examines curvature on edges and faces and forms elements without exceeding the maximum element size and the curvature normal angle specified. The CFX solver calculates the curvature normal angle, growth rate, and minimum edge length based on the geometry. The remaining parameters (minimum face size, maximum face size, and maximum size) were decided from the results of the sensitivity analysis. Initially, the following values were used as shown in Table 3: A length of 0.04 m was selected as the mesh size and 0.1 s as the time step size. The detailed results from the sensitivity analysis are shown in the results section.

Basic Set-Up
A multiphase simulation was setup in ANSYS CFX. For simplicity, one compartment of the tank was simulated. The geometry was imported from Solid Works. The multiphase homogeneous flow was modeled and gravitational acceleration was applied. A mesh sensitivity analysis was performed to determine the coordinates where convergence was achieved. These coordinates were then used in further simulations to analyze different scenarios. The velocities and accelerations against each time step for each case were recorded in an excel file. The excel file was then given as input to ANSYS CFX. Properties of air and octane [17] were used as shown in Tables 4 and 5.
Air:  Pressures are measured at points beneath the vent valve. The location of the valve is slightly different for each compartment. Hence, pressures are monitored at all three locations of the vent valve. However, this is done on a single compartment. The following figure is of a 48,000 L tank with three chambers (Figure 4). The tank is comprised of three chambers each of 16,000 L capacity.  Three monitor points were selected ( Figure 5):

Acceleration/Deceleration Function
Different cases were analyzed using various acceleration/deceleration functions. The acceleration/deceleration behavior of trucks and a truck-trailer combination was studied. This was done to develop realistic acceleration/deceleration functions to come up with accurate vapor loss generated through sloshing.
The Verizon Telematics [19] describes the minimum and maximum thresholds for hard braking and hard acceleration for light-and heavy-duty vehicles. For light-duty vehicles, it gives a minimum and maximum acceleration threshold of 0.36 g and 0.8 g respectively. For heavy-duty vehicles, it gives a minimum and maximum braking threshold of 0.22 g and 1 g respectively.
The Highway Safety Research Institute described the braking performance of trucks and tractor-trailer combinations by performing several tests [20]. It describes a test named 'brake effectiveness test'. In one such test, the tractor-trailer was moving at 30 mph and the brakes were Three monitor points were selected ( Figure 5):

1.
This point was selected to be directly below the valve of the third compartment.

Acceleration/Deceleration Function
Different cases were analyzed using various acceleration/deceleration functions. The acceleration/deceleration behavior of trucks and a truck-trailer combination was studied. This was done to develop realistic acceleration/deceleration functions to come up with accurate vapor loss generated through sloshing.
The Verizon Telematics [19] describes the minimum and maximum thresholds for hard braking and hard acceleration for light-and heavy-duty vehicles. For light-duty vehicles, it gives a minimum and maximum acceleration threshold of 0.36 g and 0.8 g respectively. For heavy-duty vehicles, it gives a minimum and maximum braking threshold of 0.22 g and 1 g respectively.
The Highway Safety Research Institute described the braking performance of trucks and tractor-trailer combinations by performing several tests [20]. It describes a test named 'brake effectiveness test'. In one such test, the tractor-trailer was moving at 30 mph and the brakes were

Acceleration/Deceleration Function
Different cases were analyzed using various acceleration/deceleration functions. The acceleration/deceleration behavior of trucks and a truck-trailer combination was studied. This was done to develop realistic acceleration/deceleration functions to come up with accurate vapor loss generated through sloshing.
The Verizon Telematics [19] describes the minimum and maximum thresholds for hard braking and hard acceleration for light-and heavy-duty vehicles. For light-duty vehicles, it gives a minimum and maximum acceleration threshold of 0.36 g and 0.8 g respectively. For heavy-duty vehicles, it gives a minimum and maximum braking threshold of 0.22 g and 1 g respectively.
The Highway Safety Research Institute described the braking performance of trucks and tractor-trailer combinations by performing several tests [20]. It describes a test named 'brake effectiveness test'. In one such test, the tractor-trailer was moving at 30 mph and the brakes were applied until the vehicle stopped. The deceleration was found to be −0.38 g and the stopping distance was 100 feet. The values from this test were used to form the normal acceleration/deceleration function because the deceleration and the cruising speed of our tractor-trailer are in close comparison to the parameters in this test.
The maximum values of deceleration from a test just before jack-knifing occurred were used to develop the abnormal case of deceleration/acceleration. The initial velocity was 30 mph and the maximum deceleration was −0.55 g with a stopping distance of 77 feet. The values from this test were used to form the abnormal acceleration/deceleration function because no wheel lock occurred at this value. The input functions were derived considering all the above information. Four cases were analyzed: 1.
Acceleration-normal case: The values (acceleration, distance covered) were kept similar to the values of one of the tests of the brake effectiveness test [20]. applied until the vehicle stopped. The deceleration was found to be −0.38 g and the stopping distance was 100 feet. The values from this test were used to form the normal acceleration/deceleration function because the deceleration and the cruising speed of our tractor-trailer are in close comparison to the parameters in this test.
The maximum values of deceleration from a test just before jack-knifing occurred were used to develop the abnormal case of deceleration/acceleration. The initial velocity was 30 mph and the maximum deceleration was −0.55 g with a stopping distance of 77 feet. The values from this test were used to form the abnormal acceleration/deceleration function because no wheel lock occurred at this value. The input functions were derived considering all the above information. Four cases were analyzed: 1. Acceleration-normal case: The values (acceleration, distance covered) were kept similar to the values of one of the tests of the brake effectiveness test [20].   Energies 2020, 13, x FOR PEER REVIEW 8 of 32 applied until the vehicle stopped. The deceleration was found to be −0.38 g and the stopping distance was 100 feet. The values from this test were used to form the normal acceleration/deceleration function because the deceleration and the cruising speed of our tractor-trailer are in close comparison to the parameters in this test.
The maximum values of deceleration from a test just before jack-knifing occurred were used to develop the abnormal case of deceleration/acceleration. The initial velocity was 30 mph and the maximum deceleration was −0.55 g with a stopping distance of 77 feet. The values from this test were used to form the abnormal acceleration/deceleration function because no wheel lock occurred at this value. The input functions were derived considering all the above information. Four cases were analyzed: 1. Acceleration-normal case: The values (acceleration, distance covered) were kept similar to the values of one of the tests of the brake effectiveness test [20].   2. Acceleration-abnormal case: The values (acceleration, distance covered) were kept similar to the maximum values of one of the tests of the brake effectiveness test [20]. The vehicle started with an initial velocity of 0.2 m/s and rose to 13.24 m/s in 3.9 s after covering a distance of 27.9 m. The maximum acceleration was kept at 0.55 g. The acceleration and velocity graphs are shown below (Figures 8 and 9):   [20]. The vehicle started with an initial velocity of 0.2 m/s and rose to 13.24 m/s in 3.9 s after covering a distance of 27.9 m. The maximum acceleration was kept at 0.55 g. The acceleration and velocity graphs are shown below (Figures 8 and 9): Energies 2020, 13, x FOR PEER REVIEW 9 of 32

3.
Deceleration-normal case: The values (acceleration, distance covered) were kept similar to the values of one of the tests of the brake effectiveness test [18]. The vehicle started from 13.5 m/s and came to rest in 5 s after covering a distance of 31 m. The maximum acceleration was −0.4 g. The velocity and acceleration curves are shown below (Figures 10 and 11): 3. Deceleration-normal case: The values (acceleration, distance covered) were kept similar to the values of one of the tests of the brake effectiveness test [18]. The vehicle started from 13.5 m/s and came to rest in 5 s after covering a distance of 31 m. The maximum acceleration was −0.4 g. The velocity and acceleration curves are shown below (Figures 10 and 11):        The lowest mode slosh frequencies were calculated using the following formula [21]:

Slosh Calculation Validation
Next, the frequency and time period were calculated: Here, a is the length of the tank, and h is the height of fluid in the tank. Hence, four different frequencies were calculated for each fluid height. Corresponding periods were calculated for each height. The sloshing simulation was set up. Exactly the same rectangular tank was used as in the theoretical calculation. A gentle Gaussian acceleration pulse of 0.1 g was applied for 1 s as shown in Figure 14 and the simulation was run for 10 s. It was run for a long enough duration so as to determine adequate wave periods to establish confidence in the simulation results.
From sloshing simulation, periods of the following three parameters were compared with the theoretical calculation for each fluid height:

1.
The pressure at a fixed monitor point at the top of the tank just beneath the valve 2.
Maximum height of the free surface 3.
Minimum height of the free surface The periods of these three parameters were compared with the calculated period. In the case of pressures, the period was measured conventionally i.e., the difference between two consecutive peaks/troughs was taken because it was measured at a fixed point. However, in the case of heights, the period completes between three consecutive peaks/troughs as it was a half-wave i.e., mode 1 [22] as shown in Figure 15. This was because heights were not measured at a fixed point, rather, the software showed the maximum/minimum height of the free surface.
Here, a is the length of the tank, and h is the height of fluid in the tank. Hence, four different frequencies were calculated for each fluid height. Corresponding periods were calculated for each height.
The sloshing simulation was set up. Exactly the same rectangular tank was used as in the theoretical calculation. A gentle Gaussian acceleration pulse of 0.1 g was applied for 1 s as shown in Figure 14 and the simulation was run for 10 s. It was run for a long enough duration so as to determine adequate wave periods to establish confidence in the simulation results. The periods of these three parameters were compared with the calculated period. In the case of pressures, the period was measured conventionally i.e., the difference between two consecutive peaks/troughs was taken because it was measured at a fixed point. However, in the case of heights, the period completes between three consecutive peaks/troughs as it was a half-wave i.e., mode 1 [22] as shown in Figure 15. This was because heights were not measured at a fixed point, rather, the software showed the maximum/minimum height of the free surface. The numerical results were found to be in close comparison to the theoretical calculation in a ee cases-pressure, maximum, and minimum height (Figures 16-18). The simulation an lculation were performed for four different fluid heights of 0.25 m, 0.5 m, 1.5 m, and 1.75 m. I ures 16-18, the calculated periods at each fluid height almost overlie the periods obtained from simulation. These results were used to build confidence in the slosh simulation setup. The sam up was then used on the cases to be analyzed. Screenshots from the simulation are shown i pendix A.  The numerical results were found to be in close comparison to the theoretical calculation in all three cases-pressure, maximum, and minimum height (Figures 16-18). The simulation and calculation were performed for four different fluid heights of 0.25 m, 0.5 m, 1.5 m, and 1.75 m. In Figures 16-18, the calculated periods at each fluid height almost overlie the periods obtained from the simulation. These results were used to build confidence in the slosh simulation setup. The same setup was then used on the cases to be analyzed. Screenshots from the simulation are shown in Appendix A.
three cases-pressure, maximum, and minimum height (Figures 16-18). The simulation and calculation were performed for four different fluid heights of 0.25 m, 0.5 m, 1.5 m, and 1.75 m. In Figures 16-18, the calculated periods at each fluid height almost overlie the periods obtained from the simulation. These results were used to build confidence in the slosh simulation setup. The same setup was then used on the cases to be analyzed. Screenshots from the simulation are shown in Appendix A.

Sensitivity Analysis
Before analyzing the actual cases in the actual tank, a sensitivity analysis was performed to gain confidence in the mesh size being used along with an acceptable time step. Screenshots from the simulation can be seen in Appendix A.
Spatial: The spatial coordinate was varied by changing the maximum element size in the curvature size function while keeping everything else constant. These were varied between 0.1 m

Sensitivity Analysis
Before analyzing the actual cases in the actual tank, a sensitivity analysis was performed to gain confidence in the mesh size being used along with an acceptable time step. Screenshots from the simulation can be seen in Appendix A.

Sensitivity Analysis
Before analyzing the actual cases in the actual tank, a sensitivity analysis was performed to gain confidence in the mesh size being used along with an acceptable time step. Screenshots from the simulation can be seen in Appendix A.
Spatial: The spatial coordinate was varied by changing the maximum element size in the curvature size function while keeping everything else constant. These were varied between 0.1 m and 0.03 m. The desired convergence was achieved at 0.04 m as shown in Table 6. The maximum pressure generated for mesh sizes of 0.04 m, 0.035 m, and 0.03 m was 1.0136 × 10 5 Pa. Hence, 0.04 m was selected as the mesh size. The maximum generated pressures inside the tank were recorded for each mesh size for 1 s with a time step size of 0.1 s in each case. Pressures generated at 0.04 m, 0.035 m, and 0.03 m almost showed the same trend. Therefore, 0.04 m was selected as the mesh size. The results are shown in Figure 19 below: The maximum generated pressures inside the tank were recorded for each mesh size for 1 s with a time step size of 0.1 s in each case. Pressures generated at 0.04 m, 0.035 m, and 0.03 m almost showed the same trend. Therefore, 0.04 m was selected as the mesh size. The results are shown in Figure 19 below: Temporal: The time step was varied between 0.05 s and 0.25 s. The desired convergence was achieved at higher step sizes of 0.2 s and 0.25 s as well, but 0.1 s was selected so as not to lose too much detail. However, 0.05 s was not selected because of time constraints. The maximum generated pressures vs. time step size are summarized in Table 7.   Temporal: The time step was varied between 0.05 s and 0.25 s. The desired convergence was achieved at higher step sizes of 0.2 s and 0.25 s as well, but 0.1 s was selected so as not to lose too much detail. However, 0.05 s was not selected because of time constraints. The maximum generated pressures vs. time step size are summarized in Table 7. The maximum pressure was recorded for each time step size for the duration of 1 s. The pressure values with a step size of 0.1 s almost follow the same trend as the values with a step size of 0.05 s. Hence, 0.1 s was selected as the step size due to time and resource constraints. The results are shown in Figure 20.

Slosh
Hydrodynamic pressures were recorded. In the case of deceleration, much higher pressures were produced because the fluid hits the valve with greater pressure compared to the acceleration scenario. Also, the valve opening pressures were reached only for the first and second compartments. This is because the valve is located towards the front of the tank in these compartments.
Case 1: Normal Acceleration The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 1 are shown in Figure 21. The opening pressure of the valve i.e., 111 kPa was not reached at any monitor point. The least pressures were produced at Monitor Point 1 i.e., third compartment. The reason was its location i.e., it was placed farthest from the left and the tanker was moving towards the left. However, a similar trend of pressures was obtained for Monitor Points 2 and 3.

Slosh
Hydrodynamic pressures were recorded. In the case of deceleration, much higher pressures were produced because the fluid hits the valve with greater pressure compared to the acceleration scenario. Also, the valve opening pressures were reached only for the first and second compartments. This is because the valve is located towards the front of the tank in these compartments.
Case 1: Normal Acceleration The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 1 are shown in Figure 21. The opening pressure of the valve i.e., 111 kPa was not reached at any monitor point. The least pressures were produced at Monitor Point 1 i.e., third compartment. The reason was its location i.e., it was placed farthest from the left and the tanker was moving towards the left. However, a similar trend of pressures was obtained for Monitor Points 2 and 3.

compartments.
Case 1: Normal Acceleration The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 1 are shown in Figure 21. The opening pressure of the valve i.e., 111 kPa was not reached at any monitor point. The least pressures were produced at Monitor Point 1 i.e., third compartment. The reason was its location i.e., it was placed farthest from the left and the tanker was moving towards the left. However, a similar trend of pressures was obtained for Monitor Points 2 and 3.

Case 2: Abnormal Acceleration
The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 2 are shown in Figure 22. In this case, as well, the opening pressure of the valve did not reach any monitor point. Similar to case 1, the least pressures were generated in the third compartment. Case 2: Abnormal Acceleration The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 2 are shown in Figure 22. In this case, as well, the opening pressure of the valve did not reach any monitor point. Similar to case 1, the least pressures were generated in the third compartment.

Case 3: Normal Deceleration
The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 3 are shown in Figure 23. Opening pressures were exceeded for the first 0.2 s at point 2 and point 3. Also, an air pocket was found beneath the valve when the valve opens i.e., for the first 0.2 s (Appendix A, Figure A31). This confirms that the air/vapor mixture leaked out during that interval.

Case 3: Normal Deceleration
The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 3 are shown in Figure 23. Opening pressures were exceeded for the first 0.2 s at point 2 and point 3. Also, an air pocket was found beneath the valve when the valve opens i.e., for the first 0.2 s (Appendix A, Figure A31). This confirms that the air/vapor mixture leaked out during that interval.
Case 4: Abnormal Deceleration The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 4 are shown in Figure 24. Similar to case 3, opening pressures were exceeded for the first 0.2 s at point 2 and point 3. However, the pressures reached in Case 4 were slightly higher than Case 3 due to the abnormal case. Also, an air pocket was found beneath the valve when the valve opened i.e., for the first 0.2 s (Appendix A, Figure A35). This confirms that the air/vapor mixture leaked out during that interval.

Case 3: Normal Deceleration
The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 3 are shown in Figure 23. Opening pressures were exceeded for the first 0.2 s at point 2 and point 3. Also, an air pocket was found beneath the valve when the valve opens i.e., for the first 0.2 s (Appendix A, Figure A31). This confirms that the air/vapor mixture leaked out during that interval.

Case 4: Abnormal Deceleration
The slosh results i.e., the pressures generated with respect to time at all three monitor points for Case 4 are shown in Figure 24. Similar to case 3, opening pressures were exceeded for the first 0.2 s at point 2 and point 3. However, the pressures reached in Case 4 were slightly higher than Case 3 due to the abnormal case. Also, an air pocket was found beneath the valve when the valve opened i.e., for the first 0.2 s (Appendix A, Figure A35). This confirms that the air/vapor mixture leaked out during that interval. The slosh results showed that the valve of the last compartment did not open even in abnormal cases of acceleration and deceleration. However, the valve of the first and second compartments did open in both cases of deceleration. However, the duration in which the valve opened was less than 0.2 s. Screenshots from the simulation can be seen in Appendix A.
The pressure obtained from sloshing in the normal and abnormal cases of deceleration in compartment 1 and 2 exceeded the cracking pressure of the valve. At 0.1 s, the pressure was much higher than the cracking pressure. At 0.2 s, the pressure dropped down below atmospheric pressure. It rose to atmospheric pressure at 0.4 s and remained nearly at atmospheric pressure for the rest of the simulation time. Hence, to find the loss from the valve in the first 0.2 s, three data points were selected. One was atmospheric pressure at 0 s. The second and third were the pressures obtained from sloshing at 0.1 s and 0.2 s respectively. A curve fit was made to fit all three data points ( Figure  25 The slosh results showed that the valve of the last compartment did not open even in abnormal cases of acceleration and deceleration. However, the valve of the first and second compartments did open in both cases of deceleration. However, the duration in which the valve opened was less than 0.2 s. Screenshots from the simulation can be seen in Appendix A. The pressure obtained from sloshing in the normal and abnormal cases of deceleration in compartment 1 and 2 exceeded the cracking pressure of the valve. At 0.1 s, the pressure was much higher than the cracking pressure. At 0.2 s, the pressure dropped down below atmospheric pressure. It rose to atmospheric pressure at 0.4 s and remained nearly at atmospheric pressure for the rest of the simulation time. Hence, to find the loss from the valve in the first 0.2 s, three data points were selected. One was atmospheric pressure at 0 s. The second and third were the pressures obtained from sloshing at 0.1 s and 0.2 s respectively. A curve fit was made to fit all three data points (Figure 25). Two more lines were plotted: Y = 111,000 Y = 131,323 the simulation time. Hence, to find the loss from the valve in the first 0.2 s, three data points were selected. One was atmospheric pressure at 0 s. The second and third were the pressures obtained from sloshing at 0.1 s and 0.2 s respectively. A curve fit was made to fit all three data points ( Figure  25). Two more lines were plotted: Both these lines intersect the curve at two points. Hence, all four points were found using simultaneous equations. To do so, the quadratic equation was set equal to 111,000 and 131,323. x1 and x4 were obtained by finding the roots of the equation when the right-hand side of the equation was set equal to 111,000. Similarly, x2 and x3 were obtained by finding the roots of the equation when the right-hand side of the equation was set equal to 131,323. The quadratic equations and its roots for points 2 and 3 in the deceleration cases are summarized in Table 8. Both these lines intersect the curve at two points. Hence, all four points were found using simultaneous equations. To do so, the quadratic equation was set equal to 111,000 and 131,323. x 1 and x 4 were obtained by finding the roots of the equation when the right-hand side of the equation was set equal to 111,000. Similarly, x 2 and x 3 were obtained by finding the roots of the equation when the right-hand side of the equation was set equal to 131,323. The quadratic equations and its roots for points 2 and 3 in the deceleration cases are summarized in Table 8.  As shown in Figure 26, the loss from x1 (0.0053 s) to x2 (0.018 s) and from x3 (0.152 s) to x4 (0.165 s) was calculated using the equation of the flow rate of region 2. Additionally, from x2 (0.018 s) to x3 (0.152 s) was calculated using the equation of the flow rate of region 3.
The loss obtained from all cases is summarized in Table 9   From x 2 to x 3 , the flow rate was found using the following value of the flow rate i.e., region 3: As shown in Figure 26, the loss from x 1 (0.0053 s) to x 2 (0.018 s) and from x 3 (0.152 s) to x 4 (0.165 s) was calculated using the equation of the flow rate of region 2. Additionally, from x 2 (0.018 s) to x 3 (0.152 s) was calculated using the equation of the flow rate of region 3.
The loss obtained from all cases is summarized in Table 9 below. A maximum of 0.009 L occurred in a normal case of deceleration per application of brakes from compartment 1 and a maximum of 0.0088 L occurred in a normal case of deceleration per application of brakes from compartment 2. Hence, if the tanker is assumed to brake 500 times over the whole trip, then, a loss of 9 L is obtained due to slosh. A maximum of 0.01 L also occurred in the abnormal case of deceleration per application of brakes from compartment 1 and a maximum of 0.0096 L occurred in the abnormal case of deceleration per application of brakes from compartment 2. Hence, if the tanker is assumed to brake 500 times over the whole trip, then, a loss of 10 L is obtained due to slosh.
As shown in Table 9, no loss was obtained from compartment 3 even in cases of deceleration because of the location of the P/V vent valve. The reason was that the valve was located farthest from the front of the compartment i.e., 2491 mm from the start of compartment 3.
The valves of the first and second compartments were placed at a distance of 1720 mm and 1780 mm respectively from the start of the compartment. Hence, moving the valve slightly along the longitudinal axis may change the loss results drastically.

Conclusions
The slosh model was validated by performing various checks on a rectangular tank. All of them were found to be closely comparable with the theoretical values. After validation, slosh simulation was performed on the actual tank. A few acceleration and deceleration functions were given as inputs to the solver. These cases were constructed after reviewing the literature and actual truck driving data. The pressures were then used to calculate the vapor loss due to slosh. The valve location was found to be of considerable importance as pressures vary largely by moving the valve slightly along the longitudinal axis. No loss was found in the third compartment even in abnormal cases of acceleration and deceleration. The location of the valve of the third compartment was farthest from the front of the compartment i.e., 2491 mm from the front of the compartment.
The slosh was not found to contribute much in terms of vapor loss. However, a negligible amount of fuel is lost to the environment in deceleration cases from the first and second compartments. Hence, valve location serves as an important factor in vapor loss. Even if the tanker is assumed to brake 500 times, the tanker still only loses a total of 9 L, and if the tanker brakes 500 times abnormally, it still loses only a total of 10 L over the whole trip.