Next Article in Journal
Numerical Simulations of Molten Breakup Behaviors of a de Laval-Type Nozzle, and the Effects of Atomization Parameters on Particle Size Distribution
Next Article in Special Issue
Study on the Law of Diesel Oil Carrying Water in Lanzhou–Chengdu–Chongqing Product Oil Pipeline Based on Large Eddy Simulation
Previous Article in Journal
Comparison of Dimethylformamide with Dimethylsulfoxide for Quality Improvement of Distillate Recovered from Waste Plastic Pyrolysis Oil
Previous Article in Special Issue
Indoor Environmental Quality (IEQ) Analysis of a Two-Sided Windcatcher Integrated with Anti-Short-Circuit Device for Low Wind Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Influence of Floating-Deck Height on Oil-Vapor Migration and Emission of Internal Floating-Roof Tank Based on Numerical Simulation and Wind-Tunnel Experiment

1
Jiangsu Key Laboratory of Oil & Gas Storage and Transportation Technology, Changzhou University, Changzhou 213164, China
2
Yangzhong Environmental Sanitation Administration Agency, Yangzhong 212200, China
*
Author to whom correspondence should be addressed.
Processes 2020, 8(9), 1026; https://doi.org/10.3390/pr8091026
Submission received: 23 June 2020 / Revised: 17 August 2020 / Accepted: 19 August 2020 / Published: 21 August 2020
(This article belongs to the Special Issue CFD Applications in Energy Engineering Research and Simulation)

Abstract

:
Internal floating-roof tanks (IFRTs) are widely used to store light oil and chemical products. However, if the annular-rim gap around the floating deck becomes wider due to abrasion and aging of the sealing arrangement, the static breathing loss from the rim gap will be correspondingly aggravated. To investigate the oil-vapor migration and emissions from an IFRT, the effects of varying both the floating-deck height and wind speed on the oil-vapor diffusion were analyzed by performing numerical simulations and wind-tunnel experiments. The results demonstrate that the gas space volume and the wind speed of an IFRT greatly influence the vapor-loss rate of the IFRT. The larger the gas space volume, the weaker the airflow exchange between the inside and outside of the tank, thereby facilitating oil-vapor accumulation in the gas space of the tank. Furthermore, the loss rate of the IFRT is positively correlated with wind speed. Meanwhile, negative pressures and the vortexes formed on the leeward side of the tank. In addition, the higher concentration areas were mainly on the three vents on the downwind side of the IFRT. The results can provide important theoretical support for the design, management, and improvement of IFRTs.

1. Introduction

Internal floating-roof tanks (IFRTs) are widely used to store light oil, oil products, and chemical products. These tanks can significantly reduce the evaporating area of the stored liquids owing to the arrangement of a floating deck, thereby effectively restraining oil evaporation and reducing oil-vapor discharge from the tank. However, even if an IFRT is equipped with a sealing device on the annular-rim gap between the floating deck and the inner tank wall, the oil surface and the tank gas space cannot be completely isolated for the convenience of the floating deck moving up and down [1]. In addition, when the elasticity of the rim seal gradually decreases, and the gap of the rim seal is widened with long-term usage and abrasion, oil evaporation from the rim gap will gradually increase. Furthermore, oil evaporation from the stored liquid into the tank gas space and vapor emissions from the tank gas space into the atmosphere will cause oil loss, environmental pollution, and potential fire hazards [2]. Therefore, the interior mechanism of oil-vapor migration and emissions in IFRTs must be analyzed to ensure their safe operation.
Scientists have conducted many relevant research works, where the combination of numerical simulations and experimental measurements have been widely used as important tools for studying the oil-evaporation loss mechanism. Pasley et al. [3], and Zhao et al. [4], numerically simulated and experimentally measured the distribution of wind speeds and flow field around an external floating-roof tank and above the floating deck, and they observed that the airflow had two different flow characteristics each for the floating deck at lower positions and higher positions, respectively. Uematsuet et al. [5,6], investigated the wind-force distribution and the buckling behavior for open-topped oil-storage tanks. Wang et al. and Huang et al. [7,8,9], and Karbasian et al. [10], compared the effects of different oil-collection methods, oil-collection rates, and initial oil-vapor-mass fractions on the oil-vapor-diffusion law in the storage tank during the collection process of oil products. Hou et al. [11], analyzed the evolution of the flow field, temperature field, pressure field, oil-concentration field, and evaporation rate during the refueling process. In addition, many researchers have also studied other factors affecting oil-vapor emissions, such as vehicle and ship loading operations [12,13], the temperature-change characteristics of oil products in a storage tank [14,15,16,17], and the characteristics of oil vapors emitted from oil depots [18,19,20]. Subsequently, the researchers proposed a series of methods for assessing oil-vapor emissions [21,22]. However, these previously conducted studies have paid little attention to oil-vapor migration and emissions from an IFRT. The oil loss in a storage tank under normal operating conditions can be roughly divided into three processes: (1) heat and mass transfer between the liquid phase and the gas phase in the storage tank; (2) oil-vapor migration in the gas phase (i.e., in the gas space) of an IFRT; (3) oil-vapor emission and diffusion from the gas space into the atmosphere. In this study, the oil-evaporation rates were measured using self-made wind-tunnel experimental measurements. Subsequently, the species transfer model and the realizable k–ε model in the ANSYS Fluent software were used to simulate the oil-vapor diffusion process in an IFRT, following which the effects of floating-deck heights and ambient wind speeds were investigated.

