From f9c0f70e1d88c0373821c0d4974a6d8de27cc09b Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Wed, 13 Mar 2019 17:50:06 +0100 Subject: Cite multi-grid, multi-domain users --- content.tex | 28 ++++++++--- main.tex | 2 +- quellen.bib | 154 ++++++++++++++++++++++++++++++++++++++++++++++++++++-------- 3 files changed, 157 insertions(+), 27 deletions(-) diff --git a/content.tex b/content.tex index 6c05ed4..6ac95a2 100644 --- a/content.tex +++ b/content.tex @@ -59,7 +59,7 @@ Sei \(\tau \in \R_{\geq 0}\) eine Relaxationszeit und \(f^\text{eq}\) die von de \[ (\partial_t + \xi \cdot \nabla_x) f = -\frac{1}{\tau} (f(x,\xi,t) - f^\text{eq}(x,\xi,t)) .\] \end{Definition} -Analog zur Boltzmann-Gleichung ist auch bei deren BGK Approximation der beschriebene Ort \(x \in \R^2\) im Allgemeinen frei gewählt. Da unser Ziel jedoch gerade die Diskretisierung der Simulationsdomäne in einem Gitter ist, wollen wir \(x\) einschränken: +Analog zur Boltzmann-Gleichung ist auch bei deren Approximation in der BGK-Boltzmann-Gleichung der beschriebene Ort \(x \in \R^2\) im Allgemeinen frei gewählt. Da unser Ziel jedoch gerade die Diskretisierung der Simulationsdomäne in einem Gitter ist, wollen wir \(x\) einschränken: \begin{Definition}[Ortsdiskretisierung] \label{def:SpatialDiscretizationLBM} @@ -74,7 +74,7 @@ Für die verbleibende Herleitung der LBM können wir diese Interpretation, d. h. \bigskip -Wir bemerken nun, dass die BGK Approximation nicht nur für beliebige Orte, sondern auch für beliebige Geschwindigkeiten \(\xi \in \R^2\) definiert ist. Da wir die LBM auf einem endlichen Rechner umsetzen wollen, müssen wir auch die Menge der betrachteten Geschwindigkeiten auf eine endliche Menge diskretisieren. +Wir bemerken nun, dass die BGK Approximation nicht nur für beliebige Orte, sondern auch für beliebige Geschwindigkeiten \(\xi \in \R^2\) definiert ist. Da wir die LBM auf einem Computer umsetzen wollen, müssen wir die Menge der betrachteten Geschwindigkeiten auf eine endliche Menge diskretisieren. Eine übliche Menge diskreter Geschwindigkeiten in 2D ist \emph{D2Q9} wobei \emph{D2} die Anzahl der Dimensionen und \emph{Q9} die Anzahl der Geschwindigkeiten verschlüsselt. @@ -223,11 +223,15 @@ Diese Darstellung können wir unter Verwendung von Definition~\ref{def:ChapmanEn \newpage \subsection{Herangehensweisen an Gitterverfeinerung} -Grundsätzlich existieren mit der \emph{Multi-Grid} und \emph{Multi-Domain} Herangehenweise zwei verschiedene Ansätze für Gitterverfeinerung in LBM \cite[Kap.~3.1]{Lagrava12}. Im Wesentlichen unterscheiden die Ansätze sich in den Ausmaßen der variabel aufgelösten Teilgitter der Simulationsdomäne. Weitere Unterschiede folgen dann aus dieser grundlegenden Struktur. +Grundsätzlich existieren mit der \emph{Multi-Grid} +und \emph{Multi-Domain} Herangehenweise zwei verschiedene Ansätze für Gitterverfeinerung in LBM \cite[Kap.~3.1]{Lagrava12}. Im Wesentlichen unterscheiden die Ansätze sich in den Ausmaßen der variabel aufgelösten Teilgitter der Simulationsdomäne. Weitere Unterschiede folgen dann aus dieser grundlegenden Struktur. \subsubsection{Multi-Grid Ansatz} -Bei Verfahren des Multi-Grid Ansatzes existiert das gröbste Gitter über der gesamten Domäne. Feinere Teilgitter werden über gröberen Gittern platziert und nicht aus deren Verarbeitung ausgeschlossen. Somit existieren über der gesamten \emph{Fläche} feinerer Gitter auch die Knoten gröberer Gitter. +Bei Verfahren des Multi-Grid Ansatzes +\cite{Lin00} +\cite{Toelke02} +existiert das gröbste Gitter über der gesamten Domäne. Feinere Teilgitter werden über gröberen Gittern platziert und nicht aus deren Verarbeitung ausgeschlossen. Somit existieren über der gesamten \emph{Fläche} feinerer Gitter auch die Knoten gröberer Gitter. \begin{figure}[h] \centering @@ -239,7 +243,17 @@ Vorteil dieser Herangehensweise ist es, dass feinere Teilgitter im Simulationsve \subsubsection{Multi-Domain Ansatz} -Kern des Multi-Domain Ansatzes ist es, außerhalb von etwaigen verfahrensbedingten Übergangsbereichen, jede Position in der Simulationsdomäne durch genau ein Teilgitter abzubilden. Konkret werden also bereits durch feinere Gitter abgedeckte Bereiche aus gröberen Teilgittern ausgespart. +Kern von Multi-Domain Ansätzen +\cite{Chen06} +\cite{DupuisChopard03} +\cite{EitelAmor13} +\cite{Fakhari16} +\cite{Filippova98} +\cite{Lagrava12} +\cite{Peng06} +%\cite{Rheinlaender05} % diffusiv +\cite{Rohde06} +ist es, außerhalb von etwaigen verfahrensbedingten Übergangsbereichen, jede Position der Simulationsdomäne durch genau ein Teilgitter abzubilden. Konkret werden also bereits durch feinere Gitter abgedeckte Bereiche aus gröberen Teilgittern ausgespart. \begin{figure}[h] \centering @@ -395,7 +409,7 @@ Zusammenfassend wird die Aufgabe der im kommenden Kapitel zu erarbeitenden Skali \subsection{Komponenten der Gitterkopplung}\label{kap:Komponenten} \subsubsection{Skalierung}\label{kap:Skalierung} -Während die Skalierung räumlicher Größen durch die Festlegung des Verfahrens auf Übergänge im Verhältnis \(1:2\) definiert ist, eröffnen sich für die zeitliche Skalierung zwei Möglichkeiten: Konvektive oder diffusive Skalierung. Unterschied der beiden Ansätze ist dabei das jeweilige Verhältnis zwischen räumlicher und zeitlicher Auflösung. +Während die Skalierung räumlicher Größen durch die Festlegung des Verfahrens auf Übergänge im Verhältnis \(1:2\) definiert ist, eröffnen sich für die zeitliche Skalierung zwei Möglichkeiten: Konvektive oder diffusive Skalierung \cite[Kap.~7.2.2.1]{Krueger17}. Unterschied der beiden Ansätze ist dabei das jeweilige Verhältnis zwischen räumlicher und zeitlicher Auflösung. \begin{Definition}[Konvektive Skalierung] Sei \(\delta t > 0\) die zeitliche und \(\delta x > 0\) die räumliche Diskretisierung. Dann gilt bei konvektiver Skalierung das Verhältnis: @@ -409,7 +423,7 @@ Sei \(\delta t > 0\) die zeitliche und \(\delta x > 0\) die räumliche Diskretis Es besteht hier also eine quadratische Proportionalität. Im Vergleich zu einer konvektiven Skalierung ist die zeitliche Auflösung somit um eine Ordnung feiner. \end{Definition} -Es ist klar zu erkennen, dass diffusive Skalierung einen deutlich größeren numerischen Aufwand gegenüber der konvektiven Skalierung nach sich zieht. Vorteil der bei diffusiver Skalierung erhöhten Schrittanzahl pro Zeiteinheit sind kleinere Fehler. +Es ist klar zu erkennen, dass diffusive Skalierung einen deutlich größeren numerischen Aufwand gegenüber der konvektiven Skalierung nach sich zieht. Vorteil der bei diffusiver Skalierung erhöhten Schrittanzahl pro Zeiteinheit sind kleinere Fehler. Darüber hinaus ist nach \cite[S.~276]{Krueger17} die Konvergenz gegen die inkompressiblen Navier-Stokes Gleichungen nur bei diffusiver Skalierung gegeben. Für die Autoren des hier erörterten Gitterverfeinerungsverfahrens überwog dabei das Argument der numerischen Effizienz, weshalb auch wir hier nun die konvektive Skalierung betrachten wollen. Die Austauschbarkeit des Skalierungsverfahrens sollte jedoch bei der Implementierung in OpenLB beachtet werden, da dieser Aspekt eine weitere prinzipiell flexible Komponente des Verfahrens darstellt. diff --git a/main.tex b/main.tex index 602c508..dca44c8 100644 --- a/main.tex +++ b/main.tex @@ -110,7 +110,7 @@ Trotz anhaltendem Wachstum der verfügbaren Rechenleistung können viele reale S OpenLB ist eine in C++ geschriebene freie Bibliothek zur Implementierung von LBM basierenden Strömungssimulationen. Aktuell bietet OpenLB noch keine Unterstützung für Gitterverfeinerung. -Ziel dieser Arbeit ist es, aus der Behebung dieser Einschränkung heraus ein Gitterverfeinerungsverfahren nachzuvollziehen und dessen Stärken und Schwächen zu evaluieren. +Ziel dieser Arbeit ist es, aus der Behebung dieser Einschränkung heraus einen Gitterverfeinerungsansatz theoretisch nachzuvollziehen und dessen Einsetzbarkeit praktisch zu evaluieren. Zu diesen Zweck werden anhand einer zweidimensionalen Lattice Boltzmann Methode verschiedene Ansätze zur Verfeinerung von Gittern diskutiert. Darauf aufbauend wird ein konkretes Verfahren detailliert ausformuliert und im Rahmen der Entwicklung eines generischen Gitterverfeinerungsframeworks umgesetzt. Der übergeordneten Fragestellung nach dem tatsächlichen Nutzen und möglichen Problemen von gitterverfeinerten Lattice Boltzmann Methoden wird durch die Evaluation von Anwendungsbeispielen Rechnung getragen. In diesem Kontext findet weiterhin eine Besprechung der formal begründeten anwendungsbezogenen Auswahl von zu verfeinernden Gebieten statt. \end{abstract} \newpage diff --git a/quellen.bib b/quellen.bib index 58e7eba..6f4f4f8 100644 --- a/quellen.bib +++ b/quellen.bib @@ -1,10 +1,10 @@ @book{Krueger17, - title = {The Lattice Boltzmann Method: Principles and Practice}, - author = {T. {Krüger} and H. {Kusumaatmaja} and A. {Kuzmin} and O. {Shardt} and G. {Silva} and E.M. {Viggen}}, + title = {The Lattice Boltzmann Method: Principles and Practice}, + author = {T. {Krüger} and H. {Kusumaatmaja} and A. {Kuzmin} and O. {Shardt} and G. {Silva} and E.M. {Viggen}}, publisher = {Springer International Publishing}, - isbn = {978-3-319-44647-9}, - year = {2017}, - series = {Graduate Texts in Physics}, + isbn = {978-3-319-44647-9}, + year = {2017}, + series = {Graduate Texts in Physics}, }, @book{Haenel04, title = {Molekulare Gasdynamik}, @@ -49,8 +49,6 @@ volume = 48, pages = {387-411}, doi = {10.1016/0021-9991(82)90058-4}, - adsurl = {http://adsabs.harvard.edu/abs/1982JCoPh..48..387G}, - adsnote = {Provided by the SAO/NASA Astrophysics Data System} }, @article{DupuisChopard03, author = {Dupuis, A. and Chopard, B.}, @@ -64,17 +62,6 @@ publisher = {American Physical Society}, doi = {10.1103/PhysRevE.67.066707}, }, -@article{Filippova98, - author = {O. {Filippova} and D. {Hänel}}, - title = {Grid Refinement for Lattice-BGK Models}, - journal = {Journal of Computational Physics}, - volume = 147, - number = 1, - pages = {219-228}, - year = 1998, - issn = {0021-9991}, - doi = {10.1006/jcph.1998.6089}, -}, @book{AmannEscherI, title = {Analysis I}, author = {H. Amann and J. Escher}, @@ -173,5 +160,134 @@ and Brian Barrett and Andrew Lumsdaine and Ralph H. Castain and David J. Daniel acmid = {615542}, publisher = {IEEE Computer Society Press}, address = {Los Alamitos, CA, USA}, +}, +@article{Filippova98, + author = {O. {Filippova} and D. {Hänel}}, + title = {Grid Refinement for Lattice-BGK Models}, + journal = {Journal of Computational Physics}, + volume = 147, + number = 1, + pages = {219-228}, + year = 1998, + issn = {0021-9991}, + doi = {10.1006/jcph.1998.6089}, +}, +@article{Chen06, + title = {Grid refinement in lattice Boltzmann methods based on volumetric formulation}, + journal = {Physica A: Statistical Mechanics and its Applications}, + volume = {362}, + number = {1}, + pages = {158 - 167}, + year = {2006}, + note = {Discrete Simulation of Fluid Dynamics}, + issn = {0378-4371}, + doi = {10.1016/j.physa.2005.09.036}, + author = {H. Chen and O. Filippova and J. Hoch and K. Molvig and R. Shock and C. Teixeira and R. Zhang}, +}, +@article{Fakhari16, + title = {A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows}, + journal = {Journal of Computational Physics}, + volume = {315}, + pages = {434 - 457}, + year = {2016}, + issn = {0021-9991}, + doi = {10.1016/j.jcp.2016.03.058}, + author = {A. Fakhari and M. Geier and T. Lee}, +}, +@phdthesis{Freudinger09, + author = {Freudiger, S.}, + year = {2009}, + month = {09}, + title = {Entwicklung eines parallelen, adaptiven, komponentenbasierten Strömungskerns für hierarchische Gitter auf Basis des Lattice-Boltzmann-Verfahrens}, + note = {Technische Universität Carolo-Wilhelmina zu Braunschweig} +}, +@article{EitelAmor13, + title = {A lattice-Boltzmann method with hierarchically refined meshes}, + journal = {Computers \& Fluids}, + volume = {75}, + pages = {127 - 139}, + year = {2013}, + issn = {0045-7930}, + doi = {10.1016/j.compfluid.2013.01.013}, + author = {G. Eitel-Amor and M. Meinke and W. Schröder}, +}, +@article{DupuisChopard03, + title = {Theory and applications of an alternative lattice Boltzmann grid refinement algorithm}, + author = {Dupuis, A. and Chopard, B.}, + journal = {Phys. Rev. E}, + volume = {67}, + issue = {6}, + year = {2003}, + month = {Jun}, + publisher = {American Physical Society}, + doi = {10.1103/PhysRevE.67.066707}, +}, +@article{Peng06, + author = {Peng, Y. and Shu, C. and Chew, Y. T. and Niu, X. D. and Lu, X. Y.}, + title = {Application of Multi-block Approach in the Immersed Boundary-lattice Boltzmann Method for Viscous Fluid Flows}, + journal = {Journal of Computational Physics}, + volume = {218}, + number = {2}, + month = nov, + year = {2006}, + issn = {0021-9991}, + pages = {460--478}, + doi = {10.1016/j.jcp.2006.02.017}, + publisher = {Academic Press Professional, Inc.}, +}, +@article{Morrison18, + title = {Sediment transport in turbulent flows with the lattice Boltzmann method}, + journal = {Computers \& Fluids}, + volume = {172}, + pages = {340 - 351}, + year = {2018}, + issn = {0045-7930}, + doi = {10.1016/j.compfluid.2018.04.015}, + author = {H.E. Morrison and A. Leder}, +}, +@article{Rohde06, + author = {Rohde, M. and Kandhai, D. and Derksen, J.J. and Van den Akker, H.}, + year = {2006}, + month = {06}, + pages = {439 - 468}, + title = {A generic, mass conservative local grid refinement technique for lattice-Boltzmann schemes}, + volume = {51}, + journal = {International Journal for Numerical Methods in Fluids}, + doi = {10.1002/fld.1140} +}, +@article{Rheinlaender05, + author = {Rheinl{\"a}nder, M.}, + title = {A Consistent Grid Coupling Method for Lattice-Boltzmann Schemes}, + journal = {Journal of Statistical Physics}, + year = {2005}, + month = {Oct}, + volume = {121}, + number = {1}, + pages = {49--74}, + issn = {1572-9613}, + doi = {10.1007/s10955-005-8412-0}, +}, +@article{Lin00, % multi grid + title = {Lattice Boltzmann method on composite grids}, + author = {Lin, Ching-Long and Lai, Yong G.}, + journal = {Phys. Rev. E}, + volume = {62}, + issue = {2}, + pages = {2219--2225}, + year = {2000}, + month = {Aug}, + publisher = {American Physical Society}, + doi = {10.1103/PhysRevE.62.2219}, +}, +@article{Toelke02, % multi grid + author = {T{\"o}lke, J. and Krafczyk, M. and Rank, E.}, + title = {A Multigrid-Solver for the Discrete Boltzmann Equation}, + journal = {Journal of Statistical Physics}, + year = {2002}, + month = {Apr}, + volume = {107}, + number = {1}, + pages = {573--591}, + issn = {1572-9613}, + doi = "10.1023/A:1014551813787", } - -- cgit v1.2.3