Two-Dimensional Numerical Modelling of a Moored Floating Body under Sloping Seabed Conditions

A coupled floating body-mooring line model is developed by combining a boundary element model for a two-dimensional floating body and a catenary mooring line model. The boundary element model is formulated in the time domain by a continuous Rankine source, and a reflection potential is introduced to account for the wave reflection due to sloping seabed. This newly developed model is validated by comparisons against available data. Then, dynamic response analyses are performed for the moored body in various seabed conditions. Compared with a flat seabed, a sloping seabed causes unsymmetrical mooring line configuration and generates noticeable effects in the motion responses of the floating body.


Introduction
Many types of offshore structures, including floating islands [1], floating shelters [2], floating wind turbines [3], floating wave energy converters [4] and floating offshore fish farms [5], are moored in complicated coastal environments. These structures are exposed to ocean waves, with the possibility of wave reflection and wave shoaling from coastal seabeds. Dynamic responses of these floating bodies are the combined effects of ocean environments and mooring systems. Accurate predictions of dynamic responses are of practical importance for the design and manufacture of these offshore structures.
Compared with that of deep water offshore structures, design and construction of floating structures in coastal environments face unique challenges and the dynamics of these structures could be more complicated because of the seabed effects. Previous researches suggest that shallow waters can excite larger responses of offshore structures under hydrodynamic loading because of the flat seabed effects [6][7][8]. As numerical approaches for flat seabeds are not directly applicable to sloping seabed profiles, [9] and [10] developed second-body models to account for the sloping seabed effect within the boundary element model frame and found that a sloping seabed significantly influences the cross coupling hydrodynamic coefficients. [11] and [12] developed multi-domain approaches that divide the fluid domain into an interior domain of variable water depth and an exterior domain of constant depth. An extra term accounting for the sloping-bottom effects is introduced to correct the incident wave potential so that the sloping seabed condition is satisfied. They found that the sloping seabed ignored. The investigation of vortices on the mooring line and body motion is out of the scope of the present study.

Coupled Floating Body-Mooring Line Model
For a free-floating body moored in coastal environment, the position and velocity of the floating body vary at each time step and mooring line supplies restoring force to keep the body in acceptable positions. The coupled motion equation for a floating body with mooring lines can be expressed in the form: where M ij andẌ j (t) respectively denote the body mass matrix and body acceleration. F wave i (t), F mooring i (t) and F static i (t) denote wave force, mooring line force and hydrostatic force, respectively. F wave i (t) consists of wave radiation and wave diffraction forces. i, j = 1, 2, 3 stand for surge, heave and pitch motion, respectively. In the present study, the coupling effects between different motion modes are ignored for simplification and therefore Equation (1) can be further expressed as:

Catenary Mooring Line Model
The boundary conditions applied on the mooring lines are the constraints at the two ends of the line: the upper end is attached to the floating platform and is forced to move with the platform, and the lower end is connected to an anchor point at seafloor. For a catenary mooring configuration, a portion of catenary line should lie on the seabed to avoid large vertical forces on the anchors.
The following key assumptions are applied during the analysis of a catenary model: (1) The mooring line moves very slowly such that the drag forces and inertial forces are negligible.
(2) The environmental loads on the mooring line are insignificant and can be excluded in the model.  5) The mooring line is lying on a horizontal seabed. In the present study, the sloping seabed inclination is small and the sloping seabed effects for the mooring line mechanism are therefore ignored. This simplification justifies the application of catenary mooring line model to horizontal seabed profiles. For a catenary mooring line with length L, under its own weight of W per unit length in water, the governing equations can be analytically expressed as [29] (pp. 258-263): where L x and L z are the horizontal and vertical distances between the two ends points of a mooring line, respectively. T x and T z are the horizontal and vertical tension, respectively. A sketch of the catenary mooring line model is presented in Figure 1. At an initial time t = 0, a mooring line pretension T pre is applied to keep the floating body at a static equilibrium position. At each time step, an iterative scheme is applied to calculate the mooring line forces based on Equation (5). This force is used as the exciting mooring line force on the right-hand side of Equation (1) and then the equation is solved using a fourth-order Runge-Kutta method. The floating body displacement works as boundary conditions for the upper end of mooring line. Then, at the next time step, Equation (5) is solved again to start another loop until the end of the numerical simulation. Figure 2 illustrates a coordinate frame OXZ with its origin O at the centre of a two-dimensional floating body-mooring line system in a fluid domain of finite water depth. The fluid domain Ω is bounded by a free surface S f , a body surface S b , a seabed surface S 0 and an enclosing surface at infinity, S ∞ . The fluid flow is assumed to be irrotational and can be described by the velocity potential, φ, which satisfies the Laplace equation as follows,

Fluid Motion Equations
Based on linear theory, the velocity potential φ can be expressed as: where φ I , φ D , φ R and φ re denote the incident wave potential, diffraction potential, radiation potential and reflection potential, respectively. The radiation potential is used to describe the velocity potential generated by floating body motions. For two-dimensional problems, a floating body can undergo surge, heave and pitch motions. For different degrees of freedom, the radiation potential essentially has different values.

Diffraction, Radiation and Reflection Potential Problems
The velocity potentials φ I , φ D , φ R and φ re satisfy the Laplace equation, and φ D , φ R and φ re can be investigated by formulating the boundary integral equations accordingly based on the work of [8,14,30]. The diffraction potential, φ D (x, z, t), describes the fluid field potential generated by the incoming wave within the seabed environment. For an incoming progressive linear wave with a wave amplitude, A, in a finite water depth, h, the velocity potential, φ I and the wave elevation, η I , can be expressed as where g, k and ω denote the gravitational acceleration, wave number and wave frequency, respectively. As no fluid penetrates the rigid body and seabed surfaces, the boundary conditions applied to the body and flat seabed are: where n b and n 0 denote normal vector on body surface and seabed surface, respectively. The linearised kinematic condition and the dynamic free surface boundary condition are given as where η D refers to the wave elevation caused by wave diffraction. The body boundary condition (Equation (11)), seabed condition (Equation (12)) and free surface boundary conditions (Equations (13) and (14)) create a closed boundary value formulation to the diffraction potential problem. The radiation potential, φ R (x, z, t), accounts for the fluid field potential caused by body motions in the fluid domain. Nonpermeable boundary condition is applied to the body and seabed conditions. The free surface boundary condition is simplified as a linear form: where v b denotes the body motion velocity. For a floating body experiencing forced oscillatory motion with an amplitude a, the body motion velocity v b can be described as: and v b = (0, aωcosωt) for heave motion. (20) In the proposed model, the free surface is divided into inner and outer domains, which are discussed in the next section. In the outer domain, the source panels are arranged in an exponentially increasing manner. Using this distribution method, a very large free surface can be covered by the source panels and therefore the generated waves can move towards the far field in this enlarged computational domain before the end of numerical simulation. The radiation condition is satisfied without applying any artificial damping on the free surface.
The sloping seabed provides a source of wave reflections into the fluid domain. In this investigation, the wave reflection potential, φ re , is introduced to satisfy the Laplace Equation (7) whose solution is determined by solving the following boundary value equations Notice that this boundary value equation is analogous to the solution of the wave diffraction potential. In this regard, the potential φ can be considered as the potential caused by the interaction of the floating body (the first body) and the sloping seabed (the second body). Equation (24) describes the sloping seabed boundary condition whereas, in Equation (12), the constant water depth refers to the deeper water depth on the left hand side of the sloping seabed. If the seabed is completely flat, the reflection potential φ re disappears. In this scenario, the reflected wave moves towards the left hand side of the fluid domain and travels into the far field. As the proposed free surface source distribution method covers a very large free surface area on both sides of the body, the reflected wave will not reflect backwards the fluid domain.

Numerical Discretisation of Velocity Potential
The velocity potential φ can be expressed as a boundary integral of Rankine sources continuously distributed on the fluid boundary surface and N 0 panels, respectively such that .., N f and k = 1, ..., N 0 . The discretised form of the boundary integral (Equation (25)) describing the diffraction problem is given as: At any time step t during the numerical simulation, the body and the seabed boundary conditions are described by their boundary value equations. The free surface boundary conditions for diffraction potential φ D , radiation potential φ R and reflection potential φ re are all updated by a 3rd-order Adams-Bashforth scheme. The time interval ∆t = T 0 /100 is applied in the numerical simulations. Here, T 0 is the wave period.
A matrix equation is formulated by discretising Equations (26) The discretisation for the wave radiation, φ R , and the reflection potential, φ re , is exactly the same as that for the diffraction potential, φ D and is therefore omitted herein. Calculation of the Rankine source lnr integral in the boundary integral equations is computationally expensive. In this study, a continuous Rankine source panel integral is developed instead of an isolated Rankine source approximation to reduce numerical errors and increase numerical efficiency.
Note that, physically, the free surface, S f , is of infinite size and the integration over this infinite dimension is not possible by numerical techniques. For this reason, the computational domain is truncated in the numerical simulations and the far-field boundary condition is satisfied on this truncated free surface. In the present study, the Rankine source panels are distributed in an exponential manner, allowing for the coverage of a very large area of free surface. The numerical computations are completed before the incident wave and reflected wave reaches the truncated boundary.

