Open access peer-reviewed chapter

Radiative Heat Transfer for Curvilinear Surfaces

Written By

Jose Maria Cabeza Lainez, Jesus Alberto Pulido Arcas, Manuel- Viggo Castilla, Carlos Rubio Bellido and Juan Manuel Bonilla Martínez

Submitted: June 4th, 2014 Reviewed: October 30th, 2014 Published: May 6th, 2015

DOI: 10.5772/59797

Chapter metrics overview

1,645 Chapter Downloads

View Full Metrics

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.

Advertisement

2. Outline of the problem

The reciprocity principle enunciated by Lambert in 1760 and expressed in Eqn. (1), yields the following well-known integral equation (2) that acts as the theoretical basis for form factor calculation between two surfaces.

d1-2=Eb1-Eb2cosθ1*cosθ2*dA1*dA2π*r2E1
1-2=Eb1-Eb2A1 A2 cosθ1*cosθ2*dA1*dA2π*r2E2

Where the terms are depicted in Figure 1,

Ebi= radiant power emitted by the corresponding surface 1 or 2

Ai= area of surface, dAi= differential of area

r = distance radiovector

θi =angle between radiovector at differential element iand the normal to the surface

Figure 1.

The reciprocity principle and equation for arbitrary surfaces A1 and A2

The previous expression states that radiant interchange for every given form depends on its shape and its relative position in the three-dimensional 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.

Advertisement

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:

F11=1E3

Bearing this in mind, in a similar surface, for instance a hemisphere, the form factor is accordingly F11= ½. The configuration factor of a differential area to a disk of radius runder the center of the disk at precisely the distance r, provides a hint in that it is also ½ [2]. For a point of the hemisphere the factor required is ½.

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.

Figure 2.

A spherical cap of height h and radius of the base a

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  A1·F12=A2·F21; since F21=1, and there is no F22 for planar surfaces, F12=A2A1, in this particular case:

F12=a2a2+h2 E4
F11=h2a2+h2=h2*RE5

Two important laws are inferred from here, which have been defined as Cabeza-Lainez laws:

Cabeza-Lainez 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.

Cabeza-Lainez 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 Yth fraction of the total area of the sphere of radius R, and recalling from trigonometry that,

(h2+a2)=2·R·h E6

Thus,

Y·(h2+a2)=4·R2 5; Y=2·Rh 6; h=2·RY E7

Consequently,

F11=h2·R=h2·Y4·R2=1YE8

Cabeza-Lainez second law:

The configuration factor of an Yth 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 F11 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, cosθand the radius R (Figure 3).

Figure 3.

Differential surfaces in the sphere of centre C and luminance L used to find the radiative exchange

Substituting, these terms in the canonical equation (1):

1-2=Eb14·π·R2A1 A2  dA1·dA2E9

4πR2is the total area of the sphere. Thus, the radiative flux transfer is dependent on the size of the surfaces but not on their position in the sphere and for given areas it is also a constant. Trying to obtain F11=11Eb1.A1from equation (7) gives the expression:

F11=A14·π·R2=1YE10

This means that spherical surfaces present these unique properties (Eqs. 3 and 8) which are crucial for our discussion crucial for our discussion.

Now Cabeza-Lainez 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 xfrom 0 to 180 degrees:

Figure 4.

Denomination of surfaces in a sector of the sphere, 1 and 2 are planar semicircles, 3 is curved.

As has been discussed, the Y portion of the sphere is, in this case 1Y=x360and thus,

F33=x360E11

Accordingly,

F31=F32=12·1-x360E12

And introducing the areas of the semicircles, πR22

F13=F23=x90·1-x360E13

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).

Figure 5.

Two semicircles of the same radius R with a common edge forming an angle X

F12=1-x90+x232400E14

Figure 6.

Radiative exchanges between two semicircles with a common edge and forming an internal anglex

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 h2R (Eq. 5), applied to the spherical cap. Form factors between the contained surfaces are as follows:

F11=h2·R=h2h2+a2E15
F12=a2h2+a2E16
F21=1E17

If we introduce at this point the dimensionless parameter β, we can simplify equation 16 as,

β2=h2a2 E18
F12=1β2+1E19

Since this principle is more general than the second one, we can extend it to non-spherical surfaces.

Advertisement

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.

Figure 7.

Prolate spheroid

Firstly the dimensionless parameter m is introduced:

m=1-a2h2E20

By virtue of the first principle,

F12=a*ma*m+h*arcsin(m)E21
F21=1E22
F11=h*arcsin(m)a*m+h*arcsin(m)E23

And making,

β2=h2a2; m=1-1β2E24
F12=1-1β21-1β2+β*arcsin(1-1β2)E25

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

Figure 8.

Oblate spheroid

Denote the parameter m1,

m1=a2h2-1E26
F12=a*m1a*m1+h*arcsinh(m1) ;F21=1E27

By the first principle and,

F11=h*arcsinh(m1)a*m1+h*arcsinh(m1)E28

With the same procedure as before to make the expression dimensionless

m1=1β2-1E29
F12=m1m1+β*arcsinh(m1)E30

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

F12=6*a*h2a2+4*h23/2-a3 ;F21=1E31
F11=1-6*a*h2a2+4*h23/2-a3E32
β=ha ;F12=6*β21+4*β23/2-1E33

Figure 9.

Paraboloid of revolution

4.4. Right cone

1 is the surface of the cone and 2 is the circular base

F12=aa2+h2 ;F21=1E34
F11=1-aa2+h2E35

Introducing the parameter  β,

F12=11+β 2E36

It is possible to compare the performance in terms of F12, of all the figures found up to now, where the cone shows better performance followed by the paraboloid.

Figure 10.

Cone

Figure 11.

Comparison of form factors for different shapes

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.

F12=a*b*3y2*ay*by+ay*hy++by*hyy ;F21=1E37
F11=1-a*b*3y2*ay*by+ay*hy++by*hyyE38

Figure 12.

Ellipsoid

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,

Figure 13.

Sphere divided in two caps of diverse heights

As an immediate consequence of Cabeza-Lainez laws, rbeing the radius of the inner circle and hthe respective heights of the caps,

F11=F21=h12h12+r2=r2h22+r2=h1*h2h22+r2=h12+r2(h1+h1)2=h1(h1+h2)=h12*RE39
F22=F12=h22h22+r2=r2h12+r2=h1*h2h12+r2=h2(h1+h2)E40

If now the caps within the same sphere are of any size and arbitrary position,

Figure 14.

Two caps of arbitrary size

In this case by virtue of Cabeza-Lainez Law,

F11=h12h12+a2; F22=h22h22+a22E41

And now we need to apply the canonical equation 9 again, substituting the respective areas of the caps; A1=2.π.R.h1.  ; A2=2.π.R.h2

1-2=Eb14*π*R2A1 A2  dA1*dA2E42
F12=h1*h2h12+a2; F21=h1*h2h22+a22E43

In the special situation that the caps are parallel, which equates a truncated cone, the flux would be Eb1.π.h1.h2and the fraction of energy from disk 1 to disk 2 (or their surrounding caps), equates h1*h2/a2or  h1*h2/a22. In the case that the bases are of equal radius a, h1=h2=h. If the perpendicular distance between the disks, called 2b, is known (Figure 15), the height of the cap would be,

h=a2+b2-bE44

Thus, the fraction is obtained as,

F12=F21=a2+2*b2-2*b*a2+b2a2E45

Figure 15.

Surfaces defined by a cylindrical volume used to find the radiative transfer

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.

Figure 16.

Volume composed of a cylinder and a spherical cap used to find the radiative transfer among those surfaces

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.

z=r12-r224*b; R=z+b2+r22E46
l=(r1-r2)2+4*b2E47
Q=R2-z2+b2-2*R*bE48
Q1=r12-Q  ;Q2=r22-QE49
D1=h12+r12E50
D2=h22+r22E51
D3=l*(r1+r2 ) E52

Figure 17.

Volume composed by spherical cap, truncated cone and hemispheroid.

And then we would obtain the corresponding nine form factors involved,

F11=h12D1; F12=QD1; F13=Q1D1E53
F22=h22D2; F21=QD2; F23=Q2D2 E54
F31=Q1D3; F32=Q2D3; F33=1-Q1+Q2D3E55

In this simple way the problem is completely solved

4.7. Straight cone

This is a limit case of the previous discussion.

Figure 18.

Right cone with a circular base

As the former also includes the cone, by making r0=0and h1=h2=0,  Q2=0, z=r124b, R=z+b, Q=0, Q1=0,  Q2=0

l=r12+4*b2E56

If D1=r12, D2=0 then

D3=r12+4*b2 E57

Only three factors remain,

F11=1E58
F31=r1r12+4*b2E59
F33=1-r1r12+4*b2E60

F31 is obviously the ratio of areas of the cone to its base which proves that the equation is true, by virtue of Cabeza-Lainez 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,

Figure 19.

Volume composed by a paraboloid, a truncated cone and a spheroid.

F22=1-6*r2*h22r22+4*h223/2-r23E61

as in the paraboloid alone

F21=6*h22*Qr2*r22+4*h223/2-r23E62
F23=6*h22*(r22-Q)r2*r22+4*h223/2-a23E63
F11=h1*arcsinh(m1)r*m1+h1*arcsinh(m1)E64

as it were in the oblate elipsoid alone

m1 is now=r12h12-1E65
F12=m1*Qr1*(r*m1+h1*arcsinhm1)E66
F13=m1*(r12-Q)r1*(r1*m1+h1*arcsinhm1)E67

F31, F32 and F33 have the same values as before as these correspond to the truncated cone and bear only nominal relation with the surfaces of the extremes,

F31=Q1D3 E68
F32=Q2D3E69
F33=1-Q1+Q2D3E70

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.

F21 is always the unit as shown by first law

SURF.Area of the revolution surfaceF11F12
Prolate semi-spheroid with circular base2πa2arcsenm·h+ama·m
m=1-a2h2
h*arcsin(m)a*m+h*arcsin(m)a*ma*m+h*arcsin(m)
Oblate semi-spheroid with elliptic base2πa2arcsenhm·h+ama·m
m=a2h2-1
h*arcsinh(m1)a*m1+h*arcsinh(m1)a*m1a*m1+h*arcsinh(m1)
Revolution paraboloid with circular baseπa2+4h23-a36a3h21-6*a*h2a2+4*h23/2-a36*a*h2a2+4*h23/2-a3
Straight cone with circular baseπaa2+h21-aa2+h2aa2+h2
Revolution Ellipsoid4πayby+ayhy+byh31y
y=85
1-(ab3y)/2ayby+ayhy++byhyy(ab3y)/2ayby+ayhy++byhyy

Table 1.

Resume of form factors for curved surfaces.

Advertisement

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 re-irradiated and the concepts of emitters and receivers entwine.

Under such circumstance, the global balance of radiant power can be found through expression 71,

Etot=Edir+ErefE71

Edir is defined as the direct power received while Eref stands for the reflected energy. The two quantities added yield the global balance of radiant energy Etot. If the problem entails several surfaces, expression 71 is expanded for an array of equations. To resolve it, we define beforehand the matrices Fd and Fr, whose elements are described as follows in a three-dimensional fashion, (see figure 16):

Fd=F11*ρ1F12*ρ2F13*ρ3F21*ρ1 F22*ρ2F23*ρ3F31*ρ1F32*ρ2F33*ρ3E72
Fr=1-F12*ρ2-F13*ρ3-F21*ρ1 1-F23*ρ3-F31*ρ1-F32*ρ21E73

Each term in equations 72 and 73 is presented in the form Fij (F11, F12...). This stands for the configuration factors already found, from one of the surfaces ito another adjacent surface j. The term ρi is defined as the reflective quotient which corresponds to a given surface i.

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 Fdwere 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:

Fr*Eref=Fd* EdirE74
Frd =Fr-1*Fd  E75
Eref =Frd *Edir  E76

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:

Eref =E*WA*(ρ+ρ2+...ρn)E77

As,

limnρn+1-1ρ-1-1 =ρ1-ρ  E78
Eref =E*WA*ρ1-ρE79

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.

Advertisement

6. Conclusions

An ever-increasing 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 energy-conscious 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:

My name is Ozymandias, King of Kings, Look on my works ye Mighty And Despair

Advertisement

Acknowledgments

Jose Cabeza would like to thank his family in Japan and Spain for failing to understand his work.

References

  1. 1. John R. Howell, A Catalogue of Radiation Heat Transfer Configuration Factors. 3rd ed., 2010. On-line version available at: http://www.engr.uky.edu/rtl/ Catalog/.
  2. 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 D-944.
  3. 3. Cabeza-Lainez Jose M. Solar Radiation In buildings. Performance and Simulation procedures. InTech. 2012.
  4. 4. Cabeza Lainez Jose M. New Configuration Factors for Curved Surfaces. Journal of Quantitative Spectroscopy and Radiative Transfer (JQSRT). Vol. 117. March 2013.
  5. 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. 6. Cabeza Lainez Jose M. Fundamentals of Luminous Radiative Transfer. Netbiblo. 256 pg. December 2010.
  7. 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. 33-40. May 2007.
  8. 8. Cabeza-Lainez Jose M. Lighting Features in Japanese Traditional Architecture. In Lessons from Vernacular Architecture. Earthscan Routledge. 215 pp. August 2013.
  9. 9. Cabeza Lainez Jose M. The quest for light in Indian Architectural Heritage. Journal of Asian Architecture and Building Engineering. Pp. 39-46. May 2008.
  10. 10. Cabeza Lainez Jose M, Jimenez Verdejo Juan R. The Key-role of Eladio Dieste, Spain and the Americas in the Evolution from Brickwork to Architectural Form. Journal of Asian Architecture and Building Engineering (JAABE). Pp. 355-362. 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.

Written By

Jose Maria Cabeza Lainez, Jesus Alberto Pulido Arcas, Manuel- Viggo Castilla, Carlos Rubio Bellido and Juan Manuel Bonilla Martínez

Submitted: June 4th, 2014 Reviewed: October 30th, 2014 Published: May 6th, 2015