1. Introduction
Induction heating is a complex multiphysics problem where electromagnetic, thermal, mechanical, and metallurgical problems are influencing each other. Therefore, modeling and simulation require the simultaneous consideration of complex interactions. The knowledge of the process variables and their complicated relationships is fundamental and necessary for simulation and modeling but not fully possible. Therefore, the process optimization and process control of induction heating and hardening is a true challenge both for the experimentalists and the theorists in the field.
There is a significant amount of literature and progressive study on the simulation and modeling of induction heating [
1]. Numerical approaches such as the finite element (FE) and finite difference (FD) are quite popular for the analysis and modeling of coupled problems like induction heating. Since these approaches are based on discretization, the resulting differential equation system is very large, in the order of millions. Control approaches based on these methods, even if possible in principle, are challenging because they are computationally expensive and are not applicable when quick decision making is necessary. For example, the authors in Ref. [
2] have studied and examined the efficiency of optimization approaches for induction heating processes based on the FE approach. They have focused on the power density distribution in the sample and the goal of achieving a prescribed temperature at a specific location near the surface of a cylindrical magnetic sample. Such models have the advantage of (1) a general description of physical mechanisms and (2) validity over a whole domain such as a sample or workpiece, as shown by [
2] and other authors [
3]. However, in face of the complexity of the algorithm presented in [
2] and the usage of computationally expensive FE tools to solve electromagnetically-thermally coupled problems, the agreement between the optimized and targeted temperatures (Figure 12 of Ref. [
2]) might not be sufficient for process optimization and especially process control concepts. Moreover, it is often necessary to consider the transient evolution of process parameters such as the current and frequency over the course of a heating cycle instead of assumed constant values. Model reduction is also an active research area for the control design and process control. In Ref. [
4] the authors apply a reduced FE model for the control of the temperature profile of an axisymmetric sample by finding an optimal control input current. However, the comparison between the full and reduced model results shows an 80 K difference in the temperature at intermediate times. Therefore, applying a reduced model in real time and on induction heating test rigs to control the temperature profile will not perform satisfactorily.
In general, the mere application of physical approaches is not applicable to the induction heating systems because of the complexity of the power systems and the implemented control algorithms to the process control inputs. In real applications the induction heating operator is just allowed to change the power supply input control set values, while the control system, which is usually not accessible to the operator, will set the output physical parameters. To bypass this complexity and to propose a general approach applicable to any induction heating set up, and to avoid modeling of the power generator circuit, we present and apply an inverse model to control and optimize an inductive heat treatment process. The modeling strategy is not meant to replace a fully physical approach, but as a powerful alternative that is better suited for performance driven applications including interactions with the system control. With the implementation of such model types in the system control, the digital process control of induction heating may become a feasible task in the near future.
By combining a simplistic mathematical model with controlled experiments, we take into account the unknown relevant factors, (which are not explicitly included in the model, e.g., nonlinear boundary conditions, geometry imperfections, power generator complexities, etc.) that are important for the temperature evolution. In our induction heating test rig a cylindrical ferromagnetic sample is heated up above its Curie temperature in static experiments, where the inductor and sample are fixed. A mathematical model of the coupled electromagnetic-thermal problem is then calibrated with experiments to extract the temperature dependent model parameter
. The parameter
corrects the so-called modified Nagaoka coefficient which has mainly been assumed constant in the literature [
5,
6]. In principle, this parameter is temperature dependent and it changes during the heating process. However, its explicit functionality is not well known apart from some conjectures [
7]. Therefore, we allow the optimization of this parameter by calibrating the model to the measurement data. In this way, the correct magnetic field in the air gap between the inductor and the sample is approximated and the modeled temperature is closer to the experiment. We also measure the current in the inductor and extract its frequency from the signal. Recording of the power level, current, frequency, temperature, and the model parameter
provides the basic ingredients to define an inverse model that enables us to predict and conduct experiments with specific goals. In industrial setups, due to instrumentation limitations, it is usually not feasible to continuously control the temperature of the samples applied to the induction heating processes. Moreover, the knowledge of the process control set values is necessary to target specific temperatures in specific times. The empirical way of trial and error is expensive, and the FEM approaches are computationally demanding or might even be impossible in industrial scales. To demonstrate the applicability of the presented approach for optimal control of induction heating processes we show two examples, where the target temperatures are chosen randomly to be achieved in specific times. Employing the inverse model we show how quickly and efficiently the target temperatures are reached by providing calculated control input set values to the power supply of our test rig.
The paper is organized as follows: In
Section 2 we present the mathematical model and the solution of the coupled electromagnetic-thermal problem under certain assumptions.
Section 3 shows the test rig setup for the experiments and in
Section 4 the material properties of the heated sample are presented. In
Section 5 we calibrate the forward model with experimental data and present some test simulations. In
Section 6 we present an inverse model and as an application in
Section 7 we demonstrate two experiments in which we desire to hit the prescribed target temperatures in specified times. Finally, in
Section 8 we summarize and make a conclusion.
3. Induction Heating Test Rig Setup
The experimental setup is shown in
Figure 3. The power supply operates both in the temperature control and power control mode with the capacity to provide a maximum power of 30 kW. The graphical user interface provides control over the operating power, heating time, movement of the inductor, rotation of the sample, and quenching. Both the water and air cooling are possible in the unit. Alternatively, the test rig can be controlled by means of the input files, which are loaded into the system control. In this way, heat treatment procedures with a high number of individual steps can be carried out.
The control system of the test rig provides the data of relevant process parameters during heating. These include the resonance frequency, running power, current, and voltage (inside the generator before the transformer, not in the inductor), temperature of the cooling water at the inlet and outlet of the solenoid, and the temperature of three points in the sample (two with thermocouples and one with a pyrometer or alternatively all three with thermocouples (NiCr-Ni, type K)). The recording of the data is performed with a sampling rate of 10 ms.
The cylinder of length mm and diameter 22 mm is fixed between two clamps. The rectangular inductor of length = 37 mm and internal diameter mm is used to heat up the part of the sample beneath the solenoid. The setup is static where neither the sample rotates nor the inductor moves. The point where the temperature is measured is located at a radial distance of mm from the surface and the thermocouple hole has a depth of 65 mm from the top side of the cylindrical sample. The position of the inductor is fixed such that the temperature measurement point (at the end of the thermocouple hole) is located at the center of the inductor. The diameter of the thermocouple measures 2.0 mm, while the hole has a diameter of 2.2 mm.
The current in the inductor is measured by implementing a Rogowski coil right after the transformer as shown in
Figure 3. A Rogowski coil is wrapped around the straight conductor, which carries the current to the inductor. The induced voltage in a Rogowski coil is proportional to the rate of the current in the straight conductor. Therefore, the measured voltage is integrated to get the current as:
where
µH is the mutual induction of the Rogowski coil and
is the voltage measured by an Oscilloscope attached to our test rig (to measure the voltage induced in the Rogowski coil). In principle, Rogowski coils are immune to external fields and show very low positioning related sensitivity. During our experiments this was checked and two different Rogowski coils have been used. The finally calculated currents from both Rogowski coils differed by a few amperes only. The frequency of the circuit is extracted from the time series of the current.
5. Model Calibration and Testing
In this section, we check the accuracy of the mathematical model used for describing the experiments. To this end, we allow the optimization of only one parameter in the model. As described in
Section 2, the optimized parameter corrects the magnetic field in the air gap between the cylindric sample and the inductor and even takes into account other effects that are not included in the model. For example, these are nonlinearity effects, inaccuracy of the material data, geometry asymmetries, and the effects originating from the test rig setup.
We used a collection of experiments to calibrate the model.
Table 2 shows the list of power variation experiments for the calculation of the model parameter
. The minimum and maximum powers in the measurements are 0.6 and 15 kW. We ran every experiment until the Curie temperature was exceeded. The repetitively used sample for the experiments was assured to be free of any cracks or other damages. To keep the sample in a good condition, we first conducted the experiments at low powers and water quenching was only applied after convective cooling of the sample in ambient air down to 400 °C. In this way, the sample does not experience extreme thermal and mechanical stresses.
In the calibration stage we construct the tables of power, temperature, current, frequency, and the correction factor
. The only input control parameter in the experiments is the set value for the power in the control system. The current and frequency are determined by internal algorithms in the power supply, which depend on the sample temperature as well as the geometries of the sample and inductor. In the model, the current and frequency are taken from the Oscilloscope measurements. To find
in Equation (15), we use a fixed integration time step
= 0.01 s that is equal to the sampling rate of the temperature measurement. The knowledge of
for several powers and its corresponding currents and frequencies provides the basis to define an inverse model (see
Section 6) with the ability to design experiments with specific goals, such as the predefined temperature evolution.
Figure 5 illustrates the results of the model calibration for the power experiments with 13.5 and 2.1 kW and their corresponding optimized parameter
as a function of the temperature. We remind that the temperature is monitored at the measurement point of 2.1 mm below the surface. In principle, one could do this for any point in the sample. As we can see, with just one parameter the measurement data are accurately reproducible. With the increasing temperature, the evolution of
shows a transition zone corresponding to the Curie point of the material. For high powers, such as 13.5 kW, we see that
is almost constant while for low powers, such as 2.1 kW, it is decreasing while approaching the Curie temperature.
After calibration, we test the model on some experiments.
Figure 6 shows the accuracy of the calibrated model for four measurements different from the calibration set. In the upper panels of
Figure 6 we show the results for 12 and 15 kW of the power experiments using the model calibrated with the 13.5 kW measurement data. As we can see, the simulated temperatures are very close to the experimental data. In the lower panels of
Figure 6 the results are shown for 1.5 and 3 kW of the experiments using the model calibrated with the 2.1 kW measurement. Additionally, for the low power regime, the agreement is very promising. These examples show that using the
optimized for a specific power can provide a very accurate predictive model for experiments close to that power. This reflects the fact that
α does not change substantially for the neighbouring power experiments. Therefore, in practice, for the model calibration a discrete set of experiments is required. Then, interpolation of the data will provide a database accurate enough for modeling.
6. Inverse Model
In this section, we present the inverse model and its application for control of the induction heating systems. The temperature control of the induction heating systems is a daunting task and up to now not completely managable. Direct temperature control is not possible in the industrial context since thermocouples can-not be applied in continuous processes and thermographic methods are not sufficiently reliable. A prior knowledge of the control set values to reach certain target temperatures in the sample in specific times has a great value for industrial setups and applications. Furthermore, optimal temperature control paves the way to hit desired targets in material properties as a final goal. The empirical way of the temperature adjustment is based on trial and error, which is not cost effective. Moreover, methods based on FEM simulations are computationally expensive and might be even impossible on industrial scales.
In our test rig power controlled static experiments, the only inputs to the control system are the power level and the total holding time. Therefore, the temperature of the sample at time
t is a complicated function of power
P, which can be expressed as
The inverse problem is the process of finding power levels, which lead to certain temperatures in the sample, i.e.,
It is impossible to derive the explicit mathematical form of
. However, employing a mathematical model, which is incorporated with some physical parameters, could assist to reverse engineer the problem. As described in
Section 5, the unknown model parameters,
in this work, are determined by calibrating the model with the experiments. Performing the experiments with different power levels and recording the model parameters (current, frequency,
), the simulated temperature
and input set value
, provides a tool to bypass the necessity of knowing the explicit form of
. With a discrete set of experiments and the interpolation of the data, we can infer the required input powers to reach specific temperatures. The linear interpolation (by using a standard Scilab built-in function [
8]) has been done for three surfaces (
, current, frequency) of a matrix with the power level and temperature axes. The grid spacing of
kW and
°C is used in the interpolations.
To calculate the power levels with the model we start with an initial power level of = 0.6 kW to simulate the temperature for a total time of by using Equation (15), where is the index for numbering the target points ( = 1, 2, 3, ...). We notice that ) and for multiple target temperatures the initial temperature in the simulations with Equation (15) for a target point is .
The model parameters of the current, frequency and
are read from the interpolated surface data, based on the current values of the temperature and power level. The power level is changed in an iterative way to reach the target temperature, i.e.,
where
is the
-th target temperature and
is the simulated temperature after a total time
. The iteration is terminated when the target temperature is simulated within a tolerance of
.
7. Inverse Model Application: Demonstration of the Temperature Control
In principle, if the model is accurate enough, the initial powers calculated by the inverse model for reaching certain prescribed temperatures in the experiment might be good enough. However, in practice extra iteration loops with a few experiments might be required to reduce the target temperature error. In
Figure 7 we show one iteration of this scenario where
and
are the target temperatures and
and
are the calculated initial powers by the inverse model as described in
Section 6. When the model is imperfect, we reach other temperatures
and
(see
Figure 7) in the experiment, which might be slightly higher (
> ) or lower (
< ) than the target temperatures. To get closer to the target temperatures, in the next experiment, a new set of input power levels is calculated by using the temperature values realized in the current experiment as target temperatures for the inverse model, i.e.,
and
. Then, we adjust the power levels for the next experiment depending on the target temperature error sign (
> 0 or < 0). We summarize the steps and the procedure for the update of the power levels as follows:
- (1)
Initial suggested powers (
and
) found based on the inverse model to hit the target temperatures, i.e.,
and
in
Figure 7.
- (2)
Run the experiments with the suggested powers.
- (3)
Compare the temperatures ( and ) with the target temperatures. Terminate the loop if the target is achieved within the specified accuracy, || < tolerance.
- (4)
If the target is not achieved, adjust the power levels. For the adjustment, the model calculates (in an iterative way as described in
Section 6) the power levels to reach the new temperatures (
and
) in step (3) and then suggests a set of power levels for the next iteration based on the difference between the new power levels and the old ones used in step (2),
. The
is 0.03 or −0.03 depending on whether the measured temperature (
) is below or above the target temperature.
- (5)
Go to step (2).
As an application for the inverse model, we show the accuracy of the approach for designing targeted experiments in the induction heating systems. For the demonstration, we show two control experiments where we target two and three temperatures at specific times. In
Figure 8 we show the results.
Figure 8a shows the experiments to hit the temperatures
= {500, 900} °C at times
= {30, 80} s. In
Figure 9 on the left panels the corresponding power levels (inferred by the inverse model) and errors in the target temperature are shown for each experiment. As we can see in three iterations the target temperatures are achieved within a tolerance of |
| < 12.
Figure 8b shows the results for the second control experiment with iterations to hit the temperatures
= {700, 825, 1030} °C times
= {11, 34, 69} s. In addition, this control experiment converged in three iterations and we have omitted the iteration number two in the plot for better readability. As we can see in
Figure 9 the power levels are not changing significantly during the iterations and already the first suggested powers lead to promising results. However, for better optimization and reducing the error (see
Figure 9, lower panels) in the target temperatures a few iterations might be necessary. The control and design of experiments similar to the one in
Figure 8b is very appealing to material scientists. In inductive heat treatments, the sample is heated up in a short time to achieve temperatures close to the Curie point and then the heating continues at a lower rate. A preliminary setup of such experiments and a good guess of the control input set values are crucial. The inverse model can substantially reduce time and energy for the design of such experiments.
Temperature control with smaller convergence tolerances is possible at the cost of increased experiment numbers. For higher accuracy and smaller tolerances each iteration should be done with a new sample. Since material properties are not changing significantly, the repeated use of a single sample is sufficient in order to show the principles of the presented method.
The optimization problem in this work may be a nonconvex one. In the application of the inverse model to find a control input power level, we always started the iterations with the lowest power, i.e., 0.6 kW. In this way, the iterative procedure most of the time converged to a global optimum. However, more advanced and general optimization approaches like the genetic algorithms [
11] are preferable and might be necessary when the parameter search space is larger (more control input variables) or more complicated optimization constraints are imposed.
8. Summary
The inverse model is used to conduct and control experiments which are very interesting for heat treatment applications. We performed two experiments with temperature targets at specific times during the course of heating. We demonstrated that by using the inverse model the target temperatures are hit within a few iterations. This is a significant improvement in comparison to the classic way of trial and error and even complicated mathematical optimization algorithms [
2,
12]. The method was only applied to a position close to the surface. However, it can easily be extended to multiple locations within the sample.
Furthermore, the presented approach can be applied to a broad diversity of non-ideal complicated scenarios. For example, to the induction heating setups where the sample is either asymmetric (and has a complex geometry) or has non-ideal characteristics, or industrial configurations where an array of coils are combined to surface harden the samples. Most often the FEM approaches fail to simulate such large-scale processes due to the complexity of the boundary conditions and multiple uncontrolled process variables.
In this work, the inverse model and the temperature control were demonstrated for static cases without the quenching or cooling steps. In principle, the inverse model can be extended to include quenching and cooling operations and even to take into account mechanical and metallurgical problems. Last but not least, the approach is computationally very efficient and that makes it an attractive alternative to control large scale induction heating systems.
The computational burden of our approach is in the order of min (1 to 10 min) and it grows quickly when the number of target temperatures is increased. The response surface models are quite useful for use in the system parameter identification and fast process control. These surrogate models [
13] are applicable with a low computational effort to both linear and nonlinear problems. It is in the interest of authors to develop and investigate the neural network based metamodels (fast running surrogate models) to control and optimize the induction heat treatments.