Home On viscosity and weak solutions for non-homogeneous p-Laplace equations
Article Open Access

On viscosity and weak solutions for non-homogeneous p-Laplace equations

  • Maria Medina and Pablo Ochoa EMAIL logo
Published/Copyright: June 4, 2017

Abstract

In this manuscript, we study the relation between viscosity and weak solutions for non-homogeneous p-Laplace equations with lower-order term depending on x, u and u . More precisely, we prove that any locally bounded viscosity solution constitutes a weak solution, extending results presented in Juutinen, Lindqvist and Manfredi [9], and Julin and Juutinen [6]. Moreover, we provide a converse statement in the full case under extra assumptions on the data.

1 Introduction and main results

In this work, we consider the following degenerate (or singular) elliptic equations of p-Laplacian type:

(1.1) - div ( | u | p - 2 u ) = f ( x , u , u ) ,

defined in an open and bounded set Ω n and for 1 < p < . The modulus of ellipticity of the p-Laplace operator is | u | p - 2 . When p > 2 , the modulus vanishes whenever u = 0 , and the equation is called degenerate at those points where that occurs. On the other hand, for p < 2 , the modulus becomes infinite when u = 0 , and the equation is called singular at those points. Observe that the case p = 2 is just the linear case and corresponds to the Laplace operator.

Different notions of solutions have been formulated for equation (1.1). We are interested in the relation between Sobolev weak solutions and viscosity solutions. For the homogeneous p-Laplace equation, this relation has already been studied by Juutinen, Lindqvist and Manfredi in [9], via the notion of p-harmonic, p-subharmonic and p-superharmonic functions. Roughly speaking, a p-harmonic function is a continuous function which solves, weakly, the homogeneous p-Laplace equation, and a p-superharmonic (p-subharmonic) function is a lower (upper) semicontinuous function that admits comparison with p-harmonic functions from below (above).

In [9], Juutinen, Lindqvist and Manfredi showed that the notion of p-harmonic solution is equivalent to the notion of viscosity solution. Moreover, it was shown in [13] that locally bounded p-harmonic functions are weak solutions. Conversely, every weak solution to the homogeneous p-Laplace equation has a representative which is lower semicontinuous and it is p-harmonic. We refer the interested reader to [5] for further details. In this way, there is an equivalence between the notion of weak and viscosity solutions for the homogeneous framework. It is worth to mention that a different and simpler proof of this equivalence was stated by Julin and Juutinen in [6] by using inf and sup convolutions. In turn, this reasoning was extended in [10] to more general second-order differential equations.

For the non-homogeneous case, the notion of p-harmonic functions is lost and we need to study directly the link between viscosity and Sobolev weak solutions. In [6], the authors showed that viscosity solutions of (1.1) are weak solutions in the case where f is continuous and depends only on x.

Our main goal in the present manuscript is to prove the equivalence of these two notions of solutions for the general structure (1.1). The implication that viscosity solutions are weak solutions is partially based on the work [6], but the non-homogeneous nature of the equation under consideration requires some extra effort to deal with the lower-order term.

On the other hand, the converse statement relies on comparison principles for weak solutions. To the best of our knowledge, the available comparison results for the full case f = f ( x , s , η ) require additional limitations in the degenerate case which do not appear in the singular context (compare Theorem A.1 and Theorem A.2). Moreover, we believe that the assumption that weak subsolutions and weak supersolutions belong to 𝒞 1 or to the Sobolev space W loc 1 , in order to have comparison is not a strong limitation since we are interested in the equivalence of weak and viscosity solutions, and for weak solutions the 𝒞 1 , α -regularity holds (see [4, 16]). Finally, in the quasi-linear case f = f ( x , u ) there is no need to impose higher regularity than W loc 1 , p 𝒞 on the solutions. We refer the reader to [15] for a survey of maximum principles and comparison results for general structures in divergence form.

Finally, we stress that the equivalence between weak and viscosity solutions may be used to prove relevant properties on the solutions. As an example, in [7], Juutinen and Lindqvist prove a Radó’s-type theorem for p-harmonic functions. Roughly speaking, they state that if a function u solves, weakly, the homogeneous p-Laplace equation in the complement of the set where u vanishes, then it is a solution in the whole set. It is an open problem to obtain a similar result for equations like (1.1). We shall return to this issue in a subsequent paper.

We recall that the p-Laplace operator is defined as

Δ p u := div ( | u | p - 2 u ) .

Let us state the different type of solutions to (1.1) we will manage.

Definition 1.1 (Sobolev weak solution).

A function u W loc 1 , p ( Ω ) is a weak supersolution to (1.1) if

Ω | u | p - 2 u ψ Ω f ( x , u , u ) ψ

for all non-negative ψ C 0 ( Ω ) . On the other hand, u is a weak subsolution if - u is a weak supersolution of the equation - Δ p u = - f ( x , - u , - u ) . We call u a weak solution if it is both a weak subsolution and a weak supersolution to (1.1).

Due to the non-homogeneous nature of (1.1), viscosity solutions are stated as in [6], considering semicontinuous envelopes of the p-Laplace operator. More precisely, we have the following definition.

Definition 1.2.

A lower semicontinuous function u : Ω ( - , + ] is a viscosity supersolution to (1.1) if u + and for every ϕ C 2 ( Ω ) such that ϕ ( x 0 ) = u ( x 0 ) , u ( x ) ϕ ( x ) and ϕ ( x ) 0 for all x x 0 , there holds

(1.2) lim r 0 sup x B r ( x 0 ) { x 0 } ( - Δ p ϕ ( x ) ) f ( x 0 , u ( x 0 ) , ϕ ( x 0 ) ) .

A function u is a viscosity subsolution if - u is a viscosity supersolution to the equation - Δ p u = - f ( x , - u , - u ) , and it is a viscosity solution if it is both a viscosity sub- and supersolution.

Remark 1.3.

Notice that condition (1.2) is established this way to avoid the problems derived from having ϕ ( x 0 ) = 0 in the case 1 < p < 2 . If p 2 , this condition can be simply replaced by

- Δ p ϕ ( x 0 ) f ( x 0 , u ( x 0 ) , ϕ ( x 0 ) ) .

We now list the main contributions of our work. The results are stated for supersolutions, but they hold for subsolutions as well.

Theorem 1.4.

Let 1 < p < . Assume that f = f ( x , s , η ) is uniformly continuous in Ω × R × R n , non-increasing in s, and satisfies the growth condition

