The reaction-diffusion master equation (RDME) is a stochastic model for spatially heterogeneous chemical systems. Stochastic models have proved to be useful for problems from molecular biology since copy numbers of participating chemical species often are small, which gives a stochastic behaviour. The RDME is a discrete space model, in contrast to spatially continuous models based on Brownian motion. In this thesis two accuracy issues of the RDME on unstructured meshes are studied. The first concerns the rates of diffusion events. Errors due to previously used rates are evaluated, and a second order accurate finite volume method, not previously used in this context, is implemented. The new discretisation improves the accuracy considerably, but unfortunately it puts constraints on the mesh, limiting its current usability. The second issue concerns the rates of bimolecular reactions. Using the macroscopic reaction coefficients these rates become too low when the spatial resolution is high. Recently, two methods to overcome this problem by calculating mesoscopic reaction rates for Cartesian meshes have been proposed. The methods are compared and evaluated, and are found to work remarkably well. Their possible extension to unstructured meshes is discussed.
| Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:uu-145978 |
| Date | January 2011 |
| Creators | Kieri, Emil |
| Publisher | Uppsala universitet, Avdelningen för teknisk databehandling |
| Source Sets | DiVA Archive at Upsalla University |
| Language | English |
| Detected Language | English |
| Type | Student thesis, info:eu-repo/semantics/bachelorThesis, text |
| Format | application/pdf |
| Rights | info:eu-repo/semantics/openAccess |
| Relation | UPTEC F, 1401-5757 ; 11014 |
Page generated in 0.002 seconds