Properties of Solutions of Kolmogorov-Fisher Type Biological Population Task with Variable Density ()
Received 5 April 2016; accepted 20 May 2016; published 23 May 2016
1. Introduction
Population models are studied for a long time. The first such work was done by Gause G.F. and Fisher R.D., and mathematical studies were performed by Kolmogorov, Petrovskii (KPP) and Piskunov (1937) in the famous paper [1] - [4] . They were interested in the behavior of the speed of the wave solutions and the resulting estimate of the speed of wave propagation.
Then there were other models of the population [5] - [8] . In recent years, intensive study of nonlinear models was based on diffusion and revealed new properties of finite speed of propagation of diffusion waves (see [3] and the literature given there). We have proposed a population model of two competing populations with non-linear double diffusion and variable density that are described by nonlinear system of competing individuals. We identify new properties, such as finite speed of propagaton, and localization of the outbreaks in a specific area. In particular, in the critical case, the rate type CPT generalizes their result.
Statement of the Task
In this paper, we investigate the properties of solutions of biological population task of Fisher-Kolmogorov type in the case of variable density. The main research method is a self-similar approach. Considering in the field, there is a parabolic system of two quasilinear equations of reaction-diffusion
(1)
, , (2)
which describes the process of biological population of Kolmogorov-Fisher in a nonlinear two-component environment, and mutual diffusion coefficients which are respectively equal to, . Numeric parameters, , are positive real numbers, and, , ,;, is desired solutions.
We study properties of solutions to problem (1), (2) based on the self-similar analysis of solutions of a system of equations constructed by the method of nonlinear splitting and a reference equations and bringing the system (1) for radially symmetric mind. Note that replacing in (1)
,
leads to the form
(3)
,. (4)
If, choosing
, , ,
we get the following system of equations:
(5)
where,
,
If, and, , system has the form:
(6)
A significant role in the study of the Cauchy problem and boundary problems for Equations (1) has self- similar solutions. Under self-similar solution we will understand as particular solutions of Equation (1), depending on the combination of t and x. Knowledge of them plays a sometimes crucial role in the study of various properties of solutions of the original equations.
Below we describe one way of obtaining self-similar system for the system of Equations (5). It consists in the following. We find first the solution of a system of ordinary differential equations
in the form
, ,
for the case of, and,. And in the case of, and, we solve a system of ordinary differential equations
in the form
, ,
then the solution of system (5) is sought in the form
(7)
and is selected so
if.
Then for we get the system of equations:
(8)
where
(9)
If, self-similar solution of system (9) has the form
, (10)
Then substituting (10) into (8) with respect to gets a self-similar system of equations
(11)
where и.
System (11) has an approximate solution of the form
, ,
where А and В are constants and
,.
In this paper, on the basis of the aforesaid methods, we studied qualitative properties of solutions of the system (1), solved the problem of choosing the initial approximation for iterative, leading to fast convergence to the solution of the Cauchy problem (1), (2), depending on the values of numerical parameters and initial data. For this purpose, as the initial approximation was used, we found the asymptotic representation of the solution. This has allowed to perform numerical experiments and visualization of the process described by system (1), depending on the values included in the system of numeric parameters.
2. Construction of Upper Solutions
Let us build an upper solution for system (11).
Note that the functions have properties
and
We choose A and B from the system of nonlinear algebraic equations
Then functions were the solution of the Zel’dovich-Kompanees for the system (1) and in the field they satisfy the system of equations
in the classical sense.
Due to the fact that
function and the flows have the following smoothness properties
.
We choose A and B such that the inequality of inequality
(12)
Since then
It is due to the fact that
from (12) we have
Then in the field Q according to the comparison principle of solutions have
Theorem 1. Let Then for the solution of problem (1) Q is the evaluation
where―above-defined functions.
Note that the solution of system (1) when has the following representation at
where B(a, b)-Euler Beta function.
It is proved that this view is the self-similar asymptotics of solutions of systems (1).
Here
3. Slow Diffusion
Case (slow diffusion). Applying the method of [1] to solve Equation (11) will receive the following features
where, ,. It is known [1] [2] that for the global existence of solutions of problem (1) function must satisfy the following inequality:
and
.
Let’s take the function, and show that they are asymptotic finite solutions of the system (11).
Theorem 2. The finite solution of system (11) when has asymptotic.
Proof. We seek a solution of Equation (8) in the following form
, (13)
where, and at, to explore the asymptotic stability of the solution of problem (11) when. Substituting (13) into (11) for gets the following equation
(14)
where above-defined function.
Note that the study of the solution of the last equation is equivalent to examining the solution of Equation (11), each of which in a certain period satisfies the inequality:
,.
Let us show first of all that decision Equations (14) have a finite limit at. We introduce the notation
.
Then for the Equation (14) has the form
.
To analyze the last expression we introduce a new helper function
,
where―real numbers. Hence it is easy to see that each value function stores the sign on a certain interval and for all either of the inequalities
,.
And so for the function there is a limit when. From the expression for it follows that
Hence, given that
, , ,
get the following algebraic equation
The latter system gives and because (14).
Theorem 2 is proved.
4. Fast Diffusion
Case (fast diffusion case). For (11), we have
where.
Theorem 3. At vanishing at infinity solution of problem (11) has the asymptotics .
Proof. In the proof of theorem used the transform
, (15)
where, which leads (11) to the following form.
Substituting (15) into (11) for gets the following equation
(16)
where is above-defined function.
Note that the study of the solution of the last equation is equivalent to examining the solution of Equation (11), each of which in a certain period satisfies the inequality:
,.
Let us show first of all that solution of the Equation (16) has a finite limit at. We in-
troduce the notation. Then Equation (15) has the form
.
To analyze the last expression we introduce a new helper function
,
where -real number. Hence it is easy to see that each value function stores the sign on a certain interval and at all satisfied either of the inequalities
,.
Therefore for the function there is a limit when. From the expression for follows that
.
Hence, when, we get the following algebraic equation
The calculation of the last equation gives and because (15).
Theorem 3 is proved.
5. Computational Experiment
Investigation of qualitative properties of system (1) has allowed to perform numerical experiment depending on the values included in the system of numeric parameters. For this purpose, the initial approximation was used to construct asymptotic solutions. The numerical solution of the problem for the linearization of system (2) was used linearization methods of Newton and Picard. To build self-similar system of equations of biological population used the method of nonlinear splitting [1] [6] .
For the numerical solution of the problem (1) we will construct a uniform grid
and temporal grid
.
Replace the problem (1) implicit difference scheme and receive differential task with the error.
It is known that the main problem for the numerical solution of nonlinear problems is the appropriate choice of the initial approximation and the method of linearization of system (1).
Consider the function:
where иabove-defined functions,
Record means. These functions have the property of finite speed of propagation of perturbations [1] [6] . Therefore, for the numerical solution of problem (1) when as an initial approximation of the proposed function,.
Created on input language Matlab the program allows you to visually trace the evolution process for different values of the parameters and data.
Numerical calculations show that in the case of arbitrary values qualitative properties of solutions do not change. Below are the results of numerical experiments for various values of the parameters (Figures 1-4).
6. Conclusions
Thus, the proposed nonlinear mathematical model of biological populations with double nonlinearity and variable density properly describes the studied process. Numerical study of nonlinear processes described by equations with a double nonlinearity and analysis results on the basis of evaluation solutions provides a comprehensive picture of the process in two-component systems competing biological population with the preservation of localization properties in the target area and the size of the flash.
Results in future will provide an opportunity to evaluate the speed of propagation of diffusive waves.