#### Abstract

Electricity generation from oceanic wave depends on the wave dynamics and the behavior of the ocean. In this paper, a permanent magnet linear generator (PMLG) has been designed and analyzed for oceanic wave energy conversion. The proposed PMLG design is suitable for the point absorber type wave energy device. A mathematical model of ocean wave is presented to observe the output characteristics and performance of the PMLG with the variation of ocean waves. The generated voltage, current, power, applied force, magnetic flux linkage, and force components of the proposed PMLG have been presented for different sea wave conditions. The commercially available software package ANSYS/ANSOFT has been used to simulate the proposed PMLG by the finite element method. The magnetic flux lines, flux density, and field intensity of the proposed PMLG that greatly varies with time are presented for transient analysis. The simulation result shows the excellent features of the PMLG for constant and variable speeds related to wave conditions. These analyses help to select proper PMLG parameters for better utilization of sea wave to maximize output power.

#### 1. Introduction

At present, scientists and engineers are facing two major challenges in the world of energy: electrical energy production and environmental issue. These problems can be amicably solved using renewable energy resources (RERs). The vital factors which have stimulated the use of RERs are energy independence, financial viability, and mainly environmental protection [1]. The problem of oil crises (1973–1983) and environmental pollution concerns urged engineers to harvest electrical energy from the available RERs, for example, wind energy [2–5], solar energy [6–9], and hydro power [10]. As the RER is variable and unpredictable, different types of converters and controls are associated with the RER-based power plants that are connected with standalone or grid systems [11, 12]. The traditional RERs have the problem of uncertainty of availability and they require a large land area. On the other hand, the key advantages of oceanic wave energy (OWE) are as follows: (i) it has huge potential compared to the solar and wind energy, (ii) it is easy to forecast, and (iii) it does not need land area. Hence, it has been of great interest to the industrial field; particularly, in energy generation, the use of wave power has been more attractive compared to other RERs [13]. OWE is a promising environmental pollution-free energy, which would make significant contribution toward saving biochemical resources and reducing carbon emissions [14]. It is estimated that the total wave energy resource in the open sea around the world is 10 TW (10,000 GW), a comparable amount of the total power consumption in the world [15].

Different types of oceanic wave energy converters (WECs) have been invented and examined for successful conversion of wave energy into electrical energy [16–19]. The wave energy device may be of a rotational or a translational type and each of these devices has different features [20]. The power takeoff devices play a vital role in converting the irregular wave motion to a regular motion for energy conversion. Linear generators (LGs) have directly been implemented to the direct drive wave energy conversion without using medium devices which has unique advantages over all other wave energy devices [21]. An LG has two major parts, namely, translator and stator. The translator mounted to a hollow cylinder is sometimes called a float or floater as shown in Figure 1. The stator is mounted to another mechanical body in order to make it stationary with respect to the sea wave.

The floater tries to float on the surface of the ocean wave with the translator that moves relative to the stator due to wave action. The floater utilizes the rise and fall of the sea wave at a single point for energy conversion. As the translator moves linearly due to the reciprocating motion of the floater, the wave energy so extracted is converted into electrical energy. A lot of mathematicians analyzed and proposed different mathematical models for understanding the nature of oceanic wave [22, 23].

Different types of permanent magnet linear generators (PMLGs) have been designed and analyzed for improved performance. The flat and tubular flux switching permanent magnet linear generators (FSPMLGs) have been proposed [21, 24, 25]. Different analyses have shown that the FSPMLGs made of PMs and steel cores have suffered from the problem of higher leakage flux leading to the reduction of electrical power generation. Tubular PMLGs have been proposed [14, 21, 26, 27] to reduce cogging forces and also increase efficiency. The maintenance of tubular PMLG is difficult due to the presence of the coils inside the periphery. It is seen in the recent works [28] that the linear switched reluctance generator (LSRG) has been proposed for high power generation. Excessive leakage flux and complex control circuit of LSRGs are responsible for the degradation of overall efficiency and reliability.

This paper has presented the mathematical model of the wave motion to analyze the behavior of a PMLG that offers low internal resistance, low loading effect, and high output power. It is essential to consider the nature of wave motion for parameter selection of the PMLG, thus maximizing electricity generation. A relationship between the oceanic wave motions with the parameters of PMLG is established. Different significant parameters, for example, size of PMs, poles, translator length, and stroke length, of the PMLG are found from the relationship and the generated voltage, current, power, applied force, magnetic flux linkage, flux lines, flux density, and the applied force with force components of the proposed PMLG are shown for different sea wave conditions.

