The rmogravitational flows of liquid me tals in ope n or close d ducts and containe rs play a re le vant role in a varie ty of applications in mechanical, mate rials and nuclear e ngineering. Such flows are known to be very sensitive to the e ffective shape of the containe r use d to host the fluid and its thermal boundary conditions. For the case of tempe rature gradie nts having the main compone nt dire cte d along the horizontal dire ction, re late d convective phe nome na fall unde r the ge neral heading of "Hadle y flow". He re we introduce a ge neral frame work for the de te rmination of the properties of these flows in the case of domains having converging or diverging top and bottom walls. The frame work is built via a hybrid approach in which typical techniques of CFD are use d in synergy with analytical solutions of the e nergy e quation. The prope r use of initial and boundar y conditions results in algorithm convergence acce le ration. The role playe d by the top and bottom wall inclination with re spect to the horizontal is assessed through parametric inve stigation.