Updated Gallery (markdown)
[dealii.wiki.git] / Gallery.md
1 # A gallery of nice pictures generated by deal.II
2
3 This page is a collection of images put here for mainly aesthetic reasons, and a little bit in order to
4 show what can be done with deal.II. Most of the pictures also contain a brief summary of what they show, in order to give an idea of the kind of problem they are meant to solve.
5
6 Feel free to add your own pictures here, but please add the name of the person who did the simulation -- others may be so excited by it that they would want to contact you!
7
8 ## Supersonic flow around a cylinder
9
10 <a href="https://www.dealii.org/images/wiki/gallery/step-69.2nd-order.2d-2.jpg">
11 <img height="400px" src="https://www.dealii.org/images/wiki/gallery/step-69.2nd-order.2d-2.jpg" />
12 </a>
13
14 Finite element simulation with an invariant-domain preserving convex limiting scheme (9.5M gridoints). The corresponding first-order method is discussed in <a href="https://www.dealii.org/developer/doxygen/deal.II/step_69.html">step-69</a>. There are also videos available of the 2D computation and a corresponding 3D simulation with 1.8 billion gridpoints (click on the pictures below):
15
16 <a href="https://tamiko.43-1.org/developer/mach-3-cylinder-9.5M.mp4">
17 <img width="350px" src="https://www.dealii.org/images/wiki/gallery/step-69.2nd-order.2d.jpg" />
18 </a>
19 <a href="https://tamiko.43-1.org/developer/mach-3-cylinder-3D-1817M.mp4">
20 <img width="350px" src="https://www.dealii.org/images/wiki/gallery/step-69.2nd-order.3d-2.jpg" />
21 </a>
22
23 (Author: Matthias Maier, 10 April 2020)
24
25
26 ## Flow through the blades of a turbine
27
28 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/turbine.png" />
29
30 High-Reynolds number
31 computations lead to instationary solutions that are, however, often
32 enormously important in practice. This picture is by Manuel Quezada de
33 Luna (Texas A&M University) and was taken from a simulation of flow
34 between the blades of a compressor or turbine.
35
36
37 ## Complex geometries
38
39 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/yuhan-1.png" />
40 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/yuhan-2.png" />
41
42 These pictures are from a semester project by Yuhan Zhou at Texas A&M
43 University. The goal was to compute parasitic conducivities and impedances
44 from multilayer chip layouts. Other than generating data, this project also
45 produced a number of pretty pictures!
46
47
48
49 ## Thermally driven convection
50
51 <table>
52   <tr><td><img width="400px" src="https://www.dealii.org/images/wiki/gallery/Step-22.candles-1.png" /></td></tr>
53 </table>
54
55 The image shows results obtained with a variant of step-22 and step-31, tutorial programs that compute thermally driven convection with Stokes flow. The image shows three non-equidistantly spaced heat sources at the bottom and the flow field that emanates from the rising plumes of hot air.
56
57 The flow pattern becomes unstable at various times switching from one configuration to another. The full dynamics only really become visible when watching these images as a movie.
58 Several movies of various configurations are posted on Wolfgang Bangerth's webpage [https://www.math.colostate.edu/~bangerth/pictures.html#convection](https://dealii.org/developer/doxygen/deal.II/step_22.html]).
59
60
61 (Author: Wolfgang Bangerth, December 2007)
62
63
64
65
66 ## Singular and Hypersingular Source Terms
67
68 <table>
69   <tr><td><img width="400px" src="https://www.dealii.org/images/wiki/gallery/Single_layer.jpeg" /></td><td><img width="400px" src="https://www.dealii.org/images/wiki/gallery/Double_layer.jpeg" /></td></tr>
70 </table>
71
72 These images show a Laplacian with singular and hypersingular right hand sides, supported on a
73 co-dimension one circle.
74
75 (Author Luca Heltai, 16 May 2007)
76
77
78 ## Topology Optimization
79
80 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/Topopt1.png" />
81
82 These images result from a constrained topology optimization problem related to the minimization of steady state thermal compliance on a 2D square plate. The optimization is performed on a hierarchically adapted grid:
83 After completing optimization on a coarse grid level, the grid is adapted based on topology indicator derivatives and the optimization loop is restarted with interpolated topology indicator field on the (adaptively) refined grid as initial guess. The figure shows the global optimized topology and a selected zoomed part.
84
85 (Author: Rohallah Tavakoli, 9 February 2007)
86
87
88 ## Multiphase flow
89
90 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Random2d.png" />
91
92 This image shows the saturation at one time step of a simulation of how a mixture of two fluids moves through a random medium. The full description as well as movies of this problem and a 3d simulation can be found as step-21 in the tutorial of deal.II.
93
94 (Author: Wolfgang Bangerth, 2006)
95
96
97 ## Phase change/crystal growth
98
99 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-dendrite1.jpg" />
100 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-dendrite2.jpg" />
101 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-dendrite3.jpg" />
102
103 Dendritic growth of crystals.
104
105 (Author: Denis Danilov)
106
107
108 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/Gallery-Grain_growth.png" />
109
110 Phase field modeling of the normal grain growth.
111
112 -Slawa, 07:49, 6 September 2010 (UTC)
113
114
115
116
117 ## Incompressible Navier-Stokes
118
119 <img width="250px" src="https://www.dealii.org/images/wiki/gallery/Gallery-block100.jpg" />
120
121
122 <img width="250px" src="https://www.dealii.org/images/wiki/gallery/Gallery-backward_facing_step_01.jpg" />
123
124 The pressure of an incompressible flow around a backward facing step is visualized. The Reynolds number is 50. The question is, if the positive singularity (the left one, on the inflow side) is a real physical phenomenon. Indeed, it is numerically stable. (Andre Große-Wöhrmann, 2005)
125
126
127 ## Convection in the earth core
128
129 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-convection.png" />
130
131 Convection in itself is a hard problem (see the entry further up this page on thermally driven convection]](http://www.math.tamu.edu/~bangerth/pictures.html#optical].), but if it is between moving surfaces or with nonvertical gravity, it is even more challenging. This is what makes simulations of the earth's interior, both of the liquid rock earth mantle as well of the liquid metal outer core, so complicated. The picture to the right shows streamlines of a simulation of a convecting fluid in a radial gravity field between a hotter inner sphere and a cooler outer sphere, both of which are rotating and dragging the fluid along. Under such conditions, very complicated flow patterns develop after a while - among which are those that are responsible for the earth's magnetic field.
132
133 (Author: Andre Grosse-Woehrmann, 2005)
134
135
136 ## Solitons
137
138 <img width="180px" src="https://www.dealii.org/images/wiki/gallery/Gallery-step25-2.png" />
139 <img width="180px" src="https://www.dealii.org/images/wiki/gallery/Gallery-step25-3.png" />
140
141 The step-25 tutorial program demonstrates the solution of the nonlinear, wave equation-type sine-Gordon equation, from which the pictures here are taken. Full movies of these solutions can be found in the results section of that program as well.
142
143 (Author: Ivan Christov, 2006)
144
145 ## Transport
146
147 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-transport_01.jpg" />
148
149 Modelling transport via the wavelike approach, numerical waves occur.
150
151 (Author: Andre Große-Wöhrmann, 2005)
152
153
154
155 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-kpp.png" />
156
157 Modeling nonlinear transport equations is slightly more complicated since they can develop shocks even if the initial conditions are smooth. This isn't the case for the initial conditions that produced this solution of the KPP equations, but the difficulty of capturing the solution's features is still apparent.
158
159 (Author: Orhan Mamedov, Vladimir Tomov, Abner Salgado, as part of a student project, 2011)
160
161
162 ## Fictitious Domain Method
163
164 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-fdm_stokes_01.jpg" />
165
166 This is the boundary supported fictitious domain method applied for the problem of Stokes.`You can see the pressure and the velocity field. The velocity is enforced on the circle via a weak condition. The key feature of this method is that this curve may be independent of the discretization mesh.
167
168 (Author: Andre Große-Wöhrmann, 2005)
169
170 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-fdm_potential_01.jpg" />
171
172 This is the boundary supported fictitious domain method applied to the potential equation. Again the curve may be independent of the mesh. Nevertheless the mesh is refined in the vicinity of the curve.
173
174 (Author: Andre Große-Wöhrmann, 2005)
175
176
177 ## Plastic and quasistatic deformation
178
179 <img width="400px" src="https://www.dealii.org/images/wiki/step-42.displacement.png" />
180
181 Deformation can be described in many ways. step-42 is a tutorial program
182 (written by Joerg Frohne, 2013) that deals with plasticity and the picture above
183 shows the displacement by pressing a printing letter in the shape of the
184 Chinese character for "force" into a metal block. (The displacement at every
185 node is exaggerated by a factor of 100.)
186
187 <img width="300px" src="https://www.dealii.org/images/wiki/gallery/Solution-0010.0000.png" />
188 &nbsp;&nbsp;
189 <img width="300px" src="https://www.dealii.org/images/wiki/gallery/Step-18-contours.png" />
190
191 On the other hand, the images above correspond to the step-18 tutorial program. It shows the gradual deformation of a cylinder under pressure from above, until it essentially fails. The color coding indicates
192 stress levels in the material. Because it looks nice, at the very bottom, there is also an isosurface view of the x-displacement early on in the simulation. More details about this simulation can be found in the step-18
193 tutorial.
194
195 (Authors: Wolfgang Bangerth, 2005; Joerg Frohne, 2013)
196
197
198 ## Complicated domains
199
200 <img width="300px" src="https://www.dealii.org/images/wiki/gallery/Left_lung.png" />
201
202 Here is an image of a mesh of a rather complicated domain. It shows a mesh created from CT data of the left lung 
203
204 (Author: Li Pan, 2006)
205
206
207 ## Crack propagation
208
209 <img width="600px" src="https://www.dealii.org/images/wiki/gallery/crack_propagation_cts-specimen_03.png" />
210
211 Quasistatic crack propagation in brittle materials can be simulated using the energy principle, where a crack grows in the direction of maximal energy release. Here, the direction of the crack is computed by an asymptotic formula for the change of energy in two dimensions using (classical) stress intensity factors (linear elasticity). The crack is elongated step-by-step and the domain is remeshed in each step.
212
213 The pictures show crack paths in an orthotropic material with two axes of elastic symmetry in a compact-tension specimen. Both pictures show the crack under the same Mode-II shear loading, but with different orientation of the initial crack to the planes of elastic symmetry.
214
215 The third picture shows a crack propagating in an isotropic base material with an anisotropic local inhomogeneity. This numerical experiment demonstrates how a local perturbation in the material can influence a crack path.
216
217
218
219 ## Fun pictures
220
221 <img width="150px" src="https://www.dealii.org/images/wiki/gallery/Gallery-rotation.jpg" />
222
223 When changing the numbering to something regular, some of the structures in DataOutRotation were forgotten. Instead of iso-surfaces symmetric to the center, we got this nice picture looking a bit like a complicated turbine. (Guido Kanschat, 17 January 2006)
224
225 <img width="400px" src="https://www.dealii.org/images/wiki/gallery/Gallery-PacmanBenchmark.png" />
226
227 Poisson equation with homogeneous Dirichlet boundary conditions solved using Rvachev R-functions method combined with finite elements. Mesh is nonconforming, boundary conditions are enforced exactly. Pacman domain is considered as a benchmark because it contains curved parts as well as reentrant corner :)
228
229 For those who are interested to learn a bit about theory of R-functions, [here is a link](http://onlinelibrary.wiley.com/doi/10.1002/pamm.200701055/abstract) to my old 2-pages conference paper. More detailed  and extensive info in English is available at [Spatial Automation Lab](http://sal-cnc.me.wisc.edu/). (Slawa, 6 September 2010)
230
231 <img width="200px" src="https://www.dealii.org/images/wiki/gallery/Gallery-Solidification-Tree.png" />
232
233 This had to be a simulation of directional solidification/dendritic growth. But something was wrong with the model and it produced such a nice tree in the end :) (Slawa, 6 September 2010)

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.