Team:Grenoble/Modeling/Amplification/Quorum

From 2012.igem.org

Revision as of 15:48, 21 September 2012 by Eliseb (Talk | contribs)

iGEM Grenoble 2012

Project
              
The goal of this part was only to check the speed of the diffusion. We use the same model as the 2011 IGEM team of Grenoble, which is based on the Bangalore 2007 IGEM team model. We first have to add the diffusion terms in the equations thus we get:



Because of the temporal and special derivatives, we couldn’t use a classic matlab solver to solve this set of equations. The approximation we used consisted in dividing the space into a grid:



It is the finite difference method. We thus get:



with h=lx/N, where lx is the length of the grid and N is the number of points of discretization along x. By using the same approximation on yi, and assuming that ly=lx, and that we have the same number of points of discretization, we get:



In addition, we also use the finite difference method on the time scale. Thus, we get:



where p=∆t/Nt where ∆t is the interval of time we work in and Nt the number of points of our time discretization.

With these approximations, we could solve the equations.