Influence of fuel hydrogen fraction on syngas fueled SI engine: Fuel thermo-physical property analysis and in-cylinder experimental investigations

Hydrogen, either in pure form or as a gaseous fuel mixture specie enhances the fuel conversion efficiency and reduce emissions in an internal combustion engine. This is due to the reduction in combustion duration attributed to higher laminar flame speeds. Hydrogen is also expected to increase the engine convective heat flux, attributed (directly or indirectly) to parameters like higher adiabatic flame temperature, laminar flame speed, thermal conductivity and diffusivity and lower flame quenching distance. These factors (adversely) affect the thermo-kinematic response and offset some of the benefits. The current work addresses the influence of mixture hydrogen fraction in syngas on the engine energy balance and the thermo-kinematic response for close to stoichiometric operating conditions. Four different bio-derived syngas compositions with fuel calorific value varying from 3.14 MJ/kg to 7.55 MJ/kg and air fuel mixture hydrogen fraction varying from 7.1% to 14.2% by volume are used. The analysis comprises of (a) use of chemical kinetics simulation package CHEMKIN for quantifying the thermo-physical properties (b) 0-D model for engine in-cylinder analysis and (c) in-cylinder investigations on a two-cylinder engine in open loop cooling mode for quantifying the thermo-kinematic response and engine energy balance. With lower adiabatic flame temperature for Syngas, the in-cylinder heat transfer analysis suggests that temperature has little effect in terms of increasing the heat flux. For typical engine like conditions (700 K and 25 bar at CR of 10), the laminar flame speed for syngas exceeds that of methane (55.5 cm/s) beyond mixture hydrogen fraction of 11% and is attributed to the increase in H based radicals. This leads to a reduction in the effective Lewis number and laminar flame thickness, potentially inducing flame instability and cellularity. Use of a thermodynamic model to assess the isolated influence of thermal conductivity and diffusivity on heat flux suggests an increase in the peak heat flux between 2% and 15% for the lowest (0.420 MW/m(2)) and highest (0.480 MW/m(2)) hydrogen containing syngas over methane (0.415 MW/m(2)) fueled operation. Experimental investigations indicate the engine cooling load for syngas fueled engine is higher by about 7% and 12% as compared to methane fueled operation; the losses are seen to increase with increasing mixture hydrogen fraction. Increase in the gas to electricity efficiency is observed from 18% to 24% as the mixture hydrogen fraction increases from 7.1% to 9.5%. Further increase in mixture hydrogen fraction to 14.2% results in the reduction of efficiency to 23%; argued due to the changes in the initial and terminal stages of combustion. On doubling of mixture hydrogen fraction, the flame kernel development and fast burn phase duration decrease by about 7% and 10% respectively and the terminal combustion duration, corresponding to 90%-98% mass burn, increases by about 23%. This increase in combustion duration arises from the cooling of the near wall mixture in the boundary layer attributed to the presence of hydrogen. The enhancement in engine cooling load and subsequent reduction in the brake thermal efficiency with increasing hydrogen fraction is evident from the engine energy balance along with the cumulative heat release profiles. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

After transition in a soft-walled microchannel

In comparison to the flow in a rigid channel, there is a multifold reduction in the transition Reynolds number for the flow in a microchannel when one of the walls is made sufficiently soft, due to a dynamical instability induced by the fluid-wall coupling, as shown by Verma & Kumaran (J. Fluid Mech., vol. 727, 2013, pp. 407-455). The flow after transition is characterised using particle image velocimetry in the x-y plane, where x is the streamwise direction and y is the cross-stream coordinate along the small dimension of the channel of height 0.2-0.3 mm. The flow after transition is characterised by a mean velocity profile that is flatter at the centre and steeper at the walls in comparison to that for a laminar flow. The root mean square of the streamwise fluctuating velocity shows a characteristic sharp increase away from the wall and a maximum close to the wall, as observed in turbulent flows in rigid-walled channels. However, the profile is asymmetric, with a significantly higher maximum close to the soft wall in comparison to that close to the hard wall, and the Reynolds stress is found to be non-zero at the soft wall, indicating that there is a stress exerted by fluid velocity fluctuations on the wall. The maximum of the root mean square of the velocity fluctuations and the Reynolds stress (divided by the fluid density) in the soft-walled microchannel for Reynolds numbers in the range 250-400, when scaled by suitable powers of the maximum velocity, are comparable to those in a rigid channel at Reynolds numbers in the range 5000-20 000. The near-wall velocity profile shows no evidence of a viscous sublayer for (y upsilon(*)/nu) as low as two, but there is a logarithmic layer for (y upsilon(*)/nu) up to approximately 30, where the von Karman constants are very different from those for a rigid-walled channel. Here, upsilon(*) is the friction velocity, nu is the kinematic viscosity and y is the distance from the soft surface. The surface of the soft wall in contact with the fluid is marked with dye spots to monitor the deformation and motion along the fluid-wall interface. Low-frequency oscillations in the displacement of the surface are observed after transition in both the streamwise and spanwise directions, indicating that the velocity fluctuations are dynamically coupled to motion in the solid.

PREDICTION OF TRANSITIONAL AND SEPARATED BOUNDARY LAYERS IN A COMPRESSOR CASCADE

Numerical simulations were performed of experiments from a cascade of stator blades at three low Reynolds numbers representative of flight conditions. Solutions were assessed by comparing blade surface pressures, velocity and turbulence intensity along blade normals at several stations along the suction surface and in the wake. At Re = 210,000 and 380,000 the laminar boundary layer over the suction surface separates and reattaches with significant turbulence fluctuations. A new 3-equation transition model, the k-k(L)-omega model, was used to simulate this flow. Predicted locations of the separation bubble, and profiles of velocity and turbulence fluctuations on blade-normal lines at various stations along the blade were found to be quite close to measurements. Suction surface pressure distributions were not as close at the lower Re. The solution with the standard k-omega SST model showed significant differences in all quantities. At Re = 640,000 transition occurs earlier and it is a turbulent boundary layer that separates near the trailing edge. The solution with the Reynolds stress model was found to be quite close to the experiment in the separated region also, unlike the k-omega SST solution. Three-dimensional computations were performed at Re = 380,000 and 640,000. In both cases there were no significant differences between the midspan solution from 3D computations and the 2D solutions. However, the 3D solutions exhibited flow features observed in the experiments the nearly 2D structure of the flow over most of the span at 380,000 and the spanwise growth of corner vortices from the endwall at 640,000.

Quantitative In Situ Analysis of Deformation in Sliding Metals: Effect of Initial Strain State

Using in situ, high-speed imaging of a hard wedge sliding against pure aluminum, and image analysis by particle image velocimetry, the deformation field in sliding is mapped at high resolution. This model system is representative of asperity contacts on engineered surfaces and die-workpiece contacts in deformation and machining processes. It is shown that large, uniform plastic strains of 1-5 can be imposed at the Al surface, up to depths of 500 mu m, under suitable sliding conditions. The spatial strain and strain rate distributions are significantly influenced by the initial deformation state of the Al, e.g., extent of work hardening, and sliding incidence angle. Uniform straining occurs only under conditions of steady laminar flow in the metal. Large pre-strains and higher sliding angles promote breakdown in laminar flow due to surface fold formation or flow localization in the form of shear bands, thus imposing limits on uniform straining by sliding. Avoidance of unsteady sliding conditions, and selection of parameters like sliding angle, thus provides a way to control the deformation field. Key characteristics of the sliding deformation such as strain and strain rate, laminar flow, folding and prow formation are well predicted by finite element simulation. The deformation field provides a quantitative basis for interpreting wear particle formation. Implications for engineering functionally graded surfaces, sliding wear and ductile failure in metals are discussed.

Numerical Study of Counterflow Diffusion Flame and Water Spray Interaction

This paper reports numerical investigation concerning the interaction of a laminar methane-air counterflow diffusion flame with monodisperse and polydisperse water spray. Commercial code ANSYS FLUENT with reduced chemistry has been used for investigation. Effects of strain rate, Sauter mean diameter (SMD), and droplet size distribution on the temperature along stagnation streamline have been studied. Flame extinction using polydisperse water spray has also been explored. Comparison of monodisperse and polydisperse droplet distribution on flame properties reveals suitability of polydisperse spray in flame temperature reduction beyond a particular SMD. This study also provides a numerical framework to study flame-spray interaction and extinction.