1. RRWPT - Particle tracking simulations for reactive transport

This is a code for performing reactive transport simulations in one-dimensional and two-dimensional systems. The reaction is described by A+B->P, with some rate constant k. The simulation is done by means of particle tracking and reduction of the mass of particles by the reaction probability. The results manifest incomplete mixing, also known as the Ovchinnikov–Zeldovich segregation phenomena. That is, the domain becomes segregated into islands where one of the species is abundant and the other is missing. As a result, the reaction slows down, in comparison with the well-mixed rate. When this happens, the concentration tends to scale like t^(-d/4) instead of scaling like t^-1 in the well-mixed solution (d is the system dimension, d=1,2,3).

Click here for further discussion and sample results

Contributors

This code was written by Amir Paster, in collaboration with Diogo Bolster (ND) and David Benson (CSM). 

It was made possible with the financial support of NSF through grant EAR-1113704 (2011-2014), Mixing and Reactions in Highly Heterogeneous Media - Can Non-local Models handle them?

Code and Documentation

QUICK START GUIDE

ZIP FOLDER WITH REQUIRED MATLAB FILES