Next Article in Journal
A Novel Approach to Multi-Attribute Group Decision-Making with q-Rung Picture Linguistic Information
Next Article in Special Issue
Some Approaches to the Calculation of Conservation Laws for a Telegraph System and Their Comparisons
Previous Article in Journal
Inflation in Mimetic f(G) Gravity
Previous Article in Special Issue
Lie and Q-Conditional Symmetries of Reaction-Diffusion-Convection Equations with Exponential Nonlinearities and Their Application for Finding Exact Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lie Symmetries of Nonlinear Parabolic-Elliptic Systems and Their Application to a Tumour Growth Model

1
Institute of Mathematics, National Academy of Sciences of Ukraine, 3, Tereshchenkivs’ka Street, 01601 Kyiv, Ukraine
2
School of Mathematical Sciences, University of Nottingham, University Park, Nottingham NG7 2RD, UK
*
Author to whom correspondence should be addressed.
Symmetry 2018, 10(5), 171; https://doi.org/10.3390/sym10050171
Submission received: 25 April 2018 / Revised: 13 May 2018 / Accepted: 14 May 2018 / Published: 17 May 2018

Abstract

:
A generalisation of the Lie symmetry method is applied to classify a coupled system of reaction-diffusion equations wherein the nonlinearities involve arbitrary functions in the limit case in which one equation of the pair is quasi-steady but the other is not. A complete Lie symmetry classification, including a number of the cases characterised as being unlikely to be identified purely by intuition, is obtained. Notably, in addition to the symmetry analysis of the PDEs themselves, the approach is extended to allow the derivation of exact solutions to specific moving-boundary problems motivated by biological applications (tumour growth). Graphical representations of the solutions are provided and a biological interpretation is briefly addressed. The results are generalised on multi-dimensional case under the assumption of the radially symmetrical shape of the tumour.

1. Introduction

It is well-known that nonlinear reaction-diffusion (RD) systems (i.e., systems of second-order parabolic PDEs with reaction terms) are used to describe various processes in physics, chemistry, biology, ecology, etc. (see, e.g., books [1,2,3,4]). Indeed, there are many special cases of such systems describing real world processes, among the most famous being Lotka–Volterra type systems [5,6] and systems used for modelling the chemical basis of morphogenesis [7] (see, e.g., ref. [2] for details in modern terminology).
The most common are two-component RD systems of the form
U t = ( D 1 ( U ) U x ) x + F ( U , V ) , V t = ( D 2 ( V ) V x ) x + G ( U , V ) ,
where U ( t , x ) and V ( t , x ) are unknown functions, the nonnegative smooth functions D 1 ( U ) and D 2 ( V ) are coefficients of diffusivity (conductivity), and F and G are arbitrary smooth functions describing the interaction between components U ( t , x ) and V ( t , x ) . Hereafter the t and x subscripts denote differentiation with respect to these variables.
It should be noted that the search of Lie symmetries of RD systems was initiated many years ago and probably paper [8] was the first in this direction. A complete Lie symmetry classification was done much later because this problem is much more complicated in comparison with the case of a single RD equation. In particular, it was necessary to treat separately the cases of constant diffusivities (i.e., semilinear equations) and nonconstant (quasilinear) ones. All possible Lie symmetries of (1) with constant diffusivities were completely described in [9,10,11]. In the case of nonconstant diffusivities, it has been done in [12,13] (in [14], the result was extended to the multi-component RD systems).
It should be stressed that there are models used to describe real world processes that are based on singular limits of a RD system (1), namely with D 2 ( V ) = 0 or with the V t term absent. The system (1) with D 2 ( V ) = 0 has been used to describe interactions involving species (cells) V, which are not able to diffuse in space (typical examples include some models for plankton). Some particular results about Lie symmetries of such systems have been recently obtained in [15], however the problem of a complete Lie symmetry classification is still not solved.
Here, however, we examine the class of systems (1) with the V t term negligible so that V is governed by a quasi-steady equation, i.e.,
U t = ( D 1 ( U ) U x ) x + F ( U , V ) , 0 = ( D 2 ( V ) V x ) x + G ( U , V ) .
Obviously this system is reduced to the form
U t = ( D 1 ( U ) U x ) x + F * ( U , V * ) , 0 = V x x * + G * ( U , V * )
by the Kirchhoff substitution for the component V V * . Here the diffusivity D 1 ( U ) is a nonnegative and nonconstant function (at least in some domains) in what follows. We remind the reader that the case D 1 ( U ) = c o n s t was studied in [9,10]. Thus, the class of parabolic-elliptic systems (3) with D 1 ( U ) c o n s t is the main object of this study. The main motivation for this follows from paper [16], in which a model for tumour growth was derived. Under some biologically motivated assumptions (see p. 570 in [16] for details), the model was simplified to a boundary-value problem with governing equations of the form (3). System (2) is of course potentially of relevance to any two-component RD model in which one species diffuses much faster than the other; it is also of interest purely from the symmetry point of view, since its symmetries differ essentially from those of (1).
The paper is organized as follows. Section 2 is devoted to the search for form-preserving transformations for the class of systems (3) in order to establish possible relations between RD systems that admit equivalent Lie symmetry algebras. Section 3 is devoted to the identification of all possible Lie symmetries that any system of the form (3) can admit. The main results of this section are presented in the form of Table 1 and Table 2. In Section 4 and Section 5, we apply the Lie symmetries derived for the reduction of a nonlinear boundary value problem (BVP) with a free boundary modeling tumour growth in order to construct its exact solution. Moreover, possible biological interpretations of the results derived are discussed. Finally, we briefly discuss the result obtained and present some conclusions in the last section.

2. Form-Preserving Transformations for the Class of Systems (3)

