Next Article in Journal
New Security Development and Trends to Secure the SCADA Sensors Automated Transmission during Critical Sessions
Next Article in Special Issue
Centrally Extended Conformal Galilei Algebras and Invariant Nonlinear PDEs
Previous Article in Journal
Generation of Supramolecular Chirality around Twofold Rotational or Helical Axes in Crystalline Assemblies of Achiral Components
Previous Article in Special Issue
Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Application of Equivalence Transformations to Reaction Diffusion Equations

Department of Mathematics and Computer Sciences, University of Catania, Catania I 95125, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2015, 7(4), 1929-1944; https://doi.org/10.3390/sym7041929
Submission received: 15 July 2015 / Revised: 13 October 2015 / Accepted: 15 October 2015 / Published: 23 October 2015
(This article belongs to the Special Issue Lie Theory and Its Applications)

Abstract

:
In this paper, we consider a quite general class of advection reaction diffusion systems. By using an equivalence generator, derived in a previous paper, the authors apply a projection theorem to determine some special forms of the constitutive functions that allow the extension by one of the two-dimensional principal Lie algebra. As an example, a special case is discussed at the end of the paper.

1. Introduction

In this paper, we focus our attention on the following family of 2 × 2 nonlinear advection reaction diffusion systems in ( 1 + 1 ) independent variables:
u t = ( f ( u ) u x ) x + g ( u , v , u x ) , v t = h ( u , v ) ,
with f ( u ) , g ( u , v , u x ) , h ( u , v ) analytic functions. These systems, apart from their own mathematical interest, offer the possibility to be analyzed as possible biomathematical models for two interacting species u and v, where one of them, the species v, does not suffer diffusion. The dependence of the function g on the gradient u x shows advection effects; in fact, the individuals of the species u could be influenced by external stimuli as wind velocity or water currents. Of course, the absence of the advective phenomena brings to the following system:
u t = ( f ( u ) u x ) x + g ( u , v ) , v t = h ( u , v ) ,
that can describe the evolution of the Aedes aegypti mosquito population in a region where wind effects are negligible or the evolution of a Proteus mirabilis bacterial colony when the diffusion coefficient depends only on the species u, that is when the system (2) is a subclass of the following wider class:
u t = ( f ( u , v ) u x ) x + g ( u , v ) , v t = h ( u , v ) ,
considered in [1,2,3,4].
One of the most important problems in modeling the phenomena of life sciences and natural sciences is to select “good” forms of arbitrary functions (constitutive equations) that fit well with the experimental data and possess mathematical properties that allow scientists to get some solutions or much news about them.
A powerful tool of investigation in this field is given from transformation groups, in particular from equivalence transformations and symmetries.
In the framework of the group analysis, the literature concerning the systems of the type (1) is scarce. There are no papers devoted to a complete Lie symmetry analysis of PDE systems with advection (convection) terms of the form (1). In [5], it is possible to find a complete solution of this problem for a class of diffusion systems with convection terms in both equations. Moreover, the paper [6] contains some description of Lie symmetries for a class of systems, which includes cases having a structure similar to system (3). However, it is possible to find some papers devoted to the complete Lie symmetry analysis of a single advection (convection) reaction diffusion equation (see, e.g., [7,8,9]).
Following, e.g., [10], an equivalence transformation for the system (1) is a non-degenerate change of the independent and dependent variables t , x , u , v into t ^ , x ^ , u ^ , v ^ :
x = x ( x ^ , t ^ , u ^ , v ^ ) , t = t ( x ^ , t ^ , u ^ , v ^ ) , u = u ( x ^ , t ^ , u ^ , v ^ ) , v = v ( x ^ , t ^ , u ^ , v ^ ) ,
that transforms a system of the class (1) in another one of the same class. That is, an equivalence transformation brings the system of the form (1) in a system preserving the differential structure, but, in general, with:
f ^ ( u ^ ) f ( u ) , g ^ ( u ^ , v ^ , u ^ x ^ ) g ( u , v , u x ) , h ^ ( u ^ , v ^ ) h ( u , v ) .
It maps a solution of a system in a solution of the equivalent system.
It could occur that the transformed equations show still the same structure, but the arbitrary functions are depending on additional variables. In this case, the equivalence is said weak.
Of course, in the case:
f ^ ( u ^ ) = f ( u ) , g ^ ( u ^ , v ^ , u ^ x ^ ) = g ( u , v , u x ) , h ^ ( u ^ , v ^ ) = h ( u , v ) ,
an equivalence transformation becomes a symmetry (a transformation of variables that leaves invariant the transformed system).
A symmetry allows one to reduce the number of independent variables of an equation so that, for instance, a PDE in 1 + 1 independent variables can become an ODE. Once solved this last one, going back to the original variables, we get a solution that is invariant with respect to the symmetries used for the reduction. It is worthwhile to note that a symmetry transforms invariant solutions into invariant solutions that are not essentially different (see Ovsiannikov [11]), but, having a different form, they could satisfy different suitable initial/boundary conditions.
The aim of this paper is an improvement of the results that we have shown in [12], bearing in mind some generalization of the special form assumed from the constitutive functions f, g and h already used in some previous papers about [4,12,13,14]. In this paper, we use the infinitesimal generator of equivalence transformations derived in [12] for the class (1) in order to obtain some extensions of the principal Lie algebra for the following subclass:
u t = ( f ( u ) u x ) x + u r u x + Γ 1 ( u ) + Γ 2 ( v ) , v t = h ( u , v ) .
Here, we assumed:
g ( u , v , u x ) = u r u x + Γ 1 ( u ) + Γ 2 ( v )
that is a generalization of that ones used in [12,13].
In the next section, after recalling, for the sake of completeness, some elements about equivalence transformations (for additional mathematical and methodical details, the interested reader can see [12,15,16]), we write the equivalence generator derived in [12]. In Section 3, the principal Lie algebra and its extensions are discussed; moreover, a simple example that could be related to the biomathematical model of Aedes aegypti is considered. The conclusions are given in the last section.

2. On Equivalence Transformations and Their Calculation for the Class (1)

It is easy to ascertain that, often, in papers on differential equations, it is possible to find several examples of equivalence transformations and their applications. In general, in the past and now, the direct search for the most general equivalence transformations through the finite form of the transformation has been used. Quite often, this search is connected to considerable computational difficulties and does not always lead to the complete solution of the problem (e.g., [17,18]).

2.1. Elements on Equivalence Transformations

Following [11,15,16,19,20] (see also, e.g., [10,21,22,23]), we look for the infinitesimal generator of the equivalence transformations of the system (1) of the form:
Y = ξ 1 x + ξ 2 t + η 1 u + η 2 v + μ 1 f + μ 2 g + μ 3 h
where the infinitesimal components ξ 1 , ξ 2 , η 1 , η 2 are sought depending on x, t, u, v, while the infinitesimal components μ i ( i = 1 , 2 , 3 ) can also depend on u t , u x , v t , v x , f, g and h. Here, we are interested in obtaining the infinitesimal coordinates ξ i , η i and μ j ( i = 1 , 2 and j = 1 , 2 , 3 ), by applying the Lie–Ovsiannikov infinitesimal criterion [11] by requiring the invariance, with respect to a suitable prolongations Y ( 1 ) and Y ( 2 ) of generator (9), of the following equations:
u t ( f u x ) x g = 0 ,
v t h = 0 ,
without requiring the invariance of the so-called auxiliary conditions [15,16,24,25]:
f t = f x = f v = f u x = f u t = f v x = f v t = g t = g x = g u t = g v t = g v x = 0 ,
h t = h x = h u x = h u t = h v x = h v t = 0 ,
that characterize the functional dependence of f, g and h.
In this way, we obtain the weak equivalence transformations [15,16].
The main difference with respect to the classical one is that the infinitesimal operators of weak equivalence transformations can generate transformations that do not preserve the functional dependence of the arbitrary elements.
With respect to the application in biomathematical models, equivalence and weak equivalence transformations were applied not only to study tumor models [26,27], but also the population dynamics in [1,3,4].

2.2. Calculation of Weak Equivalence Transformations

We need the following prolongations Y ( 1 ) and Y ( 2 ) :
Y ( 1 ) = Y + ζ 1 1 u x + ζ 2 1 u t + ζ 1 2 v x + ζ 2 2 v t + ω u 1 f u ,
Y ( 2 ) = Y ( 1 ) + ζ x x 1 u x x ,
with (see [12] for more details),
ζ 1 1 = D x η 1 u x D x ξ 1 u t D x ξ 2 ,
ζ 2 1 = D t η 1 u x D t ξ 1 u t D t ξ 2 ,
ζ 1 2 = D x η 2 v x D x ξ 1 v t D x ξ 2 ,
ζ 2 2 = D t η 2 v x D t ξ 1 v t D t ξ 2 ,
ζ 11 1 = D x ζ 1 1 u x x D x ξ 1 u t x D x ξ 2 ,
ω u 1 = D ˜ u μ 1 f u D ˜ u η 1 ,
where D x and D t are, respectively, the total derivatives with respect to x and t, while in our case, the operator D ˜ u is defined as:
D ˜ u = u + f u f + g u g + h u h .
The invariant conditions read:
ζ 2 1 2 ζ 1 1 u x f u u x 2 ω u 1 u x x μ 1 f ζ 11 1 μ 2 = 0 ,
ζ 2 2 μ 3 = 0 ,
both under the constraints (10) and (11).
Following the usual techniques, we derive the following infinitesimal components for the weak equivalence generators:
ξ 1 = α ( x ) , ξ 2 = β ( t ) , η 1 = δ ( t , u ) , η 2 = λ ( x , t , v ) ,
μ 1 = ( 2 α β ) f , μ 2 = δ t + ( δ u β ) g + ( α u x δ u u u x 2 ) f , μ 3 = ( λ v β ) h + λ t ,
where α ( x ) , β ( t ) , δ ( t , u ) , λ ( x , t , v ) are arbitrary real functions of their arguments. The corresponding infinitesimal generator is:
Y = α ( x ) x + β ( t ) t + δ ( t , u ) u + λ ( x , t , v ) v + ( 2 α β ) f f + δ t + ( δ u β ) g + ( α u x δ u u u x 2 ) f g + ( λ v β ) h + λ t h .

