Search Results

The thermal efficiency of spark-ignition engines can be enhanced by increasing the rate of exhaust gas recirculation (EGR) such that the low temperature combustion regime could be achieved. However, there is an upper limit on the amount of EGR rate, beyond which flame speed becomes slow and unstable, and local quenching starts to hurt the combustion stability, efficiency, and emission. To resolve this issue, the concept of dedicated EGR has been proposed previously to be an effective way to enhance flame propagation under lean burn condition with even higher levels of EGR with reformate hydrogen and carbon monoxide. In this study, the effects of thermochemical fuel reforming on the reformate composition under rich conditions (1.0 < ϕ < 2.0) have been studied using detailed chemistry for iso-octane, as the representative component for gasoline.

Low vibration and noise level in internal combustion engines has become an essential part of the design process. It is well known that the piston assembly can be a major source of engine mechanical friction and cold start noise, if not designed properly. The piston secondary motion and piston-bore contact pattern are critical in piston design because they affect the skirt-to-bore impact force and therefore, how the piston impact excitation energy is damped, transmitted and eventually radiated from the engine structure as noise. An analytical method is presented in this paper for simulating piston secondary dynamics and piston-bore contact for an asymmetric half piston model. The method includes several important physical attributes such as bore distortion effects due to mechanical and thermal deformation, inertia loading, piston barrelity and ovality, piston flexibility and skirt-to-bore clearance. The method accounts for piston kinematics, rigid-body dynamics and flexibility.

Flame speed data reported in most literature are acquired in conventional apparatus such as the spherical combustion bomb and counterflow burner, and are limited to atmospheric pressure and ambient or slightly elevated unburnt temperatures. As such, these data bear little relevance to internal combustion engines and gas turbines, which operate under typical pressures of 10-50 bar and unburnt temperature up to 900K or higher. These elevated temperatures and pressures not only modify dominant flame chemistry, but more importantly, they inevitably facilitate pre-ignition reactions and hence can change the upstream thermodynamic and chemical conditions of a regular hot flame leading to modified flame properties. This study focuses on how auto-ignition chemistry affects flame propagation, especially in the negative-temperature coefficient (NTC) regime, where dimethyl ether (DME), n-heptane and iso-octane are chosen for study as typical fuels exhibiting low temperature chemistry (LTC).

Tubular hydroforming is a novel process that has recently gained much attention due to its cost-effective application in the automotive industry. Hydroformed automotive parts have high strength to weight ratio and have good repeatability with high dimensional accuracy. At this time, there is little experience in modeling the hydroforming process to better understand its application and researchers have tried using stamping simulation software to analyze the process. Unlike conventional sheet stamping which is a displacement driven process, tubular hydroforming is a force driven process and its success is governed by the nature of internal pressurization. Hence, a new three-dimensional finite element model using a computationally efficient 6-noded shell element has been developed. A simple pressure prediction model has been developed and integrated into the formulation for effective control of the process.

A new method is demonstrated for rating the “severity” of a hypoid gear set duty cycle (revolutions at torque) using the intercept of T-N curve to support gearset selection and sizing decision across vehicle programs. Historically, it has been customary to compute a cumulative damage (using Miner's Rule) for a rotating component duty cycle given a T-N curve slope and intercept for the component and failure mode of interest. The slope and intercept of a T-N curve is often proprietary to the axle manufacturer and are not published. Therefore, for upfront sizing and selection purposes representative T-N properties are used to assess relative component duty cycle severity via cumulative damage (non-dimensional quantity). A similar duty cycle severity rating can also be achieved by computing the intercept of the T-N curve instead of cumulative damage, which is the focus of this study.

Fatigue life estimation, reliability and durability are important in acquisition, maintenance and operation of vehicle systems. Fatigue life is random because of the stochastic load, the inherent variability of material properties, and the uncertainty in the definition of the S-N curve. The commonly used fatigue life estimation methods calculate the mean (not the distribution) of fatigue life under Gaussian loads using the potentially restrictive narrow-band assumption. In this paper, a general methodology is presented to calculate the statistics of fatigue life for a linear vibratory system under stationary, non-Gaussian loads considering the effects of skewness and kurtosis. The input loads are first characterized using their first four moments (mean, standard deviation, skewness and kurtosis) and a correlation structure equivalent to a given Power Spectral Density (PSD).

Real-time measurement or estimation of crank-angle-resolved engine cylinder pressure may become commonplace in the next generation of engine controllers to optimize spark, valve timing, or compression ratio. Toward the development of a real-time cylinder pressure estimator, this work presents a crank-angle-resolved engine cylinder pressure estimation model that could accept inputs such as speed, manifold pressure and throttle position, and deliver crank-angle resolved cylinder pressure in real-time, at engine speeds covering the useful operating range of most engines. The model was validated by comparing simulated cylinder pressure with thirteen sets of cylinder pressure data, from two different commercial engines from two different OEMs. Estimated pressures were compared against the actual measured pressure traces. The average relative error is about 3% while the maximum relative error is 5%. Both can be improved with further tuning.

The use of commercial finite element analysis (FEA) software to perform stamping feasibility studies of automotive components has grown extensively over the last decade. Although product and process engineers have now come to rely heavily on results from FEA simulation for manufacturability of components, the prediction of springback has still not been perfected. Springback prediction for simple geometries is found to be quite accurate while springback prediction in complex components fails to compare with experimental results. Since most forming simulation FEA software uses a dynamic explicit solution method, the choice of various input parameters greatly affects the prediction of post formed stresses in the final component. Accurate stress prediction is critical for determination of springback, therefore this study focuses on the effects of some of the simulation parameters such as, element size, tool/loading speed and loading profile.

Experimental studies of soot particles showed that the intensity ratio of amorphous and graphite layers measured by Raman spectroscopy correlates to soot oxidation reactivities, which is very important for regeneration of the diesel particulate filters and gasoline particulate filters. This physical mechanism is absent in all soot models. In the present paper, a novel two-layer soot model was proposed that considers the amorphous and graphite layers in the soot particles. The soot model considers soot inception, soot surface growth, soot oxidation by O2 and OH, and soot coagulation. It is assumed that amorphous-type soot forms from fullerene. No soot coagulation is considered in the model between the amorphous- and graphitic-types of soot. Benzene is taken as the soot precursor, which is formed from acetylene. The model was implemented into a commercial CFD software CONVERGE using user defined functions. A diesel engine case was simulated.

Spark ignition engines utilize catalytic converters to reform harmful exhaust gas emissions such as carbon monoxide, unburned hydrocarbons, and oxides of nitrogen into less harmful products. Aftertreatment devices require the use of expensive catalytic metals such as platinum, palladium, and rhodium. Meanwhile, tightening automotive emissions regulations globally necessitate the development of high-performance exhaust gas catalysts. So, automotive manufactures must balance maximizing catalyst performance while minimizing production costs. There are thousands of different recipes for catalytic converters, with each having a different effect on the various catalytic chemical reactions which impact the resultant tailpipe gas composition. In the development of catalytic converters, simulation models are often used to reduce the need for physical parts and testing, thus saving significant time and money.

Vehicle weight reduction is a significant challenge for the modern automotive industry. In recent years, the amount of vehicular components constructed from aluminum alloy has increased due to its light weighting capabilities. Automotive manufacturing processes, predominantly those utilizing various stamping applications, require a thorough understanding of aluminum fracture predictions methods, in order to accurately simulate the process using Finite Element Method (FEM) software or use it in automotive engineering manufacture. This paper presents the strain distribution of A5182 aluminum samples after punch impact under various conditions by Digital Image Correlation (DIC) system, its software also measured the complete strain history, in addition to sample curvature after it was impacted; therefore obtaining the data required to determine the amount of side-wall-curl (Aluminum sheet springback) present after formation.

Using the total probability theorem, we propose a method to calculate the failure rate of a linear vibratory system with random parameters excited by stationary Gaussian processes. The response of such a system is non-stationary because of the randomness of the input parameters. A space-filling design, such as optimal symmetric Latin hypercube sampling or maximin, is first used to sample the input parameter space. For each design point, the output process is stationary and Gaussian. We present two approaches to calculate the corresponding conditional probability of failure. A Kriging metamodel is then created between the input parameters and the output conditional probabilities allowing us to estimate the conditional probabilities for any set of input parameters. The total probability theorem is finally applied to calculate the time-dependent probability of failure and the failure rate of the dynamic system. The proposed method is demonstrated using a vibratory system.

Reliability-based design optimization accounts for variation. However, it assumes that statistical information is available in the form of fully defined probabilistic distributions. This is not true for a variety of engineering problems where uncertainty is usually given in terms of interval ranges. In this case, interval analysis or possibility theory can be used instead of probability theory. This paper shows how possibility theory can be used in design and presents a computationally efficient sequential optimization algorithm. The algorithm handles problems with only uncertain or a combination of random and uncertain design variables and parameters. It consists of a sequence of cycles composed of a deterministic design optimization followed by a set of worst-case reliability evaluation loops. A crank-slider mechanism example demonstrates the accuracy and efficiency of the proposed sequential algorithm.

Finite element analysis is a well-established methodology in structural dynamics. However, optimization and/or probabilistic studies can be prohibitively expensive because they require repeated FE analyses of large models. Various reanalysis methods have been proposed in order to calculate efficiently the dynamic response of a structure after a baseline design has been modified, without recalculating the new response. The parametric reduced-order modeling (PROM) and the combined approximation (CA) methods are two re-analysis methods, which can handle large model parameter changes in a relatively efficient manner. Although both methods are promising by themselves, they can not handle large FE models with large numbers of DOF (e.g. 100,000) with a large number of design parameters (e.g. 50), which are common in practice. In this paper, the advantages and disadvantages of the PROM and CA methods are first discussed in detail.

With increasing demand for cost-effective fuel cells, it is essential to investigate alternative materials for components of the fuel cells. The objective of this paper is to implement elastomeric materials (silicon rubber) for the use of bipolar plates in a polymer electrolyte membrane fuel cell. Two different types of conductive fillers, a graphite fiber and flake, were added at different concentrations to a two-component silicone rubber slurry. Electrical, thermal and mechanical properties of the composite material were investigated. Comparable electrical and thermal conductivities were achieved to that of commercially available plates. The silicone rubber based composite material maintained elastomeric properties for improved sealing of cell fluid reactants and products.

Finite element analysis is a standard tool for deterministic or probabilistic design optimization of dynamic systems. The optimization process requires repeated eigenvalue analyses which can be computationally expensive. Several reanalysis techniques have been proposed to reduce the computational cost including Parametric Reduced Order Modeling (PROM), Combined Approximations (CA), and the Modified Combined Approximations (MCA) method. Although the cost of reanalysis is substantially reduced, it can still be high for models with a large number of degrees of freedom and a large number of design variables. Reanalysis methods use a basis composed of eigenvectors from both the baseline and the modified designs which are in general linearly dependent. To eliminate the linear dependency and improve accuracy, Gram Schmidt orthonormalization is employed which is costly itself.

Self-piercing riveting (SPR) has been used in production to join sheet materials since the early 1990s. A large amount of experimental trial work was required in order to determine an appropriate combination of rivet and anvil design to fulfill the required joint parameters. The presented study is describing the methodology of SPR joint design based on numerical simulation and experimental methods of defining required simulation input parameters. The required inputs are the stress-strain curves of sheet materials and rivets for the range of strains taking place in the SPR joining process, parameters required for a fracture model for all involved materials, and friction parameters for all interfaces of SPR process. In the current study, the normalized Cockroft-Latham fracture criterion was used for predicting fracture. Custom hole and tube expansion tests were used for predicting fracture of the riveted materials and the rivet, respectively.

Two popular critical plane models developed by Fatemi-Socie and Smith-Watson-Topper were derived from the experimental observations of the nucleation and growth of cracks during loading. The Fatemi-Socie critical plane model is applicable for the life prediction of materials for which the dominant failure mechanism is shear crack nucleation and growth, while the Smith-Watson-Topper model, for materials that fail predominantly by crack growth on planes perpendicular to the planes of maximum tensile strain or stress. The two critical plane models have been validated primarily by in-phase and 90° out-of-phase loading, and few, on the complex, non-proportional loading paths. A successful critical plane model should be able to predict both the fatigue life and the dominant failure planes. However, some experimental studies indicate the 304 stainless steel has the two possible failure modes, shear and tensile failure dominant, depending on the loading mode and stress and strain states.

This paper explores the trade-off between reliability-based design and robustness for an automotive under-hood thermal system using the iSIGHT-FD environment. The interaction between the engine cooling system and the heating, ventilating, and air-conditioning (HVAC) system is described. The engine cooling system performance is modeled using Flowmaster and a metamodel is developed in iSIGHT. The actual HVAC system performance is characterized using test bench data. A design of experiment procedure determines the dominant factors and the statistics of the HVAC performance is obtained using Monte Carlo simulation (MCS). The MCS results are used to build an overall system response metamodel in order to reduce the computational effort. A multi-objective optimization in iSIGHT maximizes the system mean performance and simultaneously minimizes its standard deviation subject to probabilistic constraints.

The Press-in-Place (PIP) gasket is a static face seal with self-retaining feature, which is used for the mating surfaces of engine components to maintain the reliability of the closed system under various operating conditions. Its design allows it to provide enough contact pressure to seal the internal fluid as well as prevent mechanical failures. Insufficient sealing pressure will lead to fluid leakage, consequently resulting in engine failures. A test fixture was designed to simulate the clamp load and internal pressure condition on a gasket bolted joint. A sensor pad in combination with TEKSCAN equipment was used to capture the overall and local pressure distribution of the PIP gasket under various engine loading conditions. Then, the test results were compared with simulated results from computer models. Through the comparisons, it was found that gasket sealing pressure of test data and CAE data shows good correlations in all internal pressure cases when the bolt load was 500 N.