This presentation will be given by the team of researchers involved in progress on mathematical and computational models of biofilm at porescale including Azhar Alhammali (Imam Abdulrahman Bin Faisal University) and Oregon State University researchers Lisa Bigler, Kyla Jones, Choah Shin, and Naren Vohra, and Malgo Peszynska. The presentation will report on some of the results in two recently published papers. Our focus is on a new unconstrained close-to-singular diffusion model for biofilm phase which features an adaptively chosen singular diffusivity. Connections to and comparison with previously studied constrained model will be given, as well as to the plethora of osmotic pressure, phase field and indvividual based models known from the literature. The model is coupled to nutrient utilization; we will also discuss reduced models such as stationary approximations with comparisons to the literature. For the flow we will discuss our newly developed local heterogeneous Brinkman model as well as the impact of the choice of the flow and biofilm model on upscaled Darcy permeability. Our entire coupled model is monolithic and computationally robust even in complex pore-scale geometries, and extends to multiple species.