Parameterization optimal control of an unsteady partial differential equations with convection term by an improved three-term spectrum conjugate gradient algorithm

The unsteady partial differential equations (UPDE) with convection term gives a clear descriptions for the solidification process of a slab in dynamic production of continuous casting. To give a suitable setting value of secondary cooling water flow rate for the dynamic control system, this study investigates an optimal control problem (OCP) of UPDE with convection term. Firstly, control vector discretization of OCP and the solution of UPDE are given. Secondly, due to the rapidity for gradient, this paper analyzes the expression of the gradient calculation method based on Hamiltonian function costate system by approximate treatment, matrix calculation and composite trapezoidal integral method. Thirdly, an improved three-term spectrum conjugate gradient algorithm (ITSCGA) is proposed to solve the OCP of UPDE, and the global convergence of the ITSCGA is demonstrated. Lastly, the performance of ITSCGA is demonstrated by experimental simulations. The results demonstrate that the ITSCGA provides a smaller temperature fluctuations, and improves the quality of a slab.


Introduction
Partial differential equations (PDEs) are widely applied in many fields, such as physics, chemistry, biology, medicine, and so on, and it is encouraged to solve the mathematical problems in economics, finance, image processing, and other fields.Especially in metallurgy, the solidification heat transfer model can be used to represent the cooling process of a slab.Most secondary cooling water control systems [1][2][3][4] of continuous casting explore this PDE model to adjust the water flow rate setting value.Therefore, the study of optimization problem based on the solidification heat transfer model of a slab is crucial for the secondary cooling systems of continuous casting process.Partial scholars [5][6][7][8][9][10] based on the unsteady heat transfer model of a slab study the optimization method for the water flow rate setting value in the static production of continuous casting.However, in actual production, some dynamic situations of continuous casting often occur due to the early or tardy arrival of a ladle and the change of superheat, which can lead to adjustments in casting speed.Therefore, an unsteady heat transfer model with convection term (UHTMCT) of a slab attracts the attention of some scholars to deal with the dynamic situations.Based on the UHTMCT of a slab, Yu et al. 11 optimizes the water flow rate setting value in dynamic environment.Hardin et al. 12 proposes a dynamic spray cooling control model, and applies them to determine the optimal amount of water volume at each cooling zone.Wang et al. 13 introduces an openloop dynamic spray cooling control strategy, which uses a control algorithm to track the casting speed and temperature of each slice.
In the view of above researches, our previous works 14 study the nonlinear model predictive control (MPC) of the water flow rate setting value based on the UHTMCT of a slab.According to Canto and Banga 15 and Ping et al., 16 optimal control can be an effective method to obtain the control trajectory for nonlinear MPC implementation.Due to the large computation in NMPC applications, it is very essential to develop more efficient solvers for the nonlinear MPC of UHTMCT.Therefore, optimal control problem (OCP) of UHTMCT and an efficient solution method are investigated in this paper.The optimization efficiency is an important work to solve the OCP, thus a highly efficient parameterization optimal method and an improved three-term spectrum conjugate gradient algorithm (ITSCGA) are presented as a prominent solver.Furthermore, for the parameterization optimal method, in order to reduce the computational expense of gradients, the gradient of OCP for UHTMCT is analyzed based on a Hamiltonian costate system by the approximate treatments of numerical integration strategy.

Optimization problem of boundary conditions of PDEs
OCP of UHTMCT 17,18 in continuous casting belongs to the optimization of boundary conditions of PDEs.At present, optimization problems of boundary conditions of PDEs are divided into two categories: (1)  Boundary condition identification problem of PDEs.0][21][22] Duda et al. 23,24 solve the transient heat conduction equation based on the control volume finite element method and the control volume method, and estimate an unknown boundary condition through the transient temperature distribution.Sun et al. 25 use the sequential quadratic programming algorithm as an optimization technique to determine the surface heat flux of the billet.Huang et al. 26 estimate the time-varying heat flux on the chip boundary, based on three-dimensional unsteady PDE by using the inverse method; Ruan et al. 27 develop a stochastic particle swarm optimization method to determine the heat transfer coefficient in the continuous casting process.Torrijos et al. 28 study the inverse heat problem for transient heat conduction equation of the heat flux on the unknown surface of the molten salt loop.(2) The optimal control problem of boundary conditions of PDEs.0][31][32] Gavrikov et al. 33 study the optimal boundary control problem with linear quadratic cost function in the heat transfer process of a cylinder; Mechelli and Volkwein 34 study an optimal boundary and bilateral control problem of heat transfer equation with convection term; Dekhkonov 35 recently studies the boundary control problem, based on the heat transfer model, for the heat exchange process; Bollo et al. 36 solve OCPs of Neumann boundary for a n-dimensional heat equation; Abbasi and Malek 37 solve pointwise optimal control problem on and inside a tissue subject to thermal wave model with Dirichlet and Rubin boundary conditions.
The algorithms for solving the optimization problem of boundary conditions of PDEs are mainly divided into two categories.(1) Stochastic algorithms.The stochastic algorithms 38,39 have the random search's global optimization characteristic, but the disadvantages of these methods are slow convergence speed and high computational cost.(2) Gradient algorithm.There are a number of efficient numerical methods based on the gradient, such as Newton method, 40 alternating direction multiplier method, 41 and conjugate gradient algorithm (CGA).For the Newton method, 40 we have done the study on this method, which needs to calculate the second derivative of objective function, which is difficult to calculate when applied to our problem.Therefore, Newton method 40 is not suitable to this paper's problem.For the alternating direction multiplier method, we only find one literature, 41 which uses the alternating direction method of multipliers (ADMM) to solve the optimal control model constrained by one-dimensional partial differential equation.In Dawodu, 41 the PDE does not include the convection term, and this PDE can be discretized into a derived discrete convex optimization form amenable to the ADMM with the help of Crank-Nicolson and Composite Simpson's methods.However, the problem in our paper contains convection term, and has a higher dimension PDE, which leads to difficult to transform the original problem into the convex optimization form amenable to the ADMM.3][44][45][46][47][48][49][50][51] Zhou et al. 42 estimate the temperature and heat flux on inaccessible surfaces by using the conjugate gradient method with temperature and heat flux measured on back surface.Yang et al. 43 propose a modified conjugate gradient method to identify the physical parameters of transient heat conduction problems in aerospace engineering.Wang et al. 44,45 propose an improved CGA to solve the OCP based on onedimensional and two-dimensional parabolic PDEs; Razzaghi et al. 46 adopts an improved CGA to minimize cost function in the unsteady heat transfer model.They solve accurately an optimization problem by using CGA.8][49] .In view of the above research findings, our previous works 50,51 combine the regularization method with the improved conjugate gradient algorithm to identify the boundary condition of the heat transfer model.Based on the above reasons, we select the CGA to solve OCP of PDE, and develop an improved three-term spectrum conjugate gradient algorithm (ITSCGA).

Motivation and innovation
Based on the above analysis, few literature study the OCP of UHTMCT in continuous casting.Therefore, in the view of above valuable optimal control methods of unsteady PDE, this paper investigates the OCP of UHTMCT, and presents a parameterization optimal method that combines gradient analysis based on a Hamiltonian function costate system and an ITSCGA.Due to the rapidity for gradient, the control vectors of the OCP are discretized in time and space, and the expression of the gradient calculation method based on the Hamiltonian function costate system is analyzed.Then, an ITSCGA is proposed to solve the OCP of UHTMCT.Meanwhile, the global convergence of ITSCGA is demonstrated.The following are the main contributions of this study: An OCP of the UHTMCT is investigated.The gradient of OCP for UHTMCT is analyzed based on a Hamiltonian costate system by approximate treatment, matrix calculation, and composite trapezoidal integral method.An ITSCGA is proposed in this paper, and the global convergence of ITSCGA is analyzed.
The structure of this paper is as follows: in Section ''The solution of dynamic PDE and control vector discretization,'' the parameterization of the control vector and the solution of dynamic PDE are described; in Section ''The analysis of gradient calculation,'' we give the strategy of fast gradient calculation; in Section ''Improved three-term spectrum conjugate gradient algorithm and its global convergence,'' the ITSCGA is proposed, and the global convergence of the ITSCGA is analyzed; Section ''Experiment simulation'' verifies the performance of the ITSCGA by experimental simulations.

The solution of dynamic PDE and control vector discretization
The optimal control model for unsteady PDE with convective term is described as follows: where O T := f(x, t, z) 2 < 3 : x 2 ½0, l, t 2 ½0, t f , z 2 ½0, L, t f .0g, t represent time; x represent the width direction of the billet; z represent length direction of the billet; J(h) is a cost function; h(t, z) is the heat flux (W=m 2 ), which satisfies h(Á) 2 H(O T ); H(O T ) : = L 2 (O T ) is the set of admissible unknown heat fluxes of h(Á).kÁk is norm; the function f(Á) = kÁk 2 .The T subscripts x, z, and t stand for the partial derivatives with respect to x, z, and t, respectively; V cast is the casting speed (m=min); T is the temperature of steel billets (K); T c is calculated temperature of steel billets (K); T d is target temperature of steel billets (K); t f is the final time (s); k 0 is a thermal diffusivity coefficient (W=m Á K) that depends on the physical characteristics of the steel. 17a is a positive number; L and l are the length and width of the steel billet (m), respectively.T(x, 0, z) = T 0 is initial condition, T 0 is the initial temperature (K), which is a constant; The boundary conditions at positions x = 0 and x = l are as Àk 0 T x (x, t, z)j x = 0 = 0 and Àk 0 T x (x, t, z)j x = l = h(t, z); The boundary conditions at positions z = 0 and z = L are as Àk 0 T z (x, t, z)j z = 0 = 0 and Àk 0 T z (x, t, z)j z = L = 0.

The solution of dynamic PDEs
Equation ( 1) contains an UHTMCT that can be used to describe the solidification process of continuous casting.As shown in Figure 1, to solve this UHTMCT, this paper considers a life cycle model method (LCMM). 8n the LCMM, the first slice is generated when continuous casting begins.The first slice moves d distance along the direction of the casting speed, and the second one is generated.The slice is created at the surface of the mold, and dead at the outlet of the secondary cooling zone.Therefore, the whole billet is composed of many slicing elements, and the LCMM is used to solve the transient convective heat transfer model.The LCMM divides the whole billet into Q slices along the direction of casting speed, and the distance between two adjacent slices is d.The heat transfer equation of each slice is given as follows: We solve the equation ( 2) by using the finite difference method, which is shown in the following: where q is the indexes of space nodes in z; Dt is time step; Dz is spatial step size.n + 1 and n are the backward and forward time intervals, respectively.By organizing the equation (3), we can obtain where r = k 0 Dt Dx 2 .The boundary condition can be written as: where T w is cold water temperature.So we can get the upper boundary: and the lower boundary: The equations ( 2) to (7) are given to solve one slice model.Furthermore, the thermal convection term ∂T = ∂z can be replaced by the movement of the slice based on the LCMM.

Control vector discretization
According to equation (1), we can obtain the cost function of the OCP as follow: In the following, the specific discretized process of the control vector h(t, z) is given to obtain an approximate problem.Firstly, the control vector is discretized into N stages [t iÀ1 , t i ] within the time range of ½0, t f , and the space direction z is divided into M sections [z jÀ1 , z j ] respectively, which can be written as follows: Secondly, the piecewise constant x ½t iÀ1, t i ) is used to discretize the time horizon of control vector, where i = 1, 2, 3, Á Á Á , N and j = 1, 2, 3, Á Á Á , M. Thirdly, by using the matrix calculation, the control vector h(t, z) can be approximately expressed as, . . .
Generally, the control vector in a dynamic system satisfies the relation h min (t, z j )4h(t, z j )4h max (t, z j ) in the range of 0 to t f , where h min (t, z j ) and h max (t, z j ) are the given boundary limitations.Therefore, the limitation of So we can obtain  , x ½t 2 , t 3 ) . . .: According to Lin et al., 52 it has proved that the state trajectory can be calculated precisely.According to the equation ( 16), we can find that h(t, z) is an approximate function of s. Figure 2 shows the discretization of the control vector in time at the j space layer.The discrete stages N of time layer can affect the degree of curve approximation.Obviously, the higher the discrete stage, the better the approximate accuracy.
Therefore, the equation ( 1) can be approximately converted into an OCP as follows: Based on the Loxton and Lin, 53 it is easy to conclude that the formula (17) converges to formula (1), so the proof process is not detailed description here.

The analysis of gradient calculation
By analyzing the equation ( 17), it can be seen that s is the decision variable of the OCP.Generally, gradientbased solvers are used to solve the OCP due to its rapidity.Therefore, the gradient calculation for the cost function J in this paper is analyzed based on the Hamiltonian function costate system.The specific description is given as follows: The Hamiltonian function of the cost function is defined as follow: where .Meanwhile, l2R n is the adjoint vector or costate vector; The T in the upper right corner of a vector represents the transposition of the vector.Correspondingly, we can obtain the costate system 16 of the Hamiltonian function, Where l is the solution of costate system (19).
Consequently, the following theorem can be used to determine the gradient of the cost function.Theorem 3.1.If the Hamiltonian function is defined by equation ( 18), the costate system is generated by equation (19).The gradient formulate of J can be expressed as below: Proof.Assuming v(t) is an arbitrary continuous function.Then J in equation ( 17) can be rewritten as: Applied the integration by parts, So the following equation can be obtained: where T(x, 0, z) = 0.
Taking the derivative of the formula ( 23) with respect to s, we have Because v(t) is a continuous function, we choose arbitrarily and let v t ð Þ = l, we can obtain the gradient formula of Hamiltonian costate system, ∂J ∂s =a Therefore, Theorem 3.1 is proved.
In order to deal with the integral term in equation ( 25), an approximate treatment and composite trapezoidal integral method is used.So equation ( 25) can be approximated by the following equation: According to the compound gradient quadrature formula, we can obtain where In equation (27), ) is a remainder term, which can be easily calculated by the compound gradient quadrature formula. 16n where j is a constant that satisfies Meanwhile, the time is divided into N equal intervals, we have Due to the following inequality, we can obtain The remainder term R n ( H( Á )) can be expressed as Proper selection of H can effectively reduce its influence on the accuracy of gradient calculation. 16On the basis of this, the impact of fast gradient calculation on optimization stability is also Insignificant.

Improved three-term spectrum conjugate gradient algorithm and its global convergence
In some famous CGAs, the choice of conjugate parameters has different influences on the optimal control problem, such as Polak-Ribiere-Polyak (PRP), 54 Fletcher-Reeves (FR) method, 55 and Dai-Yuan (DY). 56ccording to Andrei, 57  often have better computational performance, but they not generally be convergent.Therefore, many hybrid CGAs and improved CGAs have been proposed.Significantly, Liu and Li 58 proposed a spectral CGA to solve unconstrained optimization problems; Yin et al. 59 proposed a hybrid three-term CGA by incorporating the adaptive line search based on hybrid technique; Li et al. 60 proposed two three-term spectral CGAs using Quasi-Newton equations to solve unconstrained optimization problems; Waziri et al. 61 proposed a double direction three-term spectral conjugate gradient method.
Therefore, CGA is a powerful technique for solving linear optimization problems, it is employed to deal with the OCP (17).However, the cost function in this paper contains control variable boundary constraints, which cannot be solved by restricted CGA directly.Therefore, this paper applies a trigonometric conversion function method to transform the original OCP into an unconstrained optimization problem.The temporary variable v(t, z) is defined and the transformed function is given by the following: h=0:53(h max (t,z j )Àh min (t,z j ))3½cos(v(t,z))+1+h min (t,z j ); Because of À 14cos(v(t, z))41, the bounded variable h(t, z) can substitute the unbounded variable v(t, z), we have Improved three-term spectrum conjugate gradient algorithm In this section, we present an ITSCGA based on DY algorithm.An appropriate spectrum m k is found to guarantee the descent direction and global convergence by an improved Quasi-Newton transversal equation. 60,62,63Starting from an initial guess v o 2 < 3 , this method generates an approximation sequence v k f g by the following iteration process: where k is iterative number; v k is the current iteration point, a k is the step-size in some line search and the search direction d k is expressed as: where rJ k is gradient of cost function; y k = rJ k À rJ kÀ1 .The conjugate parameter b k and new parameter u k can be determined as follows: where c is a constant, which satisfies c .0. According to Li et al., 60 we can obtain where B k is the approximation of the two-order Hessian matrix r 2 J v k ð Þ.Because of the equation ( 36), the equation ( 39) can be written as Letting Multiplying the equation ( 40) by s T k B k , we have Making use of the Quasi-Newton equation B k s k = y k , the above expression can be rewritten as: Therefore, we can obtain spectral parameter However, the optimization of OCP ( 17) by the ITSCGA is divergent.Therefore, we define m k as A suitable value b is needed to find to accelerate the convergence of the ITSCGA.By trial and error, we obtain 04b40:8, and m k 2 ½t 1 , t 2 , t 1 51, t 2 420.
The solving process of the ITSCGA for OCP (17) are as follows: Step 1: Set k = 0, select v o 2 R, termination criteria e s , maximum iteration steps k max ; Step 2: Use the LCMM to calculate the temperature T of the billet; Step 3: Use equations ( 27) and ( 36) to calculate the gradient rJ(v k ) and the search direction d k respectively; Step 4: Calculate spectrum m k according to equation (44) Step 5: Update the variable v k according to equation (35); Step 6: Use equation ( 34) to calculate the cost function J(v k ), if J(v k )4e s , or the number of iterations reaches k max ; then stop; else go to step 7; Step 7: Set k = k + 1, go to step 2.
The above steps are applicable to the optimization process when the convective term ∂T=∂z = 0.The ITSCGA based optimizer is difficult to apply to the case of ∂T=∂z 6 ¼ 0. Therefore, when the convective term ∂T=∂z 6 ¼ 0, dynamic optimization process is introduced.The whole dynamic process is divided into several static times, and heat flow at each time is obtained by the ITSCGA, which is taken as the initial heat flow at the next time.Therefore, dynamic optimization is obtained through iterative calculation.

The global convergence of improved three-term spectrum conjugate gradient algorithm
In this section, global convergence of the ITSCGA under the Wolfe search criteria 16 : where 04d4s41.Theorem 4.1.If sequence v k and d k are generated by equations ( 35) and ( 36), we have the following equation Proof.Induction can be used to demonstrate the conclusion.Since rJ T 0 d 0 =2rJ 0 2 , the equation ( 47) holds for k = 0.For k51, we suppose that Theorem 4.1 is true for the last iteration k À 1, and we can obtain rJ T kÀ1 d kÀ1 \ 0. Multiplying the equation ( 36) by rJ T k and we have Case 1: If m k .1, because of 04b40:8, we have and According to equations ( 48) and ( 50), we have Due to 1 À m k \ 0, rJ T kÀ1 d kÀ1 \ 0, the sufficient descent property holds for c .0 in this case.
Case 2: If 0:8 \ m k \ 1, we have m k = 1.We also obtain the equation ( 50), so we have Due to rJ T kÀ1 d kÀ1 \ 0 and c .0 in this case.Case 3: If m k \ 0:8, we have m k = 1.We have so we can obtain the following inequality: Due to rJ T k d kÀ1 \ 0, and the sufficient descent property holds for c .0 in this case.
According to equations ( 51), ( 52) and ( 55), we can obtain rJ T k d k 4 À rJ k 2 .Theorem 4.1 is proved.Remark: In our experiments, we conduct extensive testing and debugging on parameter c and find that the optimal value range for parameter c can be chosen as 14c47.The parameter in this range makes the algorithm have better performance.
Theorem 4.2.Based on the strong Wolfe line search, 16 we have the following result, Proof.We suppose that equation ( 56) is not true, so there exist constants g 1 , g 2 , which satisfy g 1 \ rJ k 2 \ g 2 .According to equation (36), we have So we have where m k 2 ½t 1 , t 2 .Multiplying the inequality (58) by Furthermore, we can obtain Based on Wolfe line search criteria and Lipschitz continuous, 19 we have and Where L t is the Lipschitz constant.According to equations ( 35) and ( 45), we can obtain the following inequality where C =2d sÀ1 L .Hence we can obtain the following equation So we have Obviously, equation ( 61) is contradict to equation (66), so we can obtain lim k!' inf rJ k = 0. Theorem 4.2 is proved.

Experiment simulation
We carry out the static and dynamic simulation experiments in this section.In the static situation, the simulation experiments are divided into three situations.In case 1, the parameters c and b are analyzed in the ITSCGA; In case 2, the comparison of the ITSCGA with other algorithms is given; In case 3, the control effects are tested and verified.In the dynamic situation, two continuous casting experiments with variable casting speeds are tested, and the optimized control method based on ITSCGA is compared with other algorithms.

Static experiment simulation
This section considers the convective term ∂T=∂z = 0 in formula (17), and the values of process parameters are given in Table 1.Table 2 shows specifications of the continuous caster.In addition, eight control points are select.The distance (m) between each control point and the meniscus are 0.1, 1.0, 1.5, 2.0, 2.6, 3.3, 4.0, and 4.5, respectively.We set the target temperature (8C) values as 955, 644, 599, 567, 555, 547, 550, and 553, and select the maximum number of iterations k max = 100.
Case 1: In this case, we specifically analyze the values of parameters b and c in equations ( 38) and (44), where 04b40:8, 14c47.We take different values of the parameters b and c within this limit, and compare the cost function in the simulation experiments.We obtain the values of parameter b and c in Table 3.The value of parameter h min (t, z j ) and h max (t, z j ) are 20 and 2000 respectively.According to the equation (1), we can calculate cost function J(h), and Figure 3 gives the results of the ITSCGA for different parameters in Table 3.It can be seen from Figure 3 that the convergence speed of the ITSCGA gives a better result when the parameters are chosen as c = 7 and b = 0:0.
Case 2: In this case, we take the simulation for comparison, in which the ITSCGA is compared with the CGA, 55 the modified conjugate gradient algorithm (MCGA), 50 the new conjugate gradient algorithms (NCGA-JCJ2 and NCGA-JCJ + ), 60 and the hybrid conjugate gradient algorithm (HCGA). 58The convergence of the cost function and the number of iterations are shown in Figure 4.In addition, the termination criteria are selected e s = 1:0.It can be seen from Figure 4 that ITSCGA significantly accelerates the convergence speed.Table 4 provides the number of iterations and runtime of the six algorithms.By comparing the results in Table 4, it can be seen that ITSCGA only iterates 39 times and reaches the termination criterion.Compared with other methods, ITSCGA has fewer iterations and faster convergence speed.
Case 3: The control effects with six algorithms, including CGA, 55 MCGA, 50 NCG-JCJ + , 60 NCG-JCJ2, 60 HCGA, 58 and ITSCGA, are compared.We set the target temperatures as 955°C, 644°C, 599°C, 567°C, 555°C, 547°C, 550°C, and 5538C.Figure 5(a) and (b) show the control effects of the six methods.According to Figure 5(a) and (b), we can see that ITSCGA makes the temperature reach the target value within 38 s, which is smaller than other algorithms.So the control effect of ITSCGA is better than that of other five methods.Meanwhile, to compare the convergence effect more clearly, we calculate the control error for the six methods.The error curve is shown in Figure 6(a) and (b), and the averaged errors are given in Table 5 (T a is the averaged error temperature between the control and target temperature).According to the calculation results, the error of the optimizer based on ITSCGA is smaller than that of other methods.We can see from Figures 5(a

Unsteady experiment simulation
In this section, we carry out unsteady experimental simulation on the OCP.We choose c = 1, b = 0:8, and N = 60 in formula (9) and M = 40 in formula (10).In this situation, the convective term ∂T=∂z 6 ¼ 0 is considered in formula (17).T c is obtained by calculating the average value of the temperature of each cooling section.The initial water flow rates of sections 1-8 are set as 357, 331, 406, 286, 328, 290, 236, and 89 L/min, respectively.The target temperatures are 1055°C, 1035°C, 1027°C, 1017°C, 995°C, 975°C, 965°C, and 9458C.The specifications of the caster and the values of technical parameters are given in Tables 2 and 6 respectively.Case 1: The casting speed is changed from 0.8 to 1.1 m/min.We compare the optimization performance of four algorithms, including the ITSCGA, CGA, 55 MCGA, 50 and HCGA. 58Figure 7 shows the change curve of the cost function J based on four optimizers.Obviously, the convergence speed of the cost function of ITSCGA is faster than that of other methods.Then, we give the temperature curves of eight cooling zones, as shown in Figure 8(a) and (b), and the errors between the control temperature and target temperature are given in Figure 9(a) and (b).In comparison to other algorithms, the ITSCGA gives smaller temperature fluctuations, so its control effect is better than that of the other three methods.Meanwhile, to compare the convergence effect more clearly, we calculate the average error and maximum error for the four methods, which are given in Tables 7 and 8. From these two tables, we can see that the average and maximum errors of ITSCGA are smaller than other methods.
Case 2: The casting speed is changed from 1.0 to 1.6 m/min.In this case, we consider the impact of significant changes in casting speed on control performance of the four algorithms.The control effect of the simulation experiment is shown in Figure 10(a) and (b).It can be seen from Figure 10(a) and (b) that ITSCGA has faster convergence speed and more stable temperature fluctuation than other methods.Figure 11 gives the average error for the four algorithms.From Figure 11, the average error of the ITSCGA is smaller, and the average error values of all eight cooling zones are less than 48C.
Obviously, especially in the case of variable casting speed, the dynamic optimization based on ITSCGA gives smaller temperature fluctuations and obtains a better control effect compared with other methods.

Conclusion
This paper investigates the OCP based on the UHTMCT, and presents a parameterization optimal method that combines gradient analysis based on a Hamiltonian function costate system and an ITSCGA.Firstly, we analyze the expression of the gradient calculation method based on the Hamiltonian function costate system to explore the rapidity of the gradient.Secondly, an ITSCGA is proposed to solve the OCP of the UHTMCT, and the global convergence of the ITSCGA is proved.Lastly, both static and dynamic continuous casting simulation comparisons are tested Table 7.The average values of error in case 1 of dynamic experiment.

Cooling section Algorithms
ITSCGA MCGA 50 CGA 55 HCGA 58 1 T a = 0:034 T a = 2:58 the optimization performance of ITSCGA, and the results show that the ITSCGA is better than other methods.Specifically, in the static experiment, the ITSCGA has fewer iterative numbers and a faster convergence rate.Meanwhile, the optimizer based on the ITSCGA provides more stable temperature control, and the average error of the eight cooling zones is smaller than other methods.
In the dynamic experiment, we consider the simulation of continuous casting experiments under two different cases of variable casting speed to verify the stability of temperature control.When casting speed changes from 0.8 to 1.1 m/min, the ITSCGA has a better control effect, and its average error is less than 0.58C, which is smaller than that of other methods.When casting speed changes from 1.0 to 1.6 m/min, the average error of the ITSCGA is less than 48C.Therefore, we conclude that the optimizer based on the ITSCGA can converge faster and have more stable temperature fluctuations than other methods, which can help to improve the quality of the slab.
In this paper, the optimization of secondary cooling water control under the change of casting speed is studied.In future research, optimization of secondary cooling water flow rate based on the multi-physical coupling,  as well as modeling and optimization of crack index model for a slab will be considered.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Figure 1 .
Figure 1.The solidification process and life cycle model method of the cast billet.
properties, but they may have modest practical performance due to jamming.The b PRP k may ) and (b), 6(a) and (b), the optimized control based on ITSCGA gives a more stable control temperature.

Figure 4 .
Figure 4. Comparison of computational performance of algorithms.

Figure 3 .
Figure 3.The convergence of the cost function with the number of iterations.

Figure 5 .
Figure 5. (a) Control effects comparison of control points 1 to 4 in static experimental simulation.(b) Control effects comparison of control points 5 to 8 in static experimental simulation.

Figure 6 .
Figure 6.(a) Comparison of error variation curves of control points 1 to 4 in static experimental simulation.(b) Comparison of error variation curves of control points 5-8 in static experimental simulation.

Figure 8 .
Figure 8.(a) Control effect of cooling sections 1 to 4 in case 1 of dynamic experiment.(b) Control effect of cooling sections 5-8 in case 1 of dynamic experiment.

Figure 9 .
Figure 9. (a) Errors comparison of cooling sections 1 to 4 in case 1 of dynamic experiment.(b) Errors comparison of cooling sections 5-8 in case 1 of dynamic experiment.

Figure 10 .
Figure 10.(a) Control effect of cooling sections 1 to 4 in case 2 of dynamic experiment.(b) Control effect of cooling sections 5 to 8 in case 2 of dynamic experiment.

Table 1 .
Process parameters of steel billets in the static experiment.

Table 3 .
The values of parameters b and c.

Table 4 .
The comparison of algorithms.

Table 5 .
Comparison of the averaged error of the six methods.