Open access peer-reviewed chapter

Designing Mu Robust Controller in Wind Turbine in Cold Weather Conditions

By Tahere Pourseif, Majid Taheri Andani, Hamed Pourgharibshahi, Hassan Zeynali and Arash Shams

Submitted: December 12th 2017Reviewed: January 30th 2018Published: April 2nd 2018

DOI: 10.5772/intechopen.74626

Downloaded: 107

Abstract

Due to wind turbine is in class of complex nonlinear system so the precise model of this plant is not accessible, therefore it can be categorized as an uncertain model. So, controlling of this system is a demanding topic. Many of schemes which presented for controlling of wind turbines investigate these systems in a good weather condition. However, many turbines work in severe weather condition. In this study, wind turbine is suggested in cold weather, and in ice on turbine blades which they are considered as uncertainties in the model. A robust controller is designed for the wind turbine, to control the pitch angle.

Keywords

  • ice formation on wind turbine blades
  • robust mu control
  • uncertainty
  • wind turbine

1. Introduction

The growth of the population in world causes to increase of electrical power consumption. Overpopulation also leads to depletion of fossil fuel energy, degradation of environment. These consequences motivate scholars to seeking better ways to produce electrical power [1, 2]. Therefore, alternative energies where their sources are Renewable can be an effective alternative way to produce electricity that generate from natural processes that can be refilled continually and are unpredictable. Therefore, the control of this kind of energy is challengeable. The renewable energy is practical because of its being low-priced, easily accessible, and purity of the energy. Wind turbine is a tool which is operates on a noncomplex principle that transforms the kinetic energy of wind into electrical energy. Various control methods are presented to regulate the wind turbine speed for effective power generation and to maintain the turbine elements within designed speed and torque bounds [1, 3]. In [4] performance of the turbine for working in various operation is investigated and for operating in the third are a controller is designed. In [5], a well-known robust controller which is named as sliding mode controller is utilized to control power as well as to adjust wind’s rotor speed and turbine generator. Due to exist the nonlinear in the characteristics of the controller, just two uncertainties such as spring constant and damping coefficient have been considered and employed to the wind turbine system. The main problem of this scheme is the chattering that produce by Sliding mode controller on control actions in [6], another robust controller is proposed to view all of the uncertainties. These uncertainties comprise the uncertainty is produced from the linearization procedure and minor deviations, and from damping coefficient and spring constant. These approaches are applied to control the pitch angle and enlarge electrical power generation [7]. In many cases, the perfuming of the wind turbine is examined in appropriate weather conditions. But, a robust controller is designed in this paper to keep the performing of the wind turbine in a severe weather conditions, like snow or in a situation that rain droplets freeze on the turbine blades [8]. In these scenarios, when the turbine blades freeze, the mass of the rotor would be changed and the electric power decreases. Then, the controller is given to regulate the speed of the generator after employing the uncertainties in the inertia and other factors [9, 10]. In light of the pitch angle and the speed of the rotor the kind of wind turbine applied in this research has different rotor speed and different pitch angle. The regulation of the changeable pitch angle and the changing rotor speed, respectively, lead to increasing of the electrical power and decreasing of the dynamic load of the turbine. It should be mentioned that that the modification of the variable rotor speed not only minimizes the turbine’s dynamic load, but also increases the system’s lifespan [11, 12]. The aim of designing the controller is to maximize the electrical power production at low wind speed and to keep it at high wind speed. In comparison with other publications that have done on wind turbine, it is the first time which the uncertainty that arises from icing of wind blade is considered, and for this condition then robust controller is proposed.

2. Wind turbine model

2.1. Wind turbine

Effective wind speed is a nonlinear stochastic process that is approximated by a linear model in order to simplify and satisfy the control objectives [9]. There are two terms in the wind model:

V=Vt+VmE1

In this model the wind has two elements, mean value term (Vm) and turbulent term (Vt).

In which:

Vm=1TtT2t+T2vτE2
Vt=kp1s+1p2s+1e;eN10E3

The turbulent term Vtcan be modeled by:

V̇tV¨t=011p1p2p1+p2p1p2VtV̇t+0kp1p2vE4

where, vN10and the parameters p1,p2and kdepend on the mean value of the wind speed (Vm).

2.2. Nonlinear model of wind turbine system

In this section, for deriving the model, the entire wind turbine is separated into four subsystems: aerodynamics subsystem, mechanical subsystem, generator subsystem and pitch actuator subsystem. The aerodynamic subsystem transforms wind forces into mechanical torque and thrust on the rotor. The mechanical subsystem contains of the drivetrain, tower and blades. The drivetrain converts torque of the rotor to the electrical generator. The tower keeps the nacelle and resists the thrust force and the blades transform wind forces into the aerodynamics torque and thrust. The generator subsystem transforms mechanical energy to electrical energy and finally the blade-pitch and generator-torque actuator subsystems are portion of the control model. To model the overall wind turbine, models of these subsystems are attained and at the end they are linked together. A wind model is attained and increased with the wind turbine model to be considered for wind speed estimation. Figure 1 shows the basic subsystems and their interactions.

Figure 1.

Schematic of the wind turbine.

Numerous degrees of freedom can be used to model the plant, but for control design mostly just a less important degrees of freedom are used. In this study we just use two degrees of freedom, namely the rotational degree of freedom (DOF) and drivetrain torsion. The subsystems nonlinear model of wind turbine is used in the following:

2.2.1. Aerodynamics subsystem

The following nonlinear equation is the model of the aerodynamics subsystem [13]:

Pr=12ρπR2V2CpλβE5

Where, Pris the electrical power of rotor, ρ is the density of wind (Kgm3), R is the radius of the blades (m), V is the wind speed (m/s), λis the tip speed ratio, βis the pitch angle, Cpthe power coefficient, λis the ratio between the blade tip speed and the wind speed:

λ=VRωrE6

where the ωris the angular speed of rotor. The derivative of Cpwith respect to both λ and βare therefore used extensively in the control designs. In previous work, numerical derivatives of Cpwith respect to λ and θ have been successfully obtained and utilize in linear control designs. However, the nonlinear controllers have proven sensitive towards the noise introduced by the numerical derivations, making it difficult to validate the responses obtained by e.g. a feedback linearizing controller. Consequently, it has been Chosen to base the aerodynamic model on an analytic expression of Cp. The relation between Cp, β and λ can be written as [13]:

Cpλβ=0.22116λt0.6β5exp12.5λtE7

where,

1λt=1λ1+0.12β0.0351.5β2+1E8

2.2.2. Generator subsystem

The synchronous generator is supposed to be idea, so the power of generator is calculated in [14]:

Pe=TgωgE9

The Pedenoted electrical power of generator, the generator’s angular speed is shown by ωgand the Tgrepresents generator torque which is controllable, although, it is not practical to change immediately. The response of generator’s dynamic is modeled with a first order linear model [6]:

Ṫg=1τTTg+1τTTg,rE10

where Tg,rrepresents the generator torque reference value and the time constant is denoted by τT

2.2.3. Pitch actuator subsystem

The blade’s pitch begins to change by a hydraulic/mechanical actuator to push the angle of blade β to tracks the reference of blade angle βin, The pitch of the blades is changed by a hydraulic/mechanical actuator Eq. (11) is a first order linear model which is a simplified of the dynamic model [6]:

β̇=1τββ+1τββinE11

where τβis the time constant.

2.2.4. Mechanical Subsystem

Figure 2 shows a schematic of the wind turbine mechanics. The turbine is split into two parts, separated by the transmission: The rotor side and the generator side [15].

Figure 2.

Schematic of the wind turbine mechanics.

Jrdenotes the inertia which is on the rotor side and generator side is represented by Jgwhich is on the leftmost and the rightmost disc respectively. The shaft link the rotor to the transmission is subject to huge torques that leads it to twist, so the shaft is suitably modeled as a damped spring. Tr, which denotes rotor torque, excites the model on the left Tgthat represents the generator torque excites the model on the right. The torques and Tsgare the torques at each side of the transmission part and has relation by Ngrepresents the gear ratio which is a relation between Tsrand Tsgthat are named the torque of each side of the transmission part.

Tsg=TsrNgE12

The equations describing the dynamics are obtained using Newton’s second law for rotating bodies. This results in two equations: one for the rotor side and one for the generator side.

ω̇r=TrTsrE13
ω̇g=TgTsgE14

Introducing a variable δ [rad] describing the twist of the shaft, leads to the following equation describing the twist of the flexible shaft [16]:

