The Study of Cooling Process on Runout Table by Simulation Method

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.


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.
This research is an extension of the previous research [1].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 [2] [4,5].The mathematical model [4] 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 [5] 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 [6] 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 [7].Reference [8] 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 [9] 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.

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 [5]: where  is the density,  is the temperature dependent specific heat, and  is the temperature dependent thermal conductivity of the material.  is strip velocity. and  are temperature and time, respectively; q is the heat generation term representing the internal heat source released due to phase transformation.

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): where  is the normal direction of strip's surface.ℎ  is the coefficient of convection heat transfer.  is the temperature of water. is the strip temperature.Boundary condition at inlet region is Boundary condition at outlet region is Initial condition is set by (5) as where () is the temperature distribution in  direction. is the runout table length in  direction.

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: 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): Time-dependent problem is solved by using implicit time stepping scheme.The step time (Δ) 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 ,  coordinates in [12,13]: The final discretization approximated by using implicit time stepping method can be expressed by ( 9) by neglecting the source term q ; then, 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   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 / is used, and the backward derivative for /,  2 / 2 and  2 / 2 is used.This results in upwind scheme:  The formulation of each numerical method is shown below.Δ and Δ are the space between two space grid points or the space step size, and  is used to represent Δ, the time step.The above equation ( 10) is modified by evaluating the space backward derivative at time step  + 1 instead of at time step ; this results in ( 12) By rearranging that discretization equation, we finally get the final result in (13): Rearranging new form, it can be represented by ( 14): where where  = / is the constant of material; this will show stable if   /Δ ≤ 1, /Δ 2 ≤ 1, and /Δ 2 ≤ 1. Writing it in matrix form it results in ( 16): . . .
where  is the number of space grid points.The above equation ( 16) is written in linear algebraic equation form as expressed in (17): To determining x, which represents the solution at time step  + 1 or at time  = ( + 1).b represents the current solution at time step , and A is the coefficients of matrix.

Simulation Model
In simulation, we set boundary and condition of 2D heat transfer model as illustrated in Figure 2. Details are listed as follows: (i) 2D heat conduction transfer in transient analysis, 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 (  ), external fluid velocity (  ) and temperature of cooling water (  ) respectively.These affect strip temperature control at entrance and exit ROT, and it can show in Table 1 for case study in simulation.
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 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.

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.

Table 1 :
Case of simulation studied on changing three variables.

Table 2 :
Output of surface temperature at ROT.