First of all we simplify notations in system (3) in a natural way:
U t = ( D ( U ) U x ) x + F ( U , V ) , 0 = V x x + G ( U , V ) .
The class of systems (4) contains three arbitrary functions and the Lie symmetry of its different representatives depends essentially on the form of the triplet ( D , F , G ) . Thus, the problem of a complete description of all possible Lie symmetries (the so called group classification problem) arises. In order to solve this, we can apply the Lie–Ovsiannikov approach (the name of Ovsiannikov arises because he published a remarkable paper in this direction, [17]) of the Lie symmetry classification, which is based on the classical Lie scheme and a set of equivalence transformations of the differential equation in question. However, it is well-known that this approach leads to a very long list of equations with nontrivial Lie symmetry provided the given equation (system) contains several arbitrary functions (system (4) involves three such).
During the last two decades, new approaches for solving group classification problems were developed, which are important for obtaining the so called canonical list of inequivalent equations admitting nontrivial Lie symmetry algebras and allowing the solution of this problem in a more efficient way than a formal application of the Lie–Ovsiannikov approach. Here we use the algorithm based on so called form-preserving transformations [18,19], which were used initially for finding locally-equivalent PDEs, especially those nonlinear PDEs that are linearizable by a point transformation. Interestingly, such transformations were implicitly used much earlier in [20] in order to find all possible heat equations with nonlinear sources that admit the Lie symmetry either of the linear heat equation or the Burgers equation. In [21], these transformations are called ‘allowed transformations’ and they were used to classify the Lie symmetries of a class of variable coefficient Korteweg-de Vries equations. One may say that the above cited paper was the first one, in which allowed (form-preserving) transformations were systematically used for solving the Lie symmetry classification problem. Now we present the definition of form-preserving transformations in the case of the class of systems (4).
Definition 1.
A point non-degenerate transformation given by the formulae
t = a ( τ , y , u , v ) , x = b ( τ , y , u , v ) , U = φ ( τ , y , u , v ) , V = ψ ( τ , y , u , v ) ,
which maps at least one system of the form (4) into a system belonging to the same class, is called a form-preserving transformation.
Comparing this definition with the well-known definition for equivalence transformations, one immediately notes that each equivalence transformation is automatically a form-preserving transformation, but not vice versa.
It was noted later [22,23] that the form-preserving transformations (nowadays the terminology ‘admissible transformations’ is also widely used) allow an essential reduction of the number of cases obtained via the classical Lie–Ovsiannikov algorithm (see extensive discussions on this matter in [11,13,24]). In particular, it is shown in [13] that there are only 10 inequivalent cases when RD systems with nonconstant diffusivities and reaction terms of the form (1) are invariant with respect to the nontrivial Lie algebras (in contrast to a much longer list derived in [12] via the Lie–Ovsiannikov algorithm). Similarly, it was proved in [24] that the canonical list of inequivalent reaction-diffusion-convection equations consists of 15 equations only (not the 30 derived by the Lie–Ovsiannikov algorithm).
The reader may find more details (including specific examples) on this matter in Section 2.3 of the recent book [25].
Now we present the main result of this section, namely the theorem describing a general form of form-preserving transformations for the class of systems (4).
Theorem 1.
An arbitrary system of the form
u τ = D ˜ ( u ) u y y + F ˜ ( u , v ) , 0 = v y y + G ˜ ( u , v )
can be reduced to a system from class (4) by a form-preserving transformation of the form (5) if and only if the smooth functions a , b , φ and ψ have the forms
a = α ( τ ) , b = β ( τ , y ) , α τ β y 0 , φ = K ( τ , y ) u + P ( τ , y ) , K 0 , ψ = L ( τ , y ) v + Q ( τ , y ) , L 0 ,
where the functions α ( τ ) , β ( τ , y ) , K ( τ , y ) , P ( τ , y ) , L ( τ , y ) and Q ( τ , y ) are such that the following equalities
α τ D = β y 2 D ˜ , α τ F = K F ˜ + ( K y u + P y ) 2 K D ˜ u + K τ u + P τ D ˜ K y y u + P y y 2 K y u + P y K K y , β y 2 G = L G ˜ + 2 L y v + Q y L L y L y y v Q y y , 2 β y ( K y u + P y ) D ˜ u + D ˜ 2 β y K y β y y K + β τ K = 0 , 2 β y L y β y y L = 0 ,
hold.
The proof of Theorem 1
is quite similar to that for the class of RD systems (1) with D k = c o n s t ( k = 1 , 2 ) presented in [26].
The group of equivalence transformations can be easily extracted from Theorem 1 by assuming that D ˜ , F ˜ and G ˜ are arbitrary smooth functions. Solving the system (8) under such condition one arrives at the following statement. ☐
Consequence 1.
The group of equivalence transformations for (4) has the form
t = α 1 τ + α 2 , x = α 3 y + α 4 , U = α 5 u + α 6 , V = α 7 v + α 8 , D = α 3 2 α 1 D ˜ , F = α 5 α 1 F ˜ , G = α 7 α 3 2 G ˜ ,
where α l ( l = 1 , , 8 ) are arbitrary group parameters ( α 1 > 0 , α 2 k 1 0 , k = 2 , 3 , 4 ).
One sees that (8) is a nonlinear system of functional-differential equations and its general solution seems to be impossible without further restrictions. For example, assuming a constant diffusivity D ˜ , one may derive that β ( τ , y ) is linear with respect to (w.r.t.) the variable y; however, β ( τ , y ) can be nonlinear for the nonconstant diffusivity D ˜ (see Table 3 and Table 4 below for examples).
While the form of the transformations (7) is still quite general, it will be shown in the next section that only particular cases play important roles in solving the Lie symmetry classification problem for the nonlinear parabolic-elliptic systems (3).

3. Lie Symmetries of a Class of Parabolic-Elliptic Systems

