Next Article in Journal
Case Study: Optimizing Grading Ring Design for High Voltage Polymeric Insulators in Power Transmission Systems for Enhanced Electric Field and Voltage Distribution by Using a Finite Element Method
Previous Article in Journal
Clean Energy Stocks: Resilient Safe Havens in the Volatility of Dirty Cryptocurrencies
Previous Article in Special Issue
Development of Cost Correlations for the Economic Assessment of Power Plant Equipment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Non-Equilibrium Thermodynamic Approach for Analysis of Power Conversion Efficiency in the Wind Energy System

1
Department of Electric Mechatronics and Computer-Controlled Electromechanical Systems, Lviv Polytechnic National University, 79013 Lviv, Ukraine
2
Faculty of Electrical Engineering, Czestochowa University of Technology, 35-959 Czestochowa, Poland
*
Authors to whom correspondence should be addressed.
Energies 2023, 16(13), 5234; https://doi.org/10.3390/en16135234
Submission received: 7 June 2023 / Revised: 28 June 2023 / Accepted: 5 July 2023 / Published: 7 July 2023

Abstract

:
This article proposes an approach and develops an appropriate method of applying linear non-equilibrium thermodynamics to analyze energy processes, in particular using the example of the wind energy conversion system (WECS) with a directly connected vertical axis wind turbine (VAWT) and vector-controlled permanent magnet synchronous generator (PMSG). The main steps of the proposed approach are the description of the component subsystems as universal linear or linearized energy converters (ECs), which are characterized by several dimensionless parameters, the main one of which is the degree of coupling between their input and output. According to their value, as well as justified efficiency criteria, the optimal operating points of each ECs can be easily found. Such an approach makes it possible to abstract from physical laws of a different nature and equally assess the work of each of the subsystems. The next step is a connection of the received ECs. As shown in the paper, for the most common cascade connection of ECs, there are the best conditions for their connection, under which the newly formed equivalent EC can have maximum efficiency. This opens up an opportunity to analyze the influence of already real parameters of cascaded interconnected subsystems on the quality of their connection and justify specific solutions that would not have been seen without this approach. For example, in this study, from all parameters of the PMSG, only the selection of the optimal rated inductance of the armature winding made it possible to improve the quality of the connection of the PMSG with a specific VAWT and approximate the efficiency of the entire WECS to the maximum possible, especially in medium and high winds.

1. Introduction

The efficiency of energy conversion is a crucial factor in achieving energy savings within various industries’ technological processes. While renewable energy utilizes natural and freely available energy sources, enhancing the efficiency of energy transformation in renewable energy systems contributes to increased power generation at specific equipment capacities and helps reduce the cost of the obtained energy.
The analysis of energy efficiency in various systems relies on the fundamental science of thermodynamics. Initially, during its early development in the first half of the 19th century, initiated by the work of Sadi Carnot, thermodynamics primarily focused on describing reversible processes and isolated systems with the significant presence of thermal energy. The name of the science itself, “thermodynamics,” reflects this emphasis on thermal energy. Classical thermodynamics, developed by Rudolf Clausius in the latter half of the 19th century, played a pivotal role in establishing a fundamental science that enables the analysis of energy transformations in reversible processes. Clausius’s approach, based on a phenomenological macroscopic description of systems encompassing diverse physical and later chemical phenomena, provided a framework for studying energy transformations without considering the specific timing of their occurrence [1]. This makes it possible to outline the thermodynamic limits of process efficiency, which has not lost its relevance to this day. Linear non-equilibrium thermodynamics, which was developed from classical thermodynamics in the middle of the 20th century thanks to the works of Lars Onsager, summarized the general theoretical foundations for the energy regularities of not only isolated, but also closed and open systems with real irreversible processes, taking into account their flow speed. As a result, a clear mathematical apparatus was developed, especially for linear systems, which allows for the unification and description of energy transformations [2]. Thermodynamics found its further development in the works of F. Curzon and B. Ahlborn on modeling real systems by combining reversible and irreversible subsystems. This section of thermodynamics was named endoreversible thermodynamics [3]. It makes it possible to establish relationships between efficiency and power output, which is especially important for the study of natural biochemical processes. Currently, due to its fundamental nature, thermodynamics is used to describe energy transformations in all types of systems of various natures—physical, chemical, biological—and is even used to study regularities in informational, social systems, etc. Its comprehensiveness can be judged from the work [4], where all natural and anthropogenic processes taking place on planet Earth are subjected to thermodynamic analysis.
Renewable energy, which combines systems of different natures, between which energy exchange and energy transformation take place, has recently become an object of thermodynamics research. This also applies to wind energy [5], where the energy of the wind flow (aerodynamic subsystem) is transformed into the mechanical energy of the wind turbine (WT) and then, through an electromechanical transformation in the generator and electrical transformation in the power electronic system controlling its operation, into electrical energy. Most research on WTs is aimed at improving their effectiveness. As for small WTs, the search for new efficient configurations, which would be easy to start up and work productively in conditions with small and rapidly changing wind speeds, is constantly ongoing [6]. For high-capacity WTs, which are typically installed in conditions characterized by high and stable wind speeds, the configurations of the WTs themselves are already well established. However, there are new problems associated with the organization of wind farms, for example, wind farm arrangements, variations in WT capacities, effective control in different wind ranges [7], production interactions in the substantial offshore wind farms [8], increasing the reliability of work under specific operating conditions, and reductions of acoustic noise, in particular, the control of WTs according to the noise and power trade-offs graphics [9]. Although the vast majority of research in wind energy is devoted to many of its other important problems, thermodynamic analysis is already used in new works, and energy and exergy terminology is introduced [10]. For example, in [11], the comparison of the impacts on WT efficiencies of four meteorological variables, such as wind speed, pressure, temperature, and the humidity ratio was investigated. As a result of research, it was shown that although wind speed dominates the turbine’s efficiency performance, other meteorological variables also play important roles. More complex multi-generation hybrid renewable energy systems are also possible [12], with the further conversion of renewable energy in storage systems, for example, in the potential energy of compressed air [13], production of hydrogen or ammonia [14], novel combined cooling, heating systems using wind energy [15], and complex systems with a combination of different methods of energy generation and accumulation [16]. In similar works, thermodynamic analysis is used to assess the quality of energy and its exergy and search for optimal solutions [17].
Such studies analyze the flow of energy or exergy and its transformation in the specific subsystems of the WT, multiplicator, and generator, but do not determine ways to increase the energy efficiency of the entire wind energy conversion system (WECS). To identify the latter, a more detailed analysis is required, which is based on a mathematical description of the processes in the WECS components. However, a detailed description of interconnected processes of different physical natures that take place in WECS components—aerodynamic, mechanical, electromagnetic—complicates the understanding and does not reveal all the problems associated with a decrease in energy efficiency. An effective approach to solving the problem of increasing WECS efficiency can be the application of methods of linear thermodynamics of non-equilibrium processes, which made it possible to successfully solve many complex problems, especially with coupled effects. A feature of the linear thermodynamics of non-equilibrium processes is the ability to describe complex, including non-linear steady-state processes of various physical, chemical, and biological systems, according to the universal principle of linearization of the relationship between input and output power coordinates at specific operating points. At the same time, the system is considered an energy converter (EC) with a certain number of inputs and outputs, but most often with one input and one output. For the latter case, a universal method of describing the EC using a system of dimensionless parameters and characteristics of its performance has been developed [2,18]. The most important indicator of an EC is the degree of coupling between its input and output, which shows what part of the input power is transferred to the EC output. Another part of the input power “slips through”, which means that it is spent on organization of the process, but is not directly transferred to the output. The degree of coupling uniquely determines the maximum values of EC performance indicators, as well as the operating points at which these maximum values are reached. This makes it possible to optimize EC operation modes depending on the assumed optimization criteria [18]. This universal method was most widely used in bioenergetics to describe free energy transformations in biosystems, but it can be successfully used to analyze established processes of energy transformations in systems of any nature. For example, in [19], it was successfully applied to the analysis of energy conversion in electromechanics.
This work should be considered the first in the direction of the development of an applied method of the linear thermodynamics of non-equilibrium processes for further successful applications to complex objects of energy systems, in particular renewable energy systems. The universal approach to the mathematical description of systems of different natures, which is characteristic for this method, makes it possible to evaluate, from a common point of view, as universal ECs, each of the subsystems of a complex system with their parameters and maximum efficiency possibilities. As will be shown later in the work, significant reserves for improving the efficiency of the entire system lie in improving the conditions of the cascade connection of individual ECs, which most often takes place in different energy systems. In the work, a universal method of EC description was developed in the direction of the regularities of the cascade connection of two ECs, which makes it possible to evaluate the quality indicator of a specific connection and its influence on the efficiency indicators of the obtained equivalent cascade EC. For the successful application of the proposed method, a flow chart of the proposed approach of the investigation is given in Section 2.4. As a first test of this approach, the method of linear thermodynamics of non-equilibrium processes is applied to analyze the efficiency of energy conversion in the main WECS subsystems and to identify hidden opportunities for improving the efficiency of the entire system. In this work, the systems of linear equations were obtained for the first time for the steady-state modes of operation of two main ECs, which are the part of the WECS—a WT with a vertical axis of rotation (VAWT) and a vector controlled synchronous generator with permanent magnets (PMSG). Based on the obtained linear equations, these two ECs were described using a universal method. Taking into account the non-linear nature of steady-state processes in the VAWT and PMSG, systems of linear equations are obtained by linearizing the characteristics of these devices at operating points, and the mathematical model provides multiple linearization around the points of the optimal VAWT angular velocity at different wind speeds. In this paper, a universal method of EC description was also developed in the direction of the regularities of the cascade connection of two ECs, which makes it possible to evaluate the quality indicator of a specific connection and its influence on the efficiency indicators of the obtained equivalent cascade EC. Numerous studies conducted on the example of the studied WECS with specific VAWTs and PMSGs showed the imperfection of the cascade connection of the latter, which increases with decreasing wind speed. The influence of parameters and characteristics of VAWT and PMSG on the quality of this connection was studied, and ways to improve it were outlined.
The paper is structured as follows: in Section 2, the main provisions of linear non-equilibrium thermodynamics and the universal method of describing and evaluating the performance indicators of linear ECs are briefly highlighted; in addition, this description is developed in relation to the cascade connection of two ECs and the method of energy optimization of the entire system steady-state operation mode is outlined; a mathematical description of the aeromechanical regularities of the WT as a linear EC and an evaluation of the effectiveness of the studied VAWT is carried out in Section 3; in Section 4, a mathematical description of the regularities of the operation of the vector-controlled PMSG as a linear EC was carried out, and the efficiency of the generator for the studied WECS was evaluated; in Section 5, an analysis of the efficiency of energy conversion in the “VAWT—PMSG” complex was carried out, and the influence of its parameters and characteristics on improving the quality of the connection between the studied VAWT and PMSG was investigated; concluding remarks are given in Section 6.

