Reinke, Paul E. (Rochester Hills, MI, US)

Plaque It!

11/668888

10/07/2008

01/30/2007

Images are available in PDF form when logged in. To view PDFs, Login  or  Create Account (Free!)

Click for automatic bibliography generation

GM Global Technology Operations, Inc. (Detroit, MI, US)

123/673

123/434, 123/678

F02M1/00; F02D41/00

123/673, 123/689, 123/677, 123/436, 123/435, 123/434, 123/678

6640539	Engine control for low emission vehicle starting	November, 2003	Lewis et al.	60/284
6666021	Adaptive engine control for low emission vehicle starting	December, 2003	Lewis et al.	60/289
6715280	Method for low emission vehicle starting with improved fuel economy	April, 2004	Lewis et al.	60/274
7185487	Adaptive engine control for low emission vehicle starting	March, 2007	Lewis et al.	60/289
7293406	Engine control for low emission vehicle starting	November, 2007	Lewis et al.	60/284

Kwon, John T.

What is claimed is:

1. A control system for controlling a valve in a fuel system of engine, comprising: an engine airflow sensor that senses a mass engine airflow (M_engair); a barometric pressure sensor that senses a barometric pressure (P_baro); an ambient temperature sensor that senses an ambient temperature (T_amb); an actuator determination module that determines an effective area (A_eff) of said valve based on at least one of said M_engair, said P_baro, and said T_amb; and a duty cycle (DC) calculation module that determines a first duty cycle (DC) of said valve based on said A_eff.

2. The system of claim 1 wherein said first DC is based on said A_eff, a first slope (M₁), a first offset (B₁), and a first DC equation, $D{C}_{\mathrm{first}}=\left(\frac{\mathrm{FP}}{M_1}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_I}{{\left(T_{\mathrm{amb}}\right)}^{\frac{1}{2}}}\right)}-\frac{B_1}{M_1}.$

3. The system of claim 2 wherein said A_eff is based on said M_engair, a purge flow percentage (FP) of said valve that is based on said M_engair, a first predetermined constant (C₁), a pressure driving function δ, a first pressure (P₁), said T_amb, and an A_eff equation $A_{\mathrm{eff}}=\mathrm{FP}*\frac{M_{\mathrm{engair}}}{C_1*\delta *P_1*{\left(\frac{1}{T_{\mathrm{amb}}}\right)}^{\frac{1}{2}}}.$

4. The system of claim 3 wherein said P₁ is based on said P_baro, a second predetermined constant (K₁), an active purge flow percentage (FP_act), said M_engair, said T_amb, said P_baro, and a P₁ equation $P_1=P_{\mathrm{baro}}-K*{\left({\mathrm{FP}}_{\mathrm{act}}*M_{\mathrm{engair}}*\frac{T_{\mathrm{amb}}}{P_{\mathrm{baro}}}\right)}^2.$

5. The system of claim 2 wherein said DC calculation module computes a compensated DC that is based on said A_eff, a second slope (M₂), a second offset (B₂), and a compensated DC equation $D{C}_{\mathrm{first}}=\left(\frac{\mathrm{FP}}{M_1}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_1}{{\left(T_{\mathrm{amb}}\right)}^{\frac{1}{2}}}\right)}-\frac{B_2}{M_2}$ when said first DC one of exceeds a DC threshold and falls below said DC threshold, said M₂ and said B₂ are based on said B₁, and said first DC is set equal to said compensated DC.

6. The system of claim 2 wherein said B₁ is based on at least one of: a manifold vacuum (MV), said V supplied to said valve, and said T_amb, wherein said MV is based on a MV equation (P_baro−P₂).

7. The system of claim 5 wherein said DC calculation module commands said valve with a first delivered DC when a difference between said first DC and a second DC exceeds a predetermined comparison threshold, wherein said first delivered DC is based on said first DC and said second DC; and said DC calculation module commands said valve with a second delivered DC when said difference between said first DC and said second DC falls below said predetermined comparison threshold, wherein said first delivered DC is set equal to said second delivered DC.

8. The system of claim 7 wherein said second DC is set equal to said first DC.

9. The system of claim 4 wherein said δ is based on a pressure ratio of said P₁ to said P₂; wherein said δ is equal to a δ constant when said pressure ratio falls below a δ threshold; and wherein said δ is based on said P₁, said P₂, and a δ equation δ=3.8639(PR^1.42857−PR^1.71428)^0.5 when said pressure ratio exceeds a δ threshold.

FIELD OF THE INVENTION

The present invention relates to internal combustion engines, and more particularly to a system for controlling engine emissions.

BACKGROUND OF THE INVENTION

Internal combustion engines combust an air/fuel (A/F) mixture within cylinders to drive pistons and to provide drive torque. Air is delivered to the cylinders via a throttle and an intake manifold. A fuel injection system supplies fuel from a fuel tank to provide fuel to the cylinders based on a desired A/F mixture. To prevent release of fuel vapor, a vehicle may include an evaporative emissions system which includes a canister that absorbs fuel vapor from the fuel tank, a canister vent valve, and a purge valve. The canister vent valve allows air to flow into the canister. The purge valve supplies a combination of air and vaporized fuel from the canister to the intake system.

Closed-loop control systems adjust inputs of a system based on feedback from outputs of the system. By monitoring the amount of oxygen in the exhaust, closed-loop fuel control systems manage fuel delivery to an engine. Based on an output of oxygen sensors, an engine control module adjusts the fuel delivery to match an ideal A/F ratio (14.7 to 1). By monitoring engine speed variation at idle, closed-loop speed control systems manage engine intake airflows and spark advance.