Tsr=Drδ̇+KrδE15

where:

δ=ΩrΩgNg,δ̇=ωrωgNgE16

In the above equations, Drrepresents the damping and Krdenotes spring coefficient, angular speed of rotor is defined by ωr, the angular speed of generator represent by ωg, Ωrand Ωgare used to define the default shaft angle at the rotor and the default shaft angle at the generator, respectively.

2.3. Linearized models of wind turbine system

As it was discussed in previous section, for design a controller a linear model of the system is needed. The input of this model is wind [17]:

input=VβinTg,rTE17

And outputs of the system include:

output=ωrωgPeE18

Having all the equations, system equations become: [8]:

ω̇r=aDrJrωraDrJrNgωgKrJrδ+a14βE19
ω̇g=DrNgJgωrDrNg2Jgωg+KrJgNgδ1JgTgE20
δ̇=ωrωgNgE21
β̇=1τββ+1τββinE22
Ṫg=1τTTg+1τTTg,rE23
V̇t=VtE24
V¨t=1p1p2Vtp1+p2p1p2V̇t+kp1p2vE25

In which Jrand Jgare rotor and generator moments of inertia, τβand τTare time constants of the first order actuator models

3. Uncertainties model

Various uncertainties have been examined in the current literature. These uncertainties derive from approximated and process parameters in a nonlinear system which changes as the operating point changes, a matter causing the electrical power production to reduce. There are always discrepancies between real system and mathematical models, which lead to uncertain models. In this work, sources of uncertainties are taken to be:

  • Uncertainty in the drivetrain stiffness and damping parameters.

  • Uncertainty in the linearized model.

Although, all control design presented design for moderate temperature. Decreasing the temperature in winter has devastating properties on the wind turbine. Ice on the elements of wind turbine leads some serious problems. Even a few amount of ice on the blades worsens the aerodynamic performance of the wind turbine. It not only decreases the output power energy, but also raises the abrasion between the elements [11, 16]. In other word, the ice in cold places and the high density of air at cold climate have damaging properties on aerodynamics. Fluctuation of produced power and load are reason for such dysfunctionality. Masses of the ice on the turbine cause fluctuation on the frequency of the turbine’s elements and also the behavior of the system’s dynamic [10, 11, 12]. In addition, this condition has effect on control plant. In other word, the performance of the turbine system worsens through wrong data sending [11]. Previous article dealing with this issue have presented approaches like observation, the use of sensors and monitors, considering aerodynamic sound, etc., to recognize ice. The control schemes are then designed to eliminate the ice [9, 10, 17].

3.1. The impact of cold weather on the operation of wind turbine

In this research, for the first time a new approach is advanced to enhance the wind turbine performance in cold climate conditions and to stop the damage which cause about the shutdown of the turbine. Because of the structure of the wind turbine, when frozen blades’ masses changes, the rotor mass will be changed and lead to the inertia of rotor. These variations will effect on equations of the wind turbine and optimal power creation. Therefore, control of turbine and optimal power output could be possible by considering the inertia of the rotor as a new component of uncertainty in the plant. The value of this uncertainty differs with decreasing the temperature and the turbine’s production capacity. In [18], imbalance of the blade was simulated by scaling the density of mass of one blade, which generates an imbalance distribution of mass with respect to the rotor. Furthermore, the aerodynamic asymmetry was simulated by adjusting the pitch of one blade, which produces an imbalance torque across the rotor. In our work, the blade imbalance in blade is due to icing on the blade. Thus, the uneven of the blade is considered about 20% as uncertainty. A robust control is designed to control the turbine in the existence of uncertainty because of the blade imbalance icing [19]. The uncertainties in the wind turbine comprise the linearized model parameters which is extracted from the nonlinear plant, spring constant, and damping coefficient that alteration as the working point diverges; another uncertainty which is added to the system because of the presence of noise and disturbance in the input signal. All of these uncertainties are considered in appropriate weather conditions. In this work, cold climate condition and inertia of rotor are considered as other sources of uncertainty in the wind turbine system [20, 21].

3.2. Frequency response of wind turbine

These icy turbine blades change the rotor mass. Under the frequency response analysis of the system is shown in Figure 3, the rotor inertia uncertainty can be considered between the range of 0 and 20% and The Wind turbine system is stability:

Figure 3.

Bode diagram open loop system with 20% uncertainty of inertia rotor.

According to, the red color is nominal frequency response and blue color is system uncertainties that uncertainty is considered between the range of 0–20%. The system has positive phase margin and positive gain margin, so it is stable.

Figure 4 shows that the system is unstable by considering range of uncertainty between 20 and 30%. So system uncertainties are defined between 0 and 20%.

Figure 4.

Bode diagram open loop system with 30% uncertainty of inertia rotor.

The uncertainties in parameters can be shown as follows:

Kr=K¯r1+pKrδKrE26
Dr=D¯r1+pDrδDrE27
a=a¯1+paδaE28
a14=a¯141+pa14δa14E29
Jr=J¯r1+pJrδJrE30

Two reasons exist which explain that spring coefficient and the damping can be considered as uncertainties. One reason is, there exist divergence in the spring variables from fabrication to fabrication, and manufacturer to manufacturer and other reason, it should be consider that the spring parameter and damping can be changed during a long time because of continuous working, and aging. In previous equations, K¯r, D¯r, a¯, a¯14are nominal parameter values, resulting from the spring constant, damping coefficient, linearization process and rotor inertia, respectively. pKr, pDr, pa, pa14indicate maximum relative uncertainties that are for uncertainty parameters shown in Table 1.

ParameterDescriptionUnitTolerance
Krspring coefficientNm/rad±10%
DrDampingkg.m2rad.s±10%
ALinearization parameter±20%
a14Linearization parameter±20%
Jrrotor inertiakg.m2+20%

Table 1.

Parameter uncertainties of the wind turbine.

δKr, δDr, δa, δa14and δJrare relative changes in these parameters. Therefore:

δKr1,δDr1,δa1,δa141,δJr1E31

4. Robust design controller

4.1. Closed loop system design specifications

Figure 5 shows the block diagram of Wind turbine closed loop system, including the feedback structure, the controller, as well as the model uncertainties and performance objectives weights.

Figure 5.

Block diagram of the closed-loop system with performance specifications [22].

In Figure 5, (r) shows the reference input, (V) represents the wind speed which has disturbance, (n) denotes noise and euand eyare considered as two output costs. The system (R) represents the performance when model is ideal, to which the designed closed-loop system wants to reach [23]. The transfer function of model is chosen in a way that the time response of the reference signal has an overshoot less than around 5%. Inside the dots rectangle is the ideal model, which shows with, Gnomof the wind turbine model and the block that parameterizes uncertainties in the model. To find the wanted performance, inputs r can be obtained from the transfer function matrix, and we need to find disturbance in V and also n to outputs euand ey. Thus the infinity norm of that transfers function can be minor for the entire existing uncertainty variable. The position noise signal is attained by moving the unit-bounded signal which is shows by n through the weighting transfer matrix which denotes Wn. The transfer matrices Wpand Wurepresent the relative significance of the diverse frequency spans for which the performance is needed. So, the performance aim can be reorganize, with probable slight conservativeness, like that transfer function matrix infinity norm be less than 1. So Δ matrix is given in following form:

Δ=diagpKrpDrpapa14pJrE32

This transfer function can be written as [21]:

euey=WpSoGuKRWuSiKWpSoGvWuKSoGvWpSoGuKWnWuKSoWnrVnE33

where Si=I+KG1and So=I+GK1are the input and output sensitivities, repectively. Note that SoGis the transfer function between V and y.

4.1.1. Robust stability

By definition, the closed loop system achieves robust stability if the closed loop system is internally stable for each possible plant dynamics G=FuGnom.

4.1.2. Robust performance

The closed loop system must remain internally stable for each G=FuGnomand in addition the performance criterion should be satisfied for each G=FuGnom[21].

4.2. Matching transfer function and weighting transfer functions

In the case of mu controller optimization design, we have to define the model transfer function which is denote with R and the weighting transfer functions that are nominated with Wn, Wpand Wu.

The model transfer function is selected therefore the time response to the reference signal has an overshoot fewer than 50% and a settling time not more than 1 ms. A probable plant which please the requirements is:

Rs=10.48s2+0.95s+1E34

In Figure 6, the response of matching model to power electrical input is shown.

Figure 6.

