TwoSided A Posteriori Error Bounds for ElectroMagneto Static Problems
Abstract
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, electromagneto statics
AMS MSCClassifications 65 N 15, 78 A 30
Dedicated to the anniversary of Prof. Nina Nikolaevna Uraltseva
Contents:
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 electromagneto static boundary value problem and any approximation from the corresponding energy space. We discuss the method with the paradigm of a prototypical electromagneto 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 curlfree 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 HelmholtzWeyl 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 seminorm generated by the operator ) is equal to some norm of the socalled 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 (matrixvalued 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.
(1.1) 
By we denote the usual scalar Hilbert space of square integrable functions over and by the Hilbert space of realvalued 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 socalled 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 socalled Betti numbers of and depend only on topological properties of the domain (for a detailed presentation see [10]). 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 electromagneto 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
such that
where the homogeneous boundary conditions are to be understood in the weak sense.
This coupled problem is equivalent to an electromagneto 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
(1.2)  
(1.3)  
(1.4)  
(1.5) 
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 electromagneto static Maxwell problem in second order form: Find an electric field in , such that
(2.1)  
(2.2)  
(2.3)  
(2.4) 
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 HelmholtzWeyl 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
(2.5) 
Of course, there exist more general variants of PoincaréFriedrich’s estimate (2.5) for vector fields. Here, we refer to Remark 8.
Now, let be some vector field in satisfying the boundary condition (2.3) in the generalized sense, i.e., . The generalized solution
of (2.1)(2.4) is then defined by the relation
(2.6) 
If then by the CauchyScharz 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 LaxMilgram’s theorem.
First, we note some HelmholtzWeyl decompositions of , i.e. decompositions into solenoidal and curlfree fields, which will be used frequently throughout our analysis.
Lemma 1
can be decomposed as
(2.7)  
and  
(2.8) 
where all closures are taken in and . Moreover,
Remark 2
Let us denote the orthogonal projection onto in (2.8) by . Then we have for all
(2.9) 
and for all
The latter line can be written in a more compact and precise way as
Remark 3
The following theorem states the main existence result.
Theorem 4
Remark 5
Remark 6
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 [23] 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 HelmholtzWeyl 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 .
Remark 7
Henceforth, we assume that is given by a tangential trace of some vector field .
Remark 8
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 [18]. 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 divergencefree 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
(3.1) 
where is a linear and continuous functional over as well as over . Furthermore, does not depend on the exact solution .
Remark 9
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 HelmholtzWeyl 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:
Lemma 10
Let be the exact solution and be an approximation. Furthermore, let be as above and let exist, such that
holds for all . Then
(3.2)  
holds for all , for which the tangential trace coincides with the tangential trace of , i.e. , on the boundary . If additionally then  
(3.3) 
Proof We use the HelmholtzWeyl decomposition (2.8) and the projection from Remark 2. We consider a vector field with and define the vector field
which holds by (2.9). Hence, . Using CauchySchwarz’ inequality we obtain
and thus .
By the triangle inequality we get (3.2).
(3.3) is trivial setting .
Using the trace and extension operators from Remark 6 we obtain the following result:
Corollary 11
Let the assumptions of Lemma 10 be satisfied. Then
(3.4) 
Here is the constant in the inequality
(3.5) 
Proof
Setting in (3.2) and using
(3.5) proves (3.4).
We note that (3.3) follows directly from the corollary as well.
Theorem 12
Proof For any and any we have
(3.7) 
Combining (3.1) and (3.7), we obtain for all
(3.8) 
By CauchySchwarz’ inequality, PoincaréFriedrich’s estimate (2.5) and (1.1) we estimate the right hand side of (3.8)
(3.9)  
(3.10) 
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 seminorm 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 seminorm
(3.11) 
on , which obviously is a norm on
Remark 13
These two facts can be seen applying the HelmholtzWeyl 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.
Theorem 14
Let be as in Lemma 10 and additionally . Then
(3.12) 
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
(3.13) 
Remark 15
If satisfies the prescribed boundary condition and is solenoidal and perpendicular to Dirichlet fields, then (3.6) or (3.12), (3.13) imply for all
(3.14) 
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 socalled functional a posteriori error estimates.
Remark 16
Remark 17
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 .
Remark 18
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:
Lemma 19
If satisfies (2.6) then and .
Proof We need to show that
(4.1) 
Using from Remark 2, we obtain provided that . Thus, by (2.6) and the fact that , we get
(4.2) 
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
Remark 20
Let and some be given. Then the following three assertions are equivalent:

and .

and

and
Theorem 21
Let be an approximation. Then
where
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