Search Results

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.

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.

A 2-D CFD model was developed to describe the heat transfer, and reaction kinetics in a honeycomb structured ceramic diesel particulate trap. This model describes the steady state as well as the transient behavior of the flow and heat transfer during the trap regeneration processes. The trap temperature profile was determined by numerically solving the 2-D unsteady energy equation including the convective, heat conduction and viscous dissipation terms. The convective terms were based on a 2-D analytical flow field solution derived from the conservation of mass and momentum equations (Opris, 1997). The reaction kinetics were described using a discretized first order Arrhenius function. The 2-D term describing the reaction kinetics and particulate matter conservation of mass was added to the energy equation as a source term in order to represent the particulate matter oxidation. The filtration model describes the particulate matter accumulation in the trap.

Time-averaged temperatures at critical locations on the piston of a direct-fuel injected, two-stroke, 388 cm3, research engine were measured using an infrared telemetry device. The piston temperatures were compared to data [7] of a carbureted version of the two-stroke engine, that was operated at comparable conditions. All temperatures were obtained at wide open throttle, and varying engine speeds (2000-4500 rpm, at 500 rpm intervals). The temperatures were measured in a configuration that allowed for axial heat flux to be determined through the piston. The heat flux was compared to carbureted data [8] obtained using measured piston temperatures as boundary conditions for a computer model, and solving for the heat flux. The direct-fuel-injected piston temperatures and heat fluxes were significantly higher than the carbureted piston. On the exhaust side of the piston, the direct-fuel injected piston temperatures ranged from 33-73 °C higher than the conventional carbureted piston.

An experimental and modeling study was conducted to study the passive regeneration of a catalyzed particulate filter (CPF) by the oxidation of particulate matter (PM) via thermal and Nitrogen dioxide/temperature-assisted means. Emissions data in the exhaust of a John Deere 6.8 liter, turbocharged and after-cooled engine with a low-pressure loop EGR and a diesel oxidation catalyst (DOC) - catalyzed particulate filter (CPF) in the exhaust system was measured and used for this study. A series of experiments was conducted to evaluate the performance of the DOC, CPF and DOC+CPF configurations at various engine speeds and loads.

A Cummins ISB 5.9 liter medium-duty engine with cooled EGR has been used to study an early extrusion of an advanced ceramic uncatalyzed diesel particulate filter (DPF). Data for the advanced ceramic material (ACM) and an uncatalyzed cordierite filter of similar dimensions are presented. Pressure drop data as a function of mass loadings (0, 4, and 6 grams of particulate matter (PM) per liter of filter volume) for various flow rate/temperature combinations (0.115 - 0.187 kg/sec and 240 - 375 °C) based upon loads of 15, 25, 40 and 60% of full engine load (684 N-m) at 2300 rpm are presented. The data obtained from these experiments were used to calibrate the MTU 1-D 2-Layer computer model developed previously at MTU. Clean wall permeability determined from the model calibration for the ACM was 5.0e-13 m2 as compared to 3.0e-13 m2 for cordierite.

A quasi-steady 1-dimensional computer model of a catalyzed particulate filter (CPF) capable of simulating active regeneration of the CPF via diesel fuel injection upstream of a diesel oxidation catalyst (DOC) or other means to increase the exhaust gas temperature has been developed. This model is capable of predicting gaseous species concentrations (HC's, CO, NO and NO2) and exhaust gas temperatures within and after the CPF, for given input values of gaseous species and PM concentrations before the CPF and other inlet variables such as time-varying temperature of the exhaust gas at the inlet of the CPF and volumetric flow rate of exhaust gas.

The discrete droplet model is widely used to describe two-phase flows such as high-velocity dense sprays. The interaction between the liquid and the gas phase is modeled via appropriate source terms in the gas phase equations. This approach can lead to a strong dependence of the liquid-gas coupling on the spatial resolution of the gas phase. The liquid-gas coupling requires the computation of source terms using the gas phase properties, and, subsequently, these sources are then distributed onto the gas phase mesh. In this study, a Lagrange polynomial interpolation method has been developed to evaluate the source terms and also to distribute these source terms onto the gas mesh. The focus of this investigation has been on the momentum exchange between the two phases. The Lagrange polynomial interpolation and source term distribution methods are evaluated for non-evaporating sprays using KIVA3 as a modeling platform.

This paper presents the damping effectiveness of free-layer damping materials through standard Oberst bar testing, solid plate excitation (RTC3) testing, and prediction through numerical schemes. The main objective is to compare damping results from various industry test methods to performance in an automotive body structure. Existing literature on laboratory and vehicle testing of free-layer viscoelastic damping materials has received significant attention in recent history. This has created considerable confusion regarding the appropriateness of different test methods to measure material properties for damping materials/treatments used in vehicles. The ability to use the material properties calculated in these tests in vehicle CAE models has not been extensively examined. Existing literature regarding theory and testing for different industry standard damping measurement techniques is discussed.

The induction hardening process involves a complex interaction of electromagnetic heating, rapid cooling, metallurgical phase transformations, and mechanical behavior. Many factors including induction coil design, power, frequency, scanning velocity, workpiece geometry, material chemistry, and quench severity determine a process outcome. This paper demonstrates an effective application of a numerical analysis tool for understanding of induction hardening. First, an overview of the Caterpillar induction simulation tool is briefly discussed. Then, several important features of the model development are examined. Finally, two examples illustrating the use of the computer simulation tool for solving induction-hardening problems related to cracking and distortion are presented. These examples demonstrate the tool's ability to simulate changes in process parameters and latitude of modeling steel or cast iron.

A one-dimensional, two layer computational model was developed to predict the behavior of a clean and particulate-loaded catalyzed wall-flow diesel particulate filter (CPF). The model included the mechanisms of particle deposition inside the CPF porous wall and on the CPF wall surface, the exhaust flow field and temperature field inside the CPF, as well as the particulate catalytic oxidation mechanisms accounting for the catalyst-assisted particulate oxidation by the catalytic coating in addition to the conventional particulate thermal oxidation. The paper also develops the methodology for calibrating and validating the model with experimental data. Steady state loading experiments were performed to calibrate and validate the model.

Non-evaporating diesel sprays have been simulated utilizing the ETAB and the WAVE atomization and breakup models and have been compared with experimental data. The experimental penetrations and widths were determined from back-lit spray images and the droplet sizes have been measured by means of a Malvern particle sizer. The model evaluation criteria include the spray penetration, the spray width and the local droplet size. The comparisons have been performed for variations of the injection pressure, the gas density and the fuel viscosity. The fuel nozzle exit velocities used in the simulations have been computed with a special code that considers the effect of in-nozzle cavitation. The simulations showed good overall agreement with experimental data. However, the capabilities of the models to predict the droplet size for different fuels could be improved.

An analysis of vehicle engine cooling airflow by means of a one-dimensional, transient, compressible flow model was carried out and revealed that similarity theory could be applied to investigate the variation of the airflow with ambient and operating conditions. It was recognized that for a given vehicle engine cooling system, the cooling airflow behavior could be explained using several dimensionless parameters that involve the vehicle speed, fan speed, heat transfer rate through the radiator, ambient temperature and pressure, and the system characteristic dimension. Using the flow resistance and fan characteristics measured from a prototype cooling system and the computer simulation for the one-dimensional compressible flow model, a quantitative correlation of non-dimensional mass flow rate to three dimensionless parameters for a prototype heavy-duty truck was established. The results are presented in charts, tables, and formulas.

Objective of this work is the incorporation of the flame stretch effects in an Eulerian-Lagrangian model for premixed SI combustion in order to describe ignition and flame propagation under highly inhomogeneous flow conditions. To this end, effects of energy transfer from electrical circuit and turbulent flame propagation were fully decoupled. The first ones are taken into account by Lagrangian particles whose main purpose is to generate an initial burned field in the computational domain. Turbulent flame development is instead considered only in the Eulerian gas phase for a better description of the local flow effects. To improve the model predictive capabilities, flame stretch effects were introduced in the turbulent combustion model by using formulations coming from the asymptotic theory and recently verified by means of DNS studies. Experiments carried out at Michigan Tech University in a pressurized, constant-volume vessel were used to validate the proposed approach.

The performance of five positive k-factor injector tips has been assessed in this work by analyzing a comprehensive set of injected mass, momentum, and spray measurements. Using high speed shadowgraphs of the injected diesel plumes, the sensitivities of measured vapor penetration and dispersion to injection pressure (100-250MPa) and ambient density (20-52 kg/m3) have been compared with the Naber-Siebers empirical spray model to gain understanding of second order effects of orifice diameter. Varying in size from 137 to 353μm, the orifice diameters and corresponding injector tips are appropriate for a relatively wide range of engine cylinder sizes (from 0.5 to 5L). In this regime, decreasing the orifice exit diameter was found to reduce spray penetration sensitivity to differential injection pressure. The cone angle and k-factored orifice exit diameter were found to be uncorrelated.

This paper aims to extend the existing Volume of Fluid (VOF) model by implementing an evaporation sub-model in an open source Computational Fluid Dynamics (CFD) software, OpenFOAM. The paper applies the new model to numerically study the evaporation of spherical n-heptane droplets impinging on a hot wall at atmospheric pressure and a temperature above the Leidenfrost temperature. Volume of Fluid (VOF) method is chosen to track the liquid gas interface and the capability of VOF method implemented in interDyMFoam solver of OpenFOAM to simulate hydrodynamics during droplet-droplet interaction and droplet-film interaction is explored. Firstly, the in-built solver is used to simulate problems in isothermal conditions and the simulation results are compared qualitatively with the published results to validate the solver. A numerical method for modeling heat and mass transfer during evaporation is implemented in conjunction with the VOF.

Diesel combustion and emissions formation is largely spray and mixing controlled and hence understanding spray parameters, specifically vaporization, is key to determine the impact of fuel injector operation and nozzle design on combustion and emissions. In this study, an eight-hole common rail piezoelectric injector was tested in an optically accessible constant volume combustion vessel at charge gas conditions typical of full load boosted engine operation. Liquid penetration of the eight sprays was determined via processing of images acquired from Mie back scattering under vaporizing conditions by injecting into a charge gas at elevated temperature with 0% oxygen. Conditions investigated included a charge temperature sweep of 800 to 1300 K and injection pressure sweep of 1034 to 2000 bar at a constant charge density of 34.8 kg/m₃.

The paper describes an experimental activity on the spatial and temporal liquid- and vapor-phase distributions of diesel fuel at engine-like conditions. The influence of the k-factor (0 and 1.5) of a single-hole axial-disposed injector (0.100 mm diameter and 10 L/d ratio) has been studied by spraying fuel in an optically-accessible constant-volume combustion vessel. A high-speed imaging system, capable of acquiring Mie-scattering and Schlieren images in a near simultaneous fashion mode along the same line of sight, has been developed at the Michigan Technological University using a high-speed camera and a pulsed-wave LED system. The time resolved pair of schlieren and Mie-scattering images identifies the instantaneous position of both the vapor and liquid phases of the fuel spray, respectively. The studies have been performed at three injection pressures (70, 120 and 180 MPa), 23.9 kg/m3 ambient gas density and 900 K gas temperature in the vessel.

A 2-D computational model was developed to describe the flow and filtration processes, in a honeycomb structured ceramic diesel particulate trap. This model describes the steady state trap loading, as well as the transient behavior of the flow and filtration processes. The theoretical model includes the effect of a copper fuel additive on trap loading and transient operation. The convective terms were based on a 2-D analytical flow field solution derived from the conservation of mass and momentum equations. The filtration theory incorporated in the time dependent numerical code included the diffusion, inertia, and direct interception mechanisms. Based on a measured upstream particle size distribution, using the filtration theory, the downstream particle size distribution was calculated. The theoretical filtration efficiency, based on particle size distribution, agreed very well (within 1%) with experimental data for a number of different cases.

Caterpillar uses heat treatment to enhance the properties of a significant number of parts. Traditional heat treat process optimization is both time consuming and expensive when done by empirical methods. This paper describes a computer simulation of the heat treatment process, developed by Caterpillar, based upon finite element analysis. This approach combines thermal, microstructural, and stress analysis to accurately model material transformation during quenching. Examples are presented to illustrate the program.