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.