Two-Sided A Posteriori Error Bounds for Electro-Magneto Static Problems
This paper is concerned with the derivation of computable and guaranteed
upper and lower bounds of the difference between the exact and the approximate solution
of a boundary value problem for static Maxwell equations.
Our analysis is based upon purely functional argumentation
and does not attract specific properties of an approximation method.
Therefore, the estimates derived in the paper at hand are applicable
to any approximate solution that belongs to the corresponding energy space.
Such estimates (also called error majorants of the functional type)
have been derived earlier for elliptic problems [19, 20].
Key Words a posteriori error estimates of functional type, Maxwell’s boundary value problem, electro-magneto statics
AMS MSC-Classifications 65 N 15, 78 A 30
Dedicated to the anniversary of Prof. Nina Nikolaevna Uraltseva
1 Introduction and notation
The main goal of this paper is to derive guaranteed and computable upper and lower bounds of the difference between the exact solution of an electro-magneto static boundary value problem and any approximation from the corresponding energy space. We discuss the method with the paradigm of a prototypical electro-magneto static problem in a bounded domain. The generalized formulation is given by the integral identity (2.6). We show that (as in many other problems of mathematical physics) certain transformations of (2.6) lead to guaranteed and fully computable majorants and minorants of the approximation error. However, the case considered here has special features that make (at some points) the derivation procedure different from, e.g., that which has been earlier applied to other elliptic type problems. This happens because the corresponding differential operator has a nonzero kernel (which contains curl-free vector fields) and the set of trial functions in (2.6) is restricted to a rather special subspace. For these reasons, the derivation of the estimates is based on Helmholtz-Weyl decompositions of vector fields, orthogonal projections onto subspaces, and on a certain version of a Poincaré-Friedrich estimate for the differential operator . First, we show that the distance between the exact solution and the approximate solution (measured throughout the semi-norm generated by the operator ) is equal to some norm of the so-called residual functional (cf. (3.1)). If satisfies the boundary condition exactly, i.e. , then the latter functional vanishes if and only if coincides with . Lemma 10 shows that an error majorant can be expressed throughout a certain norm of (cf. (3.2)). However, in general, computing of this norm is hardly possible because it requires finding a supremum over an infinite number of vector fields.
Theorem 14 provides a computable form of the upper bound. The corresponding estimate (3.12) shows that the error majorant is the sum of five terms, which can be thought of as penalties for possible violations of the relations (2.1)-(2.4). It contains only known vector fields and global constants depending on geometrical properties of the domain. Moreover, it is easy to see that the upper bound vanishes if and only if coincides with the exact solution and a ’free variable’ encompassed in the estimate coincides with . Also, we show that the estimates derived are sharp in the sense that the estimates (3.13) and (3.14) have no irremovable gap between the left and right hand sides (Remark 16). Finally, in Section 4, we derive lower estimates of the difference between exact and approximate solutions. The corresponding result is presented by Theorem 21. This estimate is also computable, guaranteed and sharp provided that the approximation exactly satisfies the prescribed boundary condition.
Throughout this paper, we consider a bounded domain with Lipschitz continuous boundary and denote the corresponding outward unit normal vector by . and stand for the electric and magnetic vector fields, respectively, while and denote positive definite, symmetric matrices with measurable, bounded coefficients that describe properties of the media (dielectricity and permeability, respectively). For the sake of brevity, matrices (matrix-valued functions) with such properties are called ’admissible’. We note that the corresponding inverse matrices are admissible as well. In particular, there exists a constant , such that for a.e.
By we denote the usual scalar -Hilbert space of square integrable functions over and by the Hilbert space of real-valued -vector fields, i.e. . For the sake of simplicity we restrict our analysis to the case of real valued functions and vector fields. The generalization to complex valued spaces is straight forward.
Orthogonality and the orthogonal sum with respect to the scalar product of is denoted by and , respectively, i.e. if
where denotes Lebesgue’s measure. Moreover, by (respectively ) we indicate the orthogonality (respectively orthogonal sum) in terms of the weighted -scalar product generated by an admissible matrix .
Throughout the paper we will utilize the following functional spaces:
Analogously, we define the spaces associated with the operators and . Furthermore, we introduce the spaces (containing the so-called Dirichlet and Neumann fields)
Here and later on we write if . These are finite dimensional spaces, whose dimensions are denoted by and , respectively. In fact, these numbers are equal to the so-called Betti numbers of and depend only on topological properties of the domain (for a detailed presentation see ). A basis of shall be given by special vector fields .
Finally, we note that being equipped with the proper inner products all the above introduced functional spaces are Hilbert spaces.
The classical formulation of the electro-magneto static problem for a given vector field (driving force) and given , reads as follows: Find a magnetic field
and a corresponding electric field
such that in
In other words, the problem is to find vector fields and
where the homogeneous boundary conditions are to be understood in the weak sense.
This coupled problem is equivalent to an electro-magneto static Maxwell problem in second order form, which in classical terms reads as follows: Find an electric field in , such that belongs to and
holds in , i.e. find , such that and
Once has been found, the magnetic field is given by .
2 Variational formulation and solution theory
Henceforth, we consider (1.2)-(1.5) assuming that the boundary condition on may be inhomogeneous (physically, such a condition is motivated by the presence of electric currents on the boundary). Hence, we intend to discuss the following prototypical electro-magneto static Maxwell problem in second order form: Find an electric field in , such that
i.e., find in
satisfying (2.1) and (2.3). There are at least two methods to prove existence of the solution. One is based upon Helmholtz-Weyl decompositions (see, e.g. [9, 14, 15, 17, 10, 11]). The second method consists of introducing and studying a suitable generalized statement of the problem (2.1)-(2.4). In this paper, we use the second method because it provides a natural way of deriving error estimates. Both methods are based on Poincaré-Friedrich estimates, see Remark 8, and (if it is needed) exploit suitable extension operators for the boundary data. On this way, we also need a certain version of the Poincaré-Friedrich estimate, namely
Now, let be some vector field in satisfying the boundary condition (2.3) in the generalized sense, i.e., . The generalized solution
If then by the Cauchy-Scharz inequality the right hand side of (2.6) is a linear and continuous functional over . By (2.5) the left hand side of (2.6) is a strongly coercive bilinear form over . Thus, under these assumptions the problem (2.6) is uniquely solvable in by Lax-Milgram’s theorem.
First, we note some Helmholtz-Weyl decompositions of , i.e. decompositions into solenoidal and curl-free fields, which will be used frequently throughout our analysis.
can be decomposed as
where all closures are taken in and . Moreover,
Let us denote the -orthogonal projection onto in (2.8) by . Then we have for all
and for all
The latter line can be written in a more compact and precise way as
The following theorem states the main existence result.
The boundary data and its extension can be described in more detail. Since the papers [1, 3, 4] and the more general paper of Weck  we know that even for Lipschitz domains, where the non scalar trace business is a challenging task, there exist a bounded linear tangential trace operator and a corresponding bounded linear tangential extension operator (right inverse) mapping to special tangential vector fields on the boundary, i.e.
and vice verse. Here, denotes the surface . Using the Helmholtz-Weyl decomposition (2.8) we even get an improved extension operator. We have
Applied to smooth vector fields we have . Now, we may specify the boundary data and the extension as well as our variational formulation for : Find , such that
holds for all .
Henceforth, we assume that is given by a tangential trace of some vector field .
More general variants of the Poincaré-Friedrich estimate for vector fields (2.5) are known. For instance, we have
which holds for all . This estimate may be proved by an indirect argument using a ’Maxwell compact embedding property’ of , which holds true not only for Lipschitz domains, but also, if the homogeneous boundary condition is considered, for more irregular domains (cone properties), see . For inhomogeneous boundary conditions the Lipschitz assumption can not be weakened. Actually, it is just the continuity of the solution operator of the corresponding electro static boundary value problem, see [5, 6, 7].
3 Upper bounds for the deviation from the exact solution
Let be an approximation of . We assume that belongs to , which means that, in general, the boundary condition, the divergence-free condition, and the orthogonality to the Dirichlet fields might be violated, i.e. the approximation field may be such that
Moreover, for the subsequent analysis and then also for the numerical application, which is even more important, it is sufficient to assume just .
Our goal is to obtain upper bounds for the difference between and in terms of the weighted norm
First, we use (2.6) and get for all
where is a linear and continuous functional over as well as over . Furthermore, does not depend on the exact solution .
Obviously, vanishes if . Furthermore, if satisfies the boundary condition exactly, i.e. , then if and only if (or what is equivalent if and only if ). This holds by the following argument using the Helmholtz-Weyl decomposition: If then . Thus by . But then is a Dirichlet field and hence must vanish by orthogonality. Finally .
The second step is based upon the following result:
Let be the exact solution and be an approximation. Furthermore, let be as above and let exist, such that
holds for all . Then
|holds for all , for which the tangential trace coincides with the tangential trace of , i.e. , on the boundary . If additionally then|
which holds by (2.9). Hence, . Using Cauchy-Schwarz’ inequality we obtain
Using the trace and extension operators from Remark 6 we obtain the following result:
Let the assumptions of Lemma 10 be satisfied. Then
Here is the constant in the inequality
Let be as in Lemma 10. Then
where is an arbitrary vector field in .
Proof For any and any we have
We remark that the latter estimate is unable to measure adequately the deviation of the divergence of to (this is obvious since even does not need to have any divergence). On the other hand, even if then the semi-norm could not feel the lack of the constraint . The same holds true for the deviation of from the orthogonality to the Dirichlet fields. However, it is not difficult to transform the estimate into a form, in which the estimate is represented in terms of the semi-norm
on , which obviously is a norm on
These two facts can be seen applying the Helmholtz-Weyl decomposition (2.8) and the projection . In particular, by (2.9) replacing by in Theorem 12 would change nothing. In other words, the part of containing the eventually non vanishing divergence term can be added to any term in (3.6) without changing anything. To get the ’full’ norm (3.11) we just add the terms
and the sum of
Of course, the terms in the first equalities make sense for only.
Let be as in Lemma 10 and additionally . Then
holds for any . If even belongs to , i.e. if the approximation satisfies the boundary condition exactly, then is a norm for and we have for all
and the left hand side is a norm for . The estimates (3.6)-(3.14) show that deviations from exact solutions contain weighted residuals of basic relations with weights given by constants in the corresponding embedding inequalities. These are typical features of the so-called functional a posteriori error estimates.
In Theorem 12 and Theorem 14 we can replace the boundary term on the right hand side by or using Lemma 10 and Corollary 11. Especially for numerical applications the first choice is recommendable. Hence, we may assume that is always given by a tangential trace of some vector field , i.e. . Then and we do not have to know the constant .
If the domain is ’simple’ in terms of a vanishing second Betti number, i.e. there are no ’handles’, then there exist no Dirichtlet fields. Thus, for instance, in Theorem 14 the last summand in the respective estimates does not occur.
4 Lower bounds for the error
Now, we proceed to derive computable lower bounds of the error. First, we present the following subsidiary result:
If satisfies (2.6) then and .
Proof We need to show that
Since , we get (by approximation) and (4.1) follows. To be more precise, we select , for which converges in to , using and . Then
Lemma 19 implies
Let and some be given. Then the following three assertions are equivalent:
Let be an approximation. Then
and the supremum is taken over . This estimate is sharp if belongs to the latter space, i.e. if the approximation satisfies the boundary condition exactly.
Proof We start with the obvious identity
Thus, for all we obtain the estimate