Next Article in Journal
Design and Implementation of 3-D Measurement Method for Container Handling Target
Previous Article in Journal
Green Ports Analysis Using an End-to-End Tool Application in the Fishing Port of Vigo
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling and Simulation Hydrodynamics Processes in Liquefied Natural Gas Transportation Systems

by
Marijonas Bogdevicius
1,
Vigaile Semaskaite
1,
Tatjana Paulauskiene
2,*,
Jochen Uebe
2 and
Algimantas Danilevicius
1
1
Department of Mobile Machinery and Railway Transport, Faculty of Transport Engineering, Vilnius Gediminas Technical University, Plytines St. 27, 10105 Vilnius, Lithuania
2
Engineering Department, Faculty of Marine Technology and Natural Sciences, Klaipeda University, H. Manto 84, 92294 Klaipeda, Lithuania
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(12), 1960; https://doi.org/10.3390/jmse10121960
Submission received: 8 November 2022 / Revised: 30 November 2022 / Accepted: 7 December 2022 / Published: 9 December 2022
(This article belongs to the Section Ocean Engineering)

Abstract

:
Global natural gas resources are growing and are increasingly geographically diverse. A Floating Storage and Regasification Unit (FSRU) is 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. It is important to control system parameters for reliable technological processes such as tank hydrostatic pressure, vapor pressure, LNG density, LNG temperature, and phase changes between liquid and gas states. Additionally, 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. In this research study, a comprehensive hydrodynamic model for an LNG storage tank in a real-life regasification terminal (Floating Storage and Regasification Unit, LNG Terminal of Klaipeda City, Lithuania), operating in transportation mode to the regasification unit, was created. For this research, LNG is investigated as a compressible liquid and the speed of sound in LNG is evaluated. A complex mathematical model of the system allows the analysis of high-speed hydrodynamic and dynamic processes at cryogenic temperature (110 K) and evaluates the geometric parameters (tank geometry, electric motors and pumps, pipe geometric parameters, and roughness of internal surfaces) and the characteristics of pumps and electric motors. The complex mathematical model of the system was implemented using Fortran programing language and MATLAB R28a. It determined the parameters (pressure, velocity, liquid level of LNG in the tanks, electric motor angular velocity, torques, hydraulic energy losses, etc.) of the system during its start-up mode (until 5 s). It was found that hydraulic energy losses in all pipes contain 1.7% of the whole system power (the total power of the electric motors is 3132 kW). In case of increasing energy costs, this model could be used to control energy losses during the operation of the FSRU in various technological modes.

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 CO2 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):
(a)
The LNG tank:
p 1 , in = ( H tank1 ( t ) H 1 , 0 ) · ρ l · g + pv 1
Level of LNG in the tank:
H ˙ tank1 ( t ) = A 1 A tank1 Q 1 , in
(b)
The Suction Drum
p 2 , in = ( H tank2 ( t ) H tank2 , in ) · ρ l · g + pv 2
p 2 , out = ( H tank2 ( t ) H tank2 , out ) · ρ l · g + pv 2
Level of LNG in the Suction Drum:
H ˙ tank2 = 1 A tank2 [ Q 2 , in Q L 1 , out ]
where:   H tan k 1 ( t ) —liquid level in the tank by the time changing in the LNG tank, m; H tan k 1 , 0 —height between the pump and the bottom of the LNG tank, m; pv 1 and pv 2 —vapor pressure of BOG in the LNG tank and the Suction Drum, Pa; g—gravitational acceleration, m/s; A 1 and A 2 —cross- sectional areas of the LNG liquid surface in the LNG tank and in the Suction Drum, m2; A tan k 1 and A tan k 2 —cross-section areas of the LNG tank and the Suction Drum, m2; Q 1 , in and Q L 1 , out —LNG volume flow of the pipeline (L1) inlet and outlet, m3/s;   Q 2 , in   — LNG volume flow of the pipeline (L2) inlet, m3/s; ρ l LNG density, kg/m3, 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].
ρ l = M mix V mix · [ 1 + 406 · ( p l p s , corr ) · ( T pc T pc T ) 1 , 77 ]
where: 406 —constant and it is converted from MPa to Pa [23].
Molecular mass of the mixture is expressed by following equation:
M mix = M i · x i
where: M mix —molecular mass of mixture, kg/kmol; xi—molar fraction of constituent i.
V mix = M i · V i [ K 1 + ( k 2 k 1 ) · ( x N 2 0.0425 ) ] · x CH 4
where: V mix —molar volume of component i at the temperature of LNG, m3; V i —molar volume of constituent i at the temperature of LNG, m3; x i —molar fraction of constituent i, mol; k 1 ,   k 2 —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):
p s , corr = p s , CH 4 + x N 2 · 0.11   MPa · ( T 90   K ) x C 2 H 6 · 0.05   MPa · ( T 95   K )
where: 0.11   MPa and 0.05   MPa are determined according to [23].
Pseudo critical temperature:
T pc = i x i · T c , i
where: x i —molar fraction of constituent, mol; T c , i —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 m3/h, which can be reduced to 220 m3/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]:
M ˙ EM 1 = c EM 1 ( ω EM 10 ω EM 1 ) d EM 1 M EM 1
where c EM 1 , d EM 1 —parameters of the electric motor; ω EM 1 —angular velocity of the electric motor, rad/s.
Equation of rotation of the electric motor’s rotor in the SMR is equal to:
I EM 1 ω ˙ EM 1 = M EM 1 M pump 1 ( ω EM 1 , p pump 1 ) c 1 ω EM 1
where I EM 1 —mass inertia moment of the rotor of electric motor and pump’s rotor; c 1 —damping coefficient;   M pump 1 ( ω EM 1 , p pump 1 ) —load torque of the electric motor’s rotor.
M pump 1 ( ω EM 1 , p pump 1 ) = p pump 1 ( Q pump 1 , ω EM 1 ) Q L 1 , 1 η pump 1 ω EM 1
where Q pump 1 —total flow of SMR pump, m3/s;   η Pump 1 —coefficient of SMR pump, (COP); p pump 1 ( Q pump 1 , ω EM 1 ) —is total SMR pump characteristic;   Q L 1 , 1 —LNG flow in the first section of pipeline, m3/s.
p Pump 1 ( Q pump 1 , ω EM 1 ) = ( ω EM 1 ω EM 1 , 0 ) 2 · p pump 1 , 0 ( Q pump 1 , 0 = ω EM 10 ω EM 1 Q L 1 , 1 )
where p ( Q pump 1 , 0 )—total SMR characteristic at an angular velocity ω EM 1 , 0 .
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:
p ( Q pump 1 , 0 ) = 5546.7 Q pump 1 , 0 5 ump 1 , 0 Q pump 1 , 0 4 + 473.16 Q pump 1 , 0 3 ump 1 , Q pump 1 , 0 2 ump 1 , 0 Q pump 1 , 0 + 0.8757

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 ( G p , in ), which enters the passage between the pipe and the first stage of pump’s joint, is determined according to differential Equation (16):
G ˙ p , in = ( A p , in L p , in ) · [ p N p p , 1 Δ p p , in sign ( G p , in ) ]
where A p , in —section area between point N and V p 1 , m2; L p , in —length between point N and p p , 1 , m; p N —pressure of inlet pipe, Pa; p p , 1 —pressure of the passage between the pipe and the first stage of pump’s joint, Pa; Δ p p , in —pressure losses between point N and p p 1 , Pa; sign ( G p , in ) —sign function, which indicates LNG mass flow direction.
Pressure losses ( Δ p p , in   ) are calculated by Equation (17):
Δ p p , in = 0.5 · ξ p , in · G p , in 2 ρ l ( p N ) A p , in 2
Hydraulic losses coefficient ( ξ p , in ) is determined by Equation (18):
ξ p , in = ξ p , in , local + λ p , in ( Re , Δ p , in ) 8
To calculate hydraulic losses, the Δ p , in   could be determined according to these conditions:
{ λ p , i n = 64 R e   ,   w h e n   R e 2320   F ( λ p , i n ) = 1 λ p , i n 2 · l o g ( 2 , 51 R e · λ p , i n   + Δ p , i n 3.7 · D p , i n ) = 0 ,   w h e n   R e 2320 }
In that case, when LNG flow is turbulent, it could be solved by the Colebrook–White nonlinear algebraic equation F( λ p , i n ) = 0 [29,30].
The pressure of the passage between the pipe and the first stage of pump’s joint p p , 1 could be determined by Equation (20):
p ˙ p 1 = K ( p p , 1 ) V p , 1 · [ G p , in ρ l ( p N ) G p 12 ρ l ( p p , 1 ) ]
where K ( p p , 1 ) —bulk modulus of elasticity of LNG, Pa; G p 12 —LNG mass flow between p p , 1   and   p p , 2 , kg/s; ρ l ( p N ) —LNG density in inlet pipe, kg/m3; ρ l ( p p , 1 ) —LNG density in the passage between the pipe and the first stage of pump’s joint, kg/m3.
LNG mass flow between p p , 1   and   p p , 2 ( G p 12 ) is determined by differential Equation (21) below:
G ˙ p 12 = ( A p , 12 L p , 12 ) · [ p p 1 +   p pump ( ω EM , Q p , out ) p p , 2 Δ p p 12 sign ( G p 12 ) ]
where Δ p p 12 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:
p pump = N stage · p pump , stage ( Q pump ω EM )
Pump pressure of the one of the stages is determined by Equation (23):
p pump , stage ( Q pump ω EM ) = ( ω EM ω EM 0 ) 2 · p pump , stage , 0 ( Q pump 0 = ω EM 0 ω EM Q pout )
where Q pout = G out ρ l ( p p , 2 ) ; p pump , stage , 0 —one of the stages of pump’s characteristics, which is provided below:
p pump 2 , stage = 8.25511502 · 10 6 0.10 · Q pump 2 , 0 1.523478682 · 10 9 · Q pump 2 , 0 2
Pressure ( p p , 2 ) of the passage between the first stage of the pump’s joint and the outlet pipe is determined by:
p ˙ p 2 = K ( p p , 2 ) V p , 2 · [ G p 12 ρ ( p p , 1 ) G p , out ρ ( p p , 2 ) ]
where K ( p p , 2 ) —bulk modulus of elasticity of LNG, Pa; G p 12 —LNG mass flow between p p , 1   and   p p , 2 , kg/s; ρ l ( p p , 1 ) —LNG density of Vp,1, kg/m3; ρ l ( p p , 2 ) —LNG density of Vp,2, kg/m3, G p , out —mass flow of LNG from outlet of pump to the outlet pipe, kg/s.
The mass flow of LNG ( G p , out ) from outlet of pump to the outlet pipe is determined by Equation (26):
G ˙ p , out = ( A p , out L p , out ) · [ p p , 2 p p , 1 Δ p p , out sign ( G p , out ) ]
where A p , in —sectional area between p p 2 and point 1, m2; L p , in —length between   p p 2 and point 1, m; p p 1 —pressure of the passage between the pipe and the first stage of pump’s joint, Pa; p p 2 —pressure of the passage between the first stage of the pump’s joint and the outlet pipe, Pa; Δ p pout —pressure losses between p p 2 and point 1, Pa; sign ( G p , out ) —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]:
t [ A ( x ) ρ l ] + x [ A ( x ) ρ l υ ] = 0
where ρ l , υ −density and velocity of LNG, kg/m3 and m/s; A ( x )   −cross-sectional area of a pipeline, m2.
Equation of liquid flow impulse (momentum):
t [ A ( x ) ρ l υ ] + t [ A ( x )   ( p + ρ l υ 2 ) ] + Π ( x ) τ + A ( x ) ρ l a x = p l A x
where p—LNG pressure, Pa; Π ( x ) —the perimeter of cross-section of the pipeline, m; τ —tangential fluid stress on the inner surface of the pipeline, Pa;
τ = λ ( Re , Δ ) ρ υ | υ | 8  
where p—LNG pressure, Pa; Π ( x ) −the perimeter of a cross-section of the pipeline, m; a x —the acceleration along the x axis, m2/s.
a x = g · sin α
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:
[ A ] { u t } + [ B ( u ) ] { u x } = { f ( u ) }
where [ A ] , [ B ( u ) ] —matrices and { f ( u ) } —a vector which depends on t, x and elements of a vector { u } T = [ p , v ] Matrices are described:
[ A ] = [ 1 0 0 1 ] ;   [ B ( u ) ] = [ v c ρ l 2 1 / ρ l v ]
Vector is described:
{ f ( u ) } = { c 2 ρ l v S ( x ) A x τ Π ( x ) A ( x ) ρ l a x }
where c—speed of sound in the LNG, m/s.
By equating the determinant of the matrix [ B ( u ) ] [ A ] dx dt to zero, which allows determining the dx dt   derivative and determines characteristic direction. This equation has two various real roots: dx dt = λ i   , ( i = 1.2 )  
C + : dx dt = υ + c ;   C : dx dt = υ c
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:
c = K ( p , T ) ρ l 1 + K ( p , T ) · D in E · e
where K ( p , T ) —bulk modulus of elasticity of liquid, Pa; K ( p , T ) = ρ l dp d ρ l ; ρ l —density of the LNG; E—modulus of elasticity of a pipeline, Pa; D in —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 ( v ) and the pressure of the liquid ( p ) 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:
C + : Φ 1 = v D v D + 1 2 ( ρ D ρ L ) [ ( 1 ρ l c ) L + ( 1 ρ l c ) D ] Δ t 2 [ ( f 1 ρ l c ) L + ( f 1 ρ l c ) D ] Δ t 2 [ ( f 2 ) L + ( f 2 ) D ] = 0
C : Φ 2 = v D v R + 1 2 ( ρ D ρ R ) [ ( 1 ρ l c ) R + ( 1 ρ l c ) D ] Δ t 2 [ ( f 1 ρ l c ) R + ( f 1 ρ l c ) D ] Δ t 2 [ ( f 2 ) R + + ( f 2 ) D ] = 0
where f 1 ( p , υ ) = c S ( x ) ( ρ l υ S x ) ; f 2 ( υ , p ) = τ Π ( x )   ρ l S ( x ) .
The nonlinear algebraic Equations (36) and (37) are solved by the Newton and Raphson method [39]:
[ J ] i 1   { Δ Y } i = { Φ } i 1 ;   { Y } i = { Δ Y } i 1 + { Δ Y } i
[ J ] i 1   = [ Φ 1 , i 1 p Φ 1 , i 1 Δ v Φ 2 , i 1 p Φ 2 , i 1 v ]   ;   { Δ Y } T i = [ Δ p i ,   Δ v i ] ;   { Φ } i 1 T = [ Φ 1 , i 1 , Φ 2 , i 1 ]
where [ J ] i 1   — Jacobi matrix; { Δ Y } i — 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]:
Δ t | v + c | Δ x C r
where C r is Courant number.
At points L and R, variables p and v are defined using a system of two nonlinear algebraic equations with unknown pL, vL, and pR, vR:
Φ 3 = p L p C θ ( p A p C ) [ v L + c ( p L ) ] = 0
Φ 4 = v L v C θ ( v A v C ) [ v L + c ( p L ) ] = 0
And
Φ 5 = p R p C + θ ( p B p C ) [ v R c ( p R ) ] = 0
Φ 6 = v R v C + θ ( v B v C ) [ v R c ( p R ) ] = 0
where =   Δ t Δ x .
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 x = 0 , pressure is known variable p ( t ) = p D , then LNG flow velocity is determined by formula:
        Φ 2 ( v D ) = 0
  • Boundary condition: when the point is known x = 0 , LNG flow velocity is known variable v ( t ) = v D , then pressure is determined by formula:
    Φ 2 ( p D ) = 0
  • Boundary condition: when the point is known x = L , pressure is known variable p ( t ) = p D , then LNG flow velocity is determined by formula:
    Φ 1 ( v D ) = 0
  • Boundary condition: when the point is known x = L , NG flow velocity is known variable v ( t ) = v D , then pressure is determined by formula:
    Φ 1 ( p D ) = 0
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:
p ˙ V 3 = K ( p V 3 , T ) V 3 ( Q L 2 , N ) sign ( Q L 2 , N ) Q L 3 , 1 sign ( Q L 3 , 1 ) Q L 4 , 1 sign ( Q L 4 , 1 ) Q L 5 , 1 sign ( Q L 5 , 1 ) Q Losses sign ( p V 3 ) )
where Q L 2 , N , Q L 3 , 1 , Q L 4 , 1 , Q L 5 , 1 , Q Losses —LNG volume flow, m3/s; Q Losses ( p V 3 ) –losses of LNG volume flow, m3/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:
C + :   Φ 1 ( p V 3 , v L 2 , N ) = 0 C :   Φ 2 ( p V 3 , v L 3 , 1 ) = 0 C :   Φ 2 ( p V 3 , v L 4 , 1 ) = 0 C :   Φ 2 ( p V 3 , v L 5 , 1 ) = 0
For mathematical model construction, there are two important conditions:
(a) 
When pipes (L3, L4, L5) are closed until HP booster pump 1 starts working:
C + :   Φ 1 ( p L 3 ( x = L 3 ) , v L 3 end = 0 ) C + :   Φ 1 ( p L 4 ( x = L 4 ) , v L 4 end = 0 ) C + :   Φ 1 ( p L 5 ( x = L 5 ) , v L 5 end = 0 )
(b) 
When HP booster pump 1 starts working and pipes start opening:
C + :   Φ 1 ( p L 3 ( x = L 3 ) , v L 3 end = G p , in A p , in · ρ l   ) C + :   Φ 1 ( p L 4 ( x = L 4 ) , v L 4 end = G p , in A p , in · ρ l ) C + :   Φ 1 ( p L 5 ( x = L 5 ) , v L 5 end = G p , in A p , in · ρ l )
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:
C + :   Φ 1 ( p L 6 ( x = L 6 ) , v L 6 end ) = 0 C + :   Φ 1 ( p L 7 ( x = L 7 ) , v L 7 end ) = 0 C + :   Φ 1 ( p L 8 ( x = L 8 ) , v L 8 end ) = 0
The pressure at the end of the pipes L6, L7 and L8 is determined according to the expression below:
p Lk ( x = Lk ) = 0.5 · ξ Lk end · v Lk end 2 + p Lk , 0
where p Lk , 0 —constant pressure in the pipes.
Hydraulic losses coefficient:
ξ Lk = [ ( A Lk n channel , k · A channel , k ) 2 1 ]
where pipe index—k = 6, 7, 8 and A Lk —cross-section area of pipe k, m2; A channel , k —cross-section area of channel of k-th BOG recondenser, m2; n channel , k —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 m3/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 m3/h, compared to the real system according to the technical parameters, the value is 520 m3/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 ω EM 1 ( t ) , we obtain the pump curve p Pump 1 ( Q pump 1 , ω EM 1 ) at the corresponding motor angular velocity ω EM 1 ( t ) . 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.

Author Contributions

Conceptualization, V.S. and M.B.; methodology, V.S. and M.B.; software, V.S., M.B. and A.D.; validation, M.B. and V.S.; formal analysis, V.S., M.B. and A.D.; investigation, V.S. and M.B.; resources, V.S. and T.P.; data curation, V.S. and M.B.; writing—original draft preparation, V.S., J.U. and M.B.; writing—review and editing, T.P. and J.U.; visualization, V.S. and M.B.; supervision, V.S. and M.B.; project administration, V.S. and M.B.; funding acquisition, T.P. and V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Commission. The European Green Deal, Brussels 2019. Available online: https://ec.europa.eu/info/sites/default/files/european-green-deal-communication_en.pdf (accessed on 5 September 2022).
  2. DNV GL AS. Technology Outlook 2025. Norway 2016. Available online: https://brandcentral.dnv.com/fileroot/gallery/dnvgl/files/original/4f01c13e94af4438b971c5254b7d8664/4f01c13e94af4438b971c5254b7d8664_low.pdf (accessed on 5 September 2022).
  3. International Gas Union. World LNG Report 2022. IGU 2022. Available online: https://www.igu.org/resources/world-lng-report-2022/ (accessed on 5 September 2020).
  4. Statistical Review of World Energy. BP, London 2020. Available online: https://www.bp.com/content/dam/bp/business-sites/en/global/corporate/pdfs/energy-economics/statistical-review/bp-stats-review-2020-full-report.pdf (accessed on 5 September 2020).
  5. Wang, C.; Ju, Y.; Wang, T.; Zou, S. Transient performance study of high pressure fuel gas supply system for LNG fueled ships. Cryogenics 2022, 125, 103510. [Google Scholar] [CrossRef]
  6. Bosser, L.; Moreno, G. The TEN-T Review: A Welcome but Perfectible Proposal. Conference of Peripheral Maritime Regions, 2022. Available online: https://cpmr.org/wpdm-package/the-ten-t-review-a-welcome-but-perfectible-proposal/ (accessed on 5 September 2020).
  7. Ren, J.; Liang, H. Measuring the sustainability of marine fuels: A fuzzy group multi-criteria decision making approach. Transp. Res. Part D Transp. Environ. 2017, 54, 12–29. [Google Scholar] [CrossRef]
  8. International Gas Union. World LNG Report. IGU 2020. Available online: https://www.igu.org/resources/2020-world-lng-report/ (accessed on 5 September 2020).
  9. Semaskaite, V.; Bogdevicius, M.; Paulauskiene, T.; Uebe, J.; Filina-Dawidowicz, L. Improvement of Regasification Process Efficiency for Floating Storage Regasification Unit. J. Mar. Sci. Eng. 2022, 10, 897. [Google Scholar] [CrossRef]
  10. Khan, M.S.; Wood, A.D.; Qyyum, A.M.; Ansari, K.B.; Ali, W.; Wazwaz, A.; Dutta, A. Graphical approach for estimating and minimizing boil-off gas and compression energy consumption in LNG regasification terminals. J. Nat. Gas Sci. Eng. 2022, 101, 104539. [Google Scholar] [CrossRef]
  11. The American Bureau of Shipping. LNG Bunkering: Technical and Operational Advisory; The American Bureau of Shipping: Houston, TX, USA, 2015; pp. 8–65. Available online: https://superlng.adrioninterreg.eu/wp-content/uploads/2020/04/2_FILE.pdf (accessed on 7 September 2020).
  12. Andreev, A.M.; Baranov, A.Y.; Sokolova, E.V.; Malysheva, T.A.; Ivanov, L.V. Selection of technological solutions to reduce LNG loss in supply pipelines to LNG terminal. AIP Conf. Proc. 2021, 2412, 030023. [Google Scholar] [CrossRef]
  13. Koo, B. A novel implicit method of characteristics using pressure-referenced correction for transient flow in natural gas pipelines. J. Nat. Gas Sci. Eng. 2022, 104, 104665. [Google Scholar] [CrossRef]
  14. Li, J.; Hu, H.; Wang, H. Numerical investigation on flow pattern transformation and heat transfer characteristics of two-phase flow boiling in the shell side of LNG spiral wound heat exchanger. Int. J. Therm. Sci. 2020, 152, 106289. [Google Scholar] [CrossRef]
  15. Pan, J.; Wang, J.; Tang, L.; Bai, J.; Lia, R.; Lub, Y.; Wua, G. Numerical investigation on thermal-hydraulic performance of a printed circuit LNG vaporizer. Appl. Therm. Eng. 2020, 165, 114447. [Google Scholar] [CrossRef]
  16. Bai, J.; Pan, J.; He, X.; Wang, K.; Tang, L.; Yang, R. Numerical investigation on thermal hydraulic performance of supercritical LNG in sinusoidal wavy channel based printed circuit vaporizer. Appl. Therm. Eng. 2020, 175, 115379. [Google Scholar] [CrossRef]
  17. Zhao, Z.; Zhao, K.; Ji, D.; Jiang, P.; Shen, R. Numerical Investigation on the Flow and Heat Transfer Characteristics of Supercritical Liquefied Natural Gas in an Airfoil Fin Printed Circuit Heat Exchanger. Energies 2017, 10, 1828. [Google Scholar] [CrossRef]
  18. Karyakina, E.D.; Shammazov, I.A.; Shalygin, A.V. Main aspects of liquefied natural gas process line thermal and hydraulic calculations. IOP Conf. Ser. Earth Environ. Sci. 2021, 677, 052056. [Google Scholar] [CrossRef]
  19. Zhmakin, V.; Samoylov, A. Study of two-phase transportation mode of liquefied natural gas through a pipeline by the gravitational method. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1138, 012047. [Google Scholar] [CrossRef]
  20. Migliore, C.; Tuvilleja, C.; Vesovic, V. Weathering prediction model for stored liquefied natural gas (LNG). J. Nat. Gas Sci. Eng. 2015, 26, 570–580. [Google Scholar] [CrossRef] [Green Version]
  21. Huerta, F.; Vesovic, V. A realistic vapour phase heat transfer model for the weathering of LNG stored in large tanks. Energy 2019, 174, 280–291. [Google Scholar] [CrossRef]
  22. Krikkis, R.N. A thermodynamic and heat transfer model for LNG ageing during ship transportation. Towards an efficient boil-off gas management. Cryogenics 2018, 92, 76–83. [Google Scholar] [CrossRef]
  23. Tietz, C.; Richter, M.; Kleinrahm, R.; Span, R. Enhancement of the revised Klosek and McKinley method for density calculations of liquefied natural gas (LNG) over the temperature range from (100 to 135) K at pressures up to 10 MPa. Fuel Process. Technol. 2017, 165, 19–26. [Google Scholar] [CrossRef]
  24. GIIGNL. Custody Transfer Handbook, 5th ed.; International Group of Liquefied Natural Gas Importers: Neuilly sur Seine, France, 2017; pp. 1–179. Available online: https://giignl.org/system/files/giignl_cthb_5.0.web_.pdf (accessed on 1 October 2020).
  25. Shinko Ind. Ltd. Technical Sheet of Submerged Motor Pump, 2017. Available online: https://www.shinkohir.co.jp/pdf/catalog/SM.pdf (accessed on 7 September 2020).
  26. Wang, C.; Zhang, Y.; Hou, H.; Zhang, J.; Xu, C. Entropy production diagnostic analysis of energy consumption for cavitation flow in a two-stage LNG cryogenic submerged pump. Int. J. Heat Mass Transf. 2019, 129, 342–356. [Google Scholar] [CrossRef]
  27. Mokhatab, S.; Poe, W.; Mak, J. Handbook of Natural Gas Transmission and Processing, 3rd ed.; Elsevier: Walthan, MA, USA, 2015; pp. 10–70. [Google Scholar]
  28. Karpenko, M. Investigation of Energy Efficiency of Mobile Machinery Hydraulic Drives. Ph.D. Thesis, Vilnius Gediminas Technical University, Vilnius, Lithuania, 2021. Available online: http://dspace.vgtu.lt/bitstream/1/4269/1/M_Karpenko%20disertacija.pdf (accessed on 5 September 2020).
  29. Koo, B. Comparison of finite-volume method and method of characteristics for simulating transient flow in natural-gas pipeline. J. Nat. Gas Sci. Eng. 2022, 98, 104374. [Google Scholar] [CrossRef]
  30. Praks, P.; Brkić, D. Review of new flow friction equations: Constructing Colebrook’s explicit correlations accurately. Rev. Int. Métodos Numér. Cálc. Diseño Ing. 2020, 36, 41. Available online: https://www.scipedia.com/public/Praks_Brkic_2020a (accessed on 7 November 2022). [CrossRef]
  31. Semaškaitė, V.; Bogdevičius, M. Liquefied natural gas regasification technologies. In Transbaltica XII: Transportation Science and Technology. Transbaltica 2019. Lecture Notes in Intelligent Transportation and Infrastructure; Prentkovskis, O., Yatskiv (Jackiva), I., Skačkauskas, P., Junevičius, R., Maruschak, P., Eds.; Springer: Cham, Germany, 2022. [Google Scholar] [CrossRef]
  32. Bogdevicius, M.; Suslavicius, V. Hydrodynamic Processes of the Impulse Fire Extinguishing. Transp. Telecommun. 2006, 7, 342–349. [Google Scholar]
  33. Bajoraityte, I.; Bogdevicius, M. Dynamic processes in the oil pipelines. Part ii. The leak identification by numerical methods. Transport 2003, 18, 49–55. [Google Scholar] [CrossRef] [Green Version]
  34. Bajoraityte, I.; Bogdevicius, M. Dynamic processes in the main oil pipelines. Part i. Causes of leakage occurrence in oil pipelines and leak identification methods. Transport 2002, 17, 234–240. [Google Scholar] [CrossRef] [Green Version]
  35. Bogdevicius, M.; Janutiene, J.; Jonikas, K.; Guseinoviene, E.; Draksas, M. Mathematical modeling of oil transportation by pipelines using anti-turbulent additives. J. Vibroengineering 2013, 15, 419–427. [Google Scholar]
  36. Karpenko, M.; Bogdevicius, M. Investigation of hydrodynamic processes in the system—“axial piston pumps-pipeline-fittings”. In VII International Symposium of Young Researchers Transport Problems 2018; Silesian University of Technology, Faculty of Transport: Kartowice, Poland, 2018; pp. 832–843. [Google Scholar]
  37. Karpenko, M.; Bogdevicius, M. Investigation into the Hydrodynamic Processes of Fitting Connections for Determining Pressure Losses of Transport Hydraulic Drive. Transport 2020, 35, 108–120. [Google Scholar] [CrossRef] [Green Version]
  38. Jiang, Y.; Ren, Z.; Yang, X.; Li, Q.; Xu, Y. A steady-state energy flow analysis method for integrated natural gas and power systems based on topology decoupling. Appl. Energy 2022, 305, 118007. [Google Scholar] [CrossRef]
  39. Garrett, S.J. Chapter 13-introductory numerical methods. In Introduction to Actuarial and Finanacial Mathematical Methods; Garrett, S.J., Ed.; Elsevier: Amsterdam, The Netherlands, 2015; pp. 411–463. [Google Scholar] [CrossRef]
  40. Gnedin, N.Y.; Semenov, V.A.; Kravtsov, A.V. Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes. J. Comput. Phys. 2018, 359, 93–105. [Google Scholar] [CrossRef] [Green Version]
  41. Fedorko, G.; Molnar, V.; Dovicab, M.; Tothb, T.; Fabianova, J. Failure analysis of irreversible changes in the construction of the damaged rubber hoses. Eng. Fail. Anal. 2015, 58, 31–43. [Google Scholar] [CrossRef]
