Next Article in Journal
A Family of Analytic and Bi-Univalent Functions Associated with Srivastava-Attiya Operator
Previous Article in Journal
The Loss of Symmetry in Unilateral Bony Syngnathia: Case Report and Literature Review
Previous Article in Special Issue
Analysis of Solutions to a Free Boundary Problem with a Nonlinear Gradient Absorption
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetries and Solutions for Some Classes of Advective Reaction–Diffusion Systems

Dipartimento di Matematica e Informatica, Università Degli Studi di Catania, Viale Andrea Doria, 6, 95125 Catania, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(10), 2009; https://doi.org/10.3390/sym14102009
Submission received: 5 September 2022 / Revised: 19 September 2022 / Accepted: 23 September 2022 / Published: 25 September 2022
(This article belongs to the Special Issue Symmetry and Biomathematics: Recent Developments and Challenges)

Abstract

:
In this paper, we consider some reaction–advection–diffusion systems in order to obtain exact solutions via a symmetry approach. We write the determining system of a general class. Then, for particular subclasses, we obtain special forms of the arbitrary constitutive parameters that allow us to extend the principal Lie algebra. In some cases, we write the corresponding reduced system and we find special exact solutions.

1. Introduction

In several of our previous papers, we analysed, in the framework of symmetry methods, systems belonging to some subclasses of the following pde systems
u t = D x ( f ( u , v ) u x ) + g ( u , v , u x ) , v t = h ( u , v )
with them in mind, generally, regarding models of population dynamics. With D x , we denote the total differentiation with respect to x, while, usually, the subscripts denote the partial differentiations with respect to the indicated variables. Finally, in the following, the prime will denote the differentiation of functions with respect to their single argument.
The first equation, an advection–reaction–diffusion equation, is obtained by introducing a Fick–Fourier flux in the balance equation for the species of density u; the diffusion coefficient f could depend not only on u but also on density v [1,2], of a second species interacting with u. The function g may not always depend on the advection u x if the species u does not detect external stimuli (usually water currents or wind); this occurs in the swarming cells in the Proteus mirabilis colonies.
The second equation is concerned with the species of density v that is assumed to be subject to neither diffusion nor advection, as, for instance, Proteus mirabilis [2,3,4,5] or mosquito models (Aedes Aegypt [6,7,8,9], Anopheles [10,11]), respectively, for swinging cells or the so-called aquatic populations.
Here, we restrict ourselves to the subclass of phenomena modelled by the following partial differential equation system:
u t = D x ( f ( u ) u x ) + g ( u , v , u x ) , v t = h ( u , v ) ,
where, in general, there is a population of mosquitos of density u interacting in more or less a weak manner with an aquatic population of density v. Moreover, analogous systems could be found in several other fields, for instance, in the framework of the models dealing with chemotaxis [12,13]. The motivation of the Lie symmetry analysis approach is twofold. The first one, well known, is that the Lie symmetries offer a methodological way to derive exact solutions by reducing the number of dependent variables. The second one, not negligible, is that the discussion of the equation determining the symmetries could bring about requests for special forms of some constitutive functions in order to admit additional symmetries. This request, sometimes, creates an interesting feedback between the different constitutive parameters of the model and its sensitivity.
In (2), the arbitrary functions f ( u ) , g ( u , v , u x ) , and h ( u , v ) appear. Once we have obtained the so-called determining system, we will obtain some restrictions for the coordinates of the symmetry generator, and three classifying equations involving the unknown coordinates and the arbitrary elements of the system. One equation involves only the function f and another one only the function h. The last condition is more complex, because it involves both f and g, having each different functional dependence. Therefore, we must introduce some simplifying constitutive assumptions. In the following, we consider some systems of the subclass (2) where, for the sake of simplicity, we assume
f ( u ) f 0 , f 0 u ,
and
g = A ( u ) u x + Π ( u ) + Γ ( v ) ,
with
A ( u ) A 0 , A 0 u , Π ( u ) γ 1 u ( γ 2 u ) , γ 1 u ( γ 2 ln u ) ,
while the constitutive functions Γ ( v ) and h ( u , v ) remain free to perform a suitable classification with respect to them, which allows us to obtain additional symmetries. Once known, the classifying equations of (2) are obtained with arbitrary f, g, and h; by introducing some specializations of them, it is possible to ascertain the existence of extensions with respect to the principal Lie algebra L P [14] for some special form of the functions Γ ( v ) and h ( u , v ) .
In this work, we analyse three systems of the type (2) and several classes of exact solutions are found. To the best of our knowledge, the literature regarding group analysis applications to similar systems is quite poor, so we believe that our results are new. The outline of the paper is the following. In Section 2, we write the determining system for the class (2) and we establish the principal Lie algebra. Systems with linear advection are studied in Section 3, while, in Section 4, we consider a system with diffusion and advection coefficients linear in u. Some remarks about the constitutive functions Γ and h are given in Section 5. Finally, in the last section, the conclusions are shown.

