Preprint
Article

This version is not peer-reviewed.

New Boundary Harnack Principles for Divergence Form Elliptic Equations with Right Hand Side

Submitted:

03 June 2025

Posted:

05 June 2025

You are already at the latest version

Abstract
We establish new boundary Harnack principles for divergence form equations with right hand side. In Lipschitz domains, we prove new boundary Harnack principles when the right hand side exhibits polynomial decay. Moreover, in Hölder domains with the Hölder exponent α∈1/2,1, we prove new boundary Harnack principles when the right hand side exponential decay.
Keywords: 
;  ;  ;  

1. Introduction

1.1. Background

In this paper, we study boundary Harnack principles for divergence form equations with right hand side. Boundary Harnack principles originally state that, if u and v are two positive harmonic functions in Ω B 1 , and they both vanish on Ω B 1 , then they are comparable near the boundary Ω .
For elliptic equations without right hand side in Lipschitz domains, we introduce some related research results. Kemper initially proved boundary Harnack principles for harmonic functions by applying the uniqueness of the positive singularity principle and the properties of Kernel functions in [24], while Caffarelli, Fabes, Mortola and Salsa first proved boundary Harnack principles for divergence form operators with bounded measurable coefficients in [8]. In [18], Fabes, Garofalo, Marin-Malave, and Salsa extensively discussed the non-divergence form operators by Fatou theorems. In [23], Jerison and Kenig extended boundary Harnack principles for divergence form operators in non-tangentially accessible domains. Bass and Burdzy proved boundary Harnack principles for non-divergence form elliptic operators with bounded measurable coefficients by applying probabilistic techniques in [5]. More recently, De Silva and Savin discovered a simple yet unified proof of divergence and non-divergence form operators with bounded measurable coefficients in their work published in [13].
For elliptic equations with right hand side in Lipschitz domains, we present some relevant results. In [1], Mark Allen and Shahgolian proved boundary Harnack principles of the Laplace operators by using harmonic functions in spherical polar coordinates under polynomial decaying conditions of the right hand side. In [31], Ros-Oton and Torres-Latorre already proved boundary Harnack principles for divergence form operators with continuous coefficients and non-divergence form operators with bounded measurable coefficients. In addition, Ros-Oton and Torres-Latorre gave a counterexample to show that boundary Harnack principles do not hold in R + 2 for divergence form operators with bounded measurable coefficients, even though the L q norm of the right hand side is controlled by a sufficiently small constant c 0 .
For elliptic equations without right hand side in Hölder domains, we introduce some related background and research results. Boundary Harnack principles for non-divergence form operators in Hölder domains with the exponent α 1 2 , 1 were first proven by Banuelos, Bass and Burdzy in [5] Bass1994. Fausto Ferrari first proved boundary Harnack principles for divergence form operators in Hölder domains with exponent α 0 , 1 in [19]. De Silva and Savin proved boundary Harnack principles for divergence and non-divergence form operators in Hölder domains with exponent α 1 2 , 1 by using a simpler method in [13]. Further, in [14], De Silva and Savin, by optimizing the previous approach, generalize this theory in Hölder domains with exponent α 0 , 1 .
Based on our current investigation and understanding, boundary Harnack principles for elliptic equations with right hand side in Hölder domains has not been fully studied. In [31], Ros-Oton and Torres-Latorre gave a counterexample to show that boundary Harnack principles do not hold in Lipschitz domains with a large slope, even though the right hand side is bounded and less than a very small constant c 0 . Obviously, this example further shows that boundary Harnack principles do not hold in Hölder domains, if without any decaying condition of the right hand side.
In this paper, we apply the positive extension theorem to obtain the decaying estimation of solutions of homogeneous divergence form equations in Lipschitz domains (Hölder domains). Because of the counterexamples mentioned above, by adding a polynomial(exponential) decaying condition to the right hand side, we get boundary Harnack principles in Lipschitz domains(Hölder domains with exponent α 1 2 , 1 ) for divergence form operators with bounded measurable coefficients. At present, for the caes of boundary Harnack principles in Hölder domains with exponent α < 1 2 , we still have not reached a corresponding conclusion.
Finally, the content of this paper is divided into four parts. The first part mainly introduces the main theorem of this paper. In the second part, we introduce some mathematical tools including the A-B-P estimation, the positive extension theorem, the decay theorem, the interior Harnack inequality and an upper bound estimation. The third part is to prove boundary Harnack principles in Lipschitz domains. The content of this part includes a decay lemma of the solution in a Lipschitz cone, a key iteration theorem and a key proposition. The last part is mainly to prove boundary Harnack principles in Hölder domains with exponent α 1 2 , 1 . The arrangement of the content and the idea of the argument in the last part are basically the same as in the third part.

1.2. The Main Notation

x = x , x n = x 1 , . . . , x n 1 , x n is a point in R n , n > 1 .
B ρ y = x R n : x y < ρ is an open ball in R n with centre y and radius ρ . B ρ y is an open ball in R n 1 with centre y and radius ρ .
Let h R + . The set C η = x R n : η x < x n < h is called the Lipschitz cone of height h with the Lipschitz constant η . The set E η = x R n : η 2 x α < x n < h is called the Hölder cone of height h with the Hölder constant η 2 and the range α 0 , 1 .
d x = d i s t x , Ω is the distance between x and Ω .
Given a function u, we write u ± = max ± u , 0 .
Ω is the Lebesgue measure of the domain Ω , and d i a m Ω is the diameter of the domain Ω .
e i = ( 0 , . . . 0 i 1 , 1 , 0 , . . . 0 n i ) R n .

1.3. Setting

