A two-dimensional finite element model of the biofilm response to flow was developed. The numerical code sequentially coupled the fluid dynamics of turbulent, incompressible flow with the mechanical response of a single hemispherical biofilm cluster (similar to 100 mu m) attached to the flow boundary. A non-linear Burger material law was used to represent the viscoelasic response of a representative micro-bial biofilm. This constitutive law was incorporated into the numerical model as a Prony series representation of the biofilm's relation and modulus. Model simulations illuminated interesting details of this fluid-structure interaction. Simulations revealed that softer biofilms (characterized by lower elastic moduli) were highly suspectible to lift forces and consequently were subject to even greater drag forces and consequently were subject to even greater drag forces found higher in the velocity field. A bimodal deformation path due to the two Burger relaxation times was also observed in several simulations. This suggested that inter-facial biofilm may be most susceptible to hydrodynamically induced detachment during the initial relaxation time. This result may prove useful in developing removal strategies. Additionally, plots of lift versus drag suggested that the deformation paths taken by viselastic biofilms are largely insensitive to specific material coefficients. Softer biofilms merely seem to follow the same path (as a stiffer biofilm) at a faster rate. These relationships may be useful in estimating the hydrodynamic forces acting on an attached biofilm based on changes in scale and cataloged material properties.