2. Symmetries: Classifying Equations

We look for a symmetry generator of the form
X = ξ 1 ( x , t , u , v ) x + ξ 2 ( x , t , u , v ) t + η 1 ( x , t , u , v ) u + η 2 ( x , t , u , v ) v ,
that leaves the system (2) invariant. Then, we require that
X ( 2 ) u t + D x ( f ( u ) u x ) + g ( u , v , u x ) | ( 2 ) = 0 , X ( 2 ) v t + h ( u , v ) | ( 2 ) = 0 ,
where the operator X ( 2 ) is
X ( 2 ) = X + ζ t 1 u t + ζ x 1 u x + ζ x x 1 u x x + ζ t 2 v t .
As usual, the expressions of the coordinates ζ t 1 , ζ x 1 , ζ x x 1 , ζ t 2 are (see, e.g., [15,16,17,18,19,20,21,22])
ζ t 1 = D t ( η 1 ) u t D t ( ξ 1 ) u x D t ( ξ 2 ) ,
ζ x 1 = D x ( η 1 ) u t D x ( ξ 1 ) u x D x ( ξ 2 ) ,
ζ x x 1 = D x ( ζ x 1 ) u x t D x ( ξ 1 ) u x x D x ( ξ 2 ) ,
ζ t 2 = D t ( η 2 ) v t D t ( ξ 1 ) v x D t ( ξ 2 ) ,
where the operators D x , and D t are the total derivatives with respect to x and t, respectively.
The conditions (4) give us the system of determining equations for the unknown coordinates ξ 1 , ξ 2 , η 1 , and η 2 . By solving only the equations that do not depend on the form of the arbitrary elements f, g, and h, we obtain the following restrictions about the coordinates ξ 1 , ξ 2 , η 1 , and η 2 :
ξ 1 = α ( x ) , ξ 2 = β ( t ) , η 1 = ϕ ( x , t , u ) , η 2 = ψ ( x , t , v ) .
The generator (3) with the coordinates (10) will be a symmetry generator for the system (2) if the functions α , β , ϕ , and ψ satisfy the following remaining determining equations, where the arbitrary elements f, g, and h appear:
( 2 α β ) f ϕ f u = 0 , ( f g + u x 2 f 2 ) ϕ ϕ u u u x 2 ( α 2 ϕ x u ) u x + ϕ x x f 2 ψ g v f ϕ g u f + ( ϕ u 2 α ) g f ϕ t f ϕ u x 2 f f ( 2 ϕ x + ϕ u u x ) u x f f + ( α ϕ u ) u x ϕ x g u x f = 0 , ψ v β h ϕ h u ψ h v + ψ t = 0 .
It is a simple matter to derive that, for f, g, and h arbitrary, we obtain
α = 0 , β = 0 , ϕ = 0 , ψ = 0 .
Then, we can affirm that the principal Lie algebra L P is spanned by the time translations and space translations
X 1 = x , X 2 = t .

3. Symmetries and Solutions for Systems with Linear Advection

