Hotcat Black Body Calculations

Alan Fletcher Version 141027A Oct 28, 2014

1. Introduction

The Lugano Paper, Levi et al describes a 32-day test of a new Rossi ECAT.

The temperature of the outer surface is measured using an Infrared camera. From this, the total radiated and convective power can be calculated, assuming that it acts as a Black Body.

They supplied a few photographs of the hotcat, the central part of which is an Alumina cylinder, with fins (or possibly a spiral), which apparently emits some visible radiation from the inside:

MFMP images from the paper

The paper does not identify when these photos were taken, so the claimed temperature is unknown. The experiment was in three parts : a dummy, calibration run at 450C, an active run at an initial 1260C, which was raised to a final stage at 1400C. (Temperatures as calculated from the IR camera).

The main issues raised by these pictures are that the exact emissivity (effectiveness at radiating energy) is unknown, that the "glow" is too "red" for a body at 1260C or 1400C, and that energy may be escaping by transmission, rather than by radiation.

WARNING: all these calculations are based on "representative" parameters, not on the actual ceramic used, so they should be regarded as "qualitative" and not "quantitative").

2. Black Body Calculations

For a black body the power emitted varies as the fourth power of the temperature according to the Stefan–Boltzmann law:

Stefan Bolzman

where A is the area and "e" is the emissivity (ranging from 0 to 1).

The distribution of power density by wavelength is described by Planck's law:

Black Body Power Spectrum

The peak for 1400C is at about 1.8 um.

These calculations depend critically on the emissivity, "e", which is used by the infrared camera to calculate the temperature, and then to calculate the power from the temperature.

The authors do acknowledge that "e" varies with temperature:

e with temperature

However, they do not take into account the fact that the Alumina cylinder may be transparent at some wavelengths, and that the emissivity itself depends on the wavelength.

For example, Manara et al analyzed a sample of Alumina, and found that both the emissivity and transmittance depend on wavelength :

Manara et al

In this case the black body power cannot be calculated directly from Stefan-Bolzman. Instead, it must be calculated by integrating the Planck Spectrum e(w)*PS(w) over the range of wavelengths.

The visible range of about 0.4 to 0.7 um is off these charts, but the emissivity is likely to be very low (less than 0.1), and the transmittance at 1050K (770C) appears to be converging to a constant 0.1.

In the Infrared the Optris PI 160 Thermal Imagers work in the range 7.5 to 13 um. Manara Fig 5 indicates a significant deviation of emissivity from the averages used in the Lugano report., and Fig 6 indicates that transmission falls from about 0.1 to 0 between 3.5 and 5.5um, so it is most likely not a factor in the infrared.

Manara fig 5 also shows the variation of emissivity with temperature. It is encouraging to the Lugano report that the shape of the curve does not change much with temperature, so a calibration at 400C might still be valid at higher temperatures.

The black body calculator (described below) was used to calculate the total emissivity of Manara's sample:

e lugana manara

The calculator was run for a range of temperatures, and the average emissivity was calculated as the fraction of Planck power to Stefan-Boltzmann power. The overall shape of the Lugano curve, other than the peak at 300C, is very similar to the Manara curve, so it could also be the result of a wavelength dependency.

3. Sample Black Body Calculation

A good online calculator is at : Hyperphysics (R. Nave) (I used Nave's javascript functions to generate the tables in this section)

As an example, consider a black body at 777C/1050K (Chosen because we have examples of emissivity from both the Lugano report Plot 1 and from Manara Fig 5).

Temperature 776.85C 1050K
Area 125 cm2
Emmissivity 0.45

From these we can calculate the power radiated from the Stefan-Boltzmann formula :

Stefan-Boltzmann power 387.69103627734376 watts

However, if the emissivity varies by wavelength then the power must be calculated by integrating the Planck formula :

power = INTEGRATE( e(w)*PS(w) | for w = 0 to infinity) [ other terms omitted for clarity]

e(w) is the emissivity as a function of wavelength, and PS(w) is the power spectrum as a function of wavelength.

The integral can be approximated as a sum :

power = SUM( dw*e(w)*PS(w) | for w = 10 to 10000 by dw ) [ other terms omitted for clarity]

Using this sum of Planck's formula we can see how the power varies with wavelength :

Results for constant emissivity
StepW-length nmEmissivityPowerPercent of SB Power
01000.453.3701191393008734e-118.692796128745109e-12
13500.450.0015831370241682440.00040835017476022076
27000.450.210410591489522680.054272751186076065
310000.4531.6497157909602978.16364394051116
420000.4586.7086125654757822.36539007918841
530000.4581.1508083483559420.931824766333477
640000.4557.38166164107009614.800873961919716
750000.4538.109765676102589.829932113478085
860000.4525.321636537595356.531395923087819
970000.4517.1573850615948024.425530501386385
1080000.4511.9141678580303573.073108930356409
1190000.458.4790969785672722.1870758374979693
12100000.456.1745207653573741.592639547368923
13110000.454.5910995897219891.1842160793311816
14120000.453.47839335932283470.8972075786747067
15130000.452.68007505316263340.6912914672717323
16140000.452.09637707163305050.5407339544815678
17150000.451.6622045823643250.4287446514949158
18160000.451.3341779629459180.3441343332971678
19170000.451.08280993302646780.27929712882292557
20180000.456.1872878274826351.5959326495896633
21500000.450.320928949892805740.08277956410197262
22100000Total =387.69271928120986100.0004341095637
Peak at 2760

(Note : the total doesn't quite add up to 100% because of the simple integration algorithm used).

Note : at this temperature (and up to 1400C) the contribution of Black Body power in visible wavelengths is miniscule. So even if the Alumina were completely transparent at visible wavelengths it would not effect the calculation of Black Body power --- and the actual power would thus be larger than the reported power.

However, IF the ceramic had an Emmissivity which varies by wavelength, then (using Manara's data) we would have a different picture:

Results for Manara fig 5 emissivity (stepped)
StepW-length nmEmissivityPowerPercent of SB Power
01000.017.489153642890832e-131.9317324730544692e-13
13500.010.0000351808227592943140.000009074448328004905
27000.010.0046757909219893940.0012060611374683573
310000.010.70332701757689590.181414309789137
420000.011.92685805701057360.4970086684264093
530000.011.80335129663013170.46515166147407716
640000.056.375740182341121.6445415513244126
750000.216.9376736338233664.36885871710137
860000.6536.575697220971069.434238555571294
970000.9235.07732057037169.047751247278832
108000126.4759285734007656.829130956347568
119000118.84243773014954.860168527773267
1210000113.7211572563497293.539198994153164
13110000.44.0809774130862131.0526365149610504
14120000.53.86488151035870780.9968973096385636
15130000.52.97786117018070540.768101630301925
16140000.31.39758471442203370.36048930298771187
17150000.31.10813638824288360.2858297676632773
18160000.20.59296798353151910.15294859257651905
19170000.40.96249771824574880.24826411450926705
20180000.45.4998114022067931.418606799635258
21500000.40.28527017768249310.07358183475730877
22100000Total =179.2141909883273546.22603419185641
Peak at 2760

In other words, the total power actually emitted MIGHT be significantly different from the nominal value.

WARNING : these calculations are for example only, and do NOT apply directly to the hotcat. Do NOT quote this result as "proof" that the hotcat didn't produce the reported power.

4. IR Methodology

A lot of the talk on the intertubes revolves around the general use of an IR camera for the temperature measurement and power calculations.

In my opinion they might have been better off using a "fatcat" (as is used in the 1MW system) with water-only calorimetry.

McKubre (one of the best experimentalists in the trade) appears to have no fundamental objection to IR calorimetry, though he is very concerned with the lack of calibration at the operating temperature, and thus appears to rate the whole experiment as inconclusive.

Is there evidence of excess heat? My impression is “Yes” (but see below). Is this evidence unambiguous? Not as presented. ..

The calibration that was performed was at a single point well below the points of ultimate operation. This is bad and bad. We need to determine the shape of the performance manifold empirically, ideally bracketing all points of operation. This was not done. ...

The inability or unwillingness to do one or all of these things [calibration] is a glaring deficiency in the thermal output record. That is not to say that this record is incorrect, but we cannot be certain that it is correct and the measurement uncertainty is undefined. ..

On the whole I am encouraged. Considerably more work is obviously needed to validate the adopted mode of calorimetry and support better sampling and testing. But we are given something we can sink our teeth into both experimentally and theoretically: testable fuel to products nuclear burn at temperatures that have practical, economic and social potential.

However Gordon Docherty makes a strong case for the use of IR calorimetry:

October 9, 2014 at 11:34

In reply to “the Pathoskeptic camp”, and their complaints about using thermal imaging that I have seen elsewhere.

If you used any form of direct contact, it would immediately be claimed:

1. the direct contact was used to pump energy into the system
2. the direct contacts were measuring hotspots, and were not at all representative
3. loose contacts were giving false readings or, if in direct contact, the characteristics of the measurement devices moved “out of band” as the high heat invalidated use of device
4. results from the measurement system cannot be trusted, as the system’s characteristics change as the measurement system heats up
5. the measurements were being fiddled anyway (for example, by hidden equipment somewhere along the wires)

likewise, if a liquid was used, and measured, it would immediately be claimed:

1. all of the above for direct contact plus:
2. the assumptions about flow rate are wrong
3. not enough steam (remember that one!)
4. the probes measuring the liquid (or steam) temperature were in the wrong place
5. the liquid contained impurities
6. need to run the e-Cat at very low temperatures, to raise water temperature from ambient to, say, 60 degrees – at the same time, of course, claiming the ambient temperature was incorrectly measured, so invalidating the claimed (small) difference between ambient and system temperatures.

Thermal imaging, while not “perfect” (nothing is):

1. cannot be used to pump energy into the system
2. provides a view across the whole piece (and, indeed, tends to underestimate, not overestimate), so avoiding the charge of measuring hotspots
3. by (calibrating) and using more than one thermal camera, “single points of failure” in the measurement system are ruled out.
4. cameras are physically decoupled from the system being measured, so impacting least on the system being measured
5. with (known temperature) calibration points and one of the cameras being mobile, measurement error (uncertainty of measurements) can be much reduced, even when the background temperature fluctuates
6. means the measuring equipment remains cool (again, decoupled) – not the case otherwise
7. records raw data – that means any measurement adjustments made after the fact are transparent and well understood. Ultimately, there is always the raw data to go back to.

So, thermal imaging was definitely THE BEST CHOICE for measuring the heat being put out by the system (running at 1400 degrees) – with a (very conservative) COP of 3.2 – 3.6, “rounding errors” can also safely be discounted.

Also, the selected Optris camera may not be the best for this task. An interview with of Richard Skogfeldt of Williamson indicates that their spot pyrometer (which can also be used to scan an area) would not only have given a more accurate reading, but would have calculated the emissivity in the process. Skogfeldt points out that measuring Alumina is difficult, because the reflectivity, transmission and emissivity all vary by temperature and wavelength, and are different for different samples.

Unfortunately at the critical point they were talking over each other. My interpretation is that it would have been better to do the measurement at 2.5u, near the peak of the black-body radiation, and not in the 7-15u band used by the Optris.)

5. Other deficiencies

The Lugano team say that they measured all the electrical inputs correctly, using appropriate meters.

However, Ethan Siegel still includes a "wiring fake" as an explanation. (Though that's for a single-phase connection : what would a 3-phase fake look like?) So : who provided the cord between the building supply and the "upstream" meter? Rossi (possible fake), or the experimental team? If Rossi, they could have cut the cable at the end of the experiment (as Defkalion offered to do.) Since this fake was identified in 2011, why didn't they address it in 2011, 2013 or 2014?

They include visible-range photos, but do not indicate at what points of the experiment they were taken.

They provide a single IR image, during the calibration phase, but none from the active stages.

The published report was obviously cut down to a minimum length (45 pages) in the hopes of publication.

The situation could possibly be corrected by posting a long version, including a second-by-second account of Rossi's involvement, and details of the calorimetry at identified stages of the experiment.

6. Conclusions

Even if the Alumina were completely transparent at visible wavelengths it would not significantly effect the calculation of the temperature or of the Stefan-Bolzman Black Body power --- and the actual power would thus be larger than the reported power. Since the visual photographs of the hotcat were at unknown times, with an uncallibrated camera (white balance), discussions of "color temperature" cannot be used to deduce anything.

The power calculation of the hotcat depend strongly on the emissivity used:

In short, because the actual emissivity of the Alumina tube may differ greatly from textbook examples, the output power may not be calculated correctly.

The experiment was calibrated at 400C, but this is too far from the operating point, and they were unable to mount a thermocouple on the cylinder to get an independent reading.

Samples of the cylinder were sent for molecular analysis, confirming that it is Alumina. A sample should have been sent for calibration of its emissivity.

I therefore rate the Legano Report results as "inconclusive".

IF they had simply filed off a small flat area on the cylinder and correctly mounted a thermocouple, then they could have calibrated the whole experiment for temperature, and the results would have been more conclusive. Whether they asked Rossi to do this and he declined, we don't know.

For want of a nail the shoe was lost;
For want of a shoe the horse was lost;
For want of a horse the battle was lost;
For the failure of battle the kingdom was lost—

All for the want of a horse-shoe nail.

7. Notes

Wikipedia : "The term emissivity is generally used to describe a simple, homogeneous surface such as silver. Similar terms, emittance and thermal emittance, are used to describe thermal radiation measurements on complex surfaces such as insulation products." They are treated as synonyms here.

In another paper Manara reports that existing IR systems are inadequate for high temperature ceramic systems, which may be semi-transparent. They present results from a new system which analyses data from 2um to 20um and can be used up to 1600K. This is, of course, the operating range of the hotcat.

This document may be quoted, copied or translated, with the following copyright notice:

© 2014 Alan Fletcher lenr.qumbu.com

Images may only be copied WITH the annotation of where they came from.