Solenoidal difference quotients and their application to the regularity theory of the $p$-Stokes system

We prove existence of a solution to the divergence equation satisfying a new Bogovski-type estimate for the difference quotients. This enables us to give an alternative proof of the interior regularity of the solution to the $p$-Stokes problem, completely avoiding the pressure. Moreover, as a key preliminary result we prove boundedness of Calder\'on-Zygmnud operators with standard kernels in weighted Lebesgue and Orlicz spaces over a general domain.


Introduction
We are concerned with the question of interior regularity of the weak solution of the steady Stokes approximation for flows of shear thinning fluids. It is given by − div S(Du) + ∇π = f in Ω, div u = 0 in Ω, where Ω ⊂ R n , n ≥ 2, is a bounded domain. In here, u = (u 1 , . . . , u n ) ⊤ denotes the unknown velocity vector field, and π the unknown scalar pressure, while the external body force f = (f 1 , . . . , f n ) ⊤ is given. The extra stress tensor S depends only on Du := 1 2 (∇u + ∇u ⊤ ), the symmetric part of the velocity gradient ∇u. The relevant example we have in mind is S(Du) = µ(δ + |Du|) p−2 Du , (1.2) with p ∈ (1, 2], δ ≥ 0, and µ > 0. Notice that despite various efforts the optimal global regularity of this problem is still open (see [8,1,4] for partial results). However, the interior regularity we focus on here is well-known (see [15,19,4]). The standard proof uses localized difference quotients in each direction. Due to the localization, the corresponding test function is not solenoidal anymore, therefore appropriate properties of the pressure have to be used.
Here we modify this approach and use a solenoidal version of the localized difference quotient. Thus, we can completely avoid the pressure π. This is a completely new approach even for the classical Stokes problem, i.e., p = 2 in (1.2). To make it possible, we show that the solution of the divergence equation obtained via the Bogovski formula (see [5,6]) satisfies an additional estimate for the difference quotient. The proof of this is based on estimates of singular Calderón-Zygmund operators generated by standard kernels in arbitrary domains Ω ⊂ R n (see Theorem 3.4). This result is of independent interest since it shows that an analogy of the classical Calderón-Zygmund estimates (cf. [7,18]) in the whole space also hold in arbitrary domains. To prove it, we employ ideas from [5,6,16,10,23], where the divergence equation is treated, and modify them for our purposes. An important feature of our result is a careful tracking of the dependence of the constants on various quantities, which is missing in the literature.

Preliminaries
2.1. Notation. Throughout the text, we use the symbols C, c to denote generic constants which may change from line to line but are independent of "crucial" quantities. In many cases, the dependence of such constants on various quantities will be explicitly specified. Furthermore, we write f ∼ g if there exist constants c, C > 0 such that c f ≤ g ≤ C f .
A set Q ⊂ R n is called an (open) cube if there exist a 1 , . . . , a n ∈ R and ℓ(Q) > 0 such that Q = {x ∈ R n a i < x i < a i + ℓ(Q), i = 1, . . . , n}.
The number ℓ(Q) is the side length of Q and the point c := (a 1 + 1 2 ℓ(Q), . . . , a n + 1 2 ℓ(Q)) ⊤ is the center of the cube Q, which can be also denoted by Q = Q(c, ℓ(Q)). For α > 0 and a cube Q = Q(c, ℓ(Q)) we use the notation α Q := Q(c, αℓ(Q)). Moreover, a cube Q is said to be dyadic if there exist j i , . . . , j n , k ∈ Z such that Q = x ∈ R n j i 2 k < x i < (j i + 1)2 k , i = 1, . . . , n .
Notice that, in the above definition, we consider only cubes with faces parallel to the coordinate axes.
While working with function spaces, we do not distinguish between spaces of scalar, vectorvalued or tensor-valued functions. However, we denote vectors by boldface lower-case letters (e.g., u) and tensors by boldface upper-case letters (e.g., S). For vectors u, v ∈ R n , the standard tensor product u ⊗ v ∈ R n×n is defined as (u ⊗ v) ij := u i v j , and the symmetric tensor product as u . The scalar product of vectors is denoted by u · v = n i=1 u i v i and the scalar product of tensors is denoted by A · B := n i,j=1 A ij B ij . If E ⊂ R n is a measurable set, |E| denotes its Lebesgue measure. In the following definitions of function spaces, we always assume E to be a domain in R n , with a sufficiently smooth boundary, if needed. We use standard Lebesgue spaces (L p (E), · p ) and Sobolev spaces (W k,p (E), · k,p ). If the underlying domain E needs to be indicated, we denote the respective norms by · L p (E) and · W k,p (E) . If p ∈ (1, ∞), we denote by p ′ the conjugated exponent p ′ := p p−1 . Besides the standard L p spaces we will also consider their weighted variants. A weight is any measurable function ω : R n → [0, ∞). If p ∈ [1, ∞) and ω is a weight, the space L p ω (E) consists of all measurable functions f on E such that If p > 1, an A p weight is a weight such that ω and ω 1−p ′ are locally integrable and [ω] Ap := sup where B ⊂ R n are balls. In this case, we write ω ∈ A p . If p = 1, a weight is called an A 1 weight if 0 < ω < ∞ a.e. and [ω] In here, M is the maximal operator (with respect to non-centered balls), defined by .
The symbol supp f denotes the support of a function f . The set of all compactly supported, smooth functions defined on E is denoted by . If E is bounded, the space W 1,p 0 (E) may be equipped with the gradient norm ∇ · p , thanks to the Poincaré inequality. Next, we denote by W 1,p 0,div (E) the subspace of W 1,p 0 (E) consisting of solenoidal vector fields u, i.e., such that div u = 0. By L p ω,0 (E) we denote the subspace of L p ω (E) consisting of functions f with vanishing mean value, i.e., such that´E f (x) dx = 0.
Orlicz and Sobolev-Orlicz spaces also appear frequently in this article. We briefly present their elementary properties here, for details we refer to [20,21].
An N -function is a continuous, nonnegative, strictly increasing and convex function ψ on [0, ∞) which additionally satisfies ψ(0) = 0, lim t→∞ Furthermore, we define If ∆ 2 (ψ) < ∞, we say that ψ satisfies the ∆ 2 -condition. From now on, let us assume that ∆ 2 (ψ) < ∞ and ∆ 2 (ψ * ) < ∞. Then we denote by L ψ (E) and W 1,ψ (E) the classical Orlicz and Sobolev-Orlicz spaces, respectively. More precisely, f ∈ L ψ (E) if the modular is finite, and f ∈ W 1,ψ (E) if both f and ∇f belong to L ψ (E). When equipped with the Luxemburg norm the space L ψ (E) becomes a Banach space. The same holds for the space W 1,ψ (E) when equipped with the norm · ψ + ∇ · ψ . Notice that the dual space (L ψ (E)) * can be identified with the space L ψ * (E). Furthermore, by W 1,ψ 0 (E) we denote the closure of C ∞ 0 (E) in W 1,ψ (E). If E is bounded and sufficiently regular, the Poincaré inequality for Orlicz modulars (see [24,Lemma 3]) implies that W 1,ψ 0 (E) may be equipped with the gradient norm ∇ · ψ . By L ψ 0 (E) and C ∞ 0,0 (E) we denote the subspaces of L ψ (E) and C ∞ 0 (E), respectively, consisting of functions f such that´E f (x) dx = 0.
If ψ is an N -function satisfying ∆ 2 (ψ) < ∞ and ∆ 2 (ψ * ) < ∞, then for all ε > 0 there exists a constant c ε > 0 such that for all s, t ≥ 0 one has Let F : R n → R n×n be a measurable tensor field (or a vector field or a real-valued function) and h > 0. Then we define the difference quotients of F as follows: We will also use the notation Elementary calculations show that we have the following variant of the product rule for For F, G ∈ L 1 (E), extended by zero outside E, the partial integration formulâ holds. Moreover, for every h 0 > 0, every open set E ⊂ R n , every F ∈ W 1,ψ loc (R n ) and all The proof of this assertion works in the same way as in the classical L p setting, if one additionally uses the Jensen inequality. Also, the converse statement can be proved as in the L p setting (see [14]).
The set of all kernels satisfying the above conditions with a given constant κ 1 will be denoted by SK(E, κ 1 ).
Remark. As it is common in the literature, it is possible to replace the right-hand side in (2.10) and (2.11) by Although the results could be obtained in this generalized setting as well, we restrict ourselves to the case δ = 1 to avoid further complications.
Definition 2.12. We say that a linear operator T on holds whenever the right-hand side is well-defined. For a given kernel K and ε > 0 we define the truncated kernel K ε by and denote by T ε the operator generated by the kernel K ε .
Definition 2.14. Let K be a kernel and let E ⊂ R n be a domain. We call K a Calderón-Zygmund kernel with respect to E, if there exists a constant κ 2 ∈ (0, ∞) such that the function N : satisfies the following conditions: The set of all kernels satisfying the above conditions with a given constant κ 2 will be denoted by CZ(E, κ 2 ).
Remark. One can replace (2.17) by where q ∈ (1, ∞), and still retain the relevant properties of the operator generated by the kernel. For the sake of simplicity, we use the condition only with q = 2.
Let us define what it means that a tensor field S has a (p, δ)-structure.
For details on this matter, see [11,22]. For a tensor P ∈ R n×n we denote its symmetric part by The function ϕ satisfies, uniformly in t and independently of δ, the equivalences In case that the (one-sided) derivative ϕ ′′ (0) does not exist, we assume that ϕ ′′ (t) t is continuously extended by zero for t = 0. We define the shifted N -functions {ϕ a } a≥0 (cf. [11,12,22]) by Note that the family {ϕ a } a≥0 satisfies the ∆ 2 -condition uniformly with respect to a ≥ 0, i.e., , S(P) = S P sym whenever P ∈ R n×n , and S(0) = 0. We say that S has a (p, δ)-structure if for some p ∈ (1, ∞), δ ∈ [0, ∞), and the N -function ϕ = ϕ p,δ (cf. (2.18)) there exist constants γ 0 , γ 1 > 0 such that the inequalities n i,j,k,l=1 are satisfied for all P, Q ∈ R n×n with P sym = 0 and all i, j, k, l = 1, . . . , n. The constants γ 0 , γ 1 , and p are called the characteristics of S.
Remark. An important example of a tensor field S having a (p, δ)-structure is given by S(P) = ϕ ′ (|P sym |)|P sym | −1 P sym . In this case, the characteristics of S, namely γ 0 and γ 1 , depend only on p and are independent of δ ≥ 0.
Suppose that a tensor field S has a (p, δ)-structure. Then we define its associated tensor field F : R n×n → R n×n sym by The connection between S, F and {ϕ a } a≥0 is best explained by the following proposition (cf. [11,22]).
Proposition 2.21. Let S have a (p, δ)-structure, and let F be defined in (2.20). Then The constants depend only on the characteristics of S.
For a detailed discussion of the properties of S and F and their relation to Orlicz spaces and N -functions we refer the reader to [22,3]. In what follows, we shall work only with S(P) and F(P), where P is a symmetric tensor. Therefore, we can drop the superscript " sym " in the above formulas.
If S has a (p, δ)-structure, from Proposition 2.21 we easily obtain the following equivalences: (2.23) The equivalence constants depend here only on the characteristics of S. All assertions from this section may be proved by easy manipulations of definitions, and we omit their proofs.

