First Multi-Mirror Solar Oven Simulations

by Curlydock
Nov 27, 2007

nov1039a.jpg

Results published in previous posts increased confidence in my solar oven simulation program. Those tests were confined to the single-mirror configuration. By keeping things simple, results were intuitively judged. Now it gets more interesting as we encounter kaleidoscopic multi-mirror configurations with results that may challenge intuition but that are, nevertheless, we hope, still accurate.

The number of variables that define the problem make closed form expression prohibitive, at least for me. Thus, my use of simulation.

HERE IS A LIST OF THOSE VARIABLES:

The oven seems simple enough in that it is composed of two sets of hinged mirrors and a sphere called a “bounding sphere” that stands for the oven cavity. The oven cavity is the part of the oven that gets hot enough to cook food. It is constructed like a small green house, but we need not concern ourself with those details.

The bounding sphere is all we need to know about the oven cavity. The center of the sphere gives the CAVITY LOCATION in EXTENSION and ELEVATION from the “origin”, the lowest point on the vertical R1-R2 hinge. The cavity location is always in the bisecting plane of the oven’s bilateral symmetry.

The radius of that sphere gives the SIZE OF THE CAVITY. Now, to make results easier to compare across simulations, I have normalized all units of distance to the diameter of the oven cavity. So, instead of inches, feet, meters, etc, all lengths are in units of cavity diameters.

I have named the four mirrors of concern: R1, R2, R3 and R4. Hinged on a vertical axis are R1 and R2. R3 and R4 are hinged horizontally. So, we have the two angles: R1–R2 ANGLE and R3–R4 ANGLE.

We will also be concerned with where the sun is in respect to our oven, so we have: SOLAR ALTITUDE ANGLE or elevation; zero degrees on horizon and 90.0 degrees at zenith, and SOLAR AZIMUTH ANGLE. In all runs I presently contemplate, the azimuth will not change, confining the sun to the plane that bisects the symmetry of the oven.

We also have to be concerned with R1 WIDTH and R1 HEIGHT. R2 is always identical to R1 in width and height.

Since R3 and R4 are hinged horizontally, it is confusing to speak of their width and height, so I define the dimension R4 RATIO, the length of the horizontal hinge, as a portion of the opening of R1-R2.

R3 is always parallel to, or in the plane of, the horizon. R1-R2 sits on R3 like an open book with the bottom two corners on the axis of the R3-R4 hinge. Because of this configuration, the reflective part of R3 will always be a triangle with size and shape determined by the width of R1, R2 and the R1-R2 angle. The other mirrors are always square or rectangular.

R4 EXTENSION is how far R4 extends from the horizontal hinge.

Also an input variable is the QUANTITY OF REFLECTIVE PLANES. I can consider four possibilities: R3 only; R1 and R2; R1, R2 and R3; as well as all four mirrors.

I think that about covers the input variables. So far I have not concerned myself with losses based on reflective angles, surface imperfections, convection, conduction, mass, re-radiation, transmission distances, imperfect insulation, etc.

OUTPUT

What we get for output is the amount of idealized SOLAR FLUX GAIN on the oven cavity over what the cavity would have gotten directly from the sun without any concentration. This is not a temperature; it is a ratio. The minimum value will always be 1.0 unless the cavity is in shadow or blocked from direct rays, where the value could fall to zero. The maximum value depends on all the inputs.

Because of the large quantity of possible combinations of input variables, I must find a way to lock some of them down so they can be temporarily ignored. That will simplify the problem and help us make sense of where we are in our journey to the “perfect” design. I do not know if there can ever be a proof that a particular design is the best possible one. There might always be a combination of inputs hitherto unconsidered that will yield a greater flux gain than what we thought was the best. If there is such a proof, a bigger brain than my own will need to find it.

That is why I am willing to share my program listing with anyone who is interested and who agrees to keep my program or any derived from it in the public domain, not for profit.

I am also willing to do runs for others who would like to see particular results of inputs meaningful to them. Perhaps they have built an oven or two and want to see what the simulation says about comparative performance. If you are interested, just leave your request in the comments section and I may publish my response or run results in a future post.

Next, some multi-mirror runs:

For the time being I will keep R1-R2 angle at 60.0 degrees. The decision is arbitrarily based on the fact that an oven I built uses that angle. In the long run, there may be a better angle.

Also, for now I keep the width of R1 at 4.0 cavity diameters. I do that because prior runs with R3-only seemed to result in 4.0 as the optimum size for the side of R3, given the R1-R2 angle = 60.0 degrees. The length of a side of R3 is the same as the width of R1 (and R2). Lastly, to start, I make the height of R1 equal to its width.

nov1029a.jpg

Chart “nov1029a” shows the effect of cavity extension on flux gain where only R1 and R2 are in play. For these inputs, the best value seems to be 1.4, but it depends on the angle of the sun. By this chart, for sun angles higher than about 40.0 degrees we might want to make the cavity extension less than 1.4 cavity diameters.

Next, keeping cavity extension at 1.4 and everything else the same, we do four runs, one for each of four different heights of R1-R2. That results in chart “nov1030” where
red (1) = 4.0
green (2) = 6.0
blue (3) = 8.0
black (4) = 10.0
cavity diameters.

nov1030a.jpg

We see gains increasing at higher sun angles as the height of R1-R2 increases. The gain increment decreases so returns are diminishing. At this point I will make the compromise of R1-R2 height = 6.0 and keep that for the next set of runs.

I cannot resist the temptation to add more mirrors at this point. Chart “nov1032” run (1) is a repeat of “nov1030” run (2), for reference. It looks different because the vertical axis of the chart is re-calibrated to cover higher gains.

nov1032a.jpg

Note the gain still peaks at about 6.0. That gain makes sense when you consider that two mirrors at a 60.0 degree angle would allow the sun to see six images of the cavity. Run (1) is only with R1 and R2 in play.

For run (2) I merely add reflector R3 to the play, and the results show a significant increase in gain. If we consider that R3 would double 6 images to 12, the gain seen of about 9.0 makes sense, considering some of the images will be partially obscured.

Now we take a lesson from a prior post when we were testing with only R3 and learned that it is best to elevate the cavity somewhat over R3. We do just that in run (3) and see that the gain is even closer to 12.0. Supposedly, the reason is that the images are less obscured when the cavity is elevated.

In run (4) I put R4 in play and made R4 the same size and shape of R1 and R2. It is a handy size if the oven is collapsible and you stack the reflectors for storage. The solar flux gain reached 14.0. One might worry that the simulation is malfunctioning because the flux gain is below 1.0 for very low angles of the sun. More thought reveals that at these low solar angles the cavity and lower part of the oven are in the shadow of R4, so one’s confidence in the program rebounds.

Is this the best design? I doubt it. Only more testing can tell, I guess. At least, it is a benchmark for measuring more attempts.

In a future post I can make the input variables conform to the oven I have actually constructed and see what the simulation says. I will be able to compare a real oven with this benchmark and tell what changes might make it better or worse without getting my hands dirty (not meaning to disparage getting one’s hands dirty).