The use of lead reverberatory furnaces is a common part of the lead recycling process, wherein a solid feed material comprising various lead compounds is directly exposed to combustive and radiative heat sources to melt the feed material and produce lead product. Furnace design decisions ranging from burner inputs to burner placement and furnace shape are difficult to trial without risking costly losses in productivity. The use of modern computational fluid dynamics approaches can allow for more-accurate modeling of the reverberatory furnace process. A steady state simulation of the reverberatory furnace has been created for this purpose. The direct modeling of the melting and smelting reactions is bypassed as this effort would be very computationally expensive and require detailed kinetic data, a large portion of which is unknown. However, as these reactions are largely endothermic, their influence on the overall temperature profile of the furnace cannot be ignored and is critical to both validation and usefulness of the model. To this end, a multi-variable artificial heat sink boundary condition has been placed on several key boundary surfaces within the furnace for the purpose of representing the heat consumption caused by the melting and smelting reactions. This boundary condition is tuned through parameter adjustment in the base study until the temperature profile matches sample data taken from furnace operation measurements.