From e440d98649dfc7cf36fadd485c7762612f12d05c Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Fri, 18 Jan 2019 17:54:08 +0100 Subject: Finally start evaluation section --- content.tex | 26 ++++++++++++++++++++++++-- img/poiseuille2d_setup.tikz | 2 +- img/static/refined_poiseuille_1.png | Bin 33877 -> 32989 bytes main.tex | 2 ++ shell.nix | 1 + 5 files changed, 28 insertions(+), 3 deletions(-) diff --git a/content.tex b/content.tex index f58c802..2b78267 100644 --- a/content.tex +++ b/content.tex @@ -790,6 +790,15 @@ Hiermit sind die zentralen Bestandteile der Umsetzung des Verfahrens von Lagrava \newpage \section{Evaluierung} +Die Auswahl von Beispielen zur Evaluation der Qualität eines Gitterverfeinerungsverfahrens gestaltet sich ohne Detailkenntnis des Versuchsaufbaus eines real zu simulierenden Strömungsproblem zunächst unklarer, als man annehmen könnte. So stellt OpenLB zwar eine gute Auswahl verschiedener Simulationsbeispiele bereit, aber nur wenige von ihnen beinhalten auch analytische Lösungen oder auch nur einfach zu verwendende Vergleichsdaten realer Versuche. + +Unter diesen Einschränkungen ist es -- abseits offensichtlicher Gütekriterien wie dem Ausschluss divergierender Simulationen -- schwer zu sagen, ob zwei auf verschiedenen Weisen simulierte Lösungen nun besser oder schlechter als die jeweils andere Lösung sind. Entsprechend beschränken wir uns auf die Betrachtung der folgenden Gütekriterien: +\begin{enumerate} + \item Subjektive Qualität des Strömungsbildes + \item Erhalt von Masse im Gitterübergang + \item Stetigkeit der Momente im Gitterübergang +\end{enumerate} + \subsection{Rohrströmung} Als Einstiegspunkt wollen wir zunächst die grundsätzliche Funktion des Verfeinerungsverfahren an einem möglichst einfachen Beispiel gegen möglichst korrekte Vergleichsdaten testen. Ein solches Beispiel ist gegeben durch die laminare Strömung in einem geradenen Rohr mit kreisförmigem Querschnitt und ohne jegliche Hindernisse abseits der Wände. @@ -797,12 +806,25 @@ Als Einstiegspunkt wollen wir zunächst die grundsätzliche Funktion des Verfein Eine solche Rohrströmung stellt nicht nur eine der denkbar einfachsten Strömungssituationen dar, sondern besitzt als Poiseuille Fluss auch eine analytische Lösung, so dass wir ideale Vergleichbedingungen erhalten. Lieferte unser Verfahren hier keine guten Ergebnisse, wäre nicht davon auszugehen, dass dies sich in komplexeren Situationen verbessern würde. \begin{figure}[h] -\centering -%\includegraphics[scale=0.35]{img/static/poiseuille2d_materialzahlen.png} +\begin{adjustbox}{center} \input{img/poiseuille2d_setup.tikz} +\end{adjustbox} \caption{Gitterstruktur einer halbseitig verfeinerten Rohrströmung} +\label{fig:PoiseuilleSetup} \end{figure} +Wir wollen in einem Rohr mit Dimensionen \(1 \times 4\) in Längeneinheiten einen Poiseuille Fluss simulieren. Als Auflösung der 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 Abbilddung~\ref{fig:PoiseuilleSetup} sehen wir das resultierende Gitter zusammen mit den zugewiesenen Materialzahlen. + +Diese Materialzahlen werden dabei wie folgt an das Strömungsverhalten modellierende \class{Dynamics} gebunden: Wandzellen werden mit einfachen Bounce-Back Randbedingungen umgesetzt. Weiter wird für den Einfluss wird eine Geschwindigkeitsrandbedingung mit dem analytischen Poiseuille Geschwindigkeitsprofil gesetzt. Der Ausfluss erhält eine korrespondierende Druckrandbedingung. Die verbleibenden Fluidzellen werden mit Ausnahme der Durchführung des, allen Gitterzellen dieser Arbeit gemeinen, D2Q9 BGK Kollisions- und Strömungsschritt nicht besonders behandelt. + +\begin{figure}[h] +\includegraphics[width=1.1\textwidth,center]{img/static/refined_poiseuille_1.png} +\caption{Norm der Geschwindigkeiten im simulierten Poiseuille Fluss} +\label{fig:PoiseuilleSim} +\end{figure} + +Neben dieses knotenspezifischen Materialeigenschaften sei \(0.01\) die charateristische Geschwindigkeit in Lattice-Einheiten und \(10\) die modellierte Reynolds-Zahl. Konstruieren wir unsere grobe \class{Grid2D} Instanz mit diesen Werten und führen die Simulation aus, erblicken wir bei Darstellung in ParaView schließlich das in Abbildung~\ref{fig:PoiseuilleSim} ersichtliche Strömungsbild. + \subsubsection{Vergleich der Interpolationsverfahren} \begin{figure}[h] diff --git a/img/poiseuille2d_setup.tikz b/img/poiseuille2d_setup.tikz index c4e50f5..76b307d 100644 --- a/img/poiseuille2d_setup.tikz +++ b/img/poiseuille2d_setup.tikz @@ -1,5 +1,5 @@ \begin{tikzpicture}[ - scale=0.375, + scale=1.1*\textwidth/41cm, coarse/.style={circle,thick,inner sep=2}, fine/.style={cross out,thick,inner sep=0.75}, legend/.style={minimum height=1,minimum width=1} diff --git a/img/static/refined_poiseuille_1.png b/img/static/refined_poiseuille_1.png index fee3996..7d2453e 100644 Binary files a/img/static/refined_poiseuille_1.png and b/img/static/refined_poiseuille_1.png differ diff --git a/main.tex b/main.tex index c23bc63..a5c09c9 100644 --- a/main.tex +++ b/main.tex @@ -46,6 +46,8 @@ \tikzexternalize[prefix=tikzcache/] +\usepackage[export]{adjustbox} + \setlength{\topmargin}{-15mm} \input{commands.tex} diff --git a/shell.nix b/shell.nix index 33de630..1985f45 100644 --- a/shell.nix +++ b/shell.nix @@ -14,6 +14,7 @@ stdenv.mkDerivation rec { minted fvextra ifplatform framed bibtex biblatex logreq xstring pgfplots + adjustbox collectbox ; }; in [ -- cgit v1.2.3