Nicholas J. Hills

Bio: Nicholas J. Hills is an academic researcher from University of Surrey. The author has contributed to research in topics: Turbine & Computational fluid dynamics. The author has an hindex of 18, co-authored 76 publications receiving 882 citations. Previous affiliations of Nicholas J. Hills include Information Technology University & University of Sussex.

Measurement and Analysis of Ingestion Through a Turbine Rim Seal

Abstract: Experimental measurements from a new single stage turbine are presented. The turbine has 26 vanes and 59 rotating blades with a design point stage expansion ratio of 2.5 and vane exit Mach number of 0.96. A variable sealing flow is supplied to the disc cavity upstream of the rotor and then enters the annulus through a simple axial clearance seal situated on the hub between the stator and rotor. Measurements at the annulus hub wall just downstream of the vanes show the degree of circumferential pressure variation. Further pressure measurements in the disc cavity indicate the strength of the swirling flow in the cavity, and show the effects of mainstream gas ingestion at low sealing flows. Ingestion is further quantified through seeding of the sealing air with nitrous oxide or carbon dioxide and measurement of gas concentrations in the cavity. Interpretation of the measurements is aided by steady and unsteady computational fluid dynamics solutions, and comparison with an elementary model of ingestion.

Computational and Mathematical Modelling of Turbine Rim Seal Ingestion

Abstract: Understanding and modelling of main annulus gas ingestion through turbine rim seals is considered and advanced in this paper. Unsteady 3-dimensional computational fluid dynamics (CFD) calculations and results from a more elementary model are presented and compared with experimental data previously published by Hills et al (1997). The most complete CFD model presented includes both stator and rotor in the main annulus and the inter-disc cavity. The k-e model of turbulence with standard wall function approximations is assumed in the model which was constructed in a commercial CFD code employing a pressure correction solution algorithm. It is shown that considerable care is needed to ensure convergence of the CFD model to a periodic solution. Compared to previous models, results from the CFD model show encouraging agreement with pressure and gas concentration measurements. The annulus gas ingestion is shown to result from a combination of the stationary and rotating circumferential pressure asymmetries in the annulus. Inertial effects associated with the circumferential velocity component of the flow have an important effect on the degree of ingestion. The elementary model used is an extension of earlier models based on orifice theory applied locally around the rim seal circumference. The new model includes a term accounting for inertial effects. Some good qualitative and fair quantitative agreement with data is shown. Copyright © 2001 by ASME.

Efficient Finite Element Analysis/Computational Fluid Dynamics Thermal Coupling for Engineering Applications

Abstract: An efficient finite element analysis/computational fluid dynamics (FEA/CFD) thermal coupling technique has been developed and demonstrated. The thermal coupling is achieved by an iterative procedure between FEA and CFD calculations. Communication between FEA and CFD calculations ensures continuity of temperature and heat flux. In the procedure, the FEA simulation is treated as unsteady for a given transient cycle. To speed up the thermal coupling, steady CFD calculations are employed, considering that fluid flow time scales are much shorter than those for the solid heat conduction and therefore the influence of unsteadiness in fluid regions is negligible. To facilitate the thermal coupling, the procedure is designed to allow a set of CFD models to be defined at key time points/intervals in the transient cycle and to be invoked during the coupling process at specified time points. To further enhance computational efficiency, a “frozen flow” or “energy equation only” coupling option was also developed, where only the energy equation is solved, while the flow is frozen in CFD simulation during the thermal coupling process for specified time intervals. This option has proven very useful in practice, as the flow is found to be unaffected by the thermal boundary conditions over certain time intervals. The FEA solver employed is an in-house code, and the coupling has been implemented for two different CFD solvers: a commercial code and an in-house code. Test cases include an industrial low pressure (LP) turbine and a high pressure (HP) compressor, with CFD modeling of the LP turbine disk cavity and the HP compressor drive cone cavity flows, respectively. Good agreement of wall temperatures with the industrial rig test data was observed. It is shown that the coupled solutions can be obtained in sufficiently short turn-around times (typically within a week) for use in design.

Large-Eddy simulation of rim seal ingestion:

Abstract: Unsteady flow dynamics in turbine rim seals are known to be complex and attempts accurately to predict the interaction of the mainstream flow with the secondary air system cooling flows using computational fluid dynamics (CFD) with Reynolds-averaged Navier–Stokes (RANS) turbulence models have proved difficult. In particular, published results from RANS models have over-predicted the sealing effectiveness of the rim seal, although their use in this context continues to be common. Previous studies have ascribed this discrepancy to the failure to model flow structures with a scale greater than the one which can be captured in the small-sector models typically used. This article presents results from a series of Large-Eddy Simulations (LES) of a turbine stage including a rim seal and rim cavity for, it is believed by the authors, the first time. The simulations were run at a rotational Reynolds number Re ¼ 2.2 106 and a main annulus axial Reynolds number Rex ¼ 1.3 106 and with varying levels of coolant mass flow. Comparison is made with previously published experimental data and with unsteady RANS simulations. The LES models are shown to be in closer agreement with the experimental sealing effectiveness than the unsteady RANS simulations. The result indicates that the previous failure to predict rim seal effectiveness was due to turbulence model limitations in the turbine rim seal flow. Consideration is given to the flow structure in this region. K

Coupled fluid/solid heat transfer computation for turbine discs

Abstract: This paper considers the coupling of a finite element thermal conduction solver with a steady, finite volume fluid flow solver. Two methods were considered for passing boundary conditions between the two codes - transfer of metal temperatures and either convective heat fluxes or heat transfer coefficients and air temperatures. These methods have been tested on two simple rotating cavity test cases and also on a more complex real engine example. Convergence rates of the two coupling methods were compared. Passing heat transfer coefficients and air temperatures was found to give the quickest convergence. The coupled method gave agreement with the analytic solution and a conjugate solution of the simple free disc problem. The predicted heat transfer results for the real engine example showed some encouraging agreement, although some modelling issues are identified. Copyright © 2001 by ASME.

New insights into large eddy simulation

Abstract: Fluid dynamic turbulence is one of the most challenging computational physics problems because of the extremely wide range of time and space scales involved, the strong nonlinearity of the governing equations, and the many practical and important applications. While most linear fluid instabilities are well understood, the nonlinear interactions among them makes even the relatively simple limit of homogeneous isotropic turbulence difficult to treat physically, mathematically, and computationally. Turbulence is modeled computationally by a two-stage bootstrap process. The first stage, direct numerical simulation, attempts to resolve the relevant physical time and space scales but its application is limited to diffusive flows with a relatively small Reynolds number (Re). Using direct numerical simulation to provide a database, in turn, allows calibration of phenomenological turbulence models for engineering applications. Large eddy simulation incorporates a form of turbulence modeling applicable when the large-scale flows of interest are intrinsically time dependent, thus throwing common statistical models into question. A promising approach to large eddy simulation involves the use of high-resolution monotone computational fluid dynamics algorithms such as flux-corrected transport or the piecewise parabolic method which have intrinsic subgrid turbulence models coupled naturally to the resolved scales in the computed flow. The physical considerations underlying and evidence supporting this monotone integrated large eddy simulation approach are discussed.

Effects of the computational time step on numerical solutions for turbulent flow

Abstract: Effects of large computational time steps on the computed turbulence were investigated using a fully implicit method. In turbulent channel flow computations the largest computational time step in wall units which led to accurate prediction of turbulence statistics was determined. Turbulence fluctuations could not be sustained if the computational time step was near or larger than the Kolmogorov time scale.