A stationary heat conduction problem in low dimensional sets in

We study a basic linear elliptic equation on a lower dimensional rectifiable set S in RN with the Neumann boundary data. Set S is a support of a finite Borel measure μ. We will use the measure theoretic tools to interpret the equation and the Neumann boundary condition. For this purpose we recall the Sobolev-type space dependent on the measure μ. We establish existence and uniqueness of weak solutions provided that an appropriate source term is given. Mathematics Subject Classification 35J20 · 35J70 · 28A33


Introduction
We study here an old problem of determining the stationary heat distribution in a conductor S in the ambient space ⊂ R N , when the conductivity tensor A and the heat sources Q are given. We assume that the conductor S is insulated at the boundary of . Moreover, the set \ S is neither conducting nor it contains any heat source. The main objective of the paper is to investigate the case of "thin" S, i.e., the case when the dimension of S is smaller than N .
In the classical setting, S is a smooth Riemannian manifold. The stationary heat equation takes the following form, Here, denotes the Laplace-Betrami operator, i.e., where div S is a surface divergence and ∇u is the covariant derivative of u. If S is a submanifold of R N , then ∇u = P Du, where P is the projection of R N onto T x S, the tangent space of S at x. Equation (1.1) is the Euler-Lagrange equation of a variational functional, where W 1,2 (S) is the Sobolev space on the Riemannian manifold S. Our main goal is to discuss such variational functionals and to analyze their Euler-Lagrange equations in a more general context, namely when the set S is a lower dimension rectifiable set in R N . To do this we employ the measure theoretic tools and we follow the approach proposed by Bouchitté, Buttazzo and Seppecher in [1]. Namely, we consider a measure μ, supported in S, which is singular with respect to the Lebesgue measure in R N . We assume further that the conductivity tensor A and the heat sources Q are, respectively, a tensor and a scalar valued measure; both measures are absolutely continuous with respect to μ. Applying variational methods, we focus initially on minimizing problem for the functional where ν is the outward normal unit vector of ∂ . However, an effort is required to derive and give meaning to such an equation when μ is a more general measure. In this paper we consider a multijunction measure in the sense considered by Bouchitté and Fragalà in [2]. Our precise assumptions, denoted further by [S], are as follows: Such measure-oriented point of view proved to be very fruitful when dealing with variational problems considered on low dimensional structures in R N , e.g., see the anisotropic shape optimization problem considered by Bouchitté and Buttazzo in [3]. We refer the reader to [4] for an introduction to this theory.
At this point let us mention the main features of the present work. First of all, let us observe, that contrary to the approach presented in [5] and [2], we do not assume a Poincarétype inequality to hold globally on S, allowing in particular occurrence of multidimensional junctions. The issue of the Poincaré inequality is addressed further in Sect. 2 and examples are given in Sect. 5. Secondly, we consider the Neumann boundary condition in (1.3), which was not present in aforementioned papers [1,3,4]. We stress that the meaning of the boundary condition in (1.3) requires clarification. We shall see in Sect. 4 that the correct interpretation relies on the theory developed by Chen and Frid in [6], and it amounts for the Neumann boundary term to be understood as a normal trace of a measure with bounded divergence. Observe also that the ambient space is, in fact, a design region restricting the position of set S. Remark 4.1 at the end of Sect. 4 addresses the lack of influence of on the solutions to (1.3) as long as ∂ is smooth in a neighborhood of ∂ ∩ ∂ S.
And finally, we want to mention the fact that conditions [S] may be relaxed a bit. Namely, we may consider μ = θμ, thereμ satisifes [S] and the density θ is bounded and strictly positive, see (2.9) and (2.10). In this case the same results hold for such μ and their proofs do not require any changes. We explain it, while presenting the argument.
The variational approach requires introduction of the appropriate Sobolev-type space and finding the lower semicontinuous envelope of E. To this end we employ the space H 1,2 μ introduced by Bouchitté and Buttazzo, [1]. Its definition and basic properties are discussed in Sect. 2. In particular, we explain there the notion of the tangential gradient operator ∇ μ for functions from H 1,2 μ and the notion of the tangent space of μ at x, denoted by T μ (x). The relaxation of E, denoted by E μ may be calculated explicitly as in [1]-this is done in Sect. 3, see Proposition 3.1). It is defined as where A μ is given by (3.4). The Euler-Lagrange equation for E μ then takes the form (see Proposition 3.3) where f is the density of the measure Q.
Our main result reads as following: This Theorem follows from Theorem 3.1, Proposition 3.3, and Theorem 4.1. Uniqueness is addressed in Corollary 3.3. We also make a comment of the boundary conditions in Corollary 4.1. Examples of explicit boundary problems are discussed in Sect. 5. In the course of proofs, we will explain that Theorem 1.1 is valid also for μ given by (2.9), where the density θ satisfies (2.10).

Definitions
In the problem we study, is the ambient space, containing S, which is the support of a measure μ. Set S interpreted as the heat conductor. We study an elliptic problem in S, however, we are interested in the Neumann-type boundary conditions on S ∩ ∂ .
We need properly defined Sobolev spaces to study weak solvability of differential equations like (1.3). We require a definition, which is general, yet easy to use. It seems that the Sobolev space H 1, p μ introduced for any p ∈ [1, ∞) by Bouchitté-Buttazzo-Seppecher in [1] is the right tool. We will briefly recall this definition. A recent discussion on other possible definitions of Sobolev space on S is in [7].
We begin with the notion of a tangent space to a measure, there is a number of ways to introduce it. We use the approach proposed in [1] and further exposed in [2,4,8,9] in a more reader friendly way.
We assume that μ is a positive, Radon measure on R N . If we set We define the multifunction N μ : The vector subspace of R N is called the normal space to μ at x. The vector space T μ (x) ⊂ R N is defined as the orthogonal complement of N μ (x) and it is called the tangent space to μ at x. We introduce μ (x, ·) as the orthogonal projection of R N onto T μ (x). Due to measurability of T μ (·), see [1], the projection is μ-measurable.
For any u ∈ D(R n ), we define, for μ-a.e. x ∈ R n , the tangential gradient The space H 1, p μ is defined as a completion of D(R n ) in the following norm This space is a reflexive Banach space for p ∈ (1, ∞). We should stress that we could define the Sobolev space by weak derivatives, this is done in [10]. This space is denoted by W If S = supp μ ⊂ and is an open region in R N and μ(S ∩ ∂ ) = 0, then

On the validity of the Poincaré inequality
An important tool in the analysis of the well posedness of the problem (1.3) is the Poincaré inequality. In case S is a smooth, compact manifold of dimension k < N , the global Poincaré inequality holds, i.e. there is C P > 0 such that However, we are particularly interested in junctions between several thin structures of possibly different dimension, i.e. when S is a finite union of compact manifolds of possibly different Hausdorff dimensions. In such a case (2.1) may not be true. Further in this section we discuss some positive examples, but let us now concentrate on the problems associated with the lack of (2.1) and the possible remedies.
One of the consequences of the global Poincaré inequality (2.1) is the following property: In other words dim ker ∇ μ = dim{u ∈ H 1,2 μ ( ) : ∇ μ u = 0} = 1, so ker ∇ μ contains only constant functions. In a general situation of a multijunction measure this is not necessarily true, i.e. it may happen that dim ker ∇ μ > 1.
The important feature is "how big" the junction area is. Every space that supports a global Poincaré inequality is connected. We expect that removing a set of zero capacity from the space should not affect the Poincaré inequality. In particular, removing a set of zero capacity should not disconnect the space. Examples show that our expectation need not be true. We consider the situation presented on Fig. 1. Since the capacity of a point on a plane is equal to zero, the global Poincaré inequality cannot be valid for this set. A similar situation occurs for the set presented on Fig. 6. We do not want to recall the definition of a (Sobolev) capacity here (we refer the interested reader for instance to the book [11], chapter 4). However, in the situation presented on the Fig. 1, the argument is simple and straightforward. Consider the unit disk on a plane D = {(x, y) ∈ R 2 : x 2 + y 2 < 1} and the function f (x, y) = − log log(1 + (x 2 + y 2 ) −1/2 ). This is an unbounded function which belongs to the Sobolev space H 1,2 (D). Set Then ζ • f ∈ H 1,2 (D) as well. Define a sequence of functions f n as Since the H 1,2 (D) energy of f is finite, the H 1,2 (D)-energy of f n tends to zero. Now, we consider a metric space (S, μ), such that S = X ∪ P Y , where and μ X is a 2-dimensional Hausdorff measure; and μ Y is a 1-dimensional Hausdorff measure and P = (0, 0, 0).

