Showing papers in "arXiv: Fluid Dynamics in 2003"

Sensitive dependence on initial conditions in transition to turbulence in pipe flow

Abstract: The experiments by Darbyshire and Mullin (J. Fluid Mech. 289, 83 (1995)) on the transition to turbulence in pipe flow show that there is no sharp border between initial conditions that trigger turbulence and those that do not. We here relate this behaviour to the possibility that the transition to turbulence is connected with the formation of a chaotic saddle in the phase space of the system. We quantify a sensitive dependence on initial conditions and find in a statistical analysis that in the transition region the distribution of turbulent lifetimes follows an exponential law. The characteristic mean lifetime of the distribution increases rapidly with Reynolds number and becomes inaccessibly large for Reynolds numbers exceeding about 2200. Suitable experiments to further probe this concept are proposed.

Induced-Charge Electro-Osmosis

Abstract: We describe the general phenomenon of `induced-charge electro-osmosis' (ICEO) -- the nonlinear electro-osmotic slip that occurs when an applied field acts on the ionic charge it {\sl induces} around a polarizable surface. Motivated by a simple physical picture, we calculate ICEO flows around conducting cylinders in steady (DC), oscillatory (AC), and suddenly-applied electric fields. This picture, and these systems, represent perhaps the clearest example of nonlinear electrokinetic phenomena. We complement and verify this physically-motivated approach using a matched asymptotic expansion to the electrokinetic equations in the thin double-layer and low potential limits. ICEO slip velocities vary like $u_s \propto E_0^2 L$, where $E_0$ is the field strength and $L$ is a geometric length scale, and are set up on a time scale $\tau_c = \lambda_D L/D$, where $\lambda_D$ is the screening length and $D$ is the ionic diffusion constant. We propose and analyze ICEO microfluidic pumps and mixers that operate without moving parts under low applied potentials. Similar flows around metallic colloids with fixed total charge have been described in the Russian literature (largely unnoticed in the West). ICEO flows around conductors with fixed potential, on the other hand, have no colloidal analog and offer further possibilities for microfluidic applications.

Viscous heating effects in fluids with temperature-dependent viscosity: triggering of secondary flows

Abstract: Viscous heating can play an important role in the dynamics of fluids with strongly temperature-dependent viscosities because of the coupling between the energy and momentum equations. The heat generated by viscous friction produces a local temperature increase near the tube walls with a consequent decrease of the viscosity and a strong stratification in the viscosity profile. The problem of viscous heating in fluids was investigated and reviewed by Costa & Macedonio (2003) because of its important implications in the study of magma flows. Because of the strong coupling between viscosity and temperature, the temperature rise due to the viscous heating may trigger instabilities in the velocity field, which cannot be predicted by a simple isothermal Newtonian model. When viscous heating produces a pronounced peak in the temperature profile near the walls, a triggering of instabilities and a transition to secondary flows can occur because of the stratification in the viscosity profile. In this paper we focus on the thermal and mechanical effects caused by viscous heating. We will present the linear stability equations and we will show, as in certain regimes, these effects can trigger and sustain a particular class of secondary rotational flows which appear organised in coherent structures similar to roller vortices. This phenomenon can play a very important role in the dynamics of magma flows in conduits and lava flows in channels and, to our knowledge, it is the first time that it has been investigated by a direct numerical simulation.

A Bound on Mixing Efficiency for the Advection-Diffusion Equation

Abstract: An upper bound on the mixing efficiency is derived for a passive scalar under the influence of advection and diffusion with a body source. For a given stirring velocity field, the mixing efficiency is measured in terms of an equivalent diffusivity, which is the molecular diffusivity that would be required to achieve the same level of fluctuations in the scalar concentration in the absence of stirring, for the same source distribution. The bound on the equivalent diffusivity depends only on the functional "shape" of both the source and the advecting field. Direct numerical simulations performed for a simple advecting flow to test the bounds are reported.

Decay of the monochromatic capillary wave

Abstract: It was demonstrated by direct numerical simulation that, in the case of weakly nonlinear capillary waves, one can get resonant waves interaction on the discrete grid when resonant conditions are never fulfilled exactly. The waves's decay pattern was obtained. The influence of the mismatch of resonant condition was studied as well.

The Richardson's Law in Large-Eddy Simulations of Boundary Layer flows

Abstract: Relative dispersion in a neutrally stratified planetary boundary layer (PBL) is investigated by means of Large-Eddy Simulations (LES). Despite the small extension of the inertial range of scales in the simulated PBL, our Lagrangian statistics turns out to be compatible with the Richardson $t^3$ law for the average of square particle separation. This emerges from the application of nonstandard methods of analysis through which a precise measure of the Richardson constant was also possible. Its values is estimated as $C_2\sim 0.5$ in close agreement with recent experiments and three-dimensional direct numerical simulations.

Band gaps and localization of water waves over one-dimensional topographical bottoms

Abstract: In this paper, the phenomenon of band gaps and Anderson localization of water waves over one-dimensional periodic and random bottoms is investigated by the transfer matrix method. The results indicate that the range of localization in random bottoms can be coincident with the band gaps for the corresponding periodic bottoms. Inside the gap or localization regime, a collective behavior of water waves appears. The results are also compared with acoustic and optical situations.

A singular integrable equation from short capillary-gravity waves

Abstract: From a columnar approximation of the Euler equations of an incompressible fluid with surface tension, we derive in the short-wave approximation a new integrable classical 1+1 dimensional field theory for the motion of the surface. Together with a Lorentz invariance,this system has the novel feature of solutions which become multiple valued in finite time.

Quasi-geostrophic model of the instabilities of the Stewartson layer

Abstract: We study the destabilization of a shear layer, produced by differential rotation of a rotating axisymmetric container. For small forcing, this produces a shear layer, which has been studied by Stewartson and is almost invariant along the rotation axis. When the forcing increases, instabilities develop. To study the asymptotic regime (very low Ekman number $E$), we develop a quasi-geostrophic two-dimensional model, whose main original feature is to handle the mass conservation correctly, resulting in a divergent two-dimensional flow, and valid for any container provided that the top and bottom have finite slopes. We use it to derive scalings and asymptotic laws by a simple linear theory, extending the previous analyses to large slopes (as in a sphere), for which we find different scaling laws. For a flat container, the critical Rossby number for the onset of instability evolves as $E^{3/4}$ and may be understood as a Kelvin-Helmoltz shear instability. For a sloping container, the instability is a Rossby wave with a critical Rossby number proportional to $\beta E^{1/2}$, where $\beta$ is related to the slope. We also investigate the asymmetry between positive and negative differential rotation and propose corrections for finite Ekman and Rossby numbers. Implemented in a numerical code, our model allows us to study the onset over a broad range of parameters, determining the threshold but also other features such as the spatial structure. We also present a few experimental results, validating our model and showing its limits.

Effects of Bulk Viscosity in Non-linear Bubble Dynamics

Abstract: The non-linear bubble dynamics equations in a compressible liquid have been modified considering the effects of compressibility of both the liquid and the gas at the bubble interface. A new bubble boundary equation has been derived, which includes a new term resulted from the liquid bulk viscosity effects. The influence of this term has been numerically investigated considering the effects of water vapor and chemical reactions on the bubble evolution. The results clearly indicate that the new term has an important damping role at the collapse, so that its consideration decreases the amplitude of the bubble rebounds after the collapse. This damping feature is more remarkable for higher deriving pressures.

Acoustic radiation force and streaming induced by non-periodic ultrasound

Abstract: Starting from the linear flow of homogeneous fluid, five modes are defined as eigenvectors of the basic system of conservation laws. Quasi-plane geometry is considered. Projectors that separate overall perturbation of the fluid into specific modes are calculated and applied to nonlinear flow. Dynamic equations for the interacting modes are obtained in the frames of the method. A particular case of streaming caused by acoustic pulse is considered, illustrations on temporal behavior of streaming velocity and streamlines are presented.

Weak turbulence of gravity waves

Abstract: For the first time weak turbulent theory was demonstrated for the surface gravity waves. Direct numerical simulation of the dynamical equations shows Kolmogorov turbulent spectra as predicted by analytical analysis from kinetic equation.

The shell model approach to the rotating turbulence

Abstract: Applications of the shell model of turbulence to the case of rapidly rotating bodies are considered. Starting from the classical GOY model we introduce the Coriolis force and obtain a $\sim k^{-2}$ spectrum for 3D hydrodynamical turbulence for the free decay regime as well for the regime with external forcing. Additional modifications of the GOY model providing a realistic form of the helicity are proposed.

Most singular vortex structures in fully developed turbulence

Abstract: Using high Reynolds number experimental data, we search for most dissipative, most intense structures. These structures possess a scaling predicted by log-Poisson model for the dissipation field $\epsilon_r$. The probability distribution function for the exponents $\alpha$, $\epsilon_r\sim e^{\alpha a}$, has been constructed, and compared with Poisson distribution. These new experimental data suggest that the most intense structures have co-dimension less than 2. The log-Poisson statistics is compared with log-binomial which follows from the random $\beta$-model.

Analytical solution for nonlinear Schrodinger vortex reconnection

Abstract: Analysis of the nonlinear Schrodinger vortex reconnection is given in terms of coordinate-time power series. The lowest order terms in these series correspond to a solution of the linear Schrodinger equation and provide several interesting properties of the reconnection process, in particular the non-singular character of reconnections, the anti-parallel configuration of vortex filaments and a square-root law of approach just before/after reconnections. The complete infinite power series represents a fully nonlinear analytical solution in a finite volume which includes the reconnection point, and is valid for finite time provided the initial condition is an analytic function. These series solutions are free from the periodicity artifacts and discretization error of the direct computational approaches and they are easy to analyze using a computer algebra program.

The Flowing System Gasdynamics. Part 3: Saint-Venant - Wantzel formula modern form

Abstract: The modern form of the Saint-Venant - Wantzel formula for an outflow velocity of gas stream from flowing element is submitted Taking into account of contact interaction of gas stream with the streamline surface in the form of the static head law has allowed to find the spatial-energy liaison in flowing system The physically correct combination of mechanics of contact interaction and thermodynamics of fluid medium in one formula has allowed simultaneously to be liberated from the velocity coefficient and the discharge coefficient and polytropic process In the new form the formula has gained the key character for computation of parameters of motion and state of gas stream in the flowing system

A model differential equation for turbulence

Abstract: A phenomenological turbulence model in which the energy spectrum obeys a nonlinear diffusion equation is presented. This equation respects the scaling properties of the original Navier-Stokes equations and it has the Kolmogorov -5/3 cascade and the thermodynamic equilibrium spectra as exact steady state solutions. The general steady state in this model contains a nonlinear mixture of the constant-flux and thermodynamic components. Such "warm cascade" solutions describe the bottleneck phenomenon of spectrum stagnation near the dissipative scale. Self-similar solutions describing a finite-time formation of steady cascades are analysed and found to exhibit nontrivial scaling behaviour.

Using Drag to Hover

Abstract: Unlike a helicopter, an insect can, in theory, use both lift and drag to stay aloft. Here we show that a dragonfly uses mostly drag to hover by employing asymmetric up and down strokes. Computations of a family of strokes further show that using drag can be as efficient as using lift at the low Reynolds number regime appropriate for insects.

The gas equation for stream

Abstract: The equation of state for gas stream as analog of well-known equation of state of gas medium in the static conditions is submitted. The connection of parameters of state and motion, and also the volume-flow rate and the mass-flow rate of the gas stream during change of its temperature and pressure as the result of external action is shown on example of modification of the represented equation.

Transition to turbulence in strongly heated vertical natural convection boundary layers

Abstract: The mechanisms governing the transition to turbulence in natural convection boundary layers along strongly heated vertical walls remain neither very clear nor well understood, because of the lack of experiments and the difficulties of physical modelling. Our measurements bring experimental data focusing on this transition, in quiescent air along radiating and conducting plates, in the whole range of 2000 to 8000 W.m -2 heating rate. The analysis of the time series obtained by sliding window cross-correlation thermoanemometry leads us to point out coherent turbulent structures on short heights throughout the thin boundary layer, which seem to be governed by heat transfer and time microscales of turbulence through the inner sublayer. Physical interpretations are given to relate the observed heat transfer correlation and these turbulence transition structures along with radiation and conduction.

The Flowing System Gasdynamics Part 2: Euler momentum conservation equation solution

Abstract: The solution of a momentum conservation equation for the gas and liquid stream in the flowing element is obtained on the basis of the modern approach to a problem on contact interaction of bodies and mediums. A flowing element, system are: pipe, tube, orifice, mouthpiece, diffuser, etc. and its combination. The integration of the differential equation has reduced to distribution law of static head along the length of flowing element and has proved the elementary algebraic solution that introduced in the previous paper by these authors. The received solution allows to describe the motion of fluid medium in non-stationary conditions, under action of any time-varying physical factors: a roughness of streamlined surface, the area of the section of the flowing element, the heat exchange with the streamlined surface, the technical work, the additional weight flow of fluid medium.

The laminar flow instability criterion and turbulence in pipe

Abstract: The identification of stream in the straight pipe as a flexible rod has allowed to present the criterion expression for determination of transition of the laminar flow regime to the turbulent as a loss of stability of the rectilinear static structure of translational motion of stream in pipe and its transition to the flexural-vortical dynamic structure of translational motion, just as a flexible rod buckling. The introduced criterion allows to take into account an influencing of the inlet geometry, the pipe length, the flow velocity, and also of any physical factors on stability of the rectilinear flow structure. It is ascertained, that Reynolds number is the number of local hydrodynamic similarity and it is displayed that one is constituent part of the introduced stability criterion. The developed approach to a problem of stability is applicable for a problem solving on internal flow and external streamline.

A model for rapid stochastic distortions of small-scale turbulence

Abstract: We present a model describing evolution of the small-scale Navier-Stokes turbulence due to its stochastic distortions by much larger turbulent scales. This study is motivated by numerical findings (laval, 2001) that such interactions of separated scales play important role in turbulence intermittency. We introduce description of turbulence in terms of the moments of the k-space quantities using a method previously developed for the kinematic dynamo problem (Nazarenko, 2003). Working with the $k$-space moments allows to introduce new useful measures of intermittency such as the mean polarization and the spectral flatness. Our study of the 2D turbulence shows that the energy cascade is scale invariant and Gaussian whereas the enstrophy cascade is intermittent. In 3D, we show that the statistics of turbulence wavepackets deviates from gaussianity toward dominance of the plane polarizations. Such turbulence is formed by ellipsoids in the $k$-space centered at its origin and having one large, one neutral and one small axes with the velocity field pointing parallel to the smallest axis.

The Flowing System Gasdynamics Part 1: On static head in the pipe flowing element

Abstract: The solution of problem on distribution of static head along the pipe flowing element is submitted. The solution is reached on the basis of consideration of contact interaction of gas and liquid stream with the wall of the pipe flowing element. The expression for distribution of static head along the pipe flowing element is obtained on the basis of usage of three fundamental laws in fluid dynamics: Torricelli formula, Weissbach-Darcy formula and Bernoulli equation. The general solution is obtained for a gas stream. The special case of the obtained solution is retrieved for liquid stream.

The boundary and continual transfer phenomena in fluids and flows

Abstract: The clearing up of a wave nature of the energy and mass transfer phenomena in classical expressions of the molecular-kinetic theory has allowed to find a quantitative measure of intensity of processes of a thermal conductivity, viscosity and diffusion in conditions of a thermally nonequilibrium and heterogeneous composition continuum It is rotined that the appearance of a temperature drop in fluid stipulates the appearance of the continuum stratification and formation of the flowing bodies interacting among themselves It is rotined that the known expressions for a thermal conductivity, viscosity, diffusion and heat convection had been obtained for a thermally equilibrium and homogeneous continuum and produce a maximum quantity of intensity of the transfer processes The introduced expression is usable for a quantitative estimation of intensity of the transfer processes in fluid and its streams in conditions of the non-stationary heat exchange in natural conditions and technical problems

A New Damping Mechanism in Non-linear Bubble Dynamics

Abstract: Non-linear equations of radial motion of a gas bubble in a compressible viscous liquid have been modified considering effects of viscosity and compressibility more complete than all previous works. A new set of equations has been derived including new terms resulted from consideration of the viscosity and compressibility not only at the bubble interface, but also in the bulk of liquid. The new equations are two non-linear coupled equations, which can not be merged into one equation unlike all previously derived equations. Numerical calculations have been performed considering effects of heat and mass transfer at the bubble interface. The results indicate that the new terms exhibit an important damping role at the collapse, so that their consideration dramatically weakens the bubble rebounds after the collapse. Dependence of this new damping mechanism to amplitude and frequency of the deriving pressure has been investigated.

Twisting Rolls. An Heuristic Model and 3D Numerical Simulations of Vortex Patterns

Abstract: We connect an appropriate feedback loop to a model of 2D vertical eddy of airflow which unfolds a wide range of vorticity behavior. Computational fluid dynamics of the twisted roll display a class of long lifespan 3D vortices. On the one hand, the infinitely stable columnar vortex simulated describes waterspouts and tornadoes with extended lifetime. On the other hand, a light modification of the retroaction exhibits strong similarities to tropical cyclones. Moreover, we investigate the outcome of the twisting process vertically shifted. This modelisation leads to the simulation of simultaneous vortices associated to this other class of 3D vortices with short lifespan. Our heuristic dynamical systems lay the foundations of a comprehensive modelisation of vortices since it joins theory and numerical simulations.

Dynamic Point-Formation in Dielectric Fluids

Abstract: We use boundary-integral methods to compute the time-dependent deformation of a drop of dielectric fluid immersed in another dielectric fluid in a uniform electric field E. Steady state theory predicts, when the permittivity ratio, \beta, is large enough, a conical interface can exist at two cone angles, with \theta_ (\beta) unstable. Our numerical evidence instead shows a dynamical process which produces a cone-formation and a transient finite-time singularity, when E and \beta are above their critical values. Based on a scaling analysis of the electric stress and the fluid motion, we are able to apply approximate boundary conditions to compute the evolution of the tip region. We find in our non-equilibrium case where the electric stress is substantially larger than the surface tension, the ratio of the electric stress to the surface tension in the newly-grown cone region can converge to a \beta dependent value, \alpha_c(\beta)>1, while the cone angle converges to \theta_<(\beta). This new dynamical solution is self-similar.

Message Passing Fluids: molecules as processes in parallel computational fluids

Abstract: In this paper we present the concept of MPF, Message Passing Fluid, an abstract fluid where the molecules move by mean of the informations that they exchange each other, on the basis of rules and methods of a generalized Cellular Automaton. The model is intended for its simulation by mean of message passing libraries on the field of parallel computing. We present a critical analysis of the necessary computational effort in a possible implementation of such an object.

Quasi-1D Modeling of Polymer Melt Die Swell in Short Dies

Abstract: This paper describes the isothermal die swell using our recent quasi-1D model for fast (high Deborah number) contraction flows of polymers melts. Because the model analyzes the flow in several flow regions as one continuous process, it makes possible to evaluate the die swell as a qusi-1D extrudate flow in dies of various lengths. Using the asymptotic matching condition for the change in flow type at the die exit allowed us to find the swelling profile for extrudate along the flow direction. The calculations in paper performed using a multi-mode viscoelastic constitutive equation of differential type, are compared with the experimental/direct numerical data including basic rheological tests. The presented swelling model involves no fitting parameter and is applicable for calculations using any viscoelastic constitutive equation.