Available online at www.sciencedirect.com Procedia Environmental Sciences 11 (2011) 996 1000 Applying Matlab/Simulink to Study Calculation of NO x Efficiency of the SCR Yuhang Su, Jun Li, Ying Gao and Dawei Qu State key laboratory of automobile dynamic simulation Jilin University Changchun, China suyh09@mails.jlu.edu.cn, junli610@263.net Abstract An SCR system offers the opportunity for a car to get its nitrogen oxides inside the exhaust gas decreased. This can be achieved by dosing a reducing agent into the system which selectively reacts with the oxygen of the nitrogen oxides and therefore causes these emissions to be diminished. Depending on the operating point of an engine more or less NO x emissions can be reduced. Adaptation logic has the ability to find the optimum operating point at the NH 3 slip limit. The adaptation logic requires the NO x emissions downstream from the SCR which will be measured by a NO x -Sensor. The goal of the module SCRChk_NO x EffCalc is the monitoring of NO x conversion efficiency of an SCR system. The monitoring function is based on a NO x sensor downstream from the SCR and the known NO x emissions upstream from the SCR. 2011 Published by Elsevier Ltd. Open access under CC BY-NC-ND license. Selection and/or peer-review under responsibility of the Intelligent Information Technology Application Research Association. Keywords: SCR; OBD; NO x ; efficiency 1. Introduction The exhaust gas legislations of several countries regulate the maximum NOx emissions on the road and the requirements for monitoring by On Board Diagnosis (OBD). The OBD legislation regulates the monitoring for catalysts, misfire, fuel system, exhaust gas recirculation, etc. The different legislations of the countries have diverse requirements for OBD. The CARB OBDII is a very strict legislation and has the most extensive requirements. The European OBD (EOBD) is partially similar and the legislation also defines NOx emissions limits. A main requirement of all relevant exhaust gas legislations is the NOx emissions monitoring with a defined NOx emissions limit. The NOx emissions limits are always related to defined test cycles (e.g. FTP, etc.). A challenge is that the monitoring shall be active under in use conditions and not only under test cycle conditions. The goal of the module SCRChk_NOxEffCalc is the monitoring of NOx conversion efficiency of an SCR system. The monitoring function is based on a NOx sensor downstream from the SCR and the 1878-0296 2011 Published by Elsevier Ltd. Open access under CC BY-NC-ND license. Selection and/or peer-review under responsibility of the Intelligent Information Technology Application Research Association. doi:10.1016/j.proenv.2011.12.152
Yuhang Su et al. / Procedia Environmental Sciences 11 (2011) 996 1000 997 known NOx emissions upstream from the SCR. The main coordination (enabling/reset) and input signal filtering of the module SCRChk_NOxEffCalc will be done in the module SCRChk_NOxEffCo. [DStgy_swtDStgyOff] [SCRChk_dmExhFlt ] [SCRChk_dmNOxSensDsFlt ] [SCRChk_dmNOxUsFltDdT ] [SCRChk_etaEstFlt ] [SCRChk_rNOxUsFltDdT ] [SCRChk_stEtaEna ] [SCRChk_stEtaEnaPreCtl ] [SCRChk_stEtaRst ] [SCRChk_stEtaRstPreCtl ] [SCRT_tFld] SCRChk_NOxEffCalc Figure 1 Overview 2. Function The module "SCRChk_NOxEffCalc" requires the following input signals: Main enabling range SCRChk_stEtaEna; Reset signal range SCRChk_stEtaRst; SCR catalyst temperatures SCRT_tFld and other available temperatures in the class SCRT_TempSel; Filtered NOx concentration upstream from the SCR SCRChk_rNOxUsFltDdT; [SCRChk_swtSetMskRng2_C] [SCRChk_stEtaEna] [SCRChk_stEtaRst] In3 [SCRChk_stEtaEnaPreCtl] In4 Out2 [SCRChk_stEtaRstPreCtl] In5 Enable Range Settings [SCRChk_numFldTempSCR_C] SCRChk_numPosTempSCR_C] In3 [SCRChk_dmExhFlt] [SCRChk_rNOxUsFltDdT] SCRT_TempSel In4 In5 [SCRChk_stEtaCalcFld] [SCRChk_dmNOxUsFltDdT] In6 [SCRChk_dmNOxSensDsFlt] In7 Out2 [SCRChk_etaEstFlt] In8 Calculation Eta Range Bit Setting Figure 2 Overview of the module SCRChk_No x EffCalc Filtered exhaust gas mass flow SCRChk_dmExhFlt; Filtered NOx mass flow upstream from the SCR SCRChk_dmNOxUsFltDdT; Filtered NOx mass flow downstream from the SCR SCRChk_dmNOxSensDsFlt; Filtered estimated SCR efficiency range SCRChk_etaEstFlt; Currently effectively dosed ammonium equivalent mass flow SCRMod_dmNH3CurrEffSlw; Required ammonium feed rate for complete NOx reduction SCRFFC_dmNH3FdRat;
998 Yuhang Su et al. / Procedia Environmental Sciences 11 (2011) 996 1000 Author name / Procedia Environmental Sciences 00 (2011) 000 000 2.1Subfunction: "Calculation Eta Range" stetastm [stetaena ] [stetarst ] [tscr] stetaena stetarst tscr stetaena 1 stetarst 1 stetaena stetastm 3 new efficiency calculation == [rnoxusfltddt ] rnoxusfltddt stetarst [dmexhflt ] dmexhflt stcond stcond dmnoxusavrg dmnoxusavrg [dmnoxusfltddt ] [dmnoxsensdsflt ] dmnoxusfltddt Operating Point Selection dmnoxus dmnoxds tinoxitgavrg stetacalcvld tinoxitgavrg AND [stetacalcnew] [etaestflt ] dmexhflt tscr etathres dmnoxusfltddt etaestflt Calculation Threshold etathres ctcondup ctconddwn etaactavrg etathresavrg Validity Check Eta Calculation [etaactavrg ] [etathresavrg ] SCRChk_NOxEffCalc_Eta Figure3 Subfunction "Calculation Eta Range" The subfunction "Calculation Eta Range" calculates the actual average efficiency SCRChk_etaActAvrg_mp and the average efficiency threshold SCRChk_etaThresAvrg_mp. The efficiency is not calculated continuously but discretely. The NOx mass flows up/downstream from the SCR are integrated separately and after a defined NOx mass flowing through the SCR system the efficiency will be calculated. The efficiency calculation is coordinated by a state machine with the measuring point SCRChk_stEtaStM_mp. The efficiency calculation is made discretely and the status variable SCRChk_stEtaCalcNew_mp shows that a new result has been calculated. The input signals are the enabling status SCRChk_stEtaEna, reset status SCRChk_stEtaRst, SCR catalyst temperature SCRChk_tSCR_mp, exhaust gas mass flow SCRChk_dmExhFlt, NOx concentration upstream from the SCR SCRChk_rNOxUsFltDdT, NOx mass flow upstream from the SCR SCRChk_dmNOxUsFltDdT, NOx mass flow downstream from the SCR SCRChk_dmNOxSensDsFlt and estimated efficiency SCRChk_etaEstFlt. SCRChk_NOxEffCalc_Eta calculates the actual average efficiency SCRChk_etaActAvrg_mp from the NOx mass flow upstream from the SCR SCRChk_dmNOxUsFltDdT and downstream from the SCR SCRChk_dmNOxSensDsFlt. The actual average efficiency SCRChk_etaActAvrg_mp is defined as SCRChk_etaActAvrg_mp = 1 SCRChk_facEtaCor_C (1) SCRChk_dmNOxSensDsFlt_mp SCRChk_dmNOxUsFltDdT_mp The correction factor SCRChk_facEtaCor_C allows e.g. the zero adjustment of the actual efficiency without dosing. If the efficiency calculation is executed, the counter value SCRChk_ctEtaCalc_mp is raised by ONE. The two mass flows SCRChk_dmNOxUsFltDdT and SCRChk_dmNOxSensDsFlt are only integrated during special applicable conditions. The subfunction "Operating Point Selection" checks the conditions for integration and if all conditions are fulfilled the status variable SCRChk_stCond_ mp is TRUE. The status SCRChk_stCond_mp influences an up/down condition counter SCRChk_ctCond_mp. The
Yuhang Su et al. / Procedia Environmental Sciences 11 (2011) 996 1000 999 condition counter is the debouncing for controlling the integration. If the status SCRChk_stCond_mp is TRUE the up/down condition counter will be incremented by SCRChk_ctCondUp_mp otherwise decremented by SCRChk_ctCondDwn_mp. The condition counter is limited between zero and SCRChk_ctCondMax_C and starts at SCRChk_ctCondStrt_C. The increment value SCRChk_ctCondUp_mp depends on the curve SCRChk_ctCondUp_CUR dependent on the NOx mass flow SCRChk_dmNOxUsFltDdT. The decrement value SCRChk_ctCondDwn_mp depends on the curve SCRChk_ctCondDwn_CUR dependent on the NOx concentration SCRChk_rNOxUsFltDdT. The condition counter will be set imme- diately to zero in case of disabling (SCRChk_stEtaEna_mp =FALSE) or reset (SCRChk_stEtaRst_mp=TRUE). The integration (1) of the NOx mass flows SCRChk_dmNOxSensDsFlt and SCRChk_dmNOxUsFltDdT will be done with a preintegration and a main integration. Figure 4 Condition Counter At first the NOx mass flows SCRChk_dmNOxSensDsFlt and SCRChk_dmNOxUsFltDdT are integrated with two preintegrators when the condition counter starts counting. These preintegrators will be set to zero if the condition counter reaches zero. If the preintegration time reaches SCRChk_tiPreInteg_C the integration values of the pre integrators will be added to the main integrators and the main integrators continue the integration of both NOx mass flows SCRChk_dmNOxSensDsFlt and SCRChk_dmNOxUsFltDdT. The main integrators stop when the condition counter reaches zero or when the conditions of efficiency calculation are fulfilled. The main integrators will only be reset after an efficiency calculation or a reset (SCRChk_stEtaRst_ mp=true). The disabling (SCRChk_stEtaEna_mp=FALSE) only stops the main integrators but does not reset the main integrators. If the preintegration time does not reach the value SCRChk_tiPreInteg_C the integrated NOx mass flows during the preintegrations will not be taken into account for the integration (1). A reset request by SCRChk_stEtaRst_mp resets the condition counter, all pre and main integrators. The average actual efficiency SCRChk_etaActAvrg_mp (1) will be calculated after the integrated NOx mass upstream from the SCR SCRChk_mNOxUs_mp reaches the minimum NOx mass SCRChk_mNOxUsMin_C and the condition counter is zero. The efficiency can only be calculated directly after the end of a main integration. If the condition counter does not reach zero, the efficiency will be calculated after the delay time SCRChk_tiDlyMaxCalcEta_ C of reaching the NOx mass SCRChk_mNOxUsMin_C. The delay timer will be reset if the main integration is not active. The average efficiency threshold SCRChk_etaThresAvrg_mp is calculated simultaneous with the ac-
1000 Yuhang Su et al. / Procedia Environmental Sciences 11 (2011) 996 1000 Author name / Procedia Environmental Sciences 00 (2011) 000 000 tual efficiency calculation (1) and use the same condition counter. SCRChk_eta ThresAvrg_mp ( 1 SCRChk_eta Thres_mp ) SCRChk_dmN OxUsFltDdT _mp = 1 SCRChk_dmN OxUsFltDdT _mp The subfunction "Calculation Threshold" calculates the efficiency threshold SCRChk_etaThres_mp. The subfunction "Validity Check Eta Calculation" validates the operating points during the integration and the duration of the integration after the calculation of a new efficiency. 2.2Subfunction: "Bit Setting" Table 1 Table Type Styles (2) Bit Number Description 0 SCRChk_stCond_mp 1 SCRChk_stEtaCalcNew_mp The subfunction "Bit Setting" collects bit information in the bit field SCRChk _stetacalcfld. The status information of new test results SCRChk_stEtaCalcNew_mp (bit 1) are only pulses. The pulses are converted in a toggled signal in SCRChk_stEtaCalcFld. Each pulse toggles the signal. References [1] Johnson, T.V., Diesel Emission Control in Review SAE 2009-01-0121, 2009. [2] Nebergall, Hagen and Owen, Selective catalytic reduction on-board diagnostics: past and future challenges, SAE 2005-01-3603, 2005. [3] Schär, Onder, Elsener and Geering, Model-based control of an SCR system for a mobile application,presented at SAE World Congress, 2004-05-0412, 2004. [4] Chi and DaCosta, Modeling and control of a urea-scr aftertreatment system, presented at SAE World Congress, 2005-01-0966, 2005 [5] Willems, Cloudt, Van den Eijnden, Van Genderen and Verbeek, De Jager, Boomsma, Is closed-loop SCR control required to meet future emission targets?, presented at SAE World Congress, 2007-01-1574, 2007. [6] Andrew Herman, Ming-Cheng Wu, David Cabush and Mark Shost, Model Based Control of SCR Dosing and OBD Strategies with Feedback from NH 3 Sensors, SAE 2009-01-0911, 2009.