Search Results

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

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.

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.

Improving vehicle response through advanced knowledge of traffic behavior can lead to large improvements in energy consumption for the single isolated vehicle. This energy savings across multiple vehicles can even be larger if they travel together as a cohort in harmonization. Additionally, if the vehicles have enough information about their immediate path of travel, and other vehicles’ in that path (and their respective critical forward-looking information), they can safely drive close enough to each other to share aerodynamic load. These energy savings can be upwards of multiple percentage points, and are dependent on several criteria. This analysis looks at criteria that contributes to energy savings for a cohort of vehicles in synchronous motion, as well as describes a study that allows for better understanding of the potential benefits of different types of cohorted vehicles in different platoon arrangements.

This paper investigates flow and combustion in a full-cycle simulation of a four-stroke, three-valve HCCI engine by visualizing the flow with pathlines. Pathlines trace massless particles in a transient flow field. In addition to visualization, pathlines are used here to trace the history, or evolution, of flow fields and species. In this study evolution is followed from the intake port through combustion. Pathline analysis follows packets of intake charge in time and space from induction through combustion. The local scalar fields traversed by the individual packets in terms of velocity magnitude, turbulence, species concentration and temperatures are extracted from the simulation results. The results show how the intake event establishes local chemical and thermal environments in-cylinder and how the species respond (chemically react) to the local field.

The purpose of this research is to develop a prediction technique that can be used in the development of port fuel-injection (hereinafter called PFI) gasoline engines, especially for general purpose small utility engines. Utility engines have two contradictory desirable aspects: compactness and high-power at wide open throttle. Therefore, applying the port fuel injector to utility engines presents a unique intractableness that is different from application to automobiles or motorcycles. At the condition of wide open throttle, a large amount of fuel is required to output high power, and injected fuel is deposited as a wall film on the intake port wall. Despite the fuel rich condition, emissions are required to be kept under a certain level. Thus, it is significant to understand the wall film phenomenon and control film thickness in the intake ports.

The ability to make fully resolved turbulent scalar field measurements has been demonstrated in an internal combustion engine using one-dimensional fluorobenzene fluorescence measurements. Data were acquired during the intake stroke in a motored engine that had been modified such that each intake valve was fed independently, and one of the two intake streams was seeded with the fluorescent tracer. The scalar energy spectra displayed a significant inertial subrange that had a −5/3 wavenumber power dependence. The scalar dissipation spectra were found to extend in the high-wavenumber regime, to where the magnitude was more than two decades below the peak value, which indicates that for all practical purposes the measurements faithfully represent all of the scalar dissipation in the flow.

Combustion systems with advanced injection strategies have been extensively studied, but there still exists a significant fundamental knowledge gap on fuel spray interactions with the piston surface and chamber walls. This paper is meant to provide detailed data on spray-wall impingement physics and support the spray-wall model development. The experimental work of spray-wall impingement with non-vaporizing spray characterization, was carried out in a high pressure-temperature constant-volume combustion vessel. The simultaneous Mie scattering of liquid spray and schlieren of liquid and vapor spray were carried out. Diesel fuel was injected at a pressure of 1500 bar into ambient gas at a density of 22.8 kg/m3 with isothermal conditions (fuel, ambient, and plate temperatures of 423 K). A Lagrangian-Eulerian modeling approach was employed to characterize the spray-gas and spray-wall interactions in the CONVERGETM framework by means of a Reynolds-Averaged Navier-Stokes (RANS) formulation.

Fundamental simulations using DNS type procedures were used to investigate the ignition, combustion characteristics and the lift-off trends of a spatially evolving turbulent liquid fuel jet. In particular, the spatially evolving n-Heptane spray injected in a two-dimensional rectangular domain with an engine like environment was investigated. The computational results were compared to the experimental observations from an optical engine as reported in the literature. It was found that an initial fuel rich combustion downstream of the spray tip is followed by diffusion combustion. Investigations were also made to understand the effects of injection velocity, ambient temperature and the droplet radius on the lift-off length. For each of these parameters three different values in a given range were chosen. For both injection velocity and droplet radius, an increase resulted in a near linear increase in the lift-off length.

The role of mixture stratification on combustion rate has been investigated in a constant volume combustion vessel in which mixtures of different equivalence ratios can be added in a spatially and temporally controlled fashion. The experiments were performed in a regime of low fluid motion to avoid the complicating effects of turbulence generated by the injection of different masses of fluid. Different mixture combinations were investigated while maintaining a constant overall equivalence ratio and initial pressure. The results indicate that the highest combustion rate for an overall lean mixture is obtained when all of the fuel is contained in a stoichiometric mixture in the vicinity of the ignition source. This is the result of the high burning velocity of these mixtures, and the complete oxidation which releases the full chemical energy.

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.

A global optimization method has been developed for an engine simulation code and utilized in the search of optimal fuel injection strategies. This method uses a Lagrange interpolation function which interpolates engine output data generated at the vertices and the intermediate points of the input parameters. This interpolation function is then used to find a global minimum over the entire parameter set, which in turn becomes the starting point of a CFD-based optimization. The CFD optimization is based on a steepest descent method with an adaptive cost function, where the line searches are performed with a fast-converging backtracking algorithm. The adaptive cost function is based on the penalty method, where the penalty coefficient is increased after every line search. The parameter space is normalized and, thus, the optimization occurs over the unit cube in higher-dimensional space.

The lift-off length plays a significant role in spray combustion as it influences the air entrainment upstream of the lift-off location and hence the soot formation. Accurate prediction of lift-off length thus becomes a prerequisite for accurate soot prediction in lifted flames. In the present study, KIVA-3v coupled with CHEMKIN, as developed at the Engine Research Center (ERC), is used as the CFD model. Experimental data from the Sandia National Labs. is used for validating the model predictions of n-heptane lift-off lengths and soot formation details in a constant volume combustion chamber. It is seen that the model predictions, in terms of lift-off length and soot mass, agree well with the experimental results for low ambient density (14.8 kg/m3) cases with different EGR rates (21% O2 - 8% O2). However, for high density cases (30 kg/m3) with different EGR rates (15% O2 - 8% O2) disagreements were found.

Multi-zone CFD simulations with detailed kinetics were used to model iso-octane HCCI experiments performed on a single-cylinder research engine. The modeling goals were to validate the method (multi-zone combustion modeling) and the reaction mechanism (LLNL 857 species iso-octane) by comparing model results to detailed exhaust speciation data, which was obtained with gas chromatography. The model is compared to experiments run at 1200 RPM and 1.35 bar boost pressure over an equivalence ratio range from 0.08 to 0.28. Fuel was introduced far upstream to ensure fuel and air homogeneity prior to entering the 13.8:1 compression ratio, shallow-bowl combustion chamber of this 4-stroke engine. The CFD grid incorporated a very detailed representation of the crevices, including the top-land ring crevice and head-gasket crevice. The ring crevice is resolved all the way into the ring pocket volume. The detailed grid was required to capture regions where emission species are formed and retained.

Transient engine operations are modeled and simulated with a 1-D code (GT Power) using heat release and emission data computed by a 3-D CFD code (Kiva3). During each iteration step of a transient engine simulation, the 1-D code utilizes the 3-D data to interpolate the values for heat release and emissions. The 3-D CFD computations were performed for the compression and combustion stroke of strategically chosen engine operating points considering engine speed, torque and excess air. The 3-D inlet conditions were obtained from the 1-D code, which utilized 3-D heat release data from the previous 1-D unsteady computations. In most cases, only two different sets of 3-D input data are needed to interpolate the transient phase between two engine operating points. This keeps the computation time at a reasonable level. The results are demonstrated on the load response of a generator which is driven by a medium-speed diesel engine.

Fundamental simulations in a quiescent cell under adiabatic conditions were made to understand the effect of temperature, equivalence ratio and the components of the recirculated exhaust gas, viz., CO2 and H2O, on the combustion of n-Heptane. Simulations were made in single phase in which evaporated n-Heptane was uniformly distributed in the domain. Computations were made for two different temperatures and four different EGR levels. CO2 or H2O or N2was used as EGR. It was found that the initiation of the main combustion process was primarily determined by two competing factors, i.e., the amount of initial OH concentration in the domain and the specific heat of the mixture. Further, initial OH concentration can be controlled by the manipulating the ambient temperature in the domain, and the specific heat capacity of the mixture via the mixture composition. In addition to these, the pre combustion and the subsequent post combustion can also be controlled via the equivalence ratio.

The development of new high power diesel engines is continually going for increased mean effective pressures and consequently increased thermal loads on combustion chamber walls close to the limits of endurance. Therefore accurate CFD simulation of conjugate heat transfer on the walls becomes a very important part of the development. In this study the heat transfer and temperature on piston surface was studied using conjugate heat transfer model along with a variety of near wall treatments for turbulence. New wall functions that account for variable density were implemented and tested against standard wall functions and against the hybrid near wall treatment readily available in a CFD software Star-CD.

Predictive models of soot formation during Diesel combustion are of great practical interest, particularly in light of newly proposed strict regulations on particulate emissions. A modified version of the phenomenological model of soot formation developed previously has been implemented in KIVA-II CFD code. The model includes major generic processes involved in soot formation during combustion, i.e., formation of soot precursors, formation of surface growth species, soot particle nucleation, coagulation, surface growth and oxidation. The formulation of the model within the KIVA-II is fully coupled with the mass and energy balances in the system. The model performance has been tested by comparison with the results of optical in-cylinder soot measurements in a single cylinder Cummins NH Diesel engine. The predicted soot volume fraction, number density and particle size agree reasonably well with the experimental data.

In the current investigation, the authors identified conditions under which increased in-cylinder turbulence can be used to improve diesel emissions. Two separate regimes of engine operation were identified; one in which combustion was constrained by mixing and one in which it was not. These regimes were dubbed under-mixed and over-mixed, respectively. It was found that increasing mixing in the former regime had a profound effect on soot emission. Fuel injection characteristics were found to be extremely important in determining the point at which mixing became inadequate. In addition, the ratio of the fuel injection momentum flux relative to that of the gas injection was found to be important in determining how increasing mixing would effect soot emissions.