Reach Us
+44-1522-440391

**Fuat Kara ^{*}**

School of Water, Energy and Environment, Cranfield University, UK

- *Corresponding Author:
- Fuat Kara

School of Water

Energy and Environment

Cranfield University, UK

**Tel:**44 (0) 1234 754118

**E-mail:**[email protected]

**Received Date**: December 09, 2016; **Accepted Date:** January 12, 2017; **Published Date**: January 16, 2017

**Citation: ** Kara F (2017) Control of Wave Energy Converters for Maximum Power Absorption with Time Domain Analysis. J Fundam Renewable Energy Appl 7: 221. doi:10.4172/20904541.1000221

**Copyright:** © 2017 Kara F. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Fundamentals of Renewable Energy and Applications

A discrete control of latching is used to increase the bandwidth of the efficiency of the Wave Energy Converters (WEC) in regular and irregular seas. When latching control applied to WEC it increases the amplitude of the motion as well as absorbed power. It is assumed that the exciting force is known in the close future and that body is hold in position during the latching time. A heaving vertical-cylinder as a point-absorber WEC is used for the numerical prediction of the different parameters. The absorbed maximum power from the sea is achieved with a three-dimensional panel method using Neumann-Kelvin approximation in which the exact initial-boundary-value problem is linearized about a uniform flow, and recast as an integral equation using the transient free-surface Green function. The calculated response amplitude operator, absorbed power, relative capture width, and efficiency of vertical-cylinder compared with analytical results.

Time domain; Transient free-surface green function; Boundary integral equation; Absorbed power; Relative capture width; Latching control

The issue of global warming has resulted in growing pressure to switch to renewable sources (e.g. solar, wind, wave, tidal, geothermal, ocean thermal etc.) to reduce carbon emissions. Of these, ocean wave power has received much attention in Europe and other part of the world. The global power potential represented by waves that hit all coasts worldwide is estimated to be in the order of 1 TeraWatt (TW) [1]. If wave energy loss due to friction and wave breaking may be utilised, then the global wave power input is estimated to be one order of magnitude larger (˜10 TW), a quantity that is comparable with the world’s present power consumption. The idea of converting the energy of ocean waves into useful energy forms is not new. There are techniques that were first patented as early as 1799 [2]. The pioneer work for the wave energy utilisation was introduced in the late 1940s by Masuda [3], Salter [4], Budal et al. [5], McCormick [6], Mei [7]. For an axisymmetrical WEC in regular waves, it is known that the maximum power that can be absorbed equals the incident wave power associated with a wave front of width one wave length divided by 2p. This result was first derived independently by Budal et al. [8], Evans [9] and Newman [10].

Power from ocean waves can be converted useful energy (e.g. electrical energy) by the use of Oscillating Water Columns (OWC) or OverTopping Devices (OTD) exploiting power of the seawater or oscillating devices (e.g. point absorber devices) in near shore or offshore environment exploiting body motion of the devices which excited by the incident waves [11]. OWC and OTD devices have large horizontal extension and large broad bandwidth whilst point absorber devices can take incident waves at any direction and have narrow bandwidth which allows the use of optimum control strategies (e.g. latching control) effectively. It is possible to convert wave energy to some other form; mechanical, hydraulic or pneumatic at very high efficiencies. The problem of converting this primary energy to a marketable product has presented the greatest challenge to various WEC inventors. The main reason for the very difficult engineering problems is the fact that waves transfer their energy in a slow irregular fashion whereas current technology, particularly for electrical generation, is more used to fast regular movements.

