Analysis and simulation of combustion oscillations in a laboratory spray combustor

The unstable combustion that can occur in combustion chambers is a major problem for aeroengines and ground-based industrial gas turbines. Nowadays, CFD provides a flexible, low cost tool to supplement direct measurement. This paper presents simulations of combustion oscillations in a liquid-fuelled experimental rig at the University of Cambridge. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone and to develop inlet and outlet boundary conditions just upstream and downstream of the combustion region enabling the CFD calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in good agreement with experimental measurements. More details about the unstable combustion can be obtained from the simulation results. The approach developed here is expected to provide a powerful tool for the design and operation of stable combustion systems. Copyright © 2009 by ASME.

The use of object oriented modelling and simulation tools in constraint satisfaction for embodiment design

The software package Dymola, which implements the new, vendor-independent standard modelling language Modelica, exemplifies the emerging generation of object-oriented modelling and simulation tools. This paper shows how, in addition to its simulation capabilities, it may be used as an embodiment design tool, to size automatically a design assembled from a library of generic parametric components. The example used is a miniature model aircraft diesel engine. To this end, the component classes contain extra algebraic equations calculating the overload factor (or its reciprocal, the safety factor) for all the different modes of failure, such as buckling or tensile yield. Thus the simulation results contain the maximum overload or minimum safety factor for each failure mode along with the critical instant and the device state at which it occurs. The Dymola "Initial Conditions Calculation" function, controlled by a simple software script, may then be used to perform automatic component sizing. Each component is minimised in mass, subject to a chosen safety factor against failure, over a given operating cycle. Whilst the example is in the realm of mechanical design, it must be emphasised that the approach is equally applicable to the electrical or mechatronic domains, indeed to any design problem requiring numerical constraint satisfaction.

Applying quality control in building energy modelling: Comparative simulation of a high performance building

Quality control is considered from the simulator's perspective through comparative simulation of an ultra energy-efficient building with EE4-DOE2.1E and EnergyPlus. The University of Calgary's Leadership in Energy and Environmental Design Platinum Child Development Centre, with a 66% certified energy cost reduction rating, was the case study building. A Natural Resources Canada incentive program required use of EE4 interface with DOE2.1E simulation engine for energy modelling. As DOE2.1E lacks specific features to simulate advanced systems such as radiant cooling in the CDC, an EnergyPlus model was developed to further evaluate these features. The EE4-DOE2.1E model was used for quality control during development of the base EnergyPlus model and simulation results were compared. Advanced energy systems then added to the EnergyPlus model generated small difference in estimated total annual energy use. The comparative simulation process helped identify the main input errors in the draft EnergyPlus model. The comparative use of less complex simulation programs is recommended for quality control when producing more complex models. © 2009 International Building Performance Simulation Association (IBPSA).

A detailed chemistry simulation of the SI-HCCI transition

A Stochastic Reactor Model (SRM) has been used to simulate the transition from Spark Ignition (SI) mode to Homogeneous Charge Compression Ignition (HCCI) mode in a four cylinder in-line four-stroke naturally aspirated direct injection SI engine with cam profile switching. The SRM is coupled with GT-Power, a one-dimensional engine simulation tool used for modelling engine breathing during the open valve portion of the engine cycle, enabling multi-cycle simulations. The model is initially calibrated in both modes using steady state data from SI and HCCI operation. The mode change is achieved by switching the cam profiles and phasing, resulting in a Negative Valve Overlap (NVO), opening the throttle, advancing the spark timing and reducing the fuel mass as well as utilising a pilot injection. Experimental data is presented along with the simulation results. The model is used to investigate key control parameters and their effects on parameters that are difficult to measure experimentally. The effect of the spark in the first HCCI cycles is found to have a major impact on the stability of the transition. Copyright © 2010 SAE International.

Numerical simulation of post-flame oxidation of hydrocarbons in spark ignition engines

