An accelerated 3D Navier-Stokes solver for flows in turbomachines

A new three-dimensional Navier-Stokes solver for flows in turbomachines has been developed. The new solver is based on the latest version of the Denton codes, but has been implemented to run on Graphics Processing Units (GPUs) instead of the traditional Central Processing Unit (CPU). The change in processor enables an order-of-magnitude reduction in run-time due to the higher performance of the GPU. Scaling results for a 16 node GPU cluster are also presented, showing almost linear scaling for typical turbomachinery cases. For validation purposes, a test case consisting of a three-stage turbine with complete hub and casing leakage paths is described. Good agreement is obtained with previously published experimental results. The simulation runs in less than 10 minutes on a cluster with four GPUs. Copyright © 2009 by ASME.

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.

Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

Influence of surfactant solubility on the deformation and breakup of a bubble or capillary jet in a viscous fluid

In a previous study [M. Hameed, J. Fluid Mech. 594, 307 (2008)] the authors investigated the influence of insoluble surfactant on the evolution of a stretched, inviscid bubble surrounded by a viscous fluid via direct numerical simulation of the Navier-Stokes equations, and showed that the presence of surfactant can cause the bubble to contract and form a quasisteady slender thread connecting parent bubbles, instead of proceeding directly toward pinch-off as occurs for a surfactant-free bubble. Insoluble surfactant significantly retards pinch-off and the thread is stabilized by a balance between internal pressure and reduced capillary pressure due to a high concentration of surfactant that develops during the initial stage of contraction. In the present study we investigate the influence of surfactant solubility on thread formation. The adsorption-desorption kinetics for solubility is in the diffusion controlled regime. A long-wave model for the evolution of a capillary jet is also studied in the Stokes flow limit, and shows dynamics that are similar to those of the evolving bubble. With soluble surfactant, depending on parameter values, a slender thread forms but can pinch-off later due to exchange of surfactant between the interface and exterior bulk flow. © 2009 American Institute of Physics.

MEASURED AND PREDICTED LOSS GENERATION IN TRANSONIC TURBINE BLADING.

For a typical transonic turbine rotor blade, designed for use with coolant ejection, the trailing edge, or base loss is three to four times the profile boundary layer loss. The base region of such a profile is dominated by viscous effects and it seems essential to attack the problem of loss prediction by solving the compressible Navier-Stokes equations. However, such an approach is inevitably compromised by both numerical accuracy and turbulence modelling constraints. This paper describes a Navier-Stokes solver written for 2D blade-blade flows and employing a simple two-layer mixing length eddy viscosity model. Then, measured and predicted losses and base pressures are presented for two transonic rotor blades and attempts are made to assess the capabilities of the Navier-Stokes solver and to outline areas for future work.

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

Investigation of air-core vortex at hydraulic intakes

Hydrodynamic properties of the surface vortex have been investigated. Based on the Navier-Stokes equations, three sets of the new formulations for the tangential velocity distributions are derived, and verified against the experimental measurements in the literature. It is shown that one modification greatly improves the agreement with the experimental data. Physical model experiments were carried out to study the intake vortex related to the Xiluodu hydropower project. The velocity fields were measured using the Particle Tracking Velocimetry (PTV) technique. The proposed equation for tangential velocity distribution is applied to the Xiluodu project with the solid boundary being considered by the method of images. Good agreement has been observed between the formula prediction and the experimental observation. © 2010 Publishing House for Journal of Hydrodynamics.

On computing the physical sources of jet noise

We derive a closed system of equations that relates the acoustically radiating flow variables to the sources of sound for homentropic flows. We use radiating density, momentum density and modified pressure as the dependent variables which leads to simple source terms for the momentum equations. The source terms involve the non-radiating parts of the density and momentum density fields. These non-radiating components are obtained by removing the radiating wavenumbers in the Fourier domain. We demonstrate the usefulness of this new technique on an axi-symmetric jet solution of the Navier-Stokes equations, obtained by direct numerical simulation (DNS). The dominant source term is proportional to the square of the non-radiating part of the axial momentum density. We compare the sound sources to that obtained by an acoustic analogy and find that they have more realistic physical properties. Their frequency content and amplitudes are consistent with. We validate the sources by computing the radiating sound field and comparing it to the DNS solution. © 2010 by S. Sinayoko, A. Agarwal.

Comparison between empirical formulae of intake vortices

The hydrodynamic properties of free surface vortices at hydraulic intakes were investigated. Based on the axisymmetric Navier-Stokes equations and empirical assumptions, two sets of formulations for the velocity distributions and the free surface profiles are proposed and validated against measurements available in the literature. Compared with previous formulae, the modifications based on Mih's formula are found to greatly improve the agreement with the experimental data. Physical model tests were also conducted to study the intake vortex of the Xiluodu hydroelectric project in China. The proposed velocity distribution formula was applied to the solid boundary as considered by the method of images. A good agreement was again observed between the prediction and the measurements. © 2011 International Association for Hydro-Environment Engineering and Research.

Reprint of: Flow filtering and the physical sources of aerodynamic sound

The physical sources of sound are expressed in terms of the non-radiating part of the flow. The non-radiating part of the flow can be obtained from convolution filtering, as we demonstrate numerically by using an axi-symmetric jet satisfying the Navier-Stokes equations. Based on the frequency spectrum of the source, we show that the sound sources exhibit more physical behaviour than sound sources based on acoustic analogies. To validate the sources of sound, one needs to let them radiate within the non-radiating flow field. However, our results suggest that the traditional Euler operator linearized about the time-averaged part of the flow should be sufficient to compute the sound field. © 2010 Published by Elsevier Ltd.

Magnetically actuated artificial cilia: the effect of fluid inertia.

Natural cilia are hairlike microtubule-based structures that are able to move fluid on the micrometer scale using asymmetric motion. In this article, we follow a biomimetic approach to design artificial cilia lining the inner surfaces of microfluidic channels with the goal of propelling fluid. The artificial cilia consist of polymer films filled with superparamagnetic nanoparticles, which can mimic the motion of natural cilia when subjected to a rotating magnetic field. To obtain the magnetic field and associated magnetization local to the cilia, we solve the Maxwell equations, from which the magnetic body moments and forces can be deduced. To obtain the ciliary motion, we solve the dynamic equations of motion, which are then fully coupled to the Navier-Stokes equations that describe the fluid flow around the cilia, thus taking full account of fluid inertial forces. The dimensionless parameters that govern the deformation behavior of the cilia and the associated fluid flow are arrived at using the principle of virtual work. The physical response of the cilia and the fluid flow for different combinations of elastic, fluid viscous, and inertia forces are identified.

Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.