From 978a8254d6adb4121d41abcbaebb0f1435eec506 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Sun, 24 Nov 2019 15:49:47 -0500 Subject: [PATCH] step-55: doc update, add citation --- doc/doxygen/references.bib | 25 +++++++++++++++++++++++++ examples/step-55/doc/intro.dox | 17 ++++++++--------- 2 files changed, 33 insertions(+), 9 deletions(-) diff --git a/doc/doxygen/references.bib b/doc/doxygen/references.bib index 50f8639f93..9fa65721b4 100644 --- a/doc/doxygen/references.bib +++ b/doc/doxygen/references.bib @@ -307,6 +307,31 @@ MRREVIEWER = {Jose Luis Gracia}, URL = {http://dx.doi.org/10.1137/070706616}, } +% ------------------------------------ +% Step 55 +% ------------------------------------ + +@Book{elman2005, + Title = {Finite Elements and Fast Iterative Solvers with Applications in Incompressible Fluid Dynamics}, + Author = {Elman, H. C. and Silvester, D. J. and Wathen, A. J.}, + Publisher = {Oxford University Press}, + Year = {2005}, + Address = {Oxford, New York}, + Series = {Numerical Mathematics and Scientific Computation}, + ISBN = {0-19-852868-X}, +} + +@inproceedings{kovasznay1948laminar, + title={Laminar flow behind a two-dimensional grid}, + author={Kovasznay, LIG}, + booktitle={Mathematical Proceedings of the Cambridge Philosophical Society}, + volume={44}, + number={1}, + pages={58--62}, + year={1948}, + organization={Cambridge University Press} +} + % ------------------------------------ % Step 63 % ------------------------------------ diff --git a/examples/step-55/doc/intro.dox b/examples/step-55/doc/intro.dox index 6000945c85..5f70c1433c 100644 --- a/examples/step-55/doc/intro.dox +++ b/examples/step-55/doc/intro.dox @@ -2,7 +2,7 @@ This program was contributed by Timo Heister. Special thanks to Sander Rhebergen for the inspiration to finally write this tutorial. - + This material is based upon work partially supported by National Science Foundation grant DMS1522191 and the Computational Infrastructure in Geodynamics initiative (CIG), through the National Science Foundation under @@ -38,7 +38,7 @@ The learning outcomes for this tutorial are: - You understand the concept of optimal preconditioners and are able to check this for a particular problem. - + - You are able to construct manufactured solutions using the free computer algreba system SymPy (https://sympy.org). @@ -79,7 +79,8 @@ done in O(N) time. A preconditioner application can also only be O(N) at best required to solve the linear system is independent of $h$ (and therefore N), the total cost of solving the system will be O(N). It is not possible to beat this complexity, because even looking at all the entries of the right-hand -side already takes O(N) time. +side already takes O(N) time. For more information see @cite elman2005, +Chapter 2.5 (Multigrid). The preconditioner described here is even simpler than the one described in step-22 and will typically require more iterations and consequently time to @@ -122,9 +123,7 @@ where $S=-BA^{-1} B^T$ is the Schur complement. With this choice of $P$, assuming that we handle $A^{-1}$ and $S^{-1}$ exactly (which is an "idealized" situation), the preconditioned linear system has three distinct eigenvalues independent of $h$ and is therefore "optimal". See -section 6.2.1 (especially p. 292) in "Finite Elements and Fast Iterative -Solvers: with Applications in Incompressible Fluid Dynamics" by Elman, -Silvester, and Wathen (Oxford University Press (UK), 2005). For comparison, +section 6.2.1 (especially p. 292) in @cite elman2005. For comparison, using the ideal version of the upper block-triangular preconditioner in step-22 (also used in step-56) would have all eigenvalues be equal to one. @@ -152,9 +151,9 @@ size of the coarsest matrix is not predictable.

The testcase

-We will construct a manufactured solution based on the Kovasznay problem -("Laminar flow behind a two-dimensional grid" by Kovasznay, 1948). Here -is an image of the solution colored by the x velocity including +We will construct a manufactured solution based on the classical Kovasznay problem, +see @cite kovasznay1948laminar. Here +is an image of the solution colored by the x velocity including streamlines of the velocity: -- 2.39.5