aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-02-20 11:27:40 +0100
committerAdrian Kummerlaender2019-02-20 11:28:20 +0100
commit5f92c1dbea7532068edadbc64b1fcebecd6931f8 (patch)
tree8e939f8485cb8eb4f2b47a958b653f639f566e36
parent46efb766d57f76cebd10819d4bf3c3460bd17264 (diff)
downloadgrid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar.gz
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar.bz2
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar.lz
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar.xz
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.tar.zst
grid_refinement_bsc_thesis-5f92c1dbea7532068edadbc64b1fcebecd6931f8.zip
Localize decimal numbers
-rw-r--r--content.tex26
-rw-r--r--img/common/velocity_fluid.gnuplot1
-rw-r--r--img/cylinder2d_deltap_comparison.tikz5
-rw-r--r--img/cylinder2d_drag_comparison.tikz2
-rw-r--r--img/cylinder2d_high_res_comparsion.tikz21
-rw-r--r--img/cylinder2d_overview.tikz4
-rw-r--r--img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv1613
-rw-r--r--img/data/cylinder2d_generously_refined_n20_re100_drag_lift_deltap.csv1602
-rw-r--r--img/data/cylinder2d_generously_refined_n25_re100_drag_lift_deltap.csv1606
-rw-r--r--img/data/cylinder2d_generously_refined_n30_re100_drag_lift_deltap.csv1602
-rw-r--r--img/data/cylinder2d_optimized_refinement_n40_re100_drag_lift_deltap.csv349
-rw-r--r--img/massloss_interpolation_plot.tikz4
-rw-r--r--img/poiseuille2d_error_comparison.tikz8
-rw-r--r--img/poiseuille2d_grid.tikz6
-rw-r--r--img/poiseuille2d_velocity_grid.tikz9
-rw-r--r--img/poiseuille2d_velocity_outflow.tikz5
-rw-r--r--img/static/cylinder2d_improved_grid_n20_re100_16s_knudsen.pdfbin96793 -> 96793 bytes
-rw-r--r--img/static/cylinder2d_low_resolution_outflow_refine_n5_re100_16s.pdfbin95974 -> 95975 bytes
-rw-r--r--img/static/cylinder2d_optimized_refinement_n5_re100_16s.pdfbin264964 -> 264965 bytes
-rw-r--r--img/static/cylinder2d_optimized_refinement_n5_re100_16s_knudsen.pdfbin91235 -> 91235 bytes
-rw-r--r--img/static/cylinder2d_single_refinement_n20_re100_16s.pdfbin588263 -> 588263 bytes
-rw-r--r--img/static/cylinder2d_single_refinement_n20_re100_16s_knudsen.pdfbin98478 -> 98478 bytes
-rw-r--r--img/static/cylinder2d_unrefined_n12_re100_16s.pdfbin230847 -> 230847 bytes
-rw-r--r--img/static/cylinder2d_unrefined_n20_re100_16s.pdfbin305191 -> 305191 bytes
-rw-r--r--img/static/cylinder2d_unrefined_n20_re100_16s_knudsen.pdfbin80479 -> 80479 bytes
-rw-r--r--img/static/cylinder2d_unrefined_n40_re100_16s.pdfbin895747 -> 895747 bytes
-rw-r--r--img/static/cylinder2d_unrefined_n40_re100_16s_knudsen.pdfbin84684 -> 84684 bytes
-rw-r--r--main.tex7
-rw-r--r--shell.nix1
29 files changed, 6837 insertions, 34 deletions
diff --git a/content.tex b/content.tex
index 0111c9d..3a9152d 100644
--- a/content.tex
+++ b/content.tex
@@ -850,9 +850,8 @@ Mit \(u_\text{max}:=1\) vereinfacht sich damit die analytische Lösung der \(x\)
Das Geschwindigkeitsprofil des Poiseuille-Flusses ist also parabelförmig.
\end{Definition}
-Wir wollen in einem \(1 \times 4\) Längeneinheiten bemessenden Rohr einen Poiseuille-Fluss simulieren. Als Auflösung einer Längeneinheit sei dabei \(N=10\) gewählt, was in der Gitterdiskretisierung durch \(11 \times 21\) grobe und \(21 \times 43\) feine Knoten abgebildet wird.
-
-In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zusammen mit den zugewiesenen Materialzahlen. Wand- und Einflusszellen werden nach dieser Vorlage mit lokalen Geschwindigkeitsrandbedingungen umgesetzt. Während für den Einfluss dabei das Geschwindigkeitsprofil der analytischen Poiseuille-Lösung vorausgesetzt wird, erhält der Ausfluss eine Druckrandbedingung. Die noch verbleibenden gewöhnlichen Fluidzellen erfahren, abgesehen von den üblichen Kollisions- und Strömungsschritten, keine besondere Behandlung.
+Wir wollen in einem \(1 \times 4\) Meter bemessenden Rohr einen Poiseuille-Fluss simulieren. Als Auflösung einer Längeneinheit sei dabei \(N=10\) gewählt, was in der Diskretisierung durch \(11 \times 21\) grobe und \(21 \times 43\) feine Knoten abgebildet wird.
+In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zusammen mit den zugewiesenen Materialzahlen. Wand- und Einflusszellen werden nach dieser Vorlage mit lokalen Geschwindigkeitsrandbedingungen umgesetzt. Während für den Einfluss dabei das Geschwindigkeitsprofil der analytischen Poiseuille-Lösung vorausgesetzt wird, erhält der Ausfluss eine Druckrandbedingung.
\begin{figure}[H]
\begin{adjustbox}{center}
@@ -862,7 +861,7 @@ In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zu
\label{fig:PoiseuilleGridSetup}
\end{figure}
-Neben diesen knotenspezifischen Eigenschaften sei \(u=0.01\) die charateristische Geschwindigkeit in Lattice-Einheiten und \(\text{Re}=10\) die modellierte Reynolds-Zahl. Erstellen wir unsere grobe \class{Grid2D} Instanz mit diesen, die Relaxionszeit \(\tau\) fixierenden, Werten und führen die Simulation bis zur Konvergenz aus, erblicken wir bei geeigneter Aufbereitung in ParaView~\cite{paraview05} schließlich das in Abbildung~\ref{fig:PoiseuilleVelocityGrid} ersichtliche Strömungsbild. Konvergenz bedeutet in diesem Kontext, dass die durchschnittliche Energie des feinen Gitters unter einen Residuumswert, hier \(1\mathrm{e}{-5}\), gefallen ist.
+Neben diesen knotenspezifischen Eigenschaften sei \(u=\num[round-mode=off]{0.01}\) die charateristische Geschwindigkeit in Lattice-Einheiten und \(\text{Re}=10\) die modellierte Reynolds-Zahl. Erstellen wir unsere grobe \class{Grid2D} Instanz mit diesen, die Relaxionszeit \(\tau\) fixierenden, Werten und führen die Simulation bis zur Konvergenz aus, erblicken wir bei geeigneter Aufbereitung in ParaView~\cite{paraview05} schließlich das in Abbildung~\ref{fig:PoiseuilleVelocityGrid} ersichtliche Strömungsbild. Konvergenz bedeutet in diesem Kontext, dass die durchschnittliche Energie des feinen Gitters unter einen Residuumswert, hier \num{1e-5}, gefallen ist.
\begin{figure}[h]
\begin{adjustbox}{center}
@@ -872,7 +871,7 @@ Neben diesen knotenspezifischen Eigenschaften sei \(u=0.01\) die charateristisch
\label{fig:PoiseuilleVelocityGrid}
\end{figure}
-Bei erster Betrachtung lässt sich erkennen, dass die Strömung den Gitterübergang subjektiv ideal bestritten hat. Es treten also keine ungewöhnlichen Artefakte im Geschwindigkeitsbild auf und dieses setzt sich nach dem Übergang in, bis auf die neuen Zwischenwerte, unveränderter Weise fort. Tatsächlich ist bei Interpolation der Knotenzwischenbereiche zur Bildung einer geschlossenen Fläche kein Gitterübergang erkennbar.
+Bei erster Betrachtung lässt sich erkennen, dass die Strömung den Gitterübergang subjektiv ideal bestritten hat. Es treten also keine ungewöhnlichen Artefakte im Geschwindigkeitsbild auf und dieses setzt sich nach dem Übergang in, bis auf die neuen Zwischenwerte, unveränderter Weise fort. Tatsächlich ist bei Bildung einer geschlossenen Fläche durch Interpolation der Zwischenbereiche kein Gitterübergang erkennbar.
\newpage
\subsubsection{Vergleich mit der analytischen Lösung}
@@ -1143,14 +1142,14 @@ Auch im Vergleich mit diesem 145314 Knoten umfassenden Gitter bewähren sich die
& Uniform & Verfeinert & Referenzintervall \cite{SchaeferTurek96} \\
\hline
\hline
-\(\widehat{c_w}\) & 3.28433 & 3.21927 & \([3.22,\ 3.24]\) \\
-\(|\widehat{c_w}-3.23|\) & 0.05433 & 0.01073 & [0, 0.01] \\
+\(\widehat{c_w}\) & \num{3.28433} & \num{3.21927} & \([\num{3.22},\ \num{3.24}]\) \\
+\(|\widehat{c_w}-\num{3.23}|\) & \num{0.05433} & \num{0.01073} & [\num{0}, \num{0.01}] \\
\hline
-\(\widehat{c_a}\) & 1.07046 & 1.10359 & \([0.99,\ 1.01]\) \\
-\(|\widehat{c_a}-1.0|\) & 0.07046 & 0.10359 & [0, 0.01] \\
+\(\widehat{c_a}\) & \num{1.07046} & \num{1.10359} & \([\num{0.99},\ \num{1.01}]\) \\
+\(|\widehat{c_a}-1.0|\) & \num{0.07046} & \num{0.10359} & [\num{0}, \num{0.01}] \\
\hline
-\(\Delta P\) & 2.5793 & 2.44285 & \([2.46,\ 2.5]\) \\
-\(|\Delta P-2.48|\) & 0.0993 & 0.03715 & [0, 0.02] \\
+\(\Delta P\) & \num{2.5793} & \num{2.44285} & \([\num{2.46},\ \num{2.5}]\) \\
+\(|\Delta P-2.48|\) & \num{0.0993} & \num{0.03715} & [\num{0}, \num{0.02}] \\
\hline
\hline
Knotenanzahl & 145314 & 13454 & \\
@@ -1175,5 +1174,10 @@ Dazu sehen wir in Abbildung~\ref{fig:142000nodes} die charakteristischen Messwer
\label{fig:142000nodes}
\end{figure}
+\begin{figure}[H]
+\centering
+\input{img/cylinder2d_generously_refined_comparison.tikz}
+\end{figure}
+
\newpage
\section{Fazit}
diff --git a/img/common/velocity_fluid.gnuplot b/img/common/velocity_fluid.gnuplot
index 5535111..c0db2a9 100644
--- a/img/common/velocity_fluid.gnuplot
+++ b/img/common/velocity_fluid.gnuplot
@@ -7,6 +7,7 @@ set output 'tmp/'.plotname.'.png'
load 'common/moreland.pal'
set datafile separator ','
+set decimalsign ','
set size 1,1
set margin 0,0,0,0
diff --git a/img/cylinder2d_deltap_comparison.tikz b/img/cylinder2d_deltap_comparison.tikz
index 97ecd37..2d6fd99 100644
--- a/img/cylinder2d_deltap_comparison.tikz
+++ b/img/cylinder2d_deltap_comparison.tikz
@@ -15,7 +15,8 @@
xlabel={Simulierte physikalische Zeit},
ylabel={Druckdifferenz},
x unit=s,
- y unit=N/m^2
+ y unit=N/m^2,
+ y tick label style={/pgf/number format/.cd, use comma}
]
\addplot[
@@ -35,7 +36,7 @@
\addlegendentry {Problembezogen verfeinertes \(N=5\) Gitter};
\addplot[color=black]{2.48};
-\addlegendentry {Gemittelte Referenzlösung \(\Delta P := 2.48\)};
+\addlegendentry {Gemittelte Referenzlösung \(\Delta P := \num{2.48}\)};
\end{axis}
\end{tikzpicture}
diff --git a/img/cylinder2d_drag_comparison.tikz b/img/cylinder2d_drag_comparison.tikz
index 98336ec..ee0c783 100644
--- a/img/cylinder2d_drag_comparison.tikz
+++ b/img/cylinder2d_drag_comparison.tikz
@@ -34,7 +34,7 @@
\addlegendentry {Problembezogen verfeinertes \(N=5\) Gitter};
\addplot[color=black]{3.23};
-\addlegendentry {Gemittelte Referenzlösung \(c_\text{Dmax} := 3.23\)};
+\addlegendentry {Gemittelte Referenzlösung \(c_\text{Dmax} := \num{3.23}\)};
\end{axis}
\end{tikzpicture}
diff --git a/img/cylinder2d_high_res_comparsion.tikz b/img/cylinder2d_high_res_comparsion.tikz
index aa8b2a7..9bde49b 100644
--- a/img/cylinder2d_high_res_comparsion.tikz
+++ b/img/cylinder2d_high_res_comparsion.tikz
@@ -1,6 +1,7 @@
\begin{tikzpicture}
\pgfplotstableread[col sep=comma]{img/data/cylinder2d_optimized_refinement_n40_re100_drag_lift_deltap.csv}\refined
-\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n40_re100_drag_lift_deltap.csv}\uniformVeryHigh
+\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}[
@@ -15,11 +16,13 @@
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\)};
+\addlegendentry {\(\widehat{c_w} := \num[round-mode=off]{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}[
@@ -35,11 +38,13 @@
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}[
@@ -58,14 +63,16 @@
y unit=N/m^2,
ylabel absolute, every axis y label/.append style={yshift=0.4cm}
]
+
+\addplot[thick,color=black]{2.48};
+\addlegendentry {\(\Delta P := \num[round-mode=off]{2.48}\)};
\addplot[color=blue!50!white,thin ] table [x expr=\thisrow{time}, y=deltap] {\uniformVeryVeryHigh};
-\addlegendentry {Uniformes \(N=160\) Gitter mit 2298014 Knoten};
+\addlegendentry {Uniformes \(N=160\) Gitter mit \(\sim 2300000\) Knoten};
\addplot[color=red!50!white,thin ] table [x expr=\thisrow{time}, y=deltap] {\uniformVeryHigh};
-\addlegendentry {Uniformes \(N=40\) Gitter mit 576758 Knoten};
+\addlegendentry {Uniformes \(N=80\) Gitter mit \(\sim 577000\) Knoten};
\addplot[color=green!70!black,thin] table [x expr=8*\thisrow{time}, y=deltap] {\refined};
-\addlegendentry {Verfeinertes \(N=40\) Gitter mit 787592 Knoten};
-\addplot[thick,color=black]{2.48};
-\addlegendentry {\(\Delta P := 2.48\)};
+\addlegendentry {Verfeinertes \(N=40\) Gitter mit \(\sim 790000\) Knoten};
+
\end{axis}
\end{tikzpicture}
diff --git a/img/cylinder2d_overview.tikz b/img/cylinder2d_overview.tikz
index 3fa0da8..d77c1c1 100644
--- a/img/cylinder2d_overview.tikz
+++ b/img/cylinder2d_overview.tikz
@@ -10,8 +10,8 @@
\draw[very thick,fill=gray!20!white] (2,2) circle (0.5);
\draw[Latex-Latex,thin] (3,1.5) -- (3,2.5) node[pos=0.5,right] {\(D\)};
-\draw[Latex-Latex,thin] (3,0) -- (3,1.5) node[pos=0.5,right] {\(1.5 D\)};
-\draw[Latex-Latex,thin] (3,2.5) -- (3,4.1) node[pos=0.5,right] {\(1.6 D\)};
+\draw[Latex-Latex,thin] (3,0) -- (3,1.5) node[pos=0.5,right] {\(\num{1.5} D\)};
+\draw[Latex-Latex,thin] (3,2.5) -- (3,4.1) node[pos=0.5,right] {\(\num{1.6} D\)};
\draw[Latex-Latex,thin] (0,-0.5) -- (2,-0.5) node[pos=0.5,below] {\(2 D\)};
\draw[Latex-Latex,thin] (2,-0.5) -- (22,-0.5) node[pos=0.5,below] {\(20 D\)};
diff --git a/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv b/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv
new file mode 100644
index 0000000..df3ffdf
--- /dev/null
+++ b/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv
@@ -0,0 +1,1613 @@
+time,drag,lift,deltap
+0,0,0,0
+0.00124074,2.16707e-07,-1.01615e-09,1.41561e-07
+0.00248148,5.66404e-06,-2.32384e-08,3.48643e-06
+0.00372222,2.07988e-05,-4.13144e-08,1.26355e-05
+0.00496296,4.51912e-05,-7.25249e-08,2.64518e-05
+0.0062037,8.05954e-05,-1.30149e-07,4.81813e-05
+0.00744444,0.000125882,-1.51336e-07,7.29354e-05
+0.00868519,0.000179798,-2.59974e-07,0.000104258
+0.00992593,0.000249092,-2.85972e-07,0.000144944
+0.0111667,0.000322546,-3.82066e-07,0.000182332
+0.0124074,0.000413101,-5.36216e-07,0.000238647
+0.0136481,0.000513467,-4.832e-07,0.000289239
+0.0148889,0.000618901,-8.09399e-07,0.000349125
+0.0161296,0.000750258,-7.48793e-07,0.000425646
+0.0173704,0.000874775,-9.12487e-07,0.000483313
+0.0186111,0.00102628,-1.23767e-06,0.000579434
+0.0198519,0.00118623,-9.88473e-07,0.000654675
+0.0210926,0.00134556,-1.62793e-06,0.000743902
+0.0223333,0.00155243,-1.41798e-06,0.000864309
+0.0235741,0.00174183,-1.65263e-06,0.000946437
+0.0248148,0.00197321,-2.19395e-06,0.00109526
+0.0260556,0.00221027,-1.6465e-06,0.00120166
+0.0272963,0.00243669,-2.70767e-06,0.00132701
+0.028537,0.002725,-2.26102e-06,0.00149418
+0.0297778,0.00296786,-2.56064e-06,0.00158973
+0.0310185,0.00326876,-3.37514e-06,0.00178677
+0.0322593,0.00357276,-2.42635e-06,0.00191507
+0.0335,0.00385458,-3.98675e-06,0.00206945
+0.0347407,0.00422257,-3.24489e-06,0.00228287
+0.0359815,0.00451762,-3.62553e-06,0.00238808
+0.0372222,0.00489117,-4.7446e-06,0.00263737
+0.038463,0.00526355,-3.3024e-06,0.00278564
+0.0397037,0.00559963,-5.4718e-06,0.00296833
+0.0409444,0.00605,-4.35728e-06,0.00323001
+0.0421852,0.00639465,-4.81179e-06,0.00334038
+0.0434259,0.00683807,-6.30127e-06,0.00364226
+0.0446667,0.00726866,-4.2744e-06,0.00380186
+0.0459074,0.00764363,-7.12343e-06,0.00400396
+0.0471481,0.00816632,-5.58634e-06,0.00430873
+0.0483889,0.00854746,-6.1322e-06,0.00441428
+0.0496296,0.00905768,-8.02073e-06,0.00476959
+0.0508704,0.0095552,-5.3359e-06,0.00494393
+0.0521111,0.00998171,-8.96114e-06,0.00517342
+0.0533519,0.0105919,-6.93372e-06,0.00553075
+0.0545926,0.0110188,-7.57015e-06,0.00563438
+0.0558333,0.0116049,-9.90253e-06,0.00605014
+0.0570741,0.0121718,-6.4789e-06,0.00623808
+0.0583148,0.0126479,-1.09423e-05,0.00649377
+0.0595556,0.0133479,-8.36819e-06,0.00690521
+0.0607963,0.0138195,-9.11157e-06,0.00700391
+0.062037,0.014484,-1.19074e-05,0.00748331
+0.0632778,0.0151227,-7.66954e-06,0.00768405
+0.0645185,0.0156511,-1.30517e-05,0.00796762
+0.0657593,0.0164518,-9.86902e-06,0.00843995
+0.067,0.01698,-1.07163e-05,0.0085384
+0.0682407,0.0177372,-1.40013e-05,0.00909125
+0.0694815,0.0184582,-8.87881e-06,0.00930843
+0.0707222,0.0190443,-1.52312e-05,0.00962269
+0.071963,0.019944,-1.13646e-05,0.0101544
+0.0732037,0.0205145,-1.2334e-05,0.0102417
+0.0744444,0.0213522,-1.61241e-05,0.0108622
+0.0756852,0.0221451,-1.00421e-05,0.0110877
+0.0769259,0.0227795,-1.74376e-05,0.0114277
+0.0781667,0.0237761,-1.28308e-05,0.0120179
+0.0794074,0.0243875,-1.39253e-05,0.0120918
+0.0806481,0.0253069,-1.82353e-05,0.0127823
+0.0818889,0.0261708,-1.11377e-05,0.013014
+0.0831296,0.0268508,-1.96262e-05,0.0133787
+0.0843704,0.027943,-1.42167e-05,0.0140268
+0.0856111,0.0285895,-1.54563e-05,0.0140825
+0.0868519,0.0295842,-2.02904e-05,0.0148408
+0.0880926,0.0305066,-1.21208e-05,0.0150697
+0.0893333,0.0312174,-2.17588e-05,0.0154511
+0.0905741,0.032396,-1.54919e-05,0.0161524
+0.0918148,0.0330688,-1.68894e-05,0.0161839
+0.0930556,0.0341394,-2.2253e-05,0.0170124
+0.0942963,0.0351289,-1.2973e-05,0.0172427
+0.095537,0.0358826,-2.38004e-05,0.017649
+0.0967778,0.0371586,-1.66095e-05,0.01841
+0.0980185,0.0378634,-1.81946e-05,0.0184195
+0.0992593,0.0390147,-2.40723e-05,0.0193224
+0.1005,0.0400709,-1.36239e-05,0.0195521
+0.101741,0.0408651,-2.56872e-05,0.0199825
+0.102981,0.0422392,-1.75037e-05,0.0208036
+0.104222,0.042975,-1.92885e-05,0.0207897
+0.105463,0.0442098,-2.56592e-05,0.0217697
+0.106704,0.0453352,-1.40047e-05,0.0219989
+0.107944,0.0461746,-2.73359e-05,0.0224571
+0.109185,0.0476567,-1.80808e-05,0.0233438
+0.110426,0.0484344,-2.00899e-05,0.0233122
+0.111667,0.0497636,-2.69241e-05,0.0243766
+0.112907,0.0509635,-1.4011e-05,0.0246067
+0.114148,0.0518488,-2.86271e-05,0.0250936
+0.115389,0.0534327,-1.82255e-05,0.0260415
+0.11663,0.0542379,-2.04757e-05,0.0259829
+0.11787,0.0556505,-2.77208e-05,0.0271263
+0.119111,0.0569158,-1.35184e-05,0.0273509
+0.120352,0.0578412,-2.94426e-05,0.0278645
+0.121593,0.0595255,-1.78016e-05,0.0288728
+0.122833,0.0603573,-2.03241e-05,0.0287867
+0.124074,0.0618542,-2.79341e-05,0.030011
+0.125315,0.0631824,-1.23953e-05,0.0302278
+0.126556,0.0641452,-2.96405e-05,0.0307678
+0.127796,0.0659263,-1.66796e-05,0.0318342
+0.129037,0.0667778,-1.94912e-05,0.0317167
+0.130278,0.0683517,-2.73988e-05,0.0330189
+0.131519,0.0697303,-1.04963e-05,0.0332198
+0.132759,0.0707192,-2.90649e-05,0.0337816
+0.134,0.0725904,-1.46902e-05,0.0349022
+0.135241,0.0734578,-1.78219e-05,0.0347523
+0.136481,0.0751134,-2.59559e-05,0.0361368
+0.137722,0.0765515,-7.65761e-06,0.0363268
+0.138963,0.0775787,-2.75388e-05,0.0369195
+0.140204,0.0795487,-1.16557e-05,0.0380992
+0.141444,0.0804362,-1.51265e-05,0.0379199
+0.142685,0.0821772,-2.3385e-05,0.0393898
+0.143926,0.0836735,-3.66457e-06,0.0395677
+0.145167,0.084739,-2.48355e-05,0.0401931
+0.146407,0.0868092,-7.32944e-06,0.0414321
+0.147648,0.0877182,-1.11659e-05,0.0412254
+0.148889,0.0895499,-1.94374e-05,0.0427851
+0.15013,0.091109,1.73986e-06,0.0429529
+0.15137,0.0922212,-2.06875e-05,0.0436182
+0.152611,0.0944018,-1.43766e-06,0.0449221
+0.153852,0.0953435,-5.65762e-06,0.0446956
+0.155093,0.0972761,-1.3807e-05,0.0463517
+0.156333,0.0989004,8.86168e-06,0.0465103
+0.157574,0.100059,-1.47759e-05,0.0472176
+0.158815,0.102344,6.35132e-06,0.0485818
+0.160056,0.103308,1.73404e-06,0.0483308
+0.161296,0.105335,-6.15241e-06,0.0500808
+0.162537,0.107019,1.80359e-05,0.0502274
+0.163778,0.108225,-6.74843e-06,0.0509797
+0.165019,0.110615,1.63867e-05,0.0524048
+0.166259,0.111603,1.13511e-05,0.052133
+0.1675,0.113728,3.89304e-06,0.0539805
+0.168741,0.115472,2.96285e-05,0.0541149
+0.169981,0.116726,3.7627e-06,0.0549158
+0.171222,0.11922,2.90447e-05,0.0564004
+0.172463,0.120228,2.35806e-05,0.0561081
+0.173704,0.122448,1.67192e-05,0.0580527
+0.174944,0.124242,4.40192e-05,0.0581709
+0.176185,0.125541,1.71588e-05,0.0590224
+0.177426,0.128137,4.47214e-05,0.0605665
+0.178667,0.12917,3.88034e-05,0.0602601
+0.179907,0.131496,3.2729e-05,0.0623107
+0.181148,0.133354,6.16013e-05,0.0624213
+0.182389,0.134714,3.38364e-05,0.0633368
+0.18363,0.137424,6.3826e-05,0.0649464
+0.18487,0.138489,5.74355e-05,0.0646338
+0.186111,0.140929,5.23454e-05,0.0667966
+0.187352,0.142854,8.27933e-05,0.0669027
+0.188593,0.144283,5.42296e-05,0.0678908
+0.189833,0.147112,8.67805e-05,0.0695699
+0.191074,0.148219,7.98869e-05,0.0692606
+0.192315,0.150786,7.59937e-05,0.0715448
+0.193556,0.152788,0.000108,0.0716543
+0.194796,0.154302,7.87408e-05,0.072729
+0.196037,0.157267,0.000113987,0.0744867
+0.197278,0.158431,0.000106555,0.0741943
+0.198519,0.16114,0.000104067,0.0766099
+0.199759,0.163225,0.000137599,0.0767278
+0.201,0.16483,0.000107754,0.0778974
+0.202241,0.167929,0.000145808,0.0797333
+0.203481,0.169151,0.000137778,0.0794623
+0.204722,0.172004,0.000136895,0.0820138
+0.205963,0.174176,0.000171886,0.0821448
+0.207204,0.175884,0.000141544,0.0834213
+0.208444,0.179126,0.000182487,0.0853412
+0.209685,0.180415,0.000173775,0.0851037
+0.210926,0.183425,0.000174679,0.0877999
+0.212167,0.18569,0.000211018,0.0879505
+0.213407,0.18751,0.000180254,0.089345
+0.214648,0.1909,0.000224131,0.0913523
+0.215889,0.192261,0.000214617,0.0911574
+0.21713,0.195432,0.000217456,0.0940042
+0.21837,0.197792,0.000254983,0.0941787
+0.219611,0.199734,0.000223841,0.0957028
+0.220852,0.203281,0.000270635,0.0978044
+0.222093,0.204731,0.000260156,0.0976687
+0.223333,0.208086,0.000265042,0.100681
+0.224574,0.210559,0.000303528,0.100894
+0.225815,0.212648,0.000272025,0.102568
+0.227056,0.216368,0.000321668,0.104774
+0.228296,0.217923,0.000310009,0.104713
+0.229537,0.221476,0.000317009,0.107902
+0.230778,0.224074,0.000356169,0.108164
+0.232019,0.226327,0.00032427,0.110003
+0.233259,0.230235,0.000376622,0.112322
+0.2345,0.231913,0.000363516,0.112353
+0.235741,0.235686,0.000372638,0.115732
+0.236981,0.238427,0.000412111,0.116057
+0.238222,0.240868,0.000379736,0.11808
+0.239463,0.244983,0.000434588,0.120524
+0.240704,0.246806,0.00041971,0.120666
+0.241944,0.250817,0.000430902,0.124245
+0.243185,0.253711,0.000470261,0.124642
+0.244426,0.256352,0.000437279,0.126861
+0.245667,0.260677,0.000494338,0.12943
+0.246907,0.262652,0.000477314,0.129692
+0.248148,0.266909,0.00049045,0.133479
+0.249389,0.269965,0.000529189,0.133955
+0.25063,0.272823,0.000495418,0.136383
+0.25187,0.27737,0.000554312,0.139085
+0.253111,0.27951,0.000534713,0.139481
+0.254352,0.284025,0.000549607,0.14348
+0.255593,0.287249,0.00058716,0.144042
+0.256833,0.290334,0.000552374,0.146689
+0.258074,0.295105,0.000612661,0.149522
+0.259315,0.297417,0.000590017,0.150061
+0.260556,0.302194,0.000606422,0.154275
+0.261796,0.305588,0.000642169,0.154927
+0.263037,0.30891,0.000606098,0.1578
+0.264278,0.313912,0.00066727,0.16077
+0.265519,0.316411,0.000641083,0.161465
+0.266759,0.321467,0.000658693,0.165902
+0.268,0.325044,0.000691973,0.166654
+0.269241,0.328619,0.000654324,0.169767
+0.270481,0.333859,0.000715818,0.172875
+0.271722,0.336553,0.000685565,0.173737
+0.272963,0.341892,0.000704033,0.178397
+0.274204,0.345656,0.00073414,0.179256
+0.275444,0.349494,0.000694581,0.182615
+0.276685,0.354975,0.000755775,0.185863
+0.277926,0.357875,0.000720902,0.186903
+0.279167,0.363506,0.000739816,0.191793
+0.280407,0.367463,0.000765998,0.192765
+0.281648,0.371575,0.000724154,0.196381
+0.282889,0.377303,0.000784353,0.199771
+0.28413,0.380415,0.000744265,0.201
+0.28537,0.38634,0.000763142,0.206118
+0.286611,0.390485,0.000784576,0.207206
+0.287852,0.394869,0.00074001,0.211081
+0.289093,0.400832,0.000798427,0.214609
+0.290333,0.404153,0.000752458,0.216031
+0.291574,0.410366,0.000770716,0.221376
+0.292815,0.414697,0.000786494,0.222584
+0.294056,0.419355,0.000738687,0.226723
+0.295296,0.425548,0.000794425,0.230388
+0.296537,0.429078,0.000741845,0.232012
+0.297778,0.435575,0.000758805,0.237582
+0.299019,0.440085,0.000767937,0.238912
+0.300259,0.445013,0.000716304,0.243318
+0.3015,0.451426,0.000768375,0.247116
+0.302741,0.455161,0.000708402,0.248946
+0.303981,0.461933,0.000723312,0.254737
+0.305222,0.466615,0.000724756,0.256193
+0.306463,0.471815,0.000668679,0.26087
+0.307704,0.478445,0.000716055,0.264803
+0.308944,0.482392,0.000647914,0.266852
+0.310185,0.489444,0.000660016,0.272867
+0.311426,0.494302,0.000652768,0.274458
+0.312667,0.499779,0.000591684,0.279414
+0.313907,0.506624,0.00063339,0.283482
+0.315148,0.510785,0.000556418,0.285758
+0.316389,0.518115,0.000565046,0.291996
+0.31763,0.523149,0.000548242,0.293728
+0.31887,0.528908,0.000481743,0.298967
+0.320111,0.535967,0.000516971,0.303173
+0.321352,0.540351,0.000430726,0.305687
+0.322593,0.547962,0.000435427,0.31215
+0.323833,0.553178,0.000408473,0.314033
+0.325074,0.559226,0.000336431,0.319562
+0.326315,0.5665,0.000364667,0.323907
+0.327556,0.571111,0.00026906,0.326669
+0.328796,0.579,0.000269706,0.333353
+0.330037,0.584392,0.000232391,0.33539
+0.331278,0.590726,0.000155067,0.341208
+0.332519,0.598205,0.000176191,0.345689
+0.333759,0.603043,7.157e-05,0.348702
+0.335,0.611204,6.84533e-05,0.355605
+0.336241,0.616773,2.10385e-05,0.3578
+0.337481,0.623394,-6.08444e-05,0.363909
+0.338722,0.631072,-4.6476e-05,0.368525
+0.339963,0.636137,-0.000159245,0.371794
+0.341204,0.644564,-0.000165353,0.37891
+0.342444,0.650302,-0.000222078,0.381266
+0.343685,0.657205,-0.000307288,0.387663
+0.344926,0.665072,-0.000298809,0.392409
+0.346167,0.670359,-0.000418328,0.395938
+0.347407,0.679042,-0.00042616,0.403259
+0.348648,0.684946,-0.000490873,0.405779
+0.349889,0.692129,-0.00057768,0.412463
+0.35113,0.700182,-0.000573721,0.41734
+0.35237,0.705697,-0.000698079,0.421135
+0.353611,0.714635,-0.000705903,0.428659
+0.354852,0.720706,-0.000776786,0.431349
+0.356093,0.72817,-0.000863008,0.43832
+0.357333,0.736401,-0.000861764,0.443324
+0.358574,0.742144,-0.000988612,0.447389
+0.359815,0.751329,-0.000994324,0.455109
+0.361056,0.757564,-0.00106917,0.457974
+0.362296,0.765307,-0.00115229,0.465229
+0.363537,0.773713,-0.00115165,0.47036
+0.364778,0.779687,-0.00127835,0.474698
+0.366019,0.789116,-0.00127962,0.482611
+0.367259,0.795516,-0.00135599,0.485654
+0.3685,0.803538,-0.00143337,0.493191
+0.369741,0.812112,-0.00143108,0.498448
+0.370981,0.818314,-0.00155492,0.503061
+0.372222,0.827974,-0.00154943,0.511156
+0.373463,0.834531,-0.00162489,0.514377
+0.374704,0.84282,-0.00169398,0.522186
+0.375944,0.851548,-0.00168796,0.52756
+0.377185,0.857972,-0.00180641,0.532446
+0.378426,0.867852,-0.00179213,0.540714
+0.379667,0.87456,-0.00186457,0.544113
+0.380907,0.883108,-0.00192323,0.552188
+0.382148,0.891979,-0.00191183,0.557675
+0.383389,0.898619,-0.00202289,0.56283
+0.38463,0.908703,-0.00199839,0.571258
+0.38587,0.915551,-0.00206633,0.574833
+0.387111,0.924346,-0.00211317,0.583162
+0.388352,0.933345,-0.00209553,0.588754
+0.389593,0.94019,-0.00219807,0.594172
+0.390833,0.950464,-0.00216285,0.60275
+0.392074,0.957446,-0.00222577,0.606499
+0.393315,0.966478,-0.00226047,0.615073
+0.394556,0.975598,-0.00223681,0.620768
+0.395796,0.98265,-0.0023308,0.626449
+0.397037,0.993104,-0.00228559,0.635167
+0.398278,1.00022,-0.00234416,0.639091
+0.399519,1.00948,-0.00236767,0.647901
+0.400759,1.01871,-0.00233948,0.653692
+0.402,1.02596,-0.00242622,0.659631
+0.403241,1.03658,-0.00237305,0.668478
+0.404481,1.04382,-0.00242924,0.672577
+0.405722,1.0533,-0.00244384,0.681613
+0.406963,1.06263,-0.0024139,0.687497
+0.408204,1.07008,-0.00249595,0.693692
+0.409444,1.08086,-0.00243813,0.702659
+0.410685,1.08822,-0.00249508,0.706932
+0.411926,1.09792,-0.00250424,0.716184
+0.413167,1.10734,-0.00247637,0.722156
+0.414407,1.11498,-0.00255732,0.728598
+0.415648,1.1259,-0.00249905,0.737672
+0.416889,1.13337,-0.00256071,0.742113
+0.41813,1.14326,-0.00256862,0.751565
+0.41937,1.15276,-0.00254723,0.757616
+0.420611,1.16057,-0.00263113,0.764297
+0.421852,1.17162,-0.00257699,0.773464
+0.423093,1.17919,-0.00264748,0.77807
+0.424333,1.18927,-0.00265838,0.787708
+0.425574,1.19883,-0.00264774,0.793831
+0.426815,1.20681,-0.00273835,0.800739
+0.428056,1.21796,-0.00269243,0.809985
+0.429296,1.22562,-0.00277519,0.814748
+0.430537,1.23586,-0.00279253,0.824554
+0.431778,1.24547,-0.00279587,0.83074
+0.433019,1.25361,-0.00289579,0.837864
+0.434259,1.26484,-0.00286082,0.847174
+0.4355,1.27259,-0.00295774,0.852089
+0.436741,1.28298,-0.00298328,0.86205
+0.437981,1.29263,-0.00300194,0.868294
+0.439222,1.30091,-0.00311174,0.875625
+0.440463,1.31222,-0.00308828,0.884988
+0.441704,1.32004,-0.00319896,0.890051
+0.442944,1.33056,-0.00323207,0.900151
+0.444185,1.34024,-0.00326484,0.906446
+0.445426,1.34865,-0.00338255,0.913972
+0.446667,1.36003,-0.00336853,0.923376
+0.447907,1.36791,-0.00348986,0.928582
+0.449148,1.37856,-0.00352727,0.938806
+0.450389,1.38826,-0.0035703,0.945147
+0.45163,1.3968,-0.00369135,0.95286
+0.45287,1.40822,-0.00368217,0.962292
+0.454111,1.41617,-0.00380864,0.967636
+0.455352,1.42693,-0.00384482,0.97797
+0.456593,1.43663,-0.0038921,0.984349
+0.457833,1.44528,-0.00401001,0.992233
+0.459074,1.45673,-0.00399937,1.00168
+0.460315,1.46472,-0.004124,1.00715
+0.461556,1.47557,-0.00415225,1.01757
+0.462796,1.48526,-0.0041969,1.02398
+0.464037,1.49402,-0.00430468,1.03202
+0.465278,1.50548,-0.00428607,1.04147
+0.466519,1.51351,-0.00440206,1.04706
+0.467759,1.52443,-0.00441626,1.05755
+0.469,1.5341,-0.0044523,1.06398
+0.470241,1.54294,-0.0045443,1.07216
+0.471481,1.55438,-0.00451292,1.08159
+0.472722,1.56243,-0.00461558,1.08729
+0.473963,1.5734,-0.00461213,1.09783
+0.475204,1.58303,-0.00463642,1.10427
+0.476444,1.59194,-0.00471026,1.11257
+0.477685,1.60335,-0.00466487,1.12198
+0.478926,1.61142,-0.00475335,1.12777
+0.480167,1.62242,-0.00473279,1.13835
+0.481407,1.632,-0.00474653,1.14479
+0.482648,1.64096,-0.00480435,1.1532
+0.483889,1.65234,-0.00474833,1.16256
+0.48513,1.6604,-0.00482652,1.16844
+0.48637,1.67142,-0.00479411,1.17904
+0.487611,1.68094,-0.00480316,1.18547
+0.488852,1.68994,-0.00485172,1.19398
+0.490093,1.70125,-0.00479285,1.20328
+0.491333,1.70932,-0.0048688,1.20924
+0.492574,1.72033,-0.00483369,1.21984
+0.493815,1.72978,-0.00484737,1.22626
+0.495056,1.73881,-0.00489648,1.23484
+0.496296,1.75006,-0.00484506,1.24408
+0.497537,1.7581,-0.00492881,1.25012
+0.498778,1.7691,-0.00490151,1.26069
+0.500019,1.77847,-0.00492987,1.26709
+0.501259,1.78751,-0.00498939,1.27573
+0.5025,1.79866,-0.004955,1.28489
+0.503741,1.80668,-0.00505516,1.29098
+0.504981,1.81764,-0.00504402,1.30151
+0.506222,1.8269,-0.00509411,1.30789
+0.507463,1.83595,-0.00517024,1.31657
+0.508704,1.84699,-0.00515803,1.32563
+0.509944,1.85496,-0.00527813,1.33176
+0.511185,1.86586,-0.00528573,1.34224
+0.512426,1.87502,-0.00535828,1.34857
+0.513667,1.88403,-0.00545037,1.35727
+0.514907,1.89495,-0.00545822,1.36623
+0.516148,1.90286,-0.00559418,1.37239
+0.517389,1.91368,-0.0056153,1.38278
+0.51863,1.92271,-0.0057031,1.38907
+0.51987,1.93168,-0.00580265,1.39776
+0.521111,1.94245,-0.00582064,1.4066
+0.522352,1.9503,-0.00596092,1.41278
+0.523593,1.96101,-0.00598325,1.42308
+0.524833,1.96991,-0.00607252,1.4293
+0.526074,1.97883,-0.00616516,1.43798
+0.527315,1.98943,-0.0061783,1.44668
+0.528556,1.9972,-0.00630722,1.45286
+0.529796,2.0078,-0.00631535,1.46305
+0.531037,2.01655,-0.00639036,1.4692
+0.532278,2.0254,-0.00646099,1.47785
+0.533519,2.03583,-0.00645481,1.48641
+0.534759,2.04351,-0.00655858,1.49258
+0.536,2.05398,-0.00654033,1.50264
+0.537241,2.06257,-0.00658993,1.50872
+0.538481,2.07133,-0.00662941,1.51731
+0.539722,2.08158,-0.00659665,1.52572
+0.540963,2.08916,-0.00666997,1.53187
+0.542204,2.09948,-0.00662273,1.54179
+0.543444,2.10791,-0.00664636,1.54778
+0.544685,2.11657,-0.00665699,1.5563
+0.545926,2.12662,-0.00660249,1.56455
+0.547167,2.13408,-0.00665267,1.57066
+0.548407,2.14423,-0.00658671,1.58041
+0.549648,2.15248,-0.00659671,1.58631
+0.550889,2.16102,-0.00659353,1.59474
+0.55213,2.17085,-0.0065344,1.60281
+0.55337,2.17818,-0.0065803,1.60887
+0.554611,2.18815,-0.00651647,1.61844
+0.555852,2.19621,-0.00653458,1.62424
+0.557093,2.2046,-0.00654057,1.63256
+0.558333,2.2142,-0.00650013,1.64044
+0.559574,2.22139,-0.00656505,1.64643
+0.560815,2.23115,-0.00652654,1.6558
+0.562056,2.23901,-0.00657464,1.66148
+0.563296,2.24725,-0.00661068,1.66968
+0.564537,2.2566,-0.00660781,1.67736
+0.565778,2.26364,-0.00670823,1.68327
+0.567019,2.27317,-0.00670906,1.69243
+0.568259,2.28082,-0.00679767,1.69799
+0.5695,2.28889,-0.00687111,1.70604
+0.570741,2.29799,-0.0069092,1.71351
+0.571981,2.30486,-0.00704446,1.71934
+0.573222,2.31416,-0.00708008,1.72827
+0.574463,2.32159,-0.00720001,1.7337
+0.575704,2.32948,-0.00729782,1.74159
+0.576944,2.33831,-0.00735953,1.74885
+0.578185,2.34501,-0.00750832,1.75457
+0.579426,2.35406,-0.00755406,1.76327
+0.580667,2.36127,-0.00767705,1.76855
+0.581907,2.36896,-0.00776857,1.77627
+0.583148,2.37752,-0.0078208,1.78331</