3. Symmetries for a Subclass of Advection Reaction Diffusion Systems

In this section, we apply the projection theorem, introduced in [28] and successively generalized in [15,16,24], in order to carry out a symmetry classification for the following subclass of system (1):
u t = ( f ( u ) u x ) x + u r u x + Γ 1 ( u ) + Γ 2 ( v ) , v t = h ( u , v ) ,
with r 0 , Γ 1 0 and Γ 2 0 .
For the system (28), we can affirm the following:
Theorem 1. The projection of the infinitesimal weak equivalence generator Y for the system (1) on the space ( x , t , u , v ) :
X = α ( x ) x + β ( t ) t + δ ( t , u ) u + λ ( x , t , v ) v
is an infinitesimal symmetry generator of a system of the class (28) if and only if the constitutive equations, specifying the forms of f, g and h, are invariant with respect to Y.
Applying the previous theorem, in order to obtain the determining system for the subclass (28), we request the invariance with respect to Y of the following constitutive equations:
f = f ( u ) , g = u r u x + Γ 1 ( u ) + Γ 2 ( v ) , h = h ( u , v ) ,
that is
Y ( f f ( u ) ) = 0 , Y ( g u r u x Γ 1 ( u ) Γ 2 ( v ) ) = 0 , Y ( h h ( u , v ) ) = 0 ,
under the constraints (30). Then, taking into account the form (27) of generator Y, we have the following determining equations:
( 2 α β ) f δ f u = 0 ,
δ t + ( δ u β ) ( u r u x + Γ 1 + Γ 2 ) ( δ u u u x 2 α u x ) f ( δ u α ) u x u r δ ( r u r 1 u x + Γ 1 ) λ Γ 2 = 0 ,
( λ v β ) h + λ t δ h u λ h v = 0 .
We recall here that the principal Lie algebra L P [10,19] is the Lie algebra that leaves invariant the system (28) for any form of the functions f ( u ) , Γ 1 ( u ) , Γ 2 ( v ) and h ( u , v ) . Then, the principal Lie algebra is the generator (29) where the functions α, β, δ and λ are solutions of the system (32)–(34) for arbitrary functions f ( u ) , Γ 1 ( u ) , Γ 2 ( v ) and h ( u , v ) . Consequently, it is a simple matter to ascertain the following:
Corollary 2. The projection (29) of infinitesimal weak equivalence generator Y for the system (1) on the space ( x , t , u , v ) is the infinitesimal symmetry generator corresponding to the principal Lie algebra of the class (28) if and only if η i = 0 , μ j = 0 , i = 1 , 2 , j = 1 , 2 , 3 .
Then, the principal Lie algebra L P is spanned by the following generators corresponding respectively to translations in time and in the space:
X 1 = t , X 2 = x .

3.1. On the Extensions of the L P

Here, we analyze some particular cases of the extension of the principal algebra for the class (28). That is, we look for a family of particular functions f ( u ) , Γ 1 ( u ) , Γ 2 ( v ) and h ( u , v ) , such that the solution of system (32)–(34) is different from δ = λ = 0 and α and β constants, which corresponds to the generators (35).
From Equation (32), deriving with respect to x, we get:
2 α f = 0 ,
that is
α ( x ) = a 1 x + a 0 ,
with a 0 and a 1 arbitrary constants. Consequently, from Equation (33), deriving with respect to x, we get:
λ x Γ 2 = 0 ,
and taking into account that Γ 2 0 , we obtain:
λ = λ ( t , v ) .
Equation (33) becomes:
δ t + ( δ u β ) ( Γ 1 + Γ 2 ) δ Γ 1 λ Γ 2 + u x ( ( a 1 β ) u r δ r u r 1 ) + u x 2 ( δ u u f ) = 0 .
As any function does not depend on u x , from Equation (40), we derive:
δ u u f = 0 ,
( a 1 β ) u r δ r u r 1 = 0 ,
δ t + ( δ u β ) ( Γ 1 + Γ 2 ) δ Γ 1 λ Γ 2 = 0 .
Then, from Equation (41):
δ ( t , u ) = u A 1 ( t ) + A 2 ( t ) ,
with A 1 ( t ) and A 2 ( t ) arbitrary functions of t. After these partial results, for the sake of clarity, we rewrite the determining system:
( 2 a 1 β ) f ( u A 1 + A 2 ) f u = 0 ,
u ( a 1 β r A 1 ) r A 2 = 0 ,
( u A 1 + A 2 ) + ( A 1 β ) ( Γ 1 + Γ 2 ) ( u A 1 + A 2 ) Γ 1 λ Γ 2 = 0 ,
( λ v β ) h + λ t ( u A 1 + A 2 ) h u λ h v = 0 .
From Equation (46), taking into account that any function does not depend on u, we get:
A 1 = a 1 β r , A 2 = 0 ,
then the other equations become:
( 2 a 1 β ) f u a 1 β r f u = 0 ,
β r u + a 1 ( 1 + r ) β r ( Γ 1 + Γ 2 ) u a 1 β r Γ 1 λ Γ 2 = 0 ,
( λ v β ) h + λ t u a 1 β r h u λ h v = 0 .
We observe that from Equation (50), if f is arbitrary, it follows β = b 0 , a 1 = 0 , from Equation (51) λ = 0 , while the Equation (52) is satisfied. Therefore, for f arbitrary, we do not obtain the extension of the principal Lie algebra. Then, in order to look for extensions of the principal algebra, we observe that from Equation (50), the form of function f must have the following structure:
  • f = f 0 u r .
  • f = f 0 u s and s r .
We study these cases separately.
  • f = f 0 u r
    In this case, from Equation (50), we have a 1 = 0 . Moreover, by differentiating Equation (51) with respect to u, we have:
    β + β ( r Γ 1 u Γ 1 ) = 0 .
    We observe that if Γ 1 is arbitrary, we have β = b 0 , while from Equation (51), we have λ = 0 , and Equation (52) is satisfied; however, we do not obtain the extension of the principal Lie algebra. Then, in order to have extensions of the principal algebra, the following conditions must be satisfied:
    u Γ 1 r Γ 1 = γ 0 ,
    β = γ 0 β .
    We distinguish two cases: r 1 and r = 1 .
    (a)
    If r 1 , from Equation (54), we get:
    Γ 1 ( u ) = c 1 u 1 + r 1 + r γ 0 u r + c 2 .
    Consequently, from Equation (51), we obtain:
    λ ( t , v ) = ( c 2 + Γ 2 ) ( r + 1 ) β r Γ 2 ,
    while Equation (52) becomes:
    β r Γ 2 2 J 1 = 0
    with:
    J 1 h ( 1 + r ) ( c 2 + Γ 2 ) Γ 2 ( h ( 1 + 2 r ) u h u ) Γ 2 2 + ( 1 + r ) ( h v γ 0 ) ( c 2 + Γ 2 ) Γ 2 .
    We observe that if β = 0 , then λ = 0 , and we do not obtain the extension of the principal Lie algebra. Consequently, in order to have extensions of the principal algebra, the functions Γ 2 and h must satisfy the equation J 1 = 0 . In this case, we have two possible generators depending on γ 0 .
    i.
    If γ 0 0 , as from Equation (55), we have:
    β ( t ) = b 0 + b 1 e γ 0 t ,
    the additional generator is:
    X 3 = e γ 0 t t γ 0 e γ 0 t r u u ( c 2 + Γ 2 ) ( r + 1 ) γ 0 e γ 0 t r Γ 2 v .
    ii.
    If γ 0 = 0 , as from Equation (55), we have:
    β ( t ) = b 0 + b 1 t ,
    the additional generator is:
    X 3 = t t u r u ( c 2 + Γ 2 ) ( r + 1 ) r Γ 2 v .
    (b)
    If r = 1 , from Equation (54), we get:
    Γ 1 ( u ) = c 1 ln ( u ) + γ 0 u + c 2 .
    Consequently, from Equation (51), we obtain:
    λ ( t , v ) = c 1 β Γ 2 ,
    while Equation (52) becomes:
    β Γ 2 2 J 2 = 0
    with:
    J 2 h c 1 Γ 2 ( h + u h u ) Γ 2 2 + c 1 ( h v γ 0 ) Γ 2 .
    We observe that if β = 0 , then λ = 0 , and we do not obtain extension of the principal Lie algebra. Consequently, in order to have extensions of the principal algebra, the functions Γ 2 and h must satisfy the equation J 2 = 0 . In this case, we have two possible generators depending on γ 0 .
    i.
    If γ 0 0 , as from Equation (55), we have:
    β ( t ) = b 0 + b 1 e γ 0 t ,
    the additional generator is:
    X 3 = e γ 0 t t + γ 0 e γ 0 t u u c 1 γ 0 e γ 0 t Γ 2 v .
    ii.
    If γ 0 = 0 , as from Equation (55), we have:
    β ( t ) = b 0 + b 1 t ,
    the additional generator is:
    X 3 = t t u r u c 1 Γ 2 v .
  • f = f 0 u s and s r
    In this case, from Equation (50), we have:
    β ( t ) = a 1 ( 2 r s ) r s t + b 0 ,
    and Equation (51) becomes:
    a 1 u Γ 1 a 1 ( 1 + 2 r s ) ( Γ 1 + Γ 2 ) ( r s ) λ Γ 2 = 0 .
    Moreover, by differentiating with respect to u, we get:
    a 1 ( u Γ 1 + ( s 2 r ) Γ 1 ) = 0 .
    We observe that if Γ 1 is arbitrary, then we have a 1 = 0 , while from Equation (51) λ = 0 and Equation (52) is satisfied, but we do not obtain the extension of the principal Lie algebra. Then, in order to have extensions of the principal algebra, the following condition must be satisfied:
    u Γ 1 + ( s 2 r ) Γ 1 = 0 .
    We distinguish the following two cases.
    (a)
    If s 2 r + 1 , from Equation (75), we get:
    Γ 1 ( u ) = c 1 + c 2 u 1 + 2 r s .
    From Equation (51):
    λ ( t , v ) = a 1 ( s 2 r 1 ) ( c 1 + Γ 2 ) ( r s ) Γ 2 ,
    while Equation (52) becomes:
    a 1 ( r s ) Γ 2 2 J 3 = 0
    with:
    J 3 ( 1 + 2 r s ) ( c 1 + Γ 2 ) ( h Γ 2 + h v Γ 2 ) + ( u h u h ( 1 + 4 r 2 s ) ) Γ 2 2 .
    We observe that if the functions Γ 2 and h do not satisfy the equation J 3 = 0 , we do not obtain the extension of the principal Lie algebra. Then, in order to have extensions of the principal algebra, the functions Γ 2 and h must satisfy the equation J 3 = 0 . In this case, we obtain the following additional generator:
    X 3 = x x + 2 r s r s t t + 1 s r u u + ( s 2 r 1 ) ( c 1 + Γ 2 ) ( r s ) Γ 2 v .
    (b)
    If s = 2 r + 1 , from Equation (75), we get:
    Γ 1 ( u ) = c 1 ln ( u ) + c 2 ,
    and from Equation (72):
    β ( t ) = a 1 r + 1 t + b 0 .
    Consequently, from Equation (51), we obtain:
    λ ( t , v ) = c 1 a 1 ( r + 1 ) Γ 2 ,
    while Equation (52) becomes:
    a 1 ( r + 1 ) Γ 2 2 J 4 = 0
    with:
    J 4 h c 1 Γ 2 ( h + u h u ) Γ 2 2 + c 1 h v Γ 2 .
    We observe that if the functions Γ 2 and h do not satisfy the equation J 4 = 0 , we do not obtain the extension of the principal Lie algebra. Then, in order to have extensions of the principal algebra, the functions Γ 2 and h must satisfy the equation J 4 = 0 . In this case, we obtain the following additional generator:
    X 3 = x x + 1 r + 1 t t + 1 r + 1 u u c 1 ( r + 1 ) Γ 2 v .
Summarizing, we obtained six subclasses of class (28), which admit a three-dimensional Lie algebra.
  • f = f 0 u r with r 1 , Γ 1 ( u ) = c 1 u 1 + r 1 + r γ 0 u r + c 2 with γ 0 0 , the functions h and Γ 2 linked from the following relation:
    h ( 1 + r ) ( c 2 + Γ 2 ) Γ 2 ( h ( 1 + 2 r ) u h u ) Γ 2 2 + ( 1 + r ) ( h v γ 0 ) ( c 2 + Γ 2 ) Γ 2 = 0 .
  • f = f 0 u r with r 1 , Γ 1 ( u ) = c 1 u 1 + r 1 + r + c 2 and the functions h and Γ 2 linked from the following relation:
    h ( 1 + r ) ( c 2 + Γ 2 ) Γ 2 ( h ( 1 + 2 r ) u h u ) Γ 2 2 + ( 1 + r ) ( h v ) ( c 2 + Γ 2 ) Γ 2 = 0 .
  • f = f 0 u , Γ 1 ( u ) = c 1 ln ( u ) + γ 0 u + c 2 with γ 0 0 and the functions h and Γ 2 linked from the following relation:
    h c 1 Γ 2 ( h + u h u ) Γ 2 2 + c 1 ( h v γ 0 ) Γ 2 = 0 .
  • f = f 0 u , Γ 1 ( u ) = c 1 ln ( u ) + c 2 and the functions h and Γ 2 linked from the following relation:
    h c 1 Γ 2 ( h + u h u ) Γ 2 2 + c 1 ( h v ) Γ 2 = 0 .
  • f = f 0 u s with s r , 2 r + 1 , Γ 1 ( u ) = c 1 + c 2 u 1 + 2 r s and the functions h and Γ 2 linked from the following relation:
    ( 1 + 2 r s ) ( c 1 + Γ 2 ) ( h Γ 2 + h v Γ 2 ) + ( u h u h ( 1 + 4 r 2 s ) ) Γ 2 2 = 0 .
  • f = f 0 u 2 r + 1 , Γ 1 ( u ) = c 1 ln ( u ) + c 2 and the functions h and Γ 2 linked from the following relation:
    h c 1 Γ 2 ( h + u h u ) Γ 2 2 + c 1 h v Γ 2 = 0 .

3.2. A Special Case

In agreement with some news about the biological compatibility of the form of g derived from some previous papers (see, e.g., [29,30] and references insides), in this subsection, we show an example of the application of the previous results.
By selecting the case 1a from the obtained cases and assuming r = 1 and c 2 = 0 in (56), we consider f , Γ 1 , of the following form:
f = f 0 u , Γ 1 ( u ) = γ 1 u 2 γ 0 u ,
with f 0 , γ 0 , γ 1 , arbitrary constants. Moreover, we assume:
Γ 2 ( v ) = γ 2 v + γ 3 ,
with γ 2 , γ 3 , arbitrary constants.
In this case, in order to have an extension on the principal algebra, the function h ( u , v ) must satisfy the equation J 1 = 0 , that is:
( u h u 3 h ) γ 2 2 + 2 ( h v γ 0 ) ( γ 2 v + γ 3 ) γ 3 = 0 .
Solutions of this equation are functions h ( u , v ) of the form:
h ( u , v ) = u 3 H ( σ ) 2 γ 0 γ 2 ( γ 2 v + γ 3 ) ,
where H is an arbitrary function of σ = γ 2 v + γ 3 γ 2 u 2 . By assuming H ( σ ) = σ in agreement with [13,29,30], we get:
h ( u , v ) = γ 2 v + γ 3 γ 2 ( u 2 γ 0 ) .
The system (28) becomes:
u t = f 0 u x 2 + f 0 u u x x + u u x + γ 1 u 2 γ 0 u + γ 2 v + γ 3 , v t = γ 2 v + γ 3 γ 2 ( u 2 γ 0 ) .
While the third generator is obtained by specializing generator (61) and has the form:
X 3 = e γ 0 t t γ 0 e γ 0 t u u 2 γ 0 γ 2 ( γ 2 v + γ 3 ) e γ 0 t v .
By considering the generator k X 2 + X 3 , we get:
u ( t , x ) = U ( z ) e γ 0 t , v ( t , x ) = V ( z ) e 2 γ 0 t γ 3 γ 2 ,
with z = γ 0 x k + e γ 0 t , while the functions U ( z ) , V ( z ) are solutions of the reduced system:
U U f 0 γ 0 2 + U 2 f 0 γ 0 2 + k γ 0 ( U + k ) U + k 2 V + γ 1 k 2 U 2 = 0 ,
γ 0 V + U V = 0 .
This reduced system, as well as other cases of biological specializations, will be studied in later research.
Of course, the systems studied here cannot be considered, strictu sensu, as mathematical models. In fact, their constitutive parameters need to be characterized carefully from the biological point view. However, having in mind some previous models concerned with Aedes aegypti [13,14,29,30,31], we try to stress some structural features of the system (98). To this aim, we rewrite system (98) as:
u t = ( f 0 u u x ) x + u u x + γ 1 u u γ 0 γ 1 + γ 2 v + γ 3 γ 2 v t = ( u 2 γ 0 ) v + γ 3 γ 2
It is easy to ascertain a weak interaction of the equation for the aquatic population on the equation for the winged population. Moreover, in this last one appears a growth for population u having a logistic structure. By identifying γ 1 as a positive rate of maturation of the aquatic forms in winged female mosquitoes and γ 0 as the positive winged population mortality, it is possible to find a threshold value u t r s = γ 0 γ 1 . Finally, the aquatic population evolution equation shows a growth rate u 2 γ 0 ruled by the density of mosquitoes and their mortality.
Remark 1. It is a simple matter to ascertain that the system (103) admits as the special solution:
u = γ 0 γ 1 , v = γ 3 γ 2 .
Moreover, it is possible to get other solutions by assuming v = γ 3 γ 2 , while u is obtained as a solution of equation:
u t = ( f 0 u u x ) x + u u x + γ 1 u u γ 0 γ 1 .
For the interested reader, it could be worthwhile noticing that Equation (105) is a particular case of equation considered in [8,9]. Moreover, we can get the results obtained in [8,9] (see Table 1, Case 8 of both papers) by projection. Indeed, the Lie symmetry generator (99) projected in the space t , x , u , v = γ 3 γ 2 becomes the Lie symmetry T 1 of [8,9].

4. Conclusions

In this paper, we have considered a class of advection reaction diffusion systems of interest in biomathematics. After having recalled a weak equivalence generator, obtained in a previous work [12], we find some particular cases of the nonlinear system (28) admitting three-dimensional Lie algebras by using a specialization of a projection theorem [24,28]. In this subclass, the constitutive equation characterizing g is assigned as:
g = u r u x + Γ 1 ( u ) + Γ 2 ( v )
that generalizes [12]:
g = ρ u r u x s + Γ 1 u a + Γ 2 v b ,
where the constants ρ ( 0 ) , Γ 1 , Γ 2 , r , s , a and b are constitutive parameters. We derive the principal Lie algebra and the functions admitting at least an extension by one. These results are summarized at the end of Section 3.1. A special case is considered in Section 3.2.

Acknowledgments

The authors wish to thank the reviewers for their interesting observations and the editor for his very useful comments. Mariano Torrisi was supported from Gruppo Nazionale per la Fisica Matematica of Istituto Nazionale di Alta Matematica.

Author Contributions

The authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cardile, V.; Torrisi, M.; Tracinà, R. On a reaction-diffusion system arising in the development of Bacterial Colonies. In Proceedings of the 10th International Conference in Modern Group Analysis, Larnaca, Cyprus, 24–31 October 2004; Volume 32, p. 38.
  2. 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]
  3. Torrisi, M.; Tracinà, R. On a class of reaction diffusion systems: Equivalence transformations and symmetries. In Asymptotic Methods in Nonlinear Wave Phenomena; Ruggeri, T., Sammartino, M., Eds.; World Science Publishing Co. Pte. Ltd.: Singapore, 2007; pp. 207–216. [Google Scholar]
  4. Torrisi, M.; Tracinà, R. Exact solutions of a reaction-diffusion system for Proteus Mirabilis bacterial colonies. Nonlinear Anal. Real World Appl. 2011, 12, 1865–1874. [Google Scholar] [CrossRef]
  5. Cherniha, R.; Serov, M. Nonlinear systems of the burgers-type equations: Lie and Q-conditional symmetries, ansätze and solutions. J. Math. Anal. Appl. 2003, 282, 305–328. [Google Scholar] [CrossRef]
  6. Cherniha, R.; Wilhelmsson, H. Symmetry and exact solution of heat-mass transfer equations in thermonuclear plasma. Ukr. Math. J. 1996, 48, 1434–1449. [Google Scholar] [CrossRef]
  7. Cherniha, R.; Serov, M. Lie and non-Lie symmetries of nonlinear diffusion equations with convection term. Symmetry Nonlinear Math. Phys. 1997, 2, 444–449. [Google Scholar]
  8. Cherniha, R.; Serov, M. Symmetries, ansätze and exact solutions of nonlinear second-order evolution equations with convection terms. Eur. J. Appl. Math. 1998, 9, 527–542. [Google Scholar] [CrossRef]
  9. Cherniha, R.; Serov, M. Symmetries, ansätze and exact solutions of nonlinear second-order evolution equations with convection terms, II. Eur. J. Appl. Math. 2006, 17, 597–605. [Google Scholar] [CrossRef]
  10. Ibragimov, N.H.; Torrisi, M.; Valenti, A. Preliminary group classification of equation vtt = f(x, vx)vxx + g(x, vx). J. Math. Phys. 1991, 32, 2988–2995. [Google Scholar] [CrossRef]
  11. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  12. Freire, I.L.; Torrisi, M. Weak equivalence transformations for a class of models in biomathematics. Abstr. Appl. Anal. 2014. [Google Scholar] [CrossRef]
  13. Freire, I.L.; Torrisi, M. Symmetry methods in mathematical modeling Aedes aegypti dispersal dynamics. Nonlinear Anal. Real World Appl. 2013, 14, 1300–1307. [Google Scholar] [CrossRef]
  14. Freire, I.L.; Torrisi, M. Similarity solutions for systems arising from an Aedes aegypti model. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 872–879. [Google Scholar] [CrossRef]
  15. Romano, V.; Torrisi, M. Application of weak equivalence transformations to a group analysis of a drift-diffusion model. J. Phys. A Math. Gen. 1999, 32, 7953–7963. [Google Scholar] [CrossRef]
  16. Torrisi, M.; Tracinà, R. Equivalence transformations and symmetries for a heat conduction model. Int. J. Non-Linear Mech. 1998, 33, 473–487. [Google Scholar] [CrossRef]
  17. Gazeau, J.P.; Winternitz, P. Symmetries of variable-coefficient Korteweg-de Vries equations. J. Math. Phys. 1992, 33, 4087–4102. [Google Scholar] [CrossRef]
  18. Winternitz, P.; Gazeau, J.P. Allowed transformations and symmetry classes of variable coefficient Korteweg-de Vries equations. Phys. Lett. A 1992, 167, 246–250. [Google Scholar] [CrossRef]
  19. Akhatov, I.S.H.; Gazizov, R.K.; Ibragimov, N.H. Nonlocal symmetries. Heuristic approach. J. Sov. Math. 1991, 55, 1401–1450. [Google Scholar] [CrossRef]
  20. Lisle, I.G. Equivalence Transformation for Classes of Differential Equations. Ph.D. Thesis, University of British Columbia, Vancouver, BC, Canada, 1992. [Google Scholar]
  21. Khalique, C.M.; Mahomed, F.M.; Ntsime, B.P. Group classification of the generalized Emden-Fowler-type equation. Nonlinear Anal. Real World Appl. 2009, 10, 3387–3395. [Google Scholar] [CrossRef]
  22. Ibragimov, N.H. CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  23. Molati, M.; Khalique, C.M. Lie group classification of a generalized Lane–Emden Type system in two dimensions. J. Appl. Math. 2012. [Google Scholar] [CrossRef]
  24. Torrisi, M.; Tracinà, R. Equivalence transformations for systems of first order quasilinear partial differential equations. In Modern Group Analysis VI: Developments in Theory, Computation and Application; New Age International(P) Ltd.: New Delhi, India, 1996; pp. 115–135. [Google Scholar]
  25. Torrisi, M.; Tracinà, R.; Valenti, A. Group analysis approach for a non linear differential system arising in diffusion phenomena. J. Math. Phys. 1996, 37, 4758–4767. [Google Scholar] [CrossRef]
  26. Gambino, G.; Greco, A.M.; Lombardo, M.C. A group analysis via weak equivalence transformations for a model of tumour encapsulation. J. Phys. A 2004, 37, 3835–3846. [Google Scholar] [CrossRef]
  27. Ibragimov, N.H.; Säfström, N. The equivalence group and invariant solutions of a tumour growth model. Commun. Nonlinear Sci. Num. Simul. 2004, 9, 61–69. [Google Scholar] [CrossRef]
  28. Ibragimov, N.H.; Torrisi, M. A simple method for group analysis and its application to a model of detonation. J. Math. Phys. 1992, 33, 3931–3937. [Google Scholar] [CrossRef]
  29. Maidana, N.A.; Yang, H.M. Describing the geographic spread of dengue disease by traveling waves. Math. Biosci. 2008, 215, 64–77. [Google Scholar] [CrossRef] [PubMed]
  30. Takahashi, L.T.; Maidana, N.A.; Ferreira, W.C., Jr.; Pulino, P.; Yang, H.M. Mathematical models for the Aedes aegypti dispersal dynamics: Traveling waves by wing and wind. Bull. Math. Biol. 2005, 67, 509–528. [Google Scholar] [CrossRef] [PubMed]
  31. Bacani, F.; Freire, I.L.; Maidana, N.A.; Torrisi, M. Modelagem para a dinâmica populacional do Aedes aegypti via simetrias de Lie. Proc. Ser. Braz. Soc. Appl. Comput. Math. 2015, 3. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Torrisi, M.; Tracinà, R. An Application of Equivalence Transformations to Reaction Diffusion Equations. Symmetry 2015, 7, 1929-1944. https://doi.org/10.3390/sym7041929

AMA Style

Torrisi M, Tracinà R. An Application of Equivalence Transformations to Reaction Diffusion Equations. Symmetry. 2015; 7(4):1929-1944. https://doi.org/10.3390/sym7041929

Chicago/Turabian Style

Torrisi, Mariano, and Rita Tracinà. 2015. "An Application of Equivalence Transformations to Reaction Diffusion Equations" Symmetry 7, no. 4: 1929-1944. https://doi.org/10.3390/sym7041929

APA Style

Torrisi, M., & Tracinà, R. (2015). An Application of Equivalence Transformations to Reaction Diffusion Equations. Symmetry, 7(4), 1929-1944. https://doi.org/10.3390/sym7041929

Article Metrics

Back to TopTop