1. Introduction
The Type B LNG tank, classified by the International Maritime Organization, is a widely used independent tank for LNG storage and transportation. Featuring a prismatic shape and supported by a partial secondary barrier, it can be installed on various types of vessels, such as LNG carriers, floating storage and regasification units (FSRUs), and floating LNG (FLNG) platforms. The partial secondary barrier is designed to contain any potential leaks from the primary barrier, ensuring additional safety. The tank is also equipped with internal bulkheads, which help to reduce liquid sloshing, enhancing the structural integrity and stability during transport [
1,
2,
3]. The Type B LNG tank offers several advantages, including high flexibility, safety, reliability, and a low cost. However, it also faces technical challenges, particularly during the pre-cooling process. Pre-cooling is a necessary step before loading LNG into cryogenic equipment, aiming to reduce the container’s temperature to a level close to the LNG boiling point [
4]. This practice helps prevent excessive thermal stress due to rapid temperature decreases and excessive pressure rises due to violent LNG evaporation during initial filling, which can cause safety hazards and operational difficulties. If not properly managed, excessive stress can lead to tank wall crack [
5], followed by LNG leakage, posing significant risks to both the vessel and the surrounding environment, such as explosions and pool fires [
6,
7]. Additionally, such hazards can endanger the safety of crew members working on board, exposing them to potential asphyxiation and frostbite [
8]. Another aspect of LNG tank pre-cooling and other applications leveraging cold energy is the operation efficiency, which involves optimizing the refrigerant usage and cooling speed [
9].
The pre-cooling process typically involves spraying refrigerants into the tank, such as liquid nitrogen, nitrogen gas, LNG, or boil-off gas (BOG) [
10]. Unlike simple heat exchange, pre-cooling with liquid refrigerants involves complex heat and mass transfer through heat convection and phase change. Several factors affect the heat and mass transfer mechanisms during pre-cooling, including the refrigerant type, flow rate, initial temperature, and tank geometry. Most field operators manually adjust the refrigerant valve to maintain the average temperature drop, as measured by scattered thermal monitors, below 5 or 10 degrees per hour. This heuristic approach can lead to a sudden temperature drop, which creates high local temperature gradients and thermal stress. To mitigate this issue and optimize the pre-cooling process, it is essential to understand the heat transfer characteristics of various refrigerants for Type B LNG tanks.
The Type B LNG tank can be broadly categorized as a kind of cryogenic storage tank. The process of pre-cooling these cryogenic liquid storage tanks has been investigated by many researchers primarily using theoretical and numerical methods due to the difficulty of conducting cryogenic experiments. Among the models used to assess the pre-cooling process, theoretical models are the most efficient for predicting dynamic temperature and pressure variations. These models are based on the energy balance of the latent and sensible heat of the cooling medium and the heat loss by the cryogenic tank structure [
11,
12,
13]. The assumption is that the temperature and pressure of the fluid in the tank and the tank structure are uniformly distributed. However, this renders the theoretical models limited, as they cannot account for local pressure and temperature gradients, which can lead to excessive stress on the tank wall. Local temperature gradients can lead to thermal stress that exceeds the yield strength of the structure, as observed in numerous cryogenic applications [
14,
15,
16]. Even if the stress does not exceed yield strength, cyclic pre-cooling operations can cause thermal fatigue. Pressure caused by the evaporation of the refrigerant can also induce local stress concentration [
17,
18], which, while potentially lower than thermal stress, might still be significant [
19]. Three-dimensional numerical modeling, specifically Computational Fluid Dynamics (CFD), allow detailed physical information such as local temperature and pressure during the pre-cooling process to be calculated, potentially offering a higher local accuracy than theoretical models due to the absence of a uniform distribution assumption.
As a result, numerical methods have become increasingly popular with advancements in computational power. Thermal stress induced by local temperature gradients has drawn the most attention in past studies. Kang et al. [
20,
21] assessed the cooldown characteristics and distribution of thermal stress during the filling of large-scale cryogenic tanks, finding that thermal stress on the tank wall is proportional to the cooling rate and peaks as the liquid nitrogen passes by. Lu et al. [
22] conducted an LNG liquid phase pre-cooling simulation of a membrane tank and found that the temperature gradient is greatest on the innermost layer of the tank, indicating that thermal stress should be given sufficient attention to avoid tank wall damage. Chen et al. [
23] performed an LNG liquid-phase pre-cooling analysis of the re-condenser in an LNG receiving station, discovering that the fatigue life of the nozzle, where thermal stress is maximum, cannot meet fatigue life requirements if the temperature drop rate exceeds 10 °C/s.
Recent years have also seen growing interest in the optimization of the pre-cooling process, with the aim of improving safety and cost efficiency. Shin et al. [
24] introduced a model predictive control system that uses a reduced-order model and optimization approach to regulate temperature drop and pressure buildup during the pre-cooling of a membrane-type LNG tank. Lim et al. [
25] optimized the injection flow rate of cryogenic CO
2 gas used for pre-cooling CO
2 storage tanks on ships to minimize CO
2 loss. Kulitsa and Wood [
26] proposed spraying LNG into the membrane LNG tank at very low rates so that the boil-off gas (BOG) generated during the pre-cooling can be utilized by the ship’s engine system instead of being wasted.
However, most existing studies on the pre-cooling process of Type B LNG tanks and other cryogenic tanks have tested only one refrigerant, lacking a comprehensive comparison of the pre-cooling properties of different refrigerants. Therefore, there is a need for a systematic investigation into the pre-cooling process of the Type B LNG tank using various refrigerants, which can provide guidance for the optimal selection and design of the pre-cooling operation. In this paper, we present a CFD assessment of the pre-cooling process of Type B LNG tanks using several refrigerants: liquid nitrogen, nitrogen gas, LNG, BOG, and their combinations. The simulation model considers the effects of phase change, convective heat transfer, and conjugate heat exchange between the fluid in the tank and the tank structure. The simulation results show the temporal variation and distribution of temperature, as well as the heat transfer coefficient at the tank wall. The results illustrate how each refrigerant differs in pre-cooling properties, such as cooling efficiency and the threat to structural integrity. This paper provides guidance for selecting the optimal refrigerant and optimizing the pre-cooling operation for Type B LNG tanks and similar cryogenic storage tanks.
3. Results and Discussion
In this section, the results of the transient LNG B-type tank pre-cooling simulations performed using the numerical model described in
Section 2 are outlined. Six pre-cooling schemes using liquid nitrogen (LN), nitrogen gas (NG), NG + LN, LNG, BOG, and BOG + LNG are assessed. The mixed pre-cooling schemes (Case 3 and Case 6) cool the tank first with cryogenic gas before injecting the liquid refrigerant to continue the pre-cooling, a practice frequently used in LNG receiving stations. In the current simulation, the solutions in Case 2 (NG) and Case 5 (BOG) at 5.8 h were used as the initial conditions for Case 3 and Case 6, respectively, while the inlet liquid refrigerant flow velocity was kept the same as in Case 1 (LN) and Case 4 (LNG).
The flow rate and temperature of the refrigerants are listed in
Table 2. The flow rates were selected to allow the pre-cooling of all schemes to complete within roughly the same duration (when the temperature stops varying significantly), with the gas refrigerants and liquid refrigerants each maintaining consistent flow rates. In the following subsections, temperature variation, temperature distribution, and temperature gradient are compared to evaluate the different cooling schemes.
3.1. Temperature Variation
The variation in the average temperature of the inner tank layer and the bulkhead is shown in
Figure 6 and
Figure 7, respectively. The trend suggests that the cooldown rate of the tank structure gradually decreases as pre-cooling proceeds, with the bulkhead cooling more quickly than the inner layer. This can be attributed to the decreased heat flux resulting from the reduced temperature difference between the tank structure and the refrigerant as the pre-cooling progresses. Meanwhile, the variation in the average temperature of the fluid inside the tank, shown in
Figure 8, decreases at a much faster rate, reaching a steady state more quickly. This rapid fluid cooling contributes to the faster cooling rate of the bulkhead, which is immersed in the refrigerant, unlike the inner tank layer.
Liquid refrigerants, specifically liquid nitrogen (LN) in Case 1 and liquefied natural gas (LNG) in Case 4, demonstrate a superior cooling performance compared to gas refrigerants with the same chemical component. Case 1 (LN) achieves the fastest initial temperature drop of the inner tank layer, reaching below −100 °C within 8 h and stabilizing around −140 °C by 16 h. In comparison, Case 4 (LNG) shows a rapid initial cooldown similar to Case 1 but levels off at −120 °C after 10 h. The normalized mean temperature drop rate at each sampled time (calculated by trapezoidal rule) until 8 h is listed in
Table 3, which suggests that the cooling rate of LNG is only 79% of that of LN. The difference in performance between LN and LNG can be attributed to their distinct thermophysical properties. Liquid nitrogen, with its extremely low boiling point of −195 °C, provides a very high temperature difference between the refrigerant and the tank, driving faster and more substantial heat transfer. This results in a more rapid initial cooldown and a lower final temperature. In contrast, LNG, with a boiling point of around −162 °C, has a higher initial temperature compared to LN. While still effective, the smaller temperature difference between LNG and the tank results in a slower heat transfer rate, leading to a higher final temperature compared to LN.
For gas refrigerants, Case 2 (NG) exhibits a slower cooldown rate compared to liquid nitrogen, reaching only about −115 °C after 16 h. This slower cooling rate can be attributed to the higher initial temperature and absence of latent heat, which has a significant contribution to liquid pre-cooling process. Case 5 (BOG) shows a similar cooling pattern to NG but with 79.4% of its cooling rate (as calculated from results in
Table 3). The performance difference between NG and BOG can also be explained by the lower initial temperature of BOG, which leads to a smaller heat transfer rate on the tank surface.
The gas–liquid mixed cooling schemes, Case 3 (NG + LN) and Case 6 (BOG + LNG) reach steady state at approximately the same time as the pure liquid pre-cooling, despite the difference in temperature as the liquid phase cooling starts. Overall, liquid nitrogen and LNG facilitated the most rapid cooldown. The mixed schemes, while slightly slower, achieved comparable cooling rates. The gas-only refrigerants exhibited the slowest cooling rates and the highest steady-state temperatures among the cases investigated.
To further compare the efficiency of different pre-cooling schemes, the average heat transfer coefficients (HTC) on the inner tank surface were obtained and are shown in
Figure 9, and their mean values over the simulation period are listed in
Table 4. HTC is an important indicator of the effectiveness of heat transfer between the refrigerant and the tank surface. The liquid refrigerants (Case 1 and Case 4) exhibited higher HTCs than their pure gas counterparts (Case 2 and Case 5). The HTCs continued to increase, albeit at a slower rate, as cooling progressed. In contrast, the gas refrigerants (Case 2 and Case 5) reached a plateau before decreasing, showing 8% and 7% lower mean HTCs (calculated from
Table 4). This difference arises from the distinct cooling mechanisms of the two types of refrigerants. Cryogenic liquids benefit from latent heat during evaporation and both forced and natural convection of the evaporated gas and liquid. As evaporation diminishes due to lower tank temperatures, the remaining liquid still significantly contributes to heat transfer. Gas refrigerants, on the other hand, rely solely on sensible heat through forced and natural convection. Liquids achieve much larger heat transfer coefficients than gasses due to their higher thermal conductivity, leading to the continuous rise in the heat transfer coefficient for liquid refrigerants. Despite having greater thermal conductivity, LNG shows a 16% lower heat transfer coefficient than liquid nitrogen at the same inlet velocity (calculated from
Table 4). This is due to the lower temperature of its evaporated gas (i.e., BOG), which leads to weaker natural heat convection, primarily driven by the temperature gradient at the tank wall. The same rationale can account for the relative difference between the two gas refrigerants (Case 2 (NG) and Case 5 (BOG)). For the mixed strategies (Case 3 and Case 6), the heat transfer coefficient quickly matches the level of pure liquid refrigerants, aligning with the temperature variation trends shown in
Figure 6. As indicated in
Table 4, the mean heat transfer coefficient is also approximately the same as that of pure liquid refrigerants.
3.2. Temperature Distribution
To examine the temperature distribution within the tank during the pre-cooling process, temperature contour plots of the inner tank surface and bulkhead for different cooling schemes at 1 h into the pre-cooling process were extracted and are shown in
Figure 10. For all the cases assessed, the temperature distribution of the bulkhead is more uniform than that of the inner tank surface, which exhibits lower surface temperatures above the center transverse bulkhead than below it, creating a vertical temperature gradient on the side surface. The temperature of the tank’s bottom surface drops faster than that of the side surfaces, with the top surface temperature remaining the highest.
The pattern of a non-uniform temperature distribution on the inner tank surface can be explained by the velocity magnitude contour and streamline plot sampled at a YZ plane (shown in
Figure 11). The streamlines in
Figure 12 suggest that the refrigerants injected into the tank reach the bottom first, forming two nearly symmetric vortices that circulate the flow upwards along the side surfaces. As a result, the refrigerant stream arriving at the side surfaces has already lost some of its cold energy, causing slower cooling. Additionally, the center transverse bulkhead exerts a blocking effect on the inlet refrigerant flow, effectively dividing the tank into two flow regions, upper and lower. This observation helps to explain the uneven vertical temperature distribution on the inner tank surface. It is worth noting that the bulkhead, used to prevent liquid sloshing, is a key feature distinguishing the Type B LNG tank from other LNG tanks [
36,
37]. As a result, the Type B LNG tanks may be more prone to an uneven temperature distribution (and therefore thermal stress problems) compared to the other types of LNG tanks.
3.3. Temperature Gradient
In light of the uneven temperature distribution within the tank structure, we have extracted the transient variation of the average temperature gradient magnitude of the inner tank layer and bulkhead, which is an important indicator of the level of thermal stress during the pre-cooling process.
Figure 13 and
Figure 14 show a comparison of the temperature gradients for the different pre-cooling schemes. It is evident that the temperature gradient of the bulkhead is significantly smaller than that of the inner tank layer, which can be attributed to the more uniform cooling observed in
Figure 10.
For different pre-cooling schemes, the relative magnitudes of the temperature gradient are consistent between the inner tank layer and bulkhead (see
Table 5 and
Table 6) and aligned with the structure cooling rate presented in
Section 3.1. Among all the schemes, Case 5 (BOG) exhibits the lowest overall temperature gradient, followed by Case 4 (LNG) and Case 6 (BOG + LNG). The schemes using liquid nitrogen and nitrogen gas result in larger temperature gradients, with Case 1 (LN) and Case 3 (NG + LN) showing higher gradients than Case 2 (NG). It is noteworthy that the mixed cooling scheme can achieve a comparable cooling rate (as illustrated in
Section 3.1) while maintaining a smaller, if not similar, temperature gradient compared to the pure liquid cooling scheme, provided the inlet velocity of the refrigerants is properly controlled, as demonstrated in this study.
4. Conclusions
CFD simulations of the Type B LNG tank pre-cooling process using different refrigerants provide valuable insights into the thermodynamic characteristics of each cooling strategy. During the pre-cooling process, the temperature of the tank rapidly decreases before leveling off. The fluid inside the tank exhibits the fastest temperature drop, followed by the bulkheads immersed in it, while the tank hull structure shows the slowest temperature drop. Liquid nitrogen (LN) was found to be the most efficient refrigerant, achieving the highest cooling rate through both latent and sensible heat. LNG also demonstrated a relatively high cooling rate, 79% of that of the liquid nitrogen due to its higher boiling point. In contrast, gas-only pre-cooling schemes relying solely on sensible heat are less efficient, among which the cooling rate of the boil-off gas (BOG) was 79.4% of that of nitrogen gas (NG). Mixed refrigerants such as NG + LN and BOG + LNG can achieve comparable, while slightly slower, cooling rates than the pure liquid refrigerants, outperforming gas-only strategies. A further assessment of heat transfer coefficient suggests the mixed cooling schemes have almost identical heat transfer coefficient on the inner tank surface to the liquid cooling scheme, over 5% higher than the gas refrigerants.
The study also highlighted the uneven temperature distribution within the tank, particularly in the inner tank layer, due to the bulkhead’s blockage effect on refrigerant flow. This uneven distribution, characterized by high local temperature gradients, can induce significant thermal stress, potentially compromising the structural integrity of the tank. Given that the bulkhead is a key feature distinguishing the LNG B-type tank from other LNG tanks, this type of tank may be more prone to excessive thermal stress. Among the cooling schemes assessed in this study, mixed schemes exhibit higher thermal gradient than the gas schemes but lower thermal gradient than the liquid schemes while achieving the same level of high cooling efficiency if the inlet velocities of gas and liquid refrigerants are configured properly.
In summary, the selection of refrigerants for the pre-cooling process of Type B LNG tanks should consider both cooling efficiency and the potential for thermal stress. The findings in this work can hopefully provide a basis for developing more effective and safer pre-cooling procedures for Type B LNG tanks and similar cryogenic tanks. Future research could focus on evaluating additional refrigerants, quantifying the magnitude of thermal stress under different refrigerants, and validating the results through full-scale experiments.