Next Article in Journal
Order-of-Addition Orthogonal Arrays with High Strength
Next Article in Special Issue
Solitons in Neurosciences by the Laplace–Adomian Decomposition Scheme
Previous Article in Journal
3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects
Previous Article in Special Issue
Numerical Solutions to Wave Propagation and Heat Transfer Non-Linear PDEs by Using a Meshless Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regularity, Asymptotic Solutions and Travelling Waves Analysis in a Porous Medium System to Model the Interaction between Invasive and Invaded Species

by
José Luis Díaz Palencia
1,*,
Julián Roa González
2,*,
Saeed Ur Rahman
3 and
Antonio Naranjo Redondo
4
1
Escuela Politécnica Superior, Universidad Francisco de Vitoria, 28223 Madrid, Spain
2
Department of Education, Universidad a Distancia de Madrid, 28400 Madrid, Spain
3
Department of Mathematics, COMSATS University Islamabad, Abbottabad 22060, Pakistan
4
Technology Programs, Schiller International University, 28002 Madrid, Spain
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(7), 1186; https://doi.org/10.3390/math10071186
Submission received: 28 February 2022 / Revised: 29 March 2022 / Accepted: 31 March 2022 / Published: 5 April 2022
(This article belongs to the Special Issue Partial Differential Equations: Theory and Applications)

Abstract

:
This work provides an analytical approach to characterize and determine solutions to a porous medium system of equations with views in applications to invasive-invaded biological dynamics. Firstly, the existence and uniqueness of solutions are proved. Afterwards, profiles of solutions are obtained making use of the self-similar structure that permits showing the existence of a diffusive front. The solutions are then studied within the Travelling Waves (TW) domain showing the existence of potential and exponential profiles in the stable connection that converges to the stationary solutions in which the invasive species predominates. The TW profiles are shown to exist based on the geometry perturbation theory together with an analytical-topological argument in the phase plane. The finding of an exponential decaying rate (related with the advection and diffusion parameters) in the invaded species TW is not trivial in the nonlinear diffusion case and reflects the existence of a TW trajectory governed by the invaded species runaway (in the direction of the advection) and the diffusion (acting in a finite speed front or support).

1. Introduction

According to the Convention of Biological Diversity (p. 1, Ch. 1 [1]), the biological invasion is defined as “…those alien species which threaten ecosystems, habitats or species”. The problem analyzed can be understood within the perspective of an invasive species that invades a region or space previously inhabited by the invaded species.
The interaction between invaded-invasive species has been analyzed based on a nonlinear diffusion previously in [2] to analyze biological cell dynamics related to a haptotactic process in cancer. Furthermore, stability analysis for an haptotaxis dynamics in melanoma invasion is provided in [3].
The reaction–diffusion problems with advection have been a source of research to model different physical scenarios. In [4], the author analyzes the impact of an advection coefficient to understand the hydrodynamical processes in combustion. In this case, it is particularly relevant to accurately model the advection coefficient to couple different involved mechanisms: chemical kinetic, thermodynamics, and fluid dynamics. A general nonlinear time-dependent advection does not necessarily lead to purely monotonic solutions; nonetheless, in [5], the authors analyze monotone properties of bounded Travelling Waves (TW) solutions for a density dependent diffusion and a nonlinear convection.
There are some analogies between predator–prey and invasive-invaded species models. Indeed, when an invasive species reaches into a region (a living organ for example), the invaded one may extinguish. Consequently, the invasive species may vanish as there is not invaded species to interact with. This described dynamics can be regarded as a connection with predator–prey processes as the predator vanishes whenever the prey disappears. Note that predator–prey models have been analyzed recently with analytical and stability methods. In [6], different forms of functional responses are provided for modelling the predator–prey dynamic. The cited work considers the harvesting in the predator and the density-dependent mortality in the prey to assess the Hopf bifurcation in the proximity of the equilibrium point. In addition, the Hopf bifurcation method has been followed in [7] to study a delayed density-dependent predator–prey system with Beddington–DeAngelis functional response. The periodic dynamic in the solutions has been shown to exist as a consequence of stochastic disturbance for a Holling II functional response [8]. Similar stability and bifurcation methods have been used in [9] for a predator–prey model subject to the Allee effect with a discrete-time Holling type-IV functional response. The mentioned methods have been used as well in [10] to model an infected predator that consumes the prey according to a Holling type-II response.
There exist some commonalities among Evan functions, Hopf bifurcation and Geometric Perturbation Theory. Indeed, the theories mentioned focus on the searching of solutions stability and to qualitatively describe their behaviour. Solutions may be expressed as TW profiles to this end. Note that the geometric perturbation theory, as a baseline support in the search of TW profiles, is considered in this work.
The nonlinear diffusion considered in this work for modelling purposes is particularly relevant, as it drives the mathematical approaches considered. Somehow, our work is influenced by the seminal model proposed by Keller and Segel [11] to study the chemotaxis walks in cells:
u t = · ( d ( u ) u χ ( v ) u u ) x Ω , t > 0 v t = d v Δ v u v x Ω , t > 0 ,
where v is the chemical agent concentration and u is the cell density. The term d ( u ) represents the diffusivity and χ ( v ) is the distribution of chemical agents sensitive to cells. In several novel articles, the Keller and Segel model was extended to include different forms of reaction and absorption terms (refer to [12,13,14,15]). The problem discussed in [16] was proposed as an integro-differential system of equations with pressure effects to predict the behaviour of cancerous cells, as an invasive species, spreading with a nonlinear diffusion over healthy cells in a closed organ. In [17], the authors provide a complete analysis of a chemotaxis process with a nonlocal reaction source where a uniform bound of solutions are explored for different conditions in the involved coefficients.
The use of nonlinear diffusion is an oblique topic and has been used in other interesting applications where numerical and purely analytical approaches have been followed. Such nonlinear diffusion allows for accurately modelling physical phenomena in which porosity is a governing parameter. The Darcy law involving nanofluids has been considered to derive numerical solutions by the Successive Local Linearization Method [18]. The authors show the different solution profiles obtained after the numerical exercise where the exponential decay can be perceived. Such an exponential decay is immediate in the linear diffusion case; nonetheless, for the nonlinear diffusion, further analytical assessments shall be done. The effect of porosity in a partial slip for a peristaltic transport in a Jeffrey fluid has been investigated in [19]. In addition, Ref. [20] employs a nonlinear diffusion to improve the accuracy in simulating potential coagulation in an electromagnetic blood flow in annular vessel geometries. In all the cited cases, the finite speed of propagation defines a diffusive front. This is a common behaviour to all porous medium equations and shall be characterized for the problem discussed in the present work.
Other models related with perturbation under nonlocal sources and attractive and repulsion actions can be mentioned as precluding studies. In [21], the authors provide an analysis about chemotactic phenomena driven by logistic and repulsive actions. Furthermore, in [22], the authors study a porous medium equation with different reaction–absorption terms, depending in some cases on the gradient of the solution over a bounded domain in R N , N 1 .

2. Model Description and Methods

The problem P analyzed is as follows:
u t = Δ u m + c · u + v ( 1 u ) , v t = Δ v m + c · v u v , m > 1 , u 0 ( x ) , v 0 ( x ) > 0 L 1 ( R d ) L ( R d ) ,
where d > 1 or d = 1 for convenience and is later specified. The solving methods exposed in this work aim to obtain explicit analytical solutions in the form of maximal, minimal solutions and Travelling Waves. The spatial operator involves a nonlinear diffusion with advection and is considered together with nonlinear coupled reaction–absorption terms.
The problem P justification is based on studying the invasive-invaded dynamic and the mathematical properties introduced by the spatial operator. For this purpose:
Consider that the invasive species acts in accordance with the presence of the invaded species. If the quantity of invaded is high, the temporal increasing rate in the invaded ( u t ) shall be high; otherwise, the invasion will not succeed. As the invasive species proliferates over time, there shall be a limit in its concentration; then, the temporal ratio (considered positive) decreases because the invasive species progressively reaches its maximum stationary level in the medium. In addition, the invasive temporal rate depends on the invaded concentration (v). This behaviour is intended to be modelled as a reaction term of the form:
u t = v ( κ u ) ,
where κ represents the maximum invasive species concentration in the given medium. Without loss of generality, assume κ = 1 .
The invaded species time vanishing rate is considered as:
v t = v u ,
and aims to express that the decreasing rate is absolutely higher for increasing values of the invasive species u and is dependent on its own invaded species existence in the medium.
An advection term is proposed to account for possible forced movement in the media. The interaction of the advection, through the vector c, has the intention of modelling the desertion behaviour in the invaded species, so that the runaway is oriented in the direction of c and toward spatial areas not previously populated (represented by c · v ). In response to the invaded species desertion, the invasive species will follow an advection process given by the orientation c towards the same non-populated areas (represented by c · u ).
In addition, a nonlinear diffusion (Porous Medium Equation) is considered to model the random movement of both species. As an alternative to the classical linear order two diffusion, the Porous Medium Equation introduces concepts such as the finite speed of propagation in the invaded species, which is admitted to account for a further reliable model.
The study of solutions for a Porous Medium Equation with advection and no reaction has been widely analyzed in [23]. Nonetheless, the present work considers the effect of the coupled reaction and absorption for which the mixed monotone behaviour induces a different analytical treatment compared to that in [23].
This work begins by showing some regularity results together with uniqueness of solutions. The degeneracy of the diffusivity leads to considering weak solutions defined upon a test function ϕ C ( R d ) . Afterwards, profiles of maximal and minimal solutions are shown to exist and analytically obtained based on the mixed monotone properties in the reaction–absorption terms. In addition, the problem P is studied within the geometric perturbation theory to support the construction of analytical TW profiles. Finally, the finite propagation and the exponential decay for v are shown in the proximity of the null solution, i.e., when v ϵ 0 + . The finite propagation is proved with a maximal solution so that any other lower solution exhibits finite propagation as well. In addition, the finite propagation is shown to provide a compact support that is mainly dependent on the advection (invaded runaway) c. This means that the homogenizing effect of the diffusive front is lost and the advection term drives the dynamic. This behaviour can be interpreted as the invasive species concentrating efforts in the runaway direction c so that the invaded species decreases exponentially even when the diffusion is nonlinear.

3. Existence and Uniqueness of Solutions

Let us consider a test function ϕ C ( R d ) such that, for 0 < τ < t < T :
R d u ( t ) ϕ ( t ) = R d u ( τ ) ϕ ( τ ) + τ t R d [ u ϕ t + u m Δ ϕ c · ϕ u + v ( 1 u ) ϕ ] d s ,
R d v ( t ) ϕ ( t ) = R d v ( τ ) ϕ ( τ ) + τ t R d [ v ϕ t + v m Δ ϕ c · ϕ v u v ϕ ] d s .
Considering r > > r 0 > 0 , the following uniformly parabolic quasilinear set of equations is defined (named P ϕ ):
u ϕ t + u m Δ ϕ c · ϕ u + v ( 1 u ) ϕ = 0 , v ϕ t + v m Δ ϕ c · ϕ v u v ϕ = 0 ,
in B r × [ 0 , T ] , with the following set of boundary and initial conditions:
( ϕ u c ) · ν = 0 , ( ϕ v c ) · ν = 0 ,
over B r × [ 0 , T ] , where ν is the outer unitary normal vector in B r and
u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) ,
initially on B r .
The problem P ϕ has existence and uniqueness of solutions based on the mixed monotone properties of the forcing terms [24]. For this purpose, note that the forcing part ϕ is a monotone increasing (with ϕ ) function and ϕ is a monotone decreasing function.
Theorem 1.
Given u 0 ( x ) , v 0 ( x ) L 1 ( R d ) L ( R d ) , the set of solutions u ( x , t ) , v ( x , t ) are bounded for all ( x , t ) B r × [ 0 , T ] with r > > 1 .
Proof. 
Consider η error and sufficiently small, the following cut off function is defined [25]:
ψ η C 0 ( R d ) , 0 ψ η 1 , ψ η = 1 i n B r η , ψ η = 0 i n R d B r η ,
so that
| ψ η | c a η , | Δ ψ η | c a η 2 .
After multiplication of (7) by ψ η and integrating in B r × [ τ , T ] :
τ t B r u ϕ t ψ η + τ t B r u m Δ ϕ ψ η τ t B r c · ϕ u ψ η + τ t B r v ( 1 u ) ϕ ψ η = 0 , τ t B r v ϕ t ψ η + τ t B r v m Δ ϕ ψ η τ t B r c · ϕ v ψ η τ t B r u v ψ η ϕ = 0 .
The integrals for the advection term read:
τ t B r c · ϕ u ψ η = τ t B r c · ϕ u ψ η + τ t B r c · ϕ u ψ η ,
For some large r > > r 0 > 1 [25]:
τ t u m c 1 ( τ ) r 2 m m 1 ,
and
τ t u c 2 ( τ ) r 2 m 1 .
Making the integration over fixed B r domains for r > > r 0 > > 1 and admitting τ t , the functions ϕ , ψ η are considered stationary locally in time. Then, the integral term in the diffusion reads:
τ t B r u m Δ ϕ ψ η B r c 1 ( τ ) r 2 m m 1 Δ ϕ ψ η .
Given a fixed r, the following holds:
B r c 1 ( τ ) r 2 m m 1 Δ ϕ ψ η = c 1 ( τ ) r 2 m m 1 ϕ ψ η B r B r ϕ · ψ η .
Note that ( ϕ ψ η ) B r < < 1 in (17) for r > > 1 . Indeed, ψ η is a bounded non-increasing function and ϕ 0 for r > > 1 as it it will be shown afterwards for a particular choice of the test function. Now:
B r c 1 ( τ ) r 2 m m 1 ϕ · ψ η B r c 1 ( τ ) r 2 m m 1 | ϕ | c a η c 1 ( τ ) c a B r r 2 m m 1 1 | ϕ | .
Similarly for the the integrals in (13) and after arranging the spatial leading terms:
τ t B r c · ϕ u ψ η + τ t B r c · ϕ u ψ η c c 2 ( τ ) τ t r 2 m 1 1 | ϕ | + c c 2 ( τ ) c a B r r 2 m 1 1 ϕ .
Next, consider a test function ϕ of the form:
ϕ ( x , s ) = e g ( s ) ( 1 + r 2 ) β ,
where g ( s ) > 0 for 0 < s < t , g ( t ) = 0 and β shall be chosen to ensure the convergence of the integrals in (18) and (19) as r . To this end, it suffices to consider:
β = 1 m 1 .
Returning to (12):
τ t B r u ϕ t ψ η + τ t B r v ( 1 u ) ϕ ψ η = τ t B r c · ϕ u ψ η + τ t B r c · ϕ u ψ η + τ t B r u m ϕ · ψ η , τ t B r v ϕ t ψ η τ t B r u v ψ η ϕ = τ t B r c · ϕ v ψ η + τ t B r c · ϕ v ψ η , + τ t B r v m ϕ · ψ η .
The intention is to show that the left-hand integrals in (22) are finite, which is equivalent to searching for the bound on the right-hand side terms. Indeed, and for the particular expression of β in (21), the right-hand integrals are indeed bounded for r > > r 0 > 1 over B r domains. Then,
τ t B r u ϕ t ψ η + τ t B r v ( 1 u ) ϕ ψ η A ( τ , t ) , τ t B r v ϕ t ψ η τ t B r u v ψ η ϕ B ( τ , t ) ,
for some finite locally in time defined A , B .
As both integrals are finite in τ < s < t < T , it is possible to conclude on the theorem postulation about the boundness of solutions in B r × [ 0 , ) .
The next intention is to show the uniqueness of solutions:
Theorem 2.
Assuming that ( u , v ) > ( 0 , 0 ) is a minimal solution for P (2) in R d × ( 0 , T ) , then ( u , v ) coincides with the maximal solution, i.e., the solution is unique.
Proof. 
Assume ( u ^ , v ^ ) to be the maximal solution to P in R d × ( 0 , T ) , such that:
( u ^ ( x , 0 ) , v ^ ( x , 0 ) ) = ( u 0 ( x ) + ν , v 0 ( x ) + ν ) ,
with ν > 0 being arbitrarily small. In addition, define the minimal solution:
( u ( x , 0 ) , v ( x , 0 ) ) = ( u 0 ( x ) , v 0 ( x ) ) ,
such that the maximal and minimal solutions verify respectively:
u ^ t = Δ u ^ m + c · u ^ + v ^ ( 1 u ) , v ^ t = Δ v ^ m + c · v ^ u v .
u t = Δ u m + c · u + v ( 1 u ^ ) , v t = Δ v m + c · v u ^ v ^ .
Then, for every test function ϕ C ( R d ) , the following expressions hold:
0 R d ( u ^ u ) ( t ) ϕ ( t ) = 0 t R d [ ( u ^ u ) ϕ t + ( u ^ m u m ) Δ ϕ + c · ϕ ( u ^ u ) + ( v ^ ( 1 u ) v ( 1 u ^ ) ϕ ] d s ,
0 R d ( v ^ v ) ( t ) ϕ ( t ) = 0 t R d [ ( v ^ v ) ϕ t + ( v ^ m v m ) Δ ϕ + c · ϕ ( v ^ v ) + ( u ^ v ^ u v ) ϕ ] d s .
Let us define:
a 1 ( x , s ) = u ^ m u m u ^ u for u ^ u m u m 1 otherwise ,
a 2 ( x , s ) = v ^ m v m v ^ v for v ^ v m v m 1 otherwise .
Given two fixed values for x and s = t T :
0 a 1 ( x , s ) c 1 ( m , u 0 , T ) ,
0 a 2 ( x , s ) c 2 ( m , v 0 , T ) .
Consider the test function:
ϕ ( | x | , s ) = e k ( T s ) ( 1 + | x | 2 ) γ ,
for some constants k and γ .
The test function verifies:
ϕ t = k ϕ ( x , s ) , | | x | ϕ | c 3 ( γ , d ) ϕ ( x , s ) , Δ | x | ϕ c 4 ( γ , d ) ϕ ( x , s ) ,
such that:
( u ^ u ) ϕ t + ( u ^ m u m ) Δ ϕ + c · ϕ ( u ^ u ) + ( v ^ ( 1 u ) v ( 1 u ^ ) ) ϕ ( u ^ u ) k ϕ + a 1 ( u ^ u ) c 4 ϕ + c c 3 ϕ ( u ^ u ) + L 1 ( v ^ v ) ϕ ,
( v ^ v ) ϕ t + ( v ^ m v m ) Δ ϕ + c · ϕ ( v ^ v ) + ( u ^ v ^ u v ) ϕ ( v ^ v ) k ϕ + a 2 ( v ^ v ) c 4 ϕ + c c 3 ϕ ( v ^ v ) + L 2 ( u ^ u ) ϕ ,
where L 1 = max R d { ( 1 u ) } and L 2 = max R d { v ^ } are bounded as per Theorem 1.
The constant k shall be selected such that:
k + a 1 c 4 + c 3 c 0 , k + a 2 c 4 + c 3 c 0 .
It suffices to consider:
k max ( a 1 c 4 + c 3 c , a 2 c 4 + c 3 c ) ,
so that:
( k + a 1 c 4 + c 3 c ) ϕ ( u ^ u ) 0 , ( k + a 1 c 4 + c 3 c ) ϕ ( v ^ v ) 0 .
The expressions (28) and (29) read:
0 R d ( u ^ u ) ( t ) ϕ ( t ) 0 t R d L 1 ( v ^ v ) ϕ d s ,
0 R d ( v ^ v ) ( t ) ϕ ( t ) 0 t R d L 2 ( u ^ u ) ϕ d s .
Make the d / d t in the previous expressions:
0 d d t R d ( u ^ u ) ( t ) ϕ ( t ) R d L 1 ( v ^ v ) ( t ) ϕ ( t ) ,
0 d d t R d ( v ^ v ) ( t ) ϕ ( t ) R d L 2 ( u ^ u ) ( t ) ϕ ( t ) .
Letting us operate in the last two expressions with the equality, then:
R d ( u ^ u ) ( t ) ϕ ( t ) = 1 L 2 d d t R d ( v ^ v ) ( t ) ϕ ( t ) ,
which can be replaced in (43):
d 2 d t 2 R d ( v ^ v ) ( t ) ϕ ( t ) L 1 L 2 R d ( v ^ v ) ( t ) ϕ ( t ) .
Define:
g ( t ) = R d ( v ^ v ) ( t ) ϕ ( t ) ;
then, Equation (46) reads as a linear standard second order equation:
d 2 g ( t ) d t 2 = L 1 L 2 g ( t ) ,
with
g ( 0 ) = ν 0 , g ( 0 ) = ν 0 .
Then, Equation (48) converges to the null solution. For the sake of simplicity, consider g ( t ) = 0 , concluding that:
v ^ = v .
Now, returning to (41):
0 R d ( u ^ u ) ( t ) ϕ ( t ) 0 ,
to conclude that:
u ^ = u ,
showing, then, the uniqueness of any positive solution to P. □

4. Profiles of Solution

This section is devoted to the searching of upper and lower solutions considering that the invasive species u ( x , t ) shall be non-decreasing (i.e., u t 0 ) while the invaded one v ( x , t ) shall be non-increasing (i.e., v t 0 ). Note that in this chapter the advection coefficient may be considered as a 1D parameter when required by the context. This is an important hypothesis and reflects a constant invasive species behaviour when entering into a domain (and at the same time a constant behaviour in the invaded one). The constant advection hypothesis may be discussed further; for example, in [26], the advection coefficient is determined in a porous medium as the relationship between the water flow velocity and the range of volumetric water content. Probably a similar process can be followed to model the invasive-invaded dynamics involved, where water flow velocity can be regarded as the invasive species velocity of invasion and the volumetric water content as the range of total concentration.
Theorem 3.
A maximal solution v ^ ( x , t ) for the invaded species is:
v ^ ( x , t ) = v 0 ( x ) ( A t ) 1 m 1 ,
where
A δ m 1 | x | c ,
being δ = min { v } > 0 , | x | representing the spatial integration domain and 0 < t < | x | c .
In addition, v ^ propagates through the support in ( x , t ) defined by the minimal approximation:
| x | min = δ m 1 c 2 m m 1 1 / 3 t 2 / 3 .
Proof. 
Let us start with the equation in v. Considering self-similar solutions:
v ^ ( x , t ) = t α f ( ξ ) , ξ = | x | t β ,
to the upper profile equation:
v ^ t = Δ v ^ m + c · v ^ .
Then, the following holds:
α t α 1 f + β ξ t α 1 f = G ( f , f , f , m ) t α m + 2 β + c t α + β f .
Making the corresponding equalities in the time exponents:
α = 1 m 1 ; β = 1 .
The self-similar profile shall, then, be obtained as a solution to the elliptic equation:
α f + β ξ f = ( f m ) + d 1 ξ ( f m ) + c f .
Note that in the approximation:
β ξ > > c , | x | t 1 > > | c | ,
and the following elliptic equation holds for t < < | x | c :
β ξ f = ( f m ) + d 1 ξ ( f m ) + α f ,
for which a solution is available [27,28]:
f ( ξ ) = ( A B ξ 2 ) 1 m 1 ,
where A > 0 and B = ( m 1 ) 2 m . A is obtained by making ξ = 0 in (63), so that:
v ^ ( x , t ) = A 1 m 1 t α .
Note that (58) provides supersolutions under the condition that each species’ concentration is positive, hence min { u v } > 0 . As u t > 0 and assuming u 0 > 0 , then it is required that min { v } > 0 . Considering δ = min { v } > 0 , then:
v ^ ( x , t ) = A 1 m 1 t α > δ ,
so that:
A = δ m 1 t δ m 1 | x | c ,
in the inner region t < < | x | c where | x | shall be understood as the maximum spatial variable representing the integration domain.
The maximal time evolution v ^ is then given by:
v ^ ( x , t ) = v 0 ( x ) ( A t ) 1 m 1 .
Finally, v ^ is indeed an upper evolution. To show this, consider a solution v to the original equation given by the absorption term u v and define a t 0 :
0 < t 0 < | x | c ,
such that, for any τ with t 0 < τ < | x | c :
v ( x , τ ) v ^ ( x , t 0 ) .
Now, for any t with τ < t < | x | c :
v ( x , τ + t ) v ^ ( x , t 0 + t ) .
In the limit, τ 0 :
v ( x , t ) v ^ ( x , t ) ,
showing the upper behaviour of v ^ compared to v.
The propagating support is obtained making:
f ( ξ ) = 0 , A = B ξ 2 ,
so that A and B can be replaced by their corresponding values. Then, the support propagates as:
| x | δ m 1 c 2 m m 1 1 / 3 t 2 / 3 = | x | min .
Note that solutions are of the form (53) in R d × [ 0 , ] , so that, in the proximity of δ 0 + , the following time assessment holds:
t v = v 0 δ m 1 A .
The expression (74) provides an estimation in time to consider a vanishing condition for v, so that, for t > t v , the equation for u ˜ u reads:
u ˜ t = Δ u ˜ m + c · u ˜ ,
and the following theorem holds in the search for a lower solution u ˜ :
Theorem 4.
Consider
A 1 = max x R d { u 0 ( x ) } .
The lower solution u ˜ ( x , t ) for the invasive species reads:
u ˜ ( x , t ) = u 0 ( x ) + A 1 t 1 m 1 ,
for
t v < t < | x | c .
Proof. 
Again, consider self-similar solutions of the form:
u ˜ ( x , t ) = t α F ( ξ ) , ξ = | x | t β ,
to the lower profile Equation (75). Upon substitution of the self-similar solution (79) and operating analogously as in (58), the following self-similar exponents are deterred:
α = 1 m 1 ; β = 1 .
Assuming t < | x | c , the self-similar profile F is obtained as a solution to the elliptic equation:
β ξ F = ( F m ) + d 1 ξ ( F m ) + α F ,
for which the solution is, again, provided in [27,28]:
F ( ξ ) = ( A 1 B 1 ξ 2 ) 1 m 1 ,
for t v < t < | x | c and A 1 > 0 , B = ( m 1 ) 2 m .
Making ξ = 0 , the solution for u ˜ is:
u ˜ ( x , t ) = u 0 ( x ) + A 1 t 1 m 1 ,
where
A 1 = max x R d { u 0 ( x ) } .
Finally, the expression in (83) is shown to be a lower solution. For such purpose, let us define t 0 , such that t v < t 0 < | x | c , and τ such that t 0 < τ < | x | c , then:
u ˜ ( x , τ ) u ( x , t 0 ) .
This last expression is valid provided | t 0 τ | < < 1 . Considering t with τ < t < | x | c , then:
u ˜ ( x , τ + t ) u ( x , t 0 + t ) .
After a time translation t = t 0 + t :
u ˜ ( x , t ) u ( x , t ) ,
with t v < t 2 < | x | c . □
The next objective is to define a maximal solution for u and a minimal solution for v. To this end, the following theorem holds:
Theorem 5.
There exists a minimal solution v ˜ ( x , t ) and a maximal solution u ^ ( x , t ) in R d × [ 0 , T ] . Furthermore, consider:
v ˜ = γ = min x R d { v 0 ( x ) } ,
Then, a maximal solution for u is:
u ^ ( x , t ) = 1 ( 1 u 0 ( x ) ) e γ t .
Proof. 
For building the maximal solution u ^ ( x , t ) , consider the following problem in u:
u t ϵ = Δ ( u ϵ ) m + c · u ϵ + v ( 1 u ϵ ) , u ϵ ( x , 0 ) = u 0 ( x ) + ϵ , v ( x , 0 ) = v 0 ( x ) , u 0 ( x ) > 0 . v 0 ( x ) > 0 L 1 ( R d ) L ( R d ) .
Solutions to problem (90) do exist under the Lipschitz condition in the reaction term and the positive initial conditions (refer to Theorem 1).
Let us consider ϵ 2 < ϵ , so that:
u t ϵ = Δ ( u ϵ ) m + c · u ϵ + v ( 1 u ϵ ) , u t ϵ 2 = Δ ( u ϵ 2 ) m + c · u ϵ 2 + v ( 1 u ϵ ) , u ϵ ( x , 0 ) = u 0 ( x ) + ϵ > u ϵ 2 ( x , 0 ) = u 0 ( x ) + ϵ 2 > 0 .
The iteration process starts with ϵ , and, based on the monotony behaviour in the reaction terms of (91), then u ϵ 2 < u ϵ in accordance with the initial data for u ϵ . The same argument can be repeated for ϵ 3 < ϵ 2 , so that u ϵ 3 < u ϵ 2 . The sequence defined as { u ϵ } is non-negative and non-increasing. Consequently, the following maximal solution in the limit is defined:
u ^ = lim ϵ 0 u ϵ .
To construct the minimal solution for v, consider the problem:
v t ϵ = Δ ( v ϵ ) m + c · v ϵ ( lim ϵ 0 u ϵ ) v ϵ v ϵ ( x , 0 ) = v ( x , 0 ) = v 0 ( x ) > 0 .
Solutions to (93) do exist (refer to Theorem 1).
Considering ϵ 2 < ϵ , then:
u ϵ 2 < u ϵ v ϵ 2 > v ϵ ,
based on the monotony properties in the reaction term of (93). The sequence { v ϵ } is non-decreasing as ϵ 0 . The minimal solutions for v are defined as per the monotone limit:
v ˜ = lim ϵ 0 v ϵ .
Once the maximal and minimal solutions have been constructed and shown to exist, the next intention is to search for a family of flat solutions u ϵ , v ϵ via the resolution of:
u t ϵ = γ ( 1 u ϵ ) , v ϵ = γ = min x R d { v 0 ( x ) } , u ϵ ( x , 0 ) = u 0 ( x ) + ϵ , v ϵ ( x , 0 ) = v 0 ( x ) , u 0 ( x ) > 0 . v 0 ( x ) > 0 L 1 ( R d ) L ( R d ) .
This last system can be solved by separation of variables, so that:
u ϵ ( x , t ) = 1 ( 1 u 0 ( x ) + ϵ ) e γ t .
Then, for ϵ 0 ,
u ^ ( x , t ) = 1 ( 1 u 0 ( x ) ) e γ t .
On a physical basis, a value for γ shall be determined experimentally and can represent a constant homogeneous initial distribution of the invaded species in the domain. Depending on the particular application, either the maximal or minimal solution may be considered. For instance, if the intention is to model the minimal level of invaded species that can face an invasive one, then the solution v ˜ in (88) shall be taken into account together with the solution u ^ in (89). Such maximal solution would make sense in this case (minimal invaded population) as the resistance to the invasion is relatively low, leading the invasive species to prosper quickly via a maximal solution u ^ . In addition, note that the maximal solution for the invaded species (and hence any other solution below) propagates through the support given in (55) which elucidates the finite propagation feature. Note that in the inner region (inner compared to the diffusive front moving with finite propagation), the desertion will guide the invaded species to move along the preferred direction c. If c is high, then the propagating support (55), induced by the diffusion, will expand further.

5. Travelling Waves Analysis, Existence and Regularity

The TW profiles are provided by the following structured changes: u ( x , t ) = f ( ξ ) , such that ξ = x · n d a t R . Note that n d is a unitary vector in R d defining the propagating direction. In addition, a is the TW velocity and the profile f : R ( 0 , ) is such that f L ( R d ) . For the sake of simplicity, the vector n d = ( 1 , 0 , . . . , 0 ) , then u ( x , t ) = f ( ξ ) , ξ = x a t R .
Previous to the formal structure of TW solutions, let us consider some remarks on TW symmetries. Note that two TW are equivalent under translation ξ ξ + ξ 0 and symmetry ξ ξ as the D’Alambert waves coordinate x a t is invariant under the translation group. It shall be noted that such invariant symmetry does not necessarily hold when the advection coefficient, c, is time dependent or when time-periodic boundaries are given. In addition, and under these last two cases, monotonic TW solutions may not hold (the reader can refer to [29] where periodic TW-solutions are analyzed in detail for a time-dependent advection coefficient).
Considering u ( x , t ) = f ( ξ ) and v ( x , t ) = g ( ξ ) , then the problem P (2) is transformed to the TW domain:
a f = ( f m ) + c f + g ( 1 f ) , a g = ( g m ) + c g f g , f , g L ( R ) , f ( ξ ) > 0 , g ( ξ ) < 0 , f ( ) = 1 , g ( ) = 0 .
Working with the density and flux variables
X 1 = f , Y 1 = ( f m ) , X 2 = g , Y 2 = ( g m ) ,
the following system holds:
X 1 = 1 m X 1 1 m Y 1 , Y 1 = ( a + c ) 1 m X 1 1 m Y 1 X 2 ( 1 X 1 ) , X 2 = 1 m X 2 1 m Y 2 , Y 2 = X 1 X 2 ( a + c ) 1 m X 2 1 m Y 2 ,
with the critical point ( 1 , 0 , 0 , 0 ) that represents a situation in which the invasive species wins over the invaded.

Geometric Perturbation Theory

The singular geometric perturbation theory is employed in this section to show the asymptotic behaviour of a two-dimensional manifold defined to make simpler the assessment of a TW analytical profile.
For this purpose, define the two-dimensional manifold as:
M 0 = { X 1 , Y 1 , X 2 , Y 2 / X 1 = 1 m X 1 1 m Y 1 ; Y 1 = ( a + c ) 1 m X 1 1 m Y 1 } ,
with critical point ( 1 , 0 , 0 , 0 ) .
The perturbed manifold M ϵ close to M 0 is defined as:
M ϵ = { X 1 , Y 1 , X 2 , Y 2 / X 2 = ϵ ; Y 2 = 0 ; X 1 = 1 m X 1 1 m Y 1 ; Y 1 = ( a + c ) 1 m X 1 1 m Y 1 ϵ ( 1 X 1 ) } .
The intention is to use the Fenichel invariant manifold theorem [30] as formulated in [31,32]. For this purpose, it is required to show that M 0 is a normally hyperbolic manifold, i.e., the eigenvalues of M 0 in the linearized frame close to the critical point, and transversal to the tangent space, have a non-zero real part. This is shown based on the following two-dimensional equivalent flow associated with M 0 :
X 1 Y 1 = 0 1 m 0 a + c m X 1 Y 1 .
where the eigenvalues are ( 0 , ( a + c ) / m ) . For λ = 0 , the eigenvector is [ 1 , 0 ] which is tangent to M 0 to conclude that M 0 is a hyperbolic manifold. The next objective is to prove that M ϵ is locally invariant under the equations (flows) in (101). To this end, and according to [32]: for all R > 0 , for all open interval J with ( a + c ) J and for any value of i N , there exists a δ such that, for ϵ ( 0 , δ ) , the manifold M ϵ is invariant. Hence, consider i 1 and the functions:
ϕ 1 = ϵ , ϕ 2 = 0 , ϕ 3 = 1 m X 1 1 m Y 1 , ϕ 4 = ( a + c ) 1 m X 1 1 m Y 1 ϵ ( 1 X 1 ) ,
which are C i ( B R ( 0 ) ¯ × I ¯ × [ 0 , δ ] ) close to the critical condition ( 1 , 0 , 0 , 0 ) .
A value for R > 0 can be selected so as to assume that M 0 B R ( 0 ) is large enough to study the whole TW evolution. Note that δ is determined based on the computation of the distance between the flows in M ϵ and M 0 . To this end, assume that the functions in each of the flows are measurable a.e. in B R ( 0 ) :
ϕ 4 M 0 ϕ 4 M ϵ δ 1 X 1 .
The distance between the manifolds keeps the normal hyperbolic condition for δ ( 0 , ) . For simplicity, assume δ = 1 .
In the same way, given the two-dimensional manifold M 1 :
M 1 = { X 1 , Y 1 , X 2 ϵ , Y 2 / X 2 = 1 m X 2 1 m Y 2 ; Y 2 = X 2 ( a + c ) 1 m X 2 1 m Y 2 } ,
with the same critical point ( 1 , 0 , 0 , 0 ) , and the perturbed M ϵ close to M 1 :
M ϵ = { X 1 , Y 1 , X 2 , Y 2 / X 1 = 1 ; Y 1 = 0 ; X 2 = 1 m X 2 1 m Y 2 ; Y 2 = X 1 X 2 ( a + c ) 1 m X 2 1 m Y 2 } .
The Fenichel invariant manifold theorem can be applied in the same manner as for M 0 . Note that the two-dimensional equivalent flow associated with M 1 is:
X 2 Y 2 = 0 ϵ 1 m m 1 ( a + c ) ϵ 1 m m X 2 Y 2 .
The associated eigenvalues are both real ( a + c ) ϵ 1 m 2 m ± 1 2 ( ( a + c ) ϵ 1 m ) 2 m 2 + 4 ϵ 1 m m . Hence, M 1 is a hyperbolic manifold.
In the same manner, the next intention is to show that the manifold M ϵ is locally invariant under the flow (101) so that the manifold M 1 can be represented as an asymptotic approach to M ϵ . For this purpose, consider the functions:
ϕ 1 = 1 , ϕ 2 = 0 , ϕ 3 = 1 m X 2 1 m Y 2 , ϕ 4 = Y 2 = X 1 X 2 ( a + c ) 1 m X 2 1 m Y 2 ,
which are C i ( B R ( 0 ) ¯ × I ¯ × [ 0 , δ ] ) in the proximity of the critical point ( 1 , 0 , 0 , 0 ) .
In this case, δ is determined based on the following flows that are considered to be measurable a.e. in B R ( 0 ) :
ϕ 4 M 1 ϕ 4 M ϵ δ ϵ 1 X 1 = δ 1 X 1 .
The normal hyperbolic condition for both manifolds M 1 and M ϵ applies for δ ( 0 , ) . In the sake of simplicity, consider δ = 1 .
Once we have shown that M 0 and M 1 , defined as two-dimensional manifolds, remain invariant with regard to the four-dimensional manifolds M ϵ and M ϵ respectively under the flow (101), the TW profiles can be obtained operating in M 0 and M 1 .

6. Travelling Waves Profiles and Finite Propagation

Based on the normal hyperbolic condition in the manifold M 0 under the flow (101), asymptotic TW profiles can be obtained. For this purpose, consider firstly:
X 1 = 1 m X 1 1 m Y 1 ; Y 1 = ( a + c ) 1 m X 1 1 m Y 1 X 2 ( 1 X 1 ) ,
such that the following family of trajectories in the phase plane ( X 1 , Y 1 ) holds:
d Y 1 d X 1 = ( a + c ) m X 2 X 1 m 1 1 X 1 Y 1 = H ( a , c , X 1 , X 2 , Y 1 ) ,
where 0 < X 1 < 1 ; 0 < X 2 < 1 ; Y 1 > 0 .
The existence of a convergent trajectory is shown based on a comparison with subsolutions for ( a + c ) sufficiently small and supersolutions for ( a + c ) sufficiently large together with topological assessment and the continuity of H.
For this purpose, note that, when ( a + c ) 0 , then d Y 1 / d X 1 < 0 , while, when ( a + c ) < 0 with | a + c | > > 1 , d Y 1 / d X 1 > 0 . Given the continuity of H, it is thus possible to make a conclusion about the existence of a critical trajectory of the form:
( a + c ) m X 2 X 1 m 1 1 X 1 Y 1 = H ( a , c , X 1 , X 2 , Y 1 ) = 0 .
Assuming now that the selected TW speed satisfies a < < c with | a + c | > > 1 , then the same solution applies with minor effects due to the advection term, which permits assessing a TW speed only making assumptions on the bound of solutions; indeed:
a m X 2 X 1 m 1 1 X 1 Y 1 ,
for which a maximal bound is obtained as:
a < m A ,
being A = min B R ( 0 ) { Y 1 } . Now, coming to (114) and under the infinitesimal asymptotic approach X 2 ϵ , the following holds:
Y 1 = m | a + c | ϵ X 1 m 1 ( 1 X 1 ) ,
in the approach X 1 1 :
Y 1 B ( ϵ ) ( 1 X 1 ) ,
for a sufficiently small B ( ϵ ) . Equivalently, in the asymptotic approximation:
( f m ) + B f = 0 f ( m f m 2 f + B ) = 0 ,
where = d d ξ . This last equation has the solution:
f ( ξ ) = C ( α B ξ ) 1 m 1 ,
where C = m 1 m 1 m 1 and α > 0 . Note that the growing TW is obtained replacing ( ξ ) with the symmetric ( ξ ). This is a meaningful result that expresses the behaviour of the invaded species ( ξ 1 m 1 ) in the proximity of the critical point.
The same process shall be repeated for the manifold M 1 under the flow (101):
X 2 = 1 m X 2 1 m Y 2 ; Y 2 = X 2 ( a + c ) 1 m X 2 1 m Y 2 ,
such that in the phase plane ( X 2 , Y 2 ) :
d Y 2 d X 2 = a + c m + m X 2 m Y 2 1 = G ( a , c , X 2 , Y 2 ) ,
where 0 < X 2 < 1 ; Y 2 < 0 .
In the same manner, the existence of a convergent trajectory is shown based on a topological argument and the continuity of G. For ( a + c ) > 0 , d Y 2 / d X 2 < 0 , while, when ( a + c ) < 0 with | a + c | > > 1 , d Y 2 / d X 2 > 0 . Given the continuity of G, a critical trajectory is given:
a + c m + m X 2 m Y 2 1 = G ( a , c , X 2 , Y 2 ) = 0 ,
such that
Y 2 = m 2 | a + c | X 2 m ( g m ) = m 2 | a + c | g m .
The last equation can be solved with standard means:
g ( ξ ) = D e m | a + c | ξ ,
for D > 0 .
The positivity condition for f (120) permits conclusions about some regularity results in the quasilinear parabolic operator (see Section 3). In addition, the species g verifies
g ϵ 0 + i n B R T = B R ( x 0 . R ) × [ T ϵ , T + ϵ ] ,
for T > > 1 , which means the existence of a convergent tail in B R T towards the null condition in v. The next objective is to show the existence of finite propagation in the invaded species v.
Theorem 6.
Given m > 2 , a finite propagating support exists whenever:
f ϵ 0 + i n B R T = B R ( x 0 . R ) × [ T ϵ , T + ϵ ] ,
T > 0 .
Proof. 
Firstly, assume the pressure variable w:
w = m m 1 v m 1 ,
so that the equation for v in (2) reads:
w t = | w | 2 + ( m 1 ) w Δ w + c · w + m 1 m w ( m 1 m ) m 1 m w m 1 m ,
where w 0 + , then:
w t | w | 2 + c · w .
A solution to a similar equation has been provided in [27]. Consider the following function in the search of a maximal solution:
W ( x , t ) = a b t + r 1 n + , r = | x | , n N .
Both a and b > 0 are constants to be determined. In particular, for 0 τ 1 , impose:
b τ = 1 2 n .
Under this condition:
W ( x , t ) 0 f o r r < 1 2 n and 0 t τ .
Any solution to Equation (130) is bounded as per Theorem 1, then:
v ( x , t ) K 1 f o r x error , 0 t τ and K 1 ( p , u 0 ) .
The intention is to make W ( x , t ) as a maximal solution:
W ( x , t ) v ( x , t ) ,
so that
a b t + r 1 n + K 1 .
Select any r > 1 n , for example r = 2 n . Thus, for t = 0 :
a 2 n 1 n + K 1 , a n K 1 .
Note that:
W ( x , t ) v ( x , t ) ,
in r = 2 n and 0 t τ . The value of b shall be chosen in such a way that W ( x , t ) is a supersolution in 0 < r < 2 n , 0 t τ :
W t | W | 2 + c · W ,
and considering that:
W t = a b , W r = a ,
the following value for b is obtained:
b a + c .
For the values of a and b in expressions (137) and (141), respectively, the function W ( x , t ) is a supersolution locally:
W ( x , t ) w ( x , t ) , 0 < | x | < 2 n , 0 t τ .
The inequality (142) permits concluding that any local supersolution satisfies the null criteria in B R T ; then, any minimal solution w ( x , t ) satisfies such null criteria in B R T . □
The expression (125) provides the characteristic profile in relation to the diffusion front (note the parameter m) and the advection. Once the invaded species starts the desertion in the direction of c, it follows a curve in ( x , t ) given by (125). The invasive species movement concentrates in the same trajectory to reduce the invaded population that follows an exponential decay. Furthermore, the existence of a finite propagation suggests that, during the desertion, the diffusion is still relevant in the proximity of the null solution (i.e., the invaded tail). This can be interpreted as the existence of a random movement of invaded species in the tail where the invasive species influence is negligible.

7. Conclusions

The proposed problem P (2) has been discussed stressing aspects related with existence, uniqueness, behaviour of minimal and maximal solutions and Travelling Waves supported by the geometric perturbation theory. In addition, the finite speed of propagation, induced by the porous medium diffusion, has been shown and a characterization of such property has been explored. The propagation features of the species v when approaching the null solution have been shown. To this end, an exponential decreasing tail in the TW domain has been proved to exist. The invaded species trajectory, to escape from the invasive one, is given by the TW solution mentioned. Even when the invaded desertion is mainly governed by the advection, the nonlinear diffusion still acts in the tail leading to a finite speed propagating front.

Author Contributions

Formal analysis, J.L.D.P. and S.U.R.; Funding acquisition, J.L.D.P. and J.R.G.; Investigation, J.L.D.P., J.R.G. and S.U.R.; Methodology, J.L.D.P. and S.U.R.; Resources, J.R.G. and A.N.R.; Validation, A.N.R.; Writing—original draft, J.L.D.P. and A.N.R.; Writing—review & editing, J.L.D.P. and J.R.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received the support from the signing institutions and the funding from the project ID-UDIMA-2021-02.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mooney, H.; Williamson, M. The Problem of Biological Invasions; Oxford Scholarship; Oxford University Press: Oxford, UK, 2010. [Google Scholar]
  2. Perumpanani, A.J.; Sherratt, J.A.; Norbury, J.; Byrne, H.M. A two parameter family of travelling waves with a singular barrier arising from the modelling of extracellular matrix mediated cellular invasion. Physica D 1999, 12, 145–159. [Google Scholar] [CrossRef]
  3. Harley, K.; Van Heijster, R.; Marangell, R.; Pettet, G.; Roberts, T.; Wechselberger, M. Instability of Travelling Waves in a Model of Haptotaxis. SIAM J. Appl. Math. 2020, 80, 1629–1653. [Google Scholar] [CrossRef]
  4. Berestycki, H. The influence of advection on the propagation of fronts in reaction-diffusion equations. In Nonlinear PDE’s in Condensed Matter and Reactive Flows; Springer: Dordrecht, The Netherlands, 2002; pp. 11–48. [Google Scholar]
  5. Malaguti, L.; Emilia, R. Travelling Wave fronts in Reaction–Diffusion Equations with Convection Effects and Non–Regular Terms. Math. Nachr. 2002, 242, 148–164. [Google Scholar] [CrossRef]
  6. Banani, R.; Sankar, K.R.; Biswas, M.H.A. Effects on prey–predator with different functional responses. Int. J. Biomath. 2017, 10, 1750113. [Google Scholar]
  7. Li, H. Hopf Bifurcation of Delayed Density-Dependent Predator-Prey Model. Acta Math. Sci. Ser. A 2019, 39, 358–371. [Google Scholar]
  8. Zu, L.; Jiang, D.; O’Regan, D. Periodic Solution for a Stochastic Non-autonomous Predator-Prey Model with Holling II Functional Response. Acta Appl. Math. 2019, 161, 89–105. [Google Scholar] [CrossRef]
  9. Pal, S.; Sasmal, S.K.; Pal, N. Chaos control in a discrete-time predator–prey model with weak Allee effect. Int. J. Biomath. 2018, 11, 1850089. [Google Scholar] [CrossRef]
  10. Huang, C.; Zhang, H.; Cao, J.; Hu, H. Stability and Hopf Bifurcation of a Delayed Prey–Predator Model with Disease in the Predator. Int. J. Bifurc. Chaos 2019, 29, 1950091. [Google Scholar] [CrossRef]
  11. Keller, E.F.; Segel, L.A. Traveling bands of chemotactic bacteria: A theoretical analysis. J. Theor. Biol. 1971, 30, 235–248. [Google Scholar] [CrossRef]
  12. Ahn, J.; Yoon, C. Global well-posedness and stability of constant equilibria in parabolic–elliptic chemotaxis system without gradient sensing. Nonlinearity 2019, 32, 1327–1351. [Google Scholar] [CrossRef]
  13. Cho, E.; Kim, Y.J. Starvation driven diffusion as a survival strategy of biological organisms. Bull. Math. Biol. 2013, 75, 845–870. [Google Scholar] [CrossRef] [PubMed]
  14. Tao, Y.; Winkler, M. Effects of signal-dependent motilities in a keller–segel-type reactiondiffusion system. Math. Models Methods Appl. Sci. 2017, 27, 1645–1683. [Google Scholar] [CrossRef] [Green Version]
  15. Yoon, C.; Kim, Y.J. Global existence and aggregation in a keller–segel model with fokker–Planck diffusion. Acta Appl. Math. 2016, 149, 101–123. [Google Scholar] [CrossRef]
  16. Szymańska, Z.; Rodrigo, C.; Lachowicz, M.; Chaplain, M. Mathematical modelling of cancer invasion of tissue: The role and effect of nonlocal interactions. Math. Model. Methods Appl. Sci. 2009, 19, 257–281. [Google Scholar] [CrossRef] [Green Version]
  17. Li, T.; Viglialoro, G. Boundedness for a nonlocal reaction chemotaxis model even in the attraction-dominated regime. Differ. Integral Equ. 2021, 34, 315–336. [Google Scholar]
  18. Shahid, A.; Huang, H.; Bhatti, M.M.; Zhang, L.; Ellahi, R. Numerical Investigation on the Swimming of Gyrotactic Microorganisms in Nanofluids through Porous Medium over a Stretched Surface. Mathematics 2020, 8, 380. [Google Scholar] [CrossRef] [Green Version]
  19. Ellahi, R.; Hussain, F.; Ishtiaq, F.; Hussain, A. Peristaltic transport of Jeffrey fluid in a rectangular duct through a porous medium under the effect of partial slip: An application to upgrade industrial sieves/filters. Pramana-J. Phys. 2019, 93, 34. [Google Scholar] [CrossRef]
  20. Bhatti, M.; Zeeshan, A.; Ellahi, R.; Anwar Bég, O.; Kadir, A. Effects of coagulation on the two-phase peristaltic pumping of magnetized prandtl biofluid through an endoscopic annular geometry containing a porous medium. Chin. J. Phys. 2019, 58, 222–223. [Google Scholar] [CrossRef]
  21. Frassu, S.; Viglialoro, G. Boundedness in a chemotaxis system with consumed chemoattractant and produced chemorepellent. Nonlinear Anal. 2021, 123, 112505. [Google Scholar] [CrossRef]
  22. Li, T.; Pintus, N.; Viglialoro, G. Properties of solutions to porous medium problems with different sources and boundary conditions. Z. Angew. Math. Phys. 2019, 70, 86. [Google Scholar] [CrossRef] [Green Version]
  23. Zhang, Y. Degenerate Diffusions with Advection. Ph.D. Thesis, University of California, Oakland, CA, USA, 2019. [Google Scholar]
  24. Pao, C. Nonlinear Parabolic and Elliptic Equations; Springer Science + Bussiness Media: Secaucus, NJ, USA, 2012. [Google Scholar]
  25. De Pablo, A.; Vázquez, J.L. Travelling waves and finite propagation in a reaction-diffusion Equation. J. Differ. Equ. 1991, 93, 19–61. [Google Scholar] [CrossRef] [Green Version]
  26. Philip, J.R. The theory of infiltration. Soil Sci. 1957, 83, 345–357. [Google Scholar] [CrossRef]
  27. De Pablo, A. Estudio de una Ecuación de Reacción—Difusión. Ph.D. Thesis, Universidad Autónoma de Madrid, Madrid, Spain, 1989. [Google Scholar]
  28. Vázquez, J.L. The Porous Medium Equation, Mathematical Theory; Oxford Mathematical Monographs: Oxford, UK, 2006. [Google Scholar]
  29. Zhao, S.; Ruan, S. Time periodic traveling wave solutions for periodic advection–reaction–diffusion systems. J. Differ. Equ. 2014, 257, 1078–1147. [Google Scholar] [CrossRef]
  30. Fenichel, N. Persistence and smoothness of invariant manifolds for flows. Indiana Univ. Math. J. 1971, 21, 193–226. [Google Scholar] [CrossRef]
  31. Akveld, M.E.; Hulshof, J. Travelling Wave Solutions of a Fourth-Order Semilinear Diffusion Equation. Appl. Math. Lett. 1998, 11, 115–120. [Google Scholar] [CrossRef] [Green Version]
  32. Jones, C.K.R.T. Geometric Singular Perturbation Theory in Dynamical Systems; Springer: Berlín, Germany, 1995. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Díaz Palencia, J.L.; Roa González, J.; Rahman, S.U.; Redondo, A.N. Regularity, Asymptotic Solutions and Travelling Waves Analysis in a Porous Medium System to Model the Interaction between Invasive and Invaded Species. Mathematics 2022, 10, 1186. https://doi.org/10.3390/math10071186

AMA Style

Díaz Palencia JL, Roa González J, Rahman SU, Redondo AN. Regularity, Asymptotic Solutions and Travelling Waves Analysis in a Porous Medium System to Model the Interaction between Invasive and Invaded Species. Mathematics. 2022; 10(7):1186. https://doi.org/10.3390/math10071186

Chicago/Turabian Style

Díaz Palencia, José Luis, Julián Roa González, Saeed Ur Rahman, and Antonio Naranjo Redondo. 2022. "Regularity, Asymptotic Solutions and Travelling Waves Analysis in a Porous Medium System to Model the Interaction between Invasive and Invaded Species" Mathematics 10, no. 7: 1186. https://doi.org/10.3390/math10071186

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop