The perfect pan for ovenThe heat transfer in the oven includes heat conduction, heat radiation and heatconvection. We use two-dimensional Fourier heat conduction equation ∂u∂t −α(ð2u∂x2+ð2u∂y2)=f(x, y, t) to make a research on distribution of heat for the pan. Heat source heats the pan by heat radiation. The pan interacts with air in the oven in the way of natural convection, so the pan realizes heat dissipation.We calculate heat radiation based on radiation ability of heat source and heating tube area. We use heat dissipation function to show the pan's different parts' loss of heat caused by natural convection. Both of them consist in heat source function f.The area of the pan is fixed at 0.085m2in this paper. When comparing temperatures at the edges of rectangular pans with different length to width ratios ξ, we can get that the smaller ξ is, the lower the temperature of the edges is. But as long as it is still a rectangle, the amplitude of its drop won't be very big. When we make the pans with fixed area vary from square to round square to round, we find that the bigger the fillet radius is, the lower the temperature of its corners is and the extent of temperature's reducing is large.We fix the bottom area of the oven and area of the pan. Through study, we find that round square's capacity for uniform distribution of heat is far higher than other shape's (except round). The larger the fillet radius of the round square is, the larger the pan’s waste of space is. But heat distribution is more uniform. We work out the optimal solution of pan’s size under different weights p through optimizing the relationship between two conditions. Then we get several oven's width to length ratios of W/L by arranging the pans with the optimal size.I. IntroductionThe temperature of each point in the pan is different. For a rectangular pan, the corners have the highest temperature, so the food is easily overcooked. While the heat is distributed evenly over the entire outer edge and the product is not overcooked at the edges in the round pan.To illustrate the model further, the following information is worth mentioning1.1 Floor space of the panThe floor space of each pan is not the square itself necessarily. In this paper, there are 3 kinds of pans with different shapes, as rectangular pans, round pans and round rectangle pans.For rectangular pan, the floor space is the square itself, and the pans can connect closely without space.For round pans, the diagrammatic sketch of the floor space is as follows:shade stands for the round pan; square stands for the floor spaceFigure 1Round pans have the largest floor space for a certain area. The space between each pan is larger than other two kinds of pans. The coefficient of utilization for the round pans is the lowest.For round rectangle pans, the diagrammatic sketch of the floor space is as follows:shade stands for the round rectangle pan; square stands for the floor spaceFigure 2If the area of the round rectangle is the same as the other two, its floor space is between them. The coefficient of utilization of oven decreases with the radius of expansion.1.2 Introduction of ovenThe oven is usually a cube, no matter it is used in home or for business. A width to length ratio for the oven is not a certain number. There are always two racks in the oven, evenly spaced. There are one or more pans on each rack. To preserve heat for the oven, food is heated by radiation. Heating tube can be made of quartz or metal. The temperature of the tube can reach 800℃ high when the material is quartz. The heating tube is often in the top and bottom of the oven. Heating mode can be heating from top or heating from bottom, and maybe both[1].1.3 Two dimensional equation of conductionTo research the heat distribution of pan, we draw into two dimensional equation[2]of conduction:∂u ∂t −α(ð2u∂x2+ð2u∂y2)=f(x, y, t)In this equation:u- temperature of the pant- time from starting to heatx- the abscissay-ordinateα- thermal diffusivityf- heat source functionThe heat equation is a parabolic partial differential equation which describes the distribution of heat (or variation in temperature) in a given region over time. The heat equation is of fundamental importance in diverse scientific fields. In mathematics, it is the prototypical parabolic partial differential equation. In probability theory, the heat equation is connected with the study of Brownian motion via the Fokker–Planck equation. The diffusion equation, a more general version of the heat equation, arises in connection with the study of chemical diffusion and other related processes.II. The Description of the Problem2.1 The original problemWhen baking in a rectangular pan heat is concentrated in the 4 corners and the product gets overcooked at the corners (and to a lesser extent at the edges). In a round pan the heat is distributed evenly over the entire outer edge and the product is not overcooked at the edges. However, since most ovens are rectangular in shape using round pans is not efficient withrespect to using the space in an oven.Develop a model to show the distribution of heat across the outer edge of a pan for pans of different shapes -rectangular to circular and other shapes in between.Assume1. A width to length ratio of W/L for the oven which is rectangular in shape.2. Each pan must have an area of A.3. Initially two racks in the oven, evenly spaced.Develop a model that can be used to select the best type of pan (shape) under the following conditions:1. Maximize number of pans that can fit in the oven (N)2. Maximize even distribution of heat (H) for the pan3. Optimize a combination of conditions (1) and (2) where weights p and (1- p) are assigned to illustrate how the results vary with different values of W/L and p.In addition to your MCM formatted solution, prepare a one to two page advertising sheet for the new Brownie Gourmet Magazine highlighting your design and results.2.2 Problem analysisWe analyze this problem from 3 aspects, showing as follows:2.2.1 Why the edge of the pan has the highest temperature?The form of heat transfer includes heat radiation, heat conduction and heat convection. Energy of heat radiation comes from heat resource. The further from heat resource, the less energy it gets. Heat conduction happens in the interior of the pan, and heat transfers from part of high temperature to the part of low with the temperature contrast as its motivation. With the two forms of heat transfer above, we find the result is that pan center has the high temperature and the boundary has the low. D epending on that, we can’t explain why the product gets overcooked at the corners while at the edges not. We think that there is natural convection between pan and gas, because the temperature of pan is much higher than that of gas. The convection is connected with the contact area. The point in pan center has a larger contact range with air, so the energy loss from convection is more. While the point in the corners of a rectangular pan has a narrow contact area, the energy loss is less than that in the inner part. Because that above, the energy in the pan center is more than that in corner, and the corners have higher temperature.2.2.2 Analysis of heat distribution in pans in different shapesThe shape of pan includes rectangle, round rectangle and round. When these pans' area is fixed, the rectangles with different shapes can be shown with different length to width ratios. Firstly, we study temperature (maximum temperature) in the corners of rectangles with different length to width ratios. Then we study how temperature in the corners changes when the pans vary from square to round square to round. After that, we select some rectangular panand make it change from rectangle to round rectangle to study changes of temperature in the corners. To calculate distribution of heat for the pan, we would use three main equations. The first one is the Fourier equation, namely heat conduction equation, the second one is the radiation transfer equation of heat source and the third one is the equation of heat dissipation through convection. In the radiation transfer equation of heat source, we take heat source as a point. We get its radiating capacity through its absolute temperature, blackening and Stefan-Boltzmann law. We combine radiating capacity with surface area of quartz heating tube to get quantity of heat emitted by heat source per second, then we can get heat flux at each point of the pan. In the equation of heat loss through convection, Heat dissipating capacity is proportional to area of heat dissipation, its proportional coefficient can be found from the related material. Through the establishment of above three main equations, we can use pdetool in matlab to draw the figure about distribution of heat for the pan.2.2.3 How to determine the shape of the pan?To make heat distribution of the pan uniform, we must make it approach round. But under the circumstances of the pan's fixed area, the closer the pan approaches round, the larger its floor space is. In other words, the closer the pan approaches round, the lower the utilization rate of the oven is.More uniform distribution of heat for the pan is, the lower temperature in the corners of the pan is. Assuming the bottom area of the oven is fixed, the number of most pans which the oven can bear is equal to the quotient of the bottom area of the oven divided by floor space per pan. So in a certain weight P, we can get the best type of pan (shape) by optimizing the relationship between temperature in the corners and the number of most pans which the oven can bear. Then we get the oven's width to length ratio of W/L by arranging the pans with the optimal size.2.3 Practical problem parameterizationu: temperature of each point in the oven;t: heating time;x: the abscissa values;y: the ordinate value;α:thermal diffusivity;ε: degree of blackness of heat resource;E: radiating capacity of heat resource;T: absolute temperature of heat resource;T max: Highest temperature of pans’ edge;ξ: the length to width ratio for a rectangular pan;R: radius for a round pan;L: length for the oven;W: width for the oven;q: heat flux;k: coefficient of the convective heat transfer;Q: heat transfer rate;P: minimum distance from pan to the heat resource;Other definitions will be given in the specific models below2.4 Assumption of all models1. We assume the heat resource as a mass point, and it has the same radiation energy in all directions.2. The absorbtivity of pan on the radiation energy is 100%.3. The rate of heat dissipation is proportional to area of heat dissipation.4. The area of heat dissipation changes in a linear fashion from the centre of the pan to its border.5. Each pan is just a two-dimensional surface and we do not care about its thickness.6. Room temperature is 25 degrees Celsius.7.The area of the pan is a certain number.III. ModelsConsider the pan center as origin, establishing a coordinate system for pan as follows:Figure 33.1 Basic ModelIn order to explain main model better, the process of building following branch models needs to be explained specially, the explanation is as follows:3.1.1 Heat radiation modelFigure 4The proportional of energy received by B accounting for energy from A is4π(x2+y2+P2) The absolute temperature of heat resource A (the heating tube made of quartz) T= 773K when it works. The degree of blackness for quartz ε=0.94.The area of quartz heating tube is 0.0088m2.Depending on Stefan-Boltzmann law[3]E=σεT4 (σ=5.67*10-8)we can get E=18827w/m2.The radiation energy per second is 0.0088E.At last, we can get the heat flux of any point of pan. =0.0088E4π(x2+y2+P2)Synthesizing the formulas above, we can get:=13.24π(x2+y2+P2)(1) 3.1.2 Heat convection model3.1.2.1 Round panheat dissipation area of round panFigure 4When the pan is round, the coordinate of any point in the pan is (x,y). When the point is in the centre of a circle, its area of heat dissipation is dxdy. When the point is in theboundaries of round, its area of heat dissipation is 12dxdy. According to equation of heatPdissipation through convection dQ=kdxdy and assumption that the area of heat dissipation changes in a linear fashion along the radius[4], we can get the pan's equation of heat dissipation:q =k[1-0.5(x2+y2)0.5/R] (2) 3.1.2.2 Rectangular panHeat dissipation area of rectangular pan(length: M width: N)Figure 5When the pan is rectangular, the coordinate of any point in the pan is (x,y). When the point is in the centre of a rectangle, its area of heat dissipation is largest, namely dxdy. Whenthe point is in the center of the rectangular edges, its area of heat dissipation is 1dxdy. When2the point is in rectangular vertices, its area of heat dissipation is minimum, namely According to equation of heat dissipation dQ=kdxdy and assumption that the area of heat dissipation changes in a linear fashion from the centre of a rectangle to the center of the rectangular edges.The area of any other point in the pan can be regarded as the result of superposing two corresponding points' area in two lines.Heat dissipating capacity of any point is:=1−y/N(3)1−x/M3.2 Pan heat distribution Model3.2.1 Heat distribution of rectangular pansAssume that the rectangle's length is M and its width is N,the material of the pan is iron.For rectangular pans, we change its length to width ratio, establishing a model to get thetemperature of the corners (namely the highest temperature of the pan).We assume the area of pan is a certain number 0.085m2,The distance from the pan above to the top of the oven P=0.23m,By checking the data, we can know that the coefficient of the convective heat transfer k is approximately 25 if the temperature contrast between pan and oven is 100~200℃.then get a several kinds rectangular pans following:In the two dimensional equation of conduction,we get the thermal diffusivity of iron is 0.000013m2/s through checking data.Heat source function equals received thermal radiation minus loss of heat caused by heat dissipation through ly equation (1) minus equation (3).In this way, we get a more complicated partial differential equation. For example, through analyzing the NO.1 pan, we can get the following partial differential equation.∂u ∂t −α(ð2u∂x2+ð2u∂y2)=13.24π(x2+y2+0.529)−1−y/0.291551−x/0.29155It is hardly to get the analytic solutions of the partial differential equation. We utilize the method of finite element partition to analyze its numerical solution, and show it in the form of figures directly.Pdetool in matlab can solve the numerical solution to differential equation in the regular form quickly and show distribution of heat by the three-dimensional image[5]. We enter partial differential equation of two-dimensional heat conduction into it, then we get the figures about distribution of heat in different pans.When we use pdetool, we take Neumann condition as boundary condition, and we suppose that the boundary is insulated, In fact, it is not insulated, and the heat dissipation will show in the heat source function.We heat the pan for 8 minutes no matter what kind of shape the pan is. By Pdetool, the heat distribution of each pan shows as follows:Figure 6(heat distribution for NO.1 pan) Figure 7(heat distribution for NO.2 pan) Figure 8(heat distribution for NO.3 pan)This pan is just a square pan, with the lowest length to width ratio. From the figure, we can know that corners have the highest temperature which is 297.7℃.The corners have the highest temperature for the pan, which is 296℃The corners have the highest temperature for the pan, which is 294.8℃The temperature ofcorners for the pan isslightly lower than thehighest temperature, andthe highest temperature is293.7 ℃.Figure 9(heat distribution for NO.4 pan)To get a more accurate relationship between the length to width ratio(ξ) and highest temperature(T max), we make several more figures of heat distribution based on the different length to width ratio. At last, we can get its highest temperature. The specific result is as follows:ξ is the argument and T max is the dependent variable. The points in the chart are scaled out in the coordinate system by mathematical software Origin. Connecting the points by smooth curve, we can get the figure following:Figure 10(the relationship between ξ and T max )From the figure we can find that, with the length to width ratio increases, the temperature of corners will sharply fall at first, and it is namely that the heat distributes evenly . When the length to width ratio increases further, the temperature of corners drops obscurely . When the ratio reaches about 2.125, the length to width ratio of rectangular pans has little effect on the heat distribution. In contrast, the ratio is too big, it is difficult for practical application.In addition, we can also find that as long as the shape of the pan is a rectangle. The highest temperature in the corners of the pan won't change a lot whether its length to width ratio changes, From the figure 10, we can find that maximum range is about 6 degrees Celsius. So in general, it is very difficult to change high temperature in the corners of the rectangular pan.3.2.2 Heat distribution of square pans to round pans 3.2.2.1 Size definition of round squareWe need some sizes to define round square, our definition isas follows:T maxThe size of round square isdecided by l and r (l stands forthe length of the straight flange;r stands for the radius of thefillet).Figure 11We assume that the area of the pan is 0.085m2. The r of round square has its range, which is 【0，0.164488】，When the r reaches the two extremums, round square becomes square and circle.3.2.2.2 ModelWhen we make this kind of pans’ heat distribution figures, we take the heat source function as (1) and (3). When the round square becomes circle, the heat source function is (1) and (2).We get a several round squares with different r and l, and take a circle as an example. The specific examples show in the table below:The pan is still made of iron. Here we also heat the pan for 8 minutes.By Pdetool, we can obtain their heat distribution figures as follows:Figure 11(heat distribution for NO.1 pan) Figure 12(heat distribution for NO.2 pan) Figure 13(heat distribution for NO.3pan)The temperature of the corners about the pan is relatively low, and the highest temperature is 264℃. On the contrary, the heat distributes quite evenly.The highest temperature of the pan is 271.1℃The highest temperature of the pan is 274℃The highest temperatureof the pan is 279.8℃Figure 14(heat distribution for NO.4pan)The highest temperatureof the pan is 283.7℃Figure 15(heat distribution for NO.5pan)To get the different relationship between l/r and T max, we take l/r as argument, and T max as dependent variable. We change l and r of round square, and make several heat distribution figures. In this way, we can get the highest temperature T max, and the results show in the table below:Here, we alsouse the mathematical drawing software origin. We use l/r and T max to express coordinates. We use smooth curve to connect points, then we can get the trend line.Figure 16(the relationship between l/r and T max )From the figures we can know that, the value of l/r is smaller, the temperature of the corners about pan is higher. It is namely that the pan is more closely to circle, and heat distribution is more evenly. With the value of l/r increases, the temperature of corners rise very quickly at first, then the amplitude is getting smaller. When the value of l/r is infinitely great, the highest temperature of the pan go to a certain number.Analyzing in a theoretical way , the shape of pan goes to square when the value of l/r is infinitely great. At the same time, the temperature of the round square’s corners approach to that of square’s. From the figures, we can know that this function has a upper boundary ,T maxwhose value is close to the corner temperature of square. Through this, we can verify the correctness of our models.3.2.3 Heat distribution of round rectangle (except round square)From the model about heat distribution of rectangular pan, we can learn that drop of temperature in the corners of rectangular pans will be very little when its length to width ratio is bigger than 2.125. So we select the rectangle with length to width ratio of 2.125. We let the pan vary from the rectangle to round rectangle. So we can study changes of temperature in the corners of the pan.3.2.3.1 Size definition of round rectangleThe specification of the roundsquare is decided by l and r. Weset its width the same as therectangular pan before, namely0.2m.(l stands for the length ofstraight long side, and r stands forthe radius of the fillet.)Figure 17We assume that the area of the pan is 0.085m2, For round rectangle pans, with r decreases constantly, the pan finally approaches the rectangular pan before .If the r increases constantly, its shape will become that of playground. The range of r is 【0,0.1】3.2.3.2 ModelWhen we draw the figure about heat distribution of this kind of pan, heat source function equals equation (1) minus equation (2).We can get some different round rectangles by changing the value of r and l. Their detailed specifications are shown in the following table:The pan is still made of iron. Here we also heat the pan for 8 minutes.By Pdetool, we can obtain their heat distribution figures as follows:The corner temperaturewhich is 291.7℃andslightly lower than thehighest temperature ofthe pan.Figure 18(heat distribution for NO.1pan)The corner temperaturewhich is 292.54℃andslightly lower than thehighest temperature ofthe pan.Figure 19(heat distribution for NO.2pan)The corner temperaturewhich is 293.5℃andslightly lower than thehighest temperature ofthe pan.Figure 20(heat distribution for NO.3pan)To get the different relationship between l/r and T max, we take l/r as argument, and T max as dependent variable. We change l and r of round square, and make several heat distribution figures. In this way, we can get the highest temperature T max, and the results show in the table below:Here, we also use the mathematical drawing software origin. We use l/r and T max to express coordinates. We use smooth curve to connect points, then we can get the trend line.Figure 21(T max)From the figures we can get that, with l/r increases, the highest temperature of pan goes down. That is to say, the bigger radius of the fillet is, the more evenly heat distributes. In addition, with l/r increases, T max rises quickly at first, then the extent is smaller. When l/r is infinitely great, the round rectangle pans become rectangle pans, and the temperature approaches to that of the rectangle pan before.Secondly, from the figure, we can learn that change of the largest temperature is very little and its largest temperature is all very high when the pan varies from rectangle to round rectangle. compared with round square pan, round rectangle pan has much worse capacity of distributing heat.3.3 Best type of pan selection ModelFrom the model of heat distribution, we can know that the extent of heat distribution for round square pan is more than the extent of other pans with other kinds of shape( except T maxcircle). It is difficult to accept it for people because the food is easily overcook, no matter what the number of pan is. Depending on that, people will choose the round square pan.In the following models, we mainly discuss the advantages and disadvantages of round square pans with different specifications.3.3.1 Local parametersFor study's convenience, we take commercial oven of bottom area 1.21m2as an example. The distance between heating tube on the top and the nearest rack is P=0.23m.The number of pans which the oven can contain is n;The fillet radius of round square is r;The floor space of pan is S;The weight of the number of pans in the oven is P;The area of pan is still 0.085m2, the material is still iron.3.3.2 The relationship between T max and rThrough the models before, we know the relationship between l/r and T max. We transform it as the function of r and T max, and shows in the form of table below:We take r as the argument, and T max is the dependent variable. Making the dots in the coordinate system by the mathematical software Origin[6].The figure is as follows:Figurebetween r and T max )To our surprise, we can find that there is nearlylinear relation between r and T max fromthe figure. We may fit the relation with a linear function, so we can get the function of r and T max .FigureThe function we are getting is: T max =-174.5r+291.5 (4)3.3.3 The relationship between n and rWe have introduced that the floor space is not the area of the pan itself in theT maxT maxIntroduction. So we can get the formula below of the area of round square pan and r:S=r2(4−π)+0.085We have already known the floor space of the oven. So we can know the maximum number of pans that the oven can hold , in which condition the shape of pan is sure. Above all, we can get the function of n and r.n= 1.21r2(4−π)+0.085(5) 3.3.3 The optimum solutionThe weight of the number of pans which the oven can hold is P, while the weight of heat distribution is (1-P). The dimensions of the two are different. The effects are also different with the unit change of r. Depending on the message above, in order to induce the weight P. we need to eliminate their dimensions[7].Through observing the figure 23, we can know that the range of T max is 27℃with the domain of r.Then we change the value of r in its domain of definition, then we can get n's approximate range: 3.05mThen we eliminate their dimension, so they are transformed into value which could be compared.They are T max/27 and n/3.5 respectively.We hope that we can get a smaller T max and a larger n. We let T max/27 multiply by -1, then add them (T max/27 and n/3.5) together, the final result is K. K has no practical significance, and we just want to know its relative value.K=-(1-P)T max/27+P n/3.5After simplification, we can obtain:K=-(1-P)(-6.463r+10.796)+0.345Pr2(4−π)+0.085(6)How to get the pan we want with the idealized shape and its corresponding width to length ratio of the oven by using this formula?We explain it by an exampleIf some one’s ideal weight P is 0.6, the function(6) of K becomes:K=-0.4(-6.463r+10.796)+0.207r2(4−π)+0.085We can get the figure of K within the domain of r, by using the matlab. The figure is presented as follow:Figure 24(the relationship between K and r)We plugged the value of r into the equation (5), then we can get n=13.76Because the number of the pan should be an integer, we round up n, namely n=13.The largest number of pans which the oven can contain is prime number, the oven has only a width to length ratio of 1/13.So at last, we get the following conclusion:When P=0.6, n=13, W/L=1/13, r=0.058m is the best solution.To make the coefficient of oven reaches the top (namely without space), we take the radius of round square into equation (5). We should notice that n must be an integer. We adopt the method of exhaustion, and the result shows in the table below:This table will be used in the following advertizing.3.3.4 Model verificationWe can see the equation (6) , when the P tends to 0, which means the largest number of pans the oven can contain make no sense, the equation becomes: K=- (-6.463r+10.796) ; the optimum solution is r=0.164488m. which means maximize even distribution of heat for the pan is most important.On the contrary, when the P tends to 1, which means the maximize even distribution of heat for the pan make no sense, the equation becomes: K =0.345r 2(4−π)+0.085; the optimum solution is r=0m. which means the largest number of pans the oven can contain is most The value of r for thecorresponding peak valuein the figure is 0.058m。