In order to find Lie symmetry operators of a system of the form (4) via the Lie method [25,27,28,29,30,31], one needs to consider the manifold ( S 1 , S 2 )
S 1 U t ( D ( U ) U x ) x F ( U , V ) = 0 , S 2 V x x + G ( U , V ) = 0
in the space of the following variables:
t , x , U , V , U t , U x , V t , V x , U t x , U t t , U x x , V t x , V t t , V x x .
The system (4) is invariant under the transformations generated by the infinitesimal operator
X = ξ 0 ( t , x , U , V ) t + ξ 1 ( t , x , U , V ) x + η 1 ( t , x , U , V ) U + η 2 ( t , x , U , V ) V
when the following invariance conditions are satisfied:
X 2 S 1 X 2 U t ( D ( U ) U x ) x F ( U , V ) | S 1 = 0 , S 2 = 0 = 0 , X 2 S 2 X 2 V x x + G ( U , V ) | S 1 = 0 , S 2 = 0 = 0 .
The operator X 2 is the second prolongation of the operator X, i.e.,
X 2 = X + ρ t 1 U t + ρ t 1 V t + ρ x 1 U x + ρ x 2 V x + σ t t 1 U t t + σ t x 1 U t x + σ x x 1 U x x + σ t t 2 V t t + σ t x 2 V t x + σ x x 2 V x x ,
where the coefficients ρ and σ with relevant subscripts are calculated by well-known formulae (see, e.g., [29,30,31]).
Substituting (13) into (12) and eliminating the derivatives U t and V x x using (4), we can obtain the following system of determining equations (DEs):
ξ x 0 = ξ U 0 = ξ V 0 = ξ U 1 = ξ V 1 = 0 ,
η V 1 = η U 2 = η V V 2 = 0 ,
2 ξ x 1 ξ t 0 D = η 1 D U ,
2 η x V 2 = ξ x x 1 ,
ξ x x 1 2 η x U 1 D = 2 η x 1 D U + ξ t 1 ,
2 ξ x 1 ξ t 0 η U 1 D U = η U U 1 D + η 1 D U U ,
η 1 F U + η 2 F V = η t 1 D η x x 1 + F η U 1 ξ t 0 ,
η 1 G U + η 2 G V = η x x 2 + G η V 2 2 ξ x 1 .
Obviously, the general solution of the above system of DEs essentially depends on the form of the triplet ( D , F , G ) . Thus, the problem of a complete description of all possible Lie symmetries (the so called group classification problem) arises.
First of all we find the so called trivial algebra (other terminology used for this algebra is the ‘principal algebra’ and the ‘kernel of maximal invariance algebras’) of the class in question, i.e., the maximal invariance algebra (MAI) admitted by each equation of the form (4). Under the natural assumption that D , F and G are arbitrary smooth functions the above system of DEs can be easily solved and the two-dimensional Lie algebra generated by the basic operators
X 1 = t , X 2 = x
is obtained.
In order to solve the problem of Lie symmetry classification for system (4), i.e., to find all possible systems admitting three- and higher-dimensional MAI, we use the algorithm based on form-preserving transformations, which was briefly discussed in Section 2. Moreover, it will be shown below that, similarly to the case of standard nonlinear RD systems [13], such an approach is also more efficient for parabolic-elliptic systems (compared to the classical Lie–Ovsiannikov algorithm).
Let us formulate a theorem which gives complete information on the classical symmetry of the nonlinear system (4).
Theorem 2.
All possible maximal algebras of invariance (up to equivalent representations generated by transformations of the form (23)) of the system (4) for any fixed nonconstant function D and nonconstant function vectors ( F , G ) are presented in Table 1 and Table 2. Any other system of the form (6) with nontrivial Lie symmetry is reduced by a local substitution of the form
t = C 0 τ + C 1 exp ( C 2 τ ) , x = C 3 y + C 4 exp ( C 5 y ) + C 6 tan ( C 7 y ) , U = C 8 + C 9 τ + C 10 exp ( C 11 τ ) u + C 4 exp ( C 12 y ) u + C 6 cos 3 ( C 13 y ) u , V = C 14 + C 15 τ + C 16 exp ( C 17 τ ) v + C 4 exp ( C 12 y ) v + C 6 cos 1 ( C 13 y ) v + C 18 τ y 2 + C 19 exp ( C 20 τ ) y 2 ,
to one of those given in Table 1 and Table 2 (the constants C with subscripts are determined by the form of the system in question, some of them necessarily being zero in any given case).
Remark 1.
Each case in Table 1 contains semi-coupled systems (i.e., those with F V G U = 0 ) as special subcases. A system specified in such a way has the same MAI provided it is not listed in Table 2 as a special case. For example, the system arising in case 1 with an arbitrary function f ( U ) and g ( U ) = c o n s t is a semi-coupled system; however, this is not listed in Table 2, so that its MAI is still a 3-dimensional Lie algebra.
Remark 2.
In Table 1 and Table 2, α , α 0 , α 1 , α 2 , β , γ , k are arbitrary constants obeying the restriction k β ( ( α 1 ) 2 + ( α 2 ) 2 ) 0 ; D , f , g , φ and ψ are arbitrary smooth functions of the relevant arguments; the function h ( t , x ) is the general solution of the equation h x x + α 2 h = 0 , while h 0 ( x ) = 1 α 2 , i f α 2 0 , x 2 2 , i f α 2 = 0 .
The sketch of the proof of Theorem 2
Consider the determining Equations (14)–(21) under the assumption D U 0 (we remind the reader that the case D U = 0 was examined in [9,10]). Obviously, Equations (14)–(17) can be easily solved and the formulae
ξ 0 = a ( t ) , ξ 1 = b ( t , x ) , η 1 = ( 2 b x a ˙ ) D D U 1 , η 2 = r ( t , x ) V + p ( t , x ) , r x = 1 2 b x x
are obtained (here a , b , r and p are arbitrary smooth functions and the upper dot means the time derivative).
Now we can consider Equations (18)–(21) as classification equations to find the function D and the pairs of the functions ( F , G ) for which the system (4) has a nontrivial Lie symmetry, i.e., its MAI is larger than the trivial algebra (22).
It follows from (18)–(19) that there is a special case 2 b x = a t when the diffusivity D can be an arbitrary smooth function. Indeed, Equation (18) produces b t = b x x = 0 provided D is an arbitrary function. Hence, one arrives at b = b 1 x + b 0 , a = 2 b 1 t + a 0 , so that the classification Equations (20) and (21) have the form:
( r V + p ) F V = 2 b 1 F , ( r V + p ) G V + G 2 b 1 r = p x x .
The general solution of (25) depends essentially on the functions r and p ( r 2 + p 2 0 , otherwise a trivial Lie symmetry is obtained). In order to derive a complete result, one needs to examine two subcases.
Let us consider the first subcase r = 0 , p 0 , i.e., the first equation in (25) takes the form p F V = 2 b 1 F . A simple analysis says that the general solution of (25) leads to the functions F , G and coefficients (24), which produce the system and MAI arising in case 1 of Table 1 provided b 1 0 . If b 1 = 0 then case 3 of Table 2 is derived.
The second subcase r 0 leads to case 2 of Table 1 and cases 1 and 2 of Table 2.
The condition 2 b x a t 0 is a generic inequality and produces all other cases listed in Table 1 and Table 2. Using this condition and (24), one easily solves the Equation (19) and obtains
D = d ( U + C 1 ) k , d exp ( C 2 U ) ,
where d 0 , C 1 , C 2 and k 0 are arbitrary constants. Moreover, substituting (26) into (18), we obtain b = b 1 x + b 0 , b 0 , b 1 R , unless
D = d ( U + C 1 ) 4 / 3 , b = b ( x ) .
One notes that the system (4) with diffusivities (26) and (27) can be simplified to one with d = 1 , C 1 = 0 , C 2 = 1 by the appropriate equivalence transformations from (9).
Now one needs to examine other determining equations separately with power law diffusivity, exponential diffusivity and the special power law coefficient k = 4 / 3 . We present here only the last case, i.e., the power low diffusivity D = U 4 / 3 , which is the most complicated.
Let us consider system (4) with D = U 4 / 3 . It turns out that we may assume that ξ x x 1 0 (otherwise only particular cases of (4) with D = U k will be derived). The classification Equations (20) and (21) for D = U 4 / 3 take the form
3 4 a ˙ 2 b U F U + ( r V + p ) F V + a ˙ 4 + 3 b 2 F = 3 4 a ¨ U + 3 b ( 3 ) 2 U 1 / 3 ,
3 4 a ˙ 2 b U G U + ( r V + p ) G V + 2 b r G = b ( 3 ) 2 V p x x .
(because the function b depends only on x we use the standard notations for its derivatives and b ( k ) = d k b d x k , k = 3 , 4 ). Now we construct differential consequences of (28) and (29) w.r.t. x:
3 b 2 U F U + b 2 V + p x F V + 3 b 2 F = 3 b ( 4 ) 2 U 1 / 3 ,
3 b 2 U G U + b 2 V + p x G V + 3 b 2 G = b ( 4 ) 2 V p x x x .
Taking into account that b 0 one easily finds the general solution of (30) and (31),
F = U f ( ω ) + 3 b ( 4 ) 4 b U 1 / 3 , ω = U V + 2 p x b 3 ,
G = U g ( ω ) b ( 4 ) 4 b V + 2 p x b + 2 3 b p x b ( 4 ) p x x x .
where f and g are arbitrary smooth functions.
The further analysis essentially uses the fact that the functions F and G depend only on U and V (not explicitly on t and x!). This means that relevant constraints on the functions f , g , b and p must take place.
Let us consider the most interesting case when the functions f and g are still arbitrary. This happens if and only if the constraints
p x b = α 1 , b ( 4 ) b = α 2 , 2 3 b p x b ( 4 ) p x x x = α 3
(here α 1 , α 2 and α 3 are some constants) hold in (32) and (33). Because the functions (32) and (33) must satisfy Equations (28) and (29) for arbitrary functions f and g, the additional condition p = γ r springs up. Direct calculations show that γ = c o n s t arising in the system obtained can be reduced to γ = 0 by the equivalence transformation V + γ V . Thus, we arrive at system (4) with D = U 4 / 3 and
F = U f ( ω ) + 3 α 4 U 1 / 3 , ω = U V 3 ,
G = U g ( ω ) α 4 V ,
i.e., the parabolic-elliptic system
U t = ( U 4 / 3 U x ) x + U f U V 3 + 3 α 4 U 1 / 3 , 0 = V x x + U g U V 3 α 4 V ,
where f and g are arbitrary smooth functions, α is an arbitrary constant. It turns out that MAI of system (36) essentially depends on the parameter α .
Setting α = 0 and substituting (34) and (35) into (28) and (29), one obtains the system
a ˙ = 0 , r = b 2 , b ( 3 ) = 0 ,
in order to find MAI of (36). The linear ODE system (37) can be easily solved, and as a result case 8 of Table 1 is obtained.
Setting α 0 , one notes that this parameter can be easily reduced to ± 1 . If α < 0 , then the equivalence transformation t = 4 α τ , x = 2 α y reduces (36) to the system listed in case 3 of Table 3. If α > 0 , then t = 4 α τ , x = 2 α y produces the system from case 4 of Table 3. It turns out that both systems can be reduced to one listed in case 8 of Table 1. The relevant transformations are
x = tan y , U = u cos 3 y , V = v cos 1 y ,
and
x = 1 2 e 2 y , U = e 3 y u , V = e y v .
Transformations (38) and (39) are presented in the third column of Table 3 (see cases 3 and 4) and they have been found by applying Theorem 1 to the nonlinear system (36). Obviously, Formulae (38) and (39) can be derived from (23) as particular cases.
Thus, the most general system with a special form of the diffusivity D from (27) and its MAI are presented in case 8 of Table 1 and any other system with such a diffusivity admitting a nontrivial Lie symmetry is reduced to one listed in case 8 of Table 1. It should be noted that the corresponding substitutions are highly nontrivial and were found using a special case of form-preserving transformations (not equivalence transformations) derived in Theorem 1. Now we need to find all possible correctly-specified functions f and g leading to extensions of MAI of system (36). Such functions have been found by analysis of the differential consequences F x V = F t V = 0 , G x U = G t U = 0 , where the functions F and G are of the form (32) and (33). Finally, we have established only two special cases leading to the above mentioned extensions.
(i) f = α 1 U γ V 3 γ , g = α 2 U γ + 1 / 3 V 3 γ + 1 , where α 1 , α 2 and γ are arbitrary constants ( α 1 2 + α 2 2 0 ). The systems (36) with these f and g are presented in case 9 of Table 1 and cases 5–6 of Table 3.
(ii) f = c o n s t , g = c o n s t . The constant f 0 leads to case 12 of Table 2 and cases 6–7 of Table 4. The zero constant f leads to case 13 of Table 2 and cases 8–9 of Table 4.
At the final stage, we note that the form-preserving transformations (38) and (39) are still working, so that only case 9 of Table 1 and cases 12–13 of Table 2 are essentially new.
Substituting all possible functions D from (26) (with d = 1 , C 1 = 0 , C 2 = 1 ) into (20) and (21) and taking into account (24), we construct the corresponding pairs ( F , G ) and Lie algebras. The results are summarised in the form of cases 3–7, 10–14 in Table 1 and cases 4–11, 14–21 in Table 2. All other parabolic-elliptic systems of the current class with nontrivial Lie symmetry can be reduced to those in Table 1 and Table 2 by the appropriate equivalence transformations from (9) and/or form-preserving transformations and these are listed in Table 3 and Table 4.
The sketch of the proof is now completed. ☐
Consequence 2.
All the nonlinear systems of the form (4) admitting a nontrivial Lie symmetry, which are reducible by FPTs to those presented in Table 1 and Table 2, are presented in Table 3 and Table 4, respectively.
The systems listed in the second column of Table 3 and Table 4 are presented in their simplified forms taking into account the equivalence transformations (9). The form-preserving transformations listed in the third column of these tables can be united to the form (23).
The above Consequence follows from the detailed proof of Theorem 2.
We would like to conclude this section by the following observation. As it was noted in Section 2, the form-preserving transformations are used in order to solve the problem of the Lie symmetry classification for the class of systems (4). We have proved that there are 35 inequivalent systems admitting nontrivial Lie algebras (i.e., their MAI is three- and higher-dimensional) and they are listed in Table 1 and Table 2. The systems are inequivalent up to point transformations of the form (23), which are particular cases of form-preserving transformations described in Theorem 1. It is a standard routine to show that there are no other point transformations that allow the reduction of a system from Table 1 or Table 2 to another system from these tables. Thus, a canonical list of nonlinear parabolic-elliptic systems of the form (4) admitting a nontrivial Lie symmetry consists of the 35 inequivalent systems listed in Table 1 and Table 2.
A natural question arises: How many systems of the form (4) with nontrivial Lie symmetry can be derived using the Lie–Ovsiannikov algorithm? This means that we should take into account only the group of equivalence transformations (9). A formal answer is very simple: one should extend the list of the above 35 systems by those from Table 3 and Table 4, therefore 57 systems with the three- and higher-dimensional MAI will be derived (in the sense that the relevant MAI are inequivalent up to the equivalence transformations (9)).

