ISRN.INDUSTRIAL.ENGINEERING ISRN Industrial Engineering 2314-6435 Hindawi Publishing Corporation 820592 10.1155/2013/820592 820592 Research Article The Study of Cooling Process on Runout Table by Simulation Method http://orcid.org/0000-0003-4460-7480 Suebsomran Anan 1 Butdee Suthep 2 Dobmann G. Jelemensky L. 1 Department of Teacher Training in Mechanical Engineering King Mongkut’s University of Technology North Bangkok 1518 Pracharat 1 Road, Wongsawang, Bangsue, Bangkok 10800 Thailand kmutnb.ac.th 2 Department of Production Engineering King Mongkut’s University of Technology North Bangkok 1518 Pracharat 1 Road, Wongsawang, Bangsue, Bangkok 10800 Thailand kmutnb.ac.th 2013 23 10 2013 2013 29 05 2013 02 09 2013 2013 Copyright © 2013 Anan Suebsomran and Suthep Butdee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This research aims to achieving the effective cooling parameter on the runout table (ROT) of strip steel in hot rolling process. The 2-dimensional transient heat conduction is developed including the external force convection and heat source due to translational motion. The material property, boundary, and initial condition are defined and bounded to model geometry. The strip velocity, cooling water temperature, and external fluid velocity are chosen for the influent parameters during cooling process at ROT. To find the optimality of cooling operating requirement, simulation study is conducted throughout this research. To reach the objective of optimal cooling consumption at ROT, temperature distribution in the strip steel during cooling governs by the form of heat transfer equation. To solve 2-dimensional transient heat conduction by numerical methods, the backward difference formula (BDF) applies to discretization of partial differentiation equation (PDE). The parallel sparse direct linear solver (PARDISO) and conjugate gradients method are comparatively applied to computation in linear algebraic equation. The simulation studies are divided into 12 case studies with three variations subjected to cooling conditions at ROT. From simulation results, the range of such three variations can be identified in relation to economic cooling system and desired quality of products.

1. Introduction

Recently, advanced material processing technology must have become suitable for low cost of production, high productivity, and better quality of product. Manufacturing process of steel making is a long distance process. Slap products are passed to several machines to gain the desired size of product such as roughing mill machine, finishing rolling mill. After finishing rolling stand, the steel strip type is defined at this stage. Mechanical and physical properties of steel strip are controlled for the desired product quality. Temperature is one of the main parameters to control the product properties. After finishing stand process, the next process of the strip steel will arrive to runout table (ROT) as depicted in Figure 1. Runout table (ROT) employs to cool down the temperature approximately from 800–950°C at the entrance to the ROT to 450–600°C at exit from the ROT. If temperature is cooled down linearly as described, then phase transformation of strip steel will change from austenite to ferrite of metallurgy structure. The water wall uses for cooling the strip temperature by nozzle jet in both bottom and top of strip surface. To save the water consumption in cooling process at ROT, the optimal cooling of control parameters during cooling is required to conduct this research. In this study, three variations such as strip velocity, external fluid velocity, and cooling water temperature subjected to cooling conditions at ROT are used for studying the effective cooling variables. Problem of cooling at ROT process for the optimal cooling variables is derived for operational method in practical manner.

Layout of hot rolling process of strip steel production.

This research is an extension of the previous research . The extension of research is to study the influence of two kinds of linear algebra computation influent to the temperature of strip steel. Many researches attempt to analyse the cooling system at runout table. Reference  implemented the online runout table model in a hot steel mill. Coiling temperature is the main parameter to perform the online cooling at runout table such as using the inversion temperature control technique. Prediction of coiling temperature is implemented by online runout table model using PLC to control the condition of cooling in each segment with the appropriate cooling rates. This research developed only 1D conduction transient heat transfer and lacked consideration of strip velocity. Thus, prediction of temperature by online runout table model would be not exactly accurate in cooling rate control. Reference  also predicted the thermal and metallurgical behavior by using finite element model at runout table. Finite element coupled analysis of the thermal and metallurgical behavior of strip occurring at runout table for hot strip rolling process. Simulation result is compared with the thermal measurements in validation. The simulation model revealed that the effect of diverse process parameters is demonstrated in this research. And the prediction of temperature and phase transform at runout table is also described by [4, 5].

