Theory and computation of radial solutions for Neumann problems with φ-Laplacian

1 year ago

(original), ISI/JCR, Numerical Analysis, paper, PDEs (numerical)

Abstract

The paper deals with existence, localization and multiplicity of radial positive solutions in the annulus or the ball, for the Neumann problem involving a general φ $ϕ$ -Laplace operator. Our results apply in particular to the classical Laplacian and to the mean curvature operators in the Euclidean and Minkowski spaces. Numerical experiments with the MATLAB object-oriented package Chebfun are performed to obtain numerical solutions for some concrete equations.

Authors

Radu Precup
Institute of Advanced Studies in Science and Technology, Babeş-Bolyai, Cluj-Napoca, Romania
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, Cluj-Napoca, Romania

Călin-Ioan Gheorghiu
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, Cluj-Napoca, Romania

Keywords

Neumann boundary value problem; $ϕ$ -Laplace operator; Radial solution; Positive solution; Fixed point index; Harnack type inequality; Numerical solution

Paper coordinates

R. Precup, C.-I. Gheorghiu, Theory and computation of radial solutions for Neumann problems with $ϕ$ -Laplacian, Qualitative Theory of Dynamical Systems, 23 (2024), art. no. 107, https://doi.org/10.1007/s12346-024-00963-8

PDF

About this paper

Journal

Qualitative Theory of Dynamical Systems

Publisher Name

Springer International Publishing

DOI

https://doi.org/10.1007/s12346-024-00963-8

Print ISSN

1662-3592

Online ISSN

1575-5460

google scholar link

[1] Bartnik, R., Simon, L.: Spacelike hypersurfaces with prescribed boundary values and mean curvature. Commun. Math. Phys. 87, 131–152 (1982). https://doi.org/10.1007/BF01211061, Article ADS MathSciNet Google Scholar
[2] Benedikt, J., Girg, P., Kotrla, L., Takac, P.: Origin of the φ-Laplacian and A. Missbach. Electron. J. Differ. Equ. 2018(16), 1–17 (2018). http://ejde.math.txstate.edu or http://ejde.math.unt.edu
[3] Bereanu, C., Jebelean, P., Mawhin, J.: Radial solutions for Neumann problems involving mean curvature operators in Euclidean and Minkowski spaces. Math. Nachr. 283, 379–391 (2010). https://doi.org/10.1002/mana.200910083, Article MathSciNet Google Scholar
[4] Bonheure, D., Noris, B., Weth, T.: Increasing radial solutions for Neumann problems without growth restrictions. Ann. Inst. H. Poincaré Anal. Non Linéaire 29, 573–588 (2012). https://doi.org/10.1016/j.anihpc.2012.02.002, Article ADS MathSciNet Google Scholar
[5] Bonheure, D., Serra, E., Tilli, P.: Radial positive solutions of elliptic systems with Neumann boundary conditions. J. Funct. Anal. 265, 375–398 (2013). https://doi.org/10.1016/j.jfa.2013.05.027, Article MathSciNet Google Scholar
[6] Boscaggin, A., Colasuonno, F., Noris, B.: Multiple positive solutions for a class of $p$ φ-Laplacian Neumann problems without growth conditions. ESAIM Control Optim. Calc. Var. 24, 1625–1644 (2018). https://doi.org/10.1051/cocv/2017074, Article MathSciNet Google Scholar
[7] Boscaggin, A., Feltrin, G.: Pairs of positive radial solutions for a Minkowski-curvature Neumann problem with indefinite weight. Nonlinear Anal. 196, 111807 (2020). https://doi.org/10.1016/j.na.2020.111807, Article MathSciNet Google Scholar
[8] Colasuonno, F., Noris, B.: Radial positive solutions for Laplacian supercritical Neumann problems. Bruno Pini Math. Anal. Semin. 8 (1), 55–72 (2017).
https://doi.org/10.6092/issn.2240-2829/7797,

[9] Deimling, K.: Nonlinear Functional Analysis. Springer, Berlin (1985) Book Google Scholar
[10] Gheorghiu, C.I.: A third-order nonlinear BVP on the half-line, Chebfun (2020). https://www.chebfun.org/examples/ode-nonlin/GulfStream.html
[11] López-Gómez, J., Omari, P., Rivetti, S.: Positive solutions of a one-dimensional indefinite capillarity-type problem: a variational approach. J. Differ. Equ. 262, 2335–2392 (2017) Article ADS MathSciNet Google Scholar
[12] Ma, R., Chen, T., Wang, H.: Nonconstant radial positive solutions of elliptic systems with Neumann boundary conditions. J. Math. Anal. Appl. 443, 542–565 (2016). https://doi.org/10.1016/j.jmaa.2016.05.038, Article MathSciNet Google Scholar
[13] O’Regan, D., Wang, H.: Positive radial solutions for $p$ φ-Laplacian systems. Aequ. Math. 75, 43–50 (2008). https://doi.org/10.1007/s00010-007-2909-3 Article MathSciNet Google Scholar
[14] Precup, R.: Existence, localization and multiplicity results for positive radial solutions of semilinear elliptic systems. J. Math. Anal. Appl. 352, 48–56 (2009). https://doi.org/10.1016/j.jmaa.2008.01.097, Article MathSciNet Google Scholar
[15] Precup, R., Pucci, P., Varga, C.: Energy-based localization and multiplicity of radially symmetric states for the stationary $p$ φ-Laplace diffusion. Complex Var. Ellipt. Equ. 65, 1198–1209 (2020). https://doi.org/10.1080/17476933.2019.1574774 Article MathSciNet Google Scholar
[16] Precup, R., Rodriguez-Lopez, J.: Positive radial solutions for Dirichlet problems via a Harnack type inequality Math. Meth. Appl. Sci. 46(2), 2972–2985 (2023). https://doi.org/10.1002/mma.8682
[17] Trefethen, L.N.: Computing numerically with functions instead of numbers. Math. Comput. Sci. 1, 9–19 (2007). https://doi.org/10.1007/s11786-007-0001-y Article MathSciNet Google Scholar
[18] Trefethen, L.N., Birkisson, A., Driscoll, T.A.: Exploring ODEs. SIAM, Philadelphia (2018) Google Scholar
[19] Wang, H.: On the existence of positive solutions for semilinear elliptic equations in the annulus. J. Differ. Equ. 109, 1–7 (1994). https://doi.org/10.1006/jdeq.1994.1042

Paper (preprint) in HTML form

Theory and computation of radial solutions for Neumann problems with $\phi$ -Laplacian

Radu Precup, Călin-Ioan Gheorghiu R. Precup, Faculty of Mathematics and Computer Science and Institute of Advanced Studies in Science and Technology, Babeş-Bolyai University, 400084 Cluj-Napoca, Romania & Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, P.O. Box 68-1, 400110 Cluj-Napoca, Romania C.-I. Gheorghiu, Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, P.O. Box 68-1, 400110 Cluj-Napoca, Romania

Abstract.

The paper deals with existence, localization and multiplicity of radial positive solutions in the annulus or the ball, for the Neumann problem involving a general $\phi$ -Laplace operator. Our results apply in particular to the mean curvature operator in the Minkowski space. Numerical experiments with the MATLAB object-oriented package Chebfun are performed to obtain numerical solutions for some concrete equations.

Key words and phrases:

Neumann boundary value problem,

\phi

-Laplace operator, radial solution, positive solution, fixed point index, Harnack type inequality, numerical solution

1991 Mathematics Subject Classification:

35J65, 34B18, 65L10

1. Introduction

In this paper we are concerning with the existence and localization of radial positive solutions for the Neumann boundary value problem in the annulus or in the ball

(1.1)

\left\{\begin{array}[]{ll}-\text{div\ }\left(\psi\left(\left|\nabla u\right|% \right)\nabla u\right)+\varepsilon u=f\left(\left|x\right|,u\right)&\text{in }% \Omega\\ \partial_{\nu}u=0&\text{on }\partial\Omega,\end{array}\right.

where $\varepsilon>0,$ $\ \psi:\left(-a,a\right)\rightarrow\mathbb{R}_{+}\$ is such that the function

\phi:\left(-a,a\right)\rightarrow\left(-b,b\right),\ \ \ \ \phi\left(s\right)=% s\psi\left(s\right)\ \ \ \ \ \left(0<a,b\leq+\infty\right)

is an increasing homeomorphism, $f:\left[R_{0},R\right]\times\mathbb{R}_{+}\rightarrow\mathbb{R}_{+}$ is continuous and $\ \Omega=\{x\in\mathbb{R}^{n}:\ R_{0}<\left|x\right|<R\},\ n\geq 2.$ Here $\ 0\leq R_{0}<R<+\infty$ and $\nu$ is the the exterior unit normal vector to the boundary of $\Omega.$

Problems of type (1.1) arise from mathematical modeling of real processes. Thus, equations involving the $p$ -Laplacian come from fluid mechanics in porous media [2], equations with a singular homeomorphism arise from the relativistic mechanics [1], and equations involving a bounded homeomorphism intervene in capillarity problems [11].

Looking for radial solutions of (1.1), that is, functions of the form $\ u(x)=v(r)$ with $\ r=\left|x\right|,$ problem (1.1) reduces to the boundary value problem

(1.2)

\left\{\begin{array}[]{ll}-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{% \prime}+\varepsilon r^{n-1}v=r^{n-1}f(r,v)&\text{in }\left(R_{0},R\right)\\ v^{\prime}\left(R_{0}\right)=v^{\prime}(R)=0.&\end{array}\right.

Note that in the case of the ball, when $R_{0}=0,$ the equality $v^{\prime}(R)=0$ stands for the Neumann condition $\ \partial_{\nu}u=0$ on the sphere, while the additional assumption $\ v^{\prime}\left(0\right)=0$ is required as a consequence of the regularity of the radially symmetric solutions $u.$

There are many contributions to radial solutions for boundary value problems in the annulus and in the ball. For instance, in papers [5], [12], [14] and [19] there is considered the case of equations and systems with the classical Laplacian, the papers [6], [8] and [13] deal with the $p$ -Laplacian, and in [3], [7] and [16] it is considered the case of the $\phi$ -Laplacian, in particular, that of the mean curvature operators in the Euclidean and Minkowski spaces. The methods that are used are of the most spilled: fixed point principles, topological degree, upper and lower solution techniques, variational methods and shooting method. Although the problem of radial solutions returns to ordinary differential equations, the presence of a singularity at the origin makes the study more difficult. The analysis is even more difficult with the Neumann problem due mainly to the absence of an explicit expression of the solution operator (the integral type inverse of the differential operator).

In this paper, to our knowledge, the first consecrated to the localization of radial solutions for the Neumann problem involving a general $\phi$ -Laplace operator, we use the homotopy technique - already applied in [15] and [16] for the Dirichlet problem - to obtain the existence of solutions $v$ such that

\beta<\min_{\left[R_{0},R\right]}v,\ \ \ \max_{\left[R_{0},R\right]}v<\alpha,

for two given numbers $0<\beta<\alpha.$

From a physical point of view, assuming that function $v$ stands for the state of a process and $f$ is the external source, such a localization is imposed by two requirements: First, from the necessity to find the state-depending source $f\left(r,v\right)$ (feedback law) in order to guarantee that the state $v$ remains bounded between two given bounds and secondly, the state-depending source $f$ being given, to find the bounds of its corresponding state $v.$

Mathematically, such a localization immediately gives multiple solutions in the case of oscillating functions $f\left(r,s\right).$ Additionally, we show that the solutions $v$ are decreasing on $\left[R_{0},R\right]$ provided that $f\left(r,s\right)$ has suitable monotonicity properties in $r$ and $s.$ Also, a certain behavior of the decreasing solution is emphasized in terms of a Harnack type inequality which is established by a variable change meant to eliminate the first-order term of the differential operator. Our results apply particularly to the singular homeomorphisms

\phi:\left(-a,a\right)\rightarrow\mathbb{R},\ \ \phi\left(s\right)=\frac{s}{% \sqrt{a^{2}-s^{2}}}\ \ \ \ \ \text{ }(\text{here \ }a<+\infty,\ b=+\infty),

as in the mean curvature operator in the Minkowski space.

2. Solution properties

By a solution of (1.2) we mean a function $v\in C^{1}\left[R_{0},R\right]$ such that $\ v^{\prime}\in\left(-a,a\right)\$ and $\ r^{n-1}\phi\left(v^{\prime}\right)$ is differentiable and satisfies (1.2).

We look for solutions which are nonnegative on $\left[R_{0},R\right].$

First, in this section, we consider the case of external influence terms that do not depend on the state variable, namely the problem

(2.1)

\left\{\begin{array}[]{ll}L\left(v\right)\left(r\right):=-\left(r^{n-1}\phi% \left(v^{\prime}\right)\right)^{\prime}+\varepsilon r^{n-1}v=r^{n-1}h\left(r% \right)&\text{in }\left(R_{0},R\right)\\ v^{\prime}\left(R_{0}\right)=v^{\prime}(R)=0&\end{array}\right.

for a given function $h\in C\left[R_{0},R\right].$

2.1. Some basic properties

The next lemma gives a characterization of the solutions.

Lemma 2.1.

A function $v$ is a solution of (2.1) if and only if $v\in C\left[R_{0},R\right]$ and the following two conditions hold:

(2.2)

v\left(r\right)=v\left(R_{0}\right)+\int_{R_{0}}^{r}\phi^{-1}\left(s^{1-n}\int% _{s}^{R}\tau^{n-1}\left(h-\varepsilon v\right)d\tau\right)ds\ \ \ \ (r\in\left% [R_{0},R\right]);

(2.3)

\int_{R_{0}}^{R}\tau^{n-1}\left(h-\varepsilon v\right)d\tau=0.

Proof.

Let $v$ be a solution of (2.1). Then integrating from $r$ to $R$ by taking into account that $v^{\prime}\left(R\right)=0$ gives

r^{n-1}\phi\left(v^{\prime}\right)=\int_{r}^{R}\tau^{n-1}\left(h-\varepsilon v% \right)d\tau.

which immediately implies (2.2). Next (2.3) is obtained by integration in (2.1) on $\left[R_{0},R\right]$ and using $v^{\prime}\left(R_{0}\right)=v^{\prime}(R)=0.$

Conversely, if $v$ satisfies (2.2) and (2.3), then clearly $v\in C^{1}\left(R_{0},R\right)$ and by direct computation using (2.3),

	$\displaystyle v^{\prime}\left(r\right)$	$\displaystyle=$	$\displaystyle\phi^{-1}\left(r^{1-n}\int_{r}^{R}\tau^{n-1}\left(h-\varepsilon v% \right)d\tau\right)$
		$\displaystyle=$	$\displaystyle-\phi^{-1}\left(r^{1-n}\int_{R_{0}}^{r}\tau^{n-1}\left(h-% \varepsilon v\right)d\tau\right),$

whence $v^{\prime}\left(R\right)=v^{\prime}\left(R_{0}\right)=0.$ Moreover,

(2.5)

\phi\left(v^{\prime}\right)=r^{1-n}\int_{r}^{R}\tau^{n-1}\left(h-\varepsilon v% \right)d\tau,

which shows that $\phi\left(v^{\prime}\right)$ is differentiable and yields

-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}=r^{n-1}\left(h-% \varepsilon v\right)\ \ \ \text{ in }\left(R_{0},R\right).\text{ }

Thus, $v$ is a solution of (2.1). ∎

Lemma 2.2.

Let $h_{1},h_{2}\in C\left[R_{0},R\right],$ $h_{1}\leq h_{2}$ on $\left[R_{0},R\right],$ and let $v_{1},v_{2}\in C^{1}\left[R_{0},R\right]$ be such that for $i=1,2,$ one has $\ v_{i}^{\prime}\left(R_{0}\right)=v_{i}^{\prime}\left(R\right)=0$ and

L\left(v_{i}\right)\left(r\right)=r^{n-1}h_{i}\left(r\right)\ \ \ \text{for }r% \in\left(R_{0},R\right).

Then $\ v_{1}\leq v_{2}$ in $\left[R_{0},R\right].$

For each $\ h\in C\left[R_{0},R\right],$ problem (2.1) has at most one solution which is nonnegative if $h\geq 0.$ Moreover, $h/\varepsilon$ is the unique solution for any constant $h.$

Proof.

Assume the contrary. Let $I=\left(\alpha,\beta\right)$ be a maximal subinterval of $\left(R_{0},R\right)$ on which $v_{2}-v_{1}$ is strictly negative. On $I$ we then have

-\left(r^{n-1}\left(\phi\left(v_{2}^{\prime}\right)-\phi\left(v_{1}^{\prime}% \right)\right)\right)^{\prime}=r^{n-1}\left(h_{2}-h_{1}\right)-\varepsilon r^{% n-1}\left(v_{2}-v_{1}\right)>0.

Hence the function

\digamma\left(r\right):=r^{n-1}\left(\phi\left(v_{2}^{\prime}\left(r\right)% \right)-\phi\left(v_{1}^{\prime}\left(r\right)\right)\right)

is strictly decreasing on $I.$

Assume $\beta=R.$ Since $v_{i}^{\prime}(R)=0$ for both $i=1,2,$ we must have $\ \phi\left(v_{2}^{\prime}\left(r\right)\right)-\phi\left(v_{1}^{\prime}\left(% r\right)\right)>0,$ whence $\ v_{2}^{\prime}\left(r\right)-v_{1}^{\prime}\left(r\right)>0$ for $r\in I.$ This shows that the function $\ v_{2}-v_{1}$ is strictly increasing on $I.$ This together with its negativity implies that $\ v_{2}-v_{1}$ is negative on the whole interval $[R_{0},R).$ Hence $\alpha=R_{0}$ and consequently $\digamma$ is strictly decreasing on $\left[R_{0},R\right],$ which is impossible since its values at $R_{0}$ and $R$ are equal to zero. Thus $\beta<R$ and $v_{2}\left(\beta\right)-v_{1}\left(\beta\right)=0.$

Assume $\alpha=R_{0}.$ Since $v_{i}^{\prime}(R_{0})=0$ for both $i=1,2,$ we have $\phi\left(v_{2}^{\prime}\left(r\right)\right)-\phi\left(v_{1}^{\prime}\left(r% \right)\right)<0,$ whence $\ v_{2}^{\prime}\left(r\right)-v_{1}^{\prime}\left(r\right)<0\$ for $\ r\in I.$ But this is impossible in virtue of $v_{2}\left(\beta\right)-v_{1}\left(\beta\right)=0$ and the negativity of $\ v_{2}-v_{1}$ on $I.$ Therefore $R_{0}<\alpha<\beta<R$ and $\ v_{2}\left(\alpha\right)-v_{1}\left(\alpha\right)=0,\ v_{2}\left(\beta\right% )-v_{1}\left(\beta\right)=0.$ Let $\ r_{0}\in\left(\alpha,\beta\right)$ be such that $\ v_{2}\left(r_{0}\right)-v_{1}\left(r_{0}\right)=\min_{r\in\left[\alpha,\beta% \right]}\left(v_{2}\left(r\right)-v_{1}\left(r\right)\right).$ Then $\digamma\left(r_{0}\right)=0$ and since $\digamma$ is decreasing we must have $\digamma\left(r\right)>0$ on $\left(\alpha,r_{0}\right)$ and $\digamma\left(r\right)<0$ on $\left(r_{0},\beta\right).$ Consequently, the function $v_{2}-v_{1}$ should be increasing on $\left(\alpha,r_{0}\right)$ , which is impossible since $v_{2}\left(\alpha\right)-v_{1}\left(\alpha\right)=0$ and $v_{2}\left(r_{0}\right)-v_{1}\left(r_{0}\right)<0.$

Therefore, $v_{1}\leq v_{2}$ on $\left[R_{0},R\right]$ as desired.

Let $v_{1},v_{2}$ solve (2.1) for the same function $h.$ Applying the previous conclusion to $h_{1}=h_{2}=h$ gives $v_{1}\leq v_{2}$ and $v_{2}\leq v_{1}.$ Thus $v_{1}=v_{2}$ proving the uniqueness. The fact that $v=h/\varepsilon\$ for any constant function $h$ follows directly from (2.1). Finally the property that the solution $\ v$ is nonnegative if $h$ is nonnegative follows from the comparison property and the fact that zero solves the problem for $h=0.$ ∎

Lemma 2.3.

If $v$ solves the problem for a decreasing function $h\in C\left[R_{0},R\right],$ then $v$ is decreasing in $\left[R_{0},R\right]$ too.

Proof.

Assume otherwise. Then there is a maximal subinterval $\left[\alpha,\beta\right]$ of $\left[R_{0},R\right]$ on which $v:=S\left(h\right)$ is strictly increasing. If $\alpha$ is interior, i.e., $\alpha>R_{0},$ then clearly $v^{\prime}\left(\alpha\right)=0.$ Otherwise $\alpha=R_{0}$ and $v^{\prime}\left(\alpha\right)=0$ due to the Neumann condition. Similarly, $v^{\prime}\left(\beta\right)=0.$ The function $h-\varepsilon v$ being decreasing on $\left[\alpha,\beta\right],$ one has that $\ r^{1-n}\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}$ is increasing on $\left(\alpha,\beta\right).$ Hence there are only two possibilities: either (a) $\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}\geq 0$ on $\left(\alpha,\beta\right),$ or (b) there is $\gamma\in(\alpha,\beta]$ such that $\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}<0$ in $\left(\alpha,\gamma\right)$ and $\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}\geq 0$ in $\left(\gamma,\beta\right).$

