NiiloKantoniemi/
EnergySystem
(1.0.0)
current version: 1.0.1 »
Renewable Energy System Modelling
Contributed by: Niilo Kantoniemi
Installation Instructions
To install this paclet in your Wolfram Language environment,
evaluate this code:
PacletInstall["NiiloKantoniemi/EnergySystem"]
To load the code after installation, evaluate this code:
Needs["NiiloKantoniemi`EnergySystem`"]
Examples
Publisher
Niilo Kantoniemi
Compatibility
Wolfram Language Version 14
Version History
-
1.0.1
– 09 July 2024
-
1.0.0
– 09 July 2024
MIT License
Paclet Source
![Needs["NiiloKantoniemi`EnergySystem`"];
startDate = DateObject[{2022, 5, 1}, "Instant"];
endDate = DateObject[{2022, 5, 15}, "Instant"];
step = Quantity[1, "Hour"];(*Quantity[1,"Hour"]*)(*Quantity[1,"Day"]*)
datesWithStep = DateRange[startDate, endDate, step];
datesNoStep = {startDate, endDate};
paris = Entity["City", {"Paris", "IleDeFrance", "France"}];
vaasa = Entity["City", {"Vaasa", "Ostrobothnia", "Finland"}];](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/622e7c24b23e1a4b.png)
![windTS = GenerateWindSpeedTimeSeries[vaasa, datesNoStep];
solarTS = GenerateSolarIrradianceTimeSeries[paris, datesNoStep, datesWithStep];
waterTS = GenerateWaterFlowTimeSeries[Quantity[0.2, ("Meters")^3/("Seconds")],
Quantity[0.1, ("Meters")^3/("Seconds")], startDate, endDate, step];](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/0b4bb61db73a4f71.png)
![w = WindSource[windTS, Quantity[1.20457518249314, ("Kilograms")/("Meters")^3], Quantity[30, "Meters"], 0.40, Quantity[2, "Megawatts"], 2];
s = SolarSource[solarTS, Quantity[4000, ("Meters")^2], 0.20];
n = InterpretationBox[FrameBox[TagBox[TooltipBox[PaneBox[GridBox[List[List[GraphicsBox[List[Thickness[0.0025`], List[FaceForm[List[RGBColor[0.9607843137254902`, 0.5058823529411764`, 0.19607843137254902`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3]]], List[List[List[205.`, 22.863691329956055`], List[205.`, 212.31669425964355`], List[246.01799774169922`, 235.99870109558105`], List[369.0710144042969`, 307.0436840057373`], List[369.0710144042969`, 117.59068870544434`], List[205.`, 22.863691329956055`]], List[List[30.928985595703125`, 307.0436840057373`], List[153.98200225830078`, 235.99870109558105`], List[195.`, 212.31669425964355`], List[195.`, 22.863691329956055`], List[30.928985595703125`, 117.59068870544434`], List[30.928985595703125`, 307.0436840057373`]], List[List[200.`, 410.42970085144043`], List[364.0710144042969`, 315.7036876678467`], List[241.01799774169922`, 244.65868949890137`], List[200.`, 220.97669792175293`], List[158.98200225830078`, 244.65868949890137`], List[35.928985595703125`, 315.7036876678467`], List[200.`, 410.42970085144043`]], List[List[376.5710144042969`, 320.03370475769043`], List[202.5`, 420.53370475769043`], List[200.95300006866455`, 421.42667961120605`], List[199.04699993133545`, 421.42667961120605`], List[197.5`, 420.53370475769043`], List[23.428985595703125`, 320.03370475769043`], List[21.882003784179688`, 319.1406993865967`], List[20.928985595703125`, 317.4896984100342`], List[20.928985595703125`, 315.7036876678467`], List[20.928985595703125`, 114.70369529724121`], List[20.928985595703125`, 112.91769218444824`], List[21.882003784179688`, 111.26669120788574`], List[23.428985595703125`, 110.37369346618652`], List[197.5`, 9.87369155883789`], List[198.27300024032593`, 9.426692008972168`], List[199.13700008392334`, 9.203690528869629`], List[200.`, 9.203690528869629`], List[200.86299991607666`, 9.203690528869629`], List[201.72699999809265`, 9.426692008972168`], List[202.5`, 9.87369155883789`], List[376.5710144042969`, 110.37369346618652`], List[378.1179962158203`, 111.26669120788574`], List[379.0710144042969`, 112.91769218444824`], List[379.0710144042969`, 114.70369529724121`], List[379.0710144042969`, 315.7036876678467`], List[379.0710144042969`, 317.4896984100342`], List[378.1179962158203`, 319.1406993865967`], List[376.5710144042969`, 320.03370475769043`]]]]], List[FaceForm[List[RGBColor[0.5529411764705883`, 0.6745098039215687`, 0.8117647058823529`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[44.92900085449219`, 282.59088134765625`], List[181.00001525878906`, 204.0298843383789`], List[181.00001525878906`, 46.90887451171875`], List[44.92900085449219`, 125.46986389160156`], List[44.92900085449219`, 282.59088134765625`]]]]], List[FaceForm[List[RGBColor[0.6627450980392157`, 0.803921568627451`, 0.5686274509803921`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[355.0710144042969`, 282.59088134765625`], List[355.0710144042969`, 125.46986389160156`], List[219.`, 46.90887451171875`], List[219.`, 204.0298843383789`], List[355.0710144042969`, 282.59088134765625`]]]]], List[FaceForm[List[RGBColor[0.6901960784313725`, 0.5882352941176471`, 0.8117647058823529`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[200.`, 394.0606994628906`], List[336.0710144042969`, 315.4997024536133`], List[200.`, 236.93968200683594`], List[63.928985595703125`, 315.4997024536133`], List[200.`, 394.0606994628906`]]]]]], List[Rule[BaselinePosition, Scaled[0.15`]], Rule[ImageSize, 10], Rule[ImageSize, 15]]], StyleBox[RowBox[List["NuclearGeneration", " "]], Rule[ShowAutoStyles, False], Rule[ShowStringCharacters, False], Rule[FontSize, Times[0.9`, Inherited]], Rule[FontColor, GrayLevel[0.1`]]]]], Rule[GridBoxSpacings, List[Rule["Columns", List[List[0.25`]]]]]], Rule[Alignment, List[Left, Baseline]], Rule[BaselinePosition, Baseline], Rule[FrameMargins, List[List[3, 0], List[0, 0]]], Rule[BaseStyle, List[Rule[LineSpacing, List[0, 0]], Rule[LineBreakWithin, False]]]], RowBox[List["PacletSymbol", "[", RowBox[List["\"NiiloKantoniemi/EnergySystem\"", ",", "\"NiiloKantoniemi`EnergySystem`NuclearGeneration\""]], "]"]], Rule[TooltipStyle, List[Rule[ShowAutoStyles, True], Rule[ShowStringCharacters, True]]]], Function[Annotation[Slot[1], Style[Defer[PacletSymbol["NiiloKantoniemi/EnergySystem", "NiiloKantoniemi`EnergySystem`NuclearGeneration"]], Rule[ShowStringCharacters, True]], "Tooltip"]]], Rule[Background, RGBColor[0.968`, 0.976`, 0.984`]], Rule[BaselinePosition, Baseline], Rule[DefaultBaseStyle, List[]], Rule[FrameMargins, List[List[0, 0], List[1, 1]]], Rule[FrameStyle, RGBColor[0.831`, 0.847`, 0.85`]], Rule[RoundingRadius, 4]], PacletSymbol["NiiloKantoniemi/EnergySystem", "NiiloKantoniemi`EnergySystem`NuclearGeneration"], Rule[Selectable, False], Rule[SelectWithContents, True], Rule[BoxID, "PacletSymbolBox"]][
Quantity[0.3, "Megawatts"]];
vP = VariableProduction[Quantity[0.3, "Megawatts"], Quantity[0.1, "Megawatts"]];
aP = AnyProduction[powerProdTS];
c = ConstantDemand[Quantity[0.25, "Megawatts"]];
v = VariableDemand[Quantity[0.3, "Megawatts"], Quantity[0.1, "Megawatts"]];
aD = AnyDemand[powerDemTS];
b = SimpleBattery[Quantity[150 , "Hours" "Megawatts"], Quantity[5, "Megawatts"], Quantity[5, "Megawatts"], 0.98, 0.00006];
h = HydroPowerSource[Quantity[50, "Meters"], waterTS, Quantity[500000, ("Meters")^3], 0.80, Quantity[2, "Megawatts"], Quantity[9.81, ("Meters")/("Seconds")^2], Quantity[997.0, ("Kilograms")/(
"Meters")^3]];(*0.1(m)^3\[NegativeMediumSpace]/s*)
ba = BackUpPowerSource[Quantity[5, "Megawatts"], Quantity[5, "Megawatts"]];
vDR = VariableDemandResponce[Quantity[2, "Megawatts"], Quantity[0.2, "Megawatts"], Quantity[50, ("Euros")/("Hours" "Megawatts")], Quantity[1, ("Euros")/("Hours" "Megawatts")], Quantity[1, ("Hours")/("Euros")], Quantity[1, ("Hours")/("Euros")]];](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/3ac2e3b416c83d14.png)
![sim = ElectricalSystemSimulate[<|
(*the first in list gets used first*)
"WindGenerator1" -> w,
"SolarGenerator1" -> s,
"NuclearGeneration1" -> n,
"HydroPlant" -> h,
"ConstantSink" -> c,
"VariableSink" -> v,
"BatteryBank" -> b
|>,(*startdate enddate stepsize*)
{startDate, endDate, step}
];](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/4ba09313dd804f38.png)
![allComponents = ListLinePlot[Total /@ sim["NetProduction"], PlotRange -> All, TargetUnits -> "MegaWatts"];
netProd = ListLinePlot[Total[Total /@ sim["NetProduction"]], PlotRange -> All, TargetUnits -> "MegaWatts"];
storages = ListLinePlot[sim["State"], PlotRange -> All, TargetUnits -> "Megawatts"*"Hours"];
generation = ListLinePlot[sim[["NetProduction", "Sources"]], PlotRange -> All, TargetUnits -> "MegaWatts"];
demand = ListLinePlot[sim[["NetProduction", "Sinks"]], PlotRange -> All, TargetUnits -> "MegaWatts"];
storage = ListLinePlot[sim[["NetProduction", "Storage"]], PlotRange -> All, TargetUnits -> "MegaWatts"];
Grid[{{allComponents, netProd}, {generation, demand}, {storages, storage}}]](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/7a758df0a39cd516.png)
![SimpleBattery[Quantity[150 , "Hours" "Megawatts"], Quantity[5, "Megawatts"], Quantity[5, "Megawatts"], 0.98, 0.00006]](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/3de4210fa38a9eab.png){kind=small}
![ConstantDemand[Quantity[0.25, "Megawatts"]]](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/47ab9ac416d86b36.png){kind=small}
![InterpretationBox[FrameBox[TagBox[TooltipBox[PaneBox[GridBox[List[List[GraphicsBox[List[Thickness[0.0025`], List[FaceForm[List[RGBColor[0.9607843137254902`, 0.5058823529411764`, 0.19607843137254902`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]], List[List[0, 2, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3], List[0, 1, 0], List[1, 3, 3]]], List[List[List[205.`, 22.863691329956055`], List[205.`, 212.31669425964355`], List[246.01799774169922`, 235.99870109558105`], List[369.0710144042969`, 307.0436840057373`], List[369.0710144042969`, 117.59068870544434`], List[205.`, 22.863691329956055`]], List[List[30.928985595703125`, 307.0436840057373`], List[153.98200225830078`, 235.99870109558105`], List[195.`, 212.31669425964355`], List[195.`, 22.863691329956055`], List[30.928985595703125`, 117.59068870544434`], List[30.928985595703125`, 307.0436840057373`]], List[List[200.`, 410.42970085144043`], List[364.0710144042969`, 315.7036876678467`], List[241.01799774169922`, 244.65868949890137`], List[200.`, 220.97669792175293`], List[158.98200225830078`, 244.65868949890137`], List[35.928985595703125`, 315.7036876678467`], List[200.`, 410.42970085144043`]], List[List[376.5710144042969`, 320.03370475769043`], List[202.5`, 420.53370475769043`], List[200.95300006866455`, 421.42667961120605`], List[199.04699993133545`, 421.42667961120605`], List[197.5`, 420.53370475769043`], List[23.428985595703125`, 320.03370475769043`], List[21.882003784179688`, 319.1406993865967`], List[20.928985595703125`, 317.4896984100342`], List[20.928985595703125`, 315.7036876678467`], List[20.928985595703125`, 114.70369529724121`], List[20.928985595703125`, 112.91769218444824`], List[21.882003784179688`, 111.26669120788574`], List[23.428985595703125`, 110.37369346618652`], List[197.5`, 9.87369155883789`], List[198.27300024032593`, 9.426692008972168`], List[199.13700008392334`, 9.203690528869629`], List[200.`, 9.203690528869629`], List[200.86299991607666`, 9.203690528869629`], List[201.72699999809265`, 9.426692008972168`], List[202.5`, 9.87369155883789`], List[376.5710144042969`, 110.37369346618652`], List[378.1179962158203`, 111.26669120788574`], List[379.0710144042969`, 112.91769218444824`], List[379.0710144042969`, 114.70369529724121`], List[379.0710144042969`, 315.7036876678467`], List[379.0710144042969`, 317.4896984100342`], List[378.1179962158203`, 319.1406993865967`], List[376.5710144042969`, 320.03370475769043`]]]]], List[FaceForm[List[RGBColor[0.5529411764705883`, 0.6745098039215687`, 0.8117647058823529`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[44.92900085449219`, 282.59088134765625`], List[181.00001525878906`, 204.0298843383789`], List[181.00001525878906`, 46.90887451171875`], List[44.92900085449219`, 125.46986389160156`], List[44.92900085449219`, 282.59088134765625`]]]]], List[FaceForm[List[RGBColor[0.6627450980392157`, 0.803921568627451`, 0.5686274509803921`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[355.0710144042969`, 282.59088134765625`], List[355.0710144042969`, 125.46986389160156`], List[219.`, 46.90887451171875`], List[219.`, 204.0298843383789`], List[355.0710144042969`, 282.59088134765625`]]]]], List[FaceForm[List[RGBColor[0.6901960784313725`, 0.5882352941176471`, 0.8117647058823529`], Opacity[1.`]]], FilledCurveBox[List[List[List[0, 2, 0], List[0, 1, 0], List[0, 1, 0], List[0, 1, 0]]], List[List[List[200.`, 394.0606994628906`], List[336.0710144042969`, 315.4997024536133`], List[200.`, 236.93968200683594`], List[63.928985595703125`, 315.4997024536133`], List[200.`, 394.0606994628906`]]]]]], List[Rule[BaselinePosition, Scaled[0.15`]], Rule[ImageSize, 10], Rule[ImageSize, 15]]], StyleBox[RowBox[List["NuclearGeneration", " "]], Rule[ShowAutoStyles, False], Rule[ShowStringCharacters, False], Rule[FontSize, Times[0.9`, Inherited]], Rule[FontColor, GrayLevel[0.1`]]]]], Rule[GridBoxSpacings, List[Rule["Columns", List[List[0.25`]]]]]], Rule[Alignment, List[Left, Baseline]], Rule[BaselinePosition, Baseline], Rule[FrameMargins, List[List[3, 0], List[0, 0]]], Rule[BaseStyle, List[Rule[LineSpacing, List[0, 0]], Rule[LineBreakWithin, False]]]], RowBox[List["PacletSymbol", "[", RowBox[List["\"NiiloKantoniemi/EnergySystem\"", ",", "\"NiiloKantoniemi`EnergySystem`NuclearGeneration\""]], "]"]], Rule[TooltipStyle, List[Rule[ShowAutoStyles, True], Rule[ShowStringCharacters, True]]]], Function[Annotation[Slot[1], Style[Defer[PacletSymbol["NiiloKantoniemi/EnergySystem", "NiiloKantoniemi`EnergySystem`NuclearGeneration"]], Rule[ShowStringCharacters, True]], "Tooltip"]]], Rule[Background, RGBColor[0.968`, 0.976`, 0.984`]], Rule[BaselinePosition, Baseline], Rule[DefaultBaseStyle, List[]], Rule[FrameMargins, List[List[0, 0], List[1, 1]]], Rule[FrameStyle, RGBColor[0.831`, 0.847`, 0.85`]], Rule[RoundingRadius, 4]], PacletSymbol["NiiloKantoniemi/EnergySystem", "NiiloKantoniemi`EnergySystem`NuclearGeneration"], Rule[Selectable, False], Rule[SelectWithContents, True], Rule[BoxID, "PacletSymbolBox"]][Quantity[5, "Megawatts"]]](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/1-0-0/50b078ac76ff1128.png){kind=small}
