From 0b7095a37de5932ad1795b87a9dcadc78e133c96 Mon Sep 17 00:00:00 2001 From: hartmann Date: Wed, 12 Sep 2001 14:11:20 +0000 Subject: [PATCH] Small changes and small fixes. git-svn-id: https://svn.dealii.org/trunk@4978 0785d39b-7218-0410-832d-ea1e28bc413d --- .../step-10.data/intro.html | 25 ++++++++++--------- .../step-10.data/results.html | 21 ++++++++-------- 2 files changed, 23 insertions(+), 23 deletions(-) diff --git a/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/intro.html b/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/intro.html index a460d933ea..7cc6081ef1 100644 --- a/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/intro.html +++ b/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/intro.html @@ -6,7 +6,7 @@ This is a rather short example which only shows some aspects of using higher order mappings. By mapping we mean the transformation between the unit cell (i.e. the unit line, square, or cube) to the cells in real space. In all the previous examples, we have implicitly -used linear or d-linear mappings; you will not have notices this at +used linear or d-linear mappings; you will not have noticed this at all, since this is what happens if you do not do anything special. However, if your domain has curved boundaries, there are cases where the piecewise linear approximation of the boundary @@ -14,27 +14,28 @@ cases where the piecewise linear approximation of the boundary your computational domain is an approximation to the real domain using curved boundaries as well. If the boundary approximation uses piecewise quadratic parabolas to approximate the true boundary, then -we say that this is a quadratic or Q2 approximation. If we -use piecewise graphs of cubic polynomials, then this is a Q3 +we say that this is a quadratic or Q2 approximation. If we +use piecewise graphs of cubic polynomials, then this is a Q3 approximation, and so on.

For some differential equations, it is known that piecewise linear -approximations of the boundary, i.e. Q1 mappings, are not +approximations of the boundary, i.e. Q1 mappings, are not sufficient if the boundary of the domain is curved. Examples are the biharmonic equation using C1 elements, or the Euler -equation. In these cases, it is necessary to compute the integrals -using a higher order mapping. The reason, of course, is that if we do -not use a higher order mapping, the order of approximation of the -boundary dominates the order of convergence of the entire numerical -scheme, irrespective of the order of convergence of the discretization -in the interior of the domain. +equation on domains with curved reflective boundaries. In these cases, +it is necessary to compute the integrals using a higher order +mapping. The reason, of course, is that if we do not use a higher +order mapping, the order of approximation of the boundary dominates +the order of convergence of the entire numerical scheme, irrespective +of the order of convergence of the discretization in the interior of +the domain.

Rather than demonstrating the use of higher order mappings with one of -these more complicated mappings, we do only a brief computation: +these more complicated examples, we do only a brief computation: calculating the value of pi=3.141592653589793238462643... by two different methods.

@@ -45,7 +46,7 @@ unit radius and integrates the unit function over it. Of course, if the domain were the exact unit circle, then the area would be pi, but since we only use an approximation by piecewise polynomial segments, the value of the area is not exactly pi. However, it is known that as -we refine the triangulation, a Qp mapping approximates the boundary +we refine the triangulation, a Qp mapping approximates the boundary with an order hp+1, where h is the mesh width. We will check the values of the computed area of the circle and their convergence towards pi under mesh refinement for different diff --git a/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/results.html b/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/results.html index dd5ae80078..1693e85ea7 100644 --- a/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/results.html +++ b/deal.II/doc/tutorial/chapter-2.step-by-step/step-10.data/results.html @@ -72,10 +72,9 @@ with dashed lines the exact circle: -The quadratic mapping obviously quite well approximates the -boundary, while for the cubic mapping the difference between -approximated domain and true one is hardly visible already for the -coarse grid. +Obviously the quadratic mapping approximates the boundary quite well, +while for the cubic mapping the difference between approximated domain +and true one is hardly visible already for the coarse grid.

@@ -164,25 +163,25 @@ cells eval.pi error

One of the immediate observations from the output is that in all cases the values converge quickly to the true value of -pi=3.141592653589793238462643. Note that for for Q4 mapping, the last +pi=3.141592653589793238462643. Note that for the Q4 mapping, the last number is correct to 13 digits in both computations, which is already -quite a lot. However, also note that for the Q1 mapping, even on the +quite a lot. However, also note that for the Q1 mapping, even on the finest grid the accuracy is significantly worse than on the coarse -grid for a Q4 mapping! +grid for a Q4 mapping!

The last column of the output shows the convergence order, in powers of the mesh width h. In the introduction, we had stated that -the convergence order for a Qp mapping should be -hp+1. However, in the example shown, the Q2 and Q4 +the convergence order for a Qp mapping should be +hp+1. However, in the example shown, the Q2 and Q4 mappings show a convergence order of hp+2! This at first surprising fact is readily explained by the particular boundary -we have here. In fact, the circle is described by the function +we have chosen in this example. In fact, the circle is described by the function sqrt(1-x2), which has the series expansion 1-x2/2-x4/8-x6/16+... around x=0. Thus, for the quadratic mapping where the truncation error of the quadratic approximation should be cubic, there is no such term but only a quartic one, which raises the convergence -order to 4, instead of 3. The same happens for the Q4 mapping. +order to 4, instead of 3. The same happens for the Q4 mapping.

-- 2.39.5