Simulation studies on the performance of thermoacoustic prime movers and refrigerator

The prime movers and refrigerators based on thermoacoustics have gained considerable importance toward practical applications in view of the absence of moving components, reasonable efficiency, use of environmental friendly working fluids, etc. Devices such as twin Standing Wave ThermoAcoustic Prime Mover (SWTAPM), Traveling Wave ThermoAcoustic Prime Mover (TWTAPM) and thermoacoustically driven Standing Wave ThermoAcoustic Refrigerator (SWTAR) have been studied by researchers. The numerical modeling and simulation play a vital role in their development. In our efforts to build the above thermoacoustic systems, we have carried out numerical analysis using the procedures of CFD on the above systems. The results of the analysis are compared with those of DeltaEC (freeware from LANL, USA) simulations and the experimental results wherever possible. For the CFD analysis commercial code Fluent 6.3.26 has been used along with the necessary boundary conditions for different working fluids at various average pressures. The results of simulation indicate that choice of the working fluid and the average pressure are critical to the performance of the above thermoacoustic devices. Also it is observed that the predictions through the CFD analysis are closer to the experimental results in most cases, compared to those of DeltaEC simulations. (C) 2015 Elsevier Ltd. All rights reserved.

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.

``Conformational Simulation'' of Sulfamethizole by Molecular Complexation and Insights from Charge Density Analysis: Role of Intramolecular S center dot center dot center dot O Chalcogen Bonding

Intramolecular S center dot center dot center dot O chalcogen bonding and its potential to lock molecular conformation have been examined in the crystal forms of sulfamethizole, a sulfonamide antibiotic. Molecular complexes of sulfamethizole, including salts and cocrystal, have been synthesized, and their crystal structures were analyzed in order to examine the possible conformational preferences of the molecule in various ionic states and supramolecular environments (neutral/cocrystal, anionic salt, and cationic salt forms). The electrostatic potential mapped on Hirshfeld surfaces generated for these crystal forms provides insights into the possible binding modes of the drug in different environments. Further, the observed conformation locking feature has been rationalized in terms of the experimental charge density features of the intramolecular S center dot center dot O chalcogen bonding in sulfamethizole. The study quantitatively illustrates and rationalizes an intriguing case of a local minimum of molecular conformation being exclusively preferred over the global minimum, as it facilitates more efficient intermolecular interactions in a supramolecular environment.

A literature survey on the real-time computational simulation of biological organs

This paper lists some references that could in some way be relevant in the context of the real-time computational simulation of biological organs, the research area being defined in a very broad sense. This paper contains 198 references.

Simulation studies on Performance parameters of a TPC Polarimeter

X-ray polarimeters based on Time Projection Chamber (TPC) geometry are currently being studied and developed to make sensitive measurement of polarization in 2-10keV energy range. TPC soft X-ray polarimeters exploit the fact that emission direction of the photoelectron ejected via photoelectric effect in a gas proportional counter carries the information of the polarization of the incident X-ray photon. Operating parameters such as pressure, drift field and drift-gap affect the performance of a TPC polarimeter. Simulations presented here showcase the effect of these operating parameters on the modulation factor of the TPC polarimeter. Models of Garfield are used to study photoelectron interaction in gas and drift of electron cloud towards Gas Electron Multiplier (GEM). The emission direction is reconstructed from the image and modulation factor is computed. Our study has shown that Ne/DME (50/50) at lower pressure and drift field can be used for a TPC polarimeter with modulation factor of 50-65%.

Drift parameters optimization of a TPC polarimeter: a simulation study

Time Projection Chamber (TPC) based X-ray polarimeters using Gas Electron Multiplier (GEM) are currently being developed to make sensitive measurement of polarization in 2-10 keV energy range. The emission direction of the photoelectron ejected via photoelectric effect carries the information of the polarization of the incident X-ray photon. Performance of a gas based polarimeter is affected by the operating drift parameters such as gas pressure, drift field and drift-gap. We present simulation studies carried out in order to understand the effect of these operating parameters on the modulation factor of a TPC polarimeter. Models of Garfield are used to study photoelectron interaction in gas and drift of electron cloud towards GEM. Our study is aimed at achieving higher modulation factors by optimizing drift parameters. Study has shown that Ne/DME (50/50) at lower pressure and drift field can lead to desired performance of a TPC polarimeter.

MODELING AND SIMULATION OF A THREE-WHEELED MOBILE ROBOT ON UNEVEN TERRAINS WITH TWO-DEGREE-OF-FREEDOM SUSPENSION MECHANISMS