Here, we assume f = f 0 and g = A u x + Π ( u ) + Γ ( v ) , with f 0 0 and Γ 0 , so the system (2) becomes
u t = f 0 u x x + A u x + Π ( u ) + Γ ( v ) , v t = H ( u , v ) .
Remark 1.
We recall that the following change in independent and dependent variables
τ = t , z = x + A t , u = u ( τ , z ) v = v ( τ , z )
brings the system (13) in the following form
u τ = f 0 u z z + Π ( u ) + Γ ( v ) , v τ = A v z + H ( u , v ) .
We can find the symmetry generator for the class (13) by solving (11) after having input the special forms of the functions f, g, and h. Taking into account that, for the system (13),
f = f 0 , g = A u x + Π ( u ) + Γ ( v ) , h = H ( u , v ) ,
the conditions (11) become
( 2 α β ) f 0 = 0 , ϕ u u u x 2 ( α 2 ϕ x u ) u x + ϕ x x f 0 2 ψ Γ f 0 ϕ Π f 0 + ( ϕ u 2 α ) ( A u x + Π ( u ) + Γ ( v ) ) f 0 + ϕ t f 0 +
( α ϕ u ) u x ϕ x A f 0 = 0 ,
ψ v β H ϕ H u ψ H v + ψ t = 0 .
From (17), we have
α ( x ) = c 0 + c 1 x , β ( t ) = 2 c 1 t + c 2 ,
while, taking into account that no function depends on u x , from (18), we obtain
ϕ u u = 0 , c 1 A + 2 ϕ x u f 0 = 0 ,
ϕ x x f 0 + ψ Γ + ϕ Π + ( 2 c 1 ϕ u ) ( Π + Γ ) ϕ t + ϕ x A = 0 .
From conditions (21), it follows that
ϕ ( x , t , u ) = u ϕ 1 ( t ) c 1 A x 2 f 0 + ϕ 2 ( x , t ) .
In the condition (22), the arbitrary elements Π ( u ) and Γ ( v ) appear.
We study (13) by specializing Π ( u ) in these two special forms
Π ( u ) = γ 1 u ( γ 2 u ) , Π ( u ) = γ 1 u ( γ 2 ln u ) .

3.1. Logistic u- Production

A typical choice for the function Π ( u ) is given by the logistic form
Π ( u ) = γ 1 u ( γ 2 u ) ,
but taking into account (23) and this form of Π ( u ) , from the condition (22), we obtain
c 1 = 0 , ϕ 1 ( t ) = 0 , ϕ 2 ( x , t ) = 0 , ψ ( x , t , v ) = 0 ,
i.e., extensions of the principal Lie algebra do not exist.
Thus, we look for invariant solutions of the system
u t = f 0 u x x + A u x + γ 1 u ( γ 2 u ) + Γ ( v ) , v t = H ( u , v ) .
By using the generator
X = c 1 X 1 + c 2 X 2 ,
the solutions are
u = U ( σ ) , v = V ( σ ) ,
where σ = c 2 x c 1 t , and with U and V solutions of the following reduced system
c 1 U + f 0 U + c 2 A U + γ 1 U ( γ 2 U ) + Γ ( V ) = 0 , c V = H ( U , V ) .
If we specialize
Γ ( v ) = γ 2 v , h ( u , v ) = γ 2 v γ 3 γ 2 ( u 2 γ 0 ) ,
the previous reduced system becomes
c 1 U + f 0 U + c 2 A U + γ 1 U ( γ 2 U ) γ 2 V = 0 , c V = ( V γ 3 γ 2 ) ( U 2 γ 0 ) .
If U = 2 γ 0 , V must be constant too and we obtain V = 2 γ 0 γ 1 γ 2 2 γ 0 γ 2 .
Instead, if V = γ 3 γ 2 , then U must satisfy the following autonomous equation:
f 0 U + ( c 1 + c 2 A ) U + γ 1 U ( γ 2 U ) + γ 3 = 0 .
After having performed the substitution Z ( U ) = U , equation (30) becomes the following Abel equation of the second kind in the canonical form [23,24]:
f 0 Z Z + ( c 1 + c 2 A ) Z + γ 1 U ( γ 2 U ) + γ 3 = 0 .
From this, after a suitable choice of arbitrary constants c 1 and c 2 , we obtain the following separable variable equation:
f 0 Z Z + γ 1 U ( γ 2 U ) + γ 3 = 0 .

3.2. L o g -Logistic u- Production

