Next Article in Journal
Past and Future Responses of Soil Water to Climate Change in Tropical and Subtropical Rainforest Systems in South America
Next Article in Special Issue
Data-Driven Air Quality and Environmental Evaluation for Cattle Farms
Previous Article in Journal
Investigating the Diurnal Variation in Coastal Boundary Layer Winds on Hainan Island Using Three Tower Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes

1
Savannah River National Laboratory, Aiken, SC 29808, USA
2
Rochester Institute of Technology, 54 Lomb Memorial Drive, Rochester, NY 14623, USA
*
Authors to whom correspondence should be addressed.
Atmosphere 2023, 14(4), 754; https://doi.org/10.3390/atmos14040754
Submission received: 22 March 2023 / Revised: 11 April 2023 / Accepted: 18 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Remote Sensing Techniques in Air Pollution Assessment)

Abstract

:
Using a one-dimensional code, we computed the power (enthalpy discharge rate) of a twelve-cell mechanical draft cooling tower (MDCT) using over two hundred visible condensed water vapor plume volume measurements derived from images, weather data, and tower operating conditions. The plume images were simultaneously captured by multiple stationary digital cameras surrounding the cooling tower. An analysis technique combining structure from motion (SfM), a neural-network-based image segmentation algorithm, and space carving was used to quantify the volumes. Afterwards, the power output was computed using novel techniques in the one-dimensional code that included cooling tower exhaust plume adjacency effects implemented with a modified version of the entrainment function, weather data averaged from eleven stations, and fan operations at the times when plume volumes were measured. The model was then compared with the averaged observed power output, and it validated well with an average error ranging from 6 to 12%, depending on the meteorological data used in the simulations. This methodology can possibly determine power plant fuel consumption rates by applying visible imagery.

1. Introduction

Mechanical draft cooling towers (MDCT) transfer heat from industrial cooling water to the atmosphere via forced convection and evaporation driven by fans. They are typically arranged with multiple independent cells in a row, as seen in Figure 1a. Heated water is sprayed evenly into the top of the fill structure where most of the heat and mass transfer occurs, as illustrated in Figure 1b. The fill structure is made of hard plastic with a honeycomb geometry that optimizes heat and mass transfer from the water to the air. After passing through the fill, the cooling water falls into a common basin from which it returns to the industrial process. Water temperature may be regulated by changing the number of fans operating or, in some facilities, the speed of individual fans. Individual cells can operate in free convection mode (fan off) which reduces the enthalpy transfer rate.
The exhaust created by an MDCT is initially a forced jet which transitions to a buoyant plume higher up. An MDCT operating close to design specifications often produces a visible plume of water droplets that form when the saturated exhaust air mixes with cooler ambient air. Visible MDCT plume droplets grow to only 2 to 3 microns in diameter during their short existence. Drift eliminators prevent large droplets from leaving the tower, so the visible exhaust plume is almost entirely composed of very small droplets that start to form as the exhaust exits the cell. Additional entrainment of sub-saturated ambient air eventually causes the droplets to evaporate. The exhaust plume will continue to rise after all droplets have evaporated until entrainment and expansion in a stable atmosphere eliminate the plume’s buoyancy. The maximum height and volume attained by the visible plume is a strong function of ambient relative humidity, which produces a strong diurnal cycle in plume height. The droplets evaporate almost instantly when entrainment of ambient air causes the relative humidity of the air around the droplet to drop below 100%. For this reason, when the visible plume disappears, the droplets have completely evaporated.
Visible cooling tower plumes and the means of modeling them have been occasionally studied for the purposes of visual impact and the practical impacts of fogging and deposited water in relationship with the enthalpy in the atmosphere. Plume emissions from volcanic eruptions have been measured and quantified to estimate volcanic discharge, analogous to what we propose in this paper [1]. Cizek and Noicka studied how plume volumes were affected by the stack diameters as well as the temperatures of the ambient air and exhaust [2]. Three-dimensional computational fluid dynamics (CFD) studies of cooling tower plumes have appeared in recent literature. CFD modeling has been used to further investigate the plume volume, length, and width with and without a wind velocity, which has also inspired other studies on plume volume flux [3,4,5,6]. Takata et al. used a Reynolds-averaged Navier–Stokes model to predict wind effects on a rising visible plume from an MDCT in the absence of nearby buildings [3], while Fan et al. similarly modeled an MDCT near a complex building environment [7]. Chahine et al. simulated plume dispersals from four large, collocated natural draft cooling towers to determine the long-range visual impact of the facility [8].
Although completely realistic simulations of the plume originating from a multi-cell MDCT requires three-dimensional simulations (multiple plume shielding, plume merging, etc.), the difference in length scale between the size of the plume and the important features of the tower’s stacks become computationally expensive. On one hand, if the CFD model domain does not have the necessary spatial resolution, the code inaccurately simulates the entrainment and transformation from a forced jet to a buoyant plume. Conversely, CFD simulations of MDCTs with many cells whose exhaust jets are adequately resolved require millions of nodes, which is impractical for long-term monitoring of many sites due to their computationally expensive analysis. For that reason, this study used a one-dimensional model based on older literature with adjustments to account for the interactions of multiple adjacent plumes.
The one-dimensional code is based on theoretical and empirical models. The theoretical model developed by Wessel and Wisse was used to attempt to calculate the size of cooling tower plumes, followed by an additional study on the effective source height of both the cooling tower and the plume [9,10]. An empirical study of cooling tower plume lengths showed that plumes were shorter in the summer and longer in the winter, because average relative humidity is higher in the winter [11]. Based on the theoretical models and empirical equations, many one-dimensional models have simulated plumes from single or multiple cells, which was the main part of this work [12,13,14,15,16].
The application of plume modeling focuses on many of the most intensive users of fossil fuel, such as power stations, refineries, and chemical plants, which discharge significant and predictable portions of their combustion heat to cooling towers. These towers often produce visible emissions, which are easily observable from commercially available satellite imagery. With a means to quantify the plume size and its relationship to the enthalpy discharge rate, such imagery could monitor the accuracy of CO2 emission reporting from these industrial sectors. This could complement direct CO2 concentration measurements from satellites with site-specific data [17,18].
This paper describes our initial research, which used a one-dimensional entraining plume model that incorporates condensation and evaporation to simulate MDCT exhaust plumes. The one-dimensional model also includes the effects of plume overlap in multi-cell MDCTs and wind on the entrainment of ambient air. This analysis is supported by an experimental effort to obtain a large and unique collection of data on plume volumes developed by our collaborators [19]. The plumes were first captured using eight stationary ground digital cameras that surrounded a cooling tower and shot thousands of images. The plume volumes were then measured using the Region Based Convolution Neural Network (Mask R-CNN) [20], a model that is based on extending the framework Faster R-CNN and used for segmentation. It has also been previously applied in a variety of topics ranging from medical imaging [21,22], to wildfire smoke [23], and gas plume detection [24]. Therefore, this model was extended to this study by identifying the plume boundary from a series of photos with manually annotated data, followed by a space carving technique that calculated plume volumes. For validation, this study considered 289 plume volumes, plant operational data, and weather data averaged from several weather stations. The next steps in future research will investigate the potential of three-dimensional modeling to aid in our understanding of interactions between adjacent plumes when the cooling tower cells are in a row. Artificial Intelligence (AI) methods will be used to compensate for some of the limitations of one-dimensional models.

