Team:Paris Saclay/Modeling

From 2014.igem.org

(Difference between revisions)
(Modeling)
(Modeling)
Line 11: Line 11:
-
{| class="wikitable centre" width="50%"
+
{| class="wikitable centre" width="100%"
|+   
|+   
|-
|-

Revision as of 16:53, 30 September 2014


Modeling

Work in progress (Don't worry, i'll explain all the mathematic formula :D)


To realise our artwork, we use an agarose gel to obtain the shape of a lemon. To push the resemblance to the extreme, we wish to have a crust in the edge of the lemon when we seperate it. In fact, we build bacteria who produce yellow/green color in presence of oxygen. Thus we must evaluate the penetration of the oxygen in the gel !

We use the following phenomenological law suggest by Adolphe Fick in 1855:


In an homogeneous and isotropic environment, containing particles distributed inhomogeneously,appears spontaneously a volumetric flow density vector particle $\overrightarrow{J}(M,t) $. In any point $M$ in space, this vector is proportional to the gradient of the particle density $n(M,t)$. Mathematicaly, this relationship take the form: \[ \overrightarrow{J}(M,t) = - D \times \nabla n(M,t) \qquad (1) \] where $D$ is the diffusion coefficient.


Referring to the article [1], the diffusion coefficient of oxygen in agarose is $ D = 0{,}256 \times 10^{-8} m^2 s^{-1} $.

  • To simplify the problem, we consider that the diffusion of oxygen particle occurs only in one direction. So $\overrightarrow{J}(M,t) = J(x,t) \overrightarrow{e}_x $.
  • Spatial variations in the density of particles are connected to spatial variations of the vector $\overrightarrow{J}(M,t)$ by the material's equation of conservation in presence of volume distribution of particle source $\sigma (x,t)$ (device which injects or subtracted particles to the system) :

\[ \frac{\partial n}{\partial t} (x,t) = - \frac{\partial J}{\partial x} (x,t) + \sigma (x,t) \qquad (2) \]

By replacing $(2)$ in $(1)$, we obtain the following equation of diffusion \[ \forall t, \forall x, \bigg( \frac{\partial}{\partial t} - D \frac{\partial^2}{\partial x^2} \bigg) n(x,t) = \sigma (x,t) \qquad (3) .\]

As our lemon is exposed to the ambient air, we stay in steady state where the source $ \sigma (x,t) $ is equal to $N_0$ the quantity of $O_2$ in the air.

To solve this equation, we use Fourier's analysis (+ d'explications)

\[ \forall x, \forall t>0, n(x,t) = \frac{N_0}{\sqrt{4 \pi D t}} exp \bigg(- \frac{x^2}{4 D t} \bigg) + \int_{0}^{t} \underbrace{N_0 * exp \bigg(- \frac{|x|^2}{4 D \tau} \bigg)}_{= 0 \text{ by symmetry of the gaussian distribution }} \frac{d\tau}{\sqrt{4 \pi D \tau} } \]

The average dispersion particle is given by the variance $\Delta x = \sqrt{2Dt}$. Using this formula, we deduct that oxygen will penetrate $3 \times 10^{-3} m$ in $1956.522 s = 32.6082 $ minutes.

Graphical visualization

Paris Saclay oxygenGraph.jpeg

References:

[1] A.C. Hulst, H.J.H. Hens, R.M. Buitelaar and J. Tramper, Determination of the effective diffusion coefficient of oxygen in gel materials in relation to gel concentration, Biotechnology Techniques Vol 3 No 3 199-204 (1989).

[2] Vincent Renvoizé, Physique PC-PC*, Cap Prepas, Pearson Education, 2010.