ADAPTATION OF MW LEVEL NATURAL GAS LEAN BURN ENGINE FOR PRODUCER GAS OPERATION A GRID CONNECTED POWER PLANT

Availability of producer gas engines at MW being limited necessitates to adapt engine from natural gas operation. The present work focus on the development of necessary kit for adapting a 12 cylinder lean burn turbo-charged natural gas engine rated at 900 kWe (Waukesha make VHP5904LTD) to operate on producer and set up an appropriate capacity biomass gasification system for grid linked power generation in Thailand. The overall plant configuration had fuel processing, drying, reactor, cooling and cleaning system, water treatment, engine generator and power evacuation. The overall project is designed for evacuation of 1.5 MWe power to the state grid and had 2 gasification system with the above configuration and 3 engines. Two gasification system each designed for about 1100 kg/hr of woody biomass was connected to the engine using a producer gas carburetor for the necessary Air to fuel ratio control. In the use of PG to fuel IC engines, it has been recognized that the engine response will differ as compared to the response with conventional fueled operation due to the differences in the thermo-physical properties of PG. On fuelling a conventional engine with PG, power de-rating can be expected due to the lower calorific value (LCV), lower adiabatic flame temperature (AFT) and the lower than unity product to reactant more ratio. Further the A/F ratio for producer gas is about 1/10th that of natural gas and requires a different carburetor for engine operation. The research involved in developing a carburetor for varying load conditions. The patented carburetor is based on area ratio control, consisting of a zero pressure regulator and a separate gas and air line along with a mixing zone. The 95 litre engine at 1000 rpm has an electrical efficiency of 33.5 % with a heat input of 2.62 MW. Each engine had two carburetors designed for producer gas flow each capable of handling about 1200 m3/hr in order to provide similar engine heat input at a lower conversion efficiency. Cold flow studies simulating the engine carburetion system results showed that the A/F was maintained in the range of 1.3 +/- 0.1 over the entire flow range. Initially, the gasification system was tested using woody biomass and the gas composition was found to be CO 15 +/- 1.5 % H-2 22 +/- 2% CH4 2.2 +/- 0.5 CO2 11.25 +/- 1.4 % and rest N-2, with the calorific value in the range of 5.0 MJ/kg. After initial trials on the engine to fine tune the control system and adjust various engine operating parameter a peak load of 800 kWe was achieved, while a stable operating conditions was found to be at 750 kWe which is nearly 85 % of the natural gas rating. The specific fuel consumption was found to be 0.9 kg of biomass per kWh.

SWEETENING OF BIOGAS IN SPRAY TOWERS FOR POWER GENERATION

This paper presents the experience of the new design of using impinging jet spray columns for scrubbing hydrogen sulfide from biogas that has been developed by Indian Institute of Science and patented. The process uses a chelated polyvalent metal ion which oxidizes the hydrogen sulfide to sulfur as a precipitate. The sulfur generated is filtered and the scrubbing liquid recycled after oxidation. The process involves in bringing contact the sour gas with chelated liquid in the spray columns where H2S reacts with chelated Fe3+ and precipitates as sulfur, whereas Fe3+ gets reduced to Fe2+. Fe2+ is regenerated to Fe3+ by reaction of oxygen in air in a separate packed column. The regenerated liquid is recirculated. Sulfur is filtered and separated as a byproduct. The paper presents the experience in using the spray towers for hydrogen sulfide removal and further use of the clean gas for generating power using gas engines. The maximum allowable limit of H2S for the gas engine is 200 ppm (v/v) in order to prevent any corrosion of engine parts and fouling of the lubricating oil. With the current ISET process, the hydrogen sulfide from the biogas is cleaned to less than 100 ppm (v/v) and the sweet gas is used for power generation. The system is designed for 550 NM3/hr of biogas and inlet H2S concentration of 2.5 %. The inlet concentration of the H2S is about 1 - 1.5 % and average measured outlet concentration is about 30 ppm, with an average gas flow of about 300 - 350 NM3/hr, which is the current gas production rate. The sweet gas is used for power generation in a 1.2 MWe V 12 engine. The average power generation is about 650 - 750 kWe, which is the captive load of the industry. The plant is a CHP (combined heat power) unit with heat from the cylinder cooling and flue being recovered for hot water and steam generation respectively. The specific fuel consumption is 2.29 kWh/m(3) of gas. The system has been in operation for more than 13,000 hours in last one year in the industry. About 8.4 million units of electricity has been generated scrubbing about 2.1 million m3 of gas. Performance of the scrubber and the engine is discussed at daily performance level and also the overall performance with an environment sustenance by precipitating over 27 tons of sulfur.

Stability of Filaments in Star-Forming Clouds and the Formation of Prestellar Cores in Them

The exact process(es) that generate(s) dense filaments which then form prestellar cores within them is unclear. Here we study the formation of a dense filament using a relatively simple set-up of a pressure-confined, uniform-density cylinder. We examine if its propensity to form a dense filament and further, to the formation of prestellar cores along this filament, bears on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised by the approximate equality between the mass line-density for this volume and its maximum value). Though self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then fragments along its length due to the growth of a Jeans-like instability to form prestellar cores. We suggest dense filaments in typical star-forming clouds classified as gravitationally super-critical under the assumption of: (i) isothermality when in fact, they are not, and (ii) extended radial profiles as against pressure-truncated, that significantly over-estimates their mass line-density, are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent surveys of nearby star-forming regions.

An Optimizing Code Generator for a Class of Lattice-Boltzmann Computations

The Lattice-Boltzmann method (LBM), a promising new particle-based simulation technique for complex and multiscale fluid flows, has seen tremendous adoption in recent years in computational fluid dynamics. Even with a state-of-the-art LBM solver such as Palabos, a user has to still manually write the program using library-supplied primitives. We propose an automated code generator for a class of LBM computations with the objective to achieve high performance on modern architectures. Few studies have looked at time tiling for LBM codes. We exploit a key similarity between stencils and LBM to enable polyhedral optimizations and in turn time tiling for LBM. We also characterize the performance of LBM with the Roofline performance model. Experimental results for standard LBM simulations like Lid Driven Cavity, Flow Past Cylinder, and Poiseuille Flow show that our scheme consistently outperforms Palabos-on average by up to 3x while running on 16 cores of an Intel Xeon (Sandybridge). We also obtain an improvement of 2.47x on the SPEC LBM benchmark.

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.

On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Onset of plane layer magnetoconvection at low Ekman number

We study the onset of magnetoconvection between two infinite horizontal planes subject to a vertical magnetic field aligned with background rotation. In order to gain insight into the convection taking place in the Earth's tangent cylinder, we target regimes of asymptotically strong rotation. The critical Rayleigh number Ra-c and critical wavenumber a(c) are computed numerically by solving the linear stability problem in a systematic way, with either stress-free or no-slip kinematic boundary conditions. A parametric study is conducted, varying the Ekman number E (ratio of viscous to Coriolis forces) and the Elsasser number. (ratio of the Lorentz force to the Coriolis force). E is varied from 10(-9) to 10(-2) and. from 10(-3) to 1. For a wide range of thermal and magnetic Prandtl numbers, our results verify and confirm previous experimental and theoretical results showing the existence of two distinct unstable modes at low values of E-one being controlled by the magnetic field, the other being controlled by viscosity (often called the viscous mode). It is shown that oscillatory onset does not occur in the range of parameters we are interested in. Asymptotic scalings for the onset of these modes are numerically confirmed and their domain of validity is precisely quantified. We show that with no-slip boundary conditions, the asymptotic behavior is reached for E < 10(-6) and establish a map in the (E, Lambda) plane. We distinguish regions where convection sets in either through the magnetic mode or through the viscous mode. Our analysis gives the regime in which the transition between magnetic and viscous modes may be observed. We also show that within the asymptotic regime, the role played by the kinematic boundary conditions is minimal. (C) 2015 AIP Publishing LLC.

Fracture Testing at Small-Length Scales: From Plasticity in Si to Brittleness in Pt

The field of micro-/nano-mechanics of materials has been driven, on the one hand by the development of ever smaller structures in devices, and, on the other, by the need to map property variations in large systems that are microstructurally graded. Observations of `smaller is stronger' have also brought in questions of accompanying fracture property changes in the materials. In the wake of scattered articles on micro-scale fracture testing of various material classes, this review attempts to provide a holistic picture of the current state of the art. In the process, various reliable micro-scale geometries are shown, challenges with respect to instrumentation to probe ever smaller length scales are discussed and examples from recent literature are put together to exhibit the expanse of unusual fracture response of materials, from ductility in Si to brittleness in Pt. Outstanding issues related to fracture mechanics of small structures are critically examined for plausible solutions.

Primary and rechargeable zinc-air batteries using ceramic and highly stable TiCN as an oxygen reduction reaction electrocatalyst

Primary and secondary zinc-air batteries based on ceramic, stable, one dimensional titanium carbonitride (TiCN) nanostructures are reported. The optimized titanium carbonitride composition by density functional theory reveals their good activity towards the oxygen reduction reaction (ORR). Electrochemical measurements show their superior performance for the ORR in alkaline media coupled with favourable kinetics. The nanostructured TiCN lends itself amenable to be used as an air cathode material in primary and rechargeable zinc-air batteries. The battery performance and cyclability are found to be good. Further, we have demonstrated a gel-based electrolyte for rechargeable zinc-air batteries based on a TiCN cathode under ambient, atmospheric conditions without any oxygen supply from a cylinder. The present cell can work at current densities of 10-20 mA cm(2) (app. 10 000 mA g(-1) of TiCN) for several hours (63 h in the case of 10 mA cm(-2)) with a charge retention of 98%. The low cost, noble metal-free, mechanically stable and corrosion resistant TiCN is a very good alternative to Pt for metal-air battery chemistry.

Control of vortex shedding at moderate Reynolds numbers

The suppression method of vortex shedding from a circular cylinder has been studied experimentally in the Reynolds number range from 300 to 1600. The test is performed in a water channel. The model cylinder is 1 cm in diameter and 38 cm in length. A row of small rods of 0.18 cm in diameter and 1.5 cm in length are perpendicularly connected to the surface of the model cylinder and distributed along the meridian, The distance between the neighboring rods and the angle of attack of the rods can be changed so that the suppression effect on vortex shedding can be adjusted. The results show that vortex shedding can be suppressed effectively if the distance between the neighboring rods is smaller than 3 times and the cylinder diameter and the angle of attack is in the range of 30degreesless than or equal tobeta<90&DEG;.

Johnson-Cook材料模型参数的实验测定

介绍了一种通过Hopkinson拉伸实验、圆筒爆炸试验和计算机仿真来确定Johnson-Cook材料模型中相关常数的方法。

Interferometric Measurement of Density in Nonstationary Shock Wave Reflection Flow and Comparison with CFD

We present density measurements from the application of interferometry and Fourier transform fringe analysis to the problem of nonstationary shock wave reflection over a semicircular cylinder and compare our experimental measurements to theoretical results from a CFD simulation of the same problem. The experimental results demonstrate our ability to resolve detailed structure in this complex shock wave reflection problem, allowing visualization of multiple shocks in the vicinity of the triple point, plus visualization of the shear layer and an associated vortical structure. Comparison between CFD and experiment show significant discrepancies with experiment producing a double Mach Reflection when CFD predicts a transitional Mach reflection.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

狭长圆柱型燃烧室内的燃烧不稳定性

研究了狭长圆柱型燃烧室内的喷雾火焰燃烧的不稳定性，为了更清楚地了解火焰的构造，首先测量了火焰的温度场。在较大的一次风和二次风变化范围内，测量了压力的振荡特性，得出其均方根值图。结果表明，火焰的稳定是由回流区完成的，在较小的一次风燃料当量比和中等的二次风量时，振荡强度达到100Pa左右。其频率为200-230Hz，与空气/燃料比值关系不大，分析还表明燃烧室中的振荡是轴向驻波振荡。

钝体尾流的双重小波包分解

根据湍流相干结构和非相干结构不相关的特性，提出了一种钝体尾流双重小波包分解的新算法，将湍流的运动分解成相干分量和非相干分量，该算法以湍流相干分量和非相干分量的相关系数作为迭代的控制指标，减小了过去算法中的随意性，用该算法对大长宽比的钝体尾流三维超声波流速仪测量数据的分析表明：1)钝体间距与宽度之比大于4时，钝体间的相互影响可以忽略不计；2)流线型的钝体尾流紊动强度较小。