Developing of computer simulations of different engineering problems with adaptive three-dimensional computational meshes is nowadays a challenging research topic. This approach allows for simulating efficiently complex physical phenomena even on stationary computers. One of the adaptive methods of solving complex engineering problems is the adaptive finite element method. The computational domain is described by means of computational mesh consisting of so-called finite elements. The simulated physical phenomena are approximated with a linear combination of basis functions spread over the finite elements. In order to increase the accuracy of solution some of the fi nite elements are refined for example, partitioned into smaller elements. The computational problem is solved in each step of adaptation using a solver algorithm, for example, the multifrontal solver. This monograph presents a graph grammar based model describing the process of mesh generation, adaptation as well as the generation and factorization of the computational problem defined over the mesh.
Performing the adaptive three-dimensional computer simulations by means of graph grammar allows for speeding up the computations by an auto- matic construction of optimal order of numerical computations.
The graph grammar system presented in the monograph stores the history of the mesh generation and refinement. Based on this knowledge it is possible to define the optimal order of numerical computation by the multifrontal solver. This, in turn, allows to obtain faster factorizations of the systems of linear equations generated during the simulations. In case of some adaptive meshes, the resulting computational cost of the multifrontal solver is linear. The results of computations performed on three-dimensional meshes modeled by graph grammars can also be stored in a distributed manner in nodes of the graph modeling the computational mesh. This approach allows to reuse some results of computations from previous adaptive iteration in the case when only part of a mesh was refined in the next iteration. The graph grammar system used for modeling mesh generation and mesh adaptation guarantees the correctness of generated computational mesh the computational mesh will be connected, without incorrect finite elements (e.g. twisted elements or elements with duplicated faces) resulting in numerical problems during factorization of the system of linear equations generated over the computational mesh during the simulation.
Symulacje komputerowe przeprowadzane w celu rozwiązywania problemów inżynieryjnych za pomocą adaptowanych trójwymiarowych siatek obliczeniowych stanowią ważny kierunek badań naukowych, ze względu na możliwość symulowania złożonych zjawisk fizycznych. Jedną z adaptacyjnych metod rozwiązywania trudnych problemów inżynieryjnych jest adaptacyjna metoda elementów skończonych. Dziedzina obliczeniowa opisywana jest za pomocą siatki obliczeniowej, złożonej z tak zwanych elementów skończonych. Na elementach skończonych rozpina się funkcje bazowe służące do przybliżania modelowanego zjawiska. Wybrane elementy skończone są następnie adaptowane, np. łamane na mniejsze, w celu uzyskania dokładniejszego rozwiązania. W każdym kroku adaptacji rozwiązuje się problem obliczeniowy korzystając z solwera, najczęściej wielofrontalnego, gwarantującego dokładne rozwiązanie numeryczne. Obliczenia za pomocą adaptacyjnej metody elementów skończonych są więc procesem iteracyjnym, rozpoczynającym się od wygenerowanej początkowej siatki obliczeniowej, która następnie jest poprawiana w sposób iteracyjny, poprzez wykonywanie adaptacji w wybranych miejscach siatki. W pracy przedstawiam proces iteracyjnych obliczeń przeprowadzanych za pomocą adaptacyjnej metody elementów skończonych poprzez sekwencję produkcji gramatyki grafowej. Przeprowadzanie adaptacyjnych trójwymiarowych symulacji komputerowych za pomocą produkcji gramatyk grafowych umożliwia znaczne przyspieszenie obliczeń poprzez automatyczne określenie optymalnej kolejności przeprowadzania obliczeń numerycznych. Przedstawiony w pracy system gramatyk grafowych przechowuje historię generacji oraz adaptacji siatki obliczeniowej. Dzięki temu możliwe jest przyspieszanie obliczeń wykonywanych przez solwer wielofrontalny poprzez automatyczne określenie optymalnej kolejności przeprowadzania obliczeń numerycznych przez solwer wielofrontalny, w szczególnych przypadkach dającej liniowy koszt faktoryzacji. Wyniki obliczeń przeprowadzanych na trójwymiarowych siatkach obliczeniowych modelowanych za pomocą gramatyk grafowych można przechowywać w sposób rozproszony w węzłach grafu modelującego siatkę obliczeniową. Dzięki temu możliwe jest ponowne wykorzystanie części poprzednio przeprowadzonych obliczeń podczas wykonywania nowego obliczenia, w sytuacji, gdy jedynie część siatki obliczeniowej została adaptowana. Takie podejście umożliwia znaczną redukcję kosztu obliczeniowego przeprowadzanej symulacji. System gramatyk grafowych zastosowany do generacji oraz adaptacji trójwymiarowej siatki obliczeniowej gwarantuje również poprawność wygenerowanej siatki obliczeniowej - wygenerowana siatka będzie spójna, bez niepoprawnych elementów skończonych (np. przekręconych wzdłuż osi lub posiadających zduplikowane ściany) skutkujących problemami numerycznymi podczas przeprowadzania symulacji komputerowych.
Wydawnictwa nie prowadzą sprzedaży książek z serii "Rozprawy Monografie". Zainteresowanych prosimy o kontakt z ich autorami.
- Spis treści
-
Summary 9
Streszczenie 10
Index of important symbols 11
1 Introduction 13
2 State-of-the-art 16
2.1 Graph grammars 16
2.2 Graph grammar models of adaptive computations 17
2.3 Frontal and multifrontal solver algorithms 18
2.4 Orderings and elimination trees 18
2.5 Graph grammar models of multifrontal solver algorithms 19
3 Preliminaries 20
3.1 Basic definitions towards composition programmable graph grammar systems 20
3.1.1 Composition graphs 20
3.1.2 Composition graph grammar 22
3.2 Three-dimensional adaptive finite element methods 29
4 Modeling three-dimensional adaptive finite element method by graph grammar 37
4.1 Graph grammar productions for generation of the initial mesh with tetrahedral elements 37
4.2 Graph grammar productions interfacing with the solver algorithm with tetrahedral elements 47
4.3 Graph grammar productions selecting optimal refinements with tetrahedral elements 51
4.4 Graph grammar productions performing optimal refinements with tetrahedral elements 52
4.5 Formal definitions of graph grammar system modeling adaptive finite element method with tetrahedral elements 64
4.6 Extension of the graph grammar model for three-dimensional hexahedral finite elements 75
5 The utilization of graph grammar system for adaptive finite element method computations 80
5.1 Guarantee of the correctness of the mesh 80
5.2 Guarantee on the correctness of the adaptation process 83
5.3 Possibility of construction of new orderings 83
5.4 Reutilization of LU factorizations over unrefined parts of the mesh 102
6 Application of the hierarchical chromosome-based genetic algorithm for finding of quasi-optimal initial grids for three-dimensional adaptive finite element method 117
6.1 The hp-adaptive Finite Element Method Algorithm 117
6.2 Application of hierarchical chromosome-based genetic algorithm for initial mesh selection 120
6.3 Numerical verification of the HCBGA algorithm 125
7 Numerical results 129
7.1 The case study for the graph grammar modeling three-dimensional adaptive projection problem 129
7.2 The case study for the graph grammar modeling three-dimensional hp-adaptive finite element method with hexahedral elements 131
7.3 The case study for the graph grammar modeling three-dimensional h-adaptive finite element method with tetrahedral, prism and pyramid elements 140
7.3.1 Graph grammar productions expressing the mesh generation algorithm 140
7.3.2 Graph grammar productions expressing the mesh refinements algorithm 146
7.3.3 Graph grammar productions interfacing with the multifrontal direct solver 148
8 Conclusions and future work 152
8.1 Conclucions 152
8.2 Current and future work 154
Appendix A 159
Appendix B 165
Bibliography 169