WEC without control traps energy only from the wave components whose frequency lies in the band of the natural frequencies of the system. Nevertheless, this system is potentially able to catch energy in the wave spectrum outside this band provided it can be dynamically controlled to adapt its bandwidth in real time. In this way, the power output from WEC can be increased by controlling the oscillation in order to approach an optimum interaction between WEC and incident wave. Optimum control methods can be classified for WEC as reactive (continuous) and latching (discrete) controls. Applying reactive control to WEC allows absorption of the maximum power [12]. In reactive control, the applied force/moments are a continuous function of time, and generally include resistive and reactive components to approach resonance at wave frequencies lower than the natural frequency of WEC. In latching control [11] while the motion is locked, the wave exciting force continues to evolve. Latching control involves braking and releasing the oscillating body such that its velocity bears the desired phase relationship with the wave-applied exciting force.

The damping of WEC must be properly adjusted in order to obtain approximate optimum amplitude. This is accomplished by Power- Take-Off (PTO) system. PTO mechanism is idealized by damping coefficients, giving forces proportional to the velocity of WEC. In the context of mathematical modelling, since the objective with WEC is to induce a movement such that radiated waves cancel in amplitude and phase with the scattered component to a maximum extend, the accurate predictions of radiated and diffracted waves due to forced oscillation of WEC and exciting forces, respectively are crucially important. This ensures that as little energy as possible is in the resultant outgoing wave field and therefore the maximum energy is observed by WEC. The accurate and computationally efficient numerical model was developed by the author [13,14] for the prediction of the hydrodynamics and power characteristics of WEC in regular waves. In the present paper, latching control for vertical-cylinder as a wave energy converter will be used to maximize the absorbed energy from both regular and irregular ocean waves. Equation of motion is considered in section 2 while solution of the integral equation is covered in section 3. Importance of latching control strategy in order to maximise the captured power is studied in section 4 whilst the main conclusions are given in section 5.

The right-handed coordinate system is used to define the fluid action. A Cartesian coordinate system is fixed to the body. Positive
x-direction is towards the bow, positive z-direction points upwards, and the z=0 plane (or xy plane) is coincident with calm water. The body is
translating through an incident wave field with velocity U0, while it undergoes oscillatory motion about its mean body surface position. The origin
of the body-fixed coordinate system is located at the centre of the xy plane. The solution domain consists of the fluid bounded by the
free surface S_{f} (t), the body surface S_{b} (t), and the boundary surface at infinity S_{8} as in shown in **Figure 1**.

Through Newton’s law, the dynamics of a floating body’s unsteady oscillations are governed by a balance between the inertia of the floating body and the external forces acting upon it. This balance is complicated by the existence of radiated waves; as a consequence both of the floating body’s own motions and its scattering of the incident waves. This means that waves generated by the floating bodies at any given time will persist indefinitely and the waves of all frequencies will be generated on the free surface. These generated waves, in principle, affect the fluid pressure field and hence the body force of the floating bodies at all subsequent times. This situation introduces memory effects and is described mathematically by a convolution integral. Having assumed that the system is linear, the equation of motion of any floating rigid body may be written in a form which is essentially proposed by Cummins [15].

(1)

The displacement of the floating bodies from its mean position in each of its six rigid-body modes of motion is given x_{k} in Equation (1), and the
overdots indicates differentiation with respect to time. The inertia matrix of the floating body is M_{jk}, damping due to PTO system is given as PTO B and linearized hydrostatic restoring force coefficients are given by C_{jk}. The radiation impulse response functions are composed of the hydrodynamic
coefficients and the kernel of the convolution on the left-hand side of Equation (1). A radiation impulse response function is the force on the
body in the j-th direction due to an impulsive velocity in the k-th direction, with the coefficients a_{jk}, b_{jk}, c_{jk}, accounting for the instantaneous forces proportional to the acceleration, velocity, and displacement, respectively, and the memory function K_{jk}(t) accounting for the free surface effects which persist after the motion occurs. For the radiation problem the term ‘memory function’ will be used to distinguish this portion of the impulseresponse function from the instantaneous force components outside of the convolution on the left-hand side of Equation (1). The term K_{jD}(t) on the right hand side of Equation (1) are the components of the exciting force and moment’s impulse response functions due to the incident wave elevation ζ_{0}(t), defined at a prescribed reference point in the body-fixed coordinate system. Here, the kernel K_{JD}(t) is the diffraction impulse response function; the force on the body in the j-th direction due to a uni-directional impulsive wave elevation with a heading angle of ß. Once the stiffness matrix, inertia matrix, and fluid forces e.g. radiation and diffraction forces are known, the equation of motion of floating body Equation (1) may be solved using the fourth order Runge-Kutta method.

The initial boundary value problem consisting of initial condition, free surface and body boundary condition may be represented as an integral equation using a transient free surface Green’s function [16]. Applying Green’s theorem over the transient free surface Green function derives the integral equation. It is possible to show that transient free surface Green’s function satisfies the initial boundary value problem without a body [17]. Integrating Green’s theorem in terms of time from -8 to 8 using the properties of transient free surface Green’s function and potential theory, the integral equation for the source strength on the body surface may be written as in Kara [18].

(2)

And potential on the body surface

(3)

Where,

is the memory part of the transient free surface source potential, the field point, the source point, the distance between field and source point and represent the Rankine part of the source potential, the distance between field point and image point over free surface, J0 the Bessel function of zero order. The Green function represents the potential at the field point P (x, y, z) and time t due to an impulsive disturbance at source point and time τ.

The integral equation for the source strength Equation (2) is first solved, and then this source strength is used in the potential formulation Equation (3) to find potential and velocities at any point in the fluid domain. The solution of the integral equation Equation (2) is done using time marching scheme. The form of the equation is consistent for both the radiation and the diffraction potentials so that the same approach may be used for all potentials. Since the transient free surface Green function satisfies free surface boundary condition and condition at infinity automatically, in this case only the underwater surface of the body needs to be discretized using quadrilateral/triangular elements. The resultant boundary integral equation Equation (2) is discretized using quadrilateral elements over which the value of the source strength is assumed to be constant and solved using the trapezoidal rule to integrate the memory or convolution part in time. This discretization reduces the continuous singularity distribution to a finite number of unknown source strengths. The integral equation Equation (2) has satisfied at collocation points located at the null points of each panel. This gives a system of algebraic equations which are solved for the unknown source strengths. At each time step the new value of the source strength is determined on each quadrilateral panel.

The evaluation of the Rankine source type terms (e.g. 1/r, 1/r’) in Equation (2) is analytically integrated over quadrilateral panels using the method and formulas of Hess and Smith [19]. For small values of r the integrals are done exactly. For intermediate values of r a multi-pole expansion is used. For large values of r a simple monopole expansion is used.

The surface and line integrals over each quadrilateral element involving the wave term of the transient free surface Green function are solved analytically [20-22] and then integrated numerically using a coordinate mapping onto a standard region and Gaussian quadrature. For surface elements the arbitrary quadrilateral element is first mapped into a unit square. A two-dimensional Gaussian quadrature formula of any desired order is then used to numerically evaluate the integral. The line integral is evaluated by subdividing Γ(t) into a series of straight line segments. The source strength σ (t) on a line segment is assumed equal to the source strength of the panel below it.

The evaluation of the memory part of the transient free surface Green function and its derivatives with an efficient and accurate method
is one the most important elements in this problem. Depending on the values of P, Q and t five different methods are used to evaluate ; power series expansion, an asymptotic expansion, a Filon integration quadrature, Bessell function, and asymptotic expansion of complex error
function. **Figure 2** shows the memory part of transient free surface Green function. The detail analysis of radiation and diffraction
problem including Froude-Krylov force as well as numerical results that compared with analytical, experimental and other published numerical
results can be found in Kara [13,14].

Latching control, which is a discrete real time control, will be used in the present paper. Rather than adapting WEC parameters to the excitation force in order to optimize the linear body response, the latching control adapts the body response to WEC and to the excitation in a nonlinear fashion. It is a kind of parametric resonance adaptation process as can be found in nonlinear oscillatory theory, this kind of behaviour can be predicted only using time domain simulations. Latching control can magnify the amplitude of the motion whatever the frequency of the excitation force, and can improve the efficiency of WEC in terms of absorbed energy for excitation frequencies. When latching control is applied, an additional force must be introduced in the dynamic of WEC to cancel the acceleration of the controlled motion in order to lock the system temporarily. The latching control of WEC consists of locking the oscillating body in position at the instant when velocity vanishes, and releasing it after a certain delay to be determined. This latching delay has to be applied in order to maximise the response amplitude of the body. The instant of latching is imposed by the dynamics of the body itself (i.e., vanishing velocity); thus the control variable is simply the duration of the latching phase, or equivalently the instant of release [23-27]. One of the advantages of latching control is that it is passive, which means that it does not need to deliver energy to WEC while it is engaged, since the forces do no work as long as the velocity vanishes.

**Numerical results for regular waves**

A vertical-cylinder with 8 m radius R and 13 m draft Tvc is used as
an offshore floating wave energy converter in the following numerical
results. The natural period Tn of the wave energy converter is 10 s.
Geometry of vertical-cylinder is given in **Figure 3**.

**Motion of wave energy converters with and without latching
control:** Budal and Falnes [11] point out that maximizing the extracted power requires that the wave force and the velocity of the wave energy
converters be in phase. This condition is achieved when the velocity is
vanishing. This implies that the wave energy converter is at its maximum
displacement. The wave energy converter is hold until a prescribed time
before the exciting force reaches its next extreme. This prescribed time
is called the latching time which is a control variable. Our ITU-WAVE
time domain program takes into account this latching strategy for a
change in sign of the velocity which starts the latching.

In **Figure 4**, the numerical simulation in the time domain of the
body response with and without latching control is plotted. The huge
magnification of the amplitude of the motion for the controlled system
compared to the uncontrolled one is observed. The amplitude of the heave motion is increased more than by a factor of 20 in the case of
latching control.

**Figure 5** shows velocity of wave energy converter with 5 s latching
and without latching control. The device has 10 s resonance frequency
and incident wave period is 15 s. One can observe a large amplification
of the velocity of the device when latching control is applied as in
motion amplitude of the heaving vertical-cylinder in **Figure 4**.

It is pointed out earlier that maximizing the extracted power requires
that the wave force and the velocity of the wave energy converters be in
phase when latching control is applied as shown in **Figure 6**. Latching
duration can be observed from **Figure 6**.

It is assumed that the excitation force is known in the future of
the current computation time. The motion of a wave energy converter
under latching control is a succession of the time of the rest and time of
the transient. The control variable is the release (or unlatching) time or
equivalently the latching duration. We pointed out earlier that the aim
of the latching control is to maximize the amplitude of the oscillations of
the device. This implicitly means that it will also maximize the absorbed
power at the same time. **Figure 7** shows an example of numerical
simulation of the body motion in a different incident wave period with
and without latching control. It can be seen that the amplifying effect of
latching control is significant in **Figure 7**. It may be noticed from **Figure
7** that the response amplitude can be as much as by a factor of 50 in the
case of latching control in higher incident wave periods.

**Instantaneous and mean absorbed power:** A damping force
proportional to the velocity of the floating body is considered as
the action of the external Power-Take-Off (PTO) mechanism. The instantaneous power converted by the wave energy converters is
defined as:

(4)

Where, P_{e} is the power due to exciting force and P_{r} is the power due
to radiation force.

**Figure 8** shows the instantaneous power absorbed for each incident
wave period of 15s from ocean waves using a vertical-cylinder as a
wave energy converter with and without latching control. In the case
of latching control the absorbed instantaneous power is increased
significantly. It may be noticed that the unlatching results are very small
in terms of controlled latching results.

The mean power extracted by the PTO over a time range T is given by,

(5)

The averaging time T must be much larger than the characteristics
period of the wave, which is typically 5 s to 15 s. The maximum mean
power that can be absorbed is equal to the linear damping proportional
to velocity of the floating body. This is achieved when the velocity is in
phase with the exciting force [11] as it is shown in **Figure 6**.

**Figure 9** shows the absorbed mean power for the range of incident
wave frequencies using a vertical-cylinder as a wave energy converter
with and without latching control. In the case of latching control the
absorbed mean power is again increased significantly. The theoretical
maximum power absorption in regular seas is given by Budal and
Falnes [11].

(6)

Where, being the incident wave amplitude and being incident wave frequency.

**Capture width and relative capture width:** The capture width l and
maximum capture width l_{max} is defined as [11]:

(7)

Respectively where is the mean power and given by Equation (5), is the wave power in the incident wavetrain per unit crest length, being the incident wave amplitude. A good wave absorber is a body which has the ability when making waves, to concentrate the wave energy along a narrow sector rather than distribute the energy evenly over all angles. The maximum capture width equals to for an axisymmetric system in symmetric mode of motion e.g. heave. This implies that the floating body absorbs all the power in an incident wave equal to that passing a crest length of . Relative capture width is obtained by dividing capture width with width B of the wave energy converters, .

**Figure 10** shows the relative capture width for with and without
latching controlled system and the theoretical maximum in the context
of linear theory are plotted. For such a single axisymmetric system in
symmetric mode of motion e.g. heave the maximum capture width as
pointed out earlier equals to λ/(2Π). It may be noticed that it means that
for low frequency waves, the absorbed power can possibly be greater
than the incident power in the width of the body, implying capture
widths greater than 1 and this can be observed in **Figure 10**. It can
be also observed that latching control allows an improvement of the
relative capture width whatever the frequency of the incident waves are.

**Efficiency:** The efficiency η of wave energy converters is defined as , which has a maximum of 1.0 for any wavelength.

**Figure 11** shows the efficiencies for vertical-cylinder in the case of
in resonance and off-resonance periods. The efficiency converges to 1.0
(100% efficient) at resonance period T_{n}=10 s whereas off-resonance case
T=11 s shows a very low efficiency (5%).

**Figure 12** shows the efficiency plotted against the circular frequency
for a vertical-cylinder which has 10 s (0.628 rad/s) natural period. If the
natural period of vertical-cylinder equals the period of incident waves
in the case of without latching control, the device is perfectly tuned
and we expect optimal efficiency. As the difference between natural
period of device and incident wave period increases, the efficiency
of the system decreases. As can be seen in **Figure 12** latching control
increases the bandwidth of the wave energy converter for all frequency ranges. If off-resonance period is 11s (0.571 rad/s), the efficiency is
approximately 5% without latching control. However, if 1s latching
is applied to vertical-cylinder wave energy converter, it is possible to
achieve an efficiency of approximately 100%.

**Numerical results for irregular waves**

The Pierson-Moskowitz spectral formulation [28] developed
from analysis of measured data obtained in the North Atlantic by
wave recorders installed on weather ships will be used for irregular seas. Analysis for this spectrum was carried out only on selected wave
records considered to have been acquired in fully developed sea. The
Pierson-Moskowitz spectrum can be presented in terms of significant
wave height H_{S}. That is,

(8)

The Pierson-Moskowitz spectral formulation depends on a
single parameter, significant wave height, Hs. The relationship
between significant wave height Hs and peak frequency ω_{p} is given as [29]. The theoretical maximum power absorption P
in irregular seas is similar to the regular wave prediction and is given by
Budal and Falnes [11].

(9)

The benefit of the latching control depends on the spectral
characteristics of the incident wave. A wave energy converter will
work with wave environment in which peak period ω_{p} is in the range
of between 5 s and 15 s. The control variable that the body velocity
vanishes remains invariant as in regular waves as well as the moment of
unlatching in random sea environment.

**Figure 13** shows the incident wave height which is the input for the
numerical calculation of the exciting forces of diffraction and Froude-
Krylov forces, heave displacement and velocity of the wave energy
converter which is the response of the floating structure in the case
of with and without latching control, and accumulated instantaneous
absorbed power from the ocean.

**Figure 13:** Comparison of different variables with and without latching control
at natural period Tn=10 s and peak period Tp=10 s of vertical-cylinder wave
energy converter (a) incident wave elevation at the origin of the wave energy
converter (b) heave motion amplitude (c) heave velocity (d) accumulation of
instantaneous power.

**Figure 14** shows the velocity and exciting force in phase in irregular
wave case. The same strategy as in regular sea is used e.g. vanishing
velocity that starts latching duration.

**Figure 15** shows the convergence of the mean power. Prediction of
the mean power for regular and irregular seas is done differently. In the
case of regular seas, mean power is calculated for each incident period **(Figure 8)**, while the accumulation of the instantaneous power is used
for the irregular seas. It may be noticed in **Figure 15** convergence of
mean power requires more time in the case of irregular seas.

**Figure 16** shows the absorbed mean power for the range of peak
frequencies using a vertical-cylinder as a wave energy converter
with and without latching control. In the case of latching control the
absorbed mean power is again increased significantly.

Results are plotted in terms of relative capture width in **Figure 17** as a function of the peak frequencies. For each result the relative capture
width of the same device without latching control is represented for
the purpose of comparisons. As it is compared to the regular seas, the
relative capture width is reduced almost 37% and occurs in higher peak
frequencies.

**Figure 18** shows convergence of the efficiency of wave energy
converters in time at around peak frequency of 10s with and without
latching control. As compared to regular seas, the convergence of
efficiency needs more time. One of the reasons for this is that the
accumulation of instantaneous power is used for the irregular seas
rather than calculation of instantaneous power for each period.

**Figure 19** shows the efficiency of vertical-cylinder with and without
latching control. The general trend shows that latching control improves
considerably the efficiency of the wave energy converters especially as
the peak frequency of the wave environment increases. The maximum
efficiency in irregular seas is around 60% as it is compared to 100% in
regular seas.

A latching control with a transient hydrodynamic code for wave energy converters are used for the prediction of maximum power absorption from ocean waves for a point absorber heaving verticalcylinder. Results have been presented to demonstrate the convergence of the method for power prediction of the floating wave energy converter and the calculation are shown to be in satisfactory agreement with analytical results. Results show that the efficiency of wave energy converters is considerably improved by the latching control which enlarges the bandwidth of WEC especially in the low frequencies, if the exciting force is predicted in the close future of the unlatching time.

- Panicker NN (1976) Power resource potential of ocean surface waves. Proceedings of the wave and salinity gradient workshop, UAS.
- Girard (1799) French Patent (translated by AE Hidden). Queens University,Belfast, UK.
- Masuda Y (1986) An experience of wave power generator through tests and improvements. Hydrodynamics of Ocean Wave Energy Utilization.pp: 445-452.
- Salter SH (1974) Wave power. Nature 249: 720-724.
- Budal K, Falnes J (1974) Proposals for conversion of the energy in ocean waves.Internal Report, Norwegian University of Science and Technology.
- McCormick M (1974) Analysis of a wave-energy conversion buoy. J Hydronautics 8: 77- 82.
- Mei CC (1976) Power extraction from water waves. J Ship Res 20: 63-66.
- Budal K, Falnes J (1975) A resonant wave point absorber of ocean waves. Nature 256: 478-479.
- Evans D (1976) A theory for wave-power absorption by oscillating bodies. J Fluid Mech 77:1-25.
- Newman JN (1976) The interaction of stationary vessels with regular waves. Proceedings of the Eleventh Symposium on Naval Hydrodynamics.pp: 491-501.
- Budal K, Falnes J (1976) Optimum operation of wave power converter. Internal Report, Norwegian University of Science and Technology.
- Salter SH, Jeffery DC, Taylor JRM (1976) The architecture of nodding duck wave power generators. The Naval Architect.pp: 21-24.
- Kara F (2010) Time domain prediction of power absorption from ocean waves with latching control. Renew Energy 35: 423-434.
- Kara F (2016) Time domain prediction of power absorption from ocean waves with wave energy converters arrays. Renew Energy 92: 30-46.
- Cummins WE (1962) Theimpulse response function and ship motions. Shiffstechnik 9: 101-109.
- Wehausen JV, Laitone EV (1960) Surface waves. Handbuch der Physik 9: 446-778.
- Finkelstein AB (1957) Theinitial value problem for transient water waves. Commun Pure Appl Math 10: 511-522
- Kara F (2000) Time domain hydrodynamics and hydroelastic analysis of floating bodies with forward speed.
- Hess JL, Smith AM (1964) Calculation of non-lifting potential flow about arbitrary three-dimensional bodies. J Ship Res 8: 22-44.
- Liapis SJ (1985) Time domain analysis of ship motions.
- King BK, Beck RF (1987) Time domain analysis of wave exciting forces.
- Newman JN (1990) Theapproximation of free surface green functions. Massachusetts Institute Technol.pp: 107-135.
- Iversen LC (1982) Numerical method for computing the power absorbed by a phase-controlled point absorber. ApplOcean Res 4: 173-180.
- Greenhow M, White SP (1997) Optimal heave motion of some axisymmetric wave energy devices in sinusoidal waves. Appl Ocean Res 19: 141-159.
- Eidsmoen H (1998)Tight-moored amplitude limited heaving buoy wave energy converter with phase control. Appl Ocean Res 20: 157-161.
- Korde UA (2002) Latching control of deep water wave energy device using an active reference. Ocean Engineering 29: 1343-1355.
- Babarit A, Duclos G, Clement AH (2004) Comparision of latching control strategies for a heaving wave energy device in random sea. Appl Ocean Res 26: 227-238.
- Pierson WJ, Moskowitz L (1964) A proposed spectral form for fully developed wind seas based on the similarity theory of S.A.Kitaigorodskii. J Geophys Res 69: 5181-5190.
- Ochi MK (2008) Ocean waves: Thestochastic approach. Cambridge Ocean Technology Series 6 .

Select your language of interest to view the total content in your interested language

- Adsorption Process
- Advanced Photovoiltoic Cells
- Advanced Photovoiltoic Panels
- Advanced Photovoiltoic Systems
- Advanced Photovoiltoic Technologies
- Alternative Energy
- Biomass Energy
- CO2 Mitigation
- Clean Energy
- Coal Energy
- Composting
- Energy Management
- Environmental Policy
- Fossil Fuel
- Geothermal Energy
- Green Energy
- Greenhouse
- Harvesting
- Harvesting Field Greens
- Heavy Metals
- Humic Substance
- Hydro Electric Energy
- Hydro Energy
- Hydrogen Energy
- Hydrogen and Biodiesel
- Hydropower Energy
- Immobilization
- Incineration Research
- Industrial Waste
- Modern Biomass Biofuel
- Organic Matter
- Petrolium and Petrochemistry
- Photovoltoics
- Renewable Energy
- Renewable Energy and Research
- Renewable Geothermal Energy
- Resource Efficiency
- Solar Energy (Thermal and Electrical)
- Solar Radiation
- Waste Hierarchy
- Wind Energy

- Total views:
**2734** - [From(publication date):

January-2017 - Dec 13, 2019] - Breakdown by view type
- HTML page views :
**2497** - PDF downloads :
**237**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2019-20