Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts.

This paper reviews the development of computational fluid dynamics (CFD) specifically for turbomachinery simulations and with a particular focus on application to problems with complex geometry. The review is structured by considering this development as a series of paradigm shifts, followed by asymptotes. The original S1-S2 blade-blade-throughflow model is briefly described, followed by the development of two-dimensional then three-dimensional blade-blade analysis. This in turn evolved from inviscid to viscous analysis and then from steady to unsteady flow simulations. This development trajectory led over a surprisingly small number of years to an accepted approach-a 'CFD orthodoxy'. A very important current area of intense interest and activity in turbomachinery simulation is in accounting for real geometry effects, not just in the secondary air and turbine cooling systems but also associated with the primary path. The requirements here are threefold: capturing and representing these geometries in a computer model; making rapid design changes to these complex geometries; and managing the very large associated computational models on PC clusters. Accordingly, the challenges in the application of the current CFD orthodoxy to complex geometries are described in some detail. The main aim of this paper is to argue that the current CFD orthodoxy is on a new asymptote and is not in fact suited for application to complex geometries and that a paradigm shift must be sought. In particular, the new paradigm must be geometry centric and inherently parallel without serial bottlenecks. The main contribution of this paper is to describe such a potential paradigm shift, inspired by the animation industry, based on a fundamental shift in perspective from explicit to implicit geometry and then illustrate this with a number of applications to turbomachinery.

Flow of wet powder in a conical centrifugal filter-an analytical model

A one-dimensional analytical model is developed for the steady state, axisymmetric, slender flow of saturated powder in a rotating perforated cone. Both the powder and the fluid spin with the cone with negligible slip in the hoop direction. They migrate up the wall of the cone along a generator under centrifugal force, which also forces the fluid out of the cone through the powder layer and the porous wall. The flow thus evolves from an over-saturated paste at inlet into a nearly dry powder at outlet. The powder is treated as a Mohr-Coulomb granular solid of constant void fraction and permeability. The shear traction at the wall is assumed to be velocity and pressure dependent. The fluid is treated as Newtonian viscous. The model provides the position of the colour line (the transition from over- to under-saturation) and the flow velocity and thickness profiles over the cone. Surface tension effects are assumed negligible compared to the centrifugal acceleration. Two alternative conditions are considered for the flow structure at inlet: fully settled powder at inlet, and progressive settling of an initially homogeneous slurry. The position of the colour line is found to be similar for these two cases over a wide range of operating conditions. Dominant dimensionless groups are identified which control the position of the colour line in a continuous conical centrifuge. Experimental observations of centrifuges used in the sugar industry provide preliminary validation of the model. © 2011 Elsevier Ltd.

Combustion rumble prediction with integrated computational-fluid-dynamics/ low-order-model methods

The pressure oscillation within combustion chambers of aeroengines and industrial gas turbines is a major technical challenge to the development of high-performance and low-emission propulsion systems. In this paper, an approach integrating computational fluid dynamics and one-dimensional linear stability analysis is developed to predict the modes of oscillation in a combustor and their frequencies and growth rates. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone, which enables the computational fluid dynamics calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in agreement with experimental measurements. Furthermore, results from the computational fluid dynamics calculation provide the flame transfer function to describe unsteady heat release rate. Departures from ideal one-dimensional flows are described by shape factors. Combined with this information, low-order models can work out the possible oscillation modes and their initial growth rates. The approach developed here can be used in more general situations for the analysis of combustion oscillations. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

The night purging of a two-storey atrium building

We examine the fluid mechanics of night purging in a two-storey naturally ventilated atrium building. We develop a mathematical model of a simplified atrium building and focus on the rate at which warm air purges from each storey and the atrium by displacement ventilation into a still cool night environment of a constant temperature. To develop a first insight into how the geometry of the building influences the rate at which warm air purges from each storey via the atrium we neglect heat exchange with the fabric (so there is no thermal buffering) and furthermore assume that the warm air layers in each storey and the atrium are of uniform temperature. The plumes of warm air that rise from the storeys into the atrium, causing the atrium to fill with warm air, have a very strong influence on the night purge. Modelling these as axisymmetric turbulent plumes, we identify three forms of purging behaviour. Each purge is characterised by five key times identified in the progression of the night purge and physical rationale for these differing behaviours is given. An interface velocity deficit and volumetric purge deficit are introduced as measures of the efficiency of a night purge. © 2010 Elsevier Ltd.

On the two classes of global primary modal instability in laminar separation bubbles

The self-excited global instability mechanisms existing in flat-plate laminar separation bubbles are studied here, in order to shed light on the causes of unsteadiness and three- dimensionality of unforced, nominally two-dimensional separated flows. The presence of two known linear global mechanisms, namely an oscillator behavior driven by local regions of absolute inflectional instability and a centrifugal instability giving rise to a steady three- dimensionalization of the bubble, is studied in a series of model separation bubbles. Present results indicate that absolute instability, and consequently a global oscillator behavior, does not exist for two-dimensional bubbles with a peak reversed-flow velocity below 12% of the free-stream velocity. However, the three-dimensional instability becomes active for recirculation levels as low as urev ≈ 7%. These findings suggest a route to the three-dimensionality and unsteadiness observed in experiments and simulations substantially different from that usually found in the literature, in which two-dimensional vortex shedding is followed by three-dimensionalization.

Properties of two model soils stabilized with different blends and contents of GGBS, MgO, lime, and PC

This paper addresses the use of ground granulated blast furnace slag (GGBS) and reactive magnesia (MgO) blends for soil stabilization, comparing them with GGBS-lime blends and Portland cement (PC) for enhanced technical performance. A range of tests were conducted to investigate the properties of stabilized soils, including unconfined compressive strength (UCS), permeability, and microstructural analyses by using X-ray diffraction (XRD) and scanning electron microscopy (SEM). The influence of GGBS:MgO ratio, binder content, soil type, and curing period were addressed. The UCS results revealed that GGBS-MgO was more efficient than GGBS-lime as a binder for soil stabilization, with an optimum MgO content in the range of 5-20% of the blends content, varying with binder content and curing age. The 28-day UCS values of the optimum GGBS-MgO mixes were up to almost four times higher than that of corresponding PC mixes. The microstructural analyses showed the hydrotalcite was produced during the GGBS hydration activated by MgO, although the main hydration products of the GGBS-MgO stabilized soils were similar to those of PC. © 2014 American Society of Civil Engineers.