Fuel drop size measurement in a carburetted induction system with 2 throttling methods

A new throttling system far SI engines is examined. The SMD of the fuel droplets in the induction system is measured to evaluate the performance of the new device with respect to the conventional throttle plate arrangement. The measurements are conducted at steady now conditions. A forward angular scattering technique with a He-Ne laser beam is used for droplet size measurement. The experiments are carried out with different mixture strength, stream velocity and throttle positions. It is observed that A/F ratio has no effect on SMD. However, stream velocity and throttle position have a significant influence on SMD. The new throttling method is found to be more effective in reducing the SMD, particularly at low throttle opening and high stream velocity compared to the conventional throttle plate.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Vaporization and collision modeling of liquid fuel sprays in a co-axial fuel and air pre-mixer

Droplet collision occurs frequently in regions where the droplet number density is high. Even for Lean Premixed and Pre-vaporized (LPP) liquid sprays, the collision effects can be very high on the droplet size distributions, which will in turn affect the droplet vaporization process. Hence, in conjunction with vaporization modeling, collision modeling for such spray systems is also essential. The standard O'Rourke's collision model, usually implemented in CFD codes, tends to generate unphysical numerical artifact when simulations are performed on Cartesian grid and the results are not grid independent. Thus, a new collision modeling approach based on no-time-counter method (NTC) proposed by Schmidt and Rutland is implemented to replace O'Rourke's collision algorithm to solve a spray injection problem in a cylindrical coflow premixer. The so called ``four-leaf clover'' numerical artifacts are eliminated by the new collision algorithm and results from a diesel spray show very good grid independence. Next, the dispersion and vaporization processes for liquid fuel sprays are simulated in a coflow premixer. Two liquid fuels under investigation are jet-A and Rapeseed Methyl Esters (RME). Results show very good grid independence in terms of SMD distribution, droplet number distribution and fuel vapor mass flow rate. A baseline test is first established with a spray cone angle of 90 degrees and injection velocity of 3 m/s and jet-A achieves much better vaporization performance than RME due to its higher vapor pressure. To improve the vaporization performance for both fuels, a series of simulations have been done at several different combinations of spray cone angle and injection velocity. At relatively low spray cone angle and injection velocity, the collision effect on the average droplet size and the vaporization performance are very high due to relatively high coalescence rate induced by droplet collisions. Thus, at higher spray cone angle and injection velocity, the results expectedly show improvement in fuel vaporization performance since smaller droplet has a higher vaporization rate. The vaporization performance and the level of homogeneity of fuel-air mixture can be significantly improved when the dispersion level is high, which can be achieved by increasing the spray cone angle and injection velocity. (C) 2012 Elsevier Ltd. All rights reserved.

First and second law thermodynamic analysis of air and oxy-steam biomass gasification

Gasification is an energy transformation process in which solid fuel undergoes thermochemical conversion to produce gaseous fuel, and the two most important criteria involved in such process to evaluate the performance, economics and sustainability of the technology are: the total available energy (exergy) and the energy conserved (energy efficiency). Current study focuses on the energy and exergy analysis of the oxy-steam gasification and comparing with air gasification to optimize the H-2 yield, efficiency and syngas energy density. Casuarina wood is used as a fuel, and mixture of oxygen and steam in different proportion and amount is used as a gasifying media. The results are analysed with respect to varying equivalence ratio and steam to biomass ratio (SBR). Elemental mass balance technique is employed to ensure the validity of results. First and second law thermodynamic analysis is used towards time evaluation of energy and exergy analysis. Different component of energy input and output has been studied carefully to understand the influence of varying SBR on the availability of energy and irreversibility in the system to minimize the losses with change in input parameters for optimum performance. The energy and exergy losses (irreversibility) for oxy-steam gasification system are compared with the results of air gasification, and losses are found to be lower in oxy-steam thermal conversion; which has been argued and reasoned due to the presence of N-2 in the air-gasification. The maximum exergy efficiency of 85% with energy efficiency of 82% is achieved at SBR of 0.75 on the molar basis. It has been observed that increase in SBR results in lower exergy and energy efficiency, and it is argued to be due to the high energy input in steam generation and subsequent losses in the form of physical exergy of steam in the product gas, which alone accounts for over 18% in exergy input and 8.5% in exergy of product gas at SBR of 2.7. Carbon boundary point (CBP), is identified at the SBR of 1.5, and water gas shift (WGS) reaction plays a crucial role in H-2 enrichment after carbon boundary point (CBP) is reached. Effects of SBR and CBP on the H-2/CO ratio is analysed and discussed from the perspective of energy as well as the reaction chemistry. Energy density of syngas and energy efficiency is favoured at lower SBR but higher SBR favours H-2 rich gas at the expense of efficiency. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Optical diagnostics of fuel-air mixing and vortex formation in a cavity combustor

