Modeling
Physical Modelling
1. Motivation
One of the main ideas in the project is to control the flow induced by a bacterial carpet in a micro-tube. This is done by changing the tumble frequency of the bacteria when they are exposed to light. It is however not clear what the exact connection between the tumble frequency and the overall liquid flow is. This matter is complicated by the fact that the bacterial pump itself is a selforganizing system and changing the tumbling of the bacterial could therefore change the structure of the overall system completely. One way to assess this connection between tumble frequency and the flow is to create a physical model, that describes the most essential features of the system. Our work connected to the development of such a model is presented in the following. The basis for any model is an understanding of the physical system in question, therefore we will start with a short review of the hydrodynamic properties of flagella and the bacterial pumping system constructed by[artikel reference].
1
2. The real system
So the "real" system that we want to model is a bacterial pump as descriped by []. This is in principle just a small tube 15 $\mu m$ deep, 200 $\mu m$ wide and 15mm long, covered on the inside by a layer of flagellated bacteria. The bacterial layer described by [] is very dense and uniform, with a spacing between each bacterial of less that 1 $\mu m$ and 80$\%$ of the bacteria adhered to the surface as single bacteria. To get a better understanding of the origin of the created flow from this carpet, it is important to understand the structure of the bacterial flagellum.
This consists of 3 major parts, a rotary motor complex, a hook and a filament. The first part creates the rotary motion of the flagellum and the second part serves as a flexible coupling between the touque creating part and the filament. For our pourpose is the filament the most interesting part. This is responsible for the conversion of the rotary motion into a linear thust. The filament is a self-assembling polymeric structure composed of flagellin protein subunits. These are arranged in a circular way to create a hollow helical structure, with a typical width of 120-250Å and a length of 10-15 $\mu m$. A bacteria as E. coli typically has around 10 flagella.[cyber cell] These filaments is able to adopt a wide range of conformation under the induced torque. Numeric studies[2] and empiric results [1] suggests that this conformation is strongly on the hydrodynamic environment that surrounds the flagellum and the direction of rotation. When several flagella rotates counterclockwise the flagella tends to bundle together in a single helix structure, due to the hydrodynamic interactions[2]. <<random biokemi!>> at some point causes the flagella to turn clockwise instead. This induces a sequence of deformations that changes flagella structure and unravels the bundle. This is know as tumble mode.
To be able to model the flow created by a bacterial carpet it is essential to know what kind of flowfield a single flagellum/bundle will create. This has primary been investigated by nurmerical approatch, where the flagella is modeled as semiflexibel hookian systems. Several studies [][] suggests that the flow created from a single flagellum is highly non-uniform, but to some degree circular symmetric at the end of the flagellum.(see figure XX and XX) When the flagella bundle together [] suggest that this symmetry becomes less clear and flow becomes even more complicated. The total movement of the bacteria in tumble mode is zero, so to a first approximation can the field here be thought as non existing.
Now we can return to the All these results refers to flagella moving freely in aquas solution, but this is severly different from the heighly constrained microtube where the bacterial carpet is placed XX.
To summarize we have to model a very dense system of . This is indee not a simple task, and quite a few simplification assumptions have to be made. These will be startingpoint of the next partXX
3. General description
A flagellum creates propulsion by spinning around in a helical shape. Many who models a single flagellum considers this shape rigid and then calculate a flowfield from the spinning helix. Our system consists of many flagella and modelling every single flagellum in this way would take too much time.
An e. coli typically has around 8-10 flagella, but for simplicity we will consider them as one bundle.
The overall result of the spinning flagellum/flagellabundle is that the bacterium moves in an almost straight line. We will therefore consider the forces created by the flagella to be simply a pointforce on the tip of the flagellabundle pointing in the same direction as the bundle.
The size of this force can be approximated by calculating drag on a swimming bacteria. If the bacteria is considered almost spherical the drag force can be calculated by using the formula for stokes flow past a sphere:
<math> F_d = -6 \pi \eta r v </math>
where η is the viscosity of the fluid in which the bacterium is swimming, r is the radius of the bacterium and v is the velocity.
Using the following data
<math> r = 0.4*10^-6 m </math>
<math> v = 50*10^-6 m/s </math>
<math> \eta = 8.94*10^-4 Pa*s </math>
the dragforce and thereby the force created by the flagellabundle of one e. coli is
<math> F_d = 3.37*10^-13 N </math>
The system we are trying to model consists of a lot of bacteria stuck to a wall, but the part we are interested in is really the flagella and the forces they create. So we ignore the bodies of the bacteria, and instead consider the flagella as being stuck directly to the wall, with one end glued to the wall and the other pointing out into the fluid.
Next step is to figure out what kind of flowfield, such a pointforce creates. The flowfield created by a pointforce in a fluid with no walls or other obstructions near it can be calculated using the corresponding Green's function known as the stokeslet. This gives the following flowfield:
File:Indsæt billede af flowfield her
If the pointforce is placed near a wall modifications must be made, since the flowvelocity has to go to zero at the boundary (known as the no-slip condition). This is done by using the Oseen-Blake tensor. Simply described what the Oseen-Blake tensor does, is to create a mirror image of our force on the other side of the wall, thus cancelling the flow near the wall. (This is not a completely accurate description, but rather an intuitive one.)
In our case the flagella are stuck to a wall, so we'll be using the Oseen-Blake tensor. A flowfield corresponding to this is shown below.
indsæt flowfield
In our case the system we are trying to describe is a microtube. This means that the width of the tube is so small, that the forces created by the flagella are not only close to the wall to which the flagella are stuck, but also close to the opposite wall. This presents an interesting problem. Since the Oseen-Blake tensor works by creating a mirrorpoint of the real force on the opposite side of the wall, we will need a mirrorpoint behind the other wall if we are to uphold the no-slip condition. But the mirror forces also affect the flow near the other wall. In order to cancel this effect one could create another mirrorforce, corresponding to each of the mirror forces, but of course these would obstruct each other to, requirering yet more mirror forces. In the end we decided, to see how precise the system would be for one mirrorpoint behind each wall. The flowfield corresponding to this is shown below.
Indsæt flowfield
As the above figures show, there is quite a difference between the two situations. We decided to keep working with both the single-wall and the double-wall flowfields.
The next thing to be considered was whether the flagella are dependent on the flow, ie. if we place a bacterium at an angle θ with the wall will it remain at that angle or will it get pulled around by the flow in the pipe, thus changing its position. The two extreme situations would be to either keep the flagella in a stationary position or to let it be completely dependent on the flow. In order to create an intermediate situation, we could induce each flagellum with a potential, pulling it toward a favored angle. We believe the answer lies somewhere in between the two extremes, but that doesn't mean the extremes can't tell us anything. We therefore decided to create a model in which the flagella keep still, and one where they are affected by the flow and a potential. The size of the potential can always be set to zero if we want to study the flagella without it.
The basic idea now is that every flagella stuck to surface creates its own flowfield. To get the entire flowfield we add together all the flowfields created by the individual flagellum. In the case where the flagella are stationary that is basically it. For flagella that are able to move it's a bit more tricky. This is described in the next chapter.
4. Considerations about velocity
The first thing we had to figure out was how the flow created by all of the other flagella would affect one single flagellum. To do this we decided to approximate a flagellum as a string of spheres and use dragforce calculations to figure out the force with which the flowfields of the other flagella would affect the beads.
This model showed us that in less than 100ns the velocity of the flagella would be the same as the velocity of the fluid when the flagella started with a velocity of zero, after that the two velocity never diverged far from each other. Since the velocity of the flagella always went to the velocity of the fluid on such a short timescale and since these calculations took a lot of computerpower we decided that instead of force calculations we would simpy find the flowvelocity at the tip of the flagellum and convert that directly to the angular velocity of the flagellum.
5. A 2-D model of the system
First, we will make a model consisting of a 1D-grid, afterwards we will expand it to 2D.
1D-model
indsæt billede her
What we do from here is that we neglect the body of the bacteria. Instead we imagine that the flagella are glued to the surface off the pipe in such a way that they
<math>\coprod_{i=1}^N x_i</math>
Litterature
1. hest hest
"
