Multi-Objective Optimal Dispatch Considering Wind Power and Interactive Load for Power System ()

Xinxin Shi^{1,2,3}, Guangqing Bao^{1,2,3*}, Kun Ding^{4}, Liang Lu^{4}

^{1}College of Electrical and Information Engineering, Lanzhou University of Technology, Lanzhou, China.

^{2}Key Laboratory of Gansu Advanced Control for Industrial Processes, Lanzhou University of Technology, Lanzhou, China.

^{3}National Demonstration Center for Experimental Electrical and Control Engineering Education, Lanzhou University of
Technology, Lanzhou, China.

^{4}Key Laboratory of Wind Power Integration Operation and Control, Gansu Electric Power Corporation Wind Power Corpora-tion Wind Power Technology Center, Lanzhou, China.

**DOI: **10.4236/epe.2018.104B001
PDF HTML XML
978
Downloads
1,628
Views
Citations

With the rapid and large-scale development of renewable energy, the lack of new energy power transportation or consumption, and the shortage of grid peak-shifting ability have become increasingly serious. Aiming to the severe wind power curtailment issue, the characteristics of interactive load are studied upon the traditional day-ahead dispatch model to mitigate the influence of wind power fluctuation. A multi-objective optimal dispatch model with the minimum operating cost and power losses is built. Optimal power flow distribution is available when both generation and demand side participate in the resource allocation. The quantum particle swarm optimization (QPSO) algorithm is applied to convert multi-objective optimization problem into single objective optimization problem. The simulation results of IEEE 30-bus system verify that the proposed method can effectively reduce the operating cost and grid loss simultaneously enhancing the consumption of wind power.

Keywords

Wind Power, Interactive Load, Optimal Dispatch, Multi-Objective, QPSO Models

Share and Cite:

Shi, X. , Bao, G. , Ding, K. and Lu, L. (2018) Multi-Objective Optimal Dispatch Considering Wind Power and Interactive Load for Power System. *Energy and Power Engineering*, **10**, 1-10. doi: 10.4236/epe.2018.104B001.

1. Introduction

Wind power industry has been rapidly developed. Wind power generation reached 241 billion kW∙h that has been 30% year-on-year growth accounting for 4% of the total electricity generation in China. Among them, new energy installed capacity accounted for more than 30% of the total installed capacity of the local power supply in Gansu, Ningxia, Xinjiang, Qinghai, which has shown favourable prospects in reducing fossil energy consumption and pollutant emissions [1]. However, the demand of electricity is influenced by variable factors, such as weather, economy, laws, policies, electrical load conditions, and so on. These factors make electric dispatch became a task. Especially, when the wind power integrated into the main grid, the intermittent and stochastic nature of such energy brings challenges to system dispatch. Therefore, usage of wind power generation is of major importance in future power grids for economic and environmental reasons. In recent years, the researches of academic experts published which classified into two major categories, one is multi-objective dispatch modeling; another one is optimal dispatch solving technology. A stochastic multi-objective optimal reactive power dispatch model is studied concerning about load and wind power generation uncertainties, including real power losses and operation cost of wind farms [2]. The objective function is implemented for energy saving and emission-reduction [3]. The day-ahead multi-objective optimal dispatching model containing thermal power, hydro power, wind-power and pumped storage units is given to minimize the total costs and CO_{2} emission under multiple constraints. A fuzzy modeling for dynamic economic dispatch is presented in [4], which could make the dispatch to reflect the willingness of decision-makers and hereby adapt the random wind power output better. A multi-objective optimization algorithm based on the non-dominated sorting differential evolution is used to solve the economic environmental dispatch stochastic optimization model of power system connected with large scale wind farms [5]. And a solution of optimal dispatch problem with a particle swarm optimisation based on multi-agent systems is presented in [6]. Optimal scheduling strategy is considered only from the generation side above literature. In this paper, the characteristics of interaction load and integration into the traditional scheduling model is discussed. The problem is formulated as a multi-objective optimal problem through simultaneous minimizing both system operational cost and power loss. QPSO algorithm has been introduced to solve this problem. Finally, the system simulation is carried out with IEEE 30 nodes. The experimental results verify the effectiveness of the proposed method, and have some certain practical significance for power system optimal scheduling.

