Next Article in Journal
Evaluation of the Approach for the Identification of Trajectory Anomalies on CCTV Video from Road Intersections
Next Article in Special Issue
Anomalous Areas Detection in Rocks Using Time-Difference Adjoint Tomography
Previous Article in Journal
Multidimensional Fibonacci Coding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prevention of Hazards Induced by a Radiation Fireball through Computational Geometry and Parametric Design

by
Joseph M. Cabeza-Lainez
1,*,
Francisco Salguero-Andújar
2 and
Inmaculada Rodríguez-Cunill
3
1
Higher Technical School of Architecture, University of Seville, 41012 Seville, Spain
2
Science and Technology Research Center of Huelva, Faculty of Experimental Sciences, Campus de El Carmen, University of Huelva, 21007 Huelva, Spain
3
Department of Painting, Faculty of Fine Arts, University of Sevilla, 41003 Seville, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(3), 387; https://doi.org/10.3390/math10030387
Submission received: 17 December 2021 / Revised: 20 January 2022 / Accepted: 26 January 2022 / Published: 27 January 2022

Abstract

:
Radiation fireballs are singular phenomena which involve severe thermal radiation and, consequently, they need to be duly assessed and prevented. Although the radiative heat transfer produced by a sphere is relatively well known, the shadowing measures implemented to control the fireball’s devastating effects have frequently posed a difficult analytical instance, mainly due to its specific configuration. The objective of this article is to develop a parametric algorithm that provides the exact radiative configuration factors for the most general case in which the fireball is located at any distance and height above the ground, partially hidden by a protective wall over an affected area at different positions with respect to the said fireball. To this aim we use methods based on Computational Geometry and Algorithm-Aided Design; tools that, departing from the projected solid-angle principle, provide exact configuration factors, in all cases, even if they do not present a definite analytical solution. This implies dealing with spatially curved radiative sources which had not been addressed formerly in the literature due to their mathematical difficulties. Adequate application of this method may improve the safety of a significant number of facilities and reduce the number casualties among persons exposed to such risks. As a similar radiative problem appears in volcanic explosions; we hope that further extensions of the method can be adapted to the issue with advantage.

1. Introduction

A Boiling Liquid Expanding Vapor Explosion (BLEVE) is a type of explosion that occurs in tanks that store liquefied and overheated gases under pressure, in which due to rupture or leakage of the tank, the liquid inside boils and is massively incorporated into the expanding vapor. The most frequent cause of this type of explosion is an external combustion that engulfs the pressurized tank, weakens it mechanically, raises the temperature of the liquid contained and increases the pressure inside the tank. There comes a point at which the pressure reaches values that the container cannot withstand, producing a fissure or rupture of the same. This causes a sudden drop in pressure, the spontaneous nucleation process begins, and all the liquid contained changes its state to gas almost instantaneously, increasing its volume hundreds or thousands of times. BLEVE occurs even though the liquid contained is not a flammable product. The combustion of the content will take place as long as the product is combustible and flammable, but this is a second explosion that constitutes a different phenomenon named as Unconfined Vapor Cloud Explosion (UVCE), also known as fireball, and is a consequence of a BLEVE and not part of it.
The formation of a fireball, isolated or as part of a BLEVE, is characterized by the emission of intense thermal radiation, capable of causing lethal and irreversible damage even to people located at significant distances. Although knowledge about this type of accident has improved substantially, it continues to occur, so it is necessary to develop countermeasures aimed at reducing the effects of its obnoxious consequences.
The process of formation and dynamics of fireballs, including their size, duration, elevation, temperature variation and their concentration range and velocities have been studied by numerous authors [1,2,3,4,5,6,7]. Their initial phase is characterized by the typical silhouette of a mushroom, the most widely accepted model for calculating the ensuing radiative configuration factors and their effects consist of a sphere as a representative figure of the so-called fireball, sometimes at ground level and sometimes elevated.
With several caveats, a similar problem appears in some volcanic explosions that require protection for visitors and researchers.
In recent years, ways for modeling heat transfer processes have been registered in the literature [8,9,10]. However, procedures for finding radiative heat transfer through the space are scarce and inaccurate, mainly due to the difficulties posed by configuration factors.
It is true that there are some exact analytical calculations of configuration factors for the particular geometry of certain fragments of spheres [11,12,13,14,15,16], collected mainly in the catalog of Howell and Pinar Mengüç [17]. However, the problem of the effects and, therefore, of the calculation of the configuration factors, from a fireball partially hidden by a protection wall with an arbitrary elevation above the ground has been (to our best knowledge) approached only by J. M. Bonilla [18]. Said sphere affects an area located at any height with respect to its base-level although as the author himself recognizes (p. 97 of his doctoral thesis).
In the case under study, a fireball with an obstacle constituted by a flat wall, the resulting integral equation for the configuration factor F12 (Equation (1)), a concept developed below in Section 2 (Equations (2)–(6)), that measures the non-dimensional ratio of energy transferred from surface 1 to surface 2 (see Figure 1) remains analytically unsolvable [11,16]. Therefore, it is necessary to resort to numerical methods to obtain an adequate solution.
F 12 = 0 π 0 π 0 R 0 R ( r 1 2 r 2 sin θ 1 r 2 2 r 1 sin θ 2 ) d r 1 d r 2 d θ 1 d θ 2 π [ ( r 1 2 + r 2 2 + 2 r 1 r 2 cos θ 1 cos θ 2 ) ] 2
Equation (1), describes the complexity of the factor F12 (concept defined below in Equations (2)–(6), to evaluate the radiative exchanges that occur in Figure 1, as such exchanges typically appear in the problems of the fireball that we are bound to study.
Precisely, the object of this article is to develop a parametric algorithm that provides the exact radiative configuration factors for the most general case in which the fireball can be located at any position over the ground, partially hidden or not, by a protective structure in an affected area surrounding said fireball. The protection measures, possessing a length of several meters and a considerable height, become costly and for that reason they are usually built as vertical walls, However, our procedure could be applied to other shapes with which we are currently experimenting.
In order to solve the complex issues of obstructions in the way to radiation, we would use methods based on Computational Geometry and Algorithm-Aided Design; tools that, originating from the principle of projection of the solid angle, provide exact configuration factors, in all cases, even when they do not present a known analytical solution. Under such procedure, the cases with a definite solution become the results of our work.
In 1928, W. Nußelt, in an article [19] with 54 lines, 5 equations and only one figure drawn by hand, suggested that the configuration factors could be calculated graphically by (1): placing a sphere of unit radius centered on the surface element dA1; (2): obtaining the intersection between the cone with its vertex in the center of said sphere and any planar trapezoidal surface, in any orientation and position that radiates energy (A2) as directrix, and the unit sphere; (3): projecting orthogonally the previously calculated intersection onto the plane containing dA1 and, finally; (4): dividing by π the area enclosed by this last orthogonal projection. In this way we would obtain the configuration factor fdA1A2 (Figure 2). Despite the short length of Nußelt’s article, his intuition was taken at face value without proper demonstration, and consequently introduced the possibility of a geometric approach to a difficult analytical problem. It is convenient to stress that Nußelt did not mention non-planar or curved emitting sources. Unlike Nußelt, we needed to investigate precisely this possibility and, starting from his conjecture, we have developed our article by means of algorithm aided design to provide a solution for such complex problems.

2. Materials and Methods

2.1. Fundamental Equations: Form Factors and Configuration Factors

Since the fundamental equations governing radiant energy transfer are well established and systematically discussed in most articles dealing with radiant transfer, we will present here in short only those relationships strictly necessary for a complete understanding of our work. For more in-depth references on the subject, which do not address non-planar radiative emitters, see for example [20,21,22,23].
Suppose two surfaces Ai and Aj that emit diffuse radiation as shown in Figure 3. We wish to obtain the energy exchange between both surfaces. If, initially, we ignore the details of the transmission, the problem essentially becomes one of determining the amount of energy that leaves one of the surfaces and reaches the other.
If we assume for convenience that there are no inter-reflections from one surface to the other, all the incident radiation will be absorbed and the differential flux of energy gives:
d ϕ i j = ( E i E j ) cos θ i cos θ j d A i d A j π r i j 2
Being Ei and Ej the amounts of energy (in W/m2) emitted by surfaces i and j.
Next, let us define the view factors (or form factors), Fij as the fraction of energy that radiates from surface i intercepted by surface j. Then the energy that leaves the surface i and reaches the surface j will be EiAiFij and, reciprocally EjAjFji represents the energy that passes from surface j to surface i.
Then Equation (2) could be rewritten in the form:
d ϕ = E i A i F i j E j A j F j i .
Let us now consider a stationary phase of the UCVE with constant dimensions of the radiant sphere, under which there is no increase or decrease in flux and that all possible losses due to the mode of transmission are negligible. Thus = 0, which implies that EiAiFij = EjAjFji, an expression known as Lambert’s reciprocity theorem [24].
If we also assume (since it is a frequent case of heat transfer for UVCEs, the negligible emission of affected bodies) that surface Aj does not emit energy and only receives heat from Ai, the total flux will be:
ϕ = E i A i A j cos θ i cos θ j d A i d A j π r i j 2 ,
Thereupon,
E i A i F i j = E i A i A j cos θ i cos θ j d A i d A j π r i j 2 ,
and finally we get:
F i j = 1 A i A i A j cos θ i cos θ j d A i d A j π r i j 2 .
Which reduces the problem to the calculation of purely geometric expressions. A particular case of this equation for the problem of the fireball was presented in Equation1 where we expressed our concerns about the difficulties of solution that it posed.
However, the form factor that Equation (6) gives is the net radiant heat exchange from surface to surface. A surface that emits a certain amount of energy Ei successively induces a certain “average” energy E2 on an adjoining surface, which tells us very little about the detailed distribution of said energy.
When analyzing the effects of a UVCE, it is necessary to describe the flux function point by point, i.e., the radiant field, and thus be able to evaluate the possible damage to people and buildings. Therefore, if we look for the point to pint distribution of energy that reaches the surface j, instead of the final value averaged over the receiving surface, we can effectively reduce the integral of Equation (6) as Ai is condensed in a point. In other words, we only need to perform the inner part of Equation (6). Such integral will consequently be double, not quadruple, resulting in a new expression, also strictly geometric, known as pointwise view factor or configuration factor:
f d A i A j = A j cos θ i cos θ j π r i j 2 d A j

2.2. Nußelt’s Analogy or Solid Angle Proyection Law

As explained in the introduction, Nußelt described in 1928 a hypothesis to calculate form factors from planar figures, that acquired wide predicament but had not been properly demonstrated since its publication. This proposal is also known as the Solid Angle Projection Law.
His ideas follow the geometric logic explained in the segueing lines. Let us consider the planar surface element dAj and connect all the points of its contour with the center of a unit sphere. In this way we will obtain a conical surface, which in turn delimits an area dω over the sphere (Figure 4). The solid angle under which a surface is seen from a point is defined as the area of the conical projection of said surface over a sphere of unit radius centered at that same point.
The vector er is the unit vector in the direction rij and the scalar product (dot product) er·dAj = cos θjdAj represents the projection of the vector dAj in the radial direction er, that is, cos θidAj gives the projection of the element of area dAj over a perpendicular plane to the direction er. Now, a simple similarity relation between dAj at a distance rij and the subtended area in the unit sphere, allows us to calculate the solid angle element as:
d ω = cos θ j d A j r i j 2 .
Then the solid angle under a finite surface ω viewed from the center of the unit sphere will be,
ω = A j cos θ j r i j 2 d A j ,  
and its projection onto a plane that forms an angle θi with er will have a surface (with respect to the area π, of the unit circle):
ω cos θ i π ,  
Thus, substituting ω for its value we obtain:
cos θ i π A j cos θ j r i j 2 d A j ,  
Since θi is independent from Aj, the orthogonal projection of ω onto said plane (divided by the area of the unit circle π) yields the dimensionless configuration factor previously obtained (Equation (7)).
Furthermore, taking into account Lambert’s reciprocity theorem, any diffusive surface (Ak) adscribed to the cone under which Aj is seen, will have the same configuration factor since, as we have assumed, = 0 and therefore fdAiAj = fdAiAk (Figure 5). An important consequence of the above reasoning is that configuration factor problems, can be reduced in most cases to the calculation of a projected area.
We would like to outline that the values obtained for the configuration factors through the application of the projected solid-angle principle are not a mere graphic approximation but should be identical to the results of analytical methods
That said, demonstrating the former principle by integral calculus becomes, in most cases, extremely difficult [18]. Cabeza-Lainez [20,25] has exposed this hindrance of Nußelt’s work for the apparently simple case of the configuration factor of a rectangular emitter as defined in Figure 6b. The particulars of such complex integral follow.
The exact analytical formula obtained by direct integration of the configuration factor [23] is (Figure 6a):
f 12 = 1 π [ D H 1 2 + D 2 tan 1 ( M H 1 2 + D 2 ) D H 0 2 + D 2 tan 1 ( M H 0 2 + D 2 ) ] ,
Equation (11) substituting the data consigned in Figure 6b gives a result of:
f d A 1 A 2 = 1 π [ 2 5 tan 1 ( 1.5 5 ) 2 8 tan 1 ( 1.5 8 ) ] = 0.0584 .
On the other hand, and surprisingly enough, if we want to verify these values using strictly the projected solid angle principle, the only viable approach to solve the problem is to use the “generalized polar coordinates” in which: x = cosθ; y = sinθ, other more usual coordinates produced no positive results for reasons that might be related with the special nature of elliptic integrals. Those reasons will be analyzed in the discussion section.
For the said coordinates, the equation of the ellipse takes the form ρ = 1 and the unit surface element is: dA = abρdρdθ.
Thus, to calculate the area of the elliptical sector we will have to solve:
A L A M = θ 2 θ 1 0 1 a b ρ d ρ d θ = a b θ 2 θ 1 [ ρ 2 2 ] 0 1 d θ = a b ( θ 1 θ 2 ) 2 = a b θ
since, by symmetry, θ1 = θ2.
For the larger ellipse, a = 1 m, b = 0.8944 m and θ = tan−1(2/(1.5·0.8944)) = 0.9799 and its complement with respect to π/2 is (π/2) − θ = θL = 0.5908, then the larger ellipse elliptical sector area is: Al = abθL = 1 × 0.8944 × 0.5908 = 0.5284 m2.
Note the singular fact that the θ angles in these coordinates are different for each ellipse as they are colligated to the respective sizes of their major and minor axes. As a consequence, the elliptical sector of the minor ellipse area will be, a = 1 m, b = 0.7071 m and θ = tan−1(2/(1.5 ×·0.7071)) = 1.0832. Its complementary with respect to π/2 is (π/2) − θ = θm = 0.4876 therefore, Am = abθm = 1·0.7071 × 0.4876 = 0.3448 m2.
Finally, subtracting the areas of the previously calculated elliptical sectors and dividing by π, we obtain exactly the same result as in Equation (13), i.e., ((0.5284 − 0.3488)/π) = 0.0584.
Unfortunately, we spent a considerable amount research time to attain this finding which is only valid for the rectangle. Still, other simple shapes of emitters such as the circle (Figure 7) or the ellipse remain unsolved, by direct planar integration or any other method [25]. Moreover, calculations for three-dimensionally curved emitting sources, such as for instance the fragment of sphere that precisely appears in the fireball, become even more complicated.
In the case of a circular emitting source displaced from the center of the unit sphere, the intersection with the said sphere has been defined by Cabeza-Lainez as a Tomomi curve [25]. This curve (Figure 8), presents the following parametric coordinates,
x = sin α cos t  
y = sin β sin t  
z = cos 2 α cos 2 t + cos 2 β sin 2 t  
Its projection on a plane is actually a fourth-degree curve as shown in Figure 9 and can be deducted from Equation (15); therefore, it is not possible to find its enclosed area by direct integration in an exact manner.
These findings convinced us of the necessity of establishing other approaches for the issue that will be explained in detail in the ensuing section. The above results for the rectangle and the circle have been thoroughly proven by the authors with the proposed new method and they can be found in reference [25].

2.3. Computational Geometry and Algorithm-Aided Design

Amongst the different classifications of widespread methods for the determination of configuration factors, a clear line can be traced between analytical and numerical methods, with an almost general omission of experimental methods, arising in the literature.
Although it is true that at the time that Nußelt proposed the hypothesis that bears his name, the available techniques lacked sufficient accuracy to develop graphic methods (so that these methods fell into disuse). Currently, with the development of Computational Geometry and Algorithm-Aided Design, under programs such as Grasshopper® from Rhinoceros a recovery of the origins of graphical methods and their preferential use in a myriad of problems has been enforced [26]. We firmly believe that the solid angle principle after Cabeza Lainez’ generalization effected in Section 2.2, can be applied to advantage to all kinds of non-planar sources and, consequently, to the problems posed by the fireball, assuming diffuse radiation.

2.4. Calculation Algorithm

The ensuing steps necessary to obtain the code that gives the configuration factors of arbitrary points on differential surfaces, arbitrarily oriented and positioned, for a fireball of variable height are described below.

2.4.1. Fireball Position and Radius

First, we would write a simple algorithm that sets and constructs a spherical fireball of radius R (which basically depends, as we will see, of the amount and type of combustible substance) and (for convenience’s sake), center at point (0, 0, Hf) (Figure 10). Note that the R and Hf parameters are provided by two sliders of Grasshopper, so the radius of the fireball and its height above the ground are modifiable at will.
The result of conducting the previous algorithm is shown in Figure 11.

2.4.2. Position, Height, Extremes of the Protective Wall and Height of the Affected Area

The second step (Figure 12), needed for the complete definition of the problem, is to set the protective wall at a certain distance from the fireball (yw) and define its height (Hw) and its extremes. Finally, to give more generality to these types of problems, we have conceived the possibility that the fuel tanks or the vehicles that transport the combustible could be located on a hill that overlooks the affected area by means of the Hg parameter. Again, all these values are modifiable by the user depending on the specific problem to be addressed. An overview of the parameters involved in this type of situation is shown in Figure 13.

2.4.3. Calculation of the Shadow Zone Provided by the Protective Wall

To calculate the shadow area provided by the protective wall (Figure 14), we first calculate the intersection of the fireball with a vertical plane that passes through the center of the sphere and contains the midpoint of the highest line of the wall by means of the brep|plane command, obtaining a circumference perpendicular to the wall enclosed in the sphere.
Next, we draw the tangent lines to the said circumference from the midpoint of the upper line of the wall and select the highest of the two tangents found (tangent lines function). The intersection of the plane that contains the highest line of the wall and the tangent previously chosen with the ground plane (line|plane command), gives us the y coordinate of the shadow line on the exposed area, which separates the real affected zone from the area of shadow. A graphic representation of the result of the previous process is shown in Figure 15.

2.4.4. Creation of a Grid of Control Points to Calculate Configuration Factors

When systematizing and automating the process, it is convenient to create a mesh of control points on which to calculate the different configuration factors between the spherical cap, resulting from the section of the fireball produced by the shadow plane and that given by the plane perpendicular to the line connecting each point with the center of the fireball. As can be seen in the code depicted in Figure 16, the affected surface to be evaluated has the same width as the protection wall and its length can be freely modified by means of a slider. We also want to note that the first row of control points have been omitted, since in that entire row the configuration factors are null by definition (Figure 17).

2.4.5. Determination of the Radiant Heat Cones with Shadow

We will now proceed to calculate the radiant heat cones with vertex at the control point tangent to the fireball, subtracting from them the part that lies below the protective wall. To perform this, given that the radiant heat cones and the fireball are two quadrics with two planes of common symmetry—one of them vertical that passes through the center of the fireball and through the control point—we would cut the fireball by the said plane to obtain a circumference contained in the fireball and in the aforementioned vertical plane. Tracing the tangents to above-defined circumference from the control points and forming the surface of revolution that generates by rotating a 2π angle, one of the previous tangents (generatrix), around the line that joins the control point and the center of the fireball (axis), (by means of the revolution function) we obtain the complete radiant heat cone.
Next, we must cut the newly calculated cone and the fireball along the plane defined by the control point and the upper line of the protection wall (Brep|Plane command) by means of the Surface Split function, keeping the upper part of the three surfaces found by this function (Figure 18). The result is shown in Figure 19a.

2.4.6. Application of the Projected Solid Angle Principle and Determination of the Configuration Factor for an Arbitrary Point on the Surface Partially Affected by a Fireball

In order to achieve this, we would place a sphere of unit radius at the chosen point and obtain the intersection between the sector of the cone calculated in the previous section with the said sphere, obtaining an arc of circumference over the unit sphere (Figure 19b). To close the solid angle, we need to compute the intersection circumference between the unit sphere and the plane that yielded the radiating cone sector above the wall.
Of this circumference we are only interested in the smallest portion that is included between the initial and final points (End Points) of the intersection formed by the radiating cone sector and the unit sphere. To find the value of such piece, we take the midpoint of the segment that links the initial and final points and determine the location of the closest point of the unit sphere (Curve Closest Point). Thus, we will arrive to the arc that passes through the ends of the intersection circumference segment, between the unit sphere and the radiating cone sector and its nearest point.
Finally, it remains to connect the arcs previously calculated (Join Curves) and draw their orthogonal projection on any desired plane (vertical, horizontal or maximum exposure situation). In this case we have selected the orthogonal projection on the affected ground plane.
To complete the problem, Grasshopper calculates the area of said projection on the ground plane and divides it by π, yielding the configuration factor of the arbitrary point chosen (Figure 20).
By automating the process through a loop for all the points of the control grid, we obtain a matrix of n × n elements that constitute each configuration factor of the affected zone considered, in this case on a meter by meter basis. The parametric nature of the code developed, allows us to obtain them for any distance and plane (horizontal, vertical or maximum).
In the following example, the parameters that have been considered are:
  • Fireball radius (R): 5 m.
  • Fireball height (Hf): 0 m.
  • Wall height: (Hw): 2 m.
  • Wall distance (Yw): 10 m.
  • Ground level (Hg): 0 m.
  • Horizontal affected area: 40 × 40 m: 1681 points.
The output is presented in the following results section.

2.4.7. Application of the Proposed Method to Find the Excess of Thermal Radiation Due to a Fireball

Returning to the initial problem of radiative heat transfer, we must reconsider the physical parameters on which the heat exchange depends in order to control the likelihood of a fireball inducing damages on the personnel present at the facilities.
To this aim we will use an empirical formula that according to [27,28] gives Thermal Radiation Dose (DTU) that produces a small probability of lethality for an average population, with second degree burns (1% lethality), estimated at 1000 DTU:
D T U = D t w = t ( I w ) 4 3 ,
where Iw is the intensity of the radiation in kW/m2 and,
t = 0.41 M 0.34
is the time required for the fireball to reach its maximum size and M stands for the mass of fuel.
Iw must be defined as a function of the configuration factors and the emissive power:
I w = τ f m a x E p
wherein fmax is the maximum configuration factor at the point considered (given directly by the algorithm proposed in the methods section), and we need to complete the expression with,
τ = 2.02 ( P w d f ) 0.09
Which gives the atmospheric transmissivity, Pw is the partial pressure of vapor (typically 1.15 Pa for a 50% of relative humidity of the air) and df is the distance from the center of the fireball to the control point considered (also automatically provided by our algorithm).
On the other hand, Ep is the emissive power and is better estimated by the equation:
E p = η r a d M H c π D 2 t
where, as expressed in references [27,28], ηrad is the fraction of the fuel mass involved in the explosion (typically equal to 0.25). M, as we have mentioned, is the fuel mass and Hc is the calorific power of said fuel and D is provided by the expression:
D = 6.14 M 0.325
A calculation example with the relevant quantities is detailed in the results chapter.

3. Results

Following what has been explained in the previous section about our methods, we have been able to obtain an innovative representation (Figure 21) of the horizontal configuration factors (i.e., the factors projected onto the horizontal plane) for a shadowed fireball. All relevant parameters of the figure are mentioned in the caption.
In this way we can determine the radiative intensities due to a shadowed fireball of variable dimensions and positions and, at the same time, provide the disaster prevention department with an effective design tool. In Figure 22, (vide infra) we have compared these results with the case of the unobstructed fireball [25].

Applied Mathematics. Safety Checks

Although the most complex issue when calculating the damage that a fireball can produce is to obtain the configuration factors, we must not forget that it is still an issue of radiative heat transfer. Thus, we must take into account the physical magnitudes on which said transfer depends to estimate the hazards that a fireball may cause on personnel present at the site.
Using equations 16 to 21, as previously detailed, we have proceeded with the following case study. We would suppose a tank located in an airport at ground level loaded with 40,000 L of the most widely used kerosene in commercial aviation, Jet A-1, with a density of 0.8 kg/L and a calorific power Hc = 42,800 kJ/kg. We will thus have a mass of fuel M = 32,800 kg that would cause a fireball (assumed at ground level) of an estimated diameter D of = 6.14 × 32,800(0.325) ≅ 179 m.
Imagine also that a protective wall is located at 100 m from the kerosene tank and presents a height of 6 m above the ground. Substituting in Equation (18) t and Iw by the values provided by Equations (18)–(21), through processing and reworking the expressions found on Ref. [15], we would obtain a DTU of:
D T U = D t w 1.000 13.9487 [ 493.7603 f m a x ( 1.1550 d f ) 0.09 ] 4 3 ,
from where we will find a safety distance with a lethality of less than 1% (Dtw = 996.8252 (kW/m2)4/3 s) of 105.75 m measured from the protection wall, with a maximum configuration factor of fmax = 0.1528 (as provided by our method) at df = 224.17021 m.

4. Discussion

The radiant field caused by a sphere which constitutes the fireball in its most usual appearance cannot be determined by analytical methods in the presence of defensive frames that obstruct the radiation. The fragments of sphere that induce the emitting surface are of an irregular nature not suited to conventional integral calculus.
The main reason for this difficulty is, in the first place, uncertainty on the geometrical boundary of the source that can be no longer taken as a sort of disk but as complex circular or spherical fragment. In the second place, the variability of the said emitter depends on the receiving point considered. In other words, the geometry of the source is not fixed in space.
To overcome such major drawbacks, we have resorted to the principle of the projected solid angle, which was suggested by Nuβelt but was not properly demonstrated by means of area integrals. In fact, in reference [21] it is stated literally that “the integration (of solid angle projection) is generally considerably less involved than the one with general equations”. Unfortunately, we believe that such redundant adverbs suggest more wishful thinking than a fact responding to reality for most cases of configuration factors.
Thus, we have extended the validity of the principle, first by providing an analytical proof for the case of a simple rectangular source.
In this common case for instance, polar coordinates, such as x = ρcosθ; y = ρsinθ, are not well suited to the singularity of elliptic fragments. Initially they were employed for the problem since they offer the undoubted advantage that the parameter angle θ is the same for all the figures that we can represent with such coordinates. On the contrary, with generalized polar coordinates, the angle θ, must be defined for each separate figure or ellipse in this case, as it is not constant or depends on external values. This seemed a strange fact and enough reason to discard it in the initial steps, although it proved essential as we would demonstrate below. We believe that such a situation might be related to the singularities of elliptic geometry which led, for instance, to open questions about elliptic integrals or the volume of a non-revolution ellipsoid.
Subsequently, we have applied the initial steps of our proposed algorithm-aided method to a variety of configuration factors for planar shapes as described in reference [25]. In a similar fashion, and to complete the study at the demand of the editors, the unobstructed radiant field due to a fireball has been compared with our method through Figure 22, showing a total coincidence.
The main novelty of this article is that in it we have been able to extend for the first time our method to 3D spherical irregular sources and we have added the shadow systems and obstructing frames to our calculations. We believe that such findings will mitigate the casualties experienced in this kind of accident.
In this manner, as the proposed computational algorithm is parametric, we could evaluate the distance and height of the protective wall or structure to be built, for any given fuel mass that causes a previously established DTU. In addition, the effect of several fireballs acting simultaneously can be evaluated with a similar procedure. In the future we hope to be able to evaluate the dynamic performance of the fireball phenomenon and to develop some experimental validation methods for these findings.

5. Conclusions

This article proposes and develops a new graphic algorithm that gives way to new codes useful for calculating the configuration factors produced by a fireball of any diameter and height obstructed by a protective wall in an innovative and accurate manner.
The viability of the method resides firstly in the novel extension of the project solid-angle principle to non-planar emitters, as shown in Section 2.2, and secondly in identifying the hidden geometric features of the surfaces involved.
The radiant shapes in these cases are evolving in a rather unpredictable way and could not be addressed analytically.
Our procedure solves the configuration factors exactly, as was previously shown, even in the presence of dynamically evolving fireballs. The simpler cases of the phenomenon have been evaluated with our method showing a total correspondence that we are convinced that will bring us the momentum for further attempts of extension.
In a manner of conclusion our manuscript presents an innovative calculation for an exacting problem that induces many casualties and causes much economic damage each year at a worldwide scale. Besides, as the thermal radiation dose received by each point of the spatial zone defined is found and deftly determined, critical aspects of the design of protections of various types can be implemented, enhancing the safety and well-being of the personnel and industries concerned in the near future by virtue of computer aided numerical-geometric methods.

Author Contributions

Conceptualization, J.M.C.-L.; methodology, F.S.-A.; software, F.S.-A.; validation, J.M.C.-L., F.S.-A. and I.R.-C.; data provision, I.R.-C.; writing—original draft preparation, F.S.-A. and J.M.C.-L.; visualization, F.S.-A. and J.M.C.-L.; supervision, J.M.C.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

F.S.-A. would like to dedicate this work to his beloved son Enrique Salguero-Jiménez from whom he learned all he knows about Grasshopper®. J. C-L. recognizes Tomomi Odajima for strategic support and Juan Manuel Bonilla for his technical vision and contribution.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gostintsev, Y.; Solodovnik, A.; Lazarev, V. Theory of the Aerodynamics, Self-Ignition and Burnup of Turbulent Thermals, Vortex Rings, and Jets in a Free Atmosphere. Khim. Fiz. 1982, 9, 1279–1290. [Google Scholar]
  2. Prug, R.W. Quantitative Evaluation of Fireball Hazards. Process Saf. Prog. 1994, 13, 83–91. [Google Scholar] [CrossRef]
  3. Surzhikov, S.T. Thermal Radiation of Large-Scale Oxigen-Hidrogen Fireballs. Analisis of the Problem and Main Results. Teplofiz. Vysok. Temp. 1997, 35, 416–423. [Google Scholar]
  4. Makhviladze, G.M.; Roberts, J.P.; Yakush, S.E. Fireball During Combustion of Hydrocarbon Fuel Releases I. Structure and Lift Dynamics. Combust. Explos. Shock Waves 1999, 35, 359–369. [Google Scholar] [CrossRef]
  5. Makhviladze, G.M.; Roberts, J.P.; Yakush, S.E. Fireball During Combustion of Hydrocarbon Fuel Releases II. Thermal Radiation. Combust. Explos. Shock Waves 1999, 35, 219–229. [Google Scholar] [CrossRef]
  6. Roberts, T.; Gosse, A.; Hawksworth, S. Thermal Radiation from Fireballs on Failure of Liquefied Petroleum Gas Storage Vessels. Process Saf. Environ. Prot. 2000, 78, 184–192. [Google Scholar] [CrossRef]
  7. Center for Chemical Process Safety. Chapter 6. Basic Principles of BLEVEs. In Guidelines for Evaluating the Characteristics of Vapor Cloud Explosions, Flash Fires and BLEVEs, 2nd ed.; American Institute of Chemical Engineers: New York, NY, USA, 1994. [Google Scholar] [CrossRef]
  8. Yu, Y.; Luo, X.; Zhang, H.; Zhang, Q. The Solution of Backward Heat Conduction Problem with Piecewise Linear Heat Transfer Coefficient. Mathematics 2019, 7, 388. [Google Scholar] [CrossRef] [Green Version]
  9. Savovic, S.; Drljaca, B.; Djordjevich, A. A comparative study of two different finite difference methods for solving advection–diffusion reaction equation for modeling exponential traveling wave. Ric. Mat. 2021, 70, 23. [Google Scholar] [CrossRef]
  10. Ivanovic, M.; Svicevic, M.; Savovic, S. Numerical solution of Stefan problem with variable space grid method based on mixed finite element/finite difference approach. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 2682–2695. [Google Scholar] [CrossRef]
  11. Chung, B.T.F.; Sumitra, P.S. Radiation Shape Forms from Plane Point Sources. J. Heat Transf. 1972, 94, 328–330. [Google Scholar] [CrossRef]
  12. Chung, B.T.F.; Naraghi, M.H.N. Some Exact Solutions for Radiation View Factors from Spheres. AIAA J. 1981, 19, 1077–1081. [Google Scholar] [CrossRef]
  13. Juul, N.H. Diffuse Radiation View Factors from Differential Plane Sources to Spheres. J. Heat Transf. 1979, 101, 558–560. [Google Scholar] [CrossRef]
  14. Hauptmann, E.G. Angle Factors Between a Small Flat Plate and a Diffusely Radiating Sphere. AIAA J. 1968, 6, 938–939. [Google Scholar] [CrossRef]
  15. Naraghi, M.H. Radiative View Factors from Spherical Segments to Planar Surfaces. J. Thermophys. Heat Transf. 1988, 2, 373–375. [Google Scholar] [CrossRef]
  16. Cabeza-Lainez, J.M.; Pulido-Arcas, J.A.; Castilla, M.-V. New Configuration Factor between a Circle, a Sphere and a Differential Area al Random Positions. J. Quant. Spectrosc. Radiat. Transf. 2013, 129, 272–276. [Google Scholar] [CrossRef]
  17. Howell, J.R.; Mengüç, P. A Catalog of Heat Transfer Configuration Factors; Taylor & Francis/CRC: New York, NY, USA, 2010. [Google Scholar]
  18. Bonilla, J.M. Contribución al Desarrollo de Factores de Configuración de Bolas de Fuego con Presencia de Obstáculos. Ph.D. Thesis, Univesitat Politècnica de Catalunya-BarcelonaTech, Barcelona, Spain, 2017. [Google Scholar]
  19. Nußelt, W. Graphische Bestimmung des Winkelverhältnisses bei der Wärmestrahlung. Z. Ver. Dtsch. Ing. 1928, 72, 673. [Google Scholar]
  20. Moon, P.H. The Scientific Basis of Illumination Engineering; Dover Publications: New York, NY, USA, 2008. [Google Scholar]
  21. Modest, M.F. View Factors. In Radiative Heat Transfer, 3rd ed.; Academic Press: Cambridge, MA, USA, 2013; pp. 129–159. [Google Scholar]
  22. Howell, J.R.; Pinar Mengüç, M.; Daun, K.; Siegel, R. Thermal Radiation Heat Transfer, 7th ed.; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  23. Cabeza-Lainez, J.M. Fundamentos de Transferencia Radiante Luminosa o La Verdadera Naturaleza del Factor de Forma y sus Modelos de Cálculo; Netbiblio: Seville, Spain, 2010. [Google Scholar]
  24. Lambert, J.H. Photometry, Or, on the Measure and Gradiations of Light, Colors and Shade: Translation from the Latin of Photometria, Sive De Mensura et Gradibus Luminis, Colorum el Umbrae; with Introductory Monograph and Notes by David L. DiLaura; Illuminating Engineering Society of North America: New York, NY, USA, 2001. [Google Scholar]
  25. Salguero-Andújar, F.; Cabeza-Lainez, J.-M. New Computational Geometry Methods Applied to Solve Complex Problems of Radiative Transfer. Mathematics 2020, 8, 2176. [Google Scholar] [CrossRef]
  26. Tedeschi, A. AAD_Algorithms-Aided Design. Parametric Strategies Using Grasshopper®, 1st ed.; Le Penseur Publisher: Brienza, Italy, 2014. [Google Scholar]
  27. Eisenberg, N.A.; Lynch, C.J.; Breeding, R.J. Vulnerability Model: A Simulation System for Assessing Damage Resulting from Marine Spills; Final Report AD-A-015245; Enviro Control, Inc.: Rockville, MD, USA, 1975. [Google Scholar]
  28. HSE. Methods of Approximation and Determination of Human Vulnerability for Offshore Major Accident Hazard Assessment; HSE: Bootle, UK, 2010; pp. 1–55. Available online: https://www.hse.gov.uk/foi/internalops/hid_circs/technical_osd/spc_tech_osd_30/spctecosd30.pdf (accessed on 17 December 2021).
Figure 1. Radiative exchange between two semicircles forming a common angle α.
Figure 1. Radiative exchange between two semicircles forming a common angle α.
Mathematics 10 00387 g001
Figure 2. Nußelt’s conjecture about radiative emission of a trapezoid (A,B,C,D) and its respective projection onto a hemisphere of center (X0,Y0,0), radius 1 and Equation (XX0)2 + (YY0)2 + Z2 = 1 for Z > 0. Points A, B, C and D which define surface A2, are coplanar in the plane z = k. Such points must be joined with said center (X0,Y0,0) to obtain the projection.
Figure 2. Nußelt’s conjecture about radiative emission of a trapezoid (A,B,C,D) and its respective projection onto a hemisphere of center (X0,Y0,0), radius 1 and Equation (XX0)2 + (YY0)2 + Z2 = 1 for Z > 0. Points A, B, C and D which define surface A2, are coplanar in the plane z = k. Such points must be joined with said center (X0,Y0,0) to obtain the projection.
Mathematics 10 00387 g002
Figure 3. Arrangement of surface-source elements. where θi and θj represent, respectively, the angles formed by the normals to the differential elements of area dAi and dAj and rij is the vector that joins the differential surfaces dAi and dAj.
Figure 3. Arrangement of surface-source elements. where θi and θj represent, respectively, the angles formed by the normals to the differential elements of area dAi and dAj and rij is the vector that joins the differential surfaces dAi and dAj.
Mathematics 10 00387 g003
Figure 4. Solid angle and its projection.
Figure 4. Solid angle and its projection.
Mathematics 10 00387 g004
Figure 5. Analogy between the configuration factors of surfaces ascribed to the same radiant cone.
Figure 5. Analogy between the configuration factors of surfaces ascribed to the same radiant cone.
Mathematics 10 00387 g005
Figure 6. Demonstration of the solid angle projection principle for a rectangular emitter [23]. (a) definitions. (b) dimensions in the selected example.
Figure 6. Demonstration of the solid angle projection principle for a rectangular emitter [23]. (a) definitions. (b) dimensions in the selected example.
Mathematics 10 00387 g006
Figure 7. Projection on the horizontal plane of the intersection (named Tomomi curve, Equation (15)) of an oblique cone with circular basis and a unit radius sphere.
Figure 7. Projection on the horizontal plane of the intersection (named Tomomi curve, Equation (15)) of an oblique cone with circular basis and a unit radius sphere.
Mathematics 10 00387 g007
Figure 8. Axonometric depiction of Tomomi curve. Equation (15).
Figure 8. Axonometric depiction of Tomomi curve. Equation (15).
Mathematics 10 00387 g008
Figure 9. General projected intersection on a plane of Figure 7 and Figure 8, this curve is deemed to be quartic.
Figure 9. General projected intersection on a plane of Figure 7 and Figure 8, this curve is deemed to be quartic.
Mathematics 10 00387 g009
Figure 10. Initial Algorithm to represent the fireball.
Figure 10. Initial Algorithm to represent the fireball.
Mathematics 10 00387 g010
Figure 11. Fireball of radius R and height above the soil Hf.
Figure 11. Fireball of radius R and height above the soil Hf.
Mathematics 10 00387 g011
Figure 12. Other geometric parameters that define the problem in a general way.
Figure 12. Other geometric parameters that define the problem in a general way.
Mathematics 10 00387 g012
Figure 13. Graphic representation of the dependent parameters of the solution of a fireball with shadows caused by protective walls.
Figure 13. Graphic representation of the dependent parameters of the solution of a fireball with shadows caused by protective walls.
Mathematics 10 00387 g013
Figure 14. Calculation of the shadow zone provided by the protection wall.
Figure 14. Calculation of the shadow zone provided by the protection wall.
Mathematics 10 00387 g014
Figure 15. Shadow zone representation.
Figure 15. Shadow zone representation.
Mathematics 10 00387 g015
Figure 16. Algorithm to calculate an arbitrary mesh of control points.
Figure 16. Algorithm to calculate an arbitrary mesh of control points.
Mathematics 10 00387 g016
Figure 17. Control point mesh.
Figure 17. Control point mesh.
Mathematics 10 00387 g017
Figure 18. Algorithm for calculation of the radiant heat cones with shadow.
Figure 18. Algorithm for calculation of the radiant heat cones with shadow.
Mathematics 10 00387 g018
Figure 19. (a) Calculation of the radiant heat cone induced by a snowball partially behind a protective wall. (b) Detail of the previous figure to appreciate the application of the solid angle principle.
Figure 19. (a) Calculation of the radiant heat cone induced by a snowball partially behind a protective wall. (b) Detail of the previous figure to appreciate the application of the solid angle principle.
Mathematics 10 00387 g019
Figure 20. Application of the solid angle principle and output of the configuration factor for an arbitrary point of the affected area.
Figure 20. Application of the solid angle principle and output of the configuration factor for an arbitrary point of the affected area.
Mathematics 10 00387 g020
Figure 21. Newly found representation (1681 control points) of the configuration factors (non-dimensional) on the horizontal plane produced by a fireball over a 40 × 40 m2 area with the following parameters: R = 5 m; Hf = 0 m; Hw = 2 m; Yw = 10 m and Hg = 0 m.
Figure 21. Newly found representation (1681 control points) of the configuration factors (non-dimensional) on the horizontal plane produced by a fireball over a 40 × 40 m2 area with the following parameters: R = 5 m; Hf = 0 m; Hw = 2 m; Yw = 10 m and Hg = 0 m.
Mathematics 10 00387 g021
Figure 22. The radiative field obtained with the method proposed in the article from an unobstructed fireball of 5 m radius, over an affected zone of 40 × 50 m2.
Figure 22. The radiative field obtained with the method proposed in the article from an unobstructed fireball of 5 m radius, over an affected zone of 40 × 50 m2.
Mathematics 10 00387 g022
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cabeza-Lainez, J.M.; Salguero-Andújar, F.; Rodríguez-Cunill, I. Prevention of Hazards Induced by a Radiation Fireball through Computational Geometry and Parametric Design. Mathematics 2022, 10, 387. https://doi.org/10.3390/math10030387

AMA Style

Cabeza-Lainez JM, Salguero-Andújar F, Rodríguez-Cunill I. Prevention of Hazards Induced by a Radiation Fireball through Computational Geometry and Parametric Design. Mathematics. 2022; 10(3):387. https://doi.org/10.3390/math10030387

Chicago/Turabian Style

Cabeza-Lainez, Joseph M., Francisco Salguero-Andújar, and Inmaculada Rodríguez-Cunill. 2022. "Prevention of Hazards Induced by a Radiation Fireball through Computational Geometry and Parametric Design" Mathematics 10, no. 3: 387. https://doi.org/10.3390/math10030387

APA Style

Cabeza-Lainez, J. M., Salguero-Andújar, F., & Rodríguez-Cunill, I. (2022). Prevention of Hazards Induced by a Radiation Fireball through Computational Geometry and Parametric Design. Mathematics, 10(3), 387. https://doi.org/10.3390/math10030387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop