Energy Efficiency Analysis of Multi-Type Floating Bodies for a Novel Heaving Point Absorber with Application to Low-Power Unmanned Ocean Device

Long-term energy supplies hinder the application of the low-power unmanned ocean devices to the deep sea. Ocean wave energy is a renewable resource with amount stores of enormous and high density. The wave energy converter (WEC) could be miniaturized so that it can be integrated into the devices to make up the power module. In this paper, a small novel heaving point absorber of energy supply for low-power unmanned ocean devices is developed based on the counter-rotating self-adaptive mechanism. The floating body as an important part of the heaving point absorber, the geometric parameters is optimized to increase the efficiency of power production. Through constructing the constitutive relation between the geometric parameters, the wave force, the motion displacement, the motion velocity, and the capture width ratio of the floating body, the energy efficiency characteristics of the multi-type floating bodies are calculated, and the optimal shape is selected. On the other hand, in the calculation process of the wave force, the Froude-Krylov method is an effective method to accurately calculate the wave excitation force. Meanwhile, nonlinear static and dynamic Froude-Krylov force effectively overcomes the inaccuracy of the linear models and reduces the time consumed to simulate. Finally, the wave force, heaving velocity, heaving displacement, and capture width ratio of the three floating bodies are compared and analyzed, and the results show that the cylindrical floater that is vertically placed on the wave surface is more suitable for the novel heaving wave energy point absorber.


Introduction
Low-power unmanned ocean device, such as unmanned surface/under water vehicles, ocean robots, and ocean buoys are widely used in unmanned combat, deep-sea exploration, and marine communications, and so on [1][2][3][4][5].With the continuous development of deep sea and offshore strategies, the scope of work of low-power unmanned ocean devices continues to expand [6].Long-term, stable, and reliable energy supply is the basis for the rapid development of low-power unmanned ocean device [7].Currently, the common energy supply technology is a battery.However, its energy density is low and it is necessary to carry multiple batteries to meet the requirements.The device space is limited, which limits the amount of battery that is carried [8].Therefore, it is difficult for the battery to meet the requirements of the long-range and long-life capability of the current low-power unmanned ocean device.Therefore, studying the energy supply technology of high energy density and improving Energies 2018, 11, 3282 3 of 20 previous works show that the maximum blade deflection angles between 35 • and 55 • are the best efficiency characteristics [7].In order to increase the power production of the novel heaving point absorber, it is necessary to increase the elevation of the Underwater PTO in the vertical direction, that is, increase the heave amplitude of the floating body.Because the displacement of the body is performed by its immersion depth, shape, and dimension, so the geometric parameters of the body are optimized, which is of great significance for increasing the power production of the novel heaving point absorber.
The floating body dynamics of the WEC have been investigated by many researchers.For example, Black et al. [27] comparatively analyzed the wave radiation forces and scattering forces for horizontal rectangular and vertical circular cylinders using the Haskind's theorem.Mohapatra and Guedes Soares [28] studied the wave forces on a two-dimensional rectangular floating structure based on linearized Boussinesq.Rodriguez et al. [29] investigated the numerical nonlinear heave response of a rectangular box concerning the importance of the relative body dimensions.Islam et al. [30] analyzed the wave radiation of a heaving box-type floating structure based on CFD simulations with a volume of fluid method.Yeung et al. [31] and Sabuncu et al. [32] discussed the added mass and damping of a vertical cylinder in finite depth water.Calisal et al. [33] presented an efficient method of hydrodynamic coefficients calculation for vertical composite cylinders at finite depth.Mansour et al. [34] analyzed the diffraction of linear waves by a uniform vertical cylinder with cosine-type radial perturbations.Bhata et al. [35] studied scattering and radiation problems for a cylinder on nonlinear wave loading at finite depth.Kim [36] researched the hydrodynamic coefficients of the floater with elliptical cylinder and ellipsoid on a free surface.Bihs et al. [37] simulated a horizontal cylinder in heave motion and the motion of a freely floating rectangular barge in waves using the CFD model and compared the results with experimental data.Koh et al. [38] used Matched Eigenfunction Expansion Method (MEEM) to solving the radiation problem of the heaving circular cylinder in the context of linear potential theory.Wang sheng [39] discussed the added mass and damping of an ellipsoid in infinite and finite depth water.Finnegan et al. [40] determined an analytical approximation for the wave excitation forces on a floating truncated vertical cylinder in water of infinite depth and solved the appropriate boundary value using the method of separation of variables.Ghadimi et al. [41] presented a detailed analytical solution for the boundary value problem to evaluate the wave loads for the cylinder with heave and pitch motions in water of finite depth in the presence of an incident wave.Although the floating body dynamics of the WEC are more researched, the energy efficiency of multi-type floating bodies in the heaving wave energy point absorber is poorly studied.Furthermore, in the above paper, the separation variable method and the eigenfunction expansion method is usually used for the calculation of the radiation force and diffraction force of the buoy base on the potential flow theory.However, it is difficult to solve their analytical solutions and is often time consuming to calculate.
Under the Airy's wave theory, the Froude-Krylov approximation method is implemented to solve the wave force when the device dimension is considerably smaller than the wave length [42].Froude-Krylov approximation method [33][34][35][36][37][38][39][40][41][42][43][44][45] is assumed that the original wave pressure distribution of the incident wave does not change due to the presence of the floating body.Therefore, wave excitation force is the product of the force of the undisturbed incident wave pressure on the floating body (dynamic Froude-Krylov force) and the diffraction correction coefficient.The diffraction coefficients reflect the attached mass effect and diffraction effects, which is determined by the model test.Although this method is an approximation method, it is simple to calculate and it is a very practical method for estimating the wave force.Moreover, because the method is built on model tests, accurate calculations can be obtained [42].Falnes and Perlin [46,47] analyze the oscillating bodies in low-amplitude waves and obtain that the diffraction is negligible when the device dimension is considerably smaller than the wave length.Clement and Ferrant [48] described a method for the computation of free surface flows generated by submerged bodies, and obtain that radiation nonlinearities are negligible for floating bodies that are small as compared to the wave length.Merigaud et al. [43] added specific nonlinear terms to hydrodynamic models for wave energy devices, to improve the validity of such models across the full operational spectrum, showing that the response of the device is mainly affected by nonlinear FK forces, while nonlinear radiation and diffraction forces have minor effects on the system dynamics.
The focuses of the paper is on the novel heaving point absorber operating in the power production region, using the Froude-Krylov method to optimal shape design and maximize energy capture of the floating body for the power supply of the low-power unmanned ocean device.The remainder of this paper is organized as follows: Section 2 introduces the structure model of the novel heaving point absorber.Section 3 gives the mathematical model of the floating body of the novel absorber.Section 4 presents the algebraic solution of the wave excitation force in the vertical direction of floating bodies and nonlinear Froude-Krylov force integral.The numerical simulation and simulation analysis for the wave force, heaving velocity, heaving displacement, and capture width ratio of the multi-type floating bodies are in Sections 5 and 6.Some conclusions and final remarks are presented in Section 7.

Structure Model
The waves are generally generated by wind blowing across the surface of the ocean.The motion of waves is regular and periodic on the water surface, with the amplitude decreasing exponentially with depth.When the depth to be greater than half the wavelength, the wave-induced motions is only approximately 4% of those at the surface and thus could be considered to be insignificant [10].Therefore, this depth range is defined as the hydrostatic layer, in which the motion of waves is hardly perceived as the depth increases.Based on the above characteristics of wave motions, a novel heaving wave energy point absorber is designed, which can be used as a power module for the low-power unmanned ocean device [7]. Figure 1 shows the working scene of this point absorber.First, when the low-power unmanned ocean device is working on the ocean surface, the value of the battery energy is decreasing from the maximum to the lowest safe.Secondly, utilizing the release mechanism of the low-power unmanned ocean device, the novel heaving point absorber is released.The Underwater PTO subsystem of this absorber converts the captured wave energy into electricity.The value of the battery energy reaches the highest energy.Finally, utilizing the recovery mechanism of the ocean device, the absorber is recovered.The above process cycles back and forth, of which ensures the energy supply of unmanned ocean device for long-range and long-term work.showing that the response of the device is mainly affected by nonlinear FK forces, while nonlinear radiation and diffraction forces have minor effects on the system dynamics.The focuses of the paper is on the novel heaving point absorber operating in the power production region, using the Froude-Krylov method to optimal shape design and maximize energy capture of the floating body for the power supply of the low-power unmanned ocean device.The remainder of this paper is organized as follows: Section 2 introduces the structure model of the novel heaving point absorber.Section 3 gives the mathematical model of the floating body of the novel absorber.Section 4 presents the algebraic solution of the wave excitation force in the vertical direction of floating bodies and nonlinear Froude-Krylov force integral.The numerical simulation and simulation analysis for the wave force, heaving velocity, heaving displacement, and capture width ratio of the multi-type floating bodies are in Section 5 and Section 6.Some conclusions and final remarks are presented in Section 7.

Structure Model
The waves are generally generated by wind blowing across the surface of the ocean.The motion of waves is regular and periodic on the water surface, with the amplitude decreasing exponentially with depth.When the depth to be greater than half the wavelength, the wave-induced motions is only approximately 4% of those at the surface and thus could be considered to be insignificant [10].Therefore, this depth range is defined as the hydrostatic layer, in which the motion of waves is hardly perceived as the depth increases.Based on the above characteristics of wave motions, a novel heaving wave energy point absorber is designed, which can be used as a power module for the low-power unmanned ocean device [7]. Figure 1 shows the working scene of this point absorber.First, when the low-power unmanned ocean device is working on the ocean surface, the value of the battery energy is decreasing from the maximum to the lowest safe.Secondly, utilizing the release mechanism of the low-power unmanned ocean device, the novel heaving point absorber is released.The Underwater PTO subsystem of this absorber converts the captured wave energy into electricity.The value of the battery energy reaches the highest energy.Finally, utilizing the recovery mechanism of the ocean device, the absorber is recovered.The above process cycles back and forth, of which ensures the energy supply of unmanned ocean device for long-range and long-term work.The system configuration and working principle of the novel heaving point absorber are illustrated in Figure 2. The absorber mainly includes two parts, a floating body and an Underwater PTO, which are connected with each other by a steel cable.The floating body floats on the ocean surface and the Underwater PTO suspends the hydrostatic layer at a depth of about 40 m.The Underwater PTO mainly consists of a power generator, upper impeller, lower impeller, steady blade, and transmission shaft and planet-gear increaser.The impellers are mainly composed of The system configuration and working principle of the novel heaving point absorber are illustrated in Figure 2. The absorber mainly includes two parts, a floating body and an Underwater PTO, which are connected with each other by a steel cable.The floating body floats on the ocean surface and the Underwater PTO suspends the hydrostatic layer at a depth of about 40 m.The Underwater PTO mainly consists of a power generator, upper impeller, lower impeller, steady blade, and transmission shaft and planet-gear increaser.The impellers are mainly composed of blades, connecting rods, locking devices, center wheels, and external fixation rings.The center wheel and the external fixation ring are connected by eight radials arranged in connecting rods.Eight fan-shaped blades are fixed on corresponding connecting rods with locking devices and are arranged in a centrally symmetric circumferential array.In addition, the steel cable not only acts as dragging the Underwater PTO, but also transmissions electricity and control signals.The working principle of the novel heaving point absorber is shown in Figure 2. (1) During the movement of the floating body from the wave trough to the wave crest, the Underwater PTO is pulled up by the steel cable.The upper surface of the blade is impacted by the water flow and the blade adaptively swings downward, as shown in Figure 2a.Due to the limitation of the locking devices, the blade then stops swinging and is in a slanted state after reaching the maximum angle of inclination.The water flow continues to impact the slanted blades and propel the blade forward.The circumferential array of the blades enables the impeller to be subjected to circumferential thrust.Since the blades of the upper and lower impellers are arranged in opposite directions, the upper impeller is clockwise rotated by the water flow and the lower impeller is anticlockwise rotated by the water flow, and they are relatively reversed.The upper and lower impellers are fixedly connected with the stator and rotor of the generator, respectively, and then drive the generator to generate electricity.(2) During the movement of the floating body from the wave crest to the wave trough, the Underwater PTO sinks under the influence of gravity.The lower surface of the blade is impacted by the water flow and the blade adaptively swings upward, as shown in Figure 2b.Due to the limitation of the locking devices, the blade then stops swinging and it is in a slanted state after reaching the maximum angle of inclination.Since the direction of the water flow impinging on the blade does not change, the water flow continues to impact the slanted blades and propel the blade forward.The upper impeller is continuously clockwise rotated by the water flow and the lower impeller is continuously anticlockwise rotated by the water flow.Therefore, the direction of rotation of the generator does not change and it continues to generate electricity.According to the different impact directions of water flow, the impellers' blades adaptively adjust the blade deflection.The upper and lower impellers act as components that interact directly with the water flow and they provide continuous rotational motion to the generator during rising and sinking of the Underwater PTO.The working principle of the novel heaving point absorber is shown in Figure 2. (1) During the movement of the floating body from the wave trough to the wave crest, the Underwater PTO is pulled up by the steel cable.The upper surface of the blade is impacted by the water flow and the blade adaptively swings downward, as shown in Figure 2a.Due to the limitation of the locking devices, the blade then stops swinging and is in a slanted state after reaching the maximum angle of inclination.The water flow continues to impact the slanted blades and propel the blade forward.The circumferential array of the blades enables the impeller to be subjected to circumferential thrust.Since the blades of the upper and lower impellers are arranged in opposite directions, the upper impeller is clockwise rotated by the water flow and the lower impeller is anticlockwise rotated by the water flow, and they are relatively reversed.The upper and lower impellers are fixedly connected with the stator and rotor of the generator, respectively, and then drive the generator to generate electricity.(2) During the movement of the floating body from the wave crest to the wave trough, the Underwater PTO sinks under the influence of gravity.The lower surface of the blade is impacted by the water flow and the blade adaptively swings upward, as shown in Figure 2b.Due to the limitation of the locking devices, the blade then stops swinging and it is in a slanted state after reaching the maximum angle of inclination.Since the direction of the water flow impinging on the blade does not change, the water flow continues to impact the slanted blades and propel the blade forward.The upper impeller is continuously clockwise rotated by the water flow and the lower impeller is continuously anticlockwise rotated by the water flow.Therefore, the direction of rotation of the generator does not change and it continues to generate electricity.According to the different impact directions of water flow, the impellers' blades adaptively adjust the blade deflection.The upper and lower impellers act as components that interact directly with the water flow and they provide continuous rotational motion to the generator during rising and sinking of the Underwater PTO.

Mathematical Model
The floating body is an important part of the novel heaving wave energy point absorber and is in direct contact with the waves on the ocean surface.In order to increase the energy harvesting of floating body, it is necessary to increase the heave amplitude, heave velocity, and wave force in the vertical direction.It can increase the rising and sinking amplitude of the Underwater PTO and the relative rotational speed of the upper and lower impellers.Moreover, the power production of the novel heaving point absorber can also be increased.This section analysis of the force of the floating body, establishes its motion equation, and solves its wave excitation force and the Underwater PTO damping force.Also, this section derives the model of energy conversion efficiency between the floating body and the wave energy, and calculates its energy conversion efficiency.
The hydrodynamics model and energy efficiency model of the floating body is established to calculate the wave force and the conversion efficiency: Section 3.1 considers nonlinear Froude-Krylov force in the vertical direction and Underwater PTO damping force for the force analysis of the floating body in micro-wave amplitude.Section 3.2 presents the capture width ratio of the floating body for the energy efficiency analysis.

Hydrodynamics Model of the Floating Body
We consider that the floating body of the novel heaving point absorber floats freely in water of uniform depth.Under the action of linear regular waves, the floating body does micro amplitude heave motion.The fluid is assumed inviscid and the incident flow irrotational and incompressible.Figure 3 is the coordinates and force analysis of the floating body.The right-handed inertial reference frame is centered at the hydrostatic equilibrium position of the body.X is the vertical distance from the reference fluid surfaces to the body waterline (positive downwards).Y is the vertical distance from the reference fluid surface to the fluid surface (positive downwards).Y = r cos ϕt, where r is the wave amplitude and ϕ is the wave angular frequency.Z is the vertical distance from the reference fluid surface to any point in the fluid (positive downwards).D is the waterline depth of the floating body.Newton's second law can be used to describe the system dynamics, as follows: where, m is the mass of the novel heaving point absorber, X is the heaving displacement of the body from its hydrostatic equilibrium position, ..
X is the heaving acceleration of the body, G the gravity force, S the submerged surface, P the pressure, n a vector normal to the surface, Q the viscous damping force, I the inertia force of attached mass effect, and F PTO the Underwater PTO damping force.
The pressure P can be derived from the incident flow applying Bernoulli's equation: where, ρ the water density, g the acceleration of gravity, P st = −ρgz, hydrostatic pressure and φ the potential flow, which can be decomposed as the sum of the undisturbed incident flow potential φ I , the diffraction potential φ D , and the radiation potential φ R : The Airy's wave theory assumes that the motion of the floating body is a small amplitude.The solution of Equation ( 1) is solved around the equilibrium position of the buoy.Under the linear assumption, the wetted surface is constant.However, the nonlinear of Froude-Krylov forces to be considered in the actual calculation process.Thus, the wetted surface is exactly instantaneous, namely integrating the fluid pressure over the actual submerged portion of the buoy, as it moves through the water.Froude-Krylov forces include the static and dynamic forces.They depend on the instantaneous wetted surface, which depends both on the incident wave elevation and the displacement of the buoy.Froude-Frylov force can be written as: where, F FK st is the static Froude-Krylov force, given as the balance between the gravity force and the Archimedes force, and F FK dy is the dynamic Froude-Krylov force.
P dy (t)ndS (6) where, is the dynamic pressure.
Energies 2018, 11, 3282 7 of 21 through the water.Froude-Krylov forces include the static and dynamic forces.They depend on the instantaneous wetted surface, which depends both on the incident wave elevation and the displacement of the buoy.Froude-Frylov force can be written as: where, st FK F is the static Froude-Krylov force, given as the balance between the gravity force and the Archimedes force, and dy FK F is the dynamic Froude-Krylov force.The time-dependence annotation will be omitted for brevity hereafter.Equation (1) can be rewritten as [42]:

 
where, D F is the diffraction force and R F is the radiation force.
Note that, since the fluid is assumed to be inviscid, irrotational, and incompressible, no viscous force and inertia force appears in Equation (7).In addition, due to Froude-Krylov force already The time-dependence annotation will be omitted for brevity hereafter.Equation (1) can be rewritten as [42]: where, F D is the diffraction force and F R is the radiation force.Note that, since the fluid is assumed to be inviscid, irrotational, and incompressible, no viscous force and inertia force appears in Equation (7).In addition, due to Froude-Krylov force already includes the inertia force that is caused by the attached mass effect, the inertia force also does not appear in this Equation.Wave excitation force is the force of wind blows the waves to disturb the motion of the floating body on the ocean surface.This force can be written as: where, C is the diffraction correction coefficient.
The linear approach assumes that radiation and diffraction forces are linear.Therefore, the radiation and diffraction potential is negligible when the floating body dimension is considerably smaller than the wave length [43][44][45][46][47][48].When combining Equations ( 4)- (8), Equation (1) can be rewritten as: The algebraic calculation of the integral in Equation ( 9) requires the explicit definition of the dynamic pressure P dy , the infinitesimal surface element ndS, and the limits of integration.Under the Airy's wave theory for deep water waves, the dynamic pressure at any point on the floating body in the local coordinate system can be written as: where, H is wave height, x is the direction of wave propagation in the local coordinate system, z is the vertical displacement of the floating body in the local coordinate system (positive upwards), and d is water depth and ω is the wave circular frequency.Wave Number k is defined by the dispersion equation ktanh(kd) = ω 2 /g.During the motion of the novel heaving point absorber on the ocean surface, the coordinates at any point on the floating body in the global coordinate system can be written as: where, (X 0 , Y 0 , Z 0 ) is the coordinates of the center of gravity on the floating body in the global coordinate system and (x, y, z) is the coordinates at any point on the floating body in the local coordinate system.
[T] is coordinate transformation matrix between the local coordinate system and global coordinate system.Therefore, the dynamic pressure at any point on the floating body in the global coordinate system can be written as: where, x is the direction of wave propagation in the global coordinate system and z is the vertical displacement of the floating body in the global coordinate system (positive upwards).Furthermore, the force components of the wave excitation force in the vertical direction can be written as the following: where, n z is the projection of the normal for the wetted surface of the floating body in the vertical direction and C V is the diffraction correction coefficient in the vertical direction.Under the movement of the novel heaving point absorber on the ocean surface, the heave motion of the floating body does work by overcoming the damping of the water.The Underwater PTO subsystem of this absorber converts the captured wave energy of the body into electricity.In the motion model of this absorber under the Airy's wave theory, the linear damping model is used to analyze the damping force of the Underwater PTO system.Thus, the Underwater PTO damping force can be written as: where, B PTO is the Underwater PTO damping coefficient and .
X is the heaving velocity of the floating body.

Energy Efficiency Model of the Floating Body
In the researches of the WEC, the performance of power production of the device is mainly concentrated on the energy harvesting and energy conversion.Using the capture width ratio [49,50] of the floating body as an evaluation index to measure the characteristic of energy harvesting and conversion efficiency for the novel heaving point absorber.The capture width ratio is that the ratio of the average power of the energy harvesting to the wave energy of the incident wave within the width of the body, namely the efficiency of energy harvesting in the WEC.It can be written as the following: where, P WEC is the average power of the floating body to harvests wave energy and P 0 is the wave energy of the incident wave within the width of the body.
The average power of the energy harvesting is that the instantaneous power is integrated in a period and it solves the average value.It can be written as: where, F wave is the wave force, V WEC is the velocity of the floating body, and T is the wave period.The average power of the energy harvesting in the vertical direction of the floating body can be written as: where, F z PTO is the Underwater PTO damping force in the vertical direction and B z PTO is the Underwater PTO damping coefficient in the vertical direction.
The wave energy of the incident wave within the width of the floating body can be written as the following: where, B is the heading wave width of the floating body.

Algebraic Solution
The characteristic parameters, such as profile parameters, immersion depth, and wetted surface distribution, play an important role in the performance analysis, structural design, manufacture, and employment of the floating body.These parameters affect the motion amplitude, motion velocity, wave force, and wave energy harvesting of the floating body.Therefore, selects the profile parameters and solves the algebraic solution of the wave force is critical.Geometry shapes of floating body are two types of axisymmetric and unaxisymmetric.The axisymmetric body is a curved surface body formed by a generatrix rotating around a fixed vertical axis.The shape of the axisymmetric body is dependent on the shape of generatrix and relative position of the generatrix and the fixed axis.The axisymmetric bodies include cylinder, cone, sphere, and so on.The unaxisymmetric body is formed by a non-rotating curved surface.This body consists of rectangular, trapezoid, polyprism, Energies 2018, 11, 3282 10 of 20 and so on.In this paper, three kinds of shapes are selected for the floating body, such as rectangular, cylinder, and sphere.This section mainly solves the algebraic solution of the wave force of the above shape.Table 1 show that the algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.The geometric parameters of the axisymmetric body are described in cylindrical coordinates [42].The axisymmetric bodies are described in rectangular coordinates.
Table 1.The algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.

The Geometry Shapes of Floating Bodies
The Coordinates of Floating Bodies The Parameters of Profiles and Coordinates The Wave Excitation Force in the Vertical Direction of Floating Bodies Unaxisy-mmetric Horizontal rectangular

Algebraic Solution
The characteristic parameters, such as profile parameters, immersion depth, and wetted surface distribution, play an important role in the performance analysis, structural design, manufacture, and employment of the floating body.These parameters affect the motion amplitude, motion velocity, wave force, and wave energy harvesting of the floating body.Therefore, selects the profile parameters and solves the algebraic solution of the wave force is critical.Geometry shapes of floating body are two types of axisymmetric and unaxisymmetric.The axisymmetric body is a curved surface body formed by a generatrix rotating around a fixed vertical axis.The shape of the axisymmetric body is dependent on the shape of generatrix and relative position of the generatrix and the fixed axis.The axisymmetric bodies include cylinder, cone, sphere, and so on.The unaxisymmetric body is formed by a non-rotating curved surface.This body consists of rectangular, trapezoid, polyprism, and so on.In this paper, three kinds of shapes are selected for the floating body, such as rectangular, cylinder, and sphere.This section mainly solves the algebraic solution of the wave force of the above shape.Table 1 show that the algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.The geometric parameters of the axisymmetric body are described in cylindrical coordinates [42].The axisymmetric bodies are described in rectangular coordinates.

Numerical Results
Based on the above analysis of the hydrodynamic performance and energy conversion characteristic for the multi-type floating bodies of the novel heaving wave energy point absorber, the a is the length, b is the length along the wave direction, c is the height, λ(t) is the immer-sion depth in the time of t.

Algebraic Solution
The characteristic parameters, such as profile parameters, immersion depth, and wetted surface distribution, play an important role in the performance analysis, structural design, manufacture, and employment of the floating body.These parameters affect the motion amplitude, motion velocity, wave force, and wave energy harvesting of the floating body.Therefore, selects the profile parameters and solves the algebraic solution of the wave force is critical.Geometry shapes of floating body are two types of axisymmetric and unaxisymmetric.The axisymmetric body is a curved surface body formed by a generatrix rotating around a fixed vertical axis.The shape of the axisymmetric body is dependent on the shape of generatrix and relative position of the generatrix and the fixed axis.The axisymmetric bodies include cylinder, cone, sphere, and so on.The unaxisymmetric body is formed by a non-rotating curved surface.This body consists of rectangular, trapezoid, polyprism, and so on.In this paper, three kinds of shapes are selected for the floating body, such as rectangular, cylinder, and sphere.This section mainly solves the algebraic solution of the wave force of the above shape.Table 1 show that the algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.The geometric parameters of the axisymmetric body are described in cylindrical coordinates [42].The axisymmetric bodies are described in rectangular coordinates.
(See Appendix A for calculation details.)

Numerical Results
Based on the above analysis of the hydrodynamic performance and energy conversion characteristic for the multi-type floating bodies of the novel heaving wave energy point absorber, the R is the radius, l is the height, λ(t) is the immersion depth in the time of t, the cylinder equation is , then, dS = Rdαdz.

Algebraic Solution
The characteristic parameters, such as profile parameters, immersion depth, and wetted surface distribution, play an important role in the performance analysis, structural design, manufacture, and employment of the floating body.These parameters affect the motion amplitude, motion velocity, wave force, and wave energy harvesting of the floating body.Therefore, selects the profile parameters and solves the algebraic solution of the wave force is critical.Geometry shapes of floating body are two types of axisymmetric and unaxisymmetric.The axisymmetric body is a curved surface body formed by a generatrix rotating around a fixed vertical axis.The shape of the axisymmetric body is dependent on the shape of generatrix and relative position of the generatrix and the fixed axis.The axisymmetric bodies include cylinder, cone, sphere, and so on.The unaxisymmetric body is formed by a non-rotating curved surface.This body consists of rectangular, trapezoid, polyprism, and so on.In this paper, three kinds of shapes are selected for the floating body, such as rectangular, cylinder, and sphere.This section mainly solves the algebraic solution of the wave force of the above shape.Table 1 show that the algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.The geometric parameters of the axisymmetric body are described in cylindrical coordinates [42].The axisymmetric bodies are described in rectangular coordinates.R is the radius, l is the heig- ht, ) (t λ is the immersion depth in the time of t , the cylin- der equation is (See Appendix A for calculation details.)

Numerical Results
Based on the above analysis of the hydrodynamic performance and energy conversion characteristic for the multi-type floating bodies of the novel heaving wave energy point absorber, the R is the radius, λ(t) is the immersion depth in the time of t, the spherical equation is

Numerical Results
Based on the above analysis of the hydrodynamic performance and energy conversion characteristic for the multi-type floating bodies of the novel heaving wave energy point absorber, the wave force, heaving velocity, heaving displacement, and capture width ratio of the multi-type floating bodies were numerically simulated and comparatively analyzed.
During the numerical simulation, the parameters of the multi-type floating bodies are as follows: the length of horizontal rectangular a is 0.3 m, the length along the wave direction of horizontal rectangular b is 0.3 m, the height of horizontal rectangular c is 0.5 m, the radius of vertical cylinder R is 0.3 m, the height of vertical cylinder l is 0.5 m, the radius of sphere R is 0.3 m, the mean immersion depth is 0.15 m, the Underwater PTO damping coefficient in the vertical direction B z PTO is 20 KNs/m, and the mass of the novel heaving point absorber m is 10 kg.According to the model test results of previous research [51][52][53], the diffraction correction coefficient in the vertical direction C V of horizontal rectangular, vertical cylinder, and sphere can be obtained, respectively.
In this study, three reference sea states are used for evaluation the energy efficiency of the novel heaving wave energy point absorber.Sea state 1: the significant wave height H S is 0.1 m, the peak wave period T P is 2 s; Sea state 2: the significant wave height H S is 0.3 m, the peak wave period T P is 3.5 s; and, Sea state 3: the significant wave height H S is 0.5 m, the peak wave period T P is 5 s.The first sea state covers H S values from 0.0 to 0.1 m, the second sea state covers H S values from 0.1 to 0.3 m, and the third sea state covers H S values from 0.3 to 0.5 m.The combination of H S and T P is representative for the same wave tank, the wave climate date come from the simulation test.The model scales of the novel heaving point absorber, as follows: the Underwater PTO's diameter is 410 mm, the distance between the upper and lower impeller is 400 mm, and the overall height of the Underwater PTO is about 560 mm.The scale between the dimensions of the tested heaving point absorber and the full scale device is 1:5.Another, the significant wave height H S is about 2.5 m in the real wave climate with the full scale device.The considered reference water depth d is 50 m.In addition, the incident wave is the linear regular wave, the angle between the direction of the linear regular wave and the direction of the novel heaving point absorber is 0 • .According to the above parameters of the incident wave, the working scene of the novel heaving point absorber is the deep-water waves.Therefore, the other parameters of the incident wave are obtained, including wavelength, wave velocity, and velocity of water particles.These parameters are as follows: the wavelength is 6.T is representative for the same wave tank, the wave climate date come from the simulation test.The model scales of the novel heaving point absorber, as follows: the Underwater PTO's diameter is 410 mm, the distance between the upper and lower impeller is 400 mm, and the overall height of the Underwater PTO is about 560 mm.The scale between the dimensions of the tested heaving point absorber and the full scale device is 1:5.Another, the significant wave height S H is about 2.5 m in the real wave climate with the full scale device.The considered reference water depth d is 50 m.In addition, the incident wave is the linear regular wave, the angle between the direction of the linear regular wave and the direction of the novel heaving point absorber is 0°.According to the above parameters of the incident wave, the working scene of the novel heaving point absorber is the deep-water waves.Therefore, the other parameters of the incident wave are obtained, including wavelength, wave velocity, and velocity of water particles.These parameters are as follows: the wavelength is 6.

Discussion
As shown in Figures 4-6, the vibration frequency curve of the wave force, heaving velocity, and heaving displacement are less of the same for the vertical cylinder buoy and the sphere buoy.However, as to the former, the curve of the wave force and the heaving displacement is steeper, the heaving velocity is faster, wave follower is better, and the generated energy is higher.Nevertheless, the wave force and heaving displacement of the horizontal rectangular floating body are smaller than those two types of floating bodies.Therefore, the generated energy and wave follower is weaker by this one less than those two types of floater bodies.In Figure 7, the capture width ratio of the vertical cylinder floating body is higher and more stable than the sphere floating body.While, the horizontal rectangular floating body is much more less than those two types of floating bodies.
In addition, as shown in Figures 8-11, the wave force, heaving velocity, heaving displacement, and capture width ratio of the floating bodying is affected by the peak wave period and the wave Energies 2018, 11, 3282 14 of 20 height at a certain sea state.The above parameters are become higher when the sea state increases.Another, the vertical cylinder floating body outperforms the other floating body in the three reference sea states.As shown in Figure 11, the curve of the capture width ratio is steeper and the value of this parameter is higher when the sea state big changes.
In summary, with the linear regular wave, the cylindrical floater vertically placed in the wave surface is the first optional shape for the novel heaving point absorber and follower is the sphere floater, which can increase the quality of power extracting and the efficient of the WEC system design.At last, the horizontal rectangular floating body is carefully selected.

Conclusions
This paper presents a small novel heaving point absorber of energy supply for low-power unmanned ocean devices that are based on the counter-rotating self-adaptive mechanism, with the advantages of small space device, stability, and reliable energy conversion process.For improving the efficiency of this absorber's power production, the wave force and energy efficiency are analyzed by the Froude-Krylov method and the optimal floating body is selected, the following conclusions are drawn: (1) The structure model and working principle of the novel heaving point absorber are feasible.
Based on the optimal power supply strategy, the design of the release and recovery mechanism can improve the concealment, release, and recovery rate of the unmanned ocean devices.Under the different directions of water flows, the design of the upper and lower impeller with opposite rotation directions can provide continuous relative rotational motion to the generator during rising and sinking of the Underwater PTO.The design of the locking devices and impellers' blades with circumference array can achieve the adaptively adjust the blade deflection in the heaving motion and the relative rotation of the impellers.The design of the steady blade can balance the motion of Underwater PTO and makes it not limited by the motion of heaving.(2) The energy efficiency of the novel heaving point absorber is greatly affected by the geometric parameters, the wave force, the motion displacement, the motion velocity, and the capture width ratio of the floating body.The constitutive relation of the above parameters of the floating body is constructed.In order to calculate the above parameters, the Froude-Krylov method is used, which effectively overcomes the inaccuracy of the linear models and reduces the time consuming to simulate.The algebraic solution of wave excitation force for the axisymmetric and unaxisymmetric floating body is obtained and validated by the numerical simulation.(3) The wave force, heaving velocity, heaving displacement, and capture width ratio of the three floating bodies are compared and analyzed by the numerical simulation.Under the same working condition, mass and mean immersion depth, the type of vertical cylindrical floater's slopes of the wave force and the steeper amplitude, faster speed of heave motion, larger amplitude of heave motion, better follow wave, larger energy produced, as well as, higher conversion efficiency.Therefore, with the linear regular wave, the cylindrical floater vertically placed on the wave surface is the first optional shape for the novel heaving point absorber and follower is the sphere floater, which can increase the quality of power extracting and the efficient of the WEC system design.
The study that was carried out in this paper focuses on the energy generation efficiency of the novel heaving point absorber for the supply power of the low-power unmanned ocean devices.Moreover, the multi-type floating bodies are optimized base on the Froude-Krylov method.Although a series of numerical analyses above has been conducted in this paper, research on the novel absorber and floating body is still not thorough enough.Therefore, a further step toward is that the verification experiments of wave tank will be done in the power production region for the novel absorber.In addition, the geometric parameters of the vertical cylinder floater will be optimized for increasing the conversion efficiency of the floating body.

Figure 1 .
Figure 1.The working scene of the novel heaving wave energy point absorber.

Figure 1 .
Figure 1.The working scene of the novel heaving wave energy point absorber.
rods, locking devices, center wheels, and external fixation rings.The center wheel and the external fixation ring are connected by eight radials arranged in connecting rods.Eight fan-shaped blades are fixed on corresponding connecting rods with locking devices and are arranged in a centrally symmetric circumferential array.In addition, the steel cable not only acts as dragging the Underwater PTO, but also transmissions electricity and control signals.

Figure 2 .
Figure 2. The system configuration and working principle of the novel heaving point absorber.(a) Rising state; and, (b) Sinking state.

Figure 2 .
Figure 2. The system configuration and working principle of the novel heaving point absorber.(a) Rising state; and, (b) Sinking state.
the dynamic pressure.

Figure 3 .
Figure 3.The coordinates and force analysis of floating body.

Figure 3 .
Figure 3.The coordinates and force analysis of floating body.

Table 1 .a
The algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.is the length, b is the length along the wave direction, c is the height, depth in the time of t .

R
is the radius, l is the heig- ht, ) (t λ is the immersion depth in the time of t , the cylin- der equation is the time of t , the spherical equation is

Table 1 .a
The algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.is the length, b is the length along the wave direction, c is the height, depth in the time of t .

R
is the radius, l is the heig- ht, ) (t λ is the immersion depth in the time of t , the cylin- der equation is the time of t , the spherical equation is

Table 1 .a
The algebraic solutions of the wave excitation force in the vertical direction of multi-type floating bodies.is the length, b is the length along the wave direction, c is the height, depth in the time of t .
π π−arccos( R−λ(t) R ) sh(kR cos θ)shk(d + R −λ(t)) • J 0 (kR sin θ) • cos θ sin θdθ (See Appendix A for calculation details.) 2/19/39 m, the wave velocity is 3.1/5.5/7.8 m/s, and the velocity of water particles is 0.16/0.27/0.31m/s.The results of numerical simulation are shown in Figures 4-11.The wave force, heaving velocity, heaving displacement, and capture width ratio of the multi-type floating bodies in the second sea state is shown in Figures 4-7, respectively.Figures 8-11 show that the above parameters for comparison in the three reference sea states and the values obtained from the numerical simulation.S S from 0.1 to 0.3 m, and the third sea state covers S H values from 0.3 to 0.5 m.The combination of S H and P 2/19/39 m, the wave velocity is 3.1/5.5/7.8 m/s, and the velocity of water particles is 0.16/0.27/0.31m/s.The results of numerical simulation are shown in Figures 4-11.The wave force, heaving velocity, heaving displacement, and capture width ratio of the multi-type floating bodies in the second sea state is shown in Figures 4-7, respectively.Figures 8-11 show that the above parameters for comparison in the three reference sea states and the values obtained from the numerical simulation.

Figure 4 .
Figure 4.The wave force of the multi-type floating bodies (H = 0.3 m, T = 3.5 s).

Figure 4 .
Figure 4.The wave force of the multi-type floating bodies ( H = 0.3 m, T = 3.5 s).

Figure 7 .
Figure 7.The capture width ratio of the multi-type floating bodies ( H = 0.3 m, T = 3.5 s).

Figure 7 .
Figure 7.The capture width ratio of the multi-type floating bodies ( H = 0.3 m, T = 3.5 s).Figure 7. The capture width ratio of the multi-type floating bodies (H = 0.3 m, T = 3.5 s).

Figure 8 .
Figure 8.The wave forces of the multi-type floating bodies in the reference sea state.

Figure 8 .
Figure 8.The wave forces of the multi-type floating bodies in the reference sea state.

Figure 8 .Figure 9 .
Figure 8.The wave forces of the multi-type floating bodies in the reference sea state.

Figure 10 .
Figure 10.The heaving displacement of the multi-type floating bodies in the reference sea state.

Figure 11 .
Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.

Figure 9 .
Figure 9.The heaving velocity of the multi-type floating bodies in the reference sea state.

Figure 8 .Figure 9 .
Figure 8.The wave forces of the multi-type floating bodies in the reference sea state.

Figure 10 .
Figure 10.The heaving displacement of the multi-type floating bodies in the reference sea state.

Figure 11 .
Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.

Figure 10 .
Figure 10.The heaving displacement of the multi-type floating bodies in the reference sea state.

Figure 8 .Figure 9 .
Figure 8.The wave forces of the multi-type floating bodies in the reference sea state.

Figure 10 .
Figure 10.The heaving displacement of the multi-type floating bodies in the reference sea state.

Figure 11 .
Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.

Figure 11 .
Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.Figure 11.The capture width ratio of the multi-type floating bodies in the reference sea state.