Integrated CFD Modeling of Combustion, Heat Transfer, and Oxide Scale Growth in Steel Slab Reheating

In this study, a three-dimensional simulation of a walking-beam reheating furnace was developed to improve the steel slab reheating process and reduce surface oxidation kinetics using computational fluid dynamics (CFD). Combustion, heat transfer, fluid dynamics, and chemical reaction models were integrated into the numerical framework of this study. In addition, dynamic mesh remeshing was coupled through user-defined functions (UDFs), enabling the simultaneous simulation of slab movement and evolution of the involved transport phenomena. Turbulence was modeled with the realizable k-ε formulation, combustion with the Eddy Dissipation model, and radiation with the P-1 model coupled with WSGGM to include CO2 and H2O gas radiation. Scale formation was modeled using customized functions based on Arrhenius-type kinetics and Wagner’s oxidation model, evaluating its growth as a function of time, temperature, and furnace atmosphere. The predicted thermal evolution inside the furnace was validated using industrial data, yielding an average deviation of 5%. Furthermore, the proposed operating conditions led to an average slab temperature of 1289.77 °C at the exit of the homogenization zone, which was 16 °C higher than that under the current operation but still within the target range (1250 ± 50 °C). The reduction in combustion air decreased energy losses and improved product quality, lowering the molar oxygen content in the furnace atmosphere from 4.9 × 102 mol to 6.7 × 101 mol. Additionally, annual savings of 4,793,472 kg of natural gas and 13,884 tons of steel were estimated owing to reduced oxidation losses. The proposed air–fuel adjustment led to estimated annual energy savings (equivalent to 4,793,472 kg of natural gas) and a reduction in material loss due to oxidation from 4.5% to 3.75% (an absolute reduction of 0.75 percentage points; relative reduction ≈ 16.7%), which has a significant industrial impact on metal conservation and descaling cost reduction. Although there are CFD studies on plate overheating and scale growth separately, this work presents three main contributions: (1) the integration, within a single numerical framework, of combustion, radiation, species transport, oxidation kinetics, and actual plate movement using a dynamic mesh; (2) validation against continuous industrial records (16 thermocouples) and quantification of operational benefits such as fuel savings and reduced material loss; and (3) a comparative analysis between actual and optimized conditions, which standardize the air–methane ratio.