The current work reports optical diagnostic measurements of fuel-air mixing and vortex structure in a single cavity trapped vortex combustor (TVC). Specifically, the mixture fraction using acetone PLIF technique in the non-reacting flow, and PIV measurements in the reacting flow are reported for the first time in trapped vortex combustors. The fuel-air momentum flux ratio, where the air momentum corresponds to that entering the cavity through a specially-incorporated flow guide vane, is used to characterize the mixing. The acetone PLIF experiments show that at high momentum flux ratios, the fuel-air mixing in the cavity is very minimal and is enhanced as the momentum flux ratio reduces, due to a favourable vortex formation in the cavity. Stoichiometric mixture fraction surfaces show that the mixing causes the reaction surfaces to shift from non-premixed to partially-premixed stratified mixtures. PIV measurements conducted in the non-reacting flow in the cavity further reinforce this observation. The scalar dissipation rates of mixture fraction were compared with the contours of RMS of fluctuating velocity and showed very good agreement. The regions of maximum mixing are observed to be along the fuel air interface. Reacting flow Ply measurements which differ substantially from the non-reacting cases primarily because of the heat release from combustion and the resulting gas expansion show that the vortex is displaced from the centre of the cavity towards the guide vane. Overall, the measurements show interesting features of the flow including the presence of the dual cavity structure and lead to a clear understanding of the underlying physics of the cavity flow highlighting the importance of the fuel-air momentum ratio parameter. (C) 2014 Elsevier Inc. All rights reserved.

Fuel-optimal G-MPSP Guidance for Powered Descent Phase of Soft Lunar Landing

A fuel optimal nonlinear sub-optimal guidance scheme is presented in this paper for soft landing of a lunar craft during the powered descent phase. The recently developed Generalized Model Predictive Static Programming (G-MPSP) is used to compute the required magnitude and angle of the thrust vector. Both terminal position and velocity vector are imposed as hard constraints, which ensures high position accuracy and facilitates initiation of vertical descent at the end of the powered descent phase. A key feature of the G-MPSP algorithm is that it converts the nonlinear dynamic programming problem into a low-dimensional static optimization problem (of the same dimension as the output vector). The control history update is done in closed form after computing a time-varying weighting matrix through a backward integration process. This feature makes the algorithm computationally efficient, which makes it suitable for on-board applications. The effectiveness of the proposed guidance algorithm is demonstrated through promising simulation results.

Optimal Guidance for Accurate Lunar Soft Landing with Minimum Fuel Consumption using Model Predictive Static Programming

In this paper the soft lunar landing with minimum fuel expenditure is formulated as a nonlinear optimal guidance problem. The realization of pinpoint soft landing with terminal velocity and position constraints is achieved using Model Predictive Static Programming (MPSP). The high accuracy of the terminal conditions is ensured as the formulation of the MPSP inherently poses final conditions as a set of hard constraints. The computational efficiency and fast convergence make the MPSP preferable for fixed final time onboard optimal guidance algorithm. It has also been observed that the minimum fuel requirement strongly depends on the choice of the final time (a critical point that is not given due importance in many literature). Hence, to optimally select the final time, a neural network is used to learn the mapping between various initial conditions in the domain of interest and the corresponding optimal flight time. To generate the training data set, the optimal final time is computed offline using a gradient based optimization technique. The effectiveness of the proposed method is demonstrated with rigorous simulation results.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

The Head-On Colliding Process Of Binary Liquid Droplets At Low Velocity: High-Speed Photography Experiments And Modeling

The experimental and theoretical studies are reported in this paper for the head-on collisions of a liquid droplet with another of the same fluid resting on a solid substrate. The droplet on the hydrophobic polydimethylsiloxane (PDMS) substrate remains in a shape of an approximately spherical segment and is isometric to an incoming droplet. The colliding process of the binary droplets was recorded with high-speed photography. Head-on collisions saw four different types of response in our experiments: complete rebound, coalescence, partial rebound With conglutination, and coalescence accompanied by conglutination. For a complete rebound, both droplets exhibited remarkable elasticity and the contact time of the two colliding droplets was found to be in the range of 10-20 ms. With both droplets approximately considered as elastic bodies, Hertz contact theory was introduced to estimate the contact time for the complete rebound case. The estimated result Was found to be on the same order of magnitude as the experimental data, which indicates that the present model is reasonable. (C) 2008 Elsevier Inc. All rights reserved.

Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.