(1.3) | f ( x , s , η ) | γ ( | s | ) | η | p - 1 + ϕ ( x ) ,

where γ 0 is continuous, and ϕ L loc ( Ω ) . Hence, if u L loc ( Ω ) is a viscosity supersolution to (1.1), then it is a weak supersolution to (1.1).

A converse of Theorem 1.4 is given below.

Theorem 1.5.

Assume that f = f ( x , s , η ) is continuous in Ω × R × R n , non-increasing in s, and locally Lipschitz continuous with respect to η. Hence we have the following:

  1. If 1 < p 2 and if u W loc 1 , ( Ω ) is a weak supersolution to ( 1.1 ), then it is a viscosity supersolution to ( 1.1 ) in Ω.

  2. If p > 2 , f ( x , s , 0 ) = 0 for x Ω and s , and if u 𝒞 1 ( Ω ) is a weak supersolution to ( 1.1 ), then it is a viscosity supersolution to ( 1.1 ) in Ω.

  3. Finally, if p > 2 and if u W loc 1 , ( Ω ) is a weak supersolution to ( 1.1 ), with u 0 in Ω , then it is a viscosity supersolution to ( 1.1 ).

Remark 1.6.

According to recent results (see [11]), it is possible to weak the locally Lipschitz assumption in Theorem 1.5 when f takes some particular forms or it satisfies extra convexity and coercivity assumptions. For instance, as a consequence of the results in [11], if

f ( x , η ) = c | η | q + ϕ ( x ) , ϕ > 0 , c , η n ,  and  q [ 1 , p ] ,

then Theorem 1.5 holds for bounded supersolutions u in W 1 , p ( Ω ) 𝒞 ( Ω ) . It is also a consequence of [11, Theorem 1.3] that the same conclusion is obtained when

f ( x , s , η ) = h ( s ) | η | p - 1 + ϕ ( x ) ,

where h 0 is decreasing and ϕ 0 .

In view of the available regularity theory for weak solutions of (1.1), we have the following equivalence.

Corollary 1.7.

Let 1 < p < . Assume that f = f ( x , s , η ) is uniformly continuous, locally Lipschitz in η, non-increasing in s, and satisfies the growth condition (1.3). Additionally, assume that f ( x , s , 0 ) = 0 for x Ω and s R when p > 2 . Then u is a weak solution to (1.1) if and only if it is a viscosity solution to (1.1).

We point out that, in the degenerate case, it is possible to remove the assumption f ( x , s , 0 ) = 0 by imposing the non-vanishing of the gradient of the weak solution in the whole Ω. This is a straightforward consequence of Theorem 1.5 (iii).

In the particular case where f does not depend on η, we have the following converse to Theorem 1.4 which does not require the locally Lipschitz regularity of the solutions.

Theorem 1.8.

Let 1 < p < . Suppose that f = f ( x , s ) is continuous in Ω × R and non-increasing in s. If u W loc 1 , p ( Ω ) C ( Ω ) is a weak supersolution to (1.1), then it is a viscosity supersolution to (1.1).

Let us briefly discuss the above hypotheses on f. Firstly, assuming that f is non-increasing and introducing the operator

F ( x , s , η , 𝒳 ) := - | η | p - 2 ( tr ( 𝒳 ) + p - 2 | η | 2 𝒳 η η ) - f ( x , s , η ) , p 2 ,

we derive that F is proper, that is, F is non-increasing in 𝒳 and non-decreasing in s, which is a standard and useful assumption in the theory of viscosity solutions [1]. For instance, it allows to get the equivalence between classical solutions ( 𝒞 2 functions which satisfy the equations pointwise) and 𝒞 2 viscosity solutions. On the other hand, the growth property (1.3) implies the 𝒞 1 , α -regularity of weak solutions to (1.1) (see [4, 17, 16]). Moreover, under a regular Dirichlet boundary condition φ 𝒞 1 , α , the 𝒞 1 , α -regularity up to the boundary of weak solutions follows. For further details, see the reference [12]. Finally, the extra assumption f ( x , s , 0 ) = 0 appearing in Theorem 1.5 in the degenerate case is used to remove critical sets of points of the weak solution (see reference [8]). Hence, it allows the application of comparison results without assuming the non-vanishing of the gradients. We point out that other properties of f = f ( x , s , η ) , as more regularity on s and η and convexity-like conditions, may be employed to ensure comparison for weak solutions. We refer the reader to [11] and the references therein for more details.

It is worth mentioning that many equations appearing in the literature have the structure of (1.1) with the lower-order term satisfying the above assumptions on f. We refer the reader to [14, 15, 3, 2] and the references therein for examples of such f.

The paper is organized as follows: In Section 2, we provide some preliminary results concerning properties of infimal convolutions (which will be the main tool in the proof of Theorem 1.4) and a convergence result. In addition, we prove a Caccioppoli-type estimate that will provide important uniform bounds, fundamental when using approximation arguments. This result is interesting in itself.

Section 3 contains the proof of the main result of the paper, Theorem 1.4, that states under which conditions on the non-homogeneous function f in (1.1) viscosity solutions are actually weak solutions. This proof is divided into two major cases: the singular and the degenerate scenario, thus, although both cases rely on the same idea, different approximations and estimates are needed depending on the range of p.

In Section 4, we prove the reverse statement, that is, weak solutions of (1.1) are viscosity solutions. This result is based on comparison arguments, and this will determine the conditions we will need to impose on f. Finally, in Appendix A we give, for the sake of completeness, precise references and state the comparison results that we use in Section 4.

2 Preliminary results

2.1 Infimal convolution

Let us define the infimal convolution of a function u as

(2.1) u ε ( x ) := inf y Ω ( u ( y ) + | x - y | q q ε q - 1 ) ,

where q 2 and ε > 0 .

We recall some useful properties of u ε . Let u : Ω be bounded and lower semicontinuous in Ω. It is well known that u ε is an increasing sequence of semiconcave functions in Ω, which converges pointwise to u. Hence, u ε is locally Lipschitz and twice differentiable a.e. in Ω. Moreover, it is possible to write

u ε ( x ) = inf y B r ( ε ) ( x ) Ω ( u ( y ) + | x - y | q q ε q - 1 )

for r ( ε ) 0 as ε 0 . For these and further properties, see [6, Lemma A.1.] and [1].

The next lemma is the counterpart of [6, Lemma A.1 (iii)] for our setting.

Lemma 2.1.

Suppose that u : Ω R is bounded and lower semicontinuous in Ω. Let f = f ( x , s , η ) be continuous in Ω × R × R n and non-increasing in s. If u is a viscosity supersolution to

(2.2) - Δ p u = f ( x , u , u )

in Ω for 1 < p < , then u ε is a viscosity supersolution to

- Δ p u ε = f ε ( x , u ε , u ε )

in Ω ε := { x Ω : dist ( x , Ω ) > r ( ε ) } , where

f ε ( x , s , η ) := inf y B r ( ε ) ( x ) f ( y , s , η ) .

Proof.

We start by noticing that

u ε ( x ) = inf z B r ( ε ) ( 0 ) ( u ( z + x ) + | z | q q ε q - 1 ) , x Ω ε .

Let us see first that for every z B r ( ε ) ( 0 ) , the function

ϕ z ( x ) := u ( z + x ) + | z | q q ε q - 1

is a viscosity supersolution to - Δ p ϕ z = f ε in Ω ε . Indeed, let x 0 Ω ε and φ 𝒞 2 ( Ω ε ) so that

min Ω ε ( ϕ z - φ ) = ( ϕ z - φ ) ( x 0 ) = 0 .

We assume that φ ( x ) 0 for all x x 0 if 1 < p < 2 . Making y := z + x , y 0 := z + x 0 and

φ ~ ( y ) := φ ( y - z ) - | z | q q ε q - 1 ,

we derive that u - φ ~ has a local minimum at y 0 , and indeed ( u - φ ~ ) ( y 0 ) = 0 . Since u is a viscosity supersolution to (2.2), there follows

lim ρ 0 sup x B ρ ( y 0 ) { y 0 } ( - Δ p φ ~ ( x ) ) f ( y 0 , φ ~ ( y 0 ) , φ ~ ( y 0 ) ) .

Therefore,

(2.3) lim ρ 0 sup x B ρ ( x 0 ) { x 0 } ( - Δ p φ ( x ) ) = lim ρ 0 sup x B ρ ( y 0 ) { y 0 } ( - Δ p φ ~ ( x ) ) f ( y 0 , φ ~ ( y 0 ) , φ ~ ( y 0 ) ) = f ( z + x 0 , φ ~ ( z + x 0 ) , φ ( x 0 ) ) = f ( z + x 0 , φ ( x 0 ) - | z | q q ε q - 1 , φ ( x 0 ) ) f ( z + x 0 , φ ( x 0 ) , φ ( x 0 ) ) f ε ( x 0 , φ ( x 0 ) , φ ( x 0 ) ) ,

where we have used that f is non-increasing in the second variable. Let us see now that, since u ε is an infimum of supersolutions, it is itself a supersolution (observe that u ε is continuous, since it is locally Lipschitz). Let x 0 Ω ε and ϕ 𝒞 2 ( Ω ε ) so that

(2.4) min Ω ε ( u ε - ϕ ) = ( u ε - ϕ ) ( x 0 ) = 0 .

Again, ϕ ( x ) 0 for all x x 0 in the singular scenario. Moreover, we may assume that the minimum is strict. For each n, there exists z n B r ( ε ) ( 0 ) such that

(2.5) u ( z n + x 0 ) + | z n | q q ε q - 1 < u ε ( x 0 ) + 1 n .

Let x n be a sequence of points in B ¯ r ( x 0 ) Ω ε so that

u ( z n + x n ) + | z n | q q ε q - 1 - ϕ ( x n ) u ( z n + x ) + | z n | q q ε q - 1 - ϕ ( x )

for all x B ¯ r ( x 0 ) , i.e., ( ϕ z n - ϕ ) has a minimum in B ¯ r ( x 0 ) at x n . Up to a subsequence, x n y 0 as n . Furthermore, by (2.5),

(2.6) u ε ( x n ) - ϕ ( x n ) u ( z n + x n ) + | z n | q q ε q - 1 - ϕ ( x n ) u ( z n + x 0 ) + | z n | q q ε q - 1 - ϕ ( x 0 ) u ε ( x 0 ) + 1 n - ϕ ( x 0 ) .

Taking liminf and using the lower semicontinuity of u ε , we derive

u ε ( y 0 ) - ϕ ( y 0 ) u ε ( x 0 ) - ϕ ( x 0 ) .

Since the minimum in (2.4) is strict, we must have y 0 = x 0 . Moreover, taking

φ ( x ) := ϕ ( x ) + ( ϕ z n - ϕ ) ( x n )

in (2.3), we have

lim ρ 0 sup x B ρ ( x n ) { x n } ( - Δ p ϕ ( x ) ) f ( z n + x n , u ( z n + x n ) + | z n | q q ε q - 1 , ϕ ( x n ) ) .

Since f is non-increasing with respect to the second variable, by (2.6) we obtain

lim ρ 0 sup x B ρ ( x n ) { x n } ( - Δ p ϕ ( x ) ) f ( z n + x n , u ε ( x 0 ) + 1 n - ϕ ( x 0 ) + ϕ ( x n ) , ϕ ( x n ) ) .

As n , there holds

lim ρ 0 sup x B ρ ( x 0 ) { x 0 } ( - Δ p ϕ ( x ) ) f ( z + x 0 , u ε ( x 0 ) , ϕ ( x 0 ) )

for some z B r ( 0 ) ¯ . Therefore,

lim ρ 0 sup x B ρ ( x 0 ) { x 0 } ( - Δ p ϕ ( x ) ) f ε ( x 0 , ϕ ( x 0 ) , ϕ ( x 0 ) ) ,

and we conclude that u ε is a viscosity supersolution of

- Δ p u ε = f ε ( x , u ε , u ε ) in  Ω ε ,

as desired. ∎

The next lemma states the weak convergence of the lower-order terms in the particular situation of infimal convolutions.

Lemma 2.2.

Let f = f ( x , s , η ) be a uniformly continuous function, which satisfies the growth condition (1.3). Assume that u W loc 1 , p ( Ω ) is locally bounded and lower semicontinuous in Ω. For each ε > 0 define u ε as in (2.1) and f ε as in Lemma 2.1. Then, if u ε converges to u in L loc p ( Ω ) , the following holds:

lim ε 0 Ω f ε ( x , u ε , u ε ) ψ 𝑑 x = Ω f ( x , u , u ) ψ 𝑑 x

