AppliedThermalEngineering149(2019)344–353

Contents lists available at

Applied Thermal Engineering

journal homepage:

Research Paper

Improving design and operation at LNG regasification terminals through a

corrected storage tank model
Cải thiện thiết kế và vận hành tại các cảng tái khí hóa LNG thông qua một
sửa chữa mô hình bể chứa
⁎
Mohd Shariq Khana, , Surya Effendyb, I.A. Karimib, Aref Wazwaza
a
DEPARTMENT of CHEMICAL Engineering, College of Engineering, DHOFAR University, SALALAH 211, OMAN
b
DEPARTMENT of CHEMICAL & BIMOLECULAR Engineering, NATIONAL University of SINGAPORE, 117585, SINGAPORE

HIGHLIGHTS

• •Planned/unplanned
LNG tank aspect ratio
tankofpressurization
1 minimizes the boil-off gas
decreases withand average liquid
increasing compressor
level.

• Minimum recirculation that prevents 2-phase is the best operating

• Unique minimum recirculation for a given pipe dia/length/heat transfer coefficient.

NỔI BẬT

• • Tỷ lệ khung hình của bồn chứa LNG là 1 giảm thiểu lượng khí sôi và công suất máy nén trung bình.

• • Áp suất bể theo kế hoạch/không theo kế hoạch giảm khi mức chất lỏng tăng.

• • Tuần hoàn tối thiểu ngăn 2 pha là chính sách vận hành tốt nhất.

• • Tuần hoàn tối thiểu duy nhất cho một đường kính/chiều dài/hệ số truyền nhiệt nhất định

A R T I C L E I N F O BÀI
VIẾT A B S T R A C T TRỪU TƯỢNG

Keywords: The simulation of liquefied natural gas (LNG) storage tanks is often based on several problematic
LNG storage tank assumptions, for instance, estimation of boil-off gas (BOG) generation using boil-off rate, vapor-liquid
LNG regasification terminal equilibrium in the tank, static liquid level, the use of only lateral area for heat loss calculations, and etcetera.
Minimum recirculation Some of these assumptions are built into selected commercial simulators, creating further challenges in
Aspect ratio simulating the behavior of LNG tanks. The present study highlights these challenges in the context of a
Planned/unplanned shutdown
commonly used process simulator, Aspen HYSYS, and provides analytical and intuitive solutions to those
từ khóa:
problems. The resultant model is validated against an established first-principle model and then exploited for
bể chứa LNG
finding improved LNG regasification terminal design and operation strategies. Tank aspect ratio (AR) was
Thiết bị đầu cuối tái khí hóa
studied in relation to plant capacity, recirculation rate, and recirculation line length. An aspect ratio of 1
LNG Tuần hoàn tối thiểu Tỷ lệ
consistently results in minimum BOG generation, in contrast to the value of 0.5 frequently cited in the
khung hình Tắt máy theo kế
hoạch/không theo kế hoạch literature. During the planned/unplanned shutdown of re- gasification terminal, higher liquid level in the tank
decreases evaporation and thus BOG generation. Minimum recirculation rate that prevents 2-phase flow is
found to minimize BOG generation and compressor duty.
Việc mô phỏng bể chứa khí tự nhiên hóa lỏng (LNG) thường dựa trên một số giả định có vấn đề, ví dụ, ước tính
lượng khí sinh ra khi sôi (BOG) bằng cách sử dụng tốc độ sôi, cân bằng hơi-lỏng trong bể, mức chất lỏng tĩnh, việc
chỉ sử dụng diện tích bên để tính toán tổn thất nhiệt, v.v. Một số giả định này được tích hợp vào các mô phỏng
thương mại được chọn, tạo ra nhiều thách thức hơn trong việc mô phỏng hoạt động của các bể chứa LNG. Nghiên
cứu hiện tại nhấn mạnh những thách thức này trong bối cảnh của một trình giả lập quy trình thường được sử dụng,
Aspen HYSYS, và cung cấp các giải pháp phân tích và trực quan cho những vấn đề đó. Mô hình kết quả được xác
thực dựa trên mô hình nguyên tắc đầu tiên đã được thiết lập và sau đó được khai thác để tìm kiếm các chiến lược
vận hành và thiết kế thiết bị đầu cuối tái khí hóa LNG cải tiến. Tỷ lệ khung hình bể (AR) đã được nghiên cứu liên quan
đến công suất nhà máy, tốc độ tuần hoàn và chiều dài đường tuần hoàn. Tỷ lệ khung hình 1 luôn dẫn đến việc tạo
BOG tối thiểu, trái ngược với giá trị 0,5 thường được trích dẫn trong tài liệu. Trong quá trình đóng cửa thiết bị đầu
cuối tái khí hóa theo kế hoạch/ngoài kế hoạch, mức chất lỏng cao hơn trong bể làm giảm sự bay hơi và do đó tạo ra
BOG. Tỷ lệ tuần hoàn tối thiểu ngăn chặn dòng chảy 2 pha được tìm thấy để giảm thiểu nhiệm vụ tạo ra BOG và máy
nén.
1. Introduction
and requires several pre-processing step before it can be used. In
At present, natural gas (NG) accounts for roughly a quarter of par- ticular, the transport of NG from one location to another
global primary energy consumption [1]. This dominant position is typically en- tails liquefaction [5], which reduces the volume of NG
achieved by virtue of its clean-burning properties, as well as an by several orders of magnitude. After transport, the liquefied
international drive towards increased energy security amidst natural gas (LNG) needs to be regasified prior to use. Liquefaction
increasing demand for fuel. NG consumption is expected to increase [6,7] and regasification [8] are two of the most energy-intensive
by 1.6% per annum between 2016 and 2034 [2], outstripping the operations in the NG value chain, and much effort has been spent
growth in oil and coal consump- tion, and eventually accounting for the to improve its energy efficiency via opti- mization [9–11] and
majority of primary energy consumption by 2050 [3]. The prospect for process integration [12,13].
NG seems bright in the near future, given its status as a bridge fuel Liquid natural gas (LNG) storage tanks play a central role in both
towards renewable energy [4]. This motivates research into the various liquefaction and regasification. However, the design and operation of
aspects of the NG value chain, aimed at increasing its energy LNG storage tanks are often based on several problematic
efficiency and profitability. Part of this research concerns the fact that assumptions. For example, the most common method for estimating
NG is often found at remote locations, boil-off gas generation, and thus for designing relevant upstream and
downstream equipment, is the boil-off rate, a parameter ranging from
0.02 to 0.05%
1. Giới thiệu
Hiện nay, khí đốt tự nhiên (NG) chiếm khoảng một phần tư mức tiêu thụ năng lượng sơ cấp toàn cầu [1]. Vị trí thống trị
này đạt được nhờ các đặc tính đốt sạch của nó, cũng như nỗ lực quốc tế hướng tới tăng cường an ninh năng lượng
trong bối cảnh nhu cầu nhiên liệu ngày càng tăng. Tiêu thụ NG dự kiến sẽ tăng 1,6% mỗi năm từ năm 2016 đến năm
2034 [2], vượt xa tốc độ tăng tiêu thụ dầu và than, và cuối cùng chiếm phần lớn tiêu thụ năng lượng sơ cấp vào năm
2050 [3]. Triển vọng cho NG có vẻ tươi sáng trong tương lai gần, với vị thế là nhiên liệu cầu nối hướng tới năng lượng
tái tạo [4]. Điều này thúc đẩy nghiên cứu về các khía cạnh khác nhau của chuỗi giá trị NG, nhằm tăng hiệu quả sử dụng
năng lượng và lợi nhuận. Một phần của nghiên cứu này liên quan đến thực tế là NG thường được tìm thấy ở các địa
điểm xa xôi,

và yêu cầu một số bước tiền xử lý trước khi có thể sử dụng. Đặc biệt, việc vận chuyển NG từ vị trí này sang vị trí khác
thường kéo theo quá trình hóa lỏng [5], làm giảm thể tích NG xuống vài bậc độ lớn. Sau khi vận chuyển, khí tự nhiên
hóa lỏng (LNG) cần được khí hóa lại trước khi sử dụng. Hóa lỏng [6,7] và tái hóa khí [8] là hai trong số các hoạt động
sử dụng nhiều năng lượng nhất trong chuỗi giá trị NG và nhiều nỗ lực đã được dành để cải thiện hiệu quả sử dụng
năng lượng thông qua tối ưu hóa [9–11] và tích hợp quy trình [ 12,13].
Bể chứa khí thiên nhiên lỏng (LNG) đóng vai trò trung tâm trong cả quá trình hóa lỏng và tái khí hóa. Tuy nhiên, việc
thiết kế và vận hành các bể chứa LNG thường dựa trên một số giả định có vấn đề. Ví dụ, phương pháp phổ biến nhất
để ước tính lượng khí sinh ra từ quá trình sôi, và do đó để thiết kế thiết bị đầu nguồn và hạ lưu có liên quan, là tốc độ
sôi, một tham số nằm trong khoảng từ 0,02 đến 0,05%.
⁎
Corresponding author.
E-MAIL ADDRESS: [email protected] (M.S. Khan).

https://doi.org/10.1016/j.applthermaleng.2018.12.060
Received 2 August 2018; Received in revised form 27 October 2018; Accepted 12 December 2018
Availableonline13December2018
1359-4311/©2018ElsevierLtd.Allrightsreserved.
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

Nomenclature mbarg QB millibar

QS,V (gauge)
QT heat transfer from the bottom surface of the tank
AR BOG CFD
aspect
HPratio
HTCboil-off
LNG gas TL heat transfer from the shell surface to the vapor phase of
computational fluid dynamics high pressure t/h TA TV UtheUtank
heat UB
transfer from the top surface of the tank liquid temperature
heat transfer coefficient liquefied natural gas US,V+T tonnes/hour
US,V UO,S US
ambient temperature vapor-phase temperature
USGPM interface heat-transfer coefficient pipe heat transfer coefficient
bottom-surface heat-transfer coefficient
composite shell- and top-surface heat-transfer coefficient, weighted a
LL liquid level shell-surface heat-transfer coefficient outer-to-shell heat-transfer coef
MMTPA million metric tonne per annum NGnatural gas shell-to-inner-vapor heat transfer coefficient US gallon per minute

Rpm Cp,V DT
rotations
HT per minute vapor-phase heat capacity tank diameter
HL height of the tank
HE height of the liquid in tank enthalpy of vaporization evaporation rate
mE

