bluuebunny/crossref_metadata_2025 huggingface dataset:
The influence of coherent structures and microfronts on scaling laws using global and local transforms
Authors: Mahrt L. , Howell J. F.
Published: April 2006

This study examines the influence of coherent structures and attendant microfronts on scaling laws. Toward this goal, we analyse atmospheric observations of turbulence collected 45 m above a flat surface during the Lammefjord Experiment in Denmark. These observations represent more than 40 hours of nearly stationary strong wind conditions and include more than 1600 samples of the main coherent structures. These samples occupy about 40% of the total record and explain the majority of the Reynolds stress.

To study the dependence of the scaling laws on the choice of basis set, the time series of velocity fluctuations are decomposed into Fourier modes, the local Haar basis set and eigenvectors of the lagged covariance matrix. The three decompositions are compared by formulating joint projections. The decompositions are first applied to the samples of phased-locked coherent structures centred about eddy microfronts. The eigenvector decomposition is able to partially separate the small-scale variances due to the coherent eddy microfronts from that due to the small-scale structure with random phase. In the Fourier spectrum, both of these contributions to the variance appear together at the higher wavenumbers and their individual contributions cannot be separated. This effect is relatively minor for the scale distribution of energy but exerts an important influence on higher-moment statistics. Deviations from the −$\frac53$ scaling are observed to be slight and depend on choice of basis set.

The microfronts strongly influence the higher-order statistics such as the sixth-order structure function traditionally used to estimate the energy transfer variance. The intermittency of fine-scale structure, energy transfer variance and dissipation are not completely characterized by random phase, as often assumed, but are partly associated with microfronts characterized by systematic phase with respect to the main transporting eddies. These conclusions are supported by both the higher-order structure function and the higher-order Haar transform.

The Fourier and Haar spectra are also computed for the entire record. The peak of the Haar energy spectrum occurs at smaller scales than those of the Fourier spectrum. The Haar transform is local and emphasizes the width of the events. The Fourier spectrum peaks at the scale of the main periodicity, if it exists, which includes the spacing between the events.
The scattering of sound waves by a vortex: numerical simulations and analytical solutions
Authors: Colonius Tim , Lele Sanjiva K. , Moin Parviz
Published: April 2006

The scattering of plane sound waves by a vortex is investigated by solving the compressible Navier–-Stokes equations numerically, and analytically with asymptotic expansions. Numerical errors associated with discretization and boundary conditions are made small by using high-order-accurate spatial differentiation and time marching schemes along with accurate non-reflecting boundary conditions. The accuracy of computations of flow fields with acoustic waves of amplitude five orders of magnitude smaller than the hydrodynamic fluctuations is directly verified. The properties of the scattered field are examined in detail. The results reveal inadequacies in previous vortex scattering theories when the circulation of the vortex is non-zero and refraction by the slowly decaying vortex flow field is important. Approximate analytical solutions that account for the refraction effect are developed and found to be in good agreement with the computations and experiments.
Patterns and quasi-patterns in the Faraday experiment
Authors: Edwards W. S. , Fauve S.
Published: April 2006

Parametric excitation of surface waves via forced vertical oscillation of a container filled with fluid (the Faraday instability) is investigated experimentally in a small-depth large-aspect-ratio system, with a viscous fluid and with two simultaneous forcing frequencies. The asymptotic pattern observed just above the threshold for the first instability of the flat surface is found to depend strongly on the frequency ratio and the amplitudes and phases of the two sinusoidal components of the driving acceleration. Parallel lines, squares, and hexagons are observed. With viscosity 100 cS, these stable standing-wave patterns do not exhibit strong sidewall effects, and are found in containers of various shapes including an irregular shape. A ‘quasi-pattern’ of twelvefold symmetry, analogous to a two-dimensional quasi-crystal, is observed for some even/odd frequency ratios. Many of the experimental phenomena can be modelled via cubic-order amplitude equations derived from symmetry arguments.
Unsteady electrophoretic motion of a non-spherical colloidal particle in an oscillating electric field
Authors: Loewenberg Michael
Published: April 2006

The oscillatory motion of an electrically charged non-spherical colloidal particle in an oscillating electric field is investigated. The particle is immersed in an incompressible viscous fluid and assumed to have a thin electric double layer. For moderate-aspect-ration spheroids and cylinders, a simple algebraic expression is derived that accurately describes oscillatory electrophoretic particle motion in terms of the steady Stokes resistance, added mass, and Basset force. The effects of double-layer conduction and displacement currents within dielectric particles are included. The results indicate that electroacoustic measurements may be able to determine the ζ-potential, dielectric constant, surface conductivity (and microstructural information contained therein), size, density, volume fraction, and possibly shape of non-spherical particles in a dilute suspension. A simple formula is obtained for the high-frequency electrical conductivity of a dilute suspension of colloidal spheroids with arbitrary charge and dielectric constant; only the added mass and Basset force are required and the requisite parameters are given. The result is needed for electroacoustic measurements but it may also be independently useful for determining the dielectric constant, surface conductivity, volume fraction, and possibly the shape of non-spherical particles in a dilute suspension. Electroacoustic energy dissipation is described for a dilute colloidal suspension. It is shown that resistive electrical heating and viscous dissipation occur independently. Electrical and viscous dissipation coefficients that characterize the order volume fraction contributions of the suspended particles are calculated; the electrical dissipation coefficient is O(1) for all oscillation frequencies, whereas the latter vanishes at low- and high-frequencies. The fluid motion is shown to be a superposition of unsteady, viscous and potential flows past an oscillating particle with no applied electric field. The electro-osmotic flow field is insensitive to particle geometry and qualitatively different from the flow past an oscillating particle with no applied field.
Stability analysis and large-eddy simulation of rotating turbulence with organized eddies
Authors: Cambon Claude , Benoit Jean-Pierre , Shao Liang , Jacquin Laurent
Published: April 2006

Rotation strongly affects the stability of turbulent flows in the presence of large eddies. In this paper, we examine the applicability of the classic Bradshaw-Richardson criterion to flows more general than a simple combination of rotation and pure shear. Two approaches are used. Firstly the linearized theory is applied to a class of rotating two-dimensional flows having arbitrary rates of strain and vorticity and streamfunctions that are quadratic. This class includes simple shear and elliptic flows as special cases. Secondly, we describe a large-eddy simulation of initially quasi-homogeneous three-dimensional turbulence superimposed on a periodic array of two-dimensional Taylor-Green vortices in a rotating frame.

The results of both approaches indicate that, for a large structure of vorticity W and subject to rotation Ω, maximum destabilization is obtained for zero tilting vorticity (½W + 2Ω = 0) whereas stability occurs for zero absolute vorticity (2Ω = 0) These results are consistent with the Bradshaw-Richardson criterion; however the numerical results show that in other cases the Bradshaw-Richardson number $B=2\Omega(W+2\Omega)/W^2$ is not always a good indicator of the flow stability.
Numerical modelling of flux spectra formation for surface gravity waves
Authors: Polnikov V. G.
Published: April 2006

By means of direct numerical solution of the kinetic equation for surface gravity waves, it is shown that under certain conditions the constant flux spectra of nonlinear waves, first predicted by Zakharov & Filonenko (1966) for an infinite frequency domain, can be formed in a finite frequency interval. For the case of angular isotropic spectra the conditions and timescales of this flux spectra formation are evaluated.
Realizability conditions for the turbulent stress tensor in large-eddy simulation
Authors: Vreman Bert , Geurts Bernard , Kuerten Hans
Published: April 2006

The turbulent stress tensor in large-eddy simulation is examined from a theoretical point of view. Realizability conditions for the components of this tensor are derived, which hold if and only if the filter function is positive. The spectral cut-off, one of the filters frequently used in large-eddy simulation, is not positive. Consequently, the turbulent stress tensor based on spectrally filtered fields does not satisfy the realizability conditions, which leads to negative values of the generalized turbulent kinetic energyk. Positive filters, e. g. Gaussian or top-hat, always give rise to a positivek. For this reason, subgrid models which require positive values forkshould be used in conjunction with e. g. the Gaussian or top-hat filter rather than with the spectral cutoff filter. If the turbulent stress tensor satisfies the realizability conditions, it is natural to require that the subgrid model for this tensor also satisfies these conditions. With respect to this point of view several subgrid models are discussed. For eddy-viscosity models a lower bound for the generalized turbulent kinetic energy follows as a necessary condition. This result provides an inequality for the model constants appearing in a ‘Smagorinsky-type’ subgrid model for compressible flows.
New solutions for the propagation of long water waves over variable depth
Authors: Zhang Yinglong , Zhu Songping
Published: April 2006

Based on the linearized long-wave equation, two new analytical solutions are obtained respectively for the propagation of long surface gravity waves around a conical island and over a paraboloidal shoal. Having been intensively studied during the last two decades, these two problems have practical significance and are physically revealing for wave propagation over variable water depth. The newly derived analytical solutions are compared with several previously obtained numerical solutions and the accuracy of those numerical solutions is discussed. The analytical method has the potential to be used to find solutions for wave propagation over more natural bottom topographies.
A family of stochastic models for two-particle dispersion in isotropic homogeneous stationary turbulence
Authors: Borgas M. S. , Sawford B. L.
Published: April 2006

A family of Lagrangian stochastic models for the joint motion of particle pairs in isotropic homogeneous stationary turbulence is considered. The Markov assumption and well-mixed criterion of Thomson (1990) are used, and the models have quadratic-form functions of velocity for the particle accelerations. Two constraints are derived which formally require that the correct one-particle statistics are obtained by the models. These constraints involve the Eulerian expectation of the ‘acceleration’ of a fluid particle with conditioned instantaneous velocity, given either at the particle, or at some other particle's position. The Navier-Stokes equations, with Gaussian Eulerian probability distributions, are shown to give quadratic-form conditional accelerations, and models which satisfy these two constraints are found. Dispersion calculations show that the constraints do not always guarantee good one-particle statistics, but it is possible to select a constrained model that does. Thomson's model has good one-particle statistics, but is shown to have unphysical conditional accelerations. Comparisons of relative dispersion for the models are made.
An experimental study of dipolar vortex structures in a stratified fluid
Authors: Flór J. B. , Van Heijst G. J. F.
Published: April 2006

This paper describes laboratory experiments on dipolar vortex structures in a linearly stratified fluid. The dipoles are generated by a pulsed horizontal injection of a small volume of fluid, by which a localized three-dimensionally turbulent flow region is created. After the subsequent gravitational collapse the flow becomes approximately two-dimensional, and eventually a single vortex dipole emerges, as the result of the self-organizing properties of such flows. The flow evolution has been visualized both by dye and tracer particles, through which qualitative as well as quantitative information was obtained. By application of digital image analysis, the spatial distribution of the velocity ν, vorticity ω and stream function ψ were determined. It was found that dipoles in the turbulent-injection experiments are characterized by a nonlinear sinh-like relationship between ω and ψ, whereas in the case of laminar injection the (ω, ψ)-scatter plots of the dipoles reveal a linear relationship. Notwithstanding these differences, both types of dipoles show a dynamical structure that agrees very well with the Lamb–Chaplygin dipole, as was found by comparing the size, position of maximum vorticity, cross-sectional distributions of ν and ω, characteristics of the (ω, ψ)-relationship, and the translation speed of the experimental and the model dipole.
Hyperviscous vortices
Authors: Jiménez Javier
Published: April 2006

The structure of diffusing planar and axisymmetric vortices of the hyperviscous Navier-Stokes equations is studied for different orders of the dissipative operator. It is found that, except for the classical Newtonian case, the vorticity decays at large distances by means of oscillatory tails, containing circulation of alternating signs. This oscillation becomes stronger for large hyperviscosity orders, and the limit of infinite order is studied. It is argued that these solutions would become unstable for large enough Reynolds numbers, and may contribute non-trivial spurious dynamics to flow simulations using hyperviscosity.
Turbulent mixing in stratified fluids: layer formation and energetics
Authors: Park Young-Gyu , Whitehead J. A. , Gnanadeskian Anand
Published: April 2006

Water with constant initial salt stratification was mixed with a horizontally moving vertical rod. The initially linear density profile turned into a series of steps when mixing was weak, in agreement with instability theory by Phillips (1972) and Posmentier (1977). For stronger mixing no steps formed. However, in all cases mixed layers formed next to the top and bottom boundaries and expanded into the interior due to the no-flux condition at the horizontal boundaries. The critical Richardson number Rie, dividing experiments with steps and ones without, increases with Reynolds number Re as Rie ≈ exp(Re/900). Steps evolved over time, with small ones forming first and larger ones appearing later. The interior seemed to reach an equilibrium state with a collection of stationary steps. The boundary mixed layers continued to penetrate into the interior. They finally formed two mixed layers separated by a step, and ultimately acquired the same densities so the fluid became homogeneous. The length scale of the equilibrium steps, ls, is a linear function of U/Ni, where U is the speed of the stirring rod and Ni is the buoyancy frequency of the initial stratification. The mixing efficiency Rf also evolved in relation to the evolution of the density structure. During the initiation of the steps, Rf showed two completely different modes of evolution depending on the overall Richardson number of the initial state, Rio. For Rio [Gt ] Rie, Rf increased initially. However for Rio near Rie, Rf decreased. Then the steps reached an equilibrium state where Rf was constant at a value that depended on the initial stratification. The density flux was measured to be uniform in the layered interior irrespective of the interior density gradient during the equilibrium state. Thus, the density (salt) was transported from the bottom boundary mixed layer through the layered interior to the top boundary mixed layer without changing the interior density structure. The relationship between Ril and Rf was found for Ril > 1, where Ril is the Richardson number based on the thickness of the interface between the mixed layers. Rf decreases as Ril increases, consistent with the most crucial assumption of the instability theory of Phillips/Posmentier.
Instantaneous three-dimensional concentration measurements in the self-similar region of a round high-Schmidt-number jet
Authors: Yoda M. , Hesselink L. , Mungal M. G.
Published: April 2006

The virtually instantaneous three-dimensional concentration fields in the self-similar region of natural or unexcited, circularly excited and weakly buoyant round jets of Reynolds number based on nozzle diameter of 1000 to 4000 are measured experimentally at a spatial resolution of the order of the Kolmogorov length scale. Isoconcentration surfaces are extracted from the concentration field. These surfaces along with their geometrical parameters are used to deduce the structure and modal composition of the jet. The concentration gradient field is calculated, and its local topology is classified using critical-point concepts.

Large-scale structure is evident in the form of ‘clumps’ of higher-concentration jet fluid. The structure, which has a downstream extent of about the local jet diameter, is roughly axisymmetric with a conical downstream end. This structure appears to be present only in fully turbulent jets. The antisymmetric two-dimensional images previously thought to be axial slices of an expanding spiral turn out in our data to instead be slices of a simple sinusoid in three dimensions. This result suggests that the helical mode, when present, is in the form of a pair of counter-rotating spirals, or that the +1 and −1 modes are simultaneously present in the flow, with their relative phase set by initial conditions.

In terms of local structure, regions with a large magnitude in concentration gradient are shown to have a local topology which is roughly axisymmetric and compressed along the axis of symmetry. Such regions, which would be locally planar and sheet-like, may correspond to the superposition of several of the layer-like structures which are the basic structure of the fine-scale passive scalar field (Buch & Dahm 1991; Ruetsch & Maxey 1991).
Nonlinear resonance of free surface waves in a current over a sinusoidal bottom: a numerical study
Authors: Sammarco Paolo , Mei Chiang C. , Trulsen Karsten
Published: April 2006

We examine the free surface flow over a fixed bed covered by rigid sinusoidal dunes. The mean current velocity is near the critical value at which the linearized theory predicts unbounded response. By allowing transients we examine the instability of the steady and nonlinear solution of Mei (1969) and the possibility of chaos when the current has a small oscillatory component.
The teapot effect: sheet-forming flows with deflection, wetting and hysteresis
Authors: Kistler S. F. , Scriven L. E.
Published: April 2006

The flow of a two-dimensional viscous film falling from the edge of an inclined plane exhibits a distinctive set of phenomena which, in various combinations, have been referred to as the teapot effect. This paper makes plain that three basic mechanisms are at the root of these phenomena: deflection of the liquid sheet by hydrodynamic forces, contact-angle hysteresis, and multiple steady states that give rise to a purely hydrodynamic hysteresis. The evidence is drawn from Galerkin/finite-element analysis of the Navier-Stokes system, matched to a one-dimensional asymptotic approximation of the sheet flow downstream, and is corroborated experimentally by flow visualization and measurements of free-surface profiles and contact line position. The results indicate that the Gibbs inequality condition quantifies the inhibiting effect of sharp edges on spreading of static contact lines, even in the presence of flow nearby. The branchings, turning points, and isolas of families of solutions in parameter space explain abrupt flow transitions observed experimentally, and illuminate the stability of predicted flow states.
Periodically arranged co-flowing jets
Authors: Villermaux E. , Hopfinger E. J.
Published: April 2006

The problem of a periodic planar arrangement of a large number of co-flowing, interacting jets is investigated. It is shown that this interaction gives rise to strong nearfield oscillations of large-scale spatial coherence and to far-field inhomogeneities. In the experiments performed, the jets were produced behind a flat plate perforated by holes arranged in a square or triangular periodic pattern and placed perpendicular to a uniform flow. At moderate Reynolds numbers, the interaction results in a remarkable low-frequency oscillation of the merging distance of the jets downstream of the plate. A detailed description of the recirculating flow in the cavities between the jets emphasizes the role of the backflow in the cavities on the oscillatory behaviour. This description is supported by measurements of the local fluctuating velocity and pressure, two-point correlation measurements and quantitative flow visualizations. These experimental observations suggest a new formulation for the instability dynamics of such unstable recirculating flows. This formulation, based on the nonlinear delayed saturation of the jet's shear layer instability (NLDS model) predicts successfully the dependence of the oscillation of the merging distance on the jet Reynolds number and on the local geometrical features of the confinement of the jets. Furthermore, it is shown that the diffusion of mass coming from one jet, seeded with an inert dye, gives rise to an exponential diffusion front over a distance corresponding to a few mesh sizes indicating a strong local coupling of the jets. At the scale of the whole jet assembly, the oscillations are organized as large-scale travelling waves, propagating from the boundaries of the domain to its centre. This symmetry-breaking property is discussed and supplemented by a spatio-temporal simulation of an array of coupled oscillators.
Nonlinear dynamics of forced transitional jets: periodic and chaotic attractors
Authors: Broze George , Hussain Fazle
Published: April 2006

Conclusive experimental evidence is presented for the existence of a low-dimensional temporal dynamical system in an open flow, namely the near field of an axisymmetric, subsonic free jet. An initially laminar jet (4 cm air jet in the Reynolds number range 1.1 × 104[Lt ]ReD× 9.1 × 104) with a top-hat profile was studied using single-frequency, longitudinal, bulk excitation. Two non-dimensional control parameters – forcing frequencyStD(≡fexD/Ue, wherefezis the excitation frequency,Dis the jet exit diameter andUeis the exit velocity) and forcing amplitudeaf(≡u’f/Ue, whereu’fis the jet exit r.m.s. longitudinal velocity fluctuation at the excitation frequency) – were varied over the ranges 10-4<af< 0.3 and 0.3 <StD< 3.0 in order to construct a phase diagram. Periodic and chaotic states were found over large domains of the parameter space. The periodic attractors correspond to stable pairing (SP) and stable double pairing (SDP) of rolled-up vortices. One chaotic attractor, near SP in the parameter space, results from nearly periodic modulations of pairing (NPMP) of vortices. At large scales (i.e. approximately the size of the attractor) in phase space, NPMP exhibits approximately quasi-periodic behaviour, including modulation sidebands around ½fexinu-spectra, large closed loops in its Poincaré sections, correlation dimensionv∼ 2 and largest Lyapunov exponent λ1∼ 0. But investigations at smaller scales (i.e. distances greater than, but of the order of, trajectory separation) in phase space reveal chaos, as shown byv> 2 and λ1> 0. The other chaotic attractor, near SDP, results from nearly periodic modulations of the first vortex pairing but chaotic modulations of the second pairing and has a broadband spectrum, a dimension 2.5 [Lt ]v[Lt ] 3 and the largest Lyapunov exponent 0.2 [Lt ] λ1[Lt ] 0.7 bits per orbit (depending on measurement locations in physical and parameter spaces).

A definition that distinguishes between physically and dynamically open flows is proposed and justified by our experimental results. The most important conclusion of this study is that a physically open flow, even one that is apparently dynamically open due to convective instability, can exhibit dynamically closed behaviour as a result of feedback. A conceptual model for transitional jets is proposed based on twodimensional instabilities, subharmonic resonance and feedback from downstream vortical structures to the nozzle lip. Feedback was quantified and shown to affect the exit fundamental–subharmonic phase difference ϕ – a crucial variable in subharmonic resonance and, hence, vortex pairing. The effect of feedback, the sensitivity of pairings to ϕ, the phase diagram, and the documented periodic and chaotic attractors demonstrate the validity of the proposed conceptual model.
Roll waves on a shallow layer of mud modelled as a power-law fluid
Authors: Ng Chiu-On , Mei Chiang C.
Published: April 2006

We give a theory of permanent roll waves on a shallow layer of fluid mud which is modelled as a power-law fluid. Based on the long-wave approximation, Kármán's momentum integral method is applied to derive the averaged continuity and the momentum equations. Linearized instability analysis of a uniform flow shows that the growth rate of unstable disturbances increases monotonically with the wavenumber, and therefore is insufficient to suggest a preferred wavelength for the roll wave. Nonlinear roll waves are obtained next as periodic shocks connected by smooth profiles with depth increasing monotonically from the rear to the front. Among all wavelengths only those longer than a certain threshold correspond to positive energy loss across the shock, and are physically acceptable. This threshold also implies a minimum discharge, viewed in the moving system, for the roll wave to exist. These facts suggest that a roll wave developed spontaneously from infinitesimal disturbances should have the shortest wavelength corresponding to zero dissipation across the shock, though finite dissipation elsewhere. The discontinuity at the wave front is a mathematical shortcoming needing a local requirement. Predictions for the spontaneously developed roll waves in a Newtonian case are compared with available experimental data. Longer roll waves, with dissipation at the discontinuous fronts, cannot be maintained if the uniform flow is linearly stable, when the fluid is slightly non-Newtonian. However, when the fluid is highly non-Newtonian, very long roll waves may still exist even if the corresponding uniform flow is stable to infinitesimal disturbances. Numerical results are presented for the phase speed, wave height and wavenumber, and wave profiles for a representative value of the flow index of fluid mud.
Dendritic growth from a melt in an external flow: uniformly valid asymptotic solution for the steady state
Authors: Xu Jian-Jun
Published: April 2006

This work is part of an investigation dealing with the effect of external flow in a melt on dendritic growth. In this paper, we will consider steady growth with zero surface tension. Assuming that the Prandtl number Pr is large, we are able to obtain a uniformly valid asymptotic solution for the steady state in the whole growth region.
On the three-dimensional internal waves excited by topography in the flow of a stratified fluid
Authors: Hanazaki Hideshi
Published: April 2006

A numerical study of the three-dimensional internal waves excited by topography in the flow of a stratified fluid is described. In the resonant flow of a nearly two-layer fluid, it is found that the time-development of the nonlinearly excited waves agrees qualitatively with the solution of the forced KP equation or the forced extended KP equation. In this case, the upstream-advancing solitary waves become asymptotically straight crested because of abnormal reflection at the sidewall similar to Mach reflection. The same phenomenon also occurs in the subcritical flow of a nearly two-layer fluid. However, in the subcritical flow of a linearly stratified Boussinesq fluid, the two-dimensionalization of the upstream waves can be interpreted as the separation of the lateral modes due to the differences in the group velocity of the linear wave, although this does not mean in general that the generation of upstream waves is describable by the linearized equation.
Blunt-body impact on the free surface of a compressible liquid
Authors: Korobkin Alexander
Published: April 2006

The present paper deals with the plane unsteady problem of the penetration of a blunt solid contour into an ideal compressible liquid. At the initial instant of time, the solid body touches the liquid free boundary at a single point. At the initial stage, the duration of which depends on the body geometry, the displacements of liquid particles are small and the disturbed fluid motion is described within the framework of the acoustic approximation. The main feature of the problem is the existence of a contact line between the free surface of the liquid and the solid-body surface. The position of this line is not known in advance and is to be determined together with the solution of the problem. A brief description of the method that provides a solution of complicated nonlinear problems such as this is given. The pressure distribution and the velocity field in the liquid are shown to be given in quadratures and in the case of a parabolic entering contour in an explicit form. For a parabolic entering contour the pressure at the top of the contour calculated using the model of an incompressible liquid is observed to deviate from a precise value by not more than 10% of the latter after the first expansion waves have passed the contact point. The solution analysis enabled us to distinguish the regions in which the acoustic approximation fails and the liquid flow becomes essentially nonlinear.
Interfacial instability due to MHD mode coupling in aluminium reduction cells
Authors: Sneyd A. D. , Wang A.
Published: April 2006

This paper analyses instabilities on the cryolite/aluminium interface in an aluminium reduction cell. The simplified cell model is a finite rectangular tank containing the two fluid layers, and carrying a uniform normal current. The magnetic field is assumed to be a linear function of position. Several previous studies have considered waves consisting of a single Fourier component but here we consider perturbations which are a general combination of the normal gravity-wave modes. We derive a system of coupled ordinary differential equations for the time-development of the mode amplitudes, and show that instability can occur via mode interactions, the electromagnetic perturbation force due to one mode feeding energy into the other. Growth rates are determined by computing the eigenvalues of an interaction matrix, and an approximate method using only the three leading diagonals is developed. If two modes have similar frequencies they may resonate and become unstable at a very low threshold current. We consider the influence of various cell parameters and draw some general conclusions about cell design.
Large deformations of elastic cylindrical capsules in shear flows
Authors: Rao Peddada R. , Zahalak George I. , Sutera Salvatore P.
Published: April 2006

The nonlinear problem of the steady-state interaction of a closed fluid-filled cylindrical elastic membrane with a slow viscous shear flow has been solved by a series-expansion technique. The problems of successive orders were both formulated and solved by a symbolic manipulation program, and the calculations were carried to sixth order in a dimensionless parameter related to the applied shear rate. Moderately large deformations (aspect ratios approaching 3) fall within the range of this analysis, which yields the dependences of the following global variables on the system parameters: membrane deformation, orientation, and strain, as well as tank-treading frequency, and mean internal pressure. The solution for the flow field around an isolated capsule is also used to calculate the apparent viscosity of a dilute suspension of flexible cylindrical particles, which yields the paradoxical result that the apparent viscosity decreases as the internal viscosity increases.
Quantitative measurements of three-dim ensional structures in the wake of a circular cylinder
Authors: Mansy Hussein , Yang Pan-Mei , Williams David R.
Published: April 2006

The fine scale three-dimensional structures usually associated with streamwise vortices in the near wake of a circular cylinder have been studied at Reynolds numbers ranging from 170 to 2200. Spatially continuous velocity measurements along lines parallel to the cylinder axis were obtained with a scanning laser anemometer. To detect the streamwise vortices in the amplitude modulated velocity field, it was necessary to develop a spatial decomposition technique to split the total flow into a primary flow component and a secondary flow component. The primary flow is comprised of the mean flow and Strouhal vortices, while the secondary flow is the result of the three-dimensional streamwise vortices that are the essence of transition to turbulence. The three-dimensional flow amplitude increases in the primary vortex formation region, then saturates shortly after the maximum amplitude in the primary flow is reached. In the near-wake region the wavelength decreases approximately like Re−0.5, but increases with downstream distance. A discontinuous increase in wavelength occurs below Re = 300 suggesting a fundamental change in the character of the three-dimensional flow. At downstream distances (x/D = 10-20), the spanwise wavelength decreases from 1.42D to 1.03D as the Reynolds number increases from 300 to 1200.
A global stability analysis of the steady and periodic cylinder wake
Authors: Noack Bernd R. , Eckelmann Helmut
Published: April 2006

A global, three-dimensional stability analysis of the steady and the periodic cylinder wake is carried out employing a low-dimensional Galerkin method. The steady flow is found to be asymptotically stable with respect to all perturbations for Re < 54. The onset of periodicity is confirmed to be a supercritical Hopf bifurcation which can be modelled by the Landau equations. The periodic solution is observed to be only neutrally stable for 54 < Re < 170. While two-dimensional perturbations of the vortex street rapidly decay, three-dimensional perturbations with long spanwise wavelengths neither grow nor decay. The periodic solution becomes unstable at Re = 170 by a perturbation with the spanwise wavelength of 1.8 diameters. This instability is shown to be a supercritical Hopf bifurcation in the spanwise coordinate and leads to a three-dimensional periodic flow. Finally the transition scenario for higher Reynolds numbers is discussed.
Evolution of disturbances in stagnation-point flow
Authors: Criminale W. O. , Jackson T. L. , Lasseigne D. G.
Published: April 2006

The evolution of three-dimensional disturbances in an incompressible three-dimensional stagnation-point flow in an inviscid fluid is investigated. Since it is not possible to apply classical normal-mode analysis to the disturbance equations for the fully three-dimensional stagnation-point flow to obtain solutions, an initial-value problem is solved instead. The evolution of the disturbances provides the necessary information to determine stability and indeed the complete transient as well. It is found that when considering the disturbance energy, the planar stagnation-point flow, which is independent of one of the transverse coordinates, represents a neutrally stable flow whereas the fully three-dimensional flow is either stable or unstable, depending on whether the flow is away from or towards the stagnation point in the transverse direction that is neglected in the planar stagnation point.
Influence of the sampling time on the kinematics of turbulent diffusion from a continuous source
Authors: Eckman Richard M.
Published: April 2006

A kinematic description is presented of how turbulent diffusion from a continuous source varies with the sampling time in stationary, homogeneous turbulence. Unlike most previous theories, the sampling is assumed to take place at fixed downstream distances from the source. It is shown that the sampling-time effects depend on two-particle velocity statistics. Thus, time-average diffusion at fixed downstream distances is more akin to relative diffusion than to absolute diffusion. Two simple diffusion models are developed from the kinematic equations. These models are in fairly good agreement with diffusion data obtained both in a wind tunnel and in the field. Moreover, these models have significant practical implications. For example, the models indicate that care must be taken when using band-pass spectral filtering as a paradigm for turbulent diffusion. Also, the models show that the mean flow speed U has an important influence on the sampling-time effects. To account for U properly, diffusion measurements with differing sampling times Λ should be compared using the product UΛ, and not just Λ.
On reduced equations in the Hamiltonian theory of weakly nonlinear surface waves
Authors: Krasitskii Vladimir P.
Published: April 2006

Many studies of weakly nonlinear surface waves are based on so-called reduced integrodifferential equations. One of these is the widely used Zakharov four-wave equation for purely gravity waves. But the reduced equations now in use are not Hamiltonian despite the Hamiltonian structure of exact water wave equations. This is entirely due to shortcomings of their derivation. The classical method of canonical transformations, generalized to the continuous case, leads automatically to reduced equations with Hamiltonian structure. In this paper, attention is primarily paid to the Hamiltonian reduced equation describing the combined effects of four- and five-wave weakly nonlinear interactions of purely gravity waves. In this equation, for brevity called five-wave, the non-resonant quadratic, cubic and fourth-order nonlinear terms are eliminated by suitable canonical transformation. The kernels of this equation and the coefficients of the transformation are expressed in explicit form in terms of expansion coefficients of the gravity-wave Hamiltonian in integral-power series in normal variables. For capillary–gravity waves on a fluid of finite depth, expansion of the Hamiltonian in integral-power series in a normal variable with accuracy up to the fifth-order terms is also given.
The solidification of buoyancy-driven flow in a flexible-walled channel. Part 2. Continual release
Authors: Lister John R.
Published: April 2006

The model developed in Part 1 (Lister 1994) for the solidification of hot fluid flowing in a thin buoyancy-driven layer between cold solid but freely deformable boundaries is extended to study the case of continual release of fluid. In this model lubrication theory was applied to reduce the equations of mass and heat conservation to a kinematic-wave equation and an advection-diffusion equation, which were coupled by the rate of solidification. The equations allow the source flux to be specified, and the cases of constant input and of flux proportional to a power of time are considered here. The structure of the flow differs significantly from the case of constant-volume release considered in Part 1. The advective resupply of heat prevents the flow from solidifying completely at the source and, if the initial fluid temperature is greater than the melting temperature of the solid, will in fact lead to rapid melting near the source. A perturbation expansion is used to describe the development of thermal boundary layers at the flow margins and the initial self-similar extension of the zone of melting. As the flow propagates beyond its thermal entry length, the fluid temperature falls to the liquidus value and melting gives way to solidification. At large times nearly all of the fluid supplied solidifies against the margins of the flow but, provided the source flux decreases less rapidly than t−½, sufficient reaches the nose of the flow that the flow continues to increase in length indefinitely. Analytic solutions are given for this long-time regime showing, for example, that the length increases asymptotically like t1/2 for constant-flux input. The theoretical solutions, which are calculated by a combination of analytic and numerical methods, may be used to describe the propagation of a dyke fed by a large body of magma through the Earth's lithosphere or the flow of lava down the flanks of a volcano during an extensive period of eruption.
The long-time self-diffusivity in concentrated colloidal dispersions
Authors: Brady John F.
Published: April 2006

The long-time self-diffusivity in concentrated colloidal dispersions is determined from a consideration of the temporal decay of density fluctuations. For hydrodynamically interacting Brownian particles the long-time self-diffusivity, Ds∞, is shown to be expressible as the product of the hydrodynamically determined short-time self-diffusivity, Ds(ϕ), and a contribution that depends on the distortion of the equilibrium structure caused by a diffusing particle. An argument is advanced to show that as maximum packing is approached the long-time self-diffusivity scales as Ds∞(ϕ) ∼ Ds0(ϕ)/g(2; ϕ), where g(2; ϕ) is the value of the equilibrium radial-distribution function at contact and ϕ is the volume fraction of interest. This result predicts that the longtime self-diffusivity vanishes quadratically at random close packing, ϕm ≈ 0.63, i.e. Ds∞D0(1-ϕ/ϕm)2 as ϕ→ϕm, where D0 = kT/6πνa is the diffusivity of a single isolated particle of radius a in a fluid of viscosity ν. This scaling occurs because Ds0(ϕ) vanishes linearly at random close packing and the radial-distribution function at contact diverges as (1 -ϕ/ϕm)−1. A model is developed to determine the structural deformation for the entire range of volume fractions, and for hard spheres the longtime self-diffusivity can be represented by: Ds∞(ϕ) = Ds∞(ϕ)/[1 + 2ϕg(2;ϕ)]. This formula is in good agreement with experiment. For particles that interact through hard-spherelike repulsive interparticle forces characterized by a length b(> a), the same formula applies with the short-time self-diffusivity still determined by hydrodynamic interactions at the true or ‘hydrodynamic’ volume fraction ϕ, but the structural deformation and equilibrium radial-distribution function are now determined by the ‘thermodynamic’ volume fraction ϕb based on the length b. When b [Gt ] a, the long-time self-diffusivity vanishes linearly at random close packing based on the ‘thermodynamic’ volume fraction ϕbm. This change in behaviour occurs because the true or ‘hydrodynamic’ volume fraction is so low that the short-time self-diffusivity is given by its infinite-dilution value D0. It is also shown that the temporal transition from short- to long-time diffusive behaviour is inversely proportional to the dynamic viscosity and is a universal function for all volume fractions when time is nondimensionalized by a2/Ds∞(ϕ).
On Chaplygin's investigations of two-dimensional vortex structures in an inviscid fluid
Authors: Meleshko V.V. , Heijst G.J.F.van
Published: April 2006

This paper describes exact solutions of two-dimensional vortex structures that were published by Chaplygin (1899, 1903) at the turn of the last century, which seem to have escaped the attention of later investigators in this field. Chaplygin's solutions include that of an elliptical patch of uniform vorticity in an exterior field of pure shear and that of a (symmetric or non-symmetric) dipolar vortex with a continuous distribution of vorticity translating steadily along a straight path. In addition, a solution is presented for a non-symmetric vortex dipole moving along a circular trajectory. A concise account of Chaplygin's solutions is given, complemented by a more detailed analysis of some of their relevant properties.
A Hamiltonian structure with contact geometry for the semi-geostrophic equations
Authors: Roulstone I. , Norbury J.
Published: April 2006

A canonical Hamiltonian structure for the semi-geostrophic equations is presented and from this a reduced non-canonical Hamiltonian structure is derived, providing a fully nonlinear version of the approximate linearized vorticity advection representation. The structure of this model is described naturally within the framework of contact geometry. A Hamiltonian approach leading to a symplectic algorithm for calculating solutions to the equations of motion is formulated. Basic necessary functional methods are introduced and the Lagrangian and Eulerian kinematic structures are discussed, together with their relevance to symplectic integrating algorithms.
The inertial migration of non-neutrally buoyant spherical particles in two-dimensional shear flows
Authors: Hogg Andrew J.
Published: April 2006

The inertial migration of a small rigid spherical particle, suspended in a fluid flowing between two plane boundaries, is investigated theoretically to find the effect on the lateral motion. The channel Reynolds number is of order unity and thus both boundary-induced and Oseen-like inertial migration effects are important. The particle Reynolds number is small but non-zero, and singular perturbation techniques are used to calculate the component of the migration velocity which is directed perpendicular to the boundaries of the channel. The particle is non-neutrally buoyant and thus its buoyancy-induced motion may be either parallel or perpendicular to the channel boundaries, depending on the channel alignment. When the buoyancy results in motion perpendicular to the channel boundaries, the inertial migration is a first-order correction to the magnitude of this lateral motion, which significantly increases near to the boundaries. When the buoyancy produces motion parallel with the channel boundaries, the inertial migration gives the zeroth-order lateral motion either towards or away from the boundaries. It is found that those particles which have a velocity exceeding the undisturbed shear flow will migrate towards the boundaries, whereas those with velocities less than the undisturbed flow migrate towards the channel centreline. This calculation is of practical importance for various chemical engineering devices in which particles must be filtered or separated. It is useful to calculate the forces on a particle moving near to a boundary, through a shear flow. This study may also explain certain migration effects of bubbles and crystals suspended in molten rock flow flowing through volcanic conduits.