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.

Partial enclosures are commonly utilized to reduce the radiated noise from equipment. Often, enclosure openings are fitted with silencers or louvers to further reduce the noise emitted. In the past, the boundary element method (BEM) has been applied to predict the insertion loss of the airborne path with good agreement with measurement. However, an alteration at the opening requires a new model and additional computational time. In this paper, a transfer function method is proposed to reduce the time required to assess the effect of modifications to an enclosure. The proposed method requires that the impedance at openings be known. Additionally, transfer functions relating the sound pressure at one opening to the volume velocity at other openings must be measured or determined using simulation. It is assumed that openings are much smaller than an acoustic wavelength. The sound power from each opening is determined from the specific acoustic impedance and sound pressure at the opening.

It was demonstrated that a bypass duct similar to a Herschel- Quincke tube could be used to increase the transmission loss of mufflers at selected frequencies. In many cases, the duct can be short and thought of as a leak. It was shown that the optimal length and cross-sectional area could be determined by using a simple optimization technique known as the Vincent Circle. Most importantly, it was demonstrated that the attenuation at low frequencies could be improved by as much as 15 dB. To prove the concept, a muffler was designed and optimized using transfer matrix theory. Then, the optimized muffler was constructed and the transmission loss measured using the two-load method. The measured results compared well with prediction from transfer matrix theory. Boundary element simulation was then used to further study the attenuation mechanism.

High pressure diesel sprays were visualized under vaporizing and combusting conditions in a constant-volume combustion vessel. Near-simultaneous visualization of vapor and liquid phase fuel distribution were acquired using a hybrid shadowgraph/Mie-scattering imaging setup. This imaging technique used two pulsed LED's operating in an alternative manner to provide proper light sources for both shadowgraph and Mie scattering. In addition, combustion cases under the same ambient conditions were visualized through high-speed combustion luminosity measurement. Two single-hole diesel injectors with same nozzle diameters (100μm) but different k-factors (k0 and k1.5) were tested in this study. Detailed analysis based on spray penetration rate curves, rate of injection measurements, combustion indicators and 1D model comparison have been performed.

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).

Transmission loss (TL) is a good performance measure of mufflers since it represents the muffler's inherent capability of sound attenuation. There are several existing numerical methods, which have been widely used to calculate the TL from numerical simulation results, such as the four-pole and three-point methods. In this paper, a new approach is proposed to evaluate the transmission loss based on the reciprocal work identity. The proposed method does not assume plane wave propagation in the inlet and outlet ducts, and more importantly, does not explicitly apply the anechoic termination impedance at the outlet. As a result, it has the potential of extending TL computation above the plane wave cut-off frequency.

Intake, exhaust, and heating / air conditioning systems in automobiles consist of various common duct elements. Noise arises primarily due to the source and is attenuated using common elements like expansion chambers and resonators. This attenuation is straightforward to predict using plane wave simulation and more advanced numerical methods. However, flow noise is often an unexpected important noise source. Predictions require computer intensive analyses. To better understand the aeroacoustic sources in duct systems, a flow rig has been developed at the University of Kentucky. The flow rig consists of a blower, a silencer to attenuate blower noise, external noise sources, and then the test duct. The flow rig can be equipped with an anechoic termination to measure transmission loss or may be used to measure insertion loss directly. In the latter case, the sound power is measured from the pipe outlet inside of a hemi-anechoic chamber.

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.

Transmissibility is the most common metric used for isolator characterization. However, engineers are becoming increasingly concerned about energy transmission through an isolator at high frequencies and how the compliance of the machine and foundation factor into the performance. In this paper, the transfer matrix approach for isolator characterization is first reviewed. Two methods are detailed for determining the transfer matrix of an isolator using finite element simulation. This is accomplished by determining either the mobility or impedance matrix for the isolator and then converting to a transfer matrix. It is shown that results are similar using either approach. In both cases, the isolator is first pre-loaded before the transfer matrix is determined. The approach to find isolator insertion loss is demonstrated for an isolator between two plates, and the effect of making changes to the structural impedance on the machine side of the isolator by adding ribs is examined.

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.

Despite numerous research efforts, there is no reliable and widely accepted tool for the prediction of erosion prone material surfaces due to collapse of cavitation bubbles. In the present paper an Erosion Aggressiveness Index (EAI) is proposed, based on the pressure loads which develop on the material surface and the material yield stress. EAI depends on parameters of the liquid quality and includes the fourth power of the maximum bubble radius and the bubble size number density distribution. Both the newly proposed EAI and the Cavitation Aggressiveness Index (CAI), which has been previously proposed by the authors based on the total derivative of pressure at locations of bubble collapse (DP/Dt>0, Dα/Dt<0), are computed for a cavitating flow orifice, for which experimental and numerical results on material erosion have been published. The predicted surface area prone to cavitation damage, as shown by the CAI and EAI indexes, is correlated with the experiments.

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.

Microperforated panel absorbers are best considered as the combination of the perforate and the backing cavity. They are sometimes likened to Helmholtz resonators. This analogy is true in the sense that they are most effective at the resonant frequencies of the panel-cavity combination when the particle velocity is high in the perforations. However, unlike traditional Helmholtz resonators, microperforated absorbers are broader band and the attenuation mechanism is dissipative rather than reactive. It is well known that the cavity depth governs the frequency bands of high absorption. The work presented here focuses on the development, modeling and testing of novel configurations of backing constructions and materials. These configurations are aimed at both dialing in the absorption properties at specific frequencies of interest and creating broadband sound absorbers. In this work, several backing cavity strategies are considered and evaluated.

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.

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.

Autonomous vehicle operation is dependent upon accurate position estimation and thus a major concern of implementing the autonomous navigation is obtaining robust and accurate data from sensors. This is especially true, in case of Inertial Measurement Unit (IMU) sensor data. The IMU consists of a 3-axis gyro, 3-axis accelerometer, and 3-axis magnetometer. The IMU provides vehicle orientation in 3D space in terms of yaw, roll and pitch. Out of which, yaw is a major parameter to control the ground vehicle’s lateral position during navigation. The accelerometer is responsible for attitude (roll-pitch) estimates and magnetometer is responsible for yaw estimates. However, the magnetometer is prone to environmental magnetic disturbances which induce errors in the measurement.

In-cylinder ion sensing is a subject of interest due to its application in spark-ignited (SI) engines for feedback control and diagnostics including: combustion knock detection, rate and phasing of combustion, and mis-fire On Board Diagnostics (OBD). Further advancement and application is likely to continue as the result of the availability of ignition coils with integrated ion sensing circuitry making ion sensing more versatile and cost effective. In SI engines, combustion knock is controlled through closed loop feedback from sensor metrics to maintain knock near the borderline, below engine damage and NVH thresholds. Combustion knock is one of the critical applications for ion sensing in SI engines and improvement in knock detection offers the potential for increased thermal efficiency. This work analyzes and characterizes the ionization signal in reference to the cylinder pressure signal under knocking and non-knocking conditions.

In this paper, a model-based linear estimator and a non-linear control law for an Fe-zeolite urea-selective catalytic reduction (SCR) catalyst for heavy duty diesel engine applications is presented. The novel aspect of this work is that the relevant species, NO, NO2 and NH3 are estimated and controlled independently. The ability to target NH3 slip is important not only to minimize urea consumption, but also to reduce this unregulated emission. Being able to discriminate between NO and NO2 is important for two reasons. First, recent Fe-zeolite catalyst studies suggest that NOx reduction is highly favored by the NO 2 based reactions. Second, NO2 is more toxic than NO to both the environment and human health. The estimator and control law are based on a 4-state model of the urea-SCR plant. A linearized version of the model is used for state estimation while the full nonlinear model is used for control design.

Virtual Product Development (VPD) is a key enabler in CAE and depends upon accurate implementation of multibody dynamics. This paper discusses the formulation and implementation of a large-scale multibody dynamics simulation code. In the presented formulation, the joint variables are used as the generalized coordinates and spatial algebra is used to formulate the system equations of motion. Although the presented formulation utilizes the joint variables as the generalized coordinates, closed-loop mechanisms can be easily modeled using impeded constraints. Baumgart stabilization approach is used to eliminate the constraint violations without using the expensive Newton-Raphson iterations. Integrated rigid and flexible body dynamic simulation allows accurate prediction of structural loads, stress, and strains. Both modal and nodal flexible body approaches are discussed in the paper.