KINETICS-BASED SCR CONTROL MODEL IMPROVEMENT

Abstract

A method for monitoring a discrete substrate element from an ammonia-selective catalyst reduction device configured to treat an exhaust gas feedstream of an internal combustion engine includes monitoring amounts of ammonia that are adsorbed, desorbed, and oxidized and an amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream from the discrete substrate element. An amount of ammonia consumption for the discrete substrate element is determined based on the amount of ammonia that is oxidized and the amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream. The amount of ammonia that is adsorbed and the amount of ammonia that is desorbed for the discrete substrate element are compared and the amount of ammonia consumption for the discrete substrate element is adjusted when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed.

Description

TECHNICAL FIELD

This disclosure is related to exhaust aftertreatment systems for internal combustion engines.

BACKGROUND

The statements in this section merely provide background information related to the present disclosure. Accordingly, such statements are not intended to constitute an admission of prior art.

Known engine control strategies to improve fuel economy and reduce fuel consumption in internal combustion engines include operating at a lean air/fuel ratio. This includes control strategies for engines configured to operate in compression-ignition and lean-burn spark-ignition combustion modes. Engines operating at lean air/fuel ratios can have increased local combustion temperatures leading to increased NOx emissions.

A known exhaust aftertreatment system and control strategy for managing and reducing NOx emissions includes a urea injection control system and an associated ammonia-selective catalytic reduction device. The urea injection control system injects a reductant, e.g., urea into an exhaust gas feedstream upstream of the ammonia-selective catalytic reduction device. The injected urea decomposes to ammonia, which reacts with NOx in the presence of a catalyst to produce nitrogen. Some amount of ammonia can be stored on the ammonia-selective catalytic reduction device, enabling continued NOx reduction when the urea injection control system is not capable of dispensing a controlled amount of urea. Known control systems include dispensing urea at a rate that corresponds to concentrations of engine-out NOx emissions to achieve NOx reduction without using excess amounts of urea, i.e., at a urea/NOx stoichiometric ratio.

SUMMARY

A method for monitoring a discrete substrate element from an ammonia-selective catalyst reduction device configured to treat an exhaust gas feedstream of an internal combustion engine includes monitoring an amount of ammonia that is adsorbed, an amount of ammonia that is desorbed, an amount of ammonia that is oxidized and an amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream from the discrete substrate element. An amount of ammonia consumption for the discrete substrate element is determined based on the amount of ammonia that is oxidized and the amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream. The amount of ammonia that is adsorbed and the amount of ammonia that is desorbed for the discrete substrate element are compared and the amount of ammonia consumption for the discrete substrate element is adjusted when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed.

BRIEF DESCRIPTION OF THE DRAWINGS

One or more embodiments will now be described, by way of example, with reference to the accompanying drawings, in which:

FIG. 1 illustrates an ammonia-selective catalytic reduction device of an exhaust aftertreatment system, in accordance with the present disclosure;

FIG. 2 illustrates a flowchart 100 for sequentially determining ammonia storage in a stepwise fashion for each of a plurality of discrete substrate elements over a time period for an ammonia-selective catalytic reduction device, in accordance with the present disclosure; and

FIGS. 3 and 4 illustrates a model for maintaining or adjusting a determined amount of ammonia consumption for a discrete substrate element per volume of gases passing through the discrete substrate element, in accordance with the present disclosure; and

FIGS. 5 and 6 illustrate experimental and derived data of an exemplary ammonia-selective catalytic reduction (SCR) device, depicting an estimated amount of NOx at an outlet of the ammonia-SCR device and a measured amount of NOx at the outlet of the ammonia-SCR device over a period of time, in accordance with the present disclosure.

DETAILED DESCRIPTION

Referring now to the drawings, wherein the showings are for the purpose of illustrating certain exemplary embodiments only and not for the purpose of limiting the same, FIG. 1 schematically illustrates an exhaust aftertreatment device 60 for treating an exhaust gas feedstream output from an internal combustion engine, including an internal combustion engine that operates lean of stoichiometry. It is appreciated that the internal combustion engine can include an engine configured to operate in a compression-ignition combustion mode, an engine configured to operate in a lean-burn spark-ignition combustion mode, and an engine configured to operate in a homogeneous-charge compression ignition combustion mode. The exhaust aftertreatment device 60 described herein is an ammonia-selective catalytic reduction device (NH3-SCR device) including a coated substrate 52.

It is appreciated that the NH3-SCR device 60 is an element of an exhaust aftertreatment system that may include other aftertreatment devices. In one embodiment the exhaust aftertreatment system can include a plurality of aftertreatment devices each of which includes a device that employs technologies having various capabilities for treating the constituent elements of the exhaust gas feedstream. Treating the constituents of the exhaust gas feedstream may include oxidation, selective catalytic reduction using a reductant, particulate filtering, and other treatments. Design features for each aftertreatment device include total volume, space velocity, cell density, washcoat materials, loading(s) of catalytic material(s), and vehicle/engine compartment locations, which are determined for specific applications. In one embodiment, a first aftertreatment device is a three-way catalyst that is located upstream of a second aftertreatment device (i.e. NH3-SCR device 60), which is located upstream of a third aftertreatment device that includes a catalyzed particulate filter, although the concepts described herein are not so limited. The first, second, and third aftertreatment devices are fluidly connected in series using known pipes and connectors. The first, second, and third aftertreatment devices can be assembled into individual structures that are fluidly connected and assembled in an engine compartment and a vehicle underbody with one or more sensing devices placed therebetween.

The exemplary NH3-SCR device 60 is depicted using a two-dimensional schematic model with an exhaust gas feedstream flowing therethrough. In one embodiment there is a urea injection device 20 and associated urea delivery system upstream of a mixer device 25 that is upstream of the NH3-SCR device 60. The NH3-SCR device 60 includes one or more ceramic coated substrates 52 preferably fabricated from cordierite material and having a multiplicity of flowthrough passageways that are coated with washcoat and catalytic materials to store ammonia for reacting with NOx molecules present in the exhaust gas feedstream. It is appreciated that ammonia storage concentration (θ_NH3) may be unevenly distributed along a flow axis of the coated substrate(s) 52.

A control module 10 is configured to monitor and control engine operation and monitor the exhaust gas feedstream. The control module 10 monitors or otherwise determines states of parameters of the exhaust gas feedstream. The control module 10 controls operation of the urea injection device 20. The control module 10 includes a virtual sensor 55 configured to estimate ammonia storage concentration (θ_NH3) on the coated substrate 52 during ongoing operation of the engine. The virtual sensor 55 is achieved by executing routines and a plurality of predetermined calibration arrays that temporally determine the ammonia storage concentration (θ_NH3) stored on the coated substrate 52. The virtual sensor 55 is described in detail with reference to FIG. 2. The ammonia storage concentration (θ_NH3) is preferably expressed as a ratio of stored ammonia to a maximum ammonia storage capacity for the coated substrate 52 or portion thereof, and can be in terms of mass/mass, moles/moles, or other suitable measurement.

The control module 10 is configured to monitor or otherwise determine states of parameters of the exhaust gas feedstream flowing into the NH3-SCR device 60. Preferred parameters of the exhaust gas feedstream include an inlet temperature of the exhaust gas feedstream, pressure, mass flowrate, oxygen concentration, NOx concentrations, and other parameters from which concentrations of input gases including nitrogen oxide, nitrogen dioxide, nitrous oxide, oxygen, and ammonia can be determined, as is appreciated by one skilled in the art. The substrate temperature T_sub can be monitored with a temperature sensor or determined by executing a mathematical model based upon the parameters of the exhaust gas feedstream and catalytic reaction rates of the coated substrate 52.

Control module, module, control, controller, control unit, processor and similar terms mean any one or various combinations of one or more of Application Specific Integrated Circuit(s) (ASIC), electronic circuit(s), central processing unit(s) (preferably microprocessor(s)) and associated memory and storage (read only, programmable read only, random access, hard drive, etc.) executing one or more software or firmware programs or routines, combinational logic circuit(s), input/output circuit(s) and devices, appropriate signal conditioning and buffer circuitry, and other components to provide the described functionality. Software, firmware, programs, instructions, routines, code, algorithms and similar terms mean any controller executable instruction sets including calibrations and look-up tables. The control module has a set of control routines executed to provide the desired functions. Routines are executed, such as by a central processing unit, and are operable to monitor inputs from sensing devices and other networked control modules, and execute control and diagnostic routines to control operation of actuators. Routines may be executed at regular intervals, for example each 3.125, 6.25, 12.5, 25 and 100 milliseconds during ongoing engine and vehicle operation.

The coated substrate(s) 52 is analytically segmented into a plurality of discrete elements 52(i), i=1 through n, or bricks, positioned in series along a flow axis of the exhaust gas feedstream. Segmenting the coated substrate(s) 52 into the plurality of discrete elements provides an analytical framework for implementing the virtual sensor 55 to determine the ammonia storage concentration (θ_NH3) on the NH3-SCR device 60 in real-time during ongoing operation of the engine.

FIG. 1 shows a single discrete substrate element 52(i). Each discrete substrate element 52(i) is characterized in terms of substrate operating temperature (T_sub) and the ammonia storage concentration (θ_NH3). The exhaust gas feedstream passing through each discrete substrate element 52(i) is characterized in terms of concentrations of input gases 50 including nitrogen oxide [NO]in, nitrogen dioxide [NO2]in, nitrous oxide [N2O]in, oxygen [O2]in, and ammonia [NH3]in, and corresponding output gases 54 including nitrogen oxide [NO], nitrogen dioxide [NO2], nitrous oxide [N2O], oxygen [O2], and ammonia [NH3].

FIG. 2 illustrates a flowchart 100 for determining ammonia storage in a stepwise fashion for each of a plurality of discrete substrate elements over a time period for an ammonia-selective catalytic reduction device in accordance with the present disclosure. The flowchart 100 is an element of the virtual sensor 55 that is executed to determine ammonia storage concentration (θ_NH3) for the total coated substrate 52. Other elements of the virtual sensor 55 are described herein. It is appreciated that the virtual sensor 55 can be reduced to a routine and executed in the control module 10 during ongoing engine operation to determine the ammonia storage concentration (θ_NH3) for the entire coated substrate 52 in real-time. Table 1 is provided as a key to a flowchart 100 of FIG. 2 wherein the numerically labeled blocks and the corresponding functions are set forth as follows.

TABLE 1
BLOCK BLOCK CONTENTS
105   Monitoring states of input parameters 150 of exhaust gas
      feedstream including [NO]in, nitrogen dioxide [NO2]in,
      nitrous oxide [N20]in, oxygen [O2]in, ammonia and
      substrate temperature for each time period Δt.
110   Monitoring the states of input parameters for each discrete
      substrate element 52(i).
115   Determining for each discrete substrate element 52(i) an
      amount of ammonia that is adsorbed, Δ[NH3]adsorption. The
      Δ[NH3]adsorption can also be described as an ammonia
      adsorption rate.
120   Determining for each discrete substrate element 52(i) an
      amount of ammonia that is desorbed, Δ[NH3]desorption. The
      Δ[NH3]desorption can also be described as an ammonia
      desorption rate.
125   Determining for each discrete substrate element 52(i) an
      amount of ammonia that is oxidized, Δ[NH3]oxidation. The
      Δ[NH3]oxidation can also be described as an ammonia
      oxidation rate.
130   Determining for each discrete substrate element 52(i) an
      amount of ammonia that is consumed,
      Δ[NH3]NOx—conversion, during reduction of NOx in the
      exhaust gas feedstream. The Δ[NH3]NOx—conversion can also
      be described as an ammonia conversion rate.
140   Determining a change in ammonia storage concentration
      (θNH3(i)) and concentrations of other chemical species for
      each of the discrete substrate elements 52(i)
145   Has the change in ammonia storage concentration (θNH3(i))
      for the last discrete substrate element 52(i) been determined?

The flowchart 100 includes determining a change in ammonia storage for each of the discrete substrate elements 52(i) over a time period and then determining the total ammonia storage concentration (θ_NH3) on the coated substrate 52 based thereon. Determining a change in the ammonia storage concentration (θ_NH3) includes sequentially determining a change in ammonia storage in a stepwise fashion for each of the discrete substrate elements 52(i), i=1 through n, over a time period Δt based upon the concentrations of the input gases of nitrogen oxide [NO]in, nitrogen dioxide [NO2]in, nitrous oxide [N2O]in, oxygen [O2]in, and ammonia [NH3]in and substrate temperature. This includes determining, for each discrete substrate element 52(i) (110) for each time period Δt (105), an amount of ammonia that is adsorbed (115), an amount of ammonia that is desorbed (120), an amount of ammonia that is oxidized (125), and an amount of ammonia that is consumed during reduction of NOx in the exhaust gas feedstream (130). The amounts of ammonia that is adsorbed (115), desorbed (120), oxidized (125), and consumed during reduction of NOx (130) can be in any suitable units of measure, including, e.g., mass, volume, or moles.

The change in the ammonia storage concentration (θ_NH3) and concentrations of other chemical species are determined in a stepwise fashion for each of the discrete substrate elements 52(i) using the foregoing blocks 115, 120, 125 and 130 (140), which are repeated for each of the discrete substrate elements 52(i) for each time period Δt. An output for the discrete substrate element 52(i) is determined that includes corresponding concentrations of output gases 54 of nitrogen oxide [NO], nitrogen dioxide [NO2], nitrous oxide [N2O], ammonia [NH3], oxygen [O2], and a cumulative ammonia storage concentration [θ_NH3]. The control module 10 can use the information to control engine fueling and air/fuel ratio for the exemplary powertrain system of FIG. 1 based upon inputs including the ammonia storage concentration (θ_NH3) on the coated substrate 52 (NH3_storage).

The following set of relationships represents reaction chemistry occurring in each of the discrete substrate elements 52(i) of the coated substrate 52.

[1]

4NH₃+4NO+O₂=4N₂+6H₂O  (A)

2NH₃+NO+NO₂=2N₂+3H₂O  (B)

8NH₃+6NO₂=7N₂+12H₂O  (C)

4NH₃+3O₂=2N₂+6H₂O  (D)

4NH₃+5O₂=4NO+6H₂O  (E)

4NH₃+4NO+3O₂=4N₂O+6H₂O  (F)

2NH₃+2NO₂═N₂O+N₂+3H₂O  (G)

2NH₃+2O₂═N₂O+3H₂O  (H)

The kinetic reactions and ammonia adsorption and desorption occur on the catalyst surface of one of the discrete substrate elements 52(i). The ammonia storage results from dynamic balance among the adsorption, desorption and kinetic reaction rates. For each of the discrete substrate elements 52(i), i=1 to n, the inlet parameter values are the outlet parameter values of the contiguous upstream discrete substrate element 52(i−1). The kinetic reaction rates rely on the ammonia storage concentration (θ_NH3) and the other related chemical concentrations.

The amount of ammonia that is adsorbed (115), i.e., ammonia adsorption rate, can be determined in accordance with the following relationships:

$\begin{array}{cc}{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{adsorption}}=\mathrm{\eta \_adsorption}\left({\left[\mathrm{NH}3\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[\mathrm{NH}3\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{desorption}}\right)& \left[2\right]\end{array}$

wherein an adsorption efficiency term η_adsorption is preferably selected from a predetermined array F_{table—adsorp}(T_sub, ξ_adsorp) that is stored in tabular form in the control module 10. A specific value for the adsorption efficiency term η_adsorption correlates to substrate temperature T_sub and an adsorption capacity term ξ_adsorp, which are set forth in the following relationships:

$\begin{array}{cc}{\eta }_{\mathrm{adsorption}}=F_{\mathrm{table}\_\mathrm{adsorp}}\left(T_{\mathrm{sub}},{\xi }_{\mathrm{adsorp}}\right){\xi }_{\mathrm{adsorp}}=\left(1-{\theta }_{\mathrm{NH}3}\right)*\Omega *\frac{t_{\mathrm{resident}}}{\frac{t_{\mathrm{resident}}}{\Delta t}+1}& \left[3\right]\end{array}$

In the above relationships [1] and [2]:

• • [NH3]_−Δt is the NH3 concentration in the discrete substrate element 52(i) at a previous timestep;
  • [NH3]_in is the NH3 concentration at the inlet to the discrete substrate element 52(i);
  • T_sub is the substrate temperature of the discrete substrate element 52(i); Δt is the time period;
  • θ_NH3 is the ammonia storage concentration for the discrete substrate element 52(i);
  • t_resident is the gas resident time, which can be determined based upon the volume of the discrete substrate element 52(i) and the volumetric flowrate of the exhaust gas feedstream; and
  • Ω is a specific ammonia storage capacity for the discrete substrate element 52(i), which is preferably stored in the control module 10, and is considered a constant. The specific ammonia storage capacity can be in any suitable units of measure, including, e.g., mass, volume, or moles, and is preferably consistent with other measurements and estimates of ammonia storage capacity.

With known states for each of the aforementioned parameters, i.e., [NH3]_in, [NH3]_−Δt, Δ[NH3]_desorption, T_sub, θ_NH3, and t_resident the amount of ammonia that is adsorbed in the discrete substrate element 52(i), i.e., Δ[NH3]_adsorption can be determined.

The amount of ammonia that is desorbed, i.e., Δ[NH3]_desorption (120), i.e., an ammonia desorption rate, can be determined using the specific ammonia storage capacity for the discrete substrate element 52(i) Ω, the residence time t_resident, and the ammonia storage concentration (θ_NH3) for the discrete substrate element 52(i) in combination with a predetermined desorption term F_{table—desorp}(T_sub,θ_NH3) as set forth below in Eq. 4. The predetermined desorption term F_{table—desorp}(T_sub,θ_NH3) is selected from a predetermined array of values stored in a memory lookup table, and is associated with the substrate temperature T_sub and ammonia storage concentration (θ_NH3) for the discrete substrate element 52(i).

Δ[NH3]_desorption=F_{table—desorp}(T_sub,θ_NH3)*θ_NH3*Ω*t_resident  [4]

The amount of ammonia that is oxidized, i.e., Δ[NH3]_oxidation (125), i.e., an ammonia oxidation rate, can be determined in accordance with the following relationship:

Δ[NH3]_oxidation=Δ[NH3]_oxid—N2+Δ[NH3]_oxid—NO+Δ[NH3]_oxid—N2O  [5]

The terms of Eq. 5 include an amount of ammonia oxidized in forming nitrogen, i.e., Δ[NH3]_oxid—N2, an amount of ammonia oxidized in forming NO, i.e., Δ[NH3]_oxid—NO, and an amount of ammonia oxidized in forming N₂O, i.e., Δ[NH3]_oxid—N2O, which can be determined as described below. The aforementioned terms include predetermined oxidation terms F_{table—oxid—N2}, F_{table—oxid—NO} and F_{table—oxid—N2O} that are selected from corresponding predetermined arrays that are preferably stored in tabular form in the control module 10. Specific values for each of the predetermined oxidation terms correspond to the substrate temperature T_sub and ammonia storage concentration (θ_NH3) for the discrete substrate element 52(i) in accordance with the following relationships:

[6]

Δ[NH3]_oxid—N2=F_{table—oxid—N2}(T_sub,θNH3)*[O₂]*Ω*t_resident  (A)

Δ[NH3]_oxid—NO=F_{table—oxid—NO}(T_sub,θ_NH3)*[O₂]*Ω*t_resident  (B)

Δ[NH3]_oxid—N2O=F_{table—oxid—N2O}(T_sub,θ_NH3)*[O₂]Ω*t_resident  (C)

wherein [O₂] is oxygen concentration,

• • t_resident is a gas resident time in the discrete substrate element 52(i),
  • θ_NH3 is the ammonia storage concentration, and
  • Ω is the specific ammonia storage capacity for the discrete substrate element 52(i).

The three oxidation reactions set forth in Eq. 6 correspond to oxidation rates related to (A) Δ[NH3]_oxid—N2 and (B) Δ[NH3]_oxid—NO of Eq. 1, respectively, to represent an amount of ammonia that is consumed thereby.

The amount of ammonia that is consumed for NOx reduction (130), i.e., an ammonia conversion rate, can be determined in accordance with the following relationship.

$\left[7\right]$ $\begin{array}{c}{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{NOx}\_\mathrm{conversion}}=\left\{{\left[{\mathrm{NO}}_X\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_X\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{oxid}\_\mathrm{NO}}\right\}\times \\ \{\left(1-R_{{\mathrm{NO}}_2}\right)\left[1-\left(1-{\eta }_{{\mathrm{NO}}_X}\right)\left(1-{\eta }_{\mathrm{NO}}\right)\right]+\left(A\right)\\ \frac{4}{3}\left[R_{{\mathrm{NO}}_2}-{\eta }_{{\mathrm{NO}}_X}\left(1-R_{{\mathrm{NO}}_2}\right)\right]{\eta }_{{\mathrm{NO}}_2}+\left(B\right)\\ {\eta }_{{\mathrm{NO}}_X}\left(1-R_{{\mathrm{NO}}_2}\right)\}\left(C\right)\end{array}$

The terms set forth in Eq. 7 include reduction efficiency terms η_NO, η_NO2, and η_NOx that indicate efficiencies associated with consumption of ammonia in reducing NO, NO₂, and NO_x, respectively. The aforementioned efficiency terms are determined in accordance with the following relationships:

$\begin{array}{cc}{\eta }_{\mathrm{NO}}=F_{\mathrm{table}\_\mathrm{NO}}\left(T_{\mathrm{sub}},{\zeta }_{\mathrm{NO}}\right){\zeta }_{\mathrm{NO}}=f_{\mathrm{table}\_O2\_\mathrm{NO}}\left(\left[O_2\right]\right)f_{\mathrm{table}\_\theta \mathrm{abl}}\left({\theta }_{\mathrm{NH}3}\right)*\Omega *t_{\mathrm{resident}}/\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right){\eta }_{{\mathrm{NO}}_2}=F_{{\mathrm{table}\_\mathrm{NO}}_2}\left(T_{\mathrm{sub}},{\zeta }_{\mathrm{NO}2}\right){\zeta }_{\mathrm{NO}2}=f_{\mathrm{table}\_\theta \mathrm{abl}2}\left({\theta }_{\mathrm{NH}3}\right)*\Omega *t_{\mathrm{resident}}/\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right){\eta }_{\mathrm{NOx}}=F_{\mathrm{table}\_\mathrm{NOx}}\left({\eta }_{{\mathrm{NO}}_{x\_2}},R_{{\mathrm{NO}}_2}\right){\eta }_{{\mathrm{NO}}_{x\_2}}=f_{{\mathrm{table\_NO}}_{\mathrm{x\_}2}}\left({\eta }_{{\mathrm{NO}}_{x\_1}},{\left[{\mathrm{NO}}_x\right]}_{\mathrm{aver}}*\left(1-{\eta }_{\mathrm{NO}}\right)*\left(1-{\eta }_{{\mathrm{NO}}_2}\right)\right){\eta }_{{\mathrm{NO}}_{x\_1}}=f_{{\mathrm{table\_NO}}_{\mathrm{x\_}1}}\left(T_{\mathrm{sub}},{\zeta }_{\mathrm{NOx}}\right){\zeta }_{\mathrm{NOx}}=f_{\mathrm{table}\mathrm{\_\theta }\mathrm{abl}x}\left({\theta }_{\mathrm{NH}3}\right)\Omega t_{\mathrm{resident}}/\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right){\mathrm{wherein}\left[{\mathrm{NO}}_x\right]}_{\mathrm{aver}}=\frac{{\left[{\mathrm{NO}}_x\right]}_{\mathrm{aver}}\left({\left[{\mathrm{NO}}_x\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_x\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{oxid}\_\mathrm{NO}}\right)}{\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right)}R_{{\mathrm{NO}}_2}=\frac{\left({\left[{\mathrm{NO}}_2\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_2\right]}_{-\Delta t}\right)}{\left({\left[{\mathrm{NO}}_x\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_x\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{oxid}\_\mathrm{NO}}\right)}& \\ {\left[{\mathrm{NO}}_x\right]}_{in}={\left[\mathrm{NO}\right]}_{in}+{{\left[{\mathrm{NO}}_2\right]}_{in}\left[{\mathrm{NO}}_x\right]}_{-\Delta t}={\left[\mathrm{NO}\right]}_{-\Delta t}+{\left[{\mathrm{NO}}_2\right]}_{-\Delta t}& \end{array}$

The term R_NO2 denotes a ratio of NO₂/NOx in the incoming gas feedstream. The reduction efficiency terms η_NO, η_NO2, and η_NOx change with substrate temperature T_sub and consumption capacities of the substrate element 52 associated with NO, NO₂, and NO_x, respectively, which are expressed as terms ζ_NO, ζ_NO2, and ζ_NOx, respectively. Each of the consumption capacity terms associated with the discrete substrate element 52(i), i.e., ζ_NO, ζ_NO2, and ζ_NOx is calculated as a function of residence time t_resident, specific ammonia storage capacity Ω, and a related scaling term, i.e., f_{table—θ—NO2}(θ_NH3), f_{table—θ—NOx}(θ_NH3), f_{table—θ—NO}(θ_NH3) and f_{table—O2—NO}([O₂]), each of which is a function of either the ammonia storage concentration (θ_NH3) or the oxygen concentration (O₂). As is appreciated, the consumption capacities of the substrate element 52 are based upon a relationship between the ammonia storage concentration (θ_NH3), ammonia storage capacity Ω, and residence time t_resident of the exhaust gas in the discrete substrate element 52(i). The consumption capacities of the substrate element 52 associated with NO, NO₂, and NO_x, i.e., ζ_NO, ζ_NO2, and ζ_NOx, respectively, are preferably predetermined and stored as arrays in tabular form in the control module 10.

The reduction efficiency terms associated with NO, NO₂, and NO_x, i.e., η_NO, η_NO2, and η_NOx, respectively, are selected from predetermined arrays that are preferably stored in tabular form in the control module 10. The predetermined arrays are searchable in terms of the substrate temperature T_sub and the corresponding consumption capacity of the substrate, i.e., one of ζ_NO, ζ_NO2, and ζ_NOx.

Thus, the ammonia storage concentration (θ_NH3) can be determined (140) in accordance with the following relationship:

θ_NH3,t=θ_NH3,t-Δt+(Δ[NH3]_adsorption−Δ[NH3]_desorption−Δ[NH3]oxidation−Δ[NH3]NOx_conversionΔtΩtresident  [8]

• • wherein Δ[NH3]_adsorption includes an amount of ammonia adsorbed into a catalyst surface per volume of gases passing through the discrete substrate element 52(i),
  • Δ[NH3]_desorption includes an amount of ammonia desorbed from catalyst surface per volume of gases passing through the discrete substrate element 52(i),
  • Δ[NH3]_oxidation includes an amount of ammonia oxidized per volume of gases passing through the discrete substrate element 52(i), and
  • Δ[NH3]_{NOx—conversion} includes an amount of ammonia consumed for NOx reduction per volume of gases passing through the discrete substrate element 52(i).

The chemical species concentrations for the discrete substrate element 52(i) can be determined for NO, NO₂, ammonia, and N₂O concentrations in accordance with the following relationships.

$\begin{array}{cc}\left[\mathrm{NO}\right]=\frac{\begin{array}{c}\left(1-R_{{\mathrm{NO}}_2}\right)\left(1-{\eta }_{{\mathrm{NO}}_x}\right)\left(1-{\eta }_{\mathrm{NO}}\right)\\ \left({\left[{\mathrm{NO}}_x\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_x\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{oxid}\_\mathrm{NO}}\right)\end{array}}{\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right)}& \left[9\right]\\ \left[{\mathrm{NO}}_2\right]=\frac{\begin{array}{c}\left[R_{{\mathrm{NO}}_2}-{\eta }_{{\mathrm{NO}}_x}\left(1-R_{{\mathrm{NO}}_2}\right)\right]\left(1-{\eta }_{{\mathrm{NO}}_2}\right)\\ \left({\left[{\mathrm{NO}}_x\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[{\mathrm{NO}}_x\right]}_{-\Delta t}\right)\end{array}}{\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right)}& \left[10\right]\\ \left[\mathrm{NH}3\right]=\frac{\begin{array}{c}\left(1-{\eta }_{\mathrm{adsorption}}\right)\\ \left({\left[\mathrm{NH}3\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[\mathrm{NH}3\right]}_{-\Delta t}+{\Delta \left[\mathrm{NH}3\right]}_{\mathrm{desorption}}\right)\end{array}}{\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right)}& \left[11\right]\\ \left[N_2O\right]\frac{\left({\left[N_2O\right]}_{\mathrm{in}}+{\frac{t_{\mathrm{resident}}}{\Delta t}\left[N_2O\right]}_{-\Delta t}+\Delta \left[N_2O\right]\right)}{\left(1+\frac{t_{\mathrm{resident}}}{\Delta t}\right)}& \left[12\right]\end{array}$

• • wherein [NO]_−Δt, [NO₂]_−Δt, [N₂O]_−Δt and [NH3]_−Δt are the concentration values in the discrete substrate element 52(i) defined at the previous timestep for NO, NO₂, and N₂O.

Δ[N₂O]=Δ[NH3]_oxid—N2O+Δ[N₂O]_NO2+Δ[N₂O]_NO  [13]

Δ[N₂O]_NO=Δ[NO]_in(1−η_NOx)η_NOγ_{table—NO—N2O}(T_sub)  [14]

Δ[N₂O]_NO2=([NO₂]_in−η_NOx[NO]_in)η_{NO2γtable—NO2—N2O}(T_sub)  [15]

• • wherein γ_{table—NO—N2O} and γ_{table—NO2—N2O} are reaction rate terms that are selected from a predetermined array that is stored in tabular form in the control module 10. A specific value for each of the reaction rate terms is retrievable as a function of the substrate temperature T_sub.

Thus, the virtual sensor 55 can be used to determine ammonia storage concentration (θ_NH3) for the entire coated substrate 52 by sequentially determining a change in ammonia storage for each of the discrete substrate elements in a stepwise fashion for each of the discrete substrate elements 52(i), i=1 through n, over a time period, and determining the ammonia storage concentration (θ_NH3) on the ammonia-selective catalyst reduction device corresponding to the change in ammonia storage for the discrete substrate elements 52(i).

FIGS. 3 and 4 illustrate a model 300 for maintaining or adjusting a determined amount of ammonia consumption (Δ[NH3]_consumption) for the discrete substrate element 52(i) per volume of gases passing through the discrete substrate element 52(i). As aforementioned, the discrete substrate element 52(i) is one of a plurality of discrete substrate elements 52(n) analytically segmented from the ammonia-selective catalytic reduction (SCR) device 60 of FIG. 1. The Δ[NH3]_consumption can also be described as an ammonia consumption rate for the discrete substrate element 52(i) per volume of gases passing through the discrete substrate element 52(i). Discussed in greater detail below, the Δ[NH3]_consumption is based on the Δ[NH3]_oxidation (see Eq. [5]) and the Δ[NH3]_{NOx—conversion}(see Eq. [7]).

The model 300 is integrated within the virtual sensor 55 of the control module 10. The model 300 includes an ammonia adsorption module 215, an ammonia desorption module 220 and an ammonia consumption module 230. The model 300 further includes an integer unit 240, first and second comparison units 250, 260, respectively, a divider unit 270, a limiter module 280, an adjustment module 290, a maintaining module 2, a switch module 295 and a multiplication module 298.

The ammonia adsorption module 215 monitors an amount of ammonia adsorbed (Δ[NH3]_adsorbed 211) onto the surface of each discrete substrate element 52(i) per volume of gases passing through the discrete substrate element 52(i). The monitored Δ[NH3]_adsorbed 211 can be determined in block 115 of flowchart 100 utilizing Equation [2].

The ammonia desorption module 220 monitors an amount of ammonia desorbed (Δ[NH3]_desorbed) 221 from the surface of each discrete substrate element 52(i) per volume of gases passing through the discrete substrate element 52(i). The monitored Δ[NH3]_desorbed 221 can be determined in block 120 of flowchart 100 utilizing Equation [4].

The ammonia consumption module 230 monitors an amount of ammonia consumption Δ[NH3]_consumption 231 for the discrete substrate element 52(i) per volume of gases passing through the discrete substrate element 52(i). The ammonia consumption module 230 can determine the Δ[NH3]_consumption 231 based on calculating a sum of the Δ[NH3]_oxidation and the Δ[NH3]_{NOx—conversion. As aforementioned. the Δ[NH}3]_oxidation can be determined in block 125 of flowchart 100 utilizing Eq. [5] and the Δ[NH3]_{NOx—conversion} can be determined in block 130 of flowchart 100 utilizing Eq. [7]

It will be appreciated that the amount of ammonia adsorbed Δ[NH3]_adsorbed 211 can interchangeably be referred to as an ammonia adsorption rate, the amount of ammonia desorbed Δ[NH3]_desorbed 221 can interchangeably be referred to as an ammonia desorption rate, the amount of ammonia that is oxidized Δ[NH3]_oxidation can interchangeably be referred to as an ammonia oxidation rate, the amount of ammonia consumed in reducing NOx in the exhaust gas feedstream Δ[NH3]_{NOx—conversion} can interchangeable be referred to as an ammonia conversion rate and the amount of ammonia consumption Δ[NH3]_consumption 231 can interchangeably be referred to as an ammonia consumption rate.

The first comparison unit 250 compares the Δ[NH3]_adsorbed 211 and the Δ[NH3]_desorbed 221 for each discrete substrate element 52(i). As will become apparent, with reference to FIG. 4, the amount of ammonia consumption (i.e., ammonia consumption rate) Δ[NH3]_consumption 231 for each discrete substrate element 52(i) can be adjusted to determine an adjusted amount of ammonia consumption (i.e., adjusted ammonia consumption rate) 296 when the Δ[NH3]_adsorbed 211 is less than the Δ[NH3]_desorbed 221. On the other hand, with reference to FIG. 3, the Δ[NH3]_consumption 231 for each discrete element 52(i) can be maintained when the Δ[NH3]_adsorbed 211 is at least the Δ[NH3]_desorbed 221.

Situations during operation of the ammonia-SCR device 60 can arise that may result in the amount of ammonia that is adsorbed to be less than the amount of ammonia that is desorbed. In one embodiment, the amount of ammonia that is adsorbed can be less than the amount of ammonia that is desorbed when urea dosing or ammonia dosing from the dosing device 20 is momentarily stopped or terminated. Inaccuracies to the determined amount of ammonia consumption (i.e., sum of Δ[NH3]_oxidation and Δ[NH3]_{NOx—conversion}) may result when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed. Specifically, the amount of ammonia consumption Δ[NH3]_consumption 231 for each discrete substrate element 52(i) can deviate by an increased magnitude from an actual amount of ammonia consumption for each discrete substrate element 52(i) when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed. As evidenced by block 140 of flowchart 100, a determined amount of ammonia consumption that deviates from an actual amount of ammonia consumption would result in an inaccurate ammonia storage concentration (θ_NH3) in each discrete substrate element 52(i) and an inaccurate total ammonia storage concentration on the ammonia-selective catalyst reduction device 60. Thus, an inaccurate θ_NH3 can result in estimations of an amount of NOx at an outlet of each discrete substrate element 52(i) and an amount of NOx at an outlet of the ammonia-SCR device 60 deviating from measured NOx values at the respective outlets. If NOx estimations deviate from measured values, ammonia dosing or urea dosing rates cannot be accurately controlled, thereby increasing ammonia or urea consumption and diminishing efficiency of the ammonia-SCR device 60. Accordingly, a modification or adjustment to the virtual sensor 55 may be required to adjust the determined amount of ammonia consumption Δ[NH3]_consumption 231 to compensate for inaccuracies when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed.

When the Δ[NH3]_adsorbed 211 is less than the Δ[NH3]_desorbed 221, the Δ[NH3]_adsorbed 211 is input to the second comparison unit 260. In one embodiment, a determination that the Δ[NH3]_adsorbed 211 is less than the Δ[NH3]_desorbed 221 is indicative of a situation where urea dosing or ammonia dosing by the urea dosing device 20 has momentarily stopped or has been terminated.

If the second comparison unit 260 determines the Δ[NH3]_adsorbed 211 has a value greater than zero, the Δ[NH3]_adsorbed 211 is input to the divider unit 270. If the second comparison unit 260 determines the Δ[NH3]_adsorbed 211 has a value equal to zero, a predetermined non-zero integer 262 provided by the integer unit 240 is input to the divider unit 270 in place of the Δ[NH3]_adsorbed 211. In other words, when the Δ[NH3]_adsorbed 211 is equal to zero, the Δ[NH3]_adsorbed 211 is designated to be equal to the predetermined non-zero integer 262. In a non-limiting example, the predetermined non-zero integer 262 has a value equal to 1.0 e-10.

The divider unit 270 receives the Δ[NH3]_adsorbed 211 and the Δ[NH3]_consumption 231 and divides the Δ[NH3]_consumption 231 by the Δ[NH3]_adsorbed 211 to determine a ratio 272 for each discrete substrate element 52(i). As aforementioned, the Δ[NH3]_adsorbed 211 can be designated to be equal to the predetermined non-zero integer 262 when the Δ[NH3]_adsorbed 211 is equal to zero. Hence, the divider unit 270 determines a ratio 272 of the amount of ammonia consumption (i.e., ammonia consumption rate) 231 to the amount of ammonia that is adsorbed (i.e., ammonia adsorption rate) 215 for each discrete substrate element 52(i).

In the exemplary embodiment, the ratio 272 can be input to the limiter module 280. The limiter module 280 can proportionally reduce the ratio 272 within a predetermined range defined by a minimum limit and a maximum limit. In a non-limiting example, the proportionally reduced predetermined range is defined by and includes a minimum limit equal to “0” and a maximum limit equal to “20.” Accordingly, a reduced ratio 282 is provided by the limiter module 280 and is output to the adjustment module 290.

The adjustment module 290 is configured to determine an adjustment multiplier 292 based on the ratio (i.e., or the reduced ratio 282) of the amount of ammonia consumption (i.e., ammonia consumption rate) 231 to the amount of ammonia that is adsorbed (i.e., ammonia adsorption rate) 211. The adjustment multiplier 292 can be applied to the amount of ammonia consumption (i.e., ammonia consumption rate) 231 for each discrete substrate element 52(i) when the Δ[NH3]_adsorbed 211 is less than the Δ[NH3]_desorbed 221.

Referring to FIG. 4, the switch module 295 switches from a nominal state 293 to an adjustment state 294 when the first comparison unit 250 determines the Δ[NH3]_adsorbed 211 is less than the Δ[NH3]_desorbed 221. Accordingly, the adjustment multiplier 292 and the amount of ammonia consumption 231 can be input to the multiplication module 298 to determine the adjusted amount of ammonia consumption (i.e., adjusted ammonia consumption rate) 296.

Referring to FIG. 3, the nominal state 293 is selected or maintained when the first comparison unit 250 determines the Δ[NH3]_adsorbed 211 is at least the Δ[NH3]_desorbed 221. When the switch 295 is in the nominal state 293, a maintaining multiplier 1 output from the maintaining module 2 can be input to the multiplication module 298 along with the amount of ammonia consumption 231. It will be appreciated that the maintaining multiplier 1 has a value of “1.” Accordingly, the amount of ammonia consumption 231 is maintained and is not multiplied or adjusted by the adjustment multiplier 292.

The adjustment multiplier 292 has a value between zero and one. Accordingly, the adjusted amount of ammonia consumption 296 is less than the determined amount of ammonia consumption 231. Likewise, the rate of change of the adjusted ammonia consumption rate 296 is less than the determined ammonia consumption rate 231. In other words, applying the adjustment multiplier 292 to the determined ammonia consumption rate 231 decreases the rate of change of the determined ammonia consumption rate 231.

The adjustment module 290 can include look-up tables for determining the adjustment multiplier 292. In an exemplary embodiment, the look-up table is a one-dimensional table having an input corresponding to the reduced ratio 282 and an output corresponding to the determined adjustment multiplier 292. In the exemplary embodiment, the one-dimensional table utilizes an equation for slope. In a non-limiting example, the equation for slope can be expressed as follows.

$\begin{array}{cc}y=\frac{1}{2\sqrt{x}+1}& \left[16\right]\end{array}$

wherein x is an input including the reduced ratio 282, and

y is an output including the adjustment multiplier 292.

It will be appreciated that Equation [16] is merely a non-limiting example for determining the adjustment multiplier 292, and that any equation can be calibrated to determine the adjustment multiplier 292.

As aforementioned, a determined amount of ammonia consumption 231 without being adjusted when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed for each discrete substrate element 52(i) would result in an inaccurate ammonia storage concentration (θ_NH3) determination in each discrete substrate element 52(i) and an inaccurate total ammonia storage concentration on the ammonia-selective catalyst 60 (see block 140 of flowchart 100).

Accordingly, the states of parameters of the exhaust gas feedstream upstream of each discrete substrate element 52(i) (e.g., input gases 50 of FIG. 1) can be monitored and the ammonia storage concentration θ_NH3 for each discrete substrate element 52(i) can be determined based on the amount of ammonia that is adsorbed 211, the amount of ammonia that is desorbed 221 and the adjusted amount of ammonia consumption 296 when the amount of ammonia that is adsorbed 211 is less than the amount of ammonia that is desorbed 221. An amount of NOx at an outlet each discrete substrate element 52(i) can be estimated based on the determined ammonia storage concentration and the monitored states of parameters of the exhaust gas feedstream upstream of the ammonia-selective catalyst reduction device. Likewise, a total ammonia storage concentration on the ammonia-selective catalyst reduction device 60 can be determined based on the sequentially calculated ammonia adsorption rate 211, the ammonia desorption rate 221 and the adjusted ammonia consumption rate 296 when the ammonia adsorption rate 211 is less than the ammonia desorption rate 221; or the total ammonia storage concentration on the ammonia-selective catalyst reduction device 60 can be determined based on the sequentially calculated ammonia adsorption rate 211, the ammonia desorption rate 221 and the determined ammonia consumption rate 231 when the ammonia adsorption rate 211 is at least the ammonia desorption rate 221. An amount of NOx at an outlet of the ammonia-selective catalyst reduction device 60 can be estimated based on the total ammonia storage concentration on the ammonia-selective catalyst reduction device 60.

FIGS. 5 and 6 illustrate experimental and derived data of an exemplary ammonia-selective catalytic reduction (SCR) device, depicting an estimated amount of NOx at an outlet of the ammonia-SCR device and a measured amount of NOx at the outlet of the ammonia-SCR device over a period of time in accordance with the present disclosure. The horizontal axis in each of FIGS. 4 and 5 represents time in seconds. The vertical axis in each of FIGS. 4 and 5 represents NOx concentration in ppm.

With reference to FIG. 5, plot 400 illustrates a profile of an estimated amount of NOx 420 at the outlet of the ammonia-SCR device and a profile 410 of a measured amount of NOx 410 at the outlet of the ammonia-SCR device. At about 4500 seconds a condition occurs resulting in an amount of ammonia that is adsorbed on a catalyst surface being less than an amount of ammonia that is desorbed on the catalyst surface. For instance, ammonia or urea dosing may have been momentarily stopped or terminated. As a result, the profile of the estimated amount of NOx 420 sharply increases and deviates from the profile of the measured amount of NOx 410. This deviation is a result of an inaccurate amount of ammonia consumption determined when the amount of ammonia that is adsorbed on the catalyst surface is less than the amount of ammonia that is desorbed on the catalyst surface.

With reference to FIG. 6, plot 500 illustrates a profile of an estimated amount of NOx 520 at the outlet of the ammonia-SCR device and a profile 510 of a measured amount of NOx 410 at the outlet of the ammonia-SCR device. At about 4500 seconds a condition occurs resulting in an amount of ammonia that is adsorbed on a catalyst surface being less than an amount of ammonia that is desorbed on the catalyst surface. As evidenced by plot 500, when an amount of ammonia consumption is adjusted when the amount of ammonia that is adsorbed on the catalyst surface is less than the amount of ammonia that is desorbed on the catalyst surface, the profile of the estimated amount of NOx 520 closely follows the profile of the measured amount of NOx 510.

The disclosure has described certain preferred embodiments and modifications thereto. Further modifications and alterations may occur to others upon reading and understanding the specification. Therefore, it is intended that the disclosure not be limited to the particular embodiment(s) disclosed as the best mode contemplated for carrying out this disclosure, but that the disclosure will include all embodiments falling within the scope of the appended claims.

Claims

1. Method for monitoring a discrete substrate element from an ammonia-selective catalyst reduction device configured to treat an exhaust gas feedstream of an internal combustion engine, comprising: monitoring an amount of ammonia that is adsorbed, an amount of ammonia that is desorbed, an amount of ammonia that is oxidized and an amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream for the discrete substrate element;determining an amount of ammonia consumption for the discrete substrate element based on the amount of ammonia that is oxidized and the amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream;comparing the amount of ammonia that is adsorbed and the amount of ammonia that is desorbed for the discrete substrate element; andadjusting the amount of ammonia consumption for the discrete substrate element when the amount of ammonia that is adsorbed is less than the amount of ammonia that is desorbed.

2. The method of claim 1 wherein the discrete substrate element comprises one of a plurality of discrete substrate elements analytically segmented from the ammonia-selective catalytic reduction device.

3. The method of claim 1 further comprising: monitoring exhaust gas feedstream parameters upstream of the discrete substrate element;determining an ammonia storage concentration for the discrete substrate element based on the amount of ammonia that is adsorbed, the amount of ammonia that is desorbed and the adjusted amount of ammonia consumption; andestimating an amount of NOx at an outlet of the discrete substrate element based on the determined ammonia storage concentration and the exhaust gas feedstream parameters upstream of the ammonia-selective catalyst reduction device.

4. The method of claim 1 further comprising: maintaining the amount of ammonia consumption for the discrete substrate element when the amount of ammonia that is adsorbed is at least the amount of ammonia that is adsorbed.

5. The method of claim 1 wherein the amount of ammonia consumption is determined based on calculating a sum of the amount of ammonia that is oxidized and the amount of ammonia that is consumed in reducing NOx in the exhaust gas feedstream.

6. The method of claim 1 wherein adjusting the amount of ammonia consumption for the discrete substrate element is based on a ratio of the amount of ammonia consumption to the amount of ammonia that is adsorbed.

7. The method of claim 6, further comprising: when the monitored amount of ammonia that is adsorbed is equal to zero, designating the monitored amount of ammonia that is adsorbed to be equal to a predetermined non-zero integer.

8. The method of claim 6 wherein the ratio of the amount of ammonia consumption to the amount of ammonia that is adsorbed is proportionally reduced within a predetermined range.

9. The method of claim 6 wherein adjusting the amount of ammonia consumption for the discrete substrate element based on the ratio of the amount of ammonia consumption to the amount of ammonia that is adsorbed comprises: determining an adjustment multiplier based on the ratio of the amount of ammonia consumption to the amount of ammonia that is adsorbed; andapplying the adjustment multiplier to the amount of ammonia consumption for the discrete substrate element.

10. The method of claim 9 wherein the adjustment multiplier is a value between zero and one.

11. The method of claim 6 wherein the adjusted amount of ammonia consumption is less than the determined amount of ammonia consumption.

12. Method for modeling a consumption rate for an ammonia-selective catalyst reduction device configured to treat an exhaust gas feedstream of an internal combustion engine operating lean of stoichiometry, the method comprising: analytically segmenting the ammonia-selective catalytic catalyst reduction device into a plurality of discrete substrate elements;sequentially calculating a respective ammonia adsorption rate, a respective ammonia desorption rate, a respective ammonia oxidation rate, and a respective ammonia conversion rate in reducing NOx in the exhaust gas feedstream for each of the discrete substrate elements;determining for each of the discrete substrate elements a respective ammonia consumption rate based on a sum of the respective ammonia oxidation rate and the respective ammonia conversion rate;comparing for each of the discrete substrate elements the respective ammonia adsorption rate to the respective ammonia desorption rate; andbased on said comparison, for each of the discrete substrate elements, selecting one of the respective ammonia consumption rate when the respective ammonia adsorption rate is at least the respective ammonia desorption rate, anda respective adjusted ammonia consumption rate when the respective ammonia adsorption rate is less than the respective ammonia desorption rate.

13. The method of claim 12 wherein each respective adjusted ammonia consumption rate is based on a respective ratio of the respective ammonia consumption rate to the respective ammonia adsorption rate.

14. The method of claim 13 further comprising: when the respective ammonia adsorption rate is equal to zero, designating the respective ammonia adsorption rate to be equal to a predetermined non-zero integer.

15. The method of claim 13 wherein the respective ratio of the respective ammonia consumption rate to the respective ammonia adsorption rate is proportionally reduced within a predetermined range.

16. The method of claim 13 wherein each respective adjusted ammonia consumption rate is determined by: calculating a respective adjustment multiplier based on the ratio of the respective ammonia consumption rate to the respective ammonia adsorption rate; andapplying the respective adjustment multiplier to the respective ammonia consumption rate to determine the respective adjusted ammonia consumption rate.

17. The method of claim 16 wherein calculating each respective adjustment multiplier comprises: utilizing a one-dimensional table to calculate the adjustment multiplier based on the respective ratio.

18. The method of claim 16 wherein applying the respective adjustment multiplier to the respective ammonia consumption rate decreases the rate of change of the respective ammonia consumption rate.

19. The method of claim 12 further comprising: determining a total ammonia storage concentration on the ammonia-selective catalyst reduction device based on the sequentially calculated ammonia adsorption rates, the ammonia adsorption rates and the selected one of the respective ammonia consumption rate and the respective adjusted ammonia consumption rate for all the discrete substrate elements; andestimating an amount of NOx at an outlet of the ammonia-selective catalyst reduction device based on the total ammonia storage concentration on the ammonia-selective catalyst reduction device.