A Stokesian fluid in motion along a porous medium saturated by the same fluid is modelled by the Beavers—Joseph—Saffman boundary-value problem to generalized Forchheimer—Stokes—Fourier systems: what we call the Beavers—Joseph—Saffman (BJS) problem. The model has nonlinear character given by the temperature dependence of physical parameters such as the viscosity, the permeability, the thermal conductivity and the thermal expansion. The paper is concerned with the study of the steady-state and the time-dependent regimes via the Galerkin and the Faedo—Galerkin techniques, respectively.