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.

Three dimensional finite element simulation of polymer melting and flow in a single-screw extruder : optimization of screw channel geometry

Single-screw extrusion is one of the widely used processing methods in plastics industry, which was the third largest manufacturing industry in the United States in 2007 [5]. In order to optimize the single-screw extrusion process, tremendous efforts have been devoted for development of accurate models in the last fifty years, especially for polymer melting in screw extruders. This has led to a good qualitative understanding of the melting process; however, quantitative predictions of melting from various models often have a large error in comparison to the experimental data. Thus, even nowadays, process parameters and the geometry of the extruder channel for the single-screw extrusion are determined by trial and error. Since new polymers are developed frequently, finding the optimum parameters to extrude these polymers by trial and error is costly and time consuming. In order to reduce the time and experimental work required for optimizing the process parameters and the geometry of the extruder channel for a given polymer, the main goal of this research was to perform a coordinated experimental and numerical investigation of melting in screw extrusion. In this work, a full three-dimensional finite element simulation of the two-phase flow in the melting and metering zones of a single-screw extruder was performed by solving the conservation equations for mass, momentum, and energy. The only attempt for such a three-dimensional simulation of melting in screw extruder was more than twenty years back. However, that work had only a limited success because of the capability of computers and mathematical algorithms available at that time. The dramatic improvement of computational power and mathematical knowledge now make it possible to run full 3-D simulations of two-phase flow in single-screw extruders on a desktop PC. In order to verify the numerical predictions from the full 3-D simulations of two-phase flow in single-screw extruders, a detailed experimental study was performed. This experimental study included Maddock screw-freezing experiments, Screw Simulator experiments and material characterization experiments. Maddock screw-freezing experiments were performed in order to visualize the melting profile along the single-screw extruder channel with different screw geometry configurations. These melting profiles were compared with the simulation results. Screw Simulator experiments were performed to collect the shear stress and melting flux data for various polymers. Cone and plate viscometer experiments were performed to obtain the shear viscosity data which is needed in the simulations. An optimization code was developed to optimize two screw geometry parameters, namely, screw lead (pitch) and depth in the metering section of a single-screw extruder, such that the output rate of the extruder was maximized without exceeding the maximum temperature value specified at the exit of the extruder. This optimization code used a mesh partitioning technique in order to obtain the flow domain. The simulations in this flow domain was performed using the code developed to simulate the two-phase flow in single-screw extruders.

Structural response and dynamics of fluid-structure-control interaction in wind turbine blades

Reducing the uncertainties related to blade dynamics by the improvement of the quality of numerical simulations of the fluid structure interaction process is a key for a breakthrough in wind-turbine technology. A fundamental step in that direction is the implementation of aeroelastic models capable of capturing the complex features of innovative prototype blades, so they can be tested at realistic full-scale conditions with a reasonable computational cost. We make use of a code based on a combination of two advanced numerical models implemented in a parallel HPC supercomputer platform: First, a model of the structural response of heterogeneous composite blades, based on a variation of the dimensional reduction technique proposed by Hodges and Yu. This technique has the capacity of reducing the geometrical complexity of the blade section into a stiffness matrix for an equivalent beam. The reduced 1-D strain energy is equivalent to the actual 3-D strain energy in an asymptotic sense, allowing accurate modeling of the blade structure as a 1-D finite-element problem. This substantially reduces the computational effort required to model the structural dynamics at each time step. Second, a novel aerodynamic model based on an advanced implementation of the BEM(Blade ElementMomentum) Theory; where all velocities and forces are re-projected through orthogonal matrices into the instantaneous deformed configuration to fully include the effects of large displacements and rotation of the airfoil sections into the computation of aerodynamic forces. This allows the aerodynamic model to take into account the effects of the complex flexo-torsional deformation that can be captured by the more sophisticated structural model mentioned above. In this thesis we have successfully developed a powerful computational tool for the aeroelastic analysis of wind-turbine blades. Due to the particular features mentioned above in terms of a full representation of the combined modes of deformation of the blade as a complex structural part and their effects on the aerodynamic loads, it constitutes a substantial advancement ahead the state-of-the-art aeroelastic models currently available, like the FAST-Aerodyn suite. In this thesis, we also include the results of several experiments on the NREL-5MW blade, which is widely accepted today as a benchmark blade, together with some modifications intended to explore the capacities of the new code in terms of capturing features on blade-dynamic behavior, which are normally overlooked by the existing aeroelastic models.

