{'fraction_non_alphanumeric': 0.025427132573299965, 'fraction_numerical': 0.017023511601322334, 'mean_word_length': 4.891518019657808, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Evolution of Soot Particle Number, Mass and Size Distribution along the Exhaust Line of a Heavy-Duty Engine Fueled with Compressed Natural Gas

An experimental study has been conducted to provide a characterization of the transformations that particle size distributions and the number density of soot particles can encounter along the exhaust line of a modern EURO VI compliant heavy-duty engine, fueled with compressed natural gas. Being aware of the particles history in the exhausts can be of utmost importance to understand soot formation and oxidation dynamics, so that, new strategies for further reducing these emissions can be formulated and present and future regulations met. To this purpose, particle samples were collected from several points along the exhaust pipe, namely upstream and downstream of each device the exhaust gases interact with. The engine was turbocharged and equipped with a two-stage after-treatment system. The measurements were carried out in steady conditions while the engine operated in stoichiometric conditions. Particle emissions were measured using a fast-response particle size spectrometer (DMS500) so that size information was analyzed in the range between 5 and 1000 nm. Particle mass information was derived from size distribution data using a correlation available in the literature. The reported results provide more insight on the particle emission process related to natural gas engines and, in particular, point out the effects that the turbine and the after-treatment devices produce on soot particles. Furthermore, the reported observations suggest that soot particles might not derive only from the fuel, namely, external sources, such as lubricant oil, might have a relevant role in soot formation.


Introduction
Recently, many governments around the world have announced their willingness to ban internal combustion engines (ICEs) in a short amount of time. This might suggest that the era of the ICEs has come to an end. Nevertheless, several concerns have been raised about the feasibility to power the entire transportation sector with electricity, especially regarding the land and marine commercial vehicles [1][2][3][4]. Although touted as the only possible solution to zero-emissions transport, electric-based mobility does not come without its own set of side effects [5][6][7].
The strong limitations related to batteries' energy density and cost make full electrification actually convenient only for small passenger cars [2,8,9]. Moreover, the absence of tailpipe emissions brings the large part of benefits in urban areas (if electricity is generated sufficiently far from the urban center). For commercial transport, full electrification is not practical or even desirable at the moment, because it would reduce the capacity of the truck to carry goods and increase the costs [2]. Due to the high The aim of the present work is to highlight the transformations that the PSD and, consequently, the total PN and PM can undergo along the exhaust line of a turbocharged heavy-duty CNG engine. Namely, the objective is to ascertain the effects produced by the interactions between the exhaust gases and each of the devices mounted on the engine exhaust line. A fast-response differential mobility spectrometer was used to sample directly from raw, upstream and downstream of the turbine and of each of the two catalysis stages.

Experimental Apparatus
The experimental campaign was carried out on a heavy-duty Euro VI CNG engine. It was selected for the present analysis because representative of the latest generation of the heavy-duty CNG engine technology. It was a production engine and it was tested with no modifications. The main specifications of the engine are reported in Table 1. The engine had 4 cylinders and 4 valves per cylinder. It was equipped with an electronically driven high-pressure injection system, based on a multipoint sequential phased injection strategy. It was turbocharged and equipped with a two-stage aftertreatment system composed of a close-coupled pre-catalyst and an underfloor TWC. A fast-response Differential Mobility Spectrometer Cambustion DMS500 was used throughout the experimental campaign for sampling directly from the exhaust line. Its operation is based on the electrical mobility and it was used to measure PN concentrations and PSD functions, from 5 to 1000 nm, at 10 Hz. The DMS500 is equipped with a built-in two-stage dilution system specifically for direct engine exhaust measurements. The first dilution stage at the sampling point avoids condensation and agglomeration, while the second diluter is employed for sampling from a wide range of concentrations. A 1.5 m heated line was used for sampling the engine exhausts in order to avoid condensation of combustion water.

Experimental Procedure
The schematic diagram of the experimental setup is depicted in Figure 1, in which the sampling points located along the exhaust line are highlighted by means of green circles. Temperature was measured at each of these points in order to ascertain possible influences that changes of the thermodynamic conditions along the exhaust line can have on particle emissions.
The considered engine was the subject of a recent work [44], in which its entire operating range was characterized in terms of PN emissions (the engine speed was varied from 800 to 3500 rpm and the engine load from 20% to 100% of the full load). Although the results showed that particle emission levels increased with engine load and speed, the highest value did not correspond to the maximum engine speed and load. Namely, the highest PN levels were recorded for engine loads comprise between 70% and 80% and engine speeds between 65% and 80% of the corresponding maximum values. Therefore, an engine speed of 2500 rpm and 80% of the full load were considered as operating conditions, because they are deemed as the most representative in terms of soot emissions for the purpose of the present work, which extends the analysis of particle emissions to the entire exhaust line.  To ensure proper operation and reliable response of the particulate analyzers, the engine was first warmed up at 1500 rpm. The selected operating conditions were then reached, and data were recorded after a stabilization period. A sampling time of about 1.5 min was chosen for ensuring good statistics and accounting for the presence of occasional spikes in the emitted particles per minute that the present engine showed [44] and that characterize soot emissions of SI engines in general [64].

Data Reduction
The PSD functions were provided by the sampling instrument in terms of normalized concentration dN/dlogD p /cm 3 , in which dN is the number of particles in the instrument channel width D p , and dlogD p , is the difference in the log of the channel width. Throughout this work, PN data are expressed in terms of particles/kWh, for generality, which means that the dN/dlogD p /cm 3 data were post-processed by considering the values measured for the exhaust mass flow rate, the engine torque and speed, and considering that the instrument had a resolution of 16 channels per decade.
PM information was derived from the particle size and number data, and its value is reported in the present work with respect to the mass of fuel, namely in g/kg fuel . The approach proposed by Maricq and Xu [65] was considered in order to take into account the fractal structure of soot particle clusters deriving from the agglomeration of spherical carbonaceous primary particles. Given the size distribution, N D p , the corresponding total PM can be expressed in terms of the soot particle mobility volume, V D p , and effective density, ρ e D p , (i.e., the soot mass per unit mobility volume) as follows: where V D p = π 6 D 3 p ; the dependence of the particle effective density on the mobility diameter, D p , indicates that the particles exhibit a fractal-like geometry [65][66][67]. Namely, because of their dendritic structure, the mass of the soot particle agglomerates does not scale as the cube of the diameter, as in the case of convex geometries, but rather as the diameter raised to the power of the fractal dimension, d f , according to the following power law: where D 0 is the primary particle diameter and ρ 0 is the primary particle density [65][66][67].
Since the number of size bins in the differential mobility instruments employed in the present work is finite, the integral in Equation (1) can be substituted by a summation extended over the total number of the size bins, N bins , (i.e., the number of electrometers in the instrument), as: in which N i represents the number of particles in the i-th bin corresponding to a mobility diameter D p i , and the expressions relative to ρ e D p and V D p were used. Maricq and Xu [65] observed that ρ e increases steeply with the decrease of D p until about 50 nm, where it breaks from the power law dependence and then levels out to a value approximatively equal to 1 g/cm 3 , and they found that such a trend was independent of the soot source (namely, of considered fuel or vehicle type) [65,67]. It follows that the hypothesis of spherical particles having unit density, as in the "standard spherical calibration" of the DMS500, is valid below a certain particle size, while a fractal-like geometry needs to be considered above that. Therefore, as proposed in [66], total PM was evaluated as the sum of the mass of the nucleation mode (NM) particles, PM NUC , and the accumulation mode (AM) particles, PM ACC , namely: In particular, it was considered that NM comprises spherical particles (d Nuc f = 3) with D p < 50 nm and having a primary density ρ Nuc 0 of 1 g/cm 3 , and that AM comprises particles characterized by a fractal dimension d Acc f = 2.3, a primary diameter D 0 = 20 nm and a primary density a value of ρ Acc 0 of 2 g/cm 3 , as suggested in various literature works [65][66][67][68][69]. Under these hypothesis, dedicated expressions for PM NUC and PM ACC can be derived: where here N Nuc bins and N Acc bins represent the total numbers of size bins corresponding to mobility diameters smaller than 50 nm and larger than 50 nm, respectively. It follows that Equation (2) can be rewritten as: Equation (4) represents the form used in the present work for deriving PM values from PSD data provided by the sampling instrument. Figure 2 highlights the reduced slope of the particle mass curve with respect to the particle diameter after 50 nm, deriving from Equation (4) (solid red line) in comparison to the "standard spherical calibration" of the DMS500 (dashed black line). The values reported in Figure 2 are nondimensionalized with respect to the mass of a spherical particle with unit diameter and density. Figure 2 highlights the reduced slope of the particle mass curve with respect to the particle

Results
Along the exhaust line, variations of the temperature and interactions with both the turbine and the after-treatment system produce significative transformations of soot particles, which induce PN and PM variations, as shown in Figures 3 and 4, respectively. Figure 3 shows that total PN increases in the first part of the exhaust line, reaching a peak downstream of the pre-catalyst and then it is significantly reduced by the passage through the TWC. At its peak, PN is roughly two orders of magnitude higher than the value recorded upstream of the turbine. After the TWC, its value returns comparable to the initial value. Figure 4 shows that a trend similar to that of total PN is also recorded for PM. Namely, its highest value is recorded downstream of the pre-catalyst, while a significant reduction is observed after the main catalysis stage. It is worth noting that the TWC abatement action is stronger in terms of PM in comparison to PN. After the TWC, PM levels are the lowest recorded along the exhaust line, being about one order of magnitude lower than those initially observed upstream of the turbine (which are the second lowest ones). This can be explained by considering that larger particles can be more easily oxidized by the passage through the TWC, and particles with larger sizes contribute more significantly to soot mass in comparison to PN. Thus, a variation in the concentration of these particles produces more relevant variations of PM than those of PN values.

Results 192
Along the exhaust line, variations of the temperature and interactions with both the turbine and 193 the after-treatment system produce significative transformations of soot particles, which induce PN 194 and PM variations, as shown in Figures 3 and 4, respectively. Figure 3 shows that total PN increases 195 in the first part of the exhaust line, reaching a peak downstream of the pre-catalyst and then it is 196 significantly reduced by the passage through the TWC. At its peak, PN is roughly two orders of 197 magnitude higher than the value recorded upstream of the turbine. After the TWC, its value returns 198 comparable to the initial value. Figure 4 shows that a trend similar to that of total PN is also recorded 199 for PM. Namely, its highest value is recorded downstream of the pre-catalyst, while a significant       In order to highlight that particles having different sizes show different behaviors in traversing the exhaust line, total PN was decomposed in NM and AM fractions, and the result is depicted in Figure 5. This decomposition is well known in the literature, and its results are particularly useful for the analysis, because from the particle size it is possible to infer information about the nature of the particles. Namely, NM comprises particles with size smaller than 50 nm, which are generally considered formed from condensation of volatile material, primarily due to gas-to-particle conversion, whereas AM particles are solid and predominantly carbon-based products of combustion, having size comprised between 50 and 1000 nm [40,70]. Figure 5 shows that the fraction of NM particles increases along the exhaust line. If upstream of the turbine, about 60% of the particles belongs to the AM; downstream of the TWC, the NM share of PN emissions exceeds 90%. Moreover, it is possible to observe that the total PN peak (located downstream of the pre-catalyst) is characterized by the highest value of the nucleation mode fraction. Downstream of the TWC, although a slight increase of the AM fraction is observed, NM still remains markedly predominant.
Energies 2020, 13, x FOR PEER REVIEW 7 of 17 particles. Namely, NM comprises particles with size smaller than 50 nm, which are generally 216 considered formed from condensation of volatile material, primarily due to gas-to-particle 217 conversion, whereas AM particles are solid and predominantly carbon-based products of 218 combustion, having size comprised between 50 and 1000 nm [40,70].  For a more detailed analysis, the PSD functions recorded along the exhaust line are provided in Figure 6, together with the corresponding 95% confidence intervals. The comparison between the PSDs recorded upstream and downstream of the turbine, depicted in Figure 6a, highlights that the finest particles are the most affected by the passage through the turbine. Namely, NM particles increase by about one order of magnitude, while AM particles are little affected. This points out that the increase in total PN (cf. Figure 3), as well as the variation of the relative composition in terms of mode fractions (cf. Figure 5) are exclusively due to an increase of the NM particles downstream of the turbine. There are evidences suggesting that nucleation of soluble organic fraction (SOF) is promoted through the turbine [71]. A large share of the SOF is usually related to lubricant oil [46,72,73]; namely, lubricant oil components can represent up to the 80% of the SOF [74].

241
It is noteworthy that lubricant oil can cause the behavior observed downstream of the turbine in 242 more than one way. Namely, it can penetrate into the combustion chamber following many routes 243 [75,76] and contributes to increas these lubricant-related SOF components, which can subsequently 244 coagulate under the turbine action, determining a gas-to-particle conversion. On the other hand,  It is noteworthy that lubricant oil can cause the behavior observed downstream of the turbine in more than one way. Namely, it can penetrate into the combustion chamber following many routes [75,76] and contributes to increase these lubricant-related SOF components, which can subsequently coagulate under the turbine action, determining a gas-to-particle conversion. On the other hand, lubricating oil can reach the exhaust stream, for instance escaping from the turbocharger bearings lubricating circuit, and determine a direct increase of the number of the finest particles. There are many studies in which the high emission levels of very fine particles in CNG engines are indirectly associated to lubricant oil [37,40,51,52,54,77,78]. In a recent work by Amirante et al. [38], the formation of particles produced solely from lubricant oil was demonstrated by means of experimental tests designed for that very purpose. The results ascertained the concentration number and size distribution of lubricant-oil-derived particles according to the different ways by which the oil escapes from its own circuit. In particular, it was found that lubricant oil contribution is always detectable in the lowest range (never exceeding 30 nm), independently of the way it was added, confirming what was supposed and indirectly observed in previous research works.
A possible additional phenomenon that can be generated by the turbine is the fragmentation of some particles due to the interaction with the turbine blades [79]. Fragmentation mainly affects larger particles, which can become detectable by the instrument [80,81], explaining the increase of the particles with size larger than 350 nm observed in Figure 6a. The increase of the number of particles having such large sizes is also the main cause of the PM increase recorded after the turbine and shown in Figure 4. This aspect is made clearly visible in Figure 7, in which the contributions of NM and AM particles to PM are reported separately and thus their relative importance is ascertained. associated to lubricant oil [37,40,51,52,54,77,78]. In a recent work by Amirante et al. [38], the formation 249 of particles produced solely from lubricant oil was demonstrated by means of experimental tests 250 designed for that very purpose. The results ascertained the concentration number and size 251 distribution of lubricant-oil-derived particles according to the different ways by which the oil escapes 252 from its own circuit. In particular, it was found that lubricant oil contribution is always detectable in 253 the lowest range (never exceeding 30 nm), independently of the way it was added, confirming what 254 was supposed and indirectly observed in previous research works.

255
A possible additional phenomenon that can be generated by the turbine is the fragmentation of 256 some particles due to the interaction with the turbine blades [79]. Fragmentation mainly affects larger 257 particles, which can become detectable by the instrument [80,81], explaining the increase of the 258 particles with size larger than 350 nm observed in Figure 6a. The increase of the number of particles 259 having such large sizes is also the main cause of the PM increase recorded after the turbine and shown 260 in Figure 4. This aspect is made clearly visible in Figure 7, in which the contributions of NM and AM 261 particles to PM are reported separately and thus their relative importance is ascertained.

262
Once left the turbine, the exhaust gases encounter the first stage of the after-treatment system.

263
The comparison between the corresponding PSDs recorded downstream of the turbine and the pre-264 catalyst is provided in Figure 6b. It shows that NM is strongly influenced, as well as the effects on   Once left the turbine, the exhaust gases encounter the first stage of the after-treatment system. The comparison between the corresponding PSDs recorded downstream of the turbine and the pre-catalyst is provided in Figure 6b. It shows that NM is strongly influenced, as well as the effects on AM particles larger than 100 nm are not negligible. The number of particles with size ranging from 150 to 1000 nm reaches the highest level along the exhaust line. PN is increased of averagely one order of magnitude in the size range comprised between 250 and 500 nm after the passage through the pre-catalyst. The increase of both NM and AM particles explains the further increases in total PN of about one order of magnitude after the pre-catalyst shown in Figure 3. On the whole, the predominance of NM particles is further enhanced, considering that they now account for more than 95% of the total emitted particles, compared with the 85% downstream of the turbine, as shown in Figure 5. The non-negligible increase of the number of relatively large AM particles explains the further increase observed in PM as well, as shown in Figure 4.
The simultaneous concurrence of many complex phenomena is at the basis of the simultaneous and significative increase of both PN and PM due to the pre-catalyst. Fragmentation, agglomeration, oxidation and vaporization are all relevant processes occurring in the pre-catalyst [59,82]. For instance, larger particles suspended into the exhaust gases can undergo fragmentation due the interaction with the internal pre-catalyst structure, determining an increase of the number of smaller particles [60]. At the same time, the passage through the narrow sections of this device forces the particles to collide and bond together, forming larger ones. Moreover, the high temperatures can activate surface oxidation reactions that can significantly contribute to reduce the size of the particles. Simultaneously, the increased temperature can also promote the evaporation of the film that coats some particles, leading to a reduction of their sizes. Thus, a transfer from particle-to-gas phase might be possible in this catalysis stage.
The thermal conditions reigning in this device play the main role in controlling particle transformation processes. Ensuring high working temperatures in the pre-catalyst is essential for a proper operation of both catalysis stages, considering that the pre-catalyst is mainly installed for reducing the light-off time of the TWC and consequently increasing the overall efficiency of the after-treatment system [83]. Figure 8 shows that after the drastic decrease due to the energy extracted by the turbine, the exhaust gases temperature is again significantly increased in the pre-catalyst. On the contrary, a drastic fall of the temperature is recorded after the TWC, indicating that the two after-treatment devices work in different conditions and behave differently with respect to soot emissions. As a consequence, the evolution of soot PN and PM through the after-treatment system appears to be in close correlation with the temperature trend reported in Figure 8.
Energies 2020, 13, x FOR PEER REVIEW 10 of 17 some particles, leading to a reduction of their sizes. Thus, a transfer from particle-to-gas phase might 286 be possible in this catalysis stage.

287
The thermal conditions reigning in this device play the main role in controlling particle 288 transformation processes. Ensuring high working temperatures in the pre-catalyst is essential for a 289 proper operation of both catalysis stages, considering that the pre-catalyst is mainly installed for 290 reducing the light-off time of the TWC and consequently increasing the overall efficiency of the after-291 treatment system [83]. Figure 8 shows that after the drastic decrease due to the energy extracted by 292 the turbine, the exhaust gases temperature is again significantly increased in the pre-catalyst. On the 293 contrary, a drastic fall of the temperature is recorded after the TWC, indicating that the two after-294 treatment devices work in different conditions and behave differently with respect to soot emissions.

295
As a consequence, the evolution of soot PN and PM through the after-treatment system appears to 296 be in close correlation with the temperature trend reported in Figure 8. contributing to a significant decrease of PM [59,60,84]. Moreover, the passage through the TWC 304 enhances the tendency for soot particles to stick together, increasing agglomeration rates and forming 305 larger particles reducing NM particles as well [57,60]. Coagulation of liquid components can also be 306 promoted. Figure 6c indicates that the reduction concerns both NM and the AM, explaining the 307 simultaneous decrease recorded for total PN (cf. Figure 3).

308
It is worth noting that the PSD function recorded downstream of the TWC shows the largest during the regular engine operation [60]; thus, non-uniform particle release due to occasional 314 breakups of such deposits can explain the larger variability observed after the main catalysis stage.

315
The presence of particle spikes during steady-state measurements are reported in the literature in the 316 case of gasoline and CNG spark-ignition engines [64]. This aspect was highlighted in the previous 317 experimental campaign carried out on the same engine [44]. These emission spikes were observed at 318 every operating condition and it was found that they were almost completely composed of nucleation 319 mode particles; furthermore, their magnitude and frequency varied depending on the case. The The significant decrease of total PM downstream of the TWC shown in Figure 4 can be explained considering that processes having a mechanical nature (rather than a thermochemical nature) dominate in this device, due the lower temperatures. The catalyst surfaces in a TWC are extremely porous and can act like a particulate filter, trapping a significant fraction of large particles and contributing to a significant decrease of PM [59,60,84]. Moreover, the passage through the TWC enhances the tendency for soot particles to stick together, increasing agglomeration rates and forming larger particles reducing NM particles as well [57,60]. Coagulation of liquid components can also be promoted. Figure 6c indicates that the reduction concerns both NM and the AM, explaining the simultaneous decrease recorded for total PN (cf. Figure 3).
It is worth noting that the PSD function recorded downstream of the TWC shows the largest variability over time, as can be inferred from the larger width of the error bars obtained at this point of the exhaust line, in comparison to all of the other cases. Moreover, this PSD function presents a less regular trend with respect to the particle size. Several peaks can be observed, with the highest one located between 25 and 30 nm. It is known that soot particles are deposited on catalyst surfaces during the regular engine operation [60]; thus, non-uniform particle release due to occasional breakups of such deposits can explain the larger variability observed after the main catalysis stage. The presence of particle spikes during steady-state measurements are reported in the literature in the case of gasoline and CNG spark-ignition engines [64]. This aspect was highlighted in the previous experimental campaign carried out on the same engine [44]. These emission spikes were observed at every operating condition and it was found that they were almost completely composed of nucleation mode particles; furthermore, their magnitude and frequency varied depending on the case. The higher temperature in the pre-catalyst can justify the fact that this behavior is not observed for that device, which shows a much more regular and stable PSD function (cf. Figure 6c). In many research works carried out on CNG engines, ash deposits on the walls of the tailpipes are ascribed to lubricating oil [39]. A transfer from the gas to the particle phase, with nucleation of material mainly constituted by oil additives, might be enhanced after the TWC due to the lower temperatures [29,54,77]. Once formed, these deposits can be stochastically released and emitted by the engine, explaining the less regular PSD function recorded downstream of the TWC. Such a result underlines once more the primary role that lubricant oil can play in determining soot emissions in CNG engines.

Conclusions
In the present work a detailed experimental characterization of the evolution of total particle number (PN), particle mass (PM) and particle size distributions (PSDs) was provided along the exhaust line of a modern EURO VI compliant turbocharged heavy-duty engine, fueled with compressed natural gas (CNG). Namely, the influence on each of the devices mounted on the engine exhaust line was investigated. Particle samples were collected upstream and downstream of the turbine, the pre-catalyst, and the three-way catalytic converter (TWC). A fast-response differential mobility spectrometer, such as the DMS500, was employed for the analysis. The decomposition of total PN in nucleation mode (NM) and accumulation mode (AM) was also provided.
The highest value of the total PN was recorded downstream of the pre-catalyst, where it was about two orders of magnitude higher than the value recorded at the beginning of the exhaust line, namely, upstream of the turbine. A significative decrease of PN was recorded after the TWC. PM showed a similar trend, progressively increasing until the pre-catalyst, where its highest values were recorded, and then drastically reducing after the passage through the TWC. The fraction of NM particles increased along the exhaust line, from about 40% upstream of the turbine to more than 90% downstream of the TWC. The highest value of the NM fraction corresponded to the total PN peak.
It was observed that the turbine produced an increase of the total PN, and this was reasonably attributed to possible induced nucleation of soluble organic fraction components. It has been pointed out that lubricating oil significantly influences this trend in many ways. Moreover, the increase of PM detected downstream of the turbine was attributed to possible fragmentation of very large particles that became detectable by the instrument. The concurrence of various possible phenomena in the pre-catalyst justified the generalized increase recorded, independently of the particle size. Fragmentation, agglomeration, and oxidation compete with each other in the first catalyst stage. On the contrary, the TWC led to a strong reduction of both NM and AM particles, and this was explained considering the temperature reduction and its filtering action.
It has been highlighted that the temperature variations along the exhaust line can play a crucial role in determining the observed particle transformation trends. In fact, it was highlighted that the evolution of PN and PM through the exhaust line was closely correlated to the temperature variations.
Moreover, the occasional release of ash deposits, mainly due to lubricating oil, explained the larger irregularity and variability observed for the PSD function recorded downstream of the TWC. This aspect, together with the above-mentioned ones, highlights the significative role that lubricant oil play in determining soot emissions in CNG engines.
This work represents a first step toward the definition of more complete and accurate analysis procedures, which can help to better understand the formation mechanism of ultrafine soot particles. Future experimental campaigns will deal with analyses carried out in transient conditions as well as with specific tests designed to isolate the sources of ultrafine soot particles (lubricant oil first among them) along the exhaust line of heavy-duty CNG engines. A comparison with samples conditioned with a volatile particle remover system will make it possible to analyze the evolution of the particulate volatile organic fraction.
Author Contributions: All the authors (E.D., R.A., G.C., P.D.P., P.T.) performed the experiments, collected and analyzed the data, wrote, edited and revised the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflicts of interest.