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.

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.

Thermal test method for high power three-phase grid-connected inverters

Semiconductor device junction temperatures are maintained within datasheet specified limits to avoid failure in power converters. Burn-in tests are used to ensure this. In inverters, thermal time constants can be large and burn-in tests are required to be performed over long durations of time. At higher power levels, besides increased production cost, the testing requires sources and loads that can handle high power. In this study, a novel method to test a high power three-phase grid-connected inverter is proposed. The method eliminates the need for high power sources and loads. Only energy corresponding to the losses is consumed. The test is done by circulating rated current within the three legs of the inverter. All the phase legs being loaded, the method can be used to test the inverter in both cases of a common or independent cooling arrangement for the inverter phase legs. Further, the method can be used with different inverter configurations - three- or four-wire and for different pulse width modulation (PWM) techniques. The method has been experimentally validated on a 24 kVA inverter for a four-wire configuration that uses sine-triangle PWM and a three-wire configuration that uses conventional space vector PWM.

A controller design method for 3 phase 4 wire grid connected VSI with LCL filter

Closed loop control of a grid connected VSI requires line current control and dc bus voltage control. The closed loop system comprising PR current controller and grid connected VSI with LCL filter is a higher order system. Closed loop control gain expressions are therefore difficult to obtain directly for such systems. In this work a simplified approach has been adopted to find current and voltage controller gain expressions for a 3 phase 4 wire grid connected VSI with LCL filter. The closed loop system considered here utilises PR current controller in natural reference frame and PI controller for dc bus voltage control. Asymptotic frequency response plot and gain bandwidth requirements of the system have been used for current control and voltage controller design. A simplified lower order model, derived for closed loop current control, is used for the dc bus voltage controller design. The adopted design method has been verified through experiments by comparison of the time domain response.

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.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

Grid-refined LES predictions for Fan-OGV broadband noise

This work forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. In this paper, we focus on LES calculations of noise sources on and close to a blade trailing edge. We consider two test cases; one an isolated NACA0012 airfoil in flow, and the other an industry-standard rotating fan. In the first case, turbulent mean and RMS velocities and energy spectra at different locations are compared with those from experiment. 1,2The sound generated by the unsteady pressure fluctuations on the airfoil surface and by the flow turbulence will be predicted using a Ffowcs Williams Hawkings (FW-H) surface. In the second case, unsteady flow and acoustic fields around the blade passage 3 are presented for a refined mesh, and the rotor-stator tonal noise will be predicted by using the rotor-wake mean velocity profile and the methodology described in Lloyd & Peake 4. Copyright © 2009 by Qinling Li, Nigel Peake & Mark Savill.

Curvature and wrinkling of premixed flame Kernels comparison between planar laser induced fluorescence (PLIF) and 3D direct numerical simulations (DNS)

A direct comparison between time resolved PLIF measurements of OH and two dimensional slices from a full three dimensional DNS data set of turbulent premixed flame kernels in lean methane/air mixture was presented. The local flame structure and the degree of flame wrinkling were examined in response to differing turbulence intensities and turbulent Reynolds numbers. Simulations were performed using the SEGA DNS code, which is based on the solution of the compressible Navier Stokes, species, and energy equations for a lean hydrocarbon mixture. For the OH PLIF measurements, a cluster of four Nd:YAG laser was fired sequentially at high repetition rates and used to pump a dye laser. The frequency doubled laser beam was formed into a sheet of 40 mm height using a cylindrical telescope. The combination of PLIF and DNS has been demonstrated as a powerful tool for flame analysis. This research will form the basis for the development of sub-grid-scale (SGS) models for LES of lean-premixed combustion systems such as gas turbines. This is an abstract of a paper presented at the 30th International Symposium on Combustion (Chicago, IL 7/25-30/2004).

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

A Pressure-Correction Method And Its Applications On An Unstructured Chimera Grid

In this paper, an unstructured Chimera mesh method is used to compute incompressible flow around a rotating body. To implement the pressure correction algorithm on unstructured overlapping sub-grids, a novel interpolation scheme for pressure correction is proposed. This indirect interpolation scheme can ensure a tight coupling of pressure between sub-domains. A moving-mesh finite volume approach is used to treat the rotating sub-domain and the governing equations are formulated in an inertial reference frame. Since the mesh that surrounds the rotating body undergoes only solid body rotation and the background mesh remains stationary, no mesh deformation is encountered in the computation. As a benefit from the utilization of an inertial frame, tensorial transformation for velocity is not needed. Three numerical simulations are successfully performed. They include flow over a fixed circular cylinder, flow over a rotating circular cylinder and flow over a rotating elliptic cylinder. These numerical examples demonstrate the capability of the current scheme in handling moving boundaries. The numerical results are in good agreement with experimental and computational data in literature. (C) 2007 Elsevier Ltd. All rights reserved.

Sediment incipience in turbulence generated in a square tank by a vertically oscillating grid

The failure of hydraulic structures in many estuaries and coastal regions around the world has been attributed to sediment transport and local scour. The sediment incipience in homogenous turbulence generated by oscillating grid is studied in this paper. The turbulent flow is measured by particle tracer velocimetry (PTV) technique. The integral length scale and time scale of turbulence are obtained. The turbulent flow near the wall is measured by local optical magnification. The sediment incipience is described by static theory. The relationship of probability of sediment incipience and the turbulent kinetic energy were obtained experimentally and theoretically. The distribution of the turbulent kinetic energy near the wall is found to obey the power law and the turbulent energy is further identified as the dynamic mechanism of sediment incipience.