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
|
\begin{tikzpicture}
\pgfplotstableread[col sep=comma]{img/data/cylinder2d_optimized_refinement_n17_re100_drag_lift_deltap.csv}\refined
\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n40_re100_drag_lift_deltap.csv}\uniformHigh
\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n80_re100_drag_lift_deltap.csv}\uniformVeryHigh
\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n160_re100_drag_lift_deltap.csv}\uniformVeryVeryHigh
\begin{axis}[
scale only axis,
height=5cm,
width=0.9*\textwidth,
legend cell align=left,
legend pos=south east,
grid=both,
domain=0:16,
xmin=0, xmax=16,
ylabel={Widerstandskoeffizient},
ylabel absolute, every axis y label/.append style={yshift=0.4cm}
]
\addplot[thick,color=black]{3.23};
\addlegendentry {\(\widehat{c_w} := 3.23\)};
\addplot[
color=blue!50!white,
thin
] table [
x expr=\thisrow{time}, y=drag
] {\uniformVeryVeryHigh};
\addplot[
color=red!50!white,
thin
] table [
x expr=\thisrow{time}, y=drag
] {\uniformVeryHigh};
\addplot[
color=green!70!black,
thin
] table [
x expr=8*\thisrow{time}, y=drag
] {\refined};
\end{axis}
\begin{axis}[
scale only axis,
yshift=-6cm,
height=5cm,
width=0.9*\textwidth,
legend cell align=left,
legend pos=south west,
grid=both,
domain=0:16,
xmin=0, xmax=16,
ylabel={Auftriebskoeffizient},
ylabel absolute, every axis y label/.append style={yshift=0.4cm}
]
\addplot[thick,color=black]{1.0};
\addlegendentry {\(\widehat{c_a} := 1\)};
\addplot[
color=blue!50!white,
thin
] table [
x expr=\thisrow{time}, y=lift
] {\uniformVeryVeryHigh};
\addplot[
color=red!50!white,
thin
] table [
x expr=\thisrow{time}, y=lift
] {\uniformVeryHigh};
\addplot[
color=green!70!black,
thin
] table [
x expr=8*\thisrow{time}, y=lift
] {\refined};
\end{axis}
\begin{axis}[
scale only axis,
yshift=-12cm,
height=5cm,
width=0.9*\textwidth,
legend cell align=left,
legend pos=south east,
grid=both,
domain=0:16,
xmin=0, xmax=16,
xlabel={Simulierte physikalische Zeit},
ylabel={Druckdifferenz},
x unit=s,
y unit=N/m^2,
ylabel absolute, every axis y label/.append style={yshift=0.4cm}
]
\addplot[
color=blue!50!white,
thin
] table [
x expr=\thisrow{time}, y=deltap
] {\uniformVeryVeryHigh};
\addlegendentry {Uniformes \(N=40\) Gitter mit 145314 Knoten};
\addplot[
color=red!50!white,
thin
] table [
x expr=\thisrow{time}, y=deltap
] {\uniformVeryHigh};
\addlegendentry {Uniformes \(N=80\) Gitter mit 576758 Knoten};
\addplot[
color=green!70!black,
thin
] table [
x expr=8*\thisrow{time}, y=deltap
] {\refined};
\addlegendentry {Problembezogen verfeinertes \(N=17\) Gitter};
\addplot[thick,color=black]{2.48};
\addlegendentry {\(\Delta P := 2.48\)};
\end{axis}
\end{tikzpicture}
|