for every non-negative ψ C 0 ( Ω ) .

Proof.

Let ψ 𝒞 0 ( Ω ) and denote K := spt ( ψ ) . Consider ε > 0 small enough so that

K K Ω ,

where K := x K B r ( ε ) ( x ) ¯ . Since f is uniformly continuous in K × × n , for every ρ > 0 there exists δ > 0 such that

| f ( x , u ε ( x ) , u ε ( x ) ) - f ( y , u ε ( x ) , u ε ( x ) ) | < ρ if  | x - y | < δ , x , y K .

Choose ε 0 > 0 so that r ( ε ) < δ for every ε < ε 0 . Thus, from the previous inequality we get

f ( x , u ε ( x ) , u ε ( x ) ) < ρ + f ( y , u ε ( x ) , u ε ( x ) )

for every x K and y B r ( ε ) ( x ) . In particular,

f ( x , u ε ( x ) , u ε ( x ) ) < ρ + f ε ( x , u ε ( x ) , u ε ( x ) ) ,

and therefore

0 | f ( x , u ε ( x ) , u ε ( x ) ) - f ε ( x , u ε ( x ) , u ε ( x ) ) | < ρ .

Hence we arrive at the estimate

(2.7) Ω | f ( x , u ε , u ε ) - f ε ( x , u ε , u ε ) | ψ 𝑑 x ρ ψ L ( K ) | K | .

On the other hand, due to the continuity of f and the convergences of u ε and u ε ,

f ( x , u ε ( x ) , u ε ( x ) ) f ( x , u ( x ) , u ( x ) ) a.e. in  Ω .

Observe that

u ε 0 u ε u  for all  ε ε 0 .

Since u ε 0 , u belong to L loc ( Ω ) , there exists a uniform constant C > 0 so that

u ε L ( K ) C , ε ε 0 .

Thus, in view of the growth estimate on f and the continuity of γ, we have, for an appropriate positive constant C,

(2.8) | f ( x , u ε ( x ) , u ε ( x ) ) | C | u ε ( x ) | p - 1 + ϕ ( x ) .

Since | u ε | p - 1 L loc p / ( p - 1 ) ( Ω ) , Hölder’s inequality and the strong convergence of u ε imply

K | u ε | p - 1 C u ε L p ( K ) p - 1 C for all  ε .

By (2.7), (2.8) and the Lebesgue dominated convergence theorem, we conclude

lim ε 0 K f ( x , u ε , u ε ) ψ 𝑑 x = K f ( x , u , u ) ψ 𝑑 x .

2.2 A Caccioppoli’s estimate

In the next lemma we provide a Caccioppoli’s estimate for the L loc p -norm of the gradients of weak solutions.

Lemma 2.3.

Let u W 1 , p ( Ω ) be a locally bounded weak supersolution to (1.1). Assume that f is continuous in Ω × R × R n and satisfies the growth bound (1.3). Then there exists a constant C = C ( p , Ω , ϕ , γ ) > 0 such that for all test function ξ C 0 ( Ω ) , 0 ξ 1 , we have

Ω | u | p ξ p 𝑑 x C [ ( osc K u ) p Ω ( | ξ | p + 1 ) 𝑑 x + osc K u ] ,

where osc K u := sup K u - inf K u , and K := spt ( ξ ) .

Proof.

Let ξ 𝒞 0 ( Ω ) and K Ω be as in the lemma. Consider the test function

ψ ( x ) := ( sup K u - u ( x ) ) ξ p ( x ) , x Ω .

Then

Ω f ( x , u , u ) ψ 𝑑 x Ω | u | p - 2 u ψ d x = - Ω | u | p - 2 u [ ξ p u - p ξ p - 1 ξ ( sup K u - u ) ] 𝑑 x .

Therefore,

(2.9) Ω | u | p ξ p 𝑑 x p Ω ξ p - 1 | u | p - 2 u ξ ( sup K u - u ) 𝑑 x - Ω f ( x , u , u ) ψ 𝑑 x .

Observe that

Ω | | u | p - 2 u | 𝑑 x Ω | u | p - 1 𝑑 x ,

which shows that | u | p - 2 u L p / ( p - 1 ) ( Ω ) . Hence, Young’s inequality

a b δ a q + δ - 1 / ( q - 1 ) b q ,

where q and q are conjugate exponents, implies that the first integral on the right-hand side of (2.9) may be bounded by

δ Ω | u | p ξ p + δ 1 - p Ω p p | ξ | p ( osc K u ) p 𝑑 x .

Moreover, by (1.3) we have

(2.10) f ( x , u ( x ) , u ( x ) ) - γ | u ( x ) | p - 1 - ϕ L ( K )

for all x in the support of ξ, where γ := sup x K | γ ( u ( x ) ) | . Therefore, the second integral in (2.9) is estimated from above by

γ Ω | u | p - 1 ( sup K u - u ) ξ p 𝑑 x + C ( Ω , ϕ ) osc K u ,

where C ( Ω , ϕ ) is a positive constant. The assumption ξ 1 and Young’s inequality yield

γ Ω | u | p - 1 ( sup K u - u ) ξ p 𝑑 x γ Ω | u | p - 1 ( sup K u - u ) ξ p - 1 𝑑 x δ Ω | u | p ξ p 𝑑 x + δ 1 - p C ( p , Ω , γ ) ( osc K u ) p .

Therefore,

Ω | u | p ξ p 𝑑 x 2 δ Ω | u | p ξ p 𝑑 x + δ 1 - p C ( p , Ω , γ , ϕ ) [ ( osc K u ) p Ω ( | ξ | p + 1 ) 𝑑 x + osc K u ] .

Taking δ < 1 / 2 , we derive Caccioppoli’s estimate. ∎

3 Proof of Theorem 1.4

3.1 Degenerate case

We begin with the range p 2 .

Proof of Theorem 1.4.

Let u ε be the infimal convolution defined in (2.1) with q = 2 . Then

ϕ ( x ) := u ε ( x ) - C | x | 2

is concave in Ω r ( ε ) (see [6, Lemma A.2.]). By Aleksandrov’s Theorem, ϕ is twice differentiable almost everywhere in Ω r ( ε ) , and so is u ε . Therefore by Lemma 2.1,

- Δ p u ε ( x ) f ε ( x , u ε ( x ) , u ε ( x ) )

a.e. in Ω r ( ε ) . Furthermore,

Ω | u ε | p - 2 u ε ψ Ω ( - Δ p u ε ) ψ

for all non-negative test functions ψ (see the proof of [6, Theorem 3.1]). Hence, we derive

Ω f ε ( x , u ε , u ε ) ψ 𝑑 x Ω | u ε | p - 2 u ε ψ d x

for all non-negative test functions ψ and all ε > 0 . We claim that, as ε 0 , there holds

Ω f ( x , u , u ) ψ 𝑑 x Ω | u | p - 2 u ψ d x .

To prove the claim, observe first that Caccioppoli’s estimate allows us to conclude that

| u ε | p - 2 u ε

converges weakly in L loc p / ( p - 1 ) ( Ω ) . Indeed, for any compact set K Ω , choose an open set U Ω containing K and a non-negative test function 0 ξ 1 so that

K K := spt ξ U

and ξ = 1 in K. Then

(3.1) K | | u ε | p - 2 u ε | p / ( p - 1 ) 𝑑 x K | u ε | p 𝑑 x Ω | u ε | p ξ p 𝑑 x .

Observe that since f satisfies (1.3), the lower term f ε verifies the bound (2.10). Therefore, Lemma 2.3 applies and the right-hand side of (3.1) is bounded from above by

(3.2) C [ ( osc K u ε ) p Ω ( | ξ | p + 1 ) 𝑑 x + osc K u ε ] .

Moreover, since u ε is an increasing sequence and converges pointwise to u in Ω, we have

osc K u ε sup K u - inf K u ε 0

for all ε < ε 0 . Then in view of (3.1), (3.2) and the above comments, we can find a uniform bound for the integrals

K | | u ε | p - 2 u ε | p / ( p - 1 ) 𝑑 x , Ω | u ε | p ξ p 𝑑 x .

Hence | u ε | p - 2 u ε converges weakly in L loc p / ( p - 1 ) ( Ω ) , and u ε converges weakly in L loc p ( Ω ) . Since u ε converges pointwise to u, we derive that u W loc 1 , p ( Ω ) and u ε converges weakly in W loc 1 , p ( Ω ) to u.

More can be said: u ε converges strongly in L loc p ( Ω ) to u . Indeed, take

ϕ ( x ) := ( u ( x ) - u ε ( x ) ) θ ( x ) , x Ω ,

where θ is a non-negative smooth test function compactly supported in Ω. From

Ω | u ε | p - 2 u ε ϕ d x Ω f ε ( x , u ε , u ε ) ϕ 𝑑 x

we get

(3.3) Ω [ | u | p - 2 u - | u ε | p - 2 u ε ] ( u - u ε ) θ d x - Ω f ε ( x , u ε , u ε ) ϕ 𝑑 x + Ω | u | p - 2 u ( u - u ε ) θ d x .

By the weak convergence of u ε to u in W loc 1 , p ( Ω ) , the last integral in (3.3) tends to 0 as ε 0 . The left-hand side is given by

(3.4) Ω θ [ | u | p - 2 u - | u ε | p - 2 u ε ] ( u - u ε ) d x + Ω ( u - u ε ) [ | u | p - 2 u - | u ε | p - 2 u ε ] θ d x .

The second integral in (3.4) is estimated in absolute value by

θ L ( Ω ) ( spt θ | u - u ε | p 𝑑 x ) 1 / p [ ( spt θ | u | p 𝑑 x ) ( p - 1 ) / p + ( spt θ | u ε | p 𝑑 x ) ( p - 1 ) / p ] ,

which tends to 0 as ε 0 . Moreover, since

- spt θ f ε ( x , u ε , u ε ) ϕ 𝑑 x γ Ω | u ε | p - 1 ( u - u ε ) θ 𝑑 x + ϕ L ( spt ( θ ) ) Ω ( u - u ε ) θ 𝑑 x ,

with γ := sup x spt ( θ ) | γ ( u ε ( x ) ) | , does not depend on ε, it also holds that

lim sup ε 0 [ - spt θ f ε ( x , u ε , u ε ) ϕ 𝑑 x ] = 0 .

Hence

(3.5) lim ε 0 Ω θ [ | u | p - 2 u - | u ε | p - 2 u ε ] ( u - u ε ) d x = 0 ,

where we have used the fact that the integrand is always non-negative. Therefore, using the inequality

2 p - 2 | u ( x ) - u ε ( x ) | p [ | u ( x ) | p - 2 u ( x ) - | u ε ( x ) | p - 2 u ε ( x ) ] ( u ( x ) - u ε ( x ) )

valid for all p 2 , and (3.5), we conclude the strong convergence of u ε in L loc p ( Ω ) . Finally, (3.5) together with [5, Lemma 3.73]), implies

| u ε | p - 2 u ε | u | p - 2 u in  L loc p / ( p - 1 ) ( Ω ) ,

and, in turn, the strong convergence of the gradients u ε and Lemma 2.2 gives

lim ε 0 Ω f ε ( x , u ε , u ε ) ψ 𝑑 x = Ω f ( x , u , u ) ψ 𝑑 x .

This ends the proof of the claim and we deduce that u is a weak supersolution. ∎

3.2 The singular case: 1 < p < 2

Consider now the infimal convolution given in (2.1) choosing q > p / ( p - 1 ) , i.e.,

(3.6) u ε ( x ) := inf y Ω ( u ( y ) + | x - y | q q ε q - 1 ) .

Notice that q > 2 for 1 < p < 2 .

We need the following auxiliary result, which is an adaptation of [6, Lemma 4.3].

Lemma 3.1.

Suppose that u is a bounded viscosity supersolution to (1.1). If there is x ^ Ω r ( ε ) such that u ε is differentiable at x ^ and u ε ( x ^ ) = 0 , then f ε ( x ^ , u ε ( x ^ ) , u ε ( x ^ ) ) 0 .

Proof.

From [6, Lemma 4.3] we know that u ε ( x ^ ) = u ( x ^ ) , and hence

u ( y ) + | x ^ - y | q q ε q - 1 u ( x ^ ) for every  y Ω .

Define

ψ ( y ) := u ( x ^ ) - | x ^ - y | q q ε q - 1 , y Ω ,

which satisfies ψ 𝒞 2 ( Ω ) , ψ ( x ^ ) = 0 and

(3.7) lim r 0 sup y B r ( x ^ ) { x ^ } ( - Δ p ψ ( y ) ) = 0

