Pressure data analysis of formula one racing engines

The measured time-history of the cylinder pressure is the principal diagnostic in the analysis of processes within the combustion chamber. This paper defines, implements and tests a pressure analysis algorithm for a Formula One racing engine in MATLAB1. Evaluation of the software on real data is presented. The sensitivity of the model to the variability of burn parameter estimates is also discussed. Copyright © 1997 Society of Automotive Engineers, Inc.

Time resolved measurement of cold start hc concentration using the fast FID

Understanding mixture formation phenomena during the first few cycles of an engine cold start is extremely important for achieving the minimum engine-out emission levels at the time when the catalytic converter is not yet operational. Of special importance is the structure of the charge (film, droplets and vapour) which enters the cylinder during this time interval as well as its concentration profile. However, direct experimental studies of the fuel behaviour in the inlet port have so far been less than fully successful due to the brevity of the process and lack of a suitable experimental technique. We present measurements of the hydrocarbon (HC) concentration in the manifold and port of a production SI engine using the Fast Response Flame Ionisation Detector (FRFID). It has been widely reported in the past few years how the FRFID can be used to study the exhaust and in-cylinder HC concentrations with a time resolution of a few degrees of crank angle, and the device has contributed significantly to the understanding of unburned HC emissions. Using the FRFID in the inlet manifold is difficult because of the presence of liquid droplets, and the low and fluctuating pressure levels, which leads to significant changes in the response time of the instrument. However, using recently developed procedures to correct for the errors caused by these effects, the concentration at the sampling point can be reconstructed to align the FRFID signal with actual events in the engine. © 1996 Society of Automotive Engineers, Inc.

Particulate matter emission during start-up and transient operation of a spark-ignition engine

In order to understand why emissions of Particulate Matter (PM) from Spark-Ignition (SI) automobiles peak during periods of transient operation such as rapid accelerations, a study of controlled, repeatable transients was performed. Time-resolved engine-out PM emissions from a modern four-cylinder engine during transient load and air/fuel ratio operation were examined, and the results could be fit in most cases to a first order time response. The time constants for the transient response are similar to those measured for changes in intake valve temperature, reflecting the strong dependence of PM emissions on the amount of liquid fuel in the combustion chamber. In only one unrepeatable case did the time response differ from a first order function: showing an overshoot in PM emissions during transition from the initial to the final steady state PM emission level. PM emissions during controlled, motored start-up experiments show a peak at start-up followed by a period during which emissions are either relatively constant or drift somewhat. When the fuel injection and ignition are shut off, PM emissions also peak briefly, but rapidly decay to low levels. Qualitative implications on the study and modeling of PM emissions during transient engine operation are discussed. Copyright © 1999 Society of Automotive Engineers, Inc.

Liquid fuel visualization using laser-induced fluoresence during cold start

The presence of liquid fuel inside the engine cylinder is believed to be a strong contributor to the high levels of hydrocarbon emissions from spark ignition (SI) engines during the warm-up period. Quantifying and determining the fate of the liquid fuel that enters the cylinder is the first step in understanding the process of emissions formation. This work uses planar laser induced fluorescence (PLIF) to visualize the liquid fuel present in the cylinder. The fluorescing compounds in indolene, and mixtures of iso-octane with dopants of different boiling points (acetone and 3-pentanone) were used to trace the behavior of different volatility components. Images were taken of three different planes through the engine intersecting the intake valve region. A closed valve fuel injection strategy was used, as this is the strategy most commonly used in practice. Background subtraction and masking were both performed to reduce the effect of any spurious fluorescence. The images were analyzed on both a time and crank angle (CA) basis, showing the time of maximum liquid fuel present in the cylinder and the effect of engine events on the inflow of liquid fuel. The results show details of the liquid fuel distribution as it enters the engine as a function of crankangle degree, volatility and location in the cylinder. A. semi-quantitative analysis based on the integration of the image intensities provides additional information on the temporal distribution of the liquid fuel flow. © 1998 Society of Automotive Engineers, Inc.

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.

Extent of oxidation of hydrocarbons desorbing from the lubricant oil layer in spark-ignition engines

The extent of oxidation of hydrocarbons desorbing from the oil layer has been measured directly in a hydrogen-fueled, spark-ignited engine in which the lubricant oil was doped with a single component hydrocarbon. The amount of hydrocarbon desorbed and oxidized could be measured simultaneously as the dopant was only source of carbon-containing species. The fraction oxidized was strongly dependent on engine load, hydrogen fuel-air ratio and dopant chemical reactivity, but only modestly dependent on spark timing and nitrogen dilution levels below 20 percent. Fast FID measurements at the cylinder exit showed that the surviving hydrocarbons emerge late in the exhaust stroke. © Copyright 1996 Society of Automotive Engineers, Inc.

Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