2. Methodology

2.1. Experimental Protocol

To measure the oil-evaporation rate and the wind speed in the tank under different operating conditions, a model IFRT (1000 m3) was built according to the length ratio of 32:1, as depicted in Figure 1. The inner diameter, wall height, roof height, and rim gap of the tank were 360, 375, 39, and 6 mm, respectively. The size of the vents of the model tank was designed according to the selection principle of the prototype tank as follows:
B 0.06 D
where B denotes the total effective ventilation area of the model-tank vents, m2, and D is the inner diameter of the tank, m. The effective ventilation area of the model-tank vents should be greater than 0.69 m2, and the number of vents should not be less than four. Eventually, the model tank was equipped with four vents, each being 19 mm wide and 10 mm high. Further, the vents were evenly placed at the tank wall near the tank roof. In addition, four different floating-deck heights measuring 88, 176, 264, and 312 mm were set for the IFRT.
Due to the uncertainty of gasoline composition, the authors used n-hexane as the experimental oil for the convenience of research. The mass-difference method was used to measure the mass change of n-hexane in a certain period of time. The wind fields were generated using a self-made wind tunnel (DFWT-10), the size of the test section is 1.5 (H) × 1.5 (W) × 3 m (L), as depicted in Figure 1. The pitot tube anemometer is composed of a static pressure pitot tube (Kimo Instruments Co., Ltd., Bordeaux, France) and a digital micro-manometer (Yokogawa Electric Corporation, Musashino, Japan), which is used for real-time monitoring and feedback of wind speed in the wind tunnel [23]. The wind tunnel could provide wind fields with wind speeds ranging from 0.5 to 20 m·s−1. The wind speed, temperature, and humidity were measured using a hot-wire anemometer (TES-1341, TES Co., Ltd., Taiwan, China, having the following specifications: wind-speed range of 0–30 m·s−1 with a resolution of 0.01 m·s−1; temperature range of −10 to −60 °C with a resolution of 0.01 °C; humidity range of 10–95% Relative humidity (RH) with a resolution of 0.1% RH). The evaporation loss in the IFRT could be automatically measured using a high-precision electronic balance (WT-30000-1B, Xinheng Electronics Co., Ltd., Shanghai, China, having a range of 0–30 kg with a resolution of 0.01 g). The mass of n-hexane was measured before and after the experiment, and the change in the mass could be calculated as the mass loss of n-hexane during this period. Subsequently, the variation in the mass per unit time can also be calculated as the mass-loss rate of n-hexane. The experiments were arranged under the following conditions: the ambient temperature of 13 °C, and wind speeds of 4.36 and 6.36 m·s−1. The wind direction was directly opposite one of the four vents on the tank wall, and the mass-loss rates of n-hexane under different wind speeds were measured respectively.

2.2. Theoretical Models for Oil-Vapor Diffusion

2.2.1. Basic Governing Equations

The single-phase multicomponent diffusion problem without chemical reaction needs to be solved using the mass- and momentum-conservation equations. The mass-conservation equation can be written as follows:
ρ t + ( ρ v ) = S m
where Sm denotes the mass source term (in kg·m−3·s−1), and v the velocity vector (m·s−1).
The momentum-conservation equation can be written as follows:
t ( ρ v ) + · ( ρ v v ) = p + · [ μ ( v + v T ) ] + ρ g + F
where p denotes the static pressure (Pa), and ρ g and F denote the gravitational body force and external body force (both in N·m−3), respectively. Term μ denotes the dynamic viscosity (Pa·s).
The energy-conservation equation can be written as follows:
t ( ρ E ) + [ v ( ρ E + p ) ] = · ( k e f f T ) + S h
where energy E (J·kg−1) and temperature T (K) are mass average variables, keff the effective thermal conductivity (W·m−1·K−1), and Sh the energy source term (J·m−3·s−1).
In turbulent flows, the species-conservation equation can be written as follows:
t ( ρ Y i ) + ( ρ v Y i ) = ¯ · j i + S i
J i = ( ρ D i , m + μ t S c t ) Y i D T , i T T
where Yi denotes the local mass fraction of species i, and J i denotes the diffusion flux of species i (kg·m−2·s−1), which arises because of the gradients of concentration and temperature. Furthermore, Si denotes the source term (kg·m−3·s−1), Di,m the mass-diffusion coefficient for species i in the mixture, DT,i the thermal-diffusion coefficient, and Sct the turbulent Schmidt number.

2.2.2. Turbulence Equation

