Commits (2)
No preview for this file type
 ... ... @@ -170,8 +170,8 @@ Consider an electricity market with two generator types, one with the cost funct Consider the two-bus power system shown in Figure \ref{twobus}, where the two nodes represent two markets, each with different total demand $D_i$, and one generator at each node producing $P_i$. At node A the demand is $D_A = 2000 \si{\mega\watt}$, whereas at node B the demand is $D_B = 1000 \si{\mega\watt}$. Furthermore, there is a transmission line with a capacity denoted by $F_{AB}$. The marginal cost of production of the generators connected to buses A and B are given respectively by the following expressions: \begin{align*} MC_A & = 20 + 0.03 P_A \hspace{1cm}\eur/\si{\mega\watt\hour} \\ MC_B & = 15 + 0.02 P_B \hspace{1cm} \eur/\si{\mega\watt\hour} MC_A & = 20 + 0.02 P_A \hspace{1cm}\eur/\si{\mega\watt\hour} \\ MC_B & = 15 + 0.03 P_B \hspace{1cm} \eur/\si{\mega\watt\hour} \end{align*} Assume that the demands $D_A$ and $D_B$ are constant and insensitive to price, that energy is sold at its marginal cost of production and that there are no limits on the output of the generators. ... ...
 ... ... @@ -85,12 +85,13 @@ "outputs": [], "source": [ "#make the code as Python 3 compatible as possible\n", "from __future__ import print_function, division, absolute_import\n", "#from __future__ import print_function, division, absolute_import\n", "import pypsa\n", "import numpy as np\n", "import pandas as pd\n", "import os\n", "import matplotlib.pyplot as plt\n", "import cartopy.crs as ccrs\n", "plt.rcParams['figure.dpi']= 400\n", "plt.style.use('ggplot')\n", "%matplotlib inline" ... ... @@ -338,14 +339,14 @@ " n_rows_g = n_graphs_g // n_cols_g\n", "else:\n", " n_rows_g = n_graphs_g // n_cols_g + 1\n", " \n", "fig, axes = plt.subplots(nrows=n_rows_g, ncols=n_cols_g,\n", " figsize=(width_factor*n_cols_g,height_factor*n_rows_g),\n", " subplot_kw={\"projection\":ccrs.PlateCarree()})\n", "\n", "width_factor = 4\n", "height_factor = 4\n", "\n", "fig, axes = plt.subplots(nrows=n_rows_g, ncols=n_cols_g,\n", " figsize=(width_factor*n_cols_g,height_factor*n_rows_g),\n", " subplot_kw={\"projection\":ccrs.PlateCarree()})\n", "\n", "for i,tech in enumerate(techs):\n", " i_row = i // n_cols_g\n", " i_col = i % n_cols_g\n", ... ... @@ -450,7 +451,7 @@ "\n", "fig, ax = plt.subplots(figsize=(16,8))\n", "\n", "p_by_carrier.plot(\n", "p_by_carrier_gw.plot(\n", " kind=\"area\",\n", " ax=ax,\n", " linewidth=0,\n", ... ...