in view of q > p / ( p - 1 ) . Since ψ ( x ^ ) = u ( x ^ ) , ψ ( y ) u ( y ) for y Ω , ψ ( x ) 0 for all x x ^ , and u is a viscosity supersolution to (1.1),

lim r 0 sup y B r ( x ^ ) { x ^ } ( - Δ p ψ ( y ) ) f ( x ^ , ψ ( x ^ ) , ψ ( x ^ ) ) .

Noticing that ψ ( x ^ ) = u ε ( x ^ ) and ψ ( x ^ ) = u ε ( x ^ ) = 0 , by (3.7) we conclude

0 f ( x ^ , ψ ( x ^ ) , ψ ( x ^ ) ) = f ( x ^ , u ε ( x ^ ) , u ε ( x ^ ) ) f ε ( x ^ , u ε ( x ^ ) , u ε ( x ^ ) ) ,

as desired. ∎

We can prove now Theorem 1.4 in the case 1 < p < 2 .

Proof of Theorem 1.4.

Let u ε be defined in (3.6). Proceeding as in the degenerate case, by Aleksandrov’s theorem and Lemma 2.1 we obtain

- Δ p u ε f ε ( x , u ε , u ε )

a.e. in Ω r ( ε ) { u ε = 0 } . Performing the same approximation argument as in the proof of [6, Theorem 4.1], we reach that

Ω | u ε | p - 2 u ε ψ d x Ω { u ε = 0 } f ε ( x , u ε , u ε ) ψ 𝑑 x

for every ψ 𝒞 0 ( Ω ) , ψ 0 . Therefore and since by Lemma 3.1 we know f ε 0 in the set { x Ω : u ε ( x ) = 0 } , we get

(3.8) Ω | u ε | p - 2 u ε ψ d x Ω f ε ( x , u ε , u ε ) ψ 𝑑 x .

Repeating the proof for the case p 2 (and noticing that Lemma 2.3 works for every 1 < p < ), we obtain the uniform boundedness of u ε in L loc p ( Ω ) and

(3.9) lim ε 0 K [ | u | p - 2 u - | u ε | p - 2 u ε ] ( u - u ε ) d x = 0

for any compact set K Ω , and from here the convergence

(3.10) | u ε | p - 2 u ε | u | p - 2 u in  L loc p / ( p - 1 ) ( Ω ) .

Using Hölder’s inequality and the vector inequality (see [4, Chapter I])

| a - b | 2 ( | a | + | b | ) 2 - p C ( | a | p - 2 a - | b | p - 2 b ) ( a - b ) , 1 < p < 2 ,

with C = C ( n , p ) and a , b n , we obtain

K | u - u ε | p 𝑑 x ( K | u - u ε | 2 ( | u | + | u ε | ) 2 - p 𝑑 x ) p / 2 ( K ( | u | + | u ε | ) p 𝑑 x ) ( 2 - p ) / 2 C ( K | u - u ε | 2 ( | u | + | u ε | ) 2 - p 𝑑 x ) p / 2 C ( K [ | u | p - 2 u - | u ε | p - 2 u ε ] ( u - u ε ) d x ) p / 2 .

Thus, from (3.9) we deduce that u ε converges to u in L loc p ( Ω ) . By (3.10) and Lemma 2.2 we can pass to the limit in (3.8) to conclude

Ω | u | p - 1 u ψ d x Ω f ( x , u , u ) ψ 𝑑 x .

4 Proofs of Theorem 1.5 and Theorem 1.8

Proof of Theorem 1.5 (i).

Let u W loc 1 , ( Ω ) be a weak supersolution to (1.1). To reach a contradiction, assume that u is not a viscosity supersolution. By assumption, there exist x 0 Ω and φ 𝒞 2 ( Ω ) so that φ ( x ) 0 for all x x 0 ,

(4.1) u ( x 0 ) = φ ( x 0 ) , u ( x ) > φ ( x ) for all  x x 0 ,

and

lim r 0 sup x B r ( x 0 ) { x 0 } ( - Δ p φ ( x ) ) < f ( x 0 , u ( x 0 ) , φ ( x 0 ) ) .

Moreover, by the W loc 1 , -regularity of u, we may assume that u is continuous in Ω. Thus, the map

x f ( x , u ( x ) , φ ( x ) )

is continuous in Ω, and (4.1) yields

lim r 0 sup x B r ( x 0 ) { x 0 } [ - Δ p φ ( x ) - f ( x , u ( x ) , φ ( x ) ) ] < 0 .

Hence, there exists some r 0 > 0 so that

(4.2) - Δ p φ f ( x , u ( x ) , φ ( x ) ) , x B r 0 ( x 0 ) { x 0 } .

Let

m := inf B r 0 ( x 0 ) ( u - φ ) .

Then by (4.1) we have m > 0 . Consider

φ ~ ( x ) := φ ( x ) + m , x Ω .

By (4.2), φ ~ is a weak subsolution to

(4.3) - Δ p v = f ~ ( x , v )

in B r 0 ( x 0 ) , where f ~ ( x , η ) := f ( x , u ( x ) , η ) . Observe that f ~ is continuous in Ω × n and locally Lipschitz in η. Moreover, in the weak sense, we have

- Δ p u f ( x , u , u ) = f ~ ( x , u ) ,

which shows that u is a weak supersolution to (4.3). In addition, u φ ~ on B r 0 ( x 0 ) . By the comparison Theorem A.2, we conclude that u φ ~ in B r 0 ( x 0 ) . This contradicts (4.1). ∎

Proof of Theorem 1.5 (ii).

For a given weak supersolution u 𝒞 1 ( Ω ) , by following the lines above and appealing to the comparison Theorem A.1, we can show that u is a viscosity supersolution in the non-critical set

{ x Ω : u ( x ) 0 } .

By [8, Corollary 4.4], which holds true for sub- and supersolutions, u is a viscosity supersolution in the whole set Ω. ∎

Proof of Theorem 1.5 (iii).

The proof follows similarly just by using the assumption u 0 in Ω together with the comparison Theorem A.1. ∎

Remark 4.1.

Observe that the W loc 1 , -regularity of u is only needed to apply the comparison principles. In the rest of the proof, the continuity of u would be enough.

Proof of Theorem 1.8.

The result follows by reproducing the proof of Theorem 1.5 using the comparison Theorem A.3 instead of Theorems A.2 and A.1. ∎

Funding statement: The first author was supported by the grant FONDECYT Postdoctorado 2016, No. 3160077. The second author was partially supported by CONICET and grant PICT 2015-1701 AGENCIA.

A Appendix

A.1 Comparison principles for weak solutions

In this section, we provide the comparison principles for weak solutions of (1.1) that we use in the proof of Theorem 1.5. As we pointed out in Section 1, other comparison results may be employed (see [11]).

The first one is contained in [15, Corollary 3.6.3].

Theorem A.1.

Assume that f = f ( x , s , η ) is continuous in Ω × R × R n , non-increasing in s, and locally Lipschitz continuous with respect to η in Ω × R × R n . Let u W loc 1 , ( Ω ) be a weak supersolution and let v W loc 1 , ( Ω ) be a weak subsolution to (1.1) in Ω for 1 < p < . Assume that | u | + | v | > 0 in Ω. If u v on Ω , then u v in Ω.

In the singular case, the assumptions on the gradients may be removed. See [15, Corollary 3.5.2].

Theorem A.2.

Let 1 < p 2 . Assume that f = f ( x , s , η ) is continuous in Ω × R × R n , non-increasing in s, and that it is locally Lipschitz continuous in η on compact subsets of its variables. Then if u W loc 1 , ( Ω ) is a weak supersolution and if v W loc 1 , ( Ω ) is a weak subsolution to (1.1) in Ω so that u v on Ω , then u v in Ω.

Finally, in the case where f does not depend on η, we have the following result (see [15, Corollary 3.4.2]).

Theorem A.3.

Suppose that f = f ( x , s ) is continuous in Ω × R and non-increasing in s. Let u W loc 1 , p ( Ω ) C ( Ω ) be a supersolution and v W loc 1 , p ( Ω ) C ( Ω ) a subsolution so that u v on Ω . Then u v in Ω.

Remark A.4.

In [15], Theorem A.1 is stated in a more general framework of equations in divergence form as

div ( A ( x , u , u ) ) = f ( x , u , u ) ,

where the operator A = A ( x , s , η ) is assumed to be continuous in Ω × × n , continuously differentiable with respect to s and η for all s and all η 0 , and elliptic in the sense that A ( x , s , η ) is positive definite in Ω × × ( n { 0 } ) . In the particular case of the p-Laplace operators

