^{1}

^{1}

^{1}

Floating liquefied natural gas (FLNG) facility using partially filled tanks for control of pitch motion response to wave-exciting forces is investigated in this paper. The governing equations for sloshing analysis of rectangular tanks under pitch motion excitation are first established, then the spatial (boundary- value) partial derivatives are approximated by finite differences. The uncoupled pitch equation of FLNG is derived by assuming that pitch is uncoupled from other modes of vibration. By using state-space model to represent fluid-memory effect, the pitch equation can be transformed to first- order ordinary differential equations which can be solved with sloshing equations simultaneously with the given initial conditions. By using the proposed coupling model for FLNG facility and the liquefied natural gas (LNG) tanks, the performance of partially filled tanks for suppressing pitching motions of FLNG facility is numerically assessed. The parametric studies on the example FLNG show that there is a beneficial filling level by which the pitch motion of FLNG can be considerably reduced.

Floating liquefied natural gas (FLNG) facility represents a new technology with a promising future in the offshore oil industry, which combines the production, storage, transportation with a liquefied natural gas (LNG) liquefaction unit. The LNG liquefaction equipment used in FLNG facility will suffer considerable wave-induced motions which appear within a wide range of sea states and lasts for almost the whole design life of the FLNG facility. As a result, the frequent occurrence of the wave-induced motions of the FLNG facility of relatively large amplitude may cause low efficiency of liquefaction process or malfunction of the equipment. Hence, it is of practical interest to reduce excessive wave-induced motions of the FLNG facility.

The FLNG facility which is moored directly above the natural gas field possesses a number of large LNG tanks. A single point mooring (SPM) mechanism allows FLNG to rotate or weathervane freely with respect to a fixed mooring point and pitch motions play an important role in FLNG hydrodynamic responses. The wave-induced motions of the FLNG facility may cause severe sloshing in a partially filled tank and, in turn, the sloshing in the tank impacts the motions of the FLNG facility. The coupling effect between the sloshing in the tank and the FLNG facility motion is of a great concern to the FLNG operation in the production site and offloading operation of LNG carriers close to FLNG facility. Malenica et al. presented an efficient method for dynamic coupling between sloshing and seakeeping in frequency domain [

A tuned liquid damper (TLD), which usually is a rigid tank partially filled with water in it, is a kind of passive control device relying upon liquid sloshing to suppress structural vibration. TLDs have been investigated and used as economical and effective dynamic vibration absorbers to mitigate the dynamic response of structures for decades [

In this paper, a coupling model for FLNG facility and the LNG tanks is proposed. By using partially filled tanks as TLDs, the suppression of pitch motion of FLNG facility is investigated. The convolution term representing fluid-memory effect associated with the dynamics of the radiation forces is represented by the corresponding state-space model. Using finite differences to approximate the spatial partial derivatives in the continuity and Navier-Stokes equations, the Partial Differential Equations are transformed to the associated ordinary differential equations. Based on the derived formulae, the performance of partially filled tanks for suppressing pitching motion of FLNG facility is numerically assessed. It is found that if the filling level, i.e. the liquid depth, is selected suitably, the pitching motion of FLNG can be considerably reduced.

The FLNG facility considered in this study is equipped with n rectangular shaped LNG tanks. Each tank is partially filled with an incompressible and viscous liquid. The ith liquid tank is characterized by the liquid density ρ , the length 2 a i , width b i and liquid height h i . Some right-handed Cartesian coordinate systems are introduced to facilitate the formulation of the coupling problem. XYZ-system with origin at the center of gravity (COG) of the ship and Z axis pointing vertically upwards in quiescent water is used as a ship-fixed coordinate system for ship motion analysis. x i y i z i -system is the ith tank-fixed coordinate system with origin at the geometrical center of the free surface of the tank.

The continuity and Navier-Stokes equations in the ith tank-fixed coordinate system are [

∂ η i ∂ t + α i 2 u i ( η i ) ∂ η i ∂ x i + β i k i ∂ u i ( η i ) ∂ x i = 0 (1)

∂ u i ( η i ) ∂ t + α i 2 u i ( η i ) ∂ u i ( η i ) ∂ x i + [ θ ¨ e z i sin θ − θ ¨ x i + ( g − θ ¨ e x i ) cos θ ( 1 + β i k i ∂ 2 η i ∂ x i 2 ) ] ∂ η i ∂ x i = ( g − θ ¨ e x i ) sin θ − θ ¨ e z i cos θ − λ i u i ( η i ) (2)

where the subscript i denotes the ith tank, u i ( η i ) is the fluid velocity along the x i direction at the free surface, η i is the free-surface displacement, e x i and e z i are the x i - and z i -coordinates of the origin of the x i y i z i -system in the XYZ-system, θ is the pitch angle of the tank, α i = s e c h [ k i ( η i + h i ) ] , β i = tanh [ k i ( η i + h i ) ] , in which k i is the wave number; λ i = − ω l i ν i ( 1 + 2 h i / b i + S ) / ( η i + h ) , in which ν i is the kinematic viscosity coefficient and ω l i is the characteristic frequency of liquid motion, which may be taken as the natural frequency of the liquid tank, S is taken as unity.

The boundary conditions at the wall of the tank can be expressed as

u i ( − a , t ) = 0 , u i ( a , t ) = 0 , ( i = 1 , 2 , ⋯ , n ) (3)

The initial conditions for the liquid at rest initially are

η i ( x , 0 ) = 0 , u i ( x , 0 ) = 0 , ( i = 1 , 2 , ⋯ , n ) (4)

The pressure is expressed approximately as [

p i = ρ a z i ( η i − z i ) , ( i = 1 , 2 , ⋯ , n ) (5)

where

a z i = − ( g − θ ¨ e x i ) cos θ + θ ¨ x i − θ ¨ e z i sin θ (6)

The moment about the COG of the ship due to liquid sloshing is calculated as

#Math_30# (7)

The equations of motion for coupled sloshing and ship motion can be expressed as follows:

∑ k = 1 6 { [ M j k + A j k ( ∞ ) ] ξ ¨ k ( t ) + ∫ 0 t K j k ( t − τ ) ξ ˙ k ( τ ) d τ + C j k ξ k ( t ) } = F j e x t ( t ) + F j t a n k ( t ) ( j = 1 , ⋯ , 6 ) (8)

where F j e x t are the linear external wave force and moment components, F j t a n k are the force and moment components associated with sloshing; ξ k ( t ) are the displacements of the ship motion; M j k , A j k ( ∞ ) and C j k are, respectively, the components of the generalized mass, infinite-frequency added mass, and restoring matrices of the ship, K j k ( t ) are the retardation functions.

For simplicity, we assume that pitch is uncoupled from other modes of vibration although there are coupling effects between the vertical and longitudinal motions (heave, pitch, surge) of the ship. The uncoupled pitch equation can now be expressed as

#Math_39# (9)

where F 5 e x t denotes the linear external wave excitation pitch moment, M i t a n k denotes the moment associated with sloshing in the ith tank, I 55 is the ship's moment of inertia in pitch, A 55 ( ∞ ) is the infinite-frequency added mass in pitch, C 55 is the restoring coefficient in pitch, K 55 ( t ) is the retardation function in pitch. The calculation of retardation function for pitch can be given by noting the relationships [

#Math_46# (10)

A 55 ( ∞ ) = A 55 ( ω ) + 1 ω ∫ 0 ∞ K 55 ( t ) sin ( ω t ) d t (11)

in which A 55 ( ω ) and B 55 ( ω ) are added mass and damping coefficients in pitch, respectively.

The convolution integral in Equation (9) is known as fluid-memory effect and denoted by μ ( t ) , i.e.

μ ( t ) = ∫ 0 t K 55 ( t − τ ) θ ˙ ( τ ) d τ . (12)

Applying a Laplace transform to both sides of Equation (12), one obtains

μ ^ ( s ) = K ^ 55 ( s ) θ ˙ ^ ( s ) (13)

where K ^ 55 ( s ) , θ ˙ ^ ( s ) and μ ^ ( s ) are the Laplace transforms of K 55 ( t ) , θ ˙ ( t ) and μ ( t ) , respectively.

Ogilvie showed that the following frequency-domain representation holds for the retardation function [

K ^ 55 ( j ω ) = B 55 ( ω ) + j ω [ A 55 ( ω ) − A 55 ( ∞ ) ] . (14)

Hence, K ^ 55 ( s ) can be approximated by relative degree one transfer function [

K ^ 55 ( s ) = p r s r + p r − 1 s r − 1 + ⋯ + p 0 s n + q n − 1 s n − 1 + ⋯ + q 0 , r = n − 1 , n ≥ 2 (15)

As a result, μ ( t ) can be represented by the corresponding state-space model of the following form:

χ ˙ ( t ) = A μ χ ( t ) + B μ θ ˙ ( t ) (16a)

μ ( t ) = C μ χ ( t ) (16b)

where χ ( t ) is a state vector, A μ , B μ and C μ are constant matrices and given by

A μ = [ 0 1 0 ⋯ 0 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ 1 − q 0 − q 1 − q 2 ⋯ − q n − 1 ] , B μ = [ 0 0 0 ⋮ 1 ] ,#Math_71# (17)

Combining Equations (7), (9) and (16) yields

Γ ˙ ( t ) = f ( Γ , t ) (18)

where

Γ ( t ) = [ θ ψ χ ] , f ( Γ , t ) = [ ψ − g 1 ( θ ) − g 2 ( θ , t ) A μ χ + B μ ψ ] (19)

in which

g 1 ( θ ) = C μ χ ( t ) I 55 + A 55 ( ∞ ) + e 1 ( θ ) , #Math_76# (20)

e 1 ( θ ) = − 1 6 ρ ∑ i = 1 n b i σ i − ρ ∑ i = 1 n ∫ − a a b i ( η i + h i ) ( e x i + x i ) ς i d x i (21)

e 2 ( θ ) = 1 6 ρ g cos θ ∑ i = 1 n b i σ i + ρ g cos θ ∑ i = 1 n ∫ − a a b i ( η i + h i ) ( e x i + x i ) d x i (22)

σ i = ( η a i + h i ) 2 ( η a i + 3 e z i − 2 h i ) − ( η − a i + h i ) 2 ( η − a i + 3 e z i − 2 h i ) (23)

ς i = ( e x i cos θ + x i − e z i sin θ ) (24)

The spatial (boundary-value) partial derivatives in Equations (1) and (2) are approximated by finite differences. The difference equations with the given initial conditions are then solved simultaneously with Equations (18) using the Runge-Kutta method.

For application, a SPM FLNG as shown in

The irregular waves are generated using the JONSWAP spectra defined as

Parameters | Values | Parameters | Values | |
---|---|---|---|---|

Length between perpendiculars (Lpp) (m) | 200 | Volume displacement (m^{3}) | 97961 | |

Breadth (m) | 44 | Density of water (kg/m^{3}) | 1025 | |

Draught (m) | 12 | Moment of inertia, I_{55} (kg×m^{2}) | 2.9280 ´ 10^{11} | |

Mass (kg) | 1.004´10^{8} | Restoring coefficient, C_{55} (N×m) | 2.4766 ´ 10^{11} | |

Parameters | Values | Parameters | Values | |
---|---|---|---|---|

Length (m) | 56.5 | LNG density (kg/m^{3}) | 474 | |

Width (m) | 42.6 | Kinematic viscosity coefficient (m^{2}/s) | 0.29´10^{-6} | |

#Math_82# (25)

where γ = 3.3 , σ = { 0.07 ω ≤ 2 π / T p 0.09 ω > 2 π / T p , H s is the significant wave height,

and T p is the peak period.

F 5 w a v e ( ω k ) = | F 5 w a v e ( ω k ) | e j ∠ F 5 w a v e ( ω k ) (26)

where | F 5 w a v e ( ω k ) | and ∠ F 5 w a v e ( ω k ) are the pitch moment RAO amplitude and phase for frequency ω k . The amplitude and phase of pitch moment RAO are shown in

The time history of wave excitation pitch moment can be obtained from the corresponding RAO and JONSWAP spectrum using a series of trigonometric functions

#Math_93# (27)

where Δ ω is frequency interval, ε k is the corresponding random phase angle which uniformly distributed over the interval ( 0 , 2 π ) . Using equation (27), a typical simulated time history of pitch moment is calculated and shown in

estimated and plotted in

To investigate the performance of partially filled tanks for suppressing pitching motion of the FLNG, the pitch responses of the FLNG are computed against different liquid heights. Then, by comparing the results among all the pitch responses, the beneficial liquid height can be found for achieving the maximum or beneficial pitch response reduction of the FLNG with reasonable sloshing force. This study selects h of 8.09 m as a beneficial value for the hydrodynamic response analysis using the derived formulae.

Liquid sloshing inside the tanks of FLNG has been used for suppressing pitching motion of the FLNG. By using state-space model to represent the fluid-memory effect, a mathematical formulation is derived for coupled sloshing and ship motion analysis. The wave excitation pitch moment is generated using a series of trigonometric functions based on RAO amplitude and phase of the pitch moment and wave spectra. With the derived formulae, extensive parametric studies have been performed on the coupled system to seek beneficial liquid heights for maximum pitch response reductions with reasonable sloshing forces. It was

found that if liquid heights are selected appropriately, the pitch response of the FLNG could be considerably reduced.

The authors wish to acknowledge the financial support from the Project supported by No. 11572072 and National High Technology Research and Development Program of China (863 Program), No. 2014AA09A224.

Zhang, W.S., Cheng, J. and Duan, Y.Q. (2017) Liquefied Natural Gas Tanks for Suppressing Pitching Motions of FLNG Facility. Journal of Applied Mathematics and Physics, 5, 1638-1647. https://doi.org/10.4236/jamp.2017.59136