2. Interactive Load Modeling

1) Interactive load characteristics

2) Interactive load model

・ Shiftable load

The dispatching center calculates the optimal power dispatching plan according to the information provided by the intention chart which determine the optimal dispatch time of the load user and the shiftable load involved in shifting the peak. Therefore, for the rth shiftable load, the decision variable is the start variable ${U}_{sr}^{t}$ : If it starts in the t period, then ${U}_{sr}^{t}=1$ ; If it starts at other times, then ${U}_{sr}^{t}=0$ ; ${t}_{sr}$ is the start time of the rth shiftable peak load [7] [8].

The shift cost curve characterizes denotes the compensation price that the user should obtain from the grid company after providing the shift service. The mathematical description is as follows:

${C}_{sr}^{t}=\{\begin{array}{l}{m}_{sr}({t}_{s{r}_{0}}-t)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\le t<{t}_{s{r}_{0}}\\ {m}_{sr}(t-{t}_{s{r}_{0}})-{m}_{sr}{T}_{sro}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>{t}_{s{r}_{0}}+{T}_{s{r}_{0}}\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{t}_{s{r}_{0}}\le t\le {t}_{s{r}_{0}}+{T}_{s{r}_{0}}\end{array}$ (1)

where is ${C}_{sr}^{t}$ the peak cost of the rth shiftable peak load in the t period of the peak; ${m}_{sr}$ is the compensation factor for the load of the rth shiftable peak load that is determined by the load control agreement signed from the user and the power company in advance; ${t}_{s{r}_{0}}$ and ${T}_{s{r}_{0}}$ are the number of the original power load start time and load duration before the peak which belong to known parameters.

・ Interruptible load

Users through the auction to declare interruptible capacity and compensation prices as well as the dispatch center by calculating the optimal power generation scheduling program to determine the interruptible users and the optimal capacity. Dispatching the hth interruptible load of the user’s compensation cost as shown in Equation (4) [9]^{ }

${C}_{Ih}^{t}={C}_{Ih0}{P}_{Ih}^{t}{I}_{Ih}^{t}$ (2)

where ${C}_{Ih0}$ is the unit to reduce load costs of the hth interruptible load in the contract; ${P}_{Ih}^{t}$ is the load reduction of the hth interruptible load in the t period ; and ${I}_{Ih}^{t}$ is variable dispatch for interruptible load. The interruptible load was whether to dispatched basing on ${I}_{Ih}^{t}=1$ and ${I}_{Ih}^{t}=0$ .

3. Multi-Objective Optimal Statement

The problem of economic power dispatch with wind penetration consideration can be formulated as a multi-objective optimal dispatch model. The two conflicting objectives, i.e., operating cost and system power loss, should be minimized simultaneously while fulfilling certain system constraints. This problem is formulated mathematically in this section.

3.1. Problem Objectives

・ Objective 1: Minimization of operational cost

$\mathrm{min}{C}_{1}={C}_{Gi}+{C}_{Wi}$ (3)

${C}_{Gi}={\displaystyle \underset{t=1}{\overset{T}{\sum}}{\displaystyle \underset{i=1}{\overset{{N}_{c}}{\sum}}\{{U}_{Gi}^{t}[{a}_{i}+{b}_{i}{P}_{Gi}^{t}+{c}_{i}{({P}_{Gi}^{t})}^{2}]+{C}_{Ui}^{t}(1-{U}_{Gi}^{t-1})+{C}_{Ri}^{t}\}}}$ (4)

${C}_{Wi}={\displaystyle \underset{t=1}{\overset{T}{\sum}}({\displaystyle \underset{j=1}{\overset{{N}_{S}}{\sum}}{C}_{Sr}^{t}+{\displaystyle \underset{i=1}{\overset{{N}_{I}}{\sum}}{C}_{Ih}^{t}}})}$ (5)

where T is the number of hours during system dispatching; ${N}_{C}$ is the number of generating units; and ${P}_{Gi}^{t}$ , ${U}_{Gi}^{t}$ are the active output and the state variables of the ith generator in the period; indicates a shutdown state. On the contrary, ${U}_{Gi}^{t}=1$ indicates that the generator set is on; ${C}_{Ri}^{t}$ , ${C}_{Ui}^{t}$ denote the spare costs and starting cost of the ith generator in the period; ${N}_{S}$ , ${N}_{I}$ represent the number of shiftable peak load; ${C}_{Sr}^{t}$ , ${C}_{Ih}^{t}$ denote the cost of the rth removable peak load and the hth interruptible load in the t period.

・ Objective 2: Minimization of system power loss

The dispatch of interactive load will inevitably cause the change of power flow distribution, which will have some influence on the system power loss. Thus, the minimize of power loss is one goal of optimal dispatching. Here the use of B-coefficient method to calculate the power loss [10].

$\mathrm{min}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{2}={\displaystyle \underset{t=1}{\overset{T}{\sum}}({\displaystyle \underset{i=1}{\overset{K}{\sum}}{\displaystyle \underset{j=1}{\overset{K}{\sum}}{P}_{i}^{t}{B}_{i,j}^{}{P}_{j}^{t}+{\displaystyle \underset{i=1}{\overset{K}{\sum}}{B}_{i,0}^{}{P}_{i}^{t}+{B}_{0,0}}}})}$ (6)

where K is the number of system nodes; ${B}_{i,j}$ , ${B}_{i,0}$ , ${B}_{0,0}$ are the second term , which are the first term and the constant term of the coefficient; and ${P}_{i}^{t}$ , ${P}_{j}^{t}$ are active power of node i and j.

3.2. Problem Constraints

・ Constraint 1: Power balance constraint [11]

・ Constraint 2: Spare constraints

・ Constraint 3: Unit constraint

・ Constraint 4: Interactive load constraint [12]

3.3. Optimal Dispatch Modeling

The objective of this model is to minimize the operating cost and the grid loss as much as possible under all constraints. Therefore, when the operating cost and the network loss are lower, the fitness value of the fitness function is greater. Where the fitness function can be defined as [13]:

$\mu ({C}_{i})=\{\begin{array}{l}1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{i}\le {C}_{x}\\ {C}_{i}^{2}+m{C}_{i}+n\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{x}<{C}_{i}\le {C}_{x}+\Delta {C}_{i}\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{i}>{C}_{x}+\Delta {C}_{i}\end{array}$ (7)

$m=(-2{C}_{x}\Delta {C}_{i}-\Delta {C}_{i}^{2}-1)/\Delta {C}_{i}$ (8)

$n=1-{C}_{x}^{2}-{C}_{x}\cdot [(-2{C}_{x}\Delta {C}_{i}-\Delta {C}_{i}^{2}-1)/\Delta {C}_{i}]$ (9)

$\Delta {C}_{i}={C}_{i}-{C}_{x}$ (10)

where ${C}_{i}$ is the ith objective function value; ${C}_{x}$ is the ith objective function ideal value. $\Delta {C}_{i}$ is the ith objective function added value. The fitness function diagram is shown in Figure 1.

Thus, where the fitness index can be defined as:

$\mu =\mathrm{min}\left\{\mu ({C}_{1}),\mu ({C}_{2})\right\}$ (11)

Figure 1. Fitness function diagram.

where m is the minimum value for all fitness functions.

The multi-objective problem is transformed into a single objective nonlinear optimization problem that satisfies the fitness value of all constraints:

$\begin{array}{l}\mathrm{max}\text{\hspace{0.17em}}\mu \\ s.t.\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{1}+\mu \Delta {C}_{1}\le {C}_{1}+\Delta {C}_{1}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{C}_{2}+\mu \Delta {C}_{2}\le {C}_{2}+\Delta {C}_{2}\end{array}$ (12)

4. Proposed Approach

In this paper, the quantum particle swarm optimization (QPSO) [14] is used to solve the model, which is based on the particle swarm optimization (PSO) to improve the formation. All particles in the population are treated as quantum particles in the feasible solution. When updating the particle position, it mainly considers the current local optimal position and global optimal position information of each particle.

5. Simulation and Evaluation

In this study, a IEEE 30-bus system with 1-wind farm of grid-connected is used to investigate the effectiveness of the model. The system configuration is shown in Figure 2.

The system parameters including generator capacities, spare cost and fuel cost coefficients are listed in Table 1.

The interruptable capacity, compensation price, interruptable times and duration for interruptible load are listed in Table 2.

1) There are 100 wind turbines in the wind farm with a total installed capacity of 200 MW. Conventional unit positive and negative rotation standby demand for the maximum unit capacity of 15%.

2) The particle size scale of QPSO algorithm is 100 and the maximum number of iterations is 500.

