p.27 | (C) |
The reference Oden [1991] is missing on p. 343: Oden, J. T. (1991): Finite elements: An introduction. In "Handbook of numerical analysis. Volume II: Finite element methods (Part 1)." (Ciarlet, P. G. et al., eds.), pp. 3-15. North-Holland, Amsterdam The access to a previous version may be easier for some readers: Oden, J. T., Historical comments on finite elements. A history of scientific computing (Princeton, NJ, 1987), 152-166, ACM Press Hist. Ser., ACM, New York, 1990. |
p.58 | (E) |
The stiffness matrix associated to the stencil (4.9)
induces a quadratic form z'Az = Σ_{i,j} (z_{i}-z_{j})² where the sum runs over all pairs i,j that belong to neighbors on the horizontal or vertical lines of the grid in Fig. 9. |
p.61 | (E) |
5.1 Definition (3):
Shape regularity may be formulated as a condition on the angles of the
triangles in a triangulation. It follows from Oswald's re-examination
of an example by Babuška and Aziz in 1976
that convergence may fail if the triangulation is not shape regular.
Thus shape regularity or a similar condition is required.
P. Oswald, Divergence of FEM: Babuška-Aziz triangulations revisited. Applications of Mathematics 60 (2015), 473-484. |
p.72+9 | (C) | ... It is concealed in the DKT element VI.5.3. |
p.74 | (E) | Due to Problem 5.14 a piecewise polynomial vector-valued function belongs to H(div), if its normal components are continuous at interelement boundaries. Similarly due to Problm VI.4.8 it belongs to H(rot), if its tangential components are continuous. The reader is warned for a wrong conclusion on the Sobolev spaces. The relatinn H(div) ∩ H(rot) = H^{1} is not always true. |
p.103 | (C) | Clearly, we get the same stiffness matrix for the triangulations shown in Figs. 27a and 27b. |
p.109 | (E) |
The mesh-dependent norm || . ||_{h} is the broken
H^{1} semi-norm. Often a broken Sobolev norm is an appropriate mesh-dependent norm; c.f., (II.6.1). Jump terms as in the treatment of the Raviart-Thomas element in §5 are possibly added. |
p.111-3 | (C) | The remark "This corresponds with the practical observation that nonconforming elements are much more sensitive to near singularities i.e., to the appearance of large H² norms" is not true. Theorem III.5.6 says that the errors of the conforming P_{1} element and of the nonconforming P_{1} element are comparable independently of the regularity. |
p.116 | (C) |
Problem 1.11.
The difference of the bilinear forms a(.,.) and a_h(.,.) is not coercive.
Therefore a stronger assumption is required: |a(u,v)-a_{h}(u,v)| ≤ ε ||u||.||v|| Moreover assume that the two bilinear forms are coercive with the constant α. |
p.116 | (E) |
Problem 1.13.
There seems to be a natural nonconforming P_{2} element of
Crouzeix-Raviart type. Let 0 < α < ½.
The natural nodes on the unit triangle are z_{1}=(α,0), z_{2}=(1-α,0), z_{3}=(1-α,α), z_{4}=(α,1-α), z_{5}=(0,1-α), z_{6}=(0,α). Show that the interpolation by quadratic polynomials at these points, however, is not well posed, since p(z_{1}) - p(z_{2}) + p(z_{3}) - p(z_{4}) + p(z_{5}) - p(z_{6}) = 0. |
p.125+3 | (E) | The importance of Theorem 3.6 for the finite element theory was pointed out by Babuška [1971]; see also Babuška and Aziz [1972]. The theorem is also found under the name Banach-Nečas-Babuška-Theorem. It can be traced back to Nečas [1962], Nirenberg and possibly to Banach; cf. Babuška [1971]. |
p.148 | (E) |
The theorem of Prager and Synge is also called the
two energy principle, where the left hand side
of (5.5)_{v} is the complementary energy.
Aubin and Burchard pointed out that the hypercircle method
can be traced back to Friedrichs and Trefftz.
The right hand side is sometimes replaced by 2J(v)+2J^{c}(σ), where J and J^{c} denote the direct and the complementary energy. It is not only used for a posteriori error estimates, but also for a justification of plate models; cf. the addition to p. 323. K.O. Friedrichs (1929): Ein Verfahren der Variationsrechnung das Minimum eines Integrals als das Maximum eines anderen Ausdrucks darzustellen. Ges. Wiss. Göttingen, Nachrichten Math. Phys. Kl. 13-20. E. Trefftz (1928): Konvergenz und Fehlerschätzung beim Ritzschen Verfahren. Math. Ann. 100, 503-521. |
p.152 | (E) | Theorem 5.6.
An error estimate for the finite element solution
with the nonconforming P_{1} element can be included
in the comparison. It was improved by C. Carstensen and D. Petersheim. D. Braess (2009). Calcolo 46, 149-155. C. Carstensen and D. Petersheim (2011) |
p.152+17/18 | (C) |
An index h is misplaced. On the left-hand side always: σ_{h}, on the right-hand side always: σ |
p.160+1 | (C) | The inequality (6.12) is called Nečas inequality. |
p.160 (6.12) | (E) |
The inequality ||q||_{0} ≤ c ||grad q||_{-1} (and therefore the inf-sup condition for the Stokes problem) does not hold on arbitrary domains. Consider the domain with a cusp {(x,y); 0 < y < x², 0 < x < 1} and the function q(x,y) = x^{-2}. It follows from ∂q/∂x = -2 ∂/∂y (y x^{-3}), ∂q/∂y = 0 and ||y x^{-3}||_{0} < ∞ that grad q belongs to H^{-1}, but ||q||_{0} = ∞. |
p.177 (8.26) | (C) | c^{-1} h_{e}|| ... |
p.180 | (E) | Problem 8.5 An a posteriori error estimator is to be checked. Since the exact solution is not available, an (approximate) reference solution is computed by using finite elements of higher order. Is the real error overestimated or underestimated if merely the error with respect to the reference solution is evaluated? |
p.181 | (E) | The term c||f-f_{h}|| that arises from the data oscillation can even be added by the Pythagorean rule instead of the triangle inequality; see Ainsworth cited below. |
p.183 | (E) |
Theorem 9.2 and Algorithm 9.3 can be found with different notation
in Destuynder, P. and Métivet, B. (1999): Explicit error bounds in a conforming finite element method. Math. Comp. 68, 1379-1396 |
p.184 | (E) |
The a posteriori estimator in Theorem 9.4 provides not only
an estimator for the P_{1} element, but also for the
mixed method of Raviart-Thomas. It is efficient since the
error of the mixed method is not dominant; see Theorem 5.6.
The converse is required if an efficient estimator for the Raviart-Thomas element is wanted. For this purpose the construction via piecewise quadratic elements by Ainsworth proceeds in a quite different way. Ainsworth, M. (2008): A posteriori error estimation for lowest order Raviart Thomas mixed finite elements. SIAM J. Sci. Comput. 30, 189-204 |
p.185 | (E) |
Remark.
The Comparison Theorem III.5.6 for the conforming P_{1} element,
the nonconforming P_{1} element, and the Raviart-Thomas element
contains an a priori estimate. There is the question:
Why are tools from the theory of a posteriori estimates used for its proof?
The reason is that the theorem is true only modulo data oscillation,
and the latter has been introduced and understood in the framework
of a posteriori estimates.
There is another fact of a similar type. The proof of the lower bound (8.19) uses cubic (and quadratic) bubble functions. It follows that P_{4} elements yield a solution with an error that is smaller than the error for P_{1} elements multiplied by a factor smaller than 1, provided that we disregard terms arising from data oscillation. We cannot do it without this addition, since it is easy to construct a right-hand side of the elliptic equation such that the finite element solution with P_{4} elements is contained in the subset of P_{1} elements. Obviously, the distance to the P_{4} solution does not reflect the distance to the true solution in this case. The use of techniques from a posteriori estimates for the a priori analysis of plates has a different reason. |
p.185 | (E) |
Remark.
The Poincaré Lemma for the inversion of the divergence in 2-space
does not help here. To understand this determine the divergence of ½ ∫_{0}^{1} txf(tx)dt. |
p.190 | (E) |
For a first convergence proof of The Gauss-Seidel method see:
v. Mises, R. and Pollaczek-Geiringer, H. (1929): Praktische Verfahren der Gleichungsauflösung. ZAMM 9, 58-77, 152-164. |
p.216+2 | (C) | -1/a_{i} ≤ b_{i}, c_{i} < 0 |
p.263-5 | (C) |
||u_{k+1}-u||²
≤ ||u_{k}+(α_{k}²-1) v^ -u||² (v^ means hat v) |
p.293 | (C) | Usually the boundary integral in (3.1) is included with g replaced by -g, and all terms with g in this section carry the opposite sign. |
p.296 (3.14) | (C) | b(τ,v) = (div τ,v)_{0} |
p.299 | (E) |
The constant in Korn's inequality depends on the shape of the domain
if Neumann boundary conditions are given on a part of the domain.
An Example is the cantilever beam (Fig. 58) with 0 ≤ x_{1} ≤ l, 0 ≤ x_{2} ≤ t, and t/l very small. Consider the displacement v_{1} = -3x_{1}^{2}x_{2}, v_{2} = x_{1}^{3}. Elementary calculations show that |v|_{1} is large when compared with ||ε(v)|| since ε_{12}(v) = ε_{22}(v) = 0. Note that ||rot v|| is also large. |
p.299 | (E) | The inf-sup condition for the Stokes problem implies Korn's inequality. The counterexample of a domain with a cusp shows that there is no implication in the converse direction. |
p.301+10 | (C) | div σ = – f |
p.301 | (E) | The minimization of the stress energy is also denoted as principle of Castigliano. |
p.311 (3.54) | (C) | The denominator is 4 since the correct denominators in the formulas above are 12. |
p.313+5 | (C) | Equation label (3.61) is missing. |
p.313+13 | (C) | γ_{h}=t^{-2}Q_{h}(w'_{h}-θ_{h}) |
p.313 (3.62)+2^{ } | (C) |
The integrand of the second term is [Q_{h}(w_{h}'-θ_{h})]². An index h is missing in the third term. |
p.314-3 | (C) | functions v of the form (3.16) |
pp.315-325 | (C) | Identify (x,y,z) with (x_{1},x_{2},x_{3}). |
p.323 | (E) |
Hypothesis H2: Verifications of the plate models have been
performed by using the theorem of Prager and Synge (p. 148).
It was done by Morgenstern [1959] for the Kirchhoff plate
and by Braess, Sauter, and Schwab [2009] for both plate models.
While the studies above refer to the displacement model,
Alessandrini et al [1999] investigated mixed methods
for the Mindlin-Reissner plate.
Although not explicitly stated, the results show
that Hypothesis H2 makes the plates stiffer than they are.
An ansatz with a quadratic term
u_{3} = w(x,y) + z² w_{2}(x.y) is required in both models. The lack of the quadratic term was often partially compensated by shear correction factors. Shear correction factors are disregarded in the text. Morgenstern, D. (1959): Herleitung der Plattentheorie aus der dreidimensionalen Elastizitätstheorie. (German) Arch. Ration. Mech. Anal. 4, 145-152. Alessandrini, A.L., Arnold, D.N., Falk, R.S., and Madureira, A.L. (1999): Derivation and justification of plate models by variational methods. In "Plates and Shells, Quebec 1996", (M. Fortin, ed.), pp. 1-20. CRM Proceeding and Lecture Notes, vol. 21, American Mathematical Society, Providence, RI. Braess, D., Sauter, S. and Schwab, C. (2009): On the justification of plate models. J. Elasticity 103, 53-71 (2011) |
p.326+12 | (C) | The Babuska paradox holds for simply-supported Kirchhoff plates. |
p.329+12 | (C) | cf. Problem II.5.16. |
p.335+10 | (C) | ... as we did in going from (3.29) to (3.31) and ... |
p.336 | (E) |
For a more recent survey of plate elements see: Falk, R.S. (2008): Finite elements for Reissner-Mindlin plates. In "Mixed Finite Elements, Compatibility Conditions, and Applications", (D. Boffi and L. Gastaldi, eds.), pp. 195-232, Springer (2008) |
p.351+5 | (C) | Braess, D. and Schöberl, J. (2006): Equilibrated residual error estimator for edge elements. Math. Comp. 77, 651-672 (2008) |
p.355 | (C) | A missing reference to T. J. Oden is added to the corrections to page 27 above. |