This paper describes an approach to furnace charge calculation for melting processes. Some complexities of the problem, including non-homogeneous element loss, non-metal contaminations, and, melt correction have been considered in the model. The contribution of the model is that it examines a rearrangement of the non-linear complexities into an iterative standard linear programming framework. The performance of the model was evaluated and demonstrated on an industrial scale test problem. The results revealed that the non-linear phenomenological complexities can be effectively fitted in the framework of an iterative standard LP model. This feature provides a backbone for reliable and fast optimization in melting which is of significant benefit for industrial automation.