1. Introduction
Modern high-end grills are available in various styles and sizes, and can thus handle a broad range of applications [
1]. Grilling has emerged as a commonly used home-cooking method. To achieve a perfect grilled flavor and consistent food quality, the user must accurately adjust the combustion conditions and temperature distribution. Moreover, a suitable grill with controllable grilling settings must be used [
2,
3]. Consequently, the development of high-end barbecue grills has emerged as a research hotspot in recent years [
4,
5,
6,
7,
8,
9,
10].
With such research, innovative heating methods for barbecue grills [
11] such as infrared heating, have been introduced in addition to the use of traditional fuels. Through such methods, users can precisely control the cooking process for different ingredients. Nevertheless, there remains room to improve conventional grills, especially in terms of controlling the surface temperature of the grill plate to ensure a uniform temperature distribution [
3,
12,
13].
In recent years, outdoor barbecue has emerged as not only a method to prepare food, but also a popular mode of socialization [
14,
15,
16,
17,
18,
19]. The corresponding grills have evolved over time, and various fuel systems such as charcoal, natural gas, and propane, have been introduced to address different requirements [
15,
20].
Although combustion systems may have different types of energy sources, the most economical and convenient fuel systems correspond to compressed gases or liquids. The combustion system considered in this study uses propane and natural gas, which represent efficient fuel sources. Propane fuel, which is prevalent in the market, is commonly stored in pressurized tanks and available in liquid form with various weight percentages.
The process of cooking food by burning gas fuel typically involves three heat transfer modes: (1) Heat convection: The working fluid (e.g., air) around the food is heated, which then transfers energy to the food [
21]. (2) Heat conduction: An intermediate heating material, such as a heating pan, is heated and placed in direct contact with the food to transfer energy to the food [
22]. (3) Heat radiation: Many existing ovens use infrared (IR) burners to generate heat in the form of infrared rays, which deeply penetrate the surface of the food during cooking, allowing users to produce restaurant-quality barbecue food [
23].
The objective of this study was to optimize grill systems that use natural gas and propane fuels. The original design involves an iron grill pan with multiple gaps, as shown in
Figure 1a. Owing to the excellent heat transfer characteristics of the iron plate, the grid-like pores promote air convection, which can alleviate the problems associated with combustion and uneven temperature. However, iron pans are not suitable for processing several foods, such as pizza. In such instances, stone grill plates must be used. Nevertheless, stone grill plates come with their own set of problems: First, stone grill plates are prone to fracturing, as shown in
Figure 1b. Second, the combustion temperature is more uneven than that of iron plates, which may generate defects in the product. Therefore, the reasons for the fracturing of stone grill plates and uneven combustion temperatures must be investigated.
Because the combustion system and furnace structure are typically mass-produced, the product cannot be significantly modified. Therefore, optimization techniques must be used, for instance, to enhance the uniformity of the temperature distribution on the plate surface, reduce the thermal stress, and prevent breakage during use. Moreover, it is desirable to modify the product such that the manufacturing and overall cost are not significantly affected. Notably, the stress and temperature distribution of a stone grill plate cannot be visually evaluated. Other methods must be used to clarify the influence of the design of the burner tubes and heat transfer characteristics on the temperature and thermal stress across the grill surface.
Therefore, in this study, a number of simulations were performed to predict the temperature and thermal stress of bakeware, and to analyze the thermal stress model by combining a computational fluid dynamics-based combustion flow model [
24,
25,
26,
27] and heat conduction, as well as the finite element method. The temperature and thermal stress distribution of the grill plate were determined by simulating the combustion of mixed gas and heat conduction, and then the two factors were taken to optimize the barbecue design. In addition, a thermal imager was used to collect temperature data to verify the consistency of the actual and analytical models [
28,
29,
30]. The optimized design was then created to exhibit a more uniform temperature distribution pattern than the original design [
31,
32].
Many scholars improve a problem in one way, but often a non-optimal design results. The breakthrough approach of this study is that it not only improves the design of the stone grill plate, but also improves the air intake design of the fuel line. In such a way, problems such as the cracking of the stone grill plate caused by the temperature difference and thermal stress difference caused by the opening of the upper cover of the oven can be solved all at once.
2. Materials and Methods
As mentioned, the objective of this study is to explore the reasons for the fracture of stone grill plates during use, and to enhance the design of the fuel tube and the overall performance of plate and grill. To this end, Ansys CFX simulation analyses were performed, as illustrated in
Figure 2.
The first step was to clarify the causes of the failure of the stone grill plate and identify a solution. The primary purpose is to theoretically analyze the problems associated with the existing barbecue grill system, and identify solutions to alleviate the uneven temperature distribution and thermal stress on the grill surface. The second step was to examine the fuel tube design. In the third step, three stone grill plate designs were created to evaluate the failure characteristics and develop the optimal design.
At the same time, flamelet and k-Epsilon turbulence models were used in the oven combustion and flow field physics model. The stone grill plate’s heat conduction method was also taken into consideration to simulate the steady-state temperature distribution. The inlet and outlet boundary conditions of the oven chamber were set at 101,325 Pa, and the ambient temperature set at 25 °C. The coordinates were placed in the center of the stone grill plate, and a small area in the center was designated as a fixed boundary so that the stone grill plate can generate thermal expansion. With that, thermal stress could be calculated.
Fourthly, a cross-simulation analysis of the designed fuel tube and stone grill plate was performed. Additionally, an analysis system was established for the oven, measurement, and the testing equipment. Several suggestions were deduced for improving the scientific design process and solving the problems of the similar oven systems. Finally, after validating the stone grill plate and corresponding outputs, the optimal design was developed based on the simulation analysis and verification results.
2.1. Equipment
The grill structure includes a top cover, bottom housing, grill plate, tube burner, and fuel orifice, as shown in
Figure 3. The top cover has venting pores to provide combustion heat, and the exhaust gas can flow out of the oven cavity when the cover is in use. The bottom housing includes an oil collecting pan and a base cover to form a gap to provide the air required for combustion (in the bottom air inlet).
2.1.1. Model Analysis
Because the test stone grill plate typically fractures when the grill is opened for grilling, it can be assumed that when the top cover is opened, the convection of the cold air above the surface of the grill plate is intensified, resulting in a rapid temperature decrease over the pan surface. In contrast, continuous heating occurs due to the flame of the bottom combustion tube. These phenomena generate a significant temperature difference and thermal stress. At the point of opening the grill, the thermal images recorded by the thermal imager are consistent with those of the analysis model. Therefore, the image corresponding to the point at which the lid is opened is examined to analyze the stress distribution of the grill plate.
2.1.2. Boundary Conditions
The thermal power specification of the barbecue oven is 10,500 BTU/h. In the combustion simulation, the temperature is less than 25 °C, and the measured oven propane flow rate is 220 g/h, set as the boundary condition of the orifice. Gas (propane) enters the gas mixture combustion tube through the nozzle, and the gas mixture combustion tube is designed based on the Venturi principle. Specifically, ambient air is sucked into the mixed combustion tube, and it flows out through the multiple holes in the mixed combustion tube. When combustion is initiated, the control gas valve is opened, allowing the gas to pass through the gas mixture pipe for the first gas mixing. The ignition device set above the air–fuel mixture combustion hole ignites the gas.
In the first mixing process, the amount of air required for complete combustion is typically inadequate. Therefore, the gas flows outside the air-fuel mixture combustion tube from the combustion hole. At the same time, the ambient air undergoes the second aeration process to enable complete combustion. The design of such combustion tubes is relatively mature. In general, various specifications are available when designing a barbecue, and thus, no redesign is required, as shown in
Figure 4.
Moreover, the combustion system and furnace structure to be improved are typically mass-produced and must be suitable for various iron grill pans. Therefore, the modification of the furnace body and combustion tube design is highly restricted. Considering these aspects, limited modification and simulation analyses were performed in this study based on furnace body specifications.
Table 1 presents the material specifications of the stone grill plate considered to establish the analysis model.
On the other hand, the Finite Volume Method analysis grid used 0.8 cm tetrahedron elements, while the square size of the combustion tube grid was refined to 0.1 cm each to guarantee the analytical model convergency. For the stone grill plate, Finite Element analysis grid with a tetrahedron element of 0.8 cm each was adopted, but the square size of the grid in the heating area on the bottom surface of the stone grill plate was refined to 0.2 cm. The total number of squares in the grids was about 2.9 million, the flow field was 1.7 million, and the stone grill plate was 1.2 million. The same setting and same method to increase the density were applied in the different designs of the grid density analyses. The same setting and benchmark would facilitate our analysis of the differences among the testing models. Grid distribution and grid independency are shown in
Table 2.
2.2. Mathematical Model
As the stone grill plate rupture is caused by excessive thermal stress, it is also true that thermal stress is affected by many parameters. Under the existing configuration and constraints, this study proposes a new set of design parameters, including the position and method of air intake in the combustion tube, to improve the air/fuel ratio and average flame distribution of the combustion tube after conducting analyses and experiments. At the same time, the shape of the structural ribs at the bottom of the stone grill plate was modified to reduce local high temperature and thermal stress concentration. The design parameters of the stone grill plate and the combustion tube work together in a way that they not only affect the flame distribution state during combustion, but also have a direct heating impact on the stone grill plate, and the distribution of thermal stress.
A RANS-based approach with the k-ε model was selected to further exemplify this case [
33]. The continuity and momentum equations in RANS format were introduced as shown in Equations (1) and (2):
2.3. Original Design Analysis Results
The design simulation analysis was divided into two parts, corresponding to the original combustion system and stone grill plate The grill plate was changed from an iron grill plate with gaps to a stone grill plate with no openings on the surface. In actual use, the problem of uneven temperature is quite apparent in such systems.
In the simulation analysis of the combustion air intake, owing to the large shading caused by the stone grill plate, most of the air intake at the bottom air inlet is directed to the bottom of the slate and input to the combustion air mixture pipe, as shown in
Figure 5. This configuration leads to several undesirable results. First, the volume flow rate of air sucked by the gas mixing pipe is fixed. If the intake air temperature increases, the air density decreases, and thus, the air mass drawn into the combustion tube per unit time decreases for a given volume flow. In the analysis, the intake air temperature reaches 140 °C. The density decreases compared to that at an intake air temperature of 25 °C, and thus, the air mass drawn into the combustion tube per unit time decreases by 28%. This phenomenon hinders complete combustion. The second result pertains to the air intake method. The system inputs exhaust gas with a high concentration of carbon dioxide, and thus, the oxygen content is reduced.
Figure 6a shows the temperature streamlines indicating the burner distribution during the second aeration. After the streamlines flow to the upper part of the concentrated combustion, they flow to the periphery of the stone grill plate. No heat energy is transferred to the middle part of the stone grill plate. This phenomenon generates a large temperature difference between the middle of the stone grill and top of the burner when the top cover is opened. This situation is consistent with that observed in the image recorded by the thermal imager, as shown in
Figure 7b. Notably, the intake port streamline distribution shown in
Figure 6b indicates that the original design does not produce excessively many abnormal parts. Nevertheless, the air entering the combustion tube mixture does not accumulate at the air mixture intake port.
Figure 7 shows the surface temperature of the original grill plate and image recorded by the thermal imager. In
Figure 7a, the orange-red part corresponds to a temperature of 266–320 °C, and the highest temperature observed in
Figure 7b is approximately 300 °C. The high-temperature regions in
Figure 7a,b are similar, and highlight that the area with the highest temperatures pertains to the fuel orifice.
Figure 7c shows that the thermal stress on the surface of the grill pan is concentrated on the combustion pipeline. The maximum and minimum thermal stresses are 160 and 110 MPa, respectively, corresponding to a difference of 50 MPa.
Figure 7d shows the temperature distribution at the bottom of the original grill plate. The maximum temperature is approximately 400 °C, and the distribution is concentrated over the combustion line.
Figure 7e shows the thermal stress distribution at the bottom of the grill plate. The thermal stress is accumulated over the combustion pipeline, and the surface thermal stress is 250 MPa. The difference between this stress and that in the lowest thermal stress region (49 MPa) is approximately 200 MPa. In other words, the high-temperature area is consistent with the fracture area of the solid stone grill plate, and the stone grill plate fractures in the maximum stress area. The corresponding fracture mode can be observed from the thermal stress analysis results, as shown in
Figure 7f.
4. Discussion
In this study, ANSYS CFX analyses were conducted utilizing a coupled analysis model method based on combustion thermodynamics, heat transfer, and finite element methodologies. The objective is to prevent the fracturing or failure of barbecue stone grills. The process flow has been described in
Section 2 and
Figure 2. The initial goal is to understand the causes of the problem. Our analysis and experiment results indicated that the cracks on stone grilling plates were due to two main reasons. Further thermal imaging results indicated that the heat of the stone grill plate was concentrated at the outlet of the fuel nozzle. Moreover, the temperature difference on the stone grill plate was extremely large and generated local thermal stress, which caused the stone grill plate to crack in the location opposite to the fuel pipe. The second step is to improve the temperature of the mixed air entering the combustion tube to increase the flow rate and velocity of the airflow through the top of the flame to the bottom surface of the stone grill. Additionally, it is necessary to ensure that the temperature of the stone grill plate is maintained at the same level. The third step is to improve the design of the continuous ribs and quantify the improvements through simulations. A cross-analysis was performed by integrating the second and third steps, and optimal results were obtained. Finally, a physical experiment was conducted to confirm that the improved stone grill plate generated a uniform temperature across its entire surface. The proposed method can effectively prevent the fracturing of the stone grill plate caused by thermal stress and thus, its service life can be prolonged.
Based on the above, the analytical model of the stone grill plate must work hand in hand with the analysis of combustion, heat transfer, fluid flow, and thermal stress to solve the problem of thermal stress fracture. For complex nonlinear multivariate coupled models, accurate results cannot be obtained from a single analytical model. The coupling model relationship is shown in
Figure 18. Both the experimental and measurement generated the same results as the simulation analysis of this study, which proves the validity of this research.
Table 4 and
Table 5 present the quantitative results for the different stone grill plates. The temperature and stress differences represent the differences between the maximum and minimum temperatures and maximum and minimum thermal stresses, respectively. The temperature and thermal stress differences are also presented as line graphs in
Figure 17b. The broken line graph in
Figure 17b shows that the temperature difference on the surface of the stone grill plate decreases in the order CA#3 > CA#2 > CA#1. This decrease tends to decrease the surface thermal stress. Therefore, the improved stone grill plate can effectively prolong the service life.
5. Conclusions
A coupled analysis model was established by combining combustion thermodynamics with heat transfer and the finite element method. The causes of barbecue oven fracturing using a stone grill plate were studied. Simulated mixed air distributions were used to determine the temperature and thermal stress distribution of the stone grill plate after combustion, which were then used as the basis for design optimization. Thermal images were recorded to quantify and verify the accuracy of the simulation. The results demonstrated that narrowing the temperature range can reduce the thermal stress on the surface.
Moreover, the temperature across the stone grill plate could be maintained at the same level by reducing the temperature of the mixed air entering the combustion tube, and increasing the flow rate and velocity of the airflow through the top of the flame to the bottom surface. When the thermal stress was reduced, fracturing could be prevented, and the service life of the stone grilling plate could be prolonged. The proposed simulation model provided an optimized design while reducing the errors typically encountered in experiments, and could shorten the experiment duration.
The critical result of this research is that oven performance can be improved by simply adjusting the intake position and the combustion tube method under the same calorific value of propane (10500BTU/h) conditions, without the need to change the oven structure. The analysis results demonstrated that the intake air temperature of the combustion tube could be reduced from 152 °C to 67 °C, and the temperature gap is about 85 °C. Therefore, under the same airflow volume, the air mass density increased from 0.83 kg/m3 to 1.039 kg/m3, and the air mass flow could be increased by 12%. Increasing the air-to-fuel mixing ratio in the combustion tube could effectively improve the efficiency of the combustion tube. Incidentally, the new combustion tube design proposed in this study can be widely applied in the available oven models in the market. The improved combustion tube design is in itself a significant invention, and thus, a patent application has been filed.
In the future, the technology employed in this research will be applied to related research projects, and thereby the cost and time of experiments can be reduced. At the same time, it is hoped to expand the applicability of this research to more countries with different climates. Hence, a portfolio embodied different climatic conditions can be established. By then, more parameters will be verified, and the geographic limitation of the research can be overcome.