^{1}

^{1}

^{1}

We have studied the mathematical model of the growth of a micromycetes colony. The results of the numerical experiments have shown that the model given in the work quite adequately reflects the emergence of the ring structures in the micromycetes colonies.

The attempts of mathematical modeling of biological structures and processes, alongside with the search for natural and experimental models to describe and investigate biological morphogenesis, i.e. emergence of elaborate structures in the development of microorganisms colonies, have been undertaken in mathematical biology for years.

In many biological systems, the mechanism of spatial structures formation is connected with the violation of stability of a homogeneous system. Mathematical models of such processes represent systems of the equations in the derivatives of a reaction-diffusion type.

Apart from the diffusive instability, an essential factor that can promote structuring is chemotaxis, i.e. the direct movement of cells along the chemical gradient.

A number of experiments have been described in literature where the formation of spatial structures in bacterial systems is observed, wherein both moving and stationary structures are possible.

For the first time, the research results of diffusion equation related to the growth of the substance amount were applied to the study of biological problems in the works by Kolmogorov A.N., Petrovsky I.G. and Piskunova H.C. [

Moving structures were described by Adler [

The mathematical model of such a wave movement was offered by Keller and Segel in 1971 [

where—is the bacteria density at time t and position (x,y),—is the concentration of the chemoattractant,—is the specific growth rate of bacteria,—is the function that indicates the sensitivity of bacterıa to the chemoattractant,—are the coefficients of bacteria diffusion and a substratum, respectively,—are constants.

Keller Segel’s model and some of her modifications succeeded to describe the formation and distribution of hemotaksis Adler rings and various structuring processes in bacterial colonies during their interaction. It was also shown that bacterial population waves belong to a special class of autowaves.

The study of the processes of the motionless ring structures formation in the work of Polezhaev A.A. et al. [

In a number of works [5-9] there has been offered a mathematical model for the description of ash value emergence in colonies of a micromycetes, where the process of colonies formation is considered as a phenomenon of self-organization caused by the diffusive instability. As the main hypothesis of such structures formation, the assumption of the metabolism products and their diffusion impact on a micromycetes colony was considered. The following system of the differential equations describing the distribution of concentration of a mycelium (), a substratum (S) and metabolism products (m) has been studied as a mathematical model:

Here—is a large-scale coefficient of speed of metabolism products formation,

—speed of consumption of a substratum,

—maximum (resource) concentration of a mycelium,

—specific growth rate of a mycelium concentration,

—threshold concentration of the metabolism products, defining the mycelium growth,

—coefficients of diffusion of metabolism products and a substratum,

—functions of conditions1[x]—Hevisayd’s function (function of single jump),

—relative radial coordinate,

—time,

—delay time of a mycelium response to the presence of metabolism products.

The system was studied together with the boundary conditions of the second kind that correspond to a case of impenetrable external walls.

Authors of the article have carried out experimental observations of the formation of the spatial-temporal structures in the micromycetes colonies, as well as numerical experiments with the use of the given mathematical model. The results of these calculations describe the formation of such types of spatial distribution as wave (zone) structures and a continuous lawn in fungi colonies. The experimental data is well agreed with the mathematical model.

It has been, however, established long before [

We have studied the following mathematical model of the growth of a micromycetes colony:

This system describes the following model experiment. The initial colony, representing a drop with the high or low number of cells, is placed into the Petri’s cup with a certain amount of the substratum. While consuming the substratum, the cells develop metabolism products.

In the fırst equation term describes the growth of metabolism products m, which is proportional to a square concentration of the mycelium. The given dependence is obtained through the experimental data.

The term in the second equation reflects the reduction of the substratum concentration in proportion to mycelium growth taking into account the substratum exhaustion s.

The composed in the third equation is responsible for the formation of the mycelium. The mycelium growth is limited by the following factors:

- the achievement of the limit (resource) concentration- the substratum exhaustion s- the development of metabolism products to the level of the concentration determined by a condition U_{1}.

Together with the system of the equations the following entry conditions are considered:

as well as the boundary conditions of the second kind, meeting the requirements of the external walls impermeability:

To determine the operating parameters the numerical experiments with the use of the offered mathematical model were carried out. As a basic vector of parameters we chose the following:

When fixing the coordinates of a vector the variation of system parameters were made: and.

The numerical modeling of the system showed the following results (

The analysis of the numerical modeling results showed that the most significant factors of the increase in a gradient of the mycelium concentration are the increase λ (

The parameters, in the main only affect the amplitude of the oscillations of the mycelium. In the considered period of time, when the substrate is still little used up, the diffusion of the substrate due to the smallness of the diffusion coefficient also no significant effect on the dynamics of the process.

The process of colonies formation in time has also been tracked.

Upon exhaustion of the resources in the equations system under the research, the steady state will be reached, where the concentration of the mycelium, the substratum and the metabolism ceases to change. The following differential equations have been obtained:

Thus, the results of the numerical experiments have shown that the model given in the work quite adequately reflects the emergence of the ring structures in the micromycetes colonies. By means of variation of the model parameters, we succeeded to get the structures given in the work of Bulyanitsa A.L. et al. [