If we want to obtain some extensions of the principal Lie algebra but at the same time a form of Π ( u ) with the characteristics of the logistic form, we can choose
Π ( u ) = γ 1 u ( γ 2 ln u ) .
Then, we consider
u t = f 0 u x x + A u x + γ 1 u ( γ 2 ln u ) + Γ ( v ) , v t = H ( u , v ) .
We recall that for the system (13), we have already obtained
α ( x ) = c 0 + c 1 x , β ( t ) = 2 c 1 t + c 2 , ϕ ( x , t , u ) = u ϕ 1 ( t ) c 1 A x 2 f 0 + ϕ 2 ( x , t ) ,
and the following conditions
( A c 1 x + 4 c 1 f 0 2 f 0 ϕ 1 ) ( Π + Γ ) + ( 2 ϕ 2 f 0 + 2 ϕ 1 u f 0 A u c 1 x ) Π +
2 A ϕ 2 x f 0 ( A 2 c 1 + 2 ϕ 1 f 0 ) u + ϕ 2 x x f 0 2 + 2 ψ Γ f 0 2 ϕ 2 t f 0 = 0 ,
ψ v β H u ϕ 1 c 1 A x 2 f 0 + ϕ 2 H u ψ H v + ψ t = 0 .
Substituting (33), the condition (36) reads
( A c 1 γ 1 x + 4 c 1 γ 1 γ 2 f 0 A 2 c 1 2 f 0 γ 1 ϕ 1 2 ϕ 1 f 0 ) U 4 c 1 f 0 γ 1 u ln u + 2 ϕ 2 f 0 γ 1 ln u + A c 1 x Γ + 2 ϕ 2 f 0 γ 1 γ 2 + 2 A ϕ 2 x f 0 + 2 ϕ 2 x x f 0 2 + 2 ϕ 2 f 0 γ 1 2 ϕ 1 Γ f 0 + 2 ψ Γ f 0 + 4 c 1 f 0 Γ 2 ϕ 2 t f 0 = 0 .
Taking into account that in the condition (38), no function depends on u, we obtain
c 1 = 0 , γ 1 ϕ 1 ϕ 1 = 0 , ϕ 2 = 0 , ϕ 1 Γ f 0 ψ Γ f 0 + ϕ 2 t f 0 = 0 ,
i.e.,
ϕ 1 = c 3 e γ 1 t , ψ = c 3 e γ 1 Γ Γ .
Finally, by substituting these results into the condition (37), it becomes
c 3 e γ 1 t ( H u u H ) Γ 2 + ( H v + γ 1 ) Γ Γ + h Γ Γ = 0 .
We obtain the following extension of L P
X 3 = e γ 1 t u u + Γ Γ v
only when
H ( u , v ) = Φ ( ω ) γ 1 ln Γ Γ Γ
with ω = u Γ .
Remark 2.
We obtain the same result even if A = 0 , i.e., without advection.
If we use the generator
X = c 1 X 2 + c 2 X 1 + c 3 X 3 ,
with c 1 0 , the invariant solutions are
u ( t , x ) = e c 3 c 1 γ 1 e γ 1 t U ( σ ) , Γ ( v ( t , x ) ) ) = e c 3 c 1 γ 1 e γ 1 t V ( σ ) ,
where σ = c 1 x c 2 t . The functions U and V are solutions of the reduced system
c 2 U = c 1 2 f 0 U + c 1 A U + γ 1 U ( γ 2 ln U ) + V , c 2 V = Φ U V γ 1 ln V V .
In order to have a reduced system of the first order, we use the generator
X = c 2 X 1 + c 3 X 3 ,
with c 2 0 . In this case, the invariant solutions are
u ( t , x ) = e c 3 c 2 x e γ 1 t U ( t ) , Γ ( v ( t , x ) ) ) = e c 3 c 2 x e γ 1 t V ( t ) .
The functions U and V are solutions of the reduced system
U = c 3 2 c 2 2 e 2 γ 1 t f 0 U + c 3 c 2 e γ 1 t A U + γ 1 U ( γ 2 ln U ) + V , V = Φ U V γ 1 ln V V .
In this case, if we consider
Φ U V = p 0 + p 1 V U + p 2 V U + p 3 ln U V ,
we can look for solutions of type
U ( t ) = u 0 e f 1 ( t ) a n d V ( t ) = v 0 e 2 γ 1 e f 1 ( t ) .
If we take
p 0 = γ 1 ( γ 2 2 ) , p 1 = f 0 u 0 2 p 2 2 A 2 v 0 2 1 , p 3 = γ 1 ,
by choosing
c 2 = A v 0 c 3 ( u 0 ) p 2 ,
we obtain
f 1 ( t ) = γ 2 ln ( u 0 ) + u 0 v 0 p 2 t + k e γ 1 t A 2 v 0 2 + f 0 p 2 2 u 0 2 A 2 γ 1 u 0 v 0 e 2 γ 1 t ,
with k arbitrary constant. Once we have obtained U ( t ) and V ( t ) , from (48), we obtain directly u ( t , x ) . To find v ( t , x ) , we must invert the assigned function Γ (numerically or analytically).

