eduzhai > Physical Sciences > Materials Sciences >

Goemetrique induction cooking model with increasing temperature and its optimization

  • sky
  • (0) Download
  • 20211030
  • Save International Journal of Materials Engineering. 2011; 1(1): 17-20 DOI: 10.5923/j.ijme.20110101.02 Modeling and Optimization of a Goemetrique Induction Cooking for Improving the Temperature M. Rezig*, K. Srairi, L. Aloui Department of Electrical Engineering, University of Biskra, Algeria Abstract This paper presents an optimal geometric design of a method for induction heating system. The objective of this design is to obtain a desired optimal distribution temperature in the load. The dimensions of the inductors are taken as design parameters and are replaced at each iterative step, using the bisection method to minimize the error of the calculated temperature and the desired one in each node. The study of such a device requires the use of models or modeling, physical, mathematical and numerical. Representation this modeling is the basis of the understanding, design and optimization of these systems. The optimization technique is to find values of variables that maximize or minimize the objective function. Keywords Numerical Methods, Induction Heating Systems, Geometrical Optimization, Electromagnetic and Thermal Coupling 1. Introduction Through the researches achieved in the last twenty years, many industrial processes requiring heating operations make use of induction heating devices. Many industrial applications have been developed, surface heat treatment, forging, bonding by polymerization, domestic and household applications. Read phonetically However, the development of heating equipment induction is not easy. The inductors are still often dimensioned empirically, based on the experience gained by the manufacturer, and the power profile to be injected into the part (pan), obtained by trial and error. However, the tests are long and costly. Also computer-aided design, based on numerical simulation of magnetic and thermal phenomena involved, has certain advantages in terms of cost, time and optimization for the development of induction heating processes. The prediction of magneto-thermal behavior of the whole work space - inductor allows the optimization of both electrical and mechanical design of the inductor, and also the necessary power required to treat the work space. The temperature distribution is itself a key factor for predicting the effectiveness of heat treatment; physical properties of materials are strongly temperature dependent. For nonlinear problems with complex geometry, the finite volume method developed in recent years, considered as an * Corresponding author: (M. Rezig) Published online at Copyright © 2011 Scientific & Academic Publishing. All Rights Reserved important tool for modeling many branches of engineering, in fluid mechanics, heat transfer, electromagnetism etc. The mathematical analysis of the finite volume method has also helped engineers to develop the basic principles that make it an efficient tool of discretization. Numerical modeling of problems involving multi-physics coupling is now a major area of research. A good comprehension of the phenomena requires passing through an important phase which is an efficient physical modeling and numerical problems. The design optimization of induction heating process is classified within this type of approach. The methods of heating induction is by nature a good example of coupling multi-physics (electromagnetism, thermal-mechanical and metallurgical). 2. Mathematical Modelling The physical modeling of induction heating involves both magneto dynamic equations and heat equations, this first stage based mathematical formulation is useful to determine the magnetic field created by the inductor in the surrounding space, to deduce the induced currents and local power dissipated by Joule effect in the conducting part to be heated, then in the second phase the dissipated power considered as a heat source to begin the calculation the temperature. Listen Read phonetically The magneto dynamic equations are derived from Maxwell's equations and constitutive laws of materials. The equation of heat diffusion is governed by the Fourier law. These equations are written in the case of the AV formulation related to the magneto dynamic problem. A mathematical model of induction heating incorporating 18 M. Rezig et al.: Modeling and Optimization of a Goemetrique Induction Cooking for Improving the Temperature coupled multi-physics must assure the coupling of the following equations: Maxwell's equations — to access to the currents induced in the specified part; The equation of heat — to model the thermal evolution due to the power dissipated by the currents induced in the part (piece). 2.1. Electromagnetic Model The electromagnetic problem is modelled using the four Maxwell's equations (with the approximation of the quasipermanent modes), with the relation between the magnetic field induction B and the magnetic field H, and the tradi- tional law Maxwell's equations of Ohm; for an axisymmetric configuration, a judicious combination of the preceding equations leads, recall that the magneto-dynamic equation in the case or the load is motionless compared to the inductor → →→ →→ rot(ν rot A ) +σ ∂A ∂ t = J ex (1) 2.2. Thermal Model The change of the temperature is primarily controlled by the power dissipated by the induced currents. The equation of heat is thus written below : − k ∇2T + q = 0 (2) Where the source term comes from the powers dissipated by joule effect. Where f(x) = e Such as The function n ∑ e = (Ti − Td )2 (4) i =1 The suggested method is applied has induction heating with pan like numerical example to order the optimal temperature of the device, the density of current the positions of the inductor and are snuff like parameters of design and the latter will be replaced in each stages of iteration. 3.2. Optimization Algorithm The optimization problem consists in the minimization of the temperature error "e" associated to each node to satisfy a (desired) target temperature. The steps of the mechanism search are summarized as follows: Initialization of the data Choice of the initial interval of the bisection method (dimensions of inductor 2) Resolution of the equations of the model related to the electromagnetic thermal coupling Calculation of the average temperature, as well as the error on the temperature: e(T) The optimal geometry of the inductor is found after criterion convergence is satisfied, otherwise repeat step2 4. Application 3. Optimization Problem Formulation Optimization is often summarized with the mathematical techniques of resolution to which one charges the failures thereafter that one meets. However, as for the majority of the problems which one is brought to solve, optimization must be the systematic step object which involves four fundamental phases summarized in Figure 1. The principle of Cooking induction based on two inductors placed under shape of coils with various geometries (figure 2), below a plate. The inductors act as a magnet producing an electromagnetic field concentrated on the part to be heated. Safety and economy are two important advantages related to cooking by induction; the studied system comprises: Optimization Prob- lem Mathematical Model Steps of Resolution Problem Problem Evaluation Figure 1. Fundamental Phases Summarized 3.1. Problem Presentation It is an optimal design of our device (cooking with induc- tion), which takes account of the experimental model. The analysis into the sensitivity of the geometry of the model by the bisection method is applied like an algorithm of optimi- zation to solve the problems of the device with induction. The optimization problem can be formulated as the fol- lowing: Minimize f (x) = f (x1, x2, …….xn), (3) Figure 2. Model design International Journal of Materials Engineering. 2011; 1(1): 17-20 19 Part Pan Table 1. the initial shape of the Coil Characteristics Inductor1 Inductor2 Inductor3 Thickness (Exciting coils) Thickness (Pan-Exciting coils) Distance A Distance B Distance C Distance D Symbols d1 d2 d3 ee eg A B C D Dimensions 0.0175 m 0.0073 m 0.0126 m 0.0004 m 0.0006 m 0.0100 m 0.0208 m 0.0218 m 0.0005 m Part Exciting Coils Table 2. the initial shape of the pan Characteristics Pan Thickness Symbols d/2 ei Dimensions 0.0950 m 0.0012m x 10-4 module potential magnetic A(A/m2) 2 1 0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 rP(mm) Figure 3. Variation of the Module potential vector A x 10-4 3 module potential magnetic A(A/m2) 2.5 x 10-4 4 2 2 1.5 070 60 50 1 40 40 30 30 0.5 20 20 10 10 r(mm) 0 z(mm) Figure 4. Variation of the module Potential vector in the plan (r,z) 4.1. Before Optimization The device of figure 2 presents the three aligned inductors; working on a relatively high frequency: F = 30000Hz. This system is fed by a current of amplitude I = 15 A with a current density of 5.106A/m2 and at a temperature 300 K0 Figure 3: presents the radial variation of the potential module magnetic vector ‘A’. Figure 4: presents the distribution of the potential module magnetic ‘A’ in the plan (r,z). Figure 5: presents the distribution of the temperature in the plan (r,z). Temperature (K°) 312 315 310 310 308 305 306 300 295 80 60 40 20 r(mm) 00 304 40 302 30 20 10 z(mm) Figure 5. Distribution of the temperature in the plan (r,z) 4.2 Post- optimizing In this study; we want to reach a temperature of 380 K0 in some point of thepart (piece) to be heated, by applying our method to the second optimal inductor. The working conditions are the same as those indicated above in term of frequency and operating current. Figure 6: presents the radial variation of the magnitude of the potential magnetic vector A. Figure 7: shows the distribution of the module of potential magnetic ‘A’ in plan (r,z). Figure 8: presents the distribution of the temperature in the plan (r,z). For a temperature of 300 K0; the term of potential vector A, density of the induced current and the density of the induced power are carried out, from these results are then extracted from the total power average, evaluated on all the part to heat, thus the voluminal density of average power PS to be used as source with the thermal problem If one wishes to increase the temperature in the part to be heated. ; The magnetic potential vector A augment, the density of the induced currents de- creases and the induced power increases. x 10-4 6 X: 0.05122 Y: 0.0005942 5 module potential magnetic A (A/m2) 4 3 2 1 0 0 0.02 0.04 0.06 0.08 0.1 0.12 rP(mm) Figure 6. Variation of the module potential vector A 20 M. Rezig et al.: Modeling and Optimization of a Goemetrique Induction Cooking for Improving the Temperature x 10-4 x 10-4 6 module potential magnetic(A/m2) 8 5 6 4 4 2 3 040 30 20 20 10 40 60 80 z(mm) r(mm) 2 01 Figure 7. variation of the module Potential vector in plan (r,z) temperature(K°) 340 345 340 335 335 330 330 325 325 320 320 315 315 310 305 310 300 305 295 0 10 20 30 100 40 50 r(mm) z(mm) Figure 8. Distibution of temperature in the plan (r,z) C. 3 Results Table 3. Resultats before and post Optimization Before Optimization Lr2 = 0.04 m Tmax =312.85 K° Aphic max = 0.0002838 T.m Js =5.106 A/mm 2 Post- optimizing Lr2min = 0.021 m Tmax =341.35 K° Aphic max=0.0005942 T.m Js2 =1.62.107 A/mm 2 Js =5.106 A/mm 2 Although a conclusion may review the main points of the paper, do not replicate the abstract as the conclusion. A conclusion might elaborate on the importance of the work or suggests applications and extensions. List of Symbols → A The potential magnetic vector Js The current source density T The temperature µ Magnetic permeability σ Conductivity k Thermal conductivity respectively − q The heat source n Total number of nodes on the load to be heated Ti Current temperature associated to the node Td Desired temperature associated to the node REFERENCES [1] M. FELICHI, ‘‘ La Modélisation Numérique des Phénomènes Couplés ”, Congres International, ‘‘L’Induction dans les Procédés Industriels ”, Paris, France, Mai 1997 [2] K. SRAIRI, M. FELIACHI ‘‘Numerical Coupling Modes for Analyzing Dynamic Behaviors of Electromagnetic Actuators”, IEEE transaction on magnetics, Vol.34, No.5, September 1998 [3] BELLMAN.R.E, KALABA.R.E ‘’Dynamic programming, invariant imbedding and quasilinearization: Comparisons and Interconnections, Computing Methods in Optimisation Problems ’’ Edited by A.Balakrishnan and L. Neustadt, Academic Press, New York, 1964. [4] JIN-KYU BYUN, KYUNG CHOI ‘’ Optimal Design Procedure for a Practical Induction Heating Cooker”, IEEE transaction on magnetics, Vol.36, No.4, Juillet200 [5] JIN-KYU BYUN, KYUNG CHOI, SONG-YHAHN ‘’ Optimal Temperature Control for Induction Heating Devices Using Physical and Geometrical Design Sensitivity”, IEEE transaction on magnetics, Vol.33, No.5, September 1998 5. Conclusions [6] BRUNO SARENI ‘’ Méthodes d’Optimisation Multimodales Associées à la Modélisation Numérique en Electroma- gné- tisme’’ Thèse de Doctorat, L’écol Centrale de Lyon 1992 We have this work, or the use of methods of finite volumes, which allowed carrying out economic models and sufficiently precise and these resulted from our model in per- [7] L.ALLOUI ‘’ Modélisation tridimensionnelle par la méthodes des volumes finis phénomènes Eletromagnetiques et thermiques couplées dans les dispositifs de chauffage par manent mode which do not attain the point of curie. induction ‘’, Thèse de magister, Université Biskra- Algérie, In addition the method of optimization chooses; who 2003 present a good simplicity of the mathematical model and easly implemented, but that remains that the temperature desire is not attained to 100 %. In prospect; and in order to obtain a whole vision on cooking by induction: to raise the [8] FRANCOIS BAY, YANN FAVENNAC ET VALERIE LABBE ‘’ Un exemple de couplage électromagnetiquethermique-mecanique: la modélisation numérique du chauffage par induction ‘’, 16 éme congrès français de mécanique, number of the inductors A conclusion section is not required. septembre 2003, Nice. JASE

... pages left unread,continue reading

Document pages: 4 pages

Please select stars to rate!


0 comments Sign in to leave a comment.

    Data loading, please wait...