Kurs:Mathematische Modellbildung/Themen/2024-25 Wintersemester/Thema 1
Energiegewinnung und Flächennutzung: Mathematische Modellierung eines Photovoltaikmodul-Systems für den Universitätsparkplatz der RPTU Landau
[Bearbeiten]In diesem Projekt beschäftigen wir uns mit der Frage, welche Vorteile es für die Umwelt haben kann, vorhandene Parkplatzflächen mit Solarmodulen zu überdachen. Dafür betrachten wir beispielhaft den Universitätsparkplatz der RPTU am Campus Landau. Wir gehen dabei vorallem der Frage auf den Grund, wie viel nachhaltiger Strom durch solche Solarmodule erzeugt werden kann, der andernfalls durch umweltschädlichere Gewinnungsmethoden wie Atomkraft erzeugt werden müsste. Das Projekt soll als Beispiel für andere Parkplatzflächen dienen, das einen möglichst genauen Richtwert für die Menge an Energie liefert, die auf einer bestimmten Fläche innerhalb eines Jahres gewonnen werden kann.
Gruppenteilnehmer
[Bearbeiten]- Jannis Walker
- Alina Döringer
- Lukas Daubner
Zielsetzung des Modellierungsthemas
[Bearbeiten](Text)
Zuordnung zu den Nachhaltigkeitszielen
[Bearbeiten]- SDG 7: Bezahlbare und saubere Energien
Die Installation von Solarpanels trägt zur Produktion von erneuerbarer Energie bei und fördert den Zugang zu sauberer und nachhaltiger Energie
- SDG 13: Maßnahmen zum Klimaschutz
Die Nutzung von Solarenergie reduziert die Abhängigkeit von fossilen Brennstoffen und hilft den CO2-Ausstoß zu verringern, was einen Beitrag zum Klimaschutz leistet
- SDG 15:Leben an Land
Durch die Nutzung bestehender Flächen zur Energiegewinnung wird zusätzlich Flächenverbrauch vermieden, was zur Schonung natürlicher Lebensräume beitragt
Datenerhebung
[Bearbeiten](Text)
Modellierungszyklus 1(Level: Sekundarstufe 2)
[Bearbeiten]Wir beginnen unser Projekt mit der Problemstellung, wie die Sonne an einem beliebigen Jahrestag t ∈ {1,2,...,364,365} über den Landauer Himmel verläuft. Dafür betrachten wir zunächst die sogenannte Sonnenhöhe. Diese gibt an, in welchem Winkel (d.h. aus welcher "Höhe") die Sonnenstrahlen auf die Erdoberfläche treffen. Die Sonnenhöhe erreicht ihren höchsten Wert jedes Jahr am 21. Juni (t = 172), dem Tag der Sommersonnenwende, wenn die Sonne am nördlichen Wendekreis im Zenit steht. Ihren niedrigsten Wert erreicht sie stets am 21. Dezember (t = 355), dem Tag der Wintersonnenwende, wenn die Sonne am südlichen Wendekreis im Zenit steht. Zwischen diesen beiden Extremwerten wächst bzw. fällt die Sonnenhöhe und erreicht sowohl um den 21. März (t ≈ 80) als auch um den 23. September (t ≈ 266) herum ihren Mittelwert. Zu diesen beiden Zeitpunkten steht die Sonne auf dem Äquator im Zenit. Grund für die Verschiebung des Zenitstandes der Sonne und damit auch die Veränderung der Sonnenhöhe in Landau sind die Schiefstellung der Erdachse um 23,5° und die Erdrevolution.
Die Größe des Azimuth-Winkels in Abhängigkeit vom Jahrestag t lässt sich als Sinusfunktion
f: {1,...,365} ≕ D → ℝ; f(t) = a ⋅ sin(b⋅t + c) + d in °
mit passenden a,b,c,d ∈ ℝ beschreiben. Um die Parameter a,b,c,d der Sinusfunktion bestimmen zu können, legen wir eine maßstabsgetreue GeoGebra-Datei an, die uns einen Wert für den Sonnenhöhenwinkel zu unterschiedlichen Zeitpunkten ausgibt. Wir wählen dafür genau die vier bzw. drei (21. März und 23. September ergeben denselben Sonnenhöhenwinkel) Zeitpunkte aus, zu denen die Sonne an einem der Wendekreise bzw. am Äquator im Zenit steht:
https://www.geogebra.org/classic/vq9kuu7u
Aus der GeoGebra-Datei wird deutlich:
• Steht die Sonne am nördlichen Wendekreis im Zenit, beträgt der Sonnenhöhenwinkel in Landau 64,3°
• Steht die Sonne am Äquator im Zenit, beträgt der Sonnenhöhenwinkel in Landau 40,8°
• Steht die Sonne am südlichen Wendekreis im Zenit, beträgt der Sonnenhöhenwinkel in Landau 17,3°
Aus diesen Daten lassen sich die Parameter der Sinusfunktion wie folgt bestimmen:
• Die Amplitude a muss der Hälfte des Unterschiedes zwischen dem zuerst erreichten Extremwert und dem zweiten Extremwert entsprechen: a = max t ∈ D f - min t ∈ D f) / 2 = (64,3 - 17,3) / 2 = 47/2 = 23,5
• Die y-Verschiebung d muss um eine Amplitude kleiner sein als der größte Funktionswert bzw. eine Amplitude größer als der kleinste Funktionswert: d = max t ∈ D f - a = min t ∈ D f + a = 64,3 - 23,5 = 40,8
• Die Streckung in x-Richtung b lässt sich aus der Periode p = 365 der Sinusfunktion wie folgt bestimmen: p = 365 = (2⋅π) / b ⇔ b = (2⋅π) / 365
• Die Verschiebung in x-Richtung c muss nun so angepasst werden, dass die Funktionswerte die reale Situation widerspiegeln. Dafür betrachten wir den 21. Juni (t = 172). An diesem Tag soll die Funktion ihr Maximum erreichen: f(172) = 23,5 ⋅ sin((2⋅π) / 365 ⋅ 172 + c) + 40,8 = 64,3 ⇔ 23,5 ⋅ sin((2⋅π) / 365 ⋅ 172 + c) = 64,3 - 40,8 = 23,5 ⇔ sin((2⋅π) / 365 ⋅ 172 + c) = 23,5 / 23,5 = 1 ⇔ (2⋅π) / 365 ⋅ 172 + c = π/2 ⇔ c = π/2 - (344⋅π) / 365 = (365⋅π) / 730 - (688⋅π) / 730 = (-323⋅π) / 730
Somit lautet die Funktionsgleichung für f: f(t) = 23,5 ⋅ sin((2⋅π)/365 ⋅ 172 - (323⋅π)/730) + 40,8
Da die meisten Computerprogramme standardmäßig mit Werten im Bogenmaß rechnen und nicht mit Werten im Gradmaß, entscheiden wir uns, die erhaltenen Funktionswerte ins Bogenmaß umzurechnen, und erhalten somit eine neue Funktionsgleichung für f:
f(t) = (23,5⋅π)/180 ⋅ sin((2⋅π)/365 ⋅ 172 - (323⋅π)/730) + (40,8⋅π)/180
Im nächsten Schritt beschäftigen wir uns mit der Veränderung des Sonnenhöhenwinkels innerhalb eines Tages. Die "Höhe", aus der die Sonne auf den Uni-Parkplatz in Landau trifft, verändert sich natürlich nicht nur im Jahresverlauf, sondern auch im Laufe des Tages. Um den Verlauf dieses Winkels innerhalb eines Tages darstellen zu können, fertigen wir einige GeoGebra-Dateien in ℝ3 an. Unser Ziel dabei ist es, einen Zusammenhang zwischen Zeitpunkt x und Größe des Sonnenhöhenwinkels y für ein beliebiges t ∈ D zu finden. Wir betrachten für unsere Überlegungen dieselben vier Zeitpunkte, die wir im ersten Schritt auch betrachtet haben. Darüberhinaus legen wir für alle t ∈ D den Zeitpunkt, an dem die Sone aufgeht, als Zeitpunkt x=0 fest, sowie den Zeitpunkt, an dem die Sonne untergeht, als Zeitpunkt x=1.
• 21. Dezember: https://www.geogebra.org/classic/hwjw66jr
• 21. Juni: https://www.geogebra.org/classic/ar9z6qgh
• 21. März bzw. 23. September:
Wie man in der Datei sehen kann, erreicht der Sonnenhöhenwinkel für ein festes t ∈ D sein Maximum stets um die Mittagszeit, zum Zeitpunkt x = 0,5. Dabei ist zu beachten, dass dieser Zeitpunkt nicht zwangsläufig mit 12:00 Uhr mittags übereinstimmt, sondern es handelt sich um den Zeitpunkt, an dem Landau auf seiner "Kreisbahn um die Erdachse" genau die Hälfte des Weges zwischen Sonnenaufgang und Sonnenuntergang zurückgelegt hat. Gleichzeitig handelt es sich dabei um genau den Zeitpunkt, zu dem Landau (an diesem Tag t) der Sonne am nächsten ist. Wenn man die dreidimensionale GeoGebra-Datei zu diesem Zeitpunkt richtig dreht (d.h. so, dass die y-Achse zu einem Punkt wird), dann erhält man dieselbe Perspektive, die wir in unserer ersten GeoGebra-Datei eingenommen haben. Nun wird deutlich, dass dieses Maximum in x = 0,5 für alle t ∈ D mit der zuvor berechneten Sonnenhöhe f(t) übereinstimmt. Darüber hinaus beträgt der Sonnenhöhenwinkel bei Sonnenauf- und Sonnenuntergang logischerweise 0°. Diese Idee verallgemeinern wir, um eine Formel zu finden, die den Sonnenhöhenwinkel in Abhängigkeit vom Zeitpunkt x für einen beliebigen Jahrestag t beschreibt. Dafür haben wir zwei Ideen:
1) Man könnte den Verlauf des Sonnenhöhenwinkels als eine nach unten geöffnete Parabel auffassen, die ihre Nullstellen in x1 = 0 und x2 = 1 hat und ihren Scheitelpunkt S (S hat automatisch die x-Koordinate 0,5) in (0,5|f(t)). Das bedeutet, dass die Parabel die Funktionsgleichung p(x) = a ⋅ (x-1) ⋅ (x-0) mit passendem a ∈ ℝ lauten muss. Damit die y-Koordinate des Scheitelpunktes mit f(t) übereinstimmt, stellen wir folgende Gleichung auf:
p(0,5) = f(t) ⇔ a ⋅ (0,5-1) ⋅ (0,5-0) = f(t) ⇔ a ⋅ (-1/4) = f(t) ⇔ a = (-4) ⋅ f(t)
Dementsprechend lautet die Funktionsgleichung für die Parabel (nun in Abhängigkeit sowohl von t als auch von x):
p: {1,...,365} x [0;1] → ℝ; p(t,x) = (-4) ⋅ f(t) ⋅ (x-1) ⋅ (x-0)
2) Man könnte den Verlauf des Sonnenhöhenwinkels allerdings auch als eine Sinusfunktion auffassen, die Nullstellen in x1 = 0 und x2 = 1 und ihr Maximum in (0,5|f(t)) hat. Die Funktionsgleichung muss dann lauten:
g: {1,...,365} x [0;1] → ℝ; g(t,x) = f(t) ⋅ sin(π⋅x)
Um die beiden Ideen miteinander vergleichen zu können, fertigen wir eine GeoGebra-Datei an, die die beiden Funktionen im Bereich [0;1] darstellt:
https://www.geogebra.org/classic/egcyz3up
weiter gehts .....
Annäherung an den Flächeninhalt des (nicht rechteckigen) Universitätsparkplatzes durch Einheitsquadrate, Bildung von einer Art Ober- und Untersumme bezüglich der Einheitsquadrate, Verfeinerung der Einheit, sodass Flächeninhalt genauer geschätzt werden kann, Beobachtung des Konvergenzverhaltens der Anzahl an Einheiten bei immer kleineren Quadraten Beschäftigung mit der Frage, wie viele Solarmodule auf den Parkplatz passen Grobe Annäherung an die Energiemenge durch Mittelwerte, wie viel Strom ein Solarmodul erzeugen kann und Anzahl der Tage im Jahr
Voraussetzungen
[Bearbeiten](Text)
Umsetzung mit mathematischen Werkzeugen
[Bearbeiten](Text)
Ergebnisse
[Bearbeiten](...)
Diskussion
[Bearbeiten](Verbesserungsvorschläge)
Modellierungszyklus 2 (Level: Uni-Niveau)
[Bearbeiten]Betrachtung der durchschnittlichen Sonnenstunden pro Tag je nach Jahreszeit, darauf basierend durch Regression eine Sinusfunktion ermitteln, die die Sonnenstunden in Landau pro Tag je nach Jahreszeit am besten annähert, dementsprechend feinere Ergebnisse für mögliche Energiemenge der Solarmodule mithilfe von Integration der Kurve
Voraussetzungen
[Bearbeiten](Text)
- Listeneintrag
- Listeneintrag
- Listeneintrag
Umsetzung mit mathematischen Werkzeugen
[Bearbeiten](Text)
Ergebnisse
[Bearbeiten](...)
Diskussion
[Bearbeiten](Verbesserungsvorschläge)
Modellierungszyklus 3 (Level: Uni-Niveau)
[Bearbeiten]zusätzliche Betrachtung des veränderten Winkels der Sonneneinstrahlung durch Erdrevolution, annäherungsweise Berechnung dieses Winkels und erneute Regression zur Bestimmung einer Funktion, die den Winkel je nach Jahreszeit am besten annähert, daraufhin Aufstellen einer Funktion, die die Parameter Winkel, Sonnenstunden, Jahreszeit und durchschnittlicher Energiegewinn miteinander vereint und einen finalen Wert für die Gesamtenergie liefert
Voraussetzungen
[Bearbeiten](Text)
- Listeneintrag
- Listeneintrag
- Listeneintrag
Umsetzung mit mathematischen Werkzeugen
[Bearbeiten](Text)
Ergebnisse
[Bearbeiten](...)
Diskussion
[Bearbeiten](Verbesserungsvorschläge)