Modeling of diesel particulate filter filtration and regeneration for transient driving schedules

Due to their high thermal efficiency, diesel engines have excellent fuel economy and have been widely used as a power source for many vehicles. Diesel engines emit less greenhouse gases (carbon dioxide) compared with gasoline engines. However, diesel engines emit large amounts of particulate matter (PM) which can imperil human health. The best way to reduce the particulate matter is by using the Diesel Particulate Filter (DPF) system which consists of a wall-flow monolith which can trap particulates, and the DPF can be periodically regenerated to remove the collected particulates. The estimation of the PM mass accumulated in the DPF and total pressure drop across the filter are very important in order to determine when to carry out the active regeneration for the DPF. In this project, by developing a filtration model and a pressure drop model, we can estimate the PM mass and the total pressure drop, then, these two models can be linked with a regeneration model which has been developed previously to predict when to regenerate the filter. There results of this project were: 1 Reproduce a filtration model and simulate the processes of filtration. By studying the deep bed filtration and cake filtration, stages and quantity of mass accumulated in the DPF can be estimated. It was found that the filtration efficiency increases faster during the deep-bed filtration than that during the cake filtration. A “unit collector” theory was used in our filtration model which can explain the mechanism of the filtration very well. 2 Perform a parametric study on the pressure drop model for changes in engine exhaust flow rate, deposit layer thickness, and inlet temperature. It was found that there are five primary variables impacting the pressure drop in the DPF which are temperature gradient along the channel, deposit layer thickness, deposit layer permeability, wall thickness, and wall permeability. 3 Link the filtration model and the pressure drop model with the regeneration model to determine the time to carry out the regeneration of the DPF. It was found that the regeneration should be initiated when the cake layer is at a certain thickness, since a cake layer with either too big or too small an amount of particulates will need more thermal energy to reach a higher regeneration efficiency. 4 Formulate diesel particulate trap regeneration strategies for real world driving conditions to find out the best desirable conditions for DPF regeneration. It was found that the regeneration should be initiated when the vehicle’s speed is high and during which there should not be any stops from the vehicle. Moreover, the regeneration duration is about 120 seconds and the inlet temperature for the regeneration is 710K.

A comparison of ground cover and frequency estimation methods for post-harvest soil monitoring

The amount and type of ground cover is an important characteristic to measure when collecting soil disturbance monitoring data after a timber harvest. Estimates of ground cover and bare soil can be used for tracking changes in invasive species, plant growth and regeneration, woody debris loadings, and the risk of surface water runoff and soil erosion. A new method of assessing ground cover and soil disturbance was recently published by the U.S. Forest Service, the Forest Soil Disturbance Monitoring Protocol (FSDMP). This protocol uses the frequency of cover types in small circular (15cm) plots to compare ground surface in pre- and post-harvest condition. While both frequency and percent cover are common methods of describing vegetation, frequency has rarely been used to measure ground surface cover. In this study, three methods for assessing ground cover percent (step-point, 15cm dia. circular and 1x5m visual plot estimates) were compared to the FSDMP frequency method. Results show that the FSDMP method provides significantly higher estimates of ground surface condition for most soil cover types, except coarse wood. The three cover methods had similar estimates for most cover values. The FSDMP method also produced the highest value when bare soil estimates were used to model erosion risk. In a person-hour analysis, estimating ground cover percent in 15cm dia. plots required the least sampling time, and provided standard errors similar to the other cover estimates even at low sampling intensities (n=18). If ground cover estimates are desired in soil monitoring, then a small plot size (15cm dia. circle), or a step-point method can provide a more accurate estimate in less time than the current FSDMP method.

Muffler characterization with implementation of the finite element method and experimental techniques

Determining how an exhaust system will perform acoustically before a prototype muffler is built can save the designer both a substantial amount of time and resources. In order to effectively use the simulation tools available it is important to understand what is the most effective tool for the intended purpose of analysis as well as how typical elements in an exhaust system affect muffler performance. An in-depth look at the available tools and their most beneficial uses are presented in this thesis. A full parametric study was conducted using the FEM method for typical muffler elements which was also correlated to experimental results. This thesis lays out the overall ground work on how to accurately predict sound pressure levels in the free field for an exhaust system with the engine properties included. The accuracy of the model is heavily dependent on the correct temperature profile of the model in addition to the accuracy of the source properties. These factors will be discussed in detail and methods for determining them will be presented. The secondary effects of mean flow, which affects both the acoustical wave propagation and the flow noise generation, will be discussed. Effective ways for predicting these secondary effects will be described. Experimental models will be tested on a flow rig that showcases these phenomena.

STUDY OF SPARK IGNITION ENGINE COMBUSTION MODEL FOR THE ANALYSIS OF CYCLIC VARIATION AND COMBUSTION STABILITY AT LEAN OPERATING CONDITIONS

A fundamental combustion model for spark-ignition engine is studied in this report. The model is implemented in SIMULINK to simulate engine outputs (mass fraction burn and in-cylinder pressure) under various engine operation conditions. The combustion model includes a turbulent propagation and eddy burning processes based on literature [1]. The turbulence propagation and eddy burning processes are simulated by zero-dimensional method and the flame is assumed as sphere. To predict pressure, temperature and other in-cylinder variables, a two-zone thermodynamic model is used. The predicted results of this model match well with the engine test data under various engine speeds, loads, spark ignition timings and air fuel mass ratios. The developed model is used to study cyclic variation and combustion stability at lean (or diluted) combustion conditions. Several variation sources are introduced into the combustion model to simulate engine performance observed in experimental data. The relations between combustion stability and the introduced variation amount are analyzed at various lean combustion levels.

A Novel Approach to Carbon Dioxide Capture and Storage

The novel approach to carbon capture and storage (CCS) described in this dissertation is a significant departure from the conventional approach to CCS. The novel approach uses a sodium carbonate solution to first capture CO2 from post combustion flue gas streams. The captured CO2 is then reacted with an alkaline industrial waste material, at ambient conditions, to regenerate the carbonate solution and permanently store the CO2 in the form of an added value carbonate mineral. Conventional CCS makes use of a hazardous amine solution for CO2 capture, a costly thermal regeneration stage, and the underground storage of supercritical CO2. The objective of the present dissertation was to examine each individual stage (capture and storage) of the proposed approach to CCS. Study of the capture stage found that a 2% w/w sodium carbonate solution was optimal for CO2 absorption in the present system. The 2% solution yielded the best tradeoff between the CO2 absorption rate and the CO2 absorption capacity of the solutions tested. Examination of CO2 absorption in the presence of flue gas impurities (NOx and SOx) found that carbonate solutions possess a significant advantage over amine solutions, that they could be used for multi-pollutant capture. All the NOx and SOx fed to the carbonate solution was able to be captured. Optimization studies found that it was possible to increase the absorption rate of CO2 into the carbonate solution by adding a surfactant to the solution to chemically alter the gas bubble size. The absorption rate of CO2 was increased by as much as 14%. Three coal combustion fly ash materials were chosen as the alkaline industrial waste materials to study the storage CO2 and regeneration the absorbent. X-ray diffraction analysis on reacted fly ash samples confirmed that the captured CO2 reacts with the fly ash materials to form a carbonate mineral, specifically calcite. Studies found that after a five day reaction time, 75% utilization of the waste material for CO2 storage could be achieved, while regenerating the absorbent. The regenerated absorbent exhibited a nearly identical CO2 absorption capacity and CO2 absorption rate as a fresh Na2CO3 solution.

AN EXPERIMENTAL INVESTIGATION INTO THE EFFECTS OF BIODIESEL BLENDS ON PARTICULATE MATTER OXIDATION IN A CATALYZED PARTICULATE FILTER DURING ACTIVE REGENERATION

Active regeneration experiments were carried out on a production 2007 Cummins 8.9L ISL engine and associated DOC and CPF aftertreatment system. The effects of SME biodiesel blends were investigated in this study in order to determine the PM oxidation kinetics associated with active regeneration, and to determine the effect of biodiesel on them. The experimental data from this study will also be used to calibrate the MTU-1D CPF model. Accurately predicting the PM mass retained in the CPF and the oxidation characteristics will provide the basis for computation in the ECU that will minimize the fuel penalty associated with active regeneration. An active regeneration test procedure was developed based on previous experimentation at MTU. During each experiment, the PM mass in the CPF is determined by weighing the filter at various phases. In addition, DOC and CPF pressure drop, particle size distribution, gaseous emissions, temperature, and PM concentration data are collected and recorded throughout each experiment. The experiments covered a range of CPF inlet temperatures using ULSD, B10, and B20 blends of biodiesel. The majority of the tests were performed at CPF PM loading of 2.2 g/L with in-cylinder dosing, although 4.1 g/L and a post-turbo dosing injector were also used. The PM oxidation characteristics at different test conditions were studied in order to determine the effects of biodiesel on PM oxidation during active regeneration. A PM reaction rate calculation method was developed to determine the global activation energy and the corresponding pre-exponential factor for all test fuels. The changing sum of the total flow resistance of the wall, cake, and channels was also determined as part of the data analysis process in order to check on the integrity of the data and to correct input data to be consistent with the expected trends of the resistance based on the engine conditions used in the test procedure. It was determined that increasing the percent biodiesel content in the test fuel tends to increase the PM reaction rate and the regeneration efficiency of fuel dosing, i.e., at a constant CPF inlet temperature, B20 test fuel resulted in the highest PM reaction rate and regeneration efficiency of fuel dosing. Increasing the CPF inlet temperature also increases PM reaction rate and regeneration efficiency of fuel dosing. Performing active regeneration with B20 as opposed to ULSD allows for a lower CPF temperature to be used to reach the same level of regeneration efficiency, or it allows for a shorter regeneration time at a constant CPF temperature, resulting in decreased fuel consumption for the engine during active regeneration in either scenario.

Patch selection by wintering white-tailed deer : consequences for eastern hemlock regeneration at local and landscape scales

Tsuga canadensis (eastern hemlock) is a highly shade-tolerant, late-successional, and long-lived conifer species found throughout eastern North America. It is most often found in pure or nearly pure stands, because highly acidic and nutrient poor forest floor conditions are thought to favor T. canadensis regeneration while simultaneously limiting the establishment of some hardwood species with greater nutrient requirements. Once a common species, T. canadensis is currently experiencing widescale declines across its range. The hemlock woolly adelgid (Adelges tsugae) is decimating the population across its eastern distribution. Across the Upper Great Lakes region, where the adelgid is currently being held at bay by cold winter temperatures, T. canadensis has been experiencing failures in regeneration attributed, in part, to herbivory by white-tailed deer (Odocoileus virginianus). Deer utilize T. canadensis stands as winter habitat in areas of high snow depth. Tsuga canadensis, once a major component of these forests, currently exists at just a fraction of its pre-settlement abundance due to historic logging and contemporary forest management practices, and what remains is found in small remnant patches surrounded by second- and third-growth deciduous forests. The deer population across the region, however, is likely double that of pre-European settlement times. In this dissertation I explore the relationship between white-tailed deer use of T. canadensis as winter habitat and the effect this use is having on regeneration and forest succession. For this research I quantified stand composition and structure and abiotic variables of elevation and snow depth in 39 randomly selected T. canadensis stands from across the western Upper Peninsula of Michigan. I also quantified composition and the configuration of the landscapes surrounding these stands. I measured relative deer use of T. canadensis stands as pellet group piles deposited in each stand during each of three consecutive winters, 2005-06, 2006-07, and 2007-08. The results of this research suggest that deer use of T. canadensis stands as winter habitat is influenced primarily by snow depth, elevation, and the composition and configuration of the greater landscapes surrounding these stands. Specifically, stands with more heterogeneous landscapes surrounding them (i.e., a patchy mosaic of conifer, deciduous, and open cover) had higher relative deer use than stands surrounded by homogenous deciduous forest cover. Additionally, the intensity of use and the number of stands used was greater in years with higher average snow depth. Tsuga canadensis regeneration in these stands was negatively associated with deer use and Acer saccharum (sugar maple) basal area. Of the 39 stands, 17 and 22 stands had no T. canadensis regeneration in small and large sapling categories, respectively. Acer saccharum was the most common understory tree species, and the importance of A. saccharum in the understory (stems < 10 cm dbh) of the stands was positively associated with overstory A. saccharum dominance. Tsuga canadensis establishment was associated with high-decay coarse woody debris and moss, and deciduous leaf litter inputs in these stands may be limiting access to these important microsites. Furthermore, A. saccharum is more tolerant to the effects of deer herbivory than T. canadensis, giving A. saccharum a competitive advantage in stands being utilized as winter habitat by deer. My research suggests that limited microsite availability, in conjunction with deer herbivory, may be leading to an erosion in T. canadensis patch stability and an altered successional trajectory toward one of A. saccharum dominance, an alternately stable climax species.

Time-domain modeling of high-frequency electromagnetic wave propagation, overhead wires, and earth

Prediction of radiated fields from transmission lines has not previously been studied from a panoptical power system perspective. The application of BPL technologies to overhead transmission lines would benefit greatly from an ability to simulate real power system environments, not limited to the transmission lines themselves. Presently circuitbased transmission line models used by EMTP-type programs utilize Carson’s formula for a waveguide parallel to an interface. This formula is not valid for calculations at high frequencies, considering effects of earth return currents. This thesis explains the challenges of developing such improved models, explores an approach to combining circuit-based and electromagnetics modeling to predict radiated fields from transmission lines, exposes inadequacies of simulation tools, and suggests methods of extending the validity of transmission line models into very high frequency ranges. Electromagnetics programs are commonly used to study radiated fields from transmission lines. However, an approach is proposed here which is also able to incorporate the components of a power system through the combined use of EMTP-type models. Carson’s formulas address the series impedance of electrical conductors above and parallel to the earth. These equations have been analyzed to show their inherent assumptions and what the implications are. Additionally, the lack of validity into higher frequencies has been demonstrated, showing the need to replace Carson’s formulas for these types of studies. This body of work leads to several conclusions about the relatively new study of BPL. Foremost, there is a gap in modeling capabilities which has been bridged through integration of circuit-based and electromagnetics modeling, allowing more realistic prediction of BPL performance and radiated fields. The proposed approach is limited in its scope of validity due to the formulas used by EMTP-type software. To extend the range of validity, a new set of equations must be identified and implemented in the approach. Several potential methods of implementation have been explored. Though an appropriate set of equations has not yet been identified, further research in this area will benefit from a clear depiction of the next important steps and how they can be accomplished. Prediction of radiated fields from transmission lines has not previously been studied from a panoptical power system perspective. The application of BPL technologies to overhead transmission lines would benefit greatly from an ability to simulate real power system environments, not limited to the transmission lines themselves. Presently circuitbased transmission line models used by EMTP-type programs utilize Carson’s formula for a waveguide parallel to an interface. This formula is not valid for calculations at high frequencies, considering effects of earth return currents. This thesis explains the challenges of developing such improved models, explores an approach to combining circuit-based and electromagnetics modeling to predict radiated fields from transmission lines, exposes inadequacies of simulation tools, and suggests methods of extending the validity of transmission line models into very high frequency ranges. Electromagnetics programs are commonly used to study radiated fields from transmission lines. However, an approach is proposed here which is also able to incorporate the components of a power system through the combined use of EMTP-type models. Carson’s formulas address the series impedance of electrical conductors above and parallel to the earth. These equations have been analyzed to show their inherent assumptions and what the implications are. Additionally, the lack of validity into higher frequencies has been demonstrated, showing the need to replace Carson’s formulas for these types of studies. This body of work leads to several conclusions about the relatively new study of BPL. Foremost, there is a gap in modeling capabilities which has been bridged through integration of circuit-based and electromagnetics modeling, allowing more realistic prediction of BPL performance and radiated fields. The proposed approach is limited in its scope of validity due to the formulas used by EMTP-type software. To extend the range of validity, a new set of equations must be identified and implemented in the approach. Several potential methods of implementation have been explored. Though an appropriate set of equations has not yet been identified, further research in this area will benefit from a clear depiction of the next important steps and how they can be accomplished.

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.