The oil vapors diffused from the IFRT are greatly affected by the ambient wind, and the oil-diffusion process is a complex unsteady turbulent flow. Therefore, the realizable k–ε model with better turbulent accuracy is selected. The following is the modeled transport equation for k and ε in the realizable k–ε model:
t ( ρ k ) + x j ( ρ k u j ) = x j [ ( μ + μ t σ k ) k x j ] + G k + G b ρ ε Y M + S k
and
t ( ρ ε ) + x j ( ρ ε u j ) = x j [ ( μ + μ t σ ε ) ε x j ] + ρ C 1 S ε ρ C 2 ε 2 k + ν ε + C 1 ε ε k C 3 ε G b + S ε
where
C 1 = m a x [ 0.43 , η η + 5 ] ,   η = S k ε ,   S = 2 S i j S i j
In the above-mentioned three equations, k denotes the turbulence kinetic energy(m−2·s−2) and ε the dissipation rate (m2·s−3). Furthermore, Gk and Gb denote the generation of the turbulence kinetic energy due to the mean velocity gradients and buoyancy, respectively (kg·m−1·s−1). YM denotes the contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate (kg·m−1·s−1). C2 and C are constants. σk and σε are the turbulent Prandtl numbers for k and ε, respectively. Sk denotes the source of the turbulence kinetic energy (kg·m−1·s−3) and Sε the source of the dissipation rate (kg·m−1·s−4).

2.3. Physical Model and Methodology

In wind-engineering calculations, the blocking ratio is usually employed as the basis for setting the cross-sectional area of the computational domain. The blocking ratio rb is defined as follows:
r b = A m A C
where Am and Ac denote the maximum windward area of the model and the cross-sectional area of the computational domain (both in m2), respectively. If the blocking ratio is less than 5%, the simulation results of the flow field may not be affected by each boundary [24,25,26]. An experimental model that is too small will increase the difficulty of experimental data measurement and lead to greater errors. We increased the size of the tank in the wind tunnel experiment to make the measurement results more exact. The blockage rate of the wind tunnel experimental model is 6%, calculated by Equation (9). To obtain the numerical-simulation results independent of the computational domain, the computational domain was selected as a cuboid (see Figure 2) of dimensions 2 (H) × 2 (W) × 6 (L), and the size of the model tank (see Figure 3) was consistent with that of the experimental model. A three-dimensional grid-type model of the IFRT was established using ICEM CFD 14.5 software (ANSYS Inc., Pittsburgh, PA, USA).
In the numerical calculation, the locations where the physical parameters in the flow field change require particular attention. These locations include the annular-rim gap between the floating deck and the inner tank wall of the IFRT, and the four vents on the tank wall. The grid of these regions must be partly encrypted to obtain more accurate numerical solutions. In this study, a structured grid type was designed for the computational domain. The number of cells was approximately 2 million; the quality of the grid was above 0.6, and its independency was also examined. The numerical-calculation procedure was performed using the commercial software package, ANSYS Fluent 14.5. Based on the species-transfer model, the measured oil-vapor-loss rates as the initial values of the mass transfer rate in the annular rim gap. The computational-domain inlet was set as the velocity boundary, and the computational-domain outlet was set as the pressure outlet. The vents on the tank wall were set as the interior boundary. The pressure and velocity were coupled using the SIMPLE scheme, and the spatial discretization of pressure was based on the Standard algorithm. The momentum-conservation equations were discretized using a second-order upwind scheme to reduce the numerical diffusion.
In wind tunnel experiments, we chose the inner floating roof tank as the experimental model. Re was calculated as 104,495 when the characteristic length is the model diameter, the wind velocity is 4.36 m·s−1, the atmospheric density is 1.205 kg·m−3, and the dynamic viscosity is 1.81 × 10−5 Pa·s. Actually, the model is 1/32 of the actual size, which requires that the inflow speed of the wind tunnel must be 32 times the actual wind speed to achieve Re equality, it is obviously difficult to achieve. Thus, it is unrealistic to strictly achieve Re equivalence. Some researchers [27,28] believe that the speed distribution would be independent of Re if the studied model Re is over the critical Re, which is Re-independence. Based on this, a good agreement was achieved between the CFD simulation and wind-tunnel experiment.

3. Results and Analysis

3.1. Verification of the Flow Field

Before using the numerical model, the consistency of the flow field of the numerical simulation must be experimentally verified. This experiment studied the speed boundary layer in the wind tunnel when no tank was placed in the wind tunnel. When the flow speed of the wind tunnel is controlled at 4 m·s−1, the wind speed in the X direction was measured at the center line of the section 500 mm away from the entrance of the wind tunnel test section. The measuring points are 20, 40, 60, 80, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, and 600 mm from the bottom, respectively. The experimental results are compared with the simulation results, listed in Figure 4. The boundary layer is 60 mm, so the tank is located in a uniform flow. According to Figure 4, the numerical simulation can better simulate the boundary layer conditions of the wind tunnel experiment.
Four points on the tank circumference near the four vents were selected (assuming the center of the tank bottom as the origin). The ambient wind speed was set to 4.36 m·s−1, and the floating-deck height was set to 176 mm. The hot-wire anemometer (TES-1341, TES Co., Ltd., Taiwan, China) was used to measure the wind speeds at these points. Each point was measured five times, and the wind-speed results are presented in Table 1.
It can be seen from Table 1 that the average wind speeds at Points 3 and 4 are approximately equal, both of which are slightly greater than the ambient wind speed of 4.36 m·s−1 owing to the effect of the tank body. However, because the airflow was blocked by the tank wall, the average wind speed at Point 2 (which lies in the windward side of the tank) is lower than the ambient wind speed. In addition, the average wind speed at Point 1 (which lies in the back side of the tank) is the smallest. The four measured wind-speed values at Points 1 to 4 basically accord with the distribution law of the flow field around the tank. However, the measured values and the numerically simulated values are slightly different from each other. The reasons for the difference are as follows: (1) although the flow field is stabilized in general, the uncertainty of turbulence exists partly; the physical parameter (i.e., wind speed) in the flow field fluctuates within a certain range; the measured values are also slightly different at different moments; (2) the error in the measuring instrument itself and the interference of the measuring instrument with the flow field can also cause measurement errors. Therefore, the consistency of the flow field of the numerical simulation was verified using the experiment, and the numerical-simulation model can be considered to be relatively suitable for oil-vapor-diffusion simulation of an oil tank.

