In this work we consider a mathematical and computational model for biofilm growth and nutrient utilization. In particular, we are interested in a model appropriate at a scale of interface. The model is a system of two coupled nonlinear diffusion--reaction partial differential equations (PDEs). One of these PDEs is subject...