In what follows, let L be a divergence form elliptic operator with bounded measurable coefficients, i.e.
L u = D i v A x u w i t h λ I A x Λ I ,
where 0 < λ Λ .
Definition 1.
Let L be a positive constant. Let ω : B 1 R be a Lipschitz function satisfied ω 0 = 0 and ω x ω y L x y , for any x and y in B 1 . We say Ω is a Lipschitz domain with the Lipschitz constant L if Ω = x , x n B 1 × R : x n > ω x .
Definition 2.
Let constants H > 0 and α ( 0 , 1 ] . Let ω : B 1 R be a Hölder function satisfied ω 0 = 0 and ω x ω y H x y α , for any x and y in B 1 . We say Ω is a Hölder domain with exponent α and the Hölder constant H if Ω = x , x n B 1 × R : x n > ω x .
Definition 3.
(De Giorgi super-class) Let E R n be an open set. We say u D G E , C if u W l o c 1 , 2 E satisfies
B r y u k 2 d x C R r 2 B R y u k 2 d x ,
for all k R , B r y B R y E , where C is a positive constant depending on λ, Λ and the dimension n.

1.4. Main Results

Theorem 1.
There exists a large positive constant q * and small positive constants η and c 0 such that the following holds.
Let Ω be a Lipschitz domain with the Lipschitz constant L < η 2 . Let u > 0 and v > 0 be weak solutions of
L u = f i n Ω B 1 u = 0 o n Ω B 1 a n d L v = g i n Ω B 1 v = 0 o n Ω B 1 ,
with f L B r y c 0 r q * 1 and g L B r y c 0 r q * 1 , for any r 0 , 1 , y Ω B 1 . Suppose further that v e n / 2 1 and u e n / 2 1 .
Then we have
u C v i n B 1 / 2 Ω ,
where C , c 0 , η and q * depend on λ, Λ and the dimension n.
With Theorem 1, it is straightforward to obtain the following Remark. If necessary, we can take a smaller constant c 0 .
Remark 1.
Let u > 0 and v > 0 satisfy (1.2) with f and g satisfying
f x c 0 d ( x ) q * 1 g x c 0 d ( x ) q * 1 o r f x c 0 Γ h ( x ) q * 1 g x c 0 Γ h ( x ) q * 1 ,
where Γ h ( x ) = x n ω ( x ) .
Then the boundary Harnack principle (1.3) also holds.
Theorem 2.
Let Ω be a Hölder domain with α 1 2 , 1 and the Hölder constant H < 4 e 4 3 . There exists a large positive constant C 2 and a small positive constant c 0 > 0 such that the following holds.
Let u > 0 and v > 0 be weak solutions of
L u = f i n Ω B 1 u = 0 o n Ω B 1 a n d L v = g i n Ω B 1 v = 0 o n Ω B 1 ,
with f L B r y c 0 r 2 exp 4 C 2 r 1 and g L B r y c 0 r 2 exp 4 C 2 r 1 , for any r 0 , 1 / 2 , y Ω B 1 . Suppose further that v e n / 2 1 and u e n / 2 1 .
Then we have
u C v i n B 1 / 2 Ω ,
where C , c 0 , C 2 depend on α , λ , Λ and the dimension n.

2. Preliminaries

In this section, we introduce five common theorems to prepare for later proofs.
The following theorem is called the A-B-P estimation. In the proof of our main theorem, we only need to consider the case of p = n .
Theorem 3.
([22, Theorem 8.16]) Assume that Ω R n is a bounded domain. Let u C Ω ¯ and satisfy L u f with f L p Ω , p > n / 2 .
Then
sup Ω u sup Ω u + C f L p Ω ,
where C depend on d i a m Ω , p, λ, Λ and the dimension n.
The following theorem is called the positive expansion. This theorem plays a crucial role in the section 3.
Theorem 4.
([15, Theorem 9.1]) Let constants M > 0 , ϵ 0 , 1 and r 2 ( 1 2 , 1 ) . Let E R n be an open set. For any non-negative function u D G E , C . Assume that u M in B ϵ R y and B 2 R E .
Then we have
u σ ϵ q 1 M i n B R y ,
where q 1 = C ϵ r 2 2 n 1 1 r 2 2 n + 3 and σ = exp C ϵ r 2 2 n 1 1 r 2 2 n + 2 with C is a positive constant depending on λ , Λ and the dimension n.
Moreover, if ϵ e 1 , we have
u ϵ q M i n B R y ,
where the positive constant q : = 2 q 1 .
The following theorem is the decay theorem of the weak subsolution in the ball.
Theorem 5.
([31, Corollary 2.5]) Let R 0 , 1 and u satisfy
L u f i n B R u 1 i n B R ,
where f L n B R . Suppose that there exists θ 0 , 1 such that u 0 B R θ B R > 0 and f L n B R ξ θ .
Then we have
sup B R / 2 u 1 γ θ ,
where ξ θ > 0 and γ θ 0 , 1 depend on θ, λ, Λ and the dimension n.
The following theorem is the interior Harnack inequality.
Theorem 6.
([31, Lemma 3.6 and Lemma 3.7]) Let Ω be a Lipschitz domain with the Lipschitz constant L < 1 / 16 . Let u > 0 be a weak solution of
L u = f i n Ω B 1 u = 0 o n Ω B 1 ,
where f L n B 1 . We define A : = x Ω B 1 δ ¯ : Γ h x δ , where Γ h ( x ) = x n ω ( x ) , δ 0 , 1 / 3 .
Then we have
sup A u C inf A u + f L n B 1 ,
where C depend on δ, λ, Λ and the dimension n.
Moreover, if u e n / 2 1 , then there exist positive constants p and C p depending on λ, Λ, f L n B 1 and the dimension n, such that
u L p B 1 C p .
The following proposition is an upper bound estimate of the weak solution.
Proposition 1.
([31, Proposition 3.8]) Let Ω be a Lipschitz domain and u be a weak solution of
L u = f i n Ω B 1 u = 0 o n Ω B 1 ,
where f L n B 1 .
If u L p Ω B 1 for any p > 0 , then
sup B r u M n , p , λ , Λ , r u L p B 1 + f L n B 1 ,
where r 0 , 1 .