3.2. Influence of Different Floating-Deck Heights on N-Hexane-Loss Rates

It can be seen from the experimental results in Table 2 that the loss rates of n-hexane can be influenced greatly by varying both floating-deck heights and wind speeds. In general, the floating-deck heights of the IFRT were positively correlated with the evaporation rates of n-hexane. However, the loss rates increased as the floating-deck heights and wind speeds increased. For the ambient wind speeds of 4.36 and 6.36 m·s−1, setting the floating-deck height to 312 mm (the highest height), the loss rates of n-hexane reached 4.139 × 10−6 and 5.960 × 10−6 kg·s−1, respectively. However, upon setting the floating-deck height to 88 mm (the lowest height), the loss rates were only 2.189 × 10−6 and 3.773 × 10−6 kg·s−1, respectively.
The reasons for this phenomenon are as follows: the evaporation of n-hexane is driven by a diffusion process that is induced by the gradient of n-hexane-vapor concentration at the liquid surface. The larger the gas space inside the tank, the slower the gas flow, and thus the less frequent the gas exchange between the tank and the atmosphere. Therefore, the n-hexane-vapor concentration remains at a relatively stable level, and the driving force for n-hexane evaporation was reduced. Conversely, the smaller the gas space inside the tank, the more frequent the gas exchange between the tank and the atmosphere. As a result, the n-hexane vapor in the tank released into the atmosphere, resulting in fresh air entering the tank from the vents. The n-hexane-vapor concentration in the tank was maintained at a low value, increasing the n-hexane-vapor concentration difference between the gas space and the gas–liquid interface, thereby aggravating the evaporation of n-hexane.

3.3. Simulation of the Large IFRT

To investigate the concentration-distribution characteristics of the light oil products in the actual IFRT, the authors continued to simulate the oil-vapor-diffusion process in a large IFRT in the length ratio of 1:1 for a 1000 m3 IFRT. N-hexane was still chosen as the stored oil product. The inner diameter, wall height, roof height, and rim gap of the tank were 11.5, 12.0, 1.254, and 0.2 m, respectively. The outlet boundary of the flow field was set as the pressure outlet boundary condition while the gap between the floating deck and the tank wall was set as the mass-flow boundary condition. The tank bottom, tank wall, and floating deck were all set as no-slip boundaries and the ambient temperature was set at 13 °C. The mass fraction of the saturated concentration of n-hexane when it diffuses in the gap of the floating disk is 0.30, according to the calculation of the saturated vapor pressure of n-hexane at 13 °C. We took the k–ε turbulence model to describe the diffusion process of the oil vapor.

3.3.1. Flow-Field Simulation

The streamtraces inside the tank with different floating-deck heights for the ambient wind speed of 4.36 m·s−1 were simulated, as depicted in Figure 5, Figure 6, Figure 7 and Figure 8. As can be seen from Figure 5, in the larger gas space of the tank, a clockwise vortex was formed, and the streamtraces of the vortex resembled the “0” type. As the gas space became smaller, the streamtraces of the vortex in the tank gradually assumed the shape of a flat ellipse (see Figure 6 and Figure 7). Finally, the gas space became sufficiently small, following which the vortex mentioned above was divided into several small vortexes; it means that the airflow in the gas space was faster and more disorderly at this moment. Furthermore, the streamtraces outside the tank were simulated, as depicted in Figure 9. It can be seen from the figure that a part of the airflow from the inlet entered the tank and that the other part of the airflow bypassed the tank. Moreover, the airflow behind the tank did not flow in the original direction, and many vortexes were generated there.
The diffusion of the n-hexane vapors evaporating from the tank was greatly affected by the ambient wind. In the numerical simulation, the ambient wind was set along the X-axis, and the wind speeds in the inlet were set to 4.36 and 6.36 m·s−1, respectively. The velocity-contour distributions both inside and outside the IFRT in the XY plane are depicted in Figure 10 and Figure 11. The results demonstrate that the airflow was obstructed by the tank at the windward side and that the velocity of the airflow was gradually reduced from 4.36 to 0 m·s−1. Combined with the pressure-cloud diagram at the windward side of the tank wall in Figure 12, it can be concluded that the pressure on the windward side of the tank body increased because the wind speed translated into stagnation pressure. On the surface of the tank roof, the flow cross-section of the ambient wind shrunk, causing a sudden increase in the wind speed and negative pressure (see Figure 13). The wind speed near the tank roof exceeded the ambient wind speed, and the maximum wind-speed value reached 5.5 and 8.0 m·s−1, respectively. On the leeward side of the tank wall, the wind-speed iso-surface continuously shrunk toward the bottom of the tank because of the obstruction of the tank body (see Figure 10 and Figure 11). The minimum value of the wind speed was approximately 0 m·s−1 at the bottom of the leeward side of the tank wall. Compared to the windward side, the region of lower wind speed was larger in the leeward side. Combined with the wind-flow diagram of the inside and outside of the tank (see Figure 9) and the tank-wall pressure-cloud diagram on the leeward side of the tank (see Figure 13), it can be seen that negative pressure and vortexes appeared near the leeward side of the tank.
Several monitored points were set along the central axis of the tank to obtain the data of wind speed inside the tank, as depicted in Figure 14. The wind speed in the tank was low because the sizes of the vents were too small relative to the size of the gas space. It can be seen from Figure 15 that for the ambient wind speed of 4.36 m·s−1, the velocities of the airflow in the gas space of the tank were less than 0.6 m·s−1. In addition, upon lowering the floating-deck height, the speeds measured from the monitored points changed greatly. The wind speeds at the points near the roof and floating deck were higher than those at other monitored points, while the wind speed at the center of the gas space was low. The minimum wind speed in the gas space was approximately 0.05 m·s−1. Upon increasing the height of the floating deck, the low-velocity region gradually reduced. When the floating-deck height reached a certain height, the velocity distribution became irregular. This irregularity reconfirmed the conclusions drawn previously that upon changing the floating-deck height, the streamtraces of the airflow in the tank also changed.

3.3.2. Concentration Distribution of N-Hexane Vapor in the Tank

The concentration distribution at different floating-deck heights in the tank is depicted in Figure 16. It can be seen that upon increasing the wind speed under a fixed floating-deck height, the concentration of n-hexane in the gas space in the tank becomes low. The higher concentration of n-hexane vapors was mainly distributed near the rim gap, and the concentration distribution of n-hexane vapors in the upper gas space of the tank was relatively uniform. Upon setting the floating-deck height to 2.82 m, because of the resulting large gas space in the tank, the airflow inside and outside the tank is exchanged infrequently, and thus the concentration of n-hexane in the tank becomes evenly distributed. Upon increasing the floating-deck height to 5.63 m, the n-hexane-vapor concentration was affected by the air inflow, which formed a large vortex in the tank. Upon further raising the floating-deck height to 8.45 m, the gas exchange between the inside and outside of the tank was accelerated, and thus the large vortex in the tank was destroyed (wind speed = 6.36 m·s−1) or compressed (wind speed = 4.36 m·s−1) because of the reduction of the gas space in the tank. Furthermore, upon increasing the floating-deck height to 9.98 m, the n-hexane vapors were mainly concentrated in the area near the gap. In addition, n-hexane vapors at other regions were quickly discharged out of the tank along with the airflow, and it became difficult for the n-hexane vapors to stay in the tank.
The n-hexane-vapor emission from the IFRT was also investigated; the floating-deck height was at 9.98 m, wind speed 4.36 m·s−1, and the concentration distribution as depicted in Figure 17. It can be seen that the n-hexane vapors discharged from the vents on the front and rear sides of the tank were diffused along the wind direction and that the vapors discharged from the leeward-side vent spread around the vent. Irrespective of the front and rear sides or the leeward side of the tank (the vent on the windward side is the airflow inlet), the mass fraction of n-hexane vapors near the vents reached 0.002, and such high concentration of vapors discharged into the atmosphere do not meet environmental safety requirements and will cause air pollution. Therefore, if equipping the vents with an oil-vapor-recovery unit is necessary, priority should be given to the above-mentioned three vents.

4. Conclusions

The evaporative loss rates of n-hexane in the IFRT model and the wind speeds of gas space above the floating deck inside the tank were measured using wind-tunnel experiments and also simulated using ANSYS Fluent software. The conclusions drawn from the research results can be summarized as follows:
(1)
Based on numerical simulation and the wind-tunnel experiments, the oil-vapor diffusion process in the IFRT was simulated and was then verified to be relatively suitable to the oil-vapor-diffusion simulation for different sizes of IFRTs. This further revealed the law of mass transfer between the oil vapors and air for the evaporation and diffusion process in the IFRTs.
(2)
Different floating-deck heights of the IFRT corresponded to different loss rates of n-hexane. The larger the gas space inside the tank, the weaker the airflow exchange between the inside and outside of the tank became. Therefore, the gradient of the n-hexane-vapor in the tank was lower, thereby reducing the driving force for n-hexane evaporation.
(3)
Outside the tank, the direction of the ambient wind would change suddenly while bypassing the tank. In addition, negative pressure and vortexes were generated on the leeward side of the tank. Inside the tank, the n-hexane-vapor distribution was relatively uniform at lower floating-deck heights. Upon increasing the floating-deck height, a large vortex formed in the tank, intensifying the airflow disturbance in the tank. Upon further increasing the floating-deck height, the aforementioned large vortex was divided into several small vortexes, and the airflow in the gas space became faster and more disorderly than that in the previous floating-height condition.
(4)
The higher vapor-concentration regions near the vents, compared with other regions in the tank, are mainly concentrated on the front and rear sides and the leeward side of the IFRT. Therefore, if it is decided to employ an oil-vapor-recovery unit, more attention should be paid to the regions having higher vapor concentrations.

Author Contributions

Conceptualization, G.Z., W.H.; software, F.H.; formal analysis, W.H., G.Z., F.H., L.F. and J.F.; investigation, Z.Z.; data curation, W.H., G.Z., and H.J.; writing—original draft preparation, G.Z., W.H., and X.S.; writing—review and editing, G.Z., W.H., F.H., L.F. and Z.Z.; supervision, W.H.; funding acquisition, W.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51574044 and No. 51804045), the Key Research and Development Program of Jiangsu Province (Industry Foresight and Common Key Technology) (No. BE2018065), the Sci & Tech Program of Changzhou (No. CJ20180053) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (NO. SJCX19_0668).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, W.Q.; Huang, F.Y.; Fang, J.; Fu, L.P. A calculation method for the numerical simulation of oil products evaporation and vapor diffusion in an internal floating-roof tank under the unsteady operating state. J. Petrol. Sci. Eng. 2020, 188, 106867. [Google Scholar] [CrossRef]
  2. Zhu, Y.; Qian, X.M.; Liu, Z.Y.; Huang, P.; Yuan, M.Q. Analysis and assessment of the Qingdao crude oil vapor explosion accident: Lessons learnt. J. Loss Prev. Process Ind. 2015, 33, 289–303. [Google Scholar] [CrossRef]
  3. Pasley, H.; Clark, C. Computational fluid dynamics study of flow around floating-roof oil storage tanks. J. Wind Eng. Ind. Aerodyn. 2000, 86, 37–54. [Google Scholar] [CrossRef]
  4. Zhao, C.L.; Huang, W.Q.; Zhong, J.; Wang, W.J.; Xu, X.Y.; Wang, Y.X. Numerical simulation of oil vapor leakage from external floating-roof tank. CIESC J. 2014, 65, 4203–4209. [Google Scholar]
  5. Uematsu, Y.; Koo, C.; Yasunaga, J. Design wind force coefficients for open-topped oil storage tanks focusing on the wind-induced buckling. J. Wind Eng. Ind. Aerodyn. 2014, 130, 16–29. [Google Scholar] [CrossRef]
  6. Uematsu, Y.; Yasunaga, J.; Koo, C. Design wind loads for open-topped storage tanks in various arrangements. J. Wind Eng. Ind. Aerodyn. 2015, 138, 77–86. [Google Scholar] [CrossRef]
  7. Wang, Z.L.; Huang, W.Q.; Ji, H. Numerical simulation of vapor diffusion and emission in loading gasoline into dome roof tank. Acta Petrol. Sin. Petrol. Process. Sect. 2017, 33, 26–33. [Google Scholar]
  8. Huang, W.Q.; Wang, Z.L.; Ji, H.; Zhao, C.L.; Lv, A.H.; Xu, X.Y.; Wang, Y.H. Experimental determination and numerical simulation of vapor diffusion and emission in loading gasoline into tank. CIESC J. 2016, 67, 4994–5005. [Google Scholar]
  9. Jing, H.B.; Huang, W.Q.; Ji, H.; Hu, B.G. Study on similarity criteria number of wind tunnel experiment for oil vapor diffusion based on numerical simulation technology. J. Changzhou Univ. 2020, 32, 83–92. [Google Scholar]
  10. Karbasian, H.R.; Kim, D.Y.; Yoon, S.Y.; Ahn, J.H.; Kim, K.C. A new method for reducing VOCs formation during crude oil loading process. J. Mech. Sci. Technol. 2017, 31, 1701–1710. [Google Scholar] [CrossRef]
  11. Hou, Y.; Chen, J.; Zhu, L.; Sun, D.; Li, J.; Jie, P. Numerical simulation of gasoline evaporation in refueling process. Adv. Mech. Eng. 2017, 9, 1687814017714978. [Google Scholar] [CrossRef] [Green Version]
  12. Milazzo, M.F.; Ancione, G.; Lisi, R. Emissions of volatile organic compounds during the ship-loading of petroleum products: Dispersion modelling and environmental concerns. J. Environ. Manag. 2017, 204, 637–650. [Google Scholar] [CrossRef] [PubMed]
  13. Hata, H.; Yamada, H.; Kokuryo, K.; Okada, M.; Funakubo, C.; Tonokura, K. Estimation model for evaporative emissions from gasoline vehicles based on thermodynamics. Sci. Total Environ. 2018, 618, 1685–1691. [Google Scholar] [CrossRef]
  14. Sun, W.; Cheng, Q.; Zheng, A.; Zhang, T.; Liu, Y. Research on coupled characteristics of heat transfer and flow in the oil static storage process under periodic boundary conditions. Int. J. Heat Mass Transf. 2018, 122, 719–731. [Google Scholar] [CrossRef]
  15. Wang, M.; Zhang, X.; Yu, G.; Li, J.; Yu, B.; Sun, D. Numerical study on the temperature drop characteristics of waxy crude oil in a double-plate floating roof oil tank. Appl. Therm. Eng. 2017, 124, 560–570. [Google Scholar] [CrossRef]
  16. Wang, M.; Yu, B.; Zhang, X.; Yu, G.; Li, J. Experimental and numerical study on the heat transfer characteristics of waxy crude oil in a 100,000 m3 double-plate floating roof oil tank. Appl. Therm. Eng. 2018, 136, 335–348. [Google Scholar] [CrossRef]
  17. Sun, W.; Cheng, Q.; Zheng, A.; Gan, Y.; Gao, W.; Liu, Y. Heat flow coupling characteristics analysis and heating effect evaluation study of crude oil in the storage tank different structure coil heating processes. Int. J. Heat Mass Transf. 2018, 127, 89–101. [Google Scholar] [CrossRef]
  18. Jovanovic, J.; Jovanovic, M.; Jovanovic, A.; Marinovic, V. Introduction of cleaner production in the tank farm of the Pancevo Oil Refinery, Serbia. J. Clean. Prod. 2010, 18, 791–798. [Google Scholar] [CrossRef]
  19. Kountouriotis, A.; Aleiferis, P.G.; Charalambides, A.G. Numerical investigation of VOC levels in the area of petrol stations. Sci. Total Environ. 2014, 470, 1205–1224. [Google Scholar] [CrossRef] [Green Version]
  20. Wei, W.; Lv, Z.; Yang, G.; Cheng, S.; Li, Y.; Wang, L. VOCs emission rate estimate for complicated industrial area source using an inverse-dispersion calculation method: A case study on a petroleum refinery in Northern China. Environ. Pollut. 2016, 218, 681–688. [Google Scholar] [CrossRef]
  21. Sharma, Y.; Majhi, A.; Kukreti, V.; Garg, M. Stock loss studies on breathing loss of gasoline. Fuel 2010, 89, 1695–1699. [Google Scholar] [CrossRef]
  22. Tamaddoni, M.; Sotudeh-Gharebagh, R.; Nario, S.; Hajihosseinzadeh, M.; Mostoufi, N. Experimental study of the VOC emitted from crude oil tankers. Process Saf. Environ. Protect. 2014, 92, 929–937. [Google Scholar] [CrossRef]
  23. Klopfenstein, R., Jr. Air velocity and flow measurement using a Pitot tube. ISA Trans. 1998, 37, 257–263. [Google Scholar] [CrossRef]
  24. Baetke, F.; Werner, H.; Wengle, H. Numerical simulation of turbulent flow over surface-mounted obstacles with sharp edges and corners. J. Wind Eng. Ind. Aerodyn. 1990, 35, 129–147. [Google Scholar] [CrossRef]
  25. Saathoff, P.; Melbourne, W. Freestream turbulence and wind tunnel blockage effects on streamwise surface pressures. J. Wind Eng. Ind. Aerodyn. 1987, 26, 353–370. [Google Scholar] [CrossRef]
  26. Holmes, J.D. Wind Loading of Structures; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  27. Uehara, K.; Wakamatsu, S.; Ooka, R. Studies on critical Reynolds number indices for wind-tunnel experiments on flow within urban areas. Bound. Layer Meteor. 2003, 107, 353–370. [Google Scholar] [CrossRef]
  28. Cui, P.Y.; Li, Z.; TAO, W.Q. Wind-tunnel experiment and numerical analysis on flow and dispersion in street canyons. J. Eng. Thermophys. 2014, 35, 2491–2495. [Google Scholar]
Figure 1. Wind tunnel for the experiment and simulation.
Figure 1. Wind tunnel for the experiment and simulation.
Processes 08 01026 g001
Figure 2. Grid model for oil-vapor diffusion of the Internal floating-roof tank (IFRT).
Figure 2. Grid model for oil-vapor diffusion of the Internal floating-roof tank (IFRT).
Processes 08 01026 g002
Figure 3. Grid model for the IFRT.
Figure 3. Grid model for the IFRT.
Processes 08 01026 g003
Figure 4. Comparison of wind speed in the X direction between the wind-tunnel experiment and numerical simulation on the inlet of 500 mm.
Figure 4. Comparison of wind speed in the X direction between the wind-tunnel experiment and numerical simulation on the inlet of 500 mm.
Processes 08 01026 g004
Figure 5. Streamtraces inside the tank for the floating-deck height of 2.82 m.
Figure 5. Streamtraces inside the tank for the floating-deck height of 2.82 m.
Processes 08 01026 g005
Figure 6. Streamtraces inside the tank for the floating-deck height of 5.63 m.
Figure 6. Streamtraces inside the tank for the floating-deck height of 5.63 m.
Processes 08 01026 g006
Figure 7. Streamtraces inside the tank for the floating-deck height of 8.45 m.
Figure 7. Streamtraces inside the tank for the floating-deck height of 8.45 m.
Processes 08 01026 g007
Figure 8. Streamtraces inside the tank for the floating-deck height of 9.98 m.
Figure 8. Streamtraces inside the tank for the floating-deck height of 9.98 m.
Processes 08 01026 g008
Figure 9. Streamtraces inside and outside the tank.
Figure 9. Streamtraces inside and outside the tank.
Processes 08 01026 g009
Figure 10. Velocity distribution of the IFRT under the ambient wind speed of 4.36 m·s−1.
Figure 10. Velocity distribution of the IFRT under the ambient wind speed of 4.36 m·s−1.
Processes 08 01026 g010
Figure 11. Velocity distribution of the IFRT under the ambient wind speed of 6.36 m·s−1.
Figure 11. Velocity distribution of the IFRT under the ambient wind speed of 6.36 m·s−1.
Processes 08 01026 g011
Figure 12. Pressure-cloud diagram for the windward side of the tank wall (Pa).
Figure 12. Pressure-cloud diagram for the windward side of the tank wall (Pa).
Processes 08 01026 g012
Figure 13. Pressure-cloud diagram for the leeward side of the tank wall (Pa).
Figure 13. Pressure-cloud diagram for the leeward side of the tank wall (Pa).
Processes 08 01026 g013
Figure 14. Location of monitored points.
Figure 14. Location of monitored points.
Processes 08 01026 g014
Figure 15. Velocities inside the tank for different floating-deck heights.
Figure 15. Velocities inside the tank for different floating-deck heights.
Processes 08 01026 g015
Figure 16. Distribution of n-hexane vapors (mass fraction) in the IFRT at different floating-deck heights. (a) Wind speed = 4.36 m·s−1, floating-deck height = 2.82 m. (b) Wind speed = 6.36 m·s−1, floating-deck height = 2.82 m. (c) Wind speed = 6.36 m·s−1, floating-deck height = 5.63 m. (d) Wind speed = 4.36 m·s−1, floating-deck height = 5.63 m. (e) Wind speed = 6.36 m·s−1, floating-deck height = 8.45 m. (f) Wind speed = 4.36 m·s−1, floating-deck height = 8.45 m. (g) Wind speed = 6.36 m·s−1, floating-deck height = 9.98 m. (h) Wind speed = 4.36 m·s−1, floating-deck height = 9.98 m.
Figure 16. Distribution of n-hexane vapors (mass fraction) in the IFRT at different floating-deck heights. (a) Wind speed = 4.36 m·s−1, floating-deck height = 2.82 m. (b) Wind speed = 6.36 m·s−1, floating-deck height = 2.82 m. (c) Wind speed = 6.36 m·s−1, floating-deck height = 5.63 m. (d) Wind speed = 4.36 m·s−1, floating-deck height = 5.63 m. (e) Wind speed = 6.36 m·s−1, floating-deck height = 8.45 m. (f) Wind speed = 4.36 m·s−1, floating-deck height = 8.45 m. (g) Wind speed = 6.36 m·s−1, floating-deck height = 9.98 m. (h) Wind speed = 4.36 m·s−1, floating-deck height = 9.98 m.
Processes 08 01026 g016aProcesses 08 01026 g016b
Figure 17. Distribution of n-hexane-vapor concentration (mass fraction) near the vents.
Figure 17. Distribution of n-hexane-vapor concentration (mass fraction) near the vents.
Processes 08 01026 g017
Table 1. X-axis wind speeds at each measuring point of the tank.
Table 1. X-axis wind speeds at each measuring point of the tank.
Measuring Point and Its Coordinates, mmMeasured Value, m·s−1Simulated Value, m·s−1Error,%
Test 1Test 2Test 3Test 4Test 5Average
Point 1 (190, 360, 0)0.540.530.650.720.610.610.4821.3
Point 2 (−190, 360, 0)2.632.712.622.642.682.662.573.4
Point 3 (0, 360, 190)5.624.654.704.714.514.845.4011.6
Point 4 (0, 360, −190)4.764.644.624.514.664.645.3916.7
Table 2. Loss rates of n-hexane measured in model tanks *.
Table 2. Loss rates of n-hexane measured in model tanks *.
Floating-Deck Height, mmFull Coefficientr1/(10−6 kg·s−1)r2/(10−6 kg·s−1)
880.242.1893.773
1760.482.7194.453
2640.713.5605.272
3120.844.1395.960
*: r1, r2—The loss rates measured for the wind speed of 4.36 and 6.36 m·s−1, respectively.

Share and Cite

MDPI and ACS Style

Zhang, G.; Huang, F.; Huang, W.; Zhu, Z.; Fang, J.; Ji, H.; Fu, L.; Sun, X. Analysis of Influence of Floating-Deck Height on Oil-Vapor Migration and Emission of Internal Floating-Roof Tank Based on Numerical Simulation and Wind-Tunnel Experiment. Processes 2020, 8, 1026. https://doi.org/10.3390/pr8091026

AMA Style

Zhang G, Huang F, Huang W, Zhu Z, Fang J, Ji H, Fu L, Sun X. Analysis of Influence of Floating-Deck Height on Oil-Vapor Migration and Emission of Internal Floating-Roof Tank Based on Numerical Simulation and Wind-Tunnel Experiment. Processes. 2020; 8(9):1026. https://doi.org/10.3390/pr8091026

Chicago/Turabian Style

Zhang, Gao, Fengyu Huang, Weiqiu Huang, Zhongquan Zhu, Jie Fang, Hong Ji, Lipei Fu, and Xianhang Sun. 2020. "Analysis of Influence of Floating-Deck Height on Oil-Vapor Migration and Emission of Internal Floating-Roof Tank Based on Numerical Simulation and Wind-Tunnel Experiment" Processes 8, no. 9: 1026. https://doi.org/10.3390/pr8091026

APA Style

Zhang, G., Huang, F., Huang, W., Zhu, Z., Fang, J., Ji, H., Fu, L., & Sun, X. (2020). Analysis of Influence of Floating-Deck Height on Oil-Vapor Migration and Emission of Internal Floating-Roof Tank Based on Numerical Simulation and Wind-Tunnel Experiment. Processes, 8(9), 1026. https://doi.org/10.3390/pr8091026

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop