Multibody Analysis of Wave Energy Converters and a Floating Platform in the Time Domain

: Generally, new and renewable energy systems generate electricity by installing and operating multiple modules simultaneously. In the Republic of Korea, recent studies and developments have focused on asymmetric wave energy converters (hereafter referred to as rotors) suitable for marine environments off the western coast of Jeju. These rotors are arranged on a large floating truss-structure platform and designed to harness electricity from the rotors’ pitch motion. However, when multiple rotors operate on a platform, their behavior diverges from that of a single module due to hydrodynamic interactions between them. Moreover, because the rotors are connected to the floating platform, their motion is influenced by the platform’s dynamics. In this study, a time-domain multibody motion equation was established to analyze changes in the behavioral characteristics of the rotors, both with and without a floating platform. The hydrostatic and hydrodynamic coefficients were derived in the frequency domain using WAMIT, a commercial code based on linear potential flow theory for three-dimensional diffraction/radiation analyses. The motion equation was then applied under regular and irregular wave conditions using OrcaFlex version 11.3, a marine systems design and analysis program. The resulting behaviors were compared to elucidate the influence of the platform and hydrodynamic interactions on the rotors’ performance.


Introduction
Active research on power generation systems utilizing new and renewable energy sources is underway, driven by the global initiative to achieve decarbonization.Recently, many research studies have expanded from terrestrial to marine environments to overcome spatial and societal constraints [1].The ocean harbors a variety of renewable energy sources, including tides, ocean currents, waves, wind power and solar heat, as well as the energy potential from differences in seawater temperature and salinity.Consequently, this facilitates a wide range of research opportunities in various fields.Typically, renewable energy systems enhance power production within confined spaces by installing and operating multiple modules simultaneously.This process also applies to wave energy converters (WECs) that produce electrical energy from wave energy.However, when multiple WECs are operated simultaneously, their motion response differs from that of a single entity due to hydrodynamic interactions between each unit.
In 2016, the research project "Development of Fundamental Technologies for a 1 MW Class Arrayed-Buoy Wave Energy Converter in Various Sea Conditions" was undertaken in South Korea, where Salter's duck rotor, as proposed by Stephen Salter [2][3][4][5], was utilized.Subsequently, this model (hereafter referred to as a rotor) was redesigned and developed to suit the ocean environment in the western part of Jeju, South Korea (see Figure 1) [6][7][8][9].In previous studies, a one-degree-of-freedom system was constructed.An analysis of this simplified system, featuring a rotor connected to a fixed rotational axis submerged at a specific depth, was then performed.
J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 2 of 20 1) [6][7][8][9].In previous studies, a one-degree-of-freedom system was constructed.An analysis of this simplified system, featuring a rotor connected to a fixed rotational axis submerged at a specific depth, was then performed.However, in a real sea, multiple rotors are designed to extract power by connecting them to a large floating truss platform.In this design, understanding the hydrodynamic interaction between the floating platform and the multiple rotors becomes crucial.Furthermore, given that the rotors are connected to the floating platform, the dynamics of the floating platform can affect the rotors' motion response.In particular, the motion response of the rotors connected to the floating platform differs marginally from that of a single rotor connected to a fixed axis of rotation.
Numerous studies have focused on multibody interactions, primarily exploring the hydrodynamic interactions between WECs and floating structures [10][11][12][13][14][15][16].Konispoliatis et al. [17] studied the effect that a fixed vertical breakwater has on the hydrodynamic characteristics of multiple WECs parametrically.A notable study delved into the dynamic behaviors of a torus-type WEC which exhibits a heaving motion in conjunction with a spartype floating body [18].Zhou et al. [19] studied the dynamic behaviors of multiple heaving WECs and semisubmersible platforms, revealing that a prominent new power peak emerges when the platform and WECs are in a synchronized state.Furthermore, various parametric studies were conducted to optimize the system by introducing frequency-domain-coupled hydrodynamics, considering the constraints and power output through the relative motion between the integrated WECs and the semi-submersible platform [20].Kamarlouei et al. [21] explored a hybrid system comprising 12 WaveStars surrounding a semisubmersible platform, revealing that the platform's heaving and pitching motions could be effectively mitigated by WECs with appropriately tuned damping in the power take-off (PTO) system.Ghafari et al. [22] examined the effect of the number of WECs on a hybrid system's performance.They indicated that a large number of WECs contributed to system stabilization in terms of pitch.Moreover, the power output improved when the WECs' movements aligned in phase with those of the platform.Cheng et al. [23] validated the feasibility of different types of WECs integrated into floating platforms and showed their synergy.Yazdi et al. [24] conducted a parametric study on the dynamic responses of multiple Salter's duck WECs and a braceless floating semisubmersible, focusing on varying numbers of WECs, wave periods and wave heights.
In this study, a time-domain coupled equation of motion for a multibody analysis was formulated to ascertain differences in motion response.In contrast to the majority of However, in a real sea, multiple rotors are designed to extract power by connecting them to a large floating truss platform.In this design, understanding the hydrodynamic interaction between the floating platform and the multiple rotors becomes crucial.Furthermore, given that the rotors are connected to the floating platform, the dynamics of the floating platform can affect the rotors' motion response.In particular, the motion response of the rotors connected to the floating platform differs marginally from that of a single rotor connected to a fixed axis of rotation.
Numerous studies have focused on multibody interactions, primarily exploring the hydrodynamic interactions between WECs and floating structures [10][11][12][13][14][15][16].Konispoliatis et al. [17] studied the effect that a fixed vertical breakwater has on the hydrodynamic characteristics of multiple WECs parametrically.A notable study delved into the dynamic behaviors of a torus-type WEC which exhibits a heaving motion in conjunction with a spar-type floating body [18].Zhou et al. [19] studied the dynamic behaviors of multiple heaving WECs and semisubmersible platforms, revealing that a prominent new power peak emerges when the platform and WECs are in a synchronized state.Furthermore, various parametric studies were conducted to optimize the system by introducing frequencydomain-coupled hydrodynamics, considering the constraints and power output through the relative motion between the integrated WECs and the semi-submersible platform [20].Kamarlouei et al. [21] explored a hybrid system comprising 12 WaveStars surrounding a semisubmersible platform, revealing that the platform's heaving and pitching motions could be effectively mitigated by WECs with appropriately tuned damping in the power take-off (PTO) system.Ghafari et al. [22] examined the effect of the number of WECs on a hybrid system's performance.They indicated that a large number of WECs contributed to system stabilization in terms of pitch.Moreover, the power output improved when the WECs' movements aligned in phase with those of the platform.Cheng et al. [23] validated the feasibility of different types of WECs integrated into floating platforms and showed their synergy.Yazdi et al. [24] conducted a parametric study on the dynamic responses of multiple Salter's duck WECs and a braceless floating semisubmersible, focusing on varying numbers of WECs, wave periods and wave heights.
In this study, a time-domain coupled equation of motion for a multibody analysis was formulated to ascertain differences in motion response.In contrast to the majority of previous studies that considered only a single degree of freedom for the rotor, this study incorporated the surge and heave effects interacting with the floating body.A numerical model was developed to include the hydrodynamic interaction effects between the rotors and the platform, adhering to necessary constraints such as connectivity.Then, the motion response was calculated through a regular wave simulation.An irregular wave simulation was also performed to calculate the motion responses and extract power at the installation site.Finally, the results were compared with those obtained from a rotor connected to a fixed rotational axis.This approach facilitated an assessment of the in-wave power generation performance of the actual wave energy structure.

System Identification
In this study, a numerical analysis was conducted on the Salter's duck model depicted in Figure 2. The rotor was designed with a distinct shape that varied between the front (facing the direction of the incident wave) and the rear sides to enhance energy absorption efficiency [25].Previous studies conducted parametric analyses focusing on the initial beak angle and ballast weight based on these geometric characteristics.Subsequently, the rotor was redesigned to optimize its performance in the ocean environment west of Jeju Island, South Korea, with a water depth of 80 m.As illustrated in Figure 2a, the rotor achieves static equilibrium by itself when floating in calm water (without waves).In this state, the rotor's rotational center is positioned 1.6 m below the water surface, with a draft of 3.6 m and a beak angle of 60 previous studies that considered only a single degree of freedom for the rotor, this study incorporated the surge and heave effects interacting with the floating body.A numerical model was developed to include the hydrodynamic interaction effects between the rotors and the platform, adhering to necessary constraints such as connectivity.Then, the motion response was calculated through a regular wave simulation.An irregular wave simulation was also performed to calculate the motion responses and extract power at the installation site.Finally, the results were compared with those obtained from a rotor connected to a fixed rotational axis.This approach facilitated an assessment of the in-wave power generation performance of the actual wave energy structure.

System Identification
In this study, a numerical analysis was conducted on the Salter's duck model depicted in Figure 2. The rotor was designed with a distinct shape that varied between the front (facing the direction of the incident wave) and the rear sides to enhance energy absorption efficiency [25].Previous studies conducted parametric analyses focusing on the initial beak angle and ballast weight based on these geometric characteristics.Subsequently, the rotor was redesigned to optimize its performance in the ocean environment west of Jeju Island, South Korea, with a water depth of 80 m.As illustrated in Figure 2a, the rotor achieves static equilibrium by itself when floating in calm water (without waves).In this state, the rotor's rotational center is positioned 1.6 m below the water surface, with a draft of 3.6 m and a beak angle of 60°.The floating platform, a semi-submersible truss structure shown in Figure 2b, also maintains a static equilibrium state with a draft of 5.35 m in calm waters.Figure 2c shows a central platform oriented toward the direction of the incident wave which can connect four rotors.The entire platform hosts a total of 20 rotors, of which 16 rotors are split into two groups and connected to the two wing platforms on the left and right sides.The angle (α) between the central and wing platforms is adjustable based on environmental conditions and is set at 150 • during operation.The spacing between each rotor is 3.5 m, and the floating platform is connected by two hawsers to a weathervaning buoy moored by multiple lines, allowing for orientation adjustments according to changes in the incident wave direction (see Figure 3).The floating platform, a semi-submersible truss structure shown in Figure 2b, also maintains a static equilibrium state with a draft of 5.35 m in calm waters.Figure 2c shows a central platform oriented toward the direction of the incident wave which can connect four rotors.The entire platform hosts a total of 20 rotors, of which 16 rotors are split into two groups and connected to the two wing platforms on the left and right sides.The angle (α) between the central and wing platforms is adjustable based on environmental conditions and is set at 150° during operation.The spacing between each rotor is 3.5 m, and the floating platform is connected by two hawsers to a weathervaning buoy moored by multiple lines, allowing for orientation adjustments according to changes in the incident wave direction (see Figure 3).In this study, to increase computational efficiency by reducing the degrees of freedom of the entire system, as shown in Figure 2c, the rotor connected to the wing platform was modeled as a single rigid body along with the floating platform, and its individual movements were not considered.Consequently, only the four rotors on the central platform were allowed to move.In conducting a multibody analysis of one floating platform and four rotors, the mass characteristics of the remaining 16 rotors were incorporated into In this study, to increase computational efficiency by reducing the degrees of freedom of the entire system, as shown in Figure 2c, the rotor connected to the wing platform was modeled as a single rigid body along with the floating platform, and its individual movements were not considered.Consequently, only the four rotors on the central platform were allowed to move.In conducting a multibody analysis of one floating platform and four rotors, the mass characteristics of the remaining 16 rotors were incorporated into the floating platform.Hereinafter, the term "platform" refers to the platform structure, including 16 rotors on the wing platform.The specifications of the single rotor and the entire platform, excluding the four central platform rotors, are detailed in Tables 1 and 2, respectively.13 (1) CoR: center of rotation.
Table 2. Specifications of the platform.

Multibody Equation of Motion
In the initial design stages of a floating platform, the system is typically assumed to be linear, and its response in the frequency domain is calculated.This approach is advantageous due to its relative simplicity and the ability to yield results quickly.However, to more accurately simulate behavior in the actual ocean environment, it becomes essential to use a time-domain analysis which can account for nonlinear effects.Specifically, the rotor in this study exhibits nonlinear restoring characteristics, with the restoring moment varying rapidly based on the rotational angle, necessitating an analysis in the time domain.
For this study, the equation of motion was established using Newton's second law, applying the Cummins equation [26], as shown in Equation (1).Notably, a multibody coupled equation of motion was formulated to express the behavior of multiple rotors connected to a floating platform.
The equation considers the hydrodynamic interactions between each floating body when analyzing the motion of two or more adjacent floating bodies.Consequently, all terms in the equation are presented in the form of a matrix or column vector, incorporating (6 × N) degrees of freedom corresponding to the number of floating bodies (N).In the next notation, "pt f m" and "rt − n" used in the superscript of each variable denote the respective floating body as the platform and the nth rotor.Furthermore, (i, ), used in the subscript, indicate the degrees of freedom of the platform and the nth rotor, respectively.
M is the sum of two matrices, M iner and A ∞ .M iner includes the mass and mass moment of inertia, and A ∞ includes the added mass and the added mass moment of inertia at infinite frequency as components.B vis and B pto denote matrices for the linear viscous damping coefficient and the PTO damping coefficient, respectively.K res represents the coefficient matrix for the restoring load of the floating body (the platform and rotor), while K m is the coefficient matrix for the restoring load induced by the mooring line connected to the platform.L w is the wave excitation load and ..

ξ,
. ξ and ξ represent the acceleration, velocity and displacement (or angular acceleration, angular velocity and angular displacement) of the platform and rotor, respectively. ..

ξ,
. ξ, ξ and L w take the form of column vectors.
In order to depict the interdependence among the degrees of freedom of the floating body, the matrix components of M were assembled as illustrated in Equation ( 2).L rd .ξ represents the radiation damping load and can be calculated using a retardation function matrix B, as depicted in Equation ( 3).The matrix B is expressible as shown in Equation ( 4).
In Equation ( 1), the coefficients of the restoring load, hydrodynamic load (A ∞ , b rad ) and wave excitation load for the floating body were derived using Wave Analysis MIT (WAMIT version 7.1 [27]), a three-dimensional diffraction/radiation solver based on the linear potential theory.The above coefficients are a function of the incident wave frequency, and in this study, the coefficients were obtained for a total of 81 frequencies at intervals of 0.02 rad/s for incident wave frequencies ranging from 0.02 rad/s to 1.70 rad/s.For WAMIT calculations, geometrical details of the submerged portion of the floating body are necessary; these were numerically modeled using Multisurf version 9.0 [28] (Figure 4).For simplification in the hydrodynamic analysis of the floating platform, the truss structure atop the pontoon was omitted. is the sum of two matrices,  and  . includes the mass and mass moment of inertia, and  includes the added mass and the added mass moment of inertia at infinite frequency as components. and  denote matrices for the linear viscous damping coefficient and the PTO damping coefficient, respectively. represents the coefficient matrix for the restoring load of the floating body (the platform and rotor), while  is the coefficient matrix for the restoring load induced by the mooring line connected to the platform. is the wave excitation load and  ,  and  represent the acceleration, velocity and displacement (or angular acceleration, angular velocity and angular displacement) of the platform and rotor, respectively. ,  ,  and  take the form of column vectors.
In order to depict the interdependence among the degrees of freedom of the floating body, the matrix components of  were assembled as illustrated in Equation ( 2).  represents the radiation damping load and can be calculated using a retardation function matrix , as depicted in Equation ( 3).The matrix  is expressible as shown in Equation (4).
In Equation ( 1), the coefficients of the restoring load, hydrodynamic load ( ,  ) and wave excitation load for the floating body were derived using Wave Analysis MIT (WAMIT version 7.1 [27]), a three-dimensional diffraction/radiation solver based on the linear potential theory.The above coefficients are a function of the incident wave frequency, and in this study, the coefficients were obtained for a total of 81 frequencies at intervals of 0.02 rad/s for incident wave frequencies ranging from 0.02 rad/s to 1.70 rad/s.For WAMIT calculations, geometrical details of the submerged portion of the floating body are necessary; these were numerically modeled using Multisurf version 9.0 [28] (Figure 4).For simplification in the hydrodynamic analysis of the floating platform, the truss structure atop the pontoon was omitted.The coupled equation of motion assumes that both the viscous damping load and the PTO damping load are linear, expressed as loads proportional to the velocity of the floating body.Therefore, they are combined as in Equation (5).
The component of  can be calculated using the damping ratio () obtained from free decay tests, with the assumption that total damping is the sum of linear viscous damping and radiation damping during the free decay test [29], as detailed in Equation (6).Free decay tests for the rotor were conducted using a 1/11 scale standalone model in a two-dimensional wave tank at Jeju National University, while platform free decay tests utilized a 1/30 scale model in a three-dimensional ocean engineering basin at KRISO (Korea Research Institute of Ships and Ocean Engineering).The damping ratios obtained from these experiments are presented in Table 3.In Equation ( 6),  , represents the mass or mass moment of inertia for translational or rotational motion, respectively, and  , denotes the undamped natural frequency of the degree of freedom .  , and  ,  , correspond to the added mass and radiation damping coefficients at the undamped natural frequency.The undamped natural period values ( , = 2/ , ) are specified in Tables 1 and 2. In general, when calculating the linear damping coefficient of a floating body, the assumption is made that the influence of each degree of freedom is more significant than the coupling effect between degrees of freedom.Consequently, only the diagonal component  , is The coupled equation of motion assumes that both the viscous damping load and the PTO damping load are linear, expressed as loads proportional to the velocity of the floating body.Therefore, they are combined as in Equation ( 5).
The component of B vis can be calculated using the damping ratio (κ) obtained from free decay tests, with the assumption that total damping is the sum of linear viscous damping and radiation damping during the free decay test [29], as detailed in Equation (6).Free decay tests for the rotor were conducted using a 1/11 scale standalone model in a two-dimensional wave tank at Jeju National University, while platform free decay tests utilized a 1/30 scale model in a three-dimensional ocean engineering basin at KRISO (Korea Research Institute of Ships and Ocean Engineering).The damping ratios obtained from these experiments are presented in Table 3.In Equation ( 6), M iner,ii represents the mass or mass moment of inertia for translational or rotational motion, respectively, and ω 0,i denotes the undamped natural frequency of the degree of freedom i.A ii (ω 0,i ) and b rad,ii (ω 0,i ) correspond to the added mass and radiation damping coefficients at the undamped natural frequency.The undamped natural period values (T 0,i = 2π/ω 0,i ) are specified in Tables 1 and 2. In general, when calculating the linear damping coefficient of a floating body, the assumption is made that the influence of each degree of freedom is more significant than the coupling effect between degrees of freedom.Consequently, only the diagonal component B vis,ii is calculated, while the off-diagonal component B vis,ij is set to zero.In the actual free decay experiment, the setup allowed for the occurrence of only one degree of freedom at a time, whether it was surge, heave or pitch.
In this study, the rotor generates power through pitch motion.Given that the rotor and platform share the same axis of rotation, the PTO mechanism must be expressed in terms of the relative rotational velocity between the rotor and the platform.The PTO damping coefficient thus appears in individual pitch terms for both the rotor and platform, as well as in their pitch coupling terms.The PTO damping coefficient utilizes the optimal value derived from previous research.
Floating bodies typically exhibit stability against heave, roll and pitch.Nevertheless, when it comes to horizontal motions like surge, sway and yaw, stability needs to be ensured by employing supplementary devices, such as a mooring system.To achieve this, a term representing the restoring load was formulated, as illustrated in Equation (7).
A linear hydrostatic coefficient was applied to the K res of the floating platform.However, applying a linear coefficient to K res corresponding to the rotor's pitch motion poses challenges for obtaining an accurate response due to rapid changes in the volume and center of buoyancy of the rotor's submerged portion depending on its rotational angle.The restoring moment of the rotor exhibits significant nonlinearity with changing rotational displacement, as illustrated in Figure 5. Therefore, the restoring moment coefficient of the rotor is set to change nonlinearly in accordance with its angular displacement.In the time-domain simulation, the rotor's rotational displacement is evaluated at each time step, and the corresponding restoring moment coefficient, as determined from Figure 5, is applied to the equation of motion.
J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 8 of 20 calculated, while the off-diagonal component  , is set to zero.In the actual free decay experiment, the setup allowed for the occurrence of only one degree of freedom at a time, whether it was surge, heave or pitch.In this study, the rotor generates power through pitch motion.Given that the rotor and platform share the same axis of rotation, the PTO mechanism must be expressed in terms of the relative rotational velocity between the rotor and the platform.The PTO damping coefficient thus appears in individual pitch terms for both the rotor and platform, as well as in their pitch coupling terms.The PTO damping coefficient utilizes the optimal value derived from previous research.
Floating bodies typically exhibit stability against heave, roll and pitch.Nevertheless, when it comes to horizontal motions like surge, sway and yaw, stability needs to be ensured by employing supplementary devices, such as a mooring system.To achieve this, a term representing the restoring load was formulated, as illustrated in Equation (7).
A linear hydrostatic coefficient was applied to the   of the floating platform.However, applying a linear coefficient to   corresponding to the rotor's pitch motion poses challenges for obtaining an accurate response due to rapid changes in the volume and center of buoyancy of the rotor's submerged portion depending on its rotational angle.The restoring moment of the rotor exhibits significant nonlinearity with changing rotational displacement, as illustrated in Figure 5. Therefore, the restoring moment coefficient of the rotor is set to change nonlinearly in accordance with its angular displacement.In the time-domain simulation, the rotor's rotational displacement is evaluated at each time step, and the corresponding restoring moment coefficient, as determined from Figure 5, is applied to the equation of motion.The influence of the mooring line is typically addressed through a floating-bodymooring coupled analysis method using a finite element model wherein the inertia and drag forces acting on the mooring line must be considered.However, given that this study specifically focused on the coupled motion between the floating platform and the rotors, the mooring system was simplified to an equivalent spring system independently connected for each degree of freedom of the platform to represent its fundamental roles, maintaining the targeted natural period and ensuring platform stability.The equivalent spring stiffness coefficient (  ) of the platform was determined from Equation (8), thereby enabling the platform to achieve the undamped natural period specified in Table 2.The influence of the mooring line is typically addressed through a floating-bodymooring coupled analysis method using a finite element model wherein the inertia and drag forces acting on the mooring line must be considered.However, given that this study specifically focused on the coupled motion between the floating platform and the rotors, the mooring system was simplified to an equivalent spring system independently connected for each degree of freedom of the platform to represent its fundamental roles, maintaining the targeted natural period and ensuring platform stability.The equivalent spring stiffness coefficient (K m ) of the platform was determined from Equation (8), thereby enabling the platform to achieve the undamped natural period specified in Table 2.
Equation ( 8) derives from the undamped natural frequency ω 0,i = K ii /(M ii + A ii (ω 0,i )), with K res subtracted from the total restoring coefficient K ii to prevent double-counting the hydrostatic restoring stiffness of the platform in the mooring stiffness K m .
Furthermore, the rotors are linked to their corresponding rotational axes on the floating platform.As a result of this assembly, the pitch motion of the rotor was allowed freedom, while the other degrees of freedom are bound to the floating platform, moving in unison.In order to facilitate the numerical implementation of this structure, constraint conditions for the floating body were formulated, as detailed in Equation (9).
Utilizing the constraint coefficient C link indicates that identical degrees of freedom between the floating platform and rotor are connected by a robust linear spring.As pitch movement needs to be allowed, the term associated with pitch was assigned a value of zero.
The wave excitation load, denoted as L w , is represented as a column vector of values corresponding to each degree of freedom, as shown in Equation (10).However, the values take into account the diffraction effect between them when multiple floating bodies exist simultaneously.
Based on the described equation of motion, a numerical analysis simulation was performed using OrcaFlex version 11.3 [30], a three-dimensional floating multibody dynamics tool that performs simulations in the time domain.

Regular Wave Simulations
In regular wave environments, floating bodies typically exhibit periodic motion dependent on the incident wave period.Moreover, the response of the floating body to the incident wave appears minimal when the amplitude of the incident wave is considerably small.Under these conditions, the motion of the floating body can be considered a linear system response, with the frequency-domain and time-domain responses being nearly identical due to the negligible nonlinear characteristics.Therefore, a regular wave with a wave height of 0.1 m was generated under the water depth condition of the installation site (80 m) to investigate the time-domain response of the floating body.A total of 24 incident wave frequencies, encompassing surge, heave and pitch natural frequencies, were selected within a range of 0.08 rad/s to 2.00 rad/s.This range was chosen to observe the resonance phenomenon for each degree of freedom of the platform and rotors.The analysis in this study was confined to head sea conditions (the direction of the 0 • incident wave), leading to the exclusion of transverse motions such as sway, roll and yaw.
To minimize the initial transient response, the wave amplitude was gradually increased.The average motion amplitude (ξ a ) was then calculated by extracting the data segment in which the floating body's behavior reached a steady state.As shown in Equation (11), the average motion amplitude was divided by the incident wave amplitude (A wave ) to derive the response amplitude operator (RAO) in the time domain of the floating body according to the incident wave frequency.

RAO =
ξ a A wave (11) Figure 6 illustrates the results of regular wave simulations when the platform and rotor exist independently.Specifically, the amplitudes of surge, heave and pitch movements were measured on a floating platform identical to the one depicted in Figure 2b, which was connected to a linear spring, simulating a mooring role.Additionally, the rotor, separated from the floating platform, performed only pitch motion based on a fixed rotational axis underwater, as shown in Figure 2a.The motion response results (OrcaFlex) of the floating platform and rotor at a wave height of 0.1 m were compared with the RAO in the frequency domain (WAMIT).A minor discrepancy was observed in the RAO near the platform's surge natural frequency; however, the two sets of results generally align well.This indicates that the time-domain equation of motion effectively captures the hydrodynamic characteristics of the floating body calculated in the frequency domain.The noted difference in surge RAO peaks is attributed to the divergence of motion due to platform resonance.In the case of the platform, surge, heave and pitch motions are interconnected, resulting in coupling terms in the equation of motion.Consequently, peaks near the natural frequency of each motion mode are also observed in other motion modes.
Figure 6 illustrates the results of regular wave simulations when the platform and rotor exist independently.Specifically, the amplitudes of surge, heave and pitch movements were measured on a floating platform identical to the one depicted in Figure 2b, which was connected to a linear spring, simulating a mooring role.Additionally, the rotor, separated from the floating platform, performed only pitch motion based on a fixed rotational axis underwater, as shown in Figure 2a.The motion response results (OrcaFlex) of the floating platform and rotor at a wave height of 0.1 m were compared with the RAO in the frequency domain (WAMIT).A minor discrepancy was observed in the RAO near the platform's surge natural frequency; however, the two sets of results generally align well.This indicates that the time-domain equation of motion effectively captures the hydrodynamic characteristics of the floating body calculated in the frequency domain.The noted difference in surge RAO peaks is attributed to the divergence of motion due to platform resonance.In the case of the platform, surge, heave and pitch motions are interconnected, resulting in coupling terms in the equation of motion.Consequently, peaks near the natural frequency of each motion mode are also observed in other motion modes.In practical scenarios, the rotor and platform are interconnected, leading to response differences compared to their independent states, as previously mentioned.To investigate changes in response characteristics when the rotor is connected to the platform, additional analyses were conducted for three scenarios, as shown in Figure 7 (detailed in Table 4).These scenarios are as follows: (1) an arrangement of four rotors without a platform (case 1); (2) a scenario in which the rotor is connected to the platform but the platform remains immobile (case 2); and (3) a scenario in which the rotor is connected to a freely floating platform (case 3).In practical scenarios, the rotor and platform are interconnected, leading to response differences compared to their independent states, as previously mentioned.To investigate changes in response characteristics when the rotor is connected to the platform, additional analyses were conducted for three scenarios, as shown in Figure 7 (detailed in Table 4).These scenarios are as follows: (1) an arrangement of four rotors without a platform (case 1); (2) a scenario in which the rotor is connected to the platform but the platform remains immobile (case 2); and (3) a scenario in which the rotor is connected to a freely floating platform (case 3).The pitch RAO of the rotor was determined by calculating the relative rotational displacement between the platform and the rotor because the rotational axes are connected to the platform.As shown in Figure 8 and Table 5, four rotors are arranged on the central platform.The outer rotor has adjacent rotors on only one side, whereas the inner rotor has adjacent rotors on both sides.This configuration leads to slight differences in response due to hydrodynamic interactions.Consequently, it is crucial to differentiate between the outer and inner rotors and assess variations in their responses.Additionally, because the four rotors are symmetrically arranged about the platform's centerline, the responses of inner rotors 1 and 2 are identical.The two outer rotors also show identical responses to each other.The pitch RAO of the rotor was determined by calculating the relative rotational displacement between the platform and the rotor because the rotational axes are connected to the platform.As shown in Figure 8 and Table 5, four rotors are arranged on the central platform.The outer rotor has adjacent rotors on only one side, whereas the inner rotor has adjacent rotors on both sides.This configuration leads to slight differences in response due to hydrodynamic interactions.Consequently, it is crucial to differentiate between the outer and inner rotors and assess variations in their responses.Additionally, because the four rotors are symmetrically arranged about the platform's centerline, the responses of inner rotors 1 and 2 are identical.The two outer rotors also show identical responses to each other.

Multiple Rotor Effect
This study, which assumed the simultaneous operation of four rotors in an ocean environment, initially explored the impact of hydrodynamic interactions among these rotors.Simulations were conducted for both a single rotor and an array of four rotors, and their responses were compared.In both scenarios, the rotor's rotational axis was assumed to be fixed at a specific position, unchanging over time.Figure 9 shows the RAO values for a single rotor and for both the inner and outer rotors in case 1.As previously mentioned, the responses of the inner and outer rotors might differ due to hydrodynamic interactions; however, the results showed that this difference is negligible (if the spacing between rotors becomes closer than 3.5 m (current spacing), the difference due to this effect can become larger).The responses of the single rotor and four rotors exhibited no significant differences.Thus, in this study, it was concluded that the response is not significantly affected by the number of rotors when rotating around a fixed axis.

Multiple Rotor Effect
This study, which assumed the simultaneous operation of four rotors in an ocean environment, initially explored the impact of hydrodynamic interactions among these rotors.Simulations were conducted for both a single rotor and an array of four rotors, and their responses were compared.In both scenarios, the rotor's rotational axis was assumed to be fixed at a specific position, unchanging over time.Figure 9 shows the RAO values for a single rotor and for both the inner and outer rotors in case 1.As previously mentioned, the responses of the inner and outer rotors might differ due to hydrodynamic interactions; however, the results showed that this difference is negligible (if the spacing between rotors becomes closer than 3.5 m (current spacing), the difference due to this effect can become larger).The responses of the single rotor and four rotors exhibited no significant differences.Thus, in this study, it was concluded that the response is not significantly affected by the number of rotors when rotating around a fixed axis.

Fixed Platform Effect
The rotors in this study are connected to a truss platform comprising cylindrical pontoons.Hence, assessing the hydrodynamic effects when the platform structure is adjacent to the rotors is essential.The response characteristics of the rotors were compared for the four-rotor scenarios with and without the platform (Figure 10).In these scenarios, the rotor's rotational axis was assumed to remain fixed at a specific position, allowing for an isolated assessment of the hydrodynamic interaction solely attributed to the platform structure's presence.Minor differences in response were observed primarily near the resonant frequency.

Floating Platform Effect
While the previous section focused on the influence of a fixed platform, in reality, the platform is not stationary but dynamically responds to wave forces.Consequently, the behavior of the rotor connected to a moving (floating) platform was examined (Figure 11).

Fixed Platform Effect
The rotors in this study are connected to a truss platform comprising cylindrical pontoons.Hence, assessing the hydrodynamic effects when the platform structure is adjacent to the rotors is essential.The response characteristics of the rotors were compared for the four-rotor scenarios with and without the platform (Figure 10).In these scenarios, the rotor's rotational axis was assumed to remain fixed at a specific position, allowing for an isolated assessment of the hydrodynamic interaction solely attributed to the platform structure's presence.Minor differences in response were observed primarily near the resonant frequency.

Fixed Platform Effect
The rotors in this study are connected to a truss platform comprising cylindrical pontoons.Hence, assessing the hydrodynamic effects when the platform structure is adjacent to the rotors is essential.The response characteristics of the rotors were compared for the four-rotor scenarios with and without the platform (Figure 10).In these scenarios, the rotor's rotational axis was assumed to remain fixed at a specific position, allowing for an isolated assessment of the hydrodynamic interaction solely attributed to the platform structure's presence.Minor differences in response were observed primarily near the resonant frequency.

Floating Platform Effect
While the previous section focused on the influence of a fixed platform, in reality, the platform is not stationary but dynamically responds to wave forces.Consequently, the behavior of the rotor connected to a moving (floating) platform was examined (Figure 11).

Floating Platform Effect
While the previous section focused on the influence of a fixed platform, in reality, the platform is not stationary but dynamically responds to wave forces.Consequently, the behavior of the rotor connected to a moving (floating) platform was examined (Figure 11).The findings indicate a notable alteration in the rotor's motion response when the platform is freely floating.The peak frequency of the rotor in case 3 shifts marginally toward a lower frequency than that of a single rotor (1.22 rad/s).Furthermore, in case 3, two peaks emerge in the low-frequency region, accompanied by a lower response amplitude compared to case 2. This suggests an interaction effect between the platform's behavior and the rotor's behavior.To verify this, the impact of the platform's surge, heave and pitch motions on the connected rotor was analyzed.Only one degree of freedom of the platform's motion was allowed to determine the platform mode affecting the pitch of the outer rotor.As illustrated in Figure 12, the heave and pitch of the platform significantly affect the rotor's pitch RAO at specific frequencies.These frequencies align with the positive and negative response ranges observed in Figure 11.This is because the connected rotor experiences heave when the platform undergoes heave and pitch motions.Due to the characteristic shape of the rotor, the submerged volume changes rapidly and nonlinearly when it moves up and down.This results in an additional difference in moments caused by gravity and buoyancy.In other words, the heave of the rotor occurs according to the behavior of the platform, and this motion induces additional rotor pitch.Consequently, the pitch response of the rotor connected to the floating platform is the result of a combination of the pitch generated by the wave excitation load (inherent pitch) and the pitch derived from the heave (induced pitch).When considered in terms of the inherent pitch, if the induced pitch aligns with the phase of the inherent pitch, the motion will be greater.Conversely, if the phase is opposite, the motion will be less.This is the reason for the response results shown in Figure 11.This overlap positively or negatively impacts the rotor's motion response, depending on the incident wave frequency.The findings indicate a notable alteration in the rotor's motion response when the platform is freely floating.The peak frequency of the rotor in case 3 shifts marginally toward a lower frequency than that of a single rotor (1.22 rad/s).Furthermore, in case 3, two peaks emerge in the low-frequency region, accompanied by a lower response amplitude compared to case 2. This suggests an interaction effect between the platform's behavior and the rotor's behavior.To verify this, the impact of the platform's surge, heave and pitch motions on the connected rotor was analyzed.Only one degree freedom of platform's motion was allowed to determine the platform mode affecting the pitch of the outer rotor.As illustrated in Figure 12, the heave and pitch of the platform significantly affect the rotor's pitch RAO at specific frequencies.These frequencies align with the positive and negative response ranges observed in Figure 11.This is because the connected rotor experiences heave when the platform undergoes heave and pitch motions.Due to the characteristic shape of the rotor, the submerged volume changes rapidly and nonlinearly when it moves up and down.This results in an additional difference in moments caused by gravity and buoyancy.In other words, the heave of the rotor occurs according to the behavior of the platform, and this motion induces additional rotor pitch.Consequently, the pitch response of the rotor connected to the floating platform is the result of a combination of the pitch generated by the wave excitation load (inherent pitch) and the pitch derived from the heave (induced pitch).When considered in terms of the inherent pitch, if the induced pitch aligns with the phase of the inherent pitch, the motion will be greater.Conversely, if the phase is opposite, the motion will be less.This is the reason for the response results shown in Figure 11.This overlap positively or negatively impacts the rotor's motion response, depending on the incident wave frequency.

Irregular Wave Simulations
The regular wave simulation results underscore that the platform's motion significantly influences the rotor's motion characteristics.Building upon this, irregular wave simulations were conducted to examine the responses in a real sea environment.The incident wave conditions for these simulations were modeled using the operational design wave of the installation site, as detailed in Table 6 and Figure 13.

Irregular Wave Simulations
The regular wave simulation results underscore that the platform's motion significantly influences the rotor's motion characteristics.Building upon this, irregular wave simulations were conducted to examine the responses in a real sea environment.The incident wave conditions for these simulations were modeled using the operational design wave of the installation site, as detailed in Table 6 and Figure 13.To mitigate the initial transient response of the platform, the incident waves were gradually intensified.These conditions were then analyzed over a period of 10,800 s

Irregular Wave Simulations
The regular wave simulation results underscore that the platform's motion significantly influences the rotor's motion characteristics.Building upon this, irregular wave simulations were conducted to examine the responses in a real sea environment.The incident wave conditions for these simulations were modeled using the operational design wave of the installation site, as detailed in Table 6 and Figure 13.To mitigate the initial transient response of the platform, the incident waves were gradually intensified.These conditions were then analyzed over a period of 10,800 s Power spectral density of the incident wave (JONSWAP Spectrum, H s = 2.0 m, T p = 6.65 s, γ = 1.0 . To mitigate the initial transient response of the platform, the incident waves were gradually intensified.These conditions were then analyzed over a period of 10,800 s (equivalent to 3 h).Figure 14 presents a time series of the rotor's angular displacement (relative to the platform [ξ rt−n rel ] when connected to the platform) under operational conditions for each case.As depicted in Figure 15, these time series data were converted into power spectral density to facilitate a reasonable comparison of motion responses.Figure 14b reveals that the behavior of the single rotor and the behavior in case 1 (four-rotor case) are similar.However, the results for case 2 (fixed platform case) exhibit different motion patterns, even under identical incident wave conditions.As shown in Figure 7 and Table 5, the initial positions of the rotors in the x-direction differ between cases 1 and 2 on the global coordinates.In other words, even if the same wave approaches, the timing of encountering the wave varies due to the difference in location.However, as can be seen in Figure 15, the motion spectra of the two scenarios are almost similar.Since this represents the response in the regular wave simulation well (Figure 10), it is judged that the results of the 3 h simulation sufficiently reflect the characteristics of each scenario.The initial positions of the rotors in case 2 and case 3 (floating platform case) are identical, yet the trend in the rotor pitch differs between the two scenarios, influenced by the behavior of the floating platform.
The findings in Figure 15 are similar to the regular wave simulation results.While the response of multiple rotors or rotors connected to a fixed platform did not significantly differ from that of a single rotor, the peak response increased significantly when the rotors were connected to a moving (floating) platform.As observed in Figure 11, the peak frequency of the rotor in case 3 shifted marginally to a lower frequency compared to other scenarios.This shift coincides with the peak frequency area of the incident wave, resulting in rotor resonance.Due to the platform's behavior, the response around 0.6 rad/s an improvement of approximately 29%.This suggests that in our model, the hydrodynamic interaction and coupled motions between the platform and the rotors contribute positively to the overall power generation efficiency.In the integrated floating platform model, the effects of coupled motion, which are not present in a rotor connected to a fixed rotational axis, become evident.Therefore, analyzing the rotor as if it were connected to a fixed axis poses a risk of underestimating or overestimating the extraction power potential of the actual floating model.This underlines the importance of considering the dynamic interactions between the platform and rotors in assessing power generation efficiency in wave energy systems.

Conclusions
In this study, we explored the motion and power generation performance of a Salter's duck-shaped rotor and platform, modeled as a multi-degree-of-freedom system, under wave conditions.The motion responses and power generation performance of both the rotor and the floating platform were evaluated through a dynamic analysis in the time domain.The following conclusions were drawn:

•
Several scenarios were analyzed to understand the hydrodynamic interaction between the floating platform and rotors: a single rotor, an array of four rotors (case 1), rotors connected to a fixed platform (case 2) and a scenario allowing for the movement of the connected platform (case 3).The arrangement of four rotors facilitated a comparison of radiation and diffraction effects.Connecting the rotor to a fixed platform highlighted the hydrodynamic interaction between the platform structure and the rotor.Allowing for platform movement provided insights into the system's performance under operational conditions.

•
Regular wave simulations, conducted at a relatively low wave height of 0.1 m for various incident wave frequencies, aimed to characterize the rotor's behavior, assuming it operated as a linear system.The results from case 1 indicated differences in motion between inner and outer rotors due to their arrangement.Case 2 demonstrated the influence of platform pontoons on rotor behavior.However, the differences observed in these two scenarios were minor.In contrast, the motion response of the rotor changed significantly when connected to a floating platform, with the main peak frequency shifting marginally to the low-frequency region and the emergence of two peak responses in regions outside the natural frequency range.The floating platform's heave and pitch motions substantially impacted the rotor's behavior.

•
Irregular wave simulations under operational conditions were conducted to assess the rotor's performance when connected to a floating platform in real sea conditions.Time series data from a 3 h simulation were converted into power spectra for comparison.The results from the three cases (single rotor, case 1 and case 2) were generally similar.However, case 3 exhibited a higher peak value and a significant response at approximately 0.6 rad/s.The peak frequency of the rotor, approximately 1.22 rad/s, shifted marginally to the low-frequency region when connected to a floating platform.
Considering that the peak frequency of the incident wave (operating condition) was 0.945 rad/s (6.65 s), the rotor operated close to resonance.Additionally, the inherent pitch motion of the rotor caused by wave forcing was combined with the pitch response derived from the heave and pitch of the platform, increasing the response in a specified wave frequency range.

•
Compared to the power extracted by a single rotor under operating conditions, the power generated in case 1 and with rotors connected to a fixed platform increased by approximately 1.86% and 4.45%, respectively.Allowing for the platform's degree of freedom resulted in a power increase of more than 29%.Although the trend in results may vary depending on the hydrodynamic characteristics of the platform and WECs, the case in this study demonstrated a positive influence on the efficiency of power generation.In other words, examining the rotor under the assumption of a fixed axis poses the risk of either underestimating or overestimating the extraction power potential of the actual floating model.This underscores the importance of considering the dynamic interactions between the floating platform and rotors as essential factors when assessing the efficiency of power generation in real systems.

Figure 1 .
Figure 1.Location where the rotor was installed and operated (western sea of Jeju, South Korea).

Figure 1 .
Figure 1.Location where the rotor was installed and operated (western sea of Jeju, South Korea).
(a) Rotor.(b) Side view of platform.

1 (
b) Side view of platform.(c) Top view of platform.

Figure 2 .
Figure 2. Schematic of the rotor and platform.

Figure 2 .
Figure 2. Schematic of the rotor and platform.

Figure 3 .
Figure 3. Schematic of the mooring buoy and cable system.

Figure 3 .
Figure 3. Schematic of the mooring buoy and cable system.

Figure 4 .
Figure 4. Details of the submerged part of the rotor and platform.

Figure 4 .
Figure 4. Details of the submerged part of the rotor and platform.

Figure 5 .
Figure 5. Restoring moment of the rotor as a function of angular displacement.

Figure 5 .
Figure 5. Restoring moment of the rotor as a function of angular displacement.
(a) Surge of the platform.(b) Heave of the platform.(c) Pitch of the platform.(d) Pitch of a single rotor.

Figure 6 .
Figure 6.RAOs of platform and single rotor in regular waves (wave height = 0.1 m).

Figure 6 .
Figure 6.RAOs of platform and single rotor in regular waves (wave height = 0.1 m).

Figure 7 .
Figure 7. Numerical model of rotor and platform in OrcaFlex.

Figure 7 .
Figure 7. Numerical model of rotor and platform in OrcaFlex.

Figure 8 .
Figure 8. Configuration of arrayed rotors on the central platform.

Figure 9 .
Figure 9. Pitch RAO of rotors in regular waves for analyzing the effect of multiple rotors.

Figure 10 .
Figure 10.Pitch RAO of rotors in regular waves for analyzing the effect of a fixed platform.

Figure 9 .
Figure 9. Pitch RAO of rotors in regular waves for analyzing the effect of multiple rotors.

Figure 9 .
Figure 9. Pitch RAO of rotors in regular waves for analyzing the effect of multiple rotors.

Figure 10 .
Figure 10.Pitch RAO of rotors in regular waves for analyzing the effect of a fixed platform.

Figure 10 .
Figure 10.Pitch RAO of rotors in regular waves for analyzing the effect of a fixed platform.

Figure 11 .
Figure 11.Pitch RAO of rotors in regular waves for analyzing the effect of a floating platform.

Figure 11 .
Figure 11.Pitch RAO of rotors in regular waves for analyzing the effect of a floating platform.

Figure 12 .
Figure 12.Pitch RAO of the outer rotor based on the platform's limited degree of freedom conditions.

Figure 12 .
Figure 12.Pitch RAO of the outer rotor based on the platform's limited degree of freedom conditions.

Figure 14 .
Figure 14.Time series of rotor pitch angular displacement.

Table 1 .
Specifications of the single rotor.

Table 3 .
Damping ratio values of the platform and rotor.

Table 3 .
Damping ratio values of the platform and rotor.

Table 4 .
Case classification for multibody analysis.

Table 4 .
Case classification for multibody analysis.

Table 5 .
Initial position for each body coordinate system.

Table 8 .
Average of the absolute value of relative angular velocity for each case's rotors.