Search Results

Sources of unburned hydrocarbon (UHC) emissions are examined for a highly dilute (10% oxygen concentration), moderately boosted (1.5 bar), low load (3.0 bar IMEP) operating condition in a single-cylinder, light-duty, optically accessible diesel engine undergoing partially-premixed low-temperature combustion (LTC). The evolution of the in-cylinder spatial distribution of UHC is observed throughout the combustion event through measurement of liquid fuel distributions via elastic light scattering, vapor and liquid fuel distributions via laser-induced fluorescence, and velocity fields via particle image velocimetry (PIV). The measurements are complemented by and contrasted with the predictions of multi-dimensional simulations employing a realistic, though reduced, chemical mechanism to describe the combustion process.

The interaction of fuel sprays and in-cylinder flow in direct-injection engines is expected to alter kinetic energy and integral length scales at least during some portions of the engine cycle. High-speed particle image velocimetry was implemented in an optical four-valve, pent-roof spark-ignition direct-injection single-cylinder engine to quantify this effect. Non-firing motored engine tests were performed at 1300 RPM with and without fuel injection. Two fuel injection timings were investigated: injection in early intake stroke represents quasi-homogenous engine condition; and injection in mid compression stroke mimics the stratified combustion strategy. Two-dimensional crank angle resolved velocity fields were measured to examine the kinetic energy and integral length scale through critical portions of the engine cycle. Reynolds decomposition was applied on the obtained engine flow fields to extract the fluctuations as an indicator for the turbulent flow.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

In-cylinder flow and combustion processes simulated with the standard k-ε turbulence model and with an alternative model-employing a non-linear, quadratic equation for the turbulent stresses-are contrasted for both motored and fired engine operation at two loads. For motored operation, the differences observed in the predictions of mean flow development are small and do not emerge until expansion. Larger differences are found in the spatial distribution and magnitude of turbulent kinetic energy. The non-linear model generally predicts lower energy levels and larger turbulent time scales. With fuel injection and combustion, significant differences in flow structure and in the spatial distribution of soot are predicted by the two models. The models also predict considerably different combustion efficiencies and NOx emissions.

In multidimensional modeling, fuels have been represented predominantly by single components, such as octane for gasoline. Several bicomponent studies have been performed, but these are still limited in their ability to represent real fuels, which are blends of as many as 300 components. This study outlines a method by which the fuel composition is represented by a distribution function of the fuel molecular weight. This allows a much wider range of compositions to be modeled, and only requires including two additional “species” besides the fuel, namely the mean and second moment of the distribution. This approach has been previously presented but is applied here to multidimensional calculations. Results are presented for single component droplet vaporization for comparison with single component fuel predictions, as well as results for a multicomponent gasoline and a diesel droplet.

The influence of the constant C3, which multiplies the mean flow divergence term in the model equation for the turbulent kinetic energy dissipation, is examined in a motored diesel engine for three different swirl ratios and three different spatial locations. Predicted temporal histories of turbulence energy and its dissipation are compared with experimentally-derived estimates. A “best-fit” value of C3 = 1.75, with an approximate uncertainty of ±0.3 is found to minimize the error between the model predictions and the experiments. Using this best-fit value, model length scale behavior corresponds well with that of measured velocity-correlation integral scales during compression. During expansion, the model scale grows too rapidly. Restriction of the model assessment to the expansion stroke suggests that C3 = 0.9 is more appropriate during this period.

Engine development is both time consuming and economically straining. Therefore, efforts are being made to optimize the research and development process for new engine technologies. The ability to apply information gained by studying an engine of one size/application to an engine of a completely different size/application would offer savings in both time and money in engine development. In this work, a computational study of diesel engine size-scaling relationships was performed to explore engine scaling parameters and the fundamental engine operating components that should be included in valid scaling arguments. Two scaling arguments were derived and tested: a simple, equal spray penetration scaling model and an extended, equal lift-off length scaling model. The simple scaling model is based on an equation for the conservation of mass and an equation for spray tip penetration developed by Hiroyasu et al. [1].

Engine design is a time consuming process in which many costly experimental tests are usually conducted. With increasing prediction ability of engine simulation tools, engine design aided by CFD software is being given more attention by both industry and academia. It is also of much interest to be able to use design information gained from an existing engine design of one size in the design of engines of other sizes to reduce design time and costs. Therefore it is important to study size-scaling relationships for engines over wide range of operating conditions. This paper presents CFD studies conducted for two production diesel engines - a light-duty GM-Fiat engine (0.5L displacement) and a heavy-duty Caterpillar engine (2.5L displacement). Previously developed scaling arguments, including an equal spray penetration scaling model and an extended, equal flame lift-off length scaling model were employed to explore the parametric scaling connections between the two engines.

Transient two-dimensional temperature distributions in the compression stroke and in the unburned end-gas of an SI engine were measured employing laser-induced fluorescence (LIF) of a fuel marker that possesses strongly temperature-dependent spectroscopic properties. The use of two different excitation wavelengths simplifies the otherwise complicated relation between LIF signal intensity and system parameters. The temperature fields obtained in this manner can be used to correct measured tracer-LIF maps and thus help to determine fuel distributions. Averaged temperature fields are compared to model calculations based on a homogeneous reactor assumption.

Various fluorescence tracers were assessed for their applicability for simultaneously measuring fuel distributions of different volatility classes. Tracers were chosen to show significantly different boiling behaviour representing three volatility classes of non-fluorescing multi-component fuels. Fluorescence properties of the markers were investigated using a heated static high-pressure cell with respect to emission behaviour, temperature and pressure dependence and quenching influences. A combination of ketonic and aromatic tracers appeared to be ideal for simultaneous imaging purposes since fluorescence is emitted in separate spectral regions with little overlap. Simultaneous measurements of the fuel distribution of two volatility classes were performed in a port fuel injected engine showing significant differences in the fuel distributions of low and mid boiling fractions in early stages of compression.

The effect of injection rate shape on spray evolution and emission characteristics is investigated and a methodology for active control of fuel injection is proposed. Extensive validation of advanced vaporization and primary jet breakup models was performed with experimental data before studying the effects of systematic changes of injection rate shape. Excellent agreement with the experiments was obtained for liquid and vapor penetration lengths, over a broad range of gas densities and temperatures. Also the predicted flame lift-off lengths of reacting diesel fuel sprays were in good agreement with the experiments. After the validation of the models, well-defined rate shapes were used to study the effect of injection rate shape on liquid and vapor penetration, flame lift-off lengths and emission characteristics.

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

In-cylinder fluid velocity is measured in an optically accessible, fired HSDI engine at idle. The velocity field is also calculated, including the full induction stroke, using multi-dimensional fluid dynamics and combustion simulation models. A detailed comparison between the measured and calculated velocities is performed to validate the computed results and to gain a physical understanding of the flow evolution. Motored measurements are also presented, to clarify the effects of the fuel injection process and combustion on the velocity field evolution. The calculated mean in-cylinder angular momentum (swirl ratio) and mean flow structures prior to injection agree well with the measurements. Modification of the mean flow by fuel injection and combustion is also well captured.

Cycle-resolved analysis of velocity data obtained in the re-entrant bowl of a fired high-;speed, direct-injection diesel engine, demonstrates an unambiguous, approximately 100% increase in late-cycle turbulence levels over the levels measured during motored operation. Model predictions of the flow field, obtained employing RNG k-ε turbulence modeling in KIVA-3V, do not capture this increased turbulence. A combined experimental and computational approach is taken to identify the source of this turbulence. The results indicate that the dominant source of the increased turbulence is associated with the formation of an unstable distribution of mean angular momentum, characterized by a negative radial gradient. The importance of this source of flow turbulence has not previously been recognized for engine flows. The enhanced late-cycle turbulence is found to be very sensitive to the flow swirl level.

The present study uses a numerical model to investigate the effects of flow turbulence on premixed iso-octane HCCI engine combustion. Different levels of in-cylinder turbulence are generated by using different piston geometries, namely a disc-shape versus a square-shape bowl. The numerical model is based on the KIVA code which is modified to use CHEMKIN as the chemistry solver. A detailed reaction mechanism is used to simulate the fuel chemistry. It is found that turbulence has significant effects on HCCI combustion. In the current engine setup, the main effect of turbulence is to affect the wall heat transfer, and hence to change the mixture temperature which, in turn, influences the ignition timing and combustion duration. The model also predicts that the combustion duration in the square bowl case is longer than that in the disc piston case which agrees with the measurements.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.

Developing a complete understanding of the structure and behavior of the near-wall region (NWR) in reciprocating, internal combustion (IC) engines and of its interaction with the core flow is needed to support the implementation of advanced combustion and engine operation strategies, as well as predictive computational models. The NWR in IC engines is fundamentally different from the canonical steady-state turbulent boundary layers (BL), whose structure, similarity and dynamics have been thoroughly documented in the technical literature. Motivated by this need, this paper presents results from the analysis of two-component velocity data measured with particle image velocimetry near the head of a single-cylinder, optical engine. The interaction between the NWR and the core flow was quantified via statistical moments and two-point velocity correlations, determined at multiple distances from the wall and piston positions.

The present study attempted to model experimental results obtained on an optical engine at the Sandia National Laboratory. Measurements of in-cylinder unburned hydrocarbon (UHC) distributions were provided using advanced optical diagnostics on a near production type piston. Previous multidimensional modeling provided accurate pressure profiles and heat release rate (HRR) predictions. However, the experimental UHC distribution was not matched, and the model predicted UHC extending from the bowl into the squish region in the expansion stroke. To explore the causes of this discrepancy a parametric study was performed using a variety of initial conditions, boundary conditions and model constants to explore their effects on the UHC distribution. Of the initial conditions, the swirl ratio was found to have the biggest impact on the UHC distribution.

RNG k-ε closure turbulence dissipation equations are evaluated employing the CFD code KIVA-3V Release 2. The numerical evaluations start by considering simple jet flows, including incompressible air jets and compressible helium jets. The results show that the RNG closure turbulence model predicts lower jet tip penetration than the "standard" k-ε model, as well as being lower than experimental data. The reason is found to be that the turbulence kinetic energy is dissipated too slowly in the downstream region near the jet nozzle exit. In this case, the over-predicted R term in RNG model becomes a sink of dissipation in the ε-equation. As a second step, the RNG turbulence closure dissipation models are further tested in complex engine flows to compare against the measured evolution of turbulence kinetic energy, and an estimate of its dissipation rate, during both the compression and expansion processes.

This presentation is an assessment of the turbulence-stress scale-similarity in an IC engine, which is used for modeling subgrid dissipation in LES. Residual stresses and Leonard stresses were computed after applying progressively smaller spatial filters to measured and simulated velocity distributions. The velocity was measured in the TCC-II engine using planar and stereo PIV taken in three different planes and with three different spatial resolutions, thus yielding two and three velocity components, respectively. Comparisons are made between the stresses computed from the measured velocity and stress computed from the LES resolved-scale velocity from an LES simulation. The results present the degree of similarity between the residual stresses and the Leonard stresses at adjacent scales. The specified filters are systematically reduced in size to the resolution limits of the measurements and simulation.