We consider the steady, laminar natural convection heat transfer of a particulate suspension in an electrically-conducting fluid through a two-dimensional channel containing a non-Darcian porous material in the presence of a transverse magnetic field. The transport equations for both fluid and particle phases are formulated using a two-phase continuum model and a heat source term is included which simulates either absorption or generation. A set of transformations are implemented to reduce the partial differential equations for momentum and energy conservation (for both phases) from a two-dimensional coordinate system to a one-dimensional system. Finite element solutions are obtained for the transformed model. A comprehensive parametric study of the effects of the heat source parameter (E), Prandtl number (Pr), Grashof number (Gr), momentum inverse Stokes number (Skm), Darcy number (Da), Forchheimer number (Fs), particle loading parameter (PL), buoyancy parameter (B), Hartmann number (Ha), temperature inverse Stokes number (SkT), viscosity ratio [...], specific heat ratio [...], dimensionless particle-phase wall slip parameter [...] on the dimensionless fluid phase velocity (U), dimensionless particle phase velocity ( ), dimensionless fluid phase temperature [...] and the dimensionless temperature of particle phase [...] are presented graphically. In addition, we also describe numerical solutions for several special cases of the model, for example, the inviscid hydromagnetic two phase non-Darcian free convection, heat transfer [...], forced convection case (GrŽ0) etc. Fluid phase velocities are found to be strongly reduced by the magnetic field, Darcian drag and also Forchheimer drag; a lesser reduction is observed for the particle phase velocity field. The Prandtl number is shown to depress both the fluid temperature and particle phase temperature in the left hand side of the channel but to boost both temperatures at the right hand side of the channel [...]. The inverse momentum Stokes number is seen to reduce fluid phase velocities and increase particle phase velocities. The influence of other thermophysical parameters is discussed in detail and computations compared with previous studies. The model finds applications in MHD plasma accelerators, astrophysical flows, geophysics, geothermics and industrial materials processing.