#### 2. Design of the PMLG

##### 2.1. Working Principle

The translator moves vertically with the incident wave; therefore, the direction may be upward or downward. Considering the translator movement in the upward direction as shown in Figures 2 and 3 for a particular time interval, the stator and translator poles are aligned facing each other. The red and green lines are representing the north pole (N) and the south pole (S) of the permanent magnet (PM), respectively. S exists in the upper side and N exists in the lower side of the stator core for a time . The position of the translator varies with time and the direction of magnetic flux changes. The translator position changes in Figure 3 from the position shown in Figure 2. S now exists in the lower side and N exists in the upper side of the stator cores for another time .

Hence, the direction of magnetic flux according to Figure 2 is opposite for the translator position as in Figure 3 within the time interval . Therefore, the induced voltage directions across the conductor are opposite to each other and due to this reason the PMLG generates AC power.

##### 2.2. Vector Diagram and Equivalent Circuit

The winding of each phase consists of two coils of opposite phases. The induced voltages are therefore 180° phase shifted to each other. The coils of each phase are connected in series in additive polarity to get higher output voltage. The vector diagrams of induced voltages are represented in Figure 4.

The equivalent circuit diagram of the proposed PMLG for a three-phase load is shown in Figure 5. The equivalent series resistances and inductances of each phase winding coil are considered equal, denoted by and , respectively. The winding consists of two series coils; therefore, the equivalent series resistance and inductance of each coil are denoted by and , respectively. , , and are the induced voltages of phase-A, phase-B, and phase-C, respectively, for simplicity. Similarly, , , and are the induced voltages of phase-, phase-, and phase-, respectively. The magnetic excitation is fed from the translator’s PM array as shown in Figure 6. The terminal voltages, , , and , are measured across the load.

The induced voltage equation may be represented aswhere , , and , , , and −1, is the terminal voltage, is the line current, and is synchronous inductance. is the constant representing the machine construction; is the translator vertical velocity or speed; is the pole pitch; and is the vertical displacement. The terminal voltage is

##### 2.3. Construction Details

The vertical cross section of the PMLG (front side) is shown in Figure 6. The PMLG basically contains a translator which is made by some PMs with the steel core. The stator contains some copper coils wounded on the stator cores of steel situated on both sides of the translator. Phase-, phase-, and phase- are located on the right side of the translator. Phase-, phase-, and phase- are located on the left side of the translator which are 180° phase shifted from phase-, phase-, and phase-, respectively. The construction supports the translator to move in the vertical direction with respect to the stator. The orientation of the PM array should be such that N and S can be formed one after another as shown in Figure 6.

#### 3. Model of the Oceanic Wave

In most of the cases, constant speeds or sinusoidal speeds having sinusoidal shapes are common approximations for simulating PMLGs [25–27]. The typical range of vertical velocity of wave is 0–2 m/s with a time period from 4 to 6 s [28–30]. Therefore, a free oceanic wave neither forced nor dissipated on a flat seabed is presented for analysis of the LG as shown in Figure 7 following the mentioned approximations. The amplitude of waveform is . Hence, the vertical distance between the wave crest and the trough is which is equal to twice the amplitude . According to the point of view of oceanographers, it is considered as a linear wave. The sea surface is lying on plane, where components are considered zero as the wave propagates in -direction, is wavelength, is water depth of the ocean with respect to plane or water surface, seabed is at , and the ocean surface coincides with .

##### 3.1. Generalized Wave

The general description of almost any type of oceanic wave may be considered as, depending only on the interpretation of ,Here, the number of waves and is the frequency that asserts the physics and describes consideration of a water wave relating frequency and wave number. Alternatively, it can be considered as a relation between the phase speed or oceanic wave velocity, , and the wavelength. Considering gravitational acceleration, = 9.8 m/s^{2}, the phase speed which is a single basic wave that moves along the -direction can be expressed asAccording to (3) and (5), the wave velocity is along -direction only although it has the velocity along -direction. So, the velocity of oceanic wave is a vector that depends on , , and and time, . To obtain a complete description of oceanic wave, the components along , , and have to be calculated. The motion can be expressed as where represents the velocity of oceanic wave and , , and are the components along , , and . The components may be considered as zero for the oceanic wave and the and components can be expressed as follows:The and components from (7) and (8) can be simplified for individual consideration of shallow water wave, intermediate depth wave, and deep sea wave.