Multi-objective optimization of a fluid structure interaction benchmarking

The integration and application of a new multi-objective tabu search optimization algorithm for Fluid Structure Interaction (FSI) problems are presented. The aim is to enhance the computational design process for real world applications and to achieve higher performance of the whole system for the four considered objectives. The described system combines the optimizer with a well established FSI solver which is based on the fully implicit, monolithic formuFlation of the problem in the Arbitrary Lagrangian-Eulerian FEM approach. The proposed solver resolves the proposed uid-structure interaction benchmark which describes the self-induced elastic deformation of a beam attached to a cylinder in laminar channel ow. The optimized ow characteristics of the aforementioned geometrical arrangement illustrate the performance of the system in two dimensions. Special emphasis is given to the analysis of the simulation package, which is of high accuracy and is the core of application. The design process identifies the best combination of ow features for optimal system behavior and the most important objectives. In addition, the presented methodology has the potential to run in parallel, which will significantly speed-up the elapsed time. Finite Element Method (FEM), Fluid-Structure Interaction (FSI), Multi-Ojective Tabu search (MOTS2). Copyright © 2013 Tech Science Press.

Absolute and convective instability in gas turbine fuel injectors

Hydrodynamic instabilities in gas turbine fuel injectors help to mix the fuel and air but can sometimes lock into acoustic oscillations and contribute to thermoacoustic instability. This paper describes a linear stability analysis that predicts the frequencies and strengths of hydrodynamic instabilities and identifies the regions of the flow that cause them. It distinguishes between convective instabilities, which grow in time but are convected away by the flow, and absolute instabilities, which grow in time without being convected away. Convectively unstable flows amplify external perturbations, while absolutely unstable flows also oscillate at intrinsic frequencies. As an input, this analysis requires velocity and density fields, either from a steady but unstable solution to the Navier-Stokes equations, or from time-averaged numerical simulations. In the former case, the analysis is a predictive tool. In the latter case, it is a diagnostic tool. This technique is applied to three flows: a swirling wake at Re = 400, a single stream swirling fuel injector at Re - 106, and a lean premixed gas turbine injector with five swirling streams at Re - 106. Its application to the swirling wake demonstrates that this technique can correctly predict the frequency, growth rate and dominant wavemaker region of the flow. It also shows that the zone of absolute instability found from the spatio-temporal analysis is a good approximation to the wavemaker region, which is found by overlapping the direct and adjoint global modes. This approximation is used in the other two flows because it is difficult to calculate their adjoint global modes. Its application to the single stream fuel injector demonstrates that it can identify the regions of the flow that are responsible for generating the hydrodynamic oscillations seen in LES and experimental data. The frequencies predicted by this technique are within a few percent of the measured frequencies. The technique also explains why these oscillations become weaker when a central jet is injected along the centreline. This is because the absolutely unstable region that causes the oscillations becomes convectively unstable. Its application to the lean premixed gas turbine injector reveals that several regions of the flow are hydrodynamically unstable, each with a different frequency and a different strength. For example, it reveals that the central region of confined swirling flow is strongly absolutely unstable and sets up a precessing vortex core, which is likely to aid mixing throughout the injector. It also reveals that the region between the second and third streams is slightly absolutely unstable at a frequency that is likely to coincide with acoustic modes within the combustion chamber. This technique, coupled with knowledge of the acoustic modes in a combustion chamber, is likely to be a useful design tool for the passive control of mixing and combustion instability. Copyright © 2012 by ASME.

Transient response of structures with uncertain properties to nonlinear shock loading

A method is presented to predict the transient response of a structure at the driving point following an impact or a shock loading. The displacement and the contact force are calculated solving the discrete convolution between the impulse response and the contact force itself, expressed in terms of a nonlinear Hertzian contact stiffness. Application of random point process theory allows the calculation of the impulse response function from knowledge of the modal density and the geometric characteristics of the structure only. The theory is applied to a wide range of structures and results are experimentally verified for the case of a rigid object hitting a beam, a plate, a thin and a thick cylinder and for the impact between two cylinders. The modal density of the flexural modes for a thick slender cylinder is derived analytically. Good agreement is found between experimental, simulated and published results, showing the reliability of the method for a wide range of situations including impacts and pyroshock applications. © 2013 Elsevier Ltd. All rights reserved.

Angle-of-attack effects on counter-rotating propellers at take-off