Boundary Panel Distribution
In the present study, the sloping seabed z = h(x) is described mathematically as: where H is the mean water depth and θ denotes the sloping seabed angle which is defined as: where L h and L v are the projected length of the slope seabed in horizontal and vertical directions, respectively. Figure 3 illustrates the definitions of these parameters in the coupled floating body-mooring line problem under investigation. The sloping seabed is divided into three regions, h 1 , h 2 (x) and h 3 . The constants h 1 and h 3 correspond to the constant depth flat portions on the left and right hand sides of the sloping seabed whereas the variable h 2 (x) represents the depth of water in the sloping region. Both the free surface and seabed surface are divided into inner and outer domains. To approximate the inner domain, the source panel distribution method discussed by [8] is adopted here. A parameter γ = 1 + tan θ is introduced to deal with this asymmetric characteristic of the surface. For free surface and flat seabed, the angle θ = 0 and therefore γ = 1. On the right side of S 0 , starting from x = 0, the horizontal distances between neighbouring control points are expressed as: and on the left hand side of S 0 : where L b is the body surface panel length. The parameters α k define the separation distances and are decided according to the method discussed by [8]. In the outer domain the distances between the centres of neighbouring source panels are displaced in an exponentially increasing manner as aforementioned.

Model Validation
In order to validate the proposed coupled floating body-mooring line model, numerical simulations are performed for the static offset test, the free decay test without mooring line, the free decay test with mooring line and the regular wave simulation. All numerical simulations are conducted in the time domain and converted into the frequency domain for direct comparison against available published data.

Static Offset Test
The stiffness characteristics of catenary mooring lines are investigated preliminarily by the static offset tests. For a mooring line with a length of L = 100 m and a unit weight of W = 828 N/m, moored in water depth h = 25 m, Figure 4 shows the comparison of this mooring line stiffness curves with respect to horizontal force T x (a) and vertical force T z (b) between the proposed numerical method and [29] (pp. 258-263). Results from the proposed method show excellent agreement with those of [29]. It is noticed that both the horizontal force T x and vertical force T z demonstrate nonlinear characteristics as the mooring line horizontal distance L x increases and this phenomenon is also called nonlinear "hardening" behaviour by [26].

Free Decay Test
The free decay test works as a measurement to determine the natural period T and the damping ratio β of the system. In the free decay test, the floating system was offset an initial displacement and allowed to damp out until reaching the equilibrium position. A heave decay test was performed for a semi-circular cylinder with radius 10 m excluding mooring line as a preliminary study. Figure 5 shows comparison of the heave decay test for the semi-circular cylinder between the proposed method, the [31] numerical simulation and the [32] experiment. As shown, the proposed model yields a good agreement compared with the numerical prediction of [31] and experimental data of [32]. The natural period for this cylinder in heave motion is 7.18 s.
Numerical simulations are further performed for the semi-circular cylinder with different mooring line stiffness K i . The mooring line system stiffness K i is defined as: where X i denotes the body offset and is corresponding to the body displacement in Equation (1). Figure 6 shows the surge decay test for the cylinder with mooring line stiffness K 1 = 90,000 N/m, 135,000 N/m, 180,000 N/m and 225,000 N/m, respectively. It is noticed that the decay period and motion amplitude decrease as stiffness K i increases. The coupled system free decay period can be analytically evaluated by the following equation of [33] (pp. 275-278): where subscripts are omitted for simple notation. The body added mass a and damping b can be obtained from the work of [8]. Table 1 shows a comparison of the surge decay period T between the analytical method by [33] and the proposed numerical method for the mooring line stiffness of K 1 = 90,000 N/m, 135,000 N/m, 180,000 N/m and 225,000 N/m, respectively. It is observed that the proposed method agrees well with the analytical expression of [33]. The natural period of this floating body-mooring system is around 10s. Therefore, the wave frequencies dominate the system responses and low-frequency responses are expected to be small.

Regular Wave Simulation
Numerical simulations are further carried out for the floating body-mooring system subject to regular waves. As a preliminary study, numerical simulations are conducted for a rectangular cylinder with breadth B = 0.30 m and draft d = 0.20 m, floating in water depth h = 0.40 m without mooring lines attached. Figure 7 shows a comparison of the time records of surge motion X 1 (t) (a) and heave motion X 2 (t) (b) between the proposed method, the experimental data by [34] and the numerical data by [35]. The floating rectangular cylinder is subject to a regular wave with wave period T 0 = 1.2 s and wave amplitude A = 0.02 m. It is shown that the proposed model results agree favourably well with the model test data of [34] and the numerical prediction of [35]. This floating body is moored by two symmetrical mooring lines, with L = 1.6 m and a weight of W = 39.2 N/m in water. Figure 8 demonstrates the comparison of time records of surge motion X 1 (t) (a) and heave motion X 2 (t) (b) between proposed method results and [35] numerical data for the rectangular cylinder with the mooring effect included. The numerical data of [35] is shifted with T = 0 as the starting time for better comparison. It is noticed that the proposed model data agrees reasonably well with the data of [35] for both surge and heave motions. The discrepancy can be explained by the fact that the viscous effects are included in the work of [35] but excluded in the present model.  Comparison of time records of surge motion X 1 (t) (a) and heave motion X 2 (t) (b) between proposed method results and [34] experimental data and [35] numerical data without mooring effect.

Case Study
For a floating body moored in coastal areas, if the mooring line anchors are located in different water depths, then the mooring line profiles are not symmetrical. The asymmetry of seabed and mooring line profiles significantly affects fluid domain and hence the body response characteristics. By comparing the asymmetrical seabed and moorings with the flat seabed and symmetrical mooring lines, this section investigates the floating body response characteristics in sloping seabed and asymmetrical mooring lines conditions. A rectangular cylinder with breadth B = 20 m and draft d = 10 m is selected as the floating body. Such a body size is typical of floating bodies in ocean and offshore engineering. An average water depth of H = 200 m is considered for the first numerical case. Two types of seabed and mooring line configurations are used in the numerical simulations. The first type is a flat seabed with two identical mooring lines, and the second one is a sloping seabed with two different mooring line arrangement. Figure 9 shows the profiles of symmetrical and asymmetrical seabeds with mooring lines. For the flat seabed condition, the left-side mooring line length L 1 and anchor horizontal distance AL 1 equal the right-side values L 2 and AL 2 . The average water depths H for the sloping seabeds are same as the flat ones for direct comparisons. Each mooring line has a weight of W = 5 kN/m in water and two mooring lines are used on each side in the present study unless otherwise stated.  Table 2 shows parameters of the four sets of mooring lines for the symmetrical seabed condition. In the present study, 10% of the mooring line length are lying on the seabed for all the numerical simulations unless stated otherwise. The incident wave has a wave amplitude A = 1 m and the mooring line is kept in a catenary profile in all numerical simulations. Free decay test shows the heave nature period for this coupled body and mooring line system are 8.78 s. Figure 10 shows a comparison of floating body surge motion RAO 11 (ω) between without mooring line condition, mooring line set 1, mooring line set 2, mooring line set 3 and mooring line set 4. It is observed that the mooring line reduces the surge motion RAO 11 (ω) noticeably and the RAO 11 (ω) becomes smaller as the mooring line length increases. Figure 11 shows a comparison of the heave motion RAO 22 (ω) of the floating body between without mooring line condition, mooring line set 1, mooring line set 2, mooring line set 3 and mooring line set 4. It is noticed that the mooring line reduces the heave motion amplitude significantly. Longer mooring line generates larger restoring force for the floating body and therefore the heave motion amplitude reduces more. The peak value of motion amplitude is shifted towards the high wave frequency as the mooring line length increases. This finding is in line with that of [36].  In order to investigate the hydrodynamic effects of asymmetrical seabed and mooring line, a symmetrical mooring line set L 1 = L 2 = 700 m is taken as the basis. Two asymmetrical mooring line sets are used in the numerical simulations and the parameters for these two sets of mooring lines are presented in Table 3. For direct comparison, the total mooring line lengths are same for both symmetrical and asymmetrical mooring line sets. It should be emphasised that the floating body has to be at static equilibrium. This requirement is naturally satisfied for symmetrical mooring line set. For asymmetrical seabed and mooring line set, the left side mooring line is arranged similarly as the symmetrical case, i.e., 10% of the mooring line length are laid on the seabed. The right-side mooring line has to be arranged to balance the horizontal mooring line force generated by the left side. An iterative procedure is developed to calculate the anchor horizontal distance AL 2 and the anchor vertical distance AH 2 is readily available according to the seabed slope once AL 2 is obtained. In such arrangement, the difference of vertical mooring line forces generated by the left and right sides is less than 3% and has marginal effects on the pitch motion of the floating body.  Figure 12 shows a comparison of the surge motion RAO 11 (ω) of the floating body between symmetrical mooring line condition, asymmetrical mooring line set 1 and set 2. It is noticed that the surge motion RAO 11 reduces appreciably as the asymmetry level increases across the whole wave frequency range. Compared with the symmetrical case, the left-side length of asymmetrical mooring line set 2 increases by 6%, but it reduces the RAO 11 by approximately 15%. If less mooring line length is allowed to lie on the seabed, effects due to the mooring line asymmetry are expected to be more significant. It should be noticed that the difference between the left-and right-side mooring lines is in a reasonable range, otherwise the mooring line set can not keep the floating body at a static equilibrium position. The heave motion RAO 22 (ω) is marginally affected by the asymmetrical mooring line setup and therefore the data are not presented here. Numerical simulations are further performed for the floating body positioned in a water depth of H = 20 m. Both flat seabed with symmetrical mooring lines and sloping seabed with asymmetrical mooring lines are considered. The mooring lines on the sloping seabed have length L 1 = L 2 but are arranged in asymmetrical position. These mooring lines are lying on flat seabed portion whereas sloping seabed portion is located within two anchor points. Both the symmetrical and asymmetrical mooring line configurations have the same length L 1 = L 2 = 60 m for direct comparison. The 15% mooring line length is allowed to lie on the seabed for both cases. Figure 13 shows the profile of an asymmetrical seabed with mooring lines positioned in the water depth H = 20 m. The sloping seabed effect for the incident wave amplitude is first investigated, and no floating body or mooring lines are included in the fluid domain. A linear incident wave with a wave amplitude of A inc comes from the left-hand side and moves towards the right-hand side as illustrated in Figure 13. Figure 14 shows the incident wave amplitude A on the right-and left-hand sides over a sloping seabed with sloping angle tanθ = 1/8. The left-hand side has a wave amplitude of A = A inc as the incident wave comes from the left. The wave amplitude changes gradually on the sloping part due to the sloping seabed effects. The wave becomes steady when it moves into the flat part on the right-hand side. It is noticed that the incident wave amplitude increases due to the sloping seabed and the sloping seabed has largest effect for the incident wave when the wave frequency is approximately kd = 0.8. The left-side mooring line weight has a weight of W = 5 kN/m in water for the sloping seabed case, and the floating body is at a static equilibrium. A preliminary iterative procedure is developed to calculate the right-side mooring line weight per meter to reach static equilibrium. Figure 15 shows a comparison of the surge motion RAO 11 (ω) of the floating body between a flat seabed and a sloping seabed with tanθ = 1/20, tanθ = 1/10 and tanθ = 1/8, respectively. It is noticed that the surge motion RAO 11 (ω) is smaller for the flat seabed than the sloping seabed conditions and RAO 11 (ω) increases as the sloping seabed angle θ increases across the frequency range.  Figure 16 shows a comparison of the heave motion RAO 22 (ω) of the floating body between flat seabed, sloping seabed with tanθ = 1/20, tanθ = 1/10 and tanθ = 1/8, respectively. It is noticed that the peak responses of the heave motion RAO 22 (ω) are shifted towards lower wave frequency and the peak response values increase along with the sloping seabed angle. The phenomenon presented in Figures 15 and 16 implies that the sloping seabed with asymmetrical mooring line arrangement produces less mooring stiffness than does the symmetrical case. In the preliminary study to investigate the right-side mooring weight, numerical simulations show that less weight is required compared with the left side.  Figure 17 shows a comparison of the time record of the floating-body heave motion X 2 (t) for wave frequency kd = 1.0 between flat seabed, sloping seabed with tanθ = 1/20, tanθ = 1/10 and tanθ = 1/8, respectively. The heave motion X 2 (t) has the largest amplitude for flat seabed condition and the X 2 (t) motion amplitudes decrease as the sloping angles increase. These time-domain results are aligned with those frequency-domain results presented in Figure 16. In this study, the fluid domain is described by a linear boundary element model. The peak value for heave motion X 2 (t) is around 0.7 whereas the trough value is about 0.5. This nonharmonic body motion track demonstrates the nonlinear characteristics of the coupled floating body-mooring line system.   (t) has the largest amplitude for flat seabed condition and F mooring 2 (t) amplitudes decrease as the sloping angles increase. It is also noticed that the force peak values are smaller than the trough values and this phenomenon is due to the fact that the heave motion X 2 (t) peak values are bigger than the trough values. Figure 19 demonstrates a comparison of time record of floating body surge motion X 1 (t) for wave frequency kd = 1.0 between flat seabed, sloping seabed with tanθ = 1/20, tanθ = 1/10 and tanθ = 1/8, respectively. The surge motion X 1 (t) peak values are almost same for these different seabed conditions but the X 1 (t) has the least value for flat seabed condition. It is also noticed that the nonlinear characteristics become more obvious as the seabed inclination increases.   (ω) on the moored body excited by an incident wave for the slope seabed angle of tan θ = 1/8. It is noticeable that the wave force F wave 2 (ω) increases with the wave frequency for kd < 0.6. This wave force F wave 2 (ω) reaches its peak in the period 0.6 < kd < 0.7 and decreases gradually with the wave frequency for kd > 0.7. This variation trend is in line with the moored body motion as presented in Figure 16. It should be noticed that the coupled floating body and mooring line system natural frequency is within the wave frequency range and only the wave-frequency responses are investigated in this paper. The low-frequency responses are very important for large floating bodies.

Conclusions
A two-dimensional coupled floating body-mooring line model is developed to study the freely floating body motion responses in both flat and sloping seabed conditions. A continuous Rankine source based boundary element model is established to describe the fluid domain. This model comprises three time-domain boundary integral equations accounting for diffraction problem, radiation problem and reflection problem, respectively. The mooring line is formulated by a catenary mooring line model and coupled with the boundary element model at each time step.
Numerical validations were carried out for the static offset test, free decay test without mooring line, free decay test with mooring line and regular wave simulation. Reasonable agreement with published data demonstrates the accuracy of proposed model. Two numerical cases are investigated with an emphasis on the hydrodynamic effects of sloping seabed and asymmetrical mooring line. In the first numerical study, both the surge motion RAO 11 and the heave motion RAO 22 decrease as the mooring line length increases for flat seabed condition. For the sloping seabed case, the left and right sides of mooring line have different lengths and are laid on sloping seabed portion. The right-side mooring line is positioned to keep the floating body at static equilibrium. Numerical results show that the surge motion RAO 11 decreases as the asymmetrical level increases across the whole frequency range but the heave motion RAO 22 keeps constant regardless of the asymmetry level. These numerical findings indicate that in such sloping seabed conditions, the mooring line system shows better station-keeping capability than it does in flat seabed conditions.
For the second numerical case, both the left and right sides of mooring lines are laid on flat seabed portion with sloping portion between two anchor points. The left and right mooring line has same length but different weight to keep the floating body in static equilibrium position. Numerical study shows that both the surge motion RAO 11 and heave motion RAO 22 increase as the sloping seabed angle becomes larger. The coupled floating body-mooring line system demonstrates nonlinear characteristics for the body motion and mooring line force responses in the time domain. The surge motion has clearer nonlinear performance than the heave motion. An asymmetrical mooring line configuration demonstrates less station-keeping capability than an symmetrical configuration, but it also requires less mooring line weight and steel consumption. From an economical viewpoint, such a mooring line configuration could be beneficial. In actual industrial projects, the design of asymmetrical mooring line is a trade-off between station-keeping capability and economical benefit. This paper deals with two-dimensional cases and therefore has academic values, while practical engineering problems are three-dimensional. The study of three-dimensional coupled floating body-mooring line problem and an extention to wave energy converters will be part of our future work.