To reduce the fuel consumption and emissions and also enhance the molten aluminum quality,a mathematical model with user-developed melting model and burning capacity model,were established according to the features of melting process of regenerative aluminum melting furnaces.Based on validating results by heat balance test for an aluminum melting furnace,CFD (computational fluid dynamics) technique,in association with statistical experimental design were used to optimize the melting process of the aluminum melting furnace.Four important factors influencing the melting time,such as horizontal angle between burners,height-to-radius ratio,natural gas mass flow and air preheated temperature,were identified by PLACKETT-BURMAN design.A steepest descent method was undertaken to determine the optimal regions of these factors.Response surface methodology with BOX-BEHNKEN design was adopted to further investigate the mutual interactions between these variables on RSD (relative standard deviation) of aluminum temperature,RSD of furnace temperature and melting time.Multiple-response optimization by desirability function approach was used to determine the optimum melting process parameters.The results indicate that the interaction between the height-to-radius ratio and horizontal angle between burners affects the response variables significantly.The predicted results show that the minimum RSD of aluminum temperature (12.13％),RSD of furnace temperature (18.50％) and melting time (3.9 h) could be obtained trader the optimum conditions of horizontal angle between burners as 64°,height-to-radius ratio as 0.3,natural gas mass flow as 599 m3/h,and air preheated temperature as 639 °C.These predicted values were further verified by validation experiments.The excellent correlation between the predicted and experimental values confirms the validity and practicability of this statistical optimum strategy.