Monotone wave fronts for $(p, q)$-Laplacian driven reaction-diffusion equations

TL;DR

This paper investigates the existence and properties of monotone traveling wave solutions for a reaction-diffusion equation involving a $(p, q)$-Laplacian operator, providing estimates of critical speeds and numerical insights.

Contribution

It introduces new existence results for monotone heteroclinic waves in $(p, q)$-Laplacian driven equations and analyzes the influence of parameters $p$ and $q$ on wave dynamics.

Findings

Derived estimates for critical wave speeds.

Analyzed the impact of $p$ and $q$ on wave behavior.

Provided numerical simulations illustrating theoretical results.

Abstract

We study the existence of monotone heteroclinic traveling waves for the $1$-dimensional reaction-diffusion equation $$ u_t = (| u_x |^{p-2} u_x + | u_x |^{q-2} u_x)_x + f(u), $$ where the non-homogeneous operator appearing on the right-hand side is known as $(p, q)$-Laplacian. Here we assume that $2 \leq q < p$ and $f$ is a nonlinearity of Fisher type, namely it is always positive out of its zeros. We give an estimate of the critical speed and we comment on the roles of $p$ and $q$ in the dynamics, providing some numerical simulations.

Equations

$$u_t=(|u_x{|}^{p-2}{u}_x+|u_x{|}^{q-2}{u}_x{)}_x+f(u),$$

$$\left\{\begin{array}[]{l}(|v^{\prime}|^{p-2}v^{\prime}+|v^{\prime}|^{q-2}v^{\prime})^{\prime}-cv^{\prime}+f(v)=0,\qquad(v=v(z))\\ v(-\infty)=0,\;v(+\infty)=1.\end{array}\right.$$

$$(|v^{\prime }{|}^{q-2}{v}^{\prime }-|v^{\prime }{|}^{p-2}{v}^{\prime }{)}^{\prime }-c{v}^{\prime }+f(v)=0.$$

$$u_t=(|u_x{|}^{p-2}{u}_x+|u_x{|}^{q-2}{u}_x)+f(u),$$

$$\left\{\begin{array}[]{l}(|v^{\prime}|^{p-2}v^{\prime}+|v^{\prime}|^{q-2}v^{\prime})^{\prime}-cv^{\prime}+f(v)=0\\ v(-\infty)=0,\;v(+\infty)=1,\;v^{\prime}>0.\end{array}\right.$$

$$\frac{d}{dv}Q\left(\varphi (v)\right)-c\varphi (v)+f(v)=0,$$

$$Q(s)=\frac{p-1}{p}|s{|}^p+\frac{q-1}{q}|s{|}^q$$

$$\left\{\begin{array}[]{l}\displaystyle y^{\prime}=cR(y)-f(v),\\ \\ y(0)=0=y(1),\;y(v)>0\textrm{ for }v\in\,]0,1[\,.\end{array}\right.$$

$$\underset{s\to 0^{+}}{\lim }\frac{Q(s)}{s^q}=\frac{q-1}{q}\mathrm{and}\underset{s\to \infty }{\lim }\frac{Q(s)}{s^p}=\frac{p-1}{p},$$

$$\underset{s\to 0^{+}}{\lim }\frac{R(s)}{s^{1/q}}={\left(\frac{q}{q-1}\right)}^{\frac{1}{q}}\mathrm{and}\underset{s\to \infty }{\lim }\frac{R(s)}{s^{1/p}}={\left(\frac{p}{p-1}\right)}^{\frac{1}{p}}.$$

$$\underset{s\to 0^{+}}{\lim }\frac{f(s)}{s^{q^{\prime }-1}}=L_0\mathrm{and}\underset{s\in ]0,1]}{\sup }\frac{f(s)}{s^{q^{\prime }-1}}=L_{+}<+\infty$$

$$f(s)\le k(1-s),\text{ for every }s\in [0,1].$$

$$L_0^{1/q^{\prime }}{q}^{\prime 1/q^{\prime }}{q}^{1/q}\le c^{\ast }\le c_{+},$$

$$c_{+}=\left\{\begin{array}[]{ll}\displaystyle\frac{L_{+}^{\frac{1}{q^{\prime}}}\,q}{q-1}\cdot(p+q-2)^{\frac{1}{q}}&\textrm{ if }\,L_{+}\leq p+q-2\vspace{0.3cm}\\ \displaystyle\frac{p(p+q-2)}{q-1}&\textrm{ if }\,\displaystyle p+q-2<L_{+}\leq\frac{p-1}{q-1}(p+q-2)\vspace{0.3cm}\\ \displaystyle\frac{L_{+}^{\frac{1}{p^{\prime}}}\,p}{(q-1)^{\frac{1}{p}}(p-1)^{\frac{p-1}{p}}}\;(p+q-2)^{\frac{1}{p}}&\textrm{ if }\,\displaystyle L_{+}>\frac{p-1}{q-1}(p+q-2),\end{array}\right.$$

$$\frac{du}{dt}=R(y(u(t))),u(0)=1/2;$$

$$\underset{s\to 0^{+}}{\lim }\frac{f(s)}{s^{q^{\prime }-1}}=L_0,$$

$$c\ge L_0^{1/q^{\prime }}{q}^{\prime 1/q^{\prime }}{q}^{1/q}.$$

$$l=\underset{v\to 0^{+}}{\mathrm{limsup}}\frac{y(v)}{v^{q^{\prime }}}\le c\underset{v\to 0^{+}}{\mathrm{limsup}}\frac{R(y(v))}{q^{\prime }{v}^{q^{\prime }-1}}-\frac{L_0}{q^{\prime }},$$

$$l-\frac{c}{q^{\prime 1/q^{\prime }}}{l}^{1/q}+\frac{L_0}{q^{\prime }}\le 0.$$

$$\frac{L_0}{q^{\prime }}\le \underset{x\in \mathbb{R}}{\max }\frac{c}{q^{\prime 1/q^{\prime }}}{x}^{1/q}-x;$$

$$\underset{s\in ]0,1]}{\sup }\frac{f(s)}{s^{q^{\prime }-1}}=L_{+}<+\infty .$$