3. Proof of Theorem 1

The following lemma is the decay estimate of the weak supersolution in the Lipschitz cone. This lemma plays a crucial role in the proof of the following theorem.
Lemma 1.
There exists a small constant η > 0 and a large positive constant q * depending on λ, Λ and the dimension n, such that the following holds.
Let u satisfy
L u 0 i n C η u 1 i n x n 1 C η u 0 o n C η ,
where the cone C η = x R n : η x < x n < 3 .
Then for any t ( 0 , 1 3 ) , we have
u t e n t q * .
Preprints 162148 i001
Proof. 
Let ϵ = 1 3 and r 2 = 6 7 in Theorem 4. We take a i = 4 3 1 3 i 1 and construct two series of balls and denote them as B 1 3 i 0 , a i , B 1 3 i 1 0 , a i for all i N + , where 0 , a i represents a point in R n with x = 0 and x n = a i (see figure(a)). In addition, we know that i = 1 2 3 i = 1 .
Obviously, u x 1 in B 1 3 0 , a 1 and by Theorem 4, we have u x 1 3 q in B 1 0 , a 1 . Thus, u x 1 3 q in B 1 3 2 0 , a 2 and by Theorem 4 again, we have u x 1 3 2 q in B 1 3 0 , a 2 . Repeating the above process, we have u x 1 3 k + 1 q in B 1 3 k 0 , a k + 1 , for all k N .
For any t 0 , 1 3 , taking k 0 N + , such that 1 2 i = 1 k 0 + 1 1 3 i t 1 2 i = 1 k 0 1 3 i , then we have t e n B 1 3 k 0 0 , a k 0 + 1 and 1 3 k 0 + 1 t 1 3 k 0 . Hence, we have u t e n 1 3 k 0 + 1 q . Then k 0 + 1 k 0 2 implies
1 3 k 0 + 1 q 1 3 k 0 2 q t 2 q .
Taking q * : = 2 q , then we have
u t e n t 2 q = t q * , f o r a l l t 0 , 1 3 .
Finally, if we take η = 1 7 n 1 , then all the balls will be located in cone C η . □
The following theorem shows that, by constructing an iteration, the solution of L u = f is positive at some distance from the boundary.
Theorem 7.
Let Ω be a Lipschitz domain with the Lipschitz constant L < η 2 . The positive constants η and q * are determined by Lemma 1. Let u satisfy
L u = f i n Ω B 1 u = 0 o n Ω u 1 i n Ω x B 1 : Γ h x > δ u ϵ i n Ω B 1 , .
where f L B r y c 0 r q * 1 , for any r 0 , 1 , y Ω B 1 .
Then we have
u ρ q * + 1 i n Ω x B ρ : Γ h x > ρ δ u ρ q * + 1 ϵ i n Ω B ρ ,
for some small enough constants ρ, ϵ, δ, c 0 0 , 1 depending on λ, Λ and the dimension n.
Proof. 
Now that, we prove the first inequality in (3.3). For any point x 0 , x 0 n x B ρ : ρ δ < Γ h x < δ , we write x 0 ˜ : = x 0 , ω ( x 0 ) and the cone
C ρ , η : = x 0 ˜ + ρ C η = x R n : η x x 0 < x n ω x 0 < 3 ρ
and the positive upper cone
C ρ , η + : = C ρ , η x n > ω x 0 + ρ ,
where C η = x R n : η x < x n < 3 .
Let ρ = 3 δ , R 0 , ρ , ν R : = η 2 R and x ρ = ρ η + R (see Figure 2). We have
η 2 x ρ + δ = η 2 ρ η + R + ρ 3 = 5 6 ρ + ν R .
Obviously, ν R + ρ > ν R + 5 6 ρ implies that C ρ , η + Ω x B 1 : Γ h x > δ . Now we take δ small enough, then we have C ρ , η Ω B 1 . Because x 0 , x 0 n x B ρ : ρ δ < Γ h x < δ , we have
ρ δ < x 0 n ω ( x 0 ) < δ a n d t h e n ρ 3 < x 0 n ω ( x 0 ) ρ < 1 3 .
We define
u ˜ x : = u x 0 ˜ + ρ x + ϵ a n d Φ x : = L u ˜ x = ρ 2 f x 0 ˜ + ρ x .
Let u ˜ = u 1 + u 2 satisfy
L u 1 = 0 i n C η u 1 = u ˜ o n C η a n d L u 2 = Φ i n C η u 2 = 0 o n C η .
By Theorem 3, we have
u 2 L C η C Φ L n C η C ρ 2 C η f x 0 ˜ + ρ x n d x 1 n = C ρ C ρ , η f y n d y 1 n C ρ f L B ρ x 0 ˜ B ρ x 0 ˜ 1 n C C n c 0 ρ q * + 1 .
On the other hand, we apply Lemma 1 to u 1 to result that u 1 t e n t q * , for any t ( 0 , 1 3 ) . We have
u x 0 , ω x 0 + ρ t e n t q * C C n c 0 ρ q * + 1 ϵ , f o r a n y t ( 0 , 1 3 ) .
Therefore, let ρ 1 3 q * + 1 , ϵ ρ q * + 1 , c 0 1 C C n , then
u x 0 u x 0 , ω x 0 + ρ x 0 n ω ( x 0 ) ρ e n x 0 n ω ( x 0 ) ρ q * C C n c 0 ρ q * + 1 ϵ ρ 3 q * C C n c 0 ρ q * + 1 ρ q * + 1 ρ q * + 1 .
Next, we prove the second inequality in (3.3). Let x 0 B 1 2 δ and Γ h x δ (recalling δ = ρ / 3 ). We define v x : = u x ϵ 0 . Extending v by 0 below Ω , then we have L v f + ϵ in B 2 δ x 0 . Let C z be the downwards cone with slope η and vertex in z : = x 0 , ω ( x 0 ) . Then we have v = 0 in C z B 2 δ x 0 and C z B 2 δ x 0 C η B 2 δ x 0 . Taking c 0 small enough, we have
B 2 δ x 0 f + ϵ n 1 n B 3 δ z f + ϵ n 1 n 1 ϵ f L n B 3 δ z 3 δ C n ϵ f L B 3 δ z c 0 C n ( 3 δ ) q * ϵ ξ η ,
where ξ η come from Theorem 5. By Theorem 5, we have v 1 γ in B δ x 0 .
Now we define v j x : = v 1 γ j x in B δ x 0 with x 0 B 1 2 j δ and Γ h x 0 δ . By mathematical induction, it is easy to show that v 1 γ j in B 1 2 j δ , i.e. u 1 γ j ϵ in B 1 2 j δ . We take j 0 N + satisfied 1 2 j 0 + 1 δ ρ 1 2 j 0 δ (recalling ρ = 3 δ ). Hence, we have
1 2 3 j 0 + 5 3 ρ 1 2 3 j 0 + 1 .
We take ρ small enough, such that
1 γ j 0 1 2 3 j 0 + 5 3 q * + 1 ,
then 1 γ j 0 ρ q * + 1 .
As same as (3.4), we have
B 2 δ x 0 f + 1 γ j 0 ϵ n 1 n 1 1 γ j 0 ϵ f L n B 3 δ z c 0 C n ( 3 δ ) q * 1 γ j 0 ϵ ξ η .
By 1 2 j 0 δ ρ with 1 γ j 0 ρ q * + 1 and u 1 γ j 0 ϵ in B 1 2 j 0 δ , then we have u ρ q * + 1 ϵ in B ρ . □
Now, we iterate Theorem 7 to obtain the following proposition.
Proposition 2.
Let Ω be a Lipschitz domain with the Lipschitz constant L < η 2 . The positive constants η and q * are determined by Lemma 1. Let u be a weak solution of (3.2) with f L B r y c 0 r q * 1 , for any r 0 , 1 , y Ω B 1 .
Then we have
u > 0 i n Ω B 2 / 3 u t e n t q * + 1 f o r a l l t 0 , 1 .
Proof. 
Let u 0 x : = u x and f 0 x : = f x . We define
u j + 1 x : = u j ρ x ρ q * + 1 , f j + 1 x : = f j ρ x ρ q * 1 .
where x Ω j + 1 and ρ x ρ Ω j + 1 : = Ω j . Obviously, L j L j + 1 , where L j is the Lipschitz constant of Ω j .
First, we already have f 0 L B r y c 0 r q * 1 , for any r 0 , 1 , y Ω B 1 . By mathematical induction, we have
L u j = f j i n Ω j B 1 u j = 0 o n Ω j w i t h u j 1 i n Ω j x B 1 : Γ h x > δ u j ϵ i n Ω j B 1 ,
where f j L B r y c 0 r q * 1 , for any r 0 , 1 , y Ω B 1 . By Theorem 7, we have
u j ρ q * + 1 i n Ω j x B ρ : Γ h x > ρ δ u j ρ q * + 1 ϵ i n Ω j B ρ ,
then
u j + 1 1 i n Ω j + 1 x B 1 : Γ h x > δ u j + 1 ϵ i n Ω j + 1 B 1 .
Moreover, we have
f j + 1 L B r y = 1 ρ q * 1 f j ρ x L B r y c 0 ρ q * 1 ρ r q * 1 = c 0 r q * 1 .
So for all j N , we have u j t e n 1 , for every t δ , 1 .
For any t 0 , 1 , there exists j 0 N such that δ < t ρ j 0 < 1 . Thus, we have t < ρ j 0 and
u t e n = u 0 t e n = ρ q * + 1 u 1 t e n ρ = . . . = ρ j 0 q * + 1 u j 0 t e n ρ j 0 t q * + 1 .
Finally, taking ϵ , ρ and c 0 smaller, then u x 0 + x 3 also satisfy (5) for any x 0 B 2 / 3 Ω . By translating the coordinates, we get
u x 0 + t e n 3 > 0 , f o r a l l t 0 , 1 .
Then this implies
u > 0 i n Ω B 2 / 3 .

