1. Introduction
Flow measurement in stacks is of fundamental importance in the assessment of pollutant emissions, because when combined with the measurement of the concentrations, it provides the mass flux of any emitted pollutant. Inside stacks flow is not laminar in most cases and the length of straight and unobstructed pipe sections, available for measurements, is often not enough to allow the full development of a regular velocity profile. Consequently, velocity and flow measurements are often affected by large uncertainty. Small ducts (inner diameter ≤ 50 cm) are more influenced than larger ones by perturbations, even due to the presence of sensors. Ducts can include curves, scrubbers, swirlers, extraction ports, probes and other devices, which can perturb the movement of the exhaust gases. Another issue is related to the presence of pollutants in the gas stream that can produce fouling or corrosion of the measuring devices. For example, particulate matter, often present in high concentration in many plant exhausts, such as biomass boilers, can obstruct the orifices for measuring the differential pressure or clog the moving parts of an anemometer. A recent review paper [
1] addresses the origin of random and systematic errors for in-stack velocity measurements using Pitot tubes, how these problems are treated in the pertaining international standards and what effect they may have on the uncertainty of the measurements of pollutant emissions. From this survey, it emerges that, in the presence of cyclonic flows, the use of S-type Pitot tubes can result in errors of up to 12%, due to non-axial flows. A similar error can be produced by the misalignment of the Pitot tube during the measuring procedures. The errors associated with velocity measurements, using S-type Pitot tubes in cyclonic flows, have been investigated in detail also in [
2] by means of computational fluid dynamic (CFD) modelling. Different pipe configurations, producing different flow patterns have been studied and the associated errors have been analyzed. In particular, the authors pointed out that, in presence of an asymmetric velocity profile, typically with swirling flows, the maximum velocity in a cross-section of stack describes a spiral along the duct. In this case, the results of the measurements depend on the location of the measuring point, whose best position is unpredictable. Another source of error, associated with the use of S-type Pitot tubes and examined in [
2], is due to the inclined gas velocity vectors in presence of non-axial velocity components: the contribution of this error, estimated under the specific conditions simulated in [
2], is up to 2.5%. Measurements were carried out in a wind tunnel, in different experimental conditions and geometries in [
3], yielding similar results.
Several techniques are available for measuring in-stack velocity and flow rate: most of them are described in two international technical standards: EN ISO 16911-1:2013 [
4] and EN ISO 16911-2:2013, the former focuses on manual and the latter on automatic methods. The most commonly used approaches make use of Pitot tubes or vane anemometers to determine the flow velocity either on a single fixed point or on a grid of points on a measuring plane. The other methods described in [
4] include a tracer gas dilution technique and a calculation procedure, based on the energy consumption of the combustion plants. Both of them provide directly the gas flow rate and indirectly the average velocity in the stack. The tracer dilution technique follows a very basic approach, because it determines the flow rate using just two physical quantities: the mass flow rate of the injected tracer and the concentration of this tracer in the gas stream. The energy consumption approach relays on the possibility to measure the flow rate and composition of the fuel and the oxygen content in the fumes, continuously and very accurately. The Transit Time (TT) tracer gas method is another possible approach, briefly described in EN ISO 16911-1. This method allows determining the average velocity of the flowing gas within a portion of a duct having a constant cross-section. Many different solutions may be adopted to put this method into practice, which make use of a tracer gas which follows the gas flow and which can be measured with a sufficiently high time rate. For this technique to be adopted, some constraints should be fulfilled. It should be possible to inject a tracer inside the stack, as homogeneously (with respect to the transverse section of the duct) as possible. The tracer should not be naturally present in the atmosphere, should not interfere with the normal operation of the plant, should be neither poisonous or toxic, nor environmentally detrimental. It should not be a byproduct of the plant itself. For an accurate measurement the maximum response time of the detection technique must be ∼
ms, to be negligible compared to the rise time of the concentration of the tracer, which is in general ∼second. Short half-life radioactive tracers are often considered the most suitable, because they can be detected through the duct walls without the need for probes, ports or windows, but the use of these tracers is often restricted by national legislation. A detection technique that measures the average concentration of the tracer on a cross-section is to be preferred to a point measurement, because it reduces the effects of an improper mixing. Optical and laser techniques exist since long for the measurement of gas flows: Laser Two Focus [
5], Doppler effect [
6], Particle Imaging Velocimetry [
7]. Laser Two Focus investigates very small regions (0.25 mm × 0.25 mm × 0.25 mm) at time, so requiring some time, and moving optics, to obtain a raster image. Laser Doppler Velocimetry provides information along a direction at an angle with the stack axis. In order to keep this angle as small as possible, inclined beams are required. In ref. [
6] a setup with two measurement channels is shown. Particle Imaging Velocimetry requires the injection in the stack of particles with suitable dimensions, a pulsed, high power laser source, and a fast, high-resolution camera. Finally, a careful calibration procedure must be carried out. TT and Dilution techniques are much simpler, requiring low power laser sources and standard detectors, with straightforward analysis procedures. Spectroscopy is particularly suitable for the TT tracer gas method. In this case, the tracer should feature optical absorptions in wavelength regions where user-friendly laser sources are available, and these absorptions should be sufficiently strong to reduce the tracer concentration to very low levels (≤1‰).
In the framework of the EMPIR Project IMPRESS II, a laser-based device, presented in this article, has been developed, which allows us to perform both dilution and TT measurements of the in-stack velocity. The two methods have been tested and cross-validated in a stack simulator, having a small section, under strongly cyclonic flow conditions. The flow pattern inside the test rig has been simulated through CFD modeling, using Ansys® Academic Fluent, Release 15.0, ANSYS, Inc., Canonsburg, PA, USA in order to get a better understanding of the experimental conditions. In parallel, fixed point measures, using an S-type Pitot tube, have been carried out and compared under conditions which are particularly critical for using a Pitot tube, in order to evaluate the possible bias between this conventional method and the other two optical techniques.
2. Optical Detection
In order to adopt any optical technique, we had to choose the tracer molecule. Following the constrains described above, we selected acetylene (C
H
). It has a strong absorption band around 1520 nm, which is one of the regions of optical telecommunications, and both laser sources and optical devices are available off-the-shelf. It is not toxic, and there is no risk of explosion at concentrations below 1‰(Low Explosion Level for acetylene is 2.5% [
8]).
A possible drawback is that it can be produced during incomplete combustion of methane [
9]. However, most of our measurements were carried out when heating the air inside the stack by means of electrically driven resistors, so the risk of interferences was avoided. In a real case, other fuels than methane could be used and, even in the case of incomplete combustion of methane, the injection of an extra amount of acetylene can be easily detected. In order to verify which concentration should be used in the measurements, to have a sufficient Signal-to-Noise Ratio, we examined the behavior of the absorptions in the temperature range of 285–385 K.
Molecular absorptions are described by the Beer-Lambert law:
where
and
are the powers of the light beam after and before crossing the sample, respectively;
L is the length of the sample, and
where
S (cm/molecule) is the absorption strength,
(cm) the shape of the absorption (area normalized to 1), and
n (molecule/cm
) the density of the absorbing species. One of the consequences of this behavior is that the absorbance of a transition depends on temperature, as both linestrength
S and lineshape
are sensitive to temperature.
Figure 1 shows the comparison of the transmissions, in a narrow wavelength range, at 285 and 385 K, at P = 1 Bar, L = 0.6 m, of a simulated exhaust mixture including H
O 10%, CO
10%, O
4% (green line), and the same mixture with addiction of C
H
1‰ (blue) [
10].
It is evident that linestrengths decrease with increasing temperature. The best choice with respect to interference from other molecules is the line close to 1.521
m, even if it is not the most intense one. Its strength decreases by about 38% at 385 K (with respect to 285 K), nevertheless, the green line in
Figure 1 is very close to transmission 1 (i.e., unperturbed transmission) at both temperature conditions. Finally, a concentration around 1‰ is sufficient to produce absorptions of a few %, easily detectable by any optical technique. Actually, we easily worked at levels about one order of magnitude less than this. Due to all these features, and because suitable laser sources are commercially available at this wavelength, the transition at 1.521
m should be selected in real operation conditions. As a matter of fact, in many of the present measurements we used hot air, or room temperature air, rather than exhaust gases, so the concentrations of water and carbon dioxide were the atmospheric values. For this reason in our tests we used either the absorption at 1.521
m, or the more intense one at 1.520
m.
As for the detection technique to be used, our choice was between Direct Absorption (DA) and Wavelength Modulation Spectroscopy (WMS) [
11], depending on the kind of flow measurement. In DA the power of a light beam, transmitted across a sample, is converted by a detector into a voltage or current. When sweeping the laser wavelength, the absorption profile is retrieved. The analysis of this profile yields the absolute value of the concentration, provided that pressure, temperature, pathlength and molecular linestrength are known. In WMS, a modulation is applied to the laser wavelength at a period
much lower than the sweep time
(typically
). The detector signal is demodulated at a multiple of the modulation frequency. The most used choice in commercial devices is a demodulation at twice the modulation frequency. In this case, the signal is approximately the second derivative of the absorption profile, and the central peak is proportional to the density of absorbing molecules. WMS exhibits a better signal-to-noise with respect to DA [
12]. In a previous paper [
13] we demonstrated the calibration problems of WMS when the composition of the sample varies, in particular because of large percentages of carbon dioxide and water, problems which do not occur with DA [
14]. As a matter of fact, the two different optical flow measurements (TT and dilution) feature different calibration issues: TT is related to the time evolution of the concentration of the injected tracer and not to the absolute value of the mixing ratio. On the contrary, dilution requires an accurate measurement of the tracer mixing ratio. For these reasons we chose DA for dilution, and WMS for TT measurements.
5. Discussion and Perspectives
We have applied two spectroscopic detection techniques to the measurement of flow in narrow ducts, and to the calibration of standard sensors and methods. Our multipurpose device was deployed in a stack simulator, proving to fulfill all the requirements for the above task. It is already stated in EN ISO 16911-1 that the dilution method is a valuable reference method, which is intrinsically free from wall effects, or from any perturbation due to the geometry of the duct, or to any objects inserted in the duct, inducing turbulence. On the contrary, the more turbulent the flow, the more homogeneous the distribution of the tracer. TT method, in particular the one implemented in this work, based on optical techniques, compares very closely with dilution.
The conventional procedure adopted for calibrating Pitot tubes is carried out ex-situ in a standardized wind tunnel using a primary reference device. It generally provides either a single calibration factor or a set of calibration factors for the different velocity ranges, used to calculate the point velocity from the readings of the differential pressure at the Pitot ends. Calibrated Pitot tubes are then used both to measure in stack velocity directly and to periodically calibrate other automated measuring systems, installed on the stacks to measure the flow rate on line. This study shows that this approach can produce inaccurate results when cyclonic flows in small ducts are involved, in particular:
the conditions of 5 internal diameters of straight duct upstream and 3 downstream of the measuring section, without elements disturbing the flow, may not be a sufficient guarantee to provide a regular velocity profile at the measuring plain;
single-point measurements may be non-representative of the average velocity, a multipoint technique should be required even for ducts having an internal diameter smaller than 0.35 m;
manual alignment of a Pitot tube is not accurate and stable enough, large deviations of the measured values may be produced by very small variations in the placement of the probe;
techniques relying on the bulk properties of the flow, such as the concentration of a tracer injected in the gas stream, are largely insensitive to the flow pattern and the local disturbances, producing linear and accurate results.
As a consequence of these considerations, the pertaining international standards, namely UNI EN 15259:2008, EN ISO 16911-1:2013 and EN ISO 16911-2:2013, should include specific warnings for small ducts, suggesting the use of reliable and robust techniques for the in-site calibration of the automated measurement systems, such as dilution based and transit time methods. Moreover, apart from periodical checks, the necessity to repeat the calibration every time a modification occurs in the duct, upstream the sensor should be introduced.
Table 2 compares the requirements of the different techniques used in this work. Pitot tubes are undoubtedly the simplest technique for flow velocity measurements, as they require an insertion port only, and no consumables. TT requires four optical ports. Dilution can be implemented either across the stack, or after gas extraction. In the first case, two optical ports are necessary, in the second case an extraction port (downstream the Pitot, to avoid any interference) and a heated line must be used. In both the latter cases, a tracer is required, which means consumables and an injection point. Despite the higher complexity, we proved that in narrow ducts it is necessary to add ports, or injection/extraction points, suitable for the application of more complicated techniques, for intercomparison and calibration. As a final remark, the described optical techniques are not so much time consuming, as they require one workday for set-up, measurement and packing.