Home Coupled Versus Uncoupled Blow-Up Rates in Cooperative n-Species Logistic Systems
Article Open Access

Coupled Versus Uncoupled Blow-Up Rates in Cooperative n-Species Logistic Systems

  • Julián López-Gómez EMAIL logo and Luis Maire
Published/Copyright: January 26, 2017

Abstract

This paper ascertains the exact boundary blow-up rates of the large positive solutions of a class of cooperative logistic systems involving n species in a general domain of N of class 𝒞2+ν, 0<ν<1. The problem models a population divided in groups whose individuals compete with those of the same group, while simultaneously they cooperate with the members of the remaining groups. Our main result provides with the exact blow-up rates along the edges of Ω from the values of the blow-up rates of the underlying uncoupled system. Rather astonishingly, these blow-up rates are independent of the strength of the cooperative effects, which play a secondary role in the analysis carried out in this paper. No previous result of this nature is available in the specialized literature for more than n=2 species.

1 Introduction

In this paper we ascertain the blow-up rate of the classical positive solution of the singular boundary value problem

(1.1) { - Δ u i = λ i ( x ) u i + j = 1 , j i n a i j u j - 𝔞 i ( x ) u i p i in  Ω , u i = + on  Ω ,  1 i n ,

where aij>0 and pi>1 are constants, Ω is a bounded subdomain of N, N, of class 𝒞2+ν for some ν(0,1), λi𝒞ν(Ω¯), and 𝔞i𝒞ν(Ω¯) satisfy 𝔞i(x)>0 for all xΩ and 1i,jn. By a solution, u, of (1.1) it is understood any function u𝒞2+ν(Ω) satisfying the system and such that

lim dist ( x , Ω ) 0 u i ( x ) = + , 1 i n .

In this paper, we suppose that 𝔞i are power-like near 𝔞i-1(0)Ω, in the sense that there exist 𝔟i,γi𝒞(Ω), with 𝔟i(z)>0 for all zΩ and γi0 on Ω, such that

(1.2) lim x z x Ω , z Ω 𝔞 i ( x ) 𝔟 i ( z ) [ dist ( x , Ω ) ] γ i ( z ) = 1 , 1 i n .

Problem (1.1) is the most natural way of coupling the fully uncoupled problem

(1.3) { - Δ u i = λ i ( x ) u i - 𝔞 i ( x ) u i p i in  Ω , u i = + on  Ω ,  1 i n ,

in a cooperative way, i.e., with aij>0 for all 1i,jn, ij, so that the strong maximum principle holds, [30]. A simpler prototype with n=2 of the system in (1.1) goes back to [34, 35, 36]. When aij<0 for some ij, the cooperative structure of (1.1) is lost and all the available comparison techniques in the context of cooperative systems might fail. Consequently, the results of this paper might not be true.

Since (1.3) consists of n-uncoupled singular boundary value problems of logistic type, and the logistic equation is the most standard one in population dynamics and mathematical biology, [27, 37, 38], the problem of analyzing the singular problem (1.1) should deserve a significative attention in spatial ecology. Indeed, the solutions of (1.1) provide us with the asymptotic profiles of the positive solutions of wide classes of cooperative parabolic systems in the presence of spatial heterogeneities, [2, 3, 27]. A more realistic model would be to have different diffusion rates for each species, measured by di>0, 1in, but dividing the i-th equation by di, one is naturally driven to deal with (1.1).

Although there is a huge amount of literature devoted to the existence and uniqueness of large positive solutions for the single generalized logistic equation, as it becomes apparent by simply looking at [25, 23, 24, 5, 4, 43, 19, 20, 32, 33, 9, 14, 6, 7, 8, 10, 12, 13, 40, 41, 42], and the rather complete list of references in [27], and even there are some fairly astonishing multiplicity results for large positive solutions, [31], the literature on systems is very short. Moreover, the few existing references deal with either two species, n=2, as [16, 2, 3] and the very recent works [15] and [29], or with the general case of n2 species but in the radially symmetric case, where the maximum principle gives the uniqueness of the positive solution of (1.1) based on a very tricky use of the strong maximum principle, as it has been recently established by the authors, [28]. This is the first paper available in the literature where the exact blow-up rates of each of the components of a n-species system have been ascertained.

Throughout this paper, for every zΩ, we set

(1.4) μ i ( z ) := γ i ( z ) + 2 p i - 1 , 1 i n .

According to [9, 14, 22], μi(z), 1in, provide us with the blow-up rates on Ω of the positive solutions of the uncoupled problem (1.3), in the sense that for every solution of (1.3), =(1,,n), one has that

lim x z x Ω , z Ω i ( x ) [ dist ( x , Ω ) ] - μ i ( z ) = ( μ i ( z ) ( μ i ( z ) + 1 ) 𝔟 i ( z ) ) 1 p i - 1 , 1 i n .

The main result of this paper provides us with the blow-up rates of all positive solutions of (1.1) in terms of μi(z), 1in, defined in (1.4). Precisely, for any given zΩ, suppose the n equations of (1.1) have been re-ordered so that

(1.5) 0 < μ n ( z ) μ n - 1 ( z ) μ 1 ( z ) .

Then, the next result holds.

Theorem 1.1.

Let zΩ be such that (1.5) is satisfied, and consider the next partition of the subscripts set

(1.6) { I + := { i { 1 , , n } : μ i ( z ) + 2 - μ 1 ( z ) > 0 } , I 0 := { i { 1 , , n } : μ i ( z ) + 2 - μ 1 ( z ) = 0 } , I - := { i { 1 , , n } : μ i ( z ) + 2 - μ 1 ( z ) < 0 } .

Let k{1,,n} be such that

I M := { i { 1 , , n } : μ i ( z ) = μ 1 ( z ) } = { 1 , , k } .

Then, any positive solution of (1.1), u=(u1,,un), satisfies

(1.7) lim x z x Ω u i ( x ) [ dist ( x , Ω ) ] - α i ( z ) = A i ( z ) , 1 i n ,

where

(1.8) α i ( z ) := { μ i ( z ) if  i I + I 0 , μ 1 ( z ) + γ i ( z ) p i if  i I - ,

and

(1.9) A i ( z ) := { ( μ i ( z ) ( μ i ( z ) + 1 ) 𝔟 i ( z ) ) 1 p i - 1 if  i I + , [ 1 𝔟 i ( z ) j = 1 k a i j ( μ 1 ( z ) ( μ 1 ( z ) + 1 ) 𝔟 j ( z ) ) 1 p j - 1 ] 1 p i if  i I - , A 0 , i if  i I 0 ,

where A0,i stands for the unique positive solution of the equation

𝔟 i ( z ) x p i - μ i ( z ) ( μ i ( z ) + 1 ) x = j = 1 k a i j ( μ 1 ( z ) ( μ 1 ( z ) + 1 ) 𝔟 j ( z ) ) 1 p j - 1 .

Essentially, Theorem 1.1 establishes that the values of the blow-up rates at z of any solution of (1.1) only depend on the precise way the blow-up rates at z of the solution of (1.3) are interrelated, rather than on the size of the coupling coefficients of the problem, aij, ij, as one might have expected from the every beginning. Indeed, by Theorem 1.1, the coupling coefficients only alter the values of Ai(z). In particular, (1.8) provides us with the components of (1.1) whose blow-up rates equal the corresponding blow-up rates of (1.3), as well as with the blow-up rates of the remaining components. Note that, although I0 and I- might be empty, owing to (1.5), 1I+. Hence, by Theorem 1.1, the first component of (1.1) blows up with exactly the same rate as the first component of the uncoupled problem (1.3), i.e., α1(z)=μ1(z). Moreover, μ1(z) plays a significant role in the blow-up rates of αi(z) for all iI-.

Remark 1.2.

The assignations Ω[0,+), zαi(z), are uniformly continuous for all 1in.

From Theorem 1.1 the next result follows readily.

Theorem 1.3.

Problem (1.1) admits a unique positive solution.

In the special case n=3, one can easily sketch how the blow-up rates given by (1.8) depend on the blow-up rates given by (1.4). By regarding

μ := ( μ 1 , μ 2 , μ 3 ) [ 0 , + ) 3

as an independent variable, we can divide [0,+)3 into several portions, taking into account the relationships between the components of the variable μ. Figure 1 shows a partition of the set of values of the parameters, μ[0,+)3, into thirteen complementary zones according to the nature of the values of the blow-up rates of the solutions of (1.1),

α i := α i ( μ ) , i = 1 , 2 , 3 ,

depending on μ.

The central portion of Figure 1 stands for the closed hexagonal prism

μ i + 2 - μ j 0 for all  i , j { 1 , 2 , 3 } , i j ,

where, according to Theorem 1.1,

α := ( α 1 , α 2 , α 3 ) = ( μ 1 , μ 2 , μ 3 ) .

Hence, it consists of the set of values of μ where all the blow-up rates of the coupled cooperative problem (1.1) equal the corresponding blow-up rates of the uncoupled one (1.3). Set

γ i := γ i ( μ ) = μ i ( p i - 1 ) - 2 , 1 i n .

Table 1 collects all the values of the blow-up rates α in each of the remaining zones according to Theorem 1.1.

Figure 1 
          The thirteen zones where α changes its nature,depending on μ.
Figure 1

The thirteen zones where α changes its nature,depending on μ.

Table 1

Blow-up rates α versus blow-up rates μ.

Order Zone α 1 α 2 α 3
μ 2 μ 3 μ 1 1 μ 1 μ 1 + γ 2 p 2 μ 3
μ 2 μ 3 μ 1 2 μ 1 μ 1 + γ 2 p 2 μ 1 + γ 3 p 3
μ 3 μ 2 μ 1 3 μ 1 μ 2 μ 1 + γ 3 p 3
μ 3 μ 2 μ 1 4 μ 1 μ 1 + γ 2 p 2 μ 1 + γ 3 p 3
μ 3 μ 1 μ 2 5 μ 1 μ 2 μ 2 + γ 3 p 3
μ 3 μ 1 μ 2 6 μ 2 + γ 1 p 1 μ 2 μ 2 + γ 3 p 3
μ 1 μ 3 μ 2 7 μ 2 + γ 1 p 1 μ 2 μ 3
μ 1 μ 3 μ 2 8 μ 2 + γ 1 p 1 μ 2 μ 2 + γ 3 p 3
μ 1 μ 2 μ 3 9 μ 3 + γ 1 p 1 μ 2 μ 3
μ 1 μ 2 μ 3 10 μ 3 + γ 1 p 1 μ 3 + γ 2 p 2 μ 3
μ 2 μ 1 μ 3 11 μ 1 μ 3 + γ 2 p 2 μ 3
μ 2 μ 1 μ 3 12 μ 3 + γ 1 p 1 μ 3 + γ 2 p 2 μ 3

Surrounding the central prisma and labeled by the first six odd integers, we have represented the set of values of μ for which two of the blow-up rates αi equal the corresponding blow-up rates μi. Labeled by the first six even integers, we find the set of values of μ for which exactly one of the αi equals the corresponding μi. As it has been already commented after the statement of Theorem 1.1, there is not any value of μ for which the three blow-up rates αi can differ from the corresponding blow-up rates μi.

The results of this paper can be easily generalized to cover the case where the coefficients aij are replaced by positive Hölder continuous functions, aij𝒞ν(Ω), 1in. Also in this case, the blow-up rates of the solution of (1.1) at zΩ are given by (1.8), while the coefficients Ai(z) are given by (1.9), by replacing aij by aij(z) for all 1i,jn. Moreover, under the appropriate circumstances, the γi(z), introduced in (1.2), might satisfy γi(z)>-2 without affecting substantially most of the results. It should be noted that power-like nonlinearities, as those considered in this paper, always satisfy the classical conditions of Keller [17] and Osserman [39].

The structure of this paper is the following. Section 2 collects the necessary results to obtain existence of solutions for (1.1) as well as the necessary comparison results for the proofs of Theorems 1.1 and 1.3. Section 3 provides us with the construction of a supersolution for the singular problem in a ball and of a subsolution for the problem in an annulus. Finally, in Section 4 we complete the proofs of Theorems 1.1 and 1.3, and the proof of Remark 1.2.

2 Existence of Large Solutions

The main goal of this section is to sketch the proof of the existence of solution for (1.1). Some previous results were already found in [2] and [28]. Given a smooth subdomain DN, we consider the operator 𝔏:[𝒞2+ν(D)]n[𝒞ν(D)]n defined by

( 𝔏 u ) i = - Δ u i - λ i ( ) u i - j = 1 , j i n a i j u j , 1 i n .

By the main result of [30], there exists a unique σ, the principal eigenvalue of 𝔏 under Dirichlet homogeneous conditions, such that the linear eigenvalue problem

𝔏 φ = σ φ in  D , φ = 0 on  D ,

admits a positive eigenfunction φ[𝒞2+ν(D)]n, φi0, 1in, φ0. Such a value of σ will be denoted by σ[𝔏,D].

Throughout the rest of this paper, for any given u[𝒞(D)]n, it is said that u>0 in D if ui0 for all 1in but u0. Similarly, given u[𝒞1(D¯)]n, it is said that u is strongly positive in D, u0, if for any 1in,

u i ( x ) > 0 for all  x D    and    u n x ( x ) < 0 for all  x u - 1 ( 0 ) D ,

where nx is the outward unit normal vector to D at xD.

The following characterization of the strong maximum principle, going back to [30], holds.

Theorem 2.1.

The following assertions are equivalent:

  1. σ [ 𝔏 , D ] > 0 .

  2. There exists u ¯ [ 𝒞 2 ( D ) ] n [ 𝒞 ( D ¯ ) ] n such that u ¯ > 0 in D and

    𝔏 u ¯ 0 in  D ,

    and, for some 1 i 0 n , either u ¯ i 0 > 0 on D , or else

    ( 𝔏 u ¯ ) i 0 > 0 in  D .

    Should this be the case, u ¯ is said to be a positive strict supersolution of 𝔏 in D.

  3. The operator 𝔏 satisfies the strong maximum principle in D , in the sense that, for every h [ 𝒞 ν ( D ¯ ) ] n , u[𝒞2+ν(D¯)]n and w[𝒞2+ν(D)]n satisfying

    𝔏 u = h 0 in  D , u = w 0 on  D ,

    with some of these inequalities strict, one has that u 0 in D.

Using Theorem 2.1, one can easily show the monotonicity of the principal eigenvalue with respect to the potentials λi() and the coefficients aij, 1i,jn. Actually, this result was established by [30, Theorem 3.2]. As a result, if we assume that

λ ¯ i ( x ) λ ¯ i ( x ) and a ¯ i j a ¯ i j for all  x Ω ¯  and  1 i , j n ,

with some of these inequalities strict, then, setting

( 𝔏 ¯ u ) i := - Δ u i - λ ¯ i ( ) u i - j = 1 , j i n a ¯ i j u j , ( 𝔏 ¯ u ) i := - Δ u i - λ ¯ i ( ) u i - j = 1 , j i n a ¯ i j u j ,

we find that

(2.1) σ [ 𝔏 ¯ , D ] > σ [ 𝔏 ¯ , D ] .

Next, for every w[𝒞2+ν(Ω)]n, w>0, we consider the non-homogeneous Dirichlet boundary value problem

(2.2) { - Δ u i = λ i ( x ) u i + j = 1 , j i n a i j u j - 𝔞 i ( x ) u i p i in  Ω , u i = w i on  Ω ,  1 i n .

Using Theorem 2.1, the results of [35, 34, 36] can be easily adapted to obtain the next one.

Theorem 2.2.

Suppose (2.2) admits a subsolution u¯[C2+ν(Ω¯)]n and a supersolution u¯[C2+ν(Ω¯)]n satisfying u¯u¯. Then (2.2) possesses a solution u[C2+ν(Ω¯)]n such that u¯uu¯. Actually, (2.2) possesses a minimal and a maximal solution in the interval [u¯,u¯].

Using Theorems 2.1 and 2.2, the abstract results of [2, Section 3] can be easily adapted, almost mutatis mutandis, to get the next one.

Theorem 2.3.

Problem (2.2) has a unique positive solution, θ[Ω,w]. Moreover, for every positive subsolution u¯ (resp. supersolution u¯) of (2.2),

u ¯ θ [ Ω , w ] ( resp.  u ¯ θ [ Ω , w ] ) .

Proof.

Adapting [2, Theorem 3.7], it suffices to construct a supersolution u¯ such that u¯i>0 for all 1in. In the special case 𝔞i(x)>0 for all xΩ and 1in, one can take u¯=(M,,M) for some M>0 sufficiently large. In general, we may proceed as follows.

Since Ω is smooth, Ω possesses finitely many components, Γk, 1km. For each ε>0 and 1km, denote

Ω k ε := { x N : dist ( x , Γ k ) < ε } .

Let

λ := max 1 i n λ i + 1 , a := max 1 i , j n a i j ,

and let 𝔏¯ be the operator

( 𝔏 ¯ u ) i := - Δ u i - λ u i - j = 1 , j i n a u j .

Thanks to (2.1),

σ [ 𝔏 , Ω k ε ] > σ [ 𝔏 ¯ , Ω k ε ] for all  ε > 0  and  1 k m .

On the other hand, by the uniqueness of the principal eigenvalue,

σ [ 𝔏 ¯ , Ω k ε ] = σ [ - Δ - λ - ( n - 1 ) a , Ω k ε ] = σ [ - Δ , Ω k ε ] - λ - ( n - 1 ) a , 1 k m .

Thus, since

lim ε 0 | Ω k ε | = 0 , 1 k m ,

where |Ωkε| stands for the Lebesgue measure of Ωkε, the Faber–Krahn inequality, going back to [11] and [18], yields

lim ε 0 σ [ 𝔏 , Ω k ε ] = + , 1 k m

(see [21, Theorem 5.1]). Therefore, ε can be shortened, if necessary, so that min1kmσ[𝔏,Ωkε]>0. Fix ε>0 satisfying the last inequality and, for each 1km, let

φ k := ( φ k , 1 , , φ k , n )

be a principal eigenfunction associated to σ[𝔏,Ωkε]. As φk0, it is apparent that

(2.3) min { min Γ k φ k , i : 1 i n } > 0 , min { min Ω Ω k ε / 2 φ k , i : 1 i n } > 0 for all  1 k m .

Subsequently, we consider the auxiliary function Φ defined through

Φ := { φ k in  Ω ¯ m ε / 2 ,  1 k m , g in  Ω int := Ω ( k = 1 m Ω ¯ k ε / 2 ) ,

where g is any C2+ν-extension of the function φ1φm to the open set Ωint with the special requirement that infΩintg>0. Such a function exists because of (2.3). Then τΦ is a supersolution of (2.2) for sufficiently large τ>1. Indeed, by (2.3), there exists τ01 such that

τ Φ > w on  Ω  for all  τ > τ 0 .

Moreover, for every 1km and 1in, we find that, in Ωkε/2Ω,

- Δ ( τ Φ i ) - λ i ( ) τ Φ i - j = 1 , j i n a i j τ Φ j = ( 𝔏 ( τ Φ ) ) i = τ ( 𝔏 ( φ k ) ) i = τ σ [ 𝔏 , Ω k ε ] φ k , i > 0 - 𝔞 i ( ) ( τ Φ i ) p i .

Lastly, in Ωint, we have that, for every 1in,

- Δ ( τ Φ i ) = - τ Δ g i λ i ( ) g i + j = 1 , j i n a i j τ g j - 𝔞 i ( ) ( τ g i ) p i

for sufficiently large τ>1, because 𝔞i and gi are bounded away from zero in Ωint and pi>1 for all 1in. This ends the proof. ∎

Theorem 2.3 entails that the mapping

( 0 , + ) [ 𝒞 2 + ν ( Ω ¯ ) ] n , m θ [ Ω , m ] ,

where m:=(m,,m), is strongly increasing, in the sense that m1<m2 implies θ[λ,Ω,m1]θ[λ,Ω,m2]. Hence, the point-wise limit

(2.4) θ [ Ω , ] ( x ) := lim m + θ [ Ω , m ] ( x ) , x Ω ,

is well defined. In fact, the next result holds.

Theorem 2.4.

There exists a minimal and a maximal positive solution of (1.1), Lmin and Lmax, respectively, in the sense that any solution, L, of (1.1) satisfies

L min ( x ) L ( x ) L max ( x ) x Ω .

Moreover, the point-wise limit (2.4) provides us with the minimal solution

L min = θ [ Ω , ] ,

while the maximal solution is given by

L max = lim δ 0 θ [ Ω δ , ] ,

where we have denoted

Ω δ := { x Ω : d ( x , Ω ) > δ } , δ > 0 .

3 Two Pivotal Technical Results Under Radial Symmetry

In this section, for every R>0, we consider the auxiliary problem

(3.1) { - ψ i ′′ - N - 1 r ψ i = λ ψ i + j = 1 , j i n a i j ψ j - b i ( R - r ) γ i ψ i p i , 0 < r < R , ψ i ( 0 ) = 0 , ψ i ( R ) = + , 1 i n ,

where λ, γi0, pi>1, aij>0 and bi>0 for all 1i,jn. Without loss of generality, we can assume that the equations in (3.1) have been reordered so that

0 < μ n μ n - 1 μ 1 ,

where, as in (1.4), for each 1in, μi is defined by

μ i := γ i + 2 p i - 1 .

As in (1.6), we consider

(3.2) { I + := { i { 1 , , n } : μ i + 2 - μ 1 > 0 } , I 0 := { i { 1 , , n } : μ i + 2 - μ 1 = 0 } , I - := { i { 1 , , n } : μ i + 2 - μ 1 < 0 } .

Let k1 be such that

I M := { i { 1 , , n } : μ i = μ 1 } = { 1 , , k } ,

and set

(3.3) α i := { μ i if  i I + I 0 , μ 1 + γ i p i if  i I - ,

and

(3.4) A ¯ i := { ( μ i ( μ i + 1 ) b i ) 1 p i - 1 if  i I + , [ 1 b i j = 1 k a i j ( μ 1 ( μ 1 + 1 ) b j ) 1 p j - 1 ] 1 p i if  i I - , A ¯ 0 , i if  i I 0 ,

where A¯0,i stands for the unique positive solution of

b i x p i - μ i ( μ i + 1 ) x = j = 1 k a i j ( μ 1 ( μ 1 + 1 ) b j ) 1 p j - 1 .

The first result of this section reads as follows.

Lemma 3.1.

For each ε>0, there exists a constant C:=C(ε,n) such that the function ψ¯ε:=(ψ¯ε,1,,ψ¯ε,n), defined by

ψ ¯ ε , i ( r ) := ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C , 0 r < R ,  1 i n ,

provides us with a supersolution of (3.1).

Proof.

By definition, ψ¯ε is smooth and satisfies the boundary conditions. Hence, ψ¯ε is a supersolution of (3.1) if, and only if,

- ( r R ) 2 ( 1 + ε ) A ¯ i α i ( α i + 1 ) ( R - r ) - α i - 2 - ( N + 3 ) r R 2 ( 1 + ε ) A ¯ i α i ( R - r ) - α i - 1 - 2 N R 2 ( 1 + ε ) A ¯ i ( R - r ) - α i
λ [ ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C ] + j = 1 , j i n a i j [ ( 1 + ε ) A ¯ j ( R - r ) - α j ( r R ) 2 + C ]
(3.5)     - b i ( R - r ) γ i [ ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C ] p i

for every 0<r<R and 0in. Now, multiply (3.5) by (R-r)αi+2 if iI+I0, and by (R-r)-γi+αipi if iI-. Then, ψ¯ε is a supersolution if, and only if, for every 0<r<R,

- ( r R ) 2 ( 1 + ε ) A ¯ i α i ( α i + 1 ) - ( N + 3 ) r R 2 ( 1 + ε ) A ¯ i α i ( R - r ) - 2 N R 2 ( 1 + ε ) A ¯ i ( R - r ) 2
λ [ ( 1 + ε ) A ¯ i ( R - r ) 2 ( r R ) 2 + C ( R - r ) α i + 2 ] + j = 1 , j i n a i j [ ( 1 + ε ) A ¯ j ( R - r ) α i + 2 - α j ( r R ) 2 + C ( R - r ) α i + 2 ]
(3.6)     - b i ( R - r ) α i + 2 + γ i [ ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C ] p i

if iI+I0, and

- ( r R ) 2 ( 1 + ε ) A ¯ i α i ( α i + 1 ) ( R - r ) - γ i + α i p i - α i - 2 - ( N + 3 ) r R 2 ( 1 + ε ) A ¯ i α i ( R - r ) - γ i + α i p i - α i - 1
- 2 N R 2 ( 1 + ε ) A ¯ i ( R - r ) - γ i + α i p i - α i
λ [ ( 1 + ε ) A ¯ i ( R - r ) - γ i + α i p i - α i ( r R ) 2 + C ( R - r ) - γ i + α i p i ]
+ j = 1 , j i n a i j [ ( 1 + ε ) A ¯ j ( R - r ) - γ i + α i p i - α j ( r R ) 2 + C ( R - r ) - γ i + α i p i ]
(3.7) - b i ( R - r ) α i p i [ ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C ] p i

if iI-.

Let us show that

(3.8) α i + 2 - α j > 0 for all  i I +  and  1 j n , j i .

Indeed, if iI+ and jI+I0, ji, we have that

α i + 2 - α j = μ i + 2 - μ j μ i + 2 - μ 1 > 0 ,

by the definition of I+ (see (3.2)). When jI-, we may proceed as follows. According to (3.2),

γ j + 2 p j - 1 < μ 1 - 2 , j I - .

So, multiplying by pj-1 we deduce that γj<pjμ1-2pj-μ1, jI-, and dividing by pj yields

γ j p j + μ 1 p j < μ 1 - 2 , j I - .

Equivalently,

(3.9) α j < μ 1 - 2 , j I - ,

and therefore

α i + 2 - α j = μ i + 2 - α j > μ i + 2 - μ 1 + 2 > 2 > 0

for every iI+ and jI-.

Analogously, the following estimates hold:

(3.10) { α i + 2 - α j = 0 for every  i I 0 ,  1 j k , α i + 2 - α j > 0 for every  i I 0 , k + 1 j n , j i .

Indeed, by (3.3), αi=μi for all iI0. Similarly, since {1,,k}I+, we have αj=μj for all 1jk. Moreover, μj=μ1 for all 1jk. Thus,

α i + 2 - α j = μ i + 2 - μ j = μ i + 2 - μ 1 = 0

for all iI0 and 1jk, which shows the validity of the identities of (3.10). In order to check the inequalities of (3.10), we can argue as follows. Pick iI0 and j{k+1,,n}, ji. Suppose jI+I0. Then, by (3.3) and taking into account that, by construction, μj<μ1 for all k+1jn, we find that

α i + 2 - α j = μ i + 2 - μ j > μ i + 2 - μ 1 = 0 ,

because iI0. Now, suppose that jI-. Then, thanks to (3.3) and (3.9),

α i + 2 - α j = μ i + 2 - α j > μ i + 2 - μ 1 + 2 = 2 > 0 .

Therefore, (3.10) holds.

Next, we will see that

(3.11) - γ i + α i p i - α i - 2 > 0 for all  i I - .

Indeed, by the definition of μi and since iI-,

γ i + 2 p i - 1 + 2 < μ 1 , i I - .

Thus, adding γi at both sides of this inequality and taking common factor γi+2 yields

( γ i + 2 ) p i p i - 1 < μ 1 + γ i , i I - .

Hence, by (3.3),

γ i + 2 p i - 1 < μ 1 + γ i p i = α i , i I - ,

whence (3.11).

Lastly, we will establish that

(3.12) { - γ i + α i p i - α j = 0 for all  i I - ,  1 j k , - γ i + α i p i - α j > 0 for all  i I - , k + 1 j n , j i .

By (3.3), -γi+αipi=μ1 for all iI-. Thus,

- γ i + α i p i - α j = μ 1 - α j = μ 1 - μ j = 0

for all iI- and 1jk, whence the identities of (3.12). Now, pick iI- and k+1jn, ij. Suppose jI+I0. Then, due to (3.3),

- γ i + α i p i - α j = μ 1 - α j = μ 1 - μ j > 0 ,

and hence (3.12) holds in this case. Now, suppose jI-. Then, by (3.9),

- γ i + α i p i - α j = μ 1 - α j > 2 > 0 ,

and therefore (3.12) is satisfied.

By (3.3) and the definition of μi,

α i + 2 + γ i - α i p i = μ i + 2 + γ i - μ i p i = 0 for all  i I + I 0 .

Thus,

(3.13) lim r R ( b i ( R - r ) α i + 2 + γ i [ ( 1 + ε ) A ¯ i ( R - r ) - α i ( r R ) 2 + C ] p i ) = b i ( 1 + ε ) p i A ¯ i p i

for all C0 and iI+I0. Hence, thanks to (3.8), (3.10) and (3.13), we can extend (3.6) to r=R by letting rR. Similarly, (3.7) can be extended to r=R by (3.11) and (3.12). More precisely, at r=R, (3.6) provides us with

(3.14) { - ( 1 + ε ) A ¯ i α i ( α i + 1 ) - b i ( 1 + ε ) p i A ¯ i p i , i I + , - ( 1 + ε ) A ¯ i α i ( α i + 1 ) j = 1 k a i j ( 1 + ε ) A ¯ j - b i ( 1 + ε ) p i A ¯ i p i , i I 0 ,

while (3.7) at r=R becomes

(3.15) 0 j = 1 k a i j ( 1 + ε ) A ¯ j - b i ( 1 + ε ) p i A ¯ i p i , i I - .

Due to (3.4) and using that (1+ε)<(1+ε)pi (since pi>1 for all 1in), it is easily seen that (3.14) and (3.15) are true. In the derivation one should note that αj=μj=μ1 for all 1jk, by construction. Actually, all inequalities in (3.14) and (3.15) are strict, because (1+ε)<(1+ε)pi. By continuity, this entails the existence of δ:=δ(ε,n)>0 such that (3.6) and (3.7) are satisfied for all r[R-δ,R). Therefore, choosing a sufficiently large C>0, we can assume that (3.5) holds in (0,R), because, for each 1in, the function bi(R-r)γi is positive and bounded away from zero in [0,R-δ]. The proof is complete. ∎

The next result provides us with a universal subsolution on an annulus.

Lemma 3.2.

Let R2>R1>0 and consider the problem

(3.16) { - ψ i ′′ - N - 1 r ψ i = λ ψ i + j = 1 , j i n a i j ψ j - β i ( r ) ( r - R 1 ) γ i ψ i p i , R 1 < r < R 2 , ψ i ( R 1 ) = + , ψ i ( R 2 ) = 0 , 1 i n ,

where all the coefficients satisfy the same requirements as in (3.1) and, for every 1in, the function βiCν[R1,R2] satisfies βi(r)>0 for all r[R1,R2]. Then, for every ε(0,1), there exists a negative constant D:=D(ε,n)<0 such that the function

ψ ¯ ε := ( ψ ¯ ε , 1 , , ψ ¯ ε , n ) ,

defined by

ψ ¯ ε , i ( r ) := max { 0 , ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D } , R 1 < r R 2 ,  1 i n ,

provides us with a weak subsolution of (3.1), as discussed in [1], where the constants αi, 1in, are given by (3.3) and

(3.17) A ¯ i := { ( μ i ( μ i + 1 ) β i ( R 1 ) ) 1 p i - 1 if  i I + , [ 1 β i ( R 1 ) j = 1 k a i j ( μ 1 ( μ 1 + 1 ) β j ( R 1 ) ) 1 p j - 1 ] 1 p i if  i I - , A ¯ 0 , i if  i I 0 ,

where A¯0,i stands for the unique positive solution of

β i ( R 1 ) x p i - μ i ( μ i + 1 ) x = j = 1 k a i j ( μ 1 ( μ 1 + 1 ) β j ( R 1 ) ) 1 p j - 1 .

Proof.

As the maps r(1-ε)A¯i(r-R1)-αi are strictly decreasing, taking any D satisfying

D < - ( 1 - ε ) A ¯ i ( R 2 - R 1 ) - α i = - min { ( 1 - ε ) A ¯ i ( r - R 1 ) - α i : R 1 < r R 2 }

for all 1in, there exist ϱi(D)(R1,R2), 1in, such that

ψ ¯ ε , i = { ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D if  R 1 < r ϱ i ( D ) , 0 if  ϱ i ( D ) < r R 2 ,  1 i n .

Moreover, the mappings Dϱi(D), 1in, can be chosen continuous, and

(3.18) lim D - ϱ i ( D ) = R 1 , 1 i n .

Thus, ψ¯ε is a subsolution of (3.16) if, and only if, for every 1in and R1<rϱi(D),

- ( 1 - ε ) A ¯ i α i ( α i + 1 ) ( r - R 1 ) - α i - 2 + N - 1 r ( 1 - ε ) A ¯ i α i ( r - R 1 ) - α i - 1
λ [ ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D ] + j = 1 , j i n a i j [ ( 1 - ε ) A ¯ j ( r - R 1 ) - α j + D ]
(3.19)     - β i ( r ) ( r - R 1 ) γ i [ ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D ] p i .

Next, we will adapt the proof of Lemma 3.1. Multiplying (3.19) by (r-R1)αi+2 when iI+I0 and by (r-R1)-γi+αipi when iI-, it becomes apparent that ψ¯ε is a subsolution of (3.16) if, and only if,

- ( 1 - ε ) A ¯ i α i ( α i + 1 ) + N - 1 r ( 1 - ε ) A ¯ i α i ( r - R 1 )
λ [ ( 1 - ε ) A ¯ i ( r - R 1 ) 2 + D ( r - R 1 ) α i + 2 ] j = 1 , j i n a i j [ ( 1 - ε ) A ¯ j ( r - R 1 ) α i + 2 - α j + D ( r - R 1 ) α i + 2 ]
(3.20)     - β i ( r ) ( r - R 1 ) α i + 2 + γ i [ ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D ] p i

for all iI+I0 and R1<rϱi(D), and

- ( 1 - ε ) A ¯ i α i ( α i + 1 ) ( r - R 1 ) - α i - 2 - γ i + α i p i + N - 1 r ( 1 - ε ) A ¯ i α i ( r - R 1 ) - α i - 1 - γ i + α i p i
λ [ ( 1 - ε ) A ¯ i ( r - R 1 ) - α i - γ i + α i p i + D ( r - R 1 ) - γ i + α i p i ]
    j = 1 , j i n a i j [ ( 1 - ε ) A ¯ j ( r - R 1 ) - γ i + α i p i - α j + D ( r - R 1 ) - γ i + α i p i ]
(3.21)     - β i ( r ) ( r - R 1 ) α i p i [ ( 1 - ε ) A ¯ i ( r - R 1 ) - α i + D ] p i

for all iI- and R1<rϱi(D). Thanks to (3.8), (3.10), (3.11) and (3.12), letting rR1 in (3.20) and (3.21) yields

(3.22) { - ( 1 - ε ) A ¯ i α i ( α i + 1 ) - β i ( R 1 ) ( 1 - ε ) p i A ¯ i p i , i I + , - ( 1 - ε ) A ¯ i α i ( α i + 1 ) j = 1 k a i j ( 1 - ε ) A ¯ j - β i ( R 1 ) ( 1 - ε ) p i A ¯ i p i , i I 0 , 0 j = 1 k a i j ( 1 - ε ) A ¯ j - β i ( R 1 ) ( 1 - ε ) p i A ¯ i p i , i I - .

As (1-ε)>(1-ε)pi for all 1in, all inequalities in (3.22) are strict. Hence, by continuity, there exists δ=δ(ε,n)>0 such that all inequalities in (3.20) and (3.21) hold in the interval (R1,R1+δ]. Finally, by (3.18), D<0 can be taken sufficiently negative so that

max 1 i n ϱ i ( D ) R 1 + δ .

This ends the proof. ∎

4 Proofs of the Main Results

As Ω is smooth, the outward unit normal vector field to Ω is well defined at every point of Ω. We will denote it by

n : Ω N , z n z .

Since Ω𝒞2, Ω satisfies the uniform interior sphere in the strong sense onΩ (see [26, Theorem 1.9]). So, there exists R0>0 such that for every xΩ with dist(x,Ω)<R0, there is a point π(x)Ω such that

(4.1) t ( x ) := dist ( x , Ω ) = | x - π ( x ) | , B R 0 ( π ( x ) + R 0 x - π ( x ) t ( x ) ) Ω .

Moreover, R0 can be shortened so that for every zΩ,

(4.2) B ¯ R 0 ( z - R 0 n z ) Ω = { z } and B ¯ R 0 ( z + R 0 n z ) Ω ¯ = { z } .

4.1 Proof of Theorem 1.1

Fix zΩ and 0<η<1. By (1.2), there exist δ>0 such that for every 1in,

(4.3) ( 1 - η ) 𝔟 i ( z ) [ dist ( x , Ω ) ] γ i ( z ) < 𝔞 i ( x ) < ( 1 + η ) 𝔟 i ( z ) [ dist ( x , Ω ) ] γ i ( z )

for all xBδ(z)Ω. Choose R0 sufficiently small so that (4.1) and (4.2) hold, and 0<R0<δ. Set

Γ := B ¯ R 0 / 2 ( z ) Ω .

It is rather clear that there exist R>0 and ϱ0>0 such that

B R ( y - ( R + ϱ ) n y ) B δ ( z ) Ω

for all yΓ and 0ϱϱ0. Figure 2 sketches this construction scheme.

According to (4.3), for every yΓ and 0ϱϱ0,

(4.4) 𝔞 i ( x ) > ( 1 - η ) 𝔟 i ( z ) [ dist ( x , Ω ] γ i ( z ) ( 1 - η ) 𝔟 i ( z ) [ dist ( x , B R ( y - ( R + ϱ ) n y ) ) ] γ i ( z ) , 1 i n ,

for all xBR(y-(R+ϱ)ny). Set

λ ¯ := max { λ i : 1 i n } ,
y ϱ := y - ( R + ϱ ) n y , y Γ , ϱ [ 0 , ϱ 0 ] ,
b ¯ i ( z ) := ( 1 - η ) 𝔟 i ( z ) , 1 i n ,

Figure 2 
            Scheme of the construction.
Figure 2

Scheme of the construction.

and consider, for each yΓ and 0<ϱϱ0, the problem

(4.5) { - Δ u i = λ ¯ u i + j = 1 , j i n a i j u j - b ¯ i ( z ) ( R - | x - y ϱ | ) γ i ( z ) u i p i in  B R ( y ϱ ) , u i = + on  B R ( y ϱ ) ,  1 i n .

By (4.4), any positive solution of (1.1),

L = ( L 1 , , L n ) := ( u 1 , , u n ) ,

is a bounded positive subsolution of (4.5) for every yΓ and 0<ϱϱ0.

Let ε>0. Applying Lemma 3.1 to problem (3.1) with λ=λ¯, γi=γi(z), bi=b¯i(z), in, we get the functions

ψ ¯ ε , i ( r ) = ( 1 + ε ) A ¯ i ( z ) ( R - r ) - α i ( z ) ( r R ) 2 + C , 0 r < R ,  1 i n ,

where αi(z), A¯i(z), 1in, are defined through (3.3) and (3.4). By the radial symmetry of (4.5), for every yΓ and 0<ϱϱ0, the functions

L ¯ ε , i y ϱ ( x ) := ψ ¯ ε , i ( r ) , x B R ( y ϱ ) , r := | x - y ϱ | ,  1 i n ,

provide us with a supersolution of (4.5). Hence, by Theorem 2.3,

L i ( x ) L ¯ ε , i y ϱ ( x ) , x B R ( y ϱ ) ,  1 i n ,

for every yΓ and 0<ϱϱ0. Consequently, we may infer

(4.6) L i ( x ) ( 1 + ε ) A ¯ i ( z ) [ dist ( x , B R ( y 0 ) ) ] - α i ( z ) ( | x - y 0 | R ) 2 + C

for all yΓ, xBR(y0) and 1in.

On the other hand, for every x sufficiently close to Γ, we have that dist(x,Γ)=dist(x,BR(y0)), with y0=π(x)-Rnπ(x). Thus, (4.6) implies

lim sup dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) ( 1 + ε ) A ¯ i ( z ) , 1 i n .

Therefore, letting ε0 yields

(4.7) lim sup dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) A ¯ i ( z ) , 1 i n .

Figure 3 
            Scheme of the construction.
Figure 3

Scheme of the construction.

Figure 4 
            Magnification around y0{y^{0}}.
Figure 4

Magnification around y0.

Now, we will construct a subsolution of (1.1) with the appropriate growth on Γ. Due to (4.2) and taking into account that Ω is bounded, there exist R2>R1>0 and ϱ0>0 such that

Ω 0 ϱ ϱ 0 A R 1 , R 2 ( y + ( R 1 + ϱ ) n y ) and Ω ¯ A R 1 , R 2 ( y + R 1 n y ) = { y }

for all yΓ and ϱ[0,ϱ0], where, for every zN and 0<r1<r2,

A r 1 , r 2 ( z ) := { x N : r 1 < | x - z | < r 2 } .

As before, we will denote

y ϱ := y + ( R 1 + ϱ ) n y , y Γ ,  0 ϱ ϱ 0 .

Figure 3 sketches this construction.

Shortening R1 and ϱ0 if necessary, it becomes apparent that there exists ξ>0 such that

(4.8) A R 1 , R 1 + ξ ( y ϱ ) Ω ( B δ ( z ) Ω ) for all  y Γ  and  ϱ [ 0 , ϱ 0 ] .

Figure 4 shows a magnification of Figure 3 in a neighborhood of y0. Using again (4.3) and thanks to (4.8), we obtain that

(4.9) 𝔞 i ( x ) < ( 1 + η ) 𝔟 i ( z ) [ dist ( x , Ω ) ] γ i ( z ) ( 1 + η ) 𝔟 i ( z ) [ dist ( x , B R 1 ( y ϱ ) ) ] γ i ( z )

for all yΓ, ϱ[0,ϱ0], xAR1,R1+ξ(y0)Ω, 1in. For getting (4.9) in the entire Ω, we may proceed as follows. Set

K := max 1 i n { 𝔞 i ( R 1 + ξ ) γ i ( z ) + ( 1 + η ) 𝔟 i ( z ) } ,

and consider, for each 1in, the piecewise linear function

β i ( r ) = { ( 1 + η ) 𝔟 i ( z ) , R 1 r < R 1 + ξ 2 , K + K - ( 1 + η ) 𝔟 i ( z ) ξ / 2 ( r - ( R 1 + ξ ) ) , R 1 + ξ 2 r < R 1 + ξ , K , R 1 + ξ r R 2 .

Using the definition of βi and (4.9), we have that

(4.10) 𝔞 i ( x ) β i ( | x - y ϱ | ) [ | x - y ϱ | - R 1 ] γ i ( z ) , 1 i n ,

for all yΓ, 0<ϱϱ0 and xΩ.

Subsequently, we will consider the auxiliary problem

(4.11) { - Δ u i = λ ¯ u i + j = 1 , j i n a i j u j - β i ( r ) ( r - R 1 ) γ i ( z ) u i p i in  A R 1 , R 2 ( y ϱ ) , u i = + on  A R 1 , R 2 ( y ϱ ) ,  1 i n ,

where

λ ¯ := min { - λ i : 1 i n } , r := | x - y ϱ | .

Pick ε(0,1). Then, due to Lemma 3.2, for every yΓ and 0<ϱϱ0, the function Lεyϱ:=(Lε,1yϱ,,Lε,nyϱ), defined for every 1in by

L ε , i y ϱ ( x ) := ψ ¯ ε , i ( r ) = max { 0 , ( 1 - ε ) A ¯ i ( z ) ( r - R 1 ) - α i ( z ) + D } , x A R 1 , R 2 ( y ϱ ) ,

provide us with a positive subsolution for (4.11), where A¯i(z) is defined through (3.17) for every 1in.

By (4.10), the restriction Lεyϱ|Ω provides us with a bounded positive subsolution of (1.1). Hence, by Theorem 2.3,

L i ( x ) L ¯ ε , i y ϱ ( x ) , x Ω ,  1 i n ,

for all yΓ and 0<ϱϱ0. Thus, we can infer that

L i L ¯ ε , i y 0 , 1 i n ,

for all yΓ. The last inequality implies that

lim inf dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) ( 1 - ε ) A ¯ i ( z ) , 1 i n .

Hence, letting ε0 yields

(4.12) lim inf dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) A ¯ i ( z ) , 1 i n .

Consequently, owing to (4.7) and (4.12), for each zΩ and 0<η<1, there exists a compact neighborhood of zΩ, Γ, such that

A ¯ i ( z ) lim inf dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) lim sup dist ( x , Γ ) 0 L i ( x ) [ dist ( x , Γ ) ] - α i ( z ) A ¯ i ( z )

for every 1in. Therefore, letting η0 yields

lim x z x Ω L i ( x ) [ dist ( x , Ω ) ] - α i ( z ) = A i ( z ) ,

because

A i ( z ) = lim η 0 A ¯ i ( z ) = lim η 0 A ¯ i ( z ) .

This ends the proof of Theorem 1.1.

4.2 Proof of Theorem 1.3

Let L:=(L1,,Ln) and M:=(M1,,Mn) be two positive solutions of (1.1). Using (1.7) it is easily seen that the quotients

q i ( x ) := { L i ( x ) M i ( x ) x Ω , 1 x Ω ,  1 i n ,

are uniformly continuous in Ω¯. Thus, for every ε>0, there exists δ>0 such that

| q i ( x ) - q i ( π ( x ) ) | = | q i ( x ) - 1 | < ε if  | x - π ( x ) | < δ ,  1 i n .

Thus, setting

Q ξ := { x Ω ¯ : dist ( x , Ω ) ξ } , ξ > 0 ,

we find that

( 1 - ε ) M i L i ( 1 + ε ) M i in  Q δ ,  1 i n .

Moreover, L is a solution of the problem

(4.13) { - Δ u i = λ i ( x ) u i + j = 1 , j i n a i j u j - 𝔞 i ( x ) u i p i in  Ω Q δ , u i = L i on  ( Ω Q δ ) ,  1 i n ,

( 1 - ε ) M is a subsolution of (4.13) and (1+ε)M is a supersolution of (4.13). Therefore, by Theorem 2.3,

( 1 - ε ) M i L i ( 1 + ε ) M i in  Ω Q δ ,  1 i n .

Letting ε0 we find that L=M in Ω. This ends the proof.

4.3 Proof of Remark 1.2

Setting

μ max ( z ) := max 1 j n μ j ( z ) , z Ω ,

by Theorem 1.1, for every zΩ and 1in, we have

α i ( z ) = { μ i ( z ) if  μ i ( z ) + 2 - μ max ( z ) 0 , μ max ( z ) + γ i ( z ) p i if  μ i ( z ) + 2 - μ max ( z ) < 0 .

Pick zΩ and 1in. Based on the continuity of μmax(z), it easy to check that αi is continuous at z if

μ i ( z ) + 2 - μ max ( z ) 0 .

Suppose

(4.14) μ i ( z ) + 2 - μ max ( z ) = 0

for some zΩ, and let (zn)n1Ω be a sequence such that znz if n+ and

μ i ( z n ) + 2 - μ max ( z n ) < 0 for all  n 1 .

Then, invoking (4.14) and (1.4), shows that

lim n α i ( z n ) = lim n μ max ( z n ) + γ i ( z n ) p i = μ max ( z ) + γ i ( z ) p i
= μ i ( z ) + 2 + γ i ( z ) p i = γ i ( z ) + 2 p i - 1 + 2 + γ i ( z ) p i
= ( γ i ( z ) + 2 ) ( 1 p i - 1 + 1 ) p i = γ i ( z ) + 2 p i - 1
= α i ( z ) .

Therefore, αi is also continuous at zΩ if (4.14) holds. This ends the proof.


Communicated by Laurent Veron


Award Identifier / Grant number: MTM2012-30669

Award Identifier / Grant number: MTM2015-65899-P

Funding statement: Supported by the Ministry of Economy and Competitiveness of Spain under Research Grants MTM2012-30669 and MTM2015-65899-P.

References

[1] H. Amann, Maximum principle and principal eigenvalues, Ten Mathematical Essays on Approximation in Analysis and Topology, Elsevier, Amsterdam (2005), 1–60. 10.1016/B978-044451861-3/50001-XSearch in Google Scholar

[2] P. Álvarez-Caudevilla and J. López-Gómez, Metasolutions for cooperative systems, Nonlinear Anal. Real World Appl. 9 (2008), 1119–1157. 10.1016/j.nonrwa.2007.02.010Search in Google Scholar

[3] P. Álvarez-Caudevilla and J. López-Gómez, The dynamics of a class of cooperative systems, Discrete Contin. Dyn. Syst. 26 (2010), 397–415. 10.3934/dcds.2010.26.397Search in Google Scholar

[4] C. Bandle and M. Marcus, Large solutions of semilinear elliptic equations: Existence, uniqueness and asymptotic behavior, J. Anal. Math. 58 (1992), 9–24. 10.1007/BF02790355Search in Google Scholar

[5] S. Cano-Casanova and J. López-Gómez, Existence, uniqueness and blow-up rate of large solutions for a canonical class of one-dimensional problems on the half line, J. Differential Equations 244 (2008), 3180–3203. 10.1016/j.jde.2007.11.012Search in Google Scholar

[6] F. C. Cirstea and V. Radulescu, Existence and uniqueness of blow-up solutions for a class of logistic equations, Commun. Contemp. Math. 4 (2002), 559–586. 10.1142/S0219199702000737Search in Google Scholar

[7] F. C. Cirstea and V. Radulescu, Solutions with boundary blow-up for a class of nonlinear elliptic problems, Houston J. Math. 29 (2003), 821–829. Search in Google Scholar

[8] Y. Du, Order Structure and Topological Methods in Nonlinear Partial Differential Equations, World Scientific, Singapore, 2006. 10.1142/5999Search in Google Scholar

[9] Y. Du and Q. Huang, Blow-up solutions for a class of semilinear elliptic and parabolic equations, SIAM J. Math. Anal. 31 (1999), 1–18. 10.1137/S0036141099352844Search in Google Scholar

[10] S. Dumont, L. Dupaigne, O. Goubet and V. Radulescu, Back to the Keller–Osserman condition for boundary blow-up solutions, Adv. Nonlinear Stud. 7 (2007), 271–298. 10.1515/ans-2007-0205Search in Google Scholar

[11] C. Faber, Beweis daß unter allen homogenen Membranen von gleicher Fläche und gleicher Spannung die kreisförmige den tiefsten Grundton gibt, Sitzungsber. Bayer. Akad. der Wiss. Math. Phys. (1923), 169–171. Search in Google Scholar

[12] J. Garcia-Melián, Uniqueness of positive solutions for a boundary value problem, J. Math. Anal. Appl. 360 (2009), 530–536. 10.1016/j.jmaa.2009.06.077Search in Google Scholar

[13] J. Garcia-Melián, Multiplicity of positive solutions to boundary blow-up elliptic problems with sign-changing weights, J. Funct. Anal. 261 (2011), 1775–1798. 10.1016/j.jfa.2011.05.018Search in Google Scholar

[14] J. Garcia-Melián, R. Letelier and J. C. Sabina de Lis, Uniqueness and asymptotic behaviour for solutions of semilinear problems with boundary blow-up, Proc. Amer. Math. Soc. 129 (2001), 3593–3602. 10.1090/S0002-9939-01-06229-3Search in Google Scholar

[15] J. Garcia-Melián, J. Rossi and J. C. Sabina de Lis, Elliptic systems with boundary blow-up: Existence, uniqueness and applications to removability of singularities, Commun. Pure Appl. Anal. 15 (2016), no. 2, 549–562. 10.3934/cpaa.2016.15.549Search in Google Scholar

[16] J. Garcia-Melián and A. Suárez, Existence and uniqueness of positive large solutions to some cooperative elliptic systems, Adv. Nonlinear Stud. 3 (2003), no. 2, 193–206. 10.1515/ans-2003-0203Search in Google Scholar

[17] J. B. Keller, On solutions of Δu=f(u), Comm. Pure Appl. Math. 10 (1957), 503–510. 10.1002/cpa.3160100402Search in Google Scholar

[18] E. Krahn, Über eine von Rayleigh formulierte Minimaleigenschaft des Kreises, Math. Ann. 91 (1925), 97–100. 10.1007/BF01208645Search in Google Scholar

[19] A. C. Lazer and P. J. McKenna, On a problem of L. Bieberbach and H. Rademacher, Nonlinear Anal. 21 (1993), 327–335. 10.1016/0362-546X(93)90076-5Search in Google Scholar

[20] A. C. Lazer and P. J. McKenna, A singular elliptic boundary value problem, Appl. Math. Comput. 65 (1994), 183–194. 10.1016/0096-3003(94)90175-9Search in Google Scholar

[21] J. López-Gómez, The maximum principle and the existence of principal eigenvalues for some linear weighted boundary value problems, J. Differential Equations 127 (1996), 263–294. 10.1006/jdeq.1996.0070Search in Google Scholar

[22] J. López-Gómez, The boundary blow-up rate of large solutions, J. Differential Equations 195 (2003), 25–45. 10.1016/j.jde.2003.06.003Search in Google Scholar

[23] J. López-Gómez, Uniqueness of large solutions for a class of radially symmetric elliptic equations, Spectral Theory and Nonlinear Analysis with Applications to Spatial Ecology, World Scientific, Singapore (2005), 75–110. 10.1142/9789812701589_0005Search in Google Scholar

[24] J. López-Gómez, Optimal uniqueness theorems and exact blow-up rates of large solutions, J. Differential Equations 224 (2006), 385–439. 10.1016/j.jde.2005.08.008Search in Google Scholar

[25] J. López-Gómez, Uniqueness of radially symmetric large solutions, Discrete Contin. Dyn. Syst. 2007 (2007), 677–686. Search in Google Scholar

[26] J. López-Gómez, Second Order Linear Elliptic Operators, World Scientific, Singapore, 2013. 10.1142/8664Search in Google Scholar

[27] J. López-Gómez, Metasolutions of Parabolic Problems in Population Dynamics, CRC Press, Boca Raton, 2015. 10.1201/b19418Search in Google Scholar

[28] J. López-Gómez and L. Maire, Uniqueness of large positive solutions for a class of radially symmetric cooperative systems, J. Math. Anal. Appl. 435 (2016), 1738–1752. 10.1016/j.jmaa.2015.11.026Search in Google Scholar

[29] J. López-Gómez and L. Maire, Boundary blow-up rates and uniqueness of the large solution for cooperative elliptic systems of logistic type, Nonlinear Anal. Real World Appl. 33 (2017), 298–316. 10.1016/j.nonrwa.2016.07.001Search in Google Scholar

[30] J. López-Gómez and M. Molina-Meyer, The maximum principle for cooperative weakly coupled elliptic systems and some applications, Differential Integral Equations 7 (1994), 383–398. 10.57262/die/1369330435Search in Google Scholar

[31] J. López-Gómez, A. Tellini and F. Zanolin, High multiplicity and complexity of the bifurcation diagrams of large solutions for a class of superlinear indefinite problems, Commun. Pure Appl. Anal. 13 (2014), 1–73. 10.3934/cpaa.2014.13.1Search in Google Scholar

[32] M. Marcus and L. Véron, Uniqueness and asymptotic behavior of solutions with boundary blow-up for a class of nonlinear elliptic equations, Ann. Inst. H. Poincaré Anal. Non Linéaire 14 (1997), 237–274.10.1016/s0294-1449(97)80146-1Search in Google Scholar

[33] M. Marcus and L. Véron, Existence and uniqueness results for large solutions of general nonlinear elliptic equations, J. Evol. Equ. 3 (2004), 637–652. 10.1007/978-3-0348-7924-8_33Search in Google Scholar

[34] M. Molina-Meyer, Existence and uniqueness of coexistence states for some nonlinear elliptic systems, Nonlinear Anal. 25 (1995), 279–296. 10.1016/0362-546X(94)00137-7Search in Google Scholar

[35] M. Molina-Meyer, Global attractivity and singular perturbation for a class of nonlinear cooperative systems, J. Differential Equations 128 (1996), 347–378. 10.1006/jdeq.1996.0099Search in Google Scholar

[36] M. Molina-Meyer, Uniqueness and existence of positive solutions for weakly coupled general sublinear systems, Nonlinear Anal. 30 (1997), 5375–5380. 10.1016/S0362-546X(97)00380-5Search in Google Scholar

[37] J. D. Murray, Mathematical Biology, Springer, Berlin, 1989. 10.1007/978-3-662-08539-4Search in Google Scholar

[38] A. Okubo and S. A. Levin, Diffusion and Ecological Problems: Modern Perspectives, Interdiscip. Appl. Math. 14, Springer, New York, 2001. 10.1007/978-1-4757-4978-6Search in Google Scholar

[39] R. Osserman, On the inequality Δuf(u), Pacific J. Math. 7 (1957), 1641–1647. 10.2140/pjm.1957.7.1641Search in Google Scholar

[40] T. Ouyang and Z. Xie, The uniqueness of blow-up for radially symmetric semilinear elliptic equations, Nonlinear Anal. 64 (2006), 2129–2142. 10.1016/j.na.2005.08.010Search in Google Scholar

[41] T. Ouyang and Z. Xie, The exact boundary blow-up rate of large solutions for semilinear elliptic problems, Nonlinear Anal. 68 (2008), 2791–2800. 10.1016/j.na.2007.02.026Search in Google Scholar

[42] A. Ratto, M. Rigoli and L. Véron, Scalar curvature and conformal deformation of hyperbolic space, J. Funct. Anal. 121 (1994), 15–77. 10.1006/jfan.1994.1044Search in Google Scholar

[43] L. Véron, Semilinear elliptic equations with uniform blow up on the boundary, J. Anal. Math. 59 (1992), 231–250. 10.1007/BF02790229Search in Google Scholar

Received: 2016-10-19
Revised: 2017-01-04
Accepted: 2017-01-04
Published Online: 2017-01-26
Published in Print: 2017-07-01

© 2017 Walter de Gruyter GmbH, Berlin/Boston

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

Downloaded on 25.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/ans-2016-6018/html
Scroll to top button