Response of matching model to electrical power input.

The noise shaping function Wnis determined on the basis of the spectral density of the position noise signal. In the given case it is taken as the high-pass filter. In this case output has a noteworthy spectral content more than 500 Hz. For this type of filter, the position noise signal is just 0.95 V in the low-frequency values but it is 1 V in the high-frequency values that matches to a position error of around 5% without width.

Wns=95×1050.1s+10.001s+1E35

So the frequency response is shown in Figure 7.

Figure 7.

Frequency response of noise filter Wn.

The closed-loop system performance specifications are reflected by the weighting performance function Wp.

Wps=0.952s+502s+0.005E36

Figure 8 shows the frequency response of the inverses of this weighting function.

Figure 8.

Frequency response of weighting function 1Wp.

It is shown that in a selection, the objective is to obtain a minor variance between the system and outputs of the model, and a minor effect of the disturbance on outputs of the system. This will ensure nice tracking of the reference input and minor error because of the low-frequency disturbances. The weight of control function is usually selected as high-pass filters to make sure that the control action will not surpass 25°.

Wus=0.0220.8s2+10s+23.5×104s2+20×102s+2E37

Figure 9 shows the frequency response of this weighting function Wu.

Figure 9.

Frequency response of weighting function Wu.

According to the above figure, effort control is very low in low- frequencies that cause reduced the control cost.

4.3. Robust μ controller design

The objective of controller with applying the μ synthesis scheme is the stabilizing of the closed loop of the plant and pleasing all of the control demands. In the existence of measurement noise, disturbance in the wind, and uncertainties, the closed-loop plant should have the robust efficiency. Therefore, the purpose of μ controller scheme is design a controller where the wind turbine can track the reference input of generated electrical energy when noise and disturbance are exist. In Figure 10, the schematic of the closed loop model which is utilized with the defined uncertainties is depicted for μ controller model design.

Figure 10.

Closed loop system model for design μ robust controller.

From the figure, it can be seen that ycis defined as the difference between noisy output and the reference and the transfer function P(s) shows the open-loop transfer function matrix with 10 inputs and eight outputs. The upper linear fractional transformation (LFT) of the closed loop system is:

P=FuPnomΔrE38

where Pnomis the nominal transfer function matrix, rincludes five uncertainties in the wind turbine model. We assume Pis defined the structure of uncertainties block as follows:

Δp=Δr00ΔfΔr5×5ΔfC4×2E39

Controller, which is attained with this scheme, is typically of high order controller that cause to challenges in a real-world implementation. So for this purpose, it is recommended to decrease the order of the control plant until it is feasible to simplify the closed loop scheme theory and operation. The results of the singular parameter of planned after repeating five iterations of D-K procedure are demonstrated in Table 2.

Iteration numberγ value achievedMaximum μvalueController order
11201.40311.6485
25.353.83711
31.7630.74816
40.7340.73317
50.7320.73017

Table 2.

Obtained results of the robust controller (μ).

It can be seen that the maximum value of μ is 11.648 that is achieved in the first iteration. Similar to this method, next steps are done to finally value of γ is less than 1. The designed μ controller synthesis is of order 17, and it is achieved after five iterations. In the final iteration, the value of γ is reached to 0.732 and, μ reaches to 0.730 that is less than 1. In other word, the closed-loop scheme has robust performance due to the structured singular parameter is less than 1 in any frequency. Also, the gamma parameter denotes the value that the function FlPKinfinity norm is fewer than that value.

5. Simulation results of the designed robust controller

The considered wind model in this paper is as follows.

Figure 11 depicts the speed of the wind in the third operating area that is a value between 15 and 25 m/s. In this study, the model of the wind has randomly varies from 15 to 24 m/s In our research, the speed of the wind is considered be highly changeable during the times [8, 9]. As in early times the speed of wind turbine is growing, then its parameter is constant in time and its value fell in the end of the period of the time. The reason for selecting this kind of plant of the random wind is, presentation the robust controller has a good performance in various speeds. This express that a sudden variation in the wind speed, the robust controller attempts to control the pitch angle for setting the electrical power at it’s at most efficiency in the third performance area.

Figure 11.

Wind model in different speeds.