Calderón-Zygmund estimates
Our interest lies in estimates concerning Orlicz modulars. However, we are going to prove the results first in weighted L p spaces. The following known extrapolation principle (see [9,Theorem 4.15]) offers an elegant connection between the two settings.
Theorem 3.1. Let p 0 ∈ [1, ∞) and let F be a family of pairs of nonnegative measurable functions on R n . Suppose that there exists a nondecreasing function holds for all (f, g) ∈ F and all weights ω ∈ A p 0 . Then for every N -function ψ satisfying holds for all (f, g) ∈ F.
The assumption on monotonicity of C 0 with respect to [ω] Ap 0 is necessary for the proof of Theorem 3.1. However, this is rarely explicitly mentioned in the literature.
(ii) Careful tracking of the constants in the proof of [9,Theorem 4.15] reveals that the constant C in Theorem 3.1 depends on C 0 , ∆ 2 (ψ) and ∆ 2 (ψ * ). To check this, one may use [9,Proposition 3.5]) as well as the boundedness of the maximal operator in Orlicz spaces (see [17,Theorem 2.2]). The latter may be expressed by the inequalityˆR where c is the weak (1, 1)-constant of the maximal operator, which depends on n only. In fact, one can show that In here, c is a constant depending on ∆ 2 (ψ), ∆ 2 (ψ * ), the exponent α depends on ∆ 2 (ψ), and ρ ψ (M ) is the constant in the modular estimate of the maximal operator (it can be estimated by ∆ 2 (ψ * ) ∆ 2 (ψ) β , where β depends on n).
The auxiliary result below is a simple version of the Whitney covering lemma (cf. [2, p. 348]).

Proposition 3.3.
Let Ω ∈ R n be a domain. Then there exists a sequence {Q j } j∈N of dyadic cubes satisfying: Proof. Let x k k∈N be a dense sequence of points in Ω. Define Q 1 as the largest 1 dyadic cube Q such that x 1 ∈ Q and dist (Q, ∂Ω) > 4 diam Q. Next, suppose that Q 1 , . . . , Q j−1 are defined and let k ∈ N be the smallest index such that x k / ∈ j i=1 Q i . Then define Q j as the largest dyadic cube Q such that x k ∈ Q and dist (Q, ∂Ω) > 4 diam Q. The sequence {Q j } j∈N obviously satisfies (ii) and (iii). To verify (i), let x ∈ Ω be arbitrary. There exists a dyadic cube Q such that x ∈ Q and dist (Q, ∂Ω) > 4 diam Q. By density, x k ∈ Q for some k ∈ N. Then necessarily Q ⊆ Q j for some j ∈ N, hence x ∈ j∈N Q j .
The next theorem justifies the definition of Calderón-Zygmund operators with standard kernels by the Cauchy principal value of the integral (2.13) and shows their boundedness in L p ω (Ω) and L ψ (Ω). Although the result is well-known, in standard literature it appears only in the setting Ω = R n . We prove it below for any domain Ω ⊂ R n .
Proof of Theorem 3.4. We are going to prove the theorem only for the operator T . The proof for T ( * ) follows the same reasoning. At first, let us assume that where Ω := {x ∈ R n | dist (x, Ω) ≤ 4 diam Ω}.
Then the kernel satisfies the conditions (2.15) and (2.16). Moreover, it satisfies (2.17) globally, i.e., with E = R n . Thus, K ∈ CZ(R n , κ 2 ). Therefore, by [7, Theorem 2], is defined for a.e. x ∈ R n and all g ∈ C ∞ 0 (R n ) and it admits a bounded extension from L 2 (R n ) to L 2 (R n ). It follows from the proofs and comments in [7], in particular p. 295 and a remark on p. 306, that the corresponding operator norm satisfies T L 2 →L 2 ≤ c κ 2 , where c > 0 is a fixed constant depending only on the dimension n. Since T ε g(x) = T ε (gη)(x) holds for all ε > 0, g ∈ C ∞ 0 (R n ) and a.e. x ∈ R n , we may draw the same conclusions about the operator .
In particular, we get T g(x) = T (gη)(x), and thus Let us next show that K is a standard kernel with respect to R n . Suppose that x, y, z ∈ R n satisfy (2.8). Then We will distinguish two cases. At first, assume that x, z ∈ Ω and dist (y, Ω) ≤ diam Ω. (3.12) Then obviously y ∈ Ω, and since K ∈ SK( Ω, κ 1 ), we have In the second and third condition we used the estimate which holds due to (3.11) and (3.12), and the property ∇η ∞ ≤ 2(diam Ω) −1 . Now suppose that (2.8) holds but (3.12) does not. Then one of the following situations occurs: (i) dist (y, Ω) > diam Ω, in which case y / ∈ supp η.
In the next step, let ψ be an N -function with ∆ 2 (ψ) < ∞ and ∆ 2 (ψ * ) < ∞. Due to (3.13) Theorem 3.1 provides the existence of a positive constant C ψ such that for all f ∈ C ∞ 0 (R n ). To obtain this estimate we have performed the extrapolation with respect to the family The obtained constant C ψ is independent of κ 1 , κ 2 and Ω. For any f ∈ C ∞ 0 (Ω) and x ∈ Ω we have Hence, (3.14) and (3.13) yield (3.6) and (3.7), respectively, for f ∈ C ∞ 0 (Ω). Since C ∞ 0 (Ω) is dense in L p ω (Ω) as well as in L ψ (Ω), the operator T admits corresponding bounded extensions such that (3.6) and (3.7) hold for any f ∈ L p ω (Ω) and f ∈ L ψ (Ω), respectively. So far, we have proved the theorem under the stronger assumption (3.9). To prove it in full generality, let Ω ⊂ R n be an arbitrary domain and K ∈ SK(Ω, κ 1 ) ∩ CZ(Ω, κ 2 ). By Proposition 3.3 there exists a sequence {Q j } j∈N of pairwise disjoint dyadic cubes such that Ω = j∈N Q j and Hence, for any j ∈ N, we have K ∈ SK( Q j , κ 1 ) ∩ CZ( Q j , κ 2 ), and using the first part of the proof we get for all ω ∈ A p and f ∈ L p ω (Ω). In here, it is important that C p p [ω] Ap does not depend on Q j . Similarly, we obtain (3.7) for any f ∈ L ψ (Ω). Note also that Remark 3.8 is justified.

Divergence equation
We proceed with proving the new estimate concerning the Bogovski solution of the divergence equation. We prove the result in the weighted-L p setting. The variant for Orlicz modulars will be then obtained as a corollary by extrapolation.
is satisfied in Q. (ii) For every p ∈ (1, ∞) there exists a nondecreasing function C p : holds for all ω ∈ A p and all f ∈ C ∞ 0,0 (Q). The function C p is independent of Q.
(iii) For every p ∈ (1, ∞) there exists a nondecreasing function C p : holds for all h > 0, ω ∈ A p and f ∈ C ∞ 0,0 (Q). The function C p is independent of Q. (iv) For every p ∈ (1, ∞) and ω ∈ A p there exists a continuous extension of B to L p ω,0 (Q) such that (4.2), (4.3) and (4.4) hold for all f ∈ L p ω,0 (Q). In here, f and Bf are extended by zero outside Q so that their difference quotients are defined a.e. in R n .
Corollary 4.5. Let Q ⊂ R n be an open cube. Let ψ be a Young function with ∆ 2 (ψ) < ∞ and ∆ 2 (ψ * ) < ∞. Then there exists a linear operator B : L ψ 0 (Q) → W 1,ψ 0 (Q) with the following properties: holds for all h > 0 and all f ∈ L ψ 0 (Q). In here, f and Bf are extended by zero outside Q so that their difference quotients are defined a.e. in R n . The constant C ψ is independent of Q.
Proof of Theorem 4.1. At first suppose that Q = − 1 2 , 1 2 n . Choose a fixed function ̺ ∈ From [16, Chapter III, Lemma 3.1] it follows that Bf ∈ W 1,∞ 0 (Q) and it satisfies (4.2) in Q. From now on, we will use the notation p.v.ˆE g(y)K(x, y) dy := lim ε→0ˆE ∩{|x−y|>ε} g(y)K(x, y) dy for functions g ∈ C ∞ 0 (R n ) and kernels K whenever the limit on the right-hand side exists. Part (ii) is proved in [23,13] without the monotone dependence of C p on [ω] Ap . Let us sketch the proof to show how the desired relation between C p and [ω] Ap is obtained. The technique used here is the same as in part (iii), where full details will be given. Let i, j ∈ {1, . . . , n}. Since Q is convex and it contains supp ̺, by [16, p. 119] we have for any f ∈ C ∞ 0,0 (Q). In here, and U ij is a kernel such that |U ij (x, y)| ≤ c ̺ 1,∞ |x − y| 1−n , where c is a positive constant. From [10, Lemma 6.1] we get that J ij ∈ SK(Q, κ)∩ CZ(Q, κ), where κ depends only on ̺ 1,∞ and n. By Theorem 3.4 there exists a positive nondecreasing function c p such that for all f ∈ C ∞ 0,0 (Q). Furthermore, we have (see [10, p. 218]) for f ∈ C ∞ 0,0 (Q) and x ∈ Q. Hence, by [18, Theorem 7.1.9(b)] there exists a positive constant C (possibly depending on p, n) such that for all f ∈ C ∞ 0,0 (Q). Obviously, we also have Therefore, by combining (4.11), (4.12) and (4.13), we have shown that there exists a positive nondecreasing function C p such that (4.3) holds for all f ∈ C ∞ 0,0 (Q). Recall however that so far we have assumed Q = − 1 2 , 1 2 n . The result for a general cube will be obtained by rescaling at the end of the proof of part (iii). We continue with part (iii). We are going to prove it only for the difference quotient d + h,k . The proof for d − h,k is fully analogous. We are still assuming that Q = − 1 2 , 1 In here, G ij is a kernel such that |G h ij (x, y)| ≤ c ̺ 1,∞ |x − y| 1−n for all x, y ∈ Q ∪ (Q − he k ), x = y. The positive constant c here depends only on n. We have also used (4.14) here. Furthermore, the kernel K h i,j is expressed as Let us proceed by estimating the partial derivatives of a 2 (x), still under the assumption h ∈ 0, 1 4 . Since supp ̺ ⊂ − 1 4 , 1 4 n holds, the function d + h,k ̺ is supported in Q. The cube Q is convex, thus we may again use the same calculation as in [16, p. 119] to obtain In here, H h ij is a kernel satisfying |H h ij (x, y)| ≤ C|x − y| 1−n for all x, y ∈ Q, x = y, with C depending only on n and ̺ 2,∞ . Furthermore, the kernel M h ij is, for x, y ∈ R n , x = y, defined as follows: . Observe that if x, y ∈ Q and ξ > n, then Hence, the integrals over (0, ∞) in the definition of M h ij may be replaced by integrals over (0, n).
From [10, Lemma 6.1] it follows that m h,1 ij is a Calderón-Zygmund kernel with respect to Q, with a constant independent of h. Here, notice that the function d + h,k ̺ plays the role of ̺ in [10], and d + h,k ̺ 1,∞ ≤ ̺ 2,∞ . In the next step, we shall verify that m h,1 ij is a standard kernel with respect to Q with a constant independent of h. Let x, y, z ∈ Q be such that x = y and |x − z| ≤ 1 2 |x − y|. It is easy to see that Using (3.11), we obtain and ≤ (1 + 3 · 2 n (n + 1)) |x − z| |x − y| n+1 .
(4.18) Moreover, we have Now we use (4.17), (4.19), Proposition 4.6 and the inequality Analogously, we obtain the following: Therefore, m h,1 ij is a standard kernel with respect to Q with a constant independent of h. Analogously, using (3.11), (4.18) and (4.19), we show that m h,2 ij is a standard kernel with respect to Q with a constant depending only on n and ̺ 3,∞ . Hence, M h ij is a standard kernel with respect to Q with a constant independent of h. By Theorem 3.4 there exists a positive nondecreasing function c p,2 such that for all h ∈ 0, 1 4 and all f ∈ C ∞ 0,0 (Q) we have The remaining parts of (4.16) are treated analogously as their counterparts in ∂a 1 j ∂x i . Altogether, it follows that there exists a positive nondecreasing function C p,2 such that holds for all h ∈ (0, 1 4 ) and f ∈ C ∞ 0,0 (Q). Using this result and (4.15), we obtain the existence of a positive nondecreasing function C p such that the following holds for all h ∈ 0, 1 4 and Notice also that λ = ℓ(λQ). Obviously, this result remains unchanged when the cube λQ is shifted, e.g., replaced by x 0 + λQ for any fixed x 0 ∈ R n . Hence, we have now proved (i)-(iii) in full generality. Part (iv) is proven by a standard approximation argument, using density of C ∞ 0,0 (Q) in L p ω,0 (Q) for any fixed p and ω. This completes the proof.

Interior Regularity
From now on, let Ω ⊂ R n be a bounded domain. We also assume that S has a (p, δ)-structure for some p ∈ (1, ∞) and δ ≥ 0. The properties of S and the standard theory of monotone operators imply in a standard way the existence of a unique u ∈ W 1,p 0,div (Ω), satisfyinĝ Ω S(Du) · Dw dx =ˆΩ f · w dx (5.1) for all w ∈ W 1,p 0,div (Ω), i.e., u is a weak solution of (1.1). Using Proposition 2.21, the properties of S and the Poincaré, Korn and Young inequalities, we obtain that this solution satisfies the a priori estimate γ 0 (p)ˆΩ ϕ(|∇u|) dx ≤ cˆΩ ϕ * (|f |) dx.
From here on, we denote by γ i (p), i = 0, 1, the constants in Definition 2.19 for a given p. Moreover, all constants may depend on the characteristics of S, diam(Ω), |Ω|, the space dimension n and on the John constants of Ω. The dependence on these quantities will not be mentioned explicitly anymore. However, dependence on other quantities will be specified.