In case (a), the monotonicity of $r^{1-n}\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}$ on $\left(\alpha,\beta\right)$ implies that $\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}$ is increasing on $\left(\alpha,\beta\right),$ which implies the convexity on $\left(\alpha,\beta\right)$ of the function $r^{n-1}\phi\left(v^{\prime}\right).$ Since this function vanishes at $\alpha$ and $\beta$ (like $v^{\prime}$ ), we must have $r^{n-1}\phi\left(v^{\prime}\right)\leq 0$ in $\left(\alpha,\beta\right).$ But this gives $v^{\prime}\leq 0$ in $\left(\alpha,\beta\right),$ which contradicts our assumption on $v.$

Assume case (b). Then the function $r^{n-1}\phi\left(v^{\prime}\right)$ is decreasing on $\left(\alpha,\gamma\right).$ Since its value at $\alpha$ is zero, we have $\ r^{n-1}\phi\left(v^{\prime}\right)\leq 0$ in $\left(\alpha,\gamma\right),$ whence $\ v^{\prime}\leq 0$ in $\left(\alpha,\gamma\right),$ again a contradiction. ∎

2.2. A Harnack type inequality

In this section we assume that the homeomorphism $\phi$ satisfies the following condition:

(H_ϕ):

$\phi$ is $C^{1},$

(2.6)

s\phi^{\prime}\left(s\right)\leq\phi\left(s\right)\ \ \ \text{and\ \ \ }\phi^{% \prime}\left(s\right)\geq\sigma>0\ \ \ \text{for all }s\in(-a,0],

for some $\sigma>0.$

For example, such homeomorphisms are those involved by the classical Laplacian and the mean curvature operator in the Minkowski space.

Let $\ R_{1}\in\left(R_{0},R\right)$ be a fixed number and $\ v\in W^{2,\infty}\left(R_{0},R\right)\cap C^{1}\left[R_{0},R\right]$ be nonnegative on $\left[R_{0},R\right],$ decreasing on $\left[R_{1},R\right],$ with $\ v^{\prime}\left(R\right)=0$ and

(2.7)

L\left(v\right)=-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}+% \varepsilon r^{n-1}v\geq 0\ \ \ \text{a.e. in\ }\left(R_{0},R\right).

We make the change of variable $\ t=\eta\left(r\right),$ where

\eta\left(r\right)=\left\{\begin{array}[]{ll}\ln\frac{R}{r}&\text{for }n=2\\ \frac{r^{2-n}-R^{2-n}}{n-2}&\text{for }n\geq 3\end{array}\right.

by witch the interval $(R_{0},R]$ of $r$ becomes $[0,t_{1})$ for $t,$ where $\ t_{1}=\ln\left(R/R_{0}\right)$ for $n=2$ and $\ t_{1}=\left(R_{0}^{2-n}-R^{2-n}\right)/\left(n-2\right)$ for $n\geq 3.$ Note that $t_{1}=+\infty$ if $R_{0}=0.$ Also, $R_{1}$ becomes $t_{0}:=\eta\left(R_{1}\right).$ Clearly $\ 0<t_{0}<t_{1}.$ Then, letting $\ w\left(t\right)=v\left(r\right)$ and using

	$\displaystyle v^{\prime}\left(r\right)$	$\displaystyle=$	$\displaystyle-r^{1-n}w^{\prime}\left(t\right),\ \ \$
	$\displaystyle v^{\prime\prime}\left(r\right)$	$\displaystyle=$	$\displaystyle r^{2\left(1-n\right)}w^{\prime\prime}\left(t\right)+\left(n-1% \right)r^{-n}w^{\prime}\left(t\right)=r^{2\left(1-n\right)}w^{\prime\prime}% \left(t\right)-\left(n-1\right)r^{-1}v^{\prime}\left(r\right)$

and

-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}=-r^{n-1}\phi^{\prime% }\left(v^{\prime}\right)v^{\prime\prime}-\left(n-1\right)r^{n-2}\phi\left(v^{% \prime}\right),

we rewrite (2.7) as

-r^{n-1}\phi^{\prime}\left(v^{\prime}\right)\left\{r^{2\left(1-n\right)}w^{% \prime\prime}\left(t\right)-\left(n-1\right)r^{-1}v^{\prime}\left(r\right)% \right\}-\left(n-1\right)r^{n-2}\phi\left(v^{\prime}\right)+\varepsilon r^{n-1% }w\geq 0,

or equivalently

-r^{1-n}\phi^{\prime}\left(v^{\prime}\right)w^{\prime\prime}\left(t\right)+% \left(n-1\right)r^{n-2}\left\{v^{\prime}\phi^{\prime}\left(v^{\prime}\right)-% \phi\left(v^{\prime}\right)\right\}+\varepsilon r^{n-1}w\geq 0.

Since $\ v^{\prime}\leq 0$ on $\left[R_{1},R\right],$ in virtue of (2.6), one has $\ v^{\prime}\phi^{\prime}\left(v^{\prime}\right)-\phi\left(v^{\prime}\right)% \leq 0.$ It follows that

-r^{1-n}\phi^{\prime}\left(v^{\prime}\right)w^{\prime\prime}\left(t\right)+% \varepsilon r^{n-1}w\left(t\right)\geq 0\ \ \ \text{for a.e. }t\in(0,t_{0}].

Hence

-w^{\prime\prime}\left(t\right)+\varepsilon\frac{r^{2\left(n-1\right)}}{\phi^{% \prime}\left(v^{\prime}\right)}w\left(t\right)\geq 0\ \ \ \text{for a.e. }t\in% (0,t_{0}].

Since $\ w\geq 0,$ $r\leq R\$ and $\ \phi^{\prime}\left(v^{\prime}\right)\geq\sigma,\$ we deduce that

(2.8)

-w^{\prime\prime}\left(t\right)+\frac{\varepsilon R^{2\left(n-1\right)}}{% \sigma}w\left(t\right)\geq 0\ \ \ \text{for a.e. }t\in\text{ }(0,t_{0}].

Notice that in virtue of $\ v^{\prime}\left(r\right)=-r^{1-n}w^{\prime}\left(t\right),$ $\ w$ is increasing on $[0,t_{0}]$ and since $v^{\prime}\left(R\right)=0,$ one has $\ w^{\prime}\left(0\right)=0.$

Now if we first integrate in (2.8) from $0$ to $t$ $\left(t\leq t_{0}\right)$ obtaining

w^{\prime}\left(t\right)\leq\frac{\varepsilon R^{2\left(n-1\right)}}{\sigma}% \int_{0}^{t}w\left(s\right)ds\leq\frac{\varepsilon R^{2\left(n-1\right)}}{% \sigma}w\left(t_{0}\right)t\ \ \ \text{for\ \ }t\leq t_{0},

and again from $0$ to $t_{0},$ we find that

w\left(t_{0}\right)-w\left(0\right)\leq\frac{\varepsilon R^{2\left(n-1\right)}% }{2\sigma}t_{0}^{2}w\left(t_{0}\right).

Hence

w\left(0\right)\geq\left(1-\frac{\varepsilon R^{2\left(n-1\right)}}{2\sigma}t_% {0}^{2}\right)w\left(t_{0}\right).

Letting

\gamma:=1-\frac{\varepsilon R^{2\left(n-1\right)}}{2\sigma}t_{0}^{2},

assuming that $\gamma>0$ (which happens for small enough $\varepsilon$ ) and recalling that $w\left(0\right)=\min_{t\in[0,t_{0}]}w\left(t\right),$ $w\left(t_{0}\right)=\max_{t\in\left[0,t_{0}\right]}w\left(t\right),$ we have

\min_{t\in[0,t_{0}]}w\left(t\right)\geq\gamma\max_{t\in\left[0,t_{0}\right]}w% \left(t\right).

Coming back to function $v,$ we have the Harnack inequality

v\left(R\right)=\min_{r\in\left[R_{1},R\right]}v\left(r\right)\geq\gamma\max_{% r\in[R_{1},R]}v\left(r\right)=\gamma v\left(R_{1}\right).

Notice that in the case $R_{0}>0,$ we may take $R_{1}=R_{0}$ and so $t_{0}=\psi\left(R_{1}\right)$ is finite and the above reasoning remains true yielding to the better inequality

v\left(R\right)=\min_{r\in\left[R_{0},R\right]}v\left(r\right)\geq\gamma\max_{% r\in[R_{0},R]}v\left(r\right)=\gamma v\left(R_{0}\right).

Thus we have the following result.

Theorem 2.4.

Assume that condition (H_ϕ) holds. Then for every number $\ R_{1}\in[R_{0},R)$ with $\ \eta\left(R_{1}\right)<+\infty,$ there exists $\ \varepsilon_{0}=\varepsilon_{0}\left(R_{1},\sigma,R\right)>0$ such that for every $0<\varepsilon<\varepsilon_{0},$ there is a constant $\ \gamma=\gamma\left(R_{1},\sigma,R,\varepsilon\right)>0$ such that

v\left(R\right)\geq\gamma v\left(R_{1}\right)

for every $\ v\in W^{2,\infty}\left(R_{0},R\right)\cap C^{1}\left[R_{0},R\right]$ nonnegative on $\left[R_{0},R\right],$ decreasing on $\left[R_{1},R\right],$ with $\ \ v^{\prime}\left(R\right)=0$ and $\ \ L\left(v\right):=-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{\prime}% +\varepsilon r^{n-1}v\geq 0\$ a.e. in $\left(R_{0},R\right).$

3. Existence, localization and multiplicity

3.1. The solution operator

Let $K_{+}$ be the positive cone of $C\left[R_{0},R\right].$

In what follows we assume that for every function $h\in K_{+},$ problem (2.1) has a solution. This hapens for example in case of singular homeomorphisms as shows Corollary 2.4 in [3]. Then, in virtue of the previous lemmas, the solution operator $S:K_{+}\rightarrow K_{+}$ attaching to each $h$ the corresponding (unique) solution $v$ is well-defined and isotone. In addition, one has

Lemma 3.1.

The solution operator $S$ is completely continuous from $K_{+}$ to $K_{+}.$

Proof.

(a) $S\left(M\right)$ is relatively compact for every bounded set $M\subset K_{+}.$ Indeed, if $C>0$ is such that $\ \left|h\right|_{\infty}=\max_{r\in\left[R_{0},R\right]}\left|h\left(r\right)% \right|\leq C\$ for all $h\in M,$ then from $0\leq h\leq C$ one has $0\leq S\left(h\right)\leq S\left(C\right).$ Hence $\left|S\left(h\right)\right|_{\infty}\leq\left|S\left(C\right)\right|_{\infty}.$ Thus $S\left(M\right)$ is bounded in $K_{+}.$ Now from (2.1) we see that the derivatives of the functions $v$ from $S\left(M\right)$ are uniformly bounded, that is $S\left(M\right)$ is equicontinuous. Therefore $S\left(M\right)$ is relatively compact in $C\left[R_{0},R\right].$

(b) $S$ is continuous. Let $\ h_{k}\in K_{+}$ be convergent to some $\ h$ and let $\ v_{k}=S\left(h_{k}\right).$ We need to prove that $\ v_{k}\rightarrow S\left(h\right).$ According to Lemma 3.1 there is a convergent subsequence of $\left(v_{k}\right).$ Let $v$ be its limit. Passing to the limit in (2.2) and (2.3) written for $h_{k}$ and $v_{k},$ we find that $S\left(h\right)=v.$ As a result the whole sequence $\left(v_{k}\right)$ converges to $S\left(h\right).$ ∎

Now it is clear that $v$ is a nonnegative solution of (1.2) if and only if $v$ is a fixed point of the operator

T:K_{+}\rightarrow K_{+},\ \ \ T=S\circ N_{f},

where $N_{f}\left(v\right)=f\left(\cdot,v\left(\cdot\right)\right)$ is the Nemytski operator associated to $f.$ According to the previous lemmas about the solution operator, the operator $T$ is well-defined and completely continuous.

3.2. Existence and localization

Now, for any number $\alpha>0,$ consider the set

V_{\alpha}:=\left\{v\in K_{+}:\ \left|v\right|_{\infty}<\alpha\right\}.

The operator $T$ being completely continuous, the set $T\left(\overline{V}_{\alpha}\right)$ is bounded, so there is a number $\widetilde{\alpha}\geq\alpha$ such that $T\left(\overline{V}_{\alpha}\right)\subset\overline{V}_{\widetilde{\alpha}}.$ Define the modified operator $\widetilde{T}:\overline{V}_{\widetilde{\alpha}}\rightarrow\overline{V}_{% \widetilde{\alpha}}$ by

\widetilde{T}\left(v\right)=T\left(\min\left\{\frac{\alpha}{\left|v\right|_{% \infty}},1\right\}v\right).

The following two lemmas rely on the properties of the fixed point index (see, a.e., [9]).

Lemma 3.2.

(3.1)

T\left(v\right)\neq\lambda v\ \ \ \text{for }v\in K_{+}\text{ with }\left|v% \right|_{\infty}=\alpha\text{ and }\lambda\geq 1,

then the fixed point index $\ i\left(\widetilde{T},V_{\alpha},\overline{V}_{\widetilde{\alpha}}\right)=1.$

Next, denote $\left|v\right|_{0}:=\min_{r\in\left[R_{0},R\right]}v\left(r\right)\$ and for a number $\beta>0,$ consider the set

W_{\beta}:=\left\{v\in\overline{V}_{\widetilde{\alpha}}:\ \left|v\right|_{0}<% \beta\right\}.

It is clear that $W_{\beta}$ is open in $\overline{V}_{\widetilde{\alpha}}.$

Lemma 3.3.

Assume that for a function $h\in K_{+}$ such that $\left|h\right|_{\infty}=\alpha,$ $\left|h\right|_{0}>\beta,$ one has

(3.2)

\left(1-\lambda\right)T\left(v\right)+\lambda h\neq v\ \ \text{for }v\in K_{+}% \text{ with }\left|v\right|_{\infty}\leq\alpha,\ \left|v\right|_{0}=\beta\text% { and }\lambda\in\left[0,1\right].

Then $\ i\left(\widetilde{T},W_{\beta},\overline{V}_{\widetilde{\alpha}}\right)=0.$

Lemma 3.4.

Under the assumptions of Lemmas 3.2 and 3.3, the operator $T$ has a fixed point $v$ in $\ V_{\alpha}\setminus\overline{W}_{\beta},$ that is problem (1.2) has a solution $\ v$ which is nonnegative on $\left[R_{0},R\right],$ with $\ \beta<\left|v\right|_{0}\$ and $\ \left|v\right|_{\infty}<\alpha.$

Proof.

One has

	$\displaystyle 1$	$\displaystyle=$	$\displaystyle i\left(\widetilde{T},V_{\alpha},\overline{V}_{\widetilde{\alpha}% }\right)=i\left(\widetilde{T},V_{\alpha}\setminus\overline{W}_{\beta},% \overline{V}_{\widetilde{\alpha}}\right)+i\left(\widetilde{T},V_{\alpha}\cap W% _{\beta},\overline{V}_{\widetilde{\alpha}}\right),$
	$\displaystyle 0$	$\displaystyle=$	$\displaystyle i\left(\widetilde{T},W_{\beta},\overline{V}_{\widetilde{\alpha}}% \right)=i\left(\widetilde{T},W_{\beta}\setminus\overline{V}_{\alpha},\overline% {V}_{\widetilde{\alpha}}\right)+i\left(\widetilde{T},V_{\alpha}\cap W_{\beta},% \overline{V}_{\widetilde{\alpha}}\right).$

Subtracting gives

(3.3)

i\left(\widetilde{T},V_{\alpha}\setminus\overline{W}_{\beta},\overline{V}_{% \widetilde{\alpha}}\right)-i\left(\widetilde{T},W_{\beta}\setminus\overline{V}% _{\alpha},\overline{V}_{\widetilde{\alpha}}\right)=1.

Hence at least one of the numbers $i\left(\widetilde{T},V_{\alpha}\setminus\overline{W}_{\beta},\overline{V}_{% \widetilde{\alpha}}\right)$ and $i\left(\widetilde{T},W_{\beta}\setminus\overline{V}_{\alpha},\overline{V}_{% \widetilde{\alpha}}\right)$ is nonzero. We claim that the last one equals zero. Indeed, otherwise there would exist $v\in W_{\beta}\setminus\overline{V}_{\alpha}$ with $\widetilde{T}\left(v\right)=v,$ that is

T\left(\frac{\alpha}{\left|v\right|_{\infty}}v\right)=v,

or equivalently $T\left(\omega\right)=\lambda\omega,$ where $\omega=\frac{\alpha}{\left|v\right|_{\infty}}v$ and $\lambda=\frac{\left|v\right|_{\infty}}{\alpha}.$ Since $\left|\omega\right|_{\infty}=\alpha$ and $\lambda>1$ we arrived to a contradiction with (3.1). Therefore $i\left(\widetilde{T},W_{\beta}\setminus\overline{V}_{\alpha},\overline{V}_{% \widetilde{\alpha}}\right)=0$ and from (3.3) one has $i\left(\widetilde{T},V_{\alpha}\setminus\overline{W}_{\beta},\overline{V}_{% \widetilde{\alpha}}\right)=1,$ which implies our conclusion. ∎

We are now ready to state and prove our main existence and localization result.

For any numbers $0<\beta<\alpha,$ denote

	$\displaystyle m_{\alpha,\beta}$	$\displaystyle:$	$\displaystyle=\min\left\{f\left(r,s\right):\ r\in\left[R_{0},R\right],\ s\in% \left[\beta,\alpha\right]\right\},$
	$\displaystyle M_{\alpha}$	$\displaystyle:$	$\displaystyle=\max\left\{f\left(r,s\right):\ r\in\left[R_{0},R\right],\ s\in% \left[0,\alpha\right]\right\}.$

Theorem 3.5.

If for two positive numbers $\alpha,\beta$ satisfying $\alpha>\beta,$ the following conditions

(h1):: $M_{\alpha}<\varepsilon\alpha,$
(h2):: $m_{\alpha,\beta}>\varepsilon\beta$

hold, then problem (1.2) has a positive solution $\ v$ such that

\beta<\left|v\right|_{0}\ \ \ \text{and\ \ \ }\left|v\right|_{\infty}<\alpha.

Proof.

First, we remark that inequality $\alpha>\beta$ guarantees the existence of a function $h\in K_{+}$ such that $\left|h\right|_{\infty}=\alpha$ and $\left|h\right|_{0}>\beta$ as needed in Lemma 3.3. Such a function is the constant $h=\alpha.$

Assume that for some $v\in K_{+}$ with $\left|v\right|_{\infty}=\alpha$ and some $\lambda\geq 1,$ one has $T\left(v\right)=\lambda v.$ Then $v\leq\lambda v=S\left(N_{f}\left(v\right)\right)\leq S\left(M_{\alpha}\right)=% M_{\alpha}/\varepsilon,$ so $\ \alpha=\left|v\right|_{\infty}\leq M_{\alpha}/\varepsilon,$ which contradicts (h1). Hence Lemma 3.2 applies.

Assume that for some $v\in K_{+}$ with $\left|v\right|_{\infty}\leq\alpha,\ \left|v\right|_{0}=\beta$ and $\lambda\in\left[0,1\right]$ we have $\left(1-\lambda\right)T\left(v\right)+\lambda\alpha=v.$ Clearly,

T\left(v\right)=S\left(N_{f}\left(v\right)\right)\geq S\left(m_{\alpha,\beta}% \right)=\frac{m_{\alpha\beta}}{\varepsilon},

hence according to (h2),

\left|T\left(v\right)\right|_{0}\geq\frac{m_{\alpha\beta}}{\varepsilon}>\beta.

Then

\beta=\left|v\right|_{0}=\left|\left(1-\lambda\right)T\left(v\right)+\lambda% \alpha\right|_{0}>\beta.

Hence Lemma 3.3 also applies. The conclusion now follows from Lemma 3.4. ∎

Remark 3.1.

(a) If we assume that for each $r\in\left[R_{0},R\right],$ the function $f\left(r,.\right)$ is increasing on $\left[0,\alpha\right],$ then