##### 3.2. Shallow Water Wave

In shallow water, water depth is much lower than the wavelength, ; that is, or . Another property of shallow water wave is that the amplitude of wave is much smaller than wavelength; so, or . So, (7) and (8) may be written as Now, (9) and (10) may be reduced asAgain, if , can be simplified as because As , phase speed can be expressed as .

##### 3.3. Deep Water Wave

In deep water, water depth is much higher than wavelength; that is, or . So,Therefore, (7) and (8) can be expressed as follows: Again, if , can be simplified as . Phase speed, , can be expressed as . The and components, both for shallow and for deep water from (6)–(12), (15), and (16), are explained in [31]. Comparisons of relative depth, wave speed, and wavelength of shallow water wave, intermediate depth wave, and deep sea wave are shown in Table 1.

#### 4. The PMLG and Wave Model

##### 4.1. Selection of Translator Length

The translator lengths may be shorter, equal, or longer compared to the stator. The common lengths of the stator and the translator affect the power rating of the PMLG. The power rating increases with increase in the common length. Longer translator has been chosen to obtain the same common length as shown in Figure 8 for different translator positions. The vertical displacement of translator with the incident sea wave is known as stroke length, . There is a relationship between the wave amplitude, , and during the energy conversion as . is smaller than the amplitude of free oceanic wave as in (7)–(12), (15), and (16) because of power dissipation.

##### 4.2. Determination of Frequency

The translator pole pitch is the summation of the PM thickness and translator pole thickness of the PMLG as shown in Figure 9. The stator pole pitch and translator pole pitch may be the same or different depending on the design strategy. In the proposed design, the stator pole pitch is the same as the translator pole pitch and the stator pole width is the same as the translator pole width. The pitch and pole width have a vital effect on the generated voltage, shape, power, frequency, and forces. The frequency, , of the PMLG is determined by

##### 4.3. Direction of Forces

The simulation setup of the PMLG design is represented using a 3D Cartesian coordinate system. The applied force is working along -axis, as shown by bidirectional arrow in Figure 10. The width and thickness of this design are along -axis and -axis, respectively. Cogging forces between stator and translator core act along -axis. Other force components except for applied force and cogging forces work along -axis. The PMLG converts mechanical energy to electrical energy due to the force along -axis and the forces along -axis and -axis are generated which is the cause of mechanical power loss. In the oceanic wave model, the direction of wave elevation is also along -axis according to (8), (10), (12), and (16).

##### 4.4. Generated Voltage with Waves

According to Faraday’s law of electromagnetic induction, the induced voltage is found as In the wave model, the wave elevation is described with respect to -axis that represents the direction of wave propagation. The induced voltage and related parameters of the PMLG are related to time. The vertical wave displacement is assumed to be sinusoidal and the translator connected to buoy tries to follow the wave elevation; therefore, the vertical displacement, , and velocity, , of the translator can be expressed as given in (19) and (20), respectively [14]. Hence,

Here, , , and represent wave amplitude, initial phase angle, and period of oceanic wave, respectively. The translator and buoy move with the wave elevation; therefore, the flux variation with respect to time can be expressed as If is the peak voltage, combining (18) and (21), the generated EMF per phase and of the PMLG, , are obtained as [32]

#### 5. Simulation Results

Two types of speed settings have been used in this simulation; one is consideration of constant translator speed for observation of the performance of the PMLG and the other is translator motion changed with the incident sea wave. The default air gap length is 2 mm, load is 4 Ω, and the translator speed is 1 m/s for the PMLG unless otherwise specified. The input of the PMLG is mechanical thrust, , in newton. From mechanical power, , and electrical power, , efficiency, , is calculated by .

##### 5.1. Constant Translator Speed

The terminal voltage and load current for default condition are shown in Figure 11.

The load voltage and current are in the same phases due to the resistive load. The applied force, , cogging force, , and force component along -axis, , are shown in Figure 12. The induced voltage, current, and magnetic flux linkage are shown in Figure 13. The generated power is shown in Figure 14.

##### 5.2. Voltage Regulation