About 50-90 percent of the hydrocarbons that escape combustion during flame passage in spark-ignition engine operation are oxidized in the cylinder before leaving the system. The process involves the transport of unreacted fuel from cold walls towards the hotter burned gas regions and subsequent reaction. In order to understand controlling factors in the process, a transient one-dimensional reactive-diffusive model has been formulated for simulating the oxidation processes taking place in the reactive layer between hot burned gases and cold unreacted air/fuel mixture, with initial and boundary conditions provided by the emergence of hydrocarbons from the piston top land crevice. Energy and species conservation equations are solved for the entire process, using a detailed chemical kinetic mechanism for propane. Simulation results show that the post-flame oxidation process takes place within a reactive layer where intermediate hydrocarbon products are formed at temperatures above 1100-1200 K, followed by a carbon monoxide conversion region closer to the hot burned gases. Model results show that most of hydrocarbons leaving the crevice are completely oxidized inside the cylinder. The largest contribution of remaining hydrocarbons are those leaving the crevice at temperatures below 1400 K. The largest fraction of non-fuel (intermediate) hydrocarbons results from hydrocarbons leaving the crevice when core temperatures are around 1400 K Copyright © 1997 Society of Automotive Engineers, Inc.

Enhanced 10-Gbit/s link performance for directly Modulated complex-coupled DFB lasers via resonance frequency, damping rate, and chirp

This paper demonstrates the respective roles that combined index- and gain-coupling play in the overall link performance of distributed feedback (DFB) lasers. Their impacts on both static and dynamic properties such as slope efficiency, resonance frequency, damping rate, and chirp are investigated. Simulation results are compared with experimental data with good agreement. Transmission-oriented optimization is then demonstrated based on a targeted specification. The design tradeoffs are revealed, and it is shown that a modest combination of index- and gain-coupling enables optimum transmission at 10 Gbit/s.

Relay node placement for wireless sensor networks deployed in tunnels

Node placement plays a significant role in the effective and successful deployment of Wireless Sensor Networks (WSNs), i.e., meeting design goals such as cost effectiveness, coverage, connectivity, lifetime and data latency. In this paper, we propose a new strategy to assist in the placement of Relay Nodes (RNs) for a WSN monitoring underground tunnel infrastructure. By applying for the first time an accurate empirical mean path loss propagation model along with a well fitted fading distribution model specifically defined for the tunnel environment, we address the RN placement problem with guaranteed levels of radio link performance. The simulation results show that the choice of appropriate path loss model and fading distribution model for a typical environment is vital in the determination of the number and the positions of RNs. Furthermore, we adapt a two-tier clustering multi-hop framework in which the first tier of the RN placement is modelled as the minimum set cover problem, and the second tier placement is solved using the search-and-find algorithm. The implementation of the proposed scheme is evaluated by simulation, and it lays the foundations for further work in WSN planning for underground tunnel applications. © 2010 IEEE.

The ensemble statistics of the band-averaged energy of a random system

This paper is concerned with the response statistics of a dynamic system that has random properties. The frequency-band-averaged energy of the system is considered, and a closed form expression is derived for the relative variance of this quantity. The expression depends upon three parameters: the modal overlap factor m, a bandwidth parameter B, and a parameter α that defines the nature of the loading (for example single point forcing or rain-on-the-roof loading). The result is applicable to any single structural component or acoustic volume, and a comparison is made here with simulation results for a mass loaded plate. Good agreement is found between the simulations and the theory. © 2003 Published by Elsevier Ltd.

Directional algebraic reconstruction technique for electrical impedance tomography

We assume that the resistance matrix can be found in electrical impedance tomography from the assumption of linear dependence between the voltages and the currents and with the help of the resistance matrix and the transfer impedance between the electrodes, a directional algebraic reconstruction technique is proposed. The goal is to reconstruct the resistivity distribution by weighting the matrices that are obtained by calculating the orthogonal distance of the underlying mesh elements from the neighbouring port resistivity lines. These weighting matrices, which only depend on the topology of the underlying mesh, can be calculated offline and result in a computationally efficient online procedure with a reasonable image reconstruction performance. Simulation results are provided to validate this approach.

Enhanced performance of an all-optical parity checker using a single Mach-Zehnder interferometer

An all-optical parity checker is proposed that requires only a single Mach-Zehnder interferometer. Simulation results demonstrate an 8dB improvement in extinction ratio at 10 Gb/s operation. © 1999 Optical Society of America.

A novel heavily doped drift-auxiliary cathode lateral insulated gate transistor structure

A novel CMOS-compatible, heavily doped drift auxiliary cathode lateral insulated gate transistor (HDD-ACLIGT) structure is analyzed using two-dimensional device simulation techniques. Simulation results indicate that low on-resistance and a fast turn-off time of less than 50 ns can be achieved by incorporating an additional n+ region which is self-aligned to the gate between the p+ auxiliary cathode and the p well, together with an extended p buried layer in an anode-shorted modified lateral insulated gate transistor (MLIGT) structure. The on-state and its transient performance are analyzed in detail. The on-state performances of the HDD-ACLIGT and the MLIGT are compared and discussed. The results indicate that the HDD-ACLIGT structure is well suited for HVICs. The device is also well suited for integration with self-aligned digital CMOS.

Low short term timing jitter in an integrated monolithic extended cavity semiconductor mode-locked laser

Jitter measurements were performed on a monolithically integrated active/passive cavity multiple quantum well laser, actively mode-locked at 10 GHz via modulation of an absorber section. Sub-10 ps pulses were produced upon optimization of the drive conditions to the gain, distributed Bragg reflector, and absorber sections. A model was also developed using travelling wave rate equations. Simulation results suggest that spontaneous emission is the dominant cause of jitter, with carrier dynamics having a time constant of the order of 1 ns.

MOS and bipolar gated thyristor: A thyristor with IGBT switching characteristics

The IGBT has become the device of choice in many high-voltage-power electronic applications, by virtue of combining the ease of MOS gate control with an acceptable forward voltage drop. However, designers have retained an interest in MOS gated thyristor structures which have a turn-off capability. These offer low on-state losses as a result of their latching behaviour. Recently, there have been various proposals for dual-gate devices that have a thyristor on-state with IGBT-like switching. Many of these dual gated structures rely on advanced MOS technology, with inherent manufacturing difficulties. The MOS and bipolar gated thyristor offers all the advantages of dual gated performance, while employing standard IGBT processing techniques. The paper describes the MBGT in detail, and presents experimental and simulation results for devices based on realistic commercial processes. It is shown that the MBGT represents a viable power semiconductor device technology, suitable for a diverse range of applications. © IEE, 1998.

Vortex Transport and Blade Interactions in High Pressure Turbines

This paper presents a study of the three-dimensional flow field within the blade rows of a high-pressure axial flow steam turbine stage. Half-delta wings were fixed to a rotating hub to simulate an upstream rotor passage vortex. The flow field is investigated in a Low-Speed Research Turbine using pneumatic and hot-wire probes downstream of the blade row. The paper examines the impact of the delta wing vortex transport on the performance of the downstream blade row. Steady and unsteady numerical simulations were performed using structured 3D Navier-Stokes solver to further understand the flow field. The loss measurements at the exit of the stator blade showed an increase in stagnation pressure loss due to the delta wing vortex transport. The increase in loss was 21% of the datum stator loss, demonstrating the importance of this vortex interaction. The transport of the stator viscous flow through the rotor blade row is also described. The rotor exit flow was affected by the interaction between the enhanced stator passage vortex and the rotor blade row. Flow underturning near the hub and overturning towards the mid-span was observed, contrary to the classical model of overturning near the hub and underturning towards the mid-span. The unsteady numerical simulation results were further analysed to identify the entropy producing regions in the unsteady flow field.

Destruction-free parameter extraction for a physics-based circuit simulator IGCT model

This paper presents a practical destruction-free parameter extraction methodology for a new physics-based circuit simulator buffer-layer Integrated Gate Commutated Thyristor (IGCT) model. Most key parameters needed for this model can be extracted by one simple clamped inductive-load switching experiment. To validate this extraction method, a clamped inductive load switching experiment was performed, and corresponding simulations were carried out by employing the IGCT model with parameters extracted through the presented methodology. Good agreement has been obtained between the experimental data and simulation results.