Analysis of a novel waste heat recovery mechanism for an I.C. engine

Typical internal combustion engines lose about 75% of the fuel energy through the engine coolant, exhaust and surface radiation. Most of the heat generated comes from converting the chemical energy in the fuel to mechanical energy and in turn thermal energy is produced. In general, the thermal energy is unutilized and thus wasted. This report describes the analysis of a novel waste heat recovery (WHR) system that operates on a Rankine cycle. This novel WHR system consists of a second piston within the existing piston to reduce losses associated with compression and exhaust strokes in a four-cycle engine. The wasted thermal energy recovered from the coolant and exhaust systems generate a high temperature and high pressure working fluid which is used to power the modified piston assembly. Cycle simulation shows that a large, stationary natural gas spark ignition engine produces enough waste heat to operate the novel WHR system. With the use of this system, the stationary gas compression ignition engine running at 900 RPM and full load had a net increase of 177.03 kW (240.7 HP). This increase in power improved the brake fuel conversion efficiency by 4.53%.

Effect of segmented anodes on the performance and plume of a hall thruster

Development of alternative propellants for Hall thruster operation is an active area of research. Xenon is the current propellant of choice for Hall thrusters, but can be costly in large thrusters and for extended test periods. Condensible propellants may offer an alternative to xenon, as they will not require costly active pumping to remove from a test facility, and may be less expensive to purchase. A method has been developed which uses segmented electrodes in the discharge channel of a Hall thruster to divert discharge current to and from the main anode and thus control the anode temperature. By placing a propellant reservoir in the anode, the evaporation rate, and hence, mass flow of propellant can be controlled. Segmented electrodes for thermal control of a Hall thruster represent a unique strategy of thruster design, and thus the performance of the thruster must be measured to determine the effect the electrodes have on the thruster. Furthermore, the source of any changes in thruster performance due to the adjustment of discharge current between the shims and the main anode must be characterized. A Hall thruster was designed and constructed with segmented electrodes. It was then tested at anode voltages between 300 and 400 V and mass flows between 4 and 6 mg/s, as well as 100%, 75%, 50%, 25%, and <5% of the discharge current on the shim electrodes. The level of current on the shims was adjusted by changing the shim voltage. At each operating point, the thruster performance, plume divergence, ion energy, and multiply charged ion fraction were measured performance exhibited a small change with the level of discharge current on the shim electrodes. Thrust and specific impulse increased by as much as 6% and 7.7%, respectively, as discharge current was shifted from the main anode to the shims at constant anode voltage. Thruster efficiency did not change. Plume divergence was reduced by approximately 4 degrees of half-angle at high levels of current on the shims and at all combinations of mass flow and anode voltage. The fraction of singly charged xenon in the thruster plume varied between approximately 80% and 95% as the anode voltage and mass flow were changed, but did not show a significant change with shim current. Doubly and triply charged xenon made up the remainder of the ions detected. Ion energy exhibited a mixed behavior. The highest voltage present in the thruster largely dictated the most probable energy; either shim or anode voltage, depending on which was higher. The overall change in most probable ion energy was 20-30 eV, the majority of which took place while the shim voltage was higher than the anode voltage. The thrust, specific impulse, plume divergence, and ion energy all indicate that the thruster is capable of a higher performance output at high levels of discharge current on the shims. The lack of a change in efficiency and fraction of multiply charged ions indicate that the thruster can be operated at any level of current on the shims without detrimental effect, and thus a condensible propellant thruster can control the anode temperature without a decrease in efficiency or a change in the multiply charged ion fraction.

Micro combined heat and power laboratory development

Micro Combined Heat and Power (Micro-CHP) system produces both electricity and heat required for residential or small business applications. Use of Micro-CHP in a residential application not only creates energy and economic savings but also reduces the carbon foot print of the house or small business. Additionally, micro-CHP can subsidize its cost of operation by selling excess electricity produced back to the grid. Even though Micro-CHP remains attractive on paper, high initial cost and optimization issues in residential scale heat and electrical requirement has kept this technology from becoming a success. To understand and overcome all disadvantages posed my Micro-CHP system, a laboratory is developed to test different scenarios of Micro-CHP applications so that we can learn and improve the current technology. This report focuses on the development of this Micro-CHP laboratory including installation of Ecopower micro-CHP unit, developing fuel line and exhaust line for Ecopower unit, design of electrical and thermal loop, installing all the instrumentation required for data collection on the Ecopower unit and developing controls for heat load simulation using thermal loop. Also a simulation of Micro-CHP running on Syngas is done in Matlab. This work was supported through the donation of ‘Ecopower’ a Micro-CHP unit by Marathon Engine and through the support of Michigan Tech REF-IF grand.

Mechanical modeling of Intraneural Ganglion Cyst

Intraneural Ganglion Cyst is a 200 year old mystery related to nerve injury which is yet to be solved. Current treatments for the above problem are relatively simple procedures related to removal of cystic contents from the nerve. However, these treatments may result into neuropathic pain and recurrence of the cyst. The articular theory proposed by Spinner et al., (Spinner et al. 2003) takes into consideration the neurological deficit in Common Peroneal Nerve (CPN) branch of the sciatic nerve and affirms that in addition to the above treatments, ligation of articular branch results into foolproof eradication of the deficit. Mechanical Modeling of the Affected Nerve Cross Section will reinforce the articular theory (Spinner et al. 2003). As the cyst propagates, it compresses the neighboring fascicles and the nerve cross section appears like a signet ring. Hence, in order to mechanically model the affected nerve cross section; computational methods capable of modeling excessively large deformations are required. Traditional FEM produces distorted elements while modeling such deformations, resulting into inaccuracies and premature termination of the analysis. The methods described in this Master’s Thesis are effective enough to be able to simulate such deformations. The results obtained from the model adequately resemble the MRI image obtained at the same location and shows an appearance of a signet ring. This Master’s Thesis describes the neurological deficit in brief followed by detail explanation of the advanced computational methods used to simulate this problem. Finally, qualitative results show the resemblance of mechanical model to MRI images of the Nerve Cross Section at the same location validating the capability of these methods to study this neurological deficit.

Static and dynamic contact angle measurement on rough surfaces using sessile drop profile analysis with application to water management in low temperature fuel cells

Fuel Cells are a promising alternative energy technology. One of the biggest problems that exists in fuel cell is that of water management. A better understanding of wettability characteristics in the fuel cells is needed to alleviate the problem of water management. Contact angle data on gas diffusion layers (GDL) of the fuel cells can be used to characterize the wettability of GDL in fuel cells. A contact angle measurement program has been developed to measure the contact angle of sessile drops from drop images. Digitization of drop images induces pixel errors in the contact angle measurement process. The resulting uncertainty in contact angle measurement has been analyzed. An experimental apparatus has been developed for contact angle measurements at different temperature, with the feature to measure advancing and receding contact angles on gas diffusion layers of fuel cells.

Numerical modeling of the failure mechanisms in Si thin film anode for Li-ion batteries

In recent times, the demand for the storage of electrical energy has grown rapidly for both static applications and the portable electronics enforcing the substantial improvement in battery systems, and Li-ion batteries have been proven to have maximum energy storage density in all rechargeable batteries. However, major breakthroughs are required to consummate the requirement of higher energy density with lower cost to penetrate new markets. Graphite anode having limited capacity has become a bottle neck in the process of developing next generation batteries and can be replaced by higher capacity metals such as Silicon. In the present study we are focusing on the mechanical behavior of the Si-thin film anode under various operating conditions. A numerical model is developed to simulate the intercalation induced stress and the failure mechanism of the complex anode structure. Effect of the various physical phenomena such as diffusion induced stress, plasticity and the crack propagation are investigated to predict better performance parameters for improved design.

Development of an optimization model for biofuel facility size and location and a simulation model for design of a biofuel supply chain

To mitigate greenhouse gas (GHG) emissions and reduce U.S. dependence on imported oil, the United States (U.S.) is pursuing several options to create biofuels from renewable woody biomass (hereafter referred to as “biomass”). Because of the distributed nature of biomass feedstock, the cost and complexity of biomass recovery operations has significant challenges that hinder increased biomass utilization for energy production. To facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization and tapping unused forest residues, it is proposed to develop biofuel supply chain models based on optimization and simulation approaches. The biofuel supply chain is structured around four components: biofuel facility locations and sizes, biomass harvesting/forwarding, transportation, and storage. A Geographic Information System (GIS) based approach is proposed as a first step for selecting potential facility locations for biofuel production from forest biomass based on a set of evaluation criteria, such as accessibility to biomass, railway/road transportation network, water body and workforce. The development of optimization and simulation models is also proposed. The results of the models will be used to determine (1) the number, location, and size of the biofuel facilities, and (2) the amounts of biomass to be transported between the harvesting areas and the biofuel facilities over a 20-year timeframe. The multi-criteria objective is to minimize the weighted sum of the delivered feedstock cost, energy consumption, and GHG emissions simultaneously. Finally, a series of sensitivity analyses will be conducted to identify the sensitivity of the decisions, such as the optimal site selected for the biofuel facility, to changes in influential parameters, such as biomass availability and transportation fuel price. Intellectual Merit The proposed research will facilitate the exploration of a wide variety of conditions that promise profitable biomass utilization in the renewable biofuel industry. The GIS-based facility location analysis considers a series of factors which have not been considered simultaneously in previous research. Location analysis is critical to the financial success of producing biofuel. The modeling of woody biomass supply chains using both optimization and simulation, combing with the GIS-based approach as a precursor, have not been done to date. The optimization and simulation models can help to ensure the economic and environmental viability and sustainability of the entire biofuel supply chain at both the strategic design level and the operational planning level. Broader Impacts The proposed models for biorefineries can be applied to other types of manufacturing or processing operations using biomass. This is because the biomass feedstock supply chain is similar, if not the same, for biorefineries, biomass fired or co-fired power plants, or torrefaction/pelletization operations. Additionally, the research results of this research will continue to be disseminated internationally through publications in journals, such as Biomass and Bioenergy, and Renewable Energy, and presentations at conferences, such as the 2011 Industrial Engineering Research Conference. For example, part of the research work related to biofuel facility identification has been published: Zhang, Johnson and Sutherland [2011] (see Appendix A). There will also be opportunities for the Michigan Tech campus community to learn about the research through the Sustainable Future Institute.

Modeling mechanical response of heterogeneous materials

Heterogeneous materials are ubiquitous in nature and as synthetic materials. These materials provide unique combination of desirable mechanical properties emerging from its heterogeneities at different length scales. Future structural and technological applications will require the development of advanced light weight materials with superior strength and toughness. Cost effective design of the advanced high performance synthetic materials by tailoring their microstructure is the challenge facing the materials design community. Prior knowledge of structure-property relationships for these materials is imperative for optimal design. Thus, understanding such relationships for heterogeneous materials is of primary interest. Furthermore, computational burden is becoming critical concern in several areas of heterogeneous materials design. Therefore, computationally efficient and accurate predictive tools are highly essential. In the present study, we mainly focus on mechanical behavior of soft cellular materials and tough biological material such as mussel byssus thread. Cellular materials exhibit microstructural heterogeneity by interconnected network of same material phase. However, mussel byssus thread comprises of two distinct material phases. A robust numerical framework is developed to investigate the micromechanisms behind the macroscopic response of both of these materials. Using this framework, effect of microstuctural parameters has been addressed on the stress state of cellular specimens during split Hopkinson pressure bar test. A voronoi tessellation based algorithm has been developed to simulate the cellular microstructure. Micromechanisms (microinertia, microbuckling and microbending) governing macroscopic behavior of cellular solids are investigated thoroughly with respect to various microstructural and loading parameters. To understand the origin of high toughness of mussel byssus thread, a Genetic Algorithm (GA) based optimization framework has been developed. It is found that two different material phases (collagens) of mussel byssus thread are optimally distributed along the thread. These applications demonstrate that the presence of heterogeneity in the system demands high computational resources for simulation and modeling. Thus, Higher Dimensional Model Representation (HDMR) based surrogate modeling concept has been proposed to reduce computational complexity. The applicability of such methodology has been demonstrated in failure envelope construction and in multiscale finite element techniques. It is observed that surrogate based model can capture the behavior of complex material systems with sufficient accuracy. The computational algorithms presented in this thesis will further pave the way for accurate prediction of macroscopic deformation behavior of various class of advanced materials from their measurable microstructural features at a reasonable computational cost.

Development of a direct evaporation bismuth Hall thruster

Hall thrusters have been under active development around the world since the 1960’s. Thrusters using traditional propellants such as xenon have been flown on a variety of satellite orbit raising and maintenance missions with an excellent record. To expand the mission envelope, it is necessary to lower the specific impulse of the thrusters but xenon and krypton are poor performers at specific impulses below 1,200 seconds. To enhance low specific impulse performance, this dissertation examines the development of a Hall-effect thruster which uses bismuth as a propellant. Bismuth, the heaviest non-radioactive element, holds many advantages over noble gas propellants from an energetics as well as a practical economic standpoint. Low ionization energy, large electron-impact crosssection and high atomic mass make bismuth ideal for low-specific impulse applications. The primary disadvantage lies in the high temperatures which are required to generate the bismuth vapors. Previous efforts carried out in the Soviet Union relied upon the complete bismuth vaporization and gas phase delivery to the anode. While this proved successful, the power required to vaporize and maintain gas phase throughout the mass flow system quickly removed many of the efficiency gains expected from using bismuth. To solve these problems, a unique method of delivering liquid bismuth to the anode has been developed. Bismuth is contained within a hollow anode reservoir that is capped by a porous metallic disc. By utilizing the inherent waste heat generated in a Hall thruster, liquid bismuth is evaporated and the vapors pass through the porous disc into the discharge chamber. Due to the high temperatures and material compatibility requirements, the anode was fabricated out of pure molybdenum. The porous vaporizer was not available commercially so a method of creating a refractory porous plate with 40-50% open porosity was developed. Molybdenum also does not respond well to most forms of welding so a diffusion bonding process was also developed to join the molybdenum porous disc to the molybdenum anode. Operation of the direct evaporation bismuth Hall thruster revealed interesting phenomenon. By utilizing constant current mode on a discharge power supply, the discharge voltage settles out to a stable operating point which is a function of discharge current, anode face area and average pore size on the vaporizer. Oscillations with a 40 second period were also observed. Preliminary performance data suggests that the direct evaporation bismuth Hall thruster performs similar to xenon and krypton Hall thrusters. Plume interrogation with a Retarding Potential Analyzer confirmed that bismuth ions were being efficiently accelerated while Faraday probe data gave a view of the ion density in the exhausted plume.

Laser Thomson scattering measurements of electron temperature and density in a Hall-effect plasma

Hall-effect thrusters (HETs) are compact electric propulsion devices with high specific impulse used for a variety of space propulsion applications. HET technology is well developed but the electron properties in the discharge are not completely understood, mainly due to the difficulty involved in performing accurate measurements in the discharge. Measurements of electron temperature and density have been performed using electrostatic probes, but presence of the probes can significantly disrupt thruster operation, and thus alter the electron temperature and density. While fast-probe studies have expanded understanding of HET discharges, a non-invasive method of measuring the electron temperature and density in the plasma is highly desirable. An alternative to electrostatic probes is a non-perturbing laser diagnostic technique that measures Thomson scattering from the plasma. Thomson scattering is the process by which photons are elastically scattered from the free electrons in a plasma. Since the electrons have thermal energy their motion causes a Doppler shift in the scattered photons that is proportional to their velocity. Like electrostatic probes, laser Thomson scattering (LTS) can be used to determine the temperature and density of free electrons in the plasma. Since Thomson scattering measures the electron velocity distribution function directly no assumptions of the plasma conditions are required, allowing accurate measurements in anisotropic and non-Maxwellian plasmas. LTS requires a complicated measurement apparatus, but has the potential to provide accurate, non-perturbing measurements of electron temperature and density in HET discharges. In order to assess the feasibility of LTS diagnostics on HETs non-invasive measurements of electron temperature and density in the near-field plume of a Hall thruster were performed using a custom built laser Thomson scattering diagnostic. Laser measurements were processed using a maximum likelihood estimation method and results were compared to conventional electrostatic double probe measurements performed at the same thruster conditions. Electron temperature was found to range from approximately 1 – 40 eV and density ranged from approximately 1.0 x 1017 m-3 to 1.3 x 1018 m-3 over discharge voltages from 250 to 450 V and mass flow rates of 40 to 80 SCCM using xenon propellant.

Quantum dot / optical protein bio-nano hybrid system biosensing

The integration of novel nanomaterials with highly-functional biological molecules has advanced multiple fields including electronics, sensing, imaging, and energy harvesting. This work focuses on the creation of a new type of bio-nano hybrid substrate for military biosensing applications. Specifically it is shown that the nano-scale interactions of the optical protein bacteriorhodopsin and colloidal semiconductor quantum dots can be utilized as a generic sensing substrate. This work spans from the basic creation of the protein to its application in a novel biosensing system. The functionality of this sensor design originates from the unique interactions between the quantum dot and bacteriorhodopsin molecule when in nanoscale proximity. A direct energy transfer relationship has been established between coreshell quantum dots and the optical protein bacteriorhodopsin that substantially enhances the protein’s native photovoltaic capabilities. This energy transfer phenomena is largely distance dependent, in the sub-10nm realm, and is characterized experimentally at multiple separation distances. Experimental results on the energy transfer efficiency in this hybrid system correlate closely to theoretical predictions. Deposition of the hybrid system with nano-scale control has allowed for the utilization of this energy transfer phenomena as a modulation point for a functional biosensor prototype. This work reveals that quantum dots have the ability to activate the bacteriorhodopsin photocycle through both photonic and non-photonic energy transfer mechanisms. By altering the energy transferred to the bacteriorhodopsin molecule from the quantum dot, the electrical output of the protein can be modulated. A biosensing prototype was created in which the energy transfer relationship is altered upon target binding, demonstrating the applicability of a quantum dot/bacteriorhodopsin hybrid system for sensor applications. The electrical nature of this sensing substrate will allow for its efficient integration into a nanoelectronics array form, potentially leading to a small-low power sensing platform for remote toxin detection applications.

Water transport in complex, non-wetting porous layers with applications to water management in low temperature fuel cells

An experimental setup was designed to visualize water percolation inside the porous transport layer, PTL, of proton exchange membrane, PEM, fuel cells and identify the relevant characterization parameters. In parallel with the observation of the water movement, the injection pressure (pressure required to transport water through the PTL) was measured. A new scaling for the drainage in porous media has been proposed based on the ratio between the input and the dissipated energies during percolation. A proportional dependency was obtained between the energy ratio and a non-dimensional time and this relationship is not dependent on the flow regime; stable displacement or capillary fingering. Experimental results show that for different PTL samples (from different manufacturers) the proportionality is different. The identification of this proportionality allows a unique characterization of PTLs with respect to water transport. This scaling has relevance in porous media flows ranging far beyond fuel cells. In parallel with the experimental analysis, a two-dimensional numerical model was developed in order to simulate the phenomena observed in the experiments. The stochastic nature of the pore size distribution, the role of the PTL wettability and morphology properties on the water transport were analyzed. The effect of a second porous layer placed between the porous transport layer and the catalyst layer called microporous layer, MPL, was also studied. It was found that the presence of the MPL significantly reduced the water content on the PTL by enhancing fingering formation. Moreover, the presence of small defects (cracks) within the MPL was shown to enhance water management. Finally, a corroboration of the numerical simulation was carried out. A threedimensional version of the network model was developed mimicking the experimental conditions. The morphology and wettability of the PTL are tuned to the experiment data by using the new energy scaling of drainage in porous media. Once the fit between numerical and experimental data is obtained, the computational PTL structure can be used in different types of simulations where the conditions are representative of the fuel cell operating conditions.