m_{\alpha,\beta}=\min\left\{f\left(r,\beta\right):\ r\in\left[R_{0},R\right]% \right\},\ \ \ M_{\alpha}=\max\left\{f\left(r,\alpha\right):\ r\in\left[R_{0},% R\right]\right\}.

(b) If $f\left(r,s\right)=a\left(r\right)g\left(s\right),$ where $a$ is continuous and positive on $\left[R_{0},R\right]$ and $g$ is increasing on $\left[0,\alpha\right],$ then

m_{\alpha,\beta}=m_{a}g\left(\beta\right),\ \ \ M_{\alpha}=M_{a}g\left(\alpha% \right),

where $\ m_{a}=\min_{\left[R_{0},R\right]}a\left(r\right),\ M_{a}=\max_{\left[R_{0},R% \right]}a\left(r\right).$

3.3. Decreasing solutions

Here assume the following monotonicity properties of $f:$

(H_f):: $f\left(.,s\right)$ is decreasing in $\left[R_{0},R\right]$ for each $s\in\mathbb{R}_{+}$ and $f\left(r,.\right)$ is increasing in $\mathbb{R}_{+}$ for each $r\in\left[R_{0},R\right].$

Under this condition, if a nonnegative function $\ v$ is decreasing on $\left[R_{0},R\right],$ then the function $\ N_{f}\left(v\right)=f\left(.,v\left(.\right)\right)\$ is decreasing too. Thus, if we consider the sub-cone $K$ of $K_{+}$ defined by

K:=\left\{v\in K_{+}:\ v\text{ is decreasing on }\left[R_{0},R\right]\right\},

then in view of Lemma 2.3, we have $\ T\left(K\right)\subset K$ and we can apply the reasoning from the proof of Theorem 3.5, working in $K$ instead of $K_{+}.$ In this way, the existence of a decreasing solution is obtained. Using in addition Theorem 2.4, we obtain the following result.

Theorem 3.6.

Assume that conditions (H_ϕ) and (H_f) hold and that $R_{1},\ \varepsilon$ and $\ \gamma$ are as in Section 2.2. If for two numbers $\ 0<\beta<\alpha$ one has

(h1’):: $f\left(R_{0},\alpha\right)<\varepsilon\alpha;$
(h2’):: $f\left(R,\beta\right)>\varepsilon\beta,$

then problem (1.2) has a decreasing positive solution $v$ such that

\beta<v\left(R\right),\ \ \ \ v\left(R_{0}\right)<\alpha,

(3.4)

v\left(R\right)\geq\gamma v\left(R_{1}\right).

Remark 3.2.

Inequality $(\ref{rap})$ gives us the bound $1/\gamma,$ independent on $\alpha$ and $\beta,$ for the ratio $v\left(R_{1}\right)/v\left(R\right)$ between the maximum and the minimum of $v$ on the interval $\left[R_{1},R\right].$ Thus, if for such a solution $v,$ $v\left(R_{1}\right)$ is large, say $\ v\left(R_{1}\right)>k,$ then its minimum $v\left(R\right)$ is larger than $\ \gamma k;$ if its minimum $v\left(R\right)$ is small, say $\ v\left(R\right)<1/k,$ then $v\left(R_{1}\right)$ is smaller than $\ 1/\left(\gamma k\right).$ As noted above, in the case of the annulus, i.e., for $R_{0}>0,$ one may take $R_{1}=R_{0}$ and then $\ 1/\gamma$ is a bound for the ratio between the maximum and the minimum of $v$ on the whole interval $\left[R_{0},R\right].$

3.4. Multiple solutions

We first give a three-solution result.

Theorem 3.7.

Under the assumptions of Theorem 3.5, if in addition there exists $\ \alpha_{0}\in\left(0,\beta\right)$ such that

M_{\alpha_{0}}<\varepsilon\alpha_{0},

then problem (1.2) has at least three nonnegative solutions solutions $\ v_{1},v_{2},v_{3}$ such that

$\displaystyle\beta$	$\displaystyle<$	$\displaystyle\left\|v_{1}\right\|_{0},\ \ \ \left\|v_{1}\right\|_{\infty}<\alpha;$
$\displaystyle\left\|v_{2}\right\|_{\infty}$	$\displaystyle<$	$\displaystyle\alpha_{0};$
$\displaystyle\left\|v_{3}\right\|_{0}$	$\displaystyle<$	$\displaystyle\beta,\ \ \ \left\|v_{3}\right\|_{\infty}>\alpha_{0}.$

Proof.

Solution $v_{1}$ is guaranteed by Theorem 3.5. Next from $\ i\left(\widetilde{T},V_{\alpha_{0}},\overline{V}_{\widetilde{\alpha}}\right)=1$ we obtain the solution $\ v_{2}.$ Now, let us remark that $\ \overline{V}_{\alpha_{0}}\subset W_{\beta}.\$ Indeed, if $\ v\in\overline{V}_{\alpha_{0}}$ then $\ \left|v\right|_{\infty}\leq\alpha_{0}<\beta$ and so $\ \left|v\right|_{0}\leq\left|v\right|_{\infty}<\beta.$ Hence $\ v\in W_{\beta}.$ This inclusion implies

i\left(\widetilde{T},W_{\beta}\setminus\overline{V}_{\alpha_{0}},\overline{V}_% {\widetilde{\alpha}}\right)=i\left(\widetilde{T},W_{\beta},\overline{V}_{% \widetilde{\alpha}}\right)-i\left(\widetilde{T},V_{\alpha_{0}},\overline{V}_{% \widetilde{\alpha}}\right)=0-1=-1,

whence the existence of $\ v_{3}.$ ∎

Obviously, the solution $v_{2}$ can be zero. However, this is not the case if $f\left(.,0\right)\neq 0.$

Next we establish the existence of an arbitrary number of solutions, or of a sequence of solutions, by assuming a strong oscillation in $s$ of nonlinearity $f\left(r,s\right).$

Theorem 3.8.

$(1^{0})$ Let $\ \left(\alpha_{i}\right)_{1\leq i\leq k},\ \left(\beta_{i}\right)_{1\leq i% \leq k}$ $\left(k\leq+\infty\right)$ be increasing finite or infinite sequences of positive numbers with $\ \beta_{i}<\alpha_{i}\leq\beta_{i+1}$ for all $i.$ If the assumptions of Theorem 3.5 are satisfied for each couple $\ \left(\alpha_{i},\beta_{i}\right),$ then problem (1.2) has $k$ (respectively, when $k=+\infty,$ an infinite sequence of) distinct solutions $\ v_{i}$ with

(3.5)

\beta_{i}<\left|v_{i}\right|_{0},\ \ \ \left|v_{i}\right|_{\infty}<\alpha_{i}.

$\left(2^{0}\right)$ Let $\ \left(\alpha_{i}\right)_{i\geq 1},\ \left(\beta_{i}\right)_{i\geq 1}$ be decreasing infinite sequences with $\ \alpha_{i+1}\leq\beta_{i}<\alpha_{i}\$ for all $i.$ If the assumptions of Theorem 3.5 are satisfied for each couple $\ \left(\alpha_{i},\beta_{i}\right),$ then problem (1.2) has an infinite sequence of distinct solutions $\ v_{i}$ satisfying (3.5).

Proof.

Denote

K_{i}:=\left\{v\in K_{+}:\ \beta_{i}<\left|v\right|_{0},\ \ \left|v\right|_{% \infty}<\alpha_{i}\right\}.

It is sufficient to remark that $K_{i}\cap K_{i+1}=\emptyset$ for all $i.$ To prove this let us first assume that sequences $\left(\alpha_{i}\right),\left(\beta_{i}\right)$ are increasing (case ( $1^{0}$ )). Then since $\alpha_{i}\leq\beta_{i+1},$ for any $v\in K_{i}$ one has $\ \left|v\right|_{0}\leq\left|v\right|_{\infty}<\alpha_{i}\leq\beta_{i+1}.\$ Hence $\ v\notin K_{i+1}.$ Similarly, in case ( $2^{0}$ ), if $v\in K_{i},$ then $\ \left|v\right|_{\infty}\geq\left|v\right|_{0}>\beta_{i}\geq\alpha_{i+1},$ so $\ v\notin K_{i+1}.$ ∎

Remark 3.3.

Under assumptions (H_ϕ) and (H_f), if multiple decreasing solutions are obtained via Theorem 3.6, then for all of them, one has the same bound $\ 1/\gamma\$ for the ratio between their maximum and minimum on the interval $\left[R_{1},R\right].$

3.5. Existence and multiplicity under asymptotic conditions

In the situation where only the existence of solutions is of interest and not exactly their location, the asymptotic conditions on $f$ are sufficient and easier to check than the punctual conditions.

Assume here again, as in Remark 3.1 (b), the following form of $f,$

f\left(r,s\right)=a\left(r\right)g\left(s\right),

where $\ a\$ is continuous and positive on $\left[R_{0},R\right]$ and $\ g\$ is increasing on $\mathbb{R}_{+}.$

Thus the existence of two numbers $\ \alpha,$ $\beta$ with $\ \alpha>\beta\$ and satisfying (h1), (h2) obviously follows from the asymptotic conditions

\lim\inf_{\tau\rightarrow+\infty}\frac{g\left(\tau\right)}{\tau}<\frac{% \varepsilon}{M_{a}}\ \ \text{and\ }\ \lim\sup_{\tau\rightarrow 0}\frac{g\left(% \tau\right)}{\tau}>\frac{\varepsilon}{m_{a}},

respectively.

Also, two sequences $\left(\alpha_{i}\right)$ and $\left(\beta_{i}\right)$ exist as in Theorem 3.8 (1⁰) provided that

\lim\inf_{\tau\rightarrow+\infty}\frac{g\left(\tau\right)}{\tau}<\frac{% \varepsilon}{M_{a}}\ \ \text{and\ \ }\lim\sup_{\tau\rightarrow+\infty}\frac{g% \left(\tau\right)}{\tau}>\frac{\varepsilon}{m_{a}},

and as in Theorem 3.8 (2⁰) provided that

\lim\inf_{\tau\rightarrow 0}\frac{g\left(\tau\right)}{\tau}<\frac{\varepsilon}% {M_{a}},\ \ \text{\ }\lim\sup_{\tau\rightarrow 0}\frac{g\left(\tau\right)}{% \tau}>\frac{\varepsilon}{m_{a}}.

4. Numerical solutions

In order to carry out our numerical experiments we have use the MATLAB object-oriented package Chebfun. We refer only to [17] and [18] for the details on using this package although the literature on this topic is much broader.

The numerical experiments performed on a similar problem in [10], encouraged us to use this programming environment and not others. It proved to be very simple and flexible in writing a code, including in imposing the boundary conditions, an otherwise non-trivial matter. The details it provides regarding the convergence of the Newton method are extremely useful.

We present three concrete Neumann problems for which numerical solutions are obtained confirming the theoretical results.

4.1. First example

We look for a nonzero numerical solution and to confirm the theory for the Neumann boundary value problem involving the classical Laplacian

(4.1)

\left\{\begin{array}[]{ll}-\left(rv^{\prime}\right)^{\prime}+rv=r\frac{\sqrt{v% }}{r+1},&r\in\left(0,1\right)\\ v^{\prime}\left(0\right)=v^{\prime}\left(1\right).&\end{array}\right.

Here, with the notations from the previous sections, $\ n=2,\ \varepsilon=1,\ R_{0}=0,\ R=1\$ and $\ f\left(r,s\right)=\sqrt{s}/\left(r+1\right).\ \$ Notice the special form of $f,$ $f\left(r,s\right)=a\left(r\right)g\left(s\right),$ where $\ a\left(r\right)=1/\left(r+1\right)\$ is decreasing and $\ g\left(s\right)=\sqrt{s}\$ is increasing.

The theory is confirmed if a decreasing positive solution $\ v\$ and numbers $\ \alpha,$ $\beta>0,$ $\beta<\alpha\$ are found such that the following inequalities are satisfied:

	$\displaystyle m_{a}g\left(\beta\right)$	$\displaystyle>$	$\displaystyle\varepsilon\beta,\ \ \ M_{a}g\left(\alpha\right)<\varepsilon\alpha,$
	$\displaystyle\beta$	$\displaystyle<$	$\displaystyle v\left(R\right),\ \ \ v\left(R_{0}\right)<\alpha,$

which applied to the present example, for which $\varepsilon=1,$ $m_{a}=1/2,$ $M_{a}=1,$ $g\left(\alpha\right)=\sqrt{\alpha}$ and $g\left(\beta\right)=\sqrt{\beta},\$ read as:

(4.2)

\beta<0.25,\ \ \ 1<\alpha,\ \ \ \ \beta<v\left(1\right),\ \ \ \ v\left(0\right% )<\alpha.

The numerical solution $v$ is presented in Fig. 1 and the confirmation of the theory takes place, for example, with $\ \alpha=0.4$ and $\ \beta=0.35.$

Figure 1. Graph of the numerical solution of problem (4.1). The initial guess for the initialization of the Newton procedure is

v_{0}:=1.

From Fig. 2 a) we observe that the Newton method converges in five steps and has at least order two of convergence. The panel b) of the same figure shows that Chebfun use a polynomial of order $16$ whose coefficients decrease linearly to order $10^{-14}.$

Figure 2. a) The convergence of Newton method. b) The behavior of Chebyshev coefficients of solution to problem (4.1).

The residual in approximating the differential operator has been of order $10^{-11}$ and the boundary conditions have been exactly satisfied.

4.2. Second example

Here we look for a nonzero numerical solution and to confirm the theory for the Neumann boundary value problem involving the mean curvature operator in the Minkowski space,

(4.3)

\left\{\begin{array}[]{ll}-\left(r\frac{v^{\prime}}{\sqrt{1-v^{\prime 2}}}% \right)^{\prime}+rv=r\frac{\sqrt{v}}{r+1},&r\in\left(0,1\right)\\ v^{\prime}\left(0\right)=v^{\prime}\left(1\right),&\end{array}\right.

or equivalently

\left\{\begin{array}[]{ll}-v^{\prime\prime}-\frac{1}{r}\left(1-v^{\prime 2}% \right)v^{\prime}+\left(1-v^{\prime 2}\right)^{\frac{3}{2}}v=\left(1-v^{\prime 2% }\right)^{\frac{3}{2}}\frac{\sqrt{v}}{r+1},&r\in\left(0,1\right)\\ v^{\prime}\left(0\right)=v^{\prime}\left(1\right).&\end{array}\right.

Here again, $\ n=2,\ \varepsilon=1,\ R_{0}=0,\ R=1,\ \ f\left(r,s\right)=\sqrt{s}/\left(r+1\right)$ and the theory is confirmed by inequalities (4.2).

The numerical solution $\ v\$ is displayed in Fig. 3 and the confirmation of the theory takes place, for example, with $\ \alpha=1.45$ and $\ \beta=1.35.$

Figure 3. Graph of the numerical solution of problem (4.3). The initial guess for the initialization of the Newton procedure is

v_{0}:=1.

From Fig. 4a) we observe that the Newton method converges in four steps and has at least order two of convergence. The panel b) of the same figure shows that Chebfun use a polynomial of order $16$ whose coefficients decrease linearly to order $10^{-14}.$

Figure 4. a) The convergence of Newton method. b) The behavior of Chebyshev coefficients of solution to problem (4.3).

The residual in approximating the differential operator has been of order $10^{-11}$ and the boundary conditions have been exactly satisfied.

4.3. Third example

As the theory shows, the Neumann problem can have multiple positive solutions for functions $f\left(r,s\right)$ which are oscillating with respect to $s.$ To make more understandable this statement, let us first consider the simplest case of the autonomous equation (1.1), that is $f\left(r,s\right)=g\left(s\right).$ Then it is trivial to see that any constant $C$ satisfying $\ \varepsilon C=g\left(C\right)\$ is a solution of the problem

\left\{\begin{array}[]{ll}-\left(r^{n-1}\phi\left(v^{\prime}\right)\right)^{% \prime}+\varepsilon r^{n-1}v=r^{n-1}g(v)&\text{in }\left(R_{0},R\right)\\ v^{\prime}\left(R_{0}\right)=v^{\prime}(R)=0.&\end{array}\right.

Hence if the graph of $\ g\$ intersects the line of equation $\ y=\varepsilon x\$ in several points, then the problem has at least as many solutions. Therefore one obtains multiple solutions when $g$ is oscillating, here around the line $y=\varepsilon x.$ The phenomenon also occurs in the non-autonomous case, as the theory shows. Thus, for $f\left(r,s\right)=a\left(r\right)g\left(s\right),$ multiple solutions are guaranteed if $g$ oscillates up and down the lines $y=\left(\varepsilon/M_{a}\right)x$ and $y=\left(\varepsilon/m_{a}\right)x,$ respectively.

As an example of such a function, we can mention

g\left(s\right)=as+bs\sin\left(c\ln\left(s+1\right)\right),\ \ \ s\in\mathbb{R% }_{+},

where $\ a,b,c>0$ and $\ a\geq\left(c+1\right)b$ (for $g$ to be increasing). This function has a countable number of intersections with a line $\ y=\lambda x,\$ provided that $\ a-b\leq\lambda\leq a+b.\$ For numerical simulations, we choose the following values of parameters: $a=2,$ $b=c=1$ and we consider the following Neumann problem for the classical Laplacian

(4.4)

\left\{\begin{array}[]{ll}-\left(rv^{\prime}\right)^{\prime}+rv=\frac{r}{r+1}% \left(2v+v\sin\left(\ln\left(v+1\right)\right)\right),&r\in\left(0,1\right)\\ v^{\prime}\left(0\right)=v^{\prime}\left(1\right).&\end{array}\right.

Fig. 5 and 7 show two positive solutions of this problem. They were obtained by using the same procedure as before, with different initial approximations. In general, initial approximations can be suggested by the localization intervals as given by the theory.

Figure 5. Graph of a numerical solution of problem (4.4). The initial guess for the initialization of the Newton procedure is

v_{0}:=\cos\left(\pi r\right)+10^{6}.

From Fig. 6 a) we observe that the Newton method converges in six steps and the panel b) of the same figure shows that Chebfun use a polynomial of order $16$ whose coefficients decrease linearly to order $10^{-8}.$

Figure 6. a) The convergence of Newton method. b) The behavior of Chebyshev coefficients of solution to problem (4.4).

The residual in approximating the differential operator has been only of order $10^{-5}$ but the boundary conditions have been exactly satisfied.

All these last three observations lead us to state that Chebfun no longer achieves the accuracy of the previous two problems. Another solution to problem (4.4 is displayed in Fig. 7.

Figure 7. Graph of a second numerical solution of problem (4.4). The initial guess for the initialization of the Newton procedure is

v_{0}:=\cos\left(\pi r\right)+10^{2}.

To find it, Newton’s algorithm starts from a different initial guess. However, comparing these last two solutions, i.e, Fig. 5 and 7, it is observed that they have identical allures (shapes).

We can conclude that using Chebfun we have succeeded to numerically confirm, with great accuracy, some of our theoretical results regarding the existence, localization and multiplicity of positive radial solutions for Neumann problems in the ball.

References

[1] R. Bartnik, L. Simon, Spacelike hypersurfaces with prescribed boundary values and mean curvature, Commun. Math. Phys. 87 (1982) 131–152. https://doi.org/10.1007/BF01211061.
[2] J. Benedikt, P. Girg, L. Kotrla, P. Takac, Origin of the $p$ -Laplacian and A. Missbach, Electron. J. Differential Equations 2018 (2018), No. 16, 1–17. http://ejde.math.txstate.edu or http://ejde.math.unt.edu.
[3] C. Bereanu, P. Jebelean, J. Mawhin, Radial solutions for Neumann problems involving mean curvature operators in Euclidean and Minkowski spaces, Math. Nachr. 283 (2010) 379–391. https://doi.org/10.1002/mana.200910083.
[4] D. Bonheure, B. Noris, T. Weth, Increasing radial solutions for Neumann problems without growth restrictions, Ann. Inst. H. Poincaré Anal. Non Linéaire 29 (2012) 573–588. https://doi.org/10.1016/j.anihpc.2012.02.002.
[5] D. Bonheure, E. Serra, P. Tilli, Radial positive solutions of elliptic systems with Neumann boundary conditions, J. Funct. Anal. 265 (2013) 375–398. doi: 10.1016/j.jfa.2013.05.027.
[6] A. Boscaggin, F. Colasuonno, B. Noris, Multiple positive solutions for a class of $p$ -Laplacian Neumann problems without growth conditions, ESAIM Control Optim. Calc. Var. 24 (2018) 1625–1644. https://doi.org/10.1051/cocv/2017074.
[7] A. Boscaggin, G. Feltrin, Pairs of positive radial solutions for a Minkowski-curvature Neumann problem with indefinite weight, Nonlinear Anal. 196 (2020), 111807. https://doi.org/10.1016/j.na.2020.111807.
[8] F. Colasuonno, B. Noris, Radial positive solutions for $p$ -Laplacian supercritical Neumann problems, Bruno Pini Mathematical Analysis Seminar 8(1) (2017) 55–72. https://doi.org/10.6092/issn.2240-2829/7797
[9] K. Deimling, Nonlinear Functional Analysis, Springer, Berlin, 1985.
[10] C. I. Gheorghiu, A third-order nonlinear BVP on the half-line, Chebfun, January 2020. https://www.chebfun.org/examples/ode-nonlin/GulfStream.html.
[11] J. López-Gómez, P. Omari, S. Rivetti, Positive solutions of a one-dimensional indefinite capillarity-type problem: a variational approach, J. Differential Equations 262 (2017) 2335–2392.
[12] R. Ma, T. Chen, H. Wang, Nonconstant radial positive solutions of elliptic systems with Neumann boundary conditions, J. Math. Anal. Appl. 443 (2016) 542–565. https://doi.org/10.1016/j.jmaa.2016.05.038.
[13] D. O’Regan, H. Wang, Positive radial solutions for $p$ -Laplacian systems, Aequationes Math. 75 (2008) 43–50. https://doi.org/10.1007/s00010-007-2909-3.
[14] R. Precup, Existence, localization and multiplicity results for positive radial solutions of semilinear elliptic systems, J. Math. Anal. Appl. 352 (2009) 48–56. doi:10.1016/j.jmaa.2008.01.097.
[15] R. Precup, P. Pucci, C. Varga, Energy–based localization and multiplicity of radially symmetric states for the stationary $p$ -Laplace diffusion, Complex Var. Elliptic Equ. 65 (2020) 1198–1209. https://doi.org/10.1080/17476933.2019.1574774.
[16] R. Precup, J. Rodriguez-Lopez, Positive radial solutions for Dirichlet problems via a Harnack type inequality, submitted.
[17] L. N. Trefethen, Computing numerically with functions instead of numbers, Math. Comput. Sci. 1 (2007) 9–19. https://doi.org/10.1007/s11786-007-0001-y.
[18] L. N. Trefethen, A. Birkisson, T. A. Driscoll, Exploring ODEs, SIAM, Philadelphia, 2018.
[19] H. Wang, On the existence of positive solutions for semilinear elliptic equations in the annulus, J. Differential Equations 109 (1994) 1–7. https://doi.org/10.1006/jdeq.1994.1042.

2024

Recovering and posting started for the older annual reports, 1951-1989

February 8, 2024

Integral equation with maxima via fibre contraction principle

December 6, 2024

A Stancu type extension of the Cheney-Sharma Chlodovsky operators

September 9, 2024

Theory and computation of radial solutions for Neumann problems with φ-Laplacian

Abstract

Authors

Keywords

Paper coordinates

PDF

About this paper

Journal

Publisher Name

DOI

Print ISSN

Online ISSN

References

Paper (preprint) in HTML form

Theory and computation of radial solutions for Neumann problems with ϕitalic-ϕ\phiitalic_ϕ-Laplacian

Abstract.

Key words and phrases:

1991 Mathematics Subject Classification:

1. Introduction

2. Solution properties

2.1. Some basic properties

Lemma 2.1.

Proof.

Lemma 2.2.

Proof.

Lemma 2.3.

Proof.

2.2. A Harnack type inequality

Theorem 2.4.

3. Existence, localization and multiplicity

3.1. The solution operator

Lemma 3.1.

Proof.

3.2. Existence and localization

Lemma 3.2.

Lemma 3.3.

Lemma 3.4.

Proof.

Theorem 3.5.

Proof.

Remark 3.1.

3.3. Decreasing solutions

Theorem 3.6.

Remark 3.2.

3.4. Multiple solutions

Theorem 3.7.

Proof.

Theorem 3.8.

Proof.

Remark 3.3.

3.5. Existence and multiplicity under asymptotic conditions

4. Numerical solutions

4.1. First example

4.2. Second example

4.3. Third example

References

Related Posts

Theory and computation of radial solutions for Neumann problems with $\phi$ -Laplacian