A p ( η ) = | η | p - 2 η , p 2 ,
A p ( η ) = { | η | p - 2 η if  η 0 , 0 , for  η = 0 , 1 < p 2 ,

all of the assumptions above are satisfied. The positive definiteness of A p is a consequence of

i , j = 1 n A i η j ( η ) ξ i ξ j c | η | p - 2 | ξ | 2

for a positive constant c. Finally, observe that A is uniformly elliptic for 0 < | η | C if p 2 . This allows the improved comparison result in Theorem A.2.

Acknowledgements

The authors would like to thank the anonymous referee for her/his comments.

References

[1] M. Crandall, H. Ishii and P.-L. Lions, User’s guide to viscosity solutions of second order partial differential equations, Bull. Amer. Math. Soc. (N.S.) 27 (1992), no. 1, 1–67. 10.1090/S0273-0979-1992-00266-5Search in Google Scholar

[2] M. Cuesta and P. Takáč, A strong comparison principle for positive solutions of degenerate elliptic equations, Differential Integral Equations 13 (2000), no. 4–6, 721–746. 10.57262/die/1356061247Search in Google Scholar

[3] L. Damascelli and B. Sciunzi, Regularity, monotonicity and symmetry of positive solutions of m-Laplace equations, J. Differential Equations 206 (2004), 483–515. 10.1016/j.jde.2004.05.012Search in Google Scholar

[4] E. DiBenedetto, 𝒞 1 + α -local regularity of weak solutions of degenerate elliptic equations, Nonlinear Anal. 7 (1983), 827–850. 10.1016/0362-546X(83)90061-5Search in Google Scholar

[5] J. Heinonen, T. Kilpelainen and O. Martio, Non-Linear Potential Theory of Degenerate Elliptic Equations, Dover Publications, Mineola, 2006. Search in Google Scholar

[6] V. Julin and P. Juutinen, A new proof for the equivalence of weak and viscosity solutions for the p-Laplace equation, Comm. Partial Differential Equations 37 (2012), no. 5, 934–946. 10.1080/03605302.2011.615878Search in Google Scholar

[7] P. Juutinen and P. Lindqvist, A theorem of Radó’s type for the solutions of a quasi-linear equation, Math. Res. Lett. 11 (2004), 31–34. 10.4310/MRL.2004.v11.n1.a4Search in Google Scholar

[8] P. Juutinen and P. Lindqvist, Removability of a level set for solutions of quasilinear equations, Comm. Partial Differential Equations 30 (2005), no. 3, 305–321. 10.1081/PDE-200050030Search in Google Scholar

[9] P. Juutinen, P. Lindqvist and J. Manfredi, On the equivalence of viscosity solutions and weak solutions for a quasilinear equation, SIAM J. Math. Anal. 33 (2001), no. 3, 699–717. 10.1137/S0036141000372179Search in Google Scholar

[10] N. Katzourakis, Nonsmooth convex functionals and feeble viscosity solutions of singular Euler–Lagrange equations, Calc. Var. Partial Differential Equations 54 (2015), no. 1, 275–298. 10.1007/s00526-014-0786-xSearch in Google Scholar

[11] T. Leonori, A. Porretta and G. Riey, Comparison principles for p-Laplace equations with lower order terms, Ann. Mat. Pura Appl. (4) (2016), 10.1007/s10231-016-0600-9. 10.1007/s10231-016-0600-9Search in Google Scholar

[12] G. Lieberman, Boundary regularity for solutions of degenerate elliptic equations, Nonlinear Anal. 12 (1988), no. 11, 1203–1219. 10.1016/0362-546X(88)90053-3Search in Google Scholar

[13] P. Lindqvist, On the definition and properties of p-superharmonic functions, J. Reine Angew. Math. 365 (1986), 67–70. Search in Google Scholar

[14] S. Merchán, L. Montoro and B. Sciunzi, On the Harnack inequality for quasilinear elliptic equations with a first order term, preprint (2016), https://arxiv.org/abs/1601.03863. 10.1017/S030821051700035XSearch in Google Scholar

[15] P. Pucci and J. Serrin, The Maximum Principle, Progr. Nonlinear Differential Equations Appl. 73, Birkhäuser, Boston, 2007. 10.1007/978-3-7643-8145-5Search in Google Scholar

[16] P. Pucci and R. Servadei, Regularity of weak solutions of homogeneous or inhomogeneous elliptic equations, Indiana Univ. Math. J. 57 (2008), no. 7, 3329–3363. 10.1512/iumj.2008.57.3525Search in Google Scholar

[17] J. Rakotoson and W. Ziemer, Local behaviour of solutions of quasilinear elliptic equations with general structure, Trans. Amer. Math. Soc. 319 (1990), no. 2, 747–764. 10.1090/S0002-9947-1990-0998128-9Search in Google Scholar

Received: 2017-01-13
Revised: 2017-02-28
Accepted: 2017-03-02
Published Online: 2017-06-04

© 2019 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Frontmatter
  2. Asymptotic behavior of evolution systems in arbitrary Banach spaces using general almost periodic splittings
  3. Solvability of a product-type system of difference equations with six parameters
  4. On Dirichlet problem for fractional p-Laplacian with singular non-linearity
  5. Absence of Lavrentiev gap for non-autonomous functionals with (p,q)-growth
  6. On a class of fully nonlinear parabolic equations
  7. On sign-changing solutions for (p,q)-Laplace equations with two parameters
  8. Weighted Caffarelli–Kohn–Nirenberg type inequalities related to Grushin type operators
  9. On the fractional p-Laplacian equations with weight and general datum
  10. An elliptic equation with an indefinite sublinear boundary condition
  11. Liouville-type theorems for elliptic equations in half-space with mixed boundary value conditions
  12. Well/ill-posedness for the dissipative Navier–Stokes system in generalized Carleson measure spaces
  13. Hypercontractivity, supercontractivity, ultraboundedness and stability in semilinear problems
  14. Theoretical analysis of a water wave model with a nonlocal viscous dispersive term using the diffusive approach
  15. A multiplicity result for asymptotically linear Kirchhoff equations
  16. Higher-order anisotropic models in phase separation
  17. Well-posedness and maximum principles for lattice reaction-diffusion equations
  18. Existence of a bound state solution for quasilinear Schrödinger equations
  19. Existence and concentration behavior of solutions for a class of quasilinear elliptic equations with critical growth
  20. Homoclinics for strongly indefinite almost periodic second order Hamiltonian systems
  21. A new method for converting boundary value problems for impulsive fractional differential equations to integral equations and its applications
  22. Diffusive logistic equations with harvesting and heterogeneity under strong growth rate
  23. On viscosity and weak solutions for non-homogeneous p-Laplace equations
  24. Periodic impulsive fractional differential equations
  25. A result of uniqueness of solutions of the Shigesada–Kawasaki–Teramoto equations
  26. Solutions of vectorial Hamilton–Jacobi equations are rank-one absolute minimisers in L
  27. Large solutions to non-divergence structure semilinear elliptic equations with inhomogeneous term
  28. The elliptic sinh-Gordon equation in a semi-strip
  29. The Gelfand problem for the 1-homogeneous p-Laplacian
  30. Boundary layers to a singularly perturbed Klein–Gordon–Maxwell–Proca system on a compact Riemannian manifold with boundary
  31. Subharmonic solutions of Hamiltonian systems displaying some kind of sublinear growth
  32. Multiple solutions for an elliptic system with indefinite Robin boundary conditions
  33. New solutions for critical Neumann problems in ℝ2
  34. A fractional Kirchhoff problem involving a singular term and a critical nonlinearity
  35. Existence and non-existence of solutions to a Hamiltonian strongly degenerate elliptic system
  36. Characterizing the strange term in critical size homogenization: Quasilinear equations with a general microscopic boundary condition
  37. Nonlocal perturbations of the fractional Choquard equation
  38. A pathological example in nonlinear spectral theory
  39. Infinitely many solutions for cubic nonlinear Schrödinger equations in dimension four
  40. On Cauchy–Liouville-type theorems
  41. Maximal Lp -Lq regularity to the Stokes problem with Navier boundary conditions
  42. Besov regularity for solutions of p-harmonic equations
  43. The classical theory of calculus of variations for generalized functions
  44. On the Cauchy problem of a degenerate parabolic-hyperbolic PDE with Lévy noise
  45. Hölder gradient estimates for a class of singular or degenerate parabolic equations
  46. Critical and subcritical fractional Trudinger–Moser-type inequalities on
  47. Multiple nonradial solutions for a nonlinear elliptic problem with singular and decaying radial potential
  48. Quantization of energy and weakly turbulent profiles of solutions to some damped second-order evolution equations
  49. An elliptic system with logarithmic nonlinearity
  50. The Caccioppoli ultrafunctions
  51. Equilibrium of a production economy with non-compact attainable allocations set
  52. Exact behavior around isolated singularity for semilinear elliptic equations with a log-type nonlinearity
  53. The higher integrability of weak solutions of porous medium systems
  54. Classification of stable solutions for boundary value problems with nonlinear boundary conditions on Riemannian manifolds with nonnegative Ricci curvature
  55. Regularity results for p-Laplacians in pre-fractal domains
  56. Carleman estimates and null controllability of a class of singular parabolic equations
  57. Limit profiles and uniqueness of ground states to the nonlinear Choquard equations
  58. On a measure of noncompactness in the space of regulated functions and its applications
  59. p-fractional Hardy–Schrödinger–Kirchhoff systems with critical nonlinearities
  60. On the well-posedness of a multiscale mathematical model for Lithium-ion batteries
  61. Global existence of a radiative Euler system coupled to an electromagnetic field
  62. On the existence of a weak solution for some singular p ( x ) -biharmonic equation with Navier boundary conditions
  63. Choquard-type equations with Hardy–Littlewood–Sobolev upper-critical growth
  64. Clustered solutions for supercritical elliptic equations on Riemannian manifolds
  65. Ground state solutions for the Hénon prescribed mean curvature equation
  66. Quasilinear equations with indefinite nonlinearity
  67. Concentrating solutions for a planar elliptic problem with large nonlinear exponent and Robin boundary condition
  68. Retraction of: Concentrating solutions for a planar elliptic problem with large nonlinear exponent and Robin boundary condition
Downloaded on 26.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/anona-2017-0005/html
Scroll to top button