Numerical Modeling and 3d Investigation of Inwave Device

In this article, numerical studies on a tightly moored point absorber type wave energy converter called INWAVE are presented. This system consists of a buoy, subsea pulleys, and a power take off (PTO) module. The buoy is moored by three ropes that pass through the subsea pulleys to the PTO module. Owing to the counterweight in the PTO module, a constant tension, which provides a horizontal restoring force to the buoy, is constantly applied to the rope. As waves pass by, the buoy is subjected to six degrees of freedom motion, consisting of surge, heave, sway, roll, pitch, and yaw, which causes reciprocating motion in the three mooring ropes. The PTO module converts the motion of the ropes into electric power. This process is expressed as a dynamic equation based on Newtonian mechanics and the performance of the device is analyzed using time domain simulation. We introduce the concept of virtual torsion spring in order to prevent the impact error in the ratchet gear modules which convert bidirectional motion of rope drum into unidirectional rotary motion. The three-dimensional geometrical relationship between the ropes and the buoy is investigated, and the effects of the angle of the mooring rope and the direction of wave propagation are addressed to determine the interaction between the tension of the rope and the buoy. Results have shown that the mooring rope angle has a large impact on the power extraction. The simulation results present a useful starting point for future experimental work.


Introduction
Ocean wave energy is a renewable energy source with high energy density and great potential.Ocean waves have higher energy density as they are deeper offshore.The water particles here draw a circular orbital trajectory.Conversely, and the lower the depth of water, the closer the ocean is, the lower the energy density of the waves, and the water particles have the orbitals of long transverse orbits.That is, as the water depth decreases, the energy distribution in the horizontal direction becomes larger.In the shallow sea, furthermore, the irregularity of the waves becomes larger depending on the terrain of the sea floor [1,2].
Depending on the characteristics of the waves due to these water depths, the power generation mechanism from waves varies greatly depending on the place where the wave energy converter (WEC) is installed [3,4].It is advantageous to use the heave motion in the deep sea, and point absorber types such as OPT PowerBuoy [5], Wavebob [6], and CETO [7] have been studied to maximize vertical energy absorption.The attenuator-types WEC resembling the inverted oscillator such as Oyster [8,9] and WaveRoller [10] are advantageous for absorbing energy mainly in the surge direction.It is therefore suitable for areas with lower depths than the target area of the point absorber.The coastal installation of WEC is led by the fixed-structure oscillating water column (OWC) types that were built in Norway (in Toftestallen, near Bergen, 1985, [3]), Japan (Sakata, 1990, [11,12]), India (Vizhinjam, near Trivandrum, Kerala state, 1990, [13]), Portugal (Pico, Azores, 1999, [14]), UK (the LIMPET plant in Islay Island, Scotland, 2000, [15]), Spain (breakwater OWC, Mutriku, Northern Spain, 2008, [16]), and Italy (U-OWC, Civitavecchia, 2012 [17,18]).Recently, moving from the experience of the offshore over-topping devices (OTD) Wavedragon [19] and the Sea-wave Slot-cone Generator [20], an onshore OTD has been built in Italy [21].Fixed-structure OWCs and OTDs have good accessibility and maintainability, but there are topographical constraints where the wave energy density should be high and the civil construction should be possible.
The distinction by such depth or distance from the land determines the motion mechanism.In addition, the more distant the sea, the higher the cost of the facilities required for electric transmission to areas where electricity is mainly consumed.There is another alternative to electric transmission.The CETO and Oyster WECs pump the seawater to the land that drives the generator, so the sea water pipe replaces the electric transmission cable.In this case, it is very advantageous for maintenance of the generator since the generator is located on the land.That is, the way to transfer energy to the land is also a main factor in determining the mechanism of wave power generation.
Meanwhile, one of the areas where the introduction of wave power is difficult is the area where the water depth is shallow and wave irregularity is large.In this area, it is difficult to design the mechanism of conversion process from wave to mechanical energy; however, accessibility on land is comparatively good, and it is distributed all over the world.There is a need for a WEC capable of targeting these areas.It is expected that such areas will need a strategy of securing the desired capacity by installing multiple WECs with small capacity rather than a single large capacity WEC, and the cost of a single device should be cheap.
The off-grid islands have a great need for the development of these wave devices, where there is no power supply from the mainland.According to reports from the Korea Electric Power Corporation (KEPCO), electric power costs in these areas are eight times higher than those on the land [22].In many of these areas, the required amount of electric power is often small, such as several tens of kilowatts.Therefore, a small capacity wave generator is more suitable for these areas.
The wave power company INGINE has created a wave power generator called INWAVE, which is a wave generator suitable for coasts with shallow water depths and wave irregularities [23].This is a multi-tight moored moving object, sharing similar concepts with Fred Olsen's Lifesaver WEC [24].However, unlike Lifesaver, INWAVE has a reinforced mechanism to cope with the wave force in a horizontal direction as the installation depth is low.It is also advantageous for maintenance because the PTO module is located on the coast.
It mainly consists of a buoy, subsea pulleys, and a power take off (PTO) module as shown in Figure 1.Three mooring ropes are attached to three points on the bottom of the buoy and also connected to three counterweights in the PTO module via subsea pulleys on the sea-bed.The counterweight moves only vertically and this provides restoring force to the buoy, and the buoy performs six degrees of freedom movement of heave, surge, sway, roll, pitch, and yaw within a certain range.In addition, the tension prevents the rope from sagging and serves to transfer the kinetic energy of the buoy to the PTO module.As the buoy moves, the three ropes reciprocate and also the rope pulleys in the PTO module that rotate back and forth share one shaft and transmit unidirectional rotary motion to the common shaft through their ratchet gears.This rotational energy is then converted into electrical energy by the electric generator.When the buoy is pulling the rope, the tension of the rope raises the counterweight while transmitting power to the common shaft through the ratchet gear.On the other hand, when the buoy is no longer pulling the rope, the tension is provided by the counterweight, which pulls the buoy back to its equilibrium position.However, during this time, the ratchet gear cuts off the work of the counterweight, so there is no power transmitted to the common shaft.The two-dimensional modeling and simulation of INWAVE has been presented in a previous paper [25].This modeling assumes that only two ropes are connected between the buoy and the PTO module unlike the actual plant, and only surge, heave, and pitch motions of the buoy were considered while sway, roll, and yaw motions were ignored as a limitation of the two-dimensional analysis.In this paper, a three-dimensional model of the INWAVE device is presented, which considers the six degree motion of the buoy and the three mooring rope's connections in order to overcome previous assumptions.To build this dynamic model, the geometric relationship between the behavior of the buoy and the ropes in the three-dimensional space was identified, and, based on this, we derived the overall dynamics by combining the hydrodynamic coefficients for the buoy using an ANSYS AQWA simulator [26] and the PTO module based on the Newtonian mechanics and linear wave theory (LWT) [27,28].Next, we carried out simulations under regular wave conditions with varying wave periods to analyze the effects of the mooring rope angle, counterweight's mass, and wave incidence angle, which are the most important factors that influence the system output.Then, for objective efficiency comparison with other WECs, we analyzed the theoretical performance of the device by calculating the capture width ratio (CWR) that is the average power output of the device divided by wave energy density and the diameter of the buoy under irregular wave conditions based on the joint north sea wave project (JONSWAP) wave spectrum model.

Geometrical Relationship between Buoy and Ropes in Three-Dimensional Space
Consider a buoy in the sea with a depth of h, where the buoy radius and draft are R and d, respectively.The behavior of the buoy is given by X = [x, y, z, q x , q y , q z ] T corresponding to surge (x), sway (y), heave (z), roll (rx), pitch (ry), and yaw (rz), respectively, since the wave propagates in the x-direction.The top view of the buoy, subsea pulleys, and rope connection points are shown in Figure 2.An actual device is installed near the breakwater.Unfortunately, ANSYS AQWA, which is a computational fluid dynamics simulator to obtain the hydrodynamics, assumes a flat undersea surface and does not provide an effect on reflected waves.Hence, it is assumed that the buoy is installed on open sea.At the bottom of the buoy, the ropes are connected at a distance of r b , and the connection points make an angle of 120 • with each other.At the bottom surface of the sea, the three subsea pulleys are each located at a distance of r p from the center and are arranged at an angle of 120 • with respect to each other.
In the equilibrium state, since the forces exerted on the buoy are only the tensions of the ropes exerted by the counterweight and buoyancy in the heave direction, the center of the subsea pulley and the buoy are placed on the z-line, and each pulley and rope connection point are placed on the same line from the center.Since the buoy is z-axis symmetry, the direction of wave propagation only affects the the rope connection points.Hence it is assumed that the wave direction is heading to the x-direction as seen in Figure 2 and the relative angle between the wave direction and the rear pulley is defined as q w , which has a value of -60 • to 60 • , but is limited to a range of 0 • to 60 • by symmetry.
The dynamic model of the device includes the influence of the rope tensions; hence, the geometric analysis between the rope connection point and subsea pulley is necessary.First, the positions of rope connection points and the subsea pulleys must be defined in three dimensions.For convenience, let the position of ith rope connection point and subsea pulley point be X r,i and X p,i , respectively.These points can be derived using the homogeneous transformation matrices P r,i and P p,i as follows [29]: P r,i = Trans (x, y, z) Rot x (q x ) Rot y q y Rot z (q z ) Trans (r b cos (q w + q i ) , r b sin (q w + q i ) , 0) , (1) P p,i = Trans (r b cos (q w + q i ) , r b sin (q w + q i ) , 0) , where Trans, Rot x , Rot y , and Rot z are defined as cos q y 0 sin q y 0 0 1 0 0 − sin q y 0 cos q x 0 0 0 0 1 X r,i and X p,i are the last column vectors of P r,i and P p,i that can be expressed as follows: x + r b cos(q y ) {cos(q i + q w + q z ) − sin(q i + q w + q z )} y + r b cos(q i + q w + q z ) cos(q x ) + sin(q x ) sin(q y ) +sin(q i + q w + q z ) cos(q x ) − sin(q x ) sin(q y ) z + r b cos(q i + q w + q z ) sin(q x ) − cos(q x ) sin(q y ) + sin(q i + q w + q z ) sin(q x ) + cos(q x ) sin(q y ) where q i = 0, 2π/3, −2π/3 (for i = 1, 2, 3).From this relationship, the ith rope vector L i ∈ R 3×1 is given by The ith length of the rope l i , the speed of the rope dl i /dt (or li ), and the acceleration of the rope d 2 l i /dt 2 (or li ) are expressed as follows: Let Ẍr,i = A r,i (X, q i , q w ) Ẍ + B r ( Ẋ, X, q i , q w ); then, Equation (12) becomes where

Dynamics of the Buoy and Rope Tension
The equations of buoy motion in surge, sway, heave, roll, pitch, and yaw directions can be expressed in Newtonian mechanics based on the LWT [27,28].These six equations representing each axis are influenced by the excitation force due to waves, radiation damping caused by buoy motion, buoyancy, and rope tensions.The six equations are expressed as one vector equation that can be expressed as [25,27]: Here, M b ∈ R 6×6 is the mass matrix consisting of buoy mass, buoy inertia, and added mass, F r ∈ R 6×1 is the radiation damping force vector, and K h ∈ R 6×6 is the buoyancy spring coefficient matrix.H i ∈ R 6×1 is the ith matrix of the rope tension effect on the buoy and f i is the tension applied to the ith rope.F e ∈ R 6×1 is the excitation force vector.H i can be expressed as where ⊗ represents a cross product.M b and K h are given as [25] where m b , J x , J y , and J z represent the buoy's masses and inertia of the x-, y-, and z-axes, respectively.Furthermore, , and µ rz (∞) are the added masses and inertia at infinity for the buoy's surge, sway, heave, roll, pitch, and yaw motions, respectively.
ρ and g denote the sea water density and the gravitational acceleration, respectively.diag[] is a diagonal matrix.
As described above, the wave direction is limited to the x-direction and it yields the excitation forces of surge, heave, and pitch, only acting on the buoy due to the buoy's symmetry, and other components are omitted.The excitation force vector F e can be computed as a convolution of the water surface elevation ψ(t) and impulse response function from the frequency domain results.τ denotes the time shift for convolution: where f e,x , f e,z , and f e,ry denote the excitation kernel functions for surge, heave, and pitch, respectively, which can be obtained by the following formula: where fe,θ (ω) is the excitation force in the frequency domain and φ e,θ is the phase shift of the excitation force, which is π/2, 0, and π/2 for surge, heave, and pitch, respectively.This phase shift is necessary for the same reason that water particles by the waves cause the orbital motion [1].
Similarly, the radiation damping force, which is a resistance component, can be obtained by convolution of the buoy's velocity and radiation kernel function as follows [30]: where Here, fr,θ = Rθ (ω and Rθ (ω) and µ θ (ω) denote radiation damping and added mass in the frequency domain, respectively.It is assumed that there is no resistance in the yaw direction because this part is disk-shaped.

PTO Modeling
Figure 3 is a simplified representation of the configuration of the PTO module.In the PTO module, three mooring ropes hang on rope pulleys and three counterweights are attached to the end of the ropes to provide constant tension equal to their weight.The rope pulley has a ratchet gear module that only transmits torque to the common shaft when the rope is pulled by the buoy.The rotation of this shaft is increased by the gear ratio κ and transmitted to the generator.The angle and angular velocity of the ith rope pulley are Q i and Ω i , respectively, while that of the common shaft and generator are Q p and Ω p , and Q g and Ω g , respectively.Then, these variables have the following relationship: where r i is the ith pulley radius.Assuming that there is no inertia of the shaft and pulleys, the dynamics of the generator are given as follows: where J g is the generator inertia, τ i is the torque transmitted by the ith pulley, and τ g is the induced torque applied by the generator and is given as where b g is the generator damping coefficient that can be controlled using an inverter at the rear end of the generator.The instantaneous output from the generator is defined as , The ratchet gears continuously supply and block the torque from the rope pulley to the common shaft.When the torque is transmitted, the rope pulley and the common shaft are coupled and driven as one part, and when the torque is not transmitted, the two parts move independently.At the moment when the rope pulley and the common shaft are combined, the analysis is linked to their physical properties such as the strength and elastic modulus and is difficult to model using Newtonian mechanics.Therefore, to solve this problem, we introduced a virtual spring concept in the ratchet module in our previous paper [25].This spring occurs at the moment of contact between the pulley and the common shaft, momentarily stores the compressed force as potential energy, and gradually transfers this energy to the shaft without power loss.Then, as soon as negative torque is transmitted through the spring, the pulley and the common shaft are separated.This process can be expressed as follows: where k r is the spring coefficient of the virtual spring, Then, when τ i < 0, the contact between the ratchet and the shaft is released.The dynamic equation of the counterweight, including the transmitted torque to the shaft, is as follows: The power transfer efficiency of the ratchet gear to the generator during one wave peak period ηratchet is defined as follows: Substituting Equations ( 13) and (30) into Equation ( 16) yields the combined dynamics of the buoy-counterweights as Combining Equation (32) with Equations ( 26) and ( 29) completes the final dynamic system in which the wave elevation ψ(t) and constant b g in the time domain are the input states and the motions of the buoy [X, Ẋ] and the PTO [Ω 1 , Ω 2 , Ω 3 , Ω p ] are the output states.

Simulation and Analysis under Regular Wave Condition
Currently, a prototype installed on Jeju Island has a buoy diameter of 5 m and the maximum power is obtained in a period between 2.5 and 4.5 s [25].However, the wave period of the area under consideration for installation is 7 to 12 s.By the simulation of buoys with various geometries, it was found that a diameter of 12 m is suitable for the target area.Therefore, in this paper, a simulation is performed based on a buoy with a diameter of 2R = 12 m and a draft of d = 1 m.
In order to perform simulation on the WEC system, we need the hydrodynamic data such as excitation force, radiation damping, and added mass in the frequency domain.The ANSYS AQWA simulator was used to obtain the data [26].A disk-shaped buoy with a diameter of 2R = 12 m and a draft of d = 1 m was applied floating on a sea with a depth (h) of 10 m; the acquired data are shown in Figure 4.This simulator is based on LWT and assumes incompressible, irrotational flow with small amplitude motion [27,28].Therefore, the simulation only considers the case of 1 m wave height (H s ).
By substituting this data into Equations ( 21) and (23), we obtain the kernel functions, and then substituting the wave elevation ψ(t) and the excitation kernel function f e,θ into Equation (20) yields the excitation force in the time domain.As mentioned in the previous section, the overall dynamic equation consisting of Equation (32) with Equations ( 26) and ( 29) is analyzed in the time domain using the fourth-order Runge-Kutta method [31].
In the modeling of INWAVE, the influence of the rope is included in H i .Since this matrix has nonlinearities including the trigonometric function, the analysis of the overall modeling is impossible in the frequency domain.Therefore, the simulation under the regular wave condition in the time domain was carried out to analyze the exact dynamic behavior of the device with the characteristics of the power output.The following sine wave was generated for regular wave condition: where H s and w p are the significant wave height and peak wave frequency, respectively.'"( ) !"#$%& ) '.( The output power P g of the WEC is sensitive to generator damping load b g .We substitute b g between 10 3 and 10 7 Nms/rad, and denote b g,opt when the average of the maximum output power from generator Pg,opt is obtained.The data for all parameters required to perform the simulation are listed in Table 1.In order to observe the characteristics of each wave frequency, H s was fixed at 1 m and ω p was assigned at intervals of 0.1 rad/s between 0.2 and 1.4 rad/s.The total sampling time was 50 × T p where T p = 2π/ω p , and only the results within 10 × T p to 50 × T p were recorded in order to exclude the data in the transient segment.Except for ηratchet defined in Equation (31), the mean value (•) and standard deviation (σ(•)) used the data in this interval.The INWAVE device has all the states intertwined because the buoy moves in six degrees of freedom and three ropes are connected to the single generator.Therefore, it is hard to obtain the perfect periodicity of the states according to the wave period even in regular wave conditions.In the case of r p = 10 m, m c = 1000 kg, b g = 2 × 10 5 Nm•s/rad, q w = 30 • , and T p = 12.56 s (ω p = 0.5 rad/s) as presented in Table 1, while the simulation results are shown in Figure 5. Figure 5b shows that the surge, heave, and pitch motions of the buoy have periodicity with the wave period and out of these parameters, the surge has the highest amplitude.In Figure 5c-e, however, the angular speeds of the rope pulley and generator, rope tensions, and power output do not have repeated patterns.In Figure 5c, the angular velocity of the generator follows the highest value of the angular velocity of the rope drum, but not perfectly.Even then, the average efficiency of the ratchet in Figure 5f is approximately 1. Accordingly, it can be estimated that there is no power loss in the ratchet transmission.The rope tension in Figure 5d only increases when the power is transferred; otherwise, it acts on the rope only as much as the counterweight's mass.
As can be seen in Figure 5, the device has a strong nonlinearity and cannot be directly analyzed in the frequency domain.Therefore, we resolved to understand the behavior and performance of the device by taking the average data for a long period of time in the time domain while varying several parameters such as the pulley position r p , mass of the counterweight m c , and the wave direction q w .

Influence of Subsea Pulley Position
The simulation was carried out where the rope connection points at the bottom of the buoy were fixed at r b = 5 m and the position of the subsea pulley r p was selected from 5, 7, 10, 15, and 22 m.The mooring rope's angle corresponding to r p with the sea floor is given in Table 2 and Figure 6 shows the rendered image of the buoy and ropes with α according to r p .In the case of r p = 5 m, the rope tension acts on the buoy only in the vertical direction, and as r p increased, the horizontal component of the rope tension acting on the buoy also increased.Figure 7 shows simulation results of the average maximum power output Pg,opt and the corresponding damping load b g,opt at various r p values that can be obtained at each wave period.In the graph, the case of r p = 5 m showed the highest average power and the power decreased as r p increased.In addition, it can be seen that the wave period of the peak power output slightly decreased while the r p increased.The smaller the α, the larger the lateral force tension generated, which affects the resonance period.α also influences how efficiently the movement of the buoy is transmitted to the displacement of the rope.As the rope connection point moved in the tangential direction around the subsea pulley, the change in the length of the rope decreased, so that the transmission rate of the kinetic energy was reduced.Therefore, we define a rope kinematic efficiency factor ηrope in order to determine the effectiveness of transmitting the buoy's motion to the rope's displacement, which can be defined by the ratio of the absolute value of the rope connection point speed to the speed of the rope length.This can be expressed as follows: As the efficiency ηrope,i was very close to 1, the movement of the rope connection point was effectively transmitted to the change in the length of the rope.For a more accurate analysis, we compared the power output, behavior of the buoy, and displacement of the rope for variations in r p values between 5 and 22 m as seen in Table 2 at T p = 4.49, 8.98, and 12.57 s; the results are shown in Figure 8. Figure 8a shows the optimal mean power and Figure 8b-d show the normalized standard deviation σ * (•) in the direction of surge, heave and pitch, which denotes the standard deviation of the displacement divided by the amplitude of the wave.Figure 8e,f show σ * (•) of ropes 1 and 2, respectively.Since q w = 0 • , the motion of rope 3 is the same as that of rope 2, and it is omitted.Figure 8g,h are the motion transfer efficiencies of the rope.First, let us consider the results for T p = 7.85 and 15.71 s in order to analyze the concentration of power output in 7-15 s of the wave period.In Figure 8a, the power output tends to decrease as r p increased and the result is similar to σ * (x) In Figure 8b, it can be seen that σ * (x) is overwhelmingly higher than σ * (z) and σ * (ry) in Figure 8c,d, respectively.This means that the output of the device depends mostly on the surge movement.It seems to be a good choice to choose r p = 5 m in order to maximize the motion in the surge direction, but the opposite tendency can be seen from Figure 8e-h.In Figure 8e,f, when T p = 15.71s, it can be seen that σ * (rope) increased with increasing r p until r p = 15 m because the increasing rate of ηrope,i is higher than the decreasing rate of σ * (x).The dynamic modeling so far considered did not take into account resistances other than radiation damping such as friction and viscosity.If all resistive components were considered, a large σ * (x) as shown in Figure 8b cannot be expected.Therefore, the actual model should select r p considering both the effects of surge displacement and the motion transfer efficiency of the rope.In addition, the case of T p = 4.49 s is different from T p = 15.71 and 7.85 s.It should be noted that Figure 8g,h, r p showed a very high ηrope,i near 5 m.This is because the pitching motion increased as T p decreased, which affects the normal direction movement of the rope rather than the surge and heave.

Influence of Counterweight
In order to evaluate the effect of the counterweight, we carried out a simulation for the case when r p = 10 m was fixed and m c was set to 500, 1000, 2000, and 4000 kg.The other parameters remained the same as shown in Table 1. Figure 9 shows the maximum power output Pg,opt and the corresponding b g,opt .As m c increased, Pg,opt generally increased, but the extent of the increase gradually decreases.The increase in the power output by the counterweight was prominent in the long period, and the effect was small in the short period.In addition, as m c increased, the peak of the output shifts slightly to the smaller period, but this tendency is not noticeable.

Influence of Wave Direction
In order to analyze the influence of the direction of the wave, the average power output and the standard deviation were calculated for various wave directions such as q w = 10 • , 20 • , 30 • , 40 • , 50 • , when r p = 10 m, b g = 10 5 Nm•s/rad.Figure 10 shows the results of the power output with varying wave direction for T p = 15.71,7.85, and 4.49 s, where the thick green line is the average output and the red line is the standard deviation.When T p = 15.71s, the average power and standard deviation decreased slightly at approximately 20-30 • , although this is not significant.When T p = 7.85 s, the average power and standard deviation increased by approximately 30% at q w = 30 • .In the case of T p = 4.49 s, the average output and standard deviation decrease by approximately 20% as the angle increased.Synthetically, consistent trends could not be derived, so we could not identify meaningful effects of the wave direction.

Simulation and Analysis under JONSWAP Irregular Wave Condition
In this section, let us compare the optimal average power for an irregular wave condition with the theoretical wave energy density.We chose the JONSWAP wave model because the target area where the INWAVE device has been installed has γ = 2.72 of the JONSWAP spectrum [32,33].When the significant wave height is H s and the peak period is T p (= 2π/ω p ), it is expressed in the frequency domain with wave angular frequency ω as S J (ω) = where The wave power per unit crest length is defined as where and Here, the wave number k = 2π/λ (λ is the wave length) is obtained from When discussing the performance of many wave power systems, the availability and presentation of data vary greatly between sources.Therefore, many devices use a capture width ratio (CWR) to objectively compare the conversion efficiency.For a cylindrical buoy, the CWR is the average power output divided by J J and the diameter of the buoy and is expressed as [34]: For time series calculations, the spectral distribution in Equation ( 35) is discretized as the sum of a large number N of regular waves and written as Here, ω n = ω l + (n − 1)∆ω, where ω l is the lowest frequency, ∆ω is a small frequency interval, n = 1, 2, ..., N and the spectrum does not contain a significant amount of energy outside the frequency range ω l ≤ ω ≤ ω l + (N − 1)∆ω.ψ n = 2S J (ω n )∆ω and θ n are the amplitude of the wave component of order n and the initial phase randomly chosen in the interval (0, 2π), respectively.This wave elevation ψ(t) is substituted into Equation (20), and the simulation was carried out with the equivalent parameters in Table 1.
The average power output under irregular waves in Figure 11a was lower than that of the regular wave, but its pattern is similar to that of the regular wave.When r p = 5 m, the highest power output was obtained and the output power gradually decreased as r p increased.The CWR shown in Figure 11 exhibited a similar pattern, having a peak of 0.43 at r p = 5 m, and the peak gradually disappeared as the r p increased.In addition, it can be seen that the CWR was much lower at a wider spectrum than in our previous paper [25].This is a characteristic of the general point observer that can be seen as the diameter of the buoy increases.Figure 12 shows the standard deviations of the buoy's behavior in each direction, the standard deviation of the rope motion, and the motion transfer efficiency of the rope, which were also similar to that of the regular waves shown in Figure 8.At T p = 15.71 and 7.85 s, the surge showed the greatest movement, and it decreased rapidly as r p decreased.In addition, the efficiency of the rope increased as r p increased, similar to the results of the regular wave condition.

Conclusions
In this study, we derived the three-dimensional dynamics of the INWAVE device and carried out simulations in time domain under regular and irregular wave conditions while varying various parameters such as the subsea pulley position, the counterweight's mass, and the wave direction.The analysis of the results show that the distance between the subsea pulleys significantly affects the behavior of the buoy, the displacement of the rope, and the shape of the power output.The smaller the value of r p , that is, as the mooring angle α approached 90 • , the dynamic behavior of the buoy became higher; however, the motion transfer efficiency of the rope decreased.This is because a large mooring angle leads to higher movements of the rope in the tangential direction, which decreases the effect on the rope length change.Since the dynamic behavior of the actual buoy motion would be greatly reduced by the unconsidered resistance factors, it may be better to consider the motion transfer efficiency of the rope more than increasing the dynamic behavior of the buoy.It is expected that the actual rope mooring angle α would be between 45 • and 60 • .
The effect of the counterweight increased the power output as the mass increased, but the extent of the increase in the output gradually decreased.Therefore, it was concluded that this factor does not affect the output as much as the subsea pulley position.The mass of the counterweight serves to restore the buoy to its origin.Therefore, it is necessary to cope with the drift and to confine the region of the buoy motion within a certain boundary even in extreme wave conditions.In addition, since the basic of the rope is determined according to the weight of the counterweight, which also determines the specifications of mechanical components such as the rope thickness and the rope drum size, a suitable value of the counterweight cannot be selected simply by setting the performance of the power output.
In Korea, due to the monsoon effect between Eurasia and the Pacific Ocean, the Northwest and the Southwest monsoons are dominant in winter and summer, respectively.Waves also have different seasonal direction due to these winds.Thus, it is very important to analyze the influence of the direction of propagation of the wave.However, the simulation results did not yield a meaningful conclusion.This implies that when choosing the position of the subsea pulley, the influence of the traveling angle of the waves can be given less consideration.
Currently, there are three full-scale prototype devices installed on Jeju Island, Korea as a test run [23].The diameter and draft of the buoy are 5 and 0.5 m, respectively.Three mooring ropes have an angle of 60 • with the sea floor.The average sea depth is 3 m.At the end of the each rope, a 250 kg counterweight is attached.A speed increaser module in the PTO module increases the angular velocity of the rope drum by 35 times and transfers it to the 20 kW AC generator.The generated power from the generator is converted into stable AC and is transmitted to the grid through the AC/DC converter and the DC/AC inverter to be transmitted.
The most important element of this device is the rope.The rope is constantly weighed by the counterweight, and it continually contacts with the subsea pulley and the rope.Therefore, it is important to secure durability against salt water and friction.As seen in the simulation results, the maximum tension on the rope is determined by the load on the generator rather than the counterweight.Therefore, the determination of the load of the generator has a significant influence on the maximum power generation and the durability of the rope.There are currently ongoing optimization processes to ensure mechanical durability and high efficiency by testing various rope materials and various control logic in the equivalent three prototypes.The optimization process will be discussed later in the paper.
We expect that, in the near future, it will be possible to compare and verify theoretically obtained characteristics using actual data.We will also present a simulation that reflects the actual characteristics of the machine components, such as backlash and friction in the gear modules, and nonlinearity in the generator.

Figure 1 .
Figure 1.Schematics of the INWAVE device model configuration.

Figure 2 .
Figure 2. Top view of the buoy and ropes.

Figure 5 .
Figure 5. Simulation results of the INWAVE device under regular wave condition with H s = 1 m and T p = 12.56 s when b g = 2 × 10 5 Nms/rad and r p = 10 m: (a) wave elevation; (b) buoy's behavior; (c) pulley's angular speed; (d) rope tension; (e) power from generator; and (f) ratchet efficiency.

Figure 6 .
Figure 6.Rendered image of buoy and ropes for: (a) r p = 5 m; (b) r p = 10 m.

Figure 7 .
Figure 7. Graphs of: (a) Pg,opt and (b) b g,opt with varying r p for regular waves with H s = 1 m.

Figure 9 .
Figure 9. Graphs of: (a) Pg,opt and (b) b g,opt with varying counterweight mass m c for regular waves with H s = 1 m.

Figure 10 .
Figure 10.Graphs of mean and standard deviation of P g with varying wave direction q w at: (a) T p = 15.71s; (b) T p = 7.85 s; (c) T p = 4.49 s (thick green: average, red: standard deviation).

Figure 11 .
Figure 11.Graphs of: (a) Pg,opt ; (b) CWR; and (c) b g,opt with varying r p under JONSWAP irregular waves condition with H s = 1 m.

Table 1 .
Parameters and values in the simulations.

Table 2 .
The pulley position and mooring angle.