The mathematical model  is based on finite element method and the Scheil additive rule to prediction thermal and phase transformation. The mathematical and simulation model are not conflicted in the present work. The mathematical model is described in rectangular coordinate, but simulation results show the cylindrical coordinate. Reference  proposed the model of heat transfer equation of moving strip, but this proposed model studied only heat transfer in the steady state condition. From simulation result, when we increase the strip speed, the slop of cooling curve will decrease. Then, the starting point of ferrite transformation moves to the end of runout table position. This research demonstrated that the effect parameters of thermal and phase transform are affected by strip speed and initial gain austenite size of microstructure of metal changed to ferrite microstructure for a low carbon steel. Reference  also presented the modeling of temperature history and phase transforms during cooling of steel by using the finite element model. The simulation in prediction of temperature and phase transformation is compared to the experiment methods to validate the proposed prediction model. The model of deformation, temperature, and phase transformation behavior of strip steel on runout table is presented in the numerical model to simulate the followed behavior. The thermal model is formulated by using finite element and heat transfer coefficient of strip from actual mill data. The additional model included the deformation behavior into developed model proposed by . Reference  presented the prediction accuracy of strip temperature due to the effect of phase transformation latent heat in laminar cooling. The model was used to calculate the strip temperature and volume fraction of steel phase during cooling at runout table into different situations by regarding and disregarding phase transform latent heat. Comparison of prediction and experiment leaded to the errors in temperature prediction. The different steel grades showed that that the higher the carbon contents of steel, the higher the error of disregarding phase transformation latent heat. Not only the essential study on prediction temperature and phase transformation, but also  extended to describe the calculation of thermal stress affecting strip flatness change during runout table cooling in hot steel strip rolling in the effect of flatness of strip. The finite element was used to analyze the thermal stress during cooling. The commercial finite element software ABAQUS is used to analyse and calculate thermal behavior. References [10, 11] studied the cooling system. Experiment is set up to study the effect of cooling parameters at runout table. Thus control design of cooling system at runout table planned to implement to an industrial plant. From the literature, we found that most of researches are aimed to study the thermal and metallurgical behavior of strip during cooling at runout table. Method of prediction of temperature and phase transformation is mostly applied for the design of cooling process. To reach the desired product quality, temperature must be controlled according to the appropriation of the mechanical, physical, and metallurgy property. To design the cooling system at runout table, the heat transfer must be modeled in strip steel. An amount of heat flux in the strip steel must be evaluated in each condition that occurs on the actual industrial operation conditions. In this research, we have the goal of finding the effect of cooling parameters of strip steel at ROT.

2. Mathematical Model 2.1. Conduction Heat Transfer Model

Mathematical model of heat conduction in strip at runout table in rectangular coordinate can be formulated to describe the thermal behavior. The physical system is modeled in two-dimensional heat transfer equation of moving strip. The transient analysis of heat transfer of moving strip can be described by (1) as explained by : (1)x(kTx)+z(kTz)+q˙=ρcuiTx+ρcTt, where ρ is the density, c is the temperature dependent specific heat, and k is the temperature dependent thermal conductivity of the material. ui is strip velocity. T and t are temperature and time, respectively; q˙ is the heat generation term representing the internal heat source released due to phase transformation.

2.2. Boundary and Initial Conditions

The differential equation needs to be solved by numerical methods. The heat conduction occurs on the slap surface. The convection heat transfer applied the force convection and natural convection using water and air, respectively. At runout table, conduction of heat transfer in strip surface can be set to boundary conditions by (2): (2)-kTz=hw(T-Tw), where z is the normal direction of strip’s surface. hw is the coefficient of convection heat transfer. Tw is the temperature of water. T is the strip temperature.

Boundary condition at inlet region is (3)T(0,z)=T(z). Boundary condition at outlet region is (4)Tx|x=L=0. Initial condition is set by (5) as (5)T(x,z,t=0)=T0,0, where T(z) is the temperature distribution in z direction. L is the runout table length in x direction.

3. Numerical Solution 3.1. Discretization Methods

Backward difference formula (BDF) is applied to estimate the first and second order derivative equation as shown in (6). Thus, discretization can be evaluated using the continuous time derivation: (6)f(xi)  f(xi)-f(xi-1)h+O(h),f(xi)f(xi)-2f(xi-1)+f(xi-2)h2+O(h). Two-dimensional heat transfer equation in (1) can be approximated by a simple backward difference formula (BDF). In this model, partial differential equation is approximated with the proper time stepping interval in (7): (7)2Tx  Ti,j-2Ti-1,j+Ti-2,jΔx2+O(Δx),2TzTi,j-2Ti,j-1+Ti,j-2Δz2+O(Δz),TxTi,j-Ti-1,jΔx+O(Δx),TtTi,j-Ti-1,jΔt+O(Δt). Time-dependent problem is solved by using implicit time stepping scheme. The step time (Δt) must be small enough with uniform grid in spatial rectangular domain, and it can be calculated for uniform conductivity by (8) for equal space grid in spatial coordinate in reference x,  z coordinates in [12, 13]: (8)Δt<ρc(Δx)22k. The final discretization approximated by using implicit time stepping method can be expressed by (9) by neglecting the source term q˙; then, (9)Ti,j-2Ti-1,j+Ti-2,jΔx2+Ti,j-2Ti,j-1+Ti,j-2Δz2=ρckuTi,j-Ti-1,jΔx+ρckTi,j-Ti-1,jΔt.

The update time derivative of model can be calculated to the specification of time stepping interval. Based on implicit time stepping method, the iterative solution can be formulated by using the finite difference methods. Since ui is the strip velocity with unidirectional translational motion, the PDE for heat transfer referred in (1) can be solved by implicit finite difference equation by using the upwind scheme which follows to discretization in (10).

Forward time derivative for T/t is used, and the backward derivative for T/x, 2T/x2 and 2T/z2 is used. This results in upwind scheme: (10)ρckTi,jn+1-Ti,jnτ=Ti,jn-2Ti-1,jn+Ti-2,jnΔx2+Ti,jn-2Ti,j-1n+Ti,j-2nΔz2-ρcuikTi,jn-Ti-1,jnΔx,(11)Ti,jn+1-Ti,jn=kτρcΔx2(Ti,jn-2Ti-1,jn+Ti-2,jn)+kτρcΔz2(Ti,jn-2Ti,j-1n+Ti,j-2n)-kτρcρcuikΔx(Ti,jn-Ti-1,jn). The formulation of each numerical method is shown below. Δx and Δz are the space between two space grid points or the space step size, and τ is used to represent Δt, the time step. The above equation (10) is modified by evaluating the space backward derivative at time step n+1 instead of at time step n; this results in (12) (12)Ti,jn+1-kτρcΔx2(Ti,jn-2Ti-1,jn+Ti-2,jn)-kτρcΔz2(Ti,jn-2Ti,j-1n+Ti,j-2n)+kτρcρcuikΔx(Ti,jn-Ti-1n)=Ti,jn. By rearranging that discretization equation, we finally get the final result in (13): (13)Ti,jn=Ti,jn+1-τkρcΔx2Ti,jn+1-τkρcΔz2Ti,jn+1+τuiΔxTi,jn+1+τkρcΔx22Ti-1,jn+1-τuiΔxTi-1,jn+1+τkρcΔz22Ti,j-1n+1-τkρcΔx2Ti-2,jn+1-τkρcz2Ti,j-2n+1. Rearranging new form, it can be represented by (14): (14)apTi,jn+1+aw-1Ti-1,jn+1+aw-2Ti-2,jn+1+as-1Ti,j-1n+1+as-2Ti,j-2n+1=Ti,jn, where (15)ap=1-ταΔx2-ταΔz2+τuiΔx,aw-1=2ταΔx2-ujτΔxaw-2=-ταΔx2,as-1=2ταΔz2as-2=-ταΔz2, where α=k/ρc  is the constant of material; this will show stable if τui/Δx1, τα/Δx21, and τα/Δz21. Writing it in matrix form it results in (16): (16)[apas-1as-20000aw-1apas-1as-2000aw-2aw-1apas-1as-2000aw-2aw-1apas-1as-20000aw-2aw-1apas-100000aw-2aw-1ap00000000ap]×[T0,0n+1T1,1n+1T2,2n+1T3,3n+1T4,4n+1TN-1,M-1n+1TN-2,M-2n+1]=[T0,0nT1,1nT2,2nT3,3nT4,4nTN-1,M-1nTN-2,M-2n], where N is the number of space grid points. The above equation (16) is written in linear algebraic equation form as expressed in (17): (17)Ax=b. To determining x, which represents the solution at time step n+1 or at time t=(n+1)τ. b represents the current solution at time step n, and A is the coefficients of matrix.

4. Simulation Model

In simulation, we set boundary and condition of 2D heat transfer model as illustrated in Figure 2. Details are listed as follows:

2D heat conduction transfer in transient analysis,

strip dimension 0.016×47 m,

steel type AISI 4340,

convective cooling by water at total plate length 47 m with average heat transfer coefficient (hw) by external force convection,

initial value for strip temperature 1148 K,

T1=1148 K,

T4=793 K,

no heat source Q,

no heat flux generation at all boundaries,

strip velocity: Vs (m/s),

external fluid velocity: Uw (m/s),

temperature of cooling water: Tf (°K).

Boundary and condition setting.

This research utilized the Comsol Multiphysics software to build the model and solve with the numerical methods as in [14, 15]. The developed model set the configuration as shown in Figure 2. In this study, the changed variables for simulation are included by three parameters such as strip velocity (Vs), external fluid velocity (Uw) and temperature of cooling water (Tf) respectively. These affect strip temperature control at entrance and exit ROT, and it can show in Table 1 for case study in simulation.

Case of simulation studied on changing three variables.

Case study Strip velocity (Vs)(m/s) Fluid velocity (Uw)(m/s) Cooling water temperature (Tf)(°K)
1 V s = 10 U w = 20 T f = 288
2 V s = 10 U w = 20 T f = 293
3 V s = 10 U w = 30 T f = 288
4 V s = 10 U w = 30 T f = 293
5 V s = 7 U w = 20 T f = 288
6 V s = 7 U w = 20 T f = 293
7 V s = 7 U w = 30 T f = 288
8 V s = 7 U w = 30 T f = 293
9 V s = 0 U w = 20 T f = 288
10 V s = 0 U w = 20 T f = 293
11 V s = 0 U w = 30 T f = 288
12 V s = 0 U w = 30 T f = 293

At runout table of cooling process of strip steel, we specify the strip temperature at entrance ROT and exist ROT indicated in Table 2. To keep the desired product quality and requirement, cooling process requires to control the strip temperature after exist ROT as shown in Table 2. The aim of simulation study verifies the minimum strip temperature error compared to the reference strip temperature at exit ROT as shown in Table 2. Thus, percent true relative error is useful for computation by (18): (18)εt=True  errorTrue value×100. From simulation results, we found that the minimum percent relative error of surface temperature at cooling process at ROT is indicated by the case study number 1 and 2. The largest percent true relative error is in the case study number 11 and 12. The output of strip temperature from simulation program on case study number 1 by using Comsol Multiphysics software is illustrated in Figures 3 and 4.

Output of surface temperature at ROT.

Case study no. Entrance ROT temperature (°K) Reference at exit ROT temperature (°K) Actual surface temperature (°K) True error Percent relative error (εt)
1 1148 793 754.31 38.69 0.048789407
2 1148 793 720.7 72.3 0.091172762
3 1148 793 458.85 334.15 0.421374527
4 1148 793 455.25 337.75 0.42591425
5 1148 793 679.96 113.04 0.142547289
6 1148 793 455.25 337.75 0.42591425
7 1148 793 679.96 113.04 0.142547289
8 1148 793 503.6 289.4 0.364943253
9 1148 793 441.85 351.15 0.442812106
10 1148 793 440.4 352.6 0.444640605
11 1148 793 375.32 417.68 0.526708701
12 1148 793 378.7 414.3 0.522446406

Temperature of strip using the parallel sparse direct linear solver (PARDISO).

Temperature of strip using conjugate gradients method.

5. Conclusions