3.1. Proof of Theorem 1

Proof. 
Let Ω be a Lipschitz domain with the Lipschitz constant L < η / 2 . The positive constants η and q * are determined by Lemma 1. By Theorem 6 and Proposition 1, we have u M in B 3 / 4 . We consider v in the set E : = x B ¯ 3 / 4 : Γ h ( x ) 3 δ / 4 . Obviously, E x B ¯ 1 3 δ / 4 : Γ h x 3 δ / 4 . By Theorem 6 and v e n / 2 1 , we have
v C 1 f L n B 1 C 1 C n f L B 1 C 1 C n c 0 2 C 1 : = m i n E ,
where the constant c 0 is small enough such that C 1 C n c 0 2 C 1 .
We define
φ x : = 1 + ϵ m v x ϵ M u x
where ϵ > 0 is chosen later. Once we prove that φ > 0 in B 1 / 2 , by choosing C = M 1 + ϵ / m ϵ , then C v u > 0 follows.
Obviously, φ v / m 1 in E, and φ ϵ in B 3 / 4 . We have
L φ L B r y 1 + ϵ m L v L B r y + ϵ M L u L B r y 1 + ϵ m + ϵ M c 0 r q * 1 , f o r a n y r 0 , 1 , y Ω .
Let φ ˜ x : = φ 3 x / 4 . Then φ ˜ 1 in Ω x B 1 : Γ h x δ and φ ˜ ϵ in Ω B 1 . Taking ϵ , c 0 > 0 small enough and applying Proposition 2, then we get φ ˜ > 0 in B 2 / 3 . Thus, φ > 0 in B 1 / 2 . □

4. Proof of Theorem 2

