Faculty Advisor

Larsen, Christopher J.

Faculty Advisor

Vernescu, Bogdan M.


This project models the heat flow of a gas when it is forced through a porous medium in one spatial dimension. A Forchheimer conservation of momentum equation for non-Darcy flow is solved simultaneously with the conservation of mass equation using the standard Galerkin finite element method with linear elements. The conservation of energy {convection-diffusion) equation is solved with a streamline diffusion method to reduce the oscillations induced by the standard Galerkin finite element method. A Picard iteration scheme is used to handle the non-linear terms in the model and to iterate between the energy equation and the conservation of mass and momentum system. A backward Euler method is used to handle the time discretization.


Worcester Polytechnic Institute

Date Accepted

January 2001



Project Type

Major Qualifying Project


Restricted-WPI community only

Advisor Department

Mathematical Sciences