This research aims to study the effective cooling parameters at ROT of strip steel in hot rolling process. The desired quality of product must be achieved. In cooling process, the efficiency at ROT must be optimal to the customer requirement of strip steel grade. To find the optimal operating variables of ROT’s cooling process, we develop the 2-dimensional transient heat transfer of strip steel developed by using mathematical model. Boundary and initial condition are bounded valuable by considering the practical conditions. Numerical solution is applied to solve the mathematical model built using a Comsol Multiphysics software by heat transfer module. Backward difference formula (BDF) applies to discretization of partial differentiation equation (PDE). The parallel sparse direct linear solver (PARDISO) and conjugate gradients method are comparisons for computation of linear algebraic equation. There are three variables useful for each simulation of case studies such as strip velocity, external fluid velocity, and temperature of cooling water, respectively. From the simulation results, the minimum percent relative error of strip surface temperature compared to strip temperature from finishing stand is the case study number 1 as depicted in Table 2. The most suitability in operating control variables of ROT cooling system is achieved subject to industrial conditions of strip steel in hot rolling. Finally, the result has showed that the parallel sparse direct linear solver (PARDISO) is more proper than conjugate gradients method for linear algebra computation method.

Acknowledgments

The work described in this paper was partially supported by the Danieli Far East Co. LTD, Thailand, and Science and Technology Research Institute (STRI), KMUTNB, Thailand.

Suebsomran A. Butdee S. Cooling process on run out table by simulation method Case Studies in Thermal Engineering 2013 10.1016/j.csite.2013.07.002 Mukhopadhyay A. Sikdar S. Implementation of an on-line run-out table model in a hot strip mill Journal of Materials Processing Technology 2005 169 2 164 172 2-s2.0-26444478146 10.1016/j.jmatprotec.2005.04.039 Sun C. G. Han H. N. Lee J. K. Jin Y. S. Hwang S. M. A finite element model for the prediction of thermal and metallurgical behavior of strip on run-out-table in hot rolling ISIJ International 2002 42 4 392 400 2-s2.0-0036072704 Serajzadeh S. Prediction of temperature distribution and phase transformation on the run-out table in the process of hot strip rolling Applied Mathematical Modelling 2003 27 11 861 875 2-s2.0-0142088995 10.1016/S0307-904X(03)00085-4 Serajzadeh S. Prediction of temperature variations and kinetics of austenite phase change on the run-out table Materials Science and Engineering A 2006 421 1-2 260 267 2-s2.0-33645920655 10.1016/j.msea.2006.01.071 Serajzadeh S. Modelling of temperature history and phase transformations during cooling of steel Journal of Materials Processing Technology 2004 146 3 311 317 2-s2.0-1242276101 10.1016/j.jmatprotec.2003.11.010 Han H. N. Lee J. K. Kim H. J. Jin Y.-S. A model for deformation, temperature and phase transformation behavior of steels on run-out table in hot strip mill Journal of Materials Processing Technology 2002 128 1–3 216 225 2-s2.0-0037032146 10.1016/S0924-0136(02)00454-5 Edalatpour S. Saboonchi A. Hassanpour S. Effect of phase transformation latent heat on prediction accuracy of strip laminar cooling Journal of Materials Processing Technology 2011 211 11 1776 1782 2-s2.0-79961170158 10.1016/j.jmatprotec.2011.05.027 Wang X. Yang Q. He A. Calculation of thermal stress affecting strip flatness change during run-out table cooling in hot steel strip rolling Journal of Materials Processing Technology 2008 207 1–3 130 146 2-s2.0-50949131671 10.1016/j.jmatprotec.2007.12.076 Horsky J. Raudensky M. Kotrbacek P. Experimental study of long product cooling in hot rolling Journal of Materials Processing Technology 1998 80-81 337 340 Leocadio H. Cesar Passos J. da Silva A. F. C. Heat transfer behavior of a high temperature stell plate cooled by a subcooled impringing circular water jet Proceedings of the 7th ECI International Conference on Boiling Heat Transfer May 2009 Florianópolis, Brazil Chapra S. C. Canale R. P. Numerical Methods for Engineers with Programming and Software Applications 2009 McGraw-Hill Patankar S. V. Numerical Heat Transfer and Fluid Flow 1980 Taylor and Francis COMSOL Multiphysics User’s Guide, VerSion 4. 1 1998-2010 COMSOL AB COMSOL Multiphysics Reference Guide, Version 4. 1 1998-2010 COMSOL AB