The terminal voltage depends on the load for a specific PMLG as in (2), which is shown in Figure 15 to observe the terminal characteristics.

Different loads ranging from 2 Ω up to 1 MΩ are shown on a semilog graph. Load voltages measured from any one of the three phases are similar and remain almost unchanged from 25 Ω to 1 MΩ. The RMS values of currents and power of phase-A for different load conditions are shown in Figures 16 and 17, respectively. The numerical values of the current and power graphs used in the legend are tabulated in Tables 2 and 3, respectively.

##### 5.3. Variable Translator Speed

According to the oceanic wave model, the vertical position of the wave is sinusoidal, so the vertical velocity of the wave elevation is a cosine function. The vertical position and velocity of the translator are shown in Figure 18. The time period = 4 s and m which means = 0.8 m. The magnitude of the generated voltage of the PMLG depends on velocity according to (20) and is shown in Figure 19.

As = 4 s, the voltage waveform is repeated in each half cycle or s as shown in Figure 19. The voltage waveform due to the translator velocity for the time period of one cycle is shown in Figure 20. The induced voltage and current for a half-cycle time interval are shown in Figure 21. The terminal voltage lags by 90° from the magnetic flux linkage as shown in Figure 22.

The voltage magnitude is directly proportional to the translator velocity. According to (17) and (20), the frequency of voltage increases with the increase in velocity. The terminal voltage and current in a winding are shown in Figure 23. The cogging force is the reason of force ripples and is liable for abnormal operation of the PMLG, so it is maintained at a lower value as shown in Figure 24 with the applied force. The magnitude of power is related to the frequency. The power generation due to the applied force is shown in Figure 25.

Different dimensions of the PMLG active materials of the stator and translator as shown in Figure 6 are given in Table 4. The significant parameters that mainly affect the PMLG performance are given in Table 5.

##### 5.4. Selection of the PM Size

According to Figure 18, the vertical displacement of the translator is 1.6 m owing to the minimum and maximum vertical displacement of −0.8 m and 0.8 m, respectively. On the other hand, the translator length is the summation of stator length and stroke length. Combining the stroke length of 1.6 m and the stator length of 1.4 m, the translator length becomes 3 m. From Table 5, the translator and stator pole pitch is 4 cm which is the summation of the PM and the pole thickness. Therefore, 75 PMs and poles are required to make 3 m length of translator. The dimension of the translator PMs is selected as 8 × 10 × 2.4 cm.

#### 6. Analysis of Magnetic Flux

The magnetic flux density, , flux lines, and magnetic field intensity of the PMLG are analyzed to observe in the different portions. The value of for different time, , is shown in Figures 26–28. varies for different positions of the translator which can be realized from the scale given in Figure 29(a). The magnetizing curve of the steel used in the stator and translator core is shown in Figure 29(b).

**(a)**

**(b)**

The translator position varies with time and the flux lines travel through the low reluctance path. The magnetic flux lines for different times are shown in Figures 30–32.

There is a relation between magnetic flux lines and . Figures 33(a)–33(c) represent the scale of flux lines for ms, 30 ms, and 40 ms, respectively. The magnetic properties inside the PMLG are needed for analyses, proper design, and performance check.

**(a)**

**(b)**

**(c)**

The values for the magnetic field intensity at different times are shown in Figures 34–36 and the scale is shown in Figure 37. The magnetomotive force mainly exists in the air gap. When the stator and translator cores come close together, the magnetomotive force reaches a high value. Neodymium iron boron (NdFeB) permanent magnets have been used for magnetic excitation.

#### 7. Conclusions

The simulation results with synchronized translator vertical velocity with wave velocity reflect the model presented in this paper. Different voltages, currents, power, and magnetic flux linkages of phase-A are shown for the three-phase PMLG because the other two phases are just 120° phase shifted from each other and are not necessary to explain the dynamics of energy conversion from wave energy to electrical energy. The cogging force and force ripples are low compared to applied force that helps to prevent mechanical vibration and also minimize mechanical power loss. Low core loss is achieved (1.27%) which would produce less heat to prevent overheating of the PMLG and increase service life. The maximum efficiency is calculated as 83% for full-load condition. The loading effect is low and voltage regulation of the PMLG is 13.93% measured from terminal characteristics. Thick copper conductor with smaller turn number is used in the coil to minimize internal resistance.

#### Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.