Filtration and oxidation characteristics of a diesel oxidation catalyst and a catalyzed particulate filter : development of a 1-D 2-layer model

The emissions, filtration and oxidation characteristics of a diesel oxidation catalyst (DOC) and a catalyzed particulate filter (CPF) in a Johnson Matthey catalyzed continuously regenerating trap (CCRT ®) were studied by using computational models. Experimental data needed to calibrate the models were obtained by characterization experiments with raw exhaust sampling from a Cummins ISM 2002 engine with variable geometry turbocharging (VGT) and programmed exhaust gas recirculation (EGR). The experiments were performed at 20, 40, 60 and 75% of full load (1120 Nm) at rated speed (2100 rpm), with and without the DOC upstream of the CPF. This was done to study the effect of temperature and CPF-inlet NO2 concentrations on particulate matter oxidation in the CCRT ®. A previously developed computational model was used to determine the kinetic parameters describing the oxidation characteristics of HCs, CO and NO in the DOC and the pressure drop across it. The model was calibrated at five temperatures in the range of 280 – 465° C, and exhaust volumetric flow rates of 0.447 – 0.843 act-m3/sec. The downstream HCs, CO and NO concentrations were predicted by the DOC model to within ±3 ppm. The HCs and CO oxidation kinetics in the temperature range of 280 - 465°C and an exhaust volumetric flow rate of 0.447 - 0.843 act-m3/sec can be represented by one ’apparent’ activation energy and pre-exponential factor. The NO oxidation kinetics in the same temperature and exhaust flow rate range can be represented by ’apparent’ activation energies and pre-exponential factors in two regimes. The DOC pressure drop was always predicted within 0.5 kPa by the model. The MTU 1-D 2-layer CPF model was enhanced in several ways to better model the performance of the CCRT ®. A model to simulate the oxidation of particulate inside the filter wall was developed. A particulate cake layer filtration model which describes particle filtration in terms of more fundamental parameters was developed and coupled to the wall oxidation model. To better model the particulate oxidation kinetics, a model to take into account the NO2 produced in the washcoat of the CPF was developed. The overall 1-D 2-layer model can be used to predict the pressure drop of the exhaust gas across the filter, the evolution of particulate mass inside the filter, the particulate mass oxidized, the filtration efficiency and the particle number distribution downstream of the CPF. The model was used to better understand the internal performance of the CCRT®, by determining the components of the total pressure drop across the filter, by classifying the total particulate matter in layer I, layer II, the filter wall, and by the means of oxidation i.e. by O2, NO2 entering the filter and by NO2 being produced in the filter. The CPF model was calibrated at four temperatures in the range of 280 – 465 °C, and exhaust volumetric flow rates of 0.447 – 0.843 act-m3/sec, in CPF-only and CCRT ® (DOC+CPF) configurations. The clean filter wall permeability was determined to be 2.00E-13 m2, which is in agreement with values in the literature for cordierite filters. The particulate packing density in the filter wall had values between 2.92 kg/m3 - 3.95 kg/m3 for all the loads. The mean pore size of the catalyst loaded filter wall was found to be 11.0 µm. The particulate cake packing densities and permeabilities, ranged from 131 kg/m3 - 134 kg/m3, and 0.42E-14 m2 and 2.00E-14 m2 respectively, and are in agreement with the Peclet number correlations in the literature. Particulate cake layer porosities determined from the particulate cake layer filtration model ranged between 0.841 and 0.814 and decreased with load, which is about 0.1 lower than experimental and more complex discrete particle simulations in the literature. The thickness of layer I was kept constant at 20 µm. The model kinetics in the CPF-only and CCRT ® configurations, showed that no ’catalyst effect’ with O2 was present. The kinetic parameters for the NO2-assisted oxidation of particulate in the CPF were determined from the simulation of transient temperature programmed oxidation data in the literature. It was determined that the thermal and NO2 kinetic parameters do not change with temperature, exhaust flow rate or NO2 concentrations. However, different kinetic parameters are used for particulate oxidation in the wall and on the wall. Model results showed that oxidation of particulate in the pores of the filter wall can cause disproportionate decreases in the filter pressure drop with respect to particulate mass. The wall oxidation model along with the particulate cake filtration model were developed to model the sudden and rapid decreases in pressure drop across the CPF. The particulate cake and wall filtration models result in higher particulate filtration efficiencies than with just the wall filtration model, with overall filtration efficiencies of 98-99% being predicted by the model. The pre-exponential factors for oxidation by NO2 did not change with temperature or NO2 concentrations because of the NO2 wall production model. In both CPF-only and CCRT ® configurations, the model showed NO2 and layer I to be the dominant means and dominant physical location of particulate oxidation respectively. However, at temperatures of 280 °C, NO2 is not a significant oxidizer of particulate matter, which is in agreement with studies in the literature. The model showed that 8.6 and 81.6% of the CPF-inlet particulate matter was oxidized after 5 hours at 20 and 75% load in CCRT® configuration. In CPF-only configuration at the same loads, the model showed that after 5 hours, 4.4 and 64.8% of the inlet particulate matter was oxidized. The increase in NO2 concentrations across the DOC contributes significantly to the oxidation of particulate in the CPF and is supplemented by the oxidation of NO to NO2 by the catalyst in the CPF, which increases the particulate oxidation rates. From the model, it was determined that the catalyst in the CPF modeslty increases the particulate oxidation rates in the range of 4.5 – 8.3% in the CCRT® configuration. Hence, the catalyst loading in the CPF of the CCRT® could possibly be reduced without significantly decreasing particulate oxidation rates leading to catalyst cost savings and better engine performance due to lower exhaust backpressures.

Implementation of a phenomenological evaporation model into a porous network simulation for water management in low temperature fuel cells

A phenomenological transition film evaporation model was introduced to a pore network model with the consideration of pore radius, contact angle, non-isothermal interface temperature, microscale fluid flows and heat and mass transfers. This was achieved by modeling the transition film region of the menisci in each pore throughout the porous transport layer of a half-cell polymer electrolyte membrane (PEM) fuel cell. The model presented in this research is compared with the standard diffusive fuel cell modeling approach to evaporation and shown to surpass the conventional modeling approach in terms of predicting the evaporation rates in porous media. The current diffusive evaporation models used in many fuel cell transport models assumes a constant evaporation rate across the entire liquid-air interface. The transition film model was implemented into the pore network model to address this issue and create a pore size dependency on the evaporation rates. This is accomplished by evaluating the transition film evaporation rates determined by the kinetic model for every pore containing liquid water in the porous transport layer (PTL). The comparison of a transition film and diffusive evaporation model shows an increase in predicted evaporation rates for smaller pore sizes with the transition film model. This is an important parameter when considering the micro-scaled pore sizes seen in the PTL and becomes even more substantial when considering transport in fuel cells containing an MPL, or a large variance in pore size. Experimentation was performed to validate the transition film model by monitoring evaporation rates from a non-zero contact angle water droplet on a heated substrate. The substrate was a glass plate with a hydrophobic coating to reduce wettability. The tests were performed at a constant substrate temperature and relative humidity. The transition film model was able to accurately predict the drop volume as time elapsed. By implementing the transition film model to a pore network model the evaporation rates present in the PTL can be more accurately modeled. This improves the ability of a pore network model to predict the distribution of liquid water and ultimately the level of flooding exhibited in a PTL for various operating conditions.

Characterization of a Coal Fly Ash-Cement Slurry by the Absolute Foam Index

This dissertation established a standard foam index: the absolute foam index test. This test characterized a wide range of coal fly ash by the absolute volume of air-entraining admixture (AEA) necessary to produce a 15-second metastable foam in a coal fly ash-cement slurry in a specified time. The absolute foam index test was used to characterize fly ash samples having loss on ignition (LOI) values that ranged from 0.17 to 23.3 %wt. The absolute foam index characterized the fly ash samples by absolute volume of AEA, defined as the amount of undiluted AEA solution added to obtain a 15-minute endpoint signified by 15-second metastable foam. Results were compared from several foam index test time trials that used different initial test concentrations to reach termination at selected times. Based on the coefficient of variation (CV), a 15-minute endpoint, with limits of 12 to 18 minutes was chosen. Various initial test concentrations were used to accomplish consistent contact times and concentration gradients for the 15-minute test endpoint for the fly ash samples. A set of four standard concentrations for the absolute foam index test were defined by regression analyses and a procedure simplifying the test process. The set of standard concentrations for the absolute foam index test was determined by analyzing experimental results of 80 tests on coal fly ashes with loss on ignition (LOI) values ranging from 0.39 to 23.3 wt.%. A regression analysis informed selection of four concentrations (2, 6, 10, and 15 vol.% AEA) that are expected to accommodate fly ashes with 0.39 to 23.3 wt.% LOI, depending on the AEA type. Higher concentrations should be used for high-LOI fly ash when necessary. A procedure developed using these standard concentrations is expected to require only 1-3 trials to meet specified endpoint criteria for most fly ashes. The AEA solution concentration that achieved the metastable foam in the foam index test was compared to the AEA equilibrium concentration obtained from the direct adsorption isotherm test with the same fly ash. The results showed that the AEA concentration that satisfied the absolute foam index test was much less than the equilibrium concentration. This indicated that the absolute foam index test was not at or near equilibrium. Rather, it was a dynamic test where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms. Equilibrium isotherm equations obtained from direct isotherm tests were used to calculate the equilibrium concentrations and capacities of fly ash from 0.17 to 10.5% LOI. The results showed that the calculated fly ash capacity was much less than capacities obtained from isotherm tests that were conducted with higher initial concentrations. This indicated that the absolute foam index was not equilibrium. Rather, the test is dynamic where the time of the test played an important role in the results. Even though the absolute foam index was not an equilibrium condition, a correlation was made between the absolute foam index and adsorption isotherms for fly ash of 0.17 to 10.5% LOI. Several batches of mortars were mixed for the same fly ash type increasing only the AEA concentration (dosage) in each subsequent batch. Mortar air test results for each batch showed for each increase in AEA concentration, air contents increased until a point where the next increase in AEA concentration resulted in no increase in air content. This was maximum air content that could be achieved by the particular mortar system; the system reached its air capacity at the saturation limit. This concentration of AEA was compared to the critical micelle concentration (CMC) for the AEA and the absolute foam index.

Dynamic modeling of active regeneration in catalyzed and non-catalyzed diesel particulate filters

Particulate matter (PM) emissions standards set by the US Environmental Protection Agency (EPA) have become increasingly stringent over the years. The EPA regulation for PM in heavy duty diesel engines has been reduced to 0.01 g/bhp-hr for the year 2010. Heavy duty diesel engines make use of an aftertreatment filtration device, the Diesel Particulate Filter (DPF). DPFs are highly efficient in filtering PM (known as soot) and are an integral part of 2010 heavy duty diesel aftertreatment system. PM is accumulated in the DPF as the exhaust gas flows through it. This PM needs to be removed by oxidation periodically for the efficient functioning of the filter. This oxidation process is also known as regeneration. There are 2 types of regeneration processes, namely active regeneration (oxidation of PM by external means) and passive oxidation (oxidation of PM by internal means). Active regeneration occurs typically in high temperature regions, about 500 - 600 °C, which is much higher than normal diesel exhaust temperatures. Thus, the exhaust temperature has to be raised with the help of external devices like a Diesel Oxidation Catalyst (DOC) or a fuel burner. The O2 oxidizes PM producing CO2 as oxidation product. In passive oxidation, one way of regeneration is by the use of NO2. NO2 oxidizes the PM producing NO and CO2 as oxidation products. The passive oxidation process occurs at lower temperatures (200 - 400 °C) in comparison to the active regeneration temperatures. Generally, DPF substrate walls are washcoated with catalyst material to speed up the rate of PM oxidation. The catalyst washcoat is observed to increase the rate of PM oxidation. The goal of this research is to develop a simple mathematical model to simulate the PM depletion during the active regeneration process in a DPF (catalyzed and non-catalyzed). A simple, zero-dimensional kinetic model was developed in MATLAB. Experimental data required for calibration was obtained by active regeneration experiments performed on PM loaded mini DPFs in an automated flow reactor. The DPFs were loaded with PM from the exhaust of a commercial heavy duty diesel engine. The model was calibrated to the data obtained from active regeneration experiments. Numerical gradient based optimization techniques were used to estimate the kinetic parameters of the model.

Numerical investigation of effects of addition of ethanol to gasoline on Laminar Flame Speed (LFS), autoignition, and wall quenching

Since the advent of automobiles, alcohol has been considered a possible engine fuel1,2. With the recent increased concern about the high price of crude oil due to fluctuating supply and demand and environmental issues, interest in alcohol based fuels has increased2,3. However, using pure alcohols or blends with conventional fuels in high percentages requires changes to the engine and fuel system design2. This leads to the need for a simple and accurate conventional fuels-alcohol blends combustion models that can be used in developing parametric burn rate and knock combustion models for designing more efficient Spark Ignited (SI) engines. To contribute to this understanding, numerical simulations were performed to obtain detailed characteristics of Gasoline-Ethanol blends with respect to Laminar Flame Speed (LFS), autoignition and Flame-Wall interactions. The one-dimensional premixed flame code CHEMKIN® was applied to simulate the burning velocity and autoignition characteristics using the freely propagating model and closed homogeneous reactor model respectively. Computational Fluid Dynamics (CFD) was used to obtain detailed flow, temperature, and species fields for Flame-wall interactions. A semi-detailed validated chemical kinetic model for a gasoline surrogate fuel developed by Andrae and Head4 was used for the study of LFS and Autoignition. For the quenching study, a skeletal chemical kinetic mechanism of gasoline surrogate, having 50 species and 174 reactions was used. The surrogate fuel was defined as a mixture of pure n-heptane, isooctane, and toluene. For LFS study, the ethanol volume fraction was varied from 0 to 85%, initial pressure from 4 to 8 bar, initial temperature from 300 to 900K, and dilution from 0 to 32%. Whereas for Autoignition study, the ethanol volume fraction was varied between 0 to 85%, initial pressure was varied between 20 to 60 bar, initial temperature was varied between 800 to 1200K, and the dilution was varied between 0 to 32% at equivalence ratios of 0.5, 1.0 and 1.5 to represent the in-cylinder conditions of a SI engine. For quenching study three Ethanol blends, namely E0, E25 and E85 are described in detail at an initial pressure of 8 atm and 17 atm. Initial wall temperature was taken to be 400 K. Quenching thicknesses and heat fluxes to the wall were computed. The laminar flame speed was found to increase with ethanol concentration and temperature but decrease with pressure and dilution. The autoignition time was found to increase with ethanol concentration at lower temperatures but was found to decrease marginally at higher temperatures. The autoignition time was also found to decrease with pressure and equivalence ratio but increase with dilution. The average quenching thickness was found to decrease with an increase in Ethanol concentration in the blend. Heat flux to the wall increased with increase in ethanol percentage in the blend and at higher initial pressures. Whereas the wall heat flux decreased with an increase in dilution. Unburned Hydrocarbon (UHC) and CO % was also found to decrease with ethanol concentration in the blend.

Sustainable Energy Production in the United States: Life Cycle Assessment of Biofuels and Bioenergy

The United States of America is making great efforts to transform the renewable and abundant biomass resources into cost-competitive, high-performance biofuels, bioproducts, and biopower. This is the key to increase domestic production of transportation fuels and renewable energy, and reduce greenhouse gas and other pollutant emissions. This dissertation focuses specifically on assessing the life cycle environmental impacts of biofuels and bioenergy produced from renewable feedstocks, such as lignocellulosic biomass, renewable oils and fats. The first part of the dissertation presents the life cycle greenhouse gas (GHG) emissions and energy demands of renewable diesel (RD) and hydroprocessed jet fuels (HRJ). The feedstocks include soybean, camelina, field pennycress, jatropha, algae, tallow and etc. Results show that RD and HRJ produced from these feedstocks reduce GHG emissions by over 50% compared to comparably performing petroleum fuels. Fossil energy requirements are also significantly reduced. The second part of this dissertation discusses the life cycle GHG emissions, energy demands and other environmental aspects of pyrolysis oil as well as pyrolysis oil derived biofuels and bioenergy. The feedstocks include waste materials such as sawmill residues, logging residues, sugarcane bagasse and corn stover, and short rotation forestry feedstocks such as hybrid poplar and willow. These LCA results show that as much as 98% GHG emission savings is possible relative to a petroleum heavy fuel oil. Life cycle GHG savings of 77 to 99% were estimated for power generation from pyrolysis oil combustion relative to fossil fuels combustion for electricity, depending on the biomass feedstock and combustion technologies used. Transportation fuels hydroprocessed from pyrolysis oil show over 60% of GHG reductions compared to petroleum gasoline and diesel. The energy required to produce pyrolysis oil and pyrolysis oil derived biofuels and bioelectricity are mainly from renewable biomass, as opposed to fossil energy. Other environmental benefits include human health, ecosystem quality and fossil resources. The third part of the dissertation addresses the direct land use change (dLUC) impact of forest based biofuels and bioenergy. An intensive harvest of aspen in Michigan is investigated to understand the GHG mitigation with biofuels and bioenergy production. The study shows that the intensive harvest of aspen in MI compared to business as usual (BAU) harvesting can produce 18.5 billion gallons of ethanol to blend with gasoline for the transport sector over the next 250 years, or 32.2 billion gallons of bio-oil by the fast pyrolysis process, which can be combusted to generate electricity or upgraded to gasoline and diesel. Intensive harvesting of these forests can result in carbon loss initially in the aspen forest, but eventually accumulates more carbon in the ecosystem, which translates to a CO2 credit from the dLUC impact. Time required for the forest-based biofuels to reach carbon neutrality is approximately 60 years. The last part of the dissertation describes the use of depolymerization model as a tool to understand the kinetic behavior of hemicellulose hydrolysis under dilute acid conditions. Experiments are carried out to measure the concentrations of xylose and xylooligomers during dilute acid hydrolysis of aspen. The experiment data are used to fine tune the parameters of the depolymerization model. The results show that the depolymerization model successfully predicts the xylose monomer profile in the reaction, however, it overestimates the concentrations of xylooligomers.

A SCR Model based on Reactor and Engine Experimental Studies for a Cu-zeolite Catalyst

A NOx reduction efficiency higher than 95% with NH3 slip less than 30 ppm is desirable for heavy-duty diesel (HDD) engines using selective catalytic reduction (SCR) systems to meet the US EPA 2010 NOx standard and the 2014-2018 fuel consumption regulation. The SCR performance needs to be improved through experimental and modeling studies. In this research, a high fidelity global kinetic 1-dimensional 2-site SCR model with mass transfer, heat transfer and global reaction mechanisms was developed for a Cu-zeolite catalyst. The model simulates the SCR performance for the engine exhaust conditions with NH3 maldistribution and aging effects, and the details are presented. SCR experimental data were collected for the model development, calibration and validation from a reactor at Oak Ridge National Laboratory (ORNL) and an engine experimental setup at Michigan Technological University (MTU) with a Cummins 2010 ISB engine. The model was calibrated separately to the reactor and engine data. The experimental setup, test procedures including a surrogate HD-FTP cycle developed for transient studies and the model calibration process are described. Differences in the model parameters were determined between the calibrations developed from the reactor and the engine data. It was determined that the SCR inlet NH3 maldistribution is one of the reasons causing the differences. The model calibrated to the engine data served as a basis for developing a reduced order SCR estimator model. The effect of the SCR inlet NO2/NOx ratio on the SCR performance was studied through simulations using the surrogate HD-FTP cycle. The cumulative outlet NOx and the overall NOx conversion efficiency of the cycle are highest with a NO2/NOx ratio of 0.5. The outlet NH3 is lowest for the NO2/NOx ratio greater than 0.6. A combined engine experimental and simulation study was performed to quantify the NH3 maldistribution at the SCR inlet and its effects on the SCR performance and kinetics. The uniformity index (UI) of the SCR inlet NH3 and NH3/NOx ratio (ANR) was determined to be below 0.8 for the production system. The UI was improved to 0.9 after installation of a swirl mixer into the SCR inlet cone. A multi-channel model was developed to simulate the maldistribution effects. The results showed that reducing the UI of the inlet ANR from 1.0 to 0.7 caused a 5-10% decrease in NOx reduction efficiency and 10-20 ppm increase in the NH3 slip. The simulations of the steady-state engine data with the multi-channel model showed that the NH3 maldistribution is a factor causing the differences in the calibrations developed from the engine and the reactor data. The Reactor experiments were performed at ORNL using a Spaci-IR technique to study the thermal aging effects. The test results showed that the thermal aging (at 800°C for 16 hours) caused a 30% reduction in the NH3 stored on the catalyst under NH3 saturation conditions and different axial concentration profiles under SCR reaction conditions. The kinetics analysis showed that the thermal aging caused a reduction in total NH3 storage capacity (94.6 compared to 138 gmol/m3), different NH3 adsorption/desorption properties and a decrease in activation energy and the pre-exponential factor for NH3 oxidation, standard and fast SCR reactions. Both reduction in the storage capability and the change in kinetics of the major reactions contributed to the change in the axial storage and concentration profiles observed from the experiments.

Parameter identification of a copper-zeolite SCR catalyst model using reactor data

The selective catalytic reduction system is a well established technology for NOx emissions control in diesel engines. A one dimensional, single channel selective catalytic reduction (SCR) model was previously developed using Oak Ridge National Laboratory (ORNL) generated reactor data for an iron-zeolite catalyst system. Calibration of this model to fit the experimental reactor data collected at ORNL for a copper-zeolite SCR catalyst is presented. Initially a test protocol was developed in order to investigate the different phenomena responsible for the SCR system response. A SCR model with two distinct types of storage sites was used. The calibration process was started with storage capacity calculations for the catalyst sample. Then the chemical kinetics occurring at each segment of the protocol was investigated. The reactions included in this model were adsorption, desorption, standard SCR, fast SCR, slow SCR, NH3 Oxidation, NO oxidation and N2O formation. The reaction rates were identified for each temperature using a time domain optimization approach. Assuming an Arrhenius form of the reaction rates, activation energies and pre-exponential parameters were fit to the reaction rates. The results indicate that the Arrhenius form is appropriate and the reaction scheme used allows the model to fit to the experimental data and also for use in real world engine studies.

Electrospun Quaternized Chitosan Fibers for Virus Removal from Drinking Water

Membrane filtration has become an accepted technology for the removal of pathogens from drinking water. Viruses, known to contaminate water supplies, are too small to be removed by a size-exclusion mechanism without a large energy penalty. Thus, functionalized electrospun membranes that can adsorb viruses have drawn our interest. We chose a quaternized chitosan derivative (HTCC) which carries a positively-charged quaternary amine, known to bind negatively-charged virus particles, as a functionalized membrane material. The technique of electrospinning was utilized to produce nanofiber mats with large pore diameters to increase water flux and decrease membrane fouling. In this study, stable, functionalized, electrospun HTCC-PVA nanofibers that can remove 3.6 logs (99.97%) of a model virus, porcine parvovirus (PPV), from water by adsorption and filtration have been successfully produced. This technology has the potential to purify drinking water in undeveloped countries and reduce the number of deaths due to lack of sanitation.