3) The prediction curves and load forecast curves of the wind power during the last 24 hours are shown in Figure 3.

4) The willingness curve and the peak cost of the shiftable load are shown in Figure 4 and Figure 5, respectively.

It can be seen from Figure 3 that there is obvious anti-peaking characteristic of wind power and load side. The three interruptible loads of this example are at nodes 7, 19, 21 in the interactive load respectively and where three shiftable loads are at nodes 10, 12, 29.

Figure 2. Wiring diagram of IEEE 30-bus system.

Figure 3. Day-ahead wind power and load prediction curve.

Table 1. Data of conventional generator.

Table 2. Data of interruptible loads.

Table 3. Results of multi-objective optimization.

Figure 4. Shiftable load curves.

Figure 5. Cost curve of shiftable loads.

Figure 4 and Figure 5 are, respectively, the time and peak load cost curve. Through reasonable optimization calculations, interruptible load 1, 2, 3 and shiftable load 1, 2 are put into the system; although the shiftable load capacity of the load 3 is relatively large, but because of its shiftable capacity and peak surge costs will be increased system operating cost. It does not meet the economic needs of the operation of the grid, so it does not work. It can be seen from Table 3, the system operating cost of $637030 when the equal to 0.826 which the net loss is 129.47 MW. It can be seen that the power system which integrates the interactive load can provide more reserve for the power grid with wind farm and it reduce the impact of wind power randomness. While reducing the switching costs of interactive load make the power grid more economical and reliable operation.

The results are compared using the conventional optimal scheduling model and the interactive load multi-objective optimal scheduling model. Figure 6 and Figure 7 are the two models of the output curve, respectively.

It can be seen from Figure 6 and Figure 7 that the traditional optimal scheduling model usually need to reduce the output of the conventional unit in the wind power generation period or increasing the output of the conventional unit in the lesser period of wind power. This leads that the conventional unit's output curve isn’t better. Thus, the normal unit 1 and 2 must be worked to meet the system peak demand. But the interactive load multi-objective optimal scheduling model reduces the valley difference of the system. Therefore, it is not necessary that make the normal unit 2 to work.

In the following, the results of the traditional optimization scheduling model and the multi-objective optimization scheduling model considering the interactive

Figure 6. Output curve of traditional scheduling model.

Figure 7. Output curve of optimal dispatching model.

Table 4. Total operational cost of system.

Table 5. Grid loss of system.

load are compared in Table 4 and Table 5.

It can be seen from Table 4 and Table 5 that the shiftable load is better for shifting peak in power system. It is put into the system scheduling system, which can improve the operational efficiency of power systems and decrease the grid loss of system. The generating units are run at a cost-effective level by reducing the frequent shutdown of the generating units. In addition, the interruption load can also reduce the cost of spare capacity due to the random fluctuations in wind power. Compared with the traditional optimal scheduling model, the total system power generation cost and grid loss of system are reduced.

6. Conclusion

The operation cost of power grid will become increased when the capacity of wind power which have fluctuation and randomness is increased. The solution of this problem is the interactive load which is adjusted the power grid to interact with user. The interactive load decreased the operation cost because it can be a part of the reserve capacity of electric power system. In addition, on the basis that bring the interaction load into the optimalscheduling of the system which include the wind farm, a multi-objective optimization scheduling model based on the minimum operation cost and network loss are established. A quantum particle swarm optimization (QPSO) algorithm is utilized to solve the optimization objective of the model which achieve minimum of system power loss and operational cost. Simulation results show the interactive load’s ability to reduce the cost of reserve capacity due to the random fluctuations in wind power. Experimental results demonstrated that the model was successfully implemented.

