Reduced Order Model Design for Three Way Catalytic Converters
|
|
- Jasmin Fleming
- 5 years ago
- Views:
Transcription
1 Clemson University TigerPrints All Theses Theses Reduced Order Model Design for Three Way Catalytic Converters Romit Ravishankar Godi Clemson University, Follow this and additional works at: Recommended Citation Godi, Romit Ravishankar, "Reduced Order Model Design for Three Way Catalytic Converters" (2017). All Theses This Thesis is brought to you for free and open access by the Theses at TigerPrints. It has been accepted for inclusion in All Theses by an authorized administrator of TigerPrints. For more information, please contact
2 REDUCED ORDER MODEL DESIGN FOR THREE WAY CATALYTIC CONVERTERS A Thesis Presented to the Graduate School of Clemson University In Partial Fulfillment of the Requirements for the Degree Master of Science Electrical Engineering by Romit Ravishankar Godi May 2017 Accepted by: Dr. Rajendra Singh, Committee Chair Dr. Simona Onori Dr. Richard Groff
3 Abstract This research topic aims to design and simulate a reduced order model of a physics based three way catalyst (TWC) in automotive systems. A catalytic converter carries out a chemical redox reactions of pollutant gases in order to convert them to less harmful gases before emission through the vehicle exhaust. The thermal and oxygen storage dynamics of a TWC are described using a set of partial differential equations (PDE) which are coupled and non-linear in nature along the spatial axis of the catalyst. Computing the solution for a PDE based system of equations is arduous since the state of the system is distributed in an infinite domain along the spatial axis. Hence a model order reduction needs to be carried out in order to observe the behavior of PDE based systems using only a finite number of states. In order to solve for such systems, the PDEs need to be transformed into ordinary differential equations (ODE) which result in a finite, but still a large, number of states. This transformation from a set of PDEs with infinite states to an ODE system with a finite number of states is carried out using the Galerkin s projection method. Since we are dealing with a non-linear PDE, the resulting ODEs are time dependent non-linear system of equations. The model order reduction of the ODE based system is achieved using Proper Orthogonal Decomposition (POD). The POD method uses a basis representation of the system variable, and these basis functions are derived from a matrix containing the values of the variable at different instants over time and space obtained through real time measurement or simulation. The basis functions are extracted from the data matrix by performing Singular Value Decomposition (SVD). SVD computes the eigen values of the system and stores them in a matrix arranging them from the most dominant modes to the least. The number of eigen values selected represents the number of basis functions for the system. An energy estimation criterion is selected to limit the number of basis functions that will capture the entire system dynamics. Therefore we end up with a reduced number of modes that describe the order of the ODE based system. The accuracy of the reduced order model will be calculated for different number of basis functions and a selection of an optimal model order will be determined by computing the root-mean-square error ii
4 (RMSE) for each of these modes. The reduced model will be developed over the Federal Test Protocol (FTP) driving cycle and the results will be validated alongside the measurement values over the US06 driving cycle.
5 Dedication I would like to dedicate my thesis to my beloved grandparents iv
6 Acknowledgments This thesis work is a result of the help and support provided by people, departments, and organizations towards a successful completion and also in my growth as an individual both professionally and personally. I would like to thank all my committee members, starting with my committee chair Dr. Rajendra Singh who was instrumental in helping me find the research project in my final year of masters and overseeing my work by providing vital guidelines. A big thank you to my research advisor Dr. Simona Onori for being immensely helpful and patient and with whom I performed the entire research work. I would also like to thank Dr. Richard Groff for providing vital inputs through the course of my work. I greatly acknowledge the contribution made towards this research by Fiat Chrysler Automobiles LLC, Auburn Hills, MI for providing me with the measurement data required for developing the model for my project. Special thanks go to the visiting scholars Simone Gelmini from Politecnico di Milano and Umberto Sartori from University of Trento for helping me acquire data from the catalyst and identify parameters for the model equations. Most part of my thesis work was written at the Clemson University-ICAR campus in Greenville and I appreciate the help from the administration for providing me with all the resources and an office space to complete my work. I would also like to thank the Electrical and Computer Engineering Department for supporting me financially throughout the duration of the project by providing me with a teaching assistant position. Pursuing a master s degree in a different part of the world and being thousands of miles away from home is a challenging and financially difficult decision to make. I would like to thank my parents who never v
7 once doubted my decision and supported me in all ways to make this happen. Last but not least, my sincerest gratitude to all my friends who constantly motivated and supported me through the tough times and were like my family away from home. Without their influence in my life, this journey would not have been as adventurous and exciting.
8 Nomenclature TWC Three Way Catalyst ODE Ordinary Differential Equation PDE Partial Differential Equation OSC Oxygen Storage Capacity λ Normalized air fuel ratio [ mol ] m 3 t Time [s] z Spacial coordinate [m] T g Gas temperature [K] T cat TWC solid phase temperature [K] T exh Exhaust gas temperature [K] T cat0 Initial temperature [[K] ] Exhaust gas mass flow rate kg ṁ exh ρ g Exhaust gas density [ s ] kg m 3 L TWC length [m] A out TWC external surface area [m 2 ] V cat TWC volume [m 3 ] i Reactions index [ H i Reaction enthalpy difference J ] [ mol R i Reaction rate for the i th reaction mol ] [ m Q 3 s react Heat produced by reactions W ] m 3 η TWC efficiency [ c 0 Exhaust gas concentration mol ] m 3 k f i Forward reaction rate ki b Backward reaction rate P Exhaust gas pressure [[P a] ] R Universal gas constant (R = 8.314) J kg K A i Arrenhius pre-exponential factor E i Activation energy [J] K i Chemical equilibrium constant [ s v Space velocity m ] s ] Average molar mass M exh Φ State of oxygen [ kg mol vii
9 Table of Contents Title Page Abstract i ii Dedication iv Acknowledgments v Nomenclature vii 1 Introduction Emission standards and regulations Three Way Catalyst Motivation Three way catalyst model Thermal model Oxygen Storage Model Implementation Model Reduction Proper orthogonal decomposition (POD) Galerkin's Projection TWC Reduced Order Model Design Singular Value Decomposition (SVD) Applying Galerkin's projection Simulation Results Conclusions and Future Work Bibliography viii
10 Chapter 1 Introduction Over the years the need for better regulation of emission control has seen a lot of technological advancements in the automotive industry, with a focus on developing efficient after-treatment systems [1]. The US Environment Protection Agency(EPA) is the regulatory body that sets the emission limit standards for all on-road vehicles. The most common and harmful gases include carbon monoxide(co), nitrogen oxides(no x ), and hydrocarbons(hc). All automotive manufacturers must adhere to these standards. 1.1 Emission standards and regulations Figure 1.1 shows the compliance life cycle of a light duty vehicle. The compliance certification reflects the standards set for tailpipe emissions after having run the vehicle over different driving cycles. Some of the driving cycles used by the EPA are the Federal Test Procedure(FTP), Federal Urban Driving Schedule (FUDS), and US06. As far as emission standards are considered in the United States, the EPA has been following a set regulations[2]. 1
11 Figure 1.1: Compliance Life of a Light-Duty Vehicle (reproduced from [1]) Tier 1 standards These are standards that were implemented between 1994 to 1997 and it applies to all vehicles with gross vehicle weight rating(gvwr) below 8,500 pounds, [3]. The vehicles in this category are also called light-duty vehicles(ldv). The emission standards for tier 1 were measured over FTP cycle shown in Figure 1.2.
12 Figure 1.2: Tier 1 emission standard (reproduced from [3]) Tier 2 standards Tier 2 standards were phased between 2004 to 2009 and it includes medium-duty passenger vehicles(mdpv) with GVWR between 8,500 to 10,000 pounds. These standards were similar to tier 1 for LDVs but had stringent limits for heavier vehicles.
13 Figure 1.3: Tier 2 emission standard (reproduced from [2]) Tier 2 introduced a new emission standard regulation called bins wherein the bins were numbered from 1 to 11, with 1 being the cleanest and 11 the pollutant. Each bin has a pollutant gas emission level set and vehicle manufacturers had to chose which bin to classify their vehicle under. Figure 1.3 gives the emission standard over the FTP cycle Tier 3 standards The latest standards will now be decided upon to be implemented between 2017 to 2025 and will incorporate heavy-duty vehicles with GVWR upto 14,000 pounds. Tier 3 standards also use the concept of bins but they are categorized into 7 bins based on their non-methane organic gas(nmog) and nitrogen oxides(no x ) concentrations. Figure 1.4 represents the standards over FTP cycle for all vehicles using any type of fuel.
14 Figure 1.4: Tier 3 emission standard (reproduced from [4]) As seen above, the EPA set emission standards state the permissible limits of harmful gases and particulate matter for all types of vehicles. In order to meet these requirements, advances have been greatly made towards developing sophisticated engine and after-treatment systems. As the complexity in the design of automotive engines has increased, gasoline direct injection (GDI) engine manufacturers have been challenged to develop efficient after-treatment and software systems to monitor the operating conditions of various components. One such development includes the thermal oxidation of carbon monoxide and hydrocarbons, whereby air is pumped into the exhaust system to allow for a complete combustion of all gases [5]. Another technique used is the exhaust gas recirculation (EGR) [6] in which the reduction of nitrogen oxides is aided by re-introducing some of the exhaust gas back into the catalyst inlet. This diffuses the excess oxygen in the fuel, thus completing the reduction reactions. All automotive companies now use computer monitored systems on the vehicle called On-board diagnostics (OBDs) [7], which monitors and displays the operational information of exhaust related components such as the catalyst and oxygen sensor. Any malfunction in the sensor readings is brought to the drivers attention through an indicator on the dashboard display. The most widely used emission control strategy involves the use of a catalytic converter and filters to reduce the emission of harmful gases and particulate matter. 1.2 Three Way Catalyst As stated in [8], the principle under which the internal combustion engine operates is that the fuel needs to be oxidized in order to create the pressure to drive the pistons. This pressure is generated from the heat which is released after the chemical reaction occurs between the fuel and air. In an ideal combustion
15 reaction, oxygen (O 2 ) reacts with the hydrocarbons (HC) to produce carbon dioxide (CO 2 ), water (H 2 O) and heat. But since the intake is not pure oxygen but air which consists of O 2 and nitrogen (N 2 ), the reaction with the fuel results in harmful gases such as carbon monoxide (CO), oxides of nitrogen (NO x ), and unburned hydrocarbons. Hence the need arises for a catalytic converter to treat these gases before releasing them through the exhaust. A three way catalytic converter acts as a catalyst for the redox (reduction-oxidation) reactions in which the harmful pollutants are converted to less toxic emissions. It is called a three way catalyst since it aids in the reduction of nitrogen oxides to less harmful nitrogen gas (N 2 ), oxidation of carbon monoxide (CO) to carbon dioxide (CO 2 ), and oxidation of unburnt hydrocarbons (HC) to carbon dioxide (CO 2 ) and water (H 2 O). These chemical reactions occur simultaneously inside the catalyst and the efficiency with which the redox reaction occurs depends highly on the air-fuel ratio (AFR) and the temperature inside the catalyst. The AFR is the ratio of the mass of air to fuel in the mixture inside the internal combustion engine during ignition given as: where AF R = m air m fuel (1.1) m air = mass of air m fuel = mass of fuel (1.2) The AFR has a stoichiometric value, defined for gasoline and ethanol based engines, which has the exact amount of air and fuel to produce chemically complete combustion. The stoichiometric AFR is defined as: 14.7 for gasoline AFR (stoichiometric) = 9 for ethanol (1.3) AFR can be written in a normalized form by taking the ratio of the air-fuel present in the mixture at the time to combustion and the stoichiometric value. This can be written as: λ = = AFR (actual) AFR (stoichiometric) m air /m fuel AFR (stoichiometric) (1.4)
16 100 Conversion efficency [%] NOx CO HC λ [-] (a) (b) Figure 1.5: TWC conversion efficiency vs (a) AFR (reproduced from [10]) and (b) catalyst temperature (reproduced from [11]). As shown in Figure 1.5, if the AFR is within a narrow region around the stoichiometric value, the chemical reactions occur with higher efficiency, i.e. the success of complete reduction and oxidation reactions is higher. For a gasoline engine, the AFR is about 14.7:1, i.e the mixture has 14.7 parts of air and one part fuel [9]. The efficiency decreases when the AFR falls out of this range. The regions outside of the stoichiometric value are classified as lean and rich conditions. When the AFR is less than the stoichiometric value, the combustion occurs in rich conditions (λ < 1). In this case the oxygen is utilized during fuel combustion and there is not sufficient amount left to carry out the oxidation reactions inside the catalyst. When the AFR is greater than the stoichiometric value, the engine operates under lean condition (λ > 1), i.e there is excess oxygen available inside the catalyst and hence the reduction reaction does not occur completely. The catalytic conversion efficiency also depends on the temperature inside the TWC. Only when the catalyst temperature has reached a threshold level, called light-off temperature, the chemical conversions take place with an efficiency close to 50%. Therefore the catalyst temperature dynamics, plays an imperative role in the control of harmful emissions. Hence the ability to monitor and predict the catalyst temperature along the catalyst axis will help design an efficient exhaust system. The present technology does not allow for such kind of measurement in real-time application.
17 1.3 Motivation Common approaches to model development to capture the TWC dynamics rely on two types of modeling, i.e physics based models and empirical or storage based models [12], [13], [14] and [15]. The models presented in [16] and [17], develop comprehensive mathematical models consisting of non-linear PDEs which are not feasible for developing real time control application based systems due their complexity. In [18] and [19] the models are developed using numerical techniques which result in a large number of ODEs due to their discretization step size used for the species concentration and along the flow direction. The comprehensive nature os physics based models also makes it difficult to identify all of the system parameters. When it comes to empirical models as proposed in [20] and [21], neglect some of the chemical kinetics and interdependence of the thermal and oxygen storage models. This results in much simpler models, as compared to the physics based models, which are not consistent across different catalyst operating conditions. Empirical models use curve fitting techniques in order to linearize the non-linear terms in the PDEs [22]. Due to this, the model cannot be applied in real time scenarios and their accuracy is not guaranteed to remain constant across different driving cycles. This thesis proposes a reduced order control-oriented model which takes into account the gas and solid phase catalyst dynamics without being mathematically intensive by neglecting the time dynamics of the gas phase. This results in a one dimensional gas model. For simulating the oxygen storage dynamics this model will not require individual gas concentration level measurements making it computationally feasible and able to implement in real time control applications. The model design uses the POD based Galerkin's projection, which follows a basis representation of the system parameter. These basis functions are obtained by performing singular value decomposition (SVD) on the matrix containing data from the TWC thermal and oxygen storage models. The reduced order model will be developed over the FTP driving cycle and validated over the US06 cycle. To the best of our knowledge, a rigorous TWC control-oriented model design using formal reduction methods applied to an experimentally validated physics-based system has never been investigated. The rest of the thesis is divided into the following chapters. Chapter 2 discusses the physics based PDE model of the TWC. In Chapter 3, the POD-Galerkin's projection for reduced order modeling (ROM) is introduced. In Chapter 4, the ROM is developed for the catalyst system. In Chapter 5, the simulation results of the reduced model are presented. Chapter 6 gives the conclusion and scope for future work.
18 Chapter 2 Three way catalyst model The TWC model used in this work is a physics-based model borrowed from the literature [17], and consists of the thermal and oxygen storage dynamics. The catalyst is usually made up of a substrate which is designed in the shape of a honeycomb structure and covered with a washcoat material. This type of structure allows for maximum surface area for the exhaust gases to pass through. The washcoat consists of the actual catalyst material made up of Rhodium (Rh), Platinum (Pt), and Palladium (Pd). Generally, Rhodium aids the reduction reactions while Platinum and Palladium are efficient for oxidation reactions. As mentioned previously, reduction reaction occurs efficiently when the AFR mixture is rich and oxidation occurs during lean condition. Therefore in the rich case there is a need to maintain lower oxygen levels and during the lean case there needs to be higher levels of oxygen to carry out the reactions. The TWC conversion efficiency depends mainly on the air-fuel ratio (AFR) and the catalyst temperature. Therefore the dynamics of the system can be described using a thermal model and an oxygen storage model. Each of these models consists of a system of PDEs which are non-linear and coupled in nature as in [17]. The thermal model describes the behavior of the gas temperature as well as the catalyst temperature which governs the oxygen storage capacity. The inputs to the thermal model include the exhaust gas temperature, exhaust gas mass flow rate, and the ambient temperature. In case of the oxygen storage model, it is computationally difficult to monitor the concentration levels of all gases entering into the catalyst. Hence the model relies on controlling the oxygen storage levels which helps maintain the AFR near stoichiometry. 9
19 2.1 Thermal model The thermal behavior of the TWC follows a one dimensional approach along the spatial axis, and can be described by carrying out an energy balance on the exhaust gas and on the solid phase. The gas phase comprises of the exhaust gas flow coming from the engine and flowing into the catalyst and the solid phase is represented by the substrate and the washcoat lumped together [23]. The partial differential equations describing these dynamics are given as follows: Energy balance of the gas phase The PDE model describing the energy balance of the gas phase is given by [15] ρ g ɛ c pg T g t = ṁexh A cs c pg T g z + h A geo (T cat T g ) (2.1) T where ρ g ɛc g pg t, represents the time dynamics of the gas phase. ṁ exh T c g pg z, describes the convective heat transport in the axial dimension and ha geo (T cat T g ) is the heat exchange between the gas and the solid phase. Considering the gas phase dynamics to be changing rapidly in comparison to the solid phase, we assume the gas temperature T g to be at steady state. Therefore we can neglect the term ρ g ɛc pg T g t. The simplified equation can be written as: A cs ṁ exh A cs c pg T g z = h A geo (T cat T g ) (2.2) The boundary condition for the above equation is: T g (z = 0, t) = T exh (t) (2.3) where the exhaust temperature T exh, measured at the input of the catalyst, is considered to be the initial gas temperature.
20 Energy balance of the solid phase Similarly, the energy balance of the solid phase is described by: ρ s (1 ɛ) c s T cat t = (1 ɛ) λ s 2 T cat z 2 A out V cat h out (T cat T amb ) h A geo (T cat T g ) + Q react (2.4) where the term ρ s (1 ɛ)c s T cat t, represents the time dynamics of the solid phase. (1 ɛ)λ 2 T cat s z, accounts 2 for the conduction losses in the substrate. ha geo (T cat T g ) represents the heat exchange between the gas and the solid phase. Q react is the heat produced by the exothermic reactions in the washcoat. Finally the term Aout V cat h out (T cat T amb ), accounts for the radial losses in the ambient phase The term Q react describes the heat dynamics due to the chemical reactions occurring in the catalyst and is a function of the conversion efficiency. Its dependency on the exhaust mass flow rate is true only under stoichiometric conditions [15]. It can be expressed as follows: Q react = K react ṁ exh η(t cat ) (2.5) The contribution of the efficiency term in Equation (2.5) is significant only above a certain threshold temperature value called T light off and is a hyperbolic function of the catalyst temperature given as [24]: η(t cat ) = 0.5 tanh(s(t cat T light off )) (2.6) where s describes the slope of the efficiency curve. The terms A out and V cat from (2.4), representing the TWC external surface area and volume respectively, are expressed in terms of the catalyst diameter. The solid phase equation can be written as: ρ s (1 ɛ) c s T cat t = (1 ɛ) λ s 2 T cat z 2 4 D cat h out (T cat T amb ) h A geo (T cat T g ) + Q react (2.7) The initial condition for (2.4) is given as: T cat (z, t = 0) = T cat0 (z) (2.8)
21 2.2 Oxygen Storage Model The objective of a TWC is to reduce the emissions from the exhaust manifold. To achieve this, the concentrations of gases such as CO, NO x, and HC need to be constantly monitored. Estimating the concentration levels of all these gases would require a model that is complex and computationally not feasible. We know that controlling the air-fuel ratio near stoichiometry efficiently converts all the gases to less harmful emissions. Therefore a model that can predict the oxygen storage level in the TWC is used for control purposes. The use of Cerium compensates for any deviation of the AFR from stoichiometry through its capability of storing and releasing oxygen during lean and rich conditions respectively. The chemical reactions describing the oxygen storage process are presented in [13] and given as: O Ce 2 O 3 2 Ce 2 O 4 (2.9a) CO + Ce 2 O 4 CO 2 + Ce 2 O 3 (2.9b) Equation (2.9a) represents the oxygen adsorption on the catalyst surface and Equation (2.9b) represents the CO oxidation with the adsorbed oxygen on the catalyst. The terms Ce 2 O 3 and Ce 2 O 4 represent the reduced and oxidized cerium oxide concentrations respectively and together they contribute to the oxygen storage capacity in the catalyst given as: [Ce 2 O 3 ] + [Ce 2 O 4 ] = OSC (2.10) The net reaction rates due to the forward and backward reactions shown in the above equations are: R 1 = k f 1 [Ce 2O 3 ] 2 [O 2 ] k b 1 [Ce 2 O 4 ] 2 c 0 (2.11a) R 2 = k f 2 [Ce 2O 4 ] [CO] k b 2 [Ce 2 O 3 ] [CO 2 ] (2.11b) where c 0 = P R T g is the total exhaust gas concentration. The symbol [X] is used to represent the chemical concentrations for X = O 2, CO, CO 2, Ce 2 O 3, and Ce 2 O 4. The forward reaction rate terms k f 1 and kf 2 depend on the catalyst temperature and the backward reaction rate terms k b 1 and k b 2 are dependent on the forward reaction term along with the chemical equilibrium
22 constant, which is a function of the Gibbs energy difference term [25]. The concentration balance equations for the chemical species in the gas phase and washcoat can be lumped together as in [17] and written as: [O 2 ] t [CO] t [CO 2 ] t = s v [O 2] z = s v [CO] z = s v [CO 2] z R 1 R 2 + R 2 (2.12) where the space velocity given by s v is a function of the exhaust mass flow rate ṁ exh approximated as: s v = ṁ exh c 0 M exh A cs (2.13) As stated in [16], the storage term is neglected to eliminate computational stiffness involved in its implementation. The final form of the equation is written as: [O 2 ] z [CO] z [CO 2 ] z = R 1 c0 M exh A cs ṁ exh = R 2 c0 M exh A cs ṁ exh = +R 2 c0 M exh A cs ṁ exh (2.14) In order to calculate the concentration levels for each of the chemical species, their stoichiometric values are determined as [26]: [O 2 ] stoich = 0.01 c 0 [CO] stoich = 0.02 c 0 (2.15) [CO 2 ] stoich = 0.12 c 0 For the first cell in the catalyst brick, the concentration levels at stoichiometry is compared with the input signal measurement value λ pre and the maximum of the two is used for computing in further cells. The
23 equations for the first cell calculation are given as: [O 2 ] 1 st cell =max ([O 2] stoich, ((λ pre 0.5) [CO] stoich + (λ pre 1) [CO 2 ] stoich )) ( )) [CO] 1 st cell ([CO] =max [O2 ] stoich + (1 λ pre ) [CO 2 ] stoich stoich, λ pre 0.5 (2.16) [CO 2 ] 1 st cell =[CO 2] stoich = 0.12 c 0 The gas concentrations in every other cell will depend on the values from the previous cell, the reaction rates from Equation 2.11, space velocity s v, and the number of discretized cells. The equations will be further discussed in the next section. The oxygen storage level denoted by Φ represents the oxygen storage sites present in the catalyst. It is also called as the state of oxygen (SOX) and is normalized with respect to the OSC (0 < Φ < 1): Φ t = 1 OSC (2 R 1 R 2 ) (2.17) The initial condition is given as: Φ(z, t = 0) = Φ 0 (z) (2.18) The AFR term λ, measured at the center of the catalyst brick can be written in terms of the individual gas concentrations as: λ mid = 2 ([O 2] + [CO 2 ]) + [CO] 2 ([CO] + [CO 2 ]) (2.19) In order to solve the system PDEs, the parameters in it have been previously identified using a Particle Swarm Optimization algorithm [27]. Their validity has been tested against experimental data simulated over real driving conditions [26]. All the identified parameters are listed under Table Implementation The physics based thermal model of the TWC is first simulated using a numerical technique called Finite Difference Method (FDM). Although using FDM as a standalone control modeling technique does allow for temperature and oxygen storage level monitoring, the small step size for spacial and time domain results in a large system of ODEs which is not desirable for real time implementation. Decreasing the step size reduces the number of ODEs, but results in large error in comparison to measurement values.
24 0 z Figure 2.1: TWC sensor placement (reproduced from [25]). In order to carry out the FDM, Equations 2.2 and 2.7 are discretized with a step size of 40 cells over the spacial coordinate and seconds time step. The input signals to the TWC are exhaust temperature T exh, exhaust gas mass flow rate ṁ exh, and AFR λ pre measured at the input of the catalyst brick. The values of these signals have been measured with the help of sensors as shown in Figure 2.1. The simulation plots for the inputs is shown in Figure 2.2 over the FTP driving cycle. The discretized equations are written as: T j g = T j 1 g + h Ageo Acs z ṁ exh c pg 1 + h Ageo Acs z ṁ exh c pg T j cat (2.20a) dt j cat dt = λ s ρ s c s ( T j+1 cat ) 2 T j cat + T j 1 cat z 2 h A geo ρ s (1 ɛ) c s (T j cat T j g ) + Q react ρ s (1 ɛ) c s 4 h out D cat ρ s (1 ɛ) c s (T j cat T amb ) (2.20b)
25 Parameter Description Value c pg ɛ λ s ρ s c s h h out K reac Specific heat [ of the ] exhaust 1100 gas J (kg K) TWC open cross sectional area [ m 2] 0.80 TWC solid [ phase ] conductivity 1.50 W (m K) Volumetric heat capacity of the catalyst [ ] J m 3 K Convective heat transfer coefficient [ ] W m 2 K Heat transfer coefficient with ambient [ ] W m 2 K Scaling factor for Q reac term s Slope of the TWC efficiency curve function T light off Light-off temperature [K] T amb Ambient temperature [K] A cs TWC cross sectional area [ m 2] A geo TWC specific geometric area [ m 1] 3450 D cat TWC diameter [m] L cat TWC length [m] Table 2.1: Identified parameters of TWC model over FTP driving cycle. Following similar steps for the oxygen storage model equations, we get: dφ j dt = 1 OSC (2 Rj 1 Rj 2 ) (2.21a) R j 1 = kf 1 (OSC Φ)j2 [O 2 ] j k b 1 Φ j2 c 0 R j 2 = kf 2 Φj [CO] j k b 2 (OSC Φ) j [CO 2 ] j (2.21b)
26 Time [s] Figure 2.2: Inputs to the TWC model over FTP cycle (reproduced from [28]). The concentration balance equations are also discretized as follows: max ([O 2 ] stoich, ((λ pre 0.5) [CO] stoich + (λ pre 1) [CO 2 ] stoich )) for j = 1 [O 2 ] j = [O 2 ] j 1 R j 1 c0 M exh A cs z ṁ exh for 2 j Z ( ( )) [O2] max [CO] stoich, stoich +(1 λ pre) [CO 2] stoich [CO] j λ pre 0.5 for j = 1 = [CO] j 1 R j 2 c0 M exh A cs z ṁ exh for 2 j Z 0.12 c 0 for j = 1 [CO 2 ] j = [CO 2 ] j 1 + R j 2 c0 M exh A cs z for 2 j Z ṁ exh (2.22) The discretized equations specified above have been developed based on the TWC schematic shown in Figure 2.3. A 2D plot simulating the thermal dynamics at the center of the catalyst brick and AFR is shown in Figure 2.4. The plot consists of measurement values compared against FDM for 40 discretized cells along
27 Z = 40 cells T exh T cat exh mid pre 1 st j th Z th {[O 2 ], [CO], [CO 2 ]} stoich {[O 2 ], [CO], [CO 2 ]} j-1 OSC + Thermal model {[O 2 ], [CO], [CO 2 ]} j Figure 2.3: Implementation model scheme post TWC discretization [28]. the catalyst length and a time step of seconds. The simulations are carried out over the FTP driving cycle. We plot a 3D simulation in Figure 2.5 for the solid phase temperature variation using FDM. For the oxygen storage level, we plot the simulation along the discretized catalyst cells at three different time instants. In Figure 2.4, we see that Φ varies along the catalyst length as the driving cycle runs over time. Therefore we cannot lump the variations of Φ using a single state throughout the catalyst length. Therefore in this paper we look to develop an ODE solution for Φ. From the 2D plots, the root mean square (RMS) error for T cat and λ mid is found to be 3.92% and 9.89% respectively. The RMS is calculated as follows: ε Tcat = p t=1 (T cat meas T cat sim ) 2 p (2.23) ε λmid = p t=1 (λ mid meas λ mid sim ) 2 p (2.24) where p = number of samples in the driving cycle.
28 Measured Simulated Measured Simulated Time [s] (a) Time = 500 seconds Time = 1000 seconds Time = 1500 seconds 0.5 [-] Number of discretized cells along catalyst length (b) Figure 2.4: (a) Simulation results using FDM over FTP driving cycle for T cat and λ mid. (b) Simulation for Φ at various time instants.
29 (a) Figure 2.5: Simulation results using FDM over FTP driving cycle for T cat
30 Chapter 3 Model Reduction Model reduction is usually carried out to reduce the dimensionality of a complex system usually described by a set of PDEs. This makes the system computationally feasible to simulate and develop models to apply in real world scenarios. Any technique applied to achieve model order reduction must be able to capture the dynamics of the physical system and express it using fewer number of modes in comparison to the original system. Several techniques are available in literature which can be used to carry out the reduction method such as Padé approximation which approximates a function by expressing it as a ratio of two Taylor series expansions with finite degree. It can be used to achieve good approximation of systems which are divergent in nature. Another method called Balanced Truncation is used which is based on eigen value decomposition. In this case a transformation is applied to a dynamical system by first calculating the controllability and observability Gramian. These Gramians are then decomposed into a matrix containing singular values in decreasing order and their order is truncated by choosing only the most dominant values in the matrix as in [29]. For a more comprehensive review of various model reduction methods, the reader can refer to [30]. 3.1 Proper orthogonal decomposition (POD) In this thesis we choose POD and its application with Galerkin projection to achieve a reduced order model of the TWC system. POD is a numerically defined technique which captures the dynamics of a complex system by acquiring the simulated data and generating basis functions from it such that only a few modes of the basis functions can be selected too predict system performance. The input to the POD method 21
31 is a data matrix of the system, in time and spacial domain, simulated over a certain time period and the output is a set of time independent functions that represent the most energetic modes of the original system. The basis functions are generated by performing Singular Value Decomposition (SVD) on the system data and then using an energy criterion, the minimum number of modes required to represent the system are selected. It is crucial to note that the accuracy of the POD method is highly dependent on the quality of the system data generated. The data matrix, or else called a snapshot matrix [31] should consist of the system output values at various time instants across the spatial coordinate which discretized using a fixed step size. If the system output values are not recorded with high accuracy or if the data does not capture the entire system dynamics, then the basis functions calculated using POD will also not reflect the system behavior efficiently. The reduced order model is then obtained by projecting these basis functions onto the original system by using Galerkin's method POD Algorithm In order to use POD technique, we need a snapshot matrix of the PDE based system under consideration. The state variables of a TWC are T cat and Φ. Let us consider a snapshot matrix for T cat of the form, T cat (z 1, t 1 ) T cat (z 1, t 2 )... T cat (z 1, t T ) T cat (z 2, t 1 ) T cat (z 2, t 2 )... T cat (z 2, t T ) S snap = T cat (z Z, t 1 ) T cat (z Z, t 2 )... T cat (z Z, t T ) (3.1) where T cat (z, t) is the catalyst temperature. A snapshot matrix, S snap is generated which consists of the system output values taken at various time instants t (1, T ), where T is the total number of samples in a driving cycle. In this case we are simulating the model on FTP driving cycle with a time step of 0.005, which results in a total number of time samples. The spacial coordinate z (0, Z) is discretized into Z = 40 number of cells for a catalyst of length L cat = 0.068m. Since we are dealing with a complex system which is described by a set of PDEs, there may not exist an analytical solution. Hence we need to solve for an approximate solution that will satisfy the system equation such that the sum squared error (SSE) between
32 the exact and approximate solution is minimum. Let the approximate solution to the system be defined as, T cat (z, t) = N x i (t) ψ i (z) i=0 = ψ N x T (t) (3.2) where x i (t) = time dependent fourier coefficients ψ i (z) = space dependent basis functions (3.3) N = number of modes This approximate solution is substituted in the system of PDEs and solved for the time coefficients x(t) using the Galerkin's method. The solution to the system is then calculated using Equation (3.2). In order to solve for the approximate solution, we also need to have the basis functions calculated prior to solving for x(t). As mentioned earlier this is achieved by performing SVD on the snapshot matrix and then calculating the number of modes required in the basis functions based on a heuristic reasoning. The snapshot matrix S snap R zxt can be decomposed into sub-matrices of the form SV D(S snap ) = UΣV T (3.4) where ] U = [U 1 U 2... U z R zxz ] V = [V 1 V 2... V t R txt σ σ 2. Σ = R. zxt σ z (3.5) U and V are orthonormal matrices, and Σ is a diagonal matrix containing real positive singular values σ of S snap matrix arranged in decreasing order σ 1 σ 2... σ z 0. The basis function are contained in the matrix U, and for these function to accurately reconstruct the approximate solution they need to satisfy the
33 properties of orthogonality and orthonormality defined as: 1 if i=j ψi T ψ j = 0 if i j (3.6) 3.2 Galerkin's Projection The Galerkin's approach looks towards developing an approximation of the exact solution of a PDE. In general finding an exact solution of PDE is complex due to non-linearity and non-homogeneous boundary conditions. Hence finding an approximate solution with the minimum sum squared error turns out to be a very useful method. Galerkin's method follows a basis representation of the systems PDEs by splitting a continuous function (in space and time) into a discrete summation of the two quantities separated. To find the approximate solution, the inner product of the PDE and a weight function is computed over the spacial domain and equated to zero. This ensures that the approximate solution has the least error. Consider an example of a heat equation defined by a parabolic PDE of the form [32]: T t (z, t) = T zz (z, t) + f(t (z, t)) (3.7) where T t = T (z, t), T zz = 2 T (z, t) t z 2 0 < z < L, t > 0 (3.8) In the above equation, T(z,t) is the temperature variable along the z-axis and f(z,t) is a continuous function in the spacial interval. Let the approximate solution of this PDE be given by T. We now define a function called as a residual function R(z,t), which is obtained by substituting the approximate solution in the PDE and equating it to zero. Therefore we can write it as: R( T, T ) = T t (z, t) T zz (z, t) f( T (z, t)) (3.9) We use the expression for T from Equation (3.2) which is described in terms of fourier time coefficients and basis functions. As mentioned earlier, a weight function is also used for approximating the solution which is chosen
34 in a similar way as the basis function. W (z) = ψ j (z) (3.10) The final step in Galerkin's method is to solve for the approximate solution by taking the inner product of the residual and weight function and equating to zero. The final expression is written as: R( T, T ), W (z) = 0 Tt (z, t) T zz (z, t) f( T (z, t)), ψ j (z) = 0 ψn ẋ T (t), ψ j (z) = ψ Nx T (t), ψ j (z) + f(ψ N x T (t)), ψ j (z) (3.11) In order to solve for the time coefficients we multiply Equation (3.11) by ψn T to take advantage of the property from Equation (3.6). Therefore we get the expression, ψ T N ψ N ẋ T (t) = ψ T N ψ N x T (t) + ψ T Nf(ψ N x T (t)) ẋ T (t) = x T (t) + ψ T N f(ψ N x T (t)) (3.12) The above equation describes the model with all the modes of POD being considered. To get the final expression for the reduced order model, we apply the energy heuristic as in 4.2 by replacing N with n: ẋ T (t) = x T (t) + ψ T n f(ψ n x T (t)) (3.13) Equation (3.13) represents the reduced order form of the original PDE computed using POD-Galerkin method.
35 Chapter 4 TWC Reduced Order Model Design In this chapter, we develop a reduced order model of the thermal dynamics of a TWC system using POD-Galerkin projection approach. POD, which is also referred to as Karhunen-Loève decomposition (KLD) or principal component analysis (PCA) [33], is a numerical technique that generate basis functions from system output data [30]. The input to the POD method is a data matrix of the thermal model, in time and spacial domain, computed through simulations over a certain time period and the output is a set of time independent functions that represent the most energetic modes of the original system. The basis functions are generated by performing Singular Value Decomposition (SVD) on the data matrix and then using a criterion to select the minimum number of modes required to reconstruct the system. 4.1 Singular Value Decomposition (SVD) The data matrix generated from the FDM simulation is the snapshot matrix S snap, containing values from the system output in time and space domain. SVD is now performed on this matrix resulting as described in Equation (3.4). The Σ matrix contains the singular values which are arranged in decreasing order of their magnitude and U contains the POD basis vectors. The number of modes in the basis function define the order of the reduced model. Since the singular values in the Σ matrix are in decreasing order, the first few values are more dominant in the sense that they capture most of the system dynamics. Hence we need to define a 26
36 heuristic for selecting the optimal number of basis modes such that the SSE, E SSE = u(z, t) ū(z, t) 2 n (4.1) = u(z, t) x i (t)ψ i (z) 2 i=1 is minimum, where is the L2 norm. One way of selecting the number of modes is by defining an energy criteria as in [34]. We define this criteria in terms of the singular values obtained from SVD as: O n = n i=1 σ2 i N i=1 σ2 i (4.2) where O n = truncation degree N = total number of modes (4.3) n = reduced number of modes This ratio is used to determine the order of the reduced system model. The ratio O n needs to be close to 1 for the basis functions to reconstruct the approximate solution accurately Formulating POD basis vectors After having calculated the SVD matrices, we now apply the energy heuristic as described in the previous section. The plots below show the energy function for different values of mode numbers for both, the temperature and oxygen storage profiles. In order to setup a reduced order model, we set the energy criterion such that the number of modes chosen should capture at least 99% of the system dynamics for the temperature model and 90% for the oxygen storage model. From the plots we can estimate that for T cat, n = 3 number of modes satisfy this criteria for FTP the driving cycle and for Φ, it equals n = Applying Galerkin's projection The reduced form of the TWC system is derived by applying the Galerkin's projection method on the POD basis derived in the previous section. The system variable under consideration in the thermal model
37 1 T cat Figure 4.1: Truncation degree plot for T cat and Φ over FTP driving cycle. is the catalyst temperature, therefore let us define an approximate solution for T cat as: n T cat (z, t) = x i (t) ψ i (z) i=0 = ψ n x T (t) (4.4) The next step is to define the residual function. This can be written using Equation (2.20b): R( T cat, T cat ) = d T cat dt λ s d2 Tcat h A geo ρ s c s dz 2 + ( ρ s (1 ɛ) c T cat T g ) s Q react 4 h out + ( ρ s (1 ɛ) c s D cat ρ s (1 ɛ) c T cat T amb ) s (4.5) To solve for the approximate solution, we take the inner product of the residual and weight function: R( T cat, T cat ), ψ j (z) = 0 (4.6)
38 Substituting the approximate solution of T cat, we can simplify the above expression as: ψn ẋ T (t), ψ j (z) = f( T cat, u cat, t), ψ j (z) (4.7) where, f( T cat, u cat, t) = λ ( s d 2 ) ψ n ρ s c s dz 2 h A geo Q xt (ψ n x T react T g ) + ρ s (1 ɛ) c s ρ s (1 ɛ) c s 4 h out (ψ n x T T amb ) D cat ρ s (1 ɛ) c s ] u cat = [T exh ṁ exh (4.8) contains the non-linear terms of Equation (4.5). Using the property of orthogonality and orthonormality from Equation (3.6), and solving for the inner product, the reduced order non-linear model can be written as: ẋ T = ψ T f( T cat, u cat, t) (4.9) We can express the above equation in terms of the catalyst temperature by multiplying both sides with ψ n as follows: T cat = f( T cat, u cat, t) (4.10) Applying a similar approach to the oxygen storage level equation from 2.21a: m Φ(z, t) = y i (t)δ i (z) i=0 = δ m y T (t) (4.11) Applying Galerkin's projection to the POD basis of the oxygen storage model, we can write the residual function as: R( Φ, Φ) = d Φ dt 1 OSC (2 (kf 1 (OSC Φ) 2 [O 2 ] k b 1 Φ 2 c 0 ) (k f 2 Φ [CO] k b 2 (OSC Φ) [CO 2 ])) (4.12)
39 Taking the inner product of the residual and weight function, we get: R( Φ, Φ), ψ j (z) = 0 (4.13) δm ẏ T (t), δ j (z) = g( Φ, u Φ, t), δ j (z) (4.14) where, g( Φ, u Φ, t) = 1 OSC (2 (kf 1 (OSC δ my T ) 2 [O 2 ] k b 1 (δ m y T ) 2 c 0 ) (k f 2 δ my T [CO] k2 b (OSC δ m y T ) [CO 2 ])) ] u Φ = [λ pre (4.15) Solving for the inner product in the above equation, we arrive at the reduced model as follows: ẏ T = δ T g( Φ, u Φ, t) Φ = g( Φ, u Φ, t) (4.16)
40 Chapter 5 Simulation Results In order to simulate real driving conditions, the control-oriented model is validated over the US06 driving cycle for the temperature and oxygen storage dynamics. In Figure 5.1 the result of the reduced model is compared against the catalyst dynamics simulated from the center location of the TWC over US06 driving cycle. The solid phase temperature T cat simulated by the control oriented model developed using POD-Galerkin using three states is compared to the high dimensional PDE based system. Similarly the eight state oxygen storage model is also plotted against simulated results. The plot also features simulation results for λ mid value which is the AFR of the catalyst measured at the center. We see that the reduced model is able to accurately capture the deviations of the AFR from the stoichiometric value over time. In order to determine the error between the simulated and Galerkin model, we plot the root mean square error (RMSE) for the catalyst temperature and oxygen storage level for different values of POD modes based on the truncation degree plot in Figure 4.1. The RMSE for each variable is calculated as follows: ε Tcat = p t=1 (T cat sim T cat galerkin ) 2 p (5.1) ε Φ = p t=1 (Φ sim Φ galerkin ) 2 p (5.2) ε λmid = p t=1 (λ mid sim λ mid galerkin ) 2 p (5.3) 31
41 where p = number of samples in the driving cycle. As seen below, Figure 5.1 shows the RMSE plots of T cat for n = 3, 5, and 8. As we increase the number of modes, the error reduces significantly since higher number of modes are able to capture more energy. But selecting a large model order results in an increase in computational time and hence we need a balance between accuracy and execution time. Setting an energy cap of 99% as described in section 4 lets us select three modes for describing the temperature dynamics. Similarly the figure also shows the RMSE plots for Φ. In comparison to the T cat RMSE plots, the error in these plots reduces gradually as we increase the number of modes. This is also evident from the truncation degree plot since the curve for the modes of Φ has a smaller slope as compared with T cat modes. Therefore setting an energy cap of 90% for the modes, we can model the oxygen storage dynamics with eight modes. The plots also confirm the reduction in error for λ mid values simulated over US06 driving cycle.
42 Simulated Galerkin Simulated Galerkin Simulated Galerkin Time [s] Figure 5.1: T cat, Φ, and λ mid plots for simulated v/s Galerkin model at the center of the catalyst brick.
43 4 n = n = n = Time [s] (a) 6 4 n = n = n = Time [s] (b)
44 0.05 n = n = n = Time [s] (c) Figure 5.1: RMSE plot of (a) T cat, (b) Φ, and (c) λ mid for different POD modes
Available online at ScienceDirect. IFAC-PapersOnLine (2015)
Available online at www.sciencedirect.com ScienceDirect IFAC-PapersOnLine 48-5 (205) 434 440 A New Semi-Empirical Temperature Model for the Three Way Catalytic Converter Stefano Sabatini, Irfan Kil, Joseph
More informationAUTOMOTIVE EXHAUST AFTERTREATMENT
AUTOMOTIVE EXHAUST AFTERTREATMENT CATALYST FUNDAMENTLS Catalyst in its simplest term is a material that increase the rate (molecules converted by unit time) of a chemical reaction while itself not undergoing
More informationKinetic Parameters Estimation using Vehicle Data for Exhaust Aftertreatment Devices
Kinetic Parameters Estimation using Vehicle Data for Exhaust Aftertreatment Devices Karthik Ramanathan India Science Lab General Motors, Global Research and Development Center Bangalore, India Acknowledgments:
More informationA SCR Model based on Reactor and Engine Experimental Studies for a Cu-zeolite Catalyst
Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports - Open Dissertations, Master's Theses and Master's Reports 213 A SCR Model based on
More informationBalancing chemical reaction equations (stoichiometry)
Balancing chemical reaction equations (stoichiometry) This worksheet and all related files are licensed under the Creative Commons Attribution License, version 1.0. To view a copy of this license, visit
More informationINTERNAL COMBUSTION ENGINE (SKMV 3413)
INTERNAL COMBUSTION ENGINE (SKMV 3413) Dr. Mohd Farid bin Muhamad Said Room : Block P21, Level 1, Automotive Development Centre (ADC) Tel : 07-5535449 Email: mfarid@fkm.utm.my THERMOCHEMISTRY IC engine
More informationCombustion Theory and Applications in CFD
Combustion Theory and Applications in CFD Princeton Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Copyright 201 8 by Heinz Pitsch. This material is not to be sold, reproduced or distributed
More informationALE 1. Chemical Kinetics: Rates of Chemical Reactions
Name Chem 163 Section: Team Number: ALE 1. Chemical Kinetics: Rates of Chemical Reactions (Reference: Sections 16.1 16.2 + parts of 16.5 16.6 Silberberg 5 th edition) How do the surface area, concentration
More informationONE-DIMENSIONAL HIGH-FIDELITY AND REDUCED-ORDER MODELS FOR THREE- WAY CATALYTIC CONVERTER
ONE-DIMENSIONAL HIGH-FIDELITY AND REDUCED-ORDER MODELS FOR THREE- WAY CATALYTIC CONVERTER ONE-DIMENSIONAL HIGH-FIDELITY AND REDUCED-ORDER MODELS FOR THREE- WAY CATALYTIC CONVERTER By Tongrui Li, B.Eng.
More information(g) 2NH 3. (g) ΔH = 92 kj mol 1
1 The uses of catalysts have great economic and environmental importance For example, catalysts are used in ammonia production and in catalytic converters (a) Nitrogen and hydrogen react together in the
More informationINTRODUCTION TO CATALYTIC COMBUSTION
INTRODUCTION TO CATALYTIC COMBUSTION R.E. Hayes Professor of Chemical Engineering Department of Chemical and Materials Engineering University of Alberta, Canada and S.T. Kolaczkowski Professor of Chemical
More informationCopyrighted by Gabriel Tang B.Ed., B.Sc.
Unit 5: Chemical Equations and Reactions & Stoichiometry Chemistry Chapter 9: Stoichiometry 9.1: Calculating Quantities in Reactions Avogadro s Number: - a group of (6.0 10 ) molecules = 1 mole Stoichiometry:
More informationRates of Reaction HL
Name: Rates of Reaction Objectives 16. Rates of Reaction -define rate of reaction -define catalysis -monitor the rate of production of oxygen from hydrogen peroxide, using manganese dioxide as a catalyst
More informationassociation adilca FUEL COMBUSTION
FUEL COMBUSTION Where does the energy that drives the motors come from? What are hydrocarbons made of? Why does it need air to burn fuel? What is the amount of energy released by combustion? How to calculate
More informationThree Way Catalyst Control using PI-style Controller with HEGO Sensor Feedback
Three Way Catalyst Control using PI-style Controller with HEGO Sensor Feedback Per Andersson, Lars Eriksson Dept. of Vehicular Systems, Linköping University, Sweden E-mail: {peran,larer}@isy.liu.se Abstract
More informationFuel, Air, and Combustion Thermodynamics
Chapter 3 Fuel, Air, and Combustion Thermodynamics 3.1) What is the molecular weight, enthalpy (kj/kg), and entropy (kj/kg K) of a gas mixture at P = 1000 kpa and T = 500 K, if the mixture contains the
More informationKINETICS CHEMICAL CHEMIC. Unit. I. Multiple Choice Questions (Type-I)
Unit 4 CHEMICAL CHEMIC KINETICS I. Multiple Choice Questions (Type-I) 1. The role of a catalyst is to change. gibbs energy of reaction. enthalpy of reaction. activation energy of reaction. equilibrium
More informationZeolite HC Traps: Experiments and Simulations
CLEERS 2013 Zeolite HC Traps: Experiments and Simulations Presenter - Ford Motor Company April 10, 2013 p. 1 Chemical Engineering Department Outline -Introduction -Mathematical model development and simulation
More informationChapter Chemical Kinetics
CHM 51 Chapter 13.5-13.7 Chemical Kinetics Graphical Determination of the Rate Law for A Product Plots of [A] versus time, ln[a] versus time, and 1/[A] versus time allow determination of whether a reaction
More informationAAE COMBUSTION AND THERMOCHEMISTRY
5. COMBUSTIO AD THERMOCHEMISTRY Ch5 1 Overview Definition & mathematical determination of chemical equilibrium, Definition/determination of adiabatic flame temperature, Prediction of composition and temperature
More informationCombustion. Indian Institute of Science Bangalore
Combustion Indian Institute of Science Bangalore Combustion Applies to a large variety of natural and artificial processes Source of energy for most of the applications today Involves exothermic chemical
More informationLaminar Premixed Flames: Flame Structure
Laminar Premixed Flames: Flame Structure Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Course Overview Part I: Fundamentals and Laminar Flames Introduction Fundamentals and mass balances of
More informationDARS Digital Analysis of Reactive Systems
DARS Digital Analysis of Reactive Systems Introduction DARS is a complex chemical reaction analysis system, developed by DigAnaRS. Our latest version, DARS V2.0, was released in September 2008 and new
More informationName AP CHEM / / Chapter 12 Outline Chemical Kinetics
Name AP CHEM / / Chapter 12 Outline Chemical Kinetics The area of chemistry that deals with the rate at which reactions occur is called chemical kinetics. One of the goals of chemical kinetics is to understand
More informationModel Based Diagnosis of Leaks in the Air-Intake System of an SI-Engine
Model Based Diagnosis of Leaks in the Air-Intake System of an SI-Engine Mattias Nyberg, Andrej Perkovic Vehicular Systems, ISY, Linköping University S-581 83 Linköping, Sweden e-mail: matny@isy.liu.se
More informationNonlinear dynamics of three-way catalyst with microkinetics and internal diffusion
Nonlinear dynamics of three-way catalyst P. Kočí, V. Nevoral, M. Kubíček, M. Marek Center for Nonlinear Dynamics of Chemical and Biological Systems Prague Institute of Chemical Technology Technická 5,
More informationStandard enthalpies and standard Gibbs energies of formation 4. Calculating enthalpy changes and Gibbs energy changes for reactions 5
Chemical Reactions as Sources of Energy Part 1: Thermodynamics Contents Standard enthalpies and standard Gibbs energies of formation 4 Calculating enthalpy changes and Gibbs energy changes for reactions
More informationDARS overview, IISc Bangalore 18/03/2014
www.cd-adapco.com CH2O Temperatur e Air C2H4 Air DARS overview, IISc Bangalore 18/03/2014 Outline Introduction Modeling reactions in CFD CFD to DARS Introduction to DARS DARS capabilities and applications
More informationOH, is an important feedstock for the chemical industry.
1 Methanol, CH 3 OH, is an important feedstock for the chemical industry. In the manufacture of methanol, carbon dioxide and hydrogen are reacted together in the reversible reaction shown below. CO 2 (g)
More informationHigh-Pressure Volumetric Analyzer
High-Pressure Volumetric Analyzer High-Pressure Volumetric Analysis HPVA II Benefits Dual free-space measurement for accurate isotherm data Free space can be measured or entered Correction for non-ideality
More informationThermodynamic and Stochiometric Principles in Materials Balance
Thermodynamic and Stochiometric Principles in Materials Balance Typical metallurgical engineering problems based on materials and energy balance NiO is reduced in an open atmosphere furnace by excess carbon
More informationWell Stirred Reactor Stabilization of flames
Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the
More informationChapter 9 Rates of Chemical Reactions
Chapter 9 Rates of Chemical Reactions A chemical reaction occurs when collisions between molecules have sufficient energy to break the bonds in the reactants; molecules collide with the proper orientation
More informationRate of reaction refers to the amount of reactant used up or product created, per unit time. We can therefore define the rate of a reaction as:
Rates of Reaction Rate of reaction refers to the amount of reactant used up or product created, per unit time. We can therefore define the rate of a reaction as: Rate = change in concentration units: mol
More informationSTP-TS THERMOPHYSICAL PROPERTIES OF WORKING GASES USED IN WORKING GAS TURBINE APPLICATIONS
THERMOPHYSICAL PROPERTIES OF WORKING GASES USED IN WORKING GAS TURBINE APPLICATIONS THERMOPHYSICAL PROPERTIES OF WORKING GASES USED IN GAS TURBINE APPLICATIONS Prepared by: ASME Standards Technology, LLC
More informationAP* Bonding & Molecular Structure Free Response Questions page 1
AP* Bonding & Molecular Structure Free Response Questions page 1 Essay Questions 1991 a) two points ΔS will be negative. The system becomes more ordered as two gases orm a solid. b) two points ΔH must
More information* 1 K.Mohan Laxmi, 2 V.Ranjith Kumar, 3 Y.V.Hanumantha Rao
Modeling And Simulation of Gas Flow Velocity In Catalytic Converter With Porous * 1 K.Mohan Laxmi, 2 V.Ranjith Kumar, 3 Y.V.Hanumantha Rao 1 M.Tech Student, Mechanical Engineering Department, KL University,
More informationFundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur
Fundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur Lecture 09 Stoichiometric calculations for air-gas mixture Let us start this
More informationFundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur
Fundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur Lecture 10 Mixture fraction calculation for diffusion flames Let us start
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationSUPPLEMENTARY TOPIC 3 ENERGY AND CHEMICAL REACTIONS
SUPPLEMENTARY TOPIC 3 ENERGY AND CHEMICAL REACTIONS Rearranging atoms. In a chemical reaction, bonds between atoms in one or more molecules (reactants) break and new bonds are formed with other atoms to
More informationState how a catalyst speeds up a chemical reaction. ...
Q1. This question is about the use of transition metals as catalysts. (a) State how a catalyst speeds up a chemical reaction. State the characteristic property of transition metals that enables them to
More informationN L N G : C O M B U S T I O N
N L N G : C O M B U S T I O N G R A N T I N G R A M 6 T H D E C E M B E R 2 0 1 7 This short lecture provides an introduction into the basic principles of combustion and some of the background to the use
More informationMathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane
Presented at the COMSOL Conference 2008 Hannover Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane Maryam Ghadrdan Norwegian University of Science and
More informationProblem 1 (Willans Approximation)
5-0567-00 Engine Systems (HS 205) Exercise 3 Topic: Lectures 3+4 Raffi Hedinger (hraffael@ethz.ch), Norbert Zsiga (nzsiga@ethz.ch); October 9, 205 Problem (Willans Approximation) A useful simplification
More informationTHIS IS A NEW SPECIFICATION
THIS IS A NEW SPECIFICATION ADVANCED SUBSIDIARY GCE CHEMISTRY B (SALTERS) Chemistry for Life F331 * OCE / 1 091 6 * Candidates answer on the Question Paper A calculator may be used for this paper OCR Supplied
More informationDevelopment of New Diesel Oxidation and NH3 Slip Catalysts
Downloaded from orbit.dtu.dk on: Feb 11, 218 Development of New Diesel Oxidation and NH3 Slip Catalysts Hansen, Thomas Klint; Jensen, Anker Degn; Hansen, Brian Brun; Janssens, Ton V.W. Publication date:
More informationNUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL
NUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL F. S. Marra*, L. Acampora**, E. Martelli*** marra@irc.cnr.it *Istituto di Ricerche sulla Combustione CNR, Napoli, ITALY *Università
More informationChemical reactions. C2- Topic 5
Chemical reactions C2- Topic 5 What is a chemical reaction? A chemical reaction is a change that takes place when one or more substances (called reactants) form one or more new substances (called products)
More information3.2 Alkanes. Refining crude oil. N Goalby chemrevise.org 40 C 110 C 180 C. 250 C fuel oil 300 C 340 C. Fractional Distillation: Industrially
3.2 Alkanes Refining crude oil Fractional Distillation: Industrially Petroleum is a mixture consisting mainly of alkane hydrocarbons Petroleum fraction: mixture of hydrocarbons with a similar chain length
More informationThe impacts of Pdin BEA zeolite on decreasing cold start HC emission of an E85 vehicle
CLEERS presentation October, 2017 The impacts of Pdin BEA zeolite on decreasing cold start HC emission of an E85 vehicle Lifeng Xu*, Jason Lupescu, Jeffery Hepburn, Giovanni Cavataio, Kevin Guo, Paul Laing,
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationSCR Catalyst Layer Addition Specification and Management
SCR Catalyst Layer Addition Specification and Management The 13th Sulfur Dioxide, Nitrogen Oxide, Mercury and Fine Particle Pollution Control Technology & Management International Exchange Conference Wuhan,
More informationModeling for Control of HCCI Engines
Modeling for Control of HCCI Engines Gregory M. Shaver J.Christian Gerdes Matthew Roelle P.A. Caton C.F. Edwards Stanford University Dept. of Mechanical Engineering D D L ynamic esign aboratory Outline
More informationThe following gas laws describes an ideal gas, where
Alief ISD Chemistry STAAR Review Reporting Category 4: Gases and Thermochemistry C.9.A Describe and calculate the relations between volume, pressure, number of moles, and temperature for an ideal gas as
More informationChapter 6 Thermochemistry 許富銀
Chapter 6 Thermochemistry 許富銀 6.1 Chemical Hand Warmers Thermochemistry: the study of the relationships between chemistry and energy Hand warmers use the oxidation of iron as the exothermic reaction: Nature
More informationUnderstanding NH 3 emissions over a three-way catalyst in lean/rich conditions.
Understanding NH 3 emissions over a three-way catalyst in lean/rich conditions. Tiago Casinhas Ribeiro tiago.casinhas@ist.utl.pt Instituto Superior Técnico, Lisboa, Portugal December 215 Abstract Current
More informationA First Course on Kinetics and Reaction Engineering Example 1.2
Example 1.2 Problem Purpose This example shows how to determine when it is permissible to choose a basis for your calculations. It also illustrates how to use reaction progress variables and the initial
More informationExplanation: They do this by providing an alternative route or mechanism with a lower activation energy
Catalysts Definition: Catalysts increase reaction rates without getting used up. Explanation: They do this by providing an alternative route or mechanism with a lower Comparison of the activation energies
More informationMethod and process for combustion synthesized supported cobalt catalysts for fixed bed Fischer Tropsch reaction
Method and process for combustion synthesized supported cobalt catalysts for fixed bed Fischer Tropsch reaction Center for Sustainable Technologies Indian Institute of Science Bangalore IDF presentation
More information50%-50% Beam Splitters Using Transparent Substrates Coated by Single- or Double-Layer Quarter-Wave Thin Films
University of New Orleans ScholarWorks@UNO University of New Orleans Theses and Dissertations Dissertations and Theses 5-22-2006 50%-50% Beam Splitters Using Transparent Substrates Coated by Single- or
More informationAuthor Index. See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. Affiliation Index.
Author Index Downloaded via 148.251.232.83 on July 5, 2018 at 02:15:59 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. Bart, J. M., 42 Bent,
More informationSmog, Catalytic Converters and You
Smog, Catalytic Converters and You APPENDIX A Outcomes: 1. Describe a reaction mechanism and a catalyst s role in a chemical reaction. (ACC-3) 2. Describe collision theory and its connection to factors
More informationLecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory
Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where
More informationCracking. 191 minutes. 186 marks. Page 1 of 27
3.1.6.2 Cracking 191 minutes 186 marks Page 1 of 27 Q1. (a) Gas oil (diesel), kerosine (paraffin), mineral oil (lubricating oil) and petrol (gasoline) are four of the five fractions obtained by the fractional
More informationEmissions Catalyst Design using GT-SUITE. Dr. Chaitanya Sampara Viridis Chemicals Private Limited
Emissions Catalyst Design using GT-SUITE Dr. Chaitanya Sampara Viridis Chemicals Private Limited Viridis Chemicals Background Originated at KSRM College of Engineering (212) Research Focus Develop tailored
More informationConservation of mass: 44 kg on LHS and RHS one kmol of CO₂ produced by burning 1 kmol of C with one kmol of O₂
L20 Stoichiometry of Reactions Gore Combustion provides the heat source for many power cycles (cars, airplanes, power plants), for home heating, cooking applications, and for manufacturing from melting
More informationCarbon dioxide removal processes by alkanolamines in aqueous organic solvents Hamborg, Espen Steinseth
University of Groningen Carbon dioxide removal processes by alkanolamines in aqueous organic solvents Hamborg, Espen Steinseth IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's
More informationModeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell
Modeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell P. Heidebrecht 1, M. Mangold 2, M. Gundermann 1, A. Kienle 2,3, K. Sundmacher 1,2 1 Otto-von-Guericke-University Magdeburg,
More informationElectrochemistry and battery technology Contents
Electrochemistry and battery technology Contents Introduction Redox overview voltaic cells, electrolytic cells, fuel cells, Primary and secondary batteries. Other batteries; Construction, working and applications
More informationClean Air Zone (CAZ) - CLEAN VEHICLE RETROFIT CERTIFICATION INFORMATION DOCUMENT FOR CVRC LOW EMISSION ADAPTATIONS
Clean Air Zone (CAZ) - CLEAN VEHICLE RETROFIT CERTIFICATION INFORMATION DOCUMENT FOR CVRC LOW EMISSION ADAPTATIONS For use by manufacturers applying to obtain Clean Vehicle Retrofit Certification for their
More informationUnit 12: Chemical Kinetics
Unit 12: Chemical Kinetics Author: S. Michalek Introductory Resources: Zumdahl v. 5 Chapter 12 Main Ideas: Integrated rate laws Half life reactions Reaction Mechanisms Model for chemical kinetics Catalysis
More informationTOPICAL PROBLEMS OF FLUID MECHANICS 97
TOPICAL PROBLEMS OF FLUID MECHANICS 97 DOI: http://dx.doi.org/10.14311/tpfm.2016.014 DESIGN OF COMBUSTION CHAMBER FOR FLAME FRONT VISUALISATION AND FIRST NUMERICAL SIMULATION J. Kouba, J. Novotný, J. Nožička
More informationModel Based Control of Throttle, EGR and Wastegate
Master of Science Thesis in Electrical Engineering Department of Electrical Engineering, Linköping University, 7 Model Based Control of Throttle, EGR and Wastegate A System Analysis of the Gas Flows in
More informationThree Way Catalyst and Lean NO x Trap Modeling for a Lean Burn SIDI Gasoline engine
Three Way Catalyst and Lean NO x Trap Modeling for a Lean Burn SIDI Gasoline engine Jian Gong, Christopher J. Rutland Presented at CLEERS on April 12 th, 213 Engine Research Center University of Wisconsin-Madison
More informationChapter 5. Stoichiometry
Chapter 5 Stoichiometry Chapter 5 Table of Contents (5-1) Counting by weighing (5-2) Atomic masses (5-3) Learning to solve problems (5-4) The mole (5-5) Molar mass (5-6) Percent composition of compounds
More informationTHE GRADIENT PROJECTION ALGORITHM FOR ORTHOGONAL ROTATION. 2 The gradient projection algorithm
THE GRADIENT PROJECTION ALGORITHM FOR ORTHOGONAL ROTATION 1 The problem Let M be the manifold of all k by m column-wise orthonormal matrices and let f be a function defined on arbitrary k by m matrices.
More informationInstrumental Analysis
Instrumental Analysis Bonds and Infrared Energy Covalent bonds can absorb infrared energy, which causes the bonds to vibrate. Depending on the masses of the bonded atoms, different frequencies of radiation
More information(a) Reaction rates (i) Following the course of a reaction Reactions can be followed by measuring changes in concentration, mass and volume of
(a) Reaction rates (i) Following the course of a reaction Reactions can be followed by measuring changes in concentration, mass and volume of reactants or products. g Measuring a change in mass Measuring
More informationRocket Propulsion. Combustion chamber Throat Nozzle
Rocket Propulsion In the section about the rocket equation we explored some of the issues surrounding the performance of a whole rocket. What we didn t explore was the heart of the rocket, the motor. In
More information2/11/2013 INTRODUCTION. STOICHIOMETRY General Chemistry. Uses of H 2 SO 4. Stoichiometry? Big Deal.
INTRODUCTION STOICHIOMETRY General Chemistry Billions of pounds of chemicals are produced each year across the world. These chemicals help manufacture: Medicines Computer chips and electronic instruments
More informationProject-Based Inquiry Science: Air Quality Storyline. Air Quality: What s the Big Question? How Can You Improve Air Quality in Your Community?
Project-Based Inquiry Science: Air Quality Storyline Air Quality: What s the Big Question? How Can You Improve Air Quality in Your Community? In the Introduction to Air Quality, students read a parable
More informationCatalytic bead sensors are used primarily to detect
Chapter 3 Catalytic Combustible Gas Sensors Catalytic bead sensors are used primarily to detect combustible gases. They have been in use for more than 50 years. Initially, these sensors were used for monitoring
More informationSteady-State Molecular Diffusion
Steady-State Molecular Diffusion This part is an application to the general differential equation of mass transfer. The objective is to solve the differential equation of mass transfer under steady state
More informationTopic 6a Electrode Potentials Revision Notes
Topic 6a Electrode Potentials Revision Notes 1) Redox Redox reactions involve the transfer of electrons e.g. in the reaction between zinc metal and copper (II) sulphate, electrons are transferred from
More informationMARIYA INTERNATIONAL SCHOOL. Work sheet II. Term II. Level 8 Chemistry [Paper IV] Name: SULFUR AND AIR AND WATER
MARIYA INTERNATIONAL SCHOOL Work sheet II Term II Level 8 Chemistry [Paper IV] Name: SULFUR AND AIR AND WATER 1. Nitrogen dioxide and other oxides of nitrogen are formed in car engines. a) Explain how
More informationAQA A2 CHEMISTRY TOPIC 4.2 EQUILIBRIA BOOKLET OF PAST EXAMINATION QUESTIONS
AQA A2 CHEMISTRY TOPIC 4.2 EQUILIBRIA BOOKLET OF PAST EXAMINATION QUESTIONS 1 1. (a) The diagram below shows the effect of temperature and pressure on the equilibrium yield of the product in a gaseous
More informationApplication of Environmental Technology Management. (ETM) to Automobile Exhaust Emission Reduction. Meshari AL-Harbi
Application of Environmental Technology Management (ETM) to Automobile Exhaust Emission Reduction by Meshari AL-Harbi A thesis presented to the University of Waterloo in fulfillment of the thesis requirement
More informationIdentification of catalytic converter kinetic model using a genetic algorithm approach
1455 Identification of catalytic converter kinetic model using a genetic algorithm approach G N Pontikakis and A M Stamatelos* Mechanical and Industrial Engineering Department, University of Thessaly,
More informationSystems Engineering Spring Group Project #1: Process Flowsheeting Calculations for Acetic Anhydride Plant. Date: 2/25/00 Due: 3/3/00
10.551 Systems Engineering Spring 2000 Group Project #1: Process Flowsheeting Calculations for Acetic Anhydride Plant Date: 2/25/00 Due: 3/3/00 c Paul I. Barton, 14th February 2000 At our Nowhere City
More informationKinetic Parameter Identification for a DOC Catalyst Using SGB test and Advanced Optimization Algorithms
Kinetic Parameter Identification for a DOC Catalyst Using SGB test and Advanced Optimization Algorithms Mahsa RAFIGH - Politecnico di Torino Federico MILLO - Politecnico di Torino Paolo FERRERI General
More informationCOMBUSTION OF FUEL 12:57:42
COMBUSTION OF FUEL The burning of fuel in presence of air is known as combustion. It is a chemical reaction taking place between fuel and oxygen at temperature above ignition temperature. Heat is released
More informationI.C. Engine Cycles. Thermodynamic Analysis
I.C. Engine Cycles Thermodynamic Analysis AIR STANDARD CYCLES Air as a perfect gas All processes ideal and reversible Mass same throughout Constant Specific Heat. OTTO CYCLE OTTO CYCLE Efficiency is
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More informationCorrelation between LNT NH 3 and N 2 O selectivities under fast cycling conditions
Correlation between LNT and N 2 O selectivities under fast cycling conditions Jae-Soon Choi, Josh Pihl, Miyoung Kim, Bill Partridge, Stuart Daw Oak Ridge National Laboratory Petr Kočí Institute of Chemical
More informationMathematical Model of Catalytic Chemical Reactor. Sintayehu Agegnehu Matintu* School of Mathematical and Statistical Sciences, Hawassa University,
Mathematical Model of Catalytic Chemical Reactor Sintayehu Agegnehu Matintu* School of Mathematical and Statistical Sciences, Hawassa University, P.O.Box 05, Hawassa, Ethiopia Abstract This study concerns
More informationUNIVERSITY OF HOUSTON ENVIRONMENTAL HEALTH AND & LIFE SAFETY POLICIES AND PROCEDURES I. PURPOSE
UNIVERSITY OF HOUSTON ENVIRONMENTAL HEALTH AND SAFETY POLICIES AND PROCEDURES SECTION: SUBJECT: ENVIRONMENTAL HEALTH AND & LIFE SAFETY HAZARD COMMUNICATION I. PURPOSE This procedure establishes the requirements
More informationModification of Alkanes by Cracking
Modification of Alkanes by Cracking Question Paper 3 Level A Level Subject Chemistry Exam Board AQA Module 3.3 Organic Chemistry Topic 3.3.2 Alkanes Sub-Topic 3.3.2.2 Modification of Alkanes by Cracking
More informationAE 205 Materials and Energy Balances Asst. Prof. Dr. Tippabust Eksangsri. Chapter 4 Stoichiometry and MB with Reactions
AE 205 Materials and Energy Balances Asst. Prof. Dr. Tippabust Eksangsri Chapter 4 Stoichiometry and MB with Reactions Stoichiometry Stoichiometry provides a quantitative means of relating the amount of
More information