Hostname: page-component-586b7cd67f-vdxz6 Total loading time: 0 Render date: 2024-11-28T04:32:00.818Z Has data issue: false hasContentIssue false

Asymptotic profiles for positive solutions of diffusive logistic equations

Published online by Cambridge University Press:  15 February 2023

Jian-Wen Sun
Affiliation:
School of Mathematics and Statistics, Gansu Key Laboratory of Applied Mathematics and Complex Systems, Lanzhou University, Lanzhou 730000, P.R. China ([email protected])
Peng-Fei Fang
Affiliation:
School of Mathematics and Statistics, Gansu Key Laboratory of Applied Mathematics and Complex Systems, Lanzhou University, Lanzhou 730000, P.R. China ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

In this paper, we study the asymptotic profiles of positive solutions for diffusive logistic equations. The aim is to study the sharp effect of linear growth and nonlinear function. Both the classical reaction-diffusion equation and nonlocal dispersal equation are investigated. Our main results reveal that the linear and nonlinear parts of reaction term play quite different roles in the study of positive solutions.

Type
Research Article
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction and main results

In this paper, we consider the diffusive logistic equation

(1.1)\begin{align} \begin{cases} \Delta u+\lambda u-a(x)u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega, \end{cases} \end{align}

where $\Omega$ is a $C^{2+\mu }$ bounded domain in $\mathbb {R}^N$ ($N\geq 2$), $\lambda >0$ is a real parameter, $p>1$ is constant, the boundary operator $B$ is given by

\begin{align*} Bu=\alpha u_\nu+\beta u, \end{align*}

here $\nu$ is the unit outward normal to $\partial \Omega$ and either $\alpha =0$, $\beta =1$ (the Dirichlet boundary condition) or $\alpha =1$, $\beta \geq 0$ (the Neumann or Robin boundary conditions). The function $a\in C^\mu (\bar {\Omega })$ and $a(x)>0$ for $x\in \bar {\Omega }$. Problem (1.1) is a basic reaction-diffusion model used in the study of diversity phenomena in the applied sciences (see, e.g. [Reference Amann1, Reference Berestycki3, Reference Brézis and Oswald4, Reference Li, López-Gómez and Sun15]). It is also the paradigmatic model in population dynamics, the diffusive logistic model [Reference Daners and López-Gómez7, Reference Du8, Reference Henry13, Reference López-Gómez16, Reference López-Gómez and Rabinowitz17]. The function $a(x)$ measures the capacity of $\Omega$ to support the species $u(x)$. Under the above assumptions, the semilinear problem (1.1) was well studied, see [Reference Li, López-Gómez and Sun15, Reference Murray18] and references therein.

In the case of $a(x)\equiv 0$, then (1.1) reduces to the following linear eigenvalue equation

(1.2)\begin{equation} \begin{cases} \Delta u+\lambda u=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega. \end{cases} \end{equation}

We know that (1.2) admits a unique positive principal eigenvalue $\lambda _1^B(\Omega )$ associated with a positive solution $\phi (x)$. Further, (1.1) admits a unique positive solution $u(x)$ if and only if $\lambda >\lambda _1^B(\Omega )$. However, we can see that (1.2) admits positive solutions if and only if $\lambda =\lambda _1^B(\Omega )$.

In the previous work [Reference Sun20], the sharp profiles of positive solutions to (1.1) for $\lambda >\lambda _1^B(\Omega )$ have been well investigated. In this paper, we shall consider the sharp changes of positive solutions between (1.1) and (1.2). To do this, we consider the following diffusive logistic problem

(1.3)\begin{equation} \begin{cases} \Delta u+(\lambda_1^B(\Omega)+\varepsilon^\alpha) u-a_\varepsilon(x)u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega, \end{cases} \end{equation}

where $\varepsilon >0$ is a parameter, $\alpha >0$ is a given constant, $a_\varepsilon \in C^\mu (\Omega )$ is positive in $\bar \Omega$ and there exist $\beta >0$ and $a\in C^\mu (\Omega )$ such that $a(x)>0$ for $x\in \bar \Omega$ and

(1.4)\begin{equation} \lim_{\varepsilon\to0+}\frac{a_\varepsilon(x)}{\varepsilon^\beta}=a(x) \text{ uniformly in }\bar\Omega. \end{equation}

In (1.4), the constant $\beta$ is the quenching rate of nonlinear function. It follows from the classical results of reaction-diffusion equation that (1.3) admits a unique positive solution $\theta _\varepsilon \in C^{2+\mu }(\Omega )$ for every $\varepsilon >0$, see e.g. [Reference Du9, Reference Li, López-Gómez and Sun15, Reference López-Gómez16]. According to (1.3), one may think that $\theta _\varepsilon (x)$ tends to the trivial solution or the positive eigenfunction of (1.2). However, our investigations reveal that $\theta _\varepsilon (x)$ admits quite different profiles, determined by various choices of $\alpha$ and $\beta$. In the present paper, we shall investigate the sharp profiles by the classical regularity estimates and uniform estimates of solutions [Reference Henry13, Reference López-Gómez16]. More precisely, we prove the following result.

Theorem 1.1 Let $\theta _\varepsilon \in C^{2+\mu }(\Omega )$ be the unique positive solution of (1.4) for $\varepsilon >0$ and $\Omega _*$ be a compact subset of $\Omega$.

  1. (i) If $\alpha <\beta,$ then

    (1.5)\begin{equation} \lim_{\varepsilon\to0+}\theta_\varepsilon(x)=\infty \text{ uniformly in } {\Omega_*}. \end{equation}
    Further, for any $x\in \Omega _*,$ there exist positive constants $c,\,C$ such that
    (1.6)\begin{equation} c\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\theta_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\theta_\varepsilon(x)\leq C. \end{equation}
  2. (ii) If $\alpha >\beta,$ then

    (1.7)\begin{equation} \lim_{\varepsilon\to0+}\theta_\varepsilon(x)=0 \text{ uniformly in } {\Omega_*}. \end{equation}
    Further, for any $x\in \Omega _*$, there exist positive constants $c_1,\,C_1$ such that
    (1.8)\begin{equation} c_1\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\theta_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\theta_\varepsilon(x)\leq C_1. \end{equation}
  3. (iii) If $\alpha =\beta,$ subject to a subsequence, we have

    \[ \lim_{\varepsilon\to0+}\theta_\varepsilon(x)=c_0\phi(x) \text{ uniformly in }\bar{\Omega} \]
    for some positive constant $c_0$.

Remark 1.2 It follows from theorem 1.1 that the linear term and nonlinear reaction function play quite different roles in the limiting behaviour of positive solutions of (1.1). We know from (1.5) and (1.6) that the blow-up phenomenon only occurs if the nonlinear function admits a quicker quenching speed, i.e. $\alpha <\beta$. It is interesting to point out that the blow-up phenomenon appears in the diffusive logistic equation with spatial degeneracy, see [Reference Du8, Reference Li, López-Gómez and Sun15]. However, if the linear term has a quicker quenching speed to the critical value $\lambda _1^B(\Omega )$, we get from (1.7) and (1.8) that the solution will tend to the trivial solution.

Since the diffusion may take place between nonadjoint places, the research in nonlocal dispersal equation has attracted much attention in recent years. Let $J: \mathbb {R}^N\to \mathbb {R}$ be a nonnegative and symmetric function. It is known that the nonlocal dispersal equation

(1.9)\begin{equation} u_t(x,t)=\int_{\mathbb{R}^N}J(x-y)[u(y,t)-u(x,t)]\,{\rm d}y \text{ in }\mathbb{R}^N\times(0,\infty), \end{equation}

and variations of it, arise in the study of different dispersal process in material science, ecology, neurology and genetics (see, for instance, [Reference Andreu-Vaillo, Mazón, Rossi and Toledo-Melero2, Reference Cortazar, Elgueta, Rossi and Wolanski5, Reference Kao, Lou and Shen12]). As stated in [Reference Fife10], if $u(y,\,t)$ is thought of as the density at location $y$ at time $t$, and $J(x-y)$ is thought of as the probability distribution of jumping from $y$ to $x$, then $\int _{\mathbb {R}^N}J(x-y)u(y,\,t)\,{\rm d}y$ denotes the rate at which individuals are arriving to location $x$ from all other places and $\int _{\mathbb {R}^N}J(y-x)u(x,\,t)\,{\rm d}y$ is the rate at which they are leaving location $x$ to all other places. Thus the right-hand side of (1.9) is the change of density $u(x,\,t)$. There has been attracted considerable interest in the study of nonlocal dispersal equations recently, for example, the papers [Reference Chasseigne, Chaves and Rossi6, Reference Garcia-Melian and Rossi11, Reference Hutson, Martinez, Mischaikow and Vickers14, Reference Sun19, Reference Sun21Reference Zhang, Li and Sun23] and references therein.

Let us consider the nonlocal dispersal logistic equation

(1.10)\begin{equation} \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+(\lambda_p(\Omega)+\varepsilon^\alpha) u-a_\varepsilon(x)u^p(x)=0 \text{ in } \bar\Omega, \end{equation}

where $\varepsilon >0$ is a parameter, $\alpha >0$ and $a_\varepsilon \in C(\bar \Omega )$ satisfies (1.4). In (1.10), the dispersal kernel function $J\in C(\mathbb {R}^N)$ is nonnegative, symmetric such that

\[ \int_{\mathbb{R}^N}J(y)\,{\rm d}y=1 \text{ and } J(0)>0, \]

and $\lambda _p(\Omega )$ stands for the unique principal eigenvalue of

\[ \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)={-}\lambda u(x)\text{ in } \bar\Omega. \]

In the rest of paper, we denoted by $\psi (x)$ the positive eigenfunction of $\lambda _p(\Omega )$. Then for any $\varepsilon >0$, we know that (1.10) admits a unique positive solution $\omega _\varepsilon (x)$, see [Reference Garcia-Melian and Rossi11, Reference Sun, Li and Wang22].

Since the nonlocal dispersal equation shares many properties with the reaction-diffusion equation, it is interesting to investigate the sharp behaviour of positive solutions of (1.10) as $\varepsilon \to 0$. However, there is a deficiency of regularity theory and compact property for nonlocal dispersal operators, the study of sharp behaviour of (1.10) is quite different to (1.3), [Reference Amann1, Reference Chasseigne, Chaves and Rossi6, Reference Henry13]. We shall obtain the asymptotic behaviour for nonlocal dispersal problem (1.10) by the means of nonlocal estimates and comparison arguments.

In the case of nonlocal dispersal logistic equation, we have the following result.

Theorem 1.3 Let $\omega _\varepsilon \in C(\bar \Omega )$ be the unique positive solution of (1.10) for $\varepsilon >0$.

  1. (i) If $\alpha <\beta,$ then

    (1.11)\begin{equation} \lim_{\varepsilon\to0+}\omega_\varepsilon(x)=\infty \text{ uniformly in }\bar {\Omega}. \end{equation}
    Further, there exist positive constants $c,\,C$ such that
    \[ c\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\omega_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\omega_\varepsilon(x)\leq C \]
    for any $x\in \bar \Omega$.
  2. (ii) If $\alpha >\beta,$ then

    (1.12)\begin{equation} \lim_{\varepsilon\to0+}\omega_\varepsilon(x)=0 \text{ uniformly in } \bar{\Omega}. \end{equation}
    Further, there exist positive constants $c_1,\,C_1$ such that
    \[ c_1\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\omega_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}\omega_\varepsilon(x)\leq C_1 \]
    for any $x\in \bar \Omega$.
  3. (iii) If $\alpha =\beta,$ then, subject to a subsequence, we have

    \[ \lim_{\varepsilon\to0+}\omega_\varepsilon(x)=c_0\psi(x) \text{ uniformly in }\bar{\Omega} \]
    for some positive constant $c_0$.

The conclusions in theorem 1.3 provide us how the sharp profiles of positive solutions to (1.10) is determined by $\alpha$ and $\beta$. We also know that the profile for nonlocal problem is different to the classical reaction-diffusion equation. By (1.11), we obtain that the positive solution for nonlocal problem (1.10) will blow-up in the whole domain $\Omega$ when $\alpha <\beta$. Similarly, by (1.12), we know that quenching occurs for all $x\in \bar \Omega$.

The rest of this paper is organized as follows. In § 2, we investigate the profiles of reaction-diffusion equation (1.3). Section 3 is devoted to the sharp profiles of nonlocal dispersal logistic equations.

2. Profiles for reaction-diffusion equations

In this section, we investigate the limiting behaviour of positive solutions for the diffusive logistic equation (1.3). It follows from the classical results [Reference Brézis and Oswald4, Reference Henry13] that there exists a unique positive solution $\theta _\varepsilon \in C^{2+\mu }(\Omega )$ to (1.3) for every $\varepsilon >0$. Moreover, the positive solution $\theta _\varepsilon$ is continuous with respect to $\varepsilon$. In what follows, we always assume that $a_\varepsilon,\,a\in C^\mu (\Omega )$ are positive in $\bar \Omega$ and

\[ \lim_{\varepsilon\to0+}\frac{a_\varepsilon(x)}{\varepsilon^\alpha}=a(x) \text{ uniformly in }\bar\Omega. \]

We first study the following diffusive logistic equation

(2.1)\begin{equation} \begin{cases} \Delta u+\lambda u-a(x)u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega. \end{cases} \end{equation}

We can see that (2.1) admits a unique positive solution $\theta _\lambda (x)$ if and only if $\lambda >\lambda _1^B(\Omega )$. Moreover, $\theta _\lambda (x)$ is continuous with respect to $\lambda$ and

\[ \lim_{\lambda\to\lambda_1^B(\Omega)+}\theta_\lambda(x)=0 \text{ locally uniformly in }\Omega. \]

We shall give the decay estimates of $\theta _\lambda (x)$ near $\lambda _1^B(\Omega )$ as follows.

Lemma 2.1 Suppose that $\Omega _*$ is a subdomain of $\Omega$ such that $\bar \Omega _* \subset \Omega$. Let $\theta _\lambda (x)$ be the unique positive solution of (2.1) for $\lambda \in (\lambda _1^B(\Omega ),\,\lambda _1^B(\Omega )+1],$ then there exist positive constants $c$ and $C,$ independent of $\lambda$ such that

(2.2)\begin{equation} c\left[\frac{\lambda-\lambda_1^B(\Omega)}{\max_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}\leq\theta_\lambda(x)\leq C\left[\frac{\lambda-\lambda_1^B(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}} \end{equation}

for $x\in \bar \Omega _*$.

Proof. By the uniqueness of positive solution to (2.1), we can find positive constant $M$, independent of $\lambda$ such that

(2.3)\begin{equation} 0<\max_{\bar\Omega}\theta_\lambda(x)\leq M-1. \end{equation}

Let $\phi (x)$ be a positive eigenfunction of $\lambda _1^B(\Omega )$ such that $\|\phi \|_{L^\infty (\Omega )}=1$. Denote

\[ \Omega^*=\left\{x\in\bar\Omega: dist(x,\Omega_*)>\inf_{x\in\partial\Omega,y\in\partial\Omega_*}\frac{|x-y|}{2}\right\}, \]

and take $C_1>0$ such that

(2.4)\begin{equation} C_1\phi(x)>\left[\frac{\lambda-\lambda_1^B(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}} \end{equation}

for $x\in \Omega _*$. Using (2.3) and (2.4), we know that there exists smooth function $u(x)$ such that

\[ u(x)= \begin{cases} C_1\phi(x) & \text{ if }x\in\bar\Omega_*,\\ M & \text{ if }x\in\bar\Omega^*,\\ \end{cases} \]

and $u(x)$ is an upper-solution to (2.1). Since $\phi (x)$ is independent to $\lambda$, we know from the comparison principle that the right-hand side of (2.2) holds.

On the other hand, we define $v(x)=c_1\phi (x)$, where

\[ c_1=\left[\frac{\lambda-\lambda_1^B(\Omega)}{\max_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}. \]

It is easy to see that $v(x)$ is a lower-solution to (2.1) and we obtain the left-hand side of (2.2) by the uniqueness of positive solutions to (2.1). The proof is completed.

Lemma 2.2 Let $u_\varepsilon \in C^{2+\mu }(\Omega )$ be the unique positive solution of

(2.5)\begin{equation} \begin{cases} \Delta u+(\lambda_1^B(\Omega)+\varepsilon^\alpha)u-\varepsilon^\beta a(x)u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega \end{cases} \end{equation}

for $\varepsilon >0$.

  1. (i) If $\alpha <\beta$, then

    \[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=\infty \text{ locally uniformly in } {\Omega}. \]
    Further, there exist positive constants $c,\,C$ such that
    \[ c\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq C. \]
  2. (ii) If $\alpha >\beta,$ then

    \[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=0 \text{ locally uniformly in } {\Omega}. \]
    Further, there exist positive constants $c_1,\,C_1$ such that
    \[ c_1\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq C_1. \]
  3. (iii) If $\alpha =\beta,$ then, subject to a subsequence, we have

    \[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=c_0\phi(x) \text{ uniformly in }\bar{\Omega} \]
    for some positive constant $c_0$.

Proof. Set $v_\varepsilon (x)=\varepsilon ^{\frac {\beta }{p-1}}u_\varepsilon (x)$, it becomes apparent that $v_\varepsilon (x)$ is the unique positive solution of

\[ \begin{cases} \Delta u+(\lambda_1^B(\Omega)+\varepsilon^\alpha)u- a(x)u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega. \end{cases} \]

Let $\Omega _*$ be a compact subset of $\Omega$, thanks to lemma 2.1, we know that there exist $c_0,\,C_0$ such that

\[ c_0\varepsilon^{\frac{\alpha}{p-1}}\leq v_\varepsilon(x)\leq C_0\varepsilon^{\frac{\alpha}{p-1}} \]

for $x\in \bar \Omega _*$. Hence we obtain

(2.6)\begin{equation} c_0\varepsilon^{\frac{\alpha-\beta}{p-1}}\leq u_\varepsilon(x)\leq C_0\varepsilon^{\frac{\alpha-\beta}{p-1}} \end{equation}

for $x\in \bar \Omega _*$. According to (2.6), we obtain the conclusions (i) and (ii).

By standard interior estimates and (2.6), there exists a positive constant $\tilde C=\tilde C(\Omega _*)$ such that

\[ \|u_\varepsilon\|_{C^{2+\mu}(\bar\Omega_*)}\leq \tilde C. \]

Therefore, by passing to a subsequence and the diagonal argument, there exists $u\in L^2(\Omega )$ such that

\[ \lim_{\varepsilon\to0+}u_{\varepsilon}(x)=u(x) \text{ weakly in }W^{1,2}(\Omega)\;\; \text{ and strongly in } L^2(\Omega). \]

Thanks to (2.5), we know that $u(x)$ is a positive weak solution of

(2.7)\begin{equation} \begin{cases} \Delta u+\lambda_1^B(\Omega)u=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega. \end{cases} \end{equation}

By elliptic regularity, it must be a strong solution. By the uniqueness of the positive solution of (2.7), $u(x)=c_0\phi (x)$ for some positive constant $c_0$. As this argument is independent of the sequence $\varepsilon$, it is apparent from Sobolev imbedding theorem that

\[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=c_0\phi(x) \quad \hbox{ uniformly in }\bar\Omega. \]

Thus the proof is completed.

At the end of this section, we prove the main result theorem 1.1.

Proof Proof of theorem 1.1

We first take $\delta >0$ such that

\[ a(x)>\delta>0 \]

for $x\in \bar \Omega$. Then we choose $\varepsilon >0$ small, denoted by $\varepsilon <\varepsilon _0$ such that

\[ a(x)+1\geq\frac{a_\varepsilon(x)}{\varepsilon^\alpha}\geq a(x)-\delta>0 \]

for $x\in \bar \Omega$.

Now let $\hat {u}(x)$ be the unique positive solution of

\[ \begin{cases} \Delta u+(\lambda_1^B(\Omega)+\varepsilon^\alpha)u-\varepsilon^\beta[a(x)-\delta]u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega, \end{cases} \]

and $\bar {u}(x)$ be the unique positive solution of

\[ \begin{cases} \Delta u+(\lambda_1^B(\Omega)+\varepsilon^\alpha)u-\varepsilon^\beta[a(x)+1]u^p=0 & \text{ in } \Omega,\\ Bu=0 & \text{ on } \partial\Omega \end{cases} \]

for $\varepsilon >0$, respectively. A simple argument from upper–lower solutions gives

(2.8)\begin{equation} 0<\bar{u}(x)\leq\theta_\varepsilon(x)\leq \hat{u}(x) \end{equation}

for $x\in \Omega$.

Thus we know from (2.8) and lemma 2.2 that the conclusions (i)–(iii) of theorem 1.1 are true.

3. Profiles for nonlocal dispersal logistic equation

In this section, we investigate the limiting behaviour of positive solutions of (1.10) as $\varepsilon \to 0+$. It follows from [Reference Garcia-Melian and Rossi11, Reference Sun21] that there exists a unique positive solution $\omega _\varepsilon \in C(\bar \Omega )$ to (1.10) for every $\varepsilon >0$ and $\theta _\varepsilon$ is continuous with respect to $\varepsilon$. In the rest of this section, for simplicity, we always assume that $a_\varepsilon,\,a\in C(\bar \Omega )$ are positive in $\bar \Omega$ and

\[ \lim_{\varepsilon\to0+}\frac{a_\varepsilon(x)}{\varepsilon^\alpha}=a(x) \text{ uniformly in }\bar\Omega. \]

We first give some estimates for the positive solution of

(3.1)\begin{equation} \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+\lambda u-a(x)u^p(x)=0 \text{ in } \bar\Omega. \end{equation}

The positive solution problem (3.1) has been well investigated, see e.g. [Reference Garcia-Melian and Rossi11, Reference Sun20Reference Sun, Li and Wang22].

Lemma 3.1 Let $\omega _\lambda (x)$ be the unique positive solution of (1.10) for $\lambda \in (\lambda _p(\Omega ),\,\lambda _p(\Omega )+1],$ then there exist positive constants $c$ and $C,$ independent of $\lambda$ such that

\[ c\left[\frac{\lambda-\lambda_p(\Omega)}{\max_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}\leq\omega_\lambda(x)\leq C\left[\frac{\lambda-\lambda_p(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}} \]

for $x\in \bar \Omega$.

Proof. By the uniqueness of positive solution to (3.1), we can find positive constant $M$, independent of $\lambda$ such that

\[ 0<\max_{\bar\Omega}\theta_\lambda(x)\leq M. \]

Let $\psi (x)$ be a positive eigenfunction of $\lambda _p(\Omega )$ such that $\|\psi \|_{L^\infty (\Omega )}=1$. Since $\psi (x)>0$ for $x\in \bar \Omega$, we can take $C_1>0$ such that

\[ C_1\phi(x)\geq\left[\frac{\lambda-\lambda_p(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}} \]

for $x\in \bar \Omega$. Then a direct computation gives that $C_1\phi (x)$ is an upper-solution to (3.1) and we know from the uniqueness of positive solution that

\[ \omega_\lambda(x)\leq C_1\phi(x) \]

for $x\in \bar \Omega$. Hence we obtain

\[ \omega_\lambda(x)\leq C_1\phi(x) \leq C\left[\frac{\lambda-\lambda_p(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}, \]

by taking $C=[\min _{\bar \Omega }\phi (x)]^{-1}$ and

\[ C_1=[\min_{\bar\Omega}\phi(x)] \left[\frac{\lambda-\lambda_p(\Omega)}{\min_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}. \]

On the other hand, we define

\[ v(x)=\left[\frac{\lambda-\lambda_p(\Omega)}{\max_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}}\psi(x). \]

It is easy to see that $v(x)$ is a lower-solution to (3.1). But $\psi (x)$ is independent to $\lambda$, it follows from the comparison principle that there exists $c>0$ such that

\[ \omega_\lambda(x)\geq c\left[\frac{\lambda-\lambda_p(\Omega)}{\max_{\bar\Omega}a(x)}\right]^{\frac{1}{p-1}} \]

for $x\in \bar \Omega$.

Lemma 3.2 Let $u_\varepsilon \in C(\bar \Omega )$ be the unique positive solution of

(3.2)\begin{equation} \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+(\lambda_p(\Omega)+\varepsilon^\alpha) u-\varepsilon^\beta a(x)u^p(x)=0 \text{ in } \bar\Omega \end{equation}

for $\varepsilon >0$.

  1. (i) If $\alpha <\beta,$ then

    \[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=\infty \text{ uniformly in } \bar{\Omega}. \]
    Further, there exist positive constants $c,\,C$ such that
    \[ c\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq C. \]
  2. (ii) If $\alpha >\beta,$ then

    \[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=0 \text{ uniformly in } \bar{\Omega}. \]
    Further, there exist positive constants $c_1,\,C_1$ such that
    \[ c_1\leq\liminf_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq \limsup_{\varepsilon\to0+}\varepsilon^{\frac{\beta-\alpha}{p-1}}u_\varepsilon(x)\leq C_1. \]
  3. (iii) If $\alpha =\beta,$ then, subject to a subsequence, we have

    (3.3)\begin{equation} \lim_{\varepsilon\to0+}u_\varepsilon(x)=c_0\psi(x) \text{ uniformly in }\bar{\Omega} \end{equation}
    for some positive constant $c_0$.

Proof. Set $v_\varepsilon (x)=\varepsilon ^{\frac {\beta }{p-1}}u_\varepsilon (x)$, then we can see that $v_\varepsilon (x)$ is the unique positive solution of

\[ \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+(\lambda_p(\Omega)+\varepsilon^\alpha) u- a(x)u^p(x)=0 \text{ in } \bar\Omega. \]

Then we know from lemma 3.1 that there exist $c_0,\,C_0$ such that

\[ c_0\varepsilon^{\frac{\alpha}{p-1}}\leq v_\varepsilon(x)\leq C_0\varepsilon^{\frac{\alpha}{p-1}} \]

for $x\in \bar \Omega$. Hence we obtain

\[ c_0\varepsilon^{\frac{\alpha-\beta}{p-1}}\leq u_\varepsilon(x)\leq C_0\varepsilon^{\frac{\alpha-\beta}{p-1}} \]

for $x\in \bar \Omega$ and the conclusions (i) and (ii) are followed.

At last, we prove (3.3). In this case, we still have

\[ c_0\leq u_\varepsilon(x)\leq C_0 \]

for $x\in \bar \Omega$. Since $\lambda _p(\Omega )\in (0,\,1)$ and

\[ \left[1 -\lambda_p(\Omega)-\varepsilon^\alpha+{\varepsilon^\beta}a(x)(u_\varepsilon(x))^{p-1}\right]u_\varepsilon(x)=\int_{\Omega}J(x-y)u_\varepsilon(y)\,{\rm d}y \text{ in } \bar\Omega, \]

we know that there exists $\rho >0$ which is independent to $\varepsilon$ such that

(3.4)\begin{equation} 1 -\lambda_p(\Omega)-\varepsilon^\alpha+{\varepsilon^\beta}a(x)(u_\varepsilon(x))^{p-1}\geq \rho \end{equation}

for $x\in \bar \Omega$, provided $\varepsilon \in (0,\,1)$ is small. Then for any $x_1,\,x_2\in \bar \Omega$, without loss of generality, we may assume that $u_\varepsilon (x_1)>u_\varepsilon (x_2)$. A direct computation from (3.2)(3.4) shows that

\begin{align*} & (1 -\lambda_p(\Omega)-\varepsilon^\alpha+p\varepsilon^\beta a(x_2)\theta_\varepsilon^{p-1})[u_\varepsilon(x_1)-u_\varepsilon(x_2)]\\ & \quad=\int_{\Omega}(J(x_1,y)-J(x_2,y))u_\varepsilon(y)\,{\rm d}y+\varepsilon^\beta(a(x_2)-a(x_1))u_\varepsilon^p(x_1)\\ & \quad\leq C_0\int_{\Omega}|J(x_1,y)-J(x_2,y)|\,{\rm d}y+C_0^p|(a(x_2)-a(x_1))|, \end{align*}

here $\theta _\varepsilon$ is between $u_\varepsilon (x_2)$ and $u_\varepsilon (x_1)$. Thus we obtain that

(3.5)\begin{equation} |u_\varepsilon(x_1)-u_\varepsilon(x_2)|\leq \frac{C_0\int_{\Omega}|J(x_1,y)-J(x_2,y)|\,{\rm d}y+C_0^p|(a(x_2)-a(x_1))|}{\rho} \end{equation}

for $x_1,\,x_2\in \bar \Omega$. It follows from (3.5) and a compact argument that we can extract a subsequence still denoted by $\varepsilon$ and there exists positive function $V\in C(\bar \Omega )$ such that

\[ \lim_{\varepsilon\to0+}u_\varepsilon(x)=V(x) \quad \hbox{ uniformly in }\bar\Omega, \]

and

(3.6)\begin{equation} \int_{\Omega}J(x-y)V(y)\,{\rm d}y-V(x)+\lambda_p(\Omega) V(x)=0 \text{ in } \bar\Omega. \end{equation}

Note that $\lambda _p(\Omega )$ is the unique principal eigenvalue of (3.6), we know that (3.3) holds.

We are ready to prove the main result theorem 1.3.

Proof Proof of theorem 1.3

We first take $\delta >0$ such that

\[ a(x)>\delta>0 \]

for $x\in \bar \Omega$. Then we can choose $\varepsilon >0$ small, denoted by $\varepsilon <\varepsilon _0$ such that

\[ a(x)+1\geq\frac{a_\varepsilon(x)}{\varepsilon^\alpha} \geq a(x)-\delta>0 \]

for $x\in \bar \Omega$. Let $\hat {u}(x)$ be the unique positive solution of

\[ \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+(\lambda_p(\Omega)+\varepsilon^\alpha)u-\varepsilon^\beta[a(x)-\delta]u^p=0 \text{ in } \bar\Omega, \]

and $\bar {u}(x)$ be the unique positive solution of

\[ \int_{\Omega}J(x-y)u(y)\,{\rm d}y-u(x)+(\lambda_p(\Omega)+\varepsilon^\alpha)u-\varepsilon^\beta[a(x)+1]u^p=0 \text{ in } \bar\Omega \]

for $\varepsilon >0$, respectively. Thus we get from the comparison principle that

\[ 0<\bar{u}(x)\leq\omega_\varepsilon(x)\leq \hat{u}(x) \]

for $x\in \bar \Omega$.

The conclusions (i)–(iii) of theorem 1.3 are followed by lemma 3.2.

Acknowledgements

The author would like to thank the anonymous reviewer for his/her helpful comments. This work was partially supported by NSF of China (11731005), FRFCU (lzujbky-2021-52) and NSF of Gansu (21JR7RA535, 21JR7RA537).

References

Amann, H.. Existence of multiple solutions for nonlinear elliptic boundary value problems. Indiana Univ. Math. J. 21 (1971-72), 925935.CrossRefGoogle Scholar
Andreu-Vaillo, F., Mazón, J. M., Rossi, J. D. and Toledo-Melero, J.. Nonlocal diffusion problems, mathematical surveys and monographs (AMS, Providence, Rhode Island, 2010).CrossRefGoogle Scholar
Berestycki, H.. Le nombre de solutions de certains problèmes semi-linéaires elliptiques. J. Funct. Anal. 40 (1981), 129.CrossRefGoogle Scholar
Brézis, H. and Oswald, L.. Remarks on sublinear elliptic equations. Nonlinear Anal. 10 (1986), 5564.CrossRefGoogle Scholar
Cortazar, C., Elgueta, M., Rossi, J. D. and Wolanski, N.. Boundary fluxes for nonlocal diffusion. J. Differ. Equ. 234 (2007), 360390.CrossRefGoogle Scholar
Chasseigne, E., Chaves, M. and Rossi, J. D.. Asymptotic behavior for nonlocal diffusion equations. J. Math. Pures Appl. 86 (2006), 271291.CrossRefGoogle Scholar
Daners, D. and López-Gómez, J.. Global dynamics of generalized logistic equations. Adv. Nonl. Studies 18 (2018), 217236.CrossRefGoogle Scholar
Du, Y.. Spatial patterns for population models in a heterogeneous environment. Taiwanese J. Math. 8 (2004), 155182.CrossRefGoogle Scholar
Du, Y.. Order structure and topological methods in nonlinear partial differential equations. Maximum principle and applications, vol 1, Fang-Hua Lin (Courant Institute of Math. Sci., New York University) (Series ed.) (Singapore: World Scientific Publishing, 2006).CrossRefGoogle Scholar
Fife, P.. Some nonclassical trends in parabolic and parabolic-like evolutions, In: Trends in Nonlinear Analysis (Springer, Berlin, 2003), pp. 153–191.CrossRefGoogle Scholar
Garcia-Melian, J. and Rossi, J. D.. A logistic equation with refuge and nonlocal diffusion. Commun. Pure Appl. Anal. 8 (2009), 20372053.CrossRefGoogle Scholar
Kao, C. Y., Lou, Y. and Shen, W. X.. Random dispersal vs. non-local dispersal. Discrete Contin. Dyn. Syst. 26 (2010), 551596.CrossRefGoogle Scholar
Henry, D.. Geometric theory of semilinear parabolic equations, Lecture Notes Math., Vol. 840 (Springer-Verlag, 1981).CrossRefGoogle Scholar
Hutson, V., Martinez, S., Mischaikow, K. and Vickers, G. T.. The evolution of dispersal. J. Math. Biol. 47 (2003), 483517.CrossRefGoogle ScholarPubMed
Li, W. T., López-Gómez, J. and Sun, J. W.. Sharp blow-up profiles of positive solutions for a class of semilinear elliptic problems. Adv. Nonlinear Stud. 21 (2021), 751765.CrossRefGoogle Scholar
López-Gómez, J.. Metasolutions of Parabolic Equations in Population Dynamics (Boca Raton: CRC Press, 2016).Google Scholar
López-Gómez, J. and Rabinowitz, P.. The effects of spatial heterogeneities on some multiplicity results. Discrete Contin. Dyn. Syst. 127 (2016), 941952.Google Scholar
Murray, J.. Mathematical Biology, 2nd ed. (New York: Springer-verlag, 1998).Google Scholar
Sun, J. W.. Limiting solutions of nonlocal dispersal problem in inhomogeneous media. J. Dynam. Differ. Equ. 34 (2022), 14891504.CrossRefGoogle Scholar
Sun, J. W.. Asymptotic profiles in diffusive logistic equations. Z. Angew. Math. Phys. 72 (2021), 152.CrossRefGoogle Scholar
Sun, J. W.. Effects of dispersal and spatial heterogeneity on nonlocal logistic equations. Nonlinearity 34 (2021), 54345455.CrossRefGoogle Scholar
Sun, J. W., Li, W. T. and Wang, Z. C.. A nonlocal dispersal logistic equation with spatial degeneracy. Discrete Contin. Dyn. Syst. 35 (2015), 32173238.CrossRefGoogle Scholar
Zhang, G. B., Li, T. T. and Sun, Y. J.. Asymptotic behavior for nonlocal dispersal equations. Nonlinear Anal. 72 (2010), 44664474.CrossRefGoogle Scholar