Figure 1. Simplified scheme of LNG transportation system in the FSRU. Note: L1, L2, L3, L4, L5, L6, L7, L8 is presented to illustrate the different section of pipes in LNG transportation system; EM1, EM2, EM3, EM4—electric motor indexes of pumps; pv1—vapor pressure in the LNG tank; pv2—vapor pressure in the Suction Drum; TI1—LNG temperature in the LNG tank; TI2—LNG temperature in the Suction Drum.
Figure 1. Simplified scheme of LNG transportation system in the FSRU. Note: L1, L2, L3, L4, L5, L6, L7, L8 is presented to illustrate the different section of pipes in LNG transportation system; EM1, EM2, EM3, EM4—electric motor indexes of pumps; pv1—vapor pressure in the LNG tank; pv2—vapor pressure in the Suction Drum; TI1—LNG temperature in the LNG tank; TI2—LNG temperature in the Suction Drum.
Jmse 10 01960 g001
Figure 2. Principal scheme of the mathematical model.
Figure 2. Principal scheme of the mathematical model.
Jmse 10 01960 g002
Figure 3. LNG level observation in the LNG tank.
Figure 3. LNG level observation in the LNG tank.
Jmse 10 01960 g003
Figure 4. LNG level observation in the Suction Drum.
Figure 4. LNG level observation in the Suction Drum.
Jmse 10 01960 g004
Figure 5. High-pressure booster pump mathematical model. Note: N, N-1—numbers of inlet pipe sections; 1, 2—numbers of outlet pipe sections; Vp,1—the passage between the pipe and the first stage of pump’s joint; Vp,2—the passage between the pipe and the second stage of pump’s joint; Gp—mass flow of LNG.
Figure 5. High-pressure booster pump mathematical model. Note: N, N-1—numbers of inlet pipe sections; 1, 2—numbers of outlet pipe sections; Vp,1—the passage between the pipe and the first stage of pump’s joint; Vp,2—the passage between the pipe and the second stage of pump’s joint; Gp—mass flow of LNG.
Jmse 10 01960 g005
Figure 6. Method of characteristics (a) Schematic of characteristic line for MOC; (b) Calculation schemes of the first and the last points.
Figure 6. Method of characteristics (a) Schematic of characteristic line for MOC; (b) Calculation schemes of the first and the last points.
Jmse 10 01960 g006
Figure 7. Pipeline Unit V3 connection with pipes (L2, L3, L4, L5). Note: see in Figure 1.
Figure 7. Pipeline Unit V3 connection with pipes (L2, L3, L4, L5). Note: see in Figure 1.
Jmse 10 01960 g007
Figure 8. Angular velocity of electric motors: EM1, EM2, EM3, EM4.
Figure 8. Angular velocity of electric motors: EM1, EM2, EM3, EM4.
Jmse 10 01960 g008
Figure 9. Pressure and velocity in the 1st pipeline (L1) (Table 2) of the LNG transportation system. Note: (a) pressure observation when the system starts working; (b) pressure observation during simulated solution time in the first pipeline in the LNG transportation system; (c) velocity observation when the system starts working; (d) velocity observation during simulated solution time in the first pipeline in the LNG transportation system.
Figure 9. Pressure and velocity in the 1st pipeline (L1) (Table 2) of the LNG transportation system. Note: (a) pressure observation when the system starts working; (b) pressure observation during simulated solution time in the first pipeline in the LNG transportation system; (c) velocity observation when the system starts working; (d) velocity observation during simulated solution time in the first pipeline in the LNG transportation system.
Jmse 10 01960 g009aJmse 10 01960 g009b
Figure 10. LNG level observation from the beginning to the end of solution time: (a) in the LNG tank; (b) in the Suction Drum.
Figure 10. LNG level observation from the beginning to the end of solution time: (a) in the LNG tank; (b) in the Suction Drum.
Jmse 10 01960 g010
Figure 11. Pressure in the 2nd pipeline (L2) (Table 2) of the LNG transportation system: (a) the beginning of the simulation time; (b) the end of the simulation time.
Figure 11. Pressure in the 2nd pipeline (L2) (Table 2) of the LNG transportation system: (a) the beginning of the simulation time; (b) the end of the simulation time.
Jmse 10 01960 g011
Figure 12. Pressure in the Pipeline Unit of the LNG transportation system: (a) the beginning of the simulation time; (b) the end of simulation.
Figure 12. Pressure in the Pipeline Unit of the LNG transportation system: (a) the beginning of the simulation time; (b) the end of simulation.
Jmse 10 01960 g012
Figure 13. Mass flows of LNG of HP Booster Pumps: (a) the beginning of the simulation time; (b) the end of the simulation.
Figure 13. Mass flows of LNG of HP Booster Pumps: (a) the beginning of the simulation time; (b) the end of the simulation.
Jmse 10 01960 g013
Figure 14. Pressure in the 6th, 7th, and 8th pipes (Table 2) of the LNG transportation system: Note: (a) Pressure observation when the system starts working; (b) Pressure observation at the end of the solution time; (c) Velocity observation when the system starts working; (d) Velocity observation at the end of the solution time.
Figure 14. Pressure in the 6th, 7th, and 8th pipes (Table 2) of the LNG transportation system: Note: (a) Pressure observation when the system starts working; (b) Pressure observation at the end of the solution time; (c) Velocity observation when the system starts working; (d) Velocity observation at the end of the solution time.
Jmse 10 01960 g014
Figure 15. Pressure in the 6th pipe (L6) (Table 2) of the LNG transportation system: Note: (a) Pressure observation when the system starts working; (b) Pressure observation at the end of the solution time; (c) Velocity observation when the system starts working; (d) Velocity observation at the end of the solution time.
Figure 15. Pressure in the 6th pipe (L6) (Table 2) of the LNG transportation system: Note: (a) Pressure observation when the system starts working; (b) Pressure observation at the end of the solution time; (c) Velocity observation when the system starts working; (d) Velocity observation at the end of the solution time.
Jmse 10 01960 g015
Figure 16. Hydraulic energy losses. Note: (a) Hydraulic energy losses in the 1st pipeline (L1); (b) Hydraulic energy losses in all pipelines and pipes.
Figure 16. Hydraulic energy losses. Note: (a) Hydraulic energy losses in the 1st pipeline (L1); (b) Hydraulic energy losses in all pipelines and pipes.
Jmse 10 01960 g016
Figure 17. Curves at various angular velocities for the main pump (SMR pump Model SMR 200; see Figure 8) as a component of the presented model.
Figure 17. Curves at various angular velocities for the main pump (SMR pump Model SMR 200; see Figure 8) as a component of the presented model.
Jmse 10 01960 g017
Table 1. Symbols and Notations.
Table 1. Symbols and Notations.
SymbolMeaning
L1, L2, L3, L4, L5, L6, L7, L8the different section of pipes in the LNG transportation system
EM1, EM2, EM3, EM4electric motor indexes of pumps
pv1vapor pressure of BOG in the LNG tank, Pa
pv2vapor pressure of BOG in the Suction Drum, Pa
TI1LNG temperature in the LNG tank, K
TI2LNG temperature in the Suction Drum, K
H tank1 ( t ) liquid level in the tank by the time changing in the LNG tank, m
H tank1,0 height between the pump and the bottom of the LNG tank, m
H tank2 ( t ) liquid level in the tank by the time changing in the Suction Drum, m
ggravitational acceleration, m/s
A 1 and A 2 cross-sectional areas of the LNG liquid surface in the LNG tank and in the Suction Drum, m2
Q 1 , in   and Q L 1 , out LNG volume flow of the pipeline (L1) inlet and outlet, m3/s
Q 2 , in   LNG volume flow of the pipeline (L2) inlet, m3/s
ρ l LNG density, kg/m3
M mix molecular mass of mixture, kg/kmol
x i molar fraction of constituent i
V mix molar volume of component i at the temperature of LNG, m3
V i molar volume of constituent i at the temperature of LNG, m3
x i molar fraction of constituent i, mol
k 1 ,   k 2 volume correction factors, mol/l
T c , i critical temperature of constituent i, K
c EM 1 , d EM 1 parameters of the electric motor
ω EM 1 angular velocity of the electric motor, rad/s
I EM 1 mass inertia moment of the rotor of electric motor and pump’s rotor
c 1 damping coefficient
M pump 1 ( ω EM 1 , p pump 1 ) load torque of the electric motor’s rotor
Q pump 1 total flow of Submerged Motor Pump (SMR) pump, m3/s
η pump 1 coefficient of SMR pump
p pump 1 ( Q pump 1 , ω EM 1 ) total SMR pump characteristic
Q L 1 , 1 LNG flow in the first section of pipeline, m3/s
p ( Q pump 1 , 0 ) total   SMR   characteristic   at   an   angular   velocity   ω EM 1 , 0
A p , in sec tion   area   between   point   N   and   V p 1 , m2
L p , in length   between   point   N   and   p p , 1 , m
p N pressure of inlet pipe, Pa
p p , 1 pressure of the passage between the pipe and the first stage of pump’s joint, Pa
Δ p p , in pressure   losses   between   point   N   and   p p 1 , Pa
sign ( G p , in )   and   sign ( G p , out ) sign function, which indicates LNG mass flow direction
K ( p p , 1 ) bulk modulus of elasticity of LNG, Pa
G p 12 LNG   mass   flow   between   p p , 1   and   p p , 2 , kg/s
ρ l ( p N ) LNG density in inlet pipe, kg/m3
ρ l ( p p , 1 ) LNG density in the passage between the pipe and the first stage of pump’s joint, kg/m3
K ( p p , 2 ) bulk modulus of elasticity of LNG, Pa
G p 12 LNG   mass   flow   between   p p , 1   and   p p , 2 , kg/s
ρ l ( p p , 1 )   and   ρ l ( p p , 2 ) LNG density of Vp,1 and Vp,2, kg/m3
G p , out mass flow of LNG from outlet of pump to the outlet pipe, kg/s
A p , in sec tional   area   between   p p 2 and point 1, m2
L p , in length   between   p p 2 and point 1, m
p p 1 pressure of the passage between the pipe and the first stage of pump’s joint, Pa
p p 2 pressure of the passage between the first stage of pump’s joint and the outlet pipe, Pa
Δ p pout pressure   losses   between   p p 2 and point 1, Pa
ρ l density of LNG, kg/m3
υ velocity of LNG, m/s
A ( x ) cross-sectional area of a pipeline, m2
pLNG pressure, Pa
Π ( x ) the perimeter of cross-section of the pipeline, m
τ tangential   fluid   stress   on   the   inner   surface   of   the   pipeline ,   Pa  
a x the acceleration along the x axis, m2/s
α angle between horizontal axis and center line of pipe
cspeed of sound in the LNG, m/s
K ( p , T ) bulk modulus of elasticity of liquid, Pa
Emodulus of elasticity of a pipeline, Pa
D in internal diameter of a pipeline, m
ethickness of a wall of a pipeline, m
[ J ] i 1   Jacobi matrix
{ Δ Y } i increment vector of unknown variables
inumber of iterations
C r Courant number
Q L 2 , N , Q L 3 , 1 , Q L 4 , 1 , Q L 5 , 1 , Q Losses LNG volume flow, m3/s
Q Losses ( p V 3 ) losses of LNG volume flow, m3/s
p Lk , 0 constant pressure in the pipes
A Lk cross-section area of pipe k, m2
A channel , k cross-section area of channel of k-th BOG recondenser, m2
n channel , k number of channels of k-th BOG recondenser
Table 2. Specifications of FSRU equipment.
Table 2. Specifications of FSRU equipment.
EquipmentSpecification
LNG TankMaximum liquid level: 27 m.
Cross-section area: 78 m2.
Submerged Motor Cargo PumpMaximum capacity: 550 m3/h.
Electric motor power: 132 kW.
Mass inertial moment of electric motor rotor: 1.5 kg/m2.
Mass inertial moment of pump wheel: 13 kg/m2.
Suction DrumMaximum liquid level: 5 m.
Cross-section area: 27 m2.
High-Pressure Booster PumpMaximum capacity: 510 m3/h.
Electric motor power: 1000 kW.
Mass inertial moment of pump rotor: 13.
Mass inertial moment of electric motor rotor: 21 kg/m2.
Mass inertial moment of pump wheel: 13 kg/m2.
BOG Recondensers (Heat Exchangers)Number of channels (from pipe L6, 600,000.
Number of channels (from pipe L7): 650,000.
Number of channels (from pipe L8): 700,000.
Table 3. Specification of FRSU LNG transportation pipelines system.
Table 3. Specification of FRSU LNG transportation pipelines system.
Pipe NameDiameter, mLength, mTilt Angle
Pipeline from the LNG tank to the Suction Drum (L1)
L1,10.4040
L1,20.401090
L1,30.4063180
L1,40.5010180
L1,50.6011180
L1,60.50104180
L1,70.2020270
L1,80.2020
Pipeline from the Suction Drum to the Pipeline Unit (L2)
L2,10.401270
L2,20.4010
L2,30.40290
Pipe from the Pipeline Unit to the HP Booster Pump 1L30.2010
Pipe from the Pipeline Unit to the HP Booster Pump 2L40.2010
Pipe from the Pipeline Unit to the HP Booster Pump 3L50.2010
Pipe from the HP Booster Pump 1 to the BOG recondenser 1L60.2020
Pipe from the HP Booster Pump 2 to the BOG recondenser 2L70.2020
Pipe from the HP Booster Pump 3 to the BOG recondenser 3L80.2020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bogdevicius, M.; Semaskaite, V.; Paulauskiene, T.; Uebe, J.; Danilevicius, A. Modelling and Simulation Hydrodynamics Processes in Liquefied Natural Gas Transportation Systems. J. Mar. Sci. Eng. 2022, 10, 1960. https://doi.org/10.3390/jmse10121960

AMA Style

Bogdevicius M, Semaskaite V, Paulauskiene T, Uebe J, Danilevicius A. Modelling and Simulation Hydrodynamics Processes in Liquefied Natural Gas Transportation Systems. Journal of Marine Science and Engineering. 2022; 10(12):1960. https://doi.org/10.3390/jmse10121960

Chicago/Turabian Style

Bogdevicius, Marijonas, Vigaile Semaskaite, Tatjana Paulauskiene, Jochen Uebe, and Algimantas Danilevicius. 2022. "Modelling and Simulation Hydrodynamics Processes in Liquefied Natural Gas Transportation Systems" Journal of Marine Science and Engineering 10, no. 12: 1960. https://doi.org/10.3390/jmse10121960

APA Style

Bogdevicius, M., Semaskaite, V., Paulauskiene, T., Uebe, J., & Danilevicius, A. (2022). Modelling and Simulation Hydrodynamics Processes in Liquefied Natural Gas Transportation Systems. Journal of Marine Science and Engineering, 10(12), 1960. https://doi.org/10.3390/jmse10121960

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

Article Metrics

Back to TopTop