A novel microwave-based approach for monitoring the state of aftertreatment systems such as diesel particulate filters (DPFs), three-way catalytic converters (TWCs), and selective catalytic reduction (SCR) catalysts is proposed. The volume inside the metallic housing of the DPF, TWC, or SCR is considered as a wireless communication channel between two terminals of a communication system. It is shown that, depending on the transmission channel characteristics, the properties of the catalyst, such as the catalyst state, can be inferred. This is done by means of an ultra-wideband (UWB) measurement and the subsequent evaluation and processing of the waveform in the time and frequency domains.
Modern regulations for exhaust gas components like Euro 5 and Euro 6 require new methods both for the emissions control and for the monitoring and operation of the emission control systems in the diesel and gasoline engines.
During diesel engine operation, soot is captured in the pores of the DPFs (diesel particulate filters) or the SCR (selective catalytic reduction) catalysts. The ceramic filters must be periodically regenerated to avoid blocking. Hence, knowledge about the soot load is useful for the engine control in a number of ways (Feulner et al., 2013; Rose and Boger, 2009; Ochs et al., 2010; Sappok et al., 2010).
For gasoline engines, a TWC (three-way catalytic converter) is used to keep the polluting emissions low. The TWC stores and releases oxygen depending on the operating conditions. As in the DPF case, it would be of great advantage from the control engineering point of view to know the oxygen loading. However, it is the state of the art that the loading state cannot be measured directly but is derived indirectly from the oxygen concentration in the exhaust gas stream before and after the catalytic converter. The oxygen-measuring elements (so-called lambda probes) are implemented for the loading state estimation (Twigg, 2007).
There are a few known approaches for the direct measurement of the desired internal parameters (Zimmermann, 2007; Reiß et al., 2009). These and other similar ones also based on physical sensors can deliver only local information about the catalyst. However, one can learn nothing about the global state, such as the overall loading degree. Additionally, these approaches suffer from the complexities involving the mounting of the sensors in the catalyst and the communication with them.
Since the catalytic converter is typically located in a metallic housing, it
and its housing together act as a filled cavity resonator. The resonances
are influenced by the changes of material parameters in the catalyst so that
the catalyst state can be determined from the appropriate signal
characteristics of the resonances. The feasibility of this approach is well
documented in the literature (Fischerauer et al., 2008,
2010b; Eichelbaum et al., 2012; Moos et al., 2013). Typically, to observe
the cavity resonances, one or two simple thin probe feeds (short stubs) are
connected to an automatic vector network analyzer (VNA) via coaxial lines.
The changes in the resonant cavity are observed via the scattering matrix
parameters
It has also been shown that the microwave-based direct observation of TWCs enables superior engine control strategies compared to the indirect observation by way of lambda probes. This superiority consists, for instance, in lower emission levels (Schödel et al., 2014a, b).
The described laboratory equipment cannot be used in the field, i.e., inside a vehicle or in the manufacturing plant, owing to space and cost reasons. For the practical application, it should be replaced by smaller and cheaper solutions.
For solving this miniaturizing task, we have described a novel approach
towards the monitoring of exhaust gas aftertreatment systems such as DPFs,
SCRs, and TWCs (Motroniuk et al., 2014). It was proposed to consider the
interior of the catalyst metal housing as a wireless communication channel.
The effect of the catalyst state on such communication parameters as the bit
error rate (BER), the packet error rate (PER), the ratio of the energy per
bit to the noise power spectral density
In a recent conference contribution, we proposed a modification of our novel
approach. Instead of measuring certain communication channel parameters
(such as BER, PER,
In this contribution, we consider an extension of the UWB approach for the catalyst state determination. It is based on the subsequent evaluation and post-processing of the UWB waveforms in the time and frequency domains.
The rest of the paper is organized as follows. In Sect. 2, a novel wireless approach is presented, followed by a description of a hardware implementation and of an experimental verification with DPF in the time and frequency domains in Sects. 3 and 4, respectively. The conclusions are presented in Sect. 5.
It is proposed to consider the interior of the catalyst housing as a communication channel between two terminals of a UWB wireless communication system. The most important parameters during the data transmission are the characteristics of the communication channel (in our case the catalyst). From the wireless channel characteristics, the properties of the catalytic converter such as the catalyst oxygen concentration or soot loading state can be determined. The latter assertion about the measurability of the soot loading state has been demonstrated in Motroniuk et al. (2014) and will be further corroborated in the remainder of this contribution. The former assertion about the measurability of the oxygen concentration follows from the fact that the electrical properties of catalysts depend heavily on their oxidization state; hence, this state can be inferred from the transmission or reflection details of electromagnetic microwaves (Fischerauer et al., 2008; Eichelbaum et al., 2012).
Basic structure of the proposed measurement system.
The basic structure of the proposed measurement system is shown in Fig. 1. The system consists of a conductive (usually metallic) housing of the catalyst (1) and the catalyst itself (2). Two gas-permeable conductive grids (3, 4) which limit the resonating cavity are of advantage but not absolutely necessary. The grids are effectively short-circuiting the electric field, thus precisely defining the resonating cavity. Two thin-probe feeds (short stubs 5, 6) serve to connect the cavity with the system environment, i.e., with communication end devices (7, 8), by coaxial cables, which are used for signal transmission and reception. The short stubs play the same role for the communication system as the antennas in the common message transmission and the catalyst in its housing plays the same role as the wireless communication channel for the ordinary message transmission.
The signal transmitted through the wireless channel inside the catalyst suffers from such effects as attenuation, fading, multipath propagation, scattering, etc. Note, however, that the sender and receiver in a common wireless transmission link are located in the far field of each other whereas, in the present context, they are located in the near field. This, together with the resonance phenomena in the catalyst-filled cavity, introduces new propagation-path effects compared to wireless transmission links. The characteristics of the received information transferred from the transmitter to the receiver depend on the propagation medium. The effect on the signal is similar to the effect of buildings, trees, autos, rain, humidity, etc. in common wireless communication links: the channel characteristics influence the signal received, either in the time or frequency domain. Here, the closed-space propagation through the soot-loaded or otherwise chemically altered catalyst medium affects the received signal waveform.
In this paper, we specifically focus on a UWB communication system. Rather
than measuring individual channel parameters, such as BER, PER,
The goal, of course, is to estimate the internal state of a catalyst by this approach in a manner substantially simpler than possible with the system architectures demonstrated so far. This is achieved because the measurement can be performed by the existing hardware (wireless communication modules) and software (signal processing techniques for the wireless communication systems).
Functional diagram of the proposed measurement system.
The functional diagram of such a system is shown in Fig. 2. A pulse generator generates the UWB waveform, which is then transmitted through the device under test (catalyst) and distorted in this process. On the receiver side, the UWB pulse waveform is measured to indirectly estimate the catalyst state.
Exemplary hardware implementation of the proposed measurement system.
Figure 3 depicts some possible ways to implement the architecture proposed in actual hardware. The UWB pulse generation and the transmission are accomplished by a communication module (we have used a Time Domain PulsON P410 module), and the receiver is implemented by another communication module (also PulsON P410). After receiving the pulse modulated by the catalyst properties, it is a matter of experience (calibration, simulation results, etc.) to infer the catalyst state from the waveform distortion.
Devices such as Altera Stratix, Xilinx Virtex, Xilinx Spartan, ZigBee, WirelessHART, Nanonet, PulsON 410, or any other including (but not limited to) field-programmable gate-array (FPGA)-based integrated circuits and devices could be used for building an appropriate communication system for the intended application purpose. The catalyst state can be estimated by a post-processing of the received UWB waveform with the numerical computing environments and programming languages (such as MATLAB and C).
The validity of the approach was tested with DPFs loaded with different amounts of soot. The aim of the tests was to compare the effect of the soot load on the communication channel and, as a result, on the received UWB waveforms in order to prove the validity of the, as we believe, novel approach towards the monitoring and estimation of catalyst states.
Soot load of the measured DPFs (Fischerauer et al., 2010a).
One of the DPFs was clean (not used, with 0 g of soot), another one was loaded with 4.6 g of soot, and the third one was loaded with 13.6 g of soot. The filters used and their soot loading are listed in Table 1. For details of the catalyst geometry, the reader is referred to Fischerauer et al. (2010a).
UWB waveforms for three DPFs loaded with different amounts of soot.
The UWB pulse waveforms recorded during the communication between two commercial PulsON P410 modules via the interior of the DPF-loaded metal housing look substantially different for the various DPFs (Fig. 4). In particular, multiple-transit signals inside the catalyst housing (occurring after 18 ns in Fig. 4) are more attenuated the higher the soot load is. This is explained by the fact that the conducting soot introduces losses into the system. Thus, a wave packet is attenuated each time it traverses the soot-loaded DPF, and this attenuates multiple-transit signals more than the single-transit (direct input-to-output) signal. In Fig. 4 the vertical scale represents A–D (analog-to-digital) counts proportional to volts as output by the analog-to-digital converter (ADC) in the UWB chip.
The envelopes
Of course, all signals are discrete in time so that Eq. (1) had to be evaluated
in the form
Envelopes of the time-domain responses from Fig. 4.
The result is shown in Fig. 5. One clearly observes a substantial difference in the signal envelopes depending on the soot content in the DPF. In particular, this can be seen starting from the second peek at around 16 ns. The soot loading state is not observable via the first peak because of the technical characteristics of the equipment (at the beginning of the received UWB pulse, the Tx pulse coupling through the Tx switch is showing through).
To quantify the differences, we calculated the convolution of the envelope
functions for different DPF states, as always defined by
Convolution of the envelopes of the time-domain responses (Fig. 5).
The difference is observable by the width of the convolution curves as well.
Measured at an amplitude of 3.28
We also calculated the cross-correlation of the envelope functions, given by
Auto-correlation functions of the envelopes shown in Fig. 5.
For
Cross-correlation functions of the envelopes shown in Fig. 5.
Finally, Fig. 8 shows the cross-correlation functions. The amplitudes and
mean values of the peaks are among the parameters suitable for a DPF state
estimation algorithm. The maximum amplitude was 2.2
Smoothed envelopes of the time-domain responses from Fig. 4.
Owing to the rapid amplitude variation of the envelope curves in Fig. 5, it makes sense to first smooth them before further processing (Fig. 9). The amplitude of the smoothed envelopes depends on the soot load – with increasing soot load, the amplitude measured at the receiver becomes lower.
One possible way of processing is to correlate certain features of these smoothed curves with the DPF soot load. This is demonstrated in Fig. 10 for the center of gravity of the areas bounded by the curves in Fig. 9 and the time axis. It is obvious that this characteristic curve may be used to estimate the soot load from the measured time-domain response of the UWB system. It is also obvious that the sensitivity of this measurement system decreases with increasing soot load.
Time coordinate of the center of gravity of the area bounded by the smoothed envelopes and the time axis in Fig. 9.
Along the same line of reasoning, we have computed the characteristic curves linking the maximum and the average values of the smoothed envelopes with the DPF soot load (Fig. 11).
Maximum and average values of the smoothed envelopes (Fig. 9).
While the overall behavior of the characteristic curves in Figs. 10 and 11 is quite similar, there are differences in the details. It appears as if small soot loads can be better resolved (measured with a higher sensitivity) by way of the center-of-gravity feature (Fig. 10) than by the other features. This would be quite interesting, but, in our opinion, it is too early to draw such conclusions as the number of measurements we have performed so far is much too small to allow for any substantial insight into the statistics of the measurement approach.
Besides, it is also important to note that the curves will look different with different transmitting power because of receiver sensitivity. This, too, will have to be investigated in more detail.
According to the results presented in Sect. 4.1, the state of automotive aftertreatment systems can possibly be inferred from the time-domain characteristics of UWB signals. However, in some cases, it may be easier to characterize and evaluate signals in the frequency domain. This will be considered in the present section.
UWB spectra for the three waveforms from Fig. 4 are shown in Fig. 12 (the
spectral range from 3.1 to 5.3 GHz corresponds to the product specifications
of the communication units used). One clearly observes substantial
differences in the spectral amplitudes between the soot loading degrees,
especially at the resonance peaks of 3.18, 3.4, 3.57, 4.0, 4.16, 4.4 GHz,
etc. For example, the magnitude of the UWB spectrum at 3.4 GHz amounts to
7.6
UWB magnitude spectrum for three different DPF loading states corresponding to the time-domain signals in Fig. 4.
The convolution of the UWB magnitude spectra from Fig. 12 with each other is
calculated with
Convolution of the UWB magnitude spectra from Fig. 12.
Figure 14 shows the auto-correlation functions of the UWB magnitude spectra
from Fig. 12, calculated with
Auto-correlation functions of the UWB magnitude spectra in Fig. 12.
Cross-correlation functions of the UWB magnitude spectra in Fig. 12.
Figure 15 presents the cross-correlation functions of the UWB magnitude
spectra, Eq. (5) with
Furthermore, we have prepared histograms of the spectral magnitudes, measured in A–D counts, for three different DPF loading states. The results are presented in Fig. 16. The frequency range evaluated was 3.1–5.3 GHz, consistent with the product specifications of the communication units used. The area of every rectangle in the histogram represents the relative frequency of the spectral magnitudes falling into the interval marked by the rectangle basis. The height of the rectangle is the relative frequency density (the relative frequency of the spectral magnitudes associated with the rectangle interval, divided by the width of the rectangle).
An appropriate empirical probability density function (pdf) used to fit the
actually observed histogram must be one-sided and supported on the interval
[0,
Histogram of the UWB magnitude spectrum for different loading
states (bars) and best fits by exponential distributions
The exponential distributions with best-fitting parameters
One would expect that more soot load (i.e., higher attenuation) would
smooth the spectrum in the sense that the magnitude variations become
smaller; this would correspond to a pdf graph with steeper skirts and
smaller width. In numbers, this would mean that the pdf parameter
Maximum and average magnitudes of the signal spectrum at the receiver as functions of the soot load (from Fig. 12).
Finally, in analogy with the approach visualized for the time domain in
Fig. 11, Fig. 17 depicts the maximum and average values for the magnitude
spectrum in the range between 3.1 and 5.1 GHz. As before, the curves look
promising in that they are monotonic and would allow for a unique inversion of
the form “signal feature
We have reported on a novel approach for the monitoring of exhaust gas aftertreatment systems, such as DPFs, SCRs, and TWCs, the most attractive feature of which is that it is based on a cost-effective modular architecture. Hence, it can be implemented in situ, e.g., in a car.
In the approach proposed, the interior of the catalyst metal housing serves as a communication channel between two wireless communication nodes. To be more specific, we investigated a UWB-based architecture. It was demonstrated experimentally by way of various soot-loaded DPFs that the catalyst chemistry (i. e., the soot load of the DPF) affects the channel characteristics strongly enough to allow for the identification of the DPF state from the UWB signal details in the time as well as the frequency domains.
In particular, we were able to show that, in principle, certain features of the received time- and frequency-domain signals depend in a manner on the soot load which allows for the inversion of the feature vs. soot load characteristics. Among these features were the position of the center of gravity of the area bounded by the signal envelope and the time axis, the maximal value, the average value of the envelopes in the time domain and the maximum and average values in the spectra.
The authors thankfully acknowledge financial support by the German Research Foundation (DFG) under grants Fi 956/3-2 and Fi 956/5-1. Edited by: F. L. Dickert Reviewed by: two anonymous referees