This thesis covers the development of a model for fluid flow which incorporates computational fluid dynamics simulations using three-dimensional planar porous media networks. Porous media are introduced along with applications and the need for computational models is discussed. Previous experiments and models are presented as well as features of the current model. This model constructs three-dimensional planar networks from cylindrical pipes and elbows of varying length, diameter, and angle. Simulations are carried out using a finite volume based computational fluid dynamics software. A methodology is provided to discuss the source code of the executable created to automate the modeling process. This process begins with the network creation from an existing code, which generates sets of random pore networks called 'realizations' and ends with linear and polynomial regressions used to provide curve fits for Darcy's law and Forchheimer's equation. Findings of these parameters are presented for varying porosity values of Berea sandstone simulated with single phase liquid water. Results show that the model follows Forchheimer's equation for certain porosities and follows experimental results. Finally, remarks on future work are given and a closing summary is presented.