THERMAL CHARACTERIZATION OF COMBUSTION CHAMBER COMPONENTS IN A GASOLINE TURBOCHARGED DIRECT INJECTION (GTDI) ENGINE

This report is a PhD dissertation proposal to study the in-cylinder temperature and heat flux distributions within a gasoline turbocharged direct injection (GTDI) engine. Recent regulations requiring automotive manufacturers to increase the fuel efficiency of their vehicles has led to great technological achievements in internal combustion engines. These achievements have increased the power density of gasoline engines dramatically in the last two decades. Engine technologies such as variable valve timing (VVT), direct injection (DI), and turbocharging have significantly improved engine power-to-weight and power-to-displacement ratios. A popular trend for increasing vehicle fuel economy in recent years has been to downsize the engine and add VVT, DI, and turbocharging technologies so that a lighter more efficient engine can replace a larger, heavier one. With the added power density, thermal management of the engine becomes a more important issue. Engine components are being pushed to their temperature limits. Therefore it has become increasingly important to have a greater understanding of the parameters that affect in-cylinder temperatures and heat transfer. The proposed research will analyze the effects of engine speed, load, relative air-fuel ratio (AFR), and exhaust gas recirculation (EGR) on both in-cylinder and global temperature and heat transfer distributions. Additionally, the effect of knocking combustion and fuel spray impingement will be investigated. The proposed research will be conducted on a 3.5 L six cylinder GTDI engine. The research engine will be instrumented with a large number of sensors to measure in-cylinder temperatures and pressures, as well as, the temperature, pressure, and flow rates of energy streams into and out of the engine. One of the goals of this research is to create a model that will predict the energy distribution to the crankshaft, exhaust, and cooling system based on normalized values for engine speed, load, AFR, and EGR. The results could be used to aid in the engine design phase for turbocharger and cooling system sizing. Additionally, the data collected can be used for validation of engine simulation models, since in-cylinder temperature and heat flux data is not readily available in the literature..

Fabrication of nanoparticles by short-distance sputter deposition

During the past decades, tremendous research interests have been attracted to investigate nanoparticles due to their promising catalytic, magnetic, and optical properties. In this thesis, two novel methods of nanoparticle fabrication were introduced and the basic formation mechanisms were studied. Metal nanoparticles and polyurethane nanoparticles were separately fabricated by a short-distance sputter deposition technique and a reactive ion etching process. First, a sputter deposition method with a very short target-substrate distance is found to be able to generate metal nanoparticles on the glass substrate inside a RIE chamber. The distribution and morphology of nanoparticles are affected by the distance, the ion concentration and the process time. Densely-distributed nanoparticles of various compositions are deposited on the substrate surface when the target-substrate distance is smaller than 130mm. It is much less than the atoms’ mean free path, which is the threshold in previous research for nanoparticles’ formation. Island structures are formed when the distance is increased to 510mm, indicating the tendency to form continuous thin film. The trend is different from previously-reported sputtering method for nanoparticle fabrication, where longer distance between the target and the substrate facilitates the formation of nanoparticle. A mechanism based on the seeding effect of the substrate is proposed to interpret the experimental results. Secondly, in polyurethane nanoparticles’ fabrication, a mechanism is put forward based on the microphase separation phenomenon in block copolymer thin film. The synthesized polymers have formed dispersed and continuous phases because of the different properties between segments. With harder mechanical property, the dispersed phase is remained after RIE process while the continuous phase is etched away, leading to the formation of nanoparticles on the substrate. The nanoparticles distribution is found to be affected by the heating effect, the process time and the plasma power. Superhydrophilic property is found on samples with these two types of nanoparticles. The relationship between the nanostructure and the hydrophilicity is studied for further potential applications.

Finite element analysis on Titan STL3 OTR tires

Finite element tire modeling can be a challenging process, due to the overall complexities within the tire and the many variables that are required to produce capable predictive simulations. Utilizing tools from Abaqus finite element software, adequate predictive simulations that represent actual operational conditions can be made possible. Many variables that result from complex geometries and materials, multiple loading conditions, and surface contact can be incorporated into modeling simulations. This thesis outlines modeling practices used to conduct analysis on specific tire variants of the STL3 series OTR tire line, produced by Titan Tire. Finite element models were created to represent an inflated tire and rim assembly, supporting a 30,000 lb load while resting on a flat surface. Simulations were conducted with reinforcement belt cords at variable angles in order to understand how belt cord arrangement affects tire components and stiffness response.

Kalman Filter approaches on crane swing

The purpose of this study is to explore a Kalman Filter approach to estimating swing of crane-suspended loads. Measuring real-time swing is needed to implement swing damping control strategies where crane joints are used to remove energy from a swinging load. The typical solution to measuring swing uses an inertial sensor attached to the hook block. Measured hook block twist is used to resolve the other two sensed body rates into tangential and radial swing. Uncertainty in the twist measurement leads to inaccurate tangential and radial swing calculations and ineffective swing damping. A typical mitigation approach is to bandpass the inertial sensor readings to remove low frequency drift and high frequency noise. The center frequency of the bandpass filter is usually designed to track the load length and the pass band width set to trade off performance with damping loop gain. The Kalman Filter approach developed here allows all swing motions (radial, tangential and twist) to be measured without the use of a bandpass filter. This provides an alternate solution for swing damping control implementation. After developing a Kalman Filter solution for a two-dimensional swing scenario, the three-dimensional system is considered where simplifying assumptions, suggested by the two-dimensional study, are exploited. One of the interesting aspects of the three-dimensional study is the hook block twist model. Unlike the mass-independence of a pendulum's natural frequency, the twist natural frequency depends both on the pendulum length and the load’s mass distribution. The linear Kalman Filter is applied to experimental data demonstrating the ability to extract the individual swing components for complex motions. It should be noted that the three-dimensional simplifying assumptions preclude the ability to measure two "secondary" hook block rotations. The ability to segregate these motions from the primary swing degrees of freedom was illustrated in the two-dimensional study and could be included into the three-dimensional solution if they were found to be important for a particular application.

DEVELOPMENT OF A PREDICTIVE COMBUSTION MODEL OF A SPARK IGNITED ENGINE WITH GASOLINE DIRECT INJECTION, VARIABLE VALVE TIMING, DURATION AND LIFT TECHNOLOGIES

There is a need by engine manufactures for computationally efficient and accurate predictive combustion modeling tools for integration in engine simulation software for the assessment of combustion system hardware designs and early development of engine calibrations. This thesis discusses the process for the development and validation of a combustion modeling tool for Gasoline Direct Injected Spark Ignited Engine with variable valve timing, lift and duration valvetrain hardware from experimental data. Data was correlated and regressed from accepted methods for calculating the turbulent flow and flame propagation characteristics for an internal combustion engine. A non-linear regression modeling method was utilized to develop a combustion model to determine the fuel mass burn rate at multiple points during the combustion process. The computational fluid dynamic software Converge ©, was used to simulate and correlate the 3-D combustion system, port and piston geometry to the turbulent flow development within the cylinder to properly predict the experimental data turbulent flow parameters through the intake, compression and expansion processes. The engine simulation software GT-Power © is then used to determine the 1-D flow characteristics of the engine hardware being tested to correlate the regressed combustion modeling tool to experimental data to determine accuracy. The results of the combustion modeling tool show accurate trends capturing the combustion sensitivities to turbulent flow, thermodynamic and internal residual effects with changes in intake and exhaust valve timing, lift and duration.

Analysis of in-cylinder pressure transducer data quality utilizing a SIDI turbocharged engine

In-cylinder pressure transducers have been used for decades to record combustion pressure inside a running engine. However, due to the extreme operating environment, transducer design and installation must be considered in order to minimize measurement error. One such error is caused by thermal shock, where the pressure transducer experiences a high heat flux that can distort the pressure transducer diaphragm and also change the crystal sensitivity. This research focused on investigating the effects of thermal shock on in-cylinder pressure transducer data quality using a 2.0L, four-cylinder, spark-ignited, direct-injected, turbo-charged GM engine. Cylinder four was modified with five ports to accommodate pressure transducers of different manufacturers. They included an AVL GH14D, an AVL GH15D, a Kistler 6125C, and a Kistler 6054AR. The GH14D, GH15D, and 6054AR were M5 size transducers. The 6125C was a larger, 6.2mm transducer. Note that both of the AVL pressure transducers utilized a PH03 flame arrestor. Sweeps of ignition timing (spark sweep), engine speed, and engine load were performed to study the effects of thermal shock on each pressure transducer. The project consisted of two distinct phases which included experimental engine testing as well as simulation using a commercially available software package. A comparison was performed to characterize the quality of the data between the actual cylinder pressure and the simulated results. This comparison was valuable because the simulation results did not include thermal shock effects. All three sets of tests showed the peak cylinder pressure was basically unaffected by thermal shock. Comparison of the experimental data with the simulated results showed very good correlation. The spark sweep was performed at 1300 RPM and 3.3 bar NMEP and showed that the differences between the simulated results (no thermal shock) and the experimental data for the indicated mean effective pressure (IMEP) and the pumping mean effective pressure (PMEP) were significantly less than the published accuracies. All transducers had an IMEP percent difference less than 0.038% and less than 0.32% for PMEP. Kistler and AVL publish that the accuracy of their pressure transducers are within plus or minus 1% for the IMEP (AVL 2011; Kistler 2011). In addition, the difference in average exhaust absolute pressure between the simulated results and experimental data was the greatest for the two Kistler pressure transducers. The location and lack of flame arrestor are believed to be the cause of the increased error. For the engine speed sweep, the torque output was held constant at 203 Nm (150 ft-lbf) from 1500 to 4000 RPM. The difference in IMEP was less than 0.01% and the PMEP was less than 1%, except for the AVL GH14D which was 5% and the AVL GH15DK which was 2.25%. A noticeable error in PMEP appeared as the load increased during the engine speed sweeps, as expected. The load sweep was conducted at 2000 RPM over a range of NMEP from 1.1 to 14 bar. The difference in IMEP values were less 0.08% while the PMEP values were below 1% except for the AVL GH14D which was 1.8% and the AVL GH15DK which was at 1.25%. In-cylinder pressure transducer data quality was effectively analyzed using a combination of experimental data and simulation results. Several criteria can be used to investigate the impact of thermal shock on data quality as well as determine the best location and thermal protection for various transducers.

Accelerometer based measurements of combustion in an automotive turbocharged diesel engine

The capability to detect combustion in a diesel engine has the potential of being an important control feature to meet increasingly stringent emission regulations, develop alternative combustion strategies, and use of biofuels. In this dissertation, block mounted accelerometers were investigated as potential feedback sensors for detecting combustion characteristics in a high-speed, high pressure common rail (HPCR), 1.9L diesel engine. Accelerometers were positioned in multiple placements and orientations on the engine, and engine testing was conducted under motored, single and pilot-main injection conditions. Engine tests were conducted at varying injection timings, engine loads, and engine speeds to observe the resulting time and frequency domain changes of the cylinder pressure and accelerometer signals. The frequency content of the cylinder pressure based signals and the accelerometer signals between 0.5 kHz and 6 kHz indicated a strong correlation with coherence values of nearly 1. The accelerometers were used to produce estimated combustion signals using the Frequency Response Functions (FRF) measured from the frequency domain characteristics of the cylinder pressure signals and the response of the accelerometers attached to the engine block. When compared to the actual combustion signals, the estimated combustion signals produced from the accelerometer response had Root Mean Square Errors (RMSE) between 7% and 25% of the actual signals peak value. Weighting the FRF’s from multiple test conditions along their frequency axis with the coherent output power reduced the median RMSE of the estimated combustion signals and the 95th percentile of RMSE produced from each test condition. The RMSE’s of the magnitude based combustion metrics including peak cylinder pressure, MPG, peak ROHR, and work estimated from the combustion signals produced by the accelerometer responses were between 15% and 50% of their actual value. The MPG measured from the estimated pressure gradient shared a direct relationship to the actual MPG. The location based combustion metrics such as the location of peak values and burn durations were capable of RMSE measurements as low as 0.9°. Overall, accelerometer based combustion sensing system was capable of detecting combustion and providing feedback regarding the in cylinder combustion process

Sustainability study of nanomaterials including societal and occupational implications

In recent years there has been a tremendous amount of research in the area of nanotechnology. History tells us that the commercialization of technologies will always be accompanied by both positive and negative effects for society and the environment. Products containing nanomaterials are already available in the market, and yet there is still not much information regarding the potential negative effects that these products may cause. The work presented in this dissertation describes a holistic approach to address different dimensions of nanotechnology sustainability. Life cycle analysis (LCA) was used to study the potential usage of polyethylene filled with nanomaterials to manufacture automobile body panels. Results showed that the nanocomposite does not provide an environmental benefit over traditional steel panels. A new methodology based on design of experiments (DOE) techniques, coupled with LCA, was implemented to investigate the impact of inventory uncertainties. Results showed that data variability does not have a significant effect on the prediction of the environmental impacts. Material profiles for input materials did have a highly significant effect on the overall impact. Energy consumption and material characterization were identified as two mainstreams where additional research is needed in order to predict the overall impact of nanomaterials more effectively. A study was undertaken to gain insights into the behavior of small particles in contact with a surface exposed to air flow to determine particle lift-off from the surface. A mapping strategy was implemented that allows for the identification of conditions for particle liftoff based on particle size and separation distance from the wall. Main results showed that particles smaller than 0:1mm will not become airborne under shear flow unless the separation distance is greater than 15 nm. Results may be used to minimize exposure to airborne materials. Societal implications that may occur in the workplace were researched. This research task explored different topics including health, ethics, and worker perception with the aim of identifying the base knowledge available in the literature. Recommendations are given for different scenarios to describe how workers and employers could minimize the unwanted effects of nanotechnology production.

Finite element study on fracture patterns in human scaphoid wrist bone due to free fall

Scaphoid is one of the 8 carpal bones found adjacent to the thumb supported proximally by Radius bone. During the free fall, on outstretched hand, the impact load gets transferred to the scaphoid at its free anterior end. Unique arrangement of other carpal bones in the palm is also one of the reasons for the load to get transferred to scaphoid. About half of the total load acting upon carpal bone gets transferred to scaphoid at its distal pole. There are about 10 to 12 clinically observed fracture pattern in the scaphoid due to free fall. The aim of the study is to determine the orientation of the load, magnitude of the load and the corresponding fracture pattern. This study includes both static and dynamic finite element models validated by experiments. The scaphoid model has been prepared from CT scans of a 27 year old person. The 2D slices of the CT scans have been converted to 3D model by using MIMICS software. There are four cases of loading studied which are considered to occur clinically more frequently. In case (i) the load is applied at the posterior end at distal pole whereas in case (ii), (iii) and (iv), the load is applied at anterior end at different directions. The model is given a fixed boundary condition at the region which is supported by Radius bone during the impact. Same loading and boundary conditions have been used in both static and dynamic explicit finite element analysis. The site of fracture initiation and path of fracture propagation have been identified by using max principal stress / gradient and max principal strain / gradient criterion respectively in static and dynamic explicit finite element analysis. Static and dynamic impact experiments were performed on the polyurethane foam specimens to validate the finite element results. Experimental results such as load at fracture, site of fracture initiation and path of fracture propagation have been compared with the results of finite element analysis. Four different types of fracture patterns observed in clinical studies have been identified in this study.

Experimental study of physical aging effects on properties of EPON 862-DETDA

Epoxies find variety of applications and during these applications they get exposed to different conditions like elevated temperatures, hydrothermal, chemical, etc. It has been observed that properties of epoxies do get affected substantially if exposed to these conditions for extended period of time and because of the variety of applications, researchers found it necessary to study their effects on the thermal, mechanical, physical and chemical properties. However in this report the focus is on studying effects of physical aging on mechanical properties of EPON 862 with DETDA as its curing agent, where physical aging is aging is the condition which occurs due to exposure to elevated temperatures. A fair amount of computational work has been performed on EPON 862- DETDA to study the effects of physical aging, however very little known work has been done experimentally to study these effects. Young’s modulus, hardness, failure strength, strain to failure, density and glass transition are the properties which have been obtained using various experimental methods - tensile testing, nanoindentation and differential scanning calorimetry. Experimental work on other epoxies have shown no increase or very slight increase in the Young’s modulus and hardness with increased aging time, also decrease in failure strength and strain to failure and through this work on EPON 862- DETDA we can observe similar trends.

Two dimensional Lattice Boltzmann simulation of fluid flow through an idealized micro-structure of disordered media

This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.

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.

Microfluidic hydrogen generator : featuring passive on-demand gas generation and self-circulated reactant for integration with micro fuel cell

A microfluidic hydrogen generator is presented in this work. Its fabrication, characterization, and integration with a micro proton exchange membrane (PEM) fuel cell are described. Hydrogen gas is generated by the hydrolysis of aqueous ammonia borane. Gas generation, as well as the circulation of ammonia borane from a rechargeable fuel reservoir, is performed without any power consumption. To achieve this, directional growth and selective venting of hydrogen gas is maintained in the microchannels, which results in the circulation of fresh reactant from the fuel reservoir. In addition to this self-circulation mechanism, the hydrogen generator has been demonstrated to self-regulate gas generation to meet demands of a connected micro fuel cell. All of this is done without parasitic power consumption from the fuel cell. Results show its feasibility in applications of high-impedance systems. Lastly, recommendations for improvements and suggestions for future work are described

Biomechanics of meniscal horn attachments

Menisci are anchored to the tibia by means of ligament-like structures called meniscal attachments. Failure material properties of bovine meniscal attachments were obtained. There were no significant differences in the structural properties or ultimate stress between the meniscal attachments (p>0.05). Furthermore, Glycosaminoglycan (GAG) fraction and crimping frequency was obtained for each attachment using histology and differential interference contrast (DIC) respectively. Results showed that the anterior attachment’s insertion had the greatest GAG fraction when compared to the posterior attachment’s insertion. Crimp frequency of the collagen fibrils was homogeneous along the length. Moreover, Scanning Electron Microscopy (SEM) technique was used to reveal the morphology of collagen in human meniscal attachments. Its midsubstance was composed of collagen fascicles running parallel to the longitudinal axis, with a few fibrils running obliquely, and others transversely. There were no differences between attachments for crimping angle or length. Since ligamentous-type tissues are comprised mainly of water, the fluid pressure within meniscal horn attachments was measured using a Fiber Optic Microsensor (FOM). Four cadaveric human joints were subjected to 2BW compressive load (ramp) at 0-, 15-, and 30-degrees of flexion for a minute and then the load was hold for 20 minutes (equilibrium). There were significant differences between 0- and 15- (p1– c5) were obtained. Significant differences were found on the straightened collagen fibers coefficient (c5) between MP and LA attachments (p

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.