-
Notifications
You must be signed in to change notification settings - Fork 4
/
esm-lecture-16.tex
735 lines (478 loc) · 20.4 KB
/
esm-lecture-16.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
% NB: use pdflatex to compile NOT pdftex. Also make sure youngtab is
% there...
% converting eps graphics to pdf with ps2pdf generates way too much
% whitespace in the resulting pdf, so crop with pdfcrop
% cf. http://www.cora.nwra.com/~stockwel/rgspages/pdftips/pdftips.shtml
\documentclass[10pt,aspectratio=169,dvipsnames]{beamer}
\usetheme[color/block=transparent]{metropolis}
\usepackage[absolute,overlay]{textpos}
\usepackage{booktabs}
\usepackage[utf8]{inputenc}
\usepackage[scale=2]{ccicons}
\usepackage[official]{eurosym}
\usepackage{tikz}
\usetikzlibrary{arrows.meta}
\usepackage{xcolor}
\makeatletter
\def\mathcolor#1#{\@mathcolor{#1}}
\def\@mathcolor#1#2#3{%
\protect\leavevmode
\begingroup
\color#1{#2}#3%
\endgroup
}
\makeatother
%use this to add space between rows
\newcommand{\ra}[1]{\renewcommand{\arraystretch}{#1}}
\setbeamerfont{alerted text}{series=\bfseries}
\setbeamercolor{alerted text}{fg=Mahogany}
\setbeamercolor{background canvas}{bg=white}
\newcommand{\R}{\mathbb{R}}
\def\l{\lambda}
\def\m{\mu}
\def\d{\partial}
\def\cL{\mathcal{L}}
\def\co2{CO${}_2$}
% for sources http://tex.stackexchange.com/questions/48473/best-way-to-give-sources-of-images-used-in-a-beamer-presentation
\setbeamercolor{framesource}{fg=gray}
\setbeamerfont{framesource}{size=\tiny}
\newcommand{\source}[1]{\begin{textblock*}{5cm}(10.5cm,8.35cm)
\begin{beamercolorbox}[ht=0.5cm,right]{framesource}
\usebeamerfont{framesource}\usebeamercolor[fg]{framesource} Source: {#1}
\end{beamercolorbox}
\end{textblock*}}
\usepackage{hyperref}
%\usepackage[pdftex]{graphicx}
\graphicspath{{graphics/}}
\DeclareGraphicsExtensions{.pdf,.jpeg,.png,.jpg}
\def\goat#1{{\scriptsize\color{green}{[#1]}}}
\newcommand{\ubar}[1]{\text{\b{$#1$}}}
\let\olditem\item
\renewcommand{\item}{%
\olditem\vspace{5pt}}
\title{Energy System Modelling\\ Summer Semester 2020, Lecture 16}
%\subtitle{---}
\author{
{\bf Dr. Tom Brown}, \href{mailto:tom.brown@kit.edu}{tom.brown@kit.edu}, \url{https://nworbmot.org/}\\
\emph{Karlsruhe Institute of Technology (KIT), Institute for Automation and Applied Informatics (IAI)}
}
\date{}
\titlegraphic{
\vspace{0cm}
\hspace{10cm}
\includegraphics[trim=0 0cm 0 0cm,height=1.8cm,clip=true]{kit.png}
\vspace{5.1cm}
{\footnotesize
Unless otherwise stated, graphics and text are Copyright \copyright Tom Brown, 2020.
Graphics and text for which no other attribution are given are licensed under a
\href{https://creativecommons.org/licenses/by/4.0/}{Creative Commons
Attribution 4.0 International Licence}. \ccby}
}
\begin{document}
\maketitle
\begin{frame}
\frametitle{Table of Contents}
\setbeamertemplate{section in toc}[sections numbered]
\tableofcontents[hideallsubsections]
\end{frame}
\section{The World is Not a Perfect Optimization Model}
\begin{frame}
\frametitle{We should be skeptical about models and modellers}
\begin{columns}[T]
\begin{column}{7.5cm}
\vspace{.2cm}
\centering
\includegraphics[width=7cm]{2019-01-10-IEEFA-EIA-coal-all-consumption-forecasts-470-x-395-v2-768x646}
\end{column}
\begin{column}{7.5cm}
\vspace{.2cm}
\centering
\includegraphics[width=7cm]{auke.jpg}
\end{column}
\end{columns}
\end{frame}
\begin{frame}
\frametitle{We should be skeptical about models and modellers}
\begin{columns}[T]
\begin{column}{9cm}
\vspace{.2cm}
\centering
\includegraphics[width=9.5cm]{hubbert.png}
\end{column}
\begin{column}{5cm}
\vspace{.2cm}
\begin{itemize}
\item Possible scenario projected from 1956 by US geologist M. King Hubbert
\item Oil production in the US did indeed peak in the 1970s, but returned to peak height in last decade thanks to shale oil extraction with fracking
\item Nuclear expanded but plateaued
\item \alert{What might we be getting wrong in 2020?}
\end{itemize}
\end{column}
\end{columns}
\source{\href{http://www.energycrisis.com/Hubbert/1956/1956.pdf}{Hubbert, 1956}}
\end{frame}
\begin{frame}
\frametitle{We should be skeptical about models and modellers}
Models can:
\vspace{-.3cm}
\begin{itemize}
\item \alert{under- or overestimate rates of change} (e.g. under: PV uptake, over: onshore wind in UK/Germany/Netherlands)
\item \alert{underestimate social factors} (e.g. concern about nuclear / transmission / wind)
\item \alert{extrapolate based on uncertain data} (e.g. oil reserves, learning curves for PV)
\item \alert{focus on easy-to-solve rather than policy-relevant problems} (e.g. most research)
\item \alert{neglect uncertainty} (e.g. in short-term due to weather forecasts, or in long-term due to cost, political uncertainty and technological development)
\item \alert{neglect need for robustness} (e.g. securing energy system against contingencies, attack)
\item \alert{neglect complex interactions of markets and incentive structures} (e.g. abuse of market power, non-linearities not represented in models, lumpiness, etc.)
\item \alert{neglect non-linearities and non-convexities} (e.g. power flow, or also learning curves, behavioural effects, perverse local optima, many, many more)
\end{itemize}
\end{frame}
\section{Robustness to Different Weather Years}
\begin{frame}
\frametitle{Different Weather Years}
Many of the simulations we looked at in this course, and many in the literature, used single weather years to determine optimal investments.
This is problematic since:
\begin{itemize}
\item Weather changes from year to year
\item There are decadal variations of wind
\item Demand changes (particularly space heating demand during cold years)
\end{itemize}
But computing investments against 30 years of data (262,800 hours) is not feasible.
\end{frame}
\begin{frame}
\frametitle{Different Weather Years}
If we use different weather years to optimize sector-coupled European model with net-zero CO$_2$ emissions (including industry) we see broadly stable technology choices but variations in total system costs of up to 20\%. NB: In real world cannot reoptimize investment every year!
\centering
\includegraphics[width=13cm]{lin-total_costs.png}
\source{Lin Yang MA thesis}
\end{frame}
\begin{frame}
\frametitle{Different Weather Years}
Biggest changes are driven by space heating demand. Cold years (like 2010) are more expensive.
\centering
\includegraphics[width=14cm]{lin-heating.png}
\source{Lin Yang MA thesis}
\end{frame}
\begin{frame}
\frametitle{Different Weather Years}
Optimal technology investments do not change dramatically from year to year. Here we show the mean capacities with standard deviation.
\centering
\includegraphics[width=14cm]{lin-generators.png}
\source{Lin Yang MA thesis}
\end{frame}
\begin{frame}
\frametitle{Different Weather Years}
If we fix the optimal technology investments based on the weather of one year ($x$-axis), then run the dispatch over all 30 years (900 simulations in total), we can assess average curtailment and load-shedding. Using coldest year 2010 gives low load-shedding but high curtailment.
\centering
\includegraphics[width=14cm]{lin-run_30.png}
\source{Lin Yang MA thesis}
\end{frame}
\begin{frame}
\frametitle{Using 2010 investments}
Using coldest year 2010 guarantees virtually no load-shedding in entire 30 years, but leads to excess energy in most years.
Better to store excess energy from warmer years (e.g. chemically).
\centering
\includegraphics[width=14cm]{lin-use_2010.png}
\source{Lin Yang MA thesis}
\end{frame}
\section{Effects of Climate Change on Energy System}
\begin{frame}
\frametitle{Climate change}
\begin{itemize}
\item What are the consequences of climate change for highly renewable energy systems?
\item How will generation patterns for wind and solar change?
\item What will be the effects on the dimensioning of wind,
solar, storage, networks and backup generation?
\end{itemize}
\end{frame}
\begin{frame}
\frametitle{Climate change scenarios: RCP 8.5}
Take a simulated dataset of how the weather would look
between today and the year 2100 with a scenario of high
concentrations of greenhouse gases.
The scenario is called Representative Concentration Pathways 8.5
(RCP~8.5), since it estimates a radiative forcing of $\Delta P = 8.5$~
$\mathrm{W}/ \mathrm{m}^2$ (difference between insolation and energy
radiated into space) at the end of the century. It is a \alert{worst-case scenario} and extrapolates current greenhouse gas emissions
without reduction efforts (improbable given current trajectories of coal, renewables and EVs). This corresponds to a
CO$_2$-equivalent-concentration (including all forcing agents) of
approximately 1250 ppm (today around 410~pmm for CO$_2$) and an
average temperature increase of $\Delta T = 3.7 \pm 1.1$~C at the
end of the century, dependent on the model used.
Compare historical values (HIS) to begin/middle/end of the century (B/M/EOC).
\end{frame}
\begin{frame}
\frametitle{Changes to wind capacity factors}
Left: historic (HIS) wind capacity factors 1970-2005
Right: change at end of century (EOC) 2070-2100
\begin{columns}[T]
\begin{column}{9.5cm}
\includegraphics[width=10cm]{cc-wind.png}
\end{column}
\begin{column}{4.5cm}
\vspace{0.5cm}
\begin{itemize}
\item Small ($\sim 5\%$) average increase in Northern Europe
\item Small ($\sim 5\%$) average decrease in Southern Europe
\end{itemize}
\end{column}
\end{columns}
%corr_ensemble_sfcWind_abs.eps
\source{\href{https://arxiv.org/abs/1805.11673}{Schlott et al, 2018}}
\end{frame}
\begin{frame}
\frametitle{Changes to solar capacity factors}
\begin{columns}[T]
\begin{column}{9.5cm}
\includegraphics[width=10cm]{cc-solar.png}
\end{column}
\begin{column}{4.5cm}
\vspace{0.5cm}
\begin{itemize}
\item Small ($\sim 5\%$) increase in
in Southern Europe around Mediterranean
\item Smallish ($\sim 10\%$) decrease in Northern Europe (due to
increased cloud cover)
\item Solar results known to be a little unreliable because of
cloud modelling etc.
\end{itemize}
\end{column}
\end{columns}
%mean_ensemble_sfcWind_abs.eps
%corr_ensemble_sfcWind_abs.eps
\source{\href{https://arxiv.org/abs/1805.11673}{Schlott et al, 2018}}
\end{frame}
\begin{frame}
\frametitle{Correlation Length}
\begin{columns}[T]
\begin{column}{5.5cm}
\includegraphics[height=7.5cm]{simeon-correlation}
\end{column}
\begin{column}{8.5cm}
The Pearson correlation coefficient of wind time series with a point in
northern Germany decays exponentially with distance.
Determine the \alert{correlation length} $L$ by fitting the
function:
\begin{equation*}
\rho \sim e^{-\frac{x}{L}}
\end{equation*}
to the radial decay with distance $x$.
\vspace{.2cm}
\centering
\includegraphics[height=4.5cm]{correlation_length.png}
\end{column}
\end{columns}
\source{\href{https://doi.org/10.1016/j.apenergy.2011.10.039}{Hagspiel et al, 2012}}
\end{frame}
\begin{frame}
\frametitle{Changes to wind speed correlation lengths}
\begin{columns}[T]
\begin{column}{9.5cm}
\includegraphics[width=10cm]{cc-correlation.png}
\end{column}
\begin{column}{5.3cm}
\begin{itemize}
\item Correlation lengths are longer in the North than
the South because of big weather systems that roll in from
the Atlantic to the North (in the South they get dissipated).
\item With global warming, correlation lengths grow
longer in the North and shorter in the South.
\item This is because weather systems have more energy and are
bigger in the North.
\end{itemize}
\end{column}
\end{columns}
\source{\href{https://arxiv.org/abs/1805.11673}{Schlott et al, 2018}}
\end{frame}
\begin{frame}
\frametitle{Effects of climate change on power system}
Conclusions from study of effects on the power system:
\begin{itemize}
\item Most effects are small ($\sim 5-10\%$); total system costs increase by only 5\%.
\item Longer correlation lengths see greater benefit from continental transmission.
\item Impact of climate change is of a similar magnitude to the uncertainty between the different weather models.
\item Not considered: Space heating and cooling demand changes may have bigger effect on overall energy system.
\item Not considered: Impact of extreme weather events (storms, fires, droughts).
\end{itemize}
For more results, see `The Impact of Climate Change on a Cost-Optimal Highly Renewable European Electricity Network,' \url{https://arxiv.org/abs/1805.11673}
\end{frame}
%\section{Myopic Foresight for Weather Uncertainty}
%\section{Myopic Foresight for Multi-Decade Investment}
\section{Cost and Political Uncertainty}
\begin{frame}
\frametitle{Power System Model: Sensitivity to Changing Solar Cost}
In 30-node European electricity system with 95\% CO$_2$ reduction, change solar
capital cost relative to default. NB: Even at zero solar cost, there
is still wind. Why? Seasonality.\\
LV 0: No cross-border grid, LV 125: compromise grid, LV Opt: optimal grid.
\centering
\includegraphics[width=11cm]{sensitivity-solar_cost.png}
\source{\href{https://arxiv.org/abs/1803.09711}{Schlachtberger et al, 2018}}
\end{frame}
\begin{frame}
\frametitle{Power System Model: Sensitivity to Onshore Wind Installable Potential}
In electricity system with 95\% CO$_2$ reduction, reduce installable
potential for onshore wind. Onshore substituted with offshore at only small
extra system cost. BUT assumes sufficient grid capacity within each country to get offshore from coast to load.
\centering
\includegraphics[width=11cm]{sensitivity-onshore_wind.png}
\source{\href{https://arxiv.org/abs/1803.09711}{Schlachtberger et al, 2018}}
\end{frame}
\begin{frame}
\frametitle{Sensitivity of Optimisation to Cost, Weather Data and Policy Constraints}
See Schlachtberger et al, `Cost optimal scenarios of a future highly renewable European
electricity system: Exploring the influence of weather data, cost
parameters and policy constraints,' 2018,
\url{https://arxiv.org/abs/1803.09711}
\end{frame}
\section{Near-Optimal Energy Systems}
\begin{frame}
\frametitle{Flat directions near optimum}
Both for changing transmission expansion AND onshore wind installable potentials, we've seen that total system costs are \alert{flat around the optimum}.
Can we explore this \alert{near-optimal space} more systematically?
\vspace{.5cm}
\begin{columns}[T]
\begin{column}{6cm}
\includegraphics[width=6cm]{costs_per_cluster_181}
\end{column}
\begin{column}{8cm}
\includegraphics[width=8cm]{sensitivity-onshore_wind.png}
\end{column}
\end{columns}
\end{frame}
% optimal
\begin{frame}
\frametitle{Large Space of Near-Optimal Energy Systems}
There is a \alert{large degeneracy} of different possible energy systems close to the optimum.
\vspace{.2cm}
\centering
\includegraphics[page=2, trim=0.3cm 16.02cm 0cm 3.05cm, clip, height=0.83\textheight]{graphics/sketch.pdf}
\source{Fabian Neumann}
\end{frame}
% epsilon constraint
\begin{frame}
\frametitle{Large Space of Near-Optimal Energy Systems}
Consider the part of the feasible space within $\varepsilon$ of the optimum $f(x^*)$.
\vspace{.2cm}
\centering
\includegraphics[page=1, trim=0.3cm 0.7cm 0cm 18.35cm, clip, height=0.83\textheight]{graphics/sketch.pdf}
\source{Fabian Neumann}
\end{frame}
% mga
\begin{frame}
\frametitle{Large Space of Near-Optimal Energy Systems}
Now within $\varepsilon$ of the optimum $f(x^*)$, try minimising or maximising $x$, to probe space.
\vspace{.2cm}
\centering
\includegraphics[page=1, trim=0.3cm 15.45cm 0cm 3.6cm, clip, height=0.83\textheight]{graphics/sketch.pdf}
\source{Fabian Neumann}
\end{frame}
% 2-dimensional
\begin{frame}
\frametitle{Large Space of Near-Optimal Energy Systems}
NB: Decision space of variables is multi-dimensional, so can probe only one direction at a time.
\vspace{.2cm}
\centering
\includegraphics[page=1, trim=8.8cm 11.5cm 9.4cm 14.8cm, clip, height=0.83\textheight]{graphics/sketch.pdf}
\source{Fabian Neumann}
\end{frame}
\begin{frame}{Application: Highly-Renewable European Electricity System}
Apply this technique to a 100-node model of the European electricity with 100\% renewable energy.
\begin{enumerate}
\item Find the \alert{least-cost power system}.
~\\
\item For $\varepsilon\in\{0.5,1,\dots,10\} \%$ \alert{minimise/maximise} investment in
\begin{itemize}
\item generation capacity {(onshore and/or offshore wind, solar)},
\item storage capacity {(hydrogen, batteries, total storage)} and
\item transmission volume {(HVAC lines and HVDC links)}
\end{itemize}
such that \alert{total annual system costs increase by less than $\varepsilon$}.
\end{enumerate}
\vspace{.5cm}
Methodology adapted from Method to Generate Alternatives (MGA) but `alternatives' are forced in politically-interesting directions.
\source{Fabian Neumann}
\end{frame}
\begin{frame}
\frametitle{Example: 100\% renewable electricity system for Europe}
\begin{columns}[T]
\begin{column}{7cm}
Capacity expansion in optimum:
\includegraphics[width=7cm]{no-optimum.png}
\end{column}
\begin{column}{7cm}
$\varepsilon =$ 10\% above optimum, minimise new grid:
\includegraphics[width=7cm]{no-trans-10.png}
\end{column}
\end{columns}
\source{\href{https://arxiv.org/abs/1910.01891}{Neumann \& Brown, 2020}}
\end{frame}
\begin{frame}
\frametitle{Example: 100\% renewable electricity system for Europe}
\begin{columns}[T]
\begin{column}{7cm}
\includegraphics[trim=0cm 0cm 0cm 0cm, clip, width=7cm]{graphics/space-00.pdf}
\end{column}
\begin{column}{6.5cm}
Within 10\% of the optimum we can:
\begin{itemize}
\item Eliminate most grid expansion
\item Exclude onshore or offshore wind or PV
\item Exclude battery or most hydrogen storage
\end{itemize}
\vspace{.2cm}
\alert{Robust conclusions}: wind, some transmission, some storage, preferably hydrogen storage, required for a cost-effective solution.
\vspace{.2cm}
This gives space to choose solutions with \alert{higher public acceptance}.
\end{column}
\end{columns}
\source{\href{https://arxiv.org/abs/1910.01891}{Neumann \& Brown, 2020}}
\end{frame}
\begin{frame}[fragile]
\frametitle{Flat directions allow society to choose based on other criteria}
\begin{columns}[T]
\begin{column}{5cm}
\centering
\includegraphics[width=5cm]{nein_zur_monstertrasse}
\end{column}
\begin{column}{6.7cm}
\vspace{.5cm}
This flatness may allow us to choose solutions with \alert{higher
public acceptance} at only \alert{small extra cost}.
\vspace{.5cm}
These trade-offs will occupy us for the next 30 years!
\vspace{.5cm}
\includegraphics[width=7cm]{Protestplakat-gegen-den-Bau-von-Windraedern-in-Hamburg-Deutschland.jpg}
\end{column}
\end{columns}
\end{frame}
\begin{frame}{\textbf{Dependencies:} Extremes cannot be achieved simultaneously}
\begin{columns}
\begin{column}{0.4\textwidth}
\begin{center}
\vspace{-0.5cm}
\footnotesize Optimal System Layout\\
\includegraphics[height=0.8\textheight]{graphics/capacity-bar.pdf}
\end{center}
\end{column}
%\vrule
\begin{column}{0.6\textwidth}
\begin{center}
\vspace{.3cm}
\includegraphics[height=0.9\textheight]{graphics/capacity-bars.pdf}
\end{center}
\end{column}
\end{columns}
\end{frame}
\begin{frame}
\frametitle{Near-Optimal Systems: Conclusions}
\begin{itemize}
\item Optimizing a single model gives a \alert{false sense of exactness}.
\item There are many uncertainties about cost assumptions and political targets.
\item There are also \alert{structural model uncertainties} since the feasible space can be very \alert{flat} near the optimum, such that the solution chosen is random within flat area.
\item We can use these techniques to probe the \alert{near-optimal space}.
\item This gives us fuzzier but \alert{more robust} conclusions (e.g. need wind, some transmission and some long-term storage for a cost-effective solution).
\item It also allows us to find cost-effective solutions with \alert{higher public acceptance}.
\end{itemize}
More details: Fabian Neumann, Tom Brown, ``The Near-Optimal Feasible Space of a Renewable Power System Model,'' 2020, accepted to PSCC 2020, \url{https://arxiv.org/abs/1910.01891}.
\end{frame}
\end{document}