# Patent

## Title
Catalyst temperature estimation device

## Abstract

In a fuel cut period, an adsorbed CO quantity calculation section calculates the amount CATco of CO adsorbed on a catalyst, and an exhaust O2 quantity calculation section calculates the amount GASo2 of O2 in exhaust gases. CO oxidation reaction heat ΔHr produced when the calculated amount CATco of adsorbed CO reacts with the calculated amount GASo2 of O2 in the exhaust gases is calculated by a ΔHr calculation section, and used in estimation of the temperature of the catalyst.

## Background
BACKGROUND OF THE INVENTION
1. Field of the Invention
This invention relates to a catalyst temperature estimation device for estimating the temperature of a catalyst of a catalytic converter for converting exhaust to less harmful substances.
2. Description of the Related Art
A catalytic converter for converting harmful substances in exhaust gases, such as HC (carbon hydrate), CO (carbon monoxide) and NOx (nitrogen oxide), to less harmful substances is provided in an exhaust system of an engine. The temperature of the catalyst of such catalytic converter varies to a great degree depending on the operating state of the engine. Thus, it is necessary to control the engine appropriately, on the basis of the estimated value of the catalyst temperature, so that the catalyst temperature will not exceed the allowable temperature limit for the catalytic converter in any operating state. There is known a catalyst temperature estimation technique in which exhaust temperature in steady operation is stored, and exhaust temperature in transient operation is estimated by subjecting the stored temperature in steady operation to filtering.
When a fuel cut is performed for vehicle deceleration, a large amount of oxidation reaction heat is produced by a reaction in which O2 (oxygen) in the exhaust gases reacts with CO and HC adsorbed on the catalyst. Meanwhile, in a fuel recovery period after the termination of the fuel cut, a large amount of conversion reaction heat is produced by a reaction in which CO and HC in the exhaust gases are converted by O2 that was adsorbed onto the catalyst during the fuel cut. In either case, a rapid rise in catalyst temperature is caused, which makes it particularly difficult to estimate the catalyst temperature. The above-mentioned catalyst temperature estimation technique only estimates the rate of change of temperature in such transient operation to be smaller than that in steady operation, through filtering. Thus, especially when the fuel cut and fuel recovery is repeated frequently, there occurs a problem that the values departing much from the actual catalyst temperature are estimated, as indicated in dashed line in the time chart of FIG. 7.
A catalyst temperature estimation technique different from the above has also been proposed (see Japanese Unexamined Patent Publication No. 2004-263606, for example). In the technique disclosed in this publication, in the operating states other than that with a fuel cut, the temperature of exhaust from the engine is considered as standard catalyst temperature; from the standard catalyst temperature, exhaust pipe wall temperature is estimated; From the standard catalyst temperature and the exhaust pipe wall temperature, provisional catalytic-converter inflow exhaust temperature is estimated; then, catalytic-converter inflow exhaust temperature taking account of a drop in exhaust temperature due to running wind is estimated from the provisional catalytic-converter inflow exhaust temperature; and on the basis of this catalytic-converter inflow exhaust temperature, catalyst temperature is estimated. Meanwhile, in the operation with a fuel cut, this estimation process is suspended, and the value of the catalyst temperature estimated in the last processing cycle is used as the present value without change.
However, the use of the last value as the present value without change in the fuel cut period, as in the catalyst temperature estimation technique disclosed in the above-mentioned publication, means that the estimation is performed on the assumption that the catalyst temperature does not vary in the fuel cut period and the fuel recovery period. This results in great departure of estimated values from the actual catalyst temperature which shows a rapid rise due to the oxidation reaction of HC and CO, like the case of the conventional technique shown in FIG. 7, for example. Consequently, it is impossible to suppress a rise in catalyst temperature in the fuel cut period and the fuel recovery period, by performing engine control on the basis of the catalyst temperature estimated by this estimation technique. Accordingly, this estimation technique has a problem that the catalyst temperature rises beyond the allowable temperature limit, which causes damage to the catalytic converter.
SUMMARY OF THE INVENTION
This invention has been made to solve the problems as mentioned above. The primary object of this invention is to provide a catalyst temperature estimation device which can accurately estimate movements of the catalyst temperature also in the fuel cut period and the fuel recovery period of the engine, and which can therefore prevent the catalytic converter from being damaged by a rise in temperature beyond the allowable temperature limit.
To achieve this object, the catalyst temperature estimation device according to the present invention comprises a catalytic converter carrying a catalyst disposed in an exhaust passage of an engine for converting exhaust gases to less harmful substances; an adsorbed unburnt-fuel-component quantity calculation element for calculating the amount of an unburnt fuel component adsorbed on the catalyst in a fuel cut period of the engine; an exhaust O2 quantity calculation element for calculating the amount of O2 in the exhaust gases in the fuel cut period; an oxidation reaction heat calculation element for calculating reaction heat produced when the adsorbed unburnt fuel component in the amount calculated by the adsorbed unburnt-fuel-component quantity calculation element reacts with O2 in the exhaust gases in the amount calculated by the exhaust O2 quantity calculation element; and a catalyst temperature estimation element for estimating the temperature of the catalyst in the fuel cut period, on the basis of the reaction heat calculated by the oxidation reaction heat calculation element.
Thus, in the fuel cut period, the amount of the unburnt fuel component adsorbed on the catalyst is calculated by the adsorbed unburnt-fuel-component quantity calculation element, the amount of O2 in the exhaust gases is calculated by the exhaust O2 quantity calculation element, the reaction heat produced when the unburnt fuel component in the calculated amount reacts with O2 in the calculated amount is calculated by the oxidation reaction heat calculation element, and the temperature of the catalyst is estimated by the catalyst temperature estimation element on the basis of the calculated reaction heat.
In the fuel cut period, a large amount of reaction heat is produced by an oxidation reaction in which unburnt fuel components such as CO and HC adsorbed on the catalyst is oxidized by a large amount of O2 in the exhaust gases, and this reaction heat causes a rapid rise in catalyst temperature. With the configuration described above, the catalyst temperature estimation device according to the present invention can estimate the catalyst temperature taking account of the reaction heat in such phenomenon. Accordingly, it can estimate movements of the catalyst temperature in the fuel cut period, accurately, and therefore prevent the catalytic converter from being damaged by a rise in temperature beyond the allowable temperature limit.
The catalyst temperature estimation device according to the present invention can further comprise an operating state detection element for detecting the operating state of the engine; an O2 storage index calculation element for calculating an O2 storage index correlating with the amount of O2 adsorbed on the catalyst; an air-fuel ratio estimation element for estimating air-fuel ratio on the catalyst; and a conversion reaction heat calculation element for calculating reaction heat produced by a conversion reaction on the catalyst in a fuel recovery period of the engine, on the basis of the operating state of the engine detected by the operating state detection element, the O2 storage index calculated by the O2 storage index calculation element, and the air-fuel ratio estimated by the air-fuel ratio estimation element, wherein the catalyst temperature estimation element estimates the temperature of the catalyst in the fuel cut period on the basis of the reaction heat calculated by the oxidation reaction heat calculation element, and in the fuel recovery period on the basis of the reaction heat calculated by the conversion reaction heat calculation element.
Thus, in the fuel recovery period, the operating state of the engine is detected by the operating state detection element, the O2 storage index correlating with the amount of O2 adsorbed on the catalyst is calculated by the O2 storage index calculation element, and the air-fuel ratio on the catalyst is estimated by the air-fuel ratio estimation element. On the basis of these engine operating state, O2 storage index and air-fuel ratio, reaction heat produced by a conversion reaction on the catalyst is calculated by the conversion reaction heat calculation element, and on the basis of the calculated reaction heat, the temperature of the catalyst is estimated by the catalyst temperature estimation element.
In the fuel recovery period, a large amount of reaction heat is produced by a conversion reaction in which CO and HC in the exhaust gases are converted by O2 that was adsorbed onto the catalyst during the fuel cut, and this reaction heat causes a rapid rise in catalyst temperature. With the configuration described above, the catalyst temperature estimation device according to the present invention can estimate the catalyst temperature taking account of the reaction heat in such phenomenon. Accordingly, it can estimate movements of the catalyst temperature accurately not only in the fuel cut period but also in the fuel recovery period.
A further scope of applicability of the present invention will become apparent from the detailed description given hereinafter. However, it should be understood that the detailed description and specific example, while indicating preferred embodiments of the invention, are given by way of illustration only, since various changes and modifications within the spirit and scope of the invention will become apparent to those skilled in the art from this detailed description.



BRIEF DESCRIPTION OF THE DRAWINGS
The nature of this invention, as well as other objects and advantages thereof, will be explained in the following with reference to the accompanying drawings, in which like reference characters designate the same or similar parts throughout the figures and wherein:
FIG. 1 is a diagram showing the entire configuration of an embodiment of the catalyst temperature estimation device;
FIG. 2 is a control block diagram showing the processing performed by a fuel-cut-period reaction heat calculation section of an ECU;
FIG. 3 is a control block diagram showing the processing performed by an exhaust-catalyst heat transfer calculation section of the ECU;
FIG. 4 is a control block diagram showing the processing performed by a fuel-recovery-period conversion heat calculation section of the ECU;
FIG. 5 is a control block diagram showing the processing performed by an O2 adsorption ratio calculation section of the ECU;
FIG. 6 is a control block diagram showing how the processing by a catalyst temperature estimation section of the ECU goes on; and
FIG. 7 shows a test result showing how the catalyst temperature is estimated when the fuel cut and fuel recovery is repeated with a short period.



DETAILED DESCRIPTION OF THE INVENTION
A catalyst temperature estimation device which is an embodiment of the present invention will be describe below.
FIG. 1 is a diagram showing the entire configuration of this embodiment of the catalyst temperature estimation device. An engine 1 is a direct-injection spark-ignition four-cylinder-in-line gasoline engine. A spark plug 2 and an electromagnetic fuel injection valve 3, provided for each cylinder, are fitted to a cylinder head of the engine 1. From the fuel injection valve 3, fuel is injected directly into a combustion chamber. The cylinder head has intake ports 4a formed for the respective cylinders, which extend almost upright. The intake ports 4a are connected to a throttle valve 5 through an intake manifold 4, and the throttle valve 5 is connected to an intake passage (not shown).
The cylinder head also has exhaust ports 6a extending almost horizontally, and the exhaust ports 6a are connected to an exhaust passage 7 through an exhaust manifold 6. An early catalytic converter 8 is disposed in the upstream part of the exhaust passage 7, while an underfloor catalytic converter 9 is disposed in the downstream part of the exhaust passage 7. The early catalytic converter 8 and the underfloor catalytic converter 9 are three-way catalytic converters carrying a precious metal such as platinum (Pt), rhodium (Rh) or the like. In the present embodiment, temperature estimation is performed for the underfloor catalytic converter 9. It is, however, to be noted that the temperature estimation is not limited to the underfloor catalytic converter 9 but can be carried out for the early catalytic converter 8 or an NOx adsorption catalytic converter for conversion of NOx (not shown), for example. Upstream of the underfloor catalytic converter 9, an air-fuel ratio sensor 10 (exhaust air-fuel ratio detection element) is provided. The air-fuel ratio sensor 10 detects the air-fuel ratio of exhaust gases flowing into the underfloor catalytic converter 9.
In a vehicle compartment, an ECU (electronic control unit) 21 including an input-output device, memory (ROM, RAM, nonvolatile RAM, etc.), a central processing unit (CPU), a timer counter, etc. is installed. To the input of the ECU 21, various sensors, such as the above-mentioned air-fuel ratio sensor 10, a throttle sensor 22 for detecting the opening of the throttle valve 5, a revolving rate sensor 23 for detecting the revolving rate of the engine 1, and a water temperature sensor 24 for detecting the temperature of a coolant, are connected. From these sensors, detection information is fed to the ECU. To the output of the ECU 21, various devices, such as the above-mentioned spark plugs 2 and fuel injection valves 3, are connected.
The ECU 21 determines ignition timing and fuel injection timing on the basis of the detection information from the sensors, and makes the engine 1 operate by drive-controlling the spark plugs 2 and fuel injection valves 3 on the basis of determined control quantities. When the conditions set for fuel cut, such that the vehicle is being decelerated with an accelerator deactivated and that the vehicle speed is greater than or equal to a determined value, are satisfied, the ECU 21 stops fuel injection by performing a fuel cut, and when the conditions set for fuel recovery, such that the accelerator is activated, are satisfied, the ECU 21 resumes fuel injection.
As stated under the subheading “Description of the Related Art”, in the fuel cut period and the fuel recovery period, the catalyst temperature Tcatf is difficult to estimate due to a rapid rise in temperature in the underfloor catalytic converter 9, which necessarily makes it difficult to appropriately perform engine control on the basis of the catalyst temperature Tcatf in order to suppress a rise in temperature. Thus, in the present embodiment, temperature estimation for the underfloor catalytic converter 9 is performed taking account of reaction heat produced in the fuel cut period and the fuel recovery period. This temperature estimation will be described below in detail.
Prior to entering the description of the temperature estimation, reactions taking place on the catalyst in the fuel cut period and the fuel recovery period and the influences of those reactions on the catalyst temperature Tcatf will be described.
As in the other operating states, also in the fuel cut period and the fuel recovery period, the rise in catalyst temperature basically depends on heat transfer ΔHt between the exhaust gases and the catalyst 9. Thus, the estimation of the catalyst temperature Tcarf needs to take account of the heat transfer ΔHt. In the fuel cut period, in addition to the rise in catalyst temperature due to the heat transfer ΔHt, the catalyst temperature Tcatf rises due to reaction heat produced when CO and HC that were adsorbed onto the catalyst 9 before the fuel cut (both are unburnt fuel components; Note that in the description below, to represent both of these CO and HC, the term “CO” will be used representatively) are oxidized by O2 in the exhaust gases. Thus, the estimation of the catalyst temperature Tcarf needs to take account of this oxidation reaction heat ΔHr.
Meanwhile, in the fuel recovery period, the catalyst temperature Tcatf rises due to reaction heat produced when CO and HC in the exhaust gases (both are unburnt fuel components; Note that in the description below, to represent both of these CO and HC, the term “CO” will be used representatively) are converted by O2 that was adsorbed onto the catalyst 9 during the fuel cut. Thus, the estimation of the catalyst temperature Tcarf needs to take account of this conversion reaction heat ΔHc. Accordingly, as will be described below, in the fuel cut period, the catalyst temperature Tcatf is estimated taking account of the oxidation reaction heat ΔHr in addition to the heat transfer ΔHt between the exhaust gases and the catalyst 9, and in the fuel recovery period, the catalyst temperature Tcatf is estimated taking account of the conversion reaction heat ΔHc in addition to the heat transfer ΔHt.
First, the calculation of the CO oxidation reaction heat ΔHr in the fuel cut period will be described.
FIG. 2 is a control block diagram showing the processing performed by a fuel-cut-period reaction heat calculation section 50 of the ECU 21. In the fuel cut period, a phenomenon that CO is desorbed from the catalyst through reaction with O2 in the exhaust gases and a phenomenon that CO is adsorbed onto the catalyst take place in parallel. Both phenomena together affect the CO adsorption ratio θco of the catalyst 9, and therefore the oxidation reaction heat ΔHr in the fuel cut period. Thus, in the calculation of the CO oxidation reaction heat ΔHr, the adsorbed CO quantity CATco, namely the amount of CO adsorbed on the catalyst 9 is calculated taking account of the desorption of CO from the catalyst and the adsorption of CO onto the catalyst, and the CO oxidation reaction heat ΔHr is calculated taking account of this adsorbed CO quantity CATco.
In FIG. 2, a desorption constant calculation section 51 of the fuel-cut-period reaction heat calculation section 50 obtains a desorption rate constant Kred, referring to a map prepared in advance, on the basis of the catalyst temperature Tcatf estimated in the last processing cycle. The desorption rate Rred at which CO is desorbed from the catalyst varies depending on the activated state of the catalyst 9. Thus, the desorption rate constant Kred is determined depending on the catalyst temperature Tcatf which correlates with the activated state.
An adsorption constant calculation section 52 obtains an adsorption rate constant Kad, referring to a map prepared in advance, on the basis of the catalyst temperature Tcatf and the exhaust air-fuel ratio A/F detected by the air-fuel ratio sensor 10. The adsorption rate Rad at which CO is adsorbed onto the catalyst 9 is affected by not only the activated state of the catalyst 9 but also the composition of the exhaust gases which varies depending on the exhaust air-fuel ratio A/F. Thus, the adsorption rate constant Kad is determined depending on the catalyst temperature Tcatf and the exhaust air-fuel ratio A/F. It is to be noted that the exhaust air-fuel ratio A/F may be estimated from the engine operating state, etc., according to a known estimation method.
An adsorption sites calculation section 53 obtains the number ρco (mol) of all the adsorption sites of the catalyst 9 capable of adsorbing CO, referring to a map prepared in advance, on the basis of the catalyst temperature Tcatf. The number ρco of the adsorption sites is proper to the catalyst but varies depending on the activated state of the catalyst 9. Thus, the number ρco of the adsorption sites is determined depending on the catalyst temperature Tcatf which correlates with the activated state.
A CO concentration calculation section 54 determines CO concentration Pco (atm) in the exhaust gases, referring to a map prepared in advance, on the basis of the exhaust air-fuel ratio A/F. Since the CO concentration Pco correlates with the exhaust air-fuel ratio A/F, the CO concentration Pco is determined depending on the exhaust air-fuel ratio A/F.
Meanwhile, in the calculation of the oxidation reaction heat ΔHr in the last processing cycle, the CO adsorption ratio θco of the catalyst 9 has been calculated, and between the CO adsorption ratio θco and the number ρco of the adsorption sites, there exists a relationship expressed by the equation (1) below:

θco=rρco/ρco  (1)

where rρco is the number of the adsorption sites of the catalyst 9 at which CO is really adsorbed.

The CO adsorption ratio θco, the desorption rate constant Kred calculated by the desorption constant calculation section 51, and the number ρco of the adsorption sites calculated by the adsorption sites calculation section 53 are fed to a desorption rate calculation section 55. Using these input values, the desorption rate calculation section 55 calculates the CO desorption rate Rred (mol/sec) at which CO is desorbed from the catalyst, according to the equation (2) below (desorption rate calculation element):

Rred=Kred·ρco·θco  (2)

The CO adsorption ratio θco is also fed to a subtraction section 56, and the subtraction section 56 subtracts the CO adsorption ratio θco from 1 (1−θco). The value obtained by the subtraction section 56, the adsorption rate constant Kad calculated by the adsorption rate constant calculation section 52, the number ρco of the adsorption sites calculated by the adsorption sites calculation section 53, and the CO concentration Pco calculated by the CO concentration calculation section 54 are fed to an adsorption rate calculation section 57, and using these input values, the adsorption rate calculation section 57 calculates the CO adsorption rate Rad (mol/sec) at which CO is adsorbed onto the catalyst, according to the equation (3) below (adsorption rate calculation element):

Rad=Kad·ρco·(1−θco)·Pco  (3)

Further, the CO concentration Pco calculated by the CO concentration calculation section 54, total pressure 1 (atm), and the total substance quantity per unit time nall (mol/sec) of the exhaust gases, namely the amount of all the substances in the exhaust gases in unit time are fed to a partial pressure calculation section 58, and the partial pressure calculation section 58 calculates the exhaust CO quantity GASco (mol/sec), namely the amount of CO in the exhaust gases, in terms of partial pressure, according to the equation (4) below. It is to be noted that the total substance quantity per unit time nall of the exhaust gases is calculated according to a known calculation method, from the intake quantity Q, the amount of air and of fuel in terms of mole, the exhaust air-fuel ratio A/F, etc., of which the detailed explanation will be omitted.

GASco=Pco/1·nall  (4)

The CO adsorption rate Rad on the catalyst and the exhaust CO quantity GASco are fed to a minimum-value selection section 59, and the minimum-value selection section 59 chooses a smaller one of these input values, thereby determining the definitive CO adsorption rate Rad on the catalyst. Specifically, the CO adsorption rate Rad represents the capacity of the catalyst 9 to adsorb CO. When the exhaust CO quantity GASco is greater than or equal to the CO adsorption rate Rad, CO is actually adsorbed onto the catalyst at the CO adsorption rate Rad. When the exhaust CO quantity GASco is less than the CO adsorption rate Rad, the rate at which CO is actually adsorbed onto the catalyst is limited to the exhaust CO quantity GASco. For this reason, the minimum-value selection section 59 chooses a smaller one of the input values.
The CO desorption rate Rred calculated by the desorption rate calculation section 55 and the CO adsorption rate Rad on the catalyst definitively determined by the minimum-value selection section 59 are fed to an adsorbed-CO-quantity-per-unit-time calculation section 60, and the adsorbed-CO-quantity-per-unit-time calculation section 60 calculates the adsorbed CO quantity per unit time Δadθco, namely the amount of CO adsorbed onto the catalyst 9 in unit time, according to the equation (5) below:

Δadθco=Rad−Rred  (5)

The adsorbed CO quantity per unit time Δadθco thus calculated, the number ρco of the adsorption sites of the catalyst 9, and the processing period f of the ECU 21 (0.1 msec, for example) are fed to an adsorption-ratio change calculation section 61, and the adsorption-ratio change calculation section 61 calculates a change Δθco in CO adsorption ratio θco in the processing period f, according to the equation (6) below:

Δθco=Δadθco/ρco·f  (6)

The change Δθco in CO adsorption ratio θco and the CO adsorption ratio θco are fed to an adsorption ratio calculation section 62, and the adsorption ratio calculation section 62 takes this CO adsorption ratio θco as a value θco(n−1) in the last processing cycle, and calculates the present CO adsorption ratio θco(n), according to the equation (7) below:

θco(n)=θco(n−1)+Δθco  (7)

Then, the CO adsorption ratio θco(n), the number ρco of the adsorption sites, and the processing period f of the ECU 21 are fed to an adsorbed CO quantity calculation section 63, and the adsorbed CO quantity calculation section 63 calculates the adsorbed CO quantity per unit time CATco (mol/sec) corresponding to the amount of CO adsorbed on the catalyst 9 at present, according to the equation (9) below (adsorbed unburnt-fuel-component quantity calculation element):

CATco=θco(n)·ρco/f  (9)

Meanwhile, an O2 concentration calculation section 64 obtains O2 concentration Po2, referring to a map prepared in advance, on the basis of the exhaust air-fuel ratio A/F. Like the above-mentioned CO concentration Pco, the O2 concentration Po2 is determined depending on the exhaust air-fuel ratio A/F, but, for the convenience of substitution in the later-mentioned equation (10) representing the reaction rate, the map is prepared to give the square root Po21/2 of the O2 concentration Po2.
A reaction rate constant calculation section 65 obtains a reaction rate constant Kr for the reaction of CO and O2, referring to a map prepared in advance, on the basis of the catalyst temperature Tcatf. Since the rate of reaction of CO and O2 varies depending on the activated state of the catalyst 9, the reaction rate constant Kr is determined depending on the catalyst temperature Tcatf.
The O2 concentration Po21/2 calculated by the O2 concentration calculation section 64, the reaction rate constant Kr calculated by the reaction rate constant calculation section 65, the CO adsorption ratio θco and the number ρco of the adsorption sites are fed to a reaction rate calculation section 66, and the reaction rate calculation section 66 calculates the maximum reaction rate r (mol/sec) at which CO and O2 can react, according to the equation (10) below (reaction rate calculation element):

r=Kr·ρco·θco(Po2)1/2  (10)

Specifically, since the maximum reaction rate r at which CO and O2 react on the catalyst is determined depending on the product of the CO quantity CATco (=ρco·θco) and the O2 concentration Po2, and since the mole ratio of CO to O2 in the reaction is “2” as seen from the equation (11) below, the above equation (10) includes the O2 concentration in the form of square root PO21/2.

CO+1/2O2→CO2  (11)

The reaction mole ratio “2” determined from the equation (11), the O2 concentration Po21/2 calculated by the O2 concentration calculation section 64, total pressure 1 (atm), and the total substance quantity per unit time nall of the exhaust gases are fed to an exhaust O2 quantity calculation section 67, and the exhaust O2 quantity calculation section 67 calculates the exhaust O2 quantity GASo2 (mol/sec), namely the amount of O2 in the exhaust gases, according to the equation (12) below (exhaust O2 quantity calculation element).

GASo2=2·Po2/1·nall  (12)

The adsorbed CO quantity CATco on the catalyst calculated by the adsorbed CO quantity calculation section 63, the reaction rate r calculated by the reaction rate calculation section 66, and the exhaust O2 quantity GASo2 calculated by the exhaust O2 quantity calculation section 67 are fed to a minimum-value selection section 68, and the minimum-value selection section 68 chooses the smallest one of these input values, thereby determining the reacting CO quantity Rct, namely the amount of CO taking part in the reaction. While the quantity of CO reacting with O2 is basically determined depending on a smaller one of the adsorbed CO quantity CATco and the exhaust O2 quantity GASo2, the reaction at a rate higher than the reaction rate r is impossible. Thus, the smallest one of the three quantities including the reaction rate r is chosen as the reacting CO quantity Rct.
The reacting CO quantity Rct calculated by the minimum-value selection section 68 and the value of a fuel-cut flag, which is at a ON-value while the fuel cut is being performed on the engine 1, are fed to a fuel-cut determination section 69. When determining that the fuel cut is not being performed from the fuel-cut flag reset to an OFF-value, the fuel-cut determination section 69 chooses “0”, and when determining that the fuel cut is being performed from the fuel-cut flag set to the ON-value, chooses the reacting CO quantity Rct. The value chosen by the fuel-cut determination section 69, reaction heat produced by fuel (283 kJ/mol, for example), and the processing period f of the ECU 21 are fed to a ΔHr calculation section 70, and using these input values, the ΔHr calculation section 70 calculates the CO oxidation reaction heat ΔHr produced by the CO oxidation reaction on the catalyst (oxidation reaction heat calculation element). Thus, while the fuel cut is not being performed, the fuel-cut-period reaction heat calculation section 50 sends out “0”, and while the fuel cut is being performed, sends out the CO oxidation reaction heat ΔHr.
The value (Rct or “0”) chosen by the fuel-cut determination section 69 and the number ρco of the adsorption sites are fed to an adsorption-ratio decrease calculation section 71, and using these input values, the adsorption-ratio decrease calculation section 71 calculates a decrease Δrθco in adsorption ratio of the catalyst caused by the reaction of CO. The decrease Δrθco in adsorption ratio thus calculated and the CO adsorption ratio θco(n) calculated by the adsorption ratio calculation section 62 in the present processing cycle are fed to an adsorption ratio update section 72, and the adsorption ratio update section 72 calculates the latest CO adsorption ratio θco, according to the equation (13) below:

θco=θco(n)−Δrθco  (13)

The CO adsorption ratio θco calculated by the adsorption ratio update section 72 is limited by an upper clipping section 73 to the upper limit 1.0, and limited by a lower clipping section 74 to the lower limit 0. The processing by the clipping sections 73 and 74 is provided for an improper CO adsorption ratio θco which may possibly be produced for some reason. The CO adsorption ratio θco after this processing is used in the next processing cycle of the ECU 21.
Next, the calculation of heat transfer ΔHt (kJ/sec) between the exhaust gases and the catalyst 9 will be described.
FIG. 3 is a control block diagram showing the processing performed by an exhaust-catalyst heat transfer calculation section 80 of the ECU 21. The heat transfer ΔHt between the exhaust gases and the catalyst 9 depends on temperature difference between the exhaust gases and the catalyst 9, contact area between the exhaust gases and the catalyst 9, and a coefficient of heat transfer. Between these quantities, there exists a relationship expressed by the equation (14) below:

ΔHt=h·Sv·(Tex−Tcatf)·Vcat  (14)

where Sv is the specific surface area of the catalyst 9, h is a coefficient of heat transfer, Tex is the temperature of the exhaust gases obtained by estimation, Tcatf is the temperature of the catalyst, and Vcat is the volume of the catalyst. It is to be noted that although here, the exhaust gas temperature Tex is assumed to be obtained by estimation, it may be detected using an exhaust temperature sensor 11 as shown in FIG. 1.

In FIG. 3, a contact area calculation section 81 of the heat transfer calculation section 80 calculates the effective contact area Scat of the catalyst 9, referring to a map prepared in advance, on the basis of the intake quantity Q. The effective contact area Scat is the area of the catalyst 9 subjected to contact with the exhaust gases, and corresponds to the product of the specific surface area Sv and the volume Vcat of the catalyst 9 in the equation (14).
A heat transfer coefficient calculation section 82 calculates a heat transfer coefficient h, referring to a map prepared in advance, on the basis of the intake quantity Q and the exhaust gas temperature Tex. A temperature difference calculation section 83 calculates a difference ΔT between the exhaust gas temperature Tex and provisional catalyst temperature Tcat (catalyst temperature before subjected to filtering).
The effective contact area Scat, the heat transfer coefficient h, the temperature difference ΔT, and the processing period f of the ECU 21 are fed to a ΔHt calculation section 84, and using these quantities together, the ΔHt calculation section 84 calculates the heat transfer ΔHt (kJ/sec), according to the equation (14) (heat transfer calculation element).
Next, the calculation of the CO conversion reaction heat ΔHc, namely the heat produced by the reaction in which CO is converted by O2 adsorbed on the catalyst 9, in the fuel recovery period.
FIG. 4 is a control block diagram showing the processing performed by a fuel-recovery-period conversion heat calculation section 110 of the ECU 21, and FIG. 5 is a control block diagram showing the processing performed by an O2 adsorption ratio calculation section 90 of the ECU 21. The CO conversion reaction heat ΔHc, namely the heat produced by the conversion of CO on the catalyst depends on the operating state of the engine 1, specifically, engine revolving rate Ne, charging efficiency Ec, and air-fuel ratio on the catalyst, and can be determined from these quantities. For the air-fuel ratio on the catalyst, basically the exhaust air-fuel ratio A/F can be used, while the air-fuel ratio on the catalyst is affected by the O2 storage function of the catalyst 9, which needs consideration. Specifically, the catalyst 9 has a function of maintaining the exhaust air-fuel ratio A/F at the stoichiometric air-fuel ratio (=14.7), by adsorbing O2 at lean exhaust air-fuel ratios and desorbing O2 at rich exhaust air-fuel ratios. Thus, when the catalyst is exhibiting this O2 storage function, the air-fuel ratio on the catalyst is maintained at the stoichiometric, although the exhaust air-fuel ratio A/F is not the stoichiometric. Meanwhile, when the catalyst is not exhibiting the O2 storage function, the air-fuel ratio on the catalyst equals the exhaust air-fuel ratio A/F.
Thus, in the calculation of the CO conversion reaction heat ΔHc, the O2 adsorption ratio θo2 which correlates with the O2 storage function of the catalyst 9 is calculated, and the calculated O2 adsorption ratio θo2 is used to calculate the CO conversion reaction heat ΔHc (kJ/sec). First, the calculation of the O2 adsorption ratio θo2 will be described.
Like the behavior of CO on the catalyst in the fuel cut period, in the fuel recovery period, a phenomenon that O2 is desorbed from the catalyst through reaction with CO in the exhaust gases and a phenomenon that O2 in the exhaust gases is adsorbed onto the catalyst take place in parallel. Both phenomena together affect the O2 adsorption ratio θo2 of the catalyst 9, and therefore the CO conversion reaction heat ΔHc in the fuel recovery period. Thus, in the calculation of the CO conversion reaction heat ΔHc, the adsorbed O2 quantity, namely the amount of O2 adsorbed on the catalyst 9 is calculated taking account of the desorption of O2 from the catalyst and the adsorption of O2 onto the catalyst, and the CO conversion reaction heat ΔHc is calculated taking account of this adsorbed O2 quantity.
In FIG. 5, an adsorption sites calculation section 91 of an O2 adsorption ratio calculation section 90 obtains the number ρo2 (mol) of all the adsorption sites of the catalyst 9 capable of adsorbing O2, referring to a map prepared in advance, on the basis of the catalyst temperature Tcatf. As in the case of the CO adsorption sites, the number ρo2 of the adsorption sites is proper to the catalyst but varies depending on the activated state of the catalyst 9. Thus, the number ρo2 of the adsorption sites is determined depending on the catalyst temperature Tcatf which correlates with the activated state. Between the number ρo2 of the adsorption sites and the O2 adsorption ratio θo2, there exists a relationship expressed by the equation (15) below:

θo2=rρo2/ρo2  (15)

where rρo2 is the number of the adsorption sites of the catalyst 9 at which O2 is really adsorbed.

The number ρo2 of the adsorption sites calculated by the adsorption sites calculation section 91 and the O2 adsorption ratio θo2 in the last processing cycle are fed to an O2 desorption capacity calculation section 92, and using these input values, the O2 desorption capacity calculation section 92 calculates the present O2 desorption capacity Cred (mol) of the catalyst 9, according to the equation (16) below:

Cred=ρo2·θo2·(−1)  (16)

It is to be noted that in the equation (16), the O2 desorption capacity Cred is obtained as a negative value by multiplying −1, considering that the desorption of O2 decreases the O2 adsorption ratio θo2.
A subtraction section 93 subtracts the O2 adsorption ratio θo2 from 1 (1−θo2). The value obtained by the subtraction and the number ρo2 of the adsorption sites calculated by the adsorption sites calculation section 91 are fed to an O2 adsorption capacity calculation section 94, and using these input values, the O2 adsorption capacity calculation section 94 calculates the present O2 adsorption capacity Cad(mol) of the catalyst 9, according to the equation (17) below:

Cad=ρo2·(1−θo2)  (17)

Meanwhile, an O2 excess/deficiency calculation section 95 obtains O2 excess/deficiency ΔO2 (vol %) in the exhaust gases, referring to a map prepared in advance, on the basis of the exhaust air-fuel ratio A/F. The O2 excess/deficiency ΔO2 means the amount by which O2 exceeds or falls below the stoichiometric quantity, and is determined depending on the exhaust air-fuel ratio A/F. The O2 excess/deficiency ΔO2 thus obtained and the total substance quantity per unit time nail (mol/sec) of the exhaust gases are fed to a conversion section 96, and the conversion section 96 converts the O2 excess/deficiency ΔO2 to the quantity per unit time (mol/sec).
To an O2 expenditure calculation section 97, the reacting CO quantity Rct (mol/sec) from the fuel-cut determination section 69 of the fuel-cut-period reaction heat calculation section 50 and the mole ratio “2” between CO and O2 in the reaction are fed, and using these input values, the O2 expenditure calculation section 97 calculates the O2 expenditure expΔO2 in the exhaust gases, namely the amount of O2 in the exhaust gases which was consumed in the reaction with CO on the catalyst during the fuel cut, according to the equation (18) below:

expΔO2=Rct/2  (18)

To a subtraction section 98, the O2 excess/deficiency ΔO2 after the conversion by the conversion section 96, and the O2 expenditure expΔO2 calculated by the O2 expenditure calculation section 97 are fed, and the subtraction section 98 calculates the stoichometric required O2 quantity Δado2, which means the amount per unit time of O2 required to make the composition of the exhaust gases stoichiometric (or in other words, to completely convert CO in the exhaust gases), according to the equation (19) below. Specifically, when the exhaust air-fuel ratio A/F is rich so that CO in the exhaust gases needs to be oxidized by O2, the stoichometric required O2 quantity Δado2 to be calculated is the amount of O2 which needs to be desorbed from the catalyst (negative value), and when the exhaust air-fuel ratio A/F is lean so that excess O2 needs to be adsorbed, the stoichometric required O2 quantity Δado2 to be calculated is the amount of O2 which needs to be adsorbed onto the catalyst (positive value). The stoichometric required O2 quantity Δado2 thus calculated is fed to a conversion section 99 and converted to the quantity Δadθo2 in the processing period f, according to the equation (20) below.

Δado2=ΔO2−expΔO2  (19)

Δadθo2=Δado2·f  (20)

The stoichometric required O2 quantity Δadθo2 after the conversion, the O2 desorption capacity Cred and the O2 adsorption capacity Cad of the catalyst 9, and the exhaust air-fuel ratio A/F are fed to an O2 desorption/adsorption quantity calculation section 100. When the exhaust air-fuel ratio is rich, namely richer than the stoichiometric air-fuel ratio, so that it is inferred that O2 is being desorbed from the catalyst, the O2 desorption/adsorption quantity calculation section 100 chooses a greater quantity between the O2 desorption capacity Cred and the stoichometric required O2 quantity Δadθo2 (namely, the quantity having a smaller absolute value, since both quantities are negative). When the exhaust air-fuel ratio A/F is lean, namely leaner than the stoichiometric air-fuel ratio, so that it is inferred that O2 is being adsorbed onto the catalyst, the O2 desorption/adsorption quantity calculation section 100 chooses a smaller quantity between the O2 adsorption capacity Cad and the stoichometric required O2 quantity Δadθo2.
Specifically, when the O2 quantity required to make the composition of the exhaust gases stoichiometric (Δadθo2) exceeds the capacity (Cred, Cad) of the catalyst, the amount of O2 actually used is limited to the capacity of the catalyst. Meanwhile, when the capacity (Cred, Cad) of the catalyst exceeds the O2 quantity required to make the composition of the exhaust gases stoichiometric (Δadθo2), the amount of O2 actually used is the required O2 quantity. The processing of the O2 desorption/adsorption quantity calculation section 100 is provided in view of this.
The value chosen by the O2 desorption/adsorption quantity calculation section 100 (Cred or Cad or Δadθo2) and the number ρo2 of the adsorption sites are fed to an adsorption-ratio decrease calculation section 101, and using these input values, the adsorption-ratio decrease calculation section 101 calculates a decrease Δrθo2 in adsorption ratio of the catalyst caused by the reaction of O2. The decrease Δrθo2 in adsorption ratio thus calculated and the O2 adsorption ratio θo2 in the last processing cycle are fed to an adsorption ratio update section 102, and the adsorption ratio update section 102 calculates the latest O2 adsorption ratio θO2 according to the equation (21) below (O2 storage index calculation element)

θO2=θO2+Δrθo2  (21)

The O2 adsorption ratio θO2 calculated by the adsorption ratio update section 102 is limited by an upper clipping section 103 to the upper limit 1.0, and limited by a lower clipping section 104 to the lower limit 0.
The O2 adsorption ratio θO2 of the catalyst 9 calculated in the above-described manner is fed to an air-fuel ratio determination section 111 of the fuel-recovery-period conversion heat calculation section 110 shown in FIG. 4. When determining that the catalyst 9 is not exhibiting the O2 storage function from the O2 adsorption ratio θO2 being 0 (O2 desorption limit) or 1.0 (O2 adsorption limit), the air-fuel ratio determination section 111 chooses the exhaust air-fuel ratio A/F for the air-fuel ratio on the catalyst (air-fuel ratio estimation element). When determining that the catalyst 9 is exhibiting the O2 storage function from the O2 adsorption ratio θO2 satisfying the condition 0<θo2<1.0, the air-fuel ratio determination section 111 chooses the stoichometric air-fuel ratio for the air-fuel ratio on the catalyst (air-fuel ratio estimation element).
The air-fuel ratio on the catalyst determined by the air-fuel ratio determination section 111, the engine revolving rate Ne detected by the revolving rate sensor 23 (operating state detection element) and the charging efficiency Ec calculated on the basis of the operating state of the engine 1 (operating state detection element) are fed to a ΔHc calculation section 112, and the ΔHc calculation section 112 calculates the CO conversion reaction heat per unit time ΔHc (kJ/sec), referring to a map prepared in advance, on the basis of the engine revolving rate Ne, the charging efficiency Ec and the air-fuel ratio on the catalyst (conversion reaction heat calculation element).
The CO conversion reaction heat ΔHc thus calculated and the value of the fuel-cut flag are fed to a fuel-cut determination section 113. When determining that the fuel cut is not being performed from the fuel-cut flag reset to the OFF-value, the fuel-cut determination section 113 chooses the CO conversion reaction heat ΔHc, and when determining that the fuel cut is being performed from the fuel-cut flag set to the ON-value, the fuel-cut determination section 113 chooses “0”. The value chosen by the fuel-cut determination section 113 and the processing period f of the ECU 21 are fed to a period conversion section 114, and using these input values, the period conversion section 114 calculates the CO conversion reaction heat ΔHc (kJ) in the processing period f. Thus, while the fuel cut is not being performed, the fuel-recovery-period conversion heat calculation section 110 sends out the CO conversion reaction heat ΔHc, and while the fuel cut is being performed, sends out “0”.
Next, the estimation of the definitive catalyst temperature Tcatf will be described.
FIG. 6 is a control block diagram showing how the processing by a catalyst temperature estimation section 120 of the ECU 21 goes on. The CO oxidation reaction heat ΔHr in the fuel cut period calculated by the fuel-cut-period reaction heat calculation section 50, the CO conversion reaction heat ΔHc in the fuel recovery period calculated by the fuel-recovery-period reaction heat calculation section 110, the heat transfer ΔHt between the exhaust gases and the catalyst calculated by the heat transfer calculation section 80 are fed to a total heat calculation section 121. As stated above, in the fuel cut period, the CO oxidation reaction heat ΔHr and the exhaust-catalyst heat transfer ΔHt are calculated as effective values, while in the fuel recovery period, the CO conversion reaction heat ΔHc and the exhaust-catalyst heat transfer ΔHt are calculated as effective values. Thus, in either case, the total heat calculation section 121 sums the two effective values, thereby obtaining total heat ΣH that causes a rise in catalyst temperature.
The total heat ΣH thus calculated is fed to a temperature rise calculation section 122, which calculates a temperature rise ΔT of the catalyst 9 from the total heat ΣH and the heat capacity of the catalyst 9 set in advance, and feeds it to a catalyst temperature calculation section 123. The catalyst temperature calculation section 123 calculates provisional catalyst temperature Tcat in the present processing cycle, by adding the temperature rise ΔT to the provisional catalyst temperature Tcat (catalyst temperature before subjected to the filtering described below) in the last processing cycle. The provisional catalyst temperature Tcat thus calculated is stored for use as a previous value in the next processing cycle, and on the basis of this provisional catalyst temperature Tcat, definitive catalyst temperature Tcatf is calculated by a three-stage filter simulating thermal conduction in the catalyst (catalyst temperature estimation element).
A gain calculation section 124 of the three-stage filter sets a filter gain Kca, referring to a map prepared in advance, on the basis of the intake quantity Q. As shown in FIG. 6, filter gain Kca increases with intake quantity Q. The filter gain Kca thus set and the provisional catalyst temperature Tcat are fed to a first filter section 125. On the basis of the provisional catalyst temperature Tcat, the filter gain Kca, and first catalyst temperature cell Tf1 (catalyst temperature after first filtering) in the last processing cycle, the first filter section 125 calculates first catalyst temperature cell Tf1 in the present processing cycle, and the first catalyst temperature cell Tf1 thus calculated is stored and also sent to a second filter section 126.
The second filter section 126 calculates second catalyst temperature cell Tf2 in the present processing cycle, on the basis of the first catalyst temperature cell Tf1, the filter gain Kca, and second catalyst temperature cell Tf2 (catalyst temperature after second filtering) in the last processing cycle, and the second catalyst temperature cell Tf2 thus calculated is stored and also sent to a third filter section 127.
The third filter section 127 calculates third catalyst temperature cell Tf3 in the present processing cycle, on the basis of the second catalyst temperature cell Tf2, the filter gain Kca, and third catalyst temperature cell Tf3 (catalyst temperature after third filtering) in the last processing cycle, and the third catalyst temperature cell Tf3 thus calculated is stored and also sent out as the definitive catalyst temperature Tcatf. When the intake quantity Q is greater and therefore the filter gain Kca is greater, the filter sections 125 to 127 each produces an output value less affected by a previous value and more affected by a current value. Consequently, the definitive catalyst temperature Tcatf taking account of the thermal conduction in the catalyst which correlates with the intake quantity Q is calculated.
As explained above, the present embodiment of the catalyst temperature estimation device calculates, in the fuel cut period, the CO oxidation reaction heat ΔHr produced when CO adsorbed on the catalyst in the quantity CATco reacts with O2 existing in the exhaust gases in the quantity GASo2, and estimates the temperature of the catalyst of the underfloor catalytic converter 9 on the basis of this CO oxidation reaction heat ΔHr and the exhaust-catalyst heat transfer ΔHt. Meanwhile, in the fuel recovery period, it calculates the CO conversion reaction heat ΔHc produced when CO in the exhaust gases is converted by O2 adsorbed on the catalyst, and estimates the temperature of the catalyst of the underfloor catalytic converter 9 on the basis of this CO conversion reaction heat ΔHc and the heat transfer ΔHt. Thus, in either of the fuel cut period and the fuel recovery period, the movements of the catalyst temperature Tcatf can be estimated accurately. By appropriately performing engine control on the basis of this catalyst temperature Tcatf, the catalytic converter 9 can be prevented from being damaged by a temperature rise beyond the allowable temperature limit.
FIG. 7 is a test result showing how the catalyst temperature is estimated by the present invention and by the conventional technique, when the fuel cut and fuel recovery is repeated with a short period. Such operating state can be produced when the driver repeats activation and deactivation of the accelerator frequently. Due to the CO oxidation reaction heat ΔHr produced in the fuel cut period and the CO conversion reaction heat ΔHc produced in the fuel recovery period, the actual catalyst temperature rises rapidly compared with the exhaust gas temperature Tex. The values estimated by the conventional technique hardly reflect this temperature rise. Meanwhile, by the estimation method of the present embodiment, estimated values Tcatf extremely close to the actual catalyst temperature are obtained. This test result confirms the positive effects of the present invention.
Further, the present embodiment of the catalyst temperature estimation device calculates the amount CATco of CO adsorbed on the catalyst 9 at present, on the basis of the CO desorption rate Rred at which CO is desorbed from the catalyst and the CO adsorption rate Rad at which CO is adsorbed onto the catalyst, which rates are obtained according to the actual process. In addition, when the amount GASco of CO in the exhaust gases is less than the CO adsorption rate Rad on the catalyst so that the CO adsorption rate Rad is limited to the CO quantity GASco, the CO quantity GASco is chosen for the CO adsorption rate Rad. This allows the adsorbed CO quantity CATco and therefore the CO oxidation reaction heat ΔHr to be calculated with improved accuracy.
Further, the calculation of the CO oxidation reaction heat ΔHr takes account of not only the amount CATco of CO adsorbed on the catalyst and the amount GASo2 of O2 in the exhaust gases, but also the rate r of reaction of CO and O2. This allows the CO oxidation reaction heat ΔHr to be calculated with more improved accuracy.
In the above, an embodiment of the present invention has been described. The present invention is, however, not limited to the described embodiment. For example, although in the above embodiment, the catalyst temperature Tcarf is estimated in both the fuel cut period and the fuel recovery period, the estimation of the catalyst temperature Tcarf in the fuel recovery period is not indispensable. For example, the embodiment may be modified to estimate the catalyst temperature Tcarf only in the fuel cut period, on the basis the CO oxidation reaction heat ΔHr. Further, although in the described embodiment, the engine 1 is a direct-injection spark-ignition four-cylinder-in-line gasoline engine, the present invention is not limited to this type of engine; the present invention is applicable to that type of gasoline engine in which fuel is supplied to the intake ports, and also to the diesel engine.

## Claims

1. A catalyst temperature estimation device, comprising:
a catalytic converter carrying a catalyst disposed in an exhaust passage of an engine for converting exhaust gases to less harmful substances;
an adsorbed unburnt-fuel-component quantity calculation element for calculating the amount of an unburnt fuel component adsorbed on the catalyst in a fuel cut period of the engine based on a number of all adsorption sites, which is determined depending on an activated state of the catalytic converter capable of adsorbing the unburnt fuel component;
an exhaust O2 quantity calculation element for calculating the amount of O2 in the exhaust gases in the fuel cut period;
an oxidation reaction heat calculation element for calculating reaction heat produced when the adsorbed unburnt fuel component in the amount calculated by the adsorbed unburnt-fuel-component quantity calculation element reacts with O2 in the exhaust gases in the amount calculated by the exhaust O2 quantity calculation element; and
a catalyst temperature estimation element for estimating the temperature of the catalyst in the fuel cut period, on the basis of the reaction heat calculated by the oxidation reaction heat calculation element.

2. The catalyst temperature estimation device according to claim 1, further comprising:
a heat transfer calculation element for calculating heat transfer between the exhaust gases from the engine and the catalyst in the fuel cut period, wherein
the catalyst temperature estimation element estimates the temperature of the catalyst from the heat transfer calculated by the heat transfer calculation element and the reaction heat calculated by the heat transfer calculation element.

3. The catalyst temperature estimation device according to claim 1, further comprising:
a reaction rate calculation element for calculating the reaction rate at which the unburnt fuel component on the catalyst reacts with O2 in the exhaust gases, on the basis of the activated state of the catalyst, wherein
the oxidation reaction heat calculation element calculates the reaction heat on the basis of the smallest one of the amount of the unburnt fuel component adsorbed, the amount of O2 in the exhaust gases, and the reaction rate calculated.

4. The catalyst temperature estimation device according to claim 1, wherein
the catalyst temperature estimation element corrects the estimated temperature of the catalyst, by means of a three-stage filter simulating thermal conduction in the catalyst.

5. The catalyst temperature estimation device according to claim 1, further comprising:
an operating state detection element for detecting the operating state of the engine;
an O2 storage index calculation element for calculating an O2 storage index correlating with the amount of O2 adsorbed on the catalyst;
an air-fuel ratio estimation element for estimating air-fuel ratio on the catalyst; and
a conversion reaction heat calculation element for calculating reaction heat produced by a conversion reaction on the catalyst in a fuel recovery period of the engine, on the basis of the operating state of the engine detected by the operating state detection element, the O2 storage index calculated by the O2 storage index calculation element, and the air-fuel ratio estimated by the air-fuel ratio estimation element, wherein
the catalyst temperature estimation element estimates the temperature of the catalyst in the fuel cut period on the basis of the reaction heat calculated by the oxidation reaction heat calculation element, and in the fuel recovery period on the basis of the reaction heat calculated by the conversion reaction heat calculation element.

6. The catalyst temperature estimation device according to claim 5, wherein
when, from the O2 storage index calculated by the O2 storage index calculation element, determining that the catalyst is at a level close to the upper or lower limit of O2 adsorption, the air-fuel ratio estimation element determines that the air-fuel ratio on the catalyst equals the exhaust air-fuel ratio detected by an exhaust air-fuel ratio detection element, and when determining that the catalyst is not at a level close to the upper or lower limit of adsorption, determines that the air-fuel ratio on the catalyst is the stoichiometric air-fuel ratio.

7. The catalyst temperature estimation device according to claim 1, further comprising:
an adsorption rate calculation element for calculating the adsorption rate at which the unburnt fuel component is adsorbed onto the catalyst in the fuel cut period; and
a desorption rate calculation element for calculating the desorption rate at which the unburnt fuel component is desorbed from the catalyst in the fuel cut period, wherein
the adsorbed unburnt-fuel-component quantity calculation element calculates the amount of the unburnt fuel component adsorbed on the catalyst, on the basis of the adsorption rate calculated by the adsorption rate calculation element and the desorption rate calculated by the desorption rate calculation element.

8. The catalyst temperature estimation device according to claim 7, wherein
when the amount of the unburnt fuel component in the exhaust gases is less than the calculated adsorption rate of the unburnt fuel component, the adsorption rate calculation element determines that the adsorption rate equals the amount of the unburnt fuel component in the exhaust gases.

9. A catalyst temperature estimation device, comprising:
a catalytic converter carrying a catalyst disposed in an exhaust passage of an engine for converting exhaust gases to less harmful substances;
an adsorbed unburnt-fuel-component quantity calculation element for calculating the amount of an unburnt fuel component adsorbed on the catalyst in a fuel cut period of the engine;
an exhaust O2 quantity calculation element for calculating the amount of O2 in the exhaust gases in the fuel cut period;
an oxidation reaction heat calculation element for calculating reaction heat produced when the adsorbed unburnt fuel component in the amount calculated by the adsorbed unburnt-fuel-component quantity calculation element reacts with O2 in the exhaust gases in the amount calculated by the exhaust O2 quantity calculation element;
a catalyst temperature estimation element for estimating the temperature of the catalyst in the fuel cut period, on the basis of the reaction heat calculated by the oxidation reaction heat calculation element;
an operating state detection element for detecting the operating state of the engine;
an O2 storage index calculation element for calculating an O2 storage index correlating with the amount of O2 adsorbed on the catalyst;
an air-fuel ratio estimation element for estimating air-fuel ratio on the catalyst; and
a conversion reaction heat calculation element for calculating reaction heat produced by a conversion reaction on the catalyst in a fuel recovery period of the engine, on the basis of the operating state of the engine detected by the operating state detection element, the O2 storage index calculated by the O2 storage index calculation element, and the air-fuel ratio estimated by the air-fuel ratio estimation element, wherein
the catalyst temperature estimation element estimates the temperature of the catalyst in the fuel cut period on the basis of the reaction heat calculated by the oxidation reaction heat calculation element, and in the fuel recovery period on the basis of the reaction heat calculated by the conversion reaction heat calculation element wherein, when, from the O2 storage index calculated by the O2 storage index calculation element, determining that the catalyst is at a level close to the upper or lower limit of O2 adsorption, the air-fuel ratio estimation element determines that the air-fuel ratio on the catalyst equals the exhaust air-fuel ratio detected by an exhaust air-fuel ratio detection element, and when determining that the catalyst is not at a level close to the upper or lower limit of adsorption, determines that the air-fuel ratio on the catalyst is the stoichiometric air-fuel ratio.

10. A catalyst temperature estimation device, comprising:
a catalytic converter carrying a catalyst disposed in an exhaust passage of an engine for converting exhaust gases to less harmful substances;
an adsorbed unburnt-fuel-component quantity calculation element for calculating the amount of an unburnt fuel component adsorbed on the catalyst in a fuel cut period of the engine;
an exhaust O2 quantity calculation element for calculating the amount of O2 in the exhaust gases in the fuel cut period;
an oxidation reaction heat calculation element for calculating reaction heat produced when the adsorbed unburnt fuel component in the amount calculated by the adsorbed unburnt-fuel-component quantity calculation element reacts with O2 in the exhaust gases in the amount calculated by the exhaust O2 quantity calculation element;
a catalyst temperature estimation element for estimating the temperature of the catalyst in the fuel cut period, on the basis of the reaction heat calculated by the oxidation reaction heat calculation element; and
a reaction rate calculation element for calculating the reaction rate at which the unburnt fuel component on the catalyst reacts with O2 in the exhaust gases, on the basis of the activated state of the catalyst, wherein
the oxidation reaction heat calculation element calculates the reaction heat on the basis of the smallest one of the amount of the unburnt fuel component adsorbed, the amount of O2 in the exhaust gases, and the reaction rate calculated.

11. A catalyst temperature estimation device, comprising:
a catalytic converter carrying a catalyst disposed in an exhaust passage of an engine for converting exhaust gases to less harmful substances;
an adsorbed unburnt-fuel-component quantity calculation element for calculating the amount of an unburnt fuel component adsorbed on the catalyst in a fuel cut period of the engine;
an exhaust O2 quantity calculation element for calculating the amount of O2 in the exhaust gases in the fuel cut period;
an oxidation reaction heat calculation element for calculating reaction heat produced when the adsorbed unburnt fuel component in the amount calculated by the adsorbed unburnt-fuel-component quantity calculation element reacts with O2 in the exhaust gases in the amount calculated by the exhaust O2 quantity calculation element;
a catalyst temperature estimation element for estimating the temperature of the catalyst in the fuel cut period, on the basis of the reaction heat calculated by the oxidation reaction heat calculation element;
an adsorption rate calculation element for calculating the adsorption rate at which the unburnt fuel component is adsorbed onto the catalyst in the fuel cut period; and
a desorption rate calculation element for calculating the desorption rate at which the unburnt fuel component is desorbed from the catalyst in the fuel cut period, wherein
the adsorbed unburnt-fuel-component quantity calculation element calculates the amount of the unburnt fuel component adsorbed on the catalyst, on the basis of the adsorption rate calculated by the adsorption rate calculation element and the desorption rate calculated by the desorption rate calculation element.

12. The catalyst temperature estimation device according to claim 11, wherein
when the amount of the unburnt fuel component in the exhaust gases is less than the calculated adsorption rate of the unburnt fuel component, the adsorption rate calculation element determines that the adsorption rate equals the amount of the unburnt fuel component in the exhaust gases.