2. Indicators of EC Performance According to Linear Non-Equilibrium Thermodynamics

2.1. Basic Provisions of Linear Non-Equilibrium Thermodynamics

For any open system that exchanges energy and matter with the external environment, the entropy change is expressed by two components: d e S —external, due to the exchange with the external environment (entropy flow); d i S —internal, caused by non-equilibrium processes within the system (production of entropy). Based on this, the generalized form of the second law of thermodynamics for non-isolated systems is expressed based on the following dependency rate of entropy production in the system [3]:
σ s = d i S d t 0 .
Since every real irreversible process is accompanied by some internal flow J k , generated by the corresponding driving force X k , which expresses disequilibrium, it has been proven [3] that in the conditions of local thermodynamic equilibrium (LTE), the following equality takes place:
J k = i L ki X i .
The principle of LTE or microscopic reversibility, which is the basis of linear non-equilibrium thermodynamics, is that in the non-equilibrium state of the body during processes that are not very rapid, the state of its very small elements can be considered equilibrium, and the equations of classical thermodynamics can be applied to it [2]. The product of σ s based on the absolute temperature of the system is called the dissipative function, because it reflects the rate of dissipation (scattering, degradation) of the free energy of the system.
A phenomenological law is also valid for LTE conditions, which expresses the intensity of any flow in the form of products of thermodynamic forces X i and the corresponding kinetic coefficients L ki (Onsager principle of linearity) [2,3]:
σ s = i J k X k .
Expression (3) shows that one flow can be caused by different thermodynamic driving forces. However, according to the Curie principle [2], different driving forces can simultaneously affect the transfer if they are tensors of the same rank, for example, only scalars or only vectors. In Expressions (2) and (3), flows and forces have a vector characteristic, which often happens in practice.
According to the second law of thermodynamics (1), the sum in Equation (2) must be greater than zero. However, this does not mean that each of the members of this sum must be individually greater than zero. Cases are possible when some members are positive, that is, they produce entropy, while others are negative, that is, and they consume an excess amount of entropy. This phenomenon, which consists of the fact that some processes can induce other processes to flow in a thermodynamically unfavorable direction (against the driving force, “downhill”), was called thermodynamic coupling [18]. Thus, for two coupling processes to which any EC belongs, since its output parameters (o—output) are determined by the input (i—input), the following system of equations can be written:
{ J i = L ii X i + L io X o J o = L oi X i + L oo X o ,
where the kinetic or “phenomenological” coefficients are from the expressions
L jk = ( J j X k ) X j = const .
By definition, X i , X o , and J i are positive, and J o is negative, which means the absorption of energy at the input and its formation at the output; therefore, L ii , L oo and − L io must be positive.
L. Onsager discovered a general principle called Onsager reciprocal relations that can be called the fourth law of thermodynamics [18]. He claims that L jk = L kj , meaning that in the vicinity of the equilibrium state, the linear dependence of any flow J j on any force X k coincides with a similar dependence of flow J k on force X j . Onsager reciprocity ratios make it possible to significantly reduce the number of independent coefficients in systems of type (4). The validity of Onsager reciprocal relations is not limited to quasi-equilibrium systems (the reciprocity relations relate to the first derivatives of (5) at the equilibrium point), but their strength is higher the closer the system is to equilibrium. For significantly non-equilibrium systems, linear dependences (4) and reciprocity relations are often verified due to the validity of the LTE principle for individual small parts into which each large system is divided, which is under the influence of gradients of thermodynamic forces [18]. According to the second law of thermodynamics (1), limitations are imposed on the permissible values of kinetic coefficients—the matrix of coefficients L must be integrally defined, i.e.,
( L io ) 2 L ii L oo .
If mass and energy are transferred inside an isolated system under the influence of some forces, then the magnitude of this force decreases, and the system approaches a state of equilibrium, at which its entropy has a maximum value. The entropy of the system in any of its states is a measure of the system’s approach to equilibrium, and the rate of its change over time (1), or the flow of entropy J S , gives a quantitative assessment of the results of the transfer, which takes place under the influence of some driving forces. The smaller these forces are, the smaller the imbalance and the smaller the energy dissipation (2) and energy losses. For non-isolated systems, according to I. Prigozhin’s minimum principle, the flow of entropy will acquire its smallest value in stationary processes, and it is clear that it will be smaller the slower the speed of this process [20].

2.2. A Universal Method of Describing and Evaluating Performance Indicators of Linear ECs

The above-described provisions of linear non-equilibrium thermodynamics are the basis of the theory of the thermodynamic analysis of ECs [2,18].
From the analysis of system (4), it is clear that the more negative the value of L io (compared to the values of the remaining L coefficients), the higher the degree of influence of the input force on the output flow. In addition, with more negative values L io , the degree of suppression of the input flow by the force at the output increases. Thus, the coefficient L io is related to the degree of coupling q between the input and output processes in EC. To avoid the effect on q values of other kinetic coefficients, the degree of coupling was defined as the coefficient L io normalized to the coefficients L ii and L oo [18]:
q = L io L ii L oo .
According to this definition, q will be negative and will acquire, taking into account (6), values from −1 to 0.
To obtain universal EC characteristics in relative units, two more indicators are introduced [18]:
  • relation of forces
χ = X o / X i ;
  • phenomenological relationship
Z = L oo / L ii .
Having normalized the input and output fluxes to the input flux at zero output force (short circuit), the system of Equation (4) can be rewritten using the introduced notations:
{ J i L ii X i = q ( Z χ ) + 1 J o Z L ii X i = ( Z χ ) + q
and the flow ratio ca be determined as follows:
j = J o J i = Z ( Z χ ) + q q ( Z χ ) + 1 .
In expressions (10) and (11), the product ( Z χ ) is a dimensionless number that expresses the reduced force ratio. The content of q also becomes clear: when q = 1 the output flow differs from the input flow always by Z times; therefore, the flows are rigidly (completely) coupled. When q = 0 the input and output flows are caused only by their own forces and flow in the direction of their decrease, which means that the flows are not coupled at all. For values of q from –1 to 0 the output flow is supported by the input flow in the “downhill” direction, but the ratio of these flows changes along with the change in ( Z χ ) , which means that the input and output processes in this case are partially coupled.
An important characteristic of ECs is the thermodynamic efficiency of their operation, i.e., the ratio of the rate at which the EC produces free output energy to the rate at which it consumes free input energy:
η = J o X o J i X i = j χ = ( Z χ ) ( Z χ ) + q q ( Z χ ) + 1 .
Graphical dependence (12) at different values of q is shown in Figure 1a. From (12), it is easy to find the optimal ratio of forces from the point of view of thermodynamic efficiency:
( Z χ ) opt η = q 1 + 1 q 2 ,
at which the thermodynamic efficiency has a maximum value equal to
η max = ( Z χ ) opt η 2 .
The operation of EC at the point with maximum thermodynamic efficiency is attractive at first glance, but at high degrees of coupling ( | q | 1 ), it may be impractical, since the output flow approaches zero, meaning the output process proceeds at a very low speed. To achieve energy conversion at a higher speed, it is necessary to partially sacrifice thermodynamic efficiency.
One of the expedient options for choosing the EC operating point can be the maximum output power mode. The normalized value of the output power can be expressed as follows:
p = J o X o L ii X i 2 = [ ( Z χ ) + q ] ( Z χ ) .
Graphic dependence (15) at different q value is shown in Figure 1b. The maximum output power is reached at
( Z χ ) P . max = q / 2
and is
p max = q 2 / 4 .
However, as can be seen from Figure 1a,b, the maximum output power is accompanied by a rather low thermodynamic efficiency and can be recommended only for ECs connected to free energy sources, for example, natural renewable ones.
Using the obtained universal dependencies (10)–(12) and (15), combining two indicators, it is possible to easily form other optimal operation criteria appropriate for specific ECs [18], for example maximum output flow at optimal thermodynamic efficiency, maximum output power at optimal thermodynamic efficiency, maximum product [ J o / ( Z L ii X i ) ] η (so-called “economically profitable output flow”), maximum product p η (so-called “economically profitable output power”), etc. The technique of ensuring the indicated optimal mode consists of changing the available kinetic coefficients of the original model (4) in such a way to obtain the required combination of parameters ( Z χ ) and q, at which the selected criterion reaches its maximum value.

2.3. Performance Indicators of Cascaded Linear ECs

The analysis shows that elementary ECs in energy conversion systems are most often connected in cascade. For a mathematical description of the steady-state mode of operation at a given point of a system of two cascaded ECs (marked by indices 1 and 2), we will use the following systems of linear equations:
{ J 1 i = L 1 ii X i + L 1 io X o J 1 o = L 1 oi X i + L 1 oo X o ;   { J 2 i = L 2 ii X i + L 2 io X o J 2 o = L 2 oi X i + L 2 oo X o .
For the cascade connection of such ECs, the following conditions apply: X 1 o = X 2 i , J 1 o = J 2 i . If such a cascade is considered a new equivalent EC, then after mathematical transformations, its parameters are found, which are determined based on the parameters of the ECs, from which it is composed:
q Σ = q 1 q 2 ( 1 q 1 2 + γ ) ( 1 q 2 2 + γ 1 ) ;   Z Σ = Z 1 Z 2 γ 1 q 2 2 + γ 1 1 q 1 2 + γ ;   χ Σ = χ 1 χ 2 ,
where γ is the connection coefficient of two ECs in the cascade, which is equal to
γ = L 2 ii L 1 oo .
Figure 2 shows the graphical dependence q ( γ ) at the given q 1 = 0.95 and q 2 = 0.86...0.98 (for a better understanding, the abscissa axis is presented using a logarithmic scale). As can be seen, for any q 1 and q 2 , there exists some optimal connection coefficient at which the total degree of conjugation is maximal. This optimal value of the connection coefficient is obtained in the form
γ opt = 1 q 1 2 1 q 2 2 .
With an optimal connection, the maximum conjugation of two ECs is characterized by the following degree of coupling:
q Σ . max = q 1 q 2 ( 1 q 1 2 ) ( 1 q 2 2 ) + 1
and the optimal reduced force ratio is equal to
( Z χ ) Σ . opt = ( Z χ ) 1 ( Z χ ) 2 .

2.4. The Energy Optimization Method of the Stady-State Operation Mode of the System

Based on the above, it is possible to propose a method of energy optimization of the steady-state operation of the system under study, which consists of the five steps, as is shown in Figure 3.

3. Thermodynamic Analysis of Energy Conversion Efficiency in WT

3.1. Mathematical Description of the Aeromechanical Regularities of WT Operation as a Linear EC

The mechanical power on the axis of WT depends on wind velocity Vw and is determined using the following equation [21]:
P WT = 0.5 ρ S C P ( λ ) V w 3 ,
where ρ is the air density, S is the swept area of WT, CP(λ) is the wind power conversion efficiency factor of WT, λ is the tip speed ratio (TSR), ω is the turbine angular velocity, and R is the radius of WT.
TSR is the linear speed of the blade edge relative to the wind speed:
λ = ω R V w .
The dimensionless aeromechanical characteristic CP(λ) uniquely determines the appearance of the characteristics PWT(ω,Vw), which are scaled along the PWT and ω axes depending on the power (dimensions) of the WT. The CP(λ) characteristic, which is obtained through complex experimental or model studies, is often represented based on a polynomial dependence of the n-th degree:
C P ( λ ) = i = 0 n a i λ i .
To provide maximum power extraction from the wind, according to Equation (24), it is necessary to maintain the maximum value of the power factor CPmax(λopt) and therefore an optimal WT angular velocity,
ω opt = λ opt R V w ,
which is achieved by the automatic regulation of the load torque at a WT shaft.
In the steady state, the torque of a WT is balanced by the load torque, which, based on (24) and (25), is described by the following expression:
T WT = P WT ω = 0.5 ρ S R C P ( λ ) λ V w 2 .
As an EC, the WT converts the power of the wind flow into the power of rotational motion on the WT shaft. At the input of this EC, there is a force in the form of the speed of the wind flow X i = V w , which determines the input flow in the form of a force acting on the WT J i = F . At the output of the EC, the force in the form of angular velocity determines the torque on the WT shaft: X o = ω , J o = T WT . Powers at the input and output of EC are respectively equal to P i = V w F and P o = ω T WT .
Similarly to (4), the system of linear equations that describes the operation of the WT as an EC will have the form
{ F = L ii V w + L io ω T WT = L oi V w + L oo ω .
As can be seen from the mathematical description (24)–(28), the mathematical model of the processes in WT is significantly nonlinear; therefore, the system of linear Equation (29) is valid only in the vicinity of some given operating point. Figure 3 shows a typical dependence T WT ( ω ) with the point of optimal operation T WT . opt ( ω opt ) , in which the maximum WT power value (maximum power point—MPP) P WT . max is provided. Shown in Figure 4, tangent 1 to this point is described based on the second equation of the system (29), if this system is intended to simulate the operation of the WT in the MPP. Let us determine the expressions for the kinetic coefficients in the system of equations (29) using the example of the WT operation at this point.
The equation of the tangent curve shown in the Figure 4 has the form
T WT = k ω + b ,
where
k = d T WT d ω | ω = ω opt = 0.5 ρ S R V w 2 d d ω ( C P ( λ ) λ ) | ω = ω opt = 0.5 ρ S R V w 2 [ d d λ ( C P ( λ ) λ ) d λ d ω ] | ω = ω opt ,
b = T WT . opt + k ω opt .
Considering (25) and introducing the notation B = 0.5 ρ S R 2 V w , the Equation (31) will have the form
k = B d d λ ( C P ( λ ) λ ) | λ = λ opt .
For ω = 0, from the second equation of system (29), we obtain
L oi = T WT | ω = 0 V w = b V w .
For T WT = 0 , ω = b / k follows from Equation (30). Then, from the second equation of system (29), we obtain
L oo = L oi V w ω = b ω = k .
For the estimated optimal operating point from the first equation of the system (29) we obtain
L ii = F L io ω opt V w ,
where F = 0.5 ρ S V w 2 is the value of the force of the wind flow pressure on the WT depending on the wind speed V w .

3.2. The Main Performance Indicators of the Experimental WT

The aeromechanical characteristic of a VAWT taken for the study is shown in Figure 5 and is described based on the following polynomial
C P ( λ ) = i = 0 n a i λ i = 0 . 046977 0 . 128468 λ + 0 . 19596 λ 2 0 . 05705 λ 3 + 0 . 00621 λ 4 0 . 000236 λ 5 .
The main parameters of the experimental VAWT are given in Table 1.
Figure 6 shows the dependences of the torque of the experimental VAWT on its angular velocity for three values of wind speed—8 m/s, 10 m/s, and 12 m/s. Similarly to the one shown in Figure 4, the linearized dependences T WT ( ω ) are shown by lines 1, 2, and 3, respectively, at the optimal points of the maximum VAWT power T WT . opt ( ω opt ) .
As can be seen from Figure 6, near the selected points, the obtained straight lines have small deviations from the linearized curves, and the parameters k and b of these straight lines are different. This is confirmed based on the results of calculations based on Expressions (33) and (32). However, regardless of this, calculations according to Expressions (34)–(36) lead to the same results regarding the main parameters of the VAWT as a linear EC for linearized characteristics at the points of optimal operation of the experimental WT at different wind speeds: degree of coupling q = 0.9192 , reduced force ratio ( Z χ ) = 0 , 4602 . The dependences of the thermodynamic efficiency and the normalized value of the output power, obtained from Expressions (12) and (15), respectively, for these parameters are shown in Figure 7. As can be seen from the figure, the value ( Z χ ) = 0.4602 is optimal for the normalized value of the output power, at which its maximum value p max = 0.2112 is reached according to (17). This corresponds to the optimal operation mode of the VAWT, designed to obtain the maximum power value from the wind at any wind speed (in the range of incomplete VAWT power—from the starting to the nominal wind speed). At the same time, given the obtained value ( Z χ ) = 0 , 4602 , the thermodynamic efficiency of converting wind power into mechanical power of the WT rotor is not maximal.
Figure 8 shows the dependences of the absolute power values of the studied VAWT on its angular velocity (curves 1), obtained based on (15), for three values of the wind speed according to the following expression:
P WT . opt ( ω ) = L ii ( V w , ω opt ) F ( V w ) 2 [ ( Z opt χ ( ω ) ) + q opt ] ( Z opt χ ( ω ) ) ,
where the index opt indicates the corresponding parameters determined for the points of optimal operation of the experimental VAWT, at which the maximum value of the output power is ensured at a specific wind speed.
As expected, at optimal angular velocity values, the output power of the studied VAWT reaches its maximum value. However, as can be seen from Figure 4, when the angular velocity deviates from the optimum, the real WT torque decreases compared to the linearized values corresponding to straight line 1. This means, for deviations from the point of optimal operation, the accuracy of VAWT modeling based on the linear EC decreases. In the future, after consideration of the entire WECS, optimization of its operation will be ensured based on the system criterion—the maximum electrical power received at the output of the generator. This will mostly lead to some deviation of the operating point of the VAWT from its optimal mode. Therefore, in order to increase the accuracy of the mathematical description, it is advisable to model the VAWT, not based on one linear EC, as was performed above, but based on a whole series of ECs with slightly different parameters obtained for the current VAWT angular velocity by linearizing the characteristic T WT ( ω ) at this point. For this purpose, a special program has been developed in the MathCad environment, which determines the dependence of all parameters and variables on the angular velocity of the VAWT in the given vicinity of its MPP, i.e., performs multiple linearization of the characteristics T WT ( ω ) for the given wind speeds of 8 m/s, 10 m/s, and 12 m/s. As a result of such modeling, the kinetic coefficients of the system of Equation (29), determined for each wind speed, will also depend on the VAWT angular velocity. Consequently, this will also lead to the dependence of the main VAWT parameters, such as EC, on the angular velocity—degree of coupling qWT and normalized force ratio ( Z χ ) WT . As a result of the conducted research, the following dependencies of the specified parameters for the studied VAWT at the three specified wind speeds are presented in Figure 9. As can be seen from the given dependencies, the trends of their change are in contradiction: with an increase in angular velocity, qWT decreases by a modulus, but ( Z χ ) WT increases and, conversely, with a decrease in angular velocity, qWT increases by a modulus, but ( Z χ ) WT decreases. In both cases, this will lead to a deviation from the points of optimal operation marked in the figures. As proof, numbers 2 in Figure 8 indicate the dependences of the output power of the studied VAWT on its angular velocity at different wind speeds, which are similar to those indicated by numbers 1 and obtained only for the MPP operation of the studied VAWT. As can be seen, curve 2 is characterized by a sharper decrease in VAWT output power than curve 1 for deviations from the MPP.

4. Thermodynamic Analysis of the PMSG Energy Conversion Efficiency

4.1. Mathematical Description of Mechano-Electrical Regularities of PMSG Operation as a Linear EC

Taking into account losses in both copper and steel, the mathematical model of the PMSG operation in a d-q orthogonal reference frame, rotating with the rotor and oriented along the rotor field, can be obtained from the corresponding two-phase substitute circuit, which is shown in Figure 10 [22]. The diagram shows (the indices d and q refer to the corresponding axes of the coordinate system) the following: u d ,   u q are the projections of the armature voltage vector at the generator output; i d ,   i q are the projections of the armature current vector at the generator output; i d 0 ,   i q 0 are the projections of the armature current vector, which forms the electromagnetic torque of the generator; R a is the active resistance of the armature winding; L d ,   L q are the armature winding inductances relative to the axes d and q; R c is the active resistance simulating losses in generator steel; i dc ,   i qc are the projections of the armature current vector, which simulates losses in the generator steel; ω is the angular speed of the generator rotor; p is the number of generator pole pairs; ψ pm is the flux-coupling amplitude created by one pair of permanent magnet poles.
Based on the schemes presented in Figure 10, the equation of the electrical balance of currents and voltages in the d-q projections under the steady-state mode of operation of the PMSG can be described using the following systems of equations [23]:
{ i q = i q 0 + i d 0 p ω L d R c p ω ψ pm R c i d = i d 0 i q 0 p ω L q R c ;
{ u q = i q R a p ω L d i d 0 ( 1 + R a R c ) + p ω ψ pm ( 1 + R a R c ) u d = i d R a + p ω L q i q 0 ( 1 + R a R c ) .
We will assume that the permanent magnets are placed on the surface of the rotor (surface-mounted PMSM—SPMSM), for which L d = L q = L a , which mostly occurs in the real WECS [24]. Then, the electromagnetic torque created by the generator is as follows:
T G = 3 2 p ψ pm i q 0 .
Control of the PMSG load torque is carried out, most often, in a vector way using an active rectifier that is implemented by forming the armature current in phase with the EMF of the generator, i.e., ensuring the condition i d = 0 [25].
As an EC, the PMSG converts the mechanical power of the WT on its shaft into electrical power at the output of the generator armature winding. At the input of this EC, there is a force in the form of the rotor angular velocity X i = ω , which determines the input flow in the form of the rotation WT torque, which in the steady-state, is equal to the braking PMSG electromagnetic torque J i = T G . At the EC output, the force in the form of the armature voltage vector determines the armature current vector: X o = u a , J o = i a . The power at the EC input is P i = ω T G , and the power at its output is the dot product of the armature voltage and current vectors. P o = ( u a i a ) . The last expression, taking into account the vector control strategy i d = 0 , will take the following form:
P o = ( u a i a ) = 3 2 ( u d i d + u q i q ) = ( 3 2 u q ) ( 3 2 i q ) ,
where the factor 3/2 is due to the transformation of the three-phase reference frame to the two-phase one.
Based on (42), the output force and flow of the EC will be as follows: X o = 3 2 u q , J o = 3 2 i q .
Similarly to (4), the system of linear equations that describes the operation of the PMSG as an EC will have the form
{ T G = L ii ω + L io ( 3 2 u q ) ( 3 2 i q ) = L oi ω + L oo ( 3 2 u q ) .
Let us determine the kinetic coefficients in the system of Equation (43) using the PMSG mathematical model (39)–(41).
Given that i d = 0 , from the second equation of system (39), we obtain
i d 0 = i q 0 p ω L a R c .
The analysis of Expression (44) shows that i d 0 has a much smaller value compared to i q 0 due to the small factor near the latter. When substituting i d 0 into the first equation of system (39), this small factor at i q 0 will already be squared, so the second terms on the right side of this equation can be neglected without the risk of losing the accuracy of the mathematical description. Having performed this and substituting i q from the first equation of system (39) into the first equation of system (40), after transformations, we obtain
i q 0 = u q A R a + p ω ψ pm A R a ( 1 + 2 R a R c ) ,
where
A = 1 + ( p ω L a ) 2 R a R c ( 1 + R a R c ) .
Substituting (45) into (41), we obtain the first equation of system (43) in the form
T G = 3 2 ( p ψ pm ) 2 A R a ( 1 + 2 R a R c ) ω 3 2 p ψ pm A R a ( 3 2 u q ) .
Substituting (45) into the first equation of system (39), we obtain
i q = p ψ pm A R a ω [ 1 + R a R c ( 2 A ) ] + 1 A R a u q .
An analysis of the expression in square brackets shows that its value differs from unity by tenths of a percent, so it can be replaced by unity. Multiplying both parts of Equation (48) by 3 / 2 , we obtain the second equation of system (43) in the form
( 3 2 i q ) = 3 2 p ψ pm A R a ω + 1 A R a ( 3 2 u q ) .
From the equations of the system composed of Equations (47) and (49), we obtain the expressions for the kinetic coefficients of the EC, which describe the operation of the vector-controlled PMSG:
L ii = 3 2 ( p ψ pm ) 2 A R a ( 1 + 2 R a R c ) ,   L io = 3 2 p ψ pm A R a ,   L oo = 1 A R a .
Having the expressions for the kinetic coefficients (50), it is possible to obtain from (5)–(7) the expressions for the main dimensionless parameters of this EC:
q G = ( 1 + 2 R a R c ) 0.5 ,   Z G = ( p ψ pm ) 1 [ 3 2 ( 1 + 2 R a R c ) ] 0.5 ,   χ G = 3 2 u q ω .
Using Equations (45) and (51), we finally obtain
( Z χ ) G = ( 1 + 2 R a R c ) 0.5 [ 1 A R a i q 0 p ω ψ pm ( 1 + 2 R a / R c ) ] .

4.2. Determination of the Parameters of the Studied PMSG

For further research, it is necessary to obtain the parameters of the studied PMSG, which will work together with the VAWT, the parameters of which are given in Table 1. Let us find, with a reasonable approximation, expressions for determining the main parameters of the generator, based on its given energy indicators—relative energy losses in copper δ Cu and steel δ Fe .
In the case of a PMSG with its specified nominal input parameters (torque T n and angular velocity ω n ), its parameters, such as p and ψ pm , determine the level of output parameters—nominal values of voltage and armature current, which in this study, are of no fundamental importance. Therefore p and ψ pm can be chosen arbitrarily, within the limits of real values.
To determine the value of the resistance of the generator armature, we will neglect the losses in the steel and take R c = . Then, the relative energy losses in copper in the case of a vector-controlled PMSG in the nominal mode will be as follows:
δ Cu = Δ P Cu P mech + Δ P Cu = 3 / 2 i q 0 . n 2 R a 3 / 2 p ψ pm i q 0 . n ω n + 3 / 2 i q 0 . n 2 R a = ( p ω n ψ pm i q 0 . n R a + 1 ) 1 ,
where Δ P Cu is the absolute power loss in copper, P mech is the mechanical power of PMSG, and i q 0 . n is the nominal value of the projection of the armature current onto the q axis.
From (53) we obtain
R a = p ω n ψ pm i q 0 . n ( δ Cu 1 1 ) .
The inductance of the PMSG armature winding and other parameters affect the magnitude of the shift angle φ between the armature voltage v and current i vectors. From the vector diagram (Figure 11) for the nominal mode of a synchronous machine with magnets placed on the surface of the rotor, taking into account that its armature current i = i q 0 . n , and EMF e = p ω n ψ pm , the following can be written:
tg φ = p ω n L a i q 0 . n p ω n ψ pm + i q 0 . n R a .
From (55), we receive
L a = ( ψ pm i q 0 . n + R a p ω n ) tg φ .
From the scheme shown in Figure 10, the relative power losses in PMSG steel can be represented by the power losses in the resistance R c . For the nominal mode of PMSG operation, the following can be written:
δ Fe = ( p ω n ψ pm ) 2 R c . n P n ,
where R c . n is the resistance R c value at the nominal operation mode of the PMSG, and P n is the rated power of the PMSG.
Since power losses due to hysteresis and eddy currents, which together make up losses in steel, depend in different ways on the angular velocity of the machine, the resistance R c value decreases as the speed decreases. As shown in [26], such a trend can be modeled using the following equation:
R c = R c . n K f K h + 1 K f K h + 1 ω * ,
where K h and K f are the relative values of power losses due to hysteresis and eddy currents for the nominal PMSG mode (according to experimental data [23], we adopt K f / K h = 0.5694).
The PMSG parameters adopted and calculated using expressions (54), (56), and (57) are given in Table 2.

4.3. Thermodynamic Performance Indicators of the Studied PMSG

To evaluate the parameters of the studied PMSG, the following computational experiment was conducted. The generator is set in motion by an external motor with such torque to ensure its given angular velocity regardless of the electrical load of the generator. Based on the obtained expressions (51) and (52), dependencies of the main parameters of the generator as the EC on the given angular velocity were calculated for three such values of the electric load of the generator, which form its electromagnetic torque at the level of 0.25TG.n, 0.5TG.n, and TG.n. According to (41), these torque values are provided based on the following values of the armature current projection iq0: 10.5 A, 21.0 A, and 42.0 A. The calculation results are presented in Figure 12.
As can be seen from Figure 12a, as well as from the expression for qG (51), the degree of coupling depends on the resistance ratio R a / R c and does not depend on the electrical load of the generator. According to (58), as the angular velocity of the latter increases, the resistance Rc increases as well, and therefore, qG increases. According to (51), dependence ( Z χ ) G is more complicated, and Figure 12b reflects this. According to the optimal dependence ( Z χ ) G . opt ( ω ) shown in Figure 12b with a dotted line, the maximum generator efficiency η G ( ω ) possible for a specific value q G ( ω ) is provided and shown in Figure 12c. However, with a change in the electrical load of the generator, as can be seen from Figure 12b, there are deviations ( Z χ ) G in one direction and the other from the optimal value, which, respectively, leads to a decrease in η G ( ω ) , as can be seen from Figure 12c. From this figure, the highest efficiency values are obtained for 0.5TG.n, because for this torque of electrical load of the generator, the values of ( Z χ ) G are as close as possible to the optimal values (Figure 12b). The lowest efficiency values are observed when the optimal values ( Z χ ) G are exceeded, which can be seen from the curve for 0.25TG.n and curves η ( Z χ ) shown in Figure 12c. It is worth noting that the biggest decrease in efficiency occurs due to the reduction in coupling and deviation from the optimal for this degree of the coupling operation point and not due to the amount of load.

5. Thermodynamic Analysis of Energy Conversion Efficiency in the “VAWT—PMSG” Complex

5.1. Thermodynamic Performance Indicators of the Studied PMSG Driven by the Studied VAWT

When the shaft of the studied VAWT (parameters in Table 1) is connected with the studied PMSG (parameters in Table 2), steady states of the WECS at wind speeds of 8 m/s, 10 m/s, and 12 m/s will be ensured by vector control of the generator armature current projections. Their dependences on the joint angular velocity of the VAWT and PMSG complex are shown in Figure 13.
Under such conditions, the dependence q G ( ω ) for the generator remains the same as shown in Figure 12a, and the dependences ( Z χ ) G and η G ( ω ) , unlike those shown in Figure 12b,c, will have the form shown in Figure 14.
As can be seen from the obtained dependences, for every wind speed, there is a significant range of variation in the VAWT angular velocity with the PMSG, in which the generator efficiency is very close to the maximum possible. This is ensured by the fact that there is equality of ( Z χ ) G = ( Z χ ) G . opt in the middle of these ranges (see Figure 14a).

5.2. Thermodynamic Parameters of the WECS as a Cascade EC

The two ECs investigated above (VAWT and PMSG) are connected into one complex of WECSs in a cascade manner, as evidenced by the obtained linear systems of Equations (29) and (43), which describe the established operation mode of these EC components. In the following systems of equations, the necessary conditions for the cascade connection of the two investigated ECs are fulfilled: X 1 o = X 2 i = ω , J 1 o = J 2 i = T WT = T G . The new complex EC obtained by the cascade connection of the two ECs studied above is characterized by several indicators. The dependence of these indicators on the angular velocity ω of the complex is shown in Figure 15.
From the point of view of the degree of coupling of the complex EC, which most affects the efficiency of the complex, the quality of the cascade connection is characterized by the connection coefficient γ. The dependences of γ(ω) on the angular velocity of the complex are calculated according to expression (20) for the three investigated wind speeds and are presented in Figure 15a. The same figure shows the dependences on the angular velocity of the optimal connection coefficients γopt(ω) of EC components calculated using Expression (21) for the same values of the wind speed. They are shown with specific dependences on the angular velocity of the complex of their degrees of coupling, which are shown, respectively, in Figure 9a and Figure 12a. As can be seen from Figure 15a, the corresponding pairs γ(ω) and γopt(ω) for each of the wind speeds differ from each other but gradually approach each other as the angular velocity of the complex increases. This indicates the suboptimality of the cascade connection of the EC components into a complex, as can be seen from the dependences for the studied wind speeds of the degrees of coupling of the complex EC q Σ ( ω ) , calculated according to Expression (19), which are lower than the corresponding maximum possible values q Σ . max ( ω ) calculated according to (22).
In addition to the degree of coupling of the complex EC, the efficiency of its operation is also affected by the operating mode point, which is characterized by the reduced force ratio of the complex ( Z χ ) Σ . Dependencies ( Z χ ) Σ ( ω ) , calculated using expression (19) for the three investigated wind speeds, are shown in Figure 15c. This figure also shows the dependences on the angular velocity for the optimal values of the reduced force ratios of the complex ( Z χ ) Σ . opt ( ω ) corresponding to the optimal operating modes of the complex EC with the optimal connection coefficients γopt(ω) of the two components of the EC calculated using Expression (23) for the same values of the wind speeds. Since a source with free energy—VAWT powered by wind energy—is included at the input of the complex EC, it is advisable to take the maximum output power (15) as a criterion for the efficiency of the complex EC. The optimal reduced force ratios for such an EC are calculated according to Expression (16) and are presented in Figure 15c based on the curves ( Z χ ) Σ . P max ( ω ) . At the intersection points of the corresponding curves ( Z χ ) Σ ( ω ) and ( Z χ ) Σ . opt ( ω ) with ( Z χ ) Σ . P max ( ω ) , the optimal values of the angular velocity ωopt of the complex are obtained for the real and ideal connection of the two EC components, respectively.
Figure 15d shows the dependences of the wind power conversion efficiency η Σ ( ω ) of the complex EC—the studied WECS—on the angular velocity of the complex at different wind speeds, calculated according to Expression (12) for the real obtained dependences q Σ ( ω ) and ( Z χ ) Σ ( ω ) . Since the power of the wind flow washing the wind rotor is constant for each of the wind speeds, the maximum values of the output power for each of the wind speeds are also provided at the maximum points of the curves η Σ ( ω ) . Also, Figure 15d shows the maximum possible efficiency dependencies on the angular velocity η Σ , max ( ω ) of the complex at the studied wind speeds. They are calculated according to Expression (12) for the maximum values of the degree of coupling q Σ . max ( ω ) and the optimal values of the reduced force ratios ( Z χ ) Σ . opt ( ω ) of the complex, corresponding to the optimal operating modes of the complex EC with optimal connection coefficients γopt(ω). As can be seen from the obtained dependences, the optimal combination of two components of the EC (VAWT and PMSG) provides a significant efficiency increase, especially at low wind speeds. Therefore, for Vw = 8 m/s, the efficiency can increase by 9.2%.
Figure 15d shows the optimal values of the angular velocity of the complex ωopt for wind speeds of 8 m/s and 12 m/s to 11.35 s−1 and 17.75 s−1, respectively (due to Mathcad limitation of displaying only two markers on one axis). The same values of ωopt are also shown in all other parts of Figure 15a–c. As can be seen from Figure 15c, at ωopt, there is an intersection of the curves ( Z χ ) Σ . opt ( ω ) with the curves ( Z χ ) Σ . P max ( ω ) , which, as mentioned above, corresponds to the ideal connection of the two components of the EC. For the intersection points of real curves ( Z χ ) Σ ( ω ) with curves ( Z χ ) Σ . P max ( ω ) , the obtained values ωopt are slightly higher, but these deviations are small.

5.3. Research on Ways to Improve the WECS Efficiency

The simulation results obtained in Section 4 and Section 5 make it possible to determine the optimal operating points of the studied VAWT and PMSG, respectively, as well as the maximum performance indicators of their operation. The latter is conditioned by the given parameters of these devices: for the VAWT, by the characteristic C P ( λ ) , and for the PMSG, by the relative losses in copper and steel (Table 2). However, having the cascade connection of the studied VAWT and PMSG, as shown by the results obtained in Section 5.2, the maximum values of the efficiency of the studied WECS are obtained as a compromise between the efficiencies of these devices. This is mainly due to the non-ideality of their cascade connection, which can be seen from the dependencies shown in Figure 15. At the same time, as can be seen from Figure 15d, there are reserves for improving the efficiency of the studied WECS. This paragraph explores ways of realizing these reserves.
For the connection coefficient γ, it can be seen from the Expression (20) that its deviation from the optimal value (21) depends on the ratio of its own kinetic coefficients—the output Loo for VAWT and the input Lii for PMSG—and this ratio is far greater than the optimal value γopt (Figure 15a).
As can be seen from (35), the value of Loo for a VAWT with a given characteristic can be changed only by changing the turbine power. However, at the same time, the power of the generator increases respectively, from which the resistance of its armature winding is determined according to (54) and (41). Studies have shown that all of these changes lead to the invariance of the ratio between γ and γopt for the optimal VAWT angular velocities ωopt at different wind speeds; however, ωopt decreases with an increasing VAWT power. Therefore, the reserves for changing the parameters of the VAWT to optimize its connection with the PMSG remain only the choice of a WT with a different characteristic C P ( λ ) , which can be the subject of a separate study.
Regarding the change in Lii for the PMSG, as can be seen from (50), the following parameters of the generator can be variable: p, ψpm, Ra, Rcn тa La. As research has shown, the values of Ra and Rcn at given relative power losses in copper and steel are completely dependent on p and ψpm, which are independent parameters of the PMSG. However, taking into account the dependencies (54) and (57), the changes in wide ranges of p and ψpm have absolutely no effect on the Lii value obtained from (50) for the PMSG, which is easy to show analytically. The value of the inductance of the PMSG armature winding La, according to (56), also depends on the parameters p and ψpm. However, as shown above, this does not lead to a change in Lii. Nevertheless, there is another independent parameter in (56)—the shift angle φ between the voltage and current vectors of the generator armature at its nominal angular velocity and nominal load. The conducted studies showed a significant influence of this parameter on the quality of the connection of the studied VAWT and PMSG, as can be seen in Figure 16. This shows how the maximum values of the overall efficiency of the studied WECS change depending on the set value of the nominal shift angle φ between the voltage and current vectors of the PMSG armature at different wind speeds (corresponding to the maximum of the curves shown in Figure 15d).
An analysis of the results shown in Figure 16 gives the possibility to choose the value φ = 40° as an option for compromise. The dependences of the main indicators of the studied WECS obtained for this angle are shown in Figure 17. Based on Figure 17a, it can be said that at the optimal VAWT angular velocities for each of the studied wind speeds (which are slightly different from those obtained in Figure 15), the values of γ are closer to γopt than in Figure 15a. This ensures that the overall efficiency of the studied WECS is close to its maximum value (Figure 17b). It is almost perfect for Vw = 12 m/s, slightly worse for Vw = 10 m/s, and even worse for Vw = 8 m/s. However, it is a better option for the studied PMSG. An even closer option to the optimal one can be sought by taking a more complex type of PMSG, for example, with hybrid magnetoelectric and electromagnetic excitation [27], in which the parameters can be changed depending on the wind speed, which can be the subject of further research. To design such complex electromechanical systems, it is advisable to use modern approaches to realize digital and smart manufacturing, such as, for example, the digital twin [28].

6. Conclusions

The approach of linear non-equilibrium thermodynamics, in particular the method of the universal description of objects and subsystems as ECs, is promising for the analysis of energy processes and improving the efficiency of various systems, in particular renewable energy systems or energy storage systems. This approach makes it possible to assess the energy quality of subsystems of different natures from one universal point of view without delving into the physical, chemical, or other features of the processes. Having a system of universal ECs, the next step can be taken—to identify the conditions for the best connection of such ECs. In this work, this is shown based on a specific example of applying the approach of linear non-equilibrium thermodynamics to the analysis of the efficiency of energy conversion in the WECS, which consists of two devices connected in cascade—the VAWT and the vector-controlled PMSG. To the best of the author’s knowledge, this is the first time that the VAWT and PMSG have been mathematically characterized as linear ECs, enabling their linearization at operating points within a specified range of common angular velocities. By further applying the universal description of these ECs using a system of dimensionless parameters and characteristics, it became possible to obtain the degrees of coupling between the input and output of these devices, choose the parameters for their optimal operating points based on specific criteria (corresponding to the maximum output power and the maximum efficiency), and evaluate their efficiency indicators. The obtained models of the VAWT and PMSG, as linearized ECs, will be used in the future to study other power-generating complexes.
The method of the universal description of ECs has been further developed to assess the patterns of cascade connections between two ECs. This advancement allows for a determination of the quality indicator of a specific connection and its impact on the performance indicators of the resulting equivalent cascaded EC. An analysis of the connection quality of the studied VAWT and PMSG showed that there are reserves in increasing the efficiency of the studied WECS built on their basis. The conducted research made it possible to identify which of the VAWT and the PMSG parameters have an impact on the energy quality of their connection. In particular, the effective factor affecting the quality of this connection turned out to be the angle of shift between the voltage and current vectors of the generator in the nominal mode. This study made it possible to substantiate the optimal value of this angle, equal to 40°. This value provides the best approximation to the optimal connection of the investigated ECs in the entire range of wind speeds, especially with medium and high winds.

Author Contributions

Conceptualization, I.S. and M.L.; methodology, I.S. and Y.B.; software, I.S. and Y.B.; validation, M.L.; formal analysis, I.S. and M.L.; investigation, I.S. and Y.B.; resources, I.S.; data curation, Y.B.; writing—original draft preparation, I.S. and Y.B.; writing—review and editing, I.S. and M.L.; visualization, Y.B.; supervision, I.S. and M.L.; project administration, I.S. and M.L.; funding acquisition, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Turns, S.R.; Pauley, L.L. Thermodynamics: Concepts and Applications, 2nd ed.; Cambridge University Press: Cambridge, UK, 2020; 854p. [Google Scholar]
  2. Demirel, Y. Non-Equilibrium Thermodynamics: Transport and Rate Processes in Physical, Chemical and Biological Systems, 2nd ed.; Elsevier Science & Technology Books: Amsterdam, The Netherlands, 2007; 754p. [Google Scholar]
  3. Christen, T. Efficiency and Power in Energy Conversion and Storage: Basic Physical Concepts; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2018; 178p. [Google Scholar] [CrossRef]
  4. Kleidon, A. Working at the Limit: A Review of Thermodynamics and Optimality of the Earth System. Earth Syst. Dyn. Discuss. 2022, preprint. [Google Scholar] [CrossRef]
  5. Perera, S.M.H.D.; Putrus, G.; Conlon, M.; Narayana, M.; Sunderland, K. Wind energy harvesting and conversion systems: A technical review. Energies 2022, 15, 9299. [Google Scholar] [CrossRef]
  6. Shchur, I.; Klymko, V.; Xie, S.; Schmidt, D. Design features and numerical investigation of counter-rotating VAWT with co-axial rotors displaced from each other along the axis of rotation. Energies 2023, 16, 4493. [Google Scholar] [CrossRef]
  7. Wu, C.; Wang, Q.; Yuan, R.; Luo, K.; Fan, J. Large eddy simulation of the layout effects on wind farm performance coupling with wind turbine control strategies. ASME. J. Energy Resour. Technol. 2022, 144, 051304. [Google Scholar] [CrossRef]
  8. Wang, Q.; Luo, K.; Wu, C.; Tan, J.; He, R.; Ye, S.; Fan, J. Inter-farm cluster interaction of the operational and planned offshore wind power base. J. Clean. Prod. 2023, 396, 136529. [Google Scholar] [CrossRef]
  9. Zhang, S.; Luo, K.; Yuan, R.; Wang, Q.; Wang, J.; Zhang, L.; Fan, J. Influences of operating parameters on the aerodynamics and aeroacoustics of a horizontal-axis wind turbine. Energy 2018, 160, 597–611. [Google Scholar] [CrossRef]
  10. Sahin, A.D.; Dincer, I.; Rosen, M.A. Thermodynamic analysis of wind energy. Int. J. Energy Res. 2006, 30, 553–566. [Google Scholar] [CrossRef]
  11. Hu, W.; Liu, Z.; Tan, J. Thermodynamic analysis of wind energy systems. In Wind Solar Hybrid Renewable Energy System; Kenneth, O.E., Tahour, A., Aissaou, A.G., Eds.; IntechOpen: London, UK, 2019; pp. 1–19. [Google Scholar] [CrossRef] [Green Version]
  12. Diyoke, C.; Ngwaka, U. Thermodynamic analysis of a hybrid wind turbine and biomass gasifier for energy supply in a rural off-grid region of Nigeria. Energy Sources Part A Recovery Util. Environ. Eff. 2021, 1–26. [Google Scholar] [CrossRef]
  13. Zhang, Y.; Yang, K.; Li, X.; Xu, J. Thermodynamic analysis of energy conversion and transfer in hybrid system consisting of wind turbine and advanced adiabatic compressed air energy storage. Energy 2014, 77, 460–477. [Google Scholar] [CrossRef]
  14. Koschwitz, P.; Bellotti, D.; Sanz, M.C.; Alcaide-Moreno, A.; Liang, C.; Epple, B. Dynamic parameter simulations for a novel small-scale power-to-ammoniac. Processes 2023, 11, 680. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Lin, Y.; Lin, F.; Yang, K. Thermodynamic analysis of a novel combined cooling, heating, and power system consisting of wind energy and transcritical compressed CO2 energy storage. Energy Convers. Manag. 2022, 260, 115609. [Google Scholar] [CrossRef]
  16. Ehyaei, M.A.; Baloochzadeh, S.; Ahmadi, A.; Abanades, S. Energy, exergy, economic, exergoenvironmental, and environmental analyses of a multigeneration system to produce electricity, cooling, potable water, hydrogen and sodium-hypochlorite. Desalination 2021, 501, 114902. [Google Scholar] [CrossRef]
  17. Vosough, A.; Vosough, S. Different kind of renewable energy and exergy concept. Int. J. Multidiscip. Sci. Eng. 2011, 2, 1–7. [Google Scholar]
  18. Westerhoff, H.V.; van Dam, K. Thermodynamics and Control of Biological Free-Energy Transduction; Elsevier: Amsterdam, The Netherlands, 1987; 298p. [Google Scholar] [CrossRef]
  19. Shchur, I.; Rusek, A.; Lis, M. Optimal frequency control of the induction electric drive based on the thermodynamics of irreversible processes. Electrotech. Comput. Syst. 2011, 3, 377–380. Available online: https://eltecs.op.edu.ua/index.php/journal/article/view/751 (accessed on 17 February 2019).
  20. Kondepudi, D.; Prigogine, I. Modern Thermodynamics: From Heat Engines to Dissipative Structures, 2nd ed.; John Wiley & Sons: New York, NY, USA, 2015; 560p. [Google Scholar] [CrossRef]
  21. Li, Y.; Yang, S.; Feng, F.; Tagawa, K. A review on numerical simulation based on CFD technology of aerodynamic characteristics of straight-bladed vertical axis wind turbines. Energy Rep. 2023, 9, 4360–4379. [Google Scholar] [CrossRef]
  22. Krause, P.; Wasynczuk, O.; Sudhoff, S.; Pekarek, S. Analysis of Electric Machinery and Drive Systems, 3rd ed.; Wiley: Hoboken, NJ, USA, 2013; 688p. [Google Scholar] [CrossRef]
  23. Majout, B.; El Alami, H.; Salime, H.; Zine Laabidine, N.; El Mourabit, Y.; Motahhir, S.; Bouderbala, M.; Karim, M.; Bossoufi, B. A review on popular control applications in wind energy conversion system based on permanent magnet generator PMSG. Energies 2022, 15, 6238. [Google Scholar] [CrossRef]
  24. Shchur, I.; Rusek, A.; Biletskyi, Y. Energy-shaping optimal load control of PMSG in a stand-alone wind turbine as a port-controlled Hamiltonian system. Przegląd Elektrotechniczny (Electr. Rev.) 2014, 5, 50–55. [Google Scholar] [CrossRef]
  25. Torki, W.; Grouz, F.; Sbita, L. Vector control of a PMSG direct-drive wind turbine. In Proceedings of the 2017 International Conference on Green Energy Conversion Systems (GECS), Hammamet, Tunisia, 23–25 March 2017; pp. 1–6. [Google Scholar] [CrossRef]
  26. Fernandez-Bernal, F.; García-Cerrada, A.; Faure, R. Determination of parameters in interior permanent-magnet synchronous motors with iron losses without torque measurement. IEEE Trans. Ind. Appl. 2001, 37, 1265–1272. [Google Scholar] [CrossRef]
  27. Golovko, V.M.; Ostroverkhov, M.Y.; Kovalenko, M.A.; Kovalenko, I.Y.; Tsyplenkov, D.V. Mathematical simulation of autonomous wind electric installation with magnetoelectric generator. Nauk. Visnyk Natsionalnoho Hirnychoho Universytetu 2022, 5, 74–79. [Google Scholar] [CrossRef]
  28. Liu, L.; Guo, Y.; Yin, W.; Lei, G.; Zhu, J. Design and optimization technologies of permanent magnet machines and drive systems based on digital twin model. Energies 2022, 15, 6186. [Google Scholar] [CrossRef]
Figure 1. Dependencies of efficiency (a) and normalized output power (b) on the reduced force ratio at different degrees of coupling q.
Figure 1. Dependencies of efficiency (a) and normalized output power (b) on the reduced force ratio at different degrees of coupling q.
Energies 16 05234 g001
Figure 2. Dependencies of the degree of coupling for an equivalent EC on the connection coefficient of two ECs at given values of the degrees of coupling of EC 1 ( q 1 = 0.95 ) and EC 2 ( q 2 = 0.86...0.98 ).
Figure 2. Dependencies of the degree of coupling for an equivalent EC on the connection coefficient of two ECs at given values of the degrees of coupling of EC 1 ( q 1 = 0.95 ) and EC 2 ( q 2 = 0.86...0.98 ).
Energies 16 05234 g002
Figure 3. Flow chart of the application of the proposed approach.
Figure 3. Flow chart of the application of the proposed approach.
Energies 16 05234 g003
Figure 4. View of the T WT ( ω ) characteristic with its linearization at the point of optimal WT operation.
Figure 4. View of the T WT ( ω ) characteristic with its linearization at the point of optimal WT operation.
Energies 16 05234 g004
Figure 5. Power coefficient vs. TSR for the studied VAWT.
Figure 5. Power coefficient vs. TSR for the studied VAWT.
Energies 16 05234 g005
Figure 6. VAWT torque vs. its angular velocity at different wind speeds.
Figure 6. VAWT torque vs. its angular velocity at different wind speeds.
Energies 16 05234 g006
Figure 7. Dependencies of p ( Z χ ) and η ( Z χ ) , obtained for the experimental VAWT as an EC linearized at the point of optimal operation of the VAWT.
Figure 7. Dependencies of p ( Z χ ) and η ( Z χ ) , obtained for the experimental VAWT as an EC linearized at the point of optimal operation of the VAWT.
Energies 16 05234 g007
Figure 8. Dependences of the output power of the studied VAWT on its angular velocity at different wind speeds.
Figure 8. Dependences of the output power of the studied VAWT on its angular velocity at different wind speeds.
Energies 16 05234 g008
Figure 9. Dependencies of the main parameters of ECs that simulate the operation of the studied VAWT at different wind speeds on the angular velocity of the VAWT: (a) the degree of coupling q, (b) the reduced force ratio ( Z χ ) .
Figure 9. Dependencies of the main parameters of ECs that simulate the operation of the studied VAWT at different wind speeds on the angular velocity of the VAWT: (a) the degree of coupling q, (b) the reduced force ratio ( Z χ ) .
Energies 16 05234 g009
Figure 10. PMSG replacement scheme taking into account losses in steel.
Figure 10. PMSG replacement scheme taking into account losses in steel.
Energies 16 05234 g010
Figure 11. Vector diagram of a vector-controlled synchronous machine with permanent magnets placed on the surface of the rotor.
Figure 11. Vector diagram of a vector-controlled synchronous machine with permanent magnets placed on the surface of the rotor.
Energies 16 05234 g011
Figure 12. Dependences of the main parameters of the PMSG as the EC on its angular velocity: (a) degree of coupling, (b) reduced force ratio, (c) efficiency.
Figure 12. Dependences of the main parameters of the PMSG as the EC on its angular velocity: (a) degree of coupling, (b) reduced force ratio, (c) efficiency.
Energies 16 05234 g012
Figure 13. Dependencies of the PMSG armature current projections on the angular velocity of the VAWT with a generator at different values of the wind speed.
Figure 13. Dependencies of the PMSG armature current projections on the angular velocity of the VAWT with a generator at different values of the wind speed.
Energies 16 05234 g013
Figure 14. Dependencies of the reduced force ratio (a) and efficiency (b) of the PMSG driven by the studied VAWT on the angular velocity.
Figure 14. Dependencies of the reduced force ratio (a) and efficiency (b) of the PMSG driven by the studied VAWT on the angular velocity.
Energies 16 05234 g014
Figure 15. Dependencies of the main indicators of the experimental WECS as a cascade EC on the angular velocity of the VAWT with PMSG at different wind speeds: (a) the connection coefficient between two ECs, (b) total degree of coupling, (c) reduced force ratio, (d) efficiency.
Figure 15. Dependencies of the main indicators of the experimental WECS as a cascade EC on the angular velocity of the VAWT with PMSG at different wind speeds: (a) the connection coefficient between two ECs, (b) total degree of coupling, (c) reduced force ratio, (d) efficiency.
Energies 16 05234 g015
Figure 16. Dependencies of the overall efficiency of the studied WECS on the given value of the nominal shift angle φ between the voltage and current vectors of the PMSG armature at different wind speeds.
Figure 16. Dependencies of the overall efficiency of the studied WECS on the given value of the nominal shift angle φ between the voltage and current vectors of the PMSG armature at different wind speeds.
Energies 16 05234 g016
Figure 17. Dependencies of the main indicators of the studied WECS on the angular velocity of the VAWT—PMSG complex with the optimally chosen nominal value of the angle φ at different wind speeds: (a) the connection coefficient between two ECs, (b) efficiency.
Figure 17. Dependencies of the main indicators of the studied WECS on the angular velocity of the VAWT—PMSG complex with the optimally chosen nominal value of the angle φ at different wind speeds: (a) the connection coefficient between two ECs, (b) efficiency.
Energies 16 05234 g017
Table 1. Parameters of the VAWT.
Table 1. Parameters of the VAWT.
ParametersValue
Rated power, PWT.n (kW)10
Rated wind speed, Vw (m/s)12
Maximum of power coefficient, Cp.max0.3661
Optimum value of TSR, λopt3.873
Rotor radius, R (m)2.65
Rotor high, H (m)1.25
Air density, ρ (kg/m3)4.78
Table 2. Parameters of the studied PMSG.
Table 2. Parameters of the studied PMSG.
ParametersValue
Rated power, PG.n (kW) 10
Rated angular velocity, ωG.n (s−1) 17.54
Rated torque, TG.n (Nm) 570.2
Number of pole pairs, p 24
PM flux linkage, ψpm (Wb) 0.41
Relative losses in copper, δCu 0.07
Relative losses in steel, δFe 0.03
Angle of shift between armature voltage and current, φ (deg) 30
Winding resistance, Ra (Ω) 0.286
Winding inductance, La (mH) 6.1
Rated equivalent iron loss resistance, Rc.n (Ω) 84.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shchur, I.; Lis, M.; Biletskyi, Y. A Non-Equilibrium Thermodynamic Approach for Analysis of Power Conversion Efficiency in the Wind Energy System. Energies 2023, 16, 5234. https://doi.org/10.3390/en16135234

AMA Style

Shchur I, Lis M, Biletskyi Y. A Non-Equilibrium Thermodynamic Approach for Analysis of Power Conversion Efficiency in the Wind Energy System. Energies. 2023; 16(13):5234. https://doi.org/10.3390/en16135234

Chicago/Turabian Style

Shchur, Ihor, Marek Lis, and Yurii Biletskyi. 2023. "A Non-Equilibrium Thermodynamic Approach for Analysis of Power Conversion Efficiency in the Wind Energy System" Energies 16, no. 13: 5234. https://doi.org/10.3390/en16135234

APA Style

Shchur, I., Lis, M., & Biletskyi, Y. (2023). A Non-Equilibrium Thermodynamic Approach for Analysis of Power Conversion Efficiency in the Wind Energy System. Energies, 16(13), 5234. https://doi.org/10.3390/en16135234

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop