Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng-Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes. © 2013.

Experimental investigation of a control method for SI-HCCI-SI transition in a multi-cylinder gasoline engine

In HCCI engines, the Air/Fuel Ratio (AFR) and Residual Gas Fraction (RGF) are difficult to control during the SI-HCCI-SI transition, and this may result in incomplete combustion and/or high pressure raise rates. As a result, there may be undesirably high engine load fluctuations. The objectives of this work are to further understand this process and develop control methods to minimize these load fluctuations. This paper presents data on instantaneous AFR and RGF measurements, both taken by novel experimental techniques. The data provides an insight into the cyclic AFR and RGF fluctuations during the switch. These results suggest that the relatively slow change in the intake Manifold Air Pressure (MAP) and actuation time of the Variable Valve Timing (VVT) are the main causes of undesired AFR and RGF fluctuations, and hence an unacceptable Net IMEP (NIMEP) fluctuation. We also found large cylinder-to-cylinder AFR variations during the transition. Therefore, besides throttle opening control and VVT shifting, cyclic and individual cylinder fuel injection control is necessary to achieve a smooth transition. The control method was developed and implemented in a test engine, and the result was a considerably reduced NIMEP fluctuation during the mode switch. The instantaneous AFR and RGF measurements could furthermore be adopted to develop more sophisticated control methods for SI-HCCI-SI transitions. © 2010 SAE International.

Development and validation of a new method for measuring friction between skin and nonwoven materials.

A new method for measuring the coefficient of friction between nonwoven materials and the curved surface of the volar forearm has been developed and validated. The method was used to measure the coefficient of static friction for three different nonwoven materials on the normal (dry) and over-hydrated volar forearms of five female volunteers (ages 18-44). The method proved simple to run and had good repeatability: the coefficient of variation (standard deviation expressed as a percentage of the mean) for triplets of repeat measurements was usually (80 per cent of the time) less than 10 per cent. Measurements involving the geometrically simpler configuration of pulling a weighted fabric sample horizontally across a quasi-planar area of volar forearm skin proved experimentally more difficult and had poorer repeatability. However, correlations between values of coefficient of static friction derived using the two methods were good (R = 0.81 for normal (dry) skin, and 0.91 for over-hydrated skin). Measurements of the coefficient of static friction for the three nonwovens for normal (dry) and for over-hydrated skin varied in the ranges of about 0.3-0.5 and 0.9-1.3, respectively. In agreement with Amontons' law, coefficients of friction were invariant with normal pressure over the entire experimental range (0.1-8.2 kPa).

Blade row interaction in a high pressure steam turbine

This paper presents a study of the three-dimensional flow field within the blade rows of a high-pressure axial flow steam turbine stage. Compound lean angles have been employed to achieve relatively low blade loading for hub and tip section and so reduce the secondary losses. The flow field is investigated in a Low-Speed Research Turbine using pneumatic and hot-wire probes downstream of the blade row. Steady and unsteady numerical simulations were performed using structured 3D Navier-Stokes solver to further understand the flow field. Agreement between the simulations and the measurements has been found. The unsteady measurements indicate that there is a significant effect of the stator flow interaction in the downstream rotor blade. The transport of the stator viscous flow through the rotor blade row is described. Unsteady numerical simulations were found to be successful in predicting accurately the flow near the secondary flow interaction regions compared to steady simulations. A method to calculate the unsteady loss generated inside the blade row was developed from the steady numerical simulations. The contribution of various regions in the blade to the unsteady loss generation was evaluated. This method can assist the designer in identifying and optimizing the features of the flow that are responsible for the majority of the unsteady loss production. An analytical model was developed to quantify this effect for the vortex transport inside the downstream blade.

The prediction of hot streak migration in a high-pressure turbine

Accurate predictions of combustor hot streak migration enable the turbine designer to identify high-temperature regions that can limit component life. It is therefore important that these predictions are achieved within the short time scales of a design process. This article compares temperature measurements of a circular hot streak through a turning duct and a research turbine with predictions using a three-dimensional Reynolds-averaged Navier-Stokes solver. It was found that the mixing length turbulence model did not predict the hot streak dissipation accurately. However, implementation of a very simple model of the free stream turbulence (FST) significantly improved the exit temperature predictions on both the duct and research turbine. One advantage of the simple FST model described over more complex alternatives is that no additional equations are solved. This makes the method attractive for design purposes, as it is not associated with any increase in computational time.

Immersed Boundary Method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples. Copyright © 2010 by ASME.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

A novel method to produce small droplets from large nozzles.

This work presents a new method to generate droplets with diameters significantly smaller than the nozzle from which they emerge. The electrical waveform used to produce the jetting consists of a single square negative pulse. The negative edge of the pressure wave pulls the meniscus in, overturning the surface in such a way that a cavity is created. This cavity is then forced to collapse under the action of the positive edge of the pressure wave. This violent collapse produces a thin jet that eventually breaks up and produces droplets. Four droplet generator prototypes that demonstrate the capabilities of this novel mechanism are described. It is also shown that the proposed mechanism extends the existing limits of the commonly accepted inkjet operating regime.

Blade Loading and Its Application in the Mean-Line Design of Low Pressure Turbines

In order to minimize the number of iterations to a turbine design, reasonable choices of the key parameters must be made at the preliminary design stage. The choice of blade loading is of particular concern in the low pressure (LP) turbine of civil aero engines, where the use of high-lift blades is widespread. This paper considers how blade loading should be measured, compares the performance of various loss correlations, and explores the impact of blade lift on performance and lapse rates. To these ends, an analytical design study is presented for a repeating-stage, axial-flow LP turbine. It is demonstrated that the long-established Zweifel lift coefficient (Zweifel, 1945, "The Spacing of Turbomachine Blading, Especially with Large Angular Deflection" Brown Boveri Rev., 32(1), pp. 436-444) is flawed because it does not account for the blade camber. As a result the Zweifel coefficient is only meaningful for a fixed set of flow angles and cannot be used as an absolute measure of blade loading. A lift coefficient based on circulation is instead proposed that accounts for the blade curvature and is independent of the flow angles. Various existing profile and secondary loss correlations are examined for their suitability to preliminary design. A largely qualitative comparison demonstrates that the loss correlations based on Ainley and Mathieson (Ainley and Mathieson, 1957, "A Method of Performance Estimation for Axial-Flow Turbines," ARC Reports and Memoranda No. 2974; Dunham and Came, 1970, "Improvements to the Ainley-Mathieson Method of Turbine Performance Prediction," Trans. ASME: J. Eng. Gas Turbines Power, July, pp. 252-256; Kacker and Okapuu, 1982, "A Mean Line Performance Method for Axial Flow Turbine Efficiency," J. Eng. Power, 104, pp. 111-119). are not realistic, while the profile loss model of Coull and Hodson (Coull and Hodson, 2011, "Predicting the Profile Loss of High-Lift Low Pressure Turbines," J. Turbomach., 134(2), pp. 021002) and the secondary loss model of (Traupel, W, 1977, Thermische Turbomaschinen, Springer-Verlag, Berlin) are arguably the most reasonable. A quantitative comparison with multistage rig data indicates that, together, these methods over-predict lapse rates by around 30%, highlighting the need for improved loss models and a better understanding of the multistage environment. By examining the influence of blade lift across the Smith efficiency chart, the analysis demonstrates that designs with higher flow turning will tend to be less sensitive to increases in blade loading. © 2013 American Society of Mechanical Engineers.