Search Results

Obtaining and maintaining a stratified charge in a practical engine is a difficult problem. Consequently, many approaches have been proposed and reported in the scientific and patent literature. In attempting to assess the most profitable approach for future development work, it is important to group together similar approaches so that one can study their performance as a group. Making such a classification has the additional advantage of helping to standardize terminology used by different investigators. With this thought in mind, a literature study was made and a proposed classification chart prepared for the different engine combustion systems reported in the literature. For the sake of completeness, the finally proposed classification chart includes homogeneous combustion engines as well as heterogeneous combustion engines. Because of their similarity of combustion, rotary engines such as the Wankel engine are considered as “reciprocating” although gas turbines are not included.

A computer program which rapidly calculates the equilibrium mole fractions and the partial derivatives of the mole fractions with respect to temperature, pressure and equivalence ratio for the products of combustion of any hydrocarbon fuel and air is described. A subroutine is also given which calculates the gas constant, enthalpy, internal energy and the partial derivatives of these with respect to temperature, pressure and equivalence ratio. Some examples of the uses of the programs are also given.

A linear parameter varying (LPV) reduced order model (ROM) is used to approximate the volume-averaged temperature of battery cells in one of the modules of the battery pack with varying mass flow rate of cooling fluid using uniform heat source as inputs. The ROM runs orders of magnitude faster than the original CFD model. To reduce the time it takes to generate training data, used in building LPV ROM, a divide-and-conquer approach is introduced. This is done by dividing the battery module into a series of mid-cell and end-cell units. A mid-cell unit is composed of a cooling channel sandwiched in between two half -cells. A half-cell has half as much heat capacity as a full-cell. An end-cell unit is composed of a cooling channel sandwiched in between full-cell and a half-cell. A mass flow rate distribution look-up-table is generated from a set of steady-state simulations obtained by running the full CFD model at different inlet manifold mass flow rate samples.

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.

Cycle-to-cycle variation in combustion phasing and combustion rate cause knock to occur differently in every cycle. This is found to be true even if the end gas thermo-chemical time history is the same. Three cycles are shown that have matched combustion phasing, combustion rate, and time of knock onset, but have knock intensity that differs by a factor of six. Thus, the prediction of knock intensity must include a stochastic component. It is shown that there is a relationship between the maximum possible knock intensity and the unburned fuel energy at the time of knock onset. Further, for a small window of unburned energy at knock onset, the probability density function of knock intensity is self similar when scaled by the 95th percentile of the cumulative distribution, and log-normal in shape.

The development of a high speed, multichannel data acquisition system is described. A precision magnetic tape recorder is used to record analog data from highly transient phenomena. Analog-to-digital data conversion is performed on a hybrid computer and the digitized data is processed using large, high speed digital computers. A detailed example of the application of the system to the measurement of rates-of-injection, rates-of-heat release, and instantaneous rates-of-heat transfer from the cylinder gases to the cylinder walls in a high speed open-chamber diesel engine is presented.

This paper presents the current state of a three-layer surface icing model for ice crystal icing risk assessment in aircraft engines, being developed jointly by Ansys and Honeywell to account for possible heat transfer from inside an engine into the flow path where ice accretion occurs. The bottom layer of the proposed model represents a thin metal sheet as a substrate surface to conductively transfer heat from an engine-internal reservoir to the ice layer. The middle layer is accretion ice with a porous structure able to hold a certain amount of liquid water. A shallow water film layer on the top receives impinged ice crystals. A mass and energy balance calculation for the film determines ice accretion rate. Water wicking and recovery is introduced to transfer liquid water between film layer and porous ice accretion layer.

Modeling ignition kernel development in spark ignition engines is crucial to capturing the sources of cyclic variability, both with RANS and LES simulations. Appropriate kernel modeling must ensure that energy transfer from the electrodes to the gas phase has the correct timing, rate and locations, until the flame surface is large enough to be represented on the mesh by the G-Equation level-set method. However, in most kernel models, geometric details driving kernel growth are missing: either because it is described as Lagrangian particles, or because its development is simplified, i.e., down to multiple spherical flames. This paper covers the geometric aspects of kernel development, which makes up the core of a Triangulated Lagrangian Ignition Kernel model. One (or multiple, if it restrikes) spark channel is initialized as a one-dimensional Lagrangian particle thread.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

Homogeneous low temperature combustion is believed to be a promising approach to resolve the conflict of goals between high efficiency and low exhaust emissions. Disadvantageously for this kind of combustion, the whole process depends on chemical kinetics and thus is hard to control. Reactivity controlled combustion can help to overcome this difficulty. In the so-called RCCI (reactivity controlled compression ignition) combustion concept a small amount of pilot diesel that is injected directly into the combustion chamber ignites a highly diluted gasoline-air mixture. As the gasoline does not ignite without the diesel, the pilot injection timing and the ratio between diesel and gasoline can be used to control the combustion process. A phenomenological multi-zone model to predict RCCI combustion has been developed and validated against experimental and 3D-CFD data. The model captures the main physics governing ignition and combustion.

An experimental investigation of non-intrusive combustion sensing was performed using a tri-axial accelerometer mounted to the engine block of a small-bore high-speed 4-cylinder compression-ignition direct-injection (CIDI) engine. This study investigates potential techniques to extract combustion features from accelerometer signals to be used for cycle-to-cycle engine control. Selection of accelerometer location and vibration axis were performed by analyzing vibration signals for three different locations along the block for all three of the accelerometer axes. A magnitude squared coherence (MSC) statistical analysis was used to select the best location and axis. Based on previous work from the literature, the vibration signal filtering was optimized, and the filtered vibration signals were analyzed. It was found that the vibration signals correlate well with the second derivative of pressure during the initial stages of combustion.

A methodology has been implemented to calculate local turbulent flame speeds for spark ignition engines accurately and on-the-fly in 3-D CFD modeling. The approach dynamically captures fuel effects, based on detailed chemistry calculations of laminar flame speeds. Accurately modeling flame propagation is critical to predicting heat release rates and emissions. Fuels used in spark ignition engines are increasingly complex, which necessitates the use of multi-component fuels or fuel surrogates for predictive simulation. Flame speeds of the individual components in these multi-component fuels may vary substantially, making it difficult to define flame speed values, especially for stratified mixtures. In addition to fuel effects, a wide range of local conditions of temperature, pressure, equivalence ratio and EGR are expected in spark ignition engines.

A digital computer is used to calculate the free motion of the valve train of a riged valve mechanism, after jump-off from the cam profile, in order to verify the correctness of the assumptions made with regard to the friction, in the dynamic analysis of the mechanism. Close correlation between the calculated and experimentally measured valve motions appears to justify the linearization of the friction present in the mechanism.

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.

Two main approaches are available when studying droplet dynamics for in-flight icing simulations: the Lagrangian approach, in which each droplet trajectory is integrated until it impacts the vehicle under study or when it leaves it behind without impact, and the Eulerian approach, where the droplet dynamics is solved as a continuum. In both cases, the same momentum equations are solved. Each approach has its advantages. In 2D, the Lagrangian approach is easy to code and it is very efficient, particularly when used in combination with a panel method flow solver. However, it is a far less practical approach for 3D simulations, particularly on complex geometries, as it is not an easy task to accurately determine the droplet seeding region without a great number of droplet trajectories, dramatically increasing the computing cost. Converting the impact locations into a water collection distribution is also a complex task, since droplet trajectories in 3D can follow convoluted paths.

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.

A CFD simulation methodology is presented to evaluate the ice that sheds from rotating components. The shedding detection is handled by coupling the ice accretion and stress analysis solvers to periodically check for the propagation of crack fronts and possible detachment. A novel approach for crack propagation is highlighted where no change in mesh topology is required. The entire computation from flow to impingement, ice accretion and crack analysis only requires a single mesh. The accretion and stress module are validated individually with published data. The analysis is extended to demonstrate potential shedding scenarios on three complex industrially-relevant 3D cases: a helicopter blade, an engine fan blade and a turboprop propeller. The largest shed fragment will be analyzed in the context of FOD damage to neighboring aircraft/component surfaces.

Deviations between transient and steady state operation of a modern light duty diesel engine were identified by comparing rapid load transitions to steady state tests at the same speeds and fueling rates. The validity of approximating transient performance by matching the transient charge air flow rate and intake manifold pressure at steady state was also assessed. Results indicate that for low load operation with low temperature combustion strategies, transient deviations of MAF and MAP from steady state values are small in magnitude or short in duration and have relatively little effect on transient engine performance. A new approximation accounting for variations in intake temperature and excess oxygen content of the EGR was more effective at capturing transient emissions trends, but significant differences in magnitudes remained in certain cases indicating that additional sources of variation between transient and steady state performance remain unaccounted for.

THIS paper is a sequel of the paper, “Photo-Electric Combustion Analysis,” presented at the 1936 Semi-Annual Meeting of the Society. The indicator described in that paper has been used to study combustion of 28 fuels and chemicals. A complete table of information of the materials used as fuels is included. The results obtained from over 1000 oscillograms show a different shape of ignition-lag curve versus injection advance angle than it is ordinarily thought to have. Even though the cetane values for these 28 fuels varied from 24 to 100, they all had nearly the same ignition lag when injected near the dead-center position. This minimum value is shown to be about 1/1000 sec. The fuels of higher-cetane value reach this minimum at an earlier injection angle than do those of low-cetane value. The paper shows how a high-cetane fuel can be just as rough as a low-cetane fuel if the injection timing is too early.

Diesel piston bowl geometry can affect turbulent mixing and therefore it impacts heat-release rates, thermal efficiency, and soot emissions. The focus of this work is on the effects of bowl geometry and injection timing on turbulent flow structure. This computational study compares engine behavior with two pistons representing competing approaches to combustion chamber design: a conventional, re-entrant piston bowl and a stepped-lip piston bowl. Three-dimensional computational fluid dynamics (CFD) simulations are performed for a part-load, conventional diesel combustion operating point with a pilot-main injection strategy under non-combusting conditions. Two injection timings are simulated based on experimental findings: an injection timing for which the stepped-lip piston enables significant efficiency and emissions benefits, and an injection timing with diminished benefits compared to the conventional, re-entrant piston.