Difference between revisions of "Team:KU Leuven/Modeling/Top"

Starting from the biological circuit described on the Research page . Two cell types A and B are interacting. Type A cells produce a repellent called leucine which causes the cells of type B to move away. At the same time type A cells also produce OHHL, which is required by the cells of type B to move. Initially, colonies of the two cell types are placed at the center of the dish. As molecule production within the type A cells kicks in, the repellent and OHHL concentrations start to increase. This triggers the type B cells to move away from the center. Movement will continue until the concentration of OHHL is insufficient for the type B cells to move further. The behaviour of the two cell types is described by the model given below:

The model has been derived while looking at ^[1] and ^[2] . The terms that appear can be grouped into four categories. Every equation has a diffusion term given by $D_x \bigtriangledown^2 X$, diffusion smoothes peaks by spreading them out in space. The two equations related to cell densities contain logistic growth terms of the form $\gamma X(1-\frac{X}{k_x})$, which model the cell growth during simulation time. Finally the second equation describing the moving cells comes with a variable coefficient Poisson term $\bigtriangledown (P \bigtriangledown X)$ which describes the cell movement. Last but not least we have the two bottom equations. These two model concentrations. Both contain linear production and degradation terms, which look like $kX$. It is important to keep in mind that even though the degradation terms appear as linear terms in the differential equation the solution will be exponential decay.
To generate the video file above the system above has been discretized using a finite volume approach in conjunction, with an explicit Euler scheme, For finite volume methods to work we rewrite our equations as conservation laws. Then each grid point is assigned the area around it such that flux of cells or molecules leaving one cell enters another one. From discretizing the integrated conservation law the following expression is obtained in one dimension:

The for the equations given above the left hand side values at the next time step depend exclusively on data of the prefious time step as illustrated in the figure below:

Figure 1

computational molecule. Click to enlarge

The image above shows the dependency of data in time and space. The computational molecule used in this case uses only data of the previous time level $t_n$ to compute data at the next time level $t_{n+1}$. A scheme with a space time dependency like the one shown above is called an explicit scheme.

The video box above shows the solution of the discretized system in one dimension. To gain additional insight into the effect of the different terms of the model, we computed simulations of different term combinations. Use the buttons to choose from the videos. The first term in each equation is a diffusion term. Diffusion smooths out edges of an initial condition, eventually it leads to an even distribution. The initial condition as rectangular the illustrate the smoothing. An other important part of the two first equations which model bacteria density is the logistic growth term. The video which visualizes logistic growth starts with a Gaussian distributed initial condition, which is more realistic then the rectangular initial condition used in the diffusion term simulation.
Zero flux and periodic boundary conditions have been implemented. The boundaries are the edges of the domain on which the equation system is solved. Here the domain ranges from zero to eight centimetres, which is the diameter of a Petri dish. With zero flux boundaries the first derivative is set to zero at the boundaries, which means that neither bacteria nor chemicals are allowed to pass trough the boundary. Periodic boundaries connect pairs of boundaries to each other, which means that cells leaving at the top of the boundary appear at the bottom, cells leaving at the left boundary reappear at the right and so on. In the continuous context these boundary conditions have been implemented to allow comparisons with the hybrid model, where these boundaries are also used. Finally simulation has been done using the parameters given in the table below:

Using the equation system as described above, the model may also be simulated in two dimensions. Once more a finite volume approach has been taken in connection with an explicit Euler scheme. All parameters have been kept constant with the one exception of the chemotactic sensitivity $K_c$. Which has been increased to $Kc = 1.5 * 10^{-1}$$cm^2/h$, which leads to earlier pattern formation. Above four simulation videos with Gaussian initial conditions can be observed. A fifth video shows a simulation using random initial data. The two last videos illustrate the effect of zero flux and periodic boundary conditions.