Due to their potential for significant fuel consumption savings, Counter-Rotating Open Rotors (CRORs) are currently being considered as an alternative to high-bypass turbofans. When CRORs are mounted on an aircraft, several 'installation effects' arise which are not present when the engine is operated in isolation. This paper investigates how flow features arising from one such effect - The angle-of-attack of the engine centre-line relative to the oncoming flow - can influence the design of CROR engines. Three-dimensional full-annulus unsteady CFD simulations are used to predict the time-varying flow field experienced by each rotor and emphasis is put on the interaction of the frontrotor wake and tip vortex with the rear-rotor. A parametric study is presented that quantifies the rotorrotor interaction as a function of the angle-of-attack. It is shown that angle-of-attack operation significantly changes the flow field and the unsteady lift on both rotors. In particular, a frequency analysis shows that the unsteady lift exhibits sidebands around the rotor-rotor interaction frequencies. Further, a non-linear increase in the total rear-rotor tip unsteadiness is observed for moderate and high angles-of-attack. The results presented in this paper demonstrate that common techniques used to mitigate CROR noise, such as modifying the rotor-rotor axial spacing and rear-rotor crop, can not be applied correctly unless angle-of-attack effects are taken into account. Copyright © 2012 by ASME.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Dynamic models of piled foundations

The vibration behavior of piled foundations is an important consideration in fields such as earthquake engineering, construction, machine-foundation design, offshore structures, nuclear energy, and road and rail development. This paper presents a review of the past 40 years' literature on modeling the frequency-dependent behavior of pile foundations. Beginning with the earliest model of a single pile, adapted from those for embedded footings, it charts the development of the four pile-modeling techniques: the "dynamic Winkler-foundation" approach that uses springs to represent the effect of the soil; elasticcontinuum-type formulations involving the analytical solutions for displacements due to a subsurface disk, cylinder, or other element; boundary element methods; and dynamic finite-element formulations with special nonreflecting boundaries. The modeling of pile groups involves accounting for pile-soil-pile interactions, and four such methods exist: interaction factors; complete pile models; the equivalent pier method; and periodic structure theory. Approaches for validating pile models are also explored. Copyright © 2013 by ASME.

On high-frequency noise scattering by aerofoils in flow

Abstract A theoretical model is developed for the sound scattered when a sound wave is incident on a cambered aerofoil at non-zero angle of attack. The model is based on the linearization of the Euler equations about a steady subsonic flow, and is an adaptation of previous work which considered incident vortical disturbances. Only high-frequency sound waves are considered. The aerofoil thickness, camber and angle of attack are restricted such that the steady flow past the aerofoil is a small perturbation to a uniform flow. The singular perturbation analysis identifies asymptotic regions around the aerofoil; local 'inner' regions, which scale on the incident wavelength, at the leading and trailing edges of the aerofoil; Fresnel regions emanating from the leading and trailing edges of the aerofoil due to the coalescence of singularities and points of stationary phase; a wake transition region downstream of the aerofoil leading and trailing edge; and an outer region far from the aerofoil and wake. An acoustic boundary layer on the aerofoil surface and within the transition region accounts for the effects of curvature. The final result is a uniformly-valid solution for the far-field sound; the effects of angle of attack, camber and thickness are investigated. © 2013 Cambridge University Press.

Large eddy simulations in low-pressure turbines: Effect of wakes at elevated free-stream turbulence

The transition of a separated shear layer over a flat plate, in the presence of periodic wakes and elevated free-stream turbulence (FST), is numerically investigated using Large Eddy Simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low-pressure turbine (LPT) blade. Two different distributions representative of a 'high-lift' and an 'ultra high-lift' turbine blade are examined. Results obtained from the current LES compare favourably with the extensive experimental data previously obtained for these configurations. The LES results are then used to further investigate the flow physics involved in the transition process.In line with experimental experience, the benefit of wakes and FST obtained by suppressing the separation bubble, is more pronounced in 'ultra high-lift' design when compared to the 'high-lift' design. Stronger 'Klebanoff streaks' are formed in the presence of wakes when compared to the streaks due to FST alone. These streaks promoted much early transition. The weak Klebanoff streaks due to FST continued to trigger transition in between the wake passing cycles.The experimental inference regarding the origin of Klebanoff streaks at the leading edge has been confirmed by the current simulations. While the wake convects at local free-stream velocity, its impression in the boundary layer in the form of streaks convects much slowly. The 'part-span' Kelvin-Helmholtz structures, which were observed in the experiments when the wake passes over the separation bubble, are also captured. The non-phase averaged space-time plots manifest that reattachment is a localized process across the span unlike the impression of global reattachment portrayed by phase averaging. © 2013 Elsevier Inc.