
At the design stage, calculation is the only way to evaluate the quality of the blackbody, to determine its optimal parameters, and to choose suitable material. But if you bought a blackbody calibration source and intend to use the value of effective emissivity specified by the manufacturer, you should be doubly careful because:
1) Values of effective emissivity might be correct only for a part of operational spectral range.
2) None of blackbody manufacturers never purchased our blackbody emissivity modeling software. Other analogous software simply does not exist.
3) If effective emissivities were measured, you may not know which measured equipment was used and what is the precision of such measurements.
Measurements performed by independent researchers are not always confirm manufacturer's values. The picture at left shows how manufacturer's emissivity of a flat-plate blackbody agrees with that measured at NIST (adopted from: S. N. Mekhontsev, V. B. Khromchenko, L. M. Hanssen, "NIST Radiance Temperature and Infrared Spectral Radiance Scales at Near-Ambient Temperatures," Int. J. Thermophys. 29, 1026-1040 (2008) - Fig. 14. Effective spectral directional emissivity of a high-temperature flat-plate BB. Manufacturer specification is 0.95 (+0.00, -0.05)).
If such disagreement takes place for flat-plate blackbody (metallic plate coated with black paint) for which emissivity measurement is not too difficult, one can imagine the situation with the cavity radiators.
It is impossible to calibrate radiometers, radiation thermometers, optical radiation detectors etc. without knowledge of radiation characteristics (primarily, effective emissivity) of blackbody calibration sources. Frequently, uncertainty in determination of effective emissivity is the predominant component in the uncertainty budget. Obviously, the best way to determine effective emissivity is measurements.
According to its definition, direct measurement of effective emissivity should be performed by comparison of radiation emitted by the blackbody source under consideration and radiation emitted by a perfect blackbody having the same temperature. However, the perfect blackbody is the physical abstraction that doesn't exist in the real world. There are two methods allowing to resolve this problem:
1) To perform absolute measurements of artificial blackbody radiation characteristics (e. g., spectral radiance in W·m-3·sr-1), assign a certain temperature to this radiator (real-world blackbody sources are always nonisothermal), then divide its measured quantity by homonymous quantity calculated for the perfect blackbody at the same temperature (e. g., using Planck's law for the case of spectral radiance).
2) To perform relative measurements by comparison of radiation emitted by the artificial blackbody under consideration and the better-quality artificial blackbody (such as fixed-point blackbodies).
Realization of both these methods is extremely complicated; it requires special instrumentation and highest qualification of operating personnel. Presently, such measurements can be conducted only in several world's leading metrological centers.
Another approach implies indirect measurement when cavity reflectance is measured instead of its emissivity. Indirect measurements can be implemented easily, for instance, using integrating sphere with the laser radiation source. However, this approach has serious disadvantages:
1) Usually, cavity reflectance is measured for the cavity temperature that significantly differs from the temperature at which the blackbody operates.
2) Often, it is difficult to provide geometry of cavity irradiation and collecting reflected radiation which correspond to geometrical conditions of collecting radiation emitted by the cavity when it is used for calibration. This violates the reciprocity principle and raises a query about correctness of effective emissivity determination.
3) Sometimes, such measurements are impossible at all, for example, due to the difference in dimensions of the blackbody and reflectometric device.
All serious researchers compute effective emissivities using Virial's software. The list of some publications concerned with the use of our old STEEP3 program can be found here. We invite our visitors to explore our programs of STEEP3 family (STEEP320, STEEP321, STEEP322, STEEP323) and INCA333 program for a cylindrical cavity with an inclined bottom in order to choose software the most suitable for solving your tasks. The Comparative Table will help you to make right choice.
Cavity of various geometrical shapes are used as blackbody radiators. The choice of a cavity shape is determined by reflection properties of cavity wall material, conditions of cavity radiation collecting, and convenience of cavity manufacturing and treatment of its internal surface. The most frequently, cavities with axial symmetry are employed. All programs of STEEP3 family allow to model cavities formed by rotation of a non-self-intersecting polygonal line around an axis. If polygonal line has a lot of segments it can reproduce complicated surfaces of revolutions. Pictures below show some cavity shapes which can be modeled using STEEP320, STEEP321, STEEP322, and STEEP323. Figures in parentheses denote the numer of segments in the cavity generatrix.
Cone with diaphragm (2) Cylinder with diaphragm (3)
Cylindro-cone with diaphragm (3) Cone-cylindro-cone (3)
Cylinder-inner-cone with diaphragm (3) Cavity with reflector (9)
Cylinder with V-grooved bottom and diaphragm (23) Sphere (399)
Among non-axisymmetrical cavities, the cylindrical cavity with an inclined bottom is the most frequently used due to simplicity of fabrication and treatment. INCA333 allows to model just this type of cavities. Pictures below show two orthogonal cross-sections of a cylindrical cavity with an inclined bottom and flat annular diaphragm.
Effective emissivities and radiance temperatures derived from them depend on viewing conditions (i.e., geometrical conditions of collecting the radiation by a measurement device). Picture below shows the most common viewing conditions used in optical radiometry, radiation thermometry, and adjacent areas.
Normal viewing conditions might be considered as a special case of directional viewing conditions. Both these types correspond to the collimated beams. Such viewing conditions are approximately implemented for the radiometric or pyrometric measurement scheme that uses very long focal-length optics. If viewing beam radius is equal to zero (hypothetical case of infinitely thin viewing beam), distribution of the normal effective emissivity across the aperture of a cavity characterizes the spatial nonuniformity of the radiance emitted by a cavity; dependence of the directional effective emissivity upon the angle of entry of a viewing beam into cavity characterizes the angular nonuniformity of the radiance emerging by the given point of cavity aperture. For right conical viewing conditions, the beam axis is parallel the cavity axis. Integrated viewing conditions reproduce collecting of radiation emerging from a cavity by a circular detector without optical system. Hemispherical viewing conditions characterize total radiation losses of a cavity and can be considered as a special case of the integrated viewing conditions when the detector radius is equal to the cavity aperture radius and the distance between cavity aperture and detector is equal to zero.
Normal Directional
Right conical (convergent viewing beam) Right conical (divergent viewing beam)
Oblique conical (convergent viewing beam) Oblique conical (divergent viewing beam)
Integrated (detector is shown in blue) Hemispherical
The temperature distribution over the radiation surface of a blackbody cavity is the crucial factor together with the optical properties of cavity walls determining the effective emissivity and, consequently, the radiance temperature of a blackbody. High-quality blackbody radiator must be almost isothermal. Fixed-point and heat-pipe blackbodies as well as cavities immersed into thermostated liquid have the highest degree of temparature uniformity. However, real-world cavities are always nonisothermal. Even if the external surface of a cavity has ideally uniform temperature, the finite thermal resistance of a cavity wall (or layers forming it) and different radiative heat losses from the aperture for different parts of the cavity internal surface unavoidably lead to nonisothermality of radiating surface. Evaluation of this effect is possible currently only by computational methods.
Temperature nonuniformity of radiating surface affects significantly on effective emissivity of a blackbody, escpecially if the area within FOV (field-of-view) of the detector, radiometer, or radiation thermometer is nonisothermal. Nonuniformity of radiating surface leads to different contributions of radiation emitted by different elements of cavity walls at the same wavelengths to the resulting cavity radiation. This, in turn, makes the cavity spectrum non-Planckian. Effective emissivity is the ratio of cavity's radiation characteristic (e.g., spectral radiance) to the analogous characteristic of a perfect blackbody at the same wavelength and the same temperature. To avoid ambiguousness in determination of temperature for nonisothermal cavities, the reference temperature Tref is used in definition of effective emissivity. Depending of Tref choice, effective emissivity of a nonisothermal cavity might be less or greater than that of an isothermal cavity, and even be greater than unity. It is recommended to choose Tref between minimal and maximal temperature of cavity's radiation surface. Very often, the temperature of cavity bottom center is a good choice for Tref. Note, that the choice of Tref do affects only effective emissivity (it is fictitious quantity in a certain sense) but does not affect on measurable quantities (e.g., radiance, radiance temperature, radiation flux etc.)




Temperature nonuniformity of the area of cavity's internal surface within FOV (that is directly "viewed" by measuring device) results in drastic changes of spectral effective emissivity. Right graph shows the spectral effective emissivity plotted vs. wavelength for the same cavity as in previous example. Tref is also equal to 1300 K for two temperature distributions:
Temperature distribution 1. Temperature linearly decreases from 1300 K at the bottom center down to 1290 K at its periphey, then decreases linearly along cylindrical wall down to 1287 K; cavity diaphragm has constant temperature of 1287 K.
Temperature distribution 2. Temperature linearly increases from 1300 K at the bottom center up to 1310 K at its periphey, then increases linearly along cylindrical wall up to 1313 K; cavity diaphragm has constant temperature of 1313 K.
Temperature of isothermal cavities doesn't affect effective emissivity only if cavity's surround is the non-emitting media (like vacuum at 0 K) of if influence of background is negligible. Since the actual characteristics of background radiation is very seldom known, usually, one can consider cavity's surround as a perfect blackbody radiator with the "background temperature" Tbg.
The influence of Tbg = 300 K on the spectral normal effective emissivity of the above-mentioned cavity having isothermal walls at T = Tref = 400...1000 K is shown at left. Black line corresponds to the case of the isothermal cavity at Tbg = 0 K.
The radiance temperature measured in conditions of non-zero background will differ from the termodynamic temperature of a cavity not only due to non-unity effective emissivity but also due to contribution of the background radiation partially reflected by the cavity.
For nonisothermal cavities, both effects (temperature nonuniformity of a cavity and radiating background) act conjointly. Besides, temperature effects cannot be considered in isolation from the angular distribution of radiation reflected by the cavity walls. In terms of the Monte Carlo ray tracing, this can be explained by the following fact: a ray undergoing multiple reflections from the cavity walls "accumulates" spectral radiances of all reflection ponts at their temperatures; for different angular distributions of reflected radiation, these points will be distributed alse differently.
All programs of STEEP3 family allow modeling of nonisothermal blackbodies with axially-symmetric temperature distributions. Temperatures can be defined over the discrete point set on the cavity generatrix; linear interpolation is used during ray tracing.
INCA333 allows to model only one-dimensional temperature distributions, along the axis of the cavity cylindrical part; the diaphragm is supposed to be isothermal at the temperature of the adjacent edge of the cylindrical part.
In exemplification of this phenomenon, the dependences of the normal spectral effective emissivity on wavelength for gray diffuse cylindrical cavity with diaphragm having wall emissivity of 0.9 are shown below. Cavty bottom is isothermal at 1300 K; three temperature distributions are considered: (Isothermal) all cavity walls have uniform temperature of 1300 K; (10% Decrease) linear decrease of temperature along cylindrical surface from 1300 K down to 1287 K, diaphragm has uniform temperature of 1287 K; (10% Increase) linear increase of temperature along cylindrical surface from 1300 K up to 1313 K, diaphragm has uniform temperature of 1313 K. In all cases, Tref = 1300 K.

Precision of the Monte Carlo calculations is in inverse proportion to the square root of a number of rays traced; usually, it is enough to trace from 105 to 107 rays to achieve the accuracy of 10-4…10-6 in effective emissivity.
The logarithmic plot at right illustrates the stochastic convergence of the Monte Carlo modeling for several calculation processes peformed for a diffuse spherical cavity, for which the effective emissivity can be calculated exactly, using well-known analytical formula.
Monte Carlo calculation can be considered as a numerical experiment in some a measure similar to the ordinary natural experiment. Therefore, its accuracy can be assessed by performing repeated calculations and consecutive statistical processing of numerical experiment results.
Adequacy of the Monte Carlo modeling depends equally to a great extent upon two factors. The first is the temperature distribution over the blackbody walls, especially over their directly viewable areas. One can evaluate the effect of temperature non-uniformity on the radiation characteristics of a blackbody by performing calculations with several feasible temperature distributions.
Direct measurements of effective emissivity of blackbody radiation sources are often extremely difficult or even impossible. Sometimes, computational methods are the only way to determine effective emissivity. Moreover, calculation of effective emissivities should be done at the design stage. Many computational methods for effective emissivities have been developed since the second half of XX century (see their overview here). They are based on the various physical and mathematical assumptions, have different areas of applicability and provide different degrees of accuracy. At present, the most comprehensive and flexible method for calculating radiation characteristics of blackbody radiators is the Monte Carlo method. It is based on the ray tracing algorithm that models radiation heat exchange among cavity walls and propagation of radiation from a blackbody to a radiation detector. The Monte Carlo method uses the ray (geometrical) optics approximation; such phenomena as polarization and diffraction are not considered.
The computational Monte Carlo algorithm employed in Virials's blackbody emissivity modeling software is based on the optical reciprocity principle and the technique of backward ray tracing. The history of a ray directed from the point of observation into a cavity is being traced until it leaves the cavity after reflections from the walls, or until its energy becomes less then the given value (flux threshold). The last point of reflection is considered as a birth point of a ray propagating in opposite direction. The spectral radiance acquired by a ray consists of the spectral radiance of own thermal radiation (computed by Planck’s law) in the point of last reflection in the viewing direction and the sum of spectral radiances in the points of successive reflectance computed along the ray trajectory and impaired by multiple reflections. Averaging of spectral radiances of a large number of rays traced allows to estimate the effective emissivity and radiance temperature of a blackbody.
Second factor crucially affected effective emissivity is the adequacy of angular dependences of the model adopted for optical properties of materials forming a blackbody cavity. The simplest (but very often insufficient) model is the diffuse model of reflection. The most powerful models must take into account bidirectional reflectance distribution function (BRDF) that describes angular distributions of reflected radiation for every direction of incident radiation. Since the measurement data required for such a model is very often incomplete or absent at all, the specular-diffuse model of reflection representing BRDF as a sum of the Lambertian (diffuse) and the perfect specular components was in common use up to date. However, BRDF for specular-diffuse model of reflection includes Dirac’s delta-function in the expression for specular component and cannot be fitted to measured BRDFs. Separation of reflectance into the perfect diffuse (Lambertian) and perfect specular components cannot be done unambiguously. Besides, specular components of reflection for real-world materials of blackbody cavities have small but finite divergence what can lead to effects which cannot be predicted within framework of specular-diffuse model of reflection.
Recently, we developed so-called Three-Component (3C) BRDF model (click here for detail) as well as software tools for fitting model’s parameters to in-plane BRDFs (i.e., BRDFs measured in the plane of incidence) and for separation of the spectral directional-hemispherical reflectance (spectral DHR) into three components - perfectly diffuse (Lambertian, or simply diffuse), glossy (having wide specular lobe), and quasi-specular (having narrow specular lobe; simply specular, for brief). For the moment, the 3C BRDF model is implemented in STEEP323 and INCA333 blackbody emissivity modeling programs.
