The use of computer simulations in biology and the medical research field has gained in popularity. These simulations are providing researchers the opportunity to better predict the behavior of biological systems before performing long and expensive physical trials. The modeling of large biological systems would benefit from a method of approximating fluid flow quickly and accurately. Currently, no analytical solution exists; instead, many different numerical methods attempt to provide accurate approximations. They are referred to as Computational Fluid Dynamic solvers (CFD). The Discrete Event System Specification (DEVS) has rarely been used for modeling the physics of fluid flow. In this thesis we show how Cell-DEVS, a derivative of the DEVS formalism that conforms to the Cellular Automata parameters, can be used to provide realistic approximations of fluid flow. The algorithms used in the CFD presented in this thesis are based on the Navier-Stokes equations for non-linear fluid flow, which are an extension to Newton's Second Law of motion. The goal of the Cell-DEVS based CFD model will be to accurately approximate the fluid flow with minimal computational effort. Furthermore, the design of the solver should be such that it can be easily adjusted for use in a wide range of biological systems.iii