PV Water Pumping With A Peak-Power Tracker Using A Simple Six-Step Square-Wave Inverter
PV Water Pumping With A Peak-Power Tracker Using A Simple Six-Step Square-Wave Inverter
3, MAY/JUNE 1997
Fig. 2. Typical VI and power characteristics of the PV array for three III. WATER PUMP, ELECTRIC MOTOR, AND INVERTER
different levels of insolation.
The water pump chosen is a centrifugal type, where the
output power is proportional to the cube of its rotor speed.
The output of the inverter is variable-frequency ac. By The output power of the water pump can be written as
operating the six-step square wave inverter to change the
frequency output, and using the inverter as the peak-power (1)
tracker, the number of switches and the switching losses can be where is the constant of the water pump in watts/(mech
minimized. The output power of the PV array, the dc voltage, rad/s) and is the rotor speed in mech rad/s.
and the dc current are fed back to the controller, to adjust the The output torque can be written as
frequency output of the inverter and keep the system operating
at a maximum power capture. (2)
Thus, to vary the output of the water pump, the frequency
supply of the electric motor driving the pump must be varied.
II. PV ARRAY AND OUTPUT POWER CHARACTERISTICS The inverter and the electric motor are represented in Fig. 3.
Fig. 2 shows a typical VI curve for a PV array as the The method of switching chosen is a six-step quasi-square-
conditions change. In the same figure, the output power wave output voltage. An example of per-phase voltage and
versus terminal voltage of the PV array is also drawn. The line-current output of the induction motor are shown in Fig. 4.
optimum operating voltage is the voltage where the PV array The operation of the inverter and electric motor is related
generates maximum power. Thus, by adjusting the voltage at to the characteristics of the PV array. The PV array, as can be
the terminals of the PV array, the output power will vary seen from Fig. 2, behaves approximately as a current source
accordingly. It is shown that the insolation level and the when the voltage output of the PV array is lower than the
temperature of the PV array affect its VI characteristics. optimum voltage, and it behaves as a voltage source when the
Therefore, the optimum output voltage of the PV array changes output voltage of the PV array is higher than the optimum
at different conditions. In Fig. 2, three different levels of voltage.
insolations are illustrated. Curves I1 and P1 correspond to the In the dc representation shown in Fig. 5, the fundamen-
highest insolation level, while curves I3 and P3 correspond to tal component of voltage and current of the motorinverter
the lowest insolation level. The dashed line shown in Fig. 2 relationship is reflected to the dc side as a variable resistor.
716 IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, VOL. 33, NO. 3, MAY/JUNE 1997
(a)
- (4)
The mechanical output power of the water pump and the output
power of the PV array are equal if the losses in the motor and
the power converter are neglected. They can be written as
(5)
- - (6)
-
(7)
-
where
- average optimum voltage of the PV array Fig. 8. Torque-speed characteristics for starting and at optimum condition.
(known);
- current of the PV array at low voltage or short
output taken consistently at the same time step, the maximum
circuit for rated conditions (known);
power can be detected.
- current of the PV array at low voltage or short
circuit for normal conditions (measured);
B. Starting Procedure
power coefficient of the water pump.
The PV array current can be measured, and - can be The starting procedure for the system begins by computing
found from rated conditions. Knowing these values, the ap- the power capability of the PV array at the present condition.
proximate frequency at the optimum operating conditions can The predicted operating frequency is computed, and the fre-
be predicted at any insolation level. Knowing the approximate quency of the inverter is increased from zero to the predicted
optimum frequency, the search can be targeted within a closer operating frequency. As the slip decreases from unity to a
range to the optimum frequency. Hence, the only variable small value, the impedance of the induction motor presents a
is the short-circuit current measured, which depends on the higher value as the slip becomes smaller and corresponds to
insolation and the temperature of the PV panel. The short- a more inductive load. In a dc representation, this condition
circuit current is the indirect sensor of the power capability of is represented as a high resistance. The voltage across the
the PV panel under any insolation level. output of the PV array increases above the optimum operating
voltage.
The output power and slip frequency increase, as the fre-
IV. COMPUTER SIMULATION
quency is increased and the induction motor presents a higher
load to the PV array. Thus, in a dc representation, the fictitious
A. Power Output resistive load appears across the PV array, as shown in Fig.
The output power of the PV array can be tracked by sensing 5, and corresponds to a lower resistance. The voltage across
the power output of the PV array. To vary the output power the PV array is lowered, approaching the optimum voltage.
of the PV array, the frequency of the inverter is varied in As the frequency keeps increasing, the slip also increases,
steps. As the frequency is increased, the speed of the water until at one point the slip reaches that for peak torque.
pump increases. Thus, the power entering the motor increases Any small perturbation around this point (X1 in Fig. 8)
approximately proportionally to the cube of the frequency will bring the system to another operating point much lower
output of the inverter. The current (dc) generated by the PV than the optimum operating point (X3). This condition can
array increases from zero to maximum as the voltage across be illustrated with the curve f1, shown in Fig. 8. To avoid
the PV array decreases from the maximum to an optimum this condition, the frequency must be quickly reduced to one
(smaller) value. step lower (refer to curve f2 in Fig. 8). The change must
Thus, referring the output power versus the PV array be done before the rotor speed decreases to a much lower
voltage, the variation of the output power moves the operating speed. The operating point will eventually settle at point X2
point from lower power (high voltage, lower current) to the on curve f2. Otherwise, if the slip is allowed to increase to
higher power (low voltage, higher current). The power is a large value, the voltage across the PV array will collapse.
varied by varying the inverter frequency. The frequency should Therefore, it is important to detect the optimum frequency
be kept constant at the optimum voltage to get peak power where the maximum power can be tapped from the PV array.
transfer. The process is repeated periodically to account for The detection of maximum power can be done by sensing
any change in insolation. the power output of the PV array at every step. If, in the
The frequency is incremented stepwise. The transition time next step, the power output becomes smaller, it means the
between the old frequency and the new operating frequency maximum power has been passed. The controller should reset
must include the mechanical transient of the system to reach the frequency to one step lower than the current frequency
its new steady state, which is determined mostly by the (from f1 to f2).
mechanical time constant of the system. The data are sampled In the starting procedure, when the system is first started,
just before the new frequency is set. With the filtered PV-array the sensor may react as if a peak value has been detected,
718 IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, VOL. 33, NO. 3, MAY/JUNE 1997
although the frequency is much smaller than the predicted Fig. 10. DC voltage and PV current as the frequency is increased.
optimum frequency. This condition occurs as the operating
point moves from point 0 to point A in the curve shown in It is important to note that moving from point 0 to point
Fig. 8. The phenomenon corresponds to a condition where the A, the power increases and then decreases, while the dc
operating point moves from the left of the peak torque (point voltage keeps increasing. The sensor must be capable of
0), passes the peak torque, and then moves to the right of the recognizing the difference between the peak power detected
peak (point A). When the steady-state condition is achieved, during starting procedure and the normal operating condition.
the operating PV voltage is stable. The movement of operating A simple algorithm in the controller can avoid a wrong
point from point 0 to point A in Fig. 8 corresponds to the detection, especially when the starting procedure is in progress,
power characteristics of the PV array. The power change is or when the optimum power has been achieved. In the control
accompanied by a change in the dc voltage from a lower value strategy implemented, the subsection added to the control
(lower than the optimum voltage) to a higher value (higher algorithm can be explained as follows.
than the optimum voltage).
1) During starting, the dc voltage is very low, the short-
The phenomenon depicted above is illustrated in Fig. 9,
circuit current can be sensed, and the predicted optimum
where the power (PFILT is the filtered output power of the PV
frequency is computed from the equation.
array) increases as the frequency is increased. It can be seen
2) The frequency is then set to start the induction motor by
from Fig. 9 that the power reaches its peak at a frequency
adding a variable
of about 90 rad/s. This peak value is not the optimum power,
but is the first peak power when the speed begins to increase
and the operating point moves from point 0 to point A. The
power decreases from the first peak, and then the operating
point reaches point A. As the frequency is increased, the power 3) The value negative of variable is increased every
begins to increase again until it reaches the second peak. The sample time by a certain frequency increment
second peak is related to point X1 in Fig. 8. Since, at this point,
the condition is unstable, the frequency must be stepped back
one step to operate the system at point X2.
The corresponding dc voltage and PV current are presented 4) Thus, varies from 0 to approximately
in Fig. 10. As the system is started from standstill, the slip is The system controller should monitor the value of
high and the operating voltage of the PV array is very low. As If output power begins to decline while the value of
the slip is reduced and the frequency is gradually increased, is negative, the system detects a false optimum. This false
the PV voltage increases. The knee point of the VI curve, optimum occurs in the process of starting (moving from points
where the optimum power is located, is passed for the first 0 to point A in Fig. 8), where the PV power increases and
time, when the operating point moves from point 0 to point then decreases after reaching the peak between points 0 and
A. The voltage keeps increasing, and the current begins to A. In this case, the controller should continue to increase
decrease. The power decreases as the operation from point 0 the frequency until the actual optimum frequency is reached.
to point A involves reducing the slip of the induction motor. Thus, the conditions for the wrong detection during the starting
As the slip begins to increase, the PV voltage decreases, and procedure can be avoided. The speed (r/min) and the load
the PV power begins to increase. After reaching the optimum (water pump) torque of the system are given in Fig. 11. It can
point, the power begins to decrease. This power reduction is be seen that some time is required for the speed to settle down
the signal used to detect the maximum power of the PV array. as the frequency is increased step by step.
MULJADI: PV WATER TRACKING WITH A PEAK-POWER TRACKER 719
(a)
Fig. 11. Motor speed and water pump load torque as the frequency is
increased step by step.
APPENDIX
INDUCTION MACHINE AND WATER PUMP DATA
(a) ACKNOWLEDGMENT
The authors wish to thank L. Flowers from the Wind
Technology Division, National Renewable Energy Laboratory,
N. Rau, Y. H. Wan, R. Taylor, and L. Pohle for valuable
discussions and for reviewing this paper.
REFERENCES
[1] S. J. Butler, D. M. Sable, F. C.Lee, and B. H. Choo, Design of a solar
array simulator for the NASA EOS testbed, in Proc. 27th Intersociety
Energy Conversion Engineering Conf. (IECEC), San Diego, CA, 1992,
pp. 1.631.66.
[2] S. Cuk and R. D. Middlebrook, Advances in switched-mode power
(b) conversion, in Optimum Power Conversion, vol. III. Pasadena, CA:
Teslaco, 1993, pp. 325345.
Fig. 15. PV array characteristics. (a) Basic characteristics. (b) Voltage and [3] P. Huynh and B. H. Cho, Design and analysis of microprocessor-
current scales to simulate the insolation level. controlled peak-power-tracking system, in Proc. 27th Intersociety En-
ergy Conversion Engineering Conf. (IECEC) San Diego, CA, 1992, pp.
1.671.72.
power condition. The data for the power, voltage, and current [4] R. M. Hilloowala and A. M. Sharaf, A rule-based fuzzy logic controller
are sampled at the end of the to assure that the steady- for a PWM inverter in photovoltaic energy conversion scheme, in Conf.
Rec. 1992 IEEE/IAS Annu. Meeting, pp. 762767.
state data are used. The data are renewed at the end of the next [5] B. K. Bose, P. M. Szczesny, and R. L. Steigerwald, Microcomputer
The voltage and current feedback can also be filtered control of a residential photovoltaic power conditioning system, IEEE
to reduce the noise before it is fed back to the controller. Trans. Ind. Applicat., vol. IA-21, pp. 852859, Sept./Oct. 1985.
[6] G. Vetter and W. Wirth, Suitability of eccentric helical pumps for turbid
water deep well pumping in photovoltaic systems, Sol. Energy, vol. 51,
V. CONCLUSION no. 3, pp. 205214, Sept. 1993.
[7] E. W. Ter Horst, J. H. Boumans, and K. Blok, The Terschelling
The system described here is a PV system for water pump- PV/Wind system after five years of operation, in Proc. European
ing. By using a six-step square-wave inverter, both as the Community Photovoltaic Solar Energy Conf., Florence, Italy, May 913,
1988, pp. 275278.
variable-frequency source and as the peak-power tracker, the [8] H. Gabler, G. J. Gerdes, and J. Luther, Wind and solar energy supply for
switching losses in the power electronics are minimized, and a sewage plantThe Fehmarn project, in Proc. European Community
the number of switches is reduced. The system is coupled with Wind Energy Conf., Madrid, Spain, Sept. 1014, 1990, pp. 692696.
[9] P. C. Krause, Analysis of Electric Machinery. New York: McGraw-
a centrifugal water pump, and the controller is set to adapt to Hill, 1986.
changing conditions. The system can be implemented by a
simple microcontroller. The only control variables needed are
the power, voltage, and current output of the PV array to be
fed back to the microcontroller. The PV-array power is used
to detect the maximum power, the dc-bus voltage is used to Eduard Muljadi (S83M84SM94) received the
B.S.E.E. degree from Surabaya Institute of Technol-
detect the operating voltage of the PV array, and the current ogy, Surabaya, Indonesia, in 1981 and the M.S. and
output of the PV array is used to compute the approximate Ph. D. degrees from the University of Wisconsin,
optimum frequency. Madison, in 1984 and 1987, respectively.
From 1988 to 1992, he was with the Department
The insolation level is simulated by changing the coef- of Electrical and Computer Engineering, California
ficients of the voltage and current at different times. The State University, Fresno, where he became an Asso-
data for the basic VI characteristics are stored in one file, ciate Professor in 1991. In June 1992, he joined the
National Renewable Energy Laboratory (NREL),
while another file stores the light-intensity characteristics. The Golden, CO. His current research interests are in the
parameters of the induction machine and the water pump fields of electric machines and power electronics, with emphasis on renewable
can be changed to simulate different systems. The controller energy applications.
Dr. Muljadi is a member of the IEEE Industry Application Society Industrial
parameters, the size of the steps, can be changed to observe Drives, Electric Machines, and Industrial Power Converter Committees. He is
different responses. also a member of Eta Kappa Nu and Sigma Xi.