Subscriber access provided by - Access paid by the | UCSB Libraries
Remediation and Control Technologies
Raytracing for Fluence Rate Simulations in Ultraviolet Photoreactors Yousra M. Ahmed, Mark Jongewaard, Mengkai Li, and Ernest R. Blatchley Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b06250 • Publication Date (Web): 29 Mar 2018 Downloaded from http://pubs.acs.org on March 31, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 29
Environmental Science & Technology
Raytracing for Fluence Rate Simulations in Ultraviolet Photoreactors
1 2
Yousra M. Ahmed†, Mark Jongewaard‡, Mengkai Li†§,, Ernest R. Blatchley III*, †∥
3 4 5
†
6
States;
7
‡
8
§
9
Environmental Sciences, Chinese Academy of Sciences, 18 Shuang-qing Road, Beijing 100085,
Lyles School of Civil Engineering, Purdue University, West Lafayette, Indiana 47907, United
LTI Optics, LLC, Westminster, Colorado 80021, United States; Laboratory of Drinking Water Science and Technology, Research Center for Eco-
10
China;
11
∥
12
Indiana 47907, United States;
Division of Environmental & Ecological Engineering, Purdue University, West Lafayette,
13 14
*
15
[email protected]; Lyles School of Civil Engineering, Purdue University, 550 Stadium
16
Mall Drive, West Lafayette, IN 47907 USA
Corresponding author: ph. 1-765-494-0316; fax 1-765-494-0395; email
17 18
ABSTRACT: The performance of photochemical reactors is governed by the spatial
19
distribution of radiant energy within the irradiated region of the reactor. Raytracing has
20
been widely used for simulation of lighting systems. The central hypothesis of this work
21
was that raytracing can provide accurate simulations of fluence rate fields within
22
ultraviolet (UV) photoreactors by accounting for the physical and optical phenomena that
23
will govern fluence rate fields in UV photoreactors. Raytracing works by simulating the
ACS Paragon Plus Environment
Environmental Science & Technology
24
behavior of a large population of rays emanating from a radiation source to describe the
25
spatial distribution of radiant energy (i.e., fluence rate) within a system. In this study,
26
fluence rate calculations were performed using commercial raytracing software for three
27
basic UV reactors, each with a single low-pressure Hg lamp. Fluence rate calculations in
28
the raytracing program were based on the formal definition of fluence rate, calculated as
29
the incident radiant power from all directions on a small spherical receptor, divided by
30
the cross-sectional area of that sphere.
31
raytracing can provide predictions of fluence rate in UV radiative systems that are close
32
to experimental measurements and the predictions of other numerical methods.
33
Key words: Raytracing; ultraviolet; fluence rate; photo-reactors.
The results of this study demonstrate that
34
35
■ INTRODUCTION
36
The Ultraviolet (UV) spectrum spans the wavelength range from 100-400 nm. This
37
spectrum has been sub-divided into four ranges that comprise vacuum UV radiation (100-
38
200 nm); UVC, often referred to as germicidal UV (200-280 nm); UVB (280-315 nm);
39
and UVA (315-400 nm). Disinfection applications generally involve UVC radiation due
40
to strong absorbance among nucleic acids and other biomolecules, especially in the
41
wavelength range of 250-270 nm,1 but UVB and UVA radiation also demonstrate
42
germicidal behavior.2,3 UV photoreactors are also commonly applied to accomplish
43
direct photolysis or as part of an advanced oxidation process.
44
The most common terms used to quantify the spatial distribution of radiant energy of
45
UV radiation in UV photoreactors are irradiance (E) and fluence rate (Eo). Irradiance is
ACS Paragon Plus Environment
Page 2 of 29
Page 3 of 29
Environmental Science & Technology
46
defined as the total radiant power on an element surface from all upward directions,
47
divided by the area of that element, expressed in units of power per unit area (e.g.,
48
mW/cm²).4 Fluence rate (Eo) is a more general term used to quantify incident UV
49
exposure and is defined as the total radiant power incident on an infinitesimally-small
50
sphere from all directions divided by its cross-sectional area, also expressed in units of
51
power per unit area.4 In the special case of a perfectly collimated UV source incident
52
perpendicular to a plane, E and Eo are identical.4
53
The local rate of a photochemical reaction is directly proportional to the local fluence
54
rate. More generally, process performance of photoreactors is governed by the “dose”
55
(fluence) of UV radiation received by the photochemical targets (expressed in energy per
56
unit area, e.g., mJ/cm2). The UV dose received by a waterborne/airborne target in a UV
57
photoreactor is influenced by exposure time and the fluence rate history of the target. In
58
turn, these parameters are influenced by the geometry of the reactor system, type and
59
location of the UV sources, and the optical characteristics of the media that comprise the
60
system.5 Integration of the fluence rate field with the fluid mechanical behavior of the
61
system allows estimation of the fluence (dose) distribution delivered by the reactor.6
62
Detailed, accurate simulations of fluid mechanics in photochemical reactors can be
63
accomplished through the application of computational fluid dynamics (CFD). Similarly,
64
accurate and detailed simulations of the fluence rate field could yield improvements in
65
the ability to predict the performance of UV photoreactors.
66
A common germicidal UV radiation source is a low-pressure (LP) mercury vapor
67
lamp, which emits most of its ultraviolet radiation at 254 nm, close to the optimal
68
wavelength for disinfection. For water treatment systems, LP lamps are usually housed
ACS Paragon Plus Environment
Environmental Science & Technology
69
in a quartz glass sleeve inside the reactor chamber to prevent direct contact of the lamp
70
and its electrical components with water, while allowing most of the UV radiation
71
generated by the lamp to be transmitted to the surrounding fluid. UV radiation emitted
72
from the lamp encounters numerous components within the system including the lamp
73
envelope, air in the annular space between the lamp and the surrounding quartz sleeve,
74
the quartz sleeve, water, other lamps, and the reactor walls (Figure S1). Each of these
75
system components is characterized by quantifiable optical attributes, all of which will
76
affect the spatial distribution of radiant energy (i.e., fluence rate field).
77
Fluence rate field predictions and measurements. Over the last three decades,
78
several numerical models have been developed to simulate Eo distributions inside UV
79
photoreactors. Among the most common models is point source summation (PSS).7 An
80
integrated form of the PSS model over an infinite number of point sources was later
81
introduced by Blatchley8 as the line source integration (LSI) model. These models
82
provide inaccurate estimates of fluence rate in the near-field (i.e., close to the radiation
83
source) because of the simplified geometry they represent, and because they do not
84
account for the effects of reflection or refraction. Bolton introduced the multiple segment
85
source summation (MSSS) model that incorporated much of the PSS model, in which the
86
lamp was modeled as a series of differential cylindrical segments where the radiation was
87
emitted normal to the cylinder surface and decreased as a cosine function with the angle
88
of beam emission.9,10 The commercial software UVCalc3D (Bolton Photosciences Inc.,
89
Edmonton, AB, Canada), which is based on the MSSS model and includes terms to
90
account for the effects of reflection and refraction, has been applied for simulating the
91
fluence rate distributions in UV reactors. The MSSS model also accounts for the spectral
ACS Paragon Plus Environment
Page 4 of 29
Page 5 of 29
Environmental Science & Technology
92
output of the lamp and absorbance of the surrounding fluid medium. Other numerical
93
approaches aimed at modeling the E0 field around Hg lamps include the Extensive Source
94
Volumetric Emission (ESVE) and Extensive Source with Superficial Diffuse Emission
95
(ESDE) models.11,12 In the ESVE model, it is assumed that point emitters are uniformly
96
distributed within the lamp volume that emit isotropically. The ESDE model assumes
97
that radiation is emitted in a diffuse manner by point sources that are uniformly
98
distributed within the lamp volume.
99
The Discrete Ordinate (DO), finite volume, and Monte Carlo methods represent three
100
different approaches for solving the radiative heat transfer equation (RTE).13 The RTE
101
describes the propagation of radiation through homogeneous or heterogeneous media by
102
accounting for absorption, emission, and scattering behavior.
103
The limitations of the numerical fluence rate models described above motivated
104
exploration of raytracing as a tool for simulation of fluence rate fields in UV
105
photoreactors.
106
from one or more sources as a large number of individual rays; the fate of radiant energy
107
along the path of each ray can be accurately simulated using the fundamental laws of
108
optics, including those that govern refraction, reflection, absorbance, and shadowing.
109
Raytracing has previously been applied by Lau et al.14 to simulate the UV fluence rate
110
field in an in-duct UVGI system of an HVAC system. Lau et al. modeled a HVAC UV
111
system and predicted UV irradiance within a test rig of three UV lamps within the
112
irradiated air medium. They defined two irradiance terms: planar irradiance and spherical
113
irradiance. Lau et al. predicted spherical irradiance by dividing the incident radiation on
114
spherical receptors (2.54 cm diameter) by the total area of sphere, and planar irradiance
Raytracing simulates the behavior of electromagnetic radiation emitted
ACS Paragon Plus Environment
Environmental Science & Technology
115
by predicting the incident radiation on a planar surface area. However, neither of these
116
quantities represent the true value of fluence rate, which is the relevant parameter for
117
most photochemical reactor systems.
118
The primary objective of this paper was to examine the ability of the raytracing
119
approach to provide reliable, accurate predictions of Eo in UV systems that are applied in
120
water treatment. The results of raytracing calculations were compared to other physical,
121
numerical, and chemical methods.
122
Direct measurements of fluence rate are necessary to validate simulation models.
123
UV radiometers and chemical actinometers are the most common experimental methods
124
for performing irradiance/fluence rate measurements.8,15 An experimental method of
125
particular interest in this work was the method sometimes referred to as spherical
126
actinometry, wherein a chemical actinometer, such as iodide/iodate, is enclosed in one or
127
more small, spherical quartz vessels.16,
128
formal definition of fluence rate, in that it accounts for radiation power imposed on a
129
small, spherical photon receptor from all directions. However, conventional radiometry
130
and actinometry both fail in close proximity to radiation sources because of geometric
131
restrictions and the variable response of radiometers to off-angle (i.e., non-perpendicular)
132
radiation, as well as the risk of saturation of the actinometer solution under high fluence
133
rate.12 In addition, the effects of reflection and refraction may be important in spherical
134
actinometry experiments.
17
This approach was developed to mimic the
135
The micro fluorescent silica detector (MFSD) of Li et al. represents an alternative to
136
conventional radiometry and actinometry.18 The small size of this detector (1 mm long
137
and 0.3 mm o.d.) and its ability to accept incident photons from virtually all angles allows
ACS Paragon Plus Environment
Page 6 of 29
Page 7 of 29
Environmental Science & Technology
138
it to accurately measure Eo at essentially any location in a UV photoreactor, including
139
near-source regions such as the gap between the lamp and the sleeve. The MFSD has
140
been used successfully to examine in-situ, real-time Eo distributions in reactors with
141
single or multiple LP lamps in different media and with a wide range of UV
142
transmittance.19 It has also been used to measure photon fluence rate distributions in
143
photoreactors with polychromatic UV sources (e.g., medium-pressure Hg lamps).20
144
For the work described herein, raytracing was applied to simulate fluence rate fields
145
for several single-lamp UV photoreactors. The results of raytracing simulations were
146
then compared with existing numerical and experimental results to evaluate its usefulness
147
as a fluence rate calculation method.
148
■ METHODS AND MATERIALS
149
Raytracing simulations in this work were conducted using the commercial software
150
Photopia (V2015, LTI Optics, Westminster, Colorado, USA).
Photopia uses a 3D
151
probabilistic raytracing approach to predict the location and direction of ray emanation,
152
and the direction of reflected and transmitted rays. The raytracing method involves
153
digital definition of the physical boundaries of a system, as well as the physical (optical)
154
characteristics of the system components.
155
population of rays, it is possible to accurately describe the distribution of radiant energy
156
within the system.
By simulating the behavior of a large
157
Photopia includes a built-in CAD system that allows the user to create and
158
manipulate source geometry, as well as the geometry of other system components with
159
which radiation will “react.” Physical (optical) properties of the source and other system
160
components are user-assigned.
ACS Paragon Plus Environment
Environmental Science & Technology
161
UV radiation modeling in Photopia. Low-pressure mercury (LP Hg) lamps are
162
produced from quartz tubes with no internal fluorescent coating. Figure 1 illustrates the
163
interactions of rays with other materials in the model that result from a single ray. As
164
shown, rays emanating from the UV arc within the Hg lamp encounter refraction/Fresnel
165
reflection reactions through the inner and outer surfaces of the lamp, the air/quartz (i.e.,
166
inner) surface of the sleeve, and through the quartz/water (i.e., outer) surface of the
167
sleeve. All lamps were assumed to emit radiation uniformly throughout the length of the
168
lamp arc; the assumption of uniform emission was based on photographic images of
169
output radiation from operating LP Hg lamps (see Figure S10-a).
170
Refractive surfaces in the simulated system (e.g., quartz, water, and air interfaces)
171
included inner and outer surfaces (i.e., front and back) and were defined by specifying the
172
indices of refraction on both sides of each refractive surface, along with the material's
173
extinction coefficient. Transmissive materials were assigned properties that allowed
174
transmission either by allowing radiation to pass through without changing its direction,
175
or scattering based on a measured bi-directional transmittance distribution function
176
(BTDF). Clear glass is an example of a transmissive layer for visible light. To allow
177
accounting for shadowing effects in the raytracing simulation, rays were blocked by
178
shadowing surfaces in the lamp model (e.g., support structures), thereby forcing these
179
rays to find a new emanation point.
180
Most fluence rate numerical models do not explicitly model the UV lamp arc; for
181
instance, the DO model does not include the lamp geometry, the air annulus, nor the
182
thickness of the quartz sleeve. Instead it simply considers the radiation from the sleeve
183
outer surface. UVCalc® approximates the UV lamp as a linear array of cylindrical
ACS Paragon Plus Environment
Page 8 of 29
Page 9 of 29
Environmental Science & Technology
184
segments emitting UV radiation from their edges. Although interaction of rays through
185
the quartz sleeve surface (i.e., air/quartz/water interface) is accounted for in the model, it
186
disregards the actual dimensions and physical features of the lamp.
187
Different approaches can be employed to simulate the emission of rays from a UV
188
source; to choose the most accurate representation of radiation emission, the effect of the
189
emission mechanism on the fluence rate field predictions was examined by the raytracing
190
model by comparing fluence rate predictions of three different emission patterns to
191
experimental measurements. Based on this examination, which is explained in more
192
detail in the SI (Section 2), it was observed that emission of rays starting from an arc
193
surface inside an actual lamp yielded the closest estimates to the measurements.
194
Spherical representation of fluence rate. Local fluence rate was calculated by
195
dividing all incident radiant power on a small receptor sphere by the cross-sectional area
196
of the sphere. The receiving sphere needs to be small enough to accurately represent the
197
physical interpretation of local fluence rate; however, small spheres require relatively
198
large numbers of rays to ensure interaction with the rays and good resolution and stability
199
of the probabilistic raytracing solution, hence requiring more computational time.
200
Therefore, sphere size was optimized to yield a stable, converged solution within a
201
reasonable computational time. Figure S4 illustrates a radial array of spheres of 1 mm
202
diameter arranged around a UV lamp. Fluence rate values at intermediate locations were
203
estimated by linear interpolation. The effect of choosing a different sphere size (i.e.,
204
radius) on Eo predictions at several locations from the raytracing simulations only caused
205
up to 10% difference, so a subjective decision was made to use a sphere diameter of 1
206
mm in the raytracing calculations (Figure S5).
ACS Paragon Plus Environment
Environmental Science & Technology
207
Raytracing parameters. Convergence on a stable fluence rate field calculation
208
depends on the use of an appropriate number of rays emanated from the source.
209
Raytracing solutions typically converge on a stable distribution as the number of rays
210
increase. The number of rays required to converge on a stable solution depends on the
211
number of lamps (radiation sources) in the model, as well as the size and geometry of the
212
computational domain. The effect of changing the number of rays on the uniformity of
213
the fluence rate calculations around a single UV lamp is illustrated in Figure 2 and
214
discussed in the SI (Section 4).
215 216
Figure 1. Ray paths resulting from a single ray that originates inside a LP Hg lamp,
217
surrounded by a quartz jacket.
218
ACS Paragon Plus Environment
Page 10 of 29
Page 11 of 29
Environmental Science & Technology
350
2
Fluence Rate (W/cm )
300
250
200
150 106 rays 100
107 rays 5 x 107 rays 108 rays
50
5 x 108 rays 0 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
219 220
Sphere Number Figure 2. Fluence rate around a single low-pressure Hg lamp as a function of number of
221
rays emitted. Each data curve was obtained at a different number of rays ranging from
222
106 to 5 x 108. The horizontal axis represents a numerical label (1-16) for sixteen
223
spherical receptors that were distributed radially around the lamp at a fixed distance of 15
224
mm from the center of the lamp.
225 226
Fluence rate calculations in water test reactor. Raytracing simulations were
227
carried out for a single-lamp UV reactor filled with water with a UV transmittance of
228
88% (based on a 1.0 cm optical path), as illustrated in Figure S6.
229
measurements and numerical modeling for this reactor were performed by Liu et al.10 as
230
part of their evaluation of commonly used fluence rate numerical models including
231
MSSS, UVCalc3D, RAD-LSI, DO, and the View Factor Models. Liu et al.10 also used a
ACS Paragon Plus Environment
Fluence rate
Environmental Science & Technology
Page 12 of 29
232
KI/KIO3 actinometer to measure fluence rate at specific locations in the water reactor by
233
encasing the actinometer solution in small, quartz spheres.
234
actinometer spheres are listed in Table S3.
The locations for the
235 (a)
(b)
236 (c)
237 238
Figure 3. 2D view of the CAD model in Photopia for (a) the water test reactor of Liu et
239
al.,10 (b) the glass reactor with single LP UV lamp of Li et al.,18 and (c) the annular
240
reactor of Duran etal.12 constructed for raytracing simulations.
241
ACS Paragon Plus Environment
Page 13 of 29
Environmental Science & Technology
242
The dimensions of the LP UV lamp used in this study were 28 cm length and 15 mm
243
outer diameter. The thickness of the lamp envelope was assumed to be 1 mm; the total
244
output power of the lamp at 254 nm was 16 watts. For the raytracing simulation,
245
transmissive spheres were positioned at the same locations as the actinometer spheres
246
(see Figure 3-a).
247
UVCalc3D, RAD-LSI, and DO models, and measurements by the KI/KIO3 actinometer.
Raytracing simulation results were compared to those of the
248
Fluence rate calculations in annular photo-reactor. Duran et al.12 applied finite
249
volume (FV) based models to predict near-field and far-field fluence rate profiles for a
250
single-phase annular reactor.
251
fluence rate field for this same reactor; the output results from the raytracing analysis
252
were compared to the experimental measurements (i.e., actinometers) and numerical
253
model simulation results based on linear and extensive source models. UV transmittance
254
of the water medium circulated inside the annular reactor was 98% (based on a 1.0 cm
255
optical path). Duran et al. applied linear source and extensive source emission models
256
with isotropic or diffuse radiation (i.e., LSDE, LSSE, ESDE, and ESVE) to solve the
257
radiative transfer equation (RTE) in connection with finite volume CFD simulations.12 In
258
addition, a modified ESVE model was developed to include reflection/refraction and
259
absorption of radiation at the air/quartz/water interface, as well as the absorption and re-
260
emission of photons by the lamp plasma (abbreviated as: ESVERA).12 Fluence rate
261
estimates from these numerical models were validated with experimental measurements
262
using the potassium ferrioxalate actinometer for the near field irradiance at the reactor
263
outer tube wall.
A raytracing simulation was conducted to predict the
ACS Paragon Plus Environment
Environmental Science & Technology
264
Figure S7 illustrates the actual geometry of the annular reactor used by Duran et al.,12
265
while Figure 3-c illustrates the CAD model that was built for the raytracing simulation.
266
Some components of the system were not included in the raytracing model for simplicity,
267
since they were assumed to have negligible impact on the fluence rate predictions. These
268
included the PVC jacket and the electrical tape. The rubber stoppers were included as a
269
black (opaque) surface ends with zero reflection.
270
Fluence rate calculations in glass reactor. The MFSD can capture radiation from
271
essentially all angles with high sensitivity and can be used to measure Eo
272
experimentally.18 The MFSD is a Ge-doped silica cylinder that is 1 mm long and 0.3 mm
273
wide. Li et al.18 used this detector to measure Eo in a single LP UV lamp glass reactor,
274
operating with air (100% UV transmittance) and water with 85% and 95% UV
275
transmittance (based on a 1.0 cm path length). Axial Eo profiles were measured along the
276
lamp with the MFSD at four radial distances (r) from the sleeve (r = 3, 8, 18, 28 cm) and
277
were numerically predicted using UVCalc3D.
278
The raytracing approach was used to re-calculate fluence rate profiles for the reactor
279
configuration described by Li et al.18 to allow comparisons with the MFSD measurements
280
and UVCalc3D predictions. The raytracing simulations included a UV lamp arc centered
281
inside a quartz envelope and housed in a quartz sleeve in the middle of a glass reactor.
282
The physical dimensions of all refractive surfaces were included; refractive indices (n) of
283
1.51, 1.38, and 1.00 at 254 nm were assigned for quartz, water, and air, respectively.21, 22
284
The reactor was open to air, so the effect of the air/water interface was accounted for.
285
Figures S8 and S3-b illustrate the actual design of the reactor and its corresponding CAD
286
schematic for the raytracing analysis, respectively.
ACS Paragon Plus Environment
Page 14 of 29
Page 15 of 29
287
Environmental Science & Technology
■ RESULTS AND DISCUSSION
288
Water test reactor. Figure 4 displays the raytracing fluence rate predictions at
289
several locations based on the DO, UVCalc3D, and RAD-LSI numerical model results,
290
along with measurements by local actinometry.
291
observed between the results of raytracing and experimental measurements at most
292
locations than the other numerical models. At most locations (1, 2, 4, 5, 7, 9, 10, 12, and
293
13) the raytracing results were closer to the actinometry measurements compared to other
294
models, except some points like 3, 6, 8, 11, and 14. Locations of 3, 6, 8, and 11 are
295
almost at the same X and Y distance from the lamp center, and close to the top and
296
bottom of the reactor. It was not clear from Liu’s study what material was present at the
297
top and the bottom of the water test reactor. For the raytracing simulation, a reflective
298
surface was assigned to the top and the bottom; therefore, there was a contribution of
299
reflection/refraction at these locations which might me be not have been present in the
300
actual system.
Generally, better agreement was
301
It was noted that points 1, 2, 4, and 5 are the closest points to the reactor wall where
302
the reflection effects tend to be most dominant; this explains why the raytracing
303
predictions and the experimental measurements were notably higher than UVCalc, RAD-
304
LSI, and DO, since these models didn’t account for wall reflection. Moreover, Points 12,
305
13, and 14 are the closest points to the lamp and oftentimes it is mentioned in literature
306
that actinometry vessels might oversaturate in the close vicinity of the lamp and in this
307
case the measured fluence rate by the actinometer is lower than it should be. 12
308
A potential source of error in the experimental measurements with the actinometer
309
could be that the actinometer spherical vessels do not perfectly mimic the fluence rate
ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 29
310
definitions. This is because actinometer spheres have a diameter that is substantially
311
larger than that of the receiving spheres used in the raytracing simulations.
312
reflection/refraction attributes of the quartz spheres will influence the behavior of UV
313
radiation at the quartz-water interfaces of the spheres. Given the errors associated with
314
the experimental measurements as were reported by Liu et al.10 which are indicated by
315
error bars in Figure 4, the results indicate potential advantages of raytracing for fluence
316
rate field modeling in UV photoreactors, relative to other modeling approaches.
Also,
9 DO UVCALC3D RAD-LSI Raytracing Actinometry
8
2
Fluence Rate (mW/cm )
7 6 5 4 3 2 1 0 1
2
3
4
5
6
7
8
9
10
11
12
13
14
317 318
Test Point Figure 4. Fluence rate values at test points, as calculated by DO, UVCalc3, RAD-LSI,
319
and raytracing, along with local (spherical) actinometry measurements. Error bars on
320
measurements represent standard deviation.
321
ACS Paragon Plus Environment
Page 17 of 29
Environmental Science & Technology
322
Annular photo-reactor. Figure 5 illustrates the axial fluence rate profiles at the
323
outer tube surface of the annular reactor. The predictions of the FV numerical models
324
were comparable to each other and relatively close to the experimental measurements,
325
with some noticeable differences towards the ends regions of the tube, which showed a
326
steep decrease in fluence rate. Raytracing predictions of Eo followed the same pattern as
327
the other numerical models and were similar to measurements reported by the authors.12
328
The raytracing simulation over-predicted the fluence rate in the region 12 cm < axial
329
position < 15 cm by about 5-25% relative to the actinometry measurements; however,
330
they still fell within the experimental error margins, as did the ESVERA model results,
331
which also accounted for reflection, refraction, and absorption of the lamp quartz
332
envelope, as well as the absorption/remission for the Hg plasma of the lamp.
333
The ability of these numerical models, including raytracing, to predict Eo (at x=0)
334
inside the annular space of the reactor was evaluated and the results are shown in Figure
335
6. Raytracing results were in agreement with the predictions of the modified ESVE
336
models, LSDE, and ESDE, which all assumed diffuse emission from the lamp. The
337
fluence rate predicted by raytracing was approximately 5% higher than the other models
338
for some locations near the central region of the annular space. This may have been
339
partly because of differences in the model geometry itself, since in the raytracing
340
approach the thickness of the lamp, quartz, and reactor tubes were included in the
341
simulation, as was the reflective steel lamp support structure.
342
differences will contribute to variations in the fluence rate predictions. Another factor
343
that could affect the actual fluence rate is absorption/re-emission of radiation within the
344
lamp plasma, which was not accounted for in the raytracing simulations.
ACS Paragon Plus Environment
These geometric
In the
Environmental Science & Technology
345
raytracing calculations, the inner cavity of the lamp arc was assigned a material property
346
of zero absorption.
347
measurements and the ESVERA results (with at most 5% relative difference), which is
348
considered the most inclusive and accurate model among the FV numerical models.12
Generally, raytracing predictions agreed with the experimental
349 350
Figure 5. Fluence rate (irradiance) as a function of axial position measured along
351
the reactor outer wall (2.5 cm away from the lamp axis) of the annular photo-reactor.
352
Raytracing (red curve) vs. experimental measurements and other numerical models.12
353
Adapted from reference [12].
ACS Paragon Plus Environment
Page 18 of 29
Page 19 of 29
Environmental Science & Technology
354 355
Figure 6. Radial fluence rate profiles within the annular region of the UV reactor (at the
356
vertical axis of the lamp), as simulated by various models,12 including raytracing.
357
Adapted from reference [12].
358 359
Glass reactor. Estimated values of fluence rate for this reactor based on raytracing
360
were 5-9.5% and 10-40% higher for air and water, respectively, than the corresponding
361
measurements by the MSFD. Figures 7a and 7b provide comparisons of raytracing
362
simulation results and MFSD measurements for the reactor filled with air and water,
363
respectively. One reason for the discrepancy between simulated and measured values
364
was that the raytracing results were based on a spherical receptor, whereas the MSFD
365
involves a cylindrical detector. This was assessed by re-calculating fluence rate by
366
raytracing using cylindrical receptors of the same size as the MSFD detector (0.3 x 1
367
mm) and comparing these results with the predictions of fluence rate using spherical
368
receptors based on the theoretical definition of fluence rate (Figure S9).
ACS Paragon Plus Environment
Spherical
Environmental Science & Technology
369
receptors are able to receive radiation equally from all angles, hence it is expected that
370
the sphere can capture more rays than a cylinder that has a similar cross-sectional area.
371
The ratio of the Eo predicted using spherical receptors to the cylindrical receptors was
372
used as a shape factor to account for the difference between the raytracing results and the
373
experimental measurements. It should also be noted that the cylindrical receptors were
374
located parallel to the lamp, similar to the detector orientation during the experiments.
375
Qiang et al.23 showed that the vertical angular response of the MFSD deviated somewhat
376
from the ideal cosine angular response and this was corrected by applying an angular
377
response correction factor of 1.08 to the MFSD readings. This factor was determined
378
from the average fluence rate ratio of the parallel to perpendicular detector orientations
379
with respect to the lamp axis.
ACS Paragon Plus Environment
Page 20 of 29
Page 21 of 29
Environmental Science & Technology
250
(a)
Ray Tracing MFSD
2
Fluence Rate (mW/cm )
200
150
100
50
0 0
5
10
15
20
25
30
35
Radial Distance from Sleeve (mm)
380 200
(b) Ray Tracing UVT = 95% MFSD UVT = 95% Ray Tracing UVT = 85% MFSD UVT = 85%
2
Fluence Rate (mW/cm )
160
120
80
40
0 0
5
10
15
20
25
30
35
Radial Distance from Sleeve (mm)
381 382
Figure 7. Radial distribution of fluence rate: (a) and (b) fluence rate in air and water
383
medium, respectively. Raytracing simulation results are compared with measurements in
384
both panels.
ACS Paragon Plus Environment
Environmental Science & Technology
385
Figures 8 and 9 illustrate the results of raytracing simulations as radial and axial
386
distributions of fluence rate, respectively for the glass reactor, after being modified by
387
receptor shape and angular response factors against UVCalc3D results and experimental
388
measurements. After considering the factors described above, it appears that raytracing
389
can provide accurate estimates of fluence rate spatial distribution, as measured by the
390
MFSD. Raytracing estimates were within 0.7-14% of the measured values for air and 0.6
391
to 13% for water, except around the vertical end regions of the lamp. Neither of the
392
corrections explained the deviations of fluence rate values at an axial distance of 12 cm
393
or more from the lamp length center which reached up to 400% over prediction,
394
especially at the closest radial location to the lamp (3 mm radial distance) as seen in
395
Figure 9. This is discussed in more detail in the SI (section 7).
ACS Paragon Plus Environment
Page 22 of 29
Page 23 of 29
Environmental Science & Technology
250
(a) Air
Measurements Raytracing with correction
2
Fluence Rate (mW/cm )
200
150
100
50
0 0
5
10
15
20
25
30
35
Radial Distance from Sleeve (mm)
250
(b) Water
Measurements T=95% Measurements T=85% Raytracing with Correction, T=95% Raytracing with Correction, T=85%
2
Fluence Rate (mW/cm )
200
150
100
50
0 0
396
5
10
15
20
25
30
35
Radial Distance from Sleeve (mm)
397
Figure 8. Radial distribution of fluence rate: (a) and (b) fluence rate in air and water,
398
respectively, for the glass reactor based on raytracing with MFSD shape and angular
399
response correction factors included.
ACS Paragon Plus Environment
Environmental Science & Technology
(a) Air
250 2
Fluence Rate (mW/cm )
Page 24 of 29
200 150 100 50 0 0
20
40
60
80
100
120
140
160
Axial Distance from Lamp Center (mm)
(b) Water, UVT = 95%
2
Fluence Rate (mW/cm )
250 200 150 100 50 0 0
20
40
60
80
100
120
140
160
140
160
MFSD 3 mm UVCalc 3 mm Ray Tracing 3 mm MFSD 8 mm UVCalc 8 mm Ray Tracing 8 mm MFSD 13 mm UVCalc 13 mm RayTracing 13 mm MFSD 28 mm UVCalc 28 mm Ray Tracing 28 mm
Axial Distance from Lamp Center (mm)
(c) Water, UVT = 85%
2
Fluence Rate (mW/cm )
250 200 150 100 50 0 0
20
40
60
80
100
120
Axial Distance from Lamp Center (mm)
400 401
Figure 9. Comparison of fluence rate field axial distributions for glass reactor in (a) air,
402
(b) water (UVT 95%), and (c) water (UVT 85%). MFSD measurements as well as
403
UVCalc (solid curves) and raytracing (dashed curves, adjusted by shape and angular
404
response corrections) predictions are included.
ACS Paragon Plus Environment
Page 25 of 29
Environmental Science & Technology
405
The results of this study illustrate the ability of raytracing models to simulate the
406
spatial distribution of UV fluence rate in simple photochemical reactors. Raytracing
407
simulations provided predictions of fluence rate in basic reactor configurations with
408
single LP Hg UV lamps that agreed with chemical actinometer measurements and the
409
micro-fluorescent silica detector. The numerical (graphical) description of source and
410
reactor geometry, along with the ability to assign measurable optical properties to
411
materials that comprise the simulation domain allow for accurate simulations of fluence
412
rate distributions within UV systems. Specific reactor components that were accounted
413
for in these simulations included the UV lamp, quartz sleeve, reactor geometry/walls, and
414
fluid media. As such, raytracing tools may be applied as an alternative to other UV
415
irradiance numerical models, including LSI, MSSS, UVCalc, ESVE, and DO.
416
UV lamps were represented in these simulations as perfect arc cylinders enclosed in
417
a cylindrical quartz envelope. UV radiation was allowed to emanate uniformly over the
418
surface of the internal lamp arc. Future applications of raytracing may account for
419
variations in the radiation output, based on physical measurements of lamp output. These
420
models could be modified to account for non-uniformity in lamp output due to the effects
421
of lamp aging and fouling.
422
All models presented herein addressed simple single-lamp reactors; however,
423
multiple-lamp UV systems can easily be incorporated since the fluence rate is spatially
424
additive.
425
■ ASSOCIATED CONTENT
426
Supporting Information. Additional text, three tables, and twelve figures, describing
427
schematic illustration of basic photoreactor, emission pattern and modeling approaches of
ACS Paragon Plus Environment
Environmental Science & Technology
428
UV radiation, variation in fluence rate predictions according to the size of the spherical
429
receptors and the initial source of rays used, schematic illustrations of the simulated
430
photoreactors with single UV lamp, comparison between fully transparent spherical
431
receptors and cylindrical receptors, radiant intensity distribution plots, and lastly
432
comparison of fluence rate predictions from fully transparent spheres versus fully
433
absorbent spheres. This material is available free of charge via the Internet at
434
http://pubs.acs.org.
435
■ AUTHOR INFORMATION
436
Corresponding Author
437
*
438
Notes
439
The authors declare no competing financial interest.
440
■ ACKNOWLEDGEMENTS
Phone: 1-765-494-0316; e-mail:
[email protected].
441
The authors gratefully acknowledge financial support from the Citizen’s Energy, the
442
Water Environment Research Foundation (WERF), and the Edward M. Curtis Visiting
443
Professorship in the Lyles School of Civil Engineering at Purdue University.
444
■ REFERENCES
445 446 447 448
(1) Wastewater Technology Fact Sheet Ultraviolet Disinfection; US Environmental Protection Agency: Washington, DC, 1999. (2) Meierhofer, R.; Wegelin, M., Solar Water Disinfection - A guide for the application of SODIS. 2002.
ACS Paragon Plus Environment
Page 26 of 29
Page 27 of 29
449 450 451 452
Environmental Science & Technology
(3) Mbonimpa, E. G.; Vadheim, B.; Blatchley III, E. R., Continuous-flow solar UVB disinfection reactor for drinking water. Water Res. 2012, 46 (7), 2344-2354. (4) IUPAC, Compendium of Chemical Terminology Gold Book. International Union of Pure and Applied Chemistry: 2014; Vol. 2.3.3.
453
(5) Ultraviolet Disinfection Guidance Manual for The Final Long Term 2 Enhanced
454
Surface Water Treatment Rule; US Environmental Protection Agency: Washington, DC,
455
2006.
456 457 458 459
(6) Chiu, K.; Lyn, D. A.; Savoye, P.; Blatchley III, E. R., Integrated UV Disinfection Model Based on Particle Tracking. J. Environ. Eng. 1999, 125 (1), 7-16. (7) Jacob, S. M.; Dranoff, J. S., Light Intensity Profiles in a Perfectly Mixed Photoreactor. AlChE J. 1970, 16(3), 359-363.
460
(8) Blatchley III, E. R., Numerical Modelling of UV Intensity: Application to
461
Collimated-Beam Reactors and Continuous-Flow Systems. Water Res. 1997, 31 (9),
462
2205-2218.
463
(9) Bolton, J. R., Calculation of Ultraviolet Fluence Rate Distributions in an Annular
464
Reactor: Significance of Refraction and Reflection. Water Res. 2000, 34 (13), 3315-3324.
465
(10) Liu, D.; Ducoste, J. J.; Jin, S. S.; Linden, K., Evaluation of alternative fluence
466
rate distribution models. J. Water Supply Res. T. 2004, 53 (6), 391-408.
467
(11) Irazoqui, H. A.; Creda, J.; Cassano, A. E., Radiation Profiles in an Empty
468
Annular Photoreactor with a Source of Finite Spatial Dimensions. AlChE J. 1973, 19(3),
469
460-467.
ACS Paragon Plus Environment
Environmental Science & Technology
470
(12) Duran, J. E.; Taghipour, F.; Mohseni, M., Irradiance modeling in annular
471
photoreactors using finite-volume method. J. Photochem. Photobiol., A: Chemistry 2010,
472
215 (1), 81-89.
473
(13) Kowalski, W. J.; Bahnfleth, W. P.; Witham, D. L.; Severin, B. F.; Whittam, T. S.,
474
Mathematical Modeling of Ultraviolet Germicidal Irradiation for Air Disinfection.
475
Quantitative Microbiology 2000, 2 (3), 249-270.
476
(14) Lau, J.; Bahnfleth, W.; Mistrick, R.; Kompare, D., Ultraviolet Irradiance
477
Measurement and Modeling for Evaluating The Effectiveness of In-Duct Ultraviolet
478
Germicidal Irradiation Devices. HVAC&R Research 2012, 18 (4), 626-642.
479 480 481 482
(15) Kuhn, H. J.; Beaslavsky, S. E.; Schmidt, R., Chemical Actinometry. IUPAC 2004, 76 (12), 2105-2146. (16) Rahn, R. O., Potassium Iodide as a Chemical Actinometer for 254 nm Radiation: Use of lodate as an Electron Scavenger. Photochem. Photobiol. 1997, 66 (4), 450-455.
483
(17) Rahn, R. O.; Stefan, M. I.; Bolton, J. R., The Iodide/Iodate Actinometer in UV
484
Disinfection: Characteristics and Use in Determination of the Fluence Rate Distribution
485
in UV Reactors.
486
Association, Denver, CO, 2000.
Water Quality Technology Conference, American Water Works
487
(18) Li, M. K.; Qiang, Z. M.; Li, T. G.; Bolton, J. R.; Liu, C. L., In situ measurement
488
of UV fluence rate distribution by use of a micro fluorescent silica detector. Environ. Sci.
489
Technol. 2011, 45 (7), 3034-3039.
490
(19) Li, M.; Li, W.; Qiang, Z.; Blatchley III, E.R., On-Site Determination and
491
Monitoring of Real-Time Fluence Delivery for an Operating UV Reactor Based on a
ACS Paragon Plus Environment
Page 28 of 29
Page 29 of 29
Environmental Science & Technology
492
True Fluence Rate Detector, Environmental Science & Technology 2017, 51 (14),
493
8094−8100.
494
(20) Li, M.; Qiang, Z.; Wang, C.; Bolton, J.R.; Blatchley III, E.R., Experimental
495
Assessment of Photon Fluence Rate Distributions in a Medium-Pressure UV
496
Photoreactor,” Environmental Science & Technology 2017, 51, 3453−3460.
497 498 499 500
(21) Hale, G. M.; Querry, M. R., Optical Constants of Water in the 200-nm to 200-µm Wavelength Region. J. Opt. Soc. Am. 1973, 12 (3), 555-563. (22) Malitson, I. H., Interspecimen Comparison of the Refractive Index of Fused Silica. J. Opt. Soc. Am. 1965, 55 (10), 1205-1209.
501
(23) Qiang, Z.; Li, M.; Bolton, J. R.; Qu, J.; Wang, C., Estimating the Fluence
502
Delivery in UV Disinfection Reactors using a 'Detector-Model' Combination Method.
503
Chem. Eng. J. 2013, 233, 39-46.
ACS Paragon Plus Environment