A simple lemma is given here to prepare for the following iteration of the ball in the Hölder cone.
Lemma 2.
Let constants ϵ 0 , 1 , η > 0 and α 1 2 , 1 . Then there exists a decreasing and positive sequence R k k = 1 satisfied R k 0 as k + and
ϵ R k α = 1 η j = 1 k R j ,
for any k N + . Moreover, there exists a constant C 1 > 0 such that
ϵ R k 1 2 C 1 ϵ 2 η k , f o r a n y k N + .
Proof. 
We apply mathematical induction to prove (4.1). First, we show that there exists R 1 satisfying ϵ R 1 α = 1 η R 1 . Let φ t : = 1 η t ϵ t α and note that φ 0 = 1 > 0 , φ 1 η = ϵ 1 η α < 0 with φ t < 0 in 0 , 1 η . By the zero point theorem, there exists R 1 such that 1 η R 1 ϵ R 1 α = 0 . By mathematical induction, we suppose that ϵ R k α = 1 η j = 1 k R j .
Next, we will show (4.1) for k + 1 . Denoting a k = 1 η j = 1 k 1 R j , then we have 1 η j = 1 k R j = 1 η j = 1 k 1 R j η R k = a k η R k . We define
Φ t : = a k η R k η t ϵ t α = ϵ R k α ϵ t α η t .
Then we have Φ 0 = ϵ R k α and Φ R k = η R k α < 0 with Φ t < 0 in 0 , R k . By the zero point theorem, there exists R k + 1 satisfied 0 < R k + 1 < R k and a k η R k η R k + 1 ϵ R k + 1 α = 0 , i.e. ϵ R k + 1 α = a k η R k + R k + 1 .
We deserve that 0 < R k + 1 < R k and
ϵ R k ϵ R k + 1 = η R k + 1 .
We have R = lim k R k = 0 . By (4.3) and letting k , we have R = 0 . By (4.1), we get η j = 1 R j = 1 .
Finally, we show (4.2). Obviously, if α = 1 2 holds true, then α [ 1 2 , 1 ] holds even more. So, we consider α = 1 2 . We have ϵ R k 1 2 ϵ R k + 1 1 2 = η R k + 1 , i.e. R k 1 2 R k + 1 1 2 = η ϵ R k + 1 : = C 0 R k + 1 with C 0 depending only on η and ϵ. Now, we have R k = R k + 1 + 2 C 0 R k + 1 3 2 + C 0 2 R k + 1 2 . For k large enough, we have ( R k R k + 1 ) / ( 2 C 0 R k + 1 3 2 ) = 1 + C 0 2 R k + 1 1 2 1 . Hence, we get R k R k + 1 2 C 0 R k + 1 3 2 and
lim k k 2 R k = lim k k 2 / ( 1 / R k ) = lim k [ ( k + 1 ) 2 k 2 ] / ( 1 / R k + 1 1 / R k ) = lim k ( 2 k + 1 ) R k + 1 2 / ( 2 C 0 R k + 1 3 2 ) = lim k ( 2 k + 1 ) ( 2 k 1 ) / [ 2 C 0 ( 1 / R k + 1 1 / R k ) ] = lim k R k R k + 1 / [ C 0 ( R k R k + 1 ) ] = lim k R k + 1 / ( C 0 2 R k + 1 ) = 1 / C 0 2 .
Hence, we have R k 1 2 1 / ( C 0 k ) = ϵ / ( η k ) .
By ϵ R k 1 2 ϵ 2 / ( η k ) , there exists a constant C 1 such that ( C 1 ϵ 2 ) / ( η k ) ϵ R k 1 2 , for any k N + .
The following lemma is the decay lemma of the weak supersolution in the Hölder cone. This lemma plays a crucial role in the proof of the following theorem.
Lemma 3.
Let α 1 2 , 1 . There exist constants C 2 > 0 and η > 0 depending on α , λ , Λ and the dimension n, such that the following holds.
Let u be a solution of
L u 0 i n E η u 1 i n x n 1 E η u 0 o n E η ,
where the E η = x R n : η 2 x α < x n < 3 .
Then for any t ( 0 , 1 3 ) , we have u ( t e n ) e C 2 t .
Preprints 162148 i002
Proof. 
We take R 1 = 1 3 and η = 1 e 1 3 3 α . By Lemma 2, we have a decreasing and positive sequence R k satisfied R k 0 as k + and η R k α = 1 ( e 1 ) j = 1 k R j for any k N + . Moreover, there exists a constant C 1 such that C 1 ϵ 2 η k ϵ R k 1 2 , for any k N + .
We constructed two series of balls and denoted them as B R i 0 , a i and B e R i 0 , a i , for all i N + (see figure (c)). Now, there exists k 0 N + such that 1 ( e 1 ) j = 1 k 0 + 1 R j t 1 ( e 1 ) j = 1 k 0 R j , i.e. η R k 0 + 1 α t η R k 0 α and t e n B 2 e R k 0 + 1 ( x k 0 + 1 ) .
As u ( x ) 1 in B R 1 ( x 1 ) , Theorem 4 implies that u ( x ) e q in B e R 1 ( x 1 ) . Hence, u ( τ e n ) e q in B R 2 ( x 2 ) . By Theorem 4 again, we have u ( x ) e 2 q in B 2 e R 2 ( x 2 ) . Now keep doing it this way, we get u ( x ) e k q in B e R k ( x k ) for any k N + . Hence, u ( t e n ) e ( k 0 + 1 ) q = e q ( e q ) k 0 and η R k 0 1 2 C 1 η 2 ( e 1 ) k 0 . As η R k 0 1 2 η R k 0 α t , we have t C 1 η 2 ( e 1 ) k 0 , i.e. k 0 C 1 η 2 ( e 1 ) t . Because t < 1 and q > 1 , we have e q ( e q ) k 0 ( e q ) ( e q ) C 1 η 2 ( e 1 ) t e q C 1 η 2 + q t : = e C 2 t for any t ( 0 , 1 3 ) . Thus, we have u ( t e n ) e C 2 t , for any t ( 0 , 1 3 ) .
The following theorem works similarly to Theorem 7.
Theorem 8.
Let Ω be a Hölder domain with α [ 1 / 2 , 1 ) and the Hölder constant H < 4 e 4 3 . Let u satisfy
L u = f i n Ω B 1 u = 0 o n Ω u 1 i n Ω x B 1 : Γ h x > δ u ϵ i n Ω B 1 ,
. f L B r y c 0 r 2 exp 4 C 2 r 1 for any r 0 , 1 / 2 , y Ω B 1 .Then we have
u exp 4 C 2 ρ 1 i n Ω x B ρ : Γ h x > ρ δ u exp 4 C 2 ρ 1 ϵ i n Ω B ρ ,
where the positive constant C 2 is determined by Lemma 3 and some small enough constants ρ, ϵ, δ, c 0 0 , 1 depend on α , λ , Λ and the dimension n.
Proof. 
Now that, we prove the first inequality in (4.6). For any x 0 , x 0 n x B ρ : ρ δ < Γ h x < δ , we write x 0 ˜ : = x 0 , ω ( x 0 ) and the Hölder cone
E ρ , η = x 0 ˜ + ρ E η = x R n : η 2 x x 0 α < x n ω x 0 < 3 ρ
and the positive upper Hölder cone
E ρ + = E ρ , η x n > ω x 0 + ρ ,
where E η = x R n : η 2 x α < x n < 3 .
We take ρ = 3 δ , it’s easy to get E ρ + Ω x B 1 : Γ h x > δ . If δ small enough, we have E ρ , η Ω B 1 . Because x 0 = x 0 , x 0 n in x B ρ : ρ δ < Γ h x < δ , we have
ρ δ < x 0 n ω ( x 0 ) < δ a n d t h e n ρ 3 < x 0 n ω ( x 0 ) ρ < 1 3 .
We define
u ˜ x : = u x 0 ˜ + ρ x + ϵ a n d Φ x : = L u ˜ x = ρ 2 f x 0 ˜ + ρ x .
Let u ˜ = u 1 + u 2 satisfy
L u 1 = 0 i n E η u 1 = u ˜ o n E η a n d L u 2 = Φ i n E η u 2 = 0 o n E η .
By Theorem 3, we have
u 2 L E η C Φ L n E η C ρ 2 E η f x 0 ˜ + ρ x n d x 1 n = C ρ E ρ , η f y n d y 1 n C ρ f L B ρ x 0 ˜ B ρ x 0 ˜ 1 n C C n c 0 exp 4 C 2 ρ 1 .
On the other hand, we apply Lemma 3 to u 1 to result that u 1 t e n exp C 2 t 1 , for any t ( 0 , 1 3 ) . Then we get
u x 0 , ω x 0 + ρ t e n exp C 2 t 1 C C n c 0 exp 4 C 2 ρ 1 ϵ , f o r a n y t ( 0 , 1 3 ) .
Therefore,
u x 0 u x 0 , ω x 0 + ρ x 0 n ω ( x 0 ) ρ e n exp C 2 x 0 n ω ( x 0 ) ρ 1 C C n c 0 exp 4 C 2 ρ 1 ϵ exp C 2 ρ 3 1 C C n c 0 exp 4 C 2 ρ 1 ϵ 3 exp 4 C 2 ρ 1 C C n c 0 exp 4 C 2 ρ 1 ϵ exp 4 C 2 ρ 1 .
Next, we prove the second inequality in (4.6). Let x 0 B 1 2 δ and Γ h x δ . We define v x : = u x ϵ 0 and extend v by 0 below Ω . Thus, L v f + ϵ in B 2 δ x 0 . We write z : = x 0 , ω ( x 0 ) . Let E z be the downwards Hölder cone with the Hölder constant η and vertex in z. Then we have v = 0 in E z B 2 δ x 0 and E z B 2 δ x 0 C η B 2 δ x 0 . Hence, we have
B 2 δ x 0 f + ϵ n 1 n B 3 δ z f + ϵ n 1 n 3 δ C n ϵ f L B 3 δ z c 0 C n ρ 2 exp 4 C 2 ρ 1 ϵ ξ η ,
where ξ η come from Theorem 5 and make c 0 small enough such that the last inequality holds. By Theorem 5, we have v 1 γ in B δ x 0 .
Now we define v j x : = v 1 γ j x in B δ x 0 with x 0 B 1 2 j δ and Γ h x 0 δ . By mathematical induction, it is easy to show that v 1 γ j in B 1 2 j δ , i.e. u 1 γ j ϵ in B 1 2 j δ . We take j 0 N + satisfied 1 2 j 0 + 1 δ ρ 1 2 j 0 δ , with ρ = 3 δ . Hence, we have
1 2 3 j 0 + 5 3 ρ 1 2 3 j 0 + 1 .
We take ρ small enough, such that
exp j 0 ln 1 γ exp 4 C 2 2 3 j 0 + 5 3 a n d t h e n 1 γ j 0 exp 4 C 2 2 3 j 0 + 5 3 .
Hence, we have 1 γ j 0 exp 4 C 2 ρ 1 .
Finally, we have
B 2 δ x 0 f + 1 γ j 0 ϵ n 1 n B 3 δ z f + 1 γ j 0 ϵ n 1 n 1 1 γ j 0 ϵ f L n B 3 δ z 3 δ C n 1 γ j 0 ϵ f L B 3 δ z c 0 C n ρ 2 exp 4 C 2 ρ 1 1 γ j 0 ϵ ξ η ,
where ξ η come from Theorem 5 and make c 0 small enough such that the last inequality holds.
By 1 2 j 0 δ ρ with 1 γ j 0 exp 4 C 2 ρ 1 , and u 1 γ j 0 ϵ in B 1 2 j 0 δ .
Hence, we have u exp 4 C 2 ρ 1 ϵ in B ρ . □
Now, by iterating Theorem 8, we get the following proposition.
Proposition 3.
Let Ω be a Hölder domain with α [ 1 / 2 , 1 ) and the Hölder constant H < 4 e 4 3 . The positive constant C 2 is determined by Lemma 3. Let u be a weak solution of (4.5) and f L B r y c 0 r 2 exp 4 C 2 r 1 , for any r 0 , 1 / 2 , y Ω B 1 .
Then we have
u > 0 i n Ω B 2 / 3 u t e n exp 4 C 2 t 1 f o r a l l t 0 , 1 .
Proof. 
Let u 0 x : = u x and f 0 x : = f x . We define
u j + 1 x : = u j ρ x exp 4 C 2 ρ 1 , f j + 1 x : = ρ 2 f j ρ x exp 4 C 2 ρ 1
with x Ω j + 1 and ρ x ρ Ω j + 1 : = Ω j , where the Hölder constant of Ω j + 1 is the same as or smaller than Ω j .
First we already have f 0 L B r y c 0 r 2 exp 4 C 2 r 1 , for any r 0 , 1 / 2 , y Ω B 1 . By mathematical induction, we have
L u j = f j i n Ω j B 1 u j = 0 o n Ω j w i t h u j 1 i n Ω j x B 1 : Γ h x > δ u j ϵ i n Ω j B 1 ,
where f j L B r y c 0 r 2 exp 4 C 2 r 1 , for any r 0 , 1 / 2 , y Ω B 1 . By Theorem 8, we have
u j exp 4 C 2 ρ 1 i n Ω j x B ρ : Γ h x > ρ δ u j exp 4 C 2 ρ 1 ϵ i n Ω j B ρ .
So, we have
u j + 1 1 i n Ω j + 1 x B 1 : Γ h x > δ u j + 1 ϵ i n Ω j + 1 B 1 .
Moreover, we have
f j + 1 L B r y = ρ 2 exp 4 C 2 ρ 1 f j ρ x L B r y c 0 ρ 2 exp 4 C 2 ρ 1 ρ r 2 exp 4 C 2 ( ρ r ) 1 c 0 r 2 exp 4 C 2 r 1 .
Then for every j N , we have u j t e n 1 , for any t δ , 1 . For any t 0 , 1 , there exists j 0 N such that δ < t ρ j 0 < 1 . Thus, t < ρ j 0 < ρ j 0 and
u t e n = u 0 t e n = exp 4 C 2 ρ 1 u 1 t e n ρ = . . . = exp 4 j 0 C 2 ρ 1 u j 0 t e n ρ j 0 exp 4 C 2 t 1 .
Finally, taking ϵ , δ and c 0 smaller, then u x 0 + x 3 also satisfy (4.5) for any x 0 B 2 / 3 Ω . By translating the coordinates, we get
u x 0 + t e n 3 > 0 , f o r a l l t 0 , 1 .
Then this implies
u > 0 i n Ω B 2 / 3 .

4.1. Proof of Theorem 2

Proof. 
By Theorem 6 and Proposition 1, we have u M in B 3 / 4 . We consider the function v in the set E : = x B ¯ 3 / 4 : Γ h ( x ) 3 δ / 4 . Obviously, E x B ¯ 1 3 δ / 4 : Γ h x 3 δ / 4 . Thus, by Theorem 6 and v e n / 2 1 , we have
v C 1 f L n B 1 C 1 C n f L B 1 C 1 C n c 0 2 C 1 : = m i n E ,
where the constant c 0 is small enough such that C 1 C n c 0 2 C 1 .
We define
φ x : = 1 + ϵ m v x ϵ M u x ,
where ϵ > 0 to be chosen later. We will prove that φ > 0 in B 1 / 2 and therefore, choosing C = M 1 + ϵ / m ϵ , we have C v u > 0 .
Obviously, φ v / m 1 in E, and φ ϵ in B 3 / 4 . We have
L φ L B r y 1 + ϵ m L v L B r y + ϵ M L u L B r y 1 + ϵ m + ϵ M c 0 r 2 exp 4 C 2 r 1 , f o r a n y r 0 , 1 / 2 , y Ω .
Where the positive constant C 2 is determined by Lemma 3.
Let φ ˜ x : = φ 3 x / 4 . Then φ ˜ 1 in Ω x B 1 : Γ h x δ and φ ˜ ϵ in Ω B 1 . Taking small enough ϵ , c 0 > 0 to apply Proposition 3. We get φ ˜ > 0 in B 2 / 3 , thus φ > 0 in B 1 / 2 . □

Author Contributions

Writing—original draft preparation, X.L. and Y.C.; writing—review and editing, X.L. and Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (Grant No.12471464).

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Allen, M.; Shahgholian, H. A new boundary Harnack principle (equations with right hand side). Arch. Ration. Mech. Anal 2019, 234, 1413–1444. [Google Scholar] [CrossRef]
  2. Ancona, A. Principe de Harnack a la frontiere et theoreme de Fatou pour un operateur elliptique dons un domaine lipschitzien. Ann. Inst 1978, Fourier 28, 169–213. [Google Scholar]
  3. Banuelos, R.; Bass, R.F.; Burdzy, K. Hölder domains and the boundary Harnack principle. Duke Math. J 1991, 64, 195–200. [Google Scholar] [CrossRef]
  4. Bass, R.F.; Burdzy, K. A boundary Harnack principle in twisted Hölder domains. Ann. Math 1991, 134, 253–276. [Google Scholar] [CrossRef]
  5. Bass, R.F.; Burdzy, K. Further regularity for the Signorini problem. J. Lond. Math. Soc 1994, 50, 157–169. [Google Scholar] [CrossRef]
  6. Caffarelli, L.A. Further regularity for the Signorini problem. Commun. Partial Differ. Equ 1979, 4, 1067–1075. [Google Scholar] [CrossRef]
  7. Caffarelli, L.A.; Fabes, E.; Mortola, S.; Salsa, S. Boundary behavior of non-negative solutions of elliptic operators in divergence form. Indiana Univ. Math. J 1981, 49, 621–640. [Google Scholar] [CrossRef]
  8. Caffarelli, L.A.; Crandall, M.; Kocan, A.; Swiech. On viscosity solutions of fully nonlinear equations with measurable ingredients. Commun. Pure Appl. Math 1996, 49, 365–398. [Google Scholar] [CrossRef]
  9. Caffarelli, L.A.; Cabré, X. The boundary Harnack principle for non-divergence form elliptic operators. Amer. Math. Soc. Colloq. Publ 1995, 43, 1067–1075. [Google Scholar]
  10. Dahlberg, B. On estimates of harmonic measure. Arch. Ration. Mech. Anal 1977, 65, 272–288. [Google Scholar] [CrossRef]
  11. De Giorgi, E. Sulla differenziabilità e l’analiticità delle estremali degli integrali multipli regolari (Italian). Mem. Accad. Sci. Torino, Cl. Sci. Fis. Mat. Nat 1957, 3, 25–43. [Google Scholar]
  12. De Silva, D.; Savin, O. Boundary Harnack estimates in slit domains and applications to thin free boundary problems. Rev. Mat. Iberoam 2016, 32, 891–912. [Google Scholar] [CrossRef]
  13. De Silva, D.; Savin, O. A short proof of boundary Harnack inequality. C J. Differ. Equ 2020, 269, 2419–2429. [Google Scholar] [CrossRef]
  14. De Silva, D.; Savin, O. On the boundary Harnack principle in Hölder domains. Mathematics in Engineering 2022, 4, 1–12. [Google Scholar] [CrossRef]
  15. DiBenedetto, E. Partial Differential Equations: Second Edition. Birkhäuser Basel, 2010. [Google Scholar]
  16. DiBenedetto, E. Some properties of De Giorgi classes. Rend. Istit. Mat. Univ. Trieste 2016, 48, 245–260. [Google Scholar]
  17. DiBenedetto, E.; Trudinger, N. Harnack inequalities for quasiminima of variational integrals. Ann. Inst. H. Poincaré Anal 1984, 1, 295–308. [Google Scholar]
  18. Fabes, E.; Garofalo, N.; Marin-Malave, S.; Salsa, S. Fatou theorems for some nonlinear elliptic equions. Rev. Mat. Iberoam 1988, 4, 227–252. [Google Scholar] [CrossRef]
  19. Ferrari, F. On Boundary Behavior of Harmonic Functions in Hölder Domains. J. Fourier. Anal. Appl 1998, 4, 447–461. [Google Scholar] [CrossRef]
  20. Giaquinta, M.; Giusti, E. On the regularity of the minima of variational integrals. Acta Math 1982, 148, 31–46. [Google Scholar] [CrossRef]
  21. Giaquinta, M.; Giusti, E. Quasiminima. Ann. Inst. H. Poincaré Anal.Non Linéaire 1984, 1, 79–107. [Google Scholar] [CrossRef]
  22. Gilbarg, D.; Trudinger, N.S. Elliptic Partial Differential Equations of Second Order, 1st ed.Fundamental Principles of Mathematical Sciences, 1998. [Google Scholar]
  23. Jerison, D.S.; Kenig, C.E. Boundary behavior of harmonic functions in non-tangentially accessible domains. Adv. Math 1982, 42, 80–147. [Google Scholar] [CrossRef]
  24. Kemper, J.T. A boundary Harnack principle for Lipschitz domains and the principle of positive singularities. Commun. Pure Appl. Math 1972, 25, 247–255. [Google Scholar] [CrossRef]
  25. Kim, H.; Safonov, M.V. Boundary Harnack principle for second order elliptic equations with unbounded drift. Problems in mathematical analysis. No. 61. J. Math. Sci. 2011, 1, 127–143. [Google Scholar] [CrossRef]
  26. Koike, S.; Swiech, A. Weak Harnack inequality for fully nonlinear uniformly elliptic PDE with unbounded ingredients. J. Math. Soc. Jpn 2009, 61, 723–755. [Google Scholar] [CrossRef]
  27. Koike, S. Local maximum principle for Lp-viscosity solutions of fully nonlinear elliptic PDEs with unbounded coefficients. Commun. Pure Appl. Anal 2012, 11, 1897–1910. [Google Scholar] [CrossRef]
  28. Krylov, N.V.; Safonov, M.V. An estimate for the probability of a diffusion process hitting a set of positive measure (Russian). Dokl. Akad. Nauk SSSR 1979, 1, 18–30. [Google Scholar]
  29. Ladyzhenskaya, O.A.; Ural’tseva, N. Linear and quasilinear elliptic equations. Academic Press: New York-London, 1968. [Google Scholar]
  30. Lian, Y.; Zhang, K. Boundary Lipschitz regularity and the Hopf lemma for fully nonlinear elliptic equations. arXiv 2018. [Google Scholar] [CrossRef]
  31. Ros-Oton, X.; Torres-Latorre, D. New boundary Harnack inequalities with right hand side. J. Differ Equations 2021, 288, 204–249. [Google Scholar] [CrossRef]
  32. Sirakov, B. Boundary Harnack estimates and quantitative strong maximum principles for uniformly elliptic PDE. Int. Math. Res. Not. IMRN 2017, 2018, 7457–7482. [Google Scholar] [CrossRef]
  33. Sirakov, B. Global integrability and boundary estimates for uniformly elliptic PDE in divergence form. arXiv 2019. [Google Scholar] [CrossRef]
  34. Trudinger, N.S. Local estimates for subsolutions and supersolutions of general second order elliptic quasilinear equations. Invent. Math. 1980, 61, 67–79. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated