Showing papers in "Annual Review of Fluid Mechanics in 1992"

Wavelet Transforms and their Applications to Turbulence

Abstract: Wavelet transforms are recent mathematical techniques, based on group theory and square integrable representations, which allows one to unfold a signal, or a field, into both space and scale, and possibly directions. They use analyzing functions, called wavelets, which are localized in space. The scale decomposition is obtained by dilating or contracting the chosen analyzing wavelet before convolving it with the signal. The limited spatial support of wavelets is important because then the behavior of the signal at infinity does not play any role. Therefore the wavelet analysis or syn­ thesis can be performed locally on the signal, as opposed to the Fourier transform which is inherently nonlocal due to the space-filling nature of the trigonometric functions. Wavelet transforms have been applied mostly to signal processing, image coding, and numerical analysis, and they are still evolving. So far there are only two complete presentations of this topic, both written in French, one for engineers (Gasquet & Witomski 1 990) and the other for mathematicians (Meyer 1 990a), and two conference proceedings, the first in English (Combes et al 1 989), the second in French (Lemarie 1 990a). In preparation are a textbook (Holschneider 199 1 ), a course (Dau­ bee hies 1 99 1), three conference procecdings (Mcyer & Paul 199 1 , Beylkin et al 199 1b, Farge et al 1 99 1), and a special issue of IEEE Transactions

Hydrodynamic Phenomena in Suspensions of Swimming Microorganisms

Abstract: Vast numbers of microorganisms are suspended in temperate aqueous environments. Oceans and rivers, puddles and droplets, the fluid interiors of animals, all host an array of splendidly varied creatures. Although their presence is usually not casua1\y obvious, they constitute the major part of the world's biomass. Their population dynamics-replication and decay, accumulation and dispersal-modulates and regulates their own life, the life of the larger creatures that feed on them, and even the climate (Charlson et al 1987). Microorganisms interact with each other and with the world, at length scales that vary upward from the size of an individual, say 10-4 cm, to the dimensions of the entire body of fluid in which they live. Fluid mechanics, in concert with the organisms' behavior, governs the dynamics of many of these interactions. However, we are concerned here with relatively localized, small-scale phenomena. We consider only single­ celled microorganisms which are motile, i .e. self-propelled, and so small that inertial effects can be ignored in describing their locomotion. Examples to be discussed include species of algae, bacteria, and protozoa, not to mention spermatozoa. The trajectories along which individual cells swim are determined by thc

Helicity in Laminar and Turbulent Flow

Abstract: The quantity h(x, t) = u" o~ is the helicity density of the flow. Both h and ~ are pseudoscalar quantities, i.e. they change sign under change from a right-handed to a left-handed frame of reference (parity transformation). It is important therefore to specify the frame that is used; we shall always use a right-handed Cartesian (or orthogonal curvilinear) frame unless otherwise stated. The simplest (prototype) helical flow u = U+1⁄2f~ x where U and f~ are constants. Then curl

Modeling of two-phase slug flow

Abstract: Slug flow belongs to a class of intermittent flows that has very distinctive features Even if the inlet conditions are stationary, the flow as seen by an observer, is an unsteady phenomenon, dispersed flow appearing alternately with separated flow These two states follow in a random-like manner, inducing pressure and velocity fluctuations In vertical flow, the large bUbbles-typically longer than, say, one pipe diameter-rise with a round shaped front followed by a cylindrical main body surrounded by an annu­ lar liquid film In the literature, these long bubbles are frequently referred to as Taylor bubbles or Dumitrescu bubbles In the film, gravity forces the

Finite element methods for navier-stokes equations

Abstract: The main goal of this article is to address the finite element solution of the Navier-Stokes equations modeling compressible and incompressible viscous flow. It is the opinion of the authors that most general and reliable incompressible viscous flow simulators arc based on finite clement method­ ologies; these simulators, which can handle complicated geometries and boundary conditions, free surfaces, and turbulence effects, are well suited to industrial applications. On the other hand, the situation is much less satisfying concerning compressible viscous flow simulation, particularly at high Reynolds and Mach numbers and much progress must still be made in order to reach the degree of achievement obtained by the incompressible flow simulations. In this article, whose scope has been voluntarily limited, we concentrate on those topics combining our own work and the work of others with which we are familiar. The paper has been divided into two parts: In the first part (Sections I to 5) we discuss the various ingredients of a solution methodology for incompressible viscous flow based on operator splitting. Via splitting one obtains, at each time step, two families of subproblems: one of advection-diffusion type and one related to the steady Stokes

Contour Dynamics Methods

Abstract: In an early paper on the stability of fluid layers with uniform vorticity Rayleigh (1880) remarks: "... In such cases, the velocity curve is composed of portions of straight lines which meet each other at finite angles. This state of things may be supposed to be slightly disturbed by bending the surfaces of transition, and the determination of the subsequent motion depends upon that of the form of these surfaces. For co retains its constant value throughout each layer unchanged in the absence of friction, and by a well-known theorem the whole motion depends upon [omega]." We can now recognize this as essentially the principal of contour dynamics (CD), where [omega] is the uniform vorticity. The theorem referred to is the Biot-Savart law. Nearly a century later Zabusky et al (1979) presented numerical CD calculations of nonlinear vortex patch evolution. Subsequently, owing to its compact form conferring a deceptive simplicity, CD has become a widely used method for the investigation of two-dimensional rotational flow of an incompressible inviscid fluid. The aim of this article is to survey the development, technical details, and vortex-dynamic applications of the CD method in an effort to assess its impact on our understanding of the mechanics of rotational flow in two dimensions at ultrahigh Reynolds numbers. The study of the dynamics of two- and three-dimensional vortex mechanics by computational methods has been an active research area for more than two decades. Quite apart from many practical applications in the aerodynamics of separated flows, the theoretical and numerical study of vortices in incompressible fluids has been stimulated by the idea that turbulent fluid motion may be viewed as comprising ensembles of more or less coherent laminar vortex structures that interact via relatively simple dynamics and by the appeal of the vorticity equation, which does not contain the fluid pressure. Two-dimensional vortex interactions have been perceived as supposedly relevant to the origins of coherent structures observed experimentally in mixing layers, jets, and wakes, and for models of large-scale atmospheric and oceanic turbulence. Interest has often focused on the limit of infinite Reynolds number, where in the absence of boundaries, the inviscid Euler equations are assumed to properly describe the flow dynamics. The numerous surveys of progress in the study of vorticity and the use of numerical methods applied to vortex mechanics include articles by Saffman & Baker (1979) and Saffman (1981) on inviscid vortex interactions and Aref (1983) on two-dimensional flows. Numerical methods have been surveyed by Chorin (1980), and Leonard (1980, 1985). Caflisch (1988) describes work on the mathematical aspects of the subject. Zabusky (1981), Aref (1983), and Melander et al (1987b) discuss various aspects of CD. The review of Dritschel (1989) gives emphasis to numerical issues in CD and to recent computations with contour surgery. This article is confined to a discussion of vortices on a two-dimensional surface. We generally follow Saffman & Baker (1979) in matters of definition. In two dimensions a vortex sheet is a line of discontinuity in velocity while a vortex jump is a line of discontinuity in vorticity. We shall, however, use filament to denote a two-dimensional ribbon of vorticity surrounded by fluid with vorticity of different magnitude (which may be zero), rather than the more usual three-dimensional idea of a vortex tube. The ambiguity is unfortunate but is already in the literature. Additionally, a vortex patch is a finite, singly connected area of uniform vorticity while a vortex strip is an infinite strip of uniform vorticity with finite thickness, or equivalently, an infinite filament. Contour Dynamics will refer to the numerical solution of initial value problems for piecewise constant vorticity distributions by the Lagrangian method of calculating the evolution of the vorticity jumps. Such flows are often related to corresponding solutions of the Euler equations that are steady in some translating or rotating frame of reference. These solutions will be called vortex equilibria, and the numerical technique for computing their shapes based on CD is often referred to as contour statics. The mathematical foundation for the study of vorticity was laid primarily by the well-known investigations of Helmholtz, Kelvin, J. J. Thomson, Love, and others. In our century, efforts to produce numerical simulations of flows governed by the Euler equations have utilized a variety of Eulerian, Lagrangian, and hybrid methods. Among the former are the class of spectral methods that now comprise the prevailing tool for large-scale two- and three-dimensional calculations (see Hussaini & Zang 1987). The Lagrangian methods for two-dimensional flows have been predominantly vortex tracking techniques based on the Helmholtz vorticity laws. The first initial value calculations were those of Rosenhead (193l) and Westwater (1935) who attempted to calculate vortex sheet evolution by the motion of O(10) point vortices. Subsequent efforts by Moore (1974) (see also Moore 1983, 1985) and others to produce more refined computations for vortex sheets have failed for reasons related to the tendency for initially smooth vortex sheet data to produce singularities (Moore 1979). Discrete vortex methods used to study the nonlinear dynamics of vortex patches and layers have included the evolution of assemblies of point vortices by direct summation (e.g. Acton 1976) and the cloud in cell method (Roberts & Christiansen 1972, Christiansen & Zabusky 1973, Aref & Siggia 1980, 1981). For reviews see Leonard (1980) and Aref (1983). These techniques have often been criticized for their lack of accuracy and numerical convergence and because they may be subject to grid scale dispersion. However, many qualitative vortex phenomena observed in nature and in experiments, such as amalgamation events and others still under active investigation (e.g. filamentation) were first simulated numerically with discrete vortices. The contour dynamics approach is attractive because it appears to allow direct access, at least for small times, to the inviscid dynamics for vorticity distributions smoother than those of either point vortices or vortex sheets, while at the same time enabling the mapping of the two-dimensional Euler equations to a one-dimensional Lagrangian form. In Section 2 we discuss the formulation and numerical implementation of contour dynamics for the Euler equations in two dimensions. Section 3 is concerned with applications to isolated and multiple vortex systems and to vortex layers. An attempt is made to relate this work to calculations of the relevant vortex equilibria and to results obtained with other methods. Axisymmetric contour dynamics and the treatment of the multi-layer model of quasigeostrophic flows are described in Section 4 while Section 5 is devoted to a discussion of the tendency shown by vorticity jumps to undergo the strange and subtle phenomenon of filamentation.

Parabolized reduced Navier-Stokes computational techniques

Abstract: A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Numerical models of mantle convection

Abstract: An overview of numerical methods describing the structure and dynamics of the mantle is presented with attention given to novel 3D modeling techniques. The paper reviews 3D spherical and Cartesian models for constant viscosity emphasizing the assumptions regarding style of convection, time dependence, and implications for the mantle. Similarly treated are 3D Cartesian models with temperature-dependent viscosities, and briefly examined are models that are based on compressibility, nonlinear viscosity, or plates. Extensive illustrations are presented detailing: (1) temperature variations from models of 3D thermal convection in spherical shells; (2) thermal anomalies in equatorial cross sections; and (3) temperature variations in a spherical shell heated from within. The discussion relates the numerical results of the models with real mantle-convection events, and the simulations are shown to yield increasingly realistic representations of material behavior.