bluuebunny/crossref_metadata_2025_split huggingface dataset:
Nonlinear modulations of solitary waves
Authors: Pedersen Geir 
Published: April 2006

The leading optical approximation to a slowly varying solitary crest on constant depth is the plane soliton solution with the local values of amplitude and orientation substituted. This leads to two nonlinear hyperbolic equations for the local amplitude and inclination of the crest that have been reported by several authors and predict the formation of progressive wave jumps, or shocks, from any initial perturbation of the crest. In comparison to numerical solutions of the Boussinesq equations we find that this optical approximation fails to reproduce essential properties of the crest dynamics, in particular that the crest modulations are damped and that well-defined wave jumps do not necessarily evolve. One purpose of the present work is to include such features in an amended optical approximation.We obtain the leading correction to the ‘local soliton’ solution by a multiple scale technique. In addition to a modification to the wave profile the perturbation expansion also yields a diffracted wave system and a celerity that depends on the curvature of the crest. The principle of energy conservation then leads us to a second-order optical approximation consisting of transport equations of mixed hyperbolic/parabolic nature. Under additional assumptions the transport equations can be reduced to the well-known Burgers equation.Numerical simulations of the Boussinesq equations are performed for modulations on otherwise straight crests and radially converging solitons. The improved optical, or ray, theory reproduces all essential features and agrees closely with the numerical solution in both cases. Contrary to purely hyperbolic optical descriptions the present theory also predicts wave jumps of finite width that are consistent with the triad solution of Miles (1977).The present work indicates that while sinusoidal waves often are appropriately described by the lowest-order physical optics, higher-order corrections must be expected to be important for single crested waves.
A new vectorial bedload formulation and its application to the time evolution of straight river channels
Authors: Kovacs Agnes , Parker Gary 
Published: April 2006

The derivation of a new vectorial bedload formulation for the transport of coarse sediment by fluid flow is presented in the first part of the paper. This relation has been developed for slopes up to the angle of repose both in the streamwise and transverse directions. The pressure distribution is assumed to be hydrostatic. The bed shear stress for the onset of particle motion and mean particle velocity are obtained from the mean force balance on a particle. A new generalized Bagnold hypothesis is introduced to calculate the sediment content of the bedload layer. The new formulation possesses two innovative features. It is fully nonlinear and vectorial in nature, in addition, it behaves smoothly up to the angle of repose.A mathematical model of the time evolution of straight river channels is presented in the second half of the paper. This study focuses on the evolution process due to bank erosion in the presence of bedload only. The bed and bank material is taken to be coarse, non-cohesive and uniform in size. The sediment continuity and the fluid momentum conservation equations describe the time evolution of the bed topography and flow field. These equations are coupled through the fluid shear stress acting on the bed. This bed shear stress distribution is predicted with the aid of a simple algebraic turbulent closure model. As regards the computation of the sediment flux, the new fully nonlinear vectorial formulation is found to perform well and renders the evolution model fully mechanistic.The formation of an erosional front in the time development of straight river channels has been so far obscured in physical experiments. Herein, with the help of the new bedload formulation, the existence and migration speed of the front of erosion are inferred from the analysis of the sediment continuity equation.The model successfully describes the time relaxation of an initially trapezoidal channel toward an equilibrium cross-sectional shape, as evidenced by comparison with experimental data. This equilibrium is characterized by a constant width, vanishing sediment transport in the transverse direction, and a small but non-vanishing streamwise transport rate of bed sediment.
Averaged equations for inviscid disperse two-phase flow
Authors: Zhang D. Z. , Prosperetti A. 
Published: April 2006

Averaged equations governing the motion of equal rigid spheres suspended in a potential flow are derived from the equation for the probability distribution. A distinctive feature of this work is the derivation of the disperse-phase momentum equation by averaging the particle equation of motion directly, rather than the microscopic equation for the particle material. This approach is more flexible than the usual one and leads to a simpler and more fundamental description of the particle phase. The model is closed in a systematic way (i.e. with no ad hoc assumptions) in the dilute limit and in the linear limit. One of the closure quantities is related to the difference between the gradient of the average pressure and the average pressure gradient, a well-known problem in the widely used two-fluid engineering models. The present result for this quantity leads to the introduction of a modified added mass coefficient (related to Wallis's ‘exertia’) that remains very nearly constant with changes in the volume fraction and densities of the phases. Statistics of this coefficient are provided and exhibit a rather strong variability of up to 20% among different numerical simulations. A detailed comparison of the present results with those of other investigators is given in § 10.As a further illustration of the flexibility of the techniques developed in the paper, in Appendix C they are applied to the calculation of the so-called ‘particle stress’ tensor. This derivation is considerably simpler than others available in the literature.
On the stability of vertical double-diffusive interfaces. Part 2. Two parallel interfaces
Authors: Eltayeb I. A. , Loper D. E. 
Published: April 2006

This is the second part of a three-part study of the stability of vertically oriented double-diffusive interfaces having an imposed vertical stable temperature gradient. In this study, flow is forced within a fluid of infinite extent by a prescribed excess of compositionally buoyant material between two parallel interfaces. Compositional diffusivity is ignored while thermal diffusivity and viscosity are finite. The stability of the interfaces is analysed first in the limit that they are close together (compared with the salt-finger lengthscale), then for general spacing. Attention is focused on whether the preferred mode of instability is varicose or sinuous and whether its wavevector is vertical or oblique.The interfaces are found to be unstable for some wavenumber for all values of the Prandtl number and interface spacing. The preferred mode of instability for closely spaced interfaces is varicose and vertical for Prandtl number less than about 9, sinuous oblique for Prandtl number between 9 and 15 and sinuous vertical for larger Prandtl number. For general spacing each of the four possible modes of instability is preferred for some range of Prandtl number and interface separation, with no clear pattern of preference, except that the sinuous oblique mode is preferred for widely separated interfaces. The growth rate of the preferred mode is largest for interfaces having separations of from 1 to 3 salt-finger lengths.
Circulation deposition on shock-accelerated planar and curved density-stratified interfaces: models and scaling laws
Authors: Samtaney Ravi , Zabusky Norman J. 
Published: April 2006

Vorticity is deposited baroclinically by shock waves on density inhomogeneities. In two dimenslons, the circulation deposited on a planar interface may be derived analytically using shock polar analysis provided the shock refraction is regular. We present analytical expressions for Γ′, the circulation deposited per unit length of the unshocked planar interface, within and beyond the regular refraction regime. To lowest order, Γ′ scales as \[ \Gamma^\prime\propto (1-\eta^{-\frac{1}{2}})(\sin\alpha)(1+M^{-1}+2M^{-2})(M-1)(\gamma^{\frac{1}{2}}/\gamma + 1)\] where M is the Mach number of the incident shock, η is the density ratio of the gases across the interface, α is the angle between the shock and the interface and γ is the ratio of specific heats for both gases. For α ≤ 30°, the error in this approximation is less than 10% for 1.0 < M ≤ 1.32 for all η > 1, and 5.8 ≤ η ≤ 32.6 for all M. We validate our results by quantification of direct numerical simulations of the compressible Euler equations with a second-order Godunov code.We generalize the results for total circulation on non-planar (sinusoidal and circular) interfaces. For the circular bubble case, we introduce a ‘near-normality’ ansatz and obtain a model for total circulation on the bubble surface that agrees well with results of direct numerical simulations. A comparison with other models in the literature is presented.
Coherent structures in rotating three-dimensional turbulence
Authors: Bartello Peter , Métais Olivier , Lesieur Marcel 
Published: April 2006

Numerical simulations investigating the formation and stability of quasi-two-dimensional coherent vortices in rotating homogeneous three-dimensional flow are described. In a numerical study of shear flows Lesieur, Yanase & Métais (1991) found that cyclones (respectively anticyclones) with |ω2D| ∼O(2Ω), where ω2Dis the vorticity and Ω is the rotation rate, are stabilized (respectively destabilized) by the rotation. A study of triply periodic pseudo-spectral simulations (643) was undertaken in order to investigate the vorticity asymmetry in homogeneous turbulence. Specifically, we examine (i) the possible three-dimensionalization of initially two-dimensional vortices and (ii) the emergence of quasi-two-dimensional structures in initially-isotropic three-dimensional turbulence. Direct numerical simulations of the Navier—Stokes equations are compared with large-eddy simulations employing a subgridscale model based on the second-order velocity structure function evaluated at the grid separation and with simulations employing hyperviscosity.Isolated coherent two-dimensional vortices, obtained from a two-dimensional decay simulation, were superposed with a low-amplitude three-dimensional perturbation, and used to initialize the first set of simulations. With Ω = 0, a three-dimensionalization of all vortices was observed. This occurred first in the small scales in conjunction with the formation of longitudinal hairpin vortices with vorticity perpendicular to that of the initial quasi-two-dimensional flow. In agreement with centrifugal stability arguments, when 2Ω = [ω2D]rmsa rapid destabilization of anticyclones was observed to occur, whereas the initial two-dimensional cyclonic vortices persisted throughout the simulation. At larger Ω, both cyclones and anticyclones remained two-dimensional, consistent with the Taylor—Proudman theorem. A second set of simulations starting from isotropic three-dimensional fields was initialized by allowing a random velocity field to evolve (Ω = 0) until maximum energy dissipation. When the simulations were continued with 2Ω = [ω · Ω]rms/Ω, the three-dimensional flow was observed to organize into two-dimensional cyclonic vortices. At larger Ω, two-dimensional anticyclones also emerged from the initially-isotropic flow. The consequences for a variety of industrial and geophysical applications are clear. For quasi-two-dimensional eddies whose characteristic circulation times are of the order ofder of Ω−1, rotation induces a complete disruption of anticyclonic vortices, while stabilizing cyclonic ones.
The spatial structure of neutral atmospheric surface-layer turbulence
Authors: Mann Jakob 
Published: April 2006

Modelling of the complete second-order structure of homogeneous, neutrally stratified atmospheric boundary-layer turbulence, including spectra of all velocity components and cross-spectra of any combination of velocity components at two arbitrarily chosen points, is attempted. Two models based on Rapid Distortion Theory (RDT) are investigated. Both models assume the velocity profile in the height interval of interest to be approximately linear. The linearized Navier–Stokes equation together with considerations of ‘eddy’ lifetimes are then used to modify the spatial second-order structure of the turbulence. The second model differs from the first by modelling the blocking by the surface in addition to the shear. The resulting models of the spectral velocity tensor contain only three adjustable parameters: a lengthscale describing the size of the largest energy-containing eddies, a non-dimensional number used in the parametrization of ‘eddy’ lifetime, and the third parameter is a measure of the energy dissipation.Two atmospheric experiments, both designed to investigate the spatial structure of turbulence and both running for approximately one year, are used to test and calibrate the models. Even though the approximations leading to the models are very crude they are capable of predicting well the two-point second-order statistics such as cross-spectra, coherences and phases, on the basis of measurements carried out at one point. The two models give very similar predictions, the largest difference being in the coherences involving vertical velocity fluctuations, where the blocking by the surface seems to have a significant effect.
On mixing and transport at a sheared density interface
Authors: Sullivan Greg D. , List E. John 
Published: April 2006

Mixing and transport of a stratifying scalar are investigated at a density interface imbedded in a turbulent shear flow. Steady-state interfacial shear flows are generated in a laboratory water channel for layer Richardson numbers, Ri, between about 1 and 10. The flow field is made optically homogeneous, enabling the use of laser-induced fluorescence with photodiode array imaging to measure the concentration field at high resolution. False-colour images of the concentration field provide valuable insight into interfacial dynamics: when the local mean shear Richardson number, Ris, is less than about 0.40–0.45, interfacial mixing appears to be dominated by Kelvin–Helmholtz (K–H) instabilities; when Ris is somewhat larger than this, interfacial mixing appears to be dominated by shear-driven wave breaking. In both cases, vertical transport of mixed fluid from the interfacial region into adjacent turbulent layers is accomplished by large-scale turbulent eddies which impinge on the interface and scour fluid from its outer edges.Motivated by the experimental findings, a model for interfacial mixing and entrainment is developed. A local equilibrium is assumed in which the rate of loss of interfacial fluid by eddy scouring is balanced by the rate of production (local mixing) by interfacial instabilities and molecular diffusion. When a single layer is turbulent and entraining, the model results are as follows: in the molecular-diffusion-dominated regime, δ/h ∼ Pe−1/2 and E ∼ Ri−1Pe−1/2; in the wave-breaking-dominated regime, δ/h ∼ Ri−1/2 and E ∼ Ri−3/2; and in the K–H-dominated regime, δ/h ∼ Ri−1 and E ∼ Ri−2, where δ is the interface thickness, h is the boundary-layer thickness, Pe is the Péclet number, and E is the normalized entrainment velocity. In all three regimes, the maximum concentration anomaly, [gcy ]m ∼ Ri−1. When both layers are turbulent and entraining, E and δ depend on combinations of parameters from both layers.
Motion of a rigid particle in a rotating viscous flow: an integral equation approach
Authors: Tanzosh John P. , Stone H. A. 
Published: April 2006

A boundary integral method is presented for analysing particle motion in a rotating fluid for flows where the Taylor number${\cal T}$is arbitrary and the Reynolds number is small. The method determines the surface traction and drag on a particle, and also the velocity field at any location in the fluid.Numerical results show that the dimensionless drag on a spherical particle translating along the rotation axis of an unbounded fluid is determined by the empirical formula$D/6\pi = 1+(4/7){\cal T}^{1/2}+(8/9\pi){\cal T}$, which incorporates known results for the low and high Taylor number limits. Streamline portraits show that a critical Taylor number${\cal T}_c\ap 50$exists at which the character of the flow changes. For${\cal T} < {\cal T}_c$the flow field appears as a perturbation of a Stokes flow with a superimposed swirling motion. For${\cal T} > {\cal T}_c$the flow field develops two detached recirculating regions of trapped fluid located fore and aft of the particle. The recirculating regions grow in size and move farther from the particle with increasing Taylor number. This recirculation functions to deflect fluid away from the translating particle, thereby generating a columnar flow structure. The flow between the recirculating regions and the particle has a plug-like velocity profile, moving slightly slower than the particle and undergoing a uniform swirling motion. The flow in this region is matched to the particle velocity in a thin Ekman layer adjacent to the particle surface.A further study examines the translation of spheroidal particles. For large Taylor numbers, the drag is determined by the equatorial radius; details of the body shape are less important.
On the mechanism of shear flow instabilities
Authors: Baines Peter G. , Mitsudera Humio 
Published: April 2006

In homogeneous and density-stratified inviscid shear flows, the mechanism for instability that is most commonly invoked and discussed is Kelvin–Helmholtz instability, as it occurs for a simple velocity discontinuity. There is a second mechanism, the wave-interaction mechanism, which is much more general, and is the subject of this paper. This mechanism depends on two free waves that propagate in opposite directions in a stratified shear flow, and which may become stationary relative to each other because of the shear. If this occurs, and their relative phase is suitably chosen, the velocity field of each wave increases the displacement of the other, and so the disturbance grows.We show that this mechanism is responsible for instability in a general class of symmetric but otherwise arbitrary velocity and density profiles, provided that the Richardson number Ri < ¼ in a central region of arbitrarily small thickness. A critical layer exists in this central region for the growing disturbance, but its role in the instability process is incidental. When Ri > ¼ everywhere, the flow is stable because the free waves described above are absorbed by the critical layer, and hence are heavily damped. The necessary criteria of Rayleigh and Fjortoft for instability in homogeneous fluid are seen to provide a suitable geometry for two interacting waves. Some specific examples are given, including a succinct explanation of Holmboe waves.
Hypersonic aerodynamics on thin bodies with interaction and upstream influence
Authors: Khorrami A. Farid , Smith Frank T. 
Published: April 2006

In the fundamental configuration studied here, a steady hypersonic free stream flows over a thin sharp aligned airfoil or flat plate with a leading-edge shock wave, and the flow field in the shock layer (containing a viscous and an inviscid layer) is steady laminar and two-dimensional, for a perfect gas without real and high-temperature gas effects. The viscous and inviscid layers are analysed and computed simultaneously in the region from the leading edge to the trailing edge, including the upstream-influence effect present, to determine the interactive flow throughout the shock layer and the positions of the shock wave and the boundary-layer edge, where matching is required. Further theoretical analysis of the shock layer helps to explain the computational results, including the nonlinear breakdown possible when forward marching against enhanced upstream influence, for example as the wall enthalpy increases towards its insulated value. Then the viscous layer is computed by sweeping methods, for higher values of wall enthalpies, to prevent this nonlinear breakdown for airfoils including the flat plate. Thin airfoils in hypersonic viscous flow are treated, for higher values of the wall enthalpies and with the upstream-influence effect, as are hypersonic inviscid flows, by modifying the computational methods used for the flat plate. Also, the behaviour of the upstream influence for bodies of relatively large thickness, and under wall velocity slip and enthalpy jump for flat plates, is discussed briefly from a theoretical point of view.Subsequent to the present work, computations based on the Navier–Stokes and on the parabolized Navier–Stokes equations have yielded excellent and good agreement respectively with the present predictions for large Mach and Reynolds numbers.
The onset of transverse recirculations during flow of gases in horizontal ducts with differentially heated lower walls
Authors: Ingle N. K. , Mountziaris T. J. 
Published: April 2006

A computational study has been performed to identify the onset of transverse buoyancy-driven recirculations during laminar flow of hydrogen and nitrogen in horizontal ducts with cool upper walls, and lower walls consisting of three sections: a cool upstream section, a heated middle section and a cool downstream section. The motivation for this work stems from the need to identify operating conditions maximizing the thickness uniformity, the interface abruptness and the precursor utilization during growth of thin films and multi-layer structures of semiconductors by metalorganic chemical vapour deposition (MOCVD). A mathematical model describing the flow and heat transfer along the vertical midplane of MOCVD reactors with the above geometry has been developed and computer simulations were performed for a variety of operating conditions using the Galerkin finite-element method. At atmospheric pressure and low inlet velocities, transverse recirculations form near the upstream and downstream edges of the heated section. These can be suppressed either by increasing the inlet velocity of the gas, so that forced convection dominates natural convection, or by decreasing the operating pressure to reduce the effects of buoyancy. The onset of transverse recirculations has been determined for Grashof (Gr) and Reynolds (Re) numbers covering the following ranges: 10−3 < Re < 100 and 1 < Gr < 106, with Gr and Re computed using fluid properties at the inlet conditions. The computations indicate that, for abrupt temperature changes along the lower wall (worst-case scenario), transverse recirculations are always absent if the following criteria are satisfied:
\[(Gr/Re) < 100\quad {\rm for}\quad 10^{-3} < Re \leqslant 4\quad {\rm and}\quad (Gr/Re^2) < 25\quad {\rm for}\quad 4 \leqslant Re < 100.\]The predicted critical values of Re, which correspond to the onset of transverse recirculations, agree well with reported experimental observations. The above criteria can be used for optimal design and operation of horizontal MOCVD reactors and may also be useful for heat transfer studies in horizontal ducts with differentially heated lower walls.