Open Access Article
Lisa Nordhausen
a,
Marie Denise Gutsch
a,
Daniel Hodonj
b and
Marion Börnhorst
*a
aInstitute of Reaction Engineering and Catalysis, TU Dortmund University, Emil-Figge-Straße 66, 44227 Dortmund, Germany. E-mail: marion.boernhorst@tu-dortmund.de
bInstitute for Chemical Technology and Polymer Chemistry (ITCP), Karlsruhe Institute of Technology (KIT), Engesserstraße 20, 76131 Karlsruhe, Germany
First published on 4th May 2026
To investigate the spatial distribution of gas-phase species within structured catalysts under conditions representative of real exhaust gas environments, an invasive capillary-based sampling method was adapted for use in a hot gas test rig. A Fourier-transform infrared (FTIR) spectrometer was employed for analysis of gas mixtures characteristic for the selective catalytic reduction (SCR) of NOx. To meet the sample volume requirements of the FTIR, a coil-based sampling extension was introduced and its influence was systematically quantified. A linear correlation between gas composition measurements with and without the coil enabled the development of a correction factor, which was generalized using a residence time model, thereby eliminating the need for extensive recalibration under different operating conditions. In addition, the impact of the capillary on local conversion behavior was assessed using an analytical model, accounting for uncertainties in capillary positioning and reduced gas velocities. These findings facilitated reliable spatially resolved species measurements for NH3-SCR reactions. Extension of the methodology to a SCR configuration including the injection of urea water solution revealed significant challenges for a precise quantification of the gas composition due to transient phenomena such as film evaporation and deposit formation.
For this purpose, the SpaciMS technology was developed, which enables spatially resolved measurements of concentrations in individual channels of a honeycomb catalyst with the use of a capillary.3 While the capillary is sequentially moved along the channel axis, a mass spectrometer (MS) analyzes the sample composition, which requires comparatively small sample volumes for frequent analysis and therefore allows for online-monitoring of the intracatalytic concentrations.3 In parallel, a thermocouple is inserted and moved within the channel with identical protocol. This technology enables in-operando spatially resolved concentration and temperature measurements, bridging integral and differential analysis of reaction kinetics. However, as a physical probe (capillary) is inserted to one or more catalyst channels, the SpaciMS technology is invasive. Previous works have investigated and quantified this influence using analytical models, computational fluid dynamics (CFD) simulations as well as experiments.4–7 The flow disturbance of the sampling capillary within the catalyst channel is highly dependent on its radial position. When positioned at the center of the channel, the capillary induces the greatest perturbation to the flow field, whereas placement near the channel corner results in minimal flow field disruption. However, the exact radial position of the capillary cannot be accurately determined during the experiments, introducing a source of uncertainty.5
SpaciMS was used to investigate different exhaust gas aftertreatment systems such as diesel oxidation catalysts,8,9 three-way-catalysts10,11 and the selective catalytic reduction (SCR) of nitrogen oxides (NOx).12,13 In the context of SCR, challenges arise due to the superposition of the molar masses of H2O and NH3. This was described by Luo et al.14 and Hou et al.15 proposing to use a Fourier transform infrared spectrometer (FTIR) as analyzer instead of a MS. The SpaciFTIR technology enables the analysis of more complex gas mixtures, making it particularly valuable for studies related to Fischer–Tropsch synthesis16 and ammonia oxidation.17 Additionally, it provides important insights into SCR technologies in both conventional diesel engines and future hydrogen-powered engines.18,19 However, in contrast to the MS, an FTIR requires a significantly larger sample volume. While this cannot be achieved through the capillary, Lou et al.14,20 propose a dilution of the capillary sample with N2. The required degree of dilution depends on the ratio of the capillary volumetric flow rate to the inert gas flow rate. For low capillary flow rates as considered in this study, high degrees of dilution would be required, which reduce the sensitivity of the FTIR measurements, thereby decreasing accuracy and limiting the detection of trace-level species. The low capillary flow rate is attributed to the use of a thin capillary, which minimizes flow distortion but leads to a low volumetric flow rate due to the high pressure drop across the capillary.
The present study applies the SpaciFTIR technology to the example of the SCR of NOx. This technology is used for NOx emission control of, e.g., diesel engines, which are harmful to both the environment and health. Using NH3 as reducing agent, NOx are reduced following one of three reactions, depending on the NO/NOx ratio:
| 4 NO + 4 NH3 + O2 → 4 N2 + 6 H2O | (R1) |
| NO + 2 NH3 + NO2 → N2 + 3 H2O | (R2) |
| 8 NH3 + 6 NO2 → 7 N2 + 12 H2O | (R3) |
The targeted conversion of urea to ammonia takes place in two steps. The first step is the urea (CO(NH2)2) thermolysis, which produces NH3 and isocyanic acid (HNCO). This is followed by the hydrolysis of HNCO, which produces NH3 and CO2:27
| CO(NH2)2 → NH3 + HNCO | (R4) |
| HNCO + H2O → NH3 + CO2 | (R5) |
Although the urea-SCR has been representing the state-of-the-art technology for NOx abatement in mobile applications for two decades, the requirements of modern combustion engines and the ever stricter regulations and standards (Euro 7, US Tier 3)32,33 demand for even higher efficiencies. Concepts aimed at reducing engine-out NOx emissions, such as exhaust gas recirculation (EGR) or the intentional lowering of combustion temperatures, effectively decrease NOx formation during the combustion process. However, these methods result in reduced exhaust gas temperatures, which adversely affect the urea thermolysis and HNCO hydrolysis. This, in turn, delays ammonia generation and slows down the reaction kinetics within the SCR catalyst. Close-coupled SCR systems, positioned near the engine, benefit from faster heat-up times and higher exhaust temperatures during cold-start and low-load operation. However, they introduce new challenges, particularly due to the shortened mixing length between UWS injection and catalyst inlet. This can lead to spatially non-uniform ammonia distribution in the exhaust stream, thereby reducing NOx conversion.24,34 For the optimization and targeted design of future exhaust gas aftertreatment systems, enabling improved efficiency and robustness under dynamic operating conditions, detailed insight into the simultaneous reaction kinetics of urea decomposition and NOx reduction over the SCR catalyst is essential. As the majority of studies on SCR use gaseous NH3 as the reductant rather than dosing urea, there is still limited information available on urea decomposition in combination with SCR reactions. SpaciFTIR enables spatially resolved analysis of the reaction progress and enables a deeper understanding of the thermolysis of urea, the subsequent catalytic hydrolysis of HNCO and the SCR reactions.
This study introduces a novel coil-based sampling configuration coupled with SpaciFTIR technology. A methodological framework is developed and generalized on the example of NH3-SCR, explicitly accounting for systematic measurement errors. The approach is subsequently applied to the urea-SCR system, providing first insights into the simultaneous reactions of urea-derived byproducts and NOx over SCR catalysts through spatially resolved kinetic data.
The synthetic exhaust gas flows into the injection section at a flow rate of approximately 1000 l min−1 and a gas temperature of approximately 280 °C. UWS is injected with a mass flow rate of 0.877 g min−1 and a pressure of 5 bar. This is accomplished through pulsed dosing at a frequency of 1 Hz with a duty cycle of 1%. Immediately downstream the urea injection zone, a 150 mm long commercial vanadia-based SCR catalyst (Umicore) is installed. It features a 600 cpsi honeycomb structure and a cross-sectional area of 30 × 60 mm. For reference measurements, a 120 mm long uncoated (catalytically inactive) honeycomb monolith with identical cell density and cross-sectional area is employed to investigate non-catalytic effects. To enable spatially resolved gas phase mole fraction measurements within the catalyst, a deactivated fused silica capillary is inserted into a single channel of the honeycomb monolith. Unless otherwise specified, the capillary is inserted to a channel in the top third of the monolith. The capillary is mounted on a motorized linear translation stage (OptoSigma), allowing precise axial positioning for spatially resolved sampling. The extracted gas is analyzed using an FTIR spectrometer (MultiGas 2030, MKS Instruments) to determine the local mole fractions of NH3, NO, and CO2 along the axial length of the catalyst. These species are selected to facilitate the investigation of urea thermolysis, HNCO hydrolysis, and the SCR reactions. As HNCO cannot be detected with spatial resolution due to a low signal-to-noise-ratio, CO2 is used as a marker for hydrolysis progress, given its stoichiometric formation during the HNCO hydrolysis. Due to the low volumetric flow rate through the capillary of approximately 0.002 l min−1, direct FTIR analysis in a sample chamber of 0.2 l would yield poor temporal resolution and increase the risk of cross-contamination. To mitigate these limitations, a coil-based configuration, analogous to that described by Tomin et al.,35 is employed. The sample is pumped from the capillary through a 6-way valve using a vacuum membrane pump (KNF). The 6-way valve connects the pump conveying the sampled gas from the capillary, the exhaust exit, the purge gas (N2), the coil, and the FTIR analytical system. In addition, it allows the coil to be operated with flow in both directions. Depending on the valve circuit, the sample flows either into the sample coil (V = 10 mL) (“load”) or into the exhaust (“sample”). The two valve circuits are shown in Fig. 2. After filling the coil with the sample, it is purged with nitrogen and the gas sample is directed to the FTIR. In addition, the mole fraction upstream and downstream the catalyst is measured directly with the FTIR, without using the coil setup. These measurements, referred to as global measurements, represent the flow duct-averaged mole fractions upstream and downstream the catalyst and allow for the additional quantification of species that cannot be spatially resolved due to their low mole fractions, such as HNCO and NO2. The upstream global mole fraction is sampled from the cover plate region, which may introduce inaccuracies due to potential liquid film formation and deposit accumulation upstream the catalyst. In contrast, the downstream global sampling is conducted in the outlet pipe, where a more homogeneous gas-phase mixture can be assumed due to enhanced mixing conditions.
Three distinct experimental configurations are employed in this study:
1. Coil characterization tests using NH3, NO, H2, and CH4 with an uncoated monolith,
2. Steady-state spatially resolved measurements of the NH3-SCR over the catalyst and the uncoated monolith,
3. Spatially resolved measurements with UWS dosing over the catalyst, conducted with and without NO dosing.
The following sections provide a detailed description of the three experimental configurations.
Prior to the spatially resolved measurements, the sampling capillary is inserted into the catalyst from the downstream side. During the experiment, it is incrementally translated in flow direction, enabling axial profiling of the gas composition within the catalyst channel. At each axial position, three replicate measurements are performed to ensure reproducibility.
Each measurement cycle of the spatially resolved measurements consists of three steps:
1. Gas sampling for a duration of 6 minutes to fill the sampling coil (valve circuit: “load”),
2. Sample transfer to the FTIR system for a duration of 1 minute (valve circuit: “sample”),
3. Rinsing of the sampling coil with nitrogen for 3 minutes to prevent cross-contamination between measurements (valve circuit: “sample”).
As a result, each measurement point requires approximately 30 minutes, including sampling, transfer, and purging cycles. Following each axial displacement of the capillary, a waiting period of 10 minutes was implemented to ensure that operation conditions were re-established prior to data acquisition. During this waiting period, the global gas composition was additionally measured to monitor the integral progress of urea thermolysis and HNCO hydrolysis. To achieve reproducible conditions during the measurement, the system is allowed to stabilize for 30 minutes after gas dosing begins and for 60 to 90 minutes following the initiation of UWS injection before the first measurement is taken.
![]() | (1) |
![]() | (2) |
N2 and Vcoil represent the volumetric flow rate of the purge gas and the internal volume of the coil, respectively. The parameters tpeak,start and tpeak,end define the temporal bounds of the peak, yi is the measured mole fraction and Nmeas the number of data points used in the integration of the peak signal. The corresponding integral used to quantify the species amount is discretized and evaluated via a rectangular (midpoint) numerical approximation, as illustrated in eqn (2). An empirical correction factor Acor is introduced to account for potential deviations due to calibration uncertainties or signal drift. The correction area is computed as the average mole fraction over the final 50 data points preceding the peak signal.
The FTIR analysis results reflect the mole fractions present in the coil and cannot be directly interpreted as gas composition present in the catalyst channel. Because the analysis is influenced by residence time distribution effects, the measured mole fractions must be corrected to accurately represent the actual mole fraction at the sampling location. The quantification of the influence of the coil setup and a corresponding correction method are described in the results section.
Furthermore, the capillary intrusion affects the measurement results. Hlavatý et al.6 demonstrated that the flow distortion of the capillary can be compensated through appropriate adjustment of the sampling flow rate. However, under the gas velocities and capillary dimensions (inner diameter 150 μm, outer diameter 375 μm) employed in the present study, such flow balancing is technically not feasible. To ensure accurate spatially resolved measurements the authors recommend to sample approximately 50% of the local channel flow rate through the capillary. Achieving this sampling fraction under the applied conditions would require a capillary flow rate of approximately 0.3 l min−1, resulting in an unacceptably high pressure drop. Potential alternatives, such as employing a monolith with lower cell density (e.g., 200 cpsi) or reducing the total volumetric flow rate, can increase the relative sampling fraction. However, these modifications would significantly alter the flow and thermal boundary conditions, thus compromising the representativeness of the exhaust gas conditions that the test rig is intended to replicate. Consequently, the approach presented by Hettel et al.5 is adopted. The proposed model enables the quantitative assessment of capillary-induced residence time distortions as a function of the axial insertion depth, by incorporating both geometric parameters (e.g., channel and capillary diameters) and operating conditions (e.g., gas velocity).
![]() | ||
| Fig. 4 Comparison of the ratio of the coil mole fraction and global mole fraction using a thick (ID 150 μm) and a thin (ID 100 μm) capillary over different sampling times. | ||
For a detailed quantification of the deviations of the measured coil mole fractions from the actual global value, the respective mole fractions for the two main components of SCR, NH3 and NO, are measured and compared for a concentration range relevant to SCR as shown in Fig. 5. A linear correlation is observed for both components, indicating that the relationship between the coil-derived and global mole fractions is independent from the absolute concentration levels. To quantify this relationship, a linear regression was applied to the data points, constrained to pass through the origin. This constraint is physically justified, as in the absence of dosing, the corresponding mole fraction must be zero. The resulting regression line thus provides a consistent basis for correcting coil measurements across varying concentrations.
![]() | ||
Fig. 5 Measured global mole fractions over coil mole fractions of a) NO and b) NH3. g = 1000 l min−1, Tg,in = 25 °C. | ||
Based on the linear regression, a correction factor (C) is introduced which corresponds to the slope of the respective function and is defined as the ratio of the coil mole fraction and the global mole fraction (yi,global) of component i:
![]() | (3) |
![]() | (4) |
Due to the high experimental effort and the need to repeat this procedure for each operating condition, an analytical model for the correction factor is applied. This model is based on the residence time sum function, which describes the cumulative probability that a fluid element introduced at the reactor inlet at time t = 0 will exit the system within a time interval [0, t]. Given that the Reynolds number in the coil remains below 1 (Re < 1), laminar tube flow conditions are assumed in the coil. Accordingly, a laminar flow residence time distribution (RTD) model is applied. Furthermore, axial dispersion is considered negligible, as the Bodenstein number satisfies Bo ≫ 100 for all investigated gas species, indicating minimal backmixing and validating the laminar tubular flow approximation. This results in the following relationship for the hydrodynamic residence time:36
![]() | (5) |
![]() | (6) |
![]() | (7) |
). By substituting the analytical velocity profile for laminar flow in a cylindrical tube, the expression for τ can be reformulated (eqn (6)). This relationship is subsequently incorporated into the cumulative residence time distribution function, F(t), yielding eqn (7).
Since the residence time sum function must remain positive for all t, eqn (7) is not valid for short times. Therefore, for
, an alternative expression is applied to ensure mathematical and physical consistency. This model allows for the analytical prediction of the residence time sum curve as a function of both elapsed time and sampling flow rate. Based on that, a theoretical coil mole fraction value is determined by establishing a time-discretized mole balance for the species (eqn (8)). Incorporating the ideal gas law, along with the definitions of mole fraction and concentration and under the assumption of constant inlet and outlet flow rates, leads to the expression given by eqn (9),
| ni,coil(t + Δt) = ni, coil(t) + ṅi,in − ṅi,out | (8) |
![]() | (9) |
cap the volume flow through the capillary. To validate the model, the theoretically calculated correction factors are compared with experimentally determined values in Fig. 6 for different gas species. The black line represents the theoretically calculated correction factor, while the cross markers denote the corresponding experimentally determined values. The species investigated include NO and NH3, both of which are directly relevant to SCR reaction mechanisms; H2, which is of growing interest in the context of future hydrogen-fueled engines; and CH4, an important reactant in various combustion and power-to-X applications. The comparison reveals that, particularly at short sampling durations, the model provides a sufficient representation of the experimental data for all species. For higher sampling times, the model exhibits a slight underestimation of the correction factor for NO, H2, and CH4, while the overall agreement remains satisfactory. A notable deviation is observed in the case of ammonia, for which the model consistently underpredicts the correction factor across all sampling intervals. This discrepancy is attributed to adsorption phenomena on the inner wall of the sampling coil,32,33 as evidenced by the delayed desorption observed during the nitrogen purge phase, where NH3 requires significantly longer flushing times to be fully displaced.
![]() | ||
| Fig. 6 Comparison of the modeled (black line) and the experimentally (crosses) determined correction factor over the sampling time for different species: a) H2, b) NO, c) CH4 and d) NH3. | ||
The extent of this flow disturbance strongly depends on the axial and radial position of the capillary. For downstream insertion, as applied here, the influence increases as the capillary tip approaches the catalyst inlet. This can be explained by the total pressure drop consisting of contributions from the empty channel section and the capillary-containing section. Since the former is significantly smaller, extending the capillary further upstream increases the overall pressure drop. Due to the increased pressure drop, the volumetric flow rate through the probed channel decreases, leading to a reduction in flow velocity. As the capillary is moved further downstream, two opposing effects occur: the reduction in flow velocity becomes less pronounced, while the length, over which the reduced velocity increases the residence time, becomes larger.
In addition, most conversion occurs in the front section of the channel, where reactant concentrations are highest, further amplifying the impact of upstream placement. The effect of radial position is commonly evaluated using two limiting cases, as the exact radial position during operation is unknown: maximum flow perturbation for a centered capillary and minimum disturbance for placement near the corner region.5,6
As previously discussed, the influence of the sampling capillary is quantitatively evaluated in this study using the analytical model developed by Hettel et al.5 The two limiting cases regarding the radial position are evaluated using the methodology of Gossler et al.37 for determining the friction factor. For both limiting cases, the mole fractions in a reference channel, unaffected by capillary-induced flow disturbances, are calculated for each axial position, representing the minimum and maximum measurement deviation. These estimates are presented in Fig. 7 alongside the capillary channel mole fraction, corrected using the developed methodology, on the example of ammonia SCR. The black dashed lines and empty green triangles represent the estimated mole fraction within a reference channel, where the upper tip corresponds to capillary placement at the radial channel center (maximum flow distortion), and the lower tip to placement at the corner (minimum flow distortion). The reference channel mole fractions of NO are presented here as a representative example; however, the observed deviations are applicable to all detected species.5
An approximately equimolar decline in NO and NH3 mole fractions is observed, consistent with the expected stoichiometry of the standard SCR reaction (eqn (R1)). As the axial shift along the catalyst increases, the concentration profiles of NH3 and NO converge progressively. This behavior is partially attributed to the presence of small quantities of NO2 (approximately 1% of the NO mole fraction), formed via equilibrium between NO and NO2. NO2 cannot be measured in spatially resolved manner due to its low concentrations.
Comparison with global measurements—conducted without the capillary setup—reveals that the mole fraction downstream the outlet of the capillary-probed channel equals the average across all channels. This demonstrates that the effect of the coil can be accurately accounted for using the correction factor. At the inlet, deviations—most pronounced for NO—can be attributed to turbulence effects upstream the honeycomb structure as well as slight variations in capillary positioning. The small error bars, derived from three-fold measurements at each point, demonstrate a high degree of reproducibility.
The influence of the capillary increases steadily over the first 60 mm, reaches a maximum at this position, and then decreases sharply. Although the flow distortion is greatest at a position of 10 mm as discussed previously, the impact on the mole fractions continues to increase up to 60 mm. This can be attributed to the fact that a longer upstream section with reduced velocity due to the capillary leads to a greater increase in residence time. It should be noted that the predicted mole fractions under conditions of maximum flow distortion are bounded by the inlet mole fractions to ensure physical consistency. However, the calculated values exceed the inlet mole fractions by up to 170 ppm. At 60 mm the initially high reaction rate, which decreases along the channel, coincides with the progressive residence time shift induced by the capillary. Beyond this point, the influence diminishes, as the overall reaction rate declines with decreasing reactant mole fractions and the mole fraction gradients become nearly uniform toward the channel outlet. Based on the observed data, it is inferred that the capillary is located close to the corner of the channel. This is evident from the observation that the values associated with the capillary effect align more closely with the global data when the capillary is positioned in the corner of the channel. For reasons of clarity, the influence of the capillary is not depicted in the following, as its precise radial position within the channel during the experiment is unknown, thereby impeding an exact quantification of its impact on the measured concentrations.
Model predictions show a sharp decrease in HNCO mole fraction in the inlet section approaching zero at approximately 120 mm. The CO2 mole fraction increases correspondingly, which can be attributed to HNCO hydrolysis. The predicted NO mole fraction decreases along the entire catalyst length due to the standard SCR reaction (eqn (R1)). The NH3 mole fraction increases slightly near the inlet and then decreases. This results from the superposition of NH3 production via HNCO hydrolysis and NH3 consumption by the SCR reaction.
The measurement uncertainties presented in Fig. 8 are significantly larger than those observed in Fig. 7. This is primarily attributed to transient effects introduced by pulsed UWS injection. The presence of a liquid phase, along with film and deposit formation upstream the catalyst, results in a spatially and temporally non-uniform reactant distribution across the flow duct, which will be addressed in more detail below.
First, the spatially resolved data are discussed, followed by a comparison with the global data. The NH3 mole fraction increases in the first section of the probed channel as a result of urea thermolysis and HNCO hydrolysis, and decreases toward the channel outlet due to the consumption of NH3 in the standard SCR reaction with NO. Similarly, the CO2 mole fraction rises within the first 45 mm of the channel, as expected from HNCO hydrolysis, and then stabilizes. This observation aligns with the expectation that HNCO hydrolysis predominantly occurs in the first half of the catalyst. However, the observed CO2 mole fraction at the catalyst inlet (0 mm) is higher than anticipated based on stoichiometric considerations. Approximately 420 ppm CO2 originate from the CO2 content in the feed air, implying that roughly 400 ppm CO2 is formed upstream of the catalyst. If this CO2 formation results from HNCO hydrolysis, approximately 800 ppm NH3 (400 ppm from urea thermolysis and 400 ppm from HNCO hydrolysis) would be expected upstream. However, measurements show only about 400 ppm NH3 at the catalyst inlet. As the mole fractions at the channel inlet are used as model input, this discrepancy is not reflected in the model predictions.
Furthermore, an offset between mole fractions immediately upstream (−4 mm) of and at the catalyst inlet (0 mm) is observed, which is likely attributable to backmixing and turbulence effects induced by the cross-sectional contraction at the catalyst entrance.
Contrary to expectations and the model predictions, the NH3 and NO mole fractions remain approximately constant throughout the second half of the catalyst. Despite the presence of both reactants, further SCR conversion appears to be negligible, suggesting a cessation or significant reduction of the reaction, which does not align with prior results. This can be attributed to the measurement protocol. It requires approximately 30 minutes per spatially resolved data point—comprising three repetitions of 6 minutes of gas sampling and 3 minutes of nitrogen purging—resulting in limited temporal resolution. Consequently, transient fluctuations, such as the entrainment of liquid droplets into the sampling capillary, may lead to local concentration spikes that are not adequately captured by the current approach. While good agreement between spatially resolved measurements and model predictions is observed in the upstream region, deviations increase toward the downstream section. This suggests that, with longer injection durations, the influence of film and deposit formation becomes more significant and cannot be adequately resolved. Additionally, simultaneous acquisition of global and spatially resolved data is not feasible as it would require the use of a second FTIR. As a result, temporal variations in the inlet gas composition cannot be tracked continuously, making it impossible to precisely correlate the inlet conditions with individual axial measurement positions.
In contrast to the spatially resolved measurements, lower conversions of NO are observed in the global measurements, which can be attributed to the increased residence time caused by flow distortion and associated stagnation upstream the capillary in the probed channel. Upstream the catalyst, approximately 450 ppm of CO2 is detected globally, which can be attributed to ambient CO2 present in the intake air. Measurement uncertainties for CO2 are higher than for NH3 and NO, primarily due to fluctuations in laboratory air composition, which affect the baseline CO2 level. Moreover, almost no NH3 is observed upstream the catalyst, suggesting that neither thermolysis nor HNCO hydrolysis occurs significantly upstream the catalyst. However, this observation may also be influenced by non-uniform NH3 distribution within the exhaust duct. Due to the pressurized injection of UWS from the top and notable liquid film formation and deposit accumulation on the bottom wall, the NH3 distribution is expected to be inhomogeneous.38 The global sampling position at the top wall of the duct upstream the catalyst is assumed to diminish representative sampling. In contrast, downstream global measurements are taken from the well-mixed exhaust gas flow, providing a more representative composition. The global measurements reveal only a marginal increase in CO2 mole fraction—approximately 30 ppm—accompanied by a significant NO slip of about 790 ppm. The measured global NH3 mole fraction downstream the catalyst is approximately 10 ppm. The observed CO2 increase implies that around 60 ppm of NH3 was generated through thermolysis and subsequent HNCO hydrolysis, which is consistent with the simultaneous reduction of NO by approximately 55 ppm (global) and the measured NH3 slip of 10 ppm. Experimental data and observations suggest that a substantial portion of the UWS forms solid deposits rather than undergoing the intended thermal decomposition and hydrolysis reactions to produce NH3. To enable a more detailed examination, the progression of the global data throughout the duration of the experiment is evaluated.
Fig. 9 presents the global mole fractions of NH3 (blue) and HNCO (red) during an experiment conducted under identical operating conditions but without NO dosing. The prominent peaks correspond to the global measurements, while the smaller, superimposed peaks in the NH3 signal reflect spatially resolved measurements. The global mole fractions exhibit significant temporal fluctuations. Initially, a continuous increase in both NH3 and HNCO mole fraction is observed, followed by a gradual decline after a transitional point. This behavior can be attributed to the evaporation of the majority of the liquid film in the early phase of the experiment, promoting urea thermolysis. As the experiment progresses and the temperature of the bottom plate decreases due to pronounced film formation and the formation of solid decomposition products becomes more dominant, less NH3 and HNCO are generated.
Closer inspection of the global data—exemplified by the measurement at 104 minutes in Fig. 9c—reveals that nearly identical mole fractions of NH3 and HNCO are present upstream the catalyst. The slight deviation is attributed to non-uniform concentration distribution across the duct and the fixed sampling position. Downstream the catalyst, however, only a few ppm of HNCO is detected, indicating efficient hydrolysis. Ideally, the sum of the upstream NH3 and HNCO mole fractions should correspond to the total amount detected downstream. The observed discrepancy reinforces the conclusion that incomplete mixing occurs upstream the catalyst, potentially affecting reaction efficiency and measurement accuracy.
To examine the non-uniform distribution of reactive species across the catalyst cross-section in detail, two experiments conducted under identical operating conditions within different measurement channels are compared in Fig. 10. In both cases, UWS was injected without the addition of NO to study urea thermolysis and HNCO hydrolysis. In the experiment labeled “top”, the sampling capillary was positioned in a channel located in the upper third of the catalyst cross-section, whereas in the “bottom” experiment, it was placed in a channel in the lower third.
Significantly reduced NH3 and CO2 mole fractions were observed in the “top” experimental configuration. This is attributed to gravitational spray settling and film formation at the bottom of the duct, acting as a reservoir for decomposition products, as well as the 30° downward injection angle of the UWS, which increases the likelihood of droplet impingement and localized concentration peaks in the lower region of the catalyst. In the upper sampling area, the NH3 mole fraction decreases over time, indicating a gradual reduction in the supply of HNCO and NH3 upstream due to deposit formation.
In contrast, the “bottom” experiment reveals significantly higher NH3 and CO2 mole fractions. The initial increase in NH3 mole fraction is consistent with the onset of HNCO hydrolysis, while the CO2 mole fraction, which also initially increases as a product of hydrolysis, subsequently slightly decreases. This decline is likely attributable to factors such as non-uniform reactant mixing, deposit formation and varying CO2 inlet concentrations. This strongly suggests the occurrence of additional reactions within the liquid film, leading to a pronounced spatial dependence of the mole fraction profiles.
Furthermore, the data indicate that HNCO hydrolysis progresses significantly faster in the lower half of the catalyst. However, further investigation is required to determine whether these observations are the result of spatial variations or are influenced by temporal effects.
These findings underscore the pronounced impact of injection strategy and potential wall film formation on the spatial distribution of reactive species within the catalyst. The results demonstrate that a precise quantification of spatially resolved mole fractions with the proposed measurement technique is strongly impeded by the transient characteristics of UWS dosing and decomposition measurements. Static mixers potentially mitigate spatial concentration inhomogeneities, however, strong turbulence and flow distortion have to be accepted.
To further generalize this correction, a laminar tube flow-based residence time model was applied, enabling the calculation of the correction factor without the need for extensive calibration experiments under different operating conditions. This model was validated against experimental data for multiple gas species and showed good agreement, with the exception of ammonia. For NH3, systematically lower mole fractions—and consequently underestimated correction factors—were observed. This deviation is attributed to adsorption effects on the inner surfaces of the coil, corroborated by the delayed desorption of NH3 during nitrogen purging.
In addition, the influence of the capillary on local catalytic conversion was evaluated using an analytical model. As this model does not account for the gas volume removed through the capillary, it serves as the worst-case estimate for the measurement deviation induced by capillary intrusion. The primary sources of uncertainty in this context include the unknown radial position of the capillary within the channel cross-section and the high gas velocities, which amplify flow disturbances. To increase the accuracy of error estimation, a method for continuous position tracking during operation must be developed. Alternatively, the approach employed by Hlavatý et al.6 to compensate for capillary-induced effects could be adopted. However, implementing such a method poses significant technical challenges under the experimental conditions applied in this study. Compared to the dilution method first proposed by Luo et al.,14 the present approach is particularly suitable for low capillary volumetric flow rates. However, due to its increased experimental complexity and potential limitations in temporal resolution, the method proposed by Luo et al. may be preferable at higher capillary flow rates. The applicability and limitations of both methods depend on the analytical technique employed, the target species, and the corresponding detection limits.
Despite the limitations of the developed methodology, it enabled reliable acquisition of spatially resolved data for the NH3-SCR reaction system. The approach was subsequently extended to the more complex urea-SCR reaction system. The obtained spatial mole fraction profiles partially deviated from expectations, largely due to the low temporal resolution of the measurement method and the transient nature of film and deposit formation during UWS dosing.
Based on these findings, adaptations to the experimental setup are recommended to improve measurement representativeness and reliability. The implementation of a static mixer in combination with optimized flow control can enhance the homogeneity of species distribution across the catalyst cross-section, thereby improving the spatial accuracy of global measurements.39,40 Alternatively, an optimized UWS dosing strategy could be developed to reduce temporal fluctuations and approach steady-state conditions. Such improvements would be particularly valuable for laboratory-scale studies aiming to investigate urea thermolysis and HNCO hydrolysis under controlled conditions, supporting the development and optimization of SCR catalysts. Additionally, to capture transient effects with greater fidelity, alternative analytical techniques with higher temporal resolution, which in this case particularly requires a smaller sample volume, ideally circumventing spectral interference between NH3 and H2O, should be explored.
| This journal is © The Royal Society of Chemistry 2026 |