2. Methods and Modeling

A one-dimensional model that integrates the conservation equations for mass, momentum, and energy along the plume axis (Figure 2) simulated a twelve-cell cooling tower where the cells are identical (Figure 1a). The cooling tower studied in this research is a GEA model 545436-12I-33-FCF that discharges 400 to 700 MW waste heat from a natural gas-fueled power plant in the southeastern United States that generates 500 to 900 MW of electricity. The cooling water line that connects the power plant to the cooling tower has a constant flow rate of 13.9 m3/s and is 6 to 11 °C warmer than the cooling water exiting the MDCT. Furthermore, the Hudson Tuf-Lite APT-32K-12 fans all run at a single RPM and the fan blades are set at a uniform pitch of 9 degrees. The twelve plumes in Figure 1a did not merge before they evaporated, but plume mergers do occur in calm conditions with high ambient relative humidity. The one-dimensional model simulated a row of individual cooling cells, where the height above the fill was 5.0 m. Each of the twelve cooling tower cells had a radius of 5.0 m and an area-averaged exhaust velocity of 6.0 m/s. The cell edges were 6.5 m from each other. Coordinate positions for each cell were used to combine individual plume simulations into a three-dimensional volume.
The model initially assumes a top-hat profile within the plume at a given elevation. A constant wind speed and direction, a flat terrain, and no obstacles such as buildings, vehicles, or trees were also assumed. The rate at which the cooling tower transfers thermal energy to the atmosphere via the cooling tower is given by:
Q l i q = ρ w c p w F c w T c w i n T c w o u t ,
where Q l i q is the rate at which heat is removed from the process cooling water and transferred to the air passing through the cooling tower, ρ w is the density of water, c p w is the specific heat capacity of water, F c w is the process cooling water flow rate, T c w i n is the process cooling water temperature entering the cooling tower, and T c w o u t is the process cooling water temperature exiting the cooling tower. The variables F c w , T c w i n , and T c w o u t cannot be measured remotely, so Equation (1) cannot be used to find Q l i q . The corresponding energy budget equation for the air passing through the cooling tower remotely estimates Q a i r , which is equal to Q l i q , aside from small radiative and conductive losses that occur as the air and water pass through the tower.
The entraining plume model requires lower boundary conditions to start the simulation. The required information is the enthalpy flux (sensible plus latent heat transfer rate to the atmosphere), MDCT exhaust volumetric airflow, and radius. The exhaust air is assumed to be saturated at the exhaust temperature. The following equation gives the starting conditions for the plume model:
Q a i r = ρ a C p a w π R c t 2 T p l T a i r + ρ a L w π R c t 2 q p l q a i r ,
where Q a i r is the enthalpy transfer rate (power), ρ a is the density of air, C p a is the specific heat of air, w is the updraft velocity, R c t is the exhaust radius, L is the latent heat of vaporization, T a i r is the ambient air temperature, T p l is the exhaust air temperature, q p l is the specific humidity of exhaust, and q a i r is the specific humidity of ambient air. Equation (2) used increasing values of the exhaust air temperature until the enthalpy transfer rate was large enough to produce a simulated plume equal in size to the measured plume size.
The standard set of the one-dimensional entraining plume model equations was also taken from Carazzo et al. Conservation equations for water vapor and liquid water were added to the Carazzo et al. equations [1]:
d V d z = 2 α e ,
d M d z = B V M ,
d B d z = V Γ a i r g T a i r ,
d T p l V d z = d V d z T a i r ,
d q p l V d z = d V d z q a i r + E v a p C o n d ,
d q l i q V d z = C o n d E v a p ,
d v p l V d z = v a i r d V d z .
In the equation set above,
V = w R 2 ,
M = w 2 R 2 ,
B = g w R 2 ,
where R is the radius of the plume, Γ a i r is the ambient air potential temperature lapse rate, v a i r is the horizontal ambient wind speed, v p l is the horizontal component of the plume air velocity, q is the specific humidity, and the subscripts pl, air, and liq are for the exhaust, ambient, and liquid phases (visible droplets), respectively. Potential temperature accounts for differences in static air pressure density with height. It is incorporated into Γ a i r by adding the dry adiabatic lapse rate (9.73 °C/km) to the actual lapse rate. The equation set was then solved from the bottom, starting with the boundary conditions at the top of the tower. The vertical discretization was 0.1 m. A predictor–corrector method solved the finite-difference equations, and the iterations were updated until the convergence criterion was obtained (plume temperature changes less than 0.001 °C in an iteration). A separate convergence criterion was applied to updraft momentum: the change in momentum in the most recent iteration relative to the previous iteration was required to be less than 0.01% of the previous momentum. The plume simulations are repeated with increasing initial exhaust air temperatures until the enthalpy discharge rate is sufficient to produce a simulated condensed water vapor (visible) plume with a volume equal to the measured volume. The air coming out of the cooling tower is assumed to be saturated.
The entrainment function, αe, models entrainment of the ambient air for conditions where there was wind from the data set [25]:
α e = R ϕ v p l 2 + w 2 v a i r c o s θ + β v a i r s i n θ .
In Equation (13), ϕ and β are entrainment coefficients with values of 0.085 [1] and 0.71 [25], respectively. The variable θ is the elevation angle of the plume axis relative to a horizontal surface. The first term in Equation (13) models the entrainment caused by velocity shear in the direction of the plume axis. The second term models the entrainment caused by velocity differences in the direction transverse to the plume axis.
The saturated plume exiting the top of the cooling tower will initially condense water vapor into small droplets with very small fall velocities, which are carried upward in the plume. As discussed above, some liquid water will be in the cooling tower updraft before it reaches the top of the stack. Unless the atmosphere is saturated, these droplets eventually evaporate. Condensation adds sensible heat to the plume and evaporation removes it. A realistic simulation of the plume includes droplet condensation/evaporation because a plume with condensed water vapor has higher temperature and greater buoyancy than the same plume without condensation. The evaporation/condensation in the plume is therefore calculated through an iteration loop until the water mass and thermal energy are conserved.
Visible plumes composed of water droplets initially increase in width as they rise but then contract before evaporating entirely. This implies that liquid water density is significantly greater in the core of the plume than in the periphery. Laboratory measurements by Kotsovinos demonstrated that the radial profiles of temperature, velocity and other plume parameters fit a normal probability (Gaussian) function [16]. This means that the inner core of an entraining turbulent plume created by MDCT exhaust will have higher temperatures, specific humidity (gas and liquid), and velocities than the edges.
The clear implication of the radial Gaussian profiles in turbulent jets is that the turbulent mixing of properties from the core of the plume to the adjacent environmental air is accomplished by randomly distributed turbulent eddies. The random fluctuations in plume properties will average out, resulting in the well-known Gaussian plume. It is expected that individual model predictions differ significantly from any individual plume, but the average model prediction for many plumes should be close to the average measured values.
The model described above produces single average values of plume properties at a given elevation. Since the actual plume variables will decrease with the distance from the center, it is expected that the predicted model average will be lower than the core values, including liquid water content, which means that the model should, on the average, under-predict the maximum visible plume elevation. Application of the Gaussian profile should correct this. The plume radial temperature and water concentration profiles are then defined by:
q p l m = q m a x q a m b exp r 2 2 σ 2 + q a m b ,
T p l m = T m a x T a m b exp r 2 2 σ 2 + T a m b .
Define:
q = q m a x q a m b ,
and
q p l m ¯ = 0 R p l m q p l m r d r R p l m = q R p l m 0 R p l m e x p r 2 2 σ 2 d r + q a m b ,
where R p l m = 2 σ and is converted to the standard error function erf r = 2 π 0 r e x p t 2 2 d t , t = r 1 2 R p l m where r = R p l m , and where q m a x = 1.67 q p l m ¯ q a m b + q a m b with corresponding equation for T a m b . As a condition, a plume is considered visible until q m a x < q s a t .
It is noted that some MDCTs, including the site of our empirical measurements, have ten or more cells in a row. Individual cell plumes in a row with many cells will not develop as they would if they were far from the other plumes in the row because as the plumes rise and spread, they will mix and interact. Relative to a single isolated cell, multiple plumes shield each other, to some degree, from entrainment of drier environmental air. The shielding will cause the maximum heights of plumes in a row, on the average, to be somewhere between the height that an equivalent area single-cell plume would reach and the height that an individual cell in the row would reach if it was isolated from the rest of the cells and their plumes.
The following approach was developed to extend the one-dimensional code to account for plume interactions: run the one-dimensional model as a single cell with an area equal to the combined areas of all cells and account for the increased entrainment experienced by the row of smaller plumes relative to a single large plume. Generally, the semi-empirical, heuristic extension of the one-dimensional model should depend only on the basic geometric features of an MDCT, such as cell number, diameter, and spacing. All 12 cells of the MDCT have radii of 5.0 m ( r 0 ), and the combined area of these 12 cells is equal to a single cell with a radius, r 00 , of 17.3 m:
r 00 = N c r 0 r 0 ,
where Nc is the numbers of cells (12). In this approach, the simulation starts with the plume radius r = r 00 . The twelve cells have a greater combined perimeter than the equivalent single cell: p 1 = N c 2 π r 0 is larger than p 2 = N c 2 π r 00 . Since the twelve plumes have a larger initial “boundary” with the surrounding atmosphere, the entrainment of ambient air into the 12 cells will be greater than the entrainment into the equivalent single cell. However, since the cells are arranged in a row, the additional entrainment will be a function of the wind direction relative to the orientation of the axis of the cell row. There exists shielding of downwind cells when the wind direction is parallel to the cell axis, and less when it is normal to the cell axis. To account for this effect, several parameters are defined. The first is a radius equal to the radius of a plume with the combined areas of the Nc plumes:
r e q = N c r s h r s h ,
where r s h is the radius of an individual cell plume.
The next parameter is θ d , which is the difference between the wind direction and the azimuth of the row of MDCT cells relative to north. A parameter, b 1 , is next defined as:
b 1 = π r s h + d l s i n θ d / π r s h + d l ,
where d l is a length scale for the row of cells:
d l = π N c 1 r s h + N c 1 C s p .
In this equation, C s p is the gap between the edges of adjacent cell shrouds as shown in Figure 3. The parameter b 1 is designed to maximize entrainment when the wind is blowing perpendicular to the row of MDCT cells and minimizes it when the wind is blowing parallel to the row of MDCT cells. The next parameter, a 1 , is designed to account for the additional entrainment that a row of plumes will experience relative to a single cell that started with the equivalent area to the Nc cells combined:
a 1 = 2 π N c r s h 2 π r e q = N c .
Equation (13) is therefore modified by incorporating a 1 and b 1 :
α e m = R a 1 ϕ v p l 2 + w 2 v a i r c o s θ + b 1 β v a i r s i n θ .
We note that inside water-cooling towers, there can exist maldistribution of water flows as demonstrated by previous authors [26,27], but at the twelve-cell water cooling tower site in the southeastern U.S., the enthalpy transfer rates for all twelve cells are nearly identical with variations of 2 to 3% except for the tenth cell, which was 8% higher than the average as seen in Figure 4. This indicates that on a cell-by-cell basis, the thermal energy transferred from the cooling water to the air flowing through the cooling tower is the same to within a few percent. Although the figure shows slightly unequal distribution of the flow in the tower, it is small and assumed uniform for the purpose of this study. The MDCT analyzed in this paper operates at a humid subtropical site in the southeastern U.S., which has very low average wind speeds (~1 m/s). This model can also be used for winter and summer applications, but it does not account for the impact of freezing conditions on plume thermodynamics and does not simulate precipitation, which can be produced by very large plumes in cold weather.

