An Analysis of Air Flow in the Baking Chamber of a Tunnel-Type Electric Oven

: The baking process in tunnel ovens can be in ﬂ uenced by many parameters. Among these, the most important can be considered as: the baking time, the volume of dough pieces, the texture and humidity of the dough, the distribution of temperature inside the oven, as well as the ﬂ ow of air currents applied in the baking chamber. In order to obtain a constant quality of bakery or pastry products, and for the e ﬃ cient operation of the oven, it is necessary that the solution made by the designers be subjected to modelling, simulation and analysis processes, before their manufacture, and in this sense it can be applied to the Computational Fluid Dynamics (CFD) numerical simulation tool. In this study, we made an analysis of the air ﬂ ow inside the baking chamber of an oven. The analyzed oven was used very frequently on the pastry lines. After performing the modelling and simulation, the temperature distribution inside the oven was obtained in the longitudinal and transverse planes. For the experimental validation of the temperatures obtained in the computer-assisted simulation, the temperatures inside the analyzed electric oven were measured. The measured temperatures validated the simulation results with a maximum error of 7.6%.


Introduction
Baking, even if it is located towards the end of the processing line, is a very complex operation through which the starch-based semi-finished products (bread, pretzels, cookies, etc.) undergo numerous physical, chemical, and biochemical changes, including volume expansion, starch gelation, protein denaturation and coagulation, and the Mayllard reaction [1][2][3][4].These changes occur depending on the baking time, the applied temperature and the speed of heat application, and if these parameters vary, the physical, chemical and bio-chemical transformations will change significantly, [5][6][7][8][9][10].The temperature applied in the baking process varies depending on the ingredients used in the preparation of the semi-finished products, and their baking time depends on both the temperature inside the oven and also shape and size for semi-finished products [11,12].Through baking, the porous structure of the dough turns into a core and the outer surface into a crust [13].In [14] it is shown that the heating of the semi-finished products takes place due to the transmission, primarily, to the outer layers of the dough of heat from the inside of the oven.When the temperature of the layers in the peripheral area of the semi-finished products reaches the water boiling point, the force due to the water vapor pressure will maintain the expansion process of the dough until the end of the process [14,15].Heat transfer from inside the oven to the dough pieces is directly influenced by the temperature of the oven hearth and the environment inside the baking chamber [14].Also, the temperatures of the semi-finished products when being fed inside the oven.Traditionally, the time and temperature of the baking process in an electric oven are controlled, but the presence of airflow affects the temperature distribution inside the oven [16].In [17], a study was performed to identify the parameters that influence finished products' quality.Among other things, parameters such as baking time, volume of dough pieces, texture and humidity of the dough, baking temperature and the flow of air currents applied inside the oven were observed.Electric ovens for baking starch-based semi-finished products can be described as thermally insulated chambers which are utilized for drying, heating, baking or cooking food, [18].The improper placement of heat sources inside the electric oven can cause heat transfer imbalances, leading to temperature variations in the baking chamber.These imbalances will result in an uneven distribution of crust color and moisture content over different areas in the depths of products, [19].During the process, as a result of water removal from the semi-finished products, the air inside the baking chamber retains moisture [4].Therefore, it is necessary to remove it through the vents because the variation of humidity is responsible for the quality defects of the finished products [20].In order to obtain an optimal baking process, in industrial practice, the semi-finished products are baked in ovens at a constant temperature [21].During baking, heat transfer takes place through radiation, conduction and convection.Among these, radiation is that which is more frequently used.Radiation is transmitted heat from the heating sources located inside the oven and from the hot metal surfaces.Convection occurs when air, coming into contact with the heating source and hot metal parts of the oven, transfers heat to the baking surface of the product, and from here a conduction transfer takes place inside the product [3,20,[22][23][24].It should be noted that the process of thermal convection can be natural or forced.
In order to obtain products of a constant quality and for an efficient operation of the oven, it is necessary that the solution made by the designers be subjected to modelling, simulation and analysis processes, before their manufacture, and in this sense the Computational Fluid Dynamics (CFD) numerical simulation tool can be applied [3,25,26].CFD is often used to simulate many processes in the food industry, such as: mixing, ventilation, refrigeration, baking, drying, etc. [27][28][29].The use of the CFD numerical simulation combined with additional models has been studied by the authors of [25] on a forced thermal convection oven.The study aimed to determine the degree of browning of the bread by applying the developed model and to improve baking performance.The objective of [25] was to optimize the design of forced convection ovens by using time-dependent 3D numerical simulations.The baking process is influenced by a number of variables in the case of electric ovens with forced thermal convection: the correct maintenance (at the values set by the operator) of the temperatures in the baking chamber, the thermal conductivity of the materials from which the electric oven is made, the baking time, the temperature distribution in the baking chamber, as well as the humidity at feeding.One of the important parameters that influences baking is the air circulation in the baking chamber.This parameter influences the temperature distribution in the volume of the baking chamber.Forced thermal convection is the forced movement of a fluid caused by an external force moving it (pump, fan, level difference, etc.) [30].
The authors of [31] also proposed to use numerical simulation to study heat transfer by natural convection.They believed that such a simulation provides helpful information for the robust design of an oven and its performance evaluation.The authors of [3] aimed to study the temperature profile inside of a household electric oven using CFD numerical simulation.Three different models of radiation (the Discrete Transfer Radiation Model, the Surface-to-Surface Radiation Model and the Discrete Ordinates Radiation Model) were proposed in the study from [3] to determine the core and crust temperatures of the finished product.The proposed CFD model was able to provide information on the optimum temperature for air inside the oven, the baking time and starch gelatinization.Similarly, by using a CFD numerical simulation, the authors of [32] proposed a model that could predict both the temperature of the air inside the oven and of the dough during the baking process.Also, in [33], the authors proposed an algorithm that aimed to reduce the baking time of semi-finished products in a forced convection oven.The data provided by computer-aided modelling can be of real use for the efficient design of furnaces and the energy improvement of the process.The results of the research in [34] regarding the temperature distribution show the need to perform modelling and temperature analysis in the baking room of an oven.After applying the modified geometry of the oven in the model proposed by the authors of the paper, an improvement in the temperature distribution in the baking room was obtained, limiting these variations to ±2 K.
This paper aims to study the temperature profiles and air currents inside a smallcapacity electric tunnel-type oven.The computer-assisted simulation was validated with the experimental data for the temperatures recorded in four points inside an electric oven.
The article has the following structure: the materials and methods section, in which the concepts of designing the computer-assisted simulation as well as the mathematical equations that were the basis of this study are presented; the results and discussions section, where the results of the analysis are presented and discussions are made on these results, as well as the experimental validation of the simulation; and the conclusions section, in which the major conclusions of the study are presented and recommendations are made.

Materials and Methods
In the case of electric ovens, forced thermal convection is performed by two fans, most often powered at 220 V, which forcefully introduce hot air into the baking chamber.In the case of electric oven Zanolly SYNTHESIS 06/40V E (made by Zanolli Srl from Verona, Italy) [35], analyzed in this paper, baking temperature is achieved by four electrical heaters arranged as follows: two electrical heaters above the movable hearth of the oven and two electrical heaters under the hearth.
The purpose of the present study was to simulate the flow of hot air inside the abovementioned electric oven, according to the methodology presented in Figure 1.The components were assembled in the "Assembly" module.The assembly was performed by introducing geometric relationships between the entities of the parts (faces, edges and points).
After making the assembly, it was introduced in the "Flow Simulation" module to perform the analysis of air flow inside the baking chamber.As only the flow of air inside the oven was concerned, all connections to the external environment (the product supply and outlet ports, as well as the air supply openings for the fans) were closed with the help of lids, as shown in Figure 5.After isolating the inner volume of the chamber, the type of analysis was chosen, as well as the physical parameters, the type of fluid, and the material from which the baking chamber was made.Also chosen were the initial thermodynamic parameters and the thermal conditions of the process; see Figures 6-9.    Figure 8 shows that in the case of external conditions, there was set a heat transfer coefficient calculated at 0.0369 W/m 2 K.The coefficient was calculated with the relation [36]: where:  is the coefficient of thermal conductivity at temperature T;  is the fluid conductivity at 273 K (for air  = 0.0234 W/m 2 K); T is the absolute temperature (for this analysis it was considered 230 °C for the temperature inside the oven; this value can be adjusted); C is the characteristic constant of each gas (for air C = 122 K) [36].
At the same time, we set the wall conditions for the computational domain.Thus, in the feeding and discharge areas of the baking chamber, as well as in the air supply areas of the fans (i.e., where the lids were placed), the normal atmospheric pressure and the ambient temperature of 25 °C were set as analysis conditions (Figure 10).Then, on the surfaces of the blades of the two parametrically shaped fans, we set the virtual characteristics of an axial fan, Comair Rotron, powered in an alternating current, model SU/A (50 Hz).This fan was chosen because the features were the same as the characteristics of the actual fan, mounted on the electric oven, on which the experimental validation was done, but which was not found in the SolidWorks program library.Both fans had been set to introduce air into the electric oven baking chamber.Figure 11 shows the arrangement of the fans, and Figure 12 shows the variation of the air flow with the pressure difference of the chosen fan.After setting the fan parameters, the heat sources inside the baking chamber were set.For all four heaters in the oven assembly, the values for the rate of heat generation were set as follows: for the electric heaters above the conveyor, a value of 1100 W was entered, for the electrical resistance under the conveyor from the supply area of the oven, we set a value of 1700 W, and for the electrical heater in the outlet area under the conveyor, a value of 1750 W was set.These values were determined experimentally by the authors for the same oven and presented in the paper [6].
Before running the analysis, the fineness of the finite volume network in the computational domain was also set. Figure 13 shows the finite volume network.The Navier-Stokes equations (Equation ( 3)), which are used in almost all CFD flow models, were the basis for this study.By solving Navier-Stokes equations, we can predict the pressure and the velocity of one type of fluid in a certain enclosure geometry [37][38][39].We used these equations together with the continuity equation (Equation ( 2)) by conserving momentum.The conservation of mass was represented by the continuity equation.
After solving the equations, values obtained could be applied to any point in the flow of the fluid, and so all flow details could be resolved anywhere in the flow area of an enclosed geometry.In the research literature, scholars often combine those equations with some additional equations (such as the energy equation).

Results
After performing the analysis, the program recorded in its own database the circulation of air and temperature inside the baking chamber.The simulation results took into account the initial conditions set and presented above.Thus, Figure 14 shows the distribution of air currents inside the oven, through several views.From Figure 14, it can be seen that in the middle area of the oven there were created air vortices which caused the temperature in the middle area of the conveyor in the longitudinal plane to decrease to a value of 170-180 °C, under the conditions in which the operating mode requires a baking temperature of 230 °C.This is mainly due to the entrainment of cold air currents, which enter the baking chamber and are entrained in the vortex movement inside.This phenomenon can also be observed in the temperature distribution at the level of the conveyor belt (Figure 15).Chhanwal, N. et al. and Kokolj, U. obtained the same results in their paper [5,23].
Also, it was found that the temperature distribution in the three longitudinal planes was different.In the air supply area (sucked in by fans from the outside environment and fed into the baking chamber), the temperatures were distributed between 20 °C and 230 °C and would increase as the air advanced in the baking chamber to the opposite wall where the temperature distribution was in the range of 160-230 °C, here, however, the temperature of 230 °C was predominant.
It was noted in Figure 16 with Sketch 1, the longitudinal line of the conveyor inside the baking chamber located on the opposite side of the air supply, with Sketch 3, the longitudinal line located in the middle area of the conveyor inside the baking chamber, and with Sketch 2, the longitudinal line located near the feeding area with air inside the baking chamber.All three lines were located on the surface of the conveyor and were used to collect data on the variations of the heat transfer coefficient (Figure 16a), air temperature (Figure 16b), air pressure (Figure 16c) and the specific heat of the air (Figure 16d) along the length of the baking chamber.Figure 17 shows the three lines used to collect the numerical values that define the variation curves in Figure 16.  Figure 16b,c shows that in the middle area of the conveyor (Sketch 3), as the air pressure increases, the temperature inside the baking chamber decreases, with the phenomenon being explained somewhat above.Also, the specific heat of the air increased from 1004 J/kgK (the value that air has at 20 °C, [41]) up to 1104 J/kgK.This phenomenon was also observed by Rek, Z et al. in [26].
Some steps can be taken to avoid heat loss in the middle area of the oven.It is recommended that the air introduced into the oven should be preheated before being fed into it (but not before the air is sucked in by the fans).And, also, to break the vortices that form inside the oven which drive the cold air to the middle zone of the oven and to the surface of the conveyor.A mechanical solution could be the use of some deflector plates, mounted in the areas of maximum intensity of air vorticity, used to breakdown the vortices.
Vorticity is defined in the literature as a mathematical concept used in fluid mechanics that reflects the intensity of the circulation or rotation of a fluid and is most often represented by the angular velocity of rotation, [42].
In the case of the analyzed electric oven, the vorticity (on the three areas presented above) is represented in Figure 18.It can be seen that the vorticity has higher values in the area opposite to the air supply because here the air current hits the vertical wall of the baking chamber and begins to create vortices.Also, it can be observed here that the maximum value of vorticity appears again in the area near the middle of the baking chamber along its length.The results are similar to those obtained by the authors of [25], in which the authors showed a temperature and vorticity gradient with variations identical to the one in this article.It should be mentioned, however, that in the study in [25], this gradient was beneficial to the working process of the oven analyzed, and while it was not in the oven used in this article, this is why we recommend mounting deflector plates in this area.Figure 20 shows the vorticity variation across the width of the baking chamber in the three areas mentioned above.It can be seen that, for the supply area, the maximum values were close to the area of the first fan.In the case of the discharge zone, the vortex variation curve shows maximums towards the left and right ends of the discharge surface, with the vortices being small in the central zone.And, in the case of the line in the middle area of the baking chamber, in width, the maximum was reached at about 0.24 m from the outer wall.The quality of the baking process is mainly influenced by two factors: the constant maintenance of the set baking temperature and the baking time (the time the product spends inside the baking chamber).The first factor can influence the second.Thus, if the temperature were evenly distributed, both longitudinally and transversely, in the baking chamber, and was kept constant at the set value, it is possible that the baking time would be reduced, which would lead to an increase in the productivity of these types of electric ovens.
The non-uniformity of the temperature and the presence of eddies of significant sizes which were constant, are similar to the results of [43], research based on the application of numerical simulation analysis methods inside a forced convection oven, in which the authors proposed adjusting the structure of the deflectors inside the baking chamber and demonstrated that they bring a significant improvement in temperature uniformity compared to the original deflector.
To guarantee an efficient distribution of air flow, thus ensuring a fast baking process, the authors of [44] identified the best furnace design or configuration by studying three furnace design geometries with different fan rotor positions.However, we consider that a preheating of the air with a low-power electric resistance would guarantee a uniformity of temperature in the baking chamber.If this is corroborated the installation of deflector plates in the areas shown above in the article, to break the vortices created, a high efficiency of baking processes can be ensured at constant temperatures.
The design of the fan cover of a bakery oven was studied by the authors of [25] through a methodological proposal and validation involving the combination of computational fluid dynamics and numerical computations to obtain accurate results.The results of the study show that the use of simulation in the design stage gives engineers the opportunity to quickly and efficiently examine new design concepts, which leads to the acceleration of the project development process.

Validation of the Air Flow imulation inside the Baking Chamber of the Analyzed Electric Oven
In order to validate the values of the temperatures obtained during the simulation of the air flow inside the baking chamber of the electric oven, a series of measurements were made with the help of the Extech MP530 MultiPro True RMS multimeter (made by Extech Instruments,Massachusetts, U.S.A.) and an ordinary laboratory thermocouple.
The measurements required a thermocouple that was introduced into the baking chamber of the oven, and the parameters of the operating mode were set according to the simulation of the air flow inside the baking chamber.Measurements were made in sets of five at four important points inside the baking chamber (according to Figure 21), and the results for each point were mediated.Thus, the following average temperature values were obtained: for the lower input area, the temperature value was 214 °C (A1); for the output area on the lower midline, the temperature value was 226 °C (A2); for the uppermiddle zone, the temperature value was 249 °C (A3); and for the lower middle area, the temperature value was 249 °C (A4).The obtained temperatures validate the simulation results with a maximum error of 7.6%, a lower value than the one obtained by Suvanjumrat and Loksupapaiboon in [45], which was 8.99% for the analysis on a batch furnace.

Conclusions
Following the simulation performed with the help of the Solid Works program, the "Flow Simulation" module, we obtained values for some physical parameters (heat transfer coefficient, specific heat, air pressure, air temperature, vorticity) on the length and width of the baking chamber for the air flow inside the electric oven.
The novelty of this article consists in the fact that a study was carried out on an electric furnace that had not been analyzed in this way before, following which recommendations were made for construction and operation.Some steps can be taken to avoid heat loss in the middle area of the oven.It is recommended that the air introduced into the oven be preheated before being fed into it (but not before the air is sucked in by the fans).And, also, to break the vortices that form inside the oven and that drive the cold air to the middle zone of the oven and to the surface of the conveyor.A mechanical solution could be the use of some deflector plates, mounted in the areas of maximum intensity of air vorticity, used to breakdown vortices.Those recommendations were determined by computer-aided modelling-simulation.
We recommend to the builders of forced-convection tunnel-type electric ovens to perform CFD analyses in the design phase of the ovens in order to avoid possible design mistakes.It was demonstrated in this paper, through experimental validation, that the temperature distribution, longitudinally and transversely in the baking chamber at the level of the conveyor, is verified by computer-assisted simulation with a maximum error of 7.6%.

Figure 1 .
Figure 1.Methodology of analyzing the air flow through the electric oven.Firstly, we made 3D geometric models of the oven Zanolly SYNTHESIS 06/40V E's components.The 3D modelling was performed with the parameterized design program Solid Works Premium 2016 S.P. 0.0 (produced by Dassault Systems, France), in the "Parts" module.

Figures 2 -
Figures 2-4 present different views of the obtained models and the interface of the Solid Works program.

Figure 5 .
Figure 5.The option to make the lids from the "Flow Simulation" module of the analysis program.

Figure 6 .
Figure 6.Choosing the type of analysis in the "Flow Simulation" module.

Figure 7 .
Figure 7.The choice of the fluid circulating inside the baking chamber, as well as the flow characteristics in the "Flow Simulation" module.

Figure 8 .
Figure 8. Choosing external parameters in the "Flow Simulation" module.

Figure 9 .
Figure 9. Choosing the initial conditions in the "Flow Simulation" module.

Figure 10 .
Figure 10.Conditions at the borders of the computational domain.

Figure 11 .
Figure 11.The way the fans are placed.

Figure 12 .
Figure 12.The variation of the air flow with the pressure difference of the chosen fan.

Figure 14 .
Figure 14.Distribution of air currents inside the baking chamber: (a) top view; (b) view from the feeding area with baking product; (c) view from the outlet area of the baked product; (d) front view; (e) isometric view.

Figure 15 Figure 15 .
Figure15shows the temperature distribution inside the oven for some characteristic sections.

Figure 16 .
Figure 16.Variation along the length of the baking chamber from some physical parameters: heat transfer coefficient (a), air temperature (b), air pressure (c) and specific heat (d).

Figure 17 .
Figure 17.The lines used to collect the numerical values defining the curves of variation of the heat transfer coefficient, the air temperature, the air pressure and the specific heat of the air along the length of the baking chamber.

Figure 18 .
Figure 18.Variation of vorticity along the length of the baking chamber.

Figure 19 Figure 19 .
Figure19shows with Line 1 the transverse line of the conveyor in the oven outlet area, with Line 3 the transverse line located in the middle area of the conveyor, and with Line 2 the transverse line of the conveyor in the supply area.This figure shows the variation data of the heat transfer coefficient, air temperature, air pressure and air-specific heat in the transverse plane of the baking chamber.

Figure 20 .
Figure 20.Projection of vorticity on the width of the baking chamber.

Figure 21 .
Figure 21.Diagram of the baking chamber of the electric oven and the areas where the temperature measurements were made.A1-down feeding area; A2-lower evacuation area; A3-upper middle area; A4-lower middle area.