diff --git a/SLEW_Case2_Electromechanical_Phenomena.ipynb b/SLEW_Case2_Electromechanical_Phenomena.ipynb index 69daee33bfbeabc0327c1d0e973d1c24f8d3a54e..6c9e14d3e29d739ce7e27bcd6a0edda6672cdb41 100644 --- a/SLEW_Case2_Electromechanical_Phenomena.ipynb +++ b/SLEW_Case2_Electromechanical_Phenomena.ipynb @@ -2,7 +2,6 @@ "cells": [ { "cell_type": "markdown", - "id": "cf1f4a01-23c9-483d-9bf3-0a80ed440ae0", "metadata": {}, "source": [ "<div>\n", @@ -29,7 +28,6 @@ }, { "cell_type": "markdown", - "id": "a620deac-bea3-4f34-841a-b9d7e41fe1ae", "metadata": { "tags": [] }, @@ -40,7 +38,6 @@ }, { "cell_type": "markdown", - "id": "40672b8b-585f-4448-9dac-523d4f33287c", "metadata": {}, "source": [ "A $200 MVA$ round-rotor generator is connected to an infinite bus system of nominal frequency of $50 Hz$ through a step-up transformer with reactance of $0.13 p.u.$ and a transmission line with reactance of $0.17 p.u.$ . The generator’s transient reactance is $x'_{d} = 0.23 p.u.$ and the inertia constant is $H = 4 s$. The synchoronous generator mechanical power is set to $0.6 p.u.$ and the steady-state emf $E = 1.1 p.u.$" @@ -48,7 +45,6 @@ }, { "cell_type": "markdown", - "id": "ab2be9d0-04eb-47cd-988d-3af3cde391bd", "metadata": { "tags": [] }, @@ -59,7 +55,6 @@ }, { "cell_type": "markdown", - "id": "4947469f-3514-4edc-a9d9-4d3157f36513", "metadata": {}, "source": [ "1. If at the given operating point the system eigen-values are as follows:\n", @@ -75,7 +70,6 @@ }, { "cell_type": "markdown", - "id": "68a26943-3791-4b08-a38e-a4a1454d24b7", "metadata": {}, "source": [ "**You can process the results from VILLASweb (as required in 1 and 2) using the prepared notebook cells below.**" @@ -83,7 +77,6 @@ }, { "cell_type": "markdown", - "id": "8b454ae8-71a1-4017-9c5a-cf3efc494e40", "metadata": { "tags": [] }, @@ -94,7 +87,6 @@ }, { "cell_type": "markdown", - "id": "fa5f07a4-23a0-4846-8fc2-d9a3eeae3b8a", "metadata": {}, "source": [ "#### First Setup your Python for post-processing" @@ -102,7 +94,6 @@ }, { "cell_type": "markdown", - "id": "032049e7-1108-4fb1-9057-3b6c522a81bd", "metadata": {}, "source": [ "Import relevant Python packages by executing the following cell:" @@ -111,7 +102,6 @@ { "cell_type": "code", "execution_count": null, - "id": "fa7d426a-0908-454d-8b0d-0dc64eb785ca", "metadata": {}, "outputs": [], "source": [ @@ -131,7 +121,6 @@ }, { "cell_type": "markdown", - "id": "9be353cf-624a-4e88-a427-f94091c651bb", "metadata": {}, "source": [ "### Subtask 1.b\n", @@ -141,7 +130,6 @@ { "cell_type": "code", "execution_count": null, - "id": "cd8c51d0-ca49-426b-8b86-af774b6717b7", "metadata": { "tags": [] }, @@ -153,7 +141,6 @@ }, { "cell_type": "markdown", - "id": "ab551c32-ef6f-47d3-b839-57e2ce506863", "metadata": {}, "source": [ "#### Download simulation results from VILLASweb" @@ -161,7 +148,6 @@ }, { "cell_type": "markdown", - "id": "e70b8ccc-8d97-4827-ba45-251be9c6ac82", "metadata": {}, "source": [ "Enter VILLASweb under https://slew.rwth-aachen.de/. \n", @@ -172,7 +158,6 @@ { "cell_type": "code", "execution_count": null, - "id": "5c6f4d63-2c1c-49aa-acde-073f1cffd7d5", "metadata": {}, "outputs": [], "source": [ @@ -182,7 +167,6 @@ }, { "cell_type": "markdown", - "id": "85d62e17-9df8-421c-8b99-aac7be09bf09", "metadata": { "tags": [] }, @@ -193,7 +177,6 @@ { "cell_type": "code", "execution_count": null, - "id": "ec357284-6407-41bb-9d54-3ca6c0532fa9", "metadata": {}, "outputs": [], "source": [ @@ -233,7 +216,6 @@ }, { "cell_type": "markdown", - "id": "87835d8e-afcd-4a23-8b36-596eaaca7b92", "metadata": { "tags": [] }, @@ -245,7 +227,6 @@ { "cell_type": "code", "execution_count": null, - "id": "56d9f0db-b158-4fc3-a9b0-5a2d8931da9b", "metadata": {}, "outputs": [], "source": [ @@ -255,7 +236,6 @@ }, { "cell_type": "markdown", - "id": "0848fc38-50ad-421a-8442-9f8b9c26a1e0", "metadata": { "tags": [] }, @@ -265,7 +245,6 @@ }, { "cell_type": "markdown", - "id": "3c0fd278-b58a-4af7-a101-4a52e1a32c94", "metadata": {}, "source": [ "#### Download simulation results from VILLASweb" @@ -273,7 +252,6 @@ }, { "cell_type": "markdown", - "id": "af0ec7ca-9cf7-4986-9ae9-ab762a792b5f", "metadata": {}, "source": [ "Enter VILLASweb under https://slew.rwth-aachen.de/. \n", @@ -284,7 +262,6 @@ { "cell_type": "code", "execution_count": null, - "id": "8331e105-8b66-4479-9caf-6ec038c1bce0", "metadata": {}, "outputs": [], "source": [ @@ -294,7 +271,6 @@ }, { "cell_type": "markdown", - "id": "f9879ce7-72ae-4306-9d65-8f8604484c43", "metadata": {}, "source": [ "#### Plot the results" @@ -303,7 +279,6 @@ { "cell_type": "code", "execution_count": null, - "id": "ccad2758-d711-46c9-a54b-a788efb563eb", "metadata": {}, "outputs": [], "source": [ @@ -328,16 +303,16 @@ "plt.subplot(2,1,1)\n", "plt.ylabel(\"Angle (degree)\")\n", "plt.xlabel(\"time (s)\")\n", - "plt.plot(ts_res1['delta_r_gen'].time, ts_res1['delta_r_gen'].values, label='Rotor Angle task 1.b', color='C0')\n", - "plt.plot(ts_res2['delta_r_gen'].time, ts_res2['delta_r_gen'].values, label='Rotor Angle task 2.b', color='C1', linestyle=':')\n", + "plt.plot(ts_res2['delta_r_gen'].time, ts_res2['delta_r_gen'].values, label='Rotor Angle task 2.b', color='C1')\n", + "plt.plot(ts_res1['delta_r_gen'].time, ts_res1['delta_r_gen'].values, label='Rotor Angle task 1.b', color='C0', linestyle=':')\n", "plt.xlim([1, 30])\n", "plt.legend(loc='upper right')\n", " \n", "plt.subplot(2,1,2)\n", "plt.ylabel(\"Power (MW)\")\n", "plt.xlabel(\"time (s)\")\n", - "plt.plot(ts_res1['P_elec'].time, ts_res1['P_elec'].values/1e6, label='Output Power task 1.b', color='C0')\n", - "plt.plot(ts_res2['P_elec'].time, ts_res2['P_elec'].values/1e6, label='Output Power task 2.b', color='C1', linestyle=':') \n", + "plt.plot(ts_res2['P_elec'].time, ts_res2['P_elec'].values/1e6, label='Output Power task 2.b', color='C1') \n", + "plt.plot(ts_res1['P_elec'].time, ts_res1['P_elec'].values/1e6, label='Output Power task 1.b', color='C0', linestyle=':')\n", "plt.xlim([1, 30])\n", "plt.legend(loc='upper right')\n", "plt.show()" @@ -345,7 +320,6 @@ }, { "cell_type": "markdown", - "id": "8e4f8736-d4a5-4d5b-8d56-fee4b7afb499", "metadata": { "tags": [] }, @@ -355,7 +329,6 @@ }, { "cell_type": "markdown", - "id": "88f14679-151b-47d3-bf9e-0c912427e504", "metadata": {}, "source": [ "#### Download simulation results from VILLASweb" @@ -363,7 +336,6 @@ }, { "cell_type": "markdown", - "id": "0740d062-5e8f-4eb9-88bc-a6eea0003391", "metadata": {}, "source": [ "Enter VILLASweb under https://slew.rwth-aachen.de/. \n", @@ -374,7 +346,6 @@ { "cell_type": "code", "execution_count": null, - "id": "e490e568-1db0-4fb6-942a-30fafecee4bd", "metadata": {}, "outputs": [], "source": [ @@ -384,7 +355,6 @@ }, { "cell_type": "markdown", - "id": "40ee9509-a952-42d8-a67d-4901ba2d95bd", "metadata": {}, "source": [ "#### Plot the results\n" @@ -393,7 +363,6 @@ { "cell_type": "code", "execution_count": null, - "id": "36505e96-7b81-41df-847c-a024269beb05", "metadata": {}, "outputs": [], "source": [ @@ -418,25 +387,32 @@ "plt.subplot(2,1,1)\n", "plt.ylabel(\"Angle (degree)\")\n", "plt.xlabel(\"time (s)\")\n", - "plt.plot(ts_res2['delta_r_gen'].time, ts_res2['delta_r_gen'].values, label='Rotor Angle task 2.b', color='C1')\n", - "plt.plot(ts_res3['delta_r_gen'].time, ts_res3['delta_r_gen'].values, label='Rotor Angle task 2.c', color='C2', linestyle=':')\n", + "plt.plot(ts_res3['delta_r_gen'].time, ts_res3['delta_r_gen'].values, label='Rotor Angle task 2.c', color='C2')\n", + "plt.plot(ts_res2['delta_r_gen'].time, ts_res2['delta_r_gen'].values, label='Rotor Angle task 2.b', color='C1', linestyle=':')\n", "plt.xlim([1, 30])\n", "plt.legend(loc='upper right')\n", " \n", "plt.subplot(2,1,2)\n", "plt.ylabel(\"Power (MW)\")\n", "plt.xlabel(\"time (s)\")\n", - "plt.plot(ts_res2['P_elec'].time, ts_res2['P_elec'].values/1e6, label='Output Power task 2.b', color='C1') \n", - "plt.plot(ts_res3['P_elec'].time, ts_res3['P_elec'].values/1e6, label='Output Power task 2.c', color='C2', linestyle=':') \n", + "plt.plot(ts_res3['P_elec'].time, ts_res3['P_elec'].values/1e6, label='Output Power task 2.c', color='C2')\n", + "plt.plot(ts_res2['P_elec'].time, ts_res2['P_elec'].values/1e6, label='Output Power task 2.b', color='C1', linestyle=':') \n", "plt.xlim([1, 30])\n", "plt.legend(loc='upper right')\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { "kernelspec": { - "display_name": "Python 3 (ipykernel)", + "display_name": "Python 3", "language": "python", "name": "python3" }, @@ -450,7 +426,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.9.6" + "version": "3.7.7" } }, "nbformat": 4,