3. Plume Volume Measurements

Plume volumes were measured using imagery from 8 stationary camera stations deployed around the cooling tower as illustrated in Figure 5. Because off-the-shelf three-dimensional reconstruction tools, based on Structure from Motion, are poorly suited to the plume’s indistinct features and changing shape over time, a novel analysis was developed to determine plume volumes. This technique, detailed by our collaborators in Connal et al., uses a neural network-based image segmentation algorithm to identify the plume in each image and a space carving technique to estimate the volume of the plume [19]. Cameras were each calibrated, and camera pose (location and orientation of the installed camera) was determined using Structure from Motion on the solid structures in the scene. A model, based on Matterport Mask R-CNN, was trained to identify the plume in each image, using manually annotated data [28]. As visualized in Figure 6 and in Connal et al.’s study, the image segmentation model determines both a boundary box for the plume and classifies the pixels as plume or not-plume within the boundary box [19]. Different results of Mask R-CNN performances were also displayed in Connal et al.’s study [19]. Finally, the plume volume was estimated by dividing the domain into voxels and “carving” away any portion of the domain on a line of sight identified as not-plume from the image segmentation in any of the 8 camera angles, also seen in Figure 5. The project has also collected imagery for plume volume measurement via unmanned aerial vehicle flights, but ground station data have been used for plume model development due to the greater number of samples and wide variety of weather conditions seen during the continuous field study. A total of 289 data points—each based on analysis of images from all of the 8 ground station cameras—were used in this study. The data yielded plume volumes from 1274 to 2,247,000 m3.
The 8 camera stations deployed around the site each contained identical hardware consisting of a Canon M200 mirrorless camera, solar charging, a Raspberry Pi-based GPS timing and communications module, and weatherproof case. The stations were distributed around the perimeter of the cooling tower site to capture all sides of the plume and distances as far as practical within the open space around the cooling tower. Prior to deployment, the field of view from each camera was modeled in Blender [29], resulting in the selection of a Canon 10–18 mm ultra-wide-angle lens to obtain an acceptable height of plume capture.

4. Weather Data and Model Inputs

The one-dimensional model used averaged weather data from 11 regional airports. Some of these stations are in urban areas (Charleston and Columbia, SC, USA; Savannah, GA, USA; and Jacksonville, FL, USA), some are associated with military bases (Beaufort, SC and Hunter Field, GA, USA) and others are smaller, more rural airfields (Walterboro, Orangeburg, Barnwell, and Hilton Head, SC; and Waycross, GA, USA). The National Oceanic and Atmospheric Administration (NOAA) has primary responsibility for weather data collection, but management of smaller airfields is often the responsibility of the Federal Aviation Administration (FAA). An onsite weather station’s data were not used because the measurements were made near a floodplain forest which affected both relative humidity and wind speed. Airport weather stations are in open areas with no nearby buildings or trees and thus are more representative of conditions at the discharge level of an MDCT, which is typically 10 to 15 m.
Simulations of the measured plume volumes require input files for the weather and MDCT layout (number of cells, cell diameter, cell spacing) as well as entrainment coefficient data (Equation (23)). The main input file contains ten variables such as date/time, ambient air temperature, dewpoint temperature, wind direction (degrees from north), wind speed (m/s), air pressure, measured plume volume, number of cells with fans on, cell-area average exhaust velocity at 6.3 m/s, and vertical temperature gradient at −6.5 °C/km (the NOAA standard atmospheric lapse rate). The model validation data are average cell exhaust air temperature and measured power. Another input file includes entrainment coefficients and cooling tower geometric parameters such as number and spacing of cells.

5. Model Simulations Results and Discussion

For each of the 289 measured plume volumes, multiple one-dimensional model runs were performed using ascending exhaust air temperatures until the simulated combined plume volume exceeded the measured value. Each of the simulated plume volumes had a corresponding enthalpy discharge rate (power) to the atmosphere (exhaust air was assumed to be saturated). The simulated plume volume greater than the measured volume and the volume from the preceding iteration were used to linearly interpolate the corresponding powers to the value expected for the measured volume.
Figure 7 compares plumes created by combining multiple simultaneous digital photos to corresponding plumes created by combining multiple simulations of one-dimensional plumes into a single volume representing twelve cells. The measured plumes are highly irregular and appear to capture the denser parts of the plume but not the fringes, when compared to Figure 1. Eight cell fans were on when the digital imagery was taken and created the plume shown in Figure 7a,b. Figure 7c,d show images for a time when seven fans were on and five were off. The cells with fans off are obvious in Figure 7b,d because buoyant convection produces a much smaller plume relative to forced convection. There is some evidence of smaller buoyant plumes at the ends of the measured volume in Figure 7a,c, but they are much less obvious. The shapes of the simulated and measured plumes can only roughly agree because the simulations are from a one-dimensional code, and because the measured plume is irregular due to turbulent eddies in the wind that constantly deform the plume from the idealized simulated plume that was created using regionally averaged weather data.
The model was primarily validated based on the averaged computed and measured powers. Figure 8 shows that on an individual plume volume basis, the one-dimensional model could greatly over or underestimate observed power output. However, when the average meteorological data for all 11 stations were used as model input, the averaged power from the calculated and measured powers were 343 MW and 378 MW, respectively, giving an error of about −10%. The correlation coefficient for measured and computed pairs was 0.60 and the average absolute error was 103 MW. The large error for individual pairs of measured and computed powers is mostly attributable to the uncertainty in the meteorological data, specifically the dewpoint temperature depression (DPTD).
The one-dimensional model assumes a Gaussian distribution of properties (updraft velocity, temperature, liquid water content) around the centerline of the plume. It is well known that time-averaged turbulent plumes have Gaussian distributions of properties. So, it is appropriate to expect that the average of many measured plume volumes will be close to the average simulated volume, even though on an individual basis the simulated volume may be very different from the measured volume. This implies that the average computed power will also be closer to the average measured power than individual pairs of measured and computed powers. The larger standard deviation of the simulated powers in Figure 8 relative to the measured power standard deviation is the consequence of using a one-dimensional model that can only simulate Gaussian (average) plume shapes and volumes. Sometimes, the simulated powers will be greater than the measured powers and sometimes they will be smaller, as can be seen in Figure 8. If the regional meteorology is representative of the average weather conditions at the cooling tower, and the enthalpy flux from the cooling tower is correct, the average measured powers should be about equal to the average simulated powers because the average measured plume will be Gaussian, like the model plume.
Figure 9 shows how the one-dimensional model’s ability to predict average power varied according to which of the 11 regional airport meteorological data sets were used as code input. Results are also shown from simulations that used average weather data from all 11 stations and from the six more urban or military sites (Charleston, Columbia and Beaufort, SC, USA; Savannah and Hunter Field, GA, USA; and Jacksonville, FL, USA). The meteorological variable that is plotted is the DPTD at the time the digital imagery of the cooling tower plumes was collected. The DPTD was chosen because it has the greatest meteorological impact on plume volume (correlation coefficient of 0.79). The average ratio of computed-to-measured power for the six urban or military sites is 0.98 and the average ratio for all eleven sites is 0.91. The five more rural sites produced an average ratio of 0.79. The reason for the difference in model performance as a function of which stations were used in the averages is not known. It could be an effect of urban heat islands which would increase air temperature and lower dewpoint temperature, both of which increase DPTD and cause the model to require more power to produce the same volumes. Alternatively, the smaller, more rural airports may be more influenced by nearby forests which would lower air temperatures and raise dewpoint temperatures relative to the more urban airports.
The computed power ratios are less than 1.0 for all three average meteorological datasets (all 11 stations, 6 NOAA/DOD stations, 5 rural stations). This bias may be attributable to model deficiencies or the average plume volume measurements may be low. When Figure 1 is compared to Figure 7, it would appear feasible that the wispy edges of the visible plumes are not captured by the space carving technique. More plume volume measurements will be available in the future, which may make it possible to determine if, for the same meteorological conditions and power, plume images with sharp edges produce larger volumes than plume images with wispy, less distinct edges. This topic is discussed further in the next section.

6. Uncertainty Analysis

Several variables contribute to the uncertainty of the one-dimensional model predictions of enthalpy transfer rate (power) from the cooling tower to the atmosphere: measured plume volume, meteorological data, entrainment coefficients, and other limitations of a one-dimensional model. The DPTD is highly correlated with the computed power as shown in Figure 9 (R2 = 0.79), and this meteorological variable is the largest source of uncertainty in the weather data inputs. Correlation of the 11 individual station power ratios with dewpoint temperature depression in Figure 9 shows that a 1 °C error in average DPTD produces a 17% error in the average computed power. The averaged meteorology produced smaller errors than 17%, so it appears that the averaging produced an average DPTD error of less than 1 °C.
Data presented by Carazzo et al. shows a range from about 0.06 (buoyant convection) to about 0.13 (forced convection) for the entrainment coefficient ϕ in Equations (13) and (23) [1]. The one-dimensional model simulations described in this paper used a value of 0.085, which is intermediate between buoyant and forced convection but closer to the buoyant convection values. The value of 0.085 was chosen because the plumes are observed to transition rapidly from forced to free convection. A survey of experimental results for the other entrainment parameter (β) in Equations (13) and (23) found a wide range of values from 0.3 to 2.0 [25,30,31,32,33,34]. A sensitivity run was conducted with the one-dimensional code using a value of ϕ of 0.12 and a value of β of 1.0. Since these values are higher than the baseline values of 0.085 and 0.71, they increased entrainment of ambient air which required higher powers to produce the observed plume volumes. The ratio of computed to observed power rose from 0.91 to 0.97 when the 11-station average station meteorology was used.
When Equation (13) was used to simulate the 289 volumes instead of Equation (23), the computed-to-measured power ratio dropped from 0.97 to 0.91 (using average meteorology from the six NOAA or DOD stations). When meteorology averaged from all 11 stations was used, the power ratio dropped from 0.91 to 0.84. The power ratio dropped in both cases because entrainment of ambient air was reduced. This suggests that the alterations in Equation (23) to model adjacency effects of multiple plumes in a row moderately improved model performance by more correctly computing entrainment.
Figure 10 is an example of a ground-based digital image of a cooling tower plume with the mask for that plume outlined in green. Small areas of the plume that are separate from the main body of the plume are not within the mask boundary. Within the mask boundary some parts of the plume are more optically thick than others.
The uncertainty estimation of the plume volume measurements lies within the plume mask generation for boundary decision and space carving voxel density with uniform volumes for reconstruction. While considering the initial boundary of three-dimensional voxels, the individual volume of all voxels within this space can affect the resulting volume. Plume volume estimates were based on space carving voxel densities where a uniform size of voxels ranged between 0.1 and 1.0 m3. For plume mask generation, the model inferences produced by Mask R-CNN provide estimates of the plume boundary. For the uncertainty procedure, 60 plume images with the model inference mask were analyzed, in which 30 locations were sampled along the border of each mask. At each location, the distance in pixels from the edge of the model inference mask to the visible plume boundary was measured. A positive distance indicated the plume boundary was outside of the mask at that sample location, a negative distance for when the boundary fell within the mask, and zero when the two were aligned. Based on 1800 pixel distance measurements, the average distance from the model inference segmentation mask to the observed plume boundary was 11.2 pixels. A negative measurement would indicate that on average the model is overcompensating, but the result was positive, which reveals that the average inference mask needed to be expanded by approximately 11 pixels uniformly along the edge of the entire boundary. Figure 11 shows examples of mask boundary adjustments that were made to more correctly track the true plume boundary.
The volume uncertainty procedure uses the morphology operation of dilation to uniformly expand the model inference masks by 11 pixels and compare against the plume volumes produced by the model. This analysis was conducted on 50 cases, where the volumes produced by the dilated masks were treated as the ground truth. The volumes were calculated and averaged from the voxel densities ranging between 0.1 and 1.0 m3 for all cases. To uniformly expand the masks by 11 pixels, a 23-pixel by 23-pixel kernel was used to implement a dilation operation, which adds pixels to the boundary of the mask. By operating on pixels centered at the origin of the kernel, the mask is expanded uniformly in all directions within a distance of 11 pixels.
The differences between the original 50 plume volumes produced by Mask R-CNN and the dilated volumes are shown in Figure 12. The differences are expressed as percent increase in volume relative to the original volumes. The large percent increases are associated with the smallest volumes. This is attributed to the rapidly decreasing percentage of an 11-pixel surface layer on a volume relative to the total plume volume as it increases. The power law fit to the data in Figure 12 represents the uncertainty in the plume volume measurements expressed as percent increase in plume volume from the original undilated measurements.
The power law fit to the data in Figure 12 was used to increase the original 289 volume estimates. The one-dimensional plume model was run with the 289 modified volumes, which produced the average power estimates for different meteorological data sets as discussed and shown in the previous section of this paper. When the one-dimensional model was run with what is believed to be the best meteorological data (average of six NOAA and DOD stations) the average power computed from the undilated volumes was 94% of the measured average. The average power computed from the dilated volumes was 98% of the measured average. Based on this uncertainty analysis, it appears that the plume volume measurement uncertainty is a small contributor to the uncertainty of the average power estimate derived from the plume volume measurements and other model input variables.
The average rate of power generation is a more reliable indicator of fuel consumption than an estimate based on a single plume volume. The effect of random errors in the meteorological data and plume volume measurements largely average out when many plume volumes are used to compute the average power. This is clearly apparent in Figure 13, which shows the results of averaging the differences between pairs of measured and computed powers. Averages were computed by randomly selecting pairs from the 289 available. The procedure started by selecting single pairs 10,000 times and computing the average. Then two pairs were selected and averaged 10,000 times and the 10,000 averages were averaged. This process was repeated for three pairs, four pairs on up to 289 pairs. Figure 13 shows that 50 pairs produce an average error of 4.5% which is only slightly larger than the average error for all 289 pairs (2.4%).
The current uncertainty analysis of the plume volume, based on ground-based imagery and the method for volume quantification, is an initial investigation that may require additional comparisons for future research. For example, three-dimensional ray tracing combined with realistic water droplet distributions generated by three-dimensional simulations may produce more accurate plume boundaries with Mask R-CNN. The inclusion of a three-dimensional plume study could use a variety of mesh adaptations/distributions (i.e., a mesh-independent study) and additional physics-based models to accurately simulate the plume boundary.

7. Conclusions

A modified one-dimensional entraining plume model predicted average MDCT enthalpy discharge rate to the atmosphere (power) to within about 9% (343 MW computed versus 378 MW measured). These results were based on 289 measured plume volumes derived from multiple simultaneous digital images of the visible plumes. A one-dimensional plume model used average meteorology from 11 regional weather stations to predict cooling tower enthalpy dissipation rate (power). When average meteorology was computed from six urban or military sites the ratio improved to 0.98. The power ratio dropped to 0.79 when the average meteorology from five rural airfields was used. The correlation coefficient for individual pairs of measured and computed powers was 0.60. Since accurate average CO2 emissions are needed for assessment of potential climate impacts, the one-dimensional model average power determinations are of primary interest. While the power output is essential to investigate emissions, issues regarding power consumption needs to be further investigated in a future project paper. The results presented here appear to justify further research: three-dimensional modeling to determine the importance of fundamentally three-dimensional phenomena such as plume mergers, and AI (artificial intelligence) applications to improve model performance and improved meteorological input.

Author Contributions

Conceptualization, C.S., A.G. and B.d.; methodology, C.S., A.G. and B.d.; software, C.S., A.G., B.d., R.C. and S.A.; validation, C.S., A.G., R.C. and S.A.; formal analysis, C.S., A.G., B.d., R.C. and S.A.; investigation, C.S., A.G., B.d. and S.A.; resources, A.G. and B.d.; data curation, C.S., A.G., B.d., R.C. and S.A.; writing—original draft preparation, C.S.; writing—review and editing, C.S., A.G., S.A., R.C. and B.d.; visualization, C.S., A.G. and B.d.; supervision, A.G. and B.d.; project administration, B.d.; funding acquisition, B.d. All authors have read and agreed to the published version of the manuscript.

Funding

This work was produced by Battelle Savannah River Alliance, LLC under Contract No. 89303321CEM000080 with the U.S. Department of Energy. U.S. Government license to provide public access under the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan, accessed on 22 February 2023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The one-dimensional input files have the required input to run the code, which will become available upon request. This includes meteorology data from public sources and plume volumes estimated from imagery.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carazzo, G.; Kaminski, E.; Tait, S. On the rise of turbulent plumes: Quantitative effects of variable entrainment for submarine hydrothermal vents, terrestrial and extra-terrestrial explosive volcanism. J. Geophys. Res. Solid Earth 2008, 113, 1–19. [Google Scholar] [CrossRef]
  2. Cizek, J.; Nozicka, J. Cooling tower plume. In Proceedings of the 35th Meeting of Departments of Fluid Mechanics and Thermomechanics (AIP), Šamorín, Slovak Republic, 23 September 2016. [Google Scholar]
  3. Takata, K.; Michioka, T.; Kurose, R. Prediction of a visible plume from a dry and wet combined cooling tower and its mechanism of abatement. Atmosphere 2016, 7, 59. [Google Scholar] [CrossRef]
  4. Li, S.; Flynn, M.R. Cooling tower plume abatement and plume modeling: A review. Environ. Fluid Mech. 2021, 21, 521–559. [Google Scholar] [CrossRef]
  5. Li, S.; Moradi, A.; Vickers, B.; Flynn, M.R. Cooling tower plume abatement using a coaxial plume structure. Int. J. Heat Mass Transf. 2018, 120, 178–193. [Google Scholar] [CrossRef]
  6. Zargar, A.; Kodkani, A.; Vickers, B.; Flynn, M.R.; Secanell, M. A hybrid cooling tower model for plume abatement and performance analysis. Appl. Therm. Eng. 2023, 219, 119593. [Google Scholar] [CrossRef]
  7. Fan, J.; Dong, H.; Xu, X.; Teng, D.; Yan, B.; Zhao, Y. Numerical investigation on the influence of mechanical draft wet-cooling towers on the cooling performance of air-cooled condenser with complex building environment. Energies 2019, 12, 4560. [Google Scholar] [CrossRef]
  8. Chahine, A.; Matharan, P.; Wendum, D.; Musson-Genon, L.; Bresson, R.; Carissimo, B. Modelling atmospheric effects on performance and plume dispersal from natural draft wet cooling towers. J. Wind. Eng. Ind. Aerodyn. 2015, 136, 151–164. [Google Scholar] [CrossRef]
  9. Wessels, H.R.A.; Wisse, J.A. A method for calculating the size of cooling tower plumes. Atmos. Environ. 1971, 5, 743–750. [Google Scholar] [CrossRef]
  10. Hanna, S.R.; Swisher, S.D. A method for calculating the size of cooling tower plumes. Atmos. Environ. 1972, 6, 587–588. [Google Scholar] [CrossRef]
  11. Barber, F.R.; Martin, A.; Shepherd, J.G.; Spurr, G. The persistence of plumes from natural draught cooling towers. Atmos. Environ. 1974, 8, 407–418. [Google Scholar] [CrossRef]
  12. Winiarski, L.D.; Frick, W.F. Field Investigations of Mechanical Draft Cooling Tower Plumes; Corvallis Environmental Research Laboratory, Office of Research and Development, US Environmental Protection Agency: Corvallis, OR, USA, 1977; pp. 6–37.
  13. Wu, F.H.; Koh, R.C. Mathematical Model for Multiple Cooling Tower Plumes; Environmental Protection Agency, Office of Research and Development, Environmental Research Laboratory: Corvallis, OR, USA, 1978; Volume 1, pp. 4–59.
  14. Winiarski, L.D.; Frick, W.F. Cooling Tower Plume Model; US Environmental Protection Agency, Office of Research and Development, Corvallis Environmental Research Laboratory, Assessment and Criteria Development Division: Corvallis, OR, USA, 1976; pp. 5–39.
  15. Policastro, A.J.; Coke, L.; Wastag, M. User’s Manual: Cooling-Tower-Plume Prediction Code; No. EPRI-CS-3403-CCM; Argonne National Lab: Chicago, IL, USA, 1984.
  16. Kotsovinos, N.E. A Study of the Entrainment and Turbulence in a Plane Buoyant Jet; Doctor of Philosophy, California Institute of Technology: Pasadena, CA, USA, 1975. [Google Scholar]
  17. Bovensmann, H.; Buchwitz, M.; Burrows, J.P.; Reuter, M.; Krings, T.; Gerilowski, K.; Schneising, O.; Heymann, J.; Tretner, A.; Erzinger, J. A remote sensing technique for global monitoring of power plant CO2 emissions from space and related applications. Atmo. Meas. Technol. 2010, 3, 781–811. [Google Scholar] [CrossRef]
  18. Nassar, R.; Mastrogiacomo, J.; Bateman-Hemphill, W.; McCracken, C.; MachDonald, C.; Hill, T.; O’Dell, C.W.; Kiel, M.; Crisp, D. Advances in quantifying power plant CO2 emissions with OCO-2. Remote Sens. Environ. 2021, 264, 112579. [Google Scholar] [CrossRef]
  19. Connal, R.; Pines, W.; Borek, M.; Bauch, T.; Raqueno, N.; d’Entremont, B.; Garrett, A.; Salvaggio, C. Utilizing Mask R-CNN for automated segmentation of condensed water vapor plumes from multi-view imagery. In Proceedings of the Volume 11843 Applications of Machine Learning 2021, San Diego, CA, USA, 1 August 2021. [Google Scholar] [CrossRef]
  20. He, K.; Gkioxari, P.; Dollár, P.; Girshick, R. Mask R-CNN. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22 October 2017. [Google Scholar]
  21. Cai, L.; Long, T.; Dai, Y.; Huang, Y. Mask R-CNN-based detection and segmentation for pulmonary nodule 3D visualization diagnosis. IEEE Access 2020, 8, 44400–44409. [Google Scholar] [CrossRef]
  22. Anantharaman, R.; Velazquez, M.; Lee, Y. Utilizing mask R-CNN for detection and segmentation of oral diseases. In Proceedings of the IEEE International Conference on Bioinformatics and Biomedicine, Madrid, Spain, 3 December 2018. [Google Scholar]
  23. Gupta, T.; Liu, H.; Bhanu, B. Early wildfire smoke detection in videos. In Proceedings of the 25th International Conference on Pattern Recognition, Milan, Italy, 10 January 2021. [Google Scholar]
  24. Peng, X.; Qin, H.; Hu, Z.; Cai, B.; Liang, J.; Ou, H. Gas plume detection in infrared image using mask R-CNN with attention mechanism. In Proceedings of the Volume 11342 AOPC 2019: AI in Optics and Photonics, Beijing, China, 18 December 2019. [Google Scholar] [CrossRef]
  25. Hewett, T.A.; Fay, J.A.; Hoult, D.P. Laboratory experiments of smokestack plumes in a stable atmosphere. Atmos. Environ. 1971, 5, 767–789. [Google Scholar] [CrossRef]
  26. Sharifullin, V.N.; Badriev, A.I. Aerodynamic Characteristics of the Cooling Tower under the Nonuniform Distribution of the Water and Air Flows. Therm. Eng. 2019, 66, 569–574. [Google Scholar] [CrossRef]
  27. Badriev, A.I.; Makarova, I.V. The water and air distribution law analysis in natural draft cooling towers. In Proceedings of the 2021 3rd International Youth Conference on Radio Electronics, Electrical and Power Engineering (REEPE), Moscow, Russia, 11–13 March 2021; pp. 1–5. [Google Scholar] [CrossRef]
  28. Abdulla, W. Mask R-CNN for Object Detection and Instance Segmentation on Keras and Tensor Flow. n.d. Matterport. Available online: https://github.com/matterport/Mask\_RCNN (accessed on 10 October 2020).
  29. Blender Foundation “Blender” Software Package, v. 2.80. 2018. Available online: http://www.blender.org (accessed on 12 June 2020).
  30. Slawson, P.R.; Csanady, G.T. On the mean path of buoyant, bent-over chimney plumes. J. Fluid Mech. 1967, 28, 311–322. [Google Scholar] [CrossRef]
  31. Weil, J.C. Plume rise. In Lectures on Air Pollution Modeling; Venkatram, A., Wyngaard, J.C., Eds.; American Meteorological Society: Boston, MA, USA, 1988; pp. 119–166. [Google Scholar] [CrossRef]
  32. Briggs, G.A. Plume rise predictions. In Lectures on Air Pollution and Environmental Impact Analyses; Haugen, D.A., Ed.; American Meteorological Society: Boston, MA, USA, 1975; pp. 59–111. [Google Scholar] [CrossRef]
  33. Hoult, D.P.; Fay, J.A.; Forney, L.J. A theory of plume rise compared with field observations. J. Air Pollut. Control Assoc. 1969, 19, 585–590. [Google Scholar] [CrossRef]
  34. Fay, J.A.; Escudier, M.; Hoult, D.P. A correlation of field observations of plume rise. J. Air Pollut. Control Assoc. 1970, 20, 391–397. [Google Scholar] [CrossRef]
Figure 1. Photo of (a) a twelve-cell MDCT and (b) an illustration of the components of a single cell.
Figure 1. Photo of (a) a twelve-cell MDCT and (b) an illustration of the components of a single cell.
Atmosphere 14 00754 g001
Figure 2. An illustration of an air parcel for which the force balance and air entrainment form the basis of the model.
Figure 2. An illustration of an air parcel for which the force balance and air entrainment form the basis of the model.
Atmosphere 14 00754 g002
Figure 3. Basic geometric parameters for a multi-celled MDCT.
Figure 3. Basic geometric parameters for a multi-celled MDCT.
Atmosphere 14 00754 g003
Figure 4. Average enthalpy transfer rates for the twelve cells that were computed with the 5 min exhaust air temperatures, regional meteorological data, and basic cooling tower data.
Figure 4. Average enthalpy transfer rates for the twelve cells that were computed with the 5 min exhaust air temperatures, regional meteorological data, and basic cooling tower data.
Atmosphere 14 00754 g004
Figure 5. Measurement methods at the cooling tower site with the bounding plume with red dots representing the ground cameras.
Figure 5. Measurement methods at the cooling tower site with the bounding plume with red dots representing the ground cameras.
Atmosphere 14 00754 g005
Figure 6. Plume identification in photos between (a) Mask R-CNN inference and (b) RGB image from one of the stationary cameras.
Figure 6. Plume identification in photos between (a) Mask R-CNN inference and (b) RGB image from one of the stationary cameras.
Atmosphere 14 00754 g006
Figure 7. Comparison of a three-dimensional plume volume created from digital images to 12 overlapping plumes simulated by the one-dimensional model. The number of on and off fans in (a,b) were 8 and 4, respectively, whereas 7 fans on and 5 fans off for (c,d).
Figure 7. Comparison of a three-dimensional plume volume created from digital images to 12 overlapping plumes simulated by the one-dimensional model. The number of on and off fans in (a,b) were 8 and 4, respectively, whereas 7 fans on and 5 fans off for (c,d).
Atmosphere 14 00754 g007
Figure 8. Comparison of power output between 289 simulated and observed data.
Figure 8. Comparison of power output between 289 simulated and observed data.
Atmosphere 14 00754 g008
Figure 9. Comparison to measured average power ratios correlated with DPTD.
Figure 9. Comparison to measured average power ratios correlated with DPTD.
Atmosphere 14 00754 g009
Figure 10. Example Mask R-CNN inference overlay on RGB image.
Figure 10. Example Mask R-CNN inference overlay on RGB image.
Atmosphere 14 00754 g010
Figure 11. Example measurements extending (a) outwards 26 pixels. and (b) inwards 83 pixels.
Figure 11. Example measurements extending (a) outwards 26 pixels. and (b) inwards 83 pixels.
Atmosphere 14 00754 g011
Figure 12. Relationship between percent increase in dilated volumes relative to original volumes.
Figure 12. Relationship between percent increase in dilated volumes relative to original volumes.
Atmosphere 14 00754 g012
Figure 13. Percent error in power estimate as a function of number of individual power estimates averaged.
Figure 13. Percent error in power estimate as a function of number of individual power estimates averaged.
Atmosphere 14 00754 g013
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sobecki, C.; Garrett, A.; d’Entremont, B.; Connal, R.; Aleman, S. Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes. Atmosphere 2023, 14, 754. https://doi.org/10.3390/atmos14040754

AMA Style

Sobecki C, Garrett A, d’Entremont B, Connal R, Aleman S. Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes. Atmosphere. 2023; 14(4):754. https://doi.org/10.3390/atmos14040754

Chicago/Turabian Style

Sobecki, Christopher, Alfred Garrett, Brian d’Entremont, Ryan Connal, and Sebastian Aleman. 2023. "Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes" Atmosphere 14, no. 4: 754. https://doi.org/10.3390/atmos14040754

APA Style

Sobecki, C., Garrett, A., d’Entremont, B., Connal, R., & Aleman, S. (2023). Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes. Atmosphere, 14(4), 754. https://doi.org/10.3390/atmos14040754

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