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.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.

A new Lagrangian-Eulerian shell-fluid coupling algorithm based on level sets

We propose a computational method for the coupled simulation of a compressible flow interacting with a thin-shell structure undergoing large deformations. An Eulerian finite volume formulation is adopted for the fluid and a Lagrangian formulation based on subdivision finite elements is adopted for the shell response. The coupling between the fluid and the solid response is achieved via a novel approach based on level sets. The basic approach furnishes a general algorithm for coupling Lagrangian shell solvers with Cartesian grid based Eulerian fluid solvers. The efficiency and robustness of the proposed approach is demonstrated with a airbag deployment simulation. It bears emphasis that in the proposed approach the solid and the fluid components as well as their coupled interaction are considered in full detail and modeled with an equivalent level of fidelity without any oversimplifying assumptions or bias towards a particular physical aspect of the problem.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

On separating propagating and non-propagating dynamics in fluid-flow equations

The ability to separate acoustically radiating and non-radiating components in fluid flow is desirable to identify the true sources of aerodynamic sound, which can be expressed in terms of the non-radiating flow dynamics. These non-radiating components are obtained by filtering the flow field. Two linear filtering strategies are investigated: one is based on a differential operator, the other employs convolution operations. Convolution filters are found to be superior at separating radiating and non-radiating components. Their ability to decompose the flow into non-radiating and radiating components is demonstrated on two different flows: one satisfying the linearized Euler and the other the Navier-Stokes equations. In the latter case, the corresponding sound sources are computed. These sources provide good insight into the sound generation process. For source localization, they are found to be superior to the commonly used sound sources computed using the steady part of the flow. Copyright © 2009 by S. Sinayoko, A. Agarwal, Z. Hu.

Fluid mechanics of pulse pressure-gain combustors

The fluid dynamic operation of a valveless pulse combustor has been studied experimentally and numerically. Through phase-locked chemiluminescence and pressure measurements it is shown that mechanical energy is created periodically in the flame surface, with an efficiency of 1.6%. This mechanical energy leaves the pulse combustor through unsteady jets at the aerovalve inlet and the tailpipe exit stations. Two thermodynamically distinct flows are identified: a flow that is transported from inlet to exit and participates in combustion along the way, and a flow that is ingested and then ejected from the combustor without undergoing combustion. It is the latter of these two flows which has the greatest quantity of net work done on it. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

Fluid-structure interaction simulation of an inflatable aerodynamic tension-cone decelerator

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. In this paper, we focus on an inflatable tension cone design that has potential advantages over other geometries. A computational fluid-structure interaction model of a tension cone is employed to investigate the behavior of the inflatable aeroshell at supersonic speeds for conditions matching recent experimental results. A parametric study is carried out to investigate the deflections of the tension cone as a function of inflation pressure of the torus at a Mach of 2.5. Comparison of the behavior of the structure, amplitude of deformations, and determined loads are reported. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Vibration of fluid loaded conical shells.

An analytical model is presented to describe the vibration of a truncated conical shell with fluid loading in the low frequency range. The solution for the dynamic response of the shell is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips which are considered to be locally cylindrical. Analytical results are presented for different boundary conditions and have been compared with the computational results from a boundary element model. Limitations of the model to the low frequency range are discussed.