The objective of this project is to investigate the behaviour of the (neutron) diffusion equation in two dimensions and implement various boundary conditions such as Dirichlet, Robin/mixed and Neumann boundary conditions. An in-depth derivation of the diffusion equation and its boundary conditions can be found in the notebook. The geometry of the simulation is a rough two-dimensional approximation of the Demon core incident from the Manhattan project during WWII.
This project was created using Jupyter Notebooks, a scientific python IDE which is availablie for all operating systems (Mac, Windows and Linux distros). For Mac and Windows users, I highly recommend downloading Jupyter Notebooks through the Anaconda Navigator. For Linux users, I recommend installing it with pip or pip3 . Linux download link