4. Boundary Value Problems for a One-Dimensional Tumour Growth Model with Negligible Cell Viscosity

In this section we apply the results of Lie symmetry classification derived in Section 2 for exact solving of a real world model proposed in [16] (see p. 570).
Consider the pressure difference with r = 1 [16] (see p. 569 therein)
Σ ( α ) = k 0 α α * ( 1 α ) q
assuming α > α * . Setting
k ( α ) = k 0 α 1 m ( 1 α ) 2 d ( α Σ ) d α
the equation (see (23) in [16]) for the tumour cell concentration, α ( t , x ) , takes the form
α t = α m α x x + S ( α , c ) .
In order to simplify the boundary conditions we set
U = α α * , V = c c
(wherein c corresponds to the level of nutrient in [16], but might instead correspond to that of a chemotherapeutic drug) and arrive at the BVP
U t = ( U + α * ) m U x x + S ( U + α * , c V ) , 0 = V x x + Q ( U + α * , c V ) , x = 0 : U x = V x = 0 , x = R ( t ) : U = V = 0 , x = R ( t ) : R ˙ = α * m 1 U x ,
where the concentrations U ( t , x ) and V ( t , x ) and the moving boundary location R ( t ) are unknown functions, while the functions S and Q are given smooth functions, and m, α * and c are nonnegative parameters.
One may note that the governing equations of this BVP with S = f ( U ) V β and Q = g ( U ) V β + 1 are a particular case of case 2 in Table 1, so that they admit the Lie symmetry operator
X = 2 β t t + β x x + 2 V V .
It turns out that the following statement can be easily proved, using the definition of invariance for BVPs with moving boundaries [32].
Theorem 3.
The nonlinear BVP (42) with S = f ( U ) V β and Q = g ( U ) V β + 1 is invariant with respect to the two-dimensional MAI generated by the Lie symmetry operators (43) and t .
Clearly the time translation operator t cannot help to find any realistic solutions of BVP (42), while the Lie symmetry (43) guarantees a reduction of the problem, which leads to some interesting results. In fact, the operator generates the ansatz
U = φ ( ω ) , ω = x t , V = t 1 β ψ ( ω ) ,
which immediately specifies the function R ( t ) = ω 0 t with the constant ω 0 > 0 to be found. The reduced ODE problem is
( φ + α * ) m φ + m ( φ + α * ) m 1 φ 2 + ω 2 φ + f ( φ ) ψ β = 0 , ψ + g ( φ ) ψ β + 1 = 0 , ω = 0 : φ = ψ = 0 , ω = ω 0 : φ = ψ = 0 , φ = α * 1 m ω 0 2 .
The nonlinear BVP (45) is still a complicated problem and we were unable to solve it analytically in general. Let us consider a special case, namely (45) with β = 1 :
( φ + α * ) m φ + m ( φ + α * ) m 1 φ 2 + ω 2 φ + f ( φ ) ψ 1 = 0 , ψ + g ( φ ) = 0 , ω = 0 : φ = ψ = 0 , ω = ω 0 : φ = ψ = 0 , φ = α * 1 m ω 0 2 .
In this case the exact solution can be found for arbitrary m in the form
φ = α * 1 m 4 ω 0 2 ω 2 , ω < ω 0 , ψ = α * 1 m q 0 48 ω 0 2 ω 2 5 ω 0 2 ω 2 , ω < ω 0
provided the functions f ( φ ) and g ( φ ) have the form
g = q 0 φ , q 0 > 0 , f = q 0 α * m 1 3 φ φ + α * 1 m ω 0 2 × ( m + 1 2 ) α * 1 m ( φ + α * ) m m α * 2 m 4 α * m ω 0 2 + 4 ( φ + α * ) m 1 φ + α * 1 m 4 ω 0 2 .
Remark 3.
A similar result can be obtained for arbitrary diffusivity in the governing equation (40), i.e., the drag coefficient k ( α ) , however, the structure of the function f ( φ ) will essentially depend on the diffusivity.
Now we present the final formulae for the original concentrations of tumour cells α ( t , x ) and nutrient (or drug) c ( t , x ) using Formulae (41), ansatz (44) and solution (47):
α ( t , x ) = α * + α * 1 m 4 ω 0 2 x 2 t , x < ω 0 t , c ( t , x ) = c α * 1 m q 0 48 t ω 0 2 x 2 t 5 ω 0 2 x 2 t , x < ω 0 t , R ( t ) = ω 0 t .
Formulae (47) give the exact solution of the nonlinear BVP with a moving boundary for the governing equations
α t = α m α x x + f ( α α * ) c c , 0 = c x x q 0 ( α α * ) ,
and the boundary conditions
x = 0 : α x = c x = 0 , x = R ( t ) : α = α * , c = c , x = R ( t ) : R ˙ = α x α * 1 m ,
where the function f is defined by (48), in particular one obtains the cubic polynomial for m = 0 :
f = q 0 3 α * ( α α * ) ( α α 1 ) ( α 2 α ) , α 1 = α * ( 1 ω 0 2 ) , α 2 = α * 6 + ω 0 2 4 ,
hence
S ( α , c ) = q 0 3 α * ( α α * ) ( α α 1 ) ( α 2 α ) c c .
This solution for some specified parameters is presented in Figure 1, Figure 2 and Figure 3. We note that the positive parameter α * < 1 can be interpreted as a natural concentration of tumour cells, while c is the concentration of nutrient (drug) in the medium surrounding the tumour. The positive parameter ω 0 may be found from additional biologically motivated conditions. Notably the exact solution (49) can be easily generalised to one with two arbitrary parameters by the time translation t t + t 0 .
It should be noted that the functions f ( φ ) and g ( φ ) in the cases m = 0 (constant diffusion) and m = 1 (as in the porous medium equation) are positive for the solution (47), meaning that the corresponding functions S ( α , c ) and Q ( α , c ) (see p. 570 in [16]) are positive (see an example in Figure 4). The cases m 0 , 1 lead to a much more complicated structure for f ( φ ) and need to be examined separately for correctly-specified values of m. For example, the functions S ( α , c ) and Q ( α , c ) are positive on the solution for m = 1 / 2 in the case of the parameters α * = 1 / 2 and ω 0 = 1 used in Figure 2.
We now turn to regimes in which the results may be susceptible to biological application, largely returning to the notation of [16]. The effective diffusivity of the cells is given by
D ( α ) = α k ( α ) ( 1 α ) 2 d d α α Σ .
The fact that cell-cell repulsion can be expected to increase at high densities implies that
d 2 d α 2 α Σ > 0 ,
but this increasing contribution to (54) may be offset by the ( 1 α ) 2 prefactor (which itself will be less significant because k ( 1 ) = 0 typically holds) that in part reflects the overall mass balance
α v c + ( 1 α ) v w = 0
between cell and water velocities, v c and v w , respectively (so when 1 α is small so typically is the cell velocity). The upshot is nevertheless that cases with m 0 are those of most relevance, with increasing m corresponding to stronger cell-cell repulsion at high densities.
S ( α , c ) will be positive when c represents a nutrient, with S an increasing function of c in both cases (corresponding to β > 0 in the above). In this case, S will be an increasing function of α for sufficiently small α α * (and linear in α ) and S = 0 at α = α * . One may note that the curves in Figure 5 satisfy these properties. In fact, all the curves represent increasing functions, the green and red curves in Figure 5 behave as linear functions, while the blue curve is linear for sufficiently small φ = α α * (the curves are built on the interval [ 0 , φ m a x ] with φ m a x prescribed by the first formula in (47)). Moreover, f ( φ ) = 0 for φ = 0 (see (48)). Notably, S may decrease again for larger α , reflecting the limited supply of water to make new cells. The function S defined by (53) possesses such property provided α > α m a x ( α * , α 2 ) , where α m a x is a local maximum of f from (52) and can be easily calculated.
In the drug case, S would typically be taken positive for small c and negative for large c (with S a decreasing function of c)—this behaviour cannot be captured by the above nonlinearities, so in interpretation of the results, c should be regarded as everywhere sufficiently large that S is negative, again with β > 0 . Various dependencies on α can then be identified, and even S = 0 at α = 0 need not hold (though additional constraints need to be imposed if α reaches zero).
Q should be positive in both interpretations (absorption, rather than generation, though the generation of toxic materials in the drug case could allow Q < 0 to be contemplated) and an increasing function of c and α . According to Theorem 3 Q = g ( U ) V β + 1 (with U and V given in (41)), therefore one needs β > 1 and an increasing function g in order to satisfy the above requirement.
Q > 0 implies by the maximum principle that c will be an increasing function of x, having a minimum at x = 0 . Such behaviour of c is pictured in Figure 1, Figure 2 and Figure 3 for different values of m. The qualitative form of α will depend on the initial data and the form of S; in the case S > 0 , the nutrient supply is highest at the edge of the tumour, but α is held at α * there by the stress balance (more water is taken up to prevent α exceeding α * due to cell division), so, perhaps counterintuitively, it is quite possible for α to be maximal at x = 0 if S remains sufficiently large, despite the decrease in c. However, such behaviour is more natural for S < 0 (the drug case) since the cell-kill rate can be expected to be highest at the edge of the tumour.

5. Higher-Dimensional Tumour Growth Model without Cell Viscosity

In Section 4, we examined the one-dimensional model (42) describing the tumour growth. Here we generalise the model to higher dimensions under the assumption that tumour cells create a radially symmetrical shape (a circular cylinder for n = 1 and a sphere for n = 2 ). In this case, the nonlinear BVP (42) takes the following form:
U t = 1 r n r n ( U + α * ) m U r r + S ( U + α * , c V ) , 0 = 1 r n r n V r r + Q ( U + α * , c V ) , r = 0 : U r = V r = 0 , r = R ( t ) : U = V = 0 , r = R ( t ) : R ˙ = α * m 1 U r ,
where n = 1 , 2 and the same notations are used as in Section 4, however, the space variable x is replaced by the variable r = x 1 2 + + x n 2 + x n + 1 2 .
One may check that the governing equations of this BVP with S = f ( U ) V β and Q = g ( U ) V β + 1 are invariant w.r.t. the Lie symmetry operator
X = 2 β t t + β r r + 2 V V ,
which has the same structure as (43).
Theorem 4.
The nonlinear BVP (57) with S = f ( U ) V β and Q = g ( U ) V β + 1 is invariant with respect to the two-dimensional MAI generated by the Lie symmetry operators (58) and t .
Thus, the Lie symmetry (58) generates the ansatz
U = φ ( ω ) , ω = r t , V = t 1 β ψ ( ω ) ,
which immediately specifies the function R ( t ) = ω 0 t and reduces the two-dimensional BVP (57) (with S = f ( U ) V β and Q = g ( U ) V β + 1 ) to the one-dimensional problem
( φ + α * ) m φ + m ( φ + α * ) m 1 φ 2 + n ω ( φ + α * ) m φ + ω 2 φ + f ( φ ) ψ β = 0 , ψ + n ω ψ + g ( φ ) ψ β + 1 = 0 , ω = 0 : φ = ψ = 0 , ω = ω 0 : φ = ψ = 0 , φ = α * 1 m ω 0 2 .
The nonlinear problem (60) is still not integrable for arbitrary given coefficients, however, one with β = 1 is again a solvable case. Direct calculations show that Formulae (47) and (48) can be generalised in this case as follows
φ = α * 1 m 4 ω 0 2 ω 2 , ω < ω 0 , ψ = α * 1 m q 0 16 ( n + 3 ) ω 0 2 ω 2 n + 5 n + 1 ω 0 2 ω 2 , ω < ω 0
and
g = q 0 φ , q 0 > 0 , f = q 0 α * m 1 n + 3 φ φ + α * 1 m ω 0 2 n + 1 × ( m + n + 1 2 ) α * 1 m ( φ + α * ) m m α * 2 m 4 α * m ω 0 2 + 4 ( φ + α * ) m 1 φ + α * 1 m 4 ω 0 2 .
In order to obtain the final formulae for the original concentrations of tumour cells α ( t , x ) and nutrient (drug) c ( t , x ) , we use Formulae (41), ansatz (59) and solution (61). The formulae obtained coincide with those presented in (49) for α ( t , x ) and R ( t ) while for the nutrient (drug) concentration we have the solution
c ( t , r ) = c α * 1 m q 0 16 ( n + 3 ) t ω 0 2 r 2 t n + 5 n + 1 ω 0 2 r 2 t , r < ω 0 t .
The biomedical interpretation of the various regimes coincides with that in Section 4 and is omitted here.

6. Conclusions

The main part of this paper is devoted to the Lie symmetry classification of the class of parabolic-elliptic systems (3). First of all, the structure of form-preserving (admissible) transformations for the systems of the form (3) is established (see Theorem 1) in order to establish possible relations between systems that admit equivalent MAIs. Theorem 2 gives a complete list of inequivalent parabolic-elliptic systems consisting of 35 systems, which are invariant under nontrivial MAIs, i.e., under the three- and higher-dimensional Lie algebras. Moreover, we have established all possible subclasses of parabolic-elliptic systems (see Table 3 and Table 4) with nontrivial Lie symmetry, which are reducible to those listed in Table 1 and Table 2 by the relevant form-preserving transformations. Notably, our result is a new confirmation of importance of finding form-preserving transformations because they allow reducing essentially (in contrast to the group of equivalence transformations) the number of inequivalent systems. In fact, the traditional Lie–Ovsiannikov algorithm, which is based on the group of equivalence transformations, would lead to 57 parabolic-elliptic systems instead of 35 systems listed in Table 1 and Table 2 because 22 systems from Table 3 and Table 4 should be added.
From the applicability point of view, the most interesting parabolic-elliptic systems occur in Table 1 because the systems from Table 2 are semi-coupled. It is worth noting that MAIs for the systems arising in all the cases of Table 1 excepting 7 and 14 are some analogs of those for single RD equations. In fact, fixing the diffusion coefficient in the first equation of the corresponding system one may easily identify its generic analog in Table 4 [24] (see cases 3–5 and 7 therein). For example, MAIs of the systems from cases 8 and 9 contain the well-known algebra s l ( 2 , R ) as a subalgebra (see the second, third and fourth operators in the last column). The same subalgebra occurs for the diffusion equation with the same diffusivity (see case 3 in Table 4 [24]). By the way, Ovsiannikov was the first to identify this special case for nonlinear diffusion equations [17]. Of course, the systems in Table 1 have more general structures (many of them contain arbitrary functions as coefficients) in contrast to their analogs listed in Table 4 [24]. However, it happens on a regular basis if one considers systems instead of single equations (see, e.g., case 10 in Table 1 [13] involving the same diffusivity).
Cases 7 and 14 in Table 1 have no analogs among single RD equations. Moreover, the relevant MAIs generate infinite-dimensional Lie algebras because they contain the operators X involving an arbitrary smooth function φ ( t ) (see the last column of Table 1). Interestingly, both MAIs from Cases 7 and 14 again contain the well-known algebra s l ( 2 , R ) as a subalgebra. In fact, if one sets φ ( t ) = t and φ ( t ) = t 2 then the operators t , X with φ ( t ) = t and X with φ ( t ) = t 2 produce nothing else but s l ( 2 , R ) with new representations (it is a simple task to check the commutator relations), which do not happen in the case of single RD equations.
In order to demonstrate applicability of our pure theoretical result, we apply the Lie symmetry corresponding to one-parameter group of scale transformation for solving a (1+1)-dimensional BVP with a free boundary modeling tumour growth [16]. We reduce the given nonlinear BVP to that governed by ODEs. As a result, its exact solution was constructed under additional restrictions on the coefficients arising in the governing equations. Moreover, possible biological interpretations of the solution are discussed. In particular, it was shown that the results obtained allow plausible interpretation in the case when the model describes growth of tumour cells consuming a nutrient. Finally, the results for the above (1+1)-dimensional BVP were generalised on multi-dimensional case under the assumption that tumour cells create a radially symmetrical shape. In a forthcoming paper, we will generalise these results using the definition of conditional invariance for BVPs introduced in [33].

Author Contributions

All the authors contributed equally to the work. All authors read and approved the final manuscript.

Acknowledgments

This research was supported by a Marie Curie International Incoming Fellowship to the first author within the 7th European Community Framework Programme (project BVPsymmetry 912563).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ames, W.F. Nonlinear Partial Differential Equations in Engineering; Academic: New York, NY, USA, 1972; ISBN 9780080955247. [Google Scholar]
  2. Murray, J.D. Mathematical Biology; Springer: Berlin, Germany, 1989. [Google Scholar]
  3. Murray, J.D. Mathematical Biology, II: Spatial Models and Biomedical Applications; Springer: Berlin, Germany, 2003. [Google Scholar]
  4. Okubo, A.; Levin, S.A. Diffusion and Ecological Problems. Modern Perspectives, 2nd ed.; Springer: Berlin, Germany, 2001. [Google Scholar]
  5. Lotka, A.J. Undamped oscillations derived from the law of mass action. J. Am. Chem. Soc. 1920, 42, 1595–1599. [Google Scholar] [CrossRef]
  6. Volterra, V. Variazionie fluttuazioni del numero d’individui in specie animali conviventi. Mem. Acad. Lincei 1926, 2, 31–113. (In Italian) [Google Scholar]
  7. Turing, A.M. The chemical basis of morphogenesis. Philos. Trans. R. Soc. B 1952, 237, 37–72. [Google Scholar] [CrossRef]
  8. Zulehner, W.; Ames, W.F. Group analysis of a semilinear vector diffusion equation. Nonlinear Anal. 1983, 7, 945–969. [Google Scholar] [CrossRef]
  9. Cherniha, R.; King, J.R. Lie symmetries of nonlinear multidimensional reaction-diffusion systems: I. J. Phys. A Math. Gen. 2000, 33, 267–282. [Google Scholar] [CrossRef]
  10. Cherniha, R.; King, J.R. Addendum: Lie symmetries of nonlinear multidimensional reaction-dfiffusion systems: I. J. Phys. A Math. Gen. 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 Math. Gen. 2003, 36, 405–425. [Google Scholar] [CrossRef]
  12. Knyazeva, I.V.; Popov, M.D. A system of two diffusion equations. In CRC Handbook of Lie Group Analysis of Differential Equations; CRC Press: Boca Raton, FL, USA, 1994; Volume 1, pp. 171–176. [Google Scholar]
  13. Cherniha, R.; King, J.R. Non-linear reaction-diffusion systems with variable diffusivities: Lie symmetries, ansatze and exact solutions. J. Math. Anal. Appl. 2005, 308, 11–35. [Google Scholar] [CrossRef]
  14. Cherniha, R.; King, J.R. Lie symmetries and conservation laws of non-linear multidimensional reaction–diffusion systems with variable diffusivities. IMA J. Appl. Math. 2006, 71, 391–408. [Google Scholar] [CrossRef]
  15. Torrisi, M.; Tracina, R. An application of equivalence transformations to reaction diffusion equations. Symmetry 2015, 7, 1929–1944. [Google Scholar] [CrossRef]
  16. Byrne, H.; King, J.R.; McElwain, D.L.S.; Preziosi, L. A two-phase model of solid tumour growth. Appl. Math. Lett. 2003, 16, 567–573. [Google Scholar] [CrossRef]
  17. Ovsiannikov, L.V. Group relations of the equation of non-linear heat conductivity. Dokl. Akad. Nauk SSSR 1959, 125, 492–495. (In Russian) [Google Scholar]
  18. Kingston, J.G. On point transformations of evolution equations. J. Phys. A Math. Gen. 1991, 24, L769–L774. [Google Scholar] [CrossRef]
  19. Kingston, J.G.; Sophocleous, C. On form-preserving point transformations of partial differential equations. J. Phys. A Math. Gen. 1998, 31, 1597–1619. [Google Scholar] [CrossRef]
  20. Niederer, U. Schrödinger invariant generalised heat equation. Helv. Phys. Acta 1978, 51, 220–239. [Google Scholar]
  21. Gazeau, J.P.; Winternitz, P. Symmetries of variable coefficient Korteweg-de Vries equations. J. Math. Phys. 1992, 33, 4087–4102. [Google Scholar] [CrossRef]
  22. Ames, W.F.; Anderson, R.L.; Dorodnitsyn, V.A.; Ferapontov, E.V.; Gazizov, R.K.; Ibragimov, N.H.; Svirshchevskiy, S.R. CRC Handbook of Lie Group Analysis of Differential Equations, vol. 1. Symmetries, Exact Solutions And Conservation Laws; CRC Press: Boca Raton, FL, USA, 1994; pp. 133–136. [Google Scholar]
  23. Cherniha, R.; Serov, M. Symmetries, ansätze and exact solutions of nonlinear second-order evolution equations with convection term. Eur. J. Appl. Math. 1998, 9, 527–542. [Google Scholar] [CrossRef]
  24. Cherniha, R.; Serov, M.; Rassokha, I. Lie symmetries and form-preserving transformations of reaction-diffusion-convection equations. J. Math. Anal. Appl. 2008, 342, 1363–1379. [Google Scholar] [CrossRef]
  25. Cherniha, R.; Serov, M.; Pliukhin, O. Nonlinear Reaction-Diffusion-Convection Equations: Lie and Conditional Symmetry, Exact Solutions and Their Applications; CRC Press Taylor and Francis Group: Boka Raton, FL, USA, 2018; ISBN 9781498776172. [Google Scholar]
  26. Cherniha, R.; Davydovych, V. Reaction-diffusion systems with constant diffusivities: Conditional symmetries and form-preserving transformations. In Algebra, Geometry and Mathematical Physics; Springer: Berlin/Heidelberg, Germany, 2014; Volume 85, pp. 533–553. [Google Scholar]
  27. Arrigo, D.J. Symmetry Analysis of Differential Equations: An Introduction; John Wiley and Sons: Hoboken, NJ, USA, 2015; ISBN 9781118721445. [Google Scholar]
  28. Bluman, G.W.; Anco, S.C. Symmetry and Integration Methods for Differential Equations; Springer: New York, NY, USA, 2002; ISBN 9780387216492. [Google Scholar]
  29. Bluman, G.W.; Kumei, S. Symmetries and Differential Equations; Springer: Berlin, Germany, 1989. [Google Scholar]
  30. Fushchych, W.I.; Shtelen, W.M.; Serov, M.I. Symmetry Analysis and Exact Solutions of Equations of Nonlinear Mathematical Physics; Kluwer: Dordrecht, The Netherlands, 1993. [Google Scholar]
  31. Olver, P.J. Applications of Lie Groups to Differential Equations; Springer: Berlin, Germany, 1986; ISBN 978-1-4684-0274-2. [Google Scholar]
  32. Cherniha, R.; Kovalenko, S. Lie symmetries and reductions of multi-dimensional boundary value problems of the Stefan type. J. Phys. A Math. Theor. 2011, 44, 485202. [Google Scholar] [CrossRef]
  33. Cherniha, R.; King, J.R. Lie and Conditional Symmetries of a Class of Nonlinear (1 + 2)—Dimensional Boundary Value Problems. Symmetry 2015, 7, 1410–1435. [Google Scholar] [CrossRef]
