| The objectives in this paper is to present a general numerical solution for temperature field in a fully developed unidirectional flow of porous medium in a duct of any cross-section area. The computation is divided in two phases. The first is to determine the velocity field using the Brinkman's Momentum Equation. The second is to determine the temperature field using the Energy Equation.; In the first phase the Galerkin Based Integral Method is used to determine velocity field. This method allows to use a polynomial function form, basis function, which is finite continuous and single valued in a duct cross-section shape as the representative form of the velocity field.; Depending on the duct cross-section shape, formation of the basis function, which represents the velocity field, may be simple or very tedious. For the purpose of simplicity, flow through a parallel plate channel and an isosceles triangle are presented in this paper. The same algorithm can be used for any duct shape.; The governing equations for computation of temperature field are the Brinkman's Momentum equation and the energy equation. For the theoretical model solutions, the algorithm presented in section 10 and 11 of Ref. 2 is used to solve the energy equation. For the mathematica solutions, the program designed by Dr. Haji, presented in Appendix A and B, is used. |