Search Results

Premixed Charge Compression Ignition (PCCI), a Low Temperature Combustion (LTC) strategy for diesel engines is of increasing interest due to its potential to simultaneously reduce soot and NOx emissions. However, the influence of mixture preparation on combustion phasing and heat release rate in LTC is not fully understood. In the present study, the influence of injection timing on mixture preparation, combustion and emissions in PCCI mode is investigated by experimental and computational methods. A sequential coupling approach of 3D CFD with a Stochastic Reactor Model (SRM) is used to simulate the PCCI engine. The SRM accounts for detailed chemical kinetics, convective heat transfer and turbulent micro-mixing. In this integrated approach, the temperature-equivalence ratio statistics obtained using KIVA 3V are mapped onto the stochastic particle ensemble used in the SRM.

Since transient vehicle HVAC computational fluids (CFD) simulations take too long to solve in a production environment, the goal of this project is to automatically create a lumped-parameter flow network from a steady-state CFD that solves nearly instantaneously. The data mining algorithm k-means is implemented to automatically discover flow features and form the network (a reduced order model). The lumped-parameter network is implemented in the commercial thermal solver MuSES to then run as a fully transient simulation. Using this network a “localized heat transfer coefficient” is shown to be an improvement over existing techniques. Also, it was found that the use of the clustering created a new flow visualization technique. Finally, fixing clusters near equipment newly demonstrates a capability to track localized temperatures near specific objects (such as equipment in vehicles).

In the development of an FAW SUV, one of the goals is to achieve a state of the art drag level. In order to achieve such an aggressive target, feedback from aerodynamics has to be included in the early stage of the design decision process. The aerodynamic performance evaluation and improvement is mostly based on CFD simulation in combination with some wind tunnel testing for verification of the simulation results. As a first step in this process, a fully detailed simulation model is built. The styling surface is combined with engine room and underbody detailed geometry from a similar size existing vehicle. From a detailed analysis of the flow field potential areas for improvement are identified and five design parameters for modifying overall shape features of the upper body are derived. In a second step, a response surface method involving design of experiments and adaptive sampling techniques are applied for characterizing the effects of the design changes.

Coatings have the potential to improve bearing tribological performance. However, every coating application process and material combination may create different residual stresses and coating microstructures, and their effect on bearing fatigue and wear performance is unclear. The aim of this work is to investigate coating induced residual stress effects on bearing failure indicators using a microstructural contact mechanics (MSCM) finite element (FE) model. The MSCM FE model consists of a two-dimensional FE model of a coated bearing surface under sliding contact where individual grains are represented by FE domains. Interactions between FE domains are represented using contact element pairs. Unique to this layered rolling contact FE model is the use of polycrystalline material models to represent realistic bearing and coating microstructural behavior. The MSCM FE model was compared to a second non-microstructural contact mechanics (non-MSCM) model.

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.

A twist beam rear suspension system is modeled, analyzed and optimized in this paper. An ADAMS model is established based on the REC (Rigid-Elastic Coupling) Theory, which is verified by FEM (Finite Element Method) approach, the effects of the geometric parameters on the twist beam suspension performance are investigated. In order to increase the calculation efficiency and improve the simulation accuracy, a neural network model and NSGA II (Non-dominated Sorting Genetic Algorithm II) are adopted to conduct a multi-objective optimization on a twist beam rear suspension system.

This paper presents a combined numerical and experimental investigation of the characteristics of spark discharge in a spark-ignition engine. The main objective of this work is to gain insights into the spark discharge process and early flame kernel development. Experiments were conducted in an inert medium within an optically accessible constant-volume combustion vessel. The cross-flow motion in the vessel was generated using a previously developed shrouded fan. Numerical modeling was based on an existing discharge model in the literature developed by Kim and Anderson. However, this model is applicable to a limited range of gas pressures and flow fields. Therefore, the original model was evaluated and improved to predict the behavior of spark discharge at pressurized conditions up to 45 bar and high-speed cross-flows up to 32 m/s. To accomplish this goal, a parametric study on the spark channel resistance was conducted.

Structure-borne inputs to hybrid FEA/SEA models could have significant effects on the model prediction accuracy. The purpose of this work was to obtain the structure-borne noise (SBN) inputs using a simplified transfer path analysis (TPA) and identify the significance of the structure-borne and airborne contributions to the spectator sound power of an engine with enclosure for future modeling references. Force inputs to the enclosure from the engine were obtained and used as inputs to a hybrid engine enclosure model for sound prediction.

Global surrogate modeling aims to build surrogate model with high accuracy in the whole design domain. A major challenge to achieve this objective is how to reduce the number of function evaluations to the original computer simulation model. To date, the most widely used approach for global surrogate modeling is the adaptive surrogate modeling method. It starts with an initial surrogate model, which is then refined adaptively using the mean square error (MSE) or maximizing the minimum distance criteria. It is observed that current methods may not be able to effectively construct a global surrogate model when the underlying black box function is highly nonlinear in only certain regions. A new surrogate modeling method which can allocate more training points in regions with high nonlinearity is needed to overcome this challenge. This article proposes an efficient global surrogate modeling method based on a multi-layer sampling scheme.

Optimizing the trade-off between reliability and cost of operating a microgrid, including vehicles as both loads and sources, can be a challenge. Optimal energy management is crucial to develop strategies to improve the efficiency and reliability of microgrids, as well as new communication networks to support optimal and reliable operation. Prior approaches modeled the grid using MATLAB, but did not include the detailed physics of loads and sources, and therefore missed the transient effects that are present in real-time operation of a microgrid. This article discusses the implementation of a physics-based detailed microgrid model including a diesel generator, wind turbine, photovoltaic array, and utility. All elements are modeled as sources in Simulink. Various loads are also implemented including an asynchronous motor. We show how a central control algorithm optimizes the microgrid by trying to maximize reliability while reducing operational cost.

Computer-aided engineering (CAE) is an important tool routinely used to simulate complex engineering systems. Virtual simulations enhance engineering insight into prospective designs and potential design issues and can limit the need for expensive engineering prototypes. For complex engineering systems, however, the effectiveness of virtual simulations is often hindered by excessive computational cost. To minimize the cost of running expensive computer simulations, approximate models of the original model (often called surrogate models or metamodels) can provide sufficient accuracy at a lower computing overhead compared to repeated runs of a full simulation. Metamodel accuracy improves if constructed using space-filling designs of experiments (DOEs). The latter provide a collection of sample points in the design space preferably covering the entire space.

Two models are presented for the long life (107 cycles) axial fatigue strength of four ferrous powder metal (PM) material series: sintered and heat-treated iron-carbon steel, iron-copper and copper steel, iron-nickel and nickel steel, and pre-alloyed steel. The materials are defined at ranges of carbon content and densities using the broad data available in the Metal Powder Industries Federation (MPIF) Standard 35 for PM structural parts. The first model evaluates 107 cycles axial fatigue strength as a function of ultimate strength and the second model as a function of hardness. For all 118 studied materials, both models are found to have a good correlation between calculated and 107 cycles axial fatigue strength with a high Pearson correlation coefficient of 0.97. The article provides details on the model development and the reasoning for selecting the ultimate strength and hardness as the best predictors for 107 cycles axial fatigue strength.

In this paper, residual stress distributions in rectangular bars due to rolling or burnishing at very high rolling or burnishing loads are investigated by roll burnishing experiments and three-dimensional finite element analyses using ABAQUS. First, roll burnishing experiments on rectangular bars at two roller burnishing loads are presented. The results indicate the higher burnishing load induces lower residual stresses and the higher burnishing load does not improve fatigue lives. Next, in the corresponding finite element analyses, the roller is modeled as rigid and the roller rolls on the flat surface of the bar with a low coefficient of friction. The bar material is modeled as an elastic-plastic strain hardening material with a nonlinear kinematic hardening rule for loading and unloading.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

Weight reduction is very important in automotive design because of stringent demand on fuel economy. Structural optimization of dynamic systems using finite element (FE) analysis plays an important role in reducing weight while simultaneously delivering a product that meets all functional requirements for durability, crash and NVH. With advancing computer technology, the demand for solving large FE models has grown. Optimization is however costly due to repeated full-order analyses. Reanalysis methods can be used in structural vibrations to reduce the analysis cost from repeated eigenvalue analyses for both deterministic and probabilistic problems. Several reanalysis techniques have been introduced over the years including Parametric Reduced Order Modeling (PROM), Combined Approximations (CA) and the Epsilon algorithm, among others.

To improve fuel economy, there is a trend in automotive industry to use light weight, high strength materials. Automotive body structures are composed of several panels which must be downsized to reduce weight. Because this affects NVH (Noise, Vibration and Harshness) performance, engineers are challenged to recover the lost panel stiffness from down-gaging in order to improve the structure borne noise transmitted through the lightweight panels in the frequency range of 100-300 Hz where most of the booming and low medium frequency noise occurs. The loss in performance can be recovered by optimized panel geometry using beading or damping treatment. Topography optimization is a special class of shape optimization for changing sheet metal shapes by introducing beads. A large number of design variables can be handled and the process is easy to setup in commercial codes. However, optimization methods are computationally intensive because of repeated full-order analyses.

Fundamental understanding of the sources of fuel-derived Unburned Hydrocarbon (UHC) emissions in heavy duty diesel engines is a key piece of knowledge that impacts engine combustion system development. Current emissions regulations for hydrocarbons can be difficult to meet in-cylinder and thus after treatment technologies such as oxidation catalysts are typically used, which can be costly. In this work, Computational Fluid Dynamics (CFD) simulations are combined with engine experiments in an effort to build an understanding of hydrocarbon sources. In the experiments, the combustion system design was varied through injector style, injector rate shape, combustion chamber geometry, and calibration, to study the impact on UHC emissions from mixing-controlled diesel combustion.

This paper investigates an optimal design of a diesel engine and after-treatment systems for a series hybrid electric forklift application. A holistic modeling approach is developed in GT-Suite® to establish a model-based hardware definition for a diesel engine and an after-treatment system to accurately predict engine performance and emissions. The used engine model is validated with the experimental data. The engine design parameters including compression ratio, boost level, air-fuel ratio (AFR), injection timing, and injection pressure are optimized at a single operating point for the series hybrid electric vehicle, together with the performance of the after-treatment components. The engine and after-treatment models are then coupled with a series hybrid electric powertrain to evaluate the performance of the forklift in the standard VDI 2198 drive cycle.

In order to extend the operability limit of the gasoline compression ignition (GCI) engine, as an avenue for low temperature combustion (LTC) regime, the effects of parametric variations of engine operating conditions on the performance of six-stroke GCI (6S-GCI) engine cycle are numerically investigated, using an in-house 3D CFD code coupled with high-fidelity physical sub-models along with the Chemkin library. The combustion and emissions were calculated using a skeletal chemical kinetics mechanism for a 14-component gasoline surrogate fuel. Authors’ previous study highlighted the effects of the variation of injection timing and split ratio on the overall performance of 6S-GCI engine and the unique mixing-controlled burning mode of the charge mixtures during the two additional strokes. As a continuing effort, the present study details the parametric studies of initial gas temperature, boost pressure, fuel injection pressure, compression ratio, and EGR ratio.

Increasing regulatory demand for efficiency has led to development of novel combustion modes such as HCCI, GCI and RCCI for gasoline light duty engines. In order to realize HCCI as a compression ignition combustion mode system, in-cylinder compression temperatures must be elevated to reach the autoignition point of the premixed fuel/air mixture. This should be co-optimized with appropriate fuel formulations that can autoignite at such temperatures. CFD combustion modeling is used to model the auto ignition of gasoline fuel under compression ignition conditions. Using the fully detailed fuel mechanism consisting of thousands of components in the CFD simulations is computationally expensive. To overcome this challenge, the real fuel is represented by few major components of create a surrogate fuel mechanism. In this study, 9 variations of gasoline fuel sets were chosen as candidates to run in HCCI combustion mode.