This year we focused our modeling efforts on the development of our device design. We worked on reaction-diffusion systems which helped us to find the weak points in our designs and to come up with new prototypes.
The model consists of two main parts; diffusion and reaction. The diffusion part predicts how the solution, once applied, will move through the paper and the reaction part models how the different enzymes and chemicals present on the paper will react with each other.
The Base Conditions
The starting conditions of the model are made by using a png image import which defines the different parts of the biosensor i.e.: where the solution is applied, hydrophobic areas, frontiers and where the enzymes are freeze dried.
A diffusion kernel, along with a diffusion matrix, control the movement of the solution on paper and form the backbone of the diffusion model. Our “world” is updated at every time step which, in turn, updates the concentration on every point of the paper. At every time step we also use a list of the reaction rates to update the model with changing concentrations. We use vector calculus throughout the code to make efficient use of resources.