A wheeled mobile robot (WMR) will move on an uneven terrain without slip if its torus-shaped wheels tilt in a lateral direction. An independent two degree-of-freedom (DOF) suspension is required to maintain contact with uneven terrain and for lateral tilting. This article deals with the modeling and simulation of a three-wheeled mobile robot with torus-shaped wheels and four novel two-DOF suspension mechanism concepts. Simulations are performed on an uneven terrain for three representative pathsa straight line, a circular, and an S'-shaped path. Simulations show that a novel concept using double four-bar mechanism performs better than the other three concepts.

Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

A mobile robot with a two-degree-of-freedom suspension for traversing uneven terrain with minimal slip: Modeling, simulation and experiments

It is known in literature that a wheeled mobile robot (WMR) with fixed length axle will slip on an uneven terrain. One way to avoid wheel slip is to use a torus-shaped wheel with lateral tilt capability which allows the distance between the wheel-ground contact points to change even with a fixed length axle. Such an arrangement needs a two degree-of-freedom (DOF) suspension for the vertical and lateral tilting motion of the wheel. In this paper modeling, simulation, design and experimentation with a three-wheeled mobile robot, with torus-shaped wheels and a novel two DOF suspension allowing independent lateral tilt and vertical motion, is presented. The suspension is based on a four-bar mechanism and is called the double four-bar (D4Bar) suspension. Numerical simulations show that the three-wheeled mobile robot can traverse uneven terrain with low wheel slip. Experiments with a prototype three-wheeled mobile robot moving on a constructed uneven terrain along a straight line, a circular arc and a path representing a lane change, also illustrate the low slip capability of the three-wheeled mobile robot with the D4Bar suspension. (C) 2015 Elsevier Ltd. All rights reserved.

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.

Lower Bound Axisymmetric Limit Analysis Using Drucker-Prager Yield Cone and Simulation with Mohr-Coulomb Pyramid

This paper presents a lower bound limit analysis approach for solving an axisymmetric stability problem by using the Drucker-Prager (D-P) yield cone in conjunction with finite elements and nonlinear optimization. In principal stress space, the tip of the yield cone has been smoothened by applying the hyperbolic approximation. The nonlinear optimization has been performed by employing an interior point method based on the logarithmic barrier function. A new proposal has also been given to simulate the D-P yield cone with the Mohr-Coulomb hexagonal yield pyramid. For the sake of illustration, bearing capacity factors N-c, N-q and N-gamma have been computed, as a function of phi, both for smooth and rough circular foundations. The results obtained from the analysis compare quite well with the solutions reported from literature.

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.

A novel index for identification of weak nodes for reactive compensation to improve voltage stability

In this study, a new reactive power loss index (RPLI) is proposed for identification of weak buses in the system. This index is further used for determining the optimal locations for placement of reactive compensation devices in the power system for additional voltage support. The new index is computed from the reactive power support and loss allocation algorithm using Y-bus method for the system under intact condition and as well as critical/severe network contingencies cases. Fuzzy logic approach is used to select the important and critical/severe line contingencies from the contingency list. The inherent characteristics of the reactive power in system operation is properly addressed while determining the reactive power loss allocation to load buses. The proposed index is tested on sample 10-bus equivalent system and 72-bus practical equivalent system of Indian southern region power grid. The validation of the weak buses identification from the proposed index with that from other existing methods in the literature is carried out to demonstrate the effectiveness of the proposed index. Simulation results show that the identification of weak buses in the system from the new RPLI is completely non-iterative, thus requires minimal computational efforts as compared with other existing methods in the literature.

Simulation of Effective Thermal Conductivity of Metal Hydride Packed Beds

Several mathematical models are available for estimation of effective thermal conductivity of nonreactive packed beds. Keeping in view the salient differences between metal hydride beds in which chemisorption of hydrogen takes place and conventional nonreactive packed beds, modified models are proposed here to predict the effective thermal conductivity. Variation in properties such as solid thermal conductivity and porosity during hydrogen absorption and desorption processes are incorporated. These extended models have been applied to simulate the effective thermal conductivity of the MmNi(4.5)Al(0.5) hydride bed and are compared with the experimental results. Applicability of the extended models for estimation of the effective thermal conductivity at different operating conditions such as pressure, temperature, and hydrogen concentration is discussed.

Computationally efficient models for simulation of non-ideal DC-DC converters operating in continuous and discontinuous conduction modes

This paper discusses dynamic modeling of non-isolated DC-DC converters (buck, boost and buck-boost) under continuous and discontinuous modes of operation. Three types of models are presented for each converter, namely, switching model, average model and harmonic model. These models include significant non-idealities of the converters. The switching model gives the instantaneous currents and voltages of the converter. The average model provides the ripple-free currents and voltages, averaged over a switching cycle. The harmonic model gives the peak to peak values of ripple in currents and voltages. The validity of all these models is established by comparing the simulation results with the experimental results from laboratory prototypes, at different steady state and transient conditions. Simulation based on a combination of average and harmonic models is shown to provide all relevant information as obtained from the switching model, while consuming less computation time than the latter.