Simulating the contamination of populations by a polluted environment

We derive two reaction-diffusion systems modeling the impact of a nuclear waste contamination on food chain population sets. We exhibit the numerical difficulties the problem features and provide second order accurate nonnegative approximate solutions in two space dimensions.