Next Article in Journal
Anonymous Multi-Receiver Identity-Based Authenticated Encryption with CCA Security
Next Article in Special Issue
An Application of Equivalence Transformations to Reaction Diffusion Equations
Previous Article in Journal
On the Continuity of the Hutchinson Operator
Previous Article in Special Issue
Lie Group Method for Solving the Generalized Burgers’, Burgers’–KdV and KdV Equations with Time-Dependent Variable Coefficients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems

Department of Mathematics, Poltava National Technical Yuriy Kondratyuk University, 24, Pershotravnevyi Prospekt, 36601 Poltava, Ukraine
Symmetry 2015, 7(4), 1841-1855; https://doi.org/10.3390/sym7041841
Submission received: 28 June 2015 / Revised: 11 September 2015 / Accepted: 1 October 2015 / Published: 16 October 2015
(This article belongs to the Special Issue Lie Theory and Its Applications)

Abstract

:
A wide range of reaction–diffusion systems with constant diffusivities that are invariant under Q-conditional operators is found. Using the symmetries obtained, the reductions of the corresponding systems to the systems of ODEs are conducted in order to find exact solutions. In particular, the solutions of some reaction–diffusion systems of the Lotka–Volterra type in an explicit form and satisfying Dirichlet boundary conditions are obtained. An biological interpretation is presented in order to show that two different types of interaction between biological species can be described.

1. Introduction

In 1952, Alan Turing published his prominent paper [1]. In this paper he proposed the Turing hypothesis of pattern formation. He used reaction–diffusion equations of the form
λ 1 u t = ( D 1 ( u ) u x ) x + F ( u , v ) , λ 2 v t = ( D 2 ( v ) v x ) x + G ( u , v )
which are central to the field of pattern formation.
In system (1), F and G are arbitrary smooth functions, u = u ( t , x ) and v = v ( t , x ) are unknown functions of the variables t and x, while the subscripts t and x denote differentiation with respect to this variable. Nonlinear system (1) generalizes many well-known nonlinear second-order models used to describe various processes in physics [2], biology [3,4,5] and ecology [6].
Here we concentrate ourselves on the most important subclass of RD systems with the form of (1), namely that with constant coefficients of diffusivity
λ 1 u t = u x x + F ( u , v ) , λ 2 v t = v x x + G ( u , v )
System (2) has been intensely studied using different mathematical methods (see, e.g., [3,4,7] and papers cited therein). All possible Lie symmetries of system (2) have been found, in [8,9,10,11]. In particular, Q-conditional symmetries of (2) were found in [12]. Reference [13] also contains some results related with system (2).
System (1) is a natural generalization of the well-known RD equation
u t = [ D ( u ) u x ] x + F ( u )
There are many papers devoted to the construction of Q-conditional symmetries for this equation [14,15,16,17,18,19,20,21], starting from the pioneering work in [22]. There is also a non-trivial generalization of these results for the case of the reaction–diffusion–convection equation ([21] and papers cited therein).
In contrast to (3), there are not many results for searching Q-conditional symmetries of system (2). Construction of the Q-conditional symmetries (non-classical symmetries) of such systems is a very difficult task. Only a few papers have been devoted to the search of such symmetries. In [23] the Q-conditional symmetries of the system
u t = ( u k u x ) x + F ( u , v ) , v t = ( v l v x ) x + G ( u , v ) , l 2 + k 2 0 ,
have been obtained; in [24] the Q-conditional symmetries of the Lotka–Volterra system
λ 1 u t = u x x + u ( a 1 + b 1 u + c 1 v ) , λ 2 v t = v x x + v ( a 2 + b 2 u + c 2 v )
were obtained.
The paper is organized as follows. In Section 2 three theorems are presented which contain the main result for Q-conditional symmetries of system (2). In Section 3, ansätze for all systems and solutions for one of the systems are derived. In Section 4, the solutions for a generalization of the Lotka–Volterra system are obtained and analyzed. Some graphs of the exact solutions are also presented. Finally, we present some conclusions.

2. Main Result

Let us consider the reaction–diffusion system with constant diffusivities: (2). We want to find Q-conditional operators of the form
Q = t + ξ ( t , x , u , v ) x + η 1 ( t , x , u , v ) u + η 2 ( t , x , u , v ) v
under which system (2) is invariant.
The most general form of the Q-conditional operators is
Q = ξ 0 ( t , x , u , v ) t + ξ 1 ( t , x , u , v ) x + η 1 ( t , x , u , v ) u + η 2 ( t , x , u , v ) v .
In the case ξ 0 ( t , x , u , v ) 0 , this operator can be reduced to that with ξ 0 ( t , x , u , v ) = 1 [25]. So we investigate operator (5).
We write down system (2) in the following form:
u x x = λ 1 u t + C 1 ( u , v ) , λ 1 0 , v x x = λ 2 v t + C 2 ( u , v ) , λ 2 0
where C 1 ( u , v ) = F ( u , v ) , C 2 ( u , v ) = G ( u , v ) .
The determining equations for finding coefficients of operator (5) and functions C 1 ( u , v ) , C 2 ( u , v ) from system (6) have the form
1 ) ξ u u = ξ v v = ξ u v = 0 , 2 ) η v v 1 = 0 , 3 ) η u u 2 = 0 , 4 ) 2 λ 1 ξ ξ u + η u u 1 2 ξ x u = 0 , 5 ) 2 λ 2 ξ ξ v + η v v 2 2 ξ x v = 0 , 6 ) ( λ 1 + λ 2 ) ξ ξ v + 2 η u v 1 2 ξ x v = 0 , 7 ) ( λ 1 + λ 2 ) ξ ξ u + 2 η u v 2 2 ξ x u = 0 , 8 ) ( λ 1 λ 2 ) ξ η v 1 + 2 η x v 1 2 ξ v C 1 2 λ 1 ξ v η 1 = 0 , 9 ) ( λ 2 λ 1 ) ξ η u 2 + 2 η x u 2 2 ξ u C 2 2 λ 2 ξ u η 2 = 0 , 10 ) λ 1 ( 2 ξ u η 1 ξ t ξ v η 2 2 ξ ξ x ) + λ 2 ξ v η 2 + 3 ξ u C 1 + ξ v C 2 2 η x u 1 + ξ x x = 0 , 11 ) λ 2 ( 2 ξ v η 2 ξ t ξ u η 1 2 ξ ξ x ) + λ 1 ξ u η 1 + 3 ξ v C 2 + ξ u C 1 2 η x v 2 + ξ x x = 0 , 12 ) λ 1 ( η t 1 + η 2 η v 1 + 2 ξ x η 1 ) λ 2 η 2 η v 1 + η 1 C u 1 + η 2 C v 1 η u 1 C 1 + 2 ξ x C 1 η v 1 C 2 η x x 1 = 0 , 13 ) λ 2 ( η t 2 + η 1 η u 2 + 2 ξ x η 2 ) λ 1 η 1 η u 2 + η 1 C u 2 + η 2 C v 2 η u 2 C 1 + 2 ξ x C 2 η v 2 C 2 η x x 2 = 0 .
System (7) is an over-determined system of partial differential equations and there are no any general method for solving of such systems [26,27]. Thus, we were not able to find the general solution of system (7), hence we have solved it with conditions
ξ = ξ ( u , v ) , η i = η i ( u , v ) , i = 1 , 2
Solving Equations 1 ) 3 ) of system (7), we obtain
ξ = a u + b v + c , η 1 = p 1 ( u ) v + q 1 ( u ) , η 2 = p 2 ( v ) u + q 2 ( v )
where a , b , c are arbitrary constants, p 1 , p 2 , q 1 , q 2 are arbitrary smooth functions. Substituting (9) into 6 ) , 7 ) from (7) and splitting the obtained equations with respect to the powers of u and v, we arrive at the system
a 2 ( λ 1 + λ 2 ) = 0 , b 2 ( λ 1 + λ 2 ) = 0 , ( λ 1 + λ 2 ) a ( b v + c ) + 2 p v 2 = 0 , ( λ 1 + λ 2 ) b ( a u + c ) + 2 p u 1 = 0
Obviously, that solutions of first pair of equations of (10) will be λ 2 = λ 1 , or a = b = 0 .
Let us consider the case λ 2 = λ 1 (the case a = b = 0 will be considered later). In this case we obtain p 1 = c o n s t = α 1 , p 2 = c o n s t = α 2 . Substituting (9) p 1 = α 1 , p 2 = α 2 into Equations 4 ) and 5 ) of system (7) and splitting with respect to the powers of u and v, we arrive at
2 λ 1 a b = 0 , q u u 1 + 2 λ 1 a ( a u + c ) = 0 , q v v 2 2 λ 1 b ( b v + c ) = 0
Since λ 1 0 , we conclude that a b = 0 . Consider the case a = 0 , b 0 (the case b = 0 , a 0 is symmetrical). From Equation 8 ) , we obtain α 2 = 0 . Substituting (9) with the specified coefficients, namely
ξ = b v + c , η 1 = α 1 v + q 1 , η 2 = q 2 ,
into Equations 10 ) , 11 ) of system (7), we arrive at
b ( C 2 2 λ 1 q 2 ) = 0 , b ( 3 C 2 2 λ 1 q 2 ) = 0
Substituting q 2 = 0 , obtained from (12), into the third equation of system (11), we obtain λ 1 b ( b v + c ) = 0 , that is λ 1 b = 0 , but that contradicts the above restrictions.
Thus, in the case λ 2 = λ 1 we do not obtain any Q-conditional operator of the form (5).
Consider the case a = b = 0 . In this case, from Equations 4 ) , 5 ) , 6 ) and 7 ) of system (7), we obtain
p i = c o n s t = α i , i = 1 , 2 , q 1 = β 1 u + γ 1 , q 2 = β 2 v + γ 2 ,
where β i , γ i , i = 1 , 2 are the arbitrary constants. Thus, expressions (9) take the form
ξ = c , η 1 = α 1 v + β 1 u + γ 1 , η 2 = α 2 u + β 2 v + γ 2
Substituting (13) into Equations 8 ) and 9 ) of system (7), we arrive at
c α 1 ( λ 2 λ 1 ) = 0 , c α 2 ( λ 2 λ 1 ) = 0
Solving the system of algebraic Equations (14), we obtain three solutions λ 2 = λ 1 , α 1 = α 2 = 0 and c = 0 , therefore we obtain three cases. Let us consider all these cases.
Theorem 1. In the cases λ 2 = λ 1 or η v 1 = η u 2 = 0 with conditions (8), the system of determining equations for finding of the Q-conditional operators of the form (5) for system (6) coincide with the system of determining equations for finding Lie operators.
Proof. Substituting (13), with λ 2 = λ 1 , into system (7) we find that Equations 1 ) 11 ) are transformed into identities, and Equations 12 ) and 13 ) take the form
η 1 C u 1 + η 2 C v 1 η u 1 C 1 η v 1 C 2 = 0 , η 1 C u 2 + η 2 C v 2 η u 2 C 1 η v 2 C 2 = 0
In [11] the determining equations for finding of Lie symmetries with condition λ 2 = λ 1 are written down in explicit form. Substituting conditions (8) into these equations, we see that the result is completely identical to Equations (15).
Substituting (13), with α 1 = α 2 = 0 into system (7), we see that Equations 1 ) 11 ) also transform into identities, and equations 12 ) and 13 ) take the form
η 1 C u 1 + η 2 C v 1 η u 1 C 1 = 0 , η 1 C u 2 + η 2 C v 2 η v 2 C 2 = 0 .
Comparing equations (16) with equations for finding of Lie symmetries of system (6) with conditions (8) from [9], we see that they are completely identical. ☐
Thus, in the following we assume that λ 1 λ 2 , α 1 2 + α 2 2 0 .
Let us consider the case c = 0 , which is on the one hand the most interesting and on the other the most difficult. In this case, (13) takes the form
ξ = 0 , η 1 = α 1 v + β 1 u + γ 1 , η 2 = α 2 u + β 2 v + γ 2
Equations 1 ) 11 ) satisfy expressions (17) and Equations 12 ) , 13 ) take the form
( α 1 v + β 1 u + γ 1 ) C u 1 + ( α 2 u + β 2 v + γ 2 ) C v 1 β 1 C 1 α 1 C 2 = α 1 ( λ 2 λ 1 ) ( α 2 u + β 2 v + γ 2 ) , ( α 1 v + β 1 u + γ 1 ) C u 2 + ( α 2 u + β 2 v + γ 2 ) C v 2 α 2 C 1 β 2 C 2 = α 2 ( λ 1 λ 2 ) ( α 1 v + β 1 u + γ 1 )
Thus, we can formulate the following theorem.
Theorem 2. The nonlinear reaction–diffusion system (6) is Q-conditionally invariant under operator (5) with coefficients (17) if and only if the nonlinearities C 1 , C 2 are the solutions of linear system (18).
To find the general solution of system (18), one need to analyze two cases α 2 = 0 and α 2 0 . In the case α 2 = 0 , system (18) takes the form
( α 1 v + β 1 u + γ 1 ) C u 1 + ( β 2 v + γ 2 ) C v 1 = β 1 C 1 + α 1 C 2 + α 1 ( λ 2 λ 1 ) ( β 2 v + γ 2 ) , ( α 1 v + β 1 u + γ 1 ) C u 2 + ( β 2 v + γ 2 ) C v 2 = β 2 C 2
Since α 1 0 , renaming C 1 α 1 C 1 , u α 1 u and γ 1 α 1 γ 1 , and taking into account that with any coefficients β 1 , β 2 we can remove the parameter γ 1 using linear substitutions of u , v , system (19) reduces to the form
( v + β 1 u ) C u 1 + ( β 2 v + γ 2 ) C v 1 = β 1 C 1 + C 2 + ( λ 2 λ 1 ) ( β 2 v + γ 2 ) , ( v + β 1 u ) C u 2 + ( β 2 v + γ 2 ) C v 2 = β 2 C 2
One notes a particular solution of system (20), of the form
C p a r t 1 = 1 2 ( λ 2 λ 1 ) ( v + β 1 u ) , C p a r t 2 = 1 2 ( λ 1 λ 2 ) ( β 2 v + γ 2 )
Now to construct the general solution of (20), we need to solve the corresponding homogeneous system, that is
( v + β 1 u ) C u 1 + ( β 2 v + γ 2 ) C v 1 = β 1 C 1 + C 2 , ( v + β 1 u ) C u 2 + ( β 2 v + γ 2 ) C v 2 = β 2 C 2
As a result, the following statement was proved.
Theorem 3. Reaction–diffusion system (6) is Q-conditionally invariant under operator (5) with conditions (8), and η u 2 = 0 , if and only if the system and corresponding operator have one of the seven following forms (moreover λ 2 λ 1 ):
u x x = λ 1 u t + g ( ω ) v ln ( v ) + h ( ω ) v + 1 2 ( λ 2 λ 1 ) ( v + β 1 u ) , v x x = λ 2 v t + β 1 g ( ω ) v + 1 2 ( λ 1 λ 2 ) β 1 v , Q = t + ( v + β 1 u ) u + β 1 v v , ω = v 1 exp ( β 1 u v ) , β 1 0
u x x = λ 1 u t + h ( ω ) v β 1 β 2 + g ( ω ) v + 1 2 ( λ 2 λ 1 ) ( v + β 1 u ) , v x x = λ 2 v t + ( β 2 β 1 ) g ( ω ) v + 1 2 ( λ 1 λ 2 ) β 2 v , Q = t + ( v + β 1 u ) u + β 2 v v , ω = v β 1 β 2 ( ( β 1 β 2 ) u + v ) , β 1 β 2 ( β 1 β 2 ) 0
u x x = λ 1 u t + u g ( v ) + h ( v ) , v x x = λ 2 v t + v g ( v ) , Q = t + v u
u x x = λ 1 u t + g ( ω ) v + h ( ω ) + 1 2 ( λ 2 λ 1 ) v , v x x = λ 2 v t + g ( ω ) + 1 2 ( λ 1 λ 2 ) , Q = t + v u + v , ω = 2 u v 2
u x x = λ 1 u t + g ( ω ) v + h ( ω ) + 1 2 ( λ 2 λ 1 ) v , v x x = λ 2 v t + β 2 g ( ω ) v + τ 2 + 1 2 β 2 ( λ 1 λ 2 ) ( v + τ 2 ) , Q = t + v u + β 2 ( v + τ 2 ) v , ω = β 2 u v + τ 2 ln v + τ 2 , β 2 0
u x x = λ 1 u t + ( v + u ) h ( v ) g ( v ) , v x x = λ 2 v t + g ( v ) , Q = t + β 1 ( v + u ) u , β 1 0
u x x = λ 1 u t + h ( ω ) exp ( β 1 γ 2 v ) g ( ω ) + 1 2 ( λ 2 λ 1 ) ( v + β 1 u ) , v x x = λ 2 v t + β 1 g ( ω ) + 1 2 γ 2 ( λ 1 λ 2 ) , Q = t + ( v + β 1 u ) u + γ 2 v , ω = exp ( β 1 γ 2 v ) β 1 2 u + β 1 v + γ 2 , β 1 γ 2 0
Proof. To prove this theorem, it is necessary and sufficient to construct the general solution of system (22) for all possible ratios between parameters β 1 , β 2 , γ 2 . To do this we need to investigate the following seven cases:
1. β 1 β 2 0 , β 1 = β 2 ;
2. β 1 β 2 0 , β 1 β 2 ;
3. β 2 = 0 , β 1 = 0 , γ 2 = 0 ;
4. β 2 = 0 , β 1 = 0 , γ 2 0 ;
5. β 2 0 , β 1 = 0 .
6. β 2 = 0 , β 1 0 , γ 2 = 0 ;
7. β 2 = 0 , β 1 0 , γ 2 0 .
These cases take into account all possibilities that arise when we solve system (22). Let us consider these cases.
Case 1. Solving the second equation of (22), we get C 2 = β 1 v g ( ω ) and ω = v 1 exp β 1 u v . So the first equation of (22) reduces to an ODE for finding of the function C 1 :
C v 1 C 1 v = g ( ω ) .
Solving it, we get that C 1 = g ( ω ) v ln ( v ) + v h ( ω ) . Taking into account the expressions for C 1 , C 2 , ω , obtained above, C p a r t 1 , C p a r t 2 from Formulas (21) and restrictions (obtained above), finally we arrive at the reaction–diffusion system and the Q-conditional operator listed in (23) of Theorem 3.
Cases 2–7. Considering similarly these cases and using simple renamings, we arrive at systems and operators (24)–(29) of Theorem 3. ☐
In the case α 2 0 we should also assume that α 1 0 , otherwise we obtain the case α 2 = 0 up to renaming. We seek a solution of system (18) of the form
C 1 = r 1 u + r 2 v + r 3 , C 2 = s 1 u + s 2 v + s 3
Substituting (30) into (18), we obtain the system of algebraic equations
α 2 r 2 α 1 s 1 + α 1 α 2 λ 1 λ 2 = 0 , α 1 r 1 + β 2 β 1 r 2 s 2 α 1 + α 1 β 2 λ 1 λ 2 = 0 , γ 1 r 1 + γ 2 r 2 β 1 r 3 α 1 s 3 + α 1 γ 2 λ 1 λ 2 = 0 , α 2 r 1 + β 2 β 1 s 1 s 2 α 2 + α 2 β 1 λ 1 λ 2 = 0 , α 2 r 3 γ 1 s 1 γ 2 s 2 + β 2 s 3 + α 2 γ 1 λ 1 λ 2 = 0
Solving system (31), we arrive at two possibilities depending on Δ = α 1 α 2 β 1 β 2 :
I) Δ = 0 ,
r 1 = β 1 α 2 s 1 + β 1 λ 2 λ 1 , r 2 = α 1 α 2 s 1 + α 1 λ 2 λ 1 , r 3 = α 2 γ 1 + β 2 γ 2 α 2 2 s 1 β 2 α 2 s 3 + γ 1 λ 2 λ 1 , s 2 = β 2 α 2 s 1 .
II) Δ 0 ,
r 1 = β 1 β 2 α 2 s 1 + s 2 + β 1 λ 2 λ 1 , r 2 = α 1 α 2 s 1 + α 1 λ 2 λ 1 , r 3 = Δ + β 2 2 γ 1 α 1 β 2 γ 2 α 2 Δ s 1 + α 1 γ 2 β 2 γ 1 Δ s 2 + γ 1 λ 2 λ 1 , s 3 = α 1 γ 2 β 2 γ 1 Δ s 1 + α 2 γ 1 β 1 γ 2 Δ s 2 .
In Case I) s 1 = s 3 = 0 , we obtain the solution of system (18)
C p a r t 1 = λ 2 λ 1 β 1 α 2 α 2 u + β 2 v + γ 1 , C p a r t 2 = 0
In Case II) s 1 = s 2 = 0 we obtain the solution of system (18)
C p a r t 1 = λ 2 λ 1 α 1 v + β 1 u + γ 1 , C p a r t 2 = 0
Furthermore, we must solve the homogeneous system
( α 1 v + β 1 u + γ 1 ) C u 1 + ( α 2 u + β 2 v + γ 2 ) C v 1 = β 1 C 1 + α 1 C 2 , ( α 1 v + β 1 u + γ 1 ) C u 2 + ( α 2 u + β 2 v + γ 2 ) C v 2 = α 2 C 1 + β 2 C 2
Let us consider Case I). Using the condition Δ = 0 for system (34), we get
( β 1 α 2 ( α 2 u + β 2 v ) + γ 1 ) C u 1 + ( α 2 u + β 2 v + γ 2 ) C v 1 = β 1 α 2 ( α 2 C 1 + β 2 C 2 ) , ( β 1 α 2 ( α 2 u + β 2 v ) + γ 1 ) C u 2 + ( α 2 u + β 2 v + γ 2 ) C v 2 = α 2 C 1 + β 2 C 2
Multiplying the second equation of (35) by β 1 α 2 , adding to the first and renaming u u γ 2 α 2 , γ 1 γ 1 + β 1 γ 2 α 2 , we arrive at
β 1 α 2 ( α 2 u + β 2 v ) + γ 1 C 1 β 1 α 2 C 2 u + ( α 2 u + β 2 v ) C 1 β 1 α 2 C 2 v = 0
Using the substitution
C 1 = S ( u , v ) + β 1 α 2 C 2
we obtain the equation
β 1 α 2 ( α 2 u + β 2 v ) + γ 1 S u + ( α 2 u + β 2 v ) S v = 0
Solving Equation (38), we arrive at three subcases:
1) β 2 = β 1 , γ 1 = 0 , S = S ( ω ) , ω = α 2 u β 1 v ;
2) β 2 = β 1 , γ 1 0 , S = S ( ω ) , ω = ( α 2 u β 1 v ) 2 2 α 2 γ 1 v ;
3) β 2 β 1 , S = S ( ω ) , ω = α 2 u β 1 v α 2 γ 1 β 1 + β 2 ln α 2 u + β 2 v + α 2 γ 1 β 1 + β 2 .
Substituting (37) together with the function S from subcase 1) into the second equation of (35), we obtain
α 2 u β 1 v α 2 C v 2 + β 1 C u 2 = α 2 2 f ( ω )
Solving (39), using (37), (32) and renaming u β 1 u , v α 2 v we obtain the system
u x x = λ 1 u t + f ( ω ) u + β 1 ( g ( ω ) λ 1 ) u v , v x x = λ 2 v t + f ( ω ) v + β 1 ( g ( ω ) λ 2 ) u v , ω = u v ,
Q-conditionally invariant under the operator
Q = t + β 1 ( u v ) ( u + v ) .
Similarly, for subcase 2), we arrive at the system
u x x = λ 1 u t + ( f ( ω ) λ 1 ) β 1 u v + γ 1 + g ( ω ) , v x x = λ 2 v t + β 1 ( f ( ω ) λ 2 ) u v + g ( ω ) , ω = ( u v ) 2 2 γ 1 β 1 v , γ 1 0
and the operator
Q = t + ( β 1 ( u v ) + γ 1 ) u + β 1 ( u v ) v .
In the subcase 3), we obtain the system
u x x = λ 1 u t + k ( f ( ω ) λ 1 ) + β 1 ( g ( ω ) λ 1 ) u + v , v x x = λ 2 v t k ( f ( ω ) λ 2 ) + β 2 ( g ( ω ) λ 2 ) u + v , ω = β 2 u β 1 v k ln | u + v | , k = γ 1 β 1 + β 2 ,
and the operator
Q = t + ( β 1 ( u + v ) + k ) u + ( β 2 ( u + v ) k ) v .
Examination of Case II) is highly nontrivial and will be reported in another paper.

3. Ansätze and Exact Solutions of the Reaction–Diffusion System

Using standard procedures, we obtain ansätze for all operators of Theorem 3. Substituting these anzätze in the corresponding reaction–diffusion systems, we obtain the reduction systems of equations. All anzätze and reduction systems are presented in Table 1.
Table 1. Ansätze and reduction systems of Theorem 3.
Table 1. Ansätze and reduction systems of Theorem 3.
No.AnsätzeSystems of ODEs
(23) u = t + ψ e β 1 t + φ φ + ( φ ) 2 β 1 g e β 1 ψ φ + λ 1 + λ 2 2 = 0
v = e β 1 t + φ ψ + 2 φ ψ + g e β 1 ψ φ β 1 ψ φ h e β 1 ψ φ λ 1 + λ 2 2 = 0
(24) u = ψ e β 2 t + φ β 1 β 2 e β 2 t + φ β 1 β 2 φ + ( φ ) 2 + g β 1 β 2 ψ β 1 β 2 β 2 λ 1 + λ 2 2 = 0
v = e β 2 t + φ ψ + 2 β 1 β 2 φ ψ + β 1 β 2 β 1 β 2 1 ψ ( φ ) 2 β 2 g β 1 β 2 ψ h β 1 β 2 ψ = 0
(25) u = φ t + ψ φ φ g ( φ ) = 0
v = φ ψ g ( φ ) ψ h ( φ ) λ 1 φ = 0
(26) u = 1 2 t 2 + φ t + ψ φ g 2 ψ φ 2 λ 1 + λ 2 2 = 0
v = t + φ ψ g 2 ψ φ 2 φ h 2 ψ φ 2 λ 1 + λ 2 2 φ = 0
(27) u = 1 β 2 e β 2 t + φ τ 2 t + ψ φ + ( φ ) 2 β 2 g τ 2 ( φ + 1 ) + β 2 ψ λ 1 + λ 2 2 β 2 = 0
v = e β 2 t + φ τ 2 ψ h τ 2 ( φ + 1 ) + β 2 ψ + τ 2 g τ 2 ( φ + 1 ) + β 2 ψ + λ 1 + λ 2 2 τ 2 = 0
(28) u = ψ e β 1 t φ φ g ( φ ) = 0
v = φ ψ h ( φ ) + β 1 λ 1 ψ = 0
(29) u = ψ e β 1 t γ 2 t β 1 φ β 1 γ 2 β 1 2 φ g β 1 2 ψ e β 1 φ γ 2 β 1 λ 1 + λ 2 2 γ 2 = 0
v = γ 2 t + φ ψ e β 1 φ γ 2 h β 1 2 ψ e β 1 φ γ 2 β 1 λ 1 + λ 2 2 ψ = 0
It is impossible to find the general solution of the systems from Table 1 for arbitrary functions g and h. However, if we correctly specify these functions we can find the solutions of these systems.
System (27) is the most interesting one from the point of view of applicability. Let us consider system (27) with τ 2 = 0 , g ( ω ) = b 1 ω + b 0 λ 1 + λ 2 2 , h ( ω ) = a 2 ω 2 + a 1 ω + a 0 . In this case, the reduction system has the form
φ + ( φ ) 2 β 2 ( b 1 β 2 ψ + b 0 ) = 0
ψ a 2 β 2 2 ψ 2 a 1 β 2 ψ a 0 = 0
The solution of Equation (41) has the form
ψ = k = c o n s t
Substituting (42) into (40), we arrive at
φ + ( φ ) 2 = A , A = β 2 b 1 β 2 k + b 0
The solutions of Equation (43) depend on the parameter A. Solving Equation (43) we get three different solutions (up to transformations x x + C 1 , C 1 = c o n s t )
φ = ln | C x | , A = 0 ; ln | C cosh ( A x ) | , A > 0 ; ln | C cos ( A x ) | , A < 0 .
Substituting φ and (42) into corresponding ansatz from Table 1, and renaming C β 2 C , we arrive at the exact solutions
u = C e β 2 t x + k , v = C β 2 e β 2 t x ,
u = C e β 2 t cosh A x + k , v = C β 2 e β 2 t cosh A x
u = C e β 2 t cos A x + k , v = C β 2 e β 2 t cos A x
of the reaction–diffusion system
u x x = λ 1 u t + v u β 2 a 1 + a 2 v a 2 β 2 u + b 0 λ 1 * b 1 v + b 1 β 2 u v + a 0 , v x x = λ 2 v t + β 2 b 0 λ 2 * b 1 v + b 1 β 2 u v
where k is the solution of the equation a 2 β 2 2 k 2 + a 1 β 2 k + a 0 = 0 .

4. Solutions and Their Properties of Some Generalization of the Lotka–Volterra System

Let us consider in detail the case A < 0 . Renaming β 2 = B 2 C 2 , b 1 = C 2 2 B 2 , b 0 = A 2 C 2 + B 2 λ 2 B 2 , a 0 = e 0 , a 1 = A 1 C 2 B 2 , a 2 = B 1 C 2 2 B 2 2 , we obtain the exact solution
u = C e B 2 C 2 t cos A x + k , v = B 2 C C 2 e B 2 C 2 t cos A x
where A = k B 2 C 2 + A 2 C 2 + B 2 λ 2 C 2 , and k is the solution of B 1 k 2 + A 1 k + e 0 = 0 , of the reaction–diffusion system
λ 1 u t = u x x + u A 1 + B 1 u + C 1 v + e 2 v 2 + e 1 v + e 0 , λ 2 v t = v x x + v A 2 + B 2 u + C 2 v
where C 1 = 2 B 1 B 2 1 C 2 , e 2 = B 1 B 2 C 2 2 B 2 2 , e 1 = A 1 A 2 C 2 B 2 + λ 1 * λ 2 * .
System (47) is the generalized Lotka–Volterra system. With e 2 = e 1 = e 0 = 0 system (47) becomes the classical Lotka–Volterra system
λ 1 u t = u x x + u A 1 + B 1 u + C 1 v , λ 2 v t = v x x + v A 2 + B 2 u + C 2 v
Note that exact solutions of the form (46) for the classical Lotka–Volterra system (48) have been found in [24].
Figure 1. Exact solution to (51).
Figure 1. Exact solution to (51).
Symmetry 07 01841 g001
System (48) can be obtained from system (47) with e 0 = 0 , B 1 = B 2 , A 1 = A 2 C 2 + λ 2 * λ 1 * B 2 C 2 . Also, the coefficients of (46) and (48) must satisfy the equation k k + A 2 B 2 + λ 2 * λ 1 * C 2 = 0 .
It is well known [3] that three main kinds of interactions between two biological species are simulated by system (48):
(i) predator u–prey v interaction,
(ii) competition of the species,
(iii) mutualism or symbiosis.
It turns out that solution (46) can describe the predator-prey interaction on the space interval [ l , l ] , (here l = π 2 A ) provided that
B 2 < 0 , C 2 < 0 , C < 0 , k > | C |
One can easily check that solution (46) is non-negative, bounded in the domain Ω = { ( t , x ) ( 0 , + ) × ( l , l ) } and satisfies the given Dirichlet boundary conditions, i.e.,
u | x = l = k , v | x = l = 0 , u | x = l = k , v | x = l = 0
Choosing the coefficients λ 1 = 2 , λ 2 = 1 , A 2 = 1 , B 2 = 1 , C 2 = 1 3 , B 1 = 0 , C = 1 3 , e 0 = 1 , gives that A 1 = 2 , C 1 = 1 3 , k = 1 2 . Thus, from solution (46) we obtain the solution
u = 1 2 1 3 e 3 t cos ( 7 2 x ) , v = e 3 t cos ( 7 2 x )
of the system
2 u t = u x x + u 2 + v 3 + v 2 9 + 1 , v t = v x x + v 1 u v 3
which can describe predator u–prey v interaction, as its coefficients satisfy the conditions for this type of the interaction [3]. System (52) is some generalization of the Lotka–Volterra system (48) with additional nonlinearity v 2 9 + 1 in the first equation.
Figure 2. Exact solution of (55).
Figure 2. Exact solution of (55).
Symmetry 07 01841 g002
Solution (51) satisfies Dirichlet boundary conditions (50) with l = π 14 , k = 1 2 .
As an example, we present solution (51) in Figure 1. This solution can describe the predator u–prey v interaction between the species u and v when population of predator u becomes 1 2 and prey eventually dies, i.e., ( u , v ) ( 1 2 , 0 ) as t + .
Choosing coefficients λ 1 = 3 2 , λ 2 = 5 4 , A 2 = 1 , B 2 = 4 , C 2 = 1 2 , B 1 = 1 , C = 1 4 , e 0 = 0 , we get A 1 = 1 , C 1 = 3 4 , k = 1 . Renaming t 1 4 t , from solution (46) we obtain the solution
u = 1 1 4 e 2 t cos 7 x , v = 2 e 2 t cos 7 x
of the system
6 u t = u x x + u 1 + u + 3 4 v + 5 64 v 2 , 5 v t = v x x + v 1 4 u v 2
which can also describe the predator u–prey v interaction, as its coefficients satisfy conditions for this type of interaction [3]. System (52) is some generalization of Lotka–Volterra system (48) with additional nonlinearity 5 v 2 64 in the first equation.
Solution (53) satisfies Dirichlet boundary conditions (50) with l = π 2 7 and k = 1 . This solution can describe the predator u–prey v interaction between the species u and v when population of predator u becomes 1 and prey eventually die, i.e., ( u , v ) ( 1 , 0 ) as t + .
If we consider system (48) with solution (46), then we obtain the solution that can describe competition of the species. Such a solution is presented in [24].
Also, system (47) can describe mutualism—or symbiosis—of two species. Choosing the coefficients λ 1 = 2 , λ 2 = 1 , A 2 = 5 , B 2 = 2 , C 2 = 1 2 , B 1 = 1 10 , k = 0 , C = 1 , e 0 = 0 , we obtain A 1 = 9 , C 1 = 11 20 . So, from solution (46) we obtain the solution
u = e 4 t cos x , v = 4 e 4 t cos x
of the system
2 u t = u x x + u 9 u 10 + 11 v 20 21 v 2 160 , v t = v x x + v 5 v 2 + 2 u
which is a generalization of Lotka–Volterra system (48) with additional nonlinearity 21 160 v 2 in the first equation.
Solution (55) satisfies Dirichlet boundary conditions (50) with l = π 2 , k = 0 . As an example, we present solution (55) in Figure 2. Solution (55) can describe the type of the interaction between the species u and v when both populations grow unboundedly, i.e., ( u , v ) ( + , + ) if t + .

5. Conclusions

In this paper, the nonlinear RD system (2) was examined in order to find the Q-conditional operators under which this system is invariant and to construct exact solutions. Because the system of differential (7) is too complicated, we were unable (and believe it is not possible) to find all the solutions of the determining system (7) and thence to find all possible Q-conditional operators. We have found the Q–conditional operators with restrictions (8) (in the case η u 2 = 0 we have found all possible systems and operators, in the case η u 2 0 we have presented some examples) with respect to which the reaction–diffusion system of equations with constant diffusion (2) is invariant. All these operators are given in Theorem 3 of Section 2. In Section 3 the ansätze for all Q-conditional operators of Theorem 3 and the reduction systems are constructed. Section 4 contains the solutions of some generalization of the Lotka–Volterra system. These solutions are analyzed in order to present of biological interpretation. Some graphs of obtained solutions are also presented. It is shown that the obtained solutions satisfy Dirichlet boundary conditions, which are typical for biological interpretation.

Acknowledgments

The author is grateful to the referees for the useful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Turing, A.M. The chemical basis of morphogenesis. Philos. Trans. Royal Soc. Lond. B Biol. Sci. 1952, 237, 37–72. [Google Scholar] [CrossRef]
  2. Ames, W.F. Nonlinear Partial Differential Equations in Engineering; Academic Press: New York, NY, USA; London, UK, 1965; p. xii+511. [Google Scholar]
  3. Murray, J.D. Mathematical Biology I: An Introduction, 3rd ed.; Interdisciplinary Applied Mathematics; Springer-Verlag: New York, NY, USA, 2002; Volume 17, p. xxiv+551. [Google Scholar]
  4. Murray, J.D. Mathematical biology II: Spatial Models and Biomedical Applications, 3rd ed.; Interdisciplinary Applied Mathematics; Springer-Verlag: New York, NY, USA, 2003; Volume 18, p. xxvi+811. [Google Scholar]
  5. Britton, N.F. Essential Mathematical Biology; Springer Undergraduate Mathematics Series; Springer-Verlag London, Ltd.: London, UK, 2003; p. xvi+335. [Google Scholar]
  6. Okubo, A.; Levin, S.A. Diffusion and Ecological Problems: Modern Perspectives, 2nd ed.; Interdisciplinary Applied Mathematics; Springer-Verlag: New York, NY, USA, 2001; Volume 14, p. xx+467. [Google Scholar]
  7. Henry, D. Geometric Theory of Semilinear Parabolic Equations; Lecture Notes in Mathematics; Springer-Verlag: Berlin, Germany; New York, NY, USA, 1981; Volume 840, p. iv+348. [Google Scholar]
  8. Cherniha, R. Lie symmetries of nonlinear two-dimensional reaction–diffusion systems. Rep. Math. Phys. 2000, 46, 63–76. [Google Scholar] [CrossRef]
  9. Cherniha, R.; King, J.R. Lie symmetries of nonlinear multidimensional reaction–diffusion systems. I. J. Phys. A 2000, 33, 267–282. [Google Scholar] [CrossRef]
  10. Cherniha, R.; King, J.R. Addendum: “Lie symmetries of nonlinear multidimensional reaction–diffusion systems. I”. J. Phys. A 2000, 33, 7839–7841. [Google Scholar] [CrossRef]
  11. Cherniha, R.; King, J.R. Lie symmetries of nonlinear multidimensional reaction–diffusion systems. II. J. Phys. A 2003, 36, 405–425. [Google Scholar] [CrossRef]
  12. Barannik, T.A. Conditional symmetry and exact solutions of a multidimensional diffusion equation. Ukr. Math. J. 2002, 54, 1416–1420. [Google Scholar]
  13. Barannyk, T. Symmetry and exact solutions for systems of nonlinear reaction–diffusion equations. Available online: http://eqworld.ipmnet.ru/en/solutions/interesting/barannyk.pdf (accessed on 13 October 2015).
  14. Serov, N.I. Conditional invariance and exact solutions of a nonlinear heat equation. Ukr. Math. J. 1990, 42, 1370–1376. [Google Scholar]
  15. Fushchych, W.; Shtelen, W.; Serov, M. Symmetry Analysis and Exact Solutions of Equations of Nonlinear Mathematical Physics; Kluwer: Dordrecht, The Netherland, 1993. [Google Scholar]
  16. Nucci, M.C. Symmetries of linear, C-integrable, S-integrable, and nonintegrable equations. In Nonlinear Evolution Equations and Dynamical Systems (Baia Verde, 1991); World Sci. Publ.: River Edge, NJ, USA, 1992; pp. 374–381. [Google Scholar]
  17. Clarkson, P.A.; Mansfield, E.L. Symmetry reductions and exact solutions of a class of nonlinear heat equations. Physica D 1994, 70, 250–288. [Google Scholar] [CrossRef]
  18. Arrigo, D.J.; Hill, J.M.; Broadbridge, P. Nonclassical symmetry reductions of the linear diffusion equation with a nonlinear source. IMA J. Appl. Math. 1994, 52, 1–24. [Google Scholar] [CrossRef]
  19. Arrigo, D.J.; Hill, J.M. Nonclassical symmetries for nonlinear diffusion and absorption. Stud. Appl. Math. 1995, 94, 21–39. [Google Scholar] [CrossRef]
  20. Pucci, E.; Saccomandi, G. Evolution equations, invariant surface conditions and functional separation of variables. Physica D 2000, 139, 28–47. [Google Scholar] [CrossRef]
  21. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Applied Mathematical Sciences; Springer: New York, NY, USA, 2010; Volume 168. [Google Scholar]
  22. Bluman, G.W.; Cole, J.D. The general similarity solution of the heat equation. J. Math. Mech. 1968/69, 18, 1025–1042. [Google Scholar]
  23. Cherniha, R.; Pliukhin, O. New conditional symmetries and exact solutions of reaction–diffusion systems with power diffusivities. J. Phys. A 2008, 41, 185208:1–185208:14. [Google Scholar] [CrossRef]
  24. Cherniha, R.; Davydovych, V. Conditional symmetries and exact solutions of the diffusive Lotka–Volterra system. Math. Comput. Model. 2011, 54, 1238–1251. [Google Scholar] [CrossRef]
  25. Cherniha, R. Conditional symmetries for systems of PDEs: New definitions and their application for reaction–diffusion systems. J. Phys. A 2010, 43, 405207:1–405207:13. [Google Scholar] [CrossRef]
  26. Sidorov, A.F.; Shapeev, V.P.; Yanenko, N.N. Metod Differentsialnykh Svyazei i Ego Prilozheniya v Gazovoi Dinamike; “Nauka” Sibirsk. Otdel.: Novosibirsk, Russia, 1984; p. 272. (In Russian) [Google Scholar]
  27. Carini, M.; Fusco, D.; Manganaro, N. Wave-like solutions for a class of parabolic models. Nonlinear Dynam. 2003, 32, 211–222. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Pliukhin, O. Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems. Symmetry 2015, 7, 1841-1855. https://doi.org/10.3390/sym7041841

AMA Style

Pliukhin O. Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems. Symmetry. 2015; 7(4):1841-1855. https://doi.org/10.3390/sym7041841

Chicago/Turabian Style

Pliukhin, Oleksii. 2015. "Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems" Symmetry 7, no. 4: 1841-1855. https://doi.org/10.3390/sym7041841

APA Style

Pliukhin, O. (2015). Q-Conditional Symmetries and Exact Solutions of Nonlinear Reaction–Diffusion Systems. Symmetry, 7(4), 1841-1855. https://doi.org/10.3390/sym7041841

Article Metrics

Back to TopTop