4. Symmetries and Solutions for Systems with Non-Linear Advection

We consider here the following class A-RD systems
u t = ( f 0 u u x ) x + A u u x + γ 1 u ( γ 2 u ) + Γ ( v ) , v t = H ( u , v ) ,
with f 0 0 and Γ 0 . In this case, we have
f = f 0 u , g = A u u x + γ 1 u ( γ 2 u ) + Γ ( v ) , h = H ( u , v ) .
As in the previous case, we can find the symmetry generator for the class (55) by solving (11) after having input the special forms of the functions f, g, and h, i.e.,
( 2 α β ) f 0 u ϕ f 0 = 0 , ( f 0 g + u x 2 f 0 2 ) ϕ ϕ u u u x 2 ( α 2 ϕ x u ) u x + ϕ x x f 0 2 u 2 ψ Γ f 0 u + ϕ ( A u x + γ 1 γ 2 2 γ 1 γ 2 u ) f 0 u ϕ t f 0 u ( 2 ϕ x + ϕ u u x ) u x f 0 2 u + ( ϕ u 2 α ) A u u x + γ 1 u ( γ 2 u ) + Γ f 0 u +
( α ϕ u ) u x ϕ x A f 0 u 2 = 0 ,
ψ v β H ϕ H u ψ H v + ψ t = 0 .
For arbitrary H ( u , v ) , we obtain only the principal Lie algebra
X 1 = t , X 2 = x .
For
H ( u , v ) = 2 γ 1 γ 2 + Γ Φ ( ω ) Γ Γ
with ω = u Γ , we obtain also the following generator
X 3 = e γ 1 γ 2 t t + γ 1 γ 2 u u + 2 γ 1 γ 2 Γ Γ v .
By using the generator
X = c 1 X 3 + c 2 X 2 ,
the invariant solutions are
u ( t , x ) = e γ 1 γ 2 t U ( σ ) , Γ ( v ( t , x ) ) ) = e 2 γ 1 γ 2 t V ( σ )
where σ = c 2 γ 1 γ 2 e γ 1 γ 2 t c 1 x . The functions U and V are solutions of the reduced system
c 2 U = c 1 2 f 0 ( U 2 + U U ) c 1 A U U γ 1 U 2 + V , c 2 V = V 3 2 Φ U V .
If we choose
Φ ( ω ) = 2 γ 1 A 2 A 2 + 4 f 0 γ 1 1 γ 1 ω 2
we find the following particular solutions for system (65)
U ( σ ) = 2 c 2 ( 2 γ 1 σ c 1 A ) c 1 2 ( A 2 + 4 f 0 γ 1 ) ,
V ( σ ) = 16 γ 1 3 c 2 2 σ 2 c 1 4 ( A 2 + 4 f 0 γ 1 ) 2 .
From (64), we obtain directly u ( t , x ) . To find v ( t , x ) , we must invert the function Γ (numerically or analytically).

5. Some Remarks about the Constitutive Functions Γ and H

We stress that, in all systems that we have considered, the extension of L P does not imply a constraint on Γ form, while it influences the functional dependence of H that must be
H = H ( Γ , Φ ( ω ) ) ,
where Φ is an arbitrary function of ω = ω ( u , Γ ) , being ω a specific function of its arguments. After having recalled that Γ ( v ) 0 , we consider the following equivalence transformation:
w = Γ ( v ) .
The system (34) with H given by (43), after (69), reads as follows:
u t = f 0 u x x + A u x + γ 1 u ( γ 2 ln u ) + w , w t = Φ u w γ 1 ln w w .
A symmetry generator for the system (70) has the form
X = ξ 1 x + ξ 2 t + η 1 u + η * w ,
where the infinitesimal η * is linked to the infinitesimal η 2 by the relation
η * = Γ η 2
obtained by requiring the invariance of the transformation (69) with respect to the generator (see, e.g., [25])
X * = ξ 1 x + ξ 2 t + η 1 u + η 2 v + η * w .
and then its extension of L P is
X 3 = e γ 1 t u u + w w .
The system (55) with H ( u , v ) given by (61) can be written in the following equivalent form:
u t = ( f 0 u u x ) x + A u u x + γ 1 u ( γ 2 u ) + w w t = 2 γ 1 γ 2 + w Φ u w w .
Analogously, for the system (75), the extension of L P assumes the following form:
X 3 = e γ 1 γ 2 t t + γ 1 γ 2 u u + 2 γ 1 γ 2 w w .
Actually, the new forms (70), (75) of the systems (34) and (55) could be more useful in the search for additional exact solutions, once we have fixed the form of arbitrary constitutive function Φ ( ω ) , or they could suggest new ways to solve the systems. However, we must not forget that when we return to the original biological variable v, we must fix the form of the still arbitrary function Γ ( v ) .

6. Conclusions

In this paper, we continue our studies dealing with the extension of the principal Lie algebra and with our search for solutions of some subclasses of systems (2), begun in [7] and recently continued in [26,27]. In all the cases analysed here, a priori, we have considered the two constitutive functions of production Γ ( v ) and h ( u , v ) arbitrary in order to obtain a classification with respect to each of them.
In the case of linear advection, we obtain, for Π ( u ) = γ 1 u ( γ 2 ln u ) , the extension by one of L P that leaves Γ ( v ) arbitrary while h is, as expected as a solution of a PDE, an arbitrary function depending on ω = u Γ . For special forms of h, wide classes of exact solutions are derived. We remark that, in this case, the function Π ( u ) has the same numerical features, at least in a suitable interval, of the classical logistic form Π ( u ) = γ 1 u ( γ 2 u ) for which we do not obtain extensions. For this form of Π ( u ) , however, looking for travelling wave solutions, in some cases, we bring the reduced system, with appropriate transformations, to an Abel equation. For the same system without advection, we obtained additional extensions in [27].
In the system analysed in Section 4, we assumed diffusion and advection coefficients both linear in u. Additionally, in this case, we obtain extensions with Γ arbitrary, while H is still dependent on an arbitrary function Φ ( ω ) with ω = u Γ . For an appropriate choice of Φ ( ω ) , classes of exact solutions are derived.

Author Contributions

M.T. and R.T. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

M.T. performed this study within the framework of G.N.F.M. of INdAM. R. Traciná acknowledges Universitá degli Studi di Catania, Piano della Ricerca 2020/2022 Linea di intervento 2 “QICT”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Medvedev, G.S.; Kaper, T.J.; Kopell, N. A reaction diffusion system with periodic front Dynamics. SIAM J. Appl. Math. 2000, 60, 1601–1638. [Google Scholar] [CrossRef]
  2. Czirok, A.; Matsushita, M.; Vicsek, T. Theory of periodic swarming of batteria: Application to Proteus mirabilis. Phys. Rev. E 2001, 63, 031915. [Google Scholar] [CrossRef]
  3. Rauprich, O.; Matsushita, M.; Weijer, K.; Siegert, F.; Esipov, S.E.; Shapiro, J.A. Periodic phenomena in Proteus mirabilis swarm colony development. J. Bacteriol. 1996, 178, 6525–6538. [Google Scholar] [CrossRef]
  4. Esipov, S.E.; Shapiro, J.A. Kinetic model of Proteus mirabilis swarm colony development. J. Math. Biol. 1998, 36, 249–268. [Google Scholar] [CrossRef]
  5. Orhan, Ö.; Torrisi, M.; Tracinà, R. Group methods applied to a reaction-diffusion system generalizing Proteus Mirabilis models. Commun. Nonlinear Sci. Numer. Simulat. 2019, 70, 223–233. [Google Scholar] [CrossRef]
  6. Takahashi, L.T.; Maidana, N.A.; Ferreira, W.C., Jr.; Pulino, P.; Yang, H.M. Mathematical models for the Aedes aegypti dispersal dynamics: Travelling waves by wing and wind. Bull. Math. Biol. 2005, 67, 509–528. [Google Scholar] [CrossRef] [PubMed]
  7. Freire I., L.; Torrisi, M. Symmetry methods in mathematical modeling of Aedes aegypti dispersal dynamics. Nonlinear Anal. Real World Appl. 2013, 14, 1300–1307. [Google Scholar] [CrossRef]
  8. Zhang, M.; Lin, Z. A reaction-diffusion-advection model for Aedes aegypti mosquitoes in a time-periodic environment. Nonlinear Anal. Real World Appl. 2019, 46, 219–237. [Google Scholar] [CrossRef]
  9. Liu, Y.; Jiao, F.; Hu, L. Modeling mosquito population control by a coupled system. J. Math. Anal. Appl. 2022, 506, 125671. [Google Scholar] [CrossRef]
  10. Anguelov, R.; Dumont, Y.; Lubuma, J.-S. Mathematical modeling of sterile insect technology for control of anopheles mosquito. Comput. Math. Appl. 2012, 64, 374–389. [Google Scholar] [CrossRef]
  11. Olaniyi, S.; Obabiyi, O.S. Mathematical model for malaria transmission dynamics in human and mosquito populations with non linear force of infection. Int. J. Pure Appl. Math. 2013, 88, 125–156. [Google Scholar] [CrossRef]
  12. Tello, J.I.; Winkler, M. Stabilization in a two-species chemotaxis system with logistic source. Nonlinearity 2012, 25, 1413–1425. [Google Scholar] [CrossRef]
  13. Cherniha, R.; Didovych, M. Exact solutions of the simplified Keller-Segel model. Commun. Nonlinear Sci. Numer. Simulat. 2013, 18, 2960–2971. [Google Scholar] [CrossRef]
  14. Akhatov, I.S.; Gazizov, R.K.; Ibragimov, N.H. Nonlocal symmetries: A heuristic approach. J. Soviet Math. 1991, 55, 1401–1450. [Google Scholar] [CrossRef]
  15. Bluman, G.W.; Cole, J.D. Similarity Methods for Differential Equations; Springer: New York, NY, USA, 1974. [Google Scholar]
  16. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  17. Ibragimov, N.H. Transformation Groups Applied to Mathematical Physics; Reidel: Dordrecht, The Netherlands, 1985. [Google Scholar]
  18. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1986. [Google Scholar]
  19. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  20. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  21. Cantwell, B.J. Introduction to Symmetry Analysis; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  22. Ibragimov, N.H. A Practical Course in Differential Equations and Mathematical Modelling; World Scientific Publishing Co Pvt Ltd.: Singapore, 2009. [Google Scholar]
  23. Zaitsev, V.F.; Polyanin, A.D. Discrete-Group Methods for Integrating Equations of Nonlinear Mechanics; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
  24. Polyanin, A.D.; Zaitsev, V.F. Handbook of Exact Solutions for Ordinary Differential Equations, 2nd ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  25. Torrisi, M.; Tracinà, R.; Valenti, A. On the linearization of semilinear wave equations. Nonlinear Dyn. 2004, 36, 97–106. [Google Scholar] [CrossRef]
  26. Torrisi, M.; Tracinà, R. An Application of Equivalence Transformations to Reaction Diffusion Equations. Symmetry 2015, 7, 1929–1944. [Google Scholar] [CrossRef]
  27. Torrisi, M.; Tracinà, R. Lie symmetries and solutions of reaction diffusion systems arising in biomathematics. Symmetry 2021, 13, 1530. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Torrisi, M.; Tracinà, R. Symmetries and Solutions for Some Classes of Advective Reaction–Diffusion Systems. Symmetry 2022, 14, 2009. https://doi.org/10.3390/sym14102009

AMA Style

Torrisi M, Tracinà R. Symmetries and Solutions for Some Classes of Advective Reaction–Diffusion Systems. Symmetry. 2022; 14(10):2009. https://doi.org/10.3390/sym14102009

Chicago/Turabian Style

Torrisi, Mariano, and Rita Tracinà. 2022. "Symmetries and Solutions for Some Classes of Advective Reaction–Diffusion Systems" Symmetry 14, no. 10: 2009. https://doi.org/10.3390/sym14102009

APA Style

Torrisi, M., & Tracinà, R. (2022). Symmetries and Solutions for Some Classes of Advective Reaction–Diffusion Systems. Symmetry, 14(10), 2009. https://doi.org/10.3390/sym14102009

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