[14]. This assumption seems unusual; for example, one would expect nhiên liệu xe buýt đối với tốc độ tạo ra BOG và sự phát triển của thành
the amount of boil-off gas to be proportional to the heat-transfer area phần LNG tại các trạm tiếp nhiên liệu LNG. Họ kết luận rằng có thể
of the tank, as opposed to the volume of the tank, as is inherent in the giảm tổng lượng nhiên liệu thất thoát bằng cách tăng tần suất tiếp
definition of the boil-off rate. This assumption is present even in rela- nhiên liệu cho xe buýt. Sangeun Roh và cộng sự. [18] đã nghiên cứu
tively recent publications. For example, the model developed by Miana sự đối lưu tự nhiên thoáng qua trong bể chứa LNG được điều áp thông
[15] requires some constant boil-off rate parameter to estimate the qua động lực học chất lỏng tính toán (CFD). Việc tạo BOG được đánh
amount of boil-off gas generated, which closes the mass and energy giá bằng cách sử dụng cân bằng năng lượng giao diện hơi-lỏng. Ảnh
balances, and thus allows the tank to be simulated. hưởng của kích thước bể, áp suất bể và cấu hình điều áp đối với việc
[14]. Giả định này có vẻ bất thường; ví dụ, người ta mong đợi lượng tạo ra BOG sau đó được định lượng. Phân tích được trình bày trong đó
khí bay hơi tỷ lệ thuận với diện tích truyền nhiệt của bể, trái ngược với được giới hạn ở mức chất lỏng tĩnh và áp suất bể gần với áp suất gần
thể tích của bể, như vốn có trong định nghĩa về tốc độ bay hơi. Giả khí quyển của mức chất lỏng trong bể LNG thông thường được giả
định này hiện diện ngay cả trong các ấn phẩm tương đối gần đây. Ví định là tĩnh. Dimopoulos và Frangopoulos [19] đã nghiên cứu động lực
dụ, mô hình được phát triển bởi Miana học của sự bay hơi LNG trong quá trình vận chuyển bằng đường biển,
[15] yêu cầu một số tham số tốc độ sôi không đổi để ước tính lượng giả sử các pha hơi-lỏng lý tưởng và trạng thái cân bằng hơi-lỏng.
khí sôi được tạo ra, đóng cân bằng khối lượng và năng lượng, và do Several studies discussed the BOG recondensation process without
đó cho phép mô phỏng bể. calculating the amount of BOG produced from the tank. Yuan et al. [20]
The dynamics of LNG regasification terminals and LNG value for instance assumed that the average daily evaporation from LNG tank
chains have been extensively studied at various levels of granularity. varies between 0.03% and 0.1% hence the BOG flowrates varies from 2
For ex- ample, Lim et al. [16] used Aspen HYSYS, a commercial t/h to 10 t/h. Li et al. [21] optimized BOG recondensation process by
process si- mulator, to simulate unstable two-phase flows in LNG assuming 6.69 t/h of BOG flow as the designed capacity of China’s Da-
pipes. They proposed a methodology for enforcing stability by peng LNG terminal. Apparently, their optimization results seem case
replacing the unit module at the point where instability is first specific and could not be generalized for a tank loading operation when
generated. Chen et al., [17] studied the feasibility of LNG-fuelled public the spikes in BOG generation occur due to volumetric displacement,
transportation (in parti- cular buses) by modeling the effect of bus transfer pipe and pumping heat gain and flashing due to the pressure
refueling frequency on the rate of BOG generation and evolution of differential between ship and storage tank.
LNG composition at LNG re- fueling stations. They concluded that total Một số nghiên cứu đã thảo luận về quá trình tái ngưng tụ BOG mà
fuel loss can be reduced by increasing bus refueling frequency. không cần
Sangeun Roh et al. [18] in- vestigated transient natural convection in a tính lượng BOG sinh ra từ bể. Nguyên và cộng sự. [20] chẳng hạn, giả
pressurized LNG storage tank via computational fluid dynamics (CFD). định rằng lượng bay hơi trung bình hàng ngày từ bể chứa LNG thay đổi
The BOG generation is evaluated using a vapour-liquid interface trong khoảng 0,03% đến 0,1% do đó lưu lượng BOG thay đổi từ 2 t/h
energy balance. The effect of tank size, tank pressure and đến 10 t/h. Lý và cộng sự. [21] quy trình tái ngưng tụ BOG được tối ưu
pressurization profile on BOG generation is then quantified. The hóa bằng cách giả sử lưu lượng BOG là 6,69 t/h theo công suất thiết kế
analysis presented therein is limited to static liquid level and tank của nhà ga Dapeng LNG của Trung Quốc. Rõ ràng, kết quả tối ưu hóa
pressures which are close to the near-atmospheric pressures of the của chúng có vẻ cụ thể theo từng trường hợp và không thể khái quát
conventional LNG tanks liquid level are assumed static. Dimopoulos hóa cho hoạt động chất hàng vào két khi sự tăng đột biến trong quá
and Frangopoulos [19] studied the dynamics of LNG eva- poration trình tạo BOG xảy ra do sự dịch chuyển thể tích, đường ống chuyển tải
during marine transportation, assuming ideal vapour-liquid phases, và bơm tăng nhiệt và nhấp nháy do chênh lệch áp suất giữa tàu và két
and vapour-liquid equilibrium. chứa.
Động lực của các thiết bị đầu cuối tái khí hóa LNG và chuỗi giá trị LNG A relevant issue is that studies on LNG storage tanks are typically
đã được nghiên cứu rộng rãi ở các mức độ chi tiết khác nhau. Ví dụ, context-free i.e. they do not consider the effect of upstream and
Lim et al. [16] đã sử dụng Aspen HYSYS, một bộ mô phỏng quy trình downstream units, and constrained to the holding mode [22] i.e. the
thương mại, để mô phỏng dòng chảy hai pha không ổn định trong normal operation, as opposed to the loading and unloading modes,
đường ống LNG. Họ đã đề xuất một phương pháp để thực thi sự ổn wherein the tank is transferring LNG to and from marine carriers,
định bằng cách thay thế mô-đun đơn vị tại điểm mà sự mất ổn định Một vấn đề liên quan là các nghiên cứu về bể chứa LNG thường không
được tạo ra lần đầu tiên. Chen và cộng sự, [17] đã nghiên cứu tính có ngữ cảnh, tức là chúng không xem xét tác động của các thiết bị
khả thi của phương tiện giao thông công cộng sử dụng nhiên liệu LNG thượng nguồn và hạ nguồn, và bị hạn chế ở chế độ giữ [22] tức là hoạt
(đặc biệt là xe buýt) bằng cách mô hình hóa tác động của tần suất tiếp động bình thường, trái ngược với chế độ nạp và dỡ tải , trong đó bồn

34
 M.S. KHAN et AL. AppliedThermAlEngineeri
đang chuyển LNG đến và từ các hãng vận tải biển, tương ứng. Đây là dỡ hàng và tuần hoàn, nhưng các ng149(2019)344–353
chi tiết liên quan đến hoạt động của
trường hợp, ví dụ, trong mô hình thời tiết được phát triển bởi Migliore bể chứa lại không có. Thay cho bể chứa LNG, một dòng LNG không
[23]. Jourda và Probert [24] đã xem xét một bể chứa LNG trên mặt đất đổi duy nhất và một dòng khí sôi không đổi tương ứng được giả định.
với mức chất lỏng tĩnh và trạng thái cân bằng hơi-lỏng giả định để tính Đây là một vấn đề vì lượng khí sôi dự kiến sẽ tăng theo tốc độ tuần
tốc độ truyền nhiệt vào bể chứa thông qua dẫn nhiệt, đối lưu và bức hoàn. Do đó, tốc độ tuần hoàn tối ưu được trình bày trong đó, được
xạ. Mô hình này sau đó được sử dụng để nghiên cứu tỷ lệ chiều cao giả định là không đổi, có thể không chính xác. Trong một ấn phẩm gần
trên đường kính tối ưu của bể chứa LNG, được tìm thấy là ca. 0,5 bất đây, Effendy et al. [26] đã xem xét động lực học của cả đường tuần
kể giá trị của các thông số truyền nhiệt khác nhau. Việc tách các bể hoàn và bể chứa LNG. Nghiên cứu này nắm bắt một phần đáng kể
chứa khỏi các đơn vị thượng nguồn và hạ nguồn là hợp lý. Các bể hoạt động của bể chứa đẳng áp khi có đường tuần hoàn nhưng sử
chứa LNG là thiết bị cao nhất trong chế độ chứa và dự kiến sẽ không dụng hành vi chất lỏng đơn giản, ví dụ: định luật khí lý tưởng trong
bị ảnh hưởng mạnh bởi hoạt động của các thiết bị phía sau nói chung. việc tính toán một số tính chất chất lỏng. Vấn đề thứ hai liên quan đến
Ngoại lệ đối với quy tắc này là dòng tuần hoàn, cần thiết để làm mát kết quả cho thấy nhiệt độ thoát ra của khí sôi có thể cao tới 16 °C
các đường vận chuyển LNG khi kho cảng đang vận hành ở chế độ trên nhiệt độ LNG. Một nghiên cứu gần đây hơn của Migliore et al. [27]
chờ. Đường tuần hoàn trải qua quá trình truyền nhiệt đáng kể với môi gợi ý rằng nhiệt độ thoát ra tối đa phải gần hơn với 8 °C và chúng tôi
trường trong chế độ giữ và truyền nhiệt này vào bể chứa. thường đồng ý với kết quả này, ngoại trừ giả thiết về pha hơi đọng lại,
điều không thể xảy ra khi có sự truyền nhiệt từ bề mặt vỏ của bể chứa.
respectively. This is the case, for example, in the weathering model
developed by Migliore [23]. Jourda and Probert [24] considered an
above-ground LNG storage tank with static liquid level and assumed
vapor-liquid equilibrium to calculate the rate of heat ingress into the
storage tank via conduction, convection, and radiation. The model is
subsequently used to study the optimal height-to-diameter ratio of LNG
storage tanks, which is found to be ca. 0.5 regardless of the value of
the various heat-transfer parameters. The uncoupling of storage tanks
from upstream and downstream units is sensible. LNG storage tanks
are the uppermost unit during holding mode and are not expected to be
strongly affected by the operation of downstream units in general. The
exception to this rule is the recirculation flow, which is needed to cool
down the LNG transfer lines when the terminal is operating in the
holding mode. The recirculation line undergoes substantial heat transfer
with the environment during the holding mode and transfers this heat to
the storage tank.
The recirculation flow has been studied in the literature, but very
rarely in conjunction with a detailed LNG storage tank model. For ex-
ample, Park [25] used a process simulator to study the recirculation
operation, with a particular focus on the operating costs of the
BOG/HP compressor and the in-tank pump. The authors suggest that
there exist a trade-off between the operating costs of these units,
which may be minimized to give an optimal recirculation flow. The
authors further suggest that the optimal recirculation flow depends on
the amount of send-out LNG. While the study is capable of capturing
the overall pic- ture of the loading, unloading, and recirculation lines,
the details re- garding the storage tank operation are absent. In lieu of
the LNG storage tank, a single constant LNG stream and a
corresponding constant boil- gas stream is assumed. This is
problematic since the amount of boil-off gas is expected to increase
with recirculation rate. Consequently, the optimal recirculation rate
presented therein, which is assumed to be constant may be incorrect.
In a recent publication, Effendy et al. [26] considered the dynamics of
both the recirculation line and the LNG storage tank. This study
captures a substantial portion of the behavior of an isobaric storage
tank in the presence of a recirculation line but uses simplistic fluid
behavior e.g. ideal gas law in the calculation of several fluid
properties. A second problem concerns the result sug- gesting that
the boil-off gas exit temperature can be as high as 16 °C
above the LNG temperature. A more recent study by Migliore et al. [27]
suggests that the maximum exit temperature should be closer to 8 °C,
and we generally concur with the result, with the exception of the as-
sumption of a stagnant vapor-phase, which is impossible in the
presence of heat transfer from the shell surface of the storage tank.
Dòng tuần hoàn đã được nghiên cứu trong tài liệu, nhưng rất hiếm
khi kết hợp với mô hình bể chứa LNG chi tiết. Ví dụ, Park [25] đã sử
dụng trình mô phỏng quy trình để nghiên cứu hoạt động tuần hoàn,
đặc biệt tập trung vào chi phí vận hành của máy nén BOG/HP và máy
bơm trong bể. Các tác giả cho rằng tồn tại sự cân bằng giữa chi phí
vận hành của các đơn vị này, có thể được giảm thiểu để tạo ra dòng
tuần hoàn tối ưu. Các tác giả đề xuất thêm rằng lưu lượng tuần hoàn
tối ưu phụ thuộc vào lượng LNG gửi đi. Mặc dù nghiên cứu có khả
năng nắm bắt được bức tranh tổng thể về các dây chuyền bốc xếp,

34
M.S. KHAN et AL. AppliedThermAlEngineeri
Section 2, we review LNG storageng149(2019)344–353
tank design and a ty- pical
regasification termination. In Section 3, the modeling of an LNG
The dynamic behavior of an LNG storage tank has been studied storage tank in Aspen HYSYS is considered. The tank module in
in detail using computational fluid dynamics (CFD) [22,28,29]. Sale et Aspen HYSYS possesses several challenges, which are addressed via
al. [22] concluded that surface evaporation is the governing boiling some simple modifications. The correctness of these modifications is
mechanism and nucleate boiling is unlikely in an LNG storage tank shown in Section 4 via validation against the model published by [26].
with proper insulation. Maeda et al. [28] reported proper tank filing The usefulness of the model is then illustrated in Section 5 through
pro- cedures using CFD simulations aimed at improving the safety. various terminal design problems and operational scenarios.
Likewise, Zakaria et al. [30] studied boil-off gas generation from a
marine storage tank via CFD and performed some limited parametric Trong nghiên cứu hiện tại, một bể chứa LNG kết hợp với dòng tuần
study on the re- lationship between boil-off gas generation and heat- hoàn được nghiên cứu chuyên sâu bằng cách sử dụng Aspen HYSYS
transfer parameter values. Rather strangely, the authors performed Dynamics, một phần mở rộng của bộ mô phỏng quy trình được sử
validation against the substantially simpler model developed by dụng rộng rãi Aspen HYSYS [32,33–36]. Nghiên cứu nắm bắt tất cả
Hasan et al. [31], as opposed to the usual practice of validating các khía cạnh tích cực của công việc nói trên, đồng thời mang lại mức
against a more detailed model. These studies highlight the inherent độ tự do cao hơn cho các sửa đổi trong tương lai, cũng như việc sử
limitations of CFD studies on LNG storage tanks: the heavy dụng các gói chất lỏng chính xác hơn [6]. Mô hình này không đẳng áp,
computational load prevents important details from being không cân bằng và dự kiến sẽ nắm bắt được hoạt động thiết yếu của
incorporated, and at the same time, the accuracy of the result does bể chứa LNG được kết hợp với đường tuần hoàn. Nghiên cứu được
not seem to be substantially different from simpler models in many cấu trúc như sau. Trong Phần 2, chúng tôi xem xét thiết kế bể chứa
significant ways. LNG và kết thúc tái hóa khí điển hình. Trong Phần 3, mô hình hóa
Hành vi động của bể chứa LNG đã được nghiên cứu chi tiết bằng LNG bể chứa trong Aspen HYSYS được xem xét. Mô-đun xe tăng
cách sử dụng động lực học chất lỏng tính toán (CFD) [22,28,29]. Bán trong Aspen HYSYS có một số thách thức, được giải quyết thông qua
và cộng sự. [22] kết luận rằng sự bay hơi bề mặt là cơ chế sôi chi một số sửa đổi đơn giản. Tính đúng đắn của những sửa đổi này được
phối và quá trình sôi mầm khó xảy ra trong bể chứa LNG có lớp cách thể hiện trong Phần 4 thông qua xác thực đối với mô hình được xuất
nhiệt thích hợp. Maeda và cộng sự. [28] đã báo cáo quy trình nộp hồ bản bởi [26]. Tính hữu ích của mô hình sau đó được minh họa trong
sơ thích hợp bằng cách sử dụng mô phỏng CFD nhằm cải thiện sự Phần 5 thông qua các vấn đề thiết kế thiết bị đầu cuối khác nhau và
an toàn. Tương tự như vậy, Zakaria et al. [30] đã nghiên cứu quá các kịch bản vận hành
trình tạo khí sôi từ bể chứa hàng hải thông qua CFD và thực hiện một
số nghiên cứu tham số hạn chế về mối quan hệ giữa quá trình tạo khí
sôi và các giá trị tham số truyền nhiệt. Khá kỳ lạ, các tác giả đã thực
hiện xác nhận đối với mô hình đơn giản hơn đáng kể được phát triển
bởi Hasan et al. [31], trái ngược với thông lệ kiểm chứng thông
thường đối với một mô hình chi tiết hơn. Các nghiên cứu này nêu bật
những hạn chế vốn có của các nghiên cứu CFD trên bể chứa LNG:
khối lượng tính toán nặng nề khiến các chi tiết quan trọng không
được tích hợp, đồng thời, độ chính xác của kết quả dường như
không khác biệt đáng kể so với các mô hình đơn giản hơn theo nhiều
cách đáng kể .
The preceding discussion indicates several deficiencies in the models
proposed in the literature. Studies on LNG storage tanks are typically
performed in the absence of the recirculation line, while stu- dies
focusing on the recirculation neglect the dynamic behavior of the
tank. Studies on LNG storage tanks often assume static behavior,
with the boil-off gas generation calculated from either vapor-liquid
equili- brium or some assumed boil-off rate parameter. The only study
in the literature that incorporates the recirculation line, the dynamic
behavior of the tank, and vapor-liquid disequilibrium [26] assume an
isobaric tank with simplistic fluid behavior. Moreover, the model
presented therein is difficult to extend beyond the original problem
statement. Các cuộc thảo luận trước chỉ ra một số thiếu sót trong các
mô hình được đề xuất trong tài liệu. Các nghiên cứu về bể chứa LNG
thường được thực hiện khi không có đường tuần hoàn, trong khi các
nghiên cứu tập trung vào tuần hoàn lại bỏ qua đặc tính động của bể.
Các nghiên cứu về bể chứa LNG thường giả định trạng thái tĩnh, với
lượng khí sinh ra được tính toán từ trạng thái cân bằng hơi-lỏng hoặc
một số thông số tốc độ sôi giả định. Nghiên cứu duy nhất trong tài
liệu kết hợp đường tuần hoàn, hành vi động của bể và sự mất cân
bằng hơi-lỏng [26] giả định một bể đẳng áp với hành vi chất lỏng đơn
giản. Hơn nữa, mô hình được trình bày trong đó khó có thể mở rộng
ra ngoài tuyên bố vấn đề ban đầu.
In the present study, an LNG storage tank coupled to a
recirculation flow is studied in-depth using Aspen HYSYS Dynamics,
an extension of the widely-used process simulator Aspen HYSYS
[32,33–36]. The study captures all the positive aspects of the
aforementioned work, while also offering a greater degree of freedom
for future modifications, as well as the use of more accurate fluid
packages [6]. The model is non-isobaric, non-equilibrium, and is
expected to capture the essential behavior of an LNG storage tank
coupled to a recirculation line. The study is structured as follows. In
34
 M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

2. Regasification terminal operation of tank operation in all three modes in order to understand, analyze,
and improve regasification terminal operation.
A typical liquefied natural gas (LNG) regasification terminal oper- A simplified diagram of an LNG regasification terminal is illustrated
ates in three different modes, namely, holding mode, in which no in Fig. 1. In the holding mode, no transfer of LNG occurs to/from the
transfer of LNG occurs to and from LNG tanks, unloading mode, in terminal tank [37]. The BOG is generated due to heat ingress via the
which tanks receive LNG from the tanker, and reloading mode, in recirculation line (in-line BOG) and via the tank surface (in-tank BOG).
which tanks send LNG to the tanker [37]. Regardless of the mode of The BOG is removed using a compressor to maintain the tank
operation, regasification of LNG at the desired rate and the cooling of pressure and cooled down for subsequent liquefaction [39]. BOG
transfer lines connecting tanks and jetties is continuously performed removal in excess of the compressor capacity is flared. At the same
[38]. The LNG used to cool the transfer lines is recycled to the tanks time, LNG is pumped from the tank using in-tank pumps (primary
[7], generating in- line boil-off gas (BOG) in addition to the in-tank pumping) and split into a recirculating and a send-out stream [40].
BOG and modifying the temperature of the vapor phase in the tank. Part of the send- out stream is mixed with the compressed BOG and
This, in turn, modifies the heat transfer from the vapor phase into the sent to the re- condenser. The resultant stream is re-merged with the
liquid phase, which affects the in-tank BOG generation. This brief remaining send- out stream and sent to the high pressure (HP) pump
description illustrates the need for a comprehensive dynamic model (secondary

Fig. 1. Schematic of a typical liquefied natural gas regasification terminal.

1. Vận hành kho tái hóa khí

Một kho cảng tái khí hóa khí tự nhiên hóa lỏng (LNG) điển hình hoạt động ở ba chế độ khác nhau, cụ thể là
chế độ giữ, trong đó không xảy ra việc chuyển LNG đến và từ các bồn chứa LNG, chế độ dỡ hàng, trong đó
các bồn chứa nhận LNG từ tàu chở dầu và nạp lại chế độ, trong đó xe tăng gửi LNG cho tàu chở dầu [37].
Bất kể chế độ vận hành nào, việc tái khí hóa LNG ở tốc độ mong muốn và làm mát các đường ống vận
chuyển nối các bồn chứa và cầu cảng được thực hiện liên tục [38]. LNG được sử dụng để làm mát các
đường truyền được tái chế vào các bồn chứa [7], tạo ra khí bốc hơi nội tuyến (BOG) ngoài BOG trong bồn
chứa và điều chỉnh nhiệt độ của pha hơi trong bồn chứa. Đến lượt mình, điều này làm thay đổi sự truyền
nhiệt từ pha hơi sang pha lỏng, ảnh hưởng đến quá trình tạo BOG trong bể. Mô tả ngắn gọn này minh họa
sự cần thiết của một mô hình vận hành bể năng động toàn diện ở cả ba chế độ để hiểu, phân tích và cải
thiện hoạt động của thiết bị đầu cuối tái hóa khí.
Sơ đồ đơn giản hóa của thiết bị đầu cuối tái khí hóa LNG được minh họa trong Hình 1. Ở chế độ giữ, không
xảy ra việc chuyển LNG đến/từ bể đầu cuối [37]. BOG được tạo ra do sự xâm nhập của nhiệt qua đường
tuần hoàn (BOG trong dòng) và qua bề mặt bể (BOG trong bể). BOG được loại bỏ bằng máy nén để duy trì
áp suất bể và làm mát cho quá trình hóa lỏng tiếp theo [39]. Loại bỏ BOG vượt quá công suất máy nén bị
bùng. Đồng thời, LNG được bơm ra khỏi bể chứa bằng máy bơm trong bể (bơm sơ cấp) và được chia thành
dòng tuần hoàn và dòng đưa ra [40]. Một phần của luồng gửi đi được trộn với BOG đã nén và được gửi đến
bộ ngưng tụ. Luồng kết quả được hợp nhất lại với luồng gửi đi còn lại và được gửi đến máy bơm áp suất
34
 M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353
cao (HP) (thứ cấp

34
 M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

pumping) [41]. In the secondary stage, the LNG is pumped through

vaporizers at a high pressure. The vaporizers convert LNG into natural Table 1
gas (NG) [42]. Two types of vaporizers are used in regasification Parameters used for the present study. The chosen composition is that of a
terminals. Open rack vaporizers use seawater to heat LNG, and sub- typical “light” LNG [22].
merged combustion vaporizers burn a portion of the vaporized NG to Parameter Value Units
heat a water tank, which acts as a heat bath for the LNG [43]. The NG
from the vaporizer is subsequently directed via a common collector to LNG STORAGE TANK
3
the regulation, measurement, and odorization systems, and finally in- Volume 200,000 m
Diameter 63.39 m
Height 63.39 m
jected to the city gas grid [37]. Top side HTC 0.02371
2
W/m -K
2
In unloading mode, a carrier arriving at the jetty unloads its cargo Shell side liquid HTC 0.02136 W/m -K
2
through the unloading arms. The unloading arms are connected to the Shell side vapour HTC 0.02121 W/m -K
2
Bottom side HTC 0.02216 W/m -K
rundown line and transport the LNG from the jetty to the storage tank 2
Vapor-liquid interface HTC 4 W/m -K
[44]. During the unloading mode, heat ingress into the LNG (via ship
Tank bottom temperature 5 °C
pumps and ambient heat leak through the pipes) causes a fraction of Feed vapor flash efficiency 0 %
the liquid to vaporize when it flashes inside the tanks. A portion of the Holdup vapor flash efficiency 0 %
flashed LNG is recycled back to the carrier to maintain carrier tank RECIRCULATION Pipe
pressure [26]. Length 1500 m
bơm) [41]. Ở giai đoạn thứ cấp, LNG được bơm qua các Diameter 40 inch
Pipe material Smooth –
thiết bị hóa hơi ở áp suất cao. Các thiết bị hóa hơi chuyển Pipe HTC 0.7092 kJ/h-m -°C
2

đổi LNG thành khí tự nhiên (NG) [42]. Hai loại bình hóa IN-TANK Pump
hơi được sử dụng trong các thiết bị đầu cuối khí hóa. Các Speed 2975 rpm
thiết bị hóa hơi giá mở sử dụng nước biển để đốt nóng Best operating point (head) 250 M
3
LNG, và các thiết bị hóa hơi đốt chìm đốt cháy một phần Best operating point (flow)
Best operating point (efficiency)
452
75
m /h
%
NG đã hóa hơi để làm nóng một bể chứa nước, hoạt động Sendout discharge pressure 11 bar
như một bể nhiệt cho LNG [43]. NG từ máy hóa hơi sau Compressor
đó được dẫn qua một bộ thu chung đến các hệ thống điều Adiabatic Efficiency 85 %
tiết, đo lường và tạo mùi, và cuối cùng được bơm vào lưới Speed 1337 rpm
3
Capacity 1045 m /h
khí đốt của thành phố [37].
VALVES
Ở chế độ dỡ hàng, người vận chuyển đến cầu cảng dỡ Recirculation valve 1 (Cv) 150 USGPM [(60 °F, 1 psi)]
hàng của mình thông qua các cánh tay dỡ hàng. Các Recirculation valve 2 (Cv) 150 USGPM [(60 °F, 1 psi)]
cánh tay dỡ hàng được kết nối với dây chuyền tải xuống Sendout valve (Cv) 2000 USGPM [(60 °F, 1 psi)]
Ambient temperature 30 °C
và vận chuyển LNG từ cầu tàu đến bể chứa [44]. Trong
chế độ không tải, nhiệt xâm nhập vào LNG (thông qua LNG Composition
Methane 0.9985 mole fraction
máy bơm tàu và rò rỉ nhiệt xung quanh qua đường ống) Nitrogen 0.0015 mole fraction
khiến một phần chất lỏng bốc hơi khi nó bốc hơi bên trong
các bể chứa. Một phần LNG đã bốc hơi được tái chế trở
lại chất mang để duy trì áp suất của bể chứa chất mang via Aspen HYSYS Dynamics, and subsequently used to study the
[26]. time- dependent behaviour of the LNG regasification terminal.
3. Aspen HYSYS tank module The proper adjustments effected on the tank module and the
effects of those adjustments are discussed in Sections 3.1–3.5. In
A schematic of the liquefied natural gas (LNG) regasification term- Section 4, the validity of these modifications is checked against an
inal developed in Aspen HYSYS is shown in Fig. 2. The terminal is as- established model
sumed to operate in the holding mode where some amount of LNG is 1. Mô-đun bồn chứa Aspen HYSYS in the literature.
withdrawn from the storage tank by in-tank pump one and in-tank
pump two. The LNG withdrawn from the tank (LNG-1) is subsequently Sơ đồ thiết bị đầu cuối khí tự nhiên hóa lỏng (LNG) được phát triển trong
split into a send-out stream (LNG-2) and a recirculation stream (Rec- Aspen HYSYS được thể hiện trong Hình 2. Thiết bị đầu cuối được giả
LNG-0). The recirculation stream passes through the loading/un- định là hoạt động ở chế độ giữ trong đó một lượng LNG được rút ra khỏi
loading/jetty lines and eventually returns to the tank. Heat ingress into bể chứa bằng cách -bơm bể một và bơm trong bể hai. LNG được rút ra
the tank and the recirculation lines, as well as the in-tank pump duty khỏi bể chứa (LNG-1) sau đó được chia thành luồng gửi đi (LNG-2) và
associated with the recirculation stream, ultimately results in the for- luồng tuần hoàn (Rec-LNG-0). Dòng tuần hoàn đi qua các đường nạp/dỡ
mation of boil-off gas (BOG-0), which is withdrawn from the tank using hàng/cầu cảng và cuối cùng quay trở lại bể chứa. Nhiệt xâm nhập vào bể
a boil-off gas compressor. và các đường tuần hoàn, cũng như nhiệm vụ bơm trong bể liên quan đến
As the name suggests, the in-tank pumps are supposed to be dòng tuần hoàn, cuối cùng dẫn đến sự hình thành khí sôi (BOG-0), khí
located inside the tank. However, since a majority of the inefficient part này được rút ra khỏi bể bằng cách sử dụng một máy nén khí sôi.
of the pump duty i.e. that which results in temperature increase as Như tên cho thấy, các máy bơm trong bể được cho là được đặt bên trong
opposed to pressure increase, exits with LNG-1 [45]. Thus, the in-tank bể. Tuy nhiên, do phần lớn phần không hiệu quả của nhiệm vụ bơm, tức
pumps may be considered as if they are located outside of the tank. All là phần dẫn đến tăng nhiệt độ trái ngược với tăng áp suất, thoát ra cùng
parameters relevant to the present study are listed in Table 1. The với LNG-1 [45]. Do đó, các máy bơm trong bể có thể được coi như thể
tank module in Aspen HYSYS V9 is used to simulate the storage tank. chúng được đặt bên ngoài bể. Tất cả các tham số liên quan đến nghiên
However, the tank module requires some theoretical adjustments prior cứu hiện tại được liệt kê trong Bảng 1. Mô-đun bể chứa trong Aspen
to use, due to inbuilt limitations and assumptions in the tank module HYSYS V9 được sử dụng để mô phỏng bể chứa. Tuy nhiên, mô-đun xe
options. After the steady-state model is developed, it is converted to a tăng yêu cầu một số điều chỉnh lý thuyết trước khi sử dụng, do các giới
dynamic model hạn và giả định sẵn có trong các tùy chọn mô-đun xe tăng. Sau khi mô
hình trạng thái ổn định được phát triển, nó được chuyển đổi sang mô hình

35
 M.S. KHAN et AL. AppliedThermAlEngineeri
động thông qua Aspen HYSYS Dynamics, và sau đó được sử dụng để ng149(2019)344–353
Các điều chỉnh thích hợp được thực hiện trên mô-đun xe tăng và tác động
nghiên cứu hành vi phụ thuộc vào thời gian của thiết bị đầu cuối tái khí của những điều chỉnh đó được thảo luận trong Phần 3.1–3.5. Trong Phần
hóa LNG. 4, tính hợp lệ của những sửa đổi này được kiểm tra dựa trên một mô hình
đã thiết lập

Fig. 2. A schematic of the setup under consideration.

35
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

3.1. HEAT TRANSFER AREA and without insulator heat capacity using the baseline parameter
values in Table 1, and the procedure outlined in Section 3.1 through
The Aspen HYSYS tank module does not automatically incorporate Section
heat-transfer from the top- and bottom surfaces of the tank. A heating 3.5. As observed from Fig. 4, there is no appreciable difference
duty is present, which can be used to account for heat-transfer from between these two cases; the total boil-off gas generation increases by
the bottom surface of the tank: merely 0.36% when the insulator heat capacity is removed. This is
3.1. KHU VỰC TRUYỀN NHIỆT reasonable agreement with the conclusion drawn by Effendy et al. [26]
via a simple energy balance. From the perspective of long-term
Mô-đun bồn chứa Aspen HYSYS không tự động kết hợp truyền behavior, the in- tegral effect of insulator heat capacity on the boil-off
nhiệt từ bề mặt trên và dưới của bồn chứa. Có một nhiệm vụ sưởi ấm, gas generation is expected to be nearly zero, since the temperature of
có thể được sử dụng để giải thích cho sự truyền nhiệt từ bề mặt đáy the insulator must
của bể: Quá trình thay thế này có phần liên quan nhiều hơn so với khái niệm
D2 đơn giản được nêu ở trên. Vấn đề nằm ở quan sát rằng Aspen
HYSYS cho phép người dùng chỉ định bốn thông số truyền nhiệt,
tương ứng với khả năng truyền nhiệt từ bên ngoài đến vỏ (UO,S), chất
cách điện (Mỹ), khả năng truyền nhiệt từ vỏ đến vỏ và không có lớp
cách nhiệt bằng cách sử dụng các giá trị tham số cơ sở trong Bảng 1
và quy trình được nêu trong Phần 3.1 đến Phần
3.5. Theo quan sát từ Hình 4, không có sự khác biệt đáng kể nào giữa
hai trường hợp này; tổng lượng khí sinh ra từ quá trình đun sôi chỉ
tăng 0,36% khi loại bỏ công suất nhiệt của chất cách nhiệt. Đây là thỏa
thuận hợp lý với kết luận được rút ra bởi Effendy et al. [26] thông qua
cân bằng năng lượng đơn giản. Từ góc độ ứng xử dài hạn, tác động
tích hợp của công suất nhiệt của chất cách điện đối với việc tạo ra khí
sôi được cho là gần như bằng không, vì nhiệt độ của chất cách điện
phải
QB = UB reset back to some value in between ambient and liquid temperatures
(TA TL)
T
(1) when the tank is re-filled with LNG.
4
Here QB is the rate of heat transfer from the bottom surface of the tank,
which is set as the heating duty, UB is the bottom-surface heat-transfer 3.3. FLASH efficiency
coefficient, DT is the diameter of the tank, TL is the temperature of the
liquid, and TA is the ambient temperature. The default values of the flash calculation efficiencies used by the
The top-surface heat transfer cannot be incorporated into the tank module are 100%. This corresponds to the case whereby 100% of
heating duty since the heating duty is applied to the liquid phase of the the incoming liquid- and vapor-streams, and 100% of the liquid- and
tank. Instead, the top-surface heat transfer is accounted for by vapor-holdups, are participating in a flash calculation at every in-
changing the shell-to-inner-vapor heat-transfer coefficient according to tegration time step. Obviously, this is equivalent to an equilibrium as-
the height of the vapor phase: sumption. However, the equilibrium between the vapor and the liquid
Ở đây QB là tốc độ truyền nhiệt từ bề mặt đáy của bể,
được đặt làm nhiệm vụ gia nhiệt, UB là hệ số truyền nhiệt trên bề
mặt đáy, DT là đường kính của bể, TL là nhiệt độ của chất lỏng và TA
là nhiệt độ môi trường.
Việc truyền nhiệt trên bề mặt không thể được kết hợp vào nhiệm
vụ gia nhiệt vì nhiệm vụ gia nhiệt được áp dụng cho pha lỏng của bể.
Thay vào đó, sự truyền nhiệt của bề mặt trên cùng được tính bằng
cách thay đổi hệ số truyền nhiệt từ lớp vỏ sang lớp hơi bên trong theo
chiều cao của pha hơi: D
T
U U 1+ phase is unlikely, since heat and mass transfer is limited by the largely-
S,V+T S,V stagnant vapor phase [27]. To reduce the problem to a vapor-liquid
4(HT HL ) (2) disequilibrium, the recycle efficiency (vapor) and the feed (vapor) flash
Here HT is the height of the tank and HL is the height of the liquid efficiency are to be set to 0%. This corresponds to the case whereby
phase. The definitions of US,V+T and US,V will be deferred to the sub- 100% of the incoming liquid stream and 100% of the liquid holdup
sequent paragraphs. For now, it suffices to state that if the top-surface participate in a flash calculation at every time step while excluding the
heat-transfer coefficient is equal to the shell-surface heat-transfer holdup and feed vapor. Instead, the holdup and feed vapor undergo a
coefficient, then
separate single flash calculation and come into equilibrium. The
Ở đây HT là chiều cao của bể và HL là chiều cao của pha lỏng. Các
recycle and feed nozzle efficiencies may be accessed from the
định nghĩa về US,V+T và US,V sẽ được chuyển sang các đoạn tiếp
storage tank
theo. Hiện tại, chỉ cần nói rằng nếu hệ số truyền nhiệt của bề mặt trên
hiệu quả được đặt thành 0%. Điều này tương ứng với trường hợp
cùng bằng hệ số truyền nhiệt của bề mặt vỏ, thì
theo đó 100% dòng chất lỏng đi vào và 100% lượng chất lỏng giữ lại
D2
tham gia vào phép tính chớp nhoáng ở mọi bước trong khi loại trừ
lượng hơi nước nạp và giữ lại. Thay vào đó, hơi giữ và hơi nạp trải
qua một phép tính chớp nhoáng riêng biệt và đi vào trạng thái cân
bằng. Hiệu quả của vòi tái chế và vòi cấp liệu có thể được truy cập từ
bể chứa
QS,V + QT = US,V DT (HT HL )(TA TV ) + T TV ) control panel, under the “Dynamics → Holdup” tab, through the “Ad-
US,V (TA vanced” button. Another control panel for the holdup will then open,
DT 4
=U 1+
D (H H )(T T ) and the efficiencies can be found under the “Efficiencies” tab. This is
S,V +T
US,V 4(HT

35
M.S. KHAN et AL. AppliedThermAlEngineeri
DT (HT
T T illustrated in Fig. ng149(2019)344–353
HL ) A 5.

HL ) TV )
(TA

3.4. INITIAL EQUILIBRATION

Here QS,V is the heat transfer from the shell surface to the vapor phase
of the tank, QT is the heat transfer from the top surface of the tank, and The Aspen HYSYS tank initialization routine enforces equilibrium
TV is the vapor-phase temperature. It is thus clear that the replacement between the vapor and the liquid phase. This is likely not
of US,V with US,V+T allows the top-surface heat transfer to be accounted representative of a real tank, since the bulk vapor temperature will, in
for via the shell-side heat transfer. general, depart from the liquid phase. To break this vapor-liquid
Ở đây QS,V là nhiệt lượng truyền từ bề mặt vỏ sang pha hơi của bể, equilibrium, we ar- tificially raise the vapor phase temperature via the
QT là nhiệt truyền từ bề mặt trên cùng của bể và TV là nhiệt độ pha following procedure, which is motivated by the separation of
hơi. Do đó, rõ ràng là việc thay thế US,V bằng US,V+T cho phép tính
timescales between enthalpy accumulation in the vapor phase and the
toán quá trình truyền nhiệt trên bề mặt thông qua quá trình truyền
change in liquid height. For a typical regasification terminal operation
nhiệt phía vỏ.
at holding mode, the former process is much faster than the latter. This
This substitution process is somewhat more involved than the
sets up an estimate for the evaporation rate:
simple concept outlined above. The issue lies with the observation that
Aspen HYSYS allows the user to specify four heat-transfer
parameters, which correspond to the outer-to-shell (UO,S), insulator 3.4. CÂN BẰNG BAN ĐẦU
(US), shell-to-
inner-vapor (US,V ), and shell-to-inner-liquid (US,L) heat-transfer Ở đây QS,V là nhiệt lượng truyền từ bề mặt vỏ sang pha hơi của bể, QT
coeffi- là nhiệt truyền từ bề mặt trên cùng của bể và TV là nhiệt độ pha hơi. Do
cient. These are schematically shown in Fig. 3. đó, rõ ràng là việc thay thế US,V bằng US,V+T cho phép tính toán quá
trình truyền nhiệt trên bề mặt thông qua quá trình truyền nhiệt phía vỏ.
Ở đây QS,V là nhiệt lượng truyền từ bề mặt vỏ sang pha hơi của bể, QT
là nhiệt truyền từ bề mặt trên cùng của bể và TV là nhiệt độ pha hơi. Do
đó, rõ ràng là việc thay thế US,V bằng US,V+T cho phép tính toán quá
trình truyền nhiệt trên bề mặt thông qua quá trình truyền nhiệt vỏ.
Quá trình thay thế này có phần liên quan nhiều hơn so với khái niệm đơn
giản được nêu ở trên. Vấn đề nằm ở quan sát rằng Aspen HYSYS cho
phép người dùng chỉ định bốn tham số truyền nhiệt, tương ứng với lớp vỏ
ngoài với vỏ (UO,S), chất cách điện (Mỹ), lớp vỏ với lớp vỏ.
hệ số truyền nhiệt từ hơi bên trong (US,V ), và hệ số truyền nhiệt từ vỏ
sang chất lỏng bên trong (US,L)
cổ đại. Chúng được thể hiện dưới dạng sơ đồ trong Hình 3.

Quy trình khởi tạo bể Aspen HYSYS đảm bảo trạng thái cân bằng giữa
pha hơi và pha lỏng. Điều này có thể không đại diện cho một bể thực, vì
nói chung, nhiệt độ hơi khối sẽ lệch khỏi pha lỏng. Để phá vỡ trạng thái
cân bằng hơi-lỏng này, chúng tôi tăng nhiệt độ pha hơi một cách giả tạo
thông qua quy trình sau, được thúc đẩy bởi sự phân tách các thang thời
gian giữa sự tích lũy entanpy trong pha hơi và sự thay đổi chiều cao của
chất lỏng. Đối với hoạt động của thiết bị đầu cuối tái khí hóa điển hình ở
chế độ giữ, quy trình trước nhanh hơn nhiều so với quy trình sau. Điều
này thiết lập một ước tính cho tốc độ bay hơi:

Thông số truyền nhiệt cách nhiệt được chia sẻ bởi cả pha hơi và pha
lỏng, và do đó không thể điều chỉnh theo cách ngụ ý của biểu thức. (2).
Để thực hiện phương trình. (2) trong Aspen HYSYS, các giá trị của US,V
và US,L
đầu tiên được đặt thành các giá trị tổng hợp của US với US,V và US,L
tương ứng.
Giá trị của US sau đó được đặt thành 0 bằng cách giảm độ dày của chất
cách điện xuống 0. Cuối cùng, giá trị của US,V+T được tính từ biểu thức.
(2) và được sử dụng thay cho US,V . Để minh họa quy trình, hãy xem xét
ba bộ
các giá trị tham số được hiển thị trong Bảng 2. Cột đầu tiên trong Bảng 2
thể hiện các giá trị tham số “đúng”. Cột tiếp theo minh họa việc thay thế
US,V và US,L bằng các giá trị tổng hợp của US với US,V
và US,L tương ứng. Giá trị của US sau đó được đặt bằng 0 và
trong cột cuối cùng, giá trị của US,V+T được sử dụng thay cho US,V .
3.2. CÔNG SUẤT CÁCH NHIỆT
Việc thực hiện phương trình. (2) loại trừ ảnh hưởng của công suất nhiệt
của chất cách điện, vì quá trình đặt độ dày lớp cách nhiệt về 0 cũng đặt
công suất nhiệt của chất cách điện về không. Độ chính xác của loại trừ
này được định lượng bằng cách thực hiện một nghiên cứu duy nhất về
hành vi của xe tăng với
35
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

The insulation heat-transfer parameter is shared by both the vapor 2 2

(3)
and liquid phase, and thus cannot be adjusted in the manner implied mE HE DT HL (TA TL) + UB DT TL) + U DT TL)
U
S,L (TA (TV
by Eq. (2). To execute Eq. (2) in Aspen HYSYS, the values of US,V and 4 4 (4)
US,L
are first set to the composite values of US with US,V and US,L
Here mE is the evaporation rate, HE is the enthalpy of vaporization,
respectively.
The value of US is then set to zero by reducing the insulator thickness and U is the interface heat-transfer coefficient. The three terms on the
to zero. Finally, the value of US,V+T is calculated from Eq. (2) and used right-hand side of Eq. (4) corresponds to heat transfer from the shell
in lieu of US,V . To illustrate the procedure, consider the three sets
of
parameter values shown in Table 2. The first column in Table 2 re-
presents the “true” parameter values. The next column illustrates the
replacement of US,V and US,L with the composite values of US with US,V
and US,L respectively. The value of US is then set to be equal to zero,
and
in the last column, the value of US,V+T is used in lieu of US,V .

3.2. INSULATOR HEAT CAPACITY

The execution of Eq. (2) precludes the effect of insulator heat ca-
Fig. 3. Heat transfer parameters in Aspen HYSYS tank module.
pacity, as the process of setting the insulation thickness to zero also
sets the insulator heat capacity to zero. The accuracy of this
preclusion is quantified by performing a single study on the behavior of
the tank with

35
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

Table 2
An illustration of the update procedure. The superscript (k) refers to the value after the k-th update procedure.
2
“True” Values (W/m -K) Step 1 Step 2 Step 3

(0) (1) (0) (2) (1) (3) (2)

O,S=
U O,S=
U O,S O,S
U= O,S O,S
U= O,S
(0) (1) (0) (2) (3) (2)
U = U =U U =0 U =U
0.04
S S S S

U (0) = 3 U (2) = U (1) (3) (2) 1+

DT
S,V 1 1 S,V S,V US,V = US,V 4(HT HL)
U (1) = +
S,V (0) (0)
(0) SU U
S,V
U = 300 (2)
U S,L (1)
= U S,L U (3) = U (2)
1 1
S,L
U (1) +
S,L S,L
=
S,L (0)
U (0)
U
S S,L

3.5. Minimum AND MAXIMUM liquid levels

The Aspen HYSYS tank module possesses an unstated maximum

(and minimum) liquid level, above (and below) which the simulation
result cannot be trusted. The minimum liquid level in the Aspen
HYSYS tank module is set by the diameter of the liquid outer piping,
which is assumed to be located at the bottom of the shell surface of
the tank. Below this level, the vapor phase starts to exit the tank. In
practice, the minimum liquid level, also called the tank heel, is dictated
by the need to avoid cavitation in the in-tank pump. Concomitantly,
simulation studies should estimate the heel level prior to the study and
stop the simulation when the liquid level drops below the heel level.
There is no theoretical maximum to the liquid level that the Aspen
HYSYS tank module can simulate. However, as the liquid level ap-
Fig. 4. Effect of insulator heat capacity on boil-off gas generation. proaches 100%, the appropriate time step required to bound the si-
mulation error approaches zero, and thus the simulation time
increases indefinitely. This relationship between the liquid level and the
surface, bottom surface, and the interface of the liquid phase. A similar
com- putational time arises due to the relaxation time for vapor-phase
estimate can be written out for the vapor phase temperature:
tem- perature and composition, which approaches zero as the vapor
bề mặt, bề mặt đáy và giao diện của pha lỏng. Một ước tính tương tự
phase
có thể được viết ra cho nhiệt độ pha hơi:
mE Cp,V (TV 2 2 vanishes. We suggest avoiding any simulation where the liquid height
(3)
TL) + U DT TL) DT TV )
(TV U
S,V (TA (5) exceeds 98% of the tank height. Fig. 6 illustrates the effect of liquid
4 4 level on BOG generation and vapor-phase temperature.
Here Cp,V is the vapor-phase heat capacity. Eq. (4) and Eq. (5) are
solved simultaneously to obtain crude estimates of the pseudo steady 3.5. Mức chất lỏng TỐI THIỂU VÀ TỐI ĐA
evaporation rate and vapor-phase temperature.
With these pseudo steady estimates, the vapor-phase temperature Mô-đun bồn chứa Aspen HYSYS sở hữu mức chất lỏng tối đa (và tối
can be raised to a value which more closely approximates a typical thiểu) không xác định, trên (và dưới) mà kết quả mô phỏng không
thể tin cậy được. Mức chất lỏng tối thiểu trong mô-đun bồn chứa
storage tank. Heat transfer to the liquid phase is stopped by setting
Aspen HYSYS được đặt theo đường kính của đường ống dẫn chất
the shell-to-inner-liquid and the interface heat-transfer coefficient to lỏng bên ngoài, đường ống này được cho là nằm ở đáy bề mặt vỏ
zero. The simulation is run until the vapor temperature reaches the của bồn chứa. Dưới mức này, pha hơi bắt đầu thoát ra khỏi bể.
pseudo- steady value obtained from Eq. (4) and Eq. (5). This is the Trong thực tế, mức chất lỏng tối thiểu, còn được gọi là gót bể,
desired initial state, which will be used in Section 5 of the present được quyết định bởi nhu cầu tránh tạo bọt trong máy bơm trong
study. bể. Đồng thời, các nghiên cứu mô phỏng nên ước tính mức gót
Ở đây Cp,V là nhiệt dung pha hơi. phương trình (4) và phương trình. chân trước khi nghiên cứu và dừng mô phỏng khi mức chất lỏng
(5) là giảm xuống dưới mức gót chân.
được giải quyết đồng thời để có được các ước tính thô về tốc độ
Không có mức tối đa theo lý thuyết đối với mức chất lỏng mà mô-
bay hơi ổn định giả và nhiệt độ pha hơi.
đun bồn chứa Aspen HYSYS có thể mô phỏng. Tuy nhiên, khi mức
Với các ước tính ổn định giả này, nhiệt độ pha hơi có thể được chất lỏng đạt tới 100%, bước thời gian thích hợp cần thiết để giới
nâng lên một giá trị gần đúng hơn với một bể chứa điển hình. Quá hạn lỗi mô phỏng tiến tới 0, và do đó thời gian mô phỏng tăng vô
trình truyền nhiệt sang pha lỏng được dừng lại bằng cách đặt hệ số hạn. Mối quan hệ này giữa mức chất lỏng và thời gian tính toán
truyền nhiệt từ vỏ sang chất lỏng bên trong và hệ số truyền nhiệt của phát sinh do thời gian thư giãn đối với nhiệt độ và thành phần của
giao diện bằng không. Mô phỏng được chạy cho đến khi nhiệt độ hơi pha hơi, tiến tới 0 khi pha hơi biến mất. Chúng tôi khuyên bạn nên
đạt đến giá trị giả ổn định thu được từ biểu thức. (4) và phương trình. tránh bất kỳ mô phỏng nào khi chiều cao chất lỏng vượt quá 98%
(5). Đây là trạng thái ban đầu mong muốn, sẽ được sử dụng trong chiều cao bể. Hình 6 minh họa ảnh hưởng của mức chất lỏng đối
với quá trình tạo BOG và nhiệt độ pha hơi.
Phần 5 của nghiên cứu này.

35
 M.S. KHAN et AL. AppliedThermAlEngineeri
cular, there is no indication inng149(2019)344–353
the Aspen HYSYS manual on the
absence of heat transfer from the top surface of the tank. The
appropriateness of the modifications proposed in Section 3 is shown
via comparison with the model published by Effendy et al. [26], which
adheres closely to the ideal formulation described in the Aspen
HYSYS manual. The para- meters used by Effendy et al. are used to
generate the vapor-phase temperature and the rate of boil-off gas
(BOG) generation, which are
4. Xác nhận mô hình

Các sửa đổi được đề xuất trong Phần 3 là giả thuyết. Đặc biệt,
không có chỉ dẫn nào trong hướng dẫn sử dụng Aspen HYSYS về việc
không có sự truyền nhiệt từ bề mặt trên cùng của bể chứa. Sự phù
hợp của các sửa đổi được đề xuất trong Phần 3 được thể hiện thông
qua so sánh với mô hình do Effendy và cộng sự công bố. [26], tuân
4. Model validation thủ chặt chẽ công thức lý tưởng được mô tả trong sổ tay Aspen
HYSYS. Các tham số được sử dụng bởi Effendy et al. được sử dụng
The modifications proposed in Section 3 are hypothetical. In parti- để tạo ra nhiệt độ pha hơi và tốc độ tạo ra khí sôi (BOG), đó là

Fig. 5. Modifying the recycle and feed nozzle efficiencies.

35
 M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

Fig. 6. Effect of the initial liquid level (LL) on (a) BOG generation and (b) vapor-phase temperature. Observe that the predicted profiles become noisier as LL
approaches 100%.

shown in Fig. 7. The vapor-phase temperature is found to agree to khí tự nhiên hóa lỏng (LNG) trong một số điều kiện nhất định. Hiện
within 0.1 K of the temperature profile reported by Effendy et al. in tượng rollover được giải thích tốt nhất bằng sự tồn tại của nhiều lớp với
Fig. 2 of their publication, and the rate BOG generation is in general nhiệt độ và mật độ riêng biệt [46].
agreement with Fig. 7 of their publication, with the present study re- Tuy nhiên, nghiên cứu hiện tại lập luận rằng mô hình được mô tả
porting a more pronounced dip in the rate at the start of the simulation. trong Phần 2 và được hiệu chỉnh trong Phần 3 có khả năng nắm bắt
In particular, the pseudo-steady rate of BOG generation is found to được nhiều chi tiết quan trọng trong hoạt động của kho cảng tái khí hóa
agree to within 5%. LNG. Sự thống nhất chặt chẽ giữa mô hình được sử dụng trong bản
The results presented in Fig. 7 also includes the cases whereby thảo và mô hình được sử dụng bởi Effendy et al. [26] gợi ý rằng mô
(a) the correction factor in Eq. (2) is not included, and (b) the hình ít nhất có thể nắm bắt được sự cân bằng năng lượng và khối
recycle and feed nozzle efficiencies are set to 100%. It is clear that lượng thiết yếu của thiết lập được đề xuất. Trong Phần 5.1 đến Phần
the proposed modifications are necessary to match the result 5.5, tiện ích của mô hình sẽ được minh họa thông qua nhiều vấn đề
reported by the afore- mentioned publication. vận hành và thiết kế liên quan đến công nghiệp.
được hiển thị trong Hình 7. Nhiệt độ pha hơi được tìm thấy là phù hợp
trong phạm vi 0,1 K của cấu hình nhiệt độ được báo cáo bởi Effendy
et al. trong Hình 2 trong ấn phẩm của họ và tốc độ tạo BOG nói chung
phù hợp với Hình 7 trong ấn phẩm của họ, với nghiên cứu hiện tại báo
cáo tốc độ giảm rõ rệt hơn khi bắt đầu mô phỏng. Đặc biệt, tốc độ tạo
BOG giả ổn định được cho là đồng ý trong phạm vi 5%.
Các kết quả được trình bày trong Hình 7 cũng bao gồm các trường
hợp theo đó (a) hệ số hiệu chỉnh trong biểu thức. (2) không được bao
gồm và (b) hiệu quả của vòi phun tái chế và cấp liệu được đặt thành
100%. Rõ ràng là những sửa đổi được đề xuất là cần thiết để phù hợp
với kết quả được báo cáo bởi ấn phẩm nói trên.

5. Tank design and terminal operation

Several limitations to the tank module cannot be overcome using

only the tools offered by Aspen HYSYS. For example, a well-mixed as-
sumption is built into the mass and energy balances of the vapor and
liquid phase of the tank. This does not correspond to reality, as evi-
denced by, among other things, the rollover of liquefied natural gas
(LNG) under certain conditions. The rollover phenomenon is best ex-
plained by the existence of multiple layers with distinct temperatures
and densities [46].
Nevertheless, the present study argues that the model described in
Section 2 and corrected in Section 3 is capable of capturing many im-
portant details in an LNG regasification terminal operation. The close
agreement between the model used in the manuscript and that used
by Effendy et al. [26] suggests that the model is at least able to
capture the essential mass and energy balances of the proposed
setup. In Section 5.1 through Section 5.5, the utility of the model will
be illustrated through a variety of industrially relevant design and
operational problems.
5. Thiết kế bồn chứa và vận hành bến bãi
Không thể khắc phục một số hạn chế đối với mô-đun xe tăng nếu
chỉ sử dụng các công cụ do Aspen HYSYS cung cấp. Ví dụ, một giả
định trộn đều được xây dựng trong cân bằng khối lượng và năng
lượng của pha hơi và pha lỏng của bể. Điều này không tương ứng với
thực tế, bằng chứng là, trong số những điều khác, việc luân chuyển

35
M.S. KHAN et AL. AppliedThermAlEngineeri
5.1. Aspect RATIO OPTIMIZATION reported) suggests that there is ng149(2019)344–353
a generally increasing trend in the
optimal aspect ratio with in- creasing tank volume.
Optimization of tank aspect ratio (AR), defined as the ratio of tank The same aspect ratio optimization is performed using the
height to tank diameter, that results in minimum boil-off gas (BOG) algorithm of Effendy et al. [26] with the set of parameter values
generation is performed using the corrected tank model. The search reported in Table 1. The optimal aspect ratio is found to be ca. 0.9,
is performed for several representative tanks’ volumes ranging from with an in- creasing trend with respect to tank volume, which is in
100,000 to 350,000 m3, in line with the current state-of-the-art LNG agreement with the result we found using the corrected tank model.
storage tank design. For each tank volume, the AR is varied at This optimal aspect ratio may be contrasted to the usual value of
intervals of 0.5 from 0.5 to 2.5. Three parameters, namely, the 0.5 reported in the lit- erature [24], which is obtained by assuming
terminal send-out capacity, the recirculation rate, and the recirculation vapor-liquid equilibrium. The difference may be attributed to the
line length are also systematically varied about the base values listed presence of a BOG stream ex- iting the tank (stream BOG-0 in Fig.
in Table 1. The average BOG generation rate is calculated from the 1), which diverts part of the heat ingress into the vapor phase of the
simulation output and used to determine the optimal aspect ratio to tank as sensible heat, instead of evaporating the liquid phase. This
within ± 0.25 for each case. The result is summarized in Fig. 8. The reduces the importance of the tank surface area in contact with the
optimum aspect ratio is found to be within the range of 0.75 to 1.25 in vapor, which allows for a taller tank.
5.1. TỐI ƯU HÓA TỶ LỆ KHÓA HÌNH Mọi trường hợp. Một phân tích riêng biệt sử dụng phép nội suy bậc hai
(không được báo cáo) cho thấy rằng nhìn chung tỷ lệ khung hình tối
Việc tối ưu hóa tỷ lệ khung hình của bể (AR), được định nghĩa là ưu có xu hướng tăng lên khi thể tích bể tăng lên.
tỷ lệ giữa chiều cao bể và đường kính bể, dẫn đến việc tạo ra khí bay Việc tối ưu hóa tỷ lệ khung hình tương tự được thực hiện bằng
hơi (BOG) tối thiểu được thực hiện bằng cách sử dụng mô hình bể thuật toán của Effendy et al. [26] với tập hợp các giá trị tham số được
đã hiệu chỉnh. Việc tìm kiếm được thực hiện đối với một số bể chứa báo cáo trong Bảng 1. Tỷ lệ khung hình tối ưu được tìm thấy là
đại diện có thể tích từ 100.000 đến 350.000 m3, phù hợp với thiết kế khoảng. 0,9, với xu hướng tăng đối với thể tích bể, phù hợp với kết
bể chứa LNG tiên tiến nhất hiện nay. Đối với mỗi thể tích bể, AR thay quả mà chúng tôi tìm thấy khi sử dụng mô hình bể đã hiệu chỉnh. Tỉ lệ
đổi trong khoảng 0,5 từ 0,5 đến 2,5. Ba tham số, cụ thể là công suất khung hình tối ưu này có thể trái ngược với giá trị thông thường là 0,5
phát của thiết bị đầu cuối, tốc độ tuần hoàn và chiều dài đường tuần được báo cáo trong tài liệu [24], giá trị này thu được bằng cách giả
hoàn cũng được thay đổi một cách có hệ thống về các giá trị cơ sở định trạng thái cân bằng hơi-lỏng. Sự khác biệt có thể là do sự hiện
được liệt kê trong Bảng 1. Tốc độ tạo BOG trung bình được tính toán diện của luồng BOG thoát ra khỏi bể (dòng BOG-0 trong Hình 1),
từ đầu ra mô phỏng và được sử dụng để xác định tỷ lệ khung hình tối luồng này chuyển một phần nhiệt đi vào pha hơi của bể dưới dạng
ưu trong khoảng ± 0,25 cho từng trường hợp. Kết quả được tóm tắt nhiệt cảm nhận được, thay vì làm bay hơi pha lỏng. Điều này làm
trong Hình 8. Tỷ lệ khung hình tối ưu nằm trong khoảng từ 0,75 đến giảm tầm quan trọng của diện tích bề mặt bể tiếp xúc với hơi, cho
1,25 in phép bể cao hơn.
Fig. 9 illustrates the effect of using an aspect ratio of 0.5, reported
all cases. A separate analysis using quadratic interpolation (not in

Fig. 7. Variation in (a) vapor-phase temperature and (b) BOG generation rate with and without the proposed modifications.

35
M.S. KHAN et AL. AppliedThermAlEngineeri
5.2. KẾ HOẠCH tắt máy ng149(2019)344–353

Mô hình này có thể được sử dụng để lên lịch tắt máy theo kế hoạch
của thiết bị đầu cuối khí hóa lại. Về mặt hình thức, mục tiêu có thể
được nêu là tối đa hóa khoảng thời gian “an toàn” của việc tắt máy theo
kế hoạch, được định nghĩa là khoảng thời gian mà áp suất bể chứa duy
trì dưới áp suất thiết kế, được coi là giá trị điển hình của 303 mbarg
[ 47]. Đây là một sự tối ưu hóa có vấn đề vì tất cả các tham số trong
Bảng 1 đều được cố định bởi thiết kế của thiết bị đầu cuối khí hóa lại.
Biến duy nhất còn lại là trường hợp bắt đầu tắt máy. Mức độ tự do này
tương đương với việc chọn chiều cao chất lỏng mà tại đó quá trình tắt
được bắt đầu do chiều cao chất lỏng giảm đơn điệu theo thời gian ở
chế độ giữ. Bể được khởi tạo với mức chất lỏng

equal to 95% of the tank height and is allowed to drain to 75%, 50%,
Fig. 8. Optimal tank aspect ratio at various terminal send-out capacities, LNG and 25% of the tank height following the normal holding mode op-
recirculation rates, and recirculation line lengths. The baseline terminal send- eration. The tank compressor and the send-out stream are eliminated
out capacity, LNG recirculation rate, and recirculation line length are 1 when the liquid level in the tank reaches the specified values, thus si-
MMTPA, 100 m3/h, and 700 m respectively. mulating the condition of a planned shutdown.
The resultant pressurization rate is shown in Fig. 10. The pressur-
ization rate of the tank with the highest liquid level is also the slowest,
bằng 95% chiều cao bể và được phép xả nước đến 75%, 50% và 25%
chiều cao bể sau khi vận hành ở chế độ giữ bình thường. Máy nén bể
chứa và dòng gửi ra được loại bỏ khi mức chất lỏng trong bể đạt đến
các giá trị được chỉ định, do đó mô phỏng điều kiện tắt máy theo kế
hoạch.
Tốc độ điều áp kết quả được thể hiện trong Hình 10. Tốc độ điều áp
của bể có mức chất lỏng cao nhất cũng là chậm nhất,

Fig. 9. Effect of suboptimal tank aspect ratio on average compressor duty as

a function of recirculation rate for various tank volumes. Tank Volume 1:
100,000 m3, Tank Volume 2: 200,000 m3, Tank Volume 3: 300,000 m3.

terms of compressor duty, which is nearly proportional to BOG gen-

eration, at tank volume 1 (100,000 m 3), tank volume 2 (200,000 m3),
and tank volume 3 (300,000 m 3). The compressor duty increases by
no more than 2.6% over a wide range of parameter values and tanks
vo- lumes, suggesting that BOG generation is insensitive to changes in
the aspect ratio around the optimal value of 1. This flexibility implies
that the thermal design of LNG storage tanks can, if necessary, be put
aside in favor of other considerations e.g. structural soundness.
về nhiệm vụ của máy nén, gần như tỷ lệ thuận với việc tạo ra BOG, ở
thể tích bể 1 (100.000 m3), thể tích bể 2 (200.000 m3) và thể tích bể 3
(300.000 m3). Nhiệm vụ của máy nén tăng không quá 2,6% trên một
loạt các giá trị tham số và thể tích bể chứa, cho thấy rằng việc tạo ra
BOG không nhạy cảm với những thay đổi về tỷ lệ khung hình xung
quanh giá trị tối ưu là 1. Tính linh hoạt này ngụ ý rằng thiết kế nhiệt
của Nếu cần thiết, các bể chứa LNG có thể được đặt sang một bên để
ưu tiên cho các cân nhắc khác, ví dụ: độ vững chắc về cấu trúc.

5.2. PLANNED shutdown

The model may be used to schedule the planned shutdown of a

regasification terminal. Formally, the objective may be stated as max-
imizing the “safe” duration of a planned shutdown, defined as the
duration for which the tank pressure remains below the design pres-
sure, which is taken to be the typical value of 303 mbarg [47]. This is a
problematic optimization since all the parameters in Table 1 are fixed
by the design of the regasification terminal. The only remaining vari-
able is the instance at which the shutdown is initiated. This degree of
freedom is equivalent to choosing the liquid height at which the shut-
down is initiated since the liquid height decreases monotonically with
time in the holding mode. The tank is initialized with a liquid level

35
M.S. KHAN et AL. AppliedThermAlEngineeri
and the pressurization rate of the tank with 25% liquid level is roughly ng149(2019)344–353
thấy trong Hình 10, ngoại trừ thời điểm đạt được áp suất thiết kế của
triple of that of the tank with 75% liquid level. This is a consequence bể hơi bị trễ một chút. Sự chậm trễ này có thể là do ảnh hưởng của
of the heat capacity of the liquid phase. As the tank pressure việc giảm mức chất lỏng, trong đó dòng chảy ra “tháo” bể, tạo ra
increases to the design pressure, the bubble point of the liquid-phase không gian hơi bổ sung, làm giảm áp suất của bể. Điều này có một ý
also increases. In the present simulation, the liquid-phase nghĩa khá bất thường – nếu việc tắt máy theo kế hoạch là cần thiết để
temperature is always equal to the bubble point temperature, thay thế máy nén khí đun sôi, thì trên thực tế, vận hành như bình
although this need not necessarily be the case depending on the thường sẽ an toàn hơn và chỉ cần đóng đường dây nối bình chứa với
chosen heat transfer coefficients. As a result, part of the heat ingress máy nén.
into the liquid-phase is diverted towards the sensible heating of the
liquid phase, with the remainder being used to evaporate the liquid.
5.4. Minimum RECIRCULATION RATE
Since the tank with the highest liquid level also possesses the highest
liquid-phase heat capacity, this tank also experi- ences the least
During a planned shutdown, the send-out LNG flow is
evaporation, and therefore the slowest rate of pressure increase.
eliminated. This prevents the use of several common strategies for
Following this line of reasoning, it seems apparent that a
reducing boil-off gas generation, such as circumventing the return
higher liquid level would maximize the safe duration of a planned
flow to the send-out line. The simplest strategy to reduce boil-off
shutdown. However, this argument is expected to fail in the limit of
gas generation appears to be, very simply, to reduce the
a vanishing vapor phase. In this limit, the heat capacity of the
recirculation flow. This minimizes the in- tank pump duty and heat
vapor phase vanishes, making it more susceptible to fluctuations
transfer with the ambient air, due to increased recirculation line
in temperature, and thus to fluctuations in pressure. This increase
temperature. On the other hand, the increased
in variability cannot be studied by the present model since this
5.4. TỶ SUẤT TUẦN HOÀN tối thiểu
model assumes a constant ambient temperature and zero
insulator heat capacity.
và tỷ lệ điều áp của bể có mức chất lỏng 25% gần gấp ba lần so với Trong quá trình ngừng hoạt động theo kế hoạch, luồng LNG gửi đi
bể có mức chất lỏng 75%. Đây là hệ quả của nhiệt dung của pha sẽ bị loại bỏ. Điều này ngăn cản việc sử dụng một số chiến lược phổ
lỏng. Khi áp suất bể tăng đến áp suất thiết kế, điểm bong bóng của biến để giảm phát sinh khí sôi, chẳng hạn như tránh dòng chảy trở lại
pha lỏng cũng tăng. Trong mô phỏng hiện tại, nhiệt độ pha lỏng luôn đường gửi đi. Chiến lược đơn giản nhất để giảm phát sinh khí sôi
bằng nhiệt độ điểm bong bóng, mặc dù điều này không nhất thiết dường như rất đơn giản là giảm lưu lượng tuần hoàn. Điều này giảm
phải như vậy tùy thuộc vào các hệ số truyền nhiệt đã chọn. Kết quả thiểu nhiệm vụ bơm bể chứa và truyền nhiệt với không khí xung
là, một phần nhiệt đi vào pha lỏng được chuyển hướng sang sự đốt quanh, do nhiệt độ đường tuần hoàn tăng lên. Mặt khác, sự gia tăng
nóng hợp lý của pha lỏng, phần còn lại được sử dụng để làm bay hơi
chất lỏng. Vì bể có mức chất lỏng cao nhất cũng sở hữu công suất
nhiệt pha lỏng cao nhất nên bể này cũng trải qua sự bay hơi ít nhất
và do đó tốc độ tăng áp suất chậm nhất.
Theo dòng lý luận này, có vẻ như rõ ràng là mức chất lỏng cao
hơn sẽ tối đa hóa khoảng thời gian an toàn cho việc tắt máy theo kế
hoạch. Tuy nhiên, lập luận này dự kiến sẽ thất bại trong giới hạn của
pha hơi biến mất. Trong giới hạn này, nhiệt dung của pha hơi biến
mất, làm cho nó dễ bị dao động về nhiệt độ và do đó dễ bị dao động
về áp suất. Mô hình hiện tại không thể nghiên cứu sự gia tăng độ
biến thiên này vì mô hình này giả định nhiệt độ môi trường không đổi
và công suất nhiệt của chất cách điện bằng không.

5.3. UNPLANNED compressor shutdown

Fig. 10. The pressurization rate of a single tank under planned shutdown for
The conditions in an unplanned compressor shutdown are nearly the set of parameter values listed in Table 1.
identical to those of the planned shutdown, except the send-out
stream is retained. In this case, we would like to study the effect of
liquid height as a key variable in determining the duration within which
the compressor shutdown must be resolved. The pressurization rate
is shown in Fig. 11 as a function of liquid height. This is very nearly
identical to the pressurization rate observed in Fig. 10, except the
point at which the tank design pressure is reached is slightly delayed.
This delay may be attributed to the effect of decreasing liquid level,
wherein the send-out stream “drains” the tank, which creates
additional vapor space, which lowers the pressure of the tank. This
has a rather unusual implication – if a planned shutdown is needed
for replacement of the boil-off gas compressor, it is, in fact, safer to
operate as usual, and simply close the line connecting the tank to
the compressor.
5.3. Tắt máy nén KHÔNG CÓ KẾ HOẠCH

Các điều kiện trong trường hợp tắt máy nén ngoài kế hoạch gần
giống với điều kiện tắt theo kế hoạch, ngoại trừ luồng gửi đi được giữ
lại. Trong trường hợp này, chúng tôi muốn nghiên cứu ảnh hưởng
của chiều cao chất lỏng như một biến số chính trong việc xác định
khoảng thời gian mà việc tắt máy nén phải được giải quyết. Tốc độ
điều áp được thể hiện trong Hình 11 dưới dạng hàm của chiều cao
chất lỏng. Tốc độ này gần giống với tốc độ điều áp được quan sát

36
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

Fig. 11. The pressurization rate of a single LNG tank under unplanned shut-
Fig. 13. Minimum recirculation rate as a function of the product of recircula-
down for the set of parameter values listed in Table 1.
tion line heat-transfer coefficient (U in kJ/h-m2-°C), recirculation line diameter
(D in m), and recirculation line length (L in m).

the recirculation rate. The vapor expansion effect is primarily a function

of send-out rate. Once the recirculation rate is fixed to that obtained
from Fig. 13, a correlation for the average compressor duty can be
constructed:

W = 52.1 + 0.149 × (VT

200) + 0.0694 × 288)
(UDL
0.0975 × (QS 200) (6)
characterized by the product of the pipe heat transfer coefficient (U
), the pipe diameter (D), and the pipe length (L). Consequently, it
should
5.5. TỔNG QUÁT KHÁC

Fig. 12. Effect of recirculation on boil-off gas generation. Bây giờ chúng ta xem xét sự khái quát hóa lập luận được trình
bày trong Phần 5.4. Như được chỉ ra bởi sự cân bằng năng lượng
dòng tuần hoàn của Effendy et al. [26] nhiệt xâm nhập vào đường
recirculation line temperature imposes the risk of having a two-
tuần hoàn có thể được đặc trưng bởi tích của hệ số truyền nhiệt của
phase flow in the recirculation line, which may damage the line and
ống (U ), đường kính ống (D) và chiều dài ống (L). Do đó, nó nên
its as- sociated valves.
Fig. 12 shows the effect of decreasing the recirculation flow. This
figure is generated by taking the integral average of the boil-off gas be possible to choose a minimum recirculation rate if this product
released after the tank design pressure is reached. It is evident that, in is known. This result is summarized in Fig. 13.
line with the argument proposed above, the boil-off gas generation Another generalization which can be obtained concerns the average
does indeed decrease with the flow rate in the recirculation line. At compressor duty. As discussed in Section 5.2 and Section 5.3, the
average compressor duty is effectively proportional to the BOG gen-
flow rates below 13 m3/h, two-phase flows are detected within the re-
eration, which can be decomposed into effects associated with in-tank
circulation line. By operating close to the two-phase flow regime, the
BOG, in-line BOG, and vapor expansion due to decreasing liquid
BOG generation can be reduced by ca. 15% compared to the baseline
height. The in-tank BOG is primarily a function of tank volume,
value. However, this occurs at the expense of process capability, since
assuming that the aspect ratio is of the order of 1 and that the heat
the system becomes completely susceptible to fluctuations in the en-
transfer coefficients of the tank are similar to the typical values reported
vironment.
in Table 1. The in- line BOG is almost entirely characterized by the
nhiệt độ của đường tuần hoàn có nguy cơ tạo ra dòng chảy hai pha
parameters UDL and
trong đường tuần hoàn, điều này có thể làm hỏng đường và các van
có thể chọn tốc độ tuần hoàn tối thiểu nếu biết sản phẩm này. Kết quả
liên kết của nó.
này được tóm tắt trong Hình 13.
Hình 12 cho thấy tác động của việc giảm lưu lượng tuần hoàn.
Một khái quát hóa khác có thể thu được liên quan đến chế độ làm
Con số này được tạo ra bằng cách lấy giá trị trung bình tích phân của
việc trung bình của máy nén. Như đã thảo luận trong Phần 5.2 và Phần
khí sôi thoát ra sau khi đạt đến áp suất thiết kế của bể. Rõ ràng là, phù
5.3, công suất trung bình của máy nén tỷ lệ thuận với việc tạo ra BOG,
hợp với lập luận được đề xuất ở trên, việc tạo ra khí sôi thực sự giảm
có thể được phân tách thành các hiệu ứng liên quan đến BOG trong
theo tốc độ dòng chảy trong đường tuần hoàn. Ở tốc độ dòng chảy
bình, BOG trong dòng và sự giãn nở hơi do chiều cao chất lỏng giảm.
dưới 13 m3/h, dòng chảy hai pha được phát hiện trong đường tuần
BOG trong bể chủ yếu là một chức năng của thể tích bể, giả sử rằng tỷ
hoàn. Bằng cách vận hành gần với chế độ dòng chảy hai pha, việc tạo
lệ khung hình ở mức 1 và hệ số truyền nhiệt của bể tương tự như các
ra BOG có thể giảm khoảng. 15% so với giá trị cơ sở. Tuy nhiên, điều
giá trị tiêu biểu được báo cáo trong Bảng 1. BOG nội tuyến là gần như
này xảy ra với cái giá phải trả là khả năng xử lý, vì hệ thống trở nên
hoàn toàn được đặc trưng bởi các tham số UDL và
hoàn toàn dễ bị ảnh hưởng bởi những biến động trong môi trường.

5.5. Other GENERALIZATIONS

We now consider a generalization of the argument presented in

Section 5.4. As indicated by the recirculation line energy balance of
Effendy et al. [26] the heat ingress into the recirculation line can be
35
 M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

known, the corresponding average compressor duty can be easily

obtained using a simple linear correlation.
6. Kết luận
Các mô-đun xe tăng của các mô phỏng quy trình thương mại
thường mang các giả định không được công bố về tính chất vật lý
của quy trình. Nghiên cứu hiện tại đã phát hiện ra một số giả định
không được nêu rõ đó và đã đề xuất một số phương pháp đơn
giản để khắc phục những giả định đó. Mô hình bể hiệu chỉnh sau
đó được mở rộng để bao gồm các khía cạnh khác của một thiết bị
đầu cuối tái khí hóa điển hình và được sử dụng để nghiên cứu thiết
kế và vận hành các thiết bị đầu cuối tái hóa khí trong một số bối
cảnh đơn giản nhưng có liên quan. Các kinh nghiệm mới lạ sau
đây đã được phát hiện:
Một tỷ lệ khung hình xe tăng của ca. 1 giảm thiểu việc tạo ra BOG
và nhiệm vụ máy nén trung bình.
Here W is the compressor duty in kW at the minimum recirculation Tốc độ điều áp của một bể chứa khi ngừng hoạt động theo kế
rate, VT is the volume of the tank in thousands of m 3, UDL is the hoạch/không theo kế hoạch sẽ giảm khi mức chất lỏng tăng lên.
characteristic parameter of the recirculation line in kJ/h-m 2-°C, and QS Tốc độ tuần hoàn phải được đặt ở giá trị tối thiểu của nó, giá trị này
is the send-out rate in m3/h. This correlation is constructed using a 3 k được quyết định bởi sự ràng buộc của dòng một pha trong đường
factorial design, setting all other parameter values to those of Table 1, tuần hoàn. Tốc độ tuần hoàn tối thiểu là duy nhất cho một UDL
and can capture more than 90% of the total variability of the nhất định trong đường tuần hoàn. Khi đã biết tốc độ tuần hoàn tối
simulation results. The linearity of the correlation suggests that these thiểu, nhiệm vụ của máy nén trung bình tương ứng có thể dễ dàng
three major parameters are weakly interacting, in line with our thu được bằng cách sử dụng một mối tương quan tuyến tính đơn
intuition of the process. giản.
Ở đây W là nhiệm vụ của máy nén tính bằng kW ở tốc độ tuần hoàn
tối thiểu, VT là thể tích của bể tính bằng nghìn m3, UDL là thông số Acknowledgment
đặc trưng của đường tuần hoàn tính bằng kJ/h-m2-°C và QS là giá trị
gửi -out tỷ lệ trong m3 / h. Mối tương quan này được xây dựng bằng The authors would also like to acknowledge the financial support
cách sử dụng thiết kế giai thừa 3k, đặt tất cả các giá trị tham số khác from the Energy Market Authority of Singapore via grant
thành giá trị của Bảng 1 và có thể thu được hơn 90% tổng độ biến NRF2014EWT-EIRP003-008, and the National University of Singapore
thiên của kết quả mô phỏng. Tính tuyến tính của mối tương quan cho via a seed grant R261-508-001-646/733 for CENGas (Centre of
thấy rằng ba tham số chính này tương tác yếu, phù hợp với trực giác Excellence for Natural Gas).
của chúng ta về quá trình.

6. Conclusions Nhìn nhận

Các tác giả cũng muốn ghi nhận sự hỗ trợ tài chính
The tank modules of commercial process simulators often carry
từ Cơ quan Thị trường Năng lượng Singapore thông
unstated assumptions regarding the physics of the process. The
present study has uncovered several of those unstated assumptions qua khoản tài trợ NRF2014EWT-EIRP003-008 và Đại
and has proposed several simple methods for overcoming those học Quốc gia Singapore thông qua khoản tài trợ hạt
assumptions. The corrected tank model was then extended to include
other aspects of a typical regasification terminal and used to study the
giống R261-508-001-646/733 cho CENGas (Trung
design and op- eration of regasification terminals under several tâm Xuất sắc đối với khí thiên nhiên)
simple but relevant contexts. The following novel heuristics were
uncovered:

A tank aspect ratio of ca. 1 minimizes the BOG generation and the •
average compressor duty.
The pressurization rate of a tank under planned/unplanned shut-
down decreases with increasing liquid level. •
The recirculation rate should be set at its minimum value, which is
dictated by the constraint of one-phase flow in the recirculation •
line. The minimum recirculation rate is unique for a given UDL in
the recirculation line. Once the minimum recirculation rate is •
35
M.S. KHAN et AL. AppliedThermAlEngineeri
ng149(2019)344–353

Appendix I. LNG storage tank roof-, shell-, and bottom-surface design parameters. Heat-transfer coefficients are calculated using the usual
harmonic sum.

2
Thickness (m) Thermal Conductivity (W/m-K) Thermal Resistance (m K/W)

Roof Shell Bottom Roof Shell Vap Shell Liq Bottom

Tankliquid
Tank vaporfilm
film –– –– –– –– – 0.3333 – 0.3333 –
0.0033 –
0.0033
Nickle-steel inner tank 0.02 0.02 0.02 18.3 0.00109 0.00109 0.00109 0.00109
Plywood – – 0.012 0.13 – – – 0.0923
Sand layer – – 0.005 0.15 – – – 0.0333
Glass wool 0.15 0.15 – 0.05 3 3 3 –
Expanded perlite 0.9 0.9 – 0.027 33.33 33.33 33.33 –
Load bearing rigid glass wool – – 1 0.04 – – – 25
Nickle-steel subfloor plate – – 0.02 18.3 – – – 0.00109
Carbon steel liner 0.005 0.005 0.005 54 1.85E−05 1.85E−05 1.85E−05 1.80E−05
Concrete – 1 2 0.1 – 10 10 20
Suspended aluminium deck 0.003 – – 205 1.40E−05 – – –
Air cavity – – – – 3 – – –
Concrete overlay 0.2032 – – 0.1 2.032 – – –
Air film outside tank – – – – 0.4768 0.4768 0.4768 –
Total resistance – – – – 42.173 47.141 46.811 45.131
2
Heat transfer coefficient (W/m K) 0.02371 0.02121 0.02136 0.02216

References storage tanks, Appl. Energy 38 (4) (1991) 263–282.

[25] C. Park, et al., Retrofit design of a boil-off gas handling process in liquefied natural
gas receiving terminals, Energy 44 (1) (2012) 69–78.
[1] I.E. Agency, Key World Energy Statistics (2017) 2017. [26] S. Effendy, et al., Dynamic modelling and optimization of an LNG storage tank in a
[2] Outlook, B.E., edition (2017). UK: BP, 2017. regasification terminal with semi-analytical solutions for N2-free LNG, Comput.
[3] GL, D. Oil and Gas Forecast to 2050, 2017. Chem. Eng. 99 (2017) 40–50.
[4] I.A. Karimi, M.S. Khan, Special Issue on PSE Advances in Natural Gas Value [27] C. Migliore, A. Salehi, V. Vesovic, A non-equilibrium approach to modelling
Chain, ACS Publications, 2018. the weathering of stored Liquefied Natural Gas (LNG), Energy 124 (2017)
[5] H.-M. Chang, A thermodynamic review of cryogenic refrigeration cycles for 684–692.
lique- faction of natural gas, Cryogenics 72 (2015) 127–147. [28] M. Maeda, Y. Shirakawa, Development of fluid dynamics simulations for LNG sto-
[6] M.S. Khan, M. Lee, Design optimization of single mixed refrigerant natural gas li- rage tanks to enhance safety, operational flexibility, and cost performance, in: Conf
quefaction process using the particle swarm paradigm with nonlinear contraints, & Exhibition on Liquefied Natural Gas (lng15), 2007.
Energy 49 (2013) 146–155. [29] K. Koyama, CFD simulation on LNG storage tank to improve safety and reduce
[7] J. Pospíšil, et al., Energy demand of liquefaction and regasification of natural gas cost, Systems Modeling and Simulation, Springer, 2007, pp. 39–43.
and the potential of LNG for operative thermal energy storage, Renew. Sustain. [30] M.S. Zakaria, et al., Computational simulation of boil-off gas formation inside li-
Energy Rev. 99 (2019) 1–15. quefied natural gas tank using evaporation model in ANSYS fluent, Applied
[8] A. Dutta, I.A. Karimi, S. Farooq, Economic feasibility of power generation by re- Mechanics and Materials, Trans Tech Publ, 2013.
covering cold energy during LNG (liquefied natural gas) regasification, ACS Sustain. [31] M.F. Hasan, A.M. Zheng, I. Karimi, Minimizing boil-off losses in liquefied natural
Chem. Eng. 6 (8) (2018) 10687–10695. gas transportation, Ind. Eng. Chem. Res. 48 (21) (2009) 9571–9580.
[9] A. Aspelund, T. Gundersen, J. Myklebust, M.P. Nowak, A. Tomasgard, An [32] HYSYS, A., Version 8.0. Aspen Technology Inc, 2013.
optimi- zation-simulation model for a simple LNG process, Comput. Chem. Eng. [33] M.A. Qyyum, et al., Innovative propane-nitrogen two-phase expander refrigeration
34 (10) (2010) 1606–1617. cycle for energy-efficient and low-global warming potential LNG production, Appl.
[10] M.F. Hasan, M.S. Razib, I.A. Karimi, Optimization of compressor networks in LNG Therm. Eng. 139 (2018) 157–165.
operations, Comput. Aided Chem. Eng. 27 (2009) 1767–1772. [34] M.S. Khan, I. Karimi, M. Lee, Evolution and optimization of the dual mixed re-
[11] M.S. Khan, S. Lee, M. Hasan, M. Lee, Process knowledge based opportunistic opti- frigerant process of natural gas liquefaction, Appl. Therm. Eng. 96 (2016) 320–
mization of the N2-CO2 expander cycle for the economic development of stranded 329.
offshore fields, J. Nat. Gas Sci. Eng. 18 (2014) 263–273. [35] J.H. Park, M.S. Khan, M. Lee, Modified coordinate descent methodology for solving
[12] M. Mehrpooya, A. Vatani, S.A. Moosavian, Introducing a new parameter for eval- process design optimization problems: application to natural gas plant, J. Nat. Gas
uating the degree of integration in cryogenic liquid recovery processes, Chem. Eng. Sci. Eng. 27 (2015) 32–41.
Process: Process Intensif 50 (9) (2011) 916–930. [36] A. Tech, Hysys Customiztion Guide 2011. Support Aspen Tech.
[13] A. Atienza-Márquez, J.C. Bruno, A. Coronas, Cold recovery from LNG-regasification [37] J. Contreras, J.M. Ferrer, Dynamic simulation: a case study, Hydrocarbon Eng. 10
for polygeneration applications, Appl. Therm. Eng. 132 (2018) 463–478. (5) (2005) 103–107.
[14] Đ. Dobrota, B. Lalić, I. Komar, Problem of boil-off in LNG supply chain, Trans. [38] I.M. Trotter, et al., Optimal LNG (liquefied natural gas) regasification scheduling
Maritime Sci. 2 (02) (2013) 91–100. for import terminals with storage, Energy 105 (2016) 80–88.
[15] M. Miana, et al., Calculation models for prediction of Liquefied Natural Gas [39] Y.M. Kurle, S. Wang, Q. Xu, Dynamic simulation of LNG loading, BOG generation,
(LNG) ageing during ship transportation, Appl. Energy 87 (5) (2010) 1687– and BOG recovery at LNG exporting terminals, Comput. Chem. Eng. 97 (2017) 47–
1700. 58.
[16] Y. Lim, et al., Methodology for stable dynamic simulation of a LNG pipe under [40] N. Jang, et al., Dynamic simulation and optimization of the operation of boil-off
two- phase-flow generation, Ind. Eng. Chem. Res. 49 (18) (2010) 8587–8592. gas compressors in a liquefied natural gas gasification plant, Korean J. Chem.
[17] Q.-S. Chen, J. Wegrzyn, V. Prasad, Analysis of temperature and pressure changes Eng. 28
in liquefied natural gas (LNG) cryogenic tanks, Cryogenics 44 (10) (2004) 701– (5) (2011) 1166–1171.
709. [41] H.N. Rao, I.A. Karimi, Optimal design of boil-off gas reliquefaction process in LNG
[18] S. Roh, et al., Numerical study of transient natural convection in a pressurized LNG regasification terminals, Comput. Chem. Eng. (2018).
storage tank, Appl. Therm. Eng. 52 (1) (2013) 209–220. [42] H.N. Rao, I.A. Karimi, S. Farooq, Optimal design of boil-off gas liquefaction in
[19] G.G. Dimopoulos, C.A. Frangopoulos, A dynamic model for liquefied natural LNG regasification terminals, Computer Aided Chemical Engineering, Elsevier,
gas evaporation during marine transportation, Int. J. Thermodyn. 11 (3) (2008) 2018, pp. 2407–2412.
123–131. [43] N. Hisada, M. Sekiguchi, Design and analysis of open rack LNG vaporizer, ASME/
[20] Y. Zongming, et al., Performance improvement of a boil-off gas re-condensation JSME 2004 Pressure Vessels and Piping Conference, American Society of
process with pre-cooling at LNG terminals, Int. J. Thermodyn. (IJoT) 18 (2) (2015) Mechanical Engineers, 2004.
74–80. [44] S.A. Srikanth, S. Narasimhan, S. Narasimhan, Modeling and simulation of unloading
[21] Y. Li, X. Chen, M.-H. Chein, Flexible and cost-effective optimization of BOG (boil-off operations in petroleum product storage terminals, Comput. Chem. Eng. 46 (2012)
gas) recondensation process at LNG receiving terminals, Chem. Eng. Res. Des. 90 59–68.
(10) (2012) 1500–1505. [45] S. Rush, Submerged motor LNG pumps in send-out system service, in: Pump-Zone,
[22] A. Saleem, et al., A CFD simulation study of boiling mechanism and BOG 2004, www.Pump-Zone.Com.
generation in a full-scale LNG storage tank, Comput. Chem. Eng. 115 (2018) 112– [46] R.G. Scurlock, Stratification, Rollover and Handling of LNG, LPG and
120. Other Cryogenic Liquid Mixtures, Springer, 2015.
[23] C. Migliore, C. Tubilleja, V. Vesovic, Weathering prediction model for stored li- [47] L.A. Pellegrini, et al., LNG technology: the weathering in above-ground storage
quefied natural gas (LNG), J. Nat. Gas Sci. Eng. 26 (2015) 570–580. tanks, Ind. Eng. Chem. Res. 53 (10) (2014) 3931–3937.
[24] P. Jourda, S. Probert, Heat-transfer considerations for large liquefied-natural-gas

35