1. A method of computational fluid dynamic modeling of production from a subterranean region comprising a wellbore and surrounding formation, the method comprising: providing a computational geometric model of the subterranean region, the computational geometric model comprising the wellbore and surrounding formation; determining one or more zones in the surrounding formation of the computational geometric model, wherein the one or more determined zones surround the wellbore; associating a viscous resistance with at least one or more of the determined zones as a computational parameter, wherein the viscous resistance associated with one or more of the determined zones is a viscous resistance calculated from a sample material obtained from the surrounding formation, wherein the calculated viscous resistance is determined by altering the sample material obtained from the surrounding formation and exposing the sample material to one or more expected conditions associated with the subterranean region in order to provide formation damage; and computationally modeling, using a computer, the production through the one or more determined zones into the wellbore using the associated viscous resistance of each of the one or more determined zones.