We consider a relatively simple model for pool boiling processes. This model involves only the temperature distribution within the heater and describes the heat exchange with the boiling fluid via a nonlinear boundary condition imposed on the fluid-heater interface. This results in a standard heat equation with a nonlinear Neumann boundary condition on part of the boundary. In this paper we analyse the qualitative structure of steady-state solutions of this heat equation. It turns out that the model allows both multiple homogeneous and multiple heterogeneous solutions in certain regimes of the parameter space. The latter solutions originate from bifurcations on a certain branch of homogeneous solutions. We present a bifurcation analysis that reveals the multiple-solution structure in this mathematical model. In the numerical analysis a continuation algorithm is combined with the method of separation-of-variables and a Fourier collocation technique. For both the continuous and discrete problem a fundamental symmetry property is derived that implies multiplicity of heterogeneous solutions. Numerical simulations of this model problem predict phenomena that are consistent with laboratory observations for pool boiling processes.