A Boundary Domain Integral Method is used for the solution of natural conve
ction in a rectangular cavity, with vertical walls maintained at different
constant temperatures and the top and bottom walls being adiabatic. Cavity
is filled with an isotropic, homogeneous, rigid, fluid-saturated porous med
ia assuming that the solid and the sigle-phase fluid are in thermal equilib
rium The Brinkman-extended Darcy momentum equation with inertial (transport
) term, Brinkman (viscous) term and Darcy law, completed with continuity an
d heat energy equation,is considered. The respective system of non-linear p
artial differential equations is written in the Velocity-Vorticity Formulat
ion.