Coordinated Optimal Control of Secondary Cooling and Final Electromagnetic Stirring for Continuous Casting Billets

Secondary cooling and final electromagnetic stirring (F-EMS) are both key technologies for continuous casting. 5ese parameters are usually optimized and controlled separately which caused internal quality fluctuations in unsteady conditions. In this paper, a coordinated optimal control strategy based on a multiobjective particle swarm optimization (MOPSO) algorithm is proposed for the parameter optimization of secondary cooling and F-EMS, which is solved based onmultiobjective particle swarm optimization (MOPSO) algorithm. 5e solidification and heat transfer model are developed for the computation of billet temperature and the solidification, and the adaptive grid method is used to improve the diversity and robustness of optimal solutions. 5e secondary cooling water and F-EMS’ stirring current are dynamically controlled based on the optimization results. 5e results of field trials showed that the maximum carbon segregation and other quality indexes of billets can be improved significantly.


Introduction
Continuous casting of steel is a process that molten steel is solidified into "semifinished" products, such as billets, blooms, and slabs, which are subsequently sent for rolling in the mills. At the beginning of continuous casting, as steel composition is qualified, molten steel is transported to the continuous caster with a ladle. From the ladle, molten steel is poured into a tundish and then into copper molds with a certain shape, which is cooled by high-speed water. In the mold, primary cooling occurs and makes a supporting solidified shell with molten steel inside. en, the strand is continuously withdrawn through a series of secondary cooling zones (SCZs) where water cooling happens with sprays. When exit from the secondary cooling zones, the strand is cooling off by radiation in the air. And finally, it is cut into fix-length strips and transported to rolling mills for further processing. e process is schematically shown in Figure 1.
In the continuous casting process, electromagnetic stirring is commonly used to improve the billet quality, including mold electromagnetic stirring (M-EMS), secondary cooling electromagnetic stirring (S-EMS), and final electromagnetic stirring (F-EMS). Among them, F-EMS is supposed to be the key to improve the central quality of products, such as central segregation and porosity [1]. Otherwise, secondary cooling plays an important role in temperature distribution, greatly influencing the billet quality, such as internal cracks, segregation, and porosity.
However, secondary cooling and F-EMS are conventionally treated as two independent technologies without considering their interaction, and their parameters are generally optimized separately [2][3][4]. In the conventional model, secondary cooling water flows are determined by steel grade and size of billet and optimized as functions of casting speed for a specific caster. In the optimization, the corresponding constraint conditions include metallurgical criteria and serial control goals, and most of them are concerning about temperature distribution. Otherwise, the parameters including the installation position of F-EMS are determined according to the comprehensive situation of the casting machine, which often remains unchanged in the production process. A billet caster often manufactures billets with different sizes and different steel grades, which results in the change of liquid core at the F-EMS position. Even for the same steel, the liquid core also changes due to factors such as production rhythm and molten steel temperature variety. According to the changing situation, if stirring parameters remain unchanged, obviously the F-EMS cannot play a proper role to keep the stability of solidification or even cause side effects. However, to keep the stability of the liquid core just by adjusting the secondary cooling water flows, that may lead to excessive temperature deviations from the target value and cause cracks and other quality problems [2]. Otherwise, if stirring parameters change according to the size area of the liquid core at F-EMS position without adjusting the water flows, as the liquid core is too big or too small, FEMS also cannot achieve its proper role. erefore, the coordinated optimal control of secondary cooling and F-EMS is practically significant to improve the quality of casting billet.
In the coordinated optimal control problem, there are some difficulties including complex solution space and nonlinear and nondifferentiable objective functions due to the solidification process, as the steel solidification process belongs to the Stefan problem. To solve this multiobjective optimization problem (MOOP), conventional optimization methods are not suitable. In the past few years, swarm intelligence techniques are becoming increasingly popular for improving the process, including genetic algorithm [5,6], particle swarm optimization (PSO) algorithms [7], and evolutionary multiobjective optimization algorithm [8,9]. Among these techniques, PSO is most attractive for its high convergence speed and better exploration results. Normally, the MOOP is converted into a single objective optimization problem with classical weighted sum method by using a single pair of fixed weights and only one point on the Pareto frontier can be obtained, which is not suitable for the situation that if the optimal solutions lie on nonconvex or discontinuous regions of function space. erefore, to obtain the global Pareto optimum, a well-spread and diverse Pareto solution frontier must be derived firstly, and it requires that the algorithms can be executed iteratively to ensure that any combination of weights is available. PSO is an intelligent evolutionary algorithm proposed by Eberhart and Kennedy [10], which is an effective metaheuristic for solving portfolio selection and then the successful applications motivated the researchers to apply it to multiobjective optimization problems [11,12]. When developing multiobjective particle swarm optimization (MOPSO) algorithm, we will focus on the following three questions: (1) taking into account the diversity and ultimately convergence of Pareto fronts; (2) choosing the globally optimal solution (gBest) and local optimal solution (pBest) to guide the swarm particles to search for new solutions; and (3) try to be convergent and avoid falling into the local Pareto front.
In this paper, aiming at improving the production quality of continuous casting billets such as macrosegregation of high-carbon steels, a coordinated multiobjective optimization and control strategy is developed based on solidification and heat transfer model. e optimal solution and research results are applied to the dynamic control of the actual billet caster and good effects have been achieved.

