PV Integration in Low-Voltage Feeders With Demand Response: Xiangkun Li, Theodor Borsche and Göran Andersson
PV Integration in Low-Voltage Feeders With Demand Response: Xiangkun Li, Theodor Borsche and Göran Andersson
PV Integration in Low-Voltage Feeders With Demand Response: Xiangkun Li, Theodor Borsche and Göran Andersson
With over 130 GW installed worldwide [1], Photo-Voltaics (PVs) are currently the third most dominant form of renewable
energy, behind hydropower and wind, in terms of capacity. The average annual growth rate of PV has been 60 % over the
last five years [2]. Due to strong policy support, PV was projected to surpass wind as the fastest growing source of electricity
generation for the first time in 2013 [3].
However, grid integration of massive amounts of distributed, stochastic generation has introduced many technical challenges
for power system operators. In the case of PV, one key issue that arose in Europe is voltage constraint violations in distribution
networks. Approximately 80 % of PV power is fed into the 400 V Low-Voltage (LV) grid [4]. When generation in-feed exceeds
local demand, voltage levels rise, causing damage to consumer and grid equipment if no action is taken. Safety concerns are
also of importance as protection equipment now need to take into account bidirectional power flows.
Staying within the allowed voltage band is a greater challenge in LV systems due to a comparatively smaller short circuit
impedance and a larger R⁄X ratio. This often causes voltage violations before line limits are reached [5]. With ambitious renewable
energy targets to meet in the next few decades, many solutions have been proposed to increase the PV hosting capacity of
LV grids. Seven key ways are outlined below. The first two methods aim to increase the system’s ability to transfer power to
higher voltage levels, while the remaining methods attempt to handle the surplus power locally [6]:
Reactive power control Consuming reactive power while injecting active power limits voltage rise, but increased reactive power
flows may create voltage issues in the medium voltage network and lead to higher losses.
On-Load Tap-Changing Transformer (OLTC) Lowering the secondary transformer voltage may not be possible if many lines
are connected to a single bus and require a constant voltage.
Correlation with load profiles While assuming zero correlation is too conservative, central Europe exhibits poor load correla-
tions due to low levels of air conditioning.
Active power curtailment Key method implemented today due to its simplicity and allowance for emergency control.
Orientation change of PV installations A smoother generation curve is obtained when peak production periods are shifted, but
this method currently lacks economic incentives as producers are only paid for energy generated.
Decentralized storage Incentives for self-supply and private storage exist, but implementing peak production storage depends
on future regulatory and costs structures.
Demand Response (DR) Power and energy curtailment are minimized and local energy use is prioritized, potentially reducing
losses and system congestion.
Grid reinforcement and expansion is also always an option, but should be avoided if possible due to high costs and a lack of
exploitation of existing resources.
This paper investigates the potential of DR to increase the PV hosting capacity of a LV network. The charging schedule
of Electric Water Heaters (EWHs) is optimized to manage PV generation peaks in order to limit the maximum voltage
rise of the line to 3 % as per DACHCZ code specifications [7]. The control strategy relies only on power measurements at
the transformer. Uncertainties associated with weather forecasts are taken into account using model predictive control. Cost
structures are investigated, as well as limited controllability due to the grouping of assets and the trade-off between switching
frequency and PV in-feed shedding. Simulations are run on real data from Swiss households and PV installations.
The paper is organized as follows: Section II gives an overview of the test system and the data used. Section III describes
the control scheme for improving the PV hosting capacity with DR. Results are presented in Section IV and discussed in
Section V, before Section VI concludes.
II. S YSTEM D ESCRIPTION
A. Grid Topology
The grid topology used is based on a representative Swiss suburban 400 V distribution network (Figure 1). Two 630 kV A
transformers connect the 400 V LV grid to the 16 kV Medium-Voltage (MV) grid. Two of the eleven feeders connected are
considered here. This double-feeder supplies a residential area with twenty nodes. PV installations are assumed to be evenly
distributed at the load buses. Voltage at the main busbar is fixed to 1.01 p.u. and cannot be changed as it might lead to voltage
limit violations in the remaining nine feeders.
B. Measured Data
PV production profiles for active and reactive power are obtained from two installations rated at 61 kWp and 30 kWp , with
an average capacity factor of 9.83 %. Load information is extrapolated from measurements taken from 53 Swiss households.
Connections are determined based on average transformer loading considerations and the maximum allowable voltage drop
during peak periods. The simulation incorporates 200 households, with 10 households per bus, where minimum voltage levels
never drop below 0.96 p.u. Measured data for load and PV production is available from January 1st, 2013 to December 31st,
2013 with a temporal resolution in minutes. Missing data is linearly interpolated when missing values span less than 15 minutes.
Days with more than 15 consecutive minutes of missing data are eliminated from the simulation.
C. Hot Water Consumption Profile
The state of charge of the EWHs are assumed to be known at all times. An initial state is chosen randomly, and later states
are computed using a simple hot water demand model that maps Swiss consumption rates onto the Becker daily water draws
profile (Figure 2) [8]. Water heating energy demands account for approximately 14 % of total household electricity use in
Switzerland [9]. This corresponds to an average hot water demand of 1.03 kW h per household per day. The EWHs considered
are rated at 4.5 kW. They are charged two to four hours each day and have a storage capacity of 24 hours. One EWH is needed
at each node to meet the hot water demands of the households1 . For simplification, it is assumed that all households follow
the same hot water usage profile and use identical amounts of hot water each day.
D. Cost Structure
SwissIX spot market prices are given for every hour of the measurement period. For simplicity, they are assumed to be
known to the controller when dispatching the flexible loads. Taking bidding strategies into account is beyond the scope of
this paper, but would most likely not significantly affect the dispatch. All PV plant owners are assumed to be paid a fixed
feed-in tariff of 0.36 e kW−1 h−1 . Because the installed PV capacity in Switzerland is too low to cause any PV curtailment,
the rules for PV shedding during peak hours are based on the German Renewable Energy Resources Act [10]: the DSO is to
compensate plant owners for 95 % of their lost income plus any additional expenses incurred, minus any expenses saved, due
to switching actions in the event of voltage instability. This penalty is necessary to provide an economic incentive for utilities
to charge EWHs during peak price periods.
1 Multiple households may be in one building, and not all water heating is done using electricity.
Waterdraw profile
1.5
0.5
0
00:00 04:00 08:00 12:00 16:00 20:00 24:00
time of day
PV forecast error
time of forecast
300 Forecast
Measured PV
200
[kW]
100
0
00:00 04:00 08:00 12:00 16:00 20:00 24:00
time of day
E. PV Forecast
PV prediction errors are introduced to study the effects of imperfect information. A random, normally distributed error factor
is generated at every time step and multiplied with the smoothed measured data curve. The mean of the error factor increases
linearly up to a maximum intraday mean of 20 % [11]. The standard deviation is set to allow for a maximum overestimation
at any given time of 40 % and a maximum underestimation of 30 % [12]. Errors for all past time steps are set to zero, and a
noise factor is added [13]. Predictions can be updated hourly with increasing accuracy.
III. I NCREASING H OSTING C APACITY WITH DR
We propose to only measure power backflows at the transformer. Limiting the backflows to a predefined level by directly
switching flexible loads, and in extreme situations also some PV generation, should be sufficient in preventing over-voltages
in the feeder as long as the distribution of loads and PV installations is known, and PV installations are not concentrated at
PV/load prediction
EWH model spot prices
MPC
16 kV:400 V
1 p.u.:1.01 p.u.
1 2 3
S
... ...
P load + B1 P EWH P1PV P2PV P3PV
Figure 4. Proposed control scheme: Power is measured at the LV transformer (blue); the controller then solves an optimization problem based on forecasts
(green) and sends control signals (purple) to flexible loads and PV installations.
a few nodes. This approach, schematically presented in Figure 4, minimizes the number of real-time measurements needed
while providing sufficient dependability.
Flexible loads and distributed generation must allow for Distribution System Operator (DSO) control, but this controllability
may be granular in power, i.e. only groups of assets can be switched. Delays in the control signal could be included in the
framework, but is subject to future research.
Flexible loads are optimally dispatched, taking into account electricity costs, asset constraints, and an estimate of the daily
energy consumption of the EWHs. Not considering exact power flows and voltage levels in the controller simplifies calculations
by eliminating non-linear power flow computations from the optimization. The linear program describing the load dispatch
is explained in III-A. Results of the dispatch are verified using historic load and PV data on the known grid topology (see
Section II), and possible voltage violations are identified. Finally, total costs are computed using spot market prices, actual
load dispatch (including passive loads), and PV shedding penalties.
A. Deterministic Optimization with Perfect Forecasts
The goal of the optimization is to minimize system costs within voltage and EWH state of charge constraints. While it is
cheaper to heat water during night hours, if excess PV production is expected, some or all of the water heating is shifted to
daytime hours to avoid costly PV shedding penalties. The optimization is done with a temporal resolution, ∆t, of 10 minutes
as that is the relevant time scale for grid operators to take control action against voltage disturbances. The deterministic
optimization problem formulation is as follows:
N
PiEWH Bi,k ∆t + cshed
X X X
min cspot
k
shed
Ps,k ∆t
B
k=1 i∈I s∈S
!
X
+csw on
(Si,k off
+ Si,k ) , (1)
i∈I
s.t. ∀k ∈ [1, N ], ∀i ∈ I, ∀s ∈ S
X
P min ≤ Pkload − PkPV + shed
Ps,k
s∈S
(2a)
PiEWH Bi,k ≤ P max
X
,
i∈I
xSoC SoC EWH
i,k+1 = xi,k + Pi Bi,k ∆t − ξk , (2b)
xSoC
i,1 = xSoC
0 , (2c)
0≤ xSoC
i,k ≤ xSoC
max , (2d)
xSoC
i,N ≥ ξ1 , (2e)
0≤ Pkshed PV
≤ Ps,k , (2f)
Bi,k = on
Bi,k−1 + Si,k − off
Si,k . (2g)
The first term in (1) minimizes the cost of water heating: cspot corresponds to spot market prices, I is the set of all EWHs
or EWH groups, PiEWH is the power rating of EWH i, and Bi,k is the binary variable that specifies whether EWH i is on or
off during time step k. The second term minimizes PV power curtailment: cshed represents the PV shedding costs, S is the set
of all PV installations or groups, and Psshed is the power shed by PV unit s during ∆t. Finally, the last term penalizes each
switching action S on or S off with a cost of csw . Si,k
on
is one if EWH i is switched from off to on during time step k, and zero
otherwise. Si,k is defined analogously.
off
Constraint (2a) sets power flow limits at the transformer: the power used by passive loads plus any electric water heating
minus the net PV power the system is able to accept must be within system limits. The positive limit is determined by the
yearly peak-load in the double feeder (i.e. 220 kV A), while the negative limit is set at a level that prevents over-voltages in
ordinary simulation scenarios.
Constraints (2b) - (2e) control the energy behaviour of the EWHs in the system: The state of charge of each EWH at time
step k+1 is equal to the state at time step k plus any water heating done minus any water drawn (2b). The initial state of
charge must be set (2c). The maximum capacity of each EWH is limited (2d). The final state of charge at the end of the day
must be large enough to supply the water needs for the first time step of the next day (2e).
Constraint (2f) limits the amount of PV curtailed by the actual level of PV generation during the current time step, and the
final constraint (2g) affects switching behaviour: whether an EWH is on or off at time step k depends on its past state and
whether or not there was any switching action.
The EWH dispatch Bi,k is then applied to the system. Voltage evolution and losses in the feeder are simulated in minute
resolution using the data described in Section II. The deterministic optimization with perfect forecasts sets the benchmark for
the lowest achievable costs and the minimum power curtailment levels in the system. Results form a basis for comparison with
other scenarios.
B. Grouping of Control Signals
Optimal system behaviour results when EWHs and PV installations are controlled individually. However, this level of control
is currently not implemented in Switzerland, where all EWH loads of a distribution company may be controlled in only two
or three groups. To investigate the effects of limited control on power curtailment and costs, I and S are varied while the
xSoC
k , w, and ξ variables are adjusted accordingly. When PV installations are individually controlled, power outputs of each
installation is ramped down to the maximum allowed limit. When control is grouped, individual installations are turned off
completely when over-voltages occur. Loads and PV installations in one group are spread out geographically so no section of
the line is heavily affected when a group is switched in or out.
C. Switching Penalty
Increased switching of flexible loads in order to accommodate additional PV capacity can have a negative impact on the
lifetime of certain devices (e.g. heat pumps and air conditioners). The virtual cost csw is varied to find the trade-off between
energy loss and switching behaviour.
D. Model Predictive Control Formulation
A model predictive control algorithm is implemented to account for limited PV forecast accuracy. The forecast uncertainty
is created artificially as described in II-E. Optimization is performed at every full hour, with a time horizon running until
midnight, thus decreasing by one hour each time the optimization is run. The horizon resets at the beginning of the next day.
The time horizon is not kept constant due to the generation patterns of PV: there is no production at night and spot market
prices are not know a day in advance.
The optimization problem is adjusted accordingly, and Equation (1) becomes:
N −1
PiEWH Bi,t+k|t ∆t + cshed ·
X X
min cspot
t+k|t
B
k=0 i∈I
!
X X
· shed
Ps,t+k|t ∆t +c sw on
Si,t+k|t + off
Si,t+k|t , (3)
s∈S i∈I
with N − 1 as the number of time steps until the end of the day, t as the time of execution and with xt+k|t denoting the
prediction of x at time t, k steps into the future. The first Bi,t+k|t is applied; then new PV prediction is available and the
optimization is solved once again. Past time steps no longer need to be considered as long as the initial state of charge of the
EWHs at time step k is set to the final state of charge at the end of the previously implemented time step,
xSoC SoC
i,t+0|t = x̂t ∀i ∈ I . (4)
IV. R ESULTS
The following definitions from [6] are used in this paper:
PV Penetration Total yearly solar energy fed into the grid divided by the yearly energy consumption of all consumers connected
to the grid.
Hosting capacity Maximum PV penetration at which the 3 % voltage rise limit is not violated.
A. Deterministic Optimization Results
Table I shows the changes in the hosting capacity of the test grid under different PV integration methods. DACHCZ regulation
assumes zero system load with all PV installations running at maximum power. Load correlations take into account actual PV
generation and consumption levels, while demand response introduces an additional flexible capacity of 360 kW h per day into
the system. Hosting capacity is found assuming a uniform distribution of PV in-feed at every node in the grid. Comparing
DR to the DACHCZ regulation gives an increase of nearly 25 % points, effectively doubling the hosting capacity. However,
it would be more fair to compare the DR approach with load correlation, as both run under similar assumptions. Even here,
over 12.5 % points can be added: a sizable gain.
Figure 5 shows how system losses behave at different PV penetration levels. Losses here refer only to the losses on the LV
feeder. Losses in the transformer and at higher voltage levels are not regarded. At a PV penetration level of 0 %, grid losses
account for 1.81 % of the total power delivered during the year. Losses are minimized at a penetration level of 24.1 %. The
Table I
PV HOSTING CAPACITIES UNDER DIFFERENT CONDITIONS
vertical lines give the hosting capacity with the three approaches. DR is only used when penetration surpasses 40 %. The purple
line is a fit of all points below this level. While one would expect increasing losses for higher penetration levels due to the
shifting of demand to peak production hours, losses can be kept at levels similar to the no-PV case. Furthermore, simulations
not shown in the graph suggest that maximum losses will never significantly exceed initial system loss levels, as backflows
are limited and PV shedding will occur before additional losses are generated in the LV feeder.
Grid Losses
6
DACHCZ regulation
3
load correlation
Additional grid losses [%]
0
−3
−6
demand response
−9
−12
−15
−18
0 10 20 30 40 50
PV penetration level [%]
0
0 4 8 12 16 20
Number of groups
Figure 6. Increases in energy shed and total costs due to load grouping
When EWH control is grouped, energy shed and total system costs increase at almost identical rates. EWH charging costs
increase as more EWHs than necessary must be turned on during peak periods to accommodate the PV in-feed. This in turn
causes state of charge levels to rise, resulting in a lower capacity to accommodate PV at other times of the day.
When the control of PV installations is grouped and units can no longer be ramped, energy shed increases much more
dramatically as whole installations must now be turned off even if current power levels are only marginally above system
Grouped PV Control
75
Energy shed Total costs
60
Percent Increase
45
30
15
0
4 8 12 16 20
Number of groups
limits. Costs increase at a lower rate because as more PV is turned off, more EWH charging can be shifted to lower cost
periods. The decrease in charging costs helps to offset the increase in energy shedding costs.
C. Effects of Switching Penalties
The switching penalty is an arbitrary value set in relation to the PV shedding costs, assumed here to be 34.2 ct/kWh.
Intuitively, as the switching penalty increases, the amount of switching decreases (Figure 8) while energy shedding and total
system costs increase (Figure 9).
For comparison purposes, each EWH is switched an average of 4.47 times a day when the system is optimized to heat water
at the lowest spot prices, with no consideration for PV shedding. The nominal switching value per EWH per day is two: once
to turn on and once to turn off. A small switching penalty has a large initial effect, as it removes any unnecessary switching
the EWHs do between periods of similar spot prices. Higher penalties must consider the trade-off between energy shed / cost
increases and the wear and tear on the equipment due to switching.
2
0 10 20 30 40 50
Switching penalty [EUR/switch]
0
2 3 4 5 6
Number of switches per EWH per day
Figure 9. Increases in energy shed and total costs due to switching penalties
D. Effects of Forecast Inaccuracies
The effects of PV forecast inaccuracies are shown in Table II. Using day-ahead forecasts, the EWH charging schedule for the
upcoming 24-hour period is set each day at 00:00. Then hourly forecasts are assumed to be available, and the model predictive
control algorithm updates the charging schedule of the EWHs at the start of each hour based on more accurate predictions.
The performance of these two methods are compared against the deterministic optimization reference case. When faced with
uncertainties, the amount of energy that is shed is highly increased. As expected, an hourly update can considerably reduce
these losses. While it is conceivable to send a daily switching plan to each EWH, it seems much more promising to have
real-time switching capabilities.
Table II
E FFECTS OF PV F ORECAST I NACCURACIES
V. D ISCUSSION
This paper proposed and investigated a simple control architecture that increases the PV hosting capacity of a LV-feeder with
DR. The controller only relies on real-time power flow measurements at the LV transformer, thus minimizing investment costs
in measuring infrastructure. Results are promising, as the hosting capacity of the test system can be substantially increased
from 28.57 % with current regulation to over 50 % when flexible loads are used to absorb power locally. As opposed to
methods relying on reactive power management for voltage compliance, no unnecessary currents are injected into the feeder.
Furthermore, losses can be slightly reduced due to better correlation between local production and consumption.
Grouping of both flexible loads and PV installations was investigated. While control performance is reduced as the number of
groups is reduced, and thus increasing the granularity of the control action, it was shown that full or continuous controllability
of assets is unnecessary. This is a relevant result, as currently installed communication channels for load switching, such as
ripple-control in Switzerland, only allow for the switching of load groups, and many existing PV converters can only be turned
on or off even though they should technically be able to continuously set power output.
Minimal shedding of PV power is achieved when loads can be switched at any time, but this may be undesirable from a
life-time perspective and may alienate customers from participating in such a DR scheme. However, most of the switching
adds little benefit to overall system performance. By adding a virtual penalty, switching can be reduced to an acceptable level
of less than four times a day, i.e., activating a load less than twice a day.
Despite these encouraging findings, there are challenges in the implementation of the approach. Firstly, backflow limits
need to be identified. This can be done heuristically by running simulations of typical years. For this, PV in-feed and demand
profiles of local customers, as well as the locations of the PV installations must be known or estimated. One may also use
analytic approaches, considering installed peak power and minimum expected load. Note, we assumed evenly distributed loads
and PV installations. If the PV production is concentrated at one node, there may be current violations within the feeder even
if power consumption and production is balanced. This cannot be handled by the presented approach.
To optimally use the available flexibility, the daily energy demand of the flexible loads must also be known. For this, dedicated
measurements of EWHs at customer premises must be available. These measurements do not need to be communicated in
real time and only need to be given with a time resolution in days, but even this is currently unavailable. It would be straight
forward to gather this data with an Advanced Metering Infrastructure (AMI), but the Smart Meters (SMs) currently rolled out
often do not have dedicated measuring capabilities for switchable loads, requiring a second meter to be installed.
Finally, PV forecasts were subject to severe uncertainty. Model Predictive Control (MPC) approaches can handle a good
part of the uncertainty, but still leave room for improvement. Techniques from robust control and stochastic control should
be investigated, and it is the authors’ belief that improvements are still achievable. Robust control would shift most flexible
loads to the day, unnecessarily increasing costs but guaranteeing minimal shedding. Stochastic control, either using probability
constraints or a scenario-based approach to minimize expected costs, might be the favorable solution: the problem as such
does not require the strict guarantees offered by robust control – PV shedding is an undesirable but always feasible strategy.
VI. C ONCLUSION AND O UTLOOK
This paper explores the technical potential of DR for integration of PV generation in LV grids and shows its capabilities
using a simple control method: only power flow at the transformer is measured in real-time; all other relevant data, such as
the daily energy consumption of EWHs, can be gathered using an AMI. Granular control of assets was shown to be sufficient,
meaning a legacy control infrastructure, such as ripple-control can be used. The hosting capacity of the feeder was increased
to over 50 % while loads were activated less than twice a day. Further studies should look at stochastic control to consider
the inherent uncertainty in the optimization inputs. A more accurate model for hot water consumption, ideally based on real
measurements, should be implemented to account for seasonal demand variations. In this context, it would be interesting to see
comparable analyses for regions with higher PV capacity factors and that use air-conditioning as flexible loads (e.g. California,
southern Europe or South-East Asia). Furthermore, to gauge whether DR can make a feasible business case, the study must
be explored from a regulatory point of view and cost incentives must be clearly identified.
VII. ACKNOWLEDGMENT
The authors gratefully received grid data and PV production time-series from project VeIN, and household load data gathered
by ewz. Theodor Borsche is financed by the project Distributed Demand Response, which is funded by the CTI (Commission
for Technology and Innovation).
R EFERENCES
[1] G. Masson, M. Latour, M. Rekinger, L.-T. Theologitis, and M. Papoutsi, “Global Market Outlook for Photovoltaics 2013-2017,” European Photovoltaic
Industry Association, Tech. Rep., 2013.
[2] REN 21 Steering Committee, “Renewables 2013 Global Status Report,” Renewable Energy Policy Network for the 21st Century, Paris, Tech. Rep., 2013.
[3] J. Isola and A. Mccrone, “Solar to Add More Megawatts than Wind in 2013 , for First Time,” New York City, Sep. 2013. [Online]. Available:
https://www.bnef.com/PressReleases/text/324
[4] B. Ernst and B. Engel, “Grid Integration of Distributed PV-Generation,” in IEEE PES General Meeting 2012. San Diego: IEEE, 2012.
[5] T. Degner, F. I. Germany, S. M. A. Solar, T. Ag, G. Arnold, M. Breede, T. Reimann, and S. A, “Increasing the Photovoltaic-system Hosting Capacity
of Low Voltage Distribution Networks,” in 21st International Conference on Electricity Distribution, no. 1243, Frankfurt, 2011, pp. 6–9.
[6] C. Bucher, G. Andersson, and L. Küng, “Increasing the PV Hosting Capacity of Distribution Power Grids – A Comparison of Seven Methods,” in 28th
European Photovoltaic Solar Energy Conference and Exhibition (PVSEC), Paris, 2013, pp. 4231 —- 4235.
[7] J. M. Gerhard Bartak, Hansjörg Holenstein, “Technical Rules for the Assessment of Network Disturbances, 2nd ed., Austria, Germany, Czech Republic,
Switzerland,” Tech. Rep., 2007.
[8] P. Fairey and D. Parker, “A Review of Hot Water Draw Profiles Used in Performance Analysis of Residential Domestic Hot Water Systems,” Florida
Solar Energy Center / University of Central Florida, Cocoa, Tech. Rep., 2004.
[9] Environment Switzerland, “State of the Environment - Households and Consumption,” Tech. Rep., 2013.
[10] EEG, “Act on Granting Priority to Renewable Energy Sources,” Tech. Rep., 2008.
[11] V. Lenzi, A. Ulbig, and G. Andersson, “Impacts of Forecast Accuracy on Grid Integration of Renewable Energy Sources,” in IEEE PES PowerTech
2013. Grenoble: IEEE, Jun. 2013.
[12] E. Lorenz, T. Scheidsteger, J. Hurka, D. Heinemann, and C. Kurz, “Regional PV Power Prediction for Improved Grid Integration,” in 25th European
Photovoltaic Solar Energy Conference and Exhibition (PVSEC), vol. 19, no. Appl. 2011, Valencia, 2010, pp. 757–771.
[13] F. Oldewurtel, A. Parisio, C. N. Jones, M. Morari, D. Gyalistras, M. Gwerder, V. Stauch, B. Lehmann, and K. Wirth, “Energy Efficient Building Climate
Control Using Stochastic Model Predictive Control and Weather Predictions,” in American Control Conference, Baltimore, 2010, pp. 5100–5105.