1. Introduction
As highlighted in the European Green Deal, transport accounts for a quarter of the European Union’s (EU) greenhouse gas emissions, and this figure is still growing. To achieve climate neutrality in the sense of the Paris Climate Agreement, current transport emissions would have to be reduced by 90% by 2050 [
1].
The use of alternative fuels, such as Liquid Natural Gas (LNG), could contribute to reducing greenhouse gas emissions. As outlined in DNV GL’s Technology Outlook 2025 [
2], global demand for gas as an environmentally friendly alternative to conventional marine fuels has more than doubled in the past 30 years and it will continue to increase [
3].
Global natural gas resources are growing, and their geographic diversity is increasing [
4]. The LNG shipping industry kept pace with this growth, adding 57 carriers and 4 Floating Storage and Regasification Units (FSRUs) to a total of 641 active vessels by the end of April 2022, including patrol vessels, car ferries, cruise vessels, bulk carriers, container vessels and other types of LNG vessels [
3,
5]. The active fleet includes 45 FSRUs and 5 Floating Storage Units (FSUs), demonstrating the continued interest in flexible solutions to enable markets to start importing LNG or increase their LNG imports as energy demand grows. This illustrates shipowners’ expectations that the LNG trade will continue to grow in line with the increase in liquefaction capacities in the coming years.
For markets to grow, LNG and electricity supply infrastructures and their services need to be open, transparent, and interoperable. By 2025, the EU should have completed the backbone of refueling infrastructure (up to EUR 945 million in the Trans-European Transport Network (TEN-T) Core Network Corridor seaports by 2025), providing full coverage of the TEN-T core network corridors [
6]. Users also need to be able to use the whole transport network in a smooth and easy way—that is why synergies with other modes of transport should be examined as well.
LNG is natural gas converted into a liquid form. Condensation occurs at −162 °C, reducing its volume by 600 times (in comparison to the gaseous state), which makes LNG cost-efficient to transport over long distances.
Natural gas consists mainly of methane and compared to other fuels has very low amounts of sulphur oxides, nitrogen oxygens, and particulate matter when combusted [
7]. LNG therefore not only easily complies with the sulphur emission caps but also meets the requirements for upcoming regulations on NOx reduction.
In addition, LNG produces less than 10% particulate matter and 50% fewer greenhouse gases than coal when used to generate electricity, 21% less than fuel oil in transportation, and has over 95% efficiency when used to heat homes [
8]. In addition, according to DNV GL, ships powered by LNG can reduce CO
2 emissions by 15% to 20%.
Floating Storage and Regasification Units (FSRUs) are one of the most commonly used vessel types in the global ship fleet due to the possibility of storage, reloading to another ship, and regasifying it for re-injection into the natural gas grid. Commonly, an FRSU consists of special isolated tanks for LNG storage and handling inside the ship, and on the deck, there is regasification equipment, which converts LNG into a gaseous state for end users [
9,
10]. Observing LNG technological processes in FSRUs, the main focus is on LNG unloading, loading, and regasification modes [
10]. So, it is important to control system parameters for reliable technological processes, such as hydrostatic pressure in the tanks, vapor pressure, LNG density, LNG temperature, and phase changes between liquid and gas states. Furthermore, pressure monitoring is important to control during transit in port and bunkering to prevent the pressure in the tanks from exceeding the tank design pressure. During loading and unloading modes, a detailed calculation for the maximum allowable LNG level in the filling tanks should be provided to prevent overfilling [
11]. One of the main challenging processes of FSRUs is regasification. For the regasification process, low pressure (8–10 bar is commonly when the selected pump starts) and high-pressure pumps (80–120 bar) are used. At first, LNG is pumped by a low-pressure pump from the storage tank to the buffer tank or the recondenser unit. From there, the LNG is pumped by high pressure to the vaporizer.
Another challenge is the many pipes and pipelines of the LNG regasification pipeline system. One of the main tasks for engineers is to design a cryogenic pipeline with minimum length, ensuring the minimum energy losses related to the heating of liquid by heat input and hydraulic losses [
12]. It is important to protect the work of the system from technical problems by technical maintenance of regasification equipment, including observation of technological parameters in pipelines. One of the solutions could be the LNG flow observation model in the regasification system.
Most mathematical models are described for natural gas flow in the pipelines. Generally, the flow of natural gas through high-pressure pipelines is governed by the partial differential equation (PDE) of mass and momentum conservation in conjunction with the equation of state (EOS) and is used with different methods to solve [
13]. Only several studies have been conducted for modeling LNG technological processes in pipelines, and more have been performed for LNG regasification equipment, for example in [
14,
15,
16,
17]. The authors of study [
18] created a model based on thermal and hydraulic calculations of an underground cryogenic pipeline. This model allowed estimating the main production parameters such as pressure drop and temperature of LNG during its flow through the pipeline underground. In addition, this model permits calculations without initial parameters such as inlet pressure and diameter. One of the most important conditions is that the temperature should be below the boiling point of LNG at the given pressure, and the operating pressure should exceed the saturation pressure of LNG by at least 0.5 MPa. Another analytical model for describing the LNG flow in the pipeline was described in [
19], which allows the determination of parameters such as pressure, temperature, density, and viscosity for a given composition as a function of pipeline length.
This paper presents the mathematical model for a real LNG transportation system at cryogenic temperature (110 K) starting from LNG pumping from the LNG tank to the Suction Drum and pumping to the LNG vaporizer of the regasification unit. These processes are the most important for reliable further regasification equipment work. The model investigates hydrodynamic and dynamic high-velocity processes in the unsteady fluid flow when the LNG is transferred to the regasification unit and allows observing pressure changes, LNG velocity and the speed of sound in the LNG at the piping points in the system during system start-up (until 5 s). The universality of this model is that there are a lot of variables in the system to determine and it allows for evaluation by the different technologic parameters of changing lengths and diameters of pipes, liquid level, volume in the tanks, and electric motors parameters, or it uses another fluid in cases when fluid parameters are known. Of special interest to this model is that characteristics, pressure losses, and hydrodynamic processes of multistage centrifugal pumps are verified when operating in different technological regimes when LNG flow enters the LNG regasification unit. Furthermore, the work of the electric motor of the transportation pump of the main LNG tank is a part of the dynamic process observation in this model. The model is validated using technical data that characterize real technological regimes of regasification in the Floating Storage and Regasification Unit.
2. Research Methodology
2.1. Symbols and Notations
Table 1 shows the symbols used in this work.
2.2. Description of LNG Transportation System
The basic components of the LNG transportation system in the FSRU are presented in
Figure 1. The technological process starts when a submerged cargo pump from the LNG Tank transports the LNG at high pressure to a vapor-liquid separator (Suction Drum) on deck. The pipeline system from the LNG tank to the Suction Drum is divided into different sections of diameters. Different diameters of pipelines are used to minimize pressure losses in the system, in that case decreasing LNG flow in separate sections. All the pipelines of the system are marked with the letter L with a number, which is essential for calculation in further mathematical modeling.
When LNG has been transported to the Suction Drum, the LNG starts mixing with the Boil of Gas (BOG). The Suction Drum has two functions to maintain constant LNG pressure, which is required to operate the compression pumps (High-Pressure Booster Pumps), and to condensate BOG from the LNG tanks. In this research, BOG mixing between LNG is not analyzed and so BOG is not marked in the simplified scheme. From the Suction Drum, the LNG is transported to the regasification equipment by a high-pressure intensifier pump. Usually, there are several regasification trains, one of which serves as a reserve. In this model, three parallel regasification trains are considered as an example (a train is a typical regasification line consisting of the same regasification equipment as the high-pressure booster pump, the BOG recondenser, the LNG vaporizer, etc.). Firstly, LNG flow goes to the Pipeline Unit to separate LNG flow, and then it is divided into three LNG flows, which are connected to regasification lines. Three regasification lines (regasification lines indicated as L3, L4, and L5) are always in working modes. After introducing the LNG transportation system, the mathematical model of the LNG transportation system including descriptions of LNG equipment is presented in the following sections.
2.3. Model Assumptions and Evaluations
As previously described, the LNG transportation system, which is operating in different technological regimes, and a full mathematical model are built, involving individual mathematical models of equipment for technological parameter observation (
Section 2.4,
Section 2.5 and
Section 2.6). For the whole system mathematical model, LNG flow in pipes is connected to system devices (tanks and pumps) using boundary conditions (
Section 2.7). The main goal of this complex model is that by changing one of the technological parameters, it allows us to observe the impact to the whole system. The principal scheme of the mathematical model of the total system is introduced below in
Figure 2.
The key assumptions used for mathematical models are as follows:
The LNG tank and the Suction Drum vapor pressures (pv1 and pv2) are assumed to be constant. From other scientific studies, it was verified that the vapor pressure changes slowly. In real systems, it changes during the day or week, etc. In our model, it is important to observe the parameters of the system when the LNG transportation system starts working to verify that the model works properly at the beginning. The solution time of the model is chosen to be 5 s, so the vapor pressures do not change drastically at this time.
LNG temperatures (Tl1 and Tl2) are assumed to be constant due to the small heat transfer between walls and LNG in the LNG tank and the Suction Drum. For this reason, the system must be fully insulated, and the heat transfer between walls and LNG in the LNG tank and the Suction Drum is therefore very low. Otherwise, the system will not work because LNG will start evaporating and will not characterize real-life regasification cases in the FSRU.
The height of the liquid level of LNG in the tank (Htank1(t)) and the Suction Drum (Htank2(t)) varies in the LNG transportation system. The height of the liquid level of LNG in the tank strongly affects hydrodynamic processes in the system, especially in the pipeline L1, when LNG starts transferring from one tank (LNG tank) to another (Suction Drum). Additionally, the height of the liquid level is a good indicator of start-process and end time.
LNG composition is assumed to be 91.798 mol% methane, 5.698 mol% ethane, 1.303 mol% propane, 0.396 mol% n-butane and 0.805 mol% nitrogen. The composition has an impact on BOG generation, density changes, and evaporation rate. So, the composition was selected from real-life regasification cases in FSRUs (LNG Terminal of Klaipeda City, Lithuania), which is important for end users.
2.4. Description and Mathematical Model of LNG Tank and Suction Drum
The model configuration starts with hydrodynamics processes such as LNG level observations in the LNG tank (Htank1(t)) and the Suction Drum (Htank2(t)) when the first pump (the submerged motor cargo pump) starts working. The FSRU of this model shall include double-insulated tanks with inner shells. The containment system functions are to contain LNG cargo at approximately a temperature of −163 °C and insulate the cargo from the hull structure. The simplified scheme of one of the LNG tanks is introduced in
Figure 3. Transportation of LNG from the LNG tank to the Suction Drum is part of the regasification process. The Suction Drum is a vertical vessel to balance pressure for further regasification steps. Additionally, the Suction Drum operates as a buffer tank for the following Higher-Pressure Booster Pumps, ensuring that there is always sufficient LNG available to the pumps when the pumps are operating at maximum capacity. The liquid-level observation is shown in
Figure 3 and
Figure 4.
Hence, the operation of the LNG tank and the pump is described by the system of differential equations for estimating the pressure and the LNG level change in the LNG tank and the Suction Drum, according to Equations (2) and (5):
Level of LNG in the tank:
- (b)
The Suction Drum
Level of LNG in the Suction Drum:
where:
—liquid level in the tank by the time changing in the LNG tank, m;
—height between the pump and the bottom of the LNG tank, m;
and
—vapor pressure of BOG in the LNG tank and the Suction Drum, Pa; g—gravitational acceleration, m/s;
and
—cross- sectional areas of the LNG liquid surface in the LNG tank and in the Suction Drum, m
2;
and
—cross-section areas of the LNG tank and the Suction Drum, m
2;
and
—LNG volume flow of the pipeline (L1) inlet and outlet, m
3/s;
— LNG volume flow of the pipeline (L2) inlet, m
3/s;
LNG density, kg/m
3, determined according to the Enhanced Revised Klosek–McKinley method (ERKM) method by Equation (6). Commonly, this method is popular in the gas industry and is used under cryogenic conditions [
20,
21,
22]. Additionally, the ERKM method is valid under pressures up to 10 MPa and the temperature varies from 115 K to 135 K [
23].
where:
—constant and it is converted from MPa to Pa [
23].
Molecular mass of the mixture is expressed by following equation:
where:
—molecular mass of mixture, kg/kmol; x
i—molar fraction of constituent i.
where:
—molar volume of component i at the temperature of LNG, m
3;
—molar volume of constituent i at the temperature of LNG, m
3;
—molar fraction of constituent i, mol;
,
—volume correction factors, mol/l; 0.0425—constant for the LNG with the content of nitrogen or butane content less than 4% [
20,
24].
Pressure correction for partial pressure (saturation) is described Equation (9):
where:
and
are determined according to [
23].
Pseudo critical temperature:
where:
—molar fraction of constituent, mol;
—critical temperature of constituent i, K.
2.5. Description and Characteristics of LNG Submerged Cargo Pump and Electric Motor
To analyze the dynamic process of the LNG transportation system, the model is performed separately for the electric motor (EM1) of the transportation pump in the LNG tank. The submerged motor cargo pump (SMR) operates in the system of this model as in the real-case operation of the FSRU. The submerged motor cargo pump [
25] casing and outer frame of the motor are combined with a flange joint as one unit. One of the main advantages of this type of pump is that it can shift from maximum flow to minimum flow smoothly following a stable head/capacity curve, even when shutting down, to ensure stable operation without noise and vibration [
26]. The pump has a submerged motor with windings, which is cooled by the pumped LNG. So, LNG is used to lubricate and cool the pump and the motor bearings [
27]. The SMR-type pump could reach the maximum flow of 520 m
3/h, which can be reduced to 220 m
3/h according to clients’ LNG regasification demands.
At the beginning of the work of the SMR pump, the torque of the electric motors is described by the differential equation [
28]:
where
—parameters of the electric motor;
—angular velocity of the electric motor, rad/s.
Equation of rotation of the electric motor’s rotor in the SMR is equal to:
where
—mass inertia moment of the rotor of electric motor and pump’s rotor;
—damping coefficient;
—load torque of the electric motor’s rotor.
where
—total flow of SMR pump, m
3/s;
—coefficient of SMR pump, (COP);
—is total SMR pump characteristic;
—LNG flow in the first section of pipeline, m
3/s.
where
(
)—total SMR characteristic at an angular velocity
.
Total SMR characteristic is described as a fifth order polynomial with a correlation coefficient of 0.9998 which is determined from points of technical FSRU sheets:
2.6. Description and Mathematical Model of High-Pressure LNG Booster Pumps
Basically, the described electric motor model (EM1) in
Section 2.5 can be used for the rest electric motors (EM2, EM3, EM4), though another type of pump is used for LNG transfer. The last step of the LNG transportation system is transferring LNG by a high-pressure booster pump to LNG vaporizers. It is the vertical, cryogenic centrifugal pump. In the FSRU, there are three trains, each consisting of two high-pressure pumps.
Observing the pressure and LNG flow changes through passages of the pump, the model is described by differential equations. These variables have an effect on operational efficiency. The model is shown in
Figure 5 below.
Mass flow of LNG (
), which enters the passage between the pipe and the first stage of pump’s joint, is determined according to differential Equation (16):
where
—section area between point N and
, m
2;
—length between point N and
, m;
—pressure of inlet pipe, Pa;
—pressure of the passage between the pipe and the first stage of pump’s joint, Pa;
—pressure losses between point N and
, Pa;
—sign function, which indicates LNG mass flow direction.
Pressure losses (
) are calculated by Equation (17):
Hydraulic losses coefficient (
) is determined by Equation (18):
To calculate hydraulic losses, the
could be determined according to these conditions:
In that case, when LNG flow is turbulent, it could be solved by the Colebrook–White nonlinear algebraic equation
F(
) = 0 [
29,
30].
The pressure of the passage between the pipe and the first stage of pump’s joint
could be determined by Equation (20):
where
—bulk modulus of elasticity of LNG, Pa;
—LNG mass flow between
and
, kg/s;
—LNG density in inlet pipe, kg/m
3;
—LNG density in the passage between the pipe and the first stage of pump’s joint, kg/m
3.
LNG mass flow between
and
(
) is determined by differential Equation (21) below:
where
is determined according to Equation (17).
Having in mind that is the multistage pump, the pump pressure is determined by verifying all stages of the pump, including the pump’s characteristics, with Equation (24) below:
Pump pressure of the one of the stages is determined by Equation (23):
where
;
—one of the stages of pump’s characteristics, which is provided below:
Pressure (
) of the passage between the first stage of the pump’s joint and the outlet pipe is determined by:
where
—bulk modulus of elasticity of LNG, Pa;
—LNG mass flow between
and
, kg/s;
—LNG density of Vp,
1, kg/m
3;
—LNG density of Vp,
2, kg/m
3,
—mass flow of LNG from outlet of pump to the outlet pipe, kg/s.
The mass flow of LNG
from outlet of pump to the outlet pipe is determined by Equation (26):
where
—sectional area between
and point 1, m
2;
—length between
and point 1, m;
—pressure of the passage between the pipe and the first stage of pump’s joint, Pa;
—pressure of the passage between the first stage of the pump’s joint and the outlet pipe, Pa;
—pressure losses between
and point 1, Pa;
—sign function, which indicates LNG mass flow direction.
2.7. System of Equations of the Pipelines Using Method of Characteristics
In the last part of the model, the one-dimensional non-stationary movement of LNG in the variable cross-sectional A(x) area of the pipeline is studied, i.e., when the LNG velocity vector is directed along the axis of a pipeline, and pressure, velocity, and temperature changes in time, and along the axis of the pipeline (coordinate). Moreover, the cross-sectional area A(x) of the pipeline varies and depends on the x coordinates. Furthermore, the flow of LNG in pipes is described by the differential equations of flow integrity and motion with partial derivatives, which are solved by the method of characteristics, and the unknown variables (velocity of LNG, pressure) are determined at every point of the pipeline.
The equation of LNG continuity can be written in a differential form as follows [
31,
32,
33,
34,
35]:
where
−density and velocity of LNG, kg/m
3 and m/s;
−cross-sectional area of a pipeline, m
2.
Equation of liquid flow impulse (momentum):
where p—LNG pressure, Pa;
—the perimeter of cross-section of the pipeline, m;
—tangential fluid stress on the inner surface of the pipeline, Pa;
where p—LNG pressure, Pa;
−the perimeter of a cross-section of the pipeline, m;
—the acceleration along the x axis, m
2/s.
where
—angle between horizontal axis and center line of pipe.
The system of Equations (27) and (28) could be expressed according to second-order quasi-linear differential equations:
where
,
—matrices and
—a vector which depends on t, x and elements of a vector
Matrices are described:
Vector is described:
where c—speed of sound in the LNG, m/s.
By equating the determinant of the matrix
to zero, which allows determining the
derivative and determines characteristic direction. This equation has two various real roots:
Such liquid movement is characterized by the wave of increased and reduced pressure, which spreads from the place of change in each pressure vibration cross-section and in the deformation of high-pressure pipeline walls [
36,
37]. The speed of sound
c in the LNG, which is stored in the pipeline, is equal to:
where
—bulk modulus of elasticity of liquid, Pa;
;
—density of the LNG; E—modulus of elasticity of a pipeline, Pa;
—internal diameter of a pipeline, m; e − thickness of a wall of a pipeline, m.
The main idea of the characteristics method (MOC) is the fact that the unknown variable velocity (
) and the pressure of the liquid (
) at an instant moment time t + ∆t are determined according to the parameters that time t + ∆t [
29,
38]. The MOC scheme is shown in
Figure 6.
Pressure and velocity at point D at that moment of time are determined from a nonlinear algebraic equation system:
where
;
.
The nonlinear algebraic Equations (36) and (37) are solved by the Newton and Raphson method [
39]:
where
— Jacobi matrix;
— increment vector of unknown variables; i—number of iterations.
To ensure the stability of the solution, the Courant–Friedrichs–Lewy condition shall be fulfilled [
40]:
where
is Courant number.
At points L and R, variables p and
are defined using a system of two nonlinear algebraic equations with unknown p
L, v
L, and p
R, v
R:
The boundary conditions are used to connect the pipeline with the equipment of the system (LNG tank, Suction Drum, pumps, etc.), which are described below:
Boundary condition: when the point is known
pressure is known variable
, then LNG flow velocity is determined by formula:
Boundary condition: when the point is known
, LNG flow velocity is known variable
, then pressure is determined by formula:
Boundary condition: when the point is known
, pressure is known variable
, then LNG flow velocity is determined by formula:
Boundary condition: when the point is known
, NG flow velocity is known variable
, then pressure is determined by formula:
When boundary conditions are set according to Equations (45) and (48), the nonlinear algebraic equations are solved by the Newton and Raphson method. Then, the Pipeline Unit (Volume—V3) is used to determine one of the last unknown variables—pressures (
Figure 7).
At the Pipeline Unit (V3), the LNG flow balance equation is written to determine the pressure:
where
,
,
,
,
—LNG volume flow, m
3/s;
–losses of LNG volume flow, m
3/s.
Velocities of LNG in Pipeline Unit sections of pipes L2, N and L3,1, L4,1, and L5,1 are determined when pressure is known in the Pipeline Unit (V3), using boundary conditions (45), (47) by:
For mathematical model construction, there are two important conditions:
- (a)
When pipes (L3, L4, L5) are closed until HP booster pump 1 starts working:
- (b)
When HP booster pump 1 starts working and pipes start opening:
Finally, the model is finished in the LNG transfer to the BOG recondensers. The velocity of LNG in Pipeline Unit sections of pipes L6, L7, and L8 are determined when pressure is known, using boundary conditions (47) by:
The pressure at the end of the pipes L6, L7 and L8 is determined according to the expression below:
where
—constant pressure in the pipes.
Hydraulic losses coefficient:
where pipe index—k = 6, 7, 8 and
—cross-section area of pipe k, m
2;
—cross-section area of channel of k-th BOG recondenser, m
2;
—number of channels of k-th BOG recondenser.
All in all, the theoretical mathematical model was built for the full LNG transportation system in the FSRU in this section. A general mathematical model of the system is released with a special software package using Fortran programing language and MATLAB R28a.
3. Results of LNG Transportation System Dynamics and Hydrodynamics Processes
A mathematical model of the FSRU LNG transportation system was created to solve total system equations with different boundary conditions and determine the final system parameters, such as pressure and velocity. In the mathematical model, the solution time step is 10
−5 s, and the discretization step of the pipes is 0.1 m. The total solution time is chosen as 5 s. The temperature of LNG in the LNG tank, the Suction Drum, and the pipelines is 110 K. The general conditions of FSRU equipment when the LNG transportation system starts working are summarized in
Table 2.
Table 3 lists different diameters and lengths or angles of inclination of pipes and pipelines considered in the simulation of LNG flow. Before starting the mathematical model simulation, it is essential to mention that in this model, that LNG flow transfers at different diameters and lengths or tilt angles of pipes and pipelines (described in
Table 2), presenting better results and the universality of the model.
The model simulation began when the SMR pump with an electric motor (EM1) started working.
Figure 8 represents the angular velocity of electric motors and how fast the orientation of pumps changes with the solution time. Additionally, it firstly shows that that EM1 of the SMR pump starts working, indicated by the blue line. Then the second electric motor of HP Booster Pump 1 (EM2) turns on after 0.5 s (red-color curve), the electric motor of HP Booster Pump 2 (EM3) turns on after 1.0 s (yellow-color curve), and the last electric motor of HP Booster Pump 3 (EM4) starts after 1.5 s (purple-color curve).
In general, the observation of pressure and velocity simulated by the mathematical model in the first LNG pipeline of the LNG transportation system agrees well with the technical data.
Figure 9a shows the initial pressure curve along the pipeline L1 (from the LNG tank to the Suction Drum) after the start of the system, and
Figure 9b shows pressure change for 5 s in the first pipeline. Furthermore, the observation of pressure and velocity was selected in the different pipeline distances, using different color lines to verify pipeline geometry. After 4 s, it can be seen that pressure and velocity start to stabilize (
Figure 9d).
After solving (1)–(5) equations of the LNG level in the LNG transportation system, after the solution time run, the results showed that level of LNG decreases in the LNG tank (
Figure 10a) and it starts slowly filling with LNG in the Suction Drum (
Figure 10b). The filling level in the Suction Drum increases by 0.040 m in the first second of the simulation time before the filling level continuously decreases by approximately 0.030 m per second, as in the main LNG tank (about 0.025 m per 1.7 s simulation time).
Figure 11 indicates the results of simulated pressure in the second pipeline (L2). Observing hydrodynamic processes verifies that the pressure and velocity of L2 depend on the start-up time moment of HP Booster Pumps.
Figure 11a, during 0.4–0.5 s, shows pressure waves in the LNG, indicating that pipes L3, L4, and L5 are closed and that the systems work only in the Suction Drum. Comparing
Figure 10b (LNG level in the Suction Drum) and
Figure 11a, it can be seen that a moment after 0.8 s that HP Booster Pump 1 starts working, LNG is transported by it in the further system (pipeline L2, Pipeline Unit V3). However, it is verified that the HP Booster Pumps 2 and 3 do not work.
Meanwhile,
Figure 11b shows pressure pulsation in the pipeline (L2) when all the HP Booster Pumps start working. The diagram of
Figure 11b estimates that pressure in pipeline 2 changes at a frequency of about 160 Hz and is influenced by HP Booster Electric Motor 1 and pressure wave propagation velocity in part of the system (L2, L3, L6, and other L7, L8 pipes are closed). The pressure pulsation indicated that the system was ready to work for further regasification process. The pulsations of pressure are important to observe because they could load pipes, i.e., the stresses in the pipe walls start increasing, which affects the creation of cracks and the occurrence of material defects inside layers of contact [
41].
Figure 12 shows pressure changes in the Pipeline Unit (Volume—V3).
Figure 12a shows when HP Booster Pumps turn on (peaks): after 0.5 s—the first HP booster pump; after 1 s—the second HP booster pump; and after 1.5 s—the third HP booster pump. The pumps’ work indicates that the further step would be the LNG regasification process.
Figure 12b shows that pressure changes are slightly similar to those in the pipeline (L2). As
Figure 12b shows, there is pressure pulsation in the whole system.
Figure 13 shows the mass flow of the LNG of HP Booster Pumps to the outlet pipes.
Figure 13a introduces the beginning of HP Booster Pump work. At the end of the solution time (5 s), the flow rate is approximately 502 m
3/h for HP Booster Pumps (
Figure 13b). The estimated mass flow is based on Equations (16), (21) and (27), which are related to pressure and velocity changes. Additionally, it is important to mention that the analyzed LNG transportation system is completely closed and constantly operating (isolated from the environment). So, the average values of some of the main technical parameters could only be observed over a certain period of time. For example, to compare the flow: in our theoretical studies of the real system, it has an average value of 502 m
3/h, compared to the real system according to the technical parameters, the value is 520 m
3/h. A slight difference in LNG flow is because we preferred to study the hydrodynamic processes occurring in pipes L6, L7 and L8 under different boundary conditions, i.e., on the number of tubes of varying heat exchangers (which leads to different pressure losses).
To finalize the LNG transportation process, the last parameters of hydrodynamics were verified (pressure and velocity).
Figure 14 shows LNG pressure and velocity changes in pipes 6, 7 and 8 when HP Booster Pumps turn on at different times and start working. After 5 s, we observed pulsations pressure and velocity (
Figure 14b). The many pulsation peaks with different amplitudes show an unsteady flow of LNG. Due to the unsteady LNG flow, LNG flow monitoring is important for further analysis of LNG vaporization in heat exchangers. Vaporization efficiency depends on pressure, temperature, mass flow, and volume fraction changes. From the research of [
31], we take that an FSRU is a complicated system with various heat transfer mechanisms and has actual demands of improving efficiency and compactness.
Figure 15 introduces pressure and velocity changes in the sixth pipe.
Figure 15a shows that the HP Booster Pump starts working after 0.5 s, so the pressure suddenly increases. Hydrodynamic processes did not stabilize because of pressure waves in the LNG. It could be seen from
Figure 15b that pressure pulsation is smaller at the end of the pipe because of the boundary conditions input, and that pressure depends on velocity (see green line in
Figure 15b). However, at other points of the pipe, pressure pulsation is bigger (other color lines) because of the pressure losses along the L6 pipe. Furthermore, it is seen from
Figure 15c that velocity varies similarly to pressure in the pipe. At the end of solution time (5 s), the LNG average velocity reaches approximately 4.37 m/s (
Figure 14d).
The variation of hydrodynamic processes (pressure and velocity) at any time is presented at the system’s characteristic points, which indicates curves of different colors in
Figure 9,
Figure 11, and
Figure 15. These curves show different pressure and velocity changes at other points and at each time. As mentioned in the introduction, these pressure and velocity changes are important for safe and reliable FRSU operations without compliances, especially in holding mode, LNG transporting between tanks, the regasification process, and the ship-to-ship procedure.
Figure 16 introduces hydraulic energy losses in the first (
Figure 16a) and the rest of the other pipes of the LNG transportation system (
Figure 16b). Hydraulic energy losses (power of hydraulic friction forces) depend on LNG velocity, length and inner diameter of the pipe, and roughness of the inside surface of the pipe. So, it can be seen from
Figure 16a that the highest energy losses are in the first pipeline due to the different geometric parameters of the pipeline—diameter, tilt angle, and length.
After 5 s, hydraulic energy losses in the pipe L1 are equal to 50 kW, respectively, and total energy equals 53 kW. These hydraulic energy losses in all pipes contain 1.7% of the whole system power (total power of the electric motors is 3132 kW). In
Figure 16b, we can observe hydraulic energy losses in different pipes or pipelines, which is important for another regasification process. This mathematical model of total system would help control energy losses in the FSRU (LNG Terminal of Klaipeda City, Lithuania) in various technological modes.
The model validation was performed using technical data in daily intervals, as is typical for FSRUs. For the individual components of the model, such as the pumps, combined with the pump and electrical characteristics from manufacturers, we were able to generate curves at various angular velocities, as shown in
Figure 17.
Figure 17 shows how the SMR pump generates pressure as a function of our angular velocity. For comparison, the lowest blue curve from
Figure 8 shows the operating points, i.e., the dependence of the pressure on the flow rate of the first line (L1,1) at the first point. The electric motor is connected to the pump, and using the manufacturer’s pump curve p(Q) at different motor angular velocity
, we obtain the pump curve
at the corresponding motor angular velocity
. This gives the pressure generated by the real pump in the system.
In the model, all devices connect to one system, so when changing pressure, speed, or other parameters, other parameters of the system change as well (the system works automatically). According to the manufacturer’s technical data, the model does not work in any mode other than the one set. Therefore, we can assume the validity of the model in the simulated time range in seconds, although naturally, for an FSRU there are no technical values in this time scale. The uniqueness of our model is that the parameters and the linking algorithm are fixed at a 5 s interval at the start time for each point of the system.
4. Conclusions
This research investigates the hydrodynamics and dynamics processes in the real LNG transportation system of FSRUs (LNG Terminal of Klaipeda City) in system start-up (until 5 s) mode. The created total mathematical model of the LNG transportation system consists of electric motors, pumps, tanks, and pipelines, enabling the evaluation of electric motors, pumps’ dynamics characteristics, and speed of sound in the LNG when LNG flow enters the LNG regasification unit.
Using this model, it is possible to study various technological regimes (for example, starting pumps at different times and varying the initial height of filling tanks with liquid). The model simulation began when the SMR pump with an electric motor (EM1) started working. The results of the study of hydrodynamic processes are obtained when the second electric motor of HP Booster Pump 1 (EM2) is started after 0.5 s, the third electric motor of HP Booster Pump 2 (EM3) is started after 1.0 s, and the fourth electric motor of HP Booster Pump 3 (EM4) is started after 1.5 s, after the first pump is started.
The hydrodynamic processes in the considered system have not yet been steady after 5 s. During the test time of 5 s, the pressures and velocities in pipe 6 change at a frequency of about 160 Hz, which depends on the rotation of electric motors, LNG physical and mechanical properties, and the materials of pipes.
We can accurately characterize the flow parameters (velocity, pressure) entering the additional regasification equipment (LNG vaporizers) with this model. For example, our theoretical studies of the real system obtained an average value of 502 m3/h, compared to the real system according to the technical parameters, with a value of 520 m3/h. The slight difference in LNG flow is because we additionally preferred to study the hydrodynamic processes occurring in pipes L6, L7 and L8 under different boundary conditions, i.e., the number of tubes of different heat exchangers (which leads to different pressure losses).
Hydraulic energy losses (power of hydraulic friction forces) depend on LNG velocity, length and inner diameter of the pipe, and roughness of the inside surface of the pipe. It was found that hydraulic energy losses in all pipes were equal to 53 kW. These hydraulic energy losses in all pipes contain 1.7% of the whole system power (the total power of the electric motors is 3132 kW). It was determined that a slight variation in LNG temperature does not significantly affect the hydrodynamic processes in the system.
In general, the pressure and velocity results of the mathematical model agree well with the technical data, which also confirms that the mathematical model is suitable to characterize real technological regimes of regasification in Floating Storage and Regasification Systems. The model found that the LNG flow in this system is unsteady and should be monitored for further regasification processes.