Proposition 2.2
Let us suppose that S and μ are defined above. Then, dim ker ∇ μ > 1.

Proof
We define a sequence of functions u n : S → R as then it follows that u n ∈ H 1,2 μ . By the construction, we also have which implies u o ∈ H 1,2 μ which in turn falsifies the Poincaré inequality (2.1) on S and it shows that (2.2) in this situation is not true. As a result, our claim holds.

Let us denote
This is a finite dimensional subspace of L 2 ( , μ). We define as the orthogonal projection onto V (with respect to L 2 scalar product). Set d = dim V . We can split the set {1, . . . , J } from the condition [S] into a final family of disjoint sets of indices The sets J l , l = 1, . . . , d are defined in the following way. The set and J l are "maximal" in the sense that for any proper subset K ⊂ J l , the characteristic function of i∈K S i does not belong to V .
Then the projection P in (2.3) can be expressed as For a multijunction measure satisfying [S] the following weaker version of the Poincaré inequality holds.

5)
where P l is defined in (2.4).
Let us mention in advance that, the regularity of solutions depends in an essential way on whether the set S supports the inequality (2.1) or (2.5). It is known that if a metric measure space (the set S may be treated as such) supports the Poincaré inequality (2.1), one can expect solutions to be continuous. The weaker version (2.5), which is sufficient for the existence result, is not enough to ensure regularity, as shown on examples presented in Sect. 5.
On the other hand, as it is observed in the vast literature concerning analysis on metricmeasure spaces, there are many examples of spaces for which validity of Poincaré-type inequality (2.1) can be asserted. It is impossible to discuss this broad topic in all the details within this paper. Various versions of the Poincaré-type inequalities are formulated in different context. However, just to give some justification for our claim that in some examples, one can indeed assert validity of the Poincaré inequality in the form (2.1), we will sketch the framework of metric-measure spaces. An interested reader may find more information and references e.g., in the papers of Semmes [12,Section B], of Hajłasz and Koskela [15], of Heinonen and Koskela [13,Section 6] and in the book [11].
We need to formulate some definitions which are necessary within the framework. Let (X , d) be a connected and complete metric space equipped with a Borel regular measure μ. We say that: for all balls B R in X of radius 0 < R < diam X . • X is linearly locally contractible if there is a constant C such that every small ball is contractible inside a ball whose radius is C times larger.
Let U be an open set in X and let u be an arbitrary real-valued function in U . We say that a Borel function g : whenever γ is a rectifiable curve joining x and y in U . In many cases (in particular in the Euclidean setting, when the function u is differentiable) one can take simply g = |∇u|. Let us observe that this is our case, due to Proposition 2.1.
We say that X admits a local weak whenever u is a bounded continuous function in a ball λB and g is its very weak gradient. The constants C p ≥ 1 and λ ≥ 1 should be independent of B and u. Inequality (2.8) is termed "weak" because we allow a larger ball on the right-hand side than on the left. In many cases this weak estimate can be used to iterate so as to yield an inequality with the same ball on both sides. If the term "weak" is omitted, we assume λ = 1.
With such preparations, we may now state the results due to Semmes (see [12, Section B]): is a q-regular metric space which is also locally linearly contractible and a connected orientable topological manifold of dimension q. Then X admits a local (1, p)-Poincaré inequality for all p ≥ 1.

Theorem 2.3
Suppose that X is a connected, finite simplicial complex of pure dimension q > 1 such that the link of each vertex is connected. Then X admits a local (1, q)-Poincaré inequality.
Observe first that both theorems assert validity of the local Poincaré inequality, not just its local weak version. This is a good starting point, but still, local version is not sufficient for our purpose, we want the global version (2.1). Then, even as Theorem 2.3 is much more general, and actually very handy in our situation (our assumptions for the set S actually make it a simplicial complex), we want to avoid introduction of yet another language of simplicial geometry. Therefore to justify passing from the local version (2.8) of the Poincaré inequality to the global one (2.1), we shall start our arguments with the Theorem 2.2.
As observed by Heinonen and Koskela [13], two spaces that support a local (weak) version of the Poincaré inequality may be glued along a sufficiently large common part to produce a new metric space that also supports a Poincaré-type inequality. To make the statement more precise, suppose that X and Y are two locally compact q-regular metric measure spaces. Suppose also that A is a closed subset of X that has an isometric copy inside Y , i.e. there is an isometric embedding i : A → Y . In the following, we understand this embedding as fixed and think of A as a closed subset of both X and Y . We form a space which is the disjoint union of X and Y with points in the two copies of A identified. This space has a natural metric which extends the metrics from X and Y : the distance between Furthermore, the measures on X and Y add up to a measure μ on X ∪ A Y that is obviously q-regular. Recall that the Hausdorff s-content of a set E in a metric space is the number where the infimum is taken over all countable covers of the set E by ball B i of radius r i .

If both X and Y admit a local weak (1, p)-Poincaré inequality, then X ∪ A Y admits a local weak (1, p)-Poincaré inequality as well.
Finally, passing from a local weak version of the Poincaré-type inequality (2.8) to the global version (2.1) can be done if the measure μ is doubling and the space satisfies the socalled chain condition. Roughly speaking, the chain conditions means that each two points of the space can be joint by a sequence of balls with controlled radius and "thick" (in the sense of measure) intersection. An elegant and elementary proof of the fact was given by Hajłasz

A more general measure
In the Introduction, we mentioned a possible generalization of our results. Namely, we now assume thatμ satisfies [S] and we consider μ = θμ, (2.9) where the density is not only in L ∞ (S), but also θ is strictly positive. We notice that such measures share the tangent space. We set μ = θμ. Then, Hence, H If we inspect the definition of N μ , then we notice that N μ and Nμ must coincide. As a result, we see that N μ = Nμ and N μ = Nμ(x) forμ a.e. x ∈ S. Finally, forμ a.e. x ∈ S (2.11) and our first claim follows. The second claim follows from (2.11) and the definition of H 1, p μ .

The elliptic problem
We would like to use the calculus of variations to establish existence of solutions to (1.3) in the space H 1,2 μ . For this purpose we recall the definition of E(u) for H 1,2 μ . We have to specify the assumptions on A and Q.
In principle, u is an element of L 2 ( , μ), so we must specify properties of Q appropriately. This is why we make the following assumption, Q is a measure of a finite total variation supported in , absolutely continuous with respect to μ, and there is f ∈ L 2 ( , μ) such that Q = f μ. This assumption in particular makes elements of H 1,2 μ measurable with respect to Q. Now, we specify our assumptions on A. We consider A ∈ L ∞ ( ; M(N × N ), μ) and for a.e. x the matrix A(x) is symmetric and T μ (x) ⊂ Im A(x). (3.2) We have to present the positivity of A in a way suitable for dealing with ∇ μ u in the tangent space to μ. We require that at μ-almost every x ∈ , we have After these preparations, we state our first result.
Proof We may assume for the sake of simplicity of notation that f = 0. Indeed, the term f u dμ is weakly continuous, hence we have to deal with the lower semicontinuity of the quadratic term.
Step 1. Let us assume first that we are given a matrix B ∈ L ∞ ( , M(N × N ), μ) and for μ-a.e. x ∈ matrix B(x) is symmetric and positive definite. Then, for μ-a.e. x ∈ , the matrix B(x) defines a scalar product in R N : Let us choose vector fields e 1 , . . . , e N −k such that for μ-a.e. x ∈ : • they are linearly independent, • they span (T μ (x)) ⊥ , • they are orthonormal with respect to the scalar product defined by the matrix B, i.e. for μ-a.e. x ∈ The μ-measurability of these vector fields is discussed in Step 4.
Observe that for μ-a.e. x ∈ , the orthogonal projection onto (T μ (x)) ⊥ with respect to the scalar product B(x) is defined as F B (x, ∇u) dμ. By [1], E B μ , the relaxation of E B , is given by The minimization problem in the definition of F B μ may be written as A simple differentiation yields optimality conditions that is It is easy to see that for any ξ ∈ R N we have B μ ξ ∈ T μ . Indeed, Moreover, B μ is strictly positively definite on T μ (x). Indeed, for ξ ∈ T μ (x) we have Step 2. We assumed that A satisfies where P(x) is the orthogonal (with respect to the euclidean product in R N ) projection onto ker A. Hence, for every > 0 and for μ-a.e. x ∈ matrix B (x) is symmetric and positive definite and formula (3.5) applies to B .
Step 3. We claim that if B is defined by (3.6), then B μ = A μ , where A μ is given by (3.4). Indeed, the summation in (3.5) may be split, the first l ≤ N − k vectors belong to Im A ∩ (T μ ) ⊥ , while vectors e i , i = l + 1, . . . , N − k span the kernel of A. We note Then, B μ takes the following form, Step 4. We claim that the mapping x → A μ (x) is μ-measurable. We know that the multivalued map x → T ⊥ μ (x) is μ-measurable, (see [9, page 476]), so is As a result, there is a measurable selection of 1 (x), which we will call e 1 .
In the next step we consider a multivalued function Of course it is measurable, we will denotes its selection by e 2 . After the finite number of steps we will end up with measurable vector fields e i (x), i = 1, . . . , l, which are necessary in the formula for A μ .
Step 5. We will check that the lower semicontinuous envelope of E is E μ . By definition Step 3 implies that E μ = E B μ and E B μ is semicontinuous. If we denote by bar the lower semicontinuous envelope, then we see, Hence, E μ (u) =Ē(u).
Interestingly, the conclusion of Proposition 3.1 does not change if we consider μ of the form (2.9). Here it is: Proof Let us notice that (3.4) is well-defined as long as the tangent space T μ (x) exists, i.e. also if μ is of the form (2.9).
Let us now introduceĀ Due to (2.10) we see thatμ,Ā andQ satisfy the assumptions of Proposition 3.1. Hence, the relaxation of functional E, is given by Eμ, whereĀμ is defined in (3.4). We notice that where we use that fact that H 1,2 μ = H 1,2 μ . It is easy to check that the structure of the formula (3.4) implies thatĀμ once we see that vectors in R N orthogonal with respect to inner product defined by A(x) are orthogonal with respect to inner product defined byĀ(x). However, this is obvious since both inner products are propotional, with a positive proportionality constant.
Thus, for any u ∈ H 1,2 μ we have, We have to describe ker div A μ ∇ μ , what is necessary, before we can specify Q. Actually, we prove: In other words, We notice that ∇ϕ(x) = ∇ μ ϕ(x) We may take a sequence ϕ n ∈ D(R N ) such that ∇ μ ϕ n converges to ∇ μ u. This yields, Combining this with ellipticity of A μ , see (3.3), yields the claim.

Remark 3.1
The same proof works in the case of μ as in Corollary 2.1.
After these preparations we can state one of our main results.
where P is the orthogonal projection onto ker ∇ μ described by (2.4).

Proof
We have already computed E μ , the relaxation of E. Actually, we study minimizers of E μ . Let us suppose that {u k } ⊂ H is a sequence minimizing E μ , i.e., We may assume that for all k ∈ N we have, Due to (3.3), we end up with, Now, the Theorem 2.1 and Pu k = 0 yield that there exists C > 0 such that Combining this with (3.7) and Young's inequality yields, In other words the minimizing sequence is bounded in H 1,2 μ . Due to the results of [1] we deduce that there is a weakly convergent subsequence (not relabeled), with limit u. Now, we invoke the lower semicontinuity results so that we deduce, Uniqueness of minimizers will be treated separately.
The claim made in Theorem 3.1 is valid also for μ as in Corollary 2.1.

Corollary 3.2 Let us assume that the hypotheses of
Proof Since u is a minimizer, then for any test function ϕ ∈ H 1,2 μ , (in particular we may take ϕ ∈ C 1 0 (R N )), we have Thus, After replacing ϕ with tϕ, where t is real, we will deduce that (3.8) holds.

Corollary 3.3
Let us suppose that the assumptions of the previous theorem hold. If u 1 and u 2 are two minimizers of E μ , which are perpendicular to ker∇ μ , then u 1 = u 2 .
Proof Let us take u = u 2 − u 1 and take the difference of (3.8) corresponding to u 1 and u 2 , Since u = u 2 − u 1 can be approximated in the H 1,2 μ -norm by a sequence of C 1 functions, we may take ϕ = u in the identity above, i.e. u ∈ ker ∇ m u. Since we assumed that u ∈ (ker ∇ μ ) ⊥ , we deduce that u = 0, i.e. u 2 = u 1 .
We do not specify any boundary conditions so we expect that the solution satisfies the so-called natural boundary condition, however, in a weak form suitable for H 1,2 μ .

Remark 3.2
The arguments used in the proofs of Proposition 3.3 and Corollary 3.3 are valid if μ is as in Corollary 3.1.

Remark 3.3
Let us remark that in the case S is a smooth compact manifold of codimension 1, μ = H N −1 S, and A is the identity matrix, the equation (1.3) takes the familiar form of where is the Laplace-Beltrami operator on S.

On the boundary condition
After [6], we introduce the notion of a vector valued measure with bounded divergence. The advantage of measures whose divergence has a finite total variation is that one can define the trace of their normal component of ∂ . We recall, see [6, Theorem 2.2],

Proposition 4.1 We assume that the boundary of is smooth. If both F and its divergence have a bounded total variation, there exists a continuous linear functional
Actually, the original statement of the Proposition is more general, in particular less smoothness of the boundary of is sufficient. Namely, it suffices if ∂ is Lipschitz deformable, see [6,Definition 2.1]. In addition, [F, ν]| ∂ could be defined as a functional on the space Lip (γ , ), γ > 1, which is larger than C 2 (∂ ), see [6]. However, such generality is not necessary here, so we choose the simplicity of exposition.
We want to show that if u is a minimizer of E, then F := A μ ∇ μ u has divergence with a finite total variation.
Let us suppose that η k ∈ C ∞ c ( ) are the cut-off functions such that η k (x) = 0 for x ∈ at a distance smaller than 1/(k + 1) from ∂ and equal to 1 for x further than 1/k from ∂ . Then, since η k ϕ ∈ C 1 c ( ) and by using of the Euler-Lagrange Eq. (3.8) with (1 − η k )ϕ as a test function, we obtain If ω is a measure with a finite variation, then because of the Lebesgue dominated convergence Theorem. If we apply this remark to ω = Q and ω = div A μ ∇ μ u, then we deduce that for all ϕ ∈ C 2 (¯ ). Our claim follows. The corollaries discussed in the next section prompt a more general observation on the role of . We notice the lack of dependence of the boundary conditions on , as long as ∂ is smooth in a neighborhood of ∂ ∩ ∂ S. The set does not enter the problem of minimization of functional E μ . Moreover, by Corollary 3.3, the solution to the minimization problem is unique. As a result, if i , i = 1, 2, satisfy [S] and u i , i = 1, 2, are the corresponding minimizers of E μ , then u 1 = u 2 . However, the boundary conditions are a bit different, because ν 1 need not be equal ν 2 . We may vary as long as conditions [S] are satisfied. In particular, at any x ∈ S ∩ ∂ the normal ν may be any vector perpendicular to T ∂ S (x). Hence, we obtain: Proof This is so, because for any smooth mapping S x → ν(x), where ν(x) ⊥ T ∂ S we can find such satisying [S] that ν(x) ⊥ ∂ at x. Hence, [A μ ∇ μ u, ν] ∂ = 0 holds in the sense of the Chen-Frid, [6].
In Sect. 3 we established results corresponding to μ given by (2.9) with the density θ satisfying (2.10). We did it by absorbing the density intoĀ,Q. It can be easily checked that the same techique works here, in this section.

Examples of the Neumann boundary problems
We present below several examples of Neumann boundary problems on various sets S. The first two examples deal with situations when the Poincaré inequality (2.1) holds (see Sect. 2.2). We further present an example (see Proposition 5.3) showing that discontinuous solutions actually occur, even for regular data, when the set S supports only the weaker inequality (2.5).
Let us suppose that is a ball centered at 0 with radius 1, T is an inscribed isosceles triangle and S 1 , S 2 , S 3 is the set of radii connecting the center of with the vertices of T (see Fig. 3). We set S = S 1 ∪ S 2 ∪ S 3 . Moreover, μ = H 1 S. We set Q = 3 i=1 a i χ S i . We can state the following fact. Here, s is the arc-length parameter counted from the center of .
Finally, the solution of (5.1-5.2) with S u dH 1 = 0 has the following form, , except for x being the center of the ball. Hence, Proposition 3.1 yields div A μ ∇ μ u = d 2 u ds 2 in the interior of S i , i = 1, 2, 3. Also the form of boundary data (5.2) follows.
The assumption a 1 + a 2 + a 3 = 0 means that Q is perpendicular to the kernel of ∇ μ so that we may apply Theorem 1.1.
The form of any solution to (5.1-5.2) follows from solving the ODE's (5.1). While doing so, we keep in mind the requirement that the solution must be continuous at s = 0, i.e. at the center of the ball . The solution is determined up to an additive constant, which can be computed due to the zero average condition, S u dH 1 = 0.
We can also study (1.3) on sums of manifolds of dimension k ≥ 1. We describe the set S depicted on Fig. 4.

Proposition 5.2
We assume that A = I d, is the unit ball B(0, 1) ⊂ R 3 and S = D 1 ∪ D 2 are two great disks of B(0, 1), where ν i are normal to D i in the plane containing D i , i = 1, 2. The condition of orthogonality to the kernel of ∇ μ is D i u dH 2 = 0, i = 1, 2.
Proof We easily see that, due to Proposition 3.1, we will have Hence, the form of relaxed equations follows. The boundary conditions are addressed in Corollary 4.1.
We notice that A μ depends upon x ∈ despite A being constant. Moreover, the kernel of A μ is two-dimensional.
Due to the generalized Poincaré inequality (2.5), we can study (1.3) on sums of manifolds of different dimensions. We analyze an example of an equation on the domain depicted on Figure 5. We define the design region to be B(0, √ 2) ⊂ R 3 . We set, In the case of this domain, we do not expect the global Poincaré inequality (2.1) to hold. Proposition 5. 3 We assume that A = I d and We assume that Q = Q 1 χ S 1 + Q 2 χ S 2 , i.e. Q| S 3 ≡ 0. We will require that S i Q i dμ = 0, i = 1, 2.
We set Then, u, the solution to (1.5) in S with S i u dH 2 = 0, i = 1, 2 and S 3 u dH 2 = 0 is discontinuous.
Thus, the equation divA μ ∇ μ u + Q = 0 becomes The radial symmetry of q implies that 0 = The above Eq. (5.3), considered on S 1 , due to the radial symmetry takes the following form u rr + 1 r u r + q = 0. In order to solve it, we have to pay a bit of attention to the boundary conditions. Since ν = 2 −1/2 (−1, x 2 , x 3 ) for (−1, x 2 , x 3 ) ∈ ∂ S 1 , then the Neumann boundary condition where n is the outer normal to S 1 in the plane x 1 = −1. We notice that due to (5.4), the above boundary conditions are automatically satisfied, provided that a = 0. Moreover, a = 0 is necessary for u to be an element of H 1,2 μ . The formula for the solution is as follows, We may choose q so that b > 0. Due to the symmetry of Q, we deduce that u(−1, x 2 , x 3 ) = −u(1, x 2 , x 3 ). Thus, in particular u(−1, 0, 0) = b = −u(1, 0, 0). Since we set Q| S 3 = 0, then we study the problem of minimizing 1 −1 on S 3 , where u ∈ H 1,2 (−1, 1) and 1 −1 u dx 1 = 0. As a result u| S 3 = 0. We conclude that we constructed a discontinuous solution when the forcing term is continuous.
The last example, on the domain depicted on Fig. 6, illustrates how the relaxation procedure affects the matrix A. Then, A μ = e 2 ⊗ e 2 + 1 2 e 1 ⊗ e 1 χ S 1 + 1 2 e 3 ⊗ e 3 χ S 2 (5.5) and the boundary conditions take the following form, The condition of orthogonality to the kernel of ∇ μ is S i u dH 2 = 0.