Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
NiiloKantoniemi/
Energy
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
Basic Examples
| In[1]:= | ![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/3b3a102d915d5966.png) |
| In[2]:= | ![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/7b9efb79cd95818e.png) |
| In[3]:= | ![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/516f5817f6fba145.png) |
| In[4]:= | ![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/5c91d769fb57476a.png) |
| In[5]:= | ![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/6b5d6e1811988b25.png) |
| Out[11]= |  |
| In[12]:= | ![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/7d6117b5a40f3b87.png){kind=small} |
| Out[12]= |  |
| In[13]:= | ![ConstantDemand[Quantity[0.25, "Megawatts"]]](https://www.wolframcloud.com/obj/resourcesystem/images/702/702ddaf1-4212-4a5d-b39c-8ba25b1bdaca/58460fc3f937f4ba.png){kind=small} |
| Out[13]= |  |
| In[14]:= | ![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/62efe9d370aaad0b.png){kind=small} |
| Out[14]= |  |
Publisher
Compatibility
Wolfram Language Version 14
Version History
- 1.0.1 – 09 July 2024
- 1.0.0 – 09 July 2024