1. Introduction
Airborne Particulate Matter (PM) has long been identified as an environmental and health hazard [
1]. Road transport is a major contributor to ambient PM, especially in urban areas [
2]. Traffic-related PM originates from both exhaust and non-exhaust sources, with brake-wear being the main contributor to non-exhaust PM
10 [
3,
4]. The technological improvements triggered by the continuously tighter regulation worldwide resulted in a significant reduction in exhaust PM emissions [
2,
5]. The relative share of brake-wear PM emissions is projected to steadily increase owing to the exhaust emission reductions brought by advanced after-treatment technologies and the forecasted electrification of the automotive sector [
6].
Brake-wear PM is a complex mixture of particles of vastly different chemistry, morphology and size. The physicochemical properties strongly depend on the brake system configuration, the brake-friction material and the operating conditions [
7]. Two distinct formation processes have been identified. Mechanical abrasion/adhesion leads to particles having a large size, with approximately 40% of their mass being respirable [
3,
8]. At sufficiently high frictional heat generation, organic constituents like phenolic resins used as binders [
8] can evaporate and subsequently nucleate, forming particles as small as 1.3 nm [
9].
Characterization of brake-wear particle emissions under real-world operation imposes several challenges, including, amongst others, the proper collection of supermicron particles [
10] and interference on the convective cooling of the brakes [
11]. Brake-dynamometers allow for an alternative approach offering tighter control of the braking and cooling conditions, allowing for the development of a standardized test procedure.
In this direction, the United Nations Particle Measurement Programme Informal Working Group (UNECE PMP IWG) has been actively working on the development of a harmonized brake-dynamometer methodology. The first step was the development of a novel brake cycle that would be representative of real-world operating conditions [
12]. The cycle (hereinafter WLTP-Brake) was based on the analysis of the database from the Worldwide harmonized Light Duty vehicles Test Procedure, containing in-use driving data with a total mileage of 743,694 km. The following step was the definition of an appropriate sampling setup to address the inherent difficulties in handling brake-wear particles. On the one hand, the large fraction of super-micron particles imposes challenges in their transportation [
10]. Additionally, their formation is highly localized, leading to non-homogeneous concentration profiles, while the large variety of available brake systems leads to complex air flows around the brakes. Accordingly, it was decided that enclosing the brake system in a dilution tunnel would be necessary for representative PM and Particle Number (PN) sampling. Furthermore, isokinetic sampling requirements necessitate the operation of the tunnel at a constant flow rate for the entire cycle.
The establishment of a thermal load on the brake system representative of real-world operation has been identified as a critical aspect of the methodology. On the one hand, this can affect the wear-rates [
13,
14] and accordingly PM emission rates. On the other hand, the developed temperatures can affect the release rate of organic compounds and accordingly the nucleation rates and thus the number concentration of volatile species [
12]. Due to the nonlinear nature of the nucleation process, a mismatch of the developed temperatures can lead to orders of magnitude increase in number concentrations unless the particles are thermally treated [
15,
16]. Hence, embedded thermocouple measurements of disc temperatures have been proposed as a means of ensuring thermal loads representative of real-world operation. Ideally, target temperature profiles should be established from real-world driving of each individual vehicle under the WLTP-Brake speed profile and then matched on the dyno by means of appropriate selection of the tunnel flow. However, this would be impractical, especially considering the varying ambient conditions. Accordingly, the use of universally applicable threshold values was proposed, subject to amendment as more vehicle data become available.
Given the emphasis of the proposed methodology on embedded thermocouple measurements, dedicated tests were conducted to assess the potential and limitations of this approach. Previous work [
17] suggested that the repeatability of thermocouple measurements can be affected by the development of hot rings, the radial location of which within the contact region of the pads changes with time [
18]. Such hot bands arising as a result of uneven distribution of the surface pressure (due to surface roughness and patch dynamics) were recorded with thermal cameras. In order to gain insight into the spatial and temporal variation of disc temperatures, dedicated tests were conducted on the same brake system. An array of four thermocouples was installed at four radial locations equally spaced along the contact area of the pads, to allow for real-time information on the spatial distribution of the disc temperatures. The experiments were conducted on a bigger dilution tunnel to assess whether the earlier findings are reproduced at enhanced convection cooling offered by a ~10-fold higher flowrates (900 to 1200 m
3/h).
A recent study on the same brake system revealed that the burnishing procedure can have a strong effect on the particle number emissions [
19]. Elevated concentrations were observed when the bedding was performed using the WLTP-Brake cycle, proposed by the PMP IWG. In order to gain insight into the nature of these particles, dedicated tests were performed with an arsenal of instrumentation. Measurements included PM
2.5, PM
10 and number concentrations of thermally treated samples following the regulated exhaust procedure, as well as number concentration and real-time electrical mobility size distribution of untreated samples. In addition to the WLTP-Brake, a short version of the Los Angeles City Traffic (3h-LACT) cycle developed within the LowBraSys project [
20] was tested. The results provide evidence on the formation of thermally stable nanosized particles over the WLTP-Brake. Volatile particles were observed over the 3h-LACT cycle, the number concentration of which was found to be strongly affected by the tunnel flow.
4. Discussion
The proposed harmonized procedure for brake-wear measurements on a brake dyno relies on the use of an embedded thermocouple to assess whether the developed temperatures are representative of real-world operation. In order to establish representative thresholds on the measured disc temperatures, it is imperative to assess the sensitivity of such measurements on the spatial and temporal evolution of the temperature. An array of four thermocouples was installed at four different radial locations across the contact area of the pads, allowing for such an assessment over two different test cycles. With respect to the temporal resolution, the measurements at all radial locations showed that the temperatures are sharply rising during braking and then drop down at a much slower rate once the brakes are released. The significantly different time constants in the two heat processes are a direct consequence of the large differences between friction power and convection heat rates. Even under the more extreme cooling conditions, the heat transfer coefficient cannot exceed 100 W/(m
2K) [
25] which corresponds to 1.8 × 10
4 W/m
2 at the highest disc temperatures anticipated over the WLTP-Brake. In comparison, the friction power is in the order of 10
6 W/m
2 [
18], i.e., two orders of magnitude higher. Consequently, adjustments of the tunnel flow mainly affect the temperatures over cruising sections, leading to inconsistent effects on the peak disc temperatures (
Figure 4, [
17]). Thermal radiation is generally considered to be insignificant under normal operating conditions [
28,
29,
30].
The measurements also revealed a spatial distribution of temperatures during the braking phases, with the maximum occurring at a radial location within the contact area. On the one hand, the work done by friction force increases with radial location owing to the increase in the linear velocity and the contact area. On the other hand, the contact area and thus mass and heat capacity increases with radius. The actual geometry of the disc (i.e., slots and or grooves to enhance heat transfer and clear water films) also contributes to the spatial temperature distribution. At the same time, the application of a brake pressure leads to a deformation of the pads, leading to non-uniform pressure, with the effect increasing with the applied pressure. The latter could explain the systematic difference observed in the radial location of the peak disc temperature over the WLTP-Brake (15 bar maximum brake pressure) and the 3h-LACT (20 bar maximum brake pressure). Changes in the local torque transmission due to the expansion of the brake caliper would require higher brake pressure (> 30 bar) [
31]. It is, however, also possible that this difference reflects temporal heat transfer effects, since braking over the 3h-LACT is more frequent and occurs at elevated start temperatures.
Superimposed in this macroscopic heat transfer process are microscopic phenomena resulting from the development of hot rings, the radial location of which changes within the contact region of the pads over time. Such hot bands are the manifestation of an uneven distribution of the surface pressure on individual contact areas caused by surface roughness and patch dynamics [
18]. The formation and time evolution of such hot rings results in test-to-test variations in individual thermocouple measurements.
When averaged over the entire duration of the cycle, the spatial nonuniformities were marginal, with all thermocouples yielding similar results. The increase in the tunnel flow also resulted in a consistent reduction in the cycle-average temperatures at all radial locations, although the effect was generally small. This behavior is anticipated, considering that the actual duration of braking is only a small fraction of the whole cycle (10% over WLTP-Brake and 20% over 3h-LACT). The cycle-average temperatures would thus correspond to an average of the temperatures during braking and during convective cooling, weighted by their relative duration in the cycle. Spatial nonuniformities are smoothened out during the cooling of the discs, leading to similar results at all radial locations. Moreover, the larger weighting of the convective cooling sections also implies that the mean temperatures better represent the effect of adjusting the tunnel flow. Similarly, the averaging over the duration of the cycle flattens out temporal variations introduced by the formation of hot-rings, leading to more repeatable measurements.
An increase in the tunnel flow from 900 to 1200 m
3/h led to a 3.3 °C (from 63.6 °C to 60.3 °C) and a 12.1 °C (114.8 °C to 102.7 °C) reduction in the cycle-average temperature over the WLTP-Brake and the 3h-LACT cycle, respectively. Some dedicated temperature tests suggested a similar change in cycle-average temperature over WLTP-Brake when reducing the tunnel flow from 900 to 400 m
3/h. Accordingly, the effect of tunnel flow is expected to be relatively small compared to reported differences between different brake systems with large differences in their heat capacity [
17]. Establishing universally applicable thresholds would therefore require careful consideration of the different commercial systems and perhaps differences in the designs between regions.
PM
2.5 emission levels were little affected by the test procedure, averaging 2.2 mg/km and 2.6 mg/km over the WLTP-Brake and the 3h-LACT, respectively. The corresponding PM
10 emission levels were 2.5 and 3.3 times higher, respectively. The same brake system was tested using the same PM sampling train on a smaller tunnel operating at 170 to 270 m
3/h [
17] over the last section of the WLTP-Brake. Emission levels following ten repetitions of trip 10 were found to be comparable, with PM
2.5 averaging at ~2.6 mg/km and PM
10 at ~7.8 mg/km. Similarly, the 3h-LACT yielded the same levels of PM
2.5 but distinctly higher PM
10 emissions, with the PM
10/PM
2.5 ratio increasing from 3 to 3.6. Tests with other brake systems suggested that PM emission varies with the size of the brakes, with PM
2.5 over the last section of WLTP ranging between 1.5 mg/km and 4.1 mg/km, for two systems having an effective disc radius of 113 mm [
32] and 144 mm [
17], respectively. The PM
10 to PM
2.5 ratios for both brake systems (~2.7) were similar to what was observed in the present study.
The consistency in the recovered PM
10 to PM
2.5 ratios for both tunnels implies that the extracted mass-weighted size distributions were similar, given that the same PM sampling system (and thus transmission efficiency) was employed in the campaigns. The adjustment of the tunnel flows had a modest effect on both tunnels, at least at the flows employed. A simple non-dimensional analysis on the relative importance of gravitational and inertial losses can verify this. The settling velocity for a 10 μm unit-density particle is 0.003 m/s at standard conditions (1 atm and 21 °C). The corresponding average convection velocities ranged between 9.5 m/s (170 m
3/h) and 15 m/s (270 m
3/h) for the small tunnel (duct diameter of 80 mm) and between 11 m/s (900 m
3/h) and 15 m/s (1200 m
3/h) for the large tunnel (duct diameter of 175 mm). Therefore, the particle settling distance per meter of horizontal duct is in the order of 0.2 to 0.3 mm, and the associated gravitational losses are calculated to be less than 1% [
33]. Accordingly, gravitational losses are not expected to be critical in brake-wear tunnels. All other relevant loss mechanisms (aspiration, flow path changes) are linked to particle inertia. Their relative importance directly depends on the magnitude of the square root of the Stokes number (Stk
0.5). The latter equals 0.025×(U
0/D)
0.5, for a 10 μm particle at ambient conditions, where U
0 is the gas velocity in m/s and D is the characteristic size (i.e., duct diameter) in meters. The Stk
0.5 number ranges between 0.27 (170 m
3/h) and 0.34 (270 m
3/h) for the small tunnel and 0.2 (900 m
3/h) to 0.23 (1200 m
3/h) for the larger one used in the present study. As a reference, the inertial losses of 10 μm particles in a 90° bend for the above duct sizes and flows would be 15% to 23% for 80 mm ducts and 8% to 11% for 175 mm ducts. In fact, a single 90° bend is the major source of losses in the small tunnel [
32], while the larger tunnel employed here incorporates two 90° bends resulting eventually in similar losses. Overall, the effect of flow on particle losses is expected to be small at the conditions tested as confirmed by the gravimetric PM measurements. It should be emphasized that the presence of such bends is desirable as they enhance mixing to ensure homogenous concentration at the extraction plane [
27,
34].
Proper PM measurements require careful consideration of the entire flow path down to the filter. The diameters of the tubing employed to transport samples extracted from the tunnel are much smaller than the tunnel diameter (typically in the 4 to 12 mm range). While the use of exchangeable nozzles offers some flexibility, isokinetic extraction implies similar velocities at the tubing and the tunnel. On the one hand, this means that settling distances remain in the order of tenths of a millimeter per meter of horizontal tubing, but now are comparable to the actual tubing diameter. Similarly, the Stk
0.5 number will increase significantly (3-fold increase for a tubing of 1/10th of the tunnel duct diameter for the same velocities). Thus, both gravitational and inertial losses are more critical on the PM sampling train and need to be carefully considered. Both loss mechanisms can be avoided by avoiding bends and using straight vertical tubing [
27]. When this is not feasible due to the layout of the tunnel, like in the present study, the sample velocities should be maintained as low as possible. The 5 lpm flow used here corresponds to an Stk
0.5 number of 0.3 for a 10 μm particle and a probe diameter of 8 mm. The losses of 10 μm particles in the PM sampling train configuration employed were calculated to be 8% due to gravitational settling in the ~0.4 m horizontal section and ~10% due to inertial impaction on the sampling probe bend.
To our knowledge, this is the first study to reveal the release of thermally stable (at the 350 °C of the catalyst) nanosized particles under non-extreme operating conditions explicitly targeted by the novel WLTP-Brake cycle. These particles were distinctly different from volatile nanoparticles reported in previous studies, which were found to be efficiently removed by the use of either thermodenuders [
16] or catalytic strippers [
15] at 300 °C, while, the disc temperature as measured with the array of four thermocouples was lower than 200 °C, the temperature on the friction interface can be significantly higher owing to the much smaller mass, thermal conductivity and thermal diffusivity of the brake pads [
13] but also due to the thermal resistance caused by wear particles [
30,
35]. Tribochemical reactions occurring at elevated temperatures and pressures are difficult to describe owing to the complex composition of brake linings [
35]. These may include degradation of phenolic resin that can be catalyzed by metal particles [
36] but also oxidation of metallic components like iron, copper, brass and tin [
37]. The composition of the thermally stable nanoparticles observed in the present study is not known. While the possibility of being organics of low volatility cannot be excluded, it is possible that they have a solid core. Monitoring the release of such thermally stable nanoparticles peaking in the 10 nm size deserves special attention.
The formation of these thermally stable nanoparticles was found to strongly depend on the bedding-in procedure employed to condition the brake system [
19].
Figure 9 compares the PN emissions over the last section of the WLTP-Brake cycle, where these particles were released, following different bedding-in procedures for the same brake system. The same brake system was also burnished following 10 repetitions of the last and more aggressive section of the WLTP-Brake, at the same facility in the context of a parallel study conducted at the same period [
19]. Subsequent repetitions of the WLTP-Brake cycle showed no excessive PN emissions over the 295th brake event, where the thermally stable nanoparticles were observed. Unfortunately, only the UCPC was available in these tests, so no information is available on their particle size. Similarly, no release of nanosized particles was observed in an earlier study [
17], where the same brake system was tested using fresh pads over consecutive repetitions of the last section of the WLTP-Brake. The different emission behavior was consistently observed through several repetitions of the different bedding-in procedures [
19]. These findings suggest that the initial transformation of the friction surfaces by the thermal, mechanical and chemical processes occurring during the actual burnishing procedure can strongly affect the particle emission performance of the brake systems. More research is necessary to better understand the nature of emitted nanoparticles and their formation processes.
Volatile nanoparticles like those reported in brake dyno [
16,
27] and chassis dyno [
15] measurements were also observed in the present study when tested under the same test cycle (3h-LACT). In contrast to thermally treated PN, the total PN concentration was found to depend on the operating tunnel flow. More specifically, the excess number of particles measured with the UCPC dropped from 185% (±40%) to 85% (±10%) when increasing the tunnel flow from 900 to 1200 m
3/h, when APC concentrations differed by less than 10%. The most probable formation mechanism of such volatile nanoparticles is the nucleation of evaporated organic compounds like phenolic resins typically employed as binders [
8]. The nucleation rates and therefore number concentration of volatile nanoparticles will ultimately depend on the developed saturation ratios in the tunnel. Given the complex flow patterns in the vicinity of the brakes, it is expected that the vapor concentrations, the gas temperatures and therefore saturation ratios will exhibit a large spatial distribution. Despite these complexities, an increase in the tunnel flow would lead to overall higher dilutions and thus reduction in vapor concentrations, suppressing nucleation rates. Furthermore, increased tunnel flows can result in reduced temperatures (
Figure 3) that could reduce the amount of released vapors, thus further reducing nucleation rates.
The established methodology constitutes a significant step towards the harmonization of brake-wear emission measurements. The comparability between the PM10, PM2.5 and their ratios for tests contacted at dilution tunnels operating at very different flows are very promising results. A potential introduction of the methodology in a regulation would require the careful consideration of some additional performance criteria, with the repeatability and reproducibility of measurements being perhaps the most important. One topic of active research pertains to the establishment of a proper degreening procedure and the definition of appropriate criteria to ensure that the emissions are stabilized. The potential formation of nanosized particles like those observed with the brake-pad tested will result in unstable PN emissions. Owing to their very small size, such nanoparticles will make an insignificant contribution to PM. PM measurements will therefore not detect unstable emissions during the bedding procedure.
While we have not observed the release of volatile particles over the WLTP-Brake procedure, it is not clear whether this is universally applicable. On the one hand, commercially available brake linings show large differences in their chemical formulations, which may also change in the future [
38]. Therefore, the volatility of organic compounds can differ between formulations. On the other hand, embedded disc thermocouple measurements do not necessarily capture the temperatures relevant for the release of volatile material, considering the higher temperatures developed on the contact area. Furthermore, homogeneous nucleation is not a binary process, and the rate of nuclei formation depends strongly on the ambient conditions and vapor concentration/chemistry. Adjustments in the tunnel flows and differences between tunnel designs can therefore affect their relative release rates. The systematic reduction observed in the relative amount of volatile PN over the 3h-LACT tests when increasing the tunnel flow is a manifestation of the nonlinear nucleation process. Significantly higher number concentrations of volatile particles were reported in tests of another brake system over the 3h-LACT both on a chassis dyno [
15] and on a brake-dyno under more appropriate ventilation conditions [
27]. If the total particle population is to be considered, it is important to carefully assess the volatile formation potential over the WLTP-Brake for different brake systems. Thermal treatment of the samples at 300 °C was found to efficiently remove volatile species [
16], indicating that the regulated methodology for automotive exhaust would be appropriate. At the same time, it would allow for the identification of the release of thermally stable nanoparticles that might deserve special attention. More research is needed, though, to better understand the nature of both the volatile and thermally stable nanoparticles in order to assess the potential need for amendments (e.g., the need for a catalytic stripper).
The mode of the number-weighted brake-wear size distribution is generally reported to lie in the submicron size range [
32,
38]. Therefore, an upper limit of 2.5 μm should be sufficient. Elevated number concentrations are generally associated with the release of nanosized particles. Given the relatively large size of brake-wear particles, it is advisable to implement a pre-classifier to reduce the risk of clogging. This risk is considerably higher compared to exhaust measurements owing to the significantly larger size of brake-wear particles. As such, the regulatory requirement for full-flow CPC designs that would allow direct monitoring of the flow entering the detector optics becomes more relevant. The flow split inside the UCPC employed in the present study was misadjusted due to contamination, despite the use of a 1.5 μm impactor at its inlet. The lowest cut-off size of full-flow CPCs is in the range of 6 to 7 nm and can be extended at smaller sizes by means of particle magnifiers [
9]. However, calibration at sizes below 10 nm becomes very challenging, while the lowest detection size will be substantially confined if a dilution stage will be required, owing to unavoidable diffusional losses.
The observed release of nanosized particles over the WLTP-Brake led to tunnel PN concentrations up to 1,000,000 cm
−3 as measured with the regulated exhaust 10 nm method at a tunnel flow of 1200 m
3/h. Such concentrations are far exceeding the measurement range of commercial full flow CPCs in single count mode (typically 50,000 to 100,000 cm
−3). Thus, a minimum dilution of at least 10:1 is required to accurately quantify the PN concentrations. An introduction of a dilution stage necessitates the characterization of size-dependent particle losses as well as an agreed approach for their compensation. The approach followed in exhaust PN regulation involving the quantification of Particle Concentration Reduction Factors (PCRF) is well established and can be transferred. Similarly, the size ranges of 15, 30, 50 and 100 nm are suitable for capturing diffusional losses in diluters, also considering the limitations in the production of monodisperse sizes. However, the relevance of an average compensation at 30, 50 and 100 nm needs to be assessed, as it was established on the basis of typical diesel exhaust number-weighted distributions peaking at approximately 60 to 70 nm [
26]. The number-weighted distribution of brake-wear particles is known to exhibit a mode above 100 nm, with a second mode in the sub-20 nm range occurring at specific events that can dominate in number.