1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
%
% waermeleitung_beispiel.tex -- Beispiel Wärmeleitung in homogenem Stab.
%%%%%%%%%%%%%%%%%%%%%%%%%%% Erster Entwurf %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
\subsection{Wärmeleitung in einem Homogenen Stab}
In diesem Abschnitt betrachten wir das Problem der Wärmeleitung in einem
homogenen Stab und wie das Sturm-Liouville-Problem bei der Beschreibung dieses
physikalischen Phänomenes auftritt.
Zunächst wird ein eindimensionaler homogener Stab der Länge $l$ und
Wärmeleitkoeffizient $\kappa$. Somit ergibt sich für das Wärmeleitungsproblem
die partielle Differentialgleichung
\[
\frac{\partial u}{\partial t} =
\kappa \frac{\partial^{2}u}{{\partial x}^{2}}
\]
wobei der Stab in diesem Fall auf der X-Achse im Intervall $[0,l]$ liegt.
Da diese Differentialgleichung das Problem allgemein für einen homogenen
Stab beschreibt, werden zusätzliche Bedingungen benötigt, um beispielsweise
die Lösung für einen Stab zu finden, bei dem die Enden auf konstanter
Tempreatur gehalten werden.
%%%%%%%%%%%%% Randbedingungen für Stab mit konstanten Endtemperaturen %%%%%%%%%
\subsubsection{Stab mit Enden auf konstanter Temperatur}
Die Enden des Stabes auf konstanter Temperatur zu halten bedeutet, dass die
Lösungsfunktion $u(t,x)$ bei $x = 0$ und $x = l$ nur die vorgegebene
Temperatur zurückgeben darf. Es folgen nun
\[
u(t,0)
=
u(t,l)
=
0
\]
als Randbedingungen.
%%%%%%%%%%%%% Randbedingungen für Stab mit isolierten Enden %%%%%%%%%%%%%%%%%%%
\subsubsection{Stab mit isolierten Enden}
Bei isolierten Enden des Stabes können belibige Temperaturen für $x = 0$ und
$x = l$ auftreten. In diesem Fall nicht erlaubt ist es, dass Wärme vom Stab
an die Umgebung oder von der Umgebung an den Stab abgegeben wird.
Aus der Physik ist bekannt, dass Wärme immer von der höheren zur tieferen
Temperatur fliesst. Um Wärmefluss zu unterdrücken, muss also dafür gesorgt
werden, dass am Rand des Stabes keine Temperaturdifferenz existiert oder
dass die partiellen Ableitungen von $u(t,x)$ nach $x$ bei $x = 0$ und $x = l$
verschwinden. Somit folgen
\[
\frac{\partial}{\partial x} u(t, 0)
=
\frac{\partial}{\partial x} u(t, l)
=
0
\]
als Randbedingungen.
%%%%%%%%%%% Lösung der Differenzialgleichung %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsubsection{Lösung der Differenzialgleichung}
% TODO: Referenz Separationsmethode
% TODO: Formeln sauber in Text einbinden.
Da die Lösungsfunktion von zwei Variablen abhängig ist, wird als Lösungsansatz
die Separationsmethode verwendet.
\[
u(t,x)
=
T(t)X(x)
\]
Dieser Ausdruck wird in die ursprüngliche Differenzialgleichung eingesetzt:
\[
T^{\prime}(t)X(x)
=
\kappa T(t)X^{\prime \prime}(x)
\]
Nun können alle von $t$ abhängigen Ausdrücke auf die linke Seite, sowie alle
von $x$ abhängigen Ausdrücke auf die rechte Seite gebracht werden und mittels
der neuen Variablen $\mu$ gekoppelt werden:
\[
\frac{T^{\prime}(t)}{\kappa T(t)}
=
\frac{X^{\prime \prime}(x)}{X(x)}
=
\mu
\]
Durch die Einführung von $\mu$ kann das Problem nun in zwei separate
Differenzialgleichungen aufgeteilt werden:
\[
T^{\prime}(t) - \kappa \mu T(t)
=
0
\]
\[
X^{\prime \prime}(x) - \mu X(x)
=
0
\]
% TODO: Rechenweg
TODO: Rechenweg... Enden auf konstanter Temperatur:
\[
u(t,x)
=
\sum_{n=1}^{\infty}a_{n}e^{-\frac{n^{2}\pi^{2}\kappa}{l^{2}}t}
\sin\left(\frac{n\pi}{l}x\right)
\]
\[
a_{n}
=
\frac{2}{l}\int_{0}^{l}u(0,x)sin\left(\frac{n\pi}{l}x\right) dx
\]
TODO: Rechenweg... Enden isoliert:
\[
u(t,x)
=
a_{0} + \sum_{n=1}^{\infty}a_{n}e^{-\frac{n^{2}\pi^{2}\kappa}{l^{2}}t}
\cos\left(\frac{n\pi}{l}x\right)
\]
\[
a_{0}
=
\frac{1}{l}\int_{0}^{l}u(0,x) dx
\]
\[
a_{n}
=
\frac{2}{l}\int_{0}^{l}u(0,x)sin\left(\frac{n\pi}{l}x\right) dx
\]
|