Overview of the Optimal Control Strategy for Quality Control of Continuous Casting Billets
e framework for optimal control of secondary cooling and F-EMS is shown in Figure 2. It includes two parts: the coordinated optimization part and the dynamic control part. Above all, the heat transfer and solidification model serve as the basis of both parts, providing calculations of solidification states such as surface temperatures and radius of the liquid core. In the optimization part, a multiobjective particle swarm optimization (MOPSO) algorithm is adopted for coordinated optimization of secondary cooling and F-EMS. In the algorithm, the heat transfer and solidification model are used as the basis and the adaptive grid method is used to improve the diversity and robustness of solutions. In the dynamic control part, secondary cooling water flows and F-EMS stirring current are dynamically controlled based on the online model's calculations to avoid fluctuation of surface temperature and flow velocity in the solidification frontier, which benefits the quality stability.

Basis for Optimal Control Strategy-Heat
Transfer and Solidification Model e heat transfer process of continuous casting billets can be generally described by the heat transport equation. And the following assumptions are considered: (1) increasing heat conduction induced by the flow of nonsolid steel can be taken into account by an increased effective thermal conductivity; (2) the latent heat released during phase transformations can be considered in the equivalent specific heat;  and (3) heat conduction in the casting direction can be ignored as it is quite small compared to the transverse crossdirections. e equation takes the following form in the follow-up coordinate system: where T is the billet temperature, t is the time, x and y are the coordinates, and ρ is the density of steel; c eff � dH/dT, in which H is the enthalpy; k eff � f s k s + m (1 − f s ) k l , in which f s is the solid fraction and k s , k l are separately the thermal conductivity of solid phase and liquid phase. e thermal properties including ρ, H, f s , k s , and k l are assumed to be functions of T and of the steel compositions and can be calculated by a special algorithm based on pseudobinary phase diagram [13].
Considering the symmetry, only one-quarter of the billet slice (transverse cross-section) is selected as the calculation domain. Each slice starts at meniscus (assumed flat) and moves along the strand through the mold, secondary cooling zones, and then the air cooling zone. e initial condition for each slice has been assumed to be the same as the casting temperature: Moving down from the meniscus, each slice sequentially experiences the boundary conditions specified as follows.
In the mold, a type of decreasing heat flow due to air gap forming is considered [14]: where t is the experienced time of slice in the mold. In SCZ, water cooling and thermal radiation are considered: where h i is the heat transfer coefficient, ε is the emissivity, σ is the Stefan-Boltzmann constant, T w is the cooling water temperature, and T air is the air temperature. ere can be two kinds of spays for different cooling zones. In the water spray cooling zones [15], where w is the water flow density (L/m 2 /s). And in the airmist spray cooling zones, In equation (3), (5), and (6), A, B, α i , and C i are machinedependent parameters. In the air cooling zone, thermal radiation is considered: In this study, the model is solved by FVM-ADI algorithm, discretized by the finite volume method, and solved by alternative direction implicit algorithm [14]. e temperatures and solidifications are computed as shown in

Optimization Objective Functions.
e billet quality depends on the behavior of the billet surface temperature and the flow of the solidification frontier. Undercooling of the strand can result in too long a liquid pool, and overcooling can lead to internal cracks. ere must be a smooth transition of the surface temperature, with minimum reheating, as the steel passes from one cooling zone to another. e secondary cooling plays a major role in the  Figure 2: Framework of the optimal control strategy.
Journal of Control Science and Engineering casting process because the intensity of the water sprays affects highly the solidification rate.

Billet Surface
Temperature. e secondary cooling zone of the billet caster is divided into several sections. In each section, cooling water is dispersed to the billet surface to control the temperature filed. As target surface temperature is specified in each section, water flow rates Q i (i � 1, 2, ..., N) of the secondary cooling zone should make billet surface temperatures match the target temperatures as closely as possible. So the cost function can be defined as follows: where N denotes the number of sections and T i and T * i are the calculated billet surface temperature and a target temperature of section i, respectively.

Liquid Core Size.
ere is also a requirement for liquid core size, namely, the radius of the liquid core at the center position of FEMS along the casting direction. To ensure the stability of liquid core size, the cost function f 2 is defined to make the radius of the liquid core as close as possible to the target value: 4.1.3. Gradients of Surface Temperature. Besides, the gradients of the strand surface temperature along the casting directors should be controlled. us, the cooling rate C N (°C/ m) and reheating rate C p (°C/m) are both considered. e cost function f 3 is defined as follows: e optimization task is to minimize f 1 , f 2 , and f 3 according to possible cooling patterns, namely, water flow rate settings. As the three objective functions are conflicting, it is reasonable to handle this optimization problem by the multiobjective algorithm.
For different steel grades, the concerned quality problems to be solved are also different. So, the weight factors in the process of optimization should be changed accordingly.

Advanced Particle Swarm Optimization Algorithm.
To solve the above optimization problem, an advanced particle swarm optimization algorithm is proposed. In the particle swarm optimization algorithm, it is supposed that the search space is d-dimensional and there are n particles in the swarm. Particle i is located at position X i � (x i1 , x i2 , . . ., x id ) and has velocity V i � (v i1 , v i2 , . . . , v id ). e update equations for velocity and position of the d-dimension and t generation are given as follows: where v i (t) is the current velocity of i particle, i � 1, ..., n; n is the population size; the superscript i is the dimension of the particle; p i is the best previous position of the i-th particle; p g is the best previous position among all the particles in the swarm; x i (t) is the current position of the i-th particle. i � 1, 2, . . ., n; w p is nonnegative (inertia weight); c 1 and c 2 are the acceleration coefficients; r 1 and r 2 are the random numbers generated uniformly in the range [0, 1]; t is the current number of iterations.
By decreasing the inertia weight (w p ) from a relatively large value to a small value, the PSO tends to have better global search ability at the beginning and better local search ability near the end [3]. So, the inertia weight is set according to the following linear expression: where N max is the maximum number of iterations and N iter is the current number of iterations. Equation (14) restricts the value w to the range w max , w min . In this study, the maximum and minimum values of the inertia weights are set as 0.9 and 0.4, respectively.
To avoid the premature problem, mutation operation is adopted when the velocity of all particles is less than the limit value. e update functions are as follows [8]: where v m is the mutation value, β is the mutation coefficient between 0 and 1 used to adjust mutation degree, r 3 is the random number between 0 and 1, and x id is the d-th dimension selected by random of i-th particle.

Pareto Optimal Solutions Based on Adaptive Grid Method.
To maintain diversity in the Pareto optimal (nondominated) solutions, a novel mechanism with an adaptive grid is For optimization problem with r objectives, borders of the grid are set to 2r and a grid can be split into several small hypercubes (HC). e specific number of divisions is determined by the evolutionary population size and number of objectives. Supposing each HC as r i , i � i 1 , i 2 , . . ., i r and i k ∈ 1 . . . d, d is a constant number of divisions for each dimension.
In Figure 4, it is supposed that d � 6 and then corresponding to the side of each r i can be expressed as follows: where ω k is the width of each region in k-dimension, Grid is a cube under the three-dimensional space for the three-dimensional objective function. e number of particles within each cube is determined by requiring a selection of the border and population size.
Let nondominated solutions set at generation t was obtained as A t . e region width range k of each dimension is set as follows: Adaptive grids in the grid boundaries are not fixed. According to the present and the previous generation distribution of the individual, each generation evolution adaptively adjusts the boundary. If the newly generated individuals are located outside the boundaries of the previous generation and nondominated, it will be added to the archive set. If the archive set is full, other individuals will be deleted to ensure the distribution of the particles in the search space.

Pareto Optimal Solutions Weighting Computation.
After obtaining many solutions that are true Pareto optimal with uniform spread coverage, a simple procedure called pseudo-weight vector approach is employed, which facilitates final decision making. In this approach, a pseudoweight vector is calculated for each obtained solution [11]. e equation is given as follows: where w ij is the weight of the i-th objective function for the jth solutions, f max

Dynamic Control of Secondary Cooling for the Stability of Temperature.
After optimization under different casting speeds, the secondary cooling water flows were determined. e water flows of secondary cooling zones were precisely controlled by a dynamic feed-forward strategy, described as follows: where Q i is the water flow of secondary cooling zone i, β i is compensation coefficients of the superheat of cooling zone i, ΔT ei is the effective superheat and ΔT std is the standard superheat, v ei is the effective casting speed [15,16], and parameters a i , b i , and c i are optimized spray cooling coefficients by MOPSO.

Dynamic Control of F-EMS' Stirring Current for the Stability of Flow Velocity of Liquid Steel.
Stirring intensity of the solidification front can be expressed by the stirring velocity of molten steel. In the rotary stirring, the stirring velocity of the molten steel is calculated as follows [1]:

Journal of Control Science and Engineering
where K is the constant related to a magnetic stirrer, B is the magnetic induction of billet center (Gs), f is the frequency of stirrer (Hz), r is the liquid core radius (two-phase region + liquid region, cm), σ is the electric conductivity of molten steel (s/cm), ρ is the molten steel density (g/cm3), and μ is the molten steel kinematics viscosity (CP, 1CP � 0.001 Pa·s). e relationship between magnetic flux density (B) at the stirrer center and electric current (I) can be obtained from the measuring data and B was given as follows: (20) In the practical application, we first choose a suitable operating frequency of final electromagnetic stirrer and then determine the current value according to the corresponding liquid core radius and the target stirring velocity of the molten steel at the stirrer center. e current was calculated as follows:

Example Results
e above optimal control strategy for secondary cooling water and electromagnetic stirring has been applied in billet caster 1# of Sanming Iron and Steel Cooperation of China. Numerical experiments in MOOP of the casting process are performed for 65# steel with the billet size of 0.16 m × 0.16 m. e radius of the caster is 10 m. e caster has four cooling sections in the secondary cooling zone. e uppermost section is equipped with water sprays and the rest three are equipped with air-mist sprays. e length of the four sections is 0.37 m, 1.85 m, 2.20 m, and 2.32 m. e effective mold length is 0.90 m. FEMS is installed 7.70 m from the meniscus. e prescribed target billet surface temperatures are adjusted for different casting speeds. For example, the prescribed target surface temperatures at the end of four zones are sequentially 924°C, 928°C, 948°C, and 962°C at a casting speed of 2 m/min. e initial parameters in PSO are set as follows: the initial population of MOPSO is set to 50, the number of generations is set to 100, c 1 � 1.85, c 2 � 1.455, w � 1.0, and the number of nondominated solutions to be found is set to 50. e optimization problems are carried out for different casting speeds. e casting speeds used are changed from 1.6 m/min to 2.2 m/min with a step size of 0.1 m/min. Steady-state operation of the casting machine is assumed and optimization is performed in an off-line manner.
The optimization procedures were tested for both the solution quality and repeatability of results. e results of 50 nondominated solutions obtained using MOPSO after 100 iterations are shown in Figures 5 and 6. Several alternatives can be chosen at various satisfaction levels of the multiple objectives. Depending on steel grade and quality requirements and by analyzing the tradeoff between the multiple objectives, the researcher can make an appropriate decision.
To evaluate the performance of the proposed AD-MOPSO algorithm, these optimization results are compared to MOPSO without adaptive grid trimming. Due to the adaptive grid trimming and maintenance of the external archive set, the Pareto optimal solutions in Figure 6 have better distribution in the search space than those in Figure 5 and the Pareto distribution is more uniform. In the Cartesian coordinate system composed of three objective functions, the overall Pareto frontier distribution is in a plane shape. It can be seen from Figures 5 and 6 that AD-MOPSO gives a better distribution of true Pareto optimal solutions.
To reduce a large number of alternatives, the adaptive grid algorithm was used to filter the large set of nondominated solutions to a few representative solutions. Table 1 shows 30 representative Pareto optimal solutions for the three objective functions. Based on equation (21), the weights that can be given for each objective are shown in Table 1, which gives the objective values for each of the representative Pareto optimal solutions and their respective weight for each alternative. According to the different kinds of steel grade and quality requirements,  decision-makers can choose a suitable solution for optimization control.
In actual casting process control, quality indexes for different steel grades are also different, so the weight of the objective functions is also different. Since the 65# steel has high carbon, high strength, hardness, and elasticity and low cold plastic deformation characteristics, it is easy to form center segregation and center shrinkage. F-EMS is usually used to improve billet quality and a higher weight of the cost function f 2 is desired.
In the optimization process, the weight coefficients choice depends on the steel grade and quality requirements. After analyzing the available alternatives, based on the analytic hierarchy process, the final decision can be made. e weights 0.25, 0.50, and 0.25 for f 1 , f 2 , and f 3 are selected which is closer to that set of weights of 20th solution from Table 1.
e corresponding fitting parameters of the secondary cooling control model with quadratic equation (26) are shown in Table 2.
e dynamic control of secondary cooling greatly reduced the fluctuation of midface temperature at the end of secondary cooling zones with the changing casting conditions. e magnetic flux density of F-EMS is the maximum value at the frequency of around 8.5 Hz. While the frequency is greater than 10 Hz, the actual current will gradually deviate from (less than) the setting value and the magnetic flux density decreases, so the operating frequency is set at 9 Hz.
Besides, equation (10) has also been applied to the dynamic control of the electric current of F-EMS. To control the stirring velocity of 0.15 m/s induced by F-EMS stirrer in the mushy zone, the F-EMS current has been set as a function of mushy zone radius at the installation position of F-EMS. Figure 7(a) presents the tuning of F-EMS current with the changing casting conditions. e radius of the mushy zone is influenced by the casting speed and superheat. As the radius of the mushy zone is increasing, the flow in the mushy zone accelerates; to keep the flow speed in the mushy zone stable, the F-EMS current is supposed to decrease accordingly. In this paper, the F-EMS current has been set to be proportional to the reciprocal of the mushy zone radius, as shown in Figure 7. After the dynamic control strategy was applied in actual billet caster, the fluctuation of surface temperatures of billet in secondary cooling zones and the flow speed in the mushy zone were reduced. e defects of macrosegregation, porosity, and shrinkage cavity have been reduced significantly. Table 3 shows the distribution of billet internal defects of 180 sample macrographs.

Conclusions
In this paper, a coordinated optimal control strategy for secondary cooling and F-EMS has been proposed and applied to an actual caster. e optimized results are applied in the practical casting machine. Some conclusions could be drawn as follows.
An optimization method based on MOPSO for secondary cooling and FEMS of continuous casting is proposed and applied to the coordinated control of the billet continuous casting process. It integrates two independent control systems by coordinating weighting factors for various steel grades. e improved Pareto multiobjective particle swarm optimization algorithm considers both convergence and diversity of Pareto optimal population. It also improved the robustness in yielding efficient Pareto frontiers. A dynamic control model is developed and applied in the billet production process based on the coordinated control of secondary cooling water and stirring current of F-EMS. is ensures the stability of the surface temperature of the billets and the radius of the liquid core and reduces the fluctuation of surface temperature and flow velocity of liquid steel in the solidification frontier. Lots of industry trial results show that the carbon level is significantly improved. e maximum central segregation index is 1.06 and other internal quality problems are improved effectively.
Data Availability e data named "billet caster 1# of Sanming Iron and Steel Cooperation of China" used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, 12 months after publication of this article, will be considered by the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.