Figure 1. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter values m = 1 , α * = 0.5 , c = 2 , ω 0 = 1 , q 0 = 0.5 .
Figure 1. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter values m = 1 , α * = 0.5 , c = 2 , ω 0 = 1 , q 0 = 0.5 .
Symmetry 10 00171 g001
Figure 2. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter value m = 0.5 (other parameters as in Figure 1).
Figure 2. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter value m = 0.5 (other parameters as in Figure 1).
Symmetry 10 00171 g002
Figure 3. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter value m = 0 (other parameters as in Figure 1).
Figure 3. Surfaces representing the concentrations α (left) and c (right) of the form (49) for the parameter value m = 0 (other parameters as in Figure 1).
Symmetry 10 00171 g003
Figure 4. Surfaces representing the function S ( α ( t , x ) , c ( t , x ) ) (left) and Q ( α ( t , x ) , c ( t , x ) ) (right) on the known concentrations given by (49) (other parameters as in Figure 3).
Figure 4. Surfaces representing the function S ( α ( t , x ) , c ( t , x ) ) (left) and Q ( α ( t , x ) , c ( t , x ) ) (right) on the known concentrations given by (49) (other parameters as in Figure 3).
Symmetry 10 00171 g004
Figure 5. Curves representing the functions f ( φ ) (48) for parameters m = 0 (green), m = 0.5 (red) and m = 1 (blue). The other parameters are as in Figure 1.
Figure 5. Curves representing the functions f ( φ ) (48) for parameters m = 0 (green), m = 0.5 (red) and m = 1 (blue). The other parameters are as in Figure 1.
Symmetry 10 00171 g005
Table 1. Lie symmetries of system (4) in the case F V G U 0 .
Table 1. Lie symmetries of system (4) in the case F V G U 0 .
CaseRD SystemBasic Operators of MAI
1. U t = D ( U ) U x x + e V f ( U ) t , x , 2 t t + x x 2 V
0 = V x x + e V g ( U )
2. U t = D ( U ) U x x + V β f ( U ) t , x , 2 β t t + β x x 2 V V
0 = V x x + V β + 1 g ( U )
3. U t = U k U x x + U γ + 1 f e V U α t , x ,
0 = V x x + U γ k g e V U α 2 γ t t + ( γ k ) x x 2 U U + 2 α V
4. U t = U k U x x + U γ + 1 f V U β t , x ,
0 = V x x + V U γ k g V U β 2 γ t t + ( γ k ) x x 2 U U + 2 β V V
5. U t = U k U x x + α 1 e V U γ + 1 t , x , 2 t t + x x 2 V ,
0 = V x x + α 2 e V U γ k k t t U U + ( γ k ) V
6. U t = U k U x x + α 1 V β U γ + 1 t , x , 2 β t t + β x x 2 V V ,
0 = V x x + α 2 V β + 1 U γ k k β t t β U U + ( γ k ) V V
7. U t = U k U x x + α 1 U k + 1 + U V t , x , k φ ( t ) t + φ ˙ ( t ) U U +
0 = V x x + α 2 U k k φ ˙ ( t ) V + ( k + 1 ) φ ¨ ( t ) V
8. U t = U 4 / 3 U x x + U f ( U V 3 ) t , x , 2 x x 3 U U + V V ,
0 = V x x + U g ( U V 3 ) x 2 x 3 x U U + x V V
9. U t = U 4 / 3 U x x + α 1 U 1 + γ V 3 γ t , x , 2 x x 3 U U + V V ,
0 = V x x + α 2 U 4 / 3 + γ V 3 γ + 1 x 2 x 3 x U U + x V V ,
4 γ t t + 3 γ U U ( 4 / 3 + γ ) V V
10. U t = e U U x x + e ( γ + 1 ) U f V + α U t , x ,
0 = V x x + e γ U g V + α U 2 ( γ + 1 ) t t + γ x x 2 U + 2 α V
11. U t = e U U x x + e ( γ + 1 ) U f V e β U t , x ,
0 = V x x + V e γ U g V e β U 2 ( γ + 1 ) t t + γ x x 2 U + 2 β V V
12. U t = e U U x x + α 1 exp ( γ + 1 ) U + V t , x , 2 t t + x x 2 V ,
0 = V x x + α 2 exp γ U + V t t U + γ V
13. U t = e U U x x + α 1 e ( γ + 1 ) U V β t , x , 2 β t t + β x x 2 V V ,
0 = V x x + α 2 e γ U V β + 1 β t t β U + γ V V
14. U t = e U U x x + α 1 e U + V x , φ ( t ) t φ ˙ ( t ) U
0 = V x x + α 2 e U φ ˙ ( t ) V + φ ¨ ( t ) V
Table 2. Lie symmetries of system (4) in the case F V G U = 0 and F V 2 + G U 2 0 .
Table 2. Lie symmetries of system (4) in the case F V G U = 0 and F V 2 + G U 2 0 .
CaseRD SystemBasic Operators of MAI
1. U t = D ( U ) U x x t , x , 2 t t + x x + 2 V V ,
0 = V x x + g ( U ) φ ( t ) x V , ψ ( t ) V
2. U t = D ( U ) U x x + f ( U ) t , x , ψ ( t ) V V
0 = V x x + V g ( U )
3. U t = D ( U ) U x x + f ( U ) t , x , h ( t , x ) V
0 = V x x + α 2 V + g ( U )
4. U t = U k U x x + α 0 U β + 1 t , x , 2 β t t + ( β k ) x x 2 U U ,
0 = V x x + α V U β k , α 0 ψ ( t ) V V
5. U t = U k U x x + U V x , k x x + 2 U U , k φ ( t ) t +
0 = V x x φ ˙ ( t ) U U + ( k φ ˙ ( t ) V + ( k + 1 ) φ ¨ ( t ) ) V
6. U t = U k U x x + α 1 U k + 1 t , x , k t t + U U + h 0 ( x ) V , h ( t , x ) V ,
0 = V x x + ln U + α 2 V
7. U t = U k U x x + α 1 U k + 1 t , x , k t t + U U + β V V ,
0 = V x x + U β + α 2 V h ( t , x ) V
8. U t = U k U x x + α U γ + 1 , α γ 0 t , x , 2 γ t t + ( γ k ) x x 2 U U +
0 = V x x + ln U 2 ( γ k ) V + x 2 V , φ ( t ) x V , ψ ( t ) V
9. U t = U k U x x + α U γ + 1 , α γ 0 t , x , 2 γ t t + ( γ k ) x x 2 U U +
0 = V x x + U β 2 ( γ β k ) V V , φ ( t ) x V , ψ ( t ) V
10. U t = U k U x x t , x , 2 k t t + 2 U U x 2 V ,
0 = V x x + ln U 2 t t + x x + 2 V V , φ ( t ) x V , ψ ( t ) V
11. U t = U k U x x t , x , k t t + U U + β V V ,
0 = V x x + U β 2 t t + x x + 2 V V , φ ( t ) x V , ψ ( t ) V
12. U t = U 4 / 3 U x x + α U , α 0 t , x , 2 x x 3 U U + V V , φ ( t ) x V
0 = V x x + U x 2 x 3 x U U + x V V , ψ ( t ) V
13. U t = U 4 / 3 U x x t , x , 2 x x 3 U U + V V ,
0 = V x x + U x 2 x 3 x U U + x V V ,
4 t t + 3 U U + V V , φ ( t ) x V , ψ ( t ) V
14. U t = e U U x x + α 0 e ( β + 1 ) U t , x , 2 ( β + 1 ) t t + β x x 2 U ,
0 = V x x + α V e β U , α 0 ψ ( t ) V V
15. U t = e U U x x + V t , x , x x + 2 U ,
0 = V x x φ ( t ) t φ ˙ ( t ) U φ ˙ ( t ) V + φ ¨ ( t ) V
16. U t = e U U x x + α 1 e U t , x , t t + U + h 0 ( x ) V , h ( t , x ) V
0 = V x x + U + α 2 V
17. U t = e U U x x + α 1 e U t , x , t t + U + β V V , h ( t , x ) V
0 = V x x + e β U + α 2 V
18. U t = e U U x x + α e γ U , α γ 0 t , x , 2 γ t t + ( γ 1 ) x x 2 U +
0 = V x x + U 2 ( γ 1 ) V + x 2 V , φ ( t ) x V , ψ ( t ) V
19. U t = e U U x x + α e γ U , α γ 0 t , x , 2 γ t t + ( γ 1 ) x x 2 U +
0 = V x x + e β U 2 ( γ β 1 ) V V , φ ( t ) x V , ψ ( t ) V
20. U t = e U U x x t , x , 2 t t 2 U + x 2 V ,
0 = V x x + U 2 t t + x x + 2 V V , φ ( t ) x V , ψ ( t ) V
21. U t = e U U x x t , x , t t U β V V ,
0 = V x x + e β U 2 t t + x x + 2 V V , φ ( t ) x V , ψ ( t ) V
Table 3. RD systems that are reduced to those in Table 1 by the form-preserving transformations.
Table 3. RD systems that are reduced to those in Table 1 by the form-preserving transformations.
RD SystemTransformation of VariablesCase of Table 1
1. u τ = u k u y y + u k + 1 f e v u α + λ u U = e λ τ u , V = v + α λ τ , 3
0 = v y y + g e v u α t = 1 k λ e k λ τ with γ = k
2. u τ = u k u y y + u k + 1 f v u β + λ u U = e λ τ u , V = e β λ τ v , 4
0 = v y y + v g v u α t = 1 k λ e k λ τ with γ = k
3. u τ = u 4 / 3 u y y + u f u v 3 3 u 1 / 3 U = u cos 3 y , V = v cos 1 y , 8
0 = v y y + u g u v 3 + v x = tan y
4. u τ = u 4 / 3 u y y + u f u v 3 + 3 u 1 / 3 U = e 3 y u , V = e y v , 8
0 = v y y + u g u v 3 v x = 1 2 e 2 y
5. u τ = u 4 / 3 u y y + α 1 u 1 + γ v 3 γ 3 u 1 / 3 U = u cos 3 y , V = v cos 1 y , 9
0 = v y y + α 2 u 4 / 3 + γ v 3 γ + 1 + v x = tan y
6. u τ = u 4 / 3 u y y + α 1 u 1 + γ v 3 γ + 3 u 1 / 3 U = e 3 y u , V = e y v , 9
0 = v y y + α 2 u 4 / 3 + γ v 3 γ + 1 v x = 1 2 e 2 y
7. u τ = e u u y y + e u f v + α u + λ U = u λ τ , V = v + α λ τ , 10
0 = v y y + g v + α u t = 1 λ e λ τ with γ = 0
8. u τ = e u u y y + e u f v e β u + λ U = u λ τ , V = e β λ τ v , 11
0 = v y y + v g v e β u t = 1 λ e λ τ with γ = 0
Table 4. RD systems that are reduced to those in Table 2 by the form-preserving transformations.
Table 4. RD systems that are reduced to those in Table 2 by the form-preserving transformations.
RD SystemTransformation of VariablesCase of Table 2
1. u τ = u k u y y + α 1 u k + 1 + λ u U = e λ τ u , t = 1 k λ e k λ τ , 6
0 = v y y + ln u + α 2 v V = v + λ α 2 τ , i f α 2 0 , v + λ 2 τ y 2 , i f α 2 = 0
2. u τ = u k u y y + α 1 u k + 1 + λ u U = e λ τ u , t = 1 k λ e k λ τ , 7
0 = v y y + u β + α 2 v + γ V = v + γ α 2 e β λ τ , i f α 2 0 , v + γ 2 y 2 e β λ τ , i f α 2 = 0
3. u τ = u k u y y + α u γ + 1 V = v + λ 2 y 2 9, if α 0
0 = v y y + u β + λ 11, if α = 0
4. u τ = u k u y y + λ u U = e λ τ u , t = 1 k λ e k λ τ , 10
0 = v y y + ln u V = v + λ 2 τ y 2
5. u τ = u k u y y + λ u U = e λ τ u , t = 1 k λ e k λ τ , 11
0 = v y y + u β + γ V = v + γ 2 y 2 e β λ τ
6. u τ = u 4 / 3 u y y + α u 3 u 1 / 3 U = u cos 3 y , V = v cos 1 y , 12
0 = v y y + u + v x = tan y
7. u τ = u 4 / 3 u y y + α u + 3 u 1 / 3 U = e 3 y u , V = e y v , 12
0 = v y y + u v x = 1 2 e 2 y
8. u τ = u 4 / 3 u y y 3 u 1 / 3 U = u cos 3 y , V = v cos 1 y , 13
0 = v y y + u + v x = tan y
9. u τ = u 4 / 3 u y y + 3 u 1 / 3 U = e 3 y u , V = e y v , 13
0 = v y y + u v x = 1 2 e 2 y
10. u τ = e u u y y + α 1 e u + λ U = u λ τ , t = 1 λ e λ τ , 16
0 = v y y + u + α 2 v V = v + λ α 2 τ , i f α 2 0 , v + λ 2 τ y 2 , i f α 2 = 0
11. u τ = e u u y y + α 1 e u + λ U = u λ τ , t = 1 λ e λ τ , 17
0 = v y y + e β u + α 2 v + γ V = v + γ α 2 e β λ τ , i f α 2 0 , v + γ 2 y 2 e β λ τ , i f α 2 = 0
12. u τ = e u u y y + α e γ u V = v + λ 2 y 2 19, if α 0
0 = v y y + e β u + λ 21, if α = 0
13. u τ = e u u y y + λ U = u λ τ , t = 1 λ e λ τ , 20
0 = v y y + u V = v + λ 2 τ y 2
14. u τ = e u u y y + λ U = u λ τ , t = 1 λ e λ τ , 21
0 = v y y + e β u + γ V = v + γ 2 y 2 e β λ τ

Share and Cite

MDPI and ACS Style

Cherniha, R.; Davydovych, V.; King, J.R. Lie Symmetries of Nonlinear Parabolic-Elliptic Systems and Their Application to a Tumour Growth Model. Symmetry 2018, 10, 171. https://doi.org/10.3390/sym10050171

AMA Style

Cherniha R, Davydovych V, King JR. Lie Symmetries of Nonlinear Parabolic-Elliptic Systems and Their Application to a Tumour Growth Model. Symmetry. 2018; 10(5):171. https://doi.org/10.3390/sym10050171

Chicago/Turabian Style

Cherniha, Roman, Vasyl’ Davydovych, and John R. King. 2018. "Lie Symmetries of Nonlinear Parabolic-Elliptic Systems and Their Application to a Tumour Growth Model" Symmetry 10, no. 5: 171. https://doi.org/10.3390/sym10050171

APA Style

Cherniha, R., Davydovych, V., & King, J. R. (2018). Lie Symmetries of Nonlinear Parabolic-Elliptic Systems and Their Application to a Tumour Growth Model. Symmetry, 10(5), 171. https://doi.org/10.3390/sym10050171

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