$$\left\{\begin{array}[]{l}\displaystyle y^{\prime}=cR(y)-f(v),\\ \\ y(1)=0,\end{array}\right.$$

$$\left\{\begin{array}[]{l}\displaystyle y^{\prime}=cR(y)-f(v),\\ \\ y(0)=0,\end{array}\right.$$

$$y(u)=Q(\beta u^{\alpha }),u\in [0,1],$$

$$\alpha \left[(p-1){\beta }^p{u}^{\alpha p-1}+(q-1){\beta }^q{u}^{\alpha q-1}\right]-c\beta u^{\alpha }+f(u)\le 0.$$

$$(q^{\prime }-1)\left[(p-1){\beta }^p{u}^{p(q^{\prime }-1)-1}+(q-1){\beta }^q{u}^{q(q^{\prime }-1)-1}\right]-c\beta u^{q^{\prime }-1}+f(u)\le 0,$$

$$(q^{\prime }-1)\left[(p-1){\beta }^p{u}^{p(q^{\prime }-1)-q^{\prime }}+(q-1){\beta }^q\right]-c\beta +L_{+}\le 0.$$

$${\mathcal{G}}_c(\beta ):=(q^{\prime }-1)\left[(p-1){\beta }^p+(q-1){\beta }^q\right]-c\beta +L_{+}=\frac{p-1}{q-1}{\beta }^p+{\beta }^q-c\beta +L_{+}\le 0.$$

$$\bar{c}=\frac{L_{+}^{\frac{q-1}{q}}q}{q-1}(p+q-2{)}^{\frac{1}{q}}$$

Full text

[TABLE]

where the non-homogeneous operator appearing on the right-hand side is known as $(p,q)$-Laplacian. Here we assume that $2\leq q<p$ and $f$ is a nonlinearity of Fisher type, namely it is always positive out of its zeros. We give an estimate of the critical speed and we comment on the roles of $p$ and $q$ in the dynamics, providing some numerical simulations.

††AMS Subject Classification: 35K55, 35K57, 35J92, 34C37.††Keywords: reaction-diffusion equations, $(p,q)$-Laplacian, admissible speeds.

1 Introduction

This paper deals with wave fronts for the $1$-dimensional reaction-diffusion equation

[TABLE]

where $2\leq q<p$ and $f$ is a regular function such that $f(0)=0=f(1)$. The operator driving the diffusion is a $(p,q)$-Laplacian type one, and it is characterized by the presence of two (cooperative) terms, the former prevailing for large gradients and the latter prevailing for small ones.

In these last years, there has been an increasing interest for problems driven by this operator, especially in the stationary (elliptic) case, both in bounded and unbounded domains (see, e.g., [2, 15] and the references therein). A $(p,q)$-Laplacian type operator may arise quite naturally, for instance, when taking into account more terms in the Taylor expansion of a quasi-linear operator, and finds application, e.g, in modeling chemical reactions and some physical phenomena (see for instance [3, 5, 13]). Of course, equations like (1.1) have been widely studied in the case $p=q$, i.e. for the $p$-Laplacian operator (we limit ourselves to mention some of the papers which are strictly related to the topic of the present investigation, namely [1, 6, 7, 8, 12]). The main difficulty in the case $p\neq q$ comes of course from the lack of homogeneity of the operator, which turns into preventing a straight application of almost any usual argument and requires some supplementary care.

By the expression wave front, we mean a heteroclinic solution of the type $u(t,x)=v(x+ct)$ connecting two equilibria $u_{-},u_{+}$ (namely, such that $v(-\infty)=u_{-}$, $v(+\infty)=u_{+}$). In view of the assumptions on $f$, for (1.1) such profiles will connect the two equilibria $u_{-}=0$ and $u_{+}=1$. Dealing with reaction-diffusion equations, the search for these solutions is the first natural step to get a better insight into the dynamics of the model. They play the role of phase transitions between two steady states of the system and, in the context of the original genetics models proposed by Fisher [9, 14], they actually were meant to describe the relative propagation of a favourable gene inside the considered population. For this reason, it is very natural to assume that $0\leq u\leq 1$.

Writing the equation for the traveling fronts by setting $u(t,x)=v(x+ct)$, equation (1.1) leads to the problem

[TABLE]

We will limit ourselves to the case when $f$ is strictly positive out of the equilibria (a Fisher-type case). To better compare with the existing literature, moreover, we will search since the very beginning for monotone fronts, by-passing some problems which may arise from the fact that uniqueness is not ensured for the differential equation in (1.2).

The first question we are interested in regards the possible speeds of propagation of the traveling waves, namely the values $c\in\mathbb{R}$ for which (1.2) has a monotone solution, which are called admissible speeds. In the classical case $p=q=2$, this problem has been extensively studied (see, e.g., the references in [4, 10]) and it has been shown that the set of the admissible speeds is an unbounded interval, whose lower endpoint is called critical speed and is denoted by $c^{*}$. The value of $c^{*}$ depends on the behavior of the reaction term $f$: if $f$ is always below its tangent line in [math], it is $c^{*}=2\sqrt{f^{\prime}(0)}$. The existence of a critical speed $c^{*}$ has been extended to the general case $p=q$ in [6, 7, 12], showing moreover that $c^{*}$ admits a variational characterization like in the linear case.

The purpose of the present investigation is to analyze when and how the estimates for the critical speed change in presence of the two different powers $p$ and $q$ in the expression of the operator. Indeed, the regimes acting for $|z|$ large involve very small gradients, since for a monotone heteroclinic solution it is $v^{\prime}\to 0$ at infinity. Hence, the dominant term for $z=x+ct$ large will be the one with exponent $q$. However, no information is a priori available on the derivative $v^{\prime}$, so that it may be that there are regions where the term with the power $p$ is prevailing. Intuitively, such effect should be more accentuated the more there is room for the derivative to grow (and this is well seen in Figure 3 in Section 3). Indeed, it turns out that the value of the constant providing a bound for the critical speed is actually influenced by this interaction when the nonlinearity is sufficiently large, see Theorem 2.1 below. On the other hand, in the case $p=q$, we recover the results proved in [6, 7] for the $p$-Laplacian operator.

After our main statement, in Section 2.2 we give some complementary remarks mentioning, in particular, the case when the two diffusive terms are in competition, namely considering the equation

[TABLE]

Avoiding a degeneration of the problem which may appear if the derivative reaches the value $1$, we also briefly discuss the picture for fronts in this situation (see Remark 2.3).

The technique exploited makes use of a change of variables allowing to lead the original problem back to a two-point first order problem on $[0,1]$, to be dealt with through a shooting technique.

The paper ends with Section 3, devoted to some numerical simulations illustrating our main results.

2 Main result and discussion

In this section, we study the admissible speeds for the $1$-dimensional partial differential equation

[TABLE]

assuming that $2\leq q<p$ and that the reaction term $f(u)$ satisfies the following hypothesis:

(f)  $f:[0,1]\to\mathbb{R}$ is a $C^{1}$-function such that $f(0)=0=f(1)$ and $f(u)>0$ for $u\in\,]0,1[\,$.

2.1 Statement and proof of the main result

As discussed above, the search for monotone heteroclinic traveling waves $u(t,x)=v(x+ct)$ solving (2.1) leads us to study the second order boundary value problem

[TABLE]

Since $v$ is strictly monotone, the map $z\mapsto v(z)$ is invertible, so that we can take $v$ as the new independent variable and write $z=z(v)$. Setting $\phi(v)=v^{\prime}(z(v))$ and proceeding similarly, e.g., as in [7, 10], the differential equation in (2.2) is rewritten as

[TABLE]

being $Q(s)$ the primitive of $(p-1)|s|^{p-2}s+(q-1)|s|^{q-2}s$ that satisfies $Q(0)=0$. In particular, it turns out that

[TABLE]

and $Q(\cdot)$ is invertible for $s\geq 0$; it is worth noticing that we can actually invert $Q(\phi(v))$ since we search for increasing fronts, such that $\phi(v)>0$. Denoting by $R(\cdot)$ the (functional) inverse of $Q(\cdot)$ and setting $y(v)=Q(\phi(v))$, problem (2.2) is thus rephrased as

[TABLE]

Moreover, since

[TABLE]

it holds

[TABLE]

We now state our main result estimating the critical speed for problem (2.2); henceforth, for a number $r>1$, we will denote by $r^{\prime}$ its conjugate exponent, namely the positive number such that $1/r+1/r^{\prime}=1$.

Theorem 2.1**.**

Let $f:[0,1]\to\mathbb{R}$ satisfy assumption (f). Moreover, let

[TABLE]

and assume that there exists $k>0$ such that

[TABLE]

Then, there exists $c^{*}$ satisfying

[TABLE]

with $c_{+}$ given by

[TABLE]

such that every $c\geq c^{*}$ is admissible for problem (2.2).

We first observe that assumption (2.5) is not used to solve problem (2.3) but ensures that the solutions of (2.3) are in a one-to-one correspondence with global strictly monotone profiles, defined on the whole $\mathbb{R}$ and solving (2.2). Such profiles are reconstructed, up to translations, by solving the problem

[TABLE]

similarly as in [6, Proposition 2.3], in view of (2.4) we have that the solution of $y^{\prime}=cR(y)-f(u)$ has a growth controlled by $u^{q^{\prime}}$ near [math], while (2.5) ensures that, for the same solution, it holds $y(u)\leq C(1-u)^{2}$. It follows that the solution of (2.8) is defined (and strictly monotone) on the whole $\mathbb{R}$ (see [6, Proposition 2.3] for further details).

In order to prove Theorem 2.1, we first give a lower bound on the critical speed $c^{*}$. This will only come from the behavior of $R(y)$ near [math], which, in view of (2.4), will be ruled completely by the smallest power inside the operator. It is indeed sufficient to assume that

[TABLE]

for a positive constant $L_{0}$, in order to state the following result, which follows similarly as in [6].

Lemma 2.1**.**

Assume (2.9) and let $y(v)$ be a solution of the differential equation in (2.3), satisfying $y(0)=0$ and $y(v)>0$ in a right neighborhood of [math]. Then,

[TABLE]

Proof.

As usual, we want to find an equation for $l=\limsup_{v\to 0^{+}}\frac{y(v)}{v^{q^{\prime}}}\geq 0$. From the differential equation in (2.3), dividing by $q^{\prime}v^{q^{\prime}-1}$ and using De l’Hopital rule, we obtain

[TABLE]

from which, recalling (2.4), it holds

[TABLE]

For this inequality to be true, at least it has to be

[TABLE]

the maximum at the right-hand side is attained at $x=c^{q^{\prime}}/(q^{\prime}q^{q^{\prime}})$ and a direct computation shows that (2.11) can be true only if $c$ satisfies (2.10). ∎

Our second result provides the upper bound for the admissible speeds. As usual in performing this estimate, a global bound on $f$ is needed. We thus assume that

[TABLE]

The strategy is the following shooting-type one. First, we observe that, for every $c$, $y\equiv 0$ is a subsolution to the backward Cauchy problem

[TABLE]

since $R(0)=0$ and $f$ has positive sign. Hence, the solution $y(v)$ of (2.13) (uniqueness holds for (2.13) since $R$ is increasing) is nonnegative. Even more, $y(v)$ never vanishes in $\,]0,1[\,$, due to the sign of $f$, so it is well defined up to $v=0$. If $y(v)=0$, then we are done; otherwise, the idea is to search for suitable values of $c$ for which there exists a positive subsolution for the forward Cauchy problem

[TABLE]

defined on the whole $[0,1]$. In this way, a standard uniqueness argument (out of $y=0$) will ensure that, for these values of $c$, we can find the desired connection between [math] and $1$.

We have the following result.

Proposition 2.1**.**

Let $f$ satisfy assumption (2.12). Then, for every $c\geq c_{+}$, where $c_{+}$ is given by (2.7), there exists a positive subsolution of (2.14).

Proof.

We will search the subsolution of (2.14) in the form

[TABLE]

with $\alpha>0$ and $\beta>0$ to be determined. By substituting into the differential equation, noticing that $y(u)>0$ for every $u\in\,]0,1]$ so that $Q(y(u))$ is always invertible in $[0,1]$, in order for $y(u)$ to be a subsolution it has to be

[TABLE]

If we now divide by $u^{\alpha}$, it is easy to check that the smallest positive $\alpha$ such that both $\alpha p-1-\alpha$ and $\alpha q-1-\alpha$ are positive is $\alpha=q^{\prime}-1$. Hence we set $\alpha=q^{\prime}-1$ and (2.16) reads as

[TABLE]

which, dividing by $u^{q^{\prime}-1}$ and using (2.12), becomes

[TABLE]

Since (2.17) has to be satisfied for every $u\in[0,1]$, we are thus led to investigate for which values of $c$ there exists $\beta>0$ such that

[TABLE]

Of course, if we are able to find $\bar{c}$ such that $\mathcal{G}_{\bar{c}}(\bar{\beta})\leq 0$ for a suitable $\bar{\beta}>0$, then the fact that $\mathcal{G}_{c}(\beta)$ is decreasing in the variable $c$ ensures that for every $c>\bar{c}$ it will be $\mathcal{G}_{c}(\bar{\beta})\leq 0$, yielding the statement. We now briefly sketch the proof in the three cases highlighted in (2.7), which we refer to as (2.7)(i), (2.7)(ii) and (2.7)(iii).

Case (2.7)(i): $L_{+}\leq p+q-2$. In this case, we fix

[TABLE]

and show that $\min\mathcal{G}_{\bar{c}}\leq 0$. To this end, we consider the function

[TABLE]

which attains its minimum at the point

[TABLE]

In view of the assumptions on $L_{+}$ and $\bar{c}$, we have that

[TABLE]

On the other hand, since $\bar{\beta}\leq 1$, it is also $\mathcal{G}_{\bar{c}}(\bar{\beta})\leq g_{\bar{c}}(\bar{\beta})=0,$ showing that the choice $c=\bar{c}$ allows to construct a positive lower solution of the form (2.15) for $\beta=\bar{\beta}$.

• 2)

Case (2.7)(ii): $p+q-2<L_{+}\leq\frac{p-1}{q-1}(p+q-2)$. Here we set

[TABLE]

and we compute the minimum of

[TABLE]

which is attained at

[TABLE]

In view of the position done, we have $h_{\bar{c}}(\bar{\beta})\leq 0$, where the equality holds if and only if $L_{+}=(p-1)(p+q-2)/(q-1)$. Since $\bar{\beta}=1$, we have that

[TABLE]

yielding the conclusion.

• 3)

Case (2.7)(iii): $L_{+}\geq\frac{p-1}{q-1}(p+q-2)$. The picture is here similar to the previous item, apart from the fact that, fixing

[TABLE]

it turns out that $\bar{\beta}$ defined in (2.19) is greater than one for $L_{+}>\frac{p-1}{q-1}(p+q-2)$, and $\mathcal{G}_{\bar{c}}(\bar{\beta})\leq h_{\bar{c}}(\bar{\beta})=0$.

∎

Theorem 2.1 now follows combining the two estimates given in Lemma 2.1 and Proposition 2.1.

Let us observe that the estimate provided in Proposition 2.1 may be not optimal: a gap in the optimality is for sure created whenever the sign of $\mathcal{G}_{c}$ is different from the one of $g_{c}$ (or $h_{c}$) and the minimum of $\mathcal{G}_{c}$ is near $1$. A numerical search for the exact value of $\beta$ would provide better bounds, but since we are not able to solve inequality (2.18) without passing through a stronger control involving only one power of $\beta$, we have preferred to state the result with an analytical proof. We postpone to the next section some remarks and corollaries of Theorem 2.1.

2.2 Some remarks and corollaries

We begin the section by observing that Proposition 2.1 is constructed through a sufficient condition for $c$ to be admissible (implying the existence of a positive lower solution), so that it may be that the provided bounds can be further improved with another strategy of proof. However, in the expression of $c_{+}$ given in Theorem 2.1, it appears a (small) range of values for $L_{+}$, given in (2.7)(ii), where the process is really undergoing the interaction of the two powers $p$ and $q$ at the same level (while we see that, in cases (2.7)(i) and (2.7)(iii), the bounds for $c^{*}$ involve in a very marginal way $p$ and $q$, respectively). This range actually appears to exist and is not a simple drawback of our estimates: in Figure 1 we plot the graph of the function $\mathcal{G}_{c}(\beta)$ for $p=4$, $q=3$, $L_{+}=6$ (such that $p+q-2<L_{+}\leq\frac{p-1}{q-1}(p+q-2)$) and $c$ given by (2.7)(i), and we see that in this case $\mathcal{G}_{c}$ is always positive, meaning that the bound for $c_{+}$ given by (2.7)(i) has to be corrected.

In fact, when $p=q$, the range in (2.7)(ii) disappears, while the expressions of $c_{+}$ in (2.7)(i) and (2.7)(iii) coincide. Observing that, in this case, $Q(s)=2\tfrac{q-1}{q}|s|^{q}$, and thus (2.4) becomes

[TABLE]

with the same proofs as in the previous section we are able to recover the results proved for the $p$-Laplacian, up to a constant which comes from the fact that we are considering the operator given by $(2|u^{\prime}|^{p-2}u^{\prime})^{\prime}$. The following result can be compared with [6, Proposition 5.5].

Corollary 2.1**.**

Let $f$ be as in Theorem 2.1 and $p=q$. Then, there exists $c^{*}$ with

[TABLE]

such that every $c\geq c^{*}$ is admissible for problem (2.2).

Incidentally, notice that in the particular case when $L_{0}=L_{+}$ this statement actually gives the exact value of the critical speed.

We continue with a series of remarks.

Remark 2.1**.**

In the special case $p=2q$, it is possible to perform some more explicit computations. Indeed, since the expression for $Q(y)$ reads as

[TABLE]

a straightforward computation leads to

[TABLE]

We notice that, on one hand, $R(y)$ can be rewritten as

[TABLE]

showing that, for $y\to 0$, the behavior of $R(y)$ is dictated by $y^{1/q}$; on the other hand,

[TABLE]

implying that $R(y)$ goes to infinity like $y^{1/2q}$ for $y\to\infty$. Hence, (2.4) is satisfied, as expected.

Remark 2.2**.**

Theorem 2.1 is extendable to the case of density-dependent diffusions, namely for an equation like

[TABLE]

where $D(s)$ is a function with strictly positive derivative. The details of the change of variables in this case can be found, for instance, in [11, Section 2]. For the sake of briefness, we avoid giving the explicit statement, which is actually a natural adaptation of the one above. Furthermore, also the case with convection is likely to be dealt with using a similar strategy as in [11].

Remark 2.3**.**

We quote another case which we find worth briefly showing, that is, the situation when the considered equation reads as

[TABLE]

namely the two terms appearing in the expression of the operator driving the diffusion are in competition. Clearly, here $p$ has to be different from $q$ and, as before, we assume $p>q$ in order for the dynamics to be comparable with the one already examined at least at infinity: namely, we want $Q(\phi(v))$ to be positive at least in the region where $\phi(v)$ is small, that is, $v\approx 0$ and $v\approx 1$. Here the operator degenerates for $|u_{x}|=1$ and a problem of invertibility arises for the function $Q(s)=\frac{q-1}{q}|s|^{q}-\frac{p-1}{p}|s|^{p}$; nevertheless, if we are able to restrict the dynamics in the neighborhood $\mathcal{I}$ of [math] where $Q$ is invertible, then we can still recover the existence of fronts. As it is expected because of the previous discussion, in this case the critical speed $c^{*}$ will depend only on the power $q$; indeed, we can prove the following result.

Proposition 2.2**.**

Let $f$ be as in Theorem 2.1; then there exists $c^{*}$, satisfying

[TABLE]

such that every

[TABLE]

is admissible (provided that such an interval is non-empty).

Of course, the non-emptiness of the interval appearing in (2.20) turns out to impose a condition on $L_{+}$, which however appears cumbersome and not particularly significant (and for this reason we omit it). Observe also that the statement yields again the exact value of the “critical speed” whenever $L_{0}=L_{+}$.

Proof.

Let us begin the proof by observing that the neighborhood $\mathcal{I}$ of [math] where $Q$ is invertible is given by

[TABLE]

Hence, in searching a subsolution of the form $Q(\beta u^{\alpha})$, we have to ensure that $\beta u^{\alpha}\leq s_{0}$ for every $u\in[0,1]$, namely $\beta\leq s_{0}$.

On one hand, the lower bound (2.10) for the critical speed still holds without any change. On the other hand, when trying to prove an upper bound analogous to the one given in (2.6)-(2.7), by using the same strategy and assuming that $Q$ is always invertible along $\beta u^{\alpha}$ ($u\in[0,1]$), we end up with

[TABLE]

Since in this case we have a negative sign in front of the term coming from the $p$-Laplacian operator and since $\inf_{u\in[0,1]}u^{p(q^{\prime}-1)-q^{\prime}}=0,$ it turns out that it is sufficient to find a value $\bar{\beta}$ such that $\bar{\beta}^{q}-c\bar{\beta}+L_{+}\leq 0.$ Exploiting the same computations as in the proof of Proposition 2.1 leads to

[TABLE]

while the upper bound in (2.20) is needed in order for $\bar{\beta}$ to be less than $s_{0}$ (and thus for the expression (2.21) to be allowed). ∎

3 Some numerical simulations

In this section, we collect a series of numerical simulations in order to better understand the role of the $p$-Laplacian in the dynamics of the traveling waves.

We have seen in Proposition 2.1 that, in spite of the fact that the operator ruling the dynamics at infinity is the $q$-Laplacian one (since the monotone wave profiles at infinity possess very small derivative), the $p$-Laplacian seems to have some influence on the value of the critical speed, especially if the reaction term is sufficiently large. Actually, this seems confirmed also by the numerical simulations we are going to present and is better understood by recalling a simple invariance-upon-rescaling property of the arguments presented in Section 2. Namely, if we assume that $f:[0,H]\to\mathbb{R}$ is such that $f(0)=0=f(H)$ and $f>0$ on $\,]0,H[\,$, for $H>0$ fixed, and we search for monotone connections between [math] and $H$, we are led to solve

[TABLE]

Rescaling the independent variable by setting $v=u/H\in[0,1]$, we can bring back (3.1) to a problem on $[0,1]$: indeed, $y_{c}(u)$ is a solution of (3.1) if and only if $w(v)=y_{c}(Hv)$ solves

[TABLE]

In the particular case of a $q$-Laplacian equation, leading to $Q(s)=\tfrac{q-1}{q}|s|^{q}$, so that $R(s)=(\tfrac{q}{q-1}|z|)^{1/q}$, the assumption on the reaction term guaranteeing the exact computation of the critical speed reads as

[TABLE]

which turns out into asking that, setting $g(v)=Hf(Hv)$, it holds

[TABLE]

However, from the computation of the critical speed done in [6] (see also Corollary 2.1) in this particular case, we have that (3.2) has a solution if and only if

[TABLE]

but this is equivalent to

[TABLE]

namely it is possible to compute the critical speed for (3.1) with the same formula as for (3.2).

Incidentally, we notice that assuming instead $f(u)\leq f^{\prime}(0)u$ for every $u\in[0,H]$, the same argument allows to infer that the critical speed for problem (3.1) with $R(s)=\sqrt{2s}$ (namely, for the classical Fisher equation) is always given by $c^{*}=2\sqrt{f^{\prime}(0)}$.

This argument fails for our equation, as we expect also in view of the discussion at the beginning of Section 2.2. We show indeed two numerical simulations (Figures 2 and 3) with a completely different behavior: in the first one, for a small $H$, we show that the critical speed appears to be the one associated with the $q$-Laplacian, given by the bound (3.3); the second one, for a larger $H$, shows that the critical speed is strictly larger than the one provided by such bound. Again, we remark that this seems in line with the intuitive observation that, the more there is room for $v^{\prime}$ to possibly become large (i.e., the larger is $H$), the more the $p$-Laplacian passes to govern the diffusion process. Indeed, in Figure 2 we have that $y(u)$ is always very small, while this is not the case for Figure 3. We refer to the captions for further comments.

The third simulation we show is instead devoted to a competitive case, namely to the situation described in Remark 2.3, where the operator is given by the difference between a $q$-Laplacian and a $p$-Laplacian. As already remarked, in this case a problem of invertibility of the function $Q(s)$ arises, and we can search for monotone connections only as long as we remain in the region where its inverse $R(\cdot)$ is well-defined. In Figure 4, we show a case where the interval appearing in (2.20) is not empty and the critical speed is computed exactly; in the situation depicted in Figure 5, instead, there are no admissible speeds (again we refer the reader to the captions for further comments).