Due to the linear dependency among the speed of the rotor and the speed of generator over the gear ratio, the results of applying the control scheme which seeks to maintain the electrical power fixed and regulate speed of the generator- can be examined after modeling of speed of the wind turbine and modeling the robust μ controller.

Considering Figure 12(a), the deflection of output about the nominal value which is changed by noise and disturbance in wind is not high that is acceptable. Maximum variation around the reference value is 2.1% and it is equal to 0.027 kW that expresses nice disturbance cancelation in areas in the matching scheme. In Figure 12(b), at most variations around nominal values are equivalent to 3.068 rad/s and it is equal to 9.2%. It is sates a disturbance cancelation and nice following of generator reference parameter in the existence of disturbance with wind speed in the good way and with the least variations.

Figure 12.

(a) The response of electric power produced by wind turbines to track input reference electrical power using μ controller. (b) The response of speed of the wind turbine generator to track input reference speed of generator using μ controller.

In Figure 13, it can be found that the control effort or the adjustment of pitch angle is between of 14 and 22.5°. Due to changed areas of wind turbine efficiency, aim in this work is, regulating the power and speed of the generator at the nominal parameter in the third performance part. So, in this scenario by growing the value of wind to cut-out, the pitch angle has been improved, and this scenario leads to the decrease of power coefficient and the power is at its nominal value [4]. Furthermore, by reducing the speed of the wind, the blade pitch angle is decreased and at this step, to regulate the power and generator speed at nominal value, robust controller is designed. This controller attempts to control the pitch angle for accessing the high electric power and adjust the speed of generator about its nominal value. So, by applying this kind of controller, the control effort remains fewer than 25°, that is the utmost pitch angle of the wind turbine, stays bounded in the third area.

Figure 13.

Pitch angle adjustment using the μ controller.

5. Conclusion

This work wanted to control the pitch angle of the wind turbine to regulate the speed of the wind turbine generator in the third area of procedure. In the third area, the generator speed, and electrical power are fix in their nominal parameter, and do constant. The existence of noise and disturbance in the model of wind are the main reason for high error rate in the generation of electrical power and generators’ speed. Therefore, this paper suggested that the generation of electrical power and adjustment of the generator speed are practical if robust controller is employed and current errors (uncertainties) in the wind turbine system are taken into account. In most previous works have done, spring constant, damping coefficient and insignificant deviations of the linearization process are listed as uncertainties. These uncertainties are deemed to be true in appropriate weather conditions. Although, cold climate leads the turbine blades to freeze that is followed by mass growth. This mass development results in the decrease of electrical power generation, incorrect model operation and wrong data sending. So, new uncertainties were employed to the system in order to work out the mentioned challenges. After employing these uncertainties, μ controller is presented. Minimal variation about the reference parameter and fewer control action for variation the blades’ angle (to obtain optimal power and adjust the speed of the generator) needs the use of μ controller. To be more accurate, the μ controller is more justifiable system in terms of disturbance cancelation.

The wind turbine that is considered in this paper has the following specifications:

Characteristic of Vestas V29 Wind turbine

Vestas V29 wind turbine characteristicValue
Electrical power of generator225 KW
Rotor diameter29 m
Rotor RPM41/30.8 RPM
Angular speed of rotor4.29 rad/s
Angular speed of generator105.6 rad/s
Frequency50–60 Hz

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Tahere Pourseif, Majid Taheri Andani, Hamed Pourgharibshahi, Hassan Zeynali and Arash Shams (April 2nd 2018). Designing Mu Robust Controller in Wind Turbine in Cold Weather Conditions, Stability Control and Reliable Performance of Wind Turbines, Kenneth Eloghene Okedu, IntechOpen, DOI: 10.5772/intechopen.74626. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="http://www.intechopen.com/embed/stability-control-and-reliable-performance-of-wind-turbines/designing-mu-robust-controller-in-wind-turbine-in-cold-weather-conditions" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

107total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Dual Robust Control of Grid-Connected DFIGs-Based Wind- Turbine-Systems under Unbalanced Grid Voltage Conditions

By Kamel Djamel Eddine Kerrouche, Lina Wang, Alex Van Den Bossche, Azzedine Draou, Abdelkader Mezouar and Larbi Boumediene

Related Book

First chapter

Potential of the Solar Energy on Mars

By Dragos Ronald Rugescu and Radu Dan Rugescu

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us