Team:ETH Zurich/Modeling/Microfluidics
From 2011.igem.org
(→Simulation) |
(→Overview) |
||
Line 16: | Line 16: | ||
| | | | ||
== Overview == | == Overview == | ||
- | In our inital model, we chose a 5 cm long cylindrical channel with a radius of 1 mm (or equivalently 2 mm diameter). This channel is connected to a reservoir filled with medium and acetaldehyde. For the purpose of modeling the degradation of acetaldehyde in the channel, we assumed that the reservoir has constant acetaldehyde concentration (or equivalently has infinite size). | + | In our inital model, we chose a 5 cm long cylindrical channel with a radius of 1 mm (or equivalently 2 mm diameter). This channel is connected to a reservoir filled with medium and acetaldehyde. For the purpose of modeling the diffusion and degradation of acetaldehyde in the channel, we assumed that the reservoir has constant acetaldehyde concentration (or equivalently has infinite size). |
[[File:ETHZ-MicrofludicsModel.png|center]] | [[File:ETHZ-MicrofludicsModel.png|center]] | ||
Revision as of 15:29, 16 September 2011
Microfluidics Model |
| ||||
In our first spatiotemporal model, we wanted to find out if an Acetaldehyde gradient would be formed at all. The idea was to generate it as an equilibrium of Acetaldehyde diffusion out of the reservoir into our channel and Acetaldehyde degradation by the cells in it. In order to answer this question, we modeled a spatiotemporal reaction-diffusion system, first in 2D with MATLAB and then in 3D with COMSOL Multiphysics. |
OverviewIn our inital model, we chose a 5 cm long cylindrical channel with a radius of 1 mm (or equivalently 2 mm diameter). This channel is connected to a reservoir filled with medium and acetaldehyde. For the purpose of modeling the diffusion and degradation of acetaldehyde in the channel, we assumed that the reservoir has constant acetaldehyde concentration (or equivalently has infinite size). |
DiffusionFirst of all, we have to model diffusion of Acetaldehyde from the reservoir to the closed end of the microfluidics channel. The relevant partial differential equation is the one for isotropic diffusion, i.e. diffusion that is uniform in all directions: For the relevant diffusion parameter, DAcAl, see the parameters page. |
DegradationCell DensityIn order to able to determine how much Acetaldehyde our smoColi cells are able to degrade, we first need to determine the cell density of the immobilized cells in agarose. For this calculation, we assume we have a cell culture with OD600 with which we fill the channel. Parameter-wise, to determine the relative cell density in the channel voume, we need to know:
A parameter we can choose suitably is the dilution of the cells in agarose, cDilution. Then we can calculate the relative cell density in the smoColi-agarose-filled microfluidics channel: Enzymatic Degradation |
SimulationSteady StateWe can see that the diameter does not influence the mean dynamics of the acetaldehyde gradient. We note that there may be stochastic behavior for very thin channels, since there will be only very few cells. For the channel diameter we chose a mean field approximation should be sufficient though as the amount of cells involved in the degradation is sufficiently high. Dynamics |