Acknowledgements

The authors would like to thank the financial support from National Natural Science Foundation of China (51267011), and partly financed by Gansu Electric Power Research Institute (271733-KJ-07).

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Sun, B., Yu, Y.X. and Qin, C. (2017) Should China Focus on the Distributed Development of Wind and Solar Photovoltaic Power Generation? A Comparative Study. Applied Energy, 185, 421-439. https://doi.org/10.1016/j.apenergy.2016.11.004 |

[2] | Mohseni-Bonab, S.M. and Rabiee, A. (2017) Optimal Reactive Power Dispatch: A Review and a New Stochastic Voltage Stability Constrained Multi-Objective Model at the Presence of Uncertain Wind Power Generation. IET Generation, Transmission & Distribution, 11, 815-829. https://doi.org/10.1049/iet-gtd.2016.1545 |

[3] | Shi, N., Zhou, S.Q., Su, X.W., et al. (2016) Unit Commitment and Multi-Objective Optimal Dispatch Model for Wind-Hydro-Thermal Power System with Pumped Storage. 2016 IEEE 8th International Power Electronics and Motion Control Conference (IPEMC-ECCE Asia). |

[4] | Chen, H.Y., Chen, J.F. and Duan, X.Z. (2009) Fuzzy Modeling and Optimization Algorithm on Dynamic Economic Dispatch in Wind Power Integrated System. Proceeding of the CSEE, 29, 13-17. |

[5] | Sun, H.J., Peng, C.H. and Yi, H.J. (2012) Multi-Objective Stochastic Optimal Dispatch of Power System. Electric Power Automation Equipment, 32, 123-128. |

[6] | Zhao, B., Guo, C. and Cao, Y. (2005) A Multiagent-Based Particle Swarm Optimization Approach for Optimal Reactive Power Pispatch. IEEE Trans. Power Syst, 20, 1070-1078. https://doi.org/10.1109/TPWRS.2005.846064 |

[7] | Graditi, G., Di Silvestre, M.L., et al. (2014) Heuristic-Based Shiftable Load Optimal Man-agement in Smart Micro-Grids. IEEE Transactions on Industrial Informatics, 11, 271-280. https://doi.org/10.1109/TII.2014.2331000 |

[8] | Cakmak, R. and Altas, I.H. (2016) Scheduling of Domestic Shiftable Loads via Cuckoo Search Optimization Algorithm. 2016 4th International Istanbul Smart Grid Congress and Fair (ICSG). |

[9] | Wang, Q.Q., Zhao, C.H., Ma, C.F., et al. (2011) Scheduling of Interruptible Load Based on Binary Particle Swarm Optimization. Journal of Nanjing Normal University Engineering and Technology Edition, 11, 19-25. |

[10] | Guo, Z.Z., Wu, J.K. and Kong, F. (2013) Multi-Objective Optimization Scheduling for Hydrothermal Power Systems Based on Electromagnetism-Like Mechanism and Data Envelopment Analysis. Proceedings of the CSEE, 33, 53-61. |

[11] | Suzuki, H. (2014) Dynamics of Load Balancing with Constraints. The European Physical Journal: Special Topics, 223, 2631-2635. https://doi.org/10.1140/epjst/e2014-02278-7 |

[12] | Wang, B.B., Liu, X.C. and Li, Y. (2013) Day-Ahead Generation Scheduling and Operation Simulation Considering Demand Response in Large-Capacity Wind Power Integrated Systems. Proceedings of the CSEE, 33, 35-44. |

[13] | Wu, J.K, and Tang, L. (2012) Stochastic Optimization Scheduling Method for Hydrothermal Power Systems with Stochastic Loads. Proceedings of the CSEE, 32, 36-43. |

[14] | Sun, J. and Xu, W.B. (2004) A Global Search Strategy of Quantum-Behaved Particle Swarm Optimization. Proceedings of IEEE Conference on Cybernetics and Intelligent Systems, Singapore, 111-116. |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2022 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.