3.1. Initial Position of HSI Occurrence with an Elevated Hot Surface Temperature
The fuel leaked in the ship engine room evaporates through heating on a hot surface, forms combustible vapor above the hot surface, and mixes with the surrounding air. The combustible vapor mixture ignites over the hot surface. In the process of heat evaporation, the leaking marine fuel also causes heat loss to air and the surface housing the equipment due to the action of the heat source. In the process of fuel his, the leaking fuel is affected by air flow organization, environmental humidity, and wall thermal feedback. As a result, the leaking marine fuel’s HSI characteristics, ignition time and ignition probability are changed. The heat loss of combustible liquid and the effect of environmental thermal feedback are affected by the specific design of the top opening in the engine room and changes in humidity in the operating environment. NFPA 921 suggests that ventilation can influence ignition. If too much air blows across a hot surface, it can hinder the ability of the combustible liquid to volatilize quickly enough to allow ignition. Conversely, if too little air blows across the hot surface, the atmosphere can become too rich. According to Equation (1) [
33], the heat received by the surface of marine fuel may be the sum of heat conducted by the hot surface, heat transferred to the liquid surface by thermal convection, and heat transferred to the liquid surface by thermal radiation.
where
Qconv is the heat generated by thermal convection, kW/m
2;
Qrad is heat transferred by thermal radiation, kW/m
2;
Qex is heating rate of external heat source per unit area of liquid surface, kW/m
2; and
Qloss is heat loss rate per unit area of combustible liquid, kW/m
2.
According to the experimental tests in the previous section, the process of the ignition of leaking marine diesel above a hot surface is accompanied by significant and complex evaporation behavior. In this study, the evaporation and ignition characteristics of marine diesel in the vertical space of a hot surface are analyzed. The data are collected by a high-speed camera, and each image is processed according to a time axis. The experimental amount of marine diesel is 10 mL and the hot surface temperature is set at 748 K as an example.
Figure 2 presents the process of diesel leaking onto a high-temperature surface, causing evaporation and initial ignition. When marine diesel leaks onto a hot surface in the initial stage, the rapid phase transition of marine diesel can be clearly observed. Due to the high temperature of the hot surface, marine diesel is almost always in a state of film boiling. This creates a gaseous layer between the diesel and high-temperature surface, which changes the mode of heat transfer from the high-temperature surface to liquid. At 3.675 s, a significant amount of white vapor mixture is produced when marine diesel comes into contact with a hot surface. The marine diesel produces a wider range of combustible vapor mixtures above hot surfaces when compared to marine lubricant. This is due to the flow characteristics of marine diesel. The evaporation product of diesel on the high-temperature surface mixes fully with fresh air, forming a combustible mixture in the vertical space. The heat continually transfers from the high-temperature surface to the liquid surface and, coupled with thermal feedback, causes the initial ignition behavior of the high-temperature surface; ignition therefore occurs instantaneously when diesel leaks at 3.725 s. Once the initial ignition forms a core in the air, it gradually grows and generates flame propagation, as depicted at 3.733 s. This reveals that the initial ignition takes place at a higher position above the high-temperature surface.
The interface between the leaked fuel and high-temperature hot surface is a solid–liquid-state heat transfer interface. At this interface, the high-temperature surface of the ship engine system transfers heat to leaking oil by means of heat conduction. The hybrid scenario shows faster movement of liquid marine fuel due to the formation of vapor underneath and the curvature of the leading edge of leaked marine fuel caused by surface tension. If the high-temperature surface is significantly elevated, the marine fuel will boil violently on the high-temperature surface, and a film of vapor will form on the surface of the liquid fuel. The film’s boiling mode can cause the pool of marine fuel to break up into smaller globules due to surface tension. When film boiling occurs in the marine fuel above the hot surface, combustible vapor is generated at different distances from the marine fuel’s surface. The film boiling regime of leaked marine fuel starts at a specific value of excess
Ts over
Tsat and can be given by Equation (2) [
34].
where
ρv is the density of fuel released into ambient air from the hot surface, kg/m
3;
ρl is the density of marine fuel, kg/m
3;
λv is the thermal conductivity of vapor in the film between the marine fuel and hot surface, W/(m·K);
hv is the latent heat of vaporization, kJ/g;
σ is the surface tension of marine fuel; and N/m;
μv is the viscosity of the vapor in the film between the surface and fuel, mm
2/s.
When marine fuel on a high-temperature surface is in a film boiling state,
hf can be presented [
35] as follows.
where
hf is the heat transfer coefficient of film boiling, W/m
2·K.
In the film boiling state, the evaporation rate of fuel on a high-temperature surface can be obtained by Equation (4) [
36], as shown:
where
mf is the evaporation rate of leaked marine fuel within unit area, g/s;
S is the surface cross-sectional area from leaking marine fuel on the hot surface, m
2; and
B is the mass transfer number.
The mass transfer number is also called the D. B. Spalding number, which represents the ratio of chemical energy to the evaporation energy of fuel.
B can be found using Equation (5) [
37]:
where
Tb is the boiling temperature, K, and
cp is the specific heat, J/(kg·K).
Figure 2 illustrates that when the marine diesel leaks and is ignited above a hot surface, the initial flame entrains air in the engine system. The experimental phenomenon demonstrates that without mechanical ventilation, ignition always occurs vertically upward from the center line of the high-temperature surface. Marine diesel HSI typically occurs between the heated surface and the ceiling of the ship engine system. The initial flame appears on the hot surface and further spreads. This shows that the accumulation of marine diesel vapor above the heating surface, combined with the hot surface temperature, affects the entrainment of the vapor/air mixture. Equation (6) shows the relationship between the entrainment velocity of the high-temperature surface ignition flame and the rising velocity of the gaseous mixture in the initial fire core [
38].
where
ue is flame velocity at which the vapor/air mixture is entrained in the horizontal direction, m/s;
ub is the rising velocity, which is attributed to vapor and buoyancy effect, m/s; and
η is the correction coefficient.
Equation (7) expresses that the rising velocity decreases as the distance from the high-temperature surface increases in the vertical space in which a gaseous mixture is formed by the marine diesel vapor and air existing in the engine room [
39].
where
qs is the heat flux density of the solid–liquid surface, W/m
2; g is the gravitational acceleration, m/s
2;
H is the marine diesel ignition height, m;
ρ0 is the ambient air density, kg/m
3;
Ta is the initial temperature of ambient air, K; and
cp is specific heat, J/(kg·K).
As a result, by obtaining the correction coefficient, we can further determine the expression of the entrainment velocity at which the initial ignition flame is on the hot surface.
Figure 3 presents the difference between the initial ignition height of marine diesel and the heated surface’s temperature. The ignition height of marine diesel exhibits a declining trend with an increase in heated surface temperature. However, this downward trend does not follow a simple linear relationship, as can be seen in the experimental results when the probability of HSI is 50%, for instance. When the heated surface temperature is close to 748 K, the ignition height of marine diesel above the hot surface is about 0.51 m. If the surface is heated and reaches a temperature of 758 K, the ignition height drops to approximately 0.42 m, with a change amplitude of 0.009 m/K. If the surface reaches a temperature of 778 K, the diesel will be ignited on the hot surface with 100% probability. At this high surface temperature, the height of the HSI for marine diesel is 0.22 m. The increase in surface temperature causes the marine diesel that leaks onto the hot surface to evaporate rapidly, forming combustible vapor that mixes with the fresh air. The high-temperature surface has a great thermal effect in a short time, leading to HSI behavior in the area near the high-temperature surface. As the heated surface’s temperature goes up, the ignition height remains around 0.18 m above the hot surface instead of fluctuating greatly. In a previous work [
40], a vapor plume model was developed. This model is presented in Equation (8) and can be used to describe the relationship between the ignition height of marine diesel and the initial specific buoyancy flux [
41] in a certain vertical space.
where
ρ0 is the density of ambient air, kg/m
3;
ρv is the density of marine diesel released into the ambient air from the hot surface, kg/m
3;
H is the ignition height of marine diesel, m;
Bi is the initial specific buoyancy flux, m
4/s
3;
a is the correction coefficient, dimensionless; and
Qv is the volumetric flow rate of marine diesel, m
3/s.
Using the experimental data, a relationship between initial ignition height and hot surface temperature is developed in the current study, as proposed in Equation (10). This model is expected to predict the initial location of ignition given various heated surface temperatures when marine diesel leaks into a ship’s engine room.
where
H is height at which the ignition of marine diesel initially occurs, m;
H0 is the surface height of marine diesel, m;
Ts is the hot surface temperature, K; and
T0 is the initial temperature of marine diesel, K.
To alleviate the issue of increasing hot surface-driven initial ignition height, the following steps can be taken:
- (a)
The on-site safety manager can regularly inspect the vicinity of pipelines, paying special attention to thermal insulation measures to protect the hot surfaces of equipment and thus preventing the ignition of accidentally leaked marine diesel by a hot surface;
- (b)
An existing monitoring gauge in ship engine room may be used. The research results of this section suggest installing an area covering the vertical space of the high-temperature surface of the equipment in the ship engine room in order to monitor the possible initial ignition of a fire core in different positions.
3.2. Temperature Field Distribution according to HSI-Driven Flame Propagation
After the marine diesel is ignited on a high-temperature surface, the size of the fire core rapidly increases, and flame propagation forms. According to the relationship between the initial temperature of combustible fuel and its flash point, flame propagation on the liquid surface can be classified into liquid-phase controlled fire spreading and gas-phase controlled fire spreading. If ignition occurs on a hot surface, the flame propagates rapidly through premixed combustion. The size of the marine diesel leakage area on the high-temperature surface has a small influence on the flame propagation.
Figure 4 shows the transmission of combustion flames after the marine diesel is ignited above the hot surface. When marine diesel leakage occurs at 3.741 s, the fire nucleus generated in the air begins to grow larger. Because the flame consumes more combustible vapor mixture, it needs more fuel. The flame gradually begins to propagate in the direction of the hot surface. This reflects the energy transfer from the flame in the chemical reaction zone to the unburned gas mixture in the preheating zone. The flame trajectory of the marine diesel vapor in the vertical space can be observed in
Figure 4. When the flame is transmitted, it is not in the form of a small core but instead in the form of a larger fireball. The color of the HSI-driven flame as it passes is a bright red–yellow, which is different from the blue core produced during initial ignition. The HIS-driven flame surface is accompanied by the consumption of components of marine diesel vapor and the generation of combustion products and heat. As the flame passes to the hot surface, the remaining marine diesel is further ignited.
During the thermal oxygen degradation process, it can be observed that the lightweight components—such as water, trace lipids, and ethers—evaporate first as the hot surface temperature increases.
Figure 5 presents temperature field variation during marine diesel HSI as the temperature of the surface continues to increase. The specified marine diesel has 10 mL leakage, and the heated surface temperature ranges from 753 K to 778 K. The thermocouples are positioned above the center of the high-temperature surface to measure temperatures at different heights ranging from 20 mm to 400 mm above the hot surface.
Figure 5 shows that as the heated surface temperature increases, the temperature field distribution in the vertical space increases accordingly. The temperature field distribution mainly changes between 3.011 s and 9.036 s after the marine diesel contacts the hot surface, as indicated by the temperatures obtained at different times.
Figure 5a indicates the temperature field changes in the vertical space as ignition occurs at 3.011 s. The thermocouple located 20 mm above the high-temperature surface exhibits the fastest response to the temperature increase. As the surface is heated to temperatures from 753 K to 778 K, the air temperature at 20 mm changes from 532.81 K to 570.89 K. This indicates that the increase in air temperature is not the highest closest to the hot surface. When the leaked marine diesel is evaporated by the hot surface, it diffuses to the top of the hot surface to form a certain distribution of concentrated liquid and gas. Meanwhile, due to the temperature transfer generated by thermal diffusion and thermal conduction of high-temperature marine diesel and gas, a certain spatial temperature field is formed. It can be seen from
Figure 5a that because marine diesel is a large-mass molecule, the heat transfer rate is greater than the mass transfer rate in the process of marine diesel evaporation when the leakage is from 5 mL to 20 mL. As a result, the temperature in the vertical space of the hot surface rises very fast. At a vertical height of 200 mm, the temperature goes up from an initial 437.07 K on the hot surface to 529.41 K at the maximum temperature, resulting in a temperature increase of 21.14%. Based on the data from the thermocouple farthest from the hot surface, the temperature varies from an initial temperature of 416.39 K to a maximum of 484.46 K, presenting a range of over 16.34%. In conditions in which the hot surface temperature is below 750 K, the temperature field’s distribution is mainly caused by the joint action of fuel vapor and hot surface heat transfer, as there is no hot surface ignition behavior at 3.011 s. However, when the surface is heated to a temperature over 751 K, the HSI behavior occurs within 3.0 s. Combustion cores and flame transfer phenomena are observed in the vertical space. During flame transfer over a hot surface, the flame continuously heats the combustible gaseous mixture on the surface, maintaining the combustion process in the vertical space. Additionally, the remaining liquid diesel on the hot surface is vaporized by the heat transmitted by the flame. As a result, the temperature of the remaining diesel goes up, causing it to change from a liquid to a gaseous fuel vapor. The vapor that becomes volatile accumulates above the hot surface, leading to incomplete combustion with air due to the consumption of auxiliary gas in the early combustion process. Once the flame spreads to the hot surface, a fuel-rich core with a higher concentration and relatively lower temperature forms at the bottom of the flame. The temperature field distribution shows significant variation in the temperature increase between 200 mm and 300 mm in vertical height, indicating the region of initial HSI. When marine diesel in only a small amount drops on a hot surface, the liquid is decomposed by the action of heat generated by the hot surface to produce reactive free radicals. Under the action of a pressure gradient and buoyance, steam containing free radicals rises and mixes with air to form high-temperature combustible mixed gas. Meanwhile, the concentration of combustible premixed gas gradually accumulates until it reaches the limit of the concentration range of combustion. As the equipment surface is heated to an elevated temperature, droplet evaporation continues. At a vertical height of 200 mm to 300 mm, when the hot surface’s temperature is lower than the critical temperature or premixed gas does not meet the combustion limit, the ignition phenomenon will not occur. Due to the action of air entrainment, the temperature at this height exhibits a declining trend. However, marine diesel continues to evaporate and accumulate, leaving the fire hazard. Once the minimum ignition temperature of the hot surface is reached, the HSI phenomenon may occur. The change in the temperature field over time reveals higher temperatures in the vertical space between 20 mm and 150 mm, leading the direction and path of flame propagation.
3.3. Heat Flux Intensity of the Marine Diesel HSI Process
The marine fuel surface radiant heat flux meter provides data on the thermal feedback of the flame to the fuel surface during flame propagation. The experimental data quantifies the total thermal feedback received by the marine fuels.
Figure 6 shows the data collected by the heat flux meter located in different areas during the HSI of marine fuels. In this experiment, the heat flux meter is arranged sequentially, with the hot surface at the center. The heat flux meter closest to the hot surface is set at a height of 0.3 m, and the horizontal distance from the hot surface center is 0.4 m. The settings of the other heat flux meters are shown in
Figure 6, which illustrates the variation of heat flux intensity with time during the marine diesel HSI. This study analyzes the condition of 15 mL of marine diesel leaked onto a surface with a temperature of 753 K.
Figure 6 shows that as the contact time between the marine diesel and the heated surface rises, the heat flux intensity shows a sudden increase within a certain period; this is followed by decay in different forms. Heat flux intensity begins to increase at 2.198 s in the area closest to the hot surface. During this period, the region distant from the hot surface does not receive heat from the gas-phase medium due to the rapid phase transformation of the fuel on the hot surface. Between 2.6 s and 2.8 s, the heat flux intensity near the hot surface reaches a maximum of approximately 8.42 kW/m
2.
Compared with the initial state, the rate of change in heat flux intensity is significant. Based on observations, HSI of marine diesel occurs within the aforementioned time range, leading to the formation of an initial fire nucleus. During flame propagation, the surrounding area receives radiant heat from the flame. At a distance of 0.7 m and 1.0 m from the center of the hot surface, the heat flux intensities are 5.88 kW/m2 and 1.91 kW/m2, respectively. When the flame spreads to the high-temperature surface, it continues to transfer heat flux to the surrounding space, resulting in stable heat flux intensity for a period of time. As the flame spreads to the hot surface and burns the remaining marine diesel, the heat flux intensity begins to decay. The position farthest from the center of heated surface presents the most significant attenuation of heat flux intensity, dropping to 0.172 kW/m2 at 3.811 s. The results indicate that the heat flux intensity during fuel HSI is transferred to surrounding area. The area in close proximity to the hot surface may receive intense radiant heat, which increases the risk of reigniting leaked fuel in adjacent areas. Due to the properties of marine diesel, the peak value of heat flux intensity is not sustained for a long time during the HSI process. The heat flux intensity begins to decay after the flame propagates to heated surface of the equipment. Observing the thermal performance of HSI-driven flame propagation, the following conclusions can be drawn:
- (a)
It is essential to pay attention to temperature change in the vertical space of equipment surfaces (in which it is easy for a heated surface to form) in order to prevent the accidental ignition of leakage. Safety supervisors in ship engine rooms must carry out monitoring tasks during inspection using infrared thermal imaging techniques.
- (b)
According to results of this section, the heat flow generated by flames is transferred to their surroundings in the horizontal direction, so it is particularly important to protect inflammable and explosive substances near areas of ship engine rooms in which incidents are likely to occur.
3.4. Effect of Hot Surface Temperature on the HSI Delay Time of Marine Diesel
The ignition delay time plays an important role in characterizing the ignition characteristics of marine fuel. Various marine fuels exhibit different ignition characteristics on high-temperature surfaces, and the time required for HSI has discrepancies. There are two basic prerequisites for the local ignition of marine fuel that has leaked onto a hot surface. One is a suitable ratio of fuel vapor/air in the ignition area, which can promote the formation of an initial ignition core. The other is a low-velocity environmental flow field, which can ensure that the initial ignition of the fire core is stable and self-sustaining in the local area, meaning the flame can achieve stable propagation. In a ship engine room, there are more complex influencing factors. Previous research suggests that a vapor layer may be formed between the fuel and the hot surface when the hot surface has a very high temperature. This vapor layer decreases the heat transfer from the heated surface to the fuel, resulting in a decreased evaporation rate of the liquid fuel. HSI delay time is attributed to the interaction between the different fuels and temperatures of the equipment’s surfaces, as shown in Equation (11) [
34].
where
Kos is the frequency factor of the combustion reaction;
cfuel is the molar concentration of combustible material, mol/m
3;
cox is the molar concentration of oxidizer, mol/m
3;
Es is the activation energy, J/mol;
Tr is the temperature of the chemical reaction, K; and
n and
m are the reaction coefficients, respectively.
It is evident that the HSI delay time of combustible fuel is inversely proportional to vapor concentration. The mathematical model of HSI delay time can be deduced using the boiling heat transfer model and vapor plume model.
Figure 7 presents an evolution of ignition delay time for marine diesel at various heated surface temperatures. The data indicate a clear decrease in the ignition time of marine diesel vapor on heated surfaces as the temperature increases. At lower heated surface temperatures, the ignition delay time for diesel is longer, thus increasing the difficulty of ignition, which is consistent with experimental observations. The reason for this is that the heated surface temperature does not reach the critical temperature required to ignite the evaporating vapor of marine diesel. Once the surface is heated to a temperature of 748 K, the ignition delay time of diesel vapor is 3.692 s. A significant amount of marine diesel/air mixture gas accumulates in the vertical space, and the probability of ignition reaches 50%. As the surface continues to be heated, the phase transition of marine diesel absorbs a significant amount of heat, resulting in an increase in the evaporation rate of fuel. This, in turn, reduces the time required to reach critical ignition conditions. The experimental data reveal that the ignition delay time of marine diesel is reduced to 2.048 s when the surface is heated to a temperature around 755 K. If the equipment surface continues to rise up 768 K, the ignition delay time drops to 0.802 s. Compared with the previous scenario, the ignition delay times decrease by 44.5% and 60.8%, respectively. The ignition delay time of marine diesel vapor decreases significantly as the heated surface temperature increases from low to high. However, further increase in the heated surface temperature does not result in a significant change in the ignition delay time. As the surface is heated to a temperature of 773 K, the ignition delay time drops below 0.58 s. After reaching this temperature, the ignition delay time remains around 0.48 s, regardless of how much the heated surface temperature increases.
This study proposes a predictive model to examine the ignition delay time of marine diesel on a hot surface in a ship engine room based on the fitting results of experimental data, as shown in Equation (12). This predictive model could be applied to determine the HSI delay time after marine diesel leakage by analyzing the change in the heated surface temperature. The extended predictive model helps to prevent the further expansion of such incidents by combining experimental data and theoretical mechanisms.
Regarding the HSI delay time and predictive model, the following conclusions can be drawn:
- (a)
The existence of ignition delay time means that the judgment of decision makers can be impaired. Therefore, once marine diesel leakage is found in a ship engine room, it is necessary to focus on the influence of surrounding heat sources on it and to strictly prevent sudden ignition accidents caused by ignition delay on hot surfaces.
- (b)
Using the predictive model proposed in this study, the ignition delay time of marine diesel on hot surfaces can be determined in advance so that initial ignition accidents in ship engine rooms may be prevented.
3.5. Overall Probability of Ignition of Marine Diesel on Hot Surfaces
Marine diesel used in machinery equipment and pipelines is a combustible liquid. If the leakage of marine diesel occurs during ship operation, it is prone to trigger further high-temperature HSI accidents. The HSI of a liquid in contact with a hot surface generally requires a temperature substantially greater than the published auto-ignition temperature of that liquid in handbooks. In this study, the types of marine diesel selected are designed for conducting experiments regarding the probability of HSI of marine diesel in ship engine rooms, and hot surface temperatures are raised to fixed values. The amount of fuel leaked in the current study is fixed from 5 mL to 20 mL. The number of experiments used to determine the HSI temperature of marine fuels is 30. The hot surface simulator can control the heating step of each experiment from 2.0 K to 3.0 K.
Figure 8 shows the distribution of the probability of HSI of marine diesel (with a leakage rate of 10 mL) given different hot surface temperatures; it also reveals the changing trend of HSI probability. It indicates that marine diesel does not always ignite after leaking onto a hot surface. According to the experimental data, when the hot surface temperature is less than 728 K, marine diesel is rarely ignited. The probability of diesel igniting is 3.33% when the heated surface temperature reaches 731 K. Although the HSI probability of marine diesel is low, it is worth noting that the probability of HSI accidents still exists. When the hot surface temperature rises to nearly 748 K, the probability of ignition is 50% (
T50%). As the equipment surface continues to be heated to a temperature of 778 K, the marine diesel will be completely ignited. According to the variation in the HSI probability of marine diesel, the probability of ignition is distributed in three phases. In the initial phase of the temperature range (
Ts < 728 K), marine diesel cannot achieve HSI. When entering the second phase of the temperature range (728 K <
Ts < 778 K), the ignition process shows instability. The temperature range of the hot surface is 778 K <
Ts, which shows that the probability of ignition is always 100%. This means that when the leaking marine diesel flows onto a machinery surface with a hot surface temperature higher than 778 K, ignition behavior inevitably occurs. It is necessary to determine whether the presence of fuel or surrounding items can intensify the risk of a fire source originating in the engine room. By fitting the experimental data, a predictive model for the probability of ignition of leaked marine diesel is proposed, as shown in Equation (13):
where
Ts is the equipment surface temperature, K, and
T50% is the surface temperature when marine diesel has a 50% probability of ignition, K.
Through the results of HSI probability, it is clear that the ignition behavior of marine diesel on hot surfaces has a certain inevitability. In the ship engine room, when the equipment surface is in a high-temperature state for a long time, a safety manager needs to remind inspectors to pay attention to safety issues within this area. Therefore, the HSI probability of accidents occurring can be reduced by increasing attention to risk. In addition, using the model of HSI probability proposed in this study, combined with big data and optimization algorithms, an intelligent predictive model of ignition accidents can be further developed in future. In this way, it can better assist ship managers and decision makers in the safe management of ship operation.