Search Results

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.

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.

This paper reports an experimental and numerical investigation on the spatial and temporal liquid- and vapor-phase distributions of diesel fuel spray under engine-like conditions. The high pressure diesel spray was investigated in an optically-accessible constant volume combustion vessel for studying 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). Measurements were carried out by a high-speed imaging system capable of acquiring Mie-scattering and schlieren in a nearly simultaneous fashion mode 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 were 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.

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.

High-resolution planar laser-induced fluorescence (PLIF) measurements were performed on the scalar field in an optical engine. The measurements were of sufficient resolution to fully resolve all of the length scales of the flow field through the full cycle. The scalar dissipation spectrum was calculated, and by fitting the results to a model turbulent spectrum the Batchelor scale of the turbulent flow was estimated. The scalar inhomogeneity was introduced by a low-momentum gas jet injection. A consistent trend was observed in all data; the Batchelor scale showed a minimum value at top dead center (TDC) and was nearly symmetric about TDC. Increasing the engine speed resulted in a decrease of the Batchelor scale, and the presence of a shroud on the intake valve, which increased the turbulence intensity, also reduced the Batchelor scale. The effect of the shrouded valve was less significant compared to the effect of engine speed.

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

A transport equation residual model incorporating refined G-equation and detailed chemical kinetics combustion models has been developed and implemented in the ERC KIVA-3V release2 code for Gasoline Direct Injection (GDI) engine simulations for better predictions of flame propagation. In the transport equation residual model a fictitious species concept is introduced to account for the residual gases in the cylinder, which have a great effect on the laminar flame speed. The residual gases include CO2, H2O and N2 remaining from the previous engine cycle or introduced using EGR. This pseudo species is described by a transport equation. The transport equation residual model differentiates between CO2 and H2O from the previous engine cycle or EGR and that which is from the combustion products of the current engine cycle.

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.

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.

A 1-D 2-layer model developed previously at MTU was used in this research to predict the pressure drop, filtration characteristics and various properties of the particulate filter and the particulate deposit layer. The model was used along with dilute emission data to characterize two catalyzed particulate filters (CPFs) having different catalyst loading and catalyst application processes. The model was calibrated and validated with data obtained from steady state experiments conducted using a 1995 Cummins M11-330E heavy-duty diesel engine with manual EGR with different fuels for the two different CPFs. The two different catalyzed particulate filters were CPF III (5 gms/ft3 Pt) and CPF V (50 gms/ft3 Pt). Both the CPFs had cordierite substrates with CPF III and CPF V had MEX and NEX catalyst type formulation respectively. The CPF III filter was catalyzed using a solution-impregnated process while the CPF V filter was catalyzed using a wash coat process.

This investigation is a continuation of a previous study by these authors in which a Lagrange polynomial interpolation method was developed to evaluate spray source terms and also to distribute the source terms onto the gas mesh; the method was applied to the liquid-gas momentum exchange. For this investigation, the method has been extended to the mass exchange between the liquid and gas phases due to evaporation. The Lagrange polynomial interpolation and source term distribution methods are applied to the liquid-gas mass and momentum exchange and are evaluated for evaporating sprays using KIVA3 as a modeling platform. These methods are compared with the standard “nearest neighbor” method of KIVA3, and experimental data are used to establish their validity. The evaluation criteria used include the liquid and vapor spray penetration, gas velocities and the computational stability.

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.

Airflow characteristics surrounding evaporating transient diesel sprays inside a constant volume chamber under temperatures around 1100 K were investigated using a 6-hole injector and a single-hole injector. Particle Image Velocimetry (PIV) was used to measure the gas velocities surrounding a spray plume as a function of space and time. A conical control surface surrounding the spray plume was chosen as a representative side entrainment surface. The normal velocities crossing the control surface toward the spray plume for single-hole injection sprays were higher than those of 6-hole injection sprays. The velocities tangential to the control surface toward the injector tip for the single-hole injection sprays were lower than those of 6-hole injection sprays. An abrupt increase in tangential velocities near the chamber wall suggests that the recirculation of surrounding gas was accelerated by the spray wall impingement, both for non-evaporating and evaporating sprays.

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.

The effects of augmented mixing through the use of an auxiliary gas injection (AGI) were investigated in a direct-injection gasoline engine operated at a 22:1 overall air-fuel ratio, but with retarded injection timing such that the combustion was occurring in a locally rich mixture as evident by the elevated CO emissions. Two AGI gas compositions, nitrogen and air, were utilized, the gas supply temperature was ambient, and a wide range of AGI timings were investigated. The injected mass was less than 10% of the total chamber mass. The injection of nitrogen during the latter portion of the heat release phase resulted in a 25% reduction in the CO emissions. This reduction is considered to be the result of the increased mixing rate of the rich combustion products with the available excess air during a time when the temperatures are high enough to promote rapid oxidation.

The effect of ambient gas density and the number of injector holes on the characteristics of airflow surrounding non-evaporating transient diesel sprays inside a constant volume chamber were investigated using a 6-hole injector. Particle Image Velocimetry (PIV) was used to measure the gas velocities surrounding a spray plume as a function of space and time. A conical control surface surrounding the spray plume was chosen as a representative side entrainment surface. The positive normal velocities across the control surface of single-hole injection sprays were higher than those of 6-hole injection sprays. An abrupt increase in velocities tangential to the control surface near the chamber wall suggests that the recirculation of surrounding gas is accelerated by spray wall impingement.

Airflow characteristics surrounding non-evaporating transient diesel sprays were investigated using a 6-hole injector. Particle Image Velocimetry (PIV) was used to measure the gas velocities surrounding a spray plume as a function of space and time. A hydraulically actuated, electronically controlled unit injector (HEUI) system was used to supply the fuel into a pressurized constant volume chamber at room temperature. The chamber gas densities in this study were 10 kg/m3, 20 kg/m3 and 30 kg/m3. The injection pressure was 96.5 MPa. Two frequency doubled (532 nm) Nd:YAG lasers were used to create coincident laser sheets to illuminate the test section at two instances after start of injection (ASI). The double exposed images of sprays and Al2O3 seed particles were developed and velocity vectors of the gas surrounding the transient diesel sprays were obtained using a numerical autocorrelation PIV method.

Gas-phase in-cylinder mixing was examined by two different methods. The first method for observing mixing involved planar Mie scattering measurements of the instantaneous number density of silicon oil droplets which were introduced to the in-cylinder flow. The local value of the number density was assumed to be representative of the local gas concentration. Because the objective was to observe the rate in which gas concentration gradients change, to provide gradients in number density, droplets were admitted into the engine through only one of the two intake ports. Air only flowed through the other port. Three different techniques were used in analyzing the droplet images to determine the spatially dependent particle number density. Direct counting, a filtering technique, and autocorrelation were used and compared. Further, numerical experiments were performed with the autocorrelation method to check its effectiveness for determination of particle number density.

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.

This paper presents some experimental results of air velocity measurements near high pressure diesel sprays. The measurements were made using a moderately high pressure (90 MPa) common rail injector in a pressurized spray chamber. The chamber was operated at ambient temperature (25°C) and was pressurized with Argon to produce a chamber gas density of about 27 kg/m3, similar to densities found in a large turbocharged diesel near TDC. The gas phase was tagged using water droplets doped with Stilbene 420, with an estimated droplet size of 18 μm. The atomized water-Stilbene droplets were illuminated with the third harmonic of a pair of Nd:YAG lasers which caused the Stilbene to fluoresce at about 420 nm. To reduce the competing fluorescence from the injected fuel, the injector was fueled with Jet-A fuel. Using the two lasers, double exposures of the small droplets were recorded on film. The laser pulse lengths were about 6 ns, and typical times between pulses were 100 μs.