Team:TU-Delft/Modeling/Diffusion

From 2012.igem.org

(Difference between revisions)
 
(153 intermediate revisions not shown)
Line 1: Line 1:
-
{{:Team:TUDelft/Scripts}}{{:Team:TUDelft/CSSLaksh}}{{:Team:TUDelft/menu}}
+
{{:Team:TUDelft/Scripts}}{{:Team:TUDelft/CSS}}{{:Team:TUDelft/menu}}
<html>
<html>
<head><title>Diffusion Model</title></head>
<head><title>Diffusion Model</title></head>
Line 6: Line 6:
<img src="https://static.igem.org/mediawiki/igem.org/7/74/Background_greenlamp.jpg" class="bg_team"></img>
<img src="https://static.igem.org/mediawiki/igem.org/7/74/Background_greenlamp.jpg" class="bg_team"></img>
<div id="logo_ed"><a href="https://2012.igem.org/Team:TU-Delft" 'onfocus=this.blur()'><img src="https://static.igem.org/mediawiki/igem.org/8/88/Logoigemklein.png" border="0" width="100" height="100"></a></div>
<div id="logo_ed"><a href="https://2012.igem.org/Team:TU-Delft" 'onfocus=this.blur()'><img src="https://static.igem.org/mediawiki/igem.org/8/88/Logoigemklein.png" border="0" width="100" height="100"></a></div>
 +
<div id="contentbox" style="text-align:justify;">
-
<img src="https://static.igem.org/mediawiki/igem.org/8/8f/Diffusion_Header.jpg" align="middle" width="100%">
+
<img src="https://static.igem.org/mediawiki/igem.org/8/81/Diffusionmodel.jpg" align="middle" width="100%">
-
</body></html>
+
 
 +
</div>
 +
</body>
 +
</html>
<div id="contentbox" style="text-align:justify;">
<div id="contentbox" style="text-align:justify;">
-
One of the main objectives of the project was to synthesize a practical device, the '''Snifferometer''' for tuberculosis detection. As a first step towards achieving this goal, we built a temporal model of the system using '''PDE's''' which was simulated in matlab. A 2D reaction-diffusion system was then implemented in COMSOL multiphysics using the knowledge obtained from single cell pathway model,ombining the behaviours of the which helped us get a better understanding of how such a device could be implemented and the response times involved in such a process.  
+
One of the main objectives of the project was to synthesize a practical device, the '''Snifferometer''' for tuberculosis detection. As a first step towards this, we built a spatio temporal 2D reaction-diffusion system which was implemented in COMSOL multiphysics. Combining the behaviours of the diffusion and the pathway model, helped us get a better understanding of how such a device could be implemented and the response times involved in such a process.  
-
= Initial Model =
+
= The device set up =
 +
The device at the bottom has a membrance preventing the yeast coming in direct contact with [[File:SnifferometerCAD.png|200px|left|thumb|'''Figure 5''': Computer aided design of the Snifferometer]]external substances, above which is a layer of agar in which yeast is placed,the nutrients for it's growth is fed in through a channel along the sides of the tube using which the yeast cells can also be replaced as and when needed. A micro-optrode is then used to sense the photons emitted by the fluoroscent proteins. A key question was the time it took for the ligand to diffuse into the layer of yeast. To answer this question, the System was modeled as a 1D diffusion problem using PDE’s.
-
[[File:Diffusion Model1.png|320px|middle|thumb|'''Figure 1''': Full structure of diffusion device.]]
+
The values for various dimensions of the model were chosen using an iterative approach of design and analysis which would yield an optimal response. The model dimensions can be obtained from the code in the codes section. Once the structure of the device was designed, a simple reaction model assuming linear kinetics was implemented in the agar layer using the results of the single cell pathway model([https://2012.igem.org/Team:TU-Delft/Modeling/SingleCellModel#Results '''Results'''])
 +
</div>
 +
<div id="contentbox" style="text-align:justify;">
 +
 
 +
= Modeling Approach =
 +
== PDEs ==
 +
The diffusion phenomenon is modelled by '''Fick's second law'''.  
 +
https://static.igem.org/mediawiki/igem.org/8/83/DiffusionFormula1.PNG
 +
 
 +
where C represents the ratio of the pressure between gasphase and atmosphere, x is the distance from the surface of petridish to the agar layer, and D is the diffusion coefficient, the value of which is calculated from [1].
 +
 
 +
In order to solve PDEs, numerical methods are used as approximation.
 +
 
 +
On the left side of the equation, the Euler forward method is taken:
 +
 
 +
https://static.igem.org/mediawiki/igem.org/1/1a/DiffusionFormula2.PNG
 +
 
 +
On the right side of the equation, the central differential method is used:
 +
 
 +
https://static.igem.org/mediawiki/igem.org/b/b6/DiffusionFormula3.PNG
 +
 
 +
It can be noticed that the equations both to the right and left of the equality are the same.
 +
 
 +
https://static.igem.org/mediawiki/igem.org/3/35/DiffusionFormula4.PNG
 +
 
 +
Therefore, the PDEs can be numerically replaced by the equation:
-
For the initial experiments a model of the diffusion in the Petridish could be useful to evaluate the effectiveness of the experimental setup.  
+
https://static.igem.org/mediawiki/igem.org/4/4c/DiffusionFormula5.PNG
-
Modeling of this setup was especially done to evaluate the speed of diffusion from the ligand phase to the yeast cells themselves and see if this behavior is slow or quick with respect to the biochemical behavior of the yeast.  
+
where https://static.igem.org/mediawiki/igem.org/b/b8/DiffusionFormula6.PNG. In the implementation by MATLAB, dx is taken as the interval of x.
-
Partial Differential Equations (PDE) are used to solve this problem.
 
== Boundary conditions ==
== Boundary conditions ==
-
Two boundary conditions are set for these two surfaces of petridish-gas and gas-agar.
+
In order to solve this equation, boundary conditions are required. Two boundary conditions are set for these two surfaces of petridish-air and air-agar.
*At x=0, the Dirichlet boundary condition was placed:  
*At x=0, the Dirichlet boundary condition was placed:  
-
https://static.igem.org/mediawiki/igem.org/b/b8/DirichletBC.PNG
+
https://static.igem.org/mediawiki/igem.org/archive/b/b8/20120926220011%21DirichletBC.PNG
-
with Psat being the saturation pressure of the specific compound being IsoAmylAcetate as an example in our model. For IsoAmylAcetate this value is 533.3 Pa at room temperature.  
+
with Csat being the pressure ratio with respect to the saturation pressure of the specific compound being IsoAmylAcetate as an example in our model. For IsoAmylAcetate this value is 533.3 Pa at room temperature. Thus Csat is equal to 0.05333.
-
*One can take Neumann boundary condition at x=l, where l is the distant between petridish and agar and assumed to be 0.01 meter.
+
*One can take Neumann boundary condition at x=l, where l is the distant between petridish and agar and assumed to be 0.01 meter in the model. The reason we can use Neumann boundary is because the molecules bounce back into the gasphase against the agar surface.
https://static.igem.org/mediawiki/igem.org/e/ea/NeumannBC.PNG
https://static.igem.org/mediawiki/igem.org/e/ea/NeumannBC.PNG
-
In this part of the system the molecules bounce against the agar surface back into the gasphase.
+
The model was the simulated in matlab, the results of which are below
-
A van Neumann boundary is the way to give a mathematical workable structure to that.
+
-
For the problem definition to be complete and be able to start programming in MATLAB
+
-
we only need one more value, being the diffusion coefficient. For this an empirical relation from the [1] was used.
+
-
== PDE ==
+
</div>
-
Now we have a complete problem that we will solve using a numerical approach.
+
<div id="contentbox" style="text-align:justify;">
-
Put this in the pressure domain.
+
= Results =
 +
[[File:output of diffusionNew.png|320px|middle|thumb|'''Figure 2''': Simulation of diffusion model.(in 3D Cartesian coordinate, X axis is time, Y axis is the distance from the petridish, Z axis is the pressure ratio)]]Figures 2 and 3 depict the results of the simulation with distance dimension 0 - 0.01 meter and the time dimension 0 - 100 seconds. Two key observations were drawn from this model
 +
*The results suggest that at approximately two minutes, there is steady state concentration at the agar layer
 +
*The input to the single cell pathway model can be considered a step function, if we assume that the Ligand is not degraded by yeast.
 +
-
And you probably guessed but my value for h is ok in the model and in my code this translates to the line:
+
*Matlab script for diffusion model can be found [https://static.igem.org/mediawiki/igem.org/c/c2/IGEM_Model_Pertidish.txt'''here'''].
 +
</div>
-
dcdt(i,1)=  a*x(i-1,1)+a*x(i+1,1) -2*a*x(i,1);
 
-
and YES a=h here and YES that IS confusing considering the other above statements.
+
<div id="contentbox" style="text-align:justify;">
-
What I basically did is grab a bunch of Ordinary differential equations and couples them
+
-
So that it becomes a PDE ( partial differential equation ).  H*time step should also not be bigger than 0.5 and that does not happen, so the solution is stable.
+
-
Analytical solutions for this problem appeared a little more troublesome than expected.
+
= Comsol Spatio-Temporal Model =
-
There are solutions ( in Sum row form ), but I was anticipating a nice function.
+
A spatio temporal 2D reaction-diffusion system which was implemented in COMSOL multiphysics. Combining the behaviours of the diffusion and the pathway model.The simulation of which is shown below.
-
I will get a book at the library to see if there are any better ones, but for now this works out fine.
+
-
I can also use a Cranck-Nicolson approach to the problem to have another numerical solution for verification. And last but not least. The results:
+
-
 
+
-
[[File:output of diffusion.png|320px|middle|thumb|'''Figure 2''': Simulation of diffusion model.(in 3D Cartesian coordinate, X axis is time, Y axis is the distance from the petridish, Z axis is the concentration)]]
+
-
 
+
-
Figure 2 shows the simulation result with distance dimension 0 - 0.01 meter and the time dimension 0 - 100 seconds. From Figure2 it can be seen after 100 seconds it approximately reaches steady state at the x=0.01 boundary.
+
 +
<html>
 +
<center>
 +
<iframe width="560" height="315" src="http://www.youtube.com/embed/qoObgkd4tsY" frameborder="0" allowfullscreen></iframe>
 +
<h6>In case the place seems to be blank, zoom into the page by using Ctrl + Mouse-wheel and hit F5.</h6>
 +
</center>
 +
</html>
</div>
</div>
 +
<div id="contentbox" style="text-align:justify;">
<div id="contentbox" style="text-align:justify;">
-
= Snifferometer =
+
= Conclusion =
 +
From the simulations and analysis it can be concluded that the response of the system can be improved by improving the kinetics of the signalling pathway and the diffusion has a limited influence on the system response.
-
One of the other goals of the diffusion modeling was to model the device which we intended to build for the project. We made use of the finite element analysis simulator Comsol Multiphysics[[Team:TU-Delft/Modeling/Diffusion#Ref3|[3]]]
+
</div>
-
for developing this model.[[File:Snifferometer.png|300px|Right|thumb|'''Figure 1''': Snifferometer - Device with the modified yeast cells for sensing]]
+
-
Different models were analyzed for their suitability before coming to the final design in Figure. The device at the bottom has a membrance preventing the yeast coming in direct contact with external substances, above which is a layer of agar in which yeast is placed, the nutrients for it's growth is fed in through a channel along the sides of the tube using which the yeast cells can also be replaced as and when needed. A micro-optrode is then used to sense the photons emitted by the fluoroscent proteins.
+
<div id="contentbox" style="text-align:justify;">
 +
= Comsol Code =
 +
=== Final Model ===
 +
[[File:TUD-Download.png|50px|link=http://thenewview.eu/igem/feromone.mph|left]]
 +
<br clear="all" />
</div>
</div>
 +
<div id="contentbox" style="text-align:justify;">
<div id="contentbox" style="text-align:justify;">
-
== Saliva Model ==
 
-
[[File:Test1.swf|300px|Right|thumb|'''Figure 1''': Snifferometer - Device with the modified yeast cells for sensing]]
 
= References =  
= References =  
-
Chen, N. H. (1962). New Generalized Equation for Gas Diffusion Coefficient. J. Chem. Eng. Data, 37–41.
+
<h6>
-
 
+
1. Chen, N. H. (1962). New Generalized Equation for Gas Diffusion Coefficient. J. Chem. Eng. Data, 37–41.
 +
</h6>
 +
<h6>
 +
2. Comsol Multiphysics.
 +
</h6>
</div>
</div>
<html><img src="https://static.igem.org/mediawiki/igem.org/3/37/Footer_2.jpg" align="middle" width="690">
<html><img src="https://static.igem.org/mediawiki/igem.org/3/37/Footer_2.jpg" align="middle" width="690">
<a href='https://2012.igem.org/Main_Page' target="_blank"><div id='logo_igem2'></div><a/></body></html>
<a href='https://2012.igem.org/Main_Page' target="_blank"><div id='logo_igem2'></div><a/></body></html>

Latest revision as of 01:46, 27 October 2012

Menu

Diffusion Model

One of the main objectives of the project was to synthesize a practical device, the Snifferometer for tuberculosis detection. As a first step towards this, we built a spatio temporal 2D reaction-diffusion system which was implemented in COMSOL multiphysics. Combining the behaviours of the diffusion and the pathway model, helped us get a better understanding of how such a device could be implemented and the response times involved in such a process.

Contents

The device set up

The device at the bottom has a membrance preventing the yeast coming in direct contact with
Figure 5: Computer aided design of the Snifferometer
external substances, above which is a layer of agar in which yeast is placed,the nutrients for it's growth is fed in through a channel along the sides of the tube using which the yeast cells can also be replaced as and when needed. A micro-optrode is then used to sense the photons emitted by the fluoroscent proteins. A key question was the time it took for the ligand to diffuse into the layer of yeast. To answer this question, the System was modeled as a 1D diffusion problem using PDE’s.

The values for various dimensions of the model were chosen using an iterative approach of design and analysis which would yield an optimal response. The model dimensions can be obtained from the code in the codes section. Once the structure of the device was designed, a simple reaction model assuming linear kinetics was implemented in the agar layer using the results of the single cell pathway model(Results)

Modeling Approach

PDEs

The diffusion phenomenon is modelled by Fick's second law. DiffusionFormula1.PNG

where C represents the ratio of the pressure between gasphase and atmosphere, x is the distance from the surface of petridish to the agar layer, and D is the diffusion coefficient, the value of which is calculated from [1].

In order to solve PDEs, numerical methods are used as approximation.

On the left side of the equation, the Euler forward method is taken:

DiffusionFormula2.PNG

On the right side of the equation, the central differential method is used:

DiffusionFormula3.PNG

It can be noticed that the equations both to the right and left of the equality are the same.

DiffusionFormula4.PNG

Therefore, the PDEs can be numerically replaced by the equation:

DiffusionFormula5.PNG

where DiffusionFormula6.PNG. In the implementation by MATLAB, dx is taken as the interval of x.

Boundary conditions

In order to solve this equation, boundary conditions are required. Two boundary conditions are set for these two surfaces of petridish-air and air-agar.

  • At x=0, the Dirichlet boundary condition was placed:

20120926220011%21DirichletBC.PNG

with Csat being the pressure ratio with respect to the saturation pressure of the specific compound being IsoAmylAcetate as an example in our model. For IsoAmylAcetate this value is 533.3 Pa at room temperature. Thus Csat is equal to 0.05333.

  • One can take Neumann boundary condition at x=l, where l is the distant between petridish and agar and assumed to be 0.01 meter in the model. The reason we can use Neumann boundary is because the molecules bounce back into the gasphase against the agar surface.

NeumannBC.PNG

The model was the simulated in matlab, the results of which are below

Results

Figure 2: Simulation of diffusion model.(in 3D Cartesian coordinate, X axis is time, Y axis is the distance from the petridish, Z axis is the pressure ratio)
Figures 2 and 3 depict the results of the simulation with distance dimension 0 - 0.01 meter and the time dimension 0 - 100 seconds. Two key observations were drawn from this model
  • The results suggest that at approximately two minutes, there is steady state concentration at the agar layer
  • The input to the single cell pathway model can be considered a step function, if we assume that the Ligand is not degraded by yeast.


  • Matlab script for diffusion model can be found here.


Comsol Spatio-Temporal Model

A spatio temporal 2D reaction-diffusion system which was implemented in COMSOL multiphysics. Combining the behaviours of the diffusion and the pathway model.The simulation of which is shown below.

In case the place seems to be blank, zoom into the page by using Ctrl + Mouse-wheel and hit F5.

Conclusion

From the simulations and analysis it can be concluded that the response of the system can be improved by improving the kinetics of the signalling pathway and the diffusion has a limited influence on the system response.

Comsol Code

Final Model

TUD-Download.png


References

1. Chen, N. H. (1962). New Generalized Equation for Gas Diffusion Coefficient. J. Chem. Eng. Data, 37–41.
2. Comsol Multiphysics.