Typically, the fuel tank stores liquid fuel such as gasoline, diesel, methanol, or other fuels. The liquid fuel may evaporate into fuel vapor which increases pressure within the fuel tank. Evaporation of fuel is caused by energy transferred to the fuel tank via radiation, convection, and/or conduction. An evaporative emissions control (EVAP) system is designed to store and dispose of fuel vapor to prevent release. More specifically, the EVAP system returns the fuel vapor from the fuel tank to the engine for combustion therein.

The EVAP system includes an evaporative emissions canister (EEC) and a purge valve. When the fuel vapor increases within the fuel tank, the fuel vapor flows into the EEC. A purge valve controls the flow of the fuel vapor from the EEC to the intake manifold. The purge valve may be modulated between open and closed positions to adjust the flow of fuel vapor to the intake manifold. Improper operation of the purge valve may cause a variety of undesirable conditions including, but not limited to, idle surge, steady throttle surge, and undesirable emission levels.

SUMMARY OF THE INVENTION

A control system for controlling a valve in a fuel system of engine including an engine airflow sensor that senses a mass engine airflow (M _engair ), a barometric pressure sensor that senses a barometric pressure (P _baro ), and an ambient temperature sensor that senses an ambient temperature (T _amb ). An actuator determination module determines an effective area (A _eff ) of the valve based on at least one of the M _engair , the P _baro , and the T _amb . A duty cycle (DC) calculation module that determines a first duty cycle (DC) of the valve based on the A _eff .

In other features, the first DC is based on the A _eff , a first slope (M ₁ ), a first offset (B ₁ ), and a first DC equation,

${\mathrm{DC}}_{\mathrm{first}}=\left(\frac{\mathrm{FP}}{M_1}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_1}{{\left(T_{\mathrm{amb}}\right)}^{1\text{/}2}}\right)}-\frac{B_1}{M_1}.$
The A _eff is based on the M _engair , a purge flow percentage (FP) of the valve that is based on the M _engair , a first predetermined constant (C ₁ ), a pressure driving function δ, a first pressure (P ₁ ), the T _amb , and an A _eff equation

$A_{\mathrm{eff}}=\mathrm{FP}*\frac{M_{\mathrm{engair}}}{C_1*\delta *P_1*{\left(\frac{1}{T_{\mathrm{amb}}}\right)}^{1\text{/}2}}.$
The P ₁ is based on the P _baro , a second predetermined constant (K ₁ ), an active purge flow percentage (FP _act ), the M _engair , the T _amb , the P _baro , and a P ₁ equation,

$P_1=P_{\mathrm{baro}}-K*{\left({\mathrm{FP}}_{\mathrm{act}}*M_{\mathrm{engair}}*\frac{T_{\mathrm{amb}}}{P_{\mathrm{baro}}}\right)}^2.$

In other features, the DC calculation module computes a compensated DC that is based on the A _eff , a second slope (M ₂ ), a second offset (B ₂ ), and a compensated DC equation

${\mathrm{DC}}_{\mathrm{first}}=\left(\frac{\mathrm{FP}}{M_2}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_1}{{\left(T_{\mathrm{amb}}\right)}^{1\text{/}2}}\right)}-\frac{B_2}{M_2}$
when the first DC one of exceeds a DC threshold and falls below the DC threshold, the M ₂ and the B ₂ are based on the B ₁ , and the first DC is set equal to the compensated DC. The B ₁ is based on at least one of: a manifold vacuum (MV), the V supplied to the valve, and the T _amb , wherein the MV is based on a MV equation (P _baro −P ₂ ).

In other features, the DC calculation module commands the valve with a first delivered DC when a difference between the first DC and a second DC exceeds a predetermined comparison threshold, wherein the first delivered DC is based on the first DC and the second DC. The DC calculation module commands the valve with a second delivered DC when the difference between the first DC and the second DC falls below the predetermined comparison threshold, wherein the first delivered DC is set equal to the second delivered DC. The second DC is set equal to the first DC.

In other features, the δ is based on a pressure ratio of the P ₁ to the P ₂ . The δ is equal to a δ constant when the pressure ratio falls below a δ threshold. δ is based on the P ₁ , the P ₂ , and a δ equation δ=3.8639(PR ^1.42857 −PR ^1.71428 ) ^0.5 when the pressure ratio exceeds a δ threshold.

A system for controlling fuel vapor in a fuel system of an engine includes an engine airflow sensor that senses mass engine airflow (M _engair ). A purge initialization module determines an amount of the fuel vapor in the fuel system based on an active flow percentage (FP _act ) and a desired flow percentage (FP _des ), wherein the FP _act and the FP _des are based on the M _engair . A purge learn module controls the amount of the fuel vapor in the engine.

In other features, The purge initialization module increments the FP _act when the FP _act is below the FP _des . The purge initialization module calculates an adjusted second array value (PLM _adj ) corresponding to the M _engair based on an active second array value (PLM _act ), an active first array value (LTM _act ), the FP _des , the FP _act , and a PLM _adj equation

$\left({\mathrm{PLM}}_{\mathrm{adj}}\right)-\left({\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_{\mathrm{act}}\right)*\left(\frac{{\mathrm{FP}}_{\mathrm{des}}}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_{\mathrm{act}}.$
The PLM _adj , the PLM _act , and the LTM _act are multipliers to a fuel rate. The PLM _adj , the PLM _act , the LTM _act , the FP _act , and the FP _des are indexed by the M _engair .

In other features, the purge learn module calculates a limited flow percentage (FP _lim ) for a predetermined second array threshold (PLM _lim ) when the PLM _lim exceeds the PLM _adj . The purge learn module calculates at least one second array value (PLM _i ) based the FP _lim , a first array (LTM _i ), the LTM _act , the PLM _act , the FP _act , and a limit PLM equation

$\left({\mathrm{PLM}}_i\right)=\left(\left(\frac{{\mathrm{LTM}}_i}{{\mathrm{LTM}}_{\mathrm{act}}}\right)*{\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_i\right)*\left(\frac{{\mathrm{FP}}_{\lim }}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_i.$
The i represents a current operating position. The purge learn module sets the FP _act equal to the FP _lim , and the learn module decreases the FP _act by a predetermined FP _act constant (Y) when the FP _des exceeds the FP _act and the PLM _lim exceeds the PLM _act , and the purge learn module increases the FP _act by the Y when the FP _des exceeds the FP _act and the PLM _act exceeds the PLM _lim .

In other features, the purge learn module sets the FP _act equal to the FP _des when the FP _act exceeds the FP _des , and wherein the purge learn module calculates at least one second array value (PLM _i ) based the LTM _i , LTM _act , PLM _act , a current position FP (FP _i ), the FP _act , the LTM _i , and second array equation

$\left({\mathrm{PLM}}_i\right)=\left(\left(\frac{{\mathrm{LTM}}_i}{{\mathrm{LTM}}_{\mathrm{act}}}\right)*{\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_i\right)*\left(\frac{{\mathrm{FP}}_i}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_i.$

In other features, the system further comprises a threshold compensation module that determines a second array value (PLM _i ) when the amount of sad fuel vapor falls below a fuel vapor threshold. The second array is a multiplier to a fuel rate. The i represents a current position. The threshold compensation module determines a threshold flow percentage (FPDC _thres ) when a duty cycle (DC) of a valve in the fuel system exceeds a DC threshold, wherein the FPDC _thres is based on a gas constant (C ₁ ), a pressure driving function (δ), a corrected pressure (P _corr ), an ambient temperature (T _amb ), an engine airflow (m _engair ), and a FP ^thres equation

$\left({\mathrm{FP}}_{100\mathrm{DC}}\right)=\frac{\left(\begin{array}{c}{C}_1*\delta *P_1\\ T_{\mathrm{amb}}^{1\text{/}2}\end{array}\right)}{M_{\mathrm{eng}}}.$
The δ is calculated based on the P _corr , a manifold absolute pressure (P ₂ ), an air constant (C ₁ ), and the δ equation δ=3.8639(PR ^1.42857 −PR ^1.71428 ) ^0.5 . The P _corr is based on a barometric pressure (P _baro ), a DC threshold, the DC, a second slope (M ₂ ), a second offset (B ₂ ), and a first pressure (P ₁ ).

In other features, the threshold compensation module determines a DC threshold second array value (PLM _thres ) based on the FP _act , an active first array value (LTM _act ), the (FP _thres ), the FP _act , and a PLM _thres equation

$\left({\mathrm{PLM}}_{100\mathrm{DC}}\right)=\left({\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_{\mathrm{act}}\right)*\left(\frac{{\mathrm{FP}}_{100\mathrm{DC}}}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_{\mathrm{act}}.$
The threshold compensation module determines at least one first array value (LTM _i ) based on the PLM _i , an active second array value (PLM _act ), an active first array value (LTM _act ), a flow percentage (FP _i ), the FP _act , and a LTM _i equation

$\left({\mathrm{LTM}}_i\right)=\frac{{\mathrm{PLM}}_i}{\left(\left(\frac{{\mathrm{PLM}}_{\mathrm{act}}}{{\mathrm{LTM}}_{\mathrm{act}}}\right)-1\right)*\left(\frac{{\mathrm{FP}}_i}{{\mathrm{FP}}_{\mathrm{act}}}\right)+1}$
when the engine is turned off.

In other features, the purge initialization module determines the amount of the fuel vapor when an off time of the engine exceeds an off time threshold.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will become more fully understood from the detailed description and the accompanying drawings, wherein:

FIG. 1 is a functional block diagram of an engine control system and a fuel system according to the present invention;

FIG. 2 is a functional block diagram of an engine control module (ECM) according to the present invention;

FIG. 3 is a flow diagram illustrating steps of a method for determining flow variables according to the present invention;

FIG. 4 is a flow diagram illustrating steps of a method for determining the actuator variables according to the present invention;

FIG. 5 is a flow diagram illustrating steps of a method for determining transient duty cycle (DC) compensation according to the present invention;

FIG. 6 is a flow diagram illustrating steps of a method for initializing purge control according to the present invention;

FIG. 7 is a flow diagram illustrating a method for initializing purge control according to the present invention; and

FIG. 8 is a flow diagram illustrating steps of a method for controlling vapor during periods of 100% DC according to the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The following description of the preferred embodiment is merely exemplary in nature and is in no way intended to limit the invention, its application, or uses. For purposes of clarity, the same reference numbers will be used in the drawings to identify similar elements. As used herein, the term module refers to an application specific integrated circuit (ASIC), an electronic circuit, a processor (shared, dedicated, or group) and memory that execute one or more software or firmware programs, a combinational logic circuit, and/or other suitable components that provide the described functionality.

Referring to FIG. 1, a vehicle 10 includes an engine system 12 and a fuel system 14 . One or more control modules 16 communicate with the engine and fuel systems 12 , 14 . The fuel system 14 selectively supplies liquid and/or fuel vapor to the engine system 12 , as will be described in further detail below.

The engine system 12 includes an engine 18 , a fuel injection system 20 , an intake manifold 22 , and an exhaust manifold 24 . Air is drawn into the intake manifold 22 through a throttle 26 . The throttle 26 regulates mass air flow into the intake manifold 22 . Air within the intake manifold 22 is distributed into cylinders 28 . The air is mixed with fuel and the air/fuel (A/F) mixture is combusted within cylinders 28 of the engine 18 . Although two cylinders 28 are illustrated, it can be appreciated that the engine 18 can include any number of cylinders 28 including, but not limited to 1, 3, 4, 5, 6, 8, and 12 cylinders. The fuel injection system 20 includes liquid injectors that inject liquid fuel into the cylinders 28 .

Exhaust flows through the exhaust manifold 24 and is treated in a catalytic converter 30 . An exhaust oxygen sensor 32 (e.g., a wide-range A/F ratio sensor) communicates exhaust A/F ratio signals to the control module 16 . A voltage (V) sensor 33 communicates an operating V supplied to a purge valve 34 to the control module 16 . A mass engine air flow (M _engair ) sensor 36 is located within an air inlet and communicates a M _engair signal based on the mass of air flowing into the engine system 12 to the control module 16 . An intake manifold absolute pressure (P ₂ ) sensor 38 senses the pressure within the intake manifold 22 of the engine and communicates a P ₂ signal to the control module 16 . Ambient temperature (T _amb ) and barometric pressure (P _baro ) signals are generated by T _amb and P _baro sensors 39 , 40 , respectively.

The control module 16 controls the fuel and air provided to the engine based on oxygen sensor signals and throttle valve position. This form of fuel control is also referred to as closed loop fuel control. Closed loop fuel control is used to maintain the A/F mixture at or close to an ideal stoichiometric A/F ratio by commanding a desired fuel delivery to match the airflow. Stoichiometry is defined as an ideal A/F ratio (e.g., 14.7 to 1 for gasoline). The engine control may command different airflow to compensate the engine speed changes during engine idle operation.

The engine system 12 operates in a lean condition (i.e., reduced fuel) when the A/F ratio is higher than the stoichiometric A/F ratio. The engine system 12 operates in a rich condition when the A/F ratio is below the stoichiometric A/F ratio. A fuel control factor helps determine whether the A/F ratio is within an ideal range (i.e., above a minimum value and below a maximum value). An exemplary fuel control factor includes a short term integrator (STI) that provides a rapid indication of fuel enrichment based on input from the oxygen sensor signals. For example, if the signals indicate an A/F ratio greater than a specified reference, the STI is increased a step. Conversely, if the signals indicate an A/F ratio less than the specified reference, the STI is decreased a step. A fuel control modifier monitors changes in the fuel control factor over a long term. An exemplary fuel control modifier includes a long term modifier (LTM). The LTM monitors STI and uses integration to produce its output.

The fuel system 14 includes a fuel tank 42 that contains liquid fuel and fuel vapor. A fuel inlet 44 extends from the fuel tank 42 to enable fuel filling. A fuel cap 46 closes the fuel inlet 44 and may include a bleed hole (not shown). A modular reservoir assembly (MRA) 48 is disposed within the fuel tank 42 and includes a fuel pump 50 . The MRA 48 includes a liquid fuel line 52 and a fuel vapor line 54 .

The fuel pump 50 pumps liquid fuel through the liquid fuel line 52 to the fuel injection system 20 of the engine 18 . A fuel vapor system includes the fuel vapor line 54 and a canister 56 . Fuel vapor flows through the fuel vapor line 54 into the canister 56 . A fuel vapor line 58 connects the purge valve 34 to the canister 56 . The control module 16 modulates the purge valve 34 to selectively enable fuel vapor to flow into the intake system of the engine 18 . The control module 16 modulates a canister vent valve 60 to selectively enable air flow from atmosphere into the canister 56 .

Referring now to FIG. 2, an engine control system 100 will be described in further detail. The engine control system 100 includes communicating with the exhaust oxygen sensor 32 , the M _engair sensor 36 , the T _amb sensor 39 , the P ₂ sensor 38 , the P _baro sensor 40 , and the V sensor 33 . A purge initialization module 70 receives signals from the exhaust oxygen sensor 32 and the M _engair sensor 36 . The purge initialization module 70 captures initialization variables and determines initial purge control conditions by determining the amount of fuel vapors present within the fuel system 14 .

A purge learn module 72 communicates with the initialization module 70 and receives signals from the exhaust oxygen sensor 32 . The purge learn module 72 determines an array of purge learn memory (PLM) values for a current engine operating condition and adaptively learns the PLM values for the remaining engine operating conditions. The PLM values correct for errors in the fuel system 14 when the purge valve 34 .

A threshold compensation module 74 receives signals from the M _engair sensor 36 , the T _amb sensor 39 , the P ₂ sensor 38 , and the P _baro sensor 40 . The threshold compensation module 74 communicates with the purge learn module 72 and the duty cycle (DC) calculation module 80 and calculates PLM values when the DC of the purge valve 34 has exceeded a DC threshold. The threshold compensation module 74 stores the PLM values in a long term memory (LTM) array (not shown) when engine 18 is off.

A flow variable initialization module 76 receives signals from the M _engair sensor 36 , the T _amb sensor 39 , the P ₂ sensor 38 , the P _baro sensor 40 , and the V sensor 33 . The flow variable initialization module 76 communicates with the purge initialization module 70 and the DC calculation module 80 . The flow variable initialization module 76 determines the initial flow variables utilized to control the operation of the purge valve 34 .

An actuator determination module 78 communicates with the flow variable initialization module 76 . The actuator determination module 78 determines the characteristics of an actuator (not shown) of the purge valve 34 based on parameters received from the flow variable initialization module 76 .

The DC calculation module 80 communicates with the actuator determination module 78 and the threshold compensation module 74 . The DC calculation module 80 determines the DC of the purge valve 34 for DC operation below the DC threshold and for DC operation that exceeds the DC threshold. The DC calculation module 80 commands the actuator of the purge valve 34 .

The engine control system 100 controls the purge valve 34 based on adjustments made to the DC of the purge valve 34 for distinct actuator characteristics of the purge valve 34 during various engine operating conditions. The engine control system 100 controls the introduction of purge flow into of the engine 18 based on the adjusted DC, thereby providing an accurate control of A/F mixture supplied to the engine 18 .

Referring now to FIG. 3, a method 300 for determining the flow variables in the engine control system 100 will be described in more detail. In the present implementation, an exemplary embodiment of the purge valve 34 will relate to the characteristics and operation of a 95 liters per minute (LPM) purge valve. Those skilled in the art can appreciate that various other embodiments of the purge valve 34 are contemplated within the scope of this invention including, but not limited to, a 71 LPM purge valve. The method 300 begins in step 302 . In step 304 , control initializes the purge control variables by processing signals provided by the M _engair sensor 36 , the P _baro sensor 40 , the P ₂ sensor 38 , the T _amb sensor 39 , and the V sensor 33 . Control also initializes DC _second . Control determines whether the purge control is enabled in step 306 . If control determines that the purge control is not enabled, control returns to step 304 . If control determines that the purge control is enabled, control proceeds to step 308 . In step 308 , control determines an active purge flow percentage (FP _act ). Control can determine the FP _act from sources including, but not limited to, a FP _act look-up table (not shown) and a predetermined purge closed loop equation (not shown). The FP _act look-up table is a function of M _engair . In an exemplary embodiment, FP _act can be approximately 3%, although it should be appreciated that other values of FP _act are anticipated.

In step 310 , control calculates a first pressure (P ₁ ) representing the pressure upstream from the inlet of the purge valve 34 . Generally, the pressure experienced by the fuel system 14 , P _baro , is subject to a drop in pressure as air flows through the purge line 58 resulting in P ₁ . P ₁ represents a calculated pressure at the inlet of the purge valve 34 . The value of P ₁ is determined according the following equation:

$\begin{array}{cc}{P}_1=P_{\mathrm{baro}}-K*{\left({\mathrm{FP}}_{\mathrm{act}}*M_{\mathrm{engair}}*\frac{T_{\mathrm{amb}}}{P_{\mathrm{baro}}}\right)}^2& \left(1\right)\end{array}$
where K represents a calibrated constant (e.g. 0.173). Although K may be 0.173 in the current example, other values for K are anticipated. In step 312 , control calculates the pressure ratio (PR) of a second pressure (P ₂ ) relative to P ₁ . P ₂ represents the pressure downstream from the purge valve 34 and is supplied by the P ₂ sensor 38 . In the current invention, P ₂ acts as a negative pressure upon the purge valve 34 . Once control has determined PR, control next determines the pressure delta function (δ) in 314 . δ represents a compressible gas characterization. In an exemplary embodiment, δ can range from 0 to 1. Control can determine δ from sources including, but not limited to, a δ look-up table (not shown) based on PR and/or a δ equation. δ can be calculated according to the following equation:
δ=1, for PR<PR _thres (2)
δ=3.8639( PR ^1.42857 −PR ^1.71428 ) ^0.5 for PR>PR _thres (3)
Where PR _thres represents a pressure ratio threshold and A ₁ represents an air constant. In an exemplary embodiment, δ is based on a PR _thres Of approximately 0.528. In step 316 , the method 300 ends.

Referring now to FIG. 4, a method 400 of determining actuator characteristics is discussed in more detail. Control begins the method 400 in step 402 . In step 403 , control calculates a voltage level supplied to the purge valve 34 . In step 404 , control calculates a manifold vacuum (MV) value. MV is used by control to correct for characteristics of distinct actuators of the purge valve 34 . (P ₁ -P ₂ ) represents a force upon the purge valve 34 that can be approximated by MV. The force results from the difference between P _baro delivered to the fuel system 14 and P ₂ acting downstream of the purge valve 34 . In various embodiments, the purge valve 34 can be pressure compensated, therefore eliminating the need to take into account MV.

In step 406 , control determines a first base purge valve offset (B ₁ ). Control can determine B ₁ for various operating conditions from sources including, but not limited to, a B ₁ look-up table (not shown). In various embodiments, the look-up table may be a function of, but is not limited to, MV and/or operating V. Generally, operation of the purge valve 34 is primarily affected by V and MV. Typically, the absolute temperature of the fuel vapor that is estimated by the intake manifold temperature, T _inlet , affects B ₁ to a lesser degree. For the sake of simplicity, T _inlet can be assumed to be T _amb . In an exemplary embodiment, the purge valve 34 operates at approximately 12V and approximately 20 degrees Celsius. B ₁ serves as a first offset to a first linear curve, detailed below, which calculates the valve effective area (A _eff ) of the purge valve 34 .

In step 408 , control corrects B ₁ for changes in temperature as T _amb of the purge valve 34 varies from 20 degrees Celsius. In step 410 , control calculates a first purge valve duty cycle (DC _first ) for purge valve 34 under operating conditions in which DC does not deviate from linearity (e.g., approximately 98%). The DC of the purge valve 34 is the ratio of the “on time,” or time the purge valve 34 is open relative to the time of a single cycle of operation. Generally, DC is expressed as a percentage value. The DC of the purge valve 34 is calculated as follows:
M _purge =FP*M _engair (4)

$\begin{array}{cc}{M}_{\mathrm{purge}}=C_1*A_{\mathrm{eff}}*\delta *\frac{P_1}{T_{\mathrm{amb}}^{1\text{/}2}}& \left(5\right)\end{array}$

In equations (4) and (5), A _eff represents the effective flow area across the purge valve 34 . The value of A _eff is calculated for each distinct actuator of the purge valve 34 that may be used in the fuel system 14 . A _eff is a normalized, dimensionless value which ranges from 0 to 1. A _eff can be characterized by two linear curves that are based on the operation of the purge valve 34 , specifically the duty cycle (DC) of the purge valve 34 .

M _purge represents the purge mass air flow of the fuel system 14 . C ₁ is a constant that considers factors including, but not limited to, the area of the purge valve 34 , discharge coefficients of the purge system, and/or thermodynamics characteristics associated with the purge gas. Each purge valve 34 corresponds to a distinct value of C ₁ . Therefore, if the purge valve 34 implemented in the fuel system 14 is changed, the value of C ₁ will also be modified. In the present implementation, C ₁ can have a value of 5.68. In various embodiments in which a 71 LPM purge valve 34 is contemplated, C ₁ can be 4.12. Combining equations (4) and (5) and isolating A _eff leads to:

$\begin{array}{cc}{A}_{\mathrm{eff}}=\mathrm{FP}*\frac{M_{\mathrm{engair}}}{C_1*\delta *P_1*{\left(\frac{1}{T_{\mathrm{amb}}}\right)}^{\frac{1}{2}}}& \left(6\right)\end{array}$

Where M ₁ represents a slope characterizing the purge valve 34 while operating below a predetermined DC threshold, which is discussed in further detail below. Generally, M ₁ is independent of operating conditions of the purge valve 34 that contribute to the variance of B ₁ . In various embodiments of the current invention, the frequency of operation of the purge valve 34 includes, but is not limited to, 8 Hz, 16 Hz, and/or 32 Hz.
A _eff =M ₁ *DC _first +B ₁ (7)

Combining equations (6) and (7) and isolating DC leads to:

$\begin{array}{cc}D{C}_{\mathrm{first}}=\left(\frac{\mathrm{FP}}{M_1}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_1}{{\left(T_{\mathrm{amb}}\right)}^{\frac{1}{2}}}\right)}-\frac{B_1}{M_1}& \left(8\right)\end{array}$
In step 412 , the method 400 ends.

Referring now to FIG. 5, a method 500 for determining transient DC compensation will be discussed in more detail. Control begins the method 500 in step 502 . In step 504 , control determines whether the DC _first calculated in step 410 of method 400 exceeds the DC threshold. In various embodiments including, but not limited to, the 95 LPM and 71 LPM purge valves 60 , the DC threshold may be 98%. If DC _first exceeds the DC threshold, control proceeds to step 506 . In step 506 , control calculates a compensated purge valve 34 slope, M ₂ , of a second linear curve that characterizes the purge valve 34 when DC _first exceeds the DC threshold.

Next, control calculates a compensated offset B ₂ for the second linear curve in step 508 . M ₂ and B ₂ are determined as functions of M ₁ and B ₁ , respectively. The following equation represents the second linear curve:
A _{comp eff} =M ₂ *DC _comp +B ₂ (9)
In step 510 , DC _comp is calculated in similar fashion to DC _first according to:

$\begin{array}{cc}D{C}_{\mathrm{comp}}-\left(\frac{\mathrm{FP}}{M_2}\right)*\frac{M_{\mathrm{engair}}}{\left(C_1*\delta *\frac{P_1}{{\left(T_{\mathrm{amb}}\right)}^{\frac{1}{2}}}\right)}-\frac{B_2}{M_2}& \left(10\right)\end{array}$
Upon calculating DC _comp , DC _first is set equal to DC _comp and control proceeds to step 512 . If DC _first did not exceed the DC threshold in step 504 , control proceeds to step 512 . In step 512 , control determines the absolute value of the difference between DC _first and DC _second . If the difference exceeds a predetermined constant, Z, control proceeds to step 514 . In step 514 , control determines a delivered DC (DC _deliv ) as a function of DC _first and DC _second . The DC _deliv represents the output DC signal commanding the operation of the purge valve 34 . If the difference between DC _first and DC _second does not exceed the predetermined constant Z in step 512 , control proceeds to step 516 . In step 516 , control sets the DC _deliv to DC _first . In step 518 , control sets DC _second equal to DC _first and transmits DC _second to 304 of method 300 . In step 520 , method 500 ends.

Now referring to FIGS. 6A and 6B, a method 600 , for initializing purge control is discussed in more detail. Control begins the method 600 in step 602 when the engine is turned on. In step 604 , control captures the purge initialization variables including the engine off time (t _engoff ), engine airflow (M _engair ), LTM values, and PLM values. In an exemplary embodiment, the LTM and PLM can include, but are not limited to, 8 cells. LTM cells are used to determine the long-term closed loop fuel multiplier which controls the fuel injectors of the fuel injection system 20 when the purge valve 34 is off. PLM values are utilized to control the fuel injectors when the purge valve 34 is on. The value of t _engoff represents the lapse of time between the previous the previous key down of the engine 12 and the current key up of the engine 12 that occurred in step 602 .

In step 606 , control determines whether the purge valve 34 control is enabled. If the purge valve 34 control is not enabled, control proceeds to step 608 . In step 608 , control activates the LTM cell (LTM _n ) corresponding to the M _engair measurement provided by the M _engair sensor 36 . The LTM, PLM, and FP are indexed by M _engair . In an exemplary embodiment, the LTM, PLM, and FP indexing is calibrated identically. Control returns to step 604 . If the purge valve 34 control is not enabled, control proceeds to 610 . In 610 , control determines whether t _engoff exceeds a predetermined t _engoff threshold. t _engoff can be, but is not limited to, a time interval equaling approximately 1 hour. If t _engoff does not exceed the predetermined t _engoff threshold, then control activates the appropriate cell of the PLM based on the M _engair value in step 612 . In step 612 , control uses values of the PLM that were obtained during a previous key cycle. Control then advances to step 740 of method 700 .

If t _engoff exceeds the predetermined t _engoff threshold in step 610 , the control proceeds to step 614 . In step 614 , control sets the cells of the PLM (PLM _n ) equal to the corresponding cells of the LTM (e.g., PLM ₁ =LTM ₁ ). In step 616 , control activates the PLM cell (PLM _act ) based on the M _engair value supplied by the M _engair sensor 36 . In various embodiments, control can reference a look-up table indexed as a function of M _engair .

In step 618 , control clears a purge flow value (x) to avoid substantial “rich” or “lean” A/F mixture excursions. An A/F ratio below 14.7 is referred to as a rich mixture while an A/F ratio above 14.7 is called a lean mixture. An A/F mixture of 14.7 corresponds to the stoichiometric or chemically correct mixture of gasoline. In step 620 , control increments the purge flow value and then proceeds to calculate actual purge flow percentage (FP _act ) in step 622 . In various embodiments, the FP _act can be computed according to the following equation:
FP _act =0.1* x ² (11)
The previous equation progressively increments FP _act in order to effectuate a gradual increase in the level of evaporated gas supplied through the purge valve 34 . In step 624 , control initiates a hold of the method 600 to allow for an accurate reading to be taken of the closed loop integrator (STM) in step 626 . The STM is a multiplier on the fuel rate that determines a switch of fuel system 14 operating conditions from a rich A/F ratio to a lean A/F ratio and vice versa based on the oxygen sensor 32 . The oxygen sensor 32 is generally disposed across the exhaust manifold 24 of engine 12 . During rich operation of the fuel system 14 , the oxygen sensor 32 transmits a first oxygen sensor voltage level to control. At lean operation, the oxygen sensor 32 signals a second oxygen sensor voltage.

In step 628 , control determines the desired purge flow percentage (FP _des ) based on M _engair . In step 630 , control determines whether the STM value exceeds a predetermined first STM threshold (Y). If the STM value does not exceed Y, the control proceeds to step 634 . If the STM value does exceed Y, then control determines whether FP _des exceeds FP _act in step 632 . If FP _des exceeds FP _act , control returns to step 620 . However, if FP _des does not exceed FP _act , control initiates a second hold of the method 600 to obtain a PLM response in step 634 after adaptively learning the STM correction into the PLM. In step 636 , control performs a second reading of the STM.

Control determines whether the second reading of the STM exceeds a second STM threshold (Z) in step 638 . If the second reading of the STM does not exceed Z, control returns to step 634 . If the second reading of the STM exceeds Z in step 640 , control calculates an adjusted PLM flow percentage (PLM _adj ) for current active cell of the PLM according to:

$\begin{array}{cc}\left({\mathrm{PLM}}_{\mathrm{adj}}\right)=\left({\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_{\mathrm{act}}\right)*\left(\frac{{\mathrm{FP}}_{\mathrm{des}}}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_{\mathrm{act}}& \left(12\right)\end{array}$
where LTM _act represents the value stored in the LTM corresponding to the current active cell. In step 642 , the method 600 ends.

Referring now to FIG. 7, a method 700 for vapor control of the PLM is discussed in more detail. Control begins the method 700 in step 702 . Next, in step 704 , control determines whether PLM _adj is below a PLM flow percentage limit (PLM _lim ). If the PLM _adj is below the PLM _lim , control proceeds to step 706 . In step 706 control calculates a limited PLM flow percentage (FP _lim ) according to:

$\begin{array}{cc}\left({\mathrm{FP}}_{\lim }\right)=\frac{\left({\mathrm{PLM}}_{\lim }-{\mathrm{LTM}}_{\mathrm{act}}\right)}{\left({\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_{\mathrm{act}}\right)}*{\mathrm{FP}}_{\mathrm{act}}& \left(13\right)\end{array}$
In step 708 , control sets a PLM control flag (flag _PLM ) and proceeds in step 710 to calculate the PLM _lim value for each cell of the PLM according to:

$\begin{array}{cc}\left({\mathrm{PLM}}_i\right)=\left(\left(\frac{{\mathrm{LTM}}_i}{{\mathrm{LTM}}_{\mathrm{act}}}\right)*{\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_i\right)*\left(\frac{{\mathrm{FP}}_{\lim }}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_i& \left(14\right)\end{array}$
where i corresponds to the current cell position reference. In step 712 , control sets the FP _act equal to the FP _lim .

In step 714 , control initiates a hold of the method 700 to obtain a PLM response. In step 716 , control determines whether the STM value exceeds Z. If the STM value does not exceed Z, control returns to step 714 . If the STM value exceeds Z, control determines the appropriate FP cell based on the M _engair in step 718 . In step 720 , control proceeds to determine whether FP _act is below FP _des . If FP _act is below FP _des , control proceeds to step 722 . In step 722 , control determines whether PLM _act is below PLM _lim . If PLM _act is below PLM _lim , control proceeds to step 724 . In step 724 , control decreases FP _act by a factor of Y. If PLM _act exceeds PLM _lim , control proceeds to step 726 . In step 726 , control increases FP _act by a factor of Y. In step 728 , control signals a second hold of the method 700 to obtain a PLM response and the proceeds to step 730 .

Referring back to step 720 , if the FP _act exceeds the FP _des , control proceeds to step 734 . In step 734 , control clears the flag _PLM and proceeds to step 736 . In step 736 , control activates a FP _act based on the FP _des . In step 730 , control determines whether the absolute STM value exceeds Z. If the value of the STM does not exceed Z, control returns to step 728 . If the value of the STM does exceed Z, control proceeds to step 732 .

Referring back to step 704 , if the PLM _adj exceeds the PLM _lim , control proceeds to step 732 . In step 732 , control computes the PLM value for each cell of the PLM according to:

$\begin{array}{cc}\left({\mathrm{PLM}}_i\right)=\left(\left(\frac{{\mathrm{LTM}}_i}{{\mathrm{LTM}}_{\mathrm{act}}}\right)*{\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_i\right)*\left(\frac{{\mathrm{FP}}_i}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_i& \left(15\right)\end{array}$
where i corresponds to the cell position within the PLM. In step 738 , control determines whether the flag _PLM is set. If the flag _PLM is set, control returns to step 718 . If the flag _PLM is not set, control proceeds to step 740 . In step 740 , control commands the vapor control based on the FP _des . The method 700 ends in step 742 .

Referring now to FIG. 8, a method 800 for vapor control during periods of 100% DC will be discussed in more detail. Control begins the method 800 in step 802 . In step 804 , control determines the LTM and PLM control variables including, but not limited to, M _engair , DC _new , the ignition key position (key _pos ), δ, P ₁ , and t _baro . In step 806 , control determines whether the purge valve 34 control is enabled. If the purge valve 34 control is not enabled, control proceeds to 808 . In step 808 , control activates the LTM cell (LTM _act ) corresponding to the M _engair supplied by the M _engair sensor 36 after which, control advances to step 822 .

If control determines that the purge valve 34 control is enabled, control proceeds to step 810 . In step 810 , control activates the PLM cell (PLM _act ) corresponding to the M _engair . In step 812 , control determines whether DC _deliv exceeds 100%. If DC _deliv is below 100%, then control proceeds to step 822 . If DC _deliv exceeds 100%, control proceeds to step 814 . In step 814 , control computes a flow corrected P ₁ (P _corr ) for the FP _act according to:

$\begin{array}{cc}\left(P_{\mathrm{corr}}\right)=P_{\mathrm{baro}}-\frac{\left(100*M_2-B_2\right)}{\left({\mathrm{DC}}_{\mathrm{new}}*M_2-B_2\right)}*\left(P_{\mathrm{baro}}-P_1\right)& \left(16\right)\end{array}$

In step 816 , control recalculates δ with the P _corr supplied by step 814 . In step 818 , control determines the FP _act experienced by the purge valve 34 at 100% DC according to the following equation:

$\begin{array}{cc}\left({\mathrm{FP}}_{100DC}\right)=\frac{\left(\frac{C_1*\delta *P_1}{T_{\mathrm{amb}}^{\frac{1}{2}}}\right)}{M_{\mathrm{eng}}}& \left(17\right)\end{array}$

In step 820 , control computes the adjusted PLM (PLM _100DC ) while the purge valve 34 is operating at 100% DC according to the following equation:

$\begin{array}{cc}\left({\mathrm{PLM}}_{100DC}\right)=\left({\mathrm{PLM}}_{\mathrm{act}}-{\mathrm{LTM}}_{\mathrm{act}}\right)*\left(\frac{{\mathrm{FP}}_{100DC}}{{\mathrm{FP}}_{\mathrm{act}}}\right)+{\mathrm{LTM}}_{\mathrm{act}}& \left(18\right)\end{array}$
In step 822 , control determines whether the key _pos is in an off state. If control determines that the key _pos is not set to off, control returns to step 804 . If the key _pos is set to off, control proceeds to step 824 . In step 824 , control initializes the cell position reference, i, to zero. In step 826 , control increments the cell position reference. In step 828 , control compares the LTM value to the PLM value for each corresponding cell position. If the value of the LTM is equivalent to the PLM for the current cell, control returns to 826 . If the value of the LTM is not equivalent to the PLM for the current cell, control proceeds to step in 830 . In step 830 , control calculates and stores the LTM cell value based on the corresponding PLM cell according to the following equation:

$\begin{array}{cc}\left({\mathrm{LTM}}_i\right)=\frac{{\mathrm{PLM}}_i}{\left(\frac{{\mathrm{PLM}}_{\mathrm{act}}}{{\mathrm{LTM}}_{\mathrm{act}}}-1\right)*\left(\frac{{\mathrm{FP}}_i}{{\mathrm{FP}}_{\mathrm{act}}}\right)+1}& \left(19\right)\end{array}$

In step 832 , control determines whether i equals the last reference position of the PLM and LTM. If i does not equal the last position of the PLM and LTM, control returns to step 826 . If i does equal the last position of the PLM and LTM, control proceeds to step 834 . In step 834 , control ends.

Those skilled in the art can now appreciate from the foregoing description that the broad teachings of the present invention can be implemented in a variety of forms. Therefore, while this invention has been described in connection with particular examples thereof, the true scope of the invention should not be so limited since other modifications will become apparent to the skilled practitioner upon a study of the drawings, the specification and the following claims.