Resume of form factors for curved surfaces.
1. Introduction
Curved surfaces have not been thoroughly considered in radiative transfer analysis mainly due to the difficulties arising from the integration process and perhaps because of the lack of spatial vision of researchers. When dealing with them, application of the iterative method or direct calculation through integration does not provide with an exact solution, so that only approximate expressions or tables are given for a very limited number of forms [1]. In this way, a vast repertoire of significant shapes remains neglected and energy waste is evident. For this reason, further research on the matter, starting from a different approach was considered worth doing.
In previous researches from the authors, form factor calculation has been undertaken for several types of emitters. In all cases, geometric properties of those, revealed as the most powerful tool that shapes radiant interchange [3,4,5,6]. This included mainly rectangular shapes, plane forms and the volumes that can be composed with such primary geometries.
Following the same approach to radiative transfer through the basic understanding of the spatial and geometric properties of volumes, in this chapter new form factors derived from a combination of curvilinear surfaces are hereby presented. Starting from the properties of the sphere and with simple calculus, new laws are devised, which enable the authors to discover a set of configuration factors for caps and various segments of the sphere. The procedure is subsequently extended to the paraboloid, the ellipsoid or the cone, useful in issues such as rocket nozzle design and organic shapes contained in human physique. Appropriate combination of the said forms with truncated cones, produces highly articulate shapes, which frequently occur in the technical domains but were not feasible for exact calculation during a number of years. The research is duly accomplished by presenting the equations needed to evaluate interreflections in curvilinear geometries. Thus, heat transfer simulation is enhanced by such results leading to create innovative software which has been expanded in turn by the authors.
2. Outline of the problem
The reciprocity principle enunciated by Lambert in 1760 and expressed in Eqn. (1), yields the following wellknown integral equation (2) that acts as the theoretical basis for form factor calculation between two surfaces.
Where the terms are depicted in Figure 1,
E_{bi}= radiant power emitted by the corresponding surface 1 or 2
A_{i}= area of surface, dA_{i}= differential of area
r = distance radiovector
θ_{i} =angle between radiovector at differential element
The previous expression states that radiant interchange for every given form depends on its shape and its relative position in the threedimensional space (Figure 1). From the times of Lambert to our days, researchers and scientists in the fields of geometric optics and radiative transfer have sought to provide solutions to the canonical equation (2) for a variety of forms [1]. This is no minor feat, since the said equation leads in most cases to a quadruple integration and the fourth degree primitive of even simple mathematical expressions often implies lengthy calculations.
Given the fact that this equation depends on geometric parameters, it is reasonable to think that there should be an easier way to approach the problem rather than dealing directly with the integral; also, with the aid of computer simulation, mathematical solutions of complex functions can be approached in a simple and friendly way. Curvilinear forms present some characteristics that make them suitable for a different treatment in terms of radiative transfer.
3. Form factors derived from the sphere
Starting from simple forms several form factors can be calculated without hardly any calculus; later, this logic can be applied to more complex configurations. Let us consider first the simplest form, a sphere that irradiates energy from its inner surface; the irradiated energy is entirely received by itself; so that, being the sphere surface 1, the only factor that has to be considered is:
Bearing this in mind, in a similar surface, for instance a hemisphere, the form factor is accordingly F_{11}= ½. The configuration factor of a differential area to a disk of radius
Stimulated by this result, volumes composed of only two surfaces, one being planar and the other spherical, were analyzed. The first case was the spherical cap which is a generalization of the hemisphere.
Extending the reciprocity principle to a spherical cap (Fig. 2) of radius R (surface 1), and its entire base (surface 2) the factor was obtained from the relation
Two important laws are inferred from here, which have been defined as CabezaLainez laws:
CabezaLainez first law:
If a volume is encircled by two surfaces preseting one of them positive of thempositive curvature, and the second being planar, the exchange factor from the curved surfaceto the other equals the inverse ratio of areas of the aforementioned figures. The notion of positive curvature of the element is introduced to foresee stagnation of radiant flux.
CabezaLainez second law:
Within a spherical surface the form factor of any given area over itself is precisely the fraction between that area and the sphere
The second law requires of more deduction as follows
Given that a spherical cap represents an Y^{th} fraction of the total area of the sphere of radius R, and recalling from trigonometry that,
Thus,
Consequently,
CabezaLainez second law:
The configuration factor of an Y^{th} part of the sphere over itself is precisely the inverse of Y.
Thus, the assumption for the hemisphere is confirmed; in the quarter of sphere F_{11} has to be 1/4 and successively for every portion of the given sphere.
This law will hold true even if we are not dealing with spherical caps but for any fragment of the surface. Taking a critical look at the canonical equation (1) adapted to the sphere, it is logical to establish a relationship between r,
Substituting, these terms in the canonical equation (1):
This means that spherical surfaces present these unique properties (Eqs. 3 and 8) which are crucial for our discussion crucial for our discussion.
Now CabezaLainez laws can be applied to more complex volumes that involve portions of the sphere. Considering a sector of the sphere comprised between to semicircles forming an internal angle
As has been discussed, the Y portion of the sphere is, in this case
Accordingly,
And introducing the areas of the semicircles,
Following the discussion, these pair of semicircles can form any angle x between 0 and 360 degrees (Fig. 5). So that, the following equation, which has not been found expressed previously in the literature, is proposed in order to obtain the energy balance between the half disks, where x represents the value of their internal angle (Figure 5).
The latter expression (Eq. 14) is a good indicator of the factor between two inclined and equal surfaces with a common edge. If they are not too dissimilar from the semicircle, a factor that is usually lengthy and cumbersome to calculate can be devised easily.
Let us now return to the first principle, the expression
If we introduce at this point the dimensionless parameter β, we can simplify equation 16 as,
Since this principle is more general than the second one, we can extend it to nonspherical surfaces.
4. Application to common surfaces
4.1. Prolate semispheroid
Surface 1 is the spheroid and surface 2 is the circular disk that works as a base to the former, h>a.
Firstly the dimensionless parameter m is introduced:
By virtue of the first principle,
And making,
4.2. Oblate semiespheroid
Surface 1 is the spheroid and surface 2 is the circular disk that works as a base to the former, h<a
Denote the parameter m_{1},
By the first principle and,
With the same procedure as before to make the expression dimensionless
4.3. Paraboloid of revolution
Surface 1 is the paraboloid and surface 2 is the circular disk that works as a base to the former
4.4. Right cone
1 is the surface of the cone and 2 is the circular base
Introducing the parameter
It is possible to compare the performance in terms of F_{12}, of all the figures found up to now, where the cone shows better performance followed by the paraboloid.
4.5. Ellipsoid
In this case, 1 is the surface of the ellipsoid and 2 is the elliptic base; y is a parameter equal to 1.6. The example shows that the first principle is not limited to surfaces of revolution.
As the area of the ellipsoid is not exact, we can expect errors on the range of 1% depending on the values of a, b and h.
This principle can be also used in other surfaces, for example, for two complementary caps within the sphere of radius r,
As an immediate consequence of CabezaLainez laws,
If now the caps within the same sphere are of any size and arbitrary position,
In this case by virtue of CabezaLainez Law,
And now we need to apply the canonical equation 9 again, substituting the respective areas of the caps;
In the special situation that the caps are parallel, which equates a truncated cone, the flux would be
Thus, the fraction is obtained as,
By virtue of equation 45 it is feasible to address radiative transfer in several figures composed of three surfaces and limited by parallel disks like truncate paraboloids, caps and especially cylinders. Appropriate equations can be easily formed in which only two values need to be found. To the circles in the extremes of the cylinder a spherical cap could be connected (fig.16) and the radiative transfer would not be altered significantly since we have previously described the performance of caps limited by circles. In the particular case that the cap is a hemisphere, the factor already determined ought to be multiplied by 0.5 and subsequently for different curvatures, bearing in mind that the unity is the circle and null would imply a “theoretical” whole sphere
Note that values under 0.5 can also be found for this relationship in a sort of globular cap with an area bigger than the hemisphere.
The space of figure 16 has been used throughout the history of buildings in cathedrals, opera houses, museums and assembly halls. If both extremes are curved, such shape is still found at bunkers, water tanks and pressure vessels of power reactors.
4.6. Two opposed spherical caps with a common axis
In order to calculate the radiative exchanges in this relatively complex figure, we need to determine beforehand the following nine geometric parameters that depend on the geometric variables shown in Figure 17.
And then we would obtain the corresponding nine form factors involved,
In this simple way the problem is completely solved
4.7. Straight cone
This is a limit case of the previous discussion.
As the former also includes the cone, by making
If
Only three factors remain,
F_{31} is obviously the ratio of areas of the cone to its base which proves that the equation is true, by virtue of CabezaLainez Law.
4.8. Paraboloid, truncated cone and spheroid
If for instance, the upper extreme of the volume is a paraboloid and the lower surface is an oblate ellipsoid (Figure 19), we can still maintain the same factors with the following simple adaptations,
as in the paraboloid alone
as it were in the oblate elipsoid alone
F_{31}, F_{32} and F_{33} have the same values as before as these correspond to the truncated cone and bear only nominal relation with the surfaces of the extremes,
Similar results will be obtained when the truncate is a paraboloid instead of a cone as it is the case in rocket nozzles.
4.9. Summary of the findings
All the aforementioned form factors have been obtained by logical deduction. In order to provide researchers and designers with all this factors in a compact format, the following table is presented, which comprises all the volume configurations presented in this chapter.
F_{21} is always the unit as shown by first law




Prolate semispheroid with circular base 



Oblate semispheroid with elliptic base 



Revolution paraboloid with circular base 



Straight cone with circular base 



Revolution Ellipsoid 



5. Interreflections amongst surfaces in a closed volume
Until this point the discussion has dealt with primary transmission of energy but, in a closed space, if the surfaces have some degree of reflectivity a significant part of the flux would be reirradiated and the concepts of emitters and receivers entwine.
Under such circumstance, the global balance of radiant power can be found through expression 71,
E_{dir} is defined as the direct power received while E_{ref} stands for the reflected energy. The two quantities added yield the global balance of radiant energy E_{tot}. If the problem entails several surfaces, expression 71 is expanded for an array of equations. To resolve it, we define beforehand the matrices F_{d} and F_{r}, whose elements are described as follows in a threedimensional fashion, (see figure 16):
Each term in equations 72 and 73 is presented in the form F_{ij} (F_{11}, F_{12}...). This stands for the configuration factors already found, from one of the surfaces
A detailed explanation for the phenomenon is given in [3]. Formerly, as volumes considered were limited by planes, all the elements in the diagonal of matrix F_{d} were equal to zero and we could not deal with the problem while, for curved surfaces, the values of the diagonal are different from null and need to be calculated with the expressions hereby presented.
Once the value of these matrices is obtained, it is easy to establish the following relationship between direct and reflected radiation:
As the value of reflected radiation is known, the problem is solved. However, we have to bear in mind that the number of surfaces should be augmented depending on the dimensions of the case study. The procedure for interreflection can be considered iterative depending on the accuracy that is required for a particular problem [3].
The simplest case of repeated reflections appears in the sphere and is wont to be employed in lieu of the former calculations with matrices.From expression 9 and successive, it was deducted that energy impinging on a point of the sphere from an emitter contained in the same surface equates the quotient between the area of the emitting surface and the total area 4 πR2, and it can be expressed under the form W/A.
After a relevant number of reflections, the total power distributed over the sphere is defined by:
As,
In the precedent discussion ρ includes the mean of all reflection quotients ρ_{i} inside the sphere, while E represents the direct power exiting from the source. Such expression would be technically applicable to all kinds of surfaces, but its accuracy dwindles when the actual volume is not akin to a sphere. If such is the case, equation 79 would be less acceptable.
Since the reflectivity of the internal surfaces can be changed on demand, the way to treat glazed elements or voids is to assign them a high absorption coefficient to ensure that those elements play a limited role in the global energy balance.
6. Conclusions
An everincreasing number of configuration factors for curved geometries, has been deducted. The authors have extracted the former in total conformity with the procedures of optical mechanics and thus the new factors can be termed as exact in contrast with other random or discretized methods.
This represents an indubitable advance of knowledge for radiative heat transfer that is already being implemented in computer models. However, the details of the simulation procedures are not discussed in this chapter in the credence that other scientists will arrive with perfect ease to the required algorithms.
Thus, this new form factors have been programmed in computer algorithms, creating a powerful tool that is able to enrich the repertoire of forms and spaces suitable for simulation. This procedure will benefit energyconscious engineering and architecture, as has been demonstrated by the authors in previous publications [7, 8,9,10] Indeed, the prototypes based on the science of heat transfer are sure to progress in their accuracy and sophistication. Radiative devices and fixtures can be conceived departing from the findings exposed previously on a more scientific basis and this will be beneficial to expand the innumerable boons of solar radiation.
Contemplating the ruins of the colossal statues of Ramses in Egypt, Shelley once wrote:
Acknowledgments
Jose Cabeza would like to thank his family in Japan and Spain for failing to understand his work.
References
 1.
John R. Howell, A Catalogue of Radiation Heat Transfer Configuration Factors. 3rd ed., 2010. Online version available at: http://www.engr.uky.edu/rtl/ Catalog/.  2.
Buschman, Albert Jr. and Pittman, Claud M., 1961, "Configuration factors for exchange of radiant energy between axisymmetrical sections of cylinders, cones, and hemispheres and their bases," NASA TN D944.  3.
CabezaLainez Jose M. Solar Radiation In buildings. Performance and Simulation procedures. InTech. 2012.  4.
Cabeza Lainez Jose M. New Configuration Factors for Curved Surfaces. Journal of Quantitative Spectroscopy and Radiative Transfer (JQSRT). Vol. 117. March 2013.  5.
Cabeza Lainez Jose M. New configuration factor between a circle, a sphere and a differential area at random positions. Journal of Quantitative Spectroscopy and Radiative Transfer (JQSRT). Vol. 133. November 2013  6.
Cabeza Lainez Jose M. Fundamentals of Luminous Radiative Transfer. Netbiblo. 256 pg. December 2010.  7.
Cabeza Lainez Jose M, Jimenez Verdejo Juan R. The Japanese Experience of Environmental Architecture through the Works of Bruno Taut and Antonin Raymond. Journal of Asian Architecture and Building Engineering (JAABE). Pp. 3340. May 2007.  8.
CabezaLainez Jose M. Lighting Features in Japanese Traditional Architecture. In Lessons from Vernacular Architecture. Earthscan Routledge. 215 pp. August 2013.  9.
Cabeza Lainez Jose M. The quest for light in Indian Architectural Heritage. Journal of Asian Architecture and Building Engineering. Pp. 3946. May 2008.  10.
Cabeza Lainez Jose M, Jimenez Verdejo Juan R. The Keyrole of Eladio Dieste, Spain and the Americas in the Evolution from Brickwork to Architectural Form. Journal of Asian Architecture and Building Engineering (JAABE). Pp. 355362. November 2009.
Notes
 Note that values under 0.5 can also be found for this relationship in a sort of globular cap with an area bigger than the hemisphere.