Add Poiseuille data, plots

10 files changed, 261 insertions, 5 deletions

diff --git a/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_bb_walls.csv b/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_bb_walls.csv
new file mode 100644
index 0000000..48eb46f
--- /dev/null
+++ b/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_bb_walls.csv
@@ -0,0 +1,20 @@
+geometry;ux;uy;p;x;y
+4;0,10501491;0;1;4;0,050000001
+4;0,3021979;0;1;4;0,1
+4;0,47633973;0;1;4;0,15000001
+4;0,6270411;0;1;4;0,2
+4;0,75486946;0;1;4;0,25
+4;0,8595866;0;1;4;0,30000001
+4;0,94110405;0;1;4;0,34999999
+4;0,99936461;0;1;4;0,40000001
+4;1,03435;0;1;4;0,44999999
+4;1,0460171;0;1;4;0,5
+4;1,03435;0;1;4;0,55000001
+4;0,99936461;0;1;4;0,60000002
+4;0,94110405;0;1;4;0,64999998
+4;0,8595866;0;1;4;0,69999999
+4;0,75486946;0;1;4;0,75
+4;0,6270411;0;1;4;0,80000001
+4;0,47633973;0;1;4;0,85000002
+4;0,3021979;0;1;4;0,89999998
+4;0,10501491;0;1;4;0,94999999
diff --git a/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_velocity_walls.csv b/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_velocity_walls.csv
new file mode 100644
index 0000000..85b0053
--- /dev/null
+++ b/img/data/poiseuille2d_refined_order2ipol_velocity_outflow_velocity_walls.csv
@@ -0,0 +1,20 @@
+geometry,ux,uy,p,x,y
+4,0.18898219,0,1,4,0.050000001
+4,0.36055607,0,1,4,0.1
+4,0.51207048,0,1,4,0.15000001
+4,0.64336717,0,1,4,0.2
+4,0.75464326,0,1,4,0.25
+4,0.84575975,0,1,4,0.30000001
+4,0.91669858,0,1,4,0.34999999
+4,0.96739316,0,1,4,0.40000001
+4,0.99782735,0,1,4,0.44999999
+4,1.0079721,0,1,4,0.5
+4,0.99782735,0,1,4,0.55000001
+4,0.96739316,0,1,4,0.60000002
+4,0.91669858,0,1,4,0.64999998
+4,0.84575975,0,1,4,0.69999999
+4,0.75464326,0,1,4,0.75
+4,0.64336717,0,1,4,0.80000001
+4,0.51207048,0,1,4,0.85000002
+4,0.36055607,0,1,4,0.89999998
+4,0.18898219,0,1,4,0.94999999
diff --git a/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_bb_walls.csv b/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_bb_walls.csv
new file mode 100644
index 0000000..37db705
--- /dev/null
+++ b/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_bb_walls.csv
@@ -0,0 +1,20 @@
+geometry;ux;uy;p;x;y
+4;0,10513373;0;1;4;0,050000001
+4;0,30253884;0;1;4;0,1
+4;0,47687662;0;1;4;0,15000001
+4;0,62774706;0;1;4;0,2
+4;0,75571918;0;1;4;0,25
+4;0,86055434;0;1;4;0,30000001
+4;0,94216371;0;1;4;0,34999999
+4;1,00049;0;1;4;0,40000001
+4;1,035515;0;1;4;0,44999999
+4;1,0471953;0;1;4;0,5
+4;1,035515;0;1;4;0,55000001
+4;1,00049;0;1;4;0,60000002
+4;0,94216371;0;1;4;0,64999998
+4;0,86055434;0;1;4;0,69999999
+4;0,75571918;0;1;4;0,75
+4;0,62774706;0;1;4;0,80000001
+4;0,47687662;0;1;4;0,85000002
+4;0,30253884;0;1;4;0,89999998
+4;0,10513373;0;1;4;0,94999999
diff --git a/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_velocity_walls.csv b/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_velocity_walls.csv
new file mode 100644
index 0000000..b0125c3
--- /dev/null
+++ b/img/data/poiseuille2d_refined_order4ipol_velocity_outflow_velocity_walls.csv
@@ -0,0 +1,20 @@
+geometry,ux,uy,p,x,y
+4,0.18905593,0,1,4,0.050000001
+4,0.36069685,0,1,4,0.1
+4,0.51227051,0,1,4,0.15000001
+4,0.6436184,0,1,4,0.2
+4,0.75493807,0,1,4,0.25
+4,0.8460902,0,1,4,0.30000001
+4,0.91705674,0,1,4,0.34999999
+4,0.96777117,0,1,4,0.40000001
+4,0.99821723,0,1,4,0.44999999
+4,1.008366,0,1,4,0.5
+4,0.99821723,0,1,4,0.55000001
+4,0.96777117,0,1,4,0.60000002
+4,0.91705674,0,1,4,0.64999998
+4,0.8460902,0,1,4,0.69999999
+4,0.75493807,0,1,4,0.75
+4,0.6436184,0,1,4,0.80000001
+4,0.51227051,0,1,4,0.85000002
+4,0.36069685,0,1,4,0.89999998
+4,0.18905593,0,1,4,0.94999999
diff --git a/img/data/poiseuille2d_unrefined_coarse_velocity_outflow_velocity_walls.csv b/img/data/poiseuille2d_unrefined_coarse_velocity_outflow_velocity_walls.csv
new file mode 100644
index 0000000..417ccd6
--- /dev/null
+++ b/img/data/poiseuille2d_unrefined_coarse_velocity_outflow_velocity_walls.csv
@@ -0,0 +1,10 @@
+geometry,ux,uy,p,x,y
+4,0.35229355,0,1,4,0.1
+4,0.63773042,0,1,4,0.2
+4,0.84351379,0,1,4,0.30000001
+4,0.96691775,0,1,4,0.40000001
+4,1.0083944,0,1,4,0.5
+4,0.96691775,0,1,4,0.60000002
+4,0.84351379,0,1,4,0.69999999
+4,0.63773042,0,1,4,0.80000001
+4,0.35229355,0,1,4,0.89999998
diff --git a/img/data/poiseuille2d_unrefined_velocity_outflow_bb_walls.csv b/img/data/poiseuille2d_unrefined_velocity_outflow_bb_walls.csv
new file mode 100644
index 0000000..bcf70bd
--- /dev/null
+++ b/img/data/poiseuille2d_unrefined_velocity_outflow_bb_walls.csv
@@ -0,0 +1,20 @@
+geometry;ux;uy;p;x;y
+4;0,10556783;0;1;4;0,94999999
+4;0,30378494;0;1;4;0,89999998
+4;0,47883838;0;1;4;0,85000002
+4;0,63032639;0;1;4;0,80000001
+4;0,75882447;0;1;4;0,75
+4;0,86409104;0;1;4;0,69999999
+4;0,94603628;0;1;4;0,64999998
+4;1,0046031;0;1;4;0,60000002
+4;1,0397729;0;1;4;0,55000001
+4;1,0515014;0;1;4;0,5
+4;1,0397729;0;1;4;0,44999999
+4;1,0046031;0;1;4;0,40000001
+4;0,94603628;0;1;4;0,34999999
+4;0,86409104;0;1;4;0,30000001
+4;0,75882447;0;1;4;0,25
+4;0,63032639;0;1;4;0,2
+4;0,47883838;0;1;4;0,15000001
+4;0,30378494;0;1;4;0,1
+4;0,10556783;0;1;4;0,050000001
diff --git a/img/data/poiseuille2d_unrefined_velocity_outflow_velocity_walls.csv b/img/data/poiseuille2d_unrefined_velocity_outflow_velocity_walls.csv
new file mode 100644
index 0000000..ee4741c
--- /dev/null
+++ b/img/data/poiseuille2d_unrefined_velocity_outflow_velocity_walls.csv
@@ -0,0 +1,20 @@
+geometry,ux,uy,p,x,y
+4,0.1880112,0,1,4,0.050000001
+4,0.35870126,0,1,4,0.1
+4,0.50943297,0,1,4,0.15000001
+4,0.64005154,0,1,4,0.2
+4,0.75075144,0,1,4,0.25
+4,0.84139568,0,1,4,0.30000001
+4,0.91196662,0,1,4,0.34999999
+4,0.96239835,0,1,4,0.40000001
+4,0.99267453,0,1,4,0.44999999
+4,1.0027666,0,1,4,0.5
+4,0.99267453,0,1,4,0.55000001
+4,0.96239835,0,1,4,0.60000002
+4,0.91196662,0,1,4,0.64999998
+4,0.84139568,0,1,4,0.69999999
+4,0.75075144,0,1,4,0.75
+4,0.64005154,0,1,4,0.80000001
+4,0.50943297,0,1,4,0.85000002
+4,0.35870126,0,1,4,0.89999998
+4,0.1880112,0,1,4,0.94999999
diff --git a/img/poiseuille2d_grid.tikz b/img/poiseuille2d_grid.tikz
index 3df52a8..8ca285f 100644
--- a/img/poiseuille2d_grid.tikz
+++ b/img/poiseuille2d_grid.tikz
@@ -7,10 +7,10 @@
 	xmin=0, xmax=4,
 	ymin=0, ymax=1,
 	scatter/classes={
-		1={color=blue},
-		2={color=black},
-		3={color=green},
-		4={color=red}
+		1={color=blue!50!white},
+		2={color=gray!70!black},
+		3={color=green!70!black},
+		4={color=red!70!black}
 	},
 	legend style={
 		at={(0.5,-0.2)},
diff --git a/img/poiseuille2d_velocity_grid.tikz b/img/poiseuille2d_velocity_grid.tikz
index 4326784..61afaeb 100644
--- a/img/poiseuille2d_velocity_grid.tikz
+++ b/img/poiseuille2d_velocity_grid.tikz
@@ -7,14 +7,25 @@
 	xmin=0, xmax=4,
 	ymin=0, ymax=1,
 	point meta=explicit,
-	colormap/bluered,
 	colorbar,
+	colormap/RdYlBu-4,
+	colormap={reverse RdYlBu-4}{
+		indices of colormap={
+			\pgfplotscolormaplastindexof{RdYlBu-4},...,0 of RdYlBu-4
+		}
+	},
 	colorbar horizontal,
 	colorbar style={
 		/pgf/number format/precision=2,
 		at={(0.5,-0.2)},
 		anchor=north,
 		width=0.5*\pgfkeysvalueof{/pgfplots/parent axis width},
+		ylabel=\(m/s\),
+		y label style={
+			at={(axis description cs:1.3,0.5)},
+			rotate=-90,
+			anchor=east
+		}
 	},
 	scatter/use mapped color={color=mapped color}
 ]
diff --git a/img/poiseuille2d_velocity_outflow.tikz b/img/poiseuille2d_velocity_outflow.tikz
new file mode 100644
index 0000000..9b39e1e
--- /dev/null
+++ b/img/poiseuille2d_velocity_outflow.tikz
@@ -0,0 +1,115 @@
+\begin{tikzpicture}
+\pgfplotstableread[col sep=comma]{img/data/poiseuille2d_refined_order4ipol_velocity_outflow_velocity_walls.csv}\refinedCubicIpol
+\pgfplotstableread[col sep=comma]{img/data/poiseuille2d_refined_order2ipol_velocity_outflow_velocity_walls.csv}\refinedLinearIpol
+\pgfplotstableread[col sep=comma]{img/data/poiseuille2d_unrefined_velocity_outflow_velocity_walls.csv}\unrefined
+\pgfplotstableread[col sep=comma]{img/data/poiseuille2d_unrefined_coarse_velocity_outflow_velocity_walls.csv}\coarse
+
+\begin{axis}[
+	scale only axis,
+	height=8cm,
+	width=0.4*\textwidth,
+	mark size=4,
+	legend cell align=left,
+	legend style={at={(0.9,-0.1)},anchor=north},
+	grid=both,
+	ylabel=\(x\)-Geschwindigkeit
+	xtick={0,0.25,0.5,0.75,1},
+	title=Geschwindigkeitsprofil,
+	every axis plot/.append style={thick}
+]
+
+\addplot[
+	only marks,
+	mark=x,
+	color=green!70!black
+] table [
+	x=y, y=ux
+] {\refinedCubicIpol};
+\addlegendentry{Halbseitig verfeinertes Gitter mit kubischer Interpolation};
+
+\addplot[
+	only marks,
+	mark=+,
+	color=red!70!black
+] table [
+	x=y, y=ux
+] {\refinedLinearIpol};
+\addlegendentry{Halbseitig verfeinertes Gitter mit linearer Interpolation};
+
+\addplot[
+	only marks,
+	mark=asterisk,
+	color=blue!50!white
+] table[
+	x=y, y=ux
+] {\unrefined};
+\addlegendentry{Uniform fein aufgelöstes Gitter};
+
+\addplot[
+	only marks,
+	mark=star,
+	color=gray
+] table[
+	x=y, y=ux
+] {\coarse};
+\addlegendentry{Uniform grob aufgelöstes Gitter};
+
+\addplot [domain=0.05:0.95, samples=100]{-4*x*(x-1)};
+\addlegendentry{Analytische Lösung}
+\end{axis}
+
+\begin{axis}[
+	scale only axis,
+	height=8cm,
+	width=0.4*\textwidth,
+	mark size=4,
+	grid=both,
+	xtick={0,0.25,0.5,0.75,1},
+	xshift=7cm,
+	scaled y ticks=false,
+	ylabel=Fehler,
+	yticklabel pos=right,
+	y tick label style={/pgf/number format/sci},
+	title=Fehler zur analytischen Lösung,
+	every axis plot/.append style={thick}
+]
+
+\addplot[
+	only marks,
+	mark=x,
+	color=green!70!black
+] table [
+	x=y,
+	y expr=\thisrow{ux}+4*\thisrow{y}*(\thisrow{y}-1)
+] {\refinedCubicIpol};
+
+\addplot[
+	only marks,
+	mark=+,
+	color=red!70!black
+] table [
+	x=y,
+	y expr=\thisrow{ux}+4*\thisrow{y}*(\thisrow{y}-1)
+] {\refinedLinearIpol};
+
+\addplot[
+	only marks,
+	mark=asterisk,
+	color=blue!50!white
+] table[
+	x=y,
+	y expr=\thisrow{ux}+4*\thisrow{y}*(\thisrow{y}-1)
+] {\unrefined};
+
+\addplot[
+	only marks,
+	mark=star,
+	color=gray
+] table[
+	x=y,
+	y expr=\thisrow{ux}+4*\thisrow{y}*(\thisrow{y}-1)
+] {\coarse};
+
+\addplot [domain=0:1, samples=100]{0};
+\end{axis}
+\end{tikzpicture}