Showing papers in "arXiv: Fluid Dynamics in 2001"

Measurement of Particle Accelerations in Fully Developed Turbulence

Abstract: We use silicon strip detectors (originally developed for the CLEO III high energy particle physics experiment) to measure fluid particle trajectories in turbulence with temporal resolution of up to 70,000 frames per second. This high frame rate allows the Kolmogorov time scale of a turbulent water flow to be fully resolved for 140 = 500. The acceleration flatness is found to increase with Reynolds number, and to exceed 60 at R_lambda = 970. The coupling of the acceleration to the large scale anisotropy is found to be large at low Reynolds number and to decrease as the Reynolds number increases, but to persist at all Reynolds numbers measured. The dependence of the acceleration variance on the size and density of the tracer particles is measured. The autocorrelation function of an acceleration component is measured, and is found to scale with the Kolmogorov time tau_eta.

Equations relating structure functions of all orders

Abstract: The hierarchy of exact equations is given that relates two-spatial-point velocity structure functions of arbitrary order with other statistics. Because no assumption is used, the exact statistical equations can apply to any flow for which the Navier-Stokes equations are accurate, and they apply no matter how small the number of samples in the ensemble. The exact statistical equations can be used to verify DNS computations and to detect their limitations. For example,if DNS data are used to evaluate the exact statistical equations, then the equations should balance to within numerical precision, otherwise a computational problem is indicated. The equations allow quantification of the approach to local homogeneity and to local isotropy. Testing the balance of the equations allows detection of scaling ranges for quantification of scaling-range exponents. The second-order equations lead to Kolmogorov's equation. All higher-order equations contain a statistic composed of one factor of the two-point difference of the pressure gradient multiplied by factors of velocity difference. Investigation of this pressure-gradient-difference statistic can reveal much about two issues: 1) whether or not different components of the velocity structure function of given order have differing exponents in the inertial range, and 2) the increasing deviation of those exponents from Kolmogorov scaling as the order increases. Full disclosure of the mathematical methods is in xxx.lanl.gov/list/physics.flu-dyn/0102055.

Non-locality and Intermittency in 3D Turbulence

Abstract: Numerical simulations are used to determine the influence of the non-local and local interactions on the intermittency corrections in the scaling properties of 3D turbulence. We show that neglect of local interactions leads to an enhanced small-scale energy spectrum and to a significantly larger number of very intense vortices (tornadoes) and stronger intermittency. On the other hand, neglect of the non-local interactions results in even stronger small-scale spectrum but significantly weaker intermittency. Based on these observations, a new model of turbulence is proposed, in which non-local (RDT-like) interactions couple large and small scale via a multiplicative process with additive noise and the local interactions are modeled by a turbulent viscosity. This model is used to derive a simple toy version of the Langevin equations for small-scale velocity increments. A Gaussian approximation for the large scale fields yields the Fokker-Planck equation for the probability distribution function of the velocity increments. Steady state solutions of this equation allows to qualitatively explain the anomalous corrections and the skewness generation along scale.

Gravitational oscillations of a liquid column

Abstract: We report gravity oscillations of a liquid column partially immersed in a bath of liquid. We stress in particular some peculiarities of this system, namely (i) the fact that the mass of this oscillator constantly changes with time; (ii) the singular character of the beginning of the rise, for which the mass of the oscillator is zero; (iii) the sources of dissipation in this system, which is found to be dominated at low viscosity by the entrance (or exit) effects, leading to a long-range damping of the oscillations. We conclude with some qualitative description of a second-order phenomenon, namely the eruption of a jet at the beginning of the rise.

A characteristic frequency of two mutually interacting gas bubbles in an acoustic field

Abstract: Transition frequencies of two spherical gas bubbles interacting in an acoustic field are discussed theoretically. In the present study, transition frequency is defined as the frequency of external sound for which the phase difference between a bubble's pulsation and the external sound is $\pi / 2$. It is shown by a linear theory that a bubble interacting with a neighboring bubble has three (or fewer) transition frequencies but only two natural frequencies. This result means that the bubble has a characteristic frequency besides the natural frequencies.

Passive scalar anisotropy in a heated turbulent wake: new observations and implications for LES

Abstract: The effects of passive scalar anisotropy on subgrid-scale (SGS) physics and modeling for Large-Eddy Simulations are studied experimentally. Measurements are performed across a moderate Reynolds number wake flow generated by a heated cylinder, using an array of four X-wire and four cold-wire probes. By varying the separation distance among probes in the array, we obtain filtered and subgrid quantities at three different filter sizes. We compute several terms that comprise the subgrid dissipation tensor of kinetic energy and scalar-variance and test for isotropic behavior, as a function of filter scale. We find that whereas the kinetic energy dissipation tensor tends towards isotropy at small scales, the SGS scalar-variance dissipation remains anisotropic independent of filter scale. The eddy-diffusion model predicts isotropic behavior, whereas the nonlinear (or tensor eddy diffusivity) model reproduces the correct trends, but overestimates the level of scalar dissipation anisotropy. These results provide some support for so-called mixed models but raise new questions about the causes of the observed anisotropy.

Time resolved tracking of a sound scatterer in a turbulent flow: non-stationary signal analysis and applications

Abstract: It is known that ultrasound techniques yield non-intrusive measurements of hydrodynamic flows. For example, the study of the echoes produced by a large number of particle insonified by pulsed wavetrains has led to a now standard velocimetry technique. In this paper, we propose to extend the method to the continuous tracking of one single particle embedded in a complex flow. This gives a Lagrangian measurement of the fluid motion, which is of importance in mixing and turbulence studies. The method relies on the ability to resolve in time the Doppler shift of the sound scattered by the continuously insonfied particle. For this signal processing problem two classes of approaches are used: time-frequency analysis and parametric high resolution methods. In the first class we consider the spectrogram and reassigned spectrogram, and we apply it to detect the motion of a small bead settling in a fluid at rest. In more non-stationary turbulent flows where methods in the second class are more robust, we have adapted an Approximated Maximum Likelihood technique coupled with a generalized Kalman filter.

Is the Critical Reynolds Number Universal

Abstract: The laminar-turbulent transition in the circular pipe flow has been tested experimentally. The critical Reynolds numbers for the flows of different gases (He, Ne, Ar, Kr, Xe, N2, CO2, SF6) and liquids (H2O, D2O, C2H5OH) have been compared. The difference up to 40% was observed. The possible reasons of non-universality of critical Reynolds number are discussed. The conclusion is the statistical approach is needed to penetrate into the nature of observed phenomenon.

Topology and Turbulence

Abstract: Over a given regular domain of independent variables {x,y,z,t}, every covariant vector field of flow can be constructed in terms a differential 1-form of Action. The associated Cartan topology permits the definition of four basic topological equivalence classes of flows based on the Pfaff dimension of the 1-form of Action. Potential flows or streamline processes are generated by an Action 1-form of Pfaff dimension 1 and 2, respectively. Chaotic flows must be associated with domains of Pfaff dimension 3 or more. Turbulent flows are associated with domains of Pfaff dimension 4. It will be demonstrated that the Navier-Stokes equations are related to Action 1-forms of Pfaff dimension 4. The Cartan Topology is a disconnected topology if the Pfaff dimension is greater than 2. This fact implies that the creation of turbulence (a state of Pfaff dimension 4 and a disconnected Cartan topology) from a streamline flow (a state of Pfaff dimension 2 and a connected topology) can take place only by discontinuous processes which induce shocks and tangential discontinuities. On the otherhand, the decay of turbulence can be described by continuous, but irreversible, processes. Numerical procedures that force continuity of slope and value cannot in principle describe the creation of turbulence, but such techniques of forced continuity can be used to describe the decay of turbulence.

Interaction of two deformable viscous drops under external temperature gradient

Abstract: The axisymmetric deformation and motion of interacting droplets in an imposed temperature gradient is considered using boundary-integral techniques for slow viscous motion. Results showing temporal drop motion, deformations and separation are presented for equal-viscosity fluids. The focus is on cases when the drops are of equal radii or when the smaller drop trails behind the larger drop. For equal-size drops, our analysis shows that the motion of a leading drop is retarded while the motion of the trailing one is enchanced compared to the undeformable case. The distance between the centers of equal-sized deformable drops decreases with time. When a small drop follows a large one, two patterns of behavior may exist. For moderate or large initial separation the drops separate. However, if the initial separation is small there is a transient period in which the separation distance initially decreases and only afterwards the drops separate. This behavior stems from the multiple time scales that exist in the system.

Critical curves of plane Poiseuille flow with slip boundary conditions

Abstract: We investigate the linear stability of plane Poiseuille flow in 2D under slip boundary conditions. The slip s is defined by the tangential velocity at the wall in units of the maximal flow velocity. As it turns out, the critical Reynolds number depends smoothly on s but increases quite rapidly.

A route to magnetic field reversals: an example of an ABC-forced non-linear dynamo

Abstract: We are investigating numerically the non-linear behavior of a space-periodic MHD system with ABC forcing. Most computations are performed for magnetic Reynolds numbers increasing from 0 to 60 and a fixed kinematic Reynolds number, small enough for the trivial solution with a zero magnetic field to be stable to velocity perturbations. At the critical magnetic Reynolds number for the onset of instability of the trivial solution the dominant eigenvalue of the kinematic dynamo problem is real. In agreement with the bifurcation theory new steady states with non-vanishing magnetic field appear in this bifurcation. Subsequent bifurcations are investigated. A regime is detected, where chaotic variations of the magnetic field orientation (analogous to magnetic field reversals) are observed in the temporal evolution of the system.

Alternative to R_lambda-scaling of Small-Scale Turbulence Statistics

Abstract: Traditionally, trends of universal turbulence statistics are presented versus R-lambda, which is the Reynolds number based on Taylor's scale, lambda, and the root-mean-squared (rms) velocity component, u'. Taylor's scale and u', and hence R-lambda, do not have the attribute of universality. The ratio of rms fluid-particle acceleration to rms viscous acceleration is an alternative to R-lambda which has the advantage of being determined by the small scales of turbulence. This ratio, denoted R-a, has the following attributes: it is a Reynolds number, it is composed of statistics of the small scales of turbulence, it can be evaluated with single-wire hot-wire anemometry, and like R-lambda, can be partially evaluated by means of flow similarity. For isotropic turbulence the relationship between R-a and R-lambda is shown. Velocity derivative flatness measured between counter-rotating blades appears monotonic when graphed versus R-a, unlike when graphed versus R-lambda.

Nonlinear Theory of Flame Front Instability

Abstract: Nonlinear non-stationary equation describing evolution of weakly curved premixed flames with arbitrary gas expansion, subject to the Landau-Darrieus instability, is derived. The new equation respects all the conservation laws to be satisfied across the flame front, as well as correctly takes into account influence of vorticity, generated in the flame, on the flame front structure and flame velocity. Analytical solutions of the derived equation are found.

Some closed form solutions to the Navier Stokes equations

Abstract: An algorithm for generating a class of closed form solutions to the Navier-Stokes equations is suggested, with examples. Of particular interest are those exact solutions that exhibit intermittency, tertiary Hopf bifurcations, flow reversal, and hysteresis.

Particles, fluids and vortices

Abstract: Classical particle mechanics on curved spaces is related to the flow of ideal fluids, by a dual interpretation of the Hamilton-Jacobi equation. As in second quantization, the procedure relates the description of a system with a finite number of degrees of freedom to one with infinitely many degrees of freedom. In some two-dimensional fluid mechanics models a duality transformation between the velocity potential and the stream function can be performed relating sources and sinks in one model to vortices in the other. The particle mechanics counterpart of the dual theory is reconstructed. In the quantum theory the strength of sources and sinks, as well as vorticity are quantized; for the duality between theories to be preserved these quantization conditions must be related.

Edge waves along a sloping beach

Abstract: We construct a family of explicit rotational solutions to the nonlinear governing equations for water waves, describing edge waves propagating over a plane-sloping beach. A detailed analysis of the edge wave dynamics and of the run-up pattern is made possible by the use of the Lagrangian approach to the water motion. A graphical representation of the edge wave is also presented.

Numerical evidence of breaking of vortex lines in an ideal fluid

Abstract: Emergence of singularity of vorticity at a single point, not related to any symmetry of the initial distribution, has been demonstrated numerically for the first time. Behavior of the maximum of vorticity near the point of collapse closely follows the dependence 1/(t0-t), where t0 is the time of collapse. This agrees with the interpretation of collapse in an ideal incompressible fluid as of the process of vortex lines breaking.

The Inverse Energy Cascade of Two-Dimensional Turbulence

Abstract: This thesis presents an experimental study of the inverse energy cascade as it occurs in an electromagnetically forced soap film It focuses on characterizing important features of the inverse cascade such as it's range, how energy is distributed over the range and how energy flows through the range The thesis also probes the assumption of scale invariance that is associated with the existence of an inverse cascade These investigations demonstrate that the extent of the inverse cascade range and the behavior of the energy distribution are in agreement with dimensional predictions The energy flow in the inverse cascade range is shown to be well described by exact mathematical predictions obtained from the Navier-Stokes equation At no time does the energy flow in the inverse cascade range produced by the e-m cell behave inertially or in a scale invariant manner Evidence that the cascade could become scale invariant should an inertial range develop is presented, as are the requirements that a system must satisfy to create such an inertial range

Remarks on unsolved basic problems of the Navier-Stokes equations

Abstract: There is renewed interest in the question of whether the Navier-Stokes equations (NSE), one of the fundamental models of classical physics and widely used in engineering applications, are actually self-consistent. After recalling the essential physical assumptions inherent in the NSE, the notion of weak solutions, possible implications for the energy conservation law, as well as existence and uniqueness in the incompressible case are discussed. Emphasis will be placed on the possibility of finite time singularities and their consequences for length scales which should be consistent with the continuum hypothesis.

Matrix Fluid Dynamics

Abstract: The paper reports the recent results on application and extension of the matrix formulation of lagrangian hydrodynamic equations. The matrix approach is based on the notion of continuous deformation of infinitesimal material elements and treats the Jacobi matrix of their derivatives with respect to lagrangian variables as the fundamental quantity completely describing fluid motion. We begin with brief review of the governing matrix equation (for the detailed discussion see Yakubovich, E.I. & Zenkovich, D.A. 2001 Matrix approach to Lagrangian fluid dynamics. J. Fluid Mech. vol. 443, 167-196). The new general relationship between two- and three-dimensional solutions of the matrix equations is then found and applied to 3-D generalization of plane 'Ptolemaic' vortices. As a result, families of non-stationary swirling vortices and of stretched vortices in an axisymmetric strain are constructed and studied. The set of matrix equations is extended to include the effect of viscosity by incorporating an evolution equation for the Cauchy invariants, which represents a lagrangian analogue of the Helmholtz equation for a viscous fluid.

Scaling of acceleration in locally isotropic turbulence

Abstract: The variances of the fluid-particle acceleration and of the pressure-gradient and viscous force are given. The scaling parameters for these variances are velocity statistics measureable with a single-wire anemometer. For both high and low Reynolds numbers, asymptotic scaling formulas are given; these agree quantitatively with DNS data. Thus, the scaling can be presumed known for all Reynolds numbers. Fluid-particle acceleration variance does not obey K41 scaling at any Reynolds number; this is consistent with recent experimental data. The non-dimensional pressure-gradient variance named lambda-sub{T} /lambda-sub{P} is shown to be obsolete.

On the thermocapillary motion of deformable droplets

Abstract: In studies on Marangoni type motion of particles the surface tension is often approximated as a linear function of temperature. For deformable particles in a linear external temperature gradient far from the reference point this approximation yields a negative surface tension which is physically unrealistic. It is shown that H.Zhou and R.H.Davis J. Colloid Interface Sci., n.181,60,1996 presented calculation where the leading deformable drop moved into a region of negative surface tension. With respect numerical studies the restriction of the migration of two deformable drops is given in terms of the drift time.

On the small-scale statistics of Lagrangian turbulence

Abstract: We provide evidence that the small-scale statistics of the acceleration of a test particle in high-Reynolds number Lagrangian turbulence is correctly described by Tsallis statistics with entropic index q=3/2. We present theoretical arguments why Tsallis statistics can naturally arise in Lagrangian turbulence and why at the smallest scales q=3/2 is relevant. A generalized Heisenberg-Yaglom formula is derived from the nonextensive model.

The Kraichnan-Kazantsev dynamo

Abstract: The problem of the dynamo effect for a Kraichnan incompressible helicity-free velocity field is considered. Exploiting a quantum formalism first introduced by Kazantsev (A.P. Kazantsev, Sov. Phys. JETP 26, 1031-1034 (1968)), we show that a critical magnetic Reynolds number exists for the presence of dynamo. The value of the Prandtl number influences the spatial distribution of the magnetic field and its growth in time. The magnetic field correlation length is always the largest between the diffusive scale and the viscous scale of the flow. In the same way the field growth is characterized by a time scale that corresponds to the largest between the diffusive and the viscous characteristic time.

Global stability of systems related to the Navier-Stokes equations

Abstract: A generalized Lyapunov method is outlined which predicts global stability of a broad class of dissipative dynamical systems. The method is applied to the complex Lorenz model and to the Navier-Stokes equations. In both cases one finds compact domains in phase space which contain the omega sets of all trajectories, in particular the fixed points, limit cycles, and strange attractors.

Temperature behavior of vortices of a 3D thermoconducting viscous fluid

Abstract: The Navier-Stokes-Fourier model for a 3D thermoconducting viscous fluid, where the evolution equation for the temperature T contains a term proportional to the rate of energy dissipation, is investigated analitically at the light of the rotational invariance property. Two cases are considered: the Couette flow and a flow with a radial velocity between two rotating impermeable and porous coaxial cylinders, respectively. In both cases, we show the existence of a maximum value of T, T_max, when the difference of temperature Delta T=T_2-T_1 on the surfaces of the cylinders is assigned. The role of T_max is discussed in the context of different physical situations.

Mathematics of Structure-Function Equations of All Orders

Abstract: Exact equations are derived that relate velocity structure functions of arbitrary order with other statistics. "Exact" means that no approximations are used except that the Navier-Stokes equation and incompressibility condition are assumed to be accurate. The exact equations are used to determine the structure-function equations of all orders for locally homogeneous but anisotropic turbulence as well as for the locally isotropic case. These equations can be used for investigating the approach to local homogeneity and to local isotropy as well as the balance of the equations and identification of scaling ranges.

Searching for vortex structures in high Reynolds number turbulence

Abstract: In experimental study of very high Reynolds number turbulence, we found evidences that there are distinguished vortex structures in the intermediate range, that is, between the Kolmogorov and Taylor microscales, where they are indeed expected to be present These structures are responsible for the intermittency, and, in the same time, they contribute into asymmetry of turbulent statistics, the latter following from the Kolmogorov law

Relativistic Mechanics of Continuous Media

Abstract: In this work we study the relativistic mechanics of continuous media on a fundamental level using a manifestly covariant proper time procedure. We formulate equations of motion and continuity (and constitutive equations) that are the starting point for any calculations regarding continuous media. In the force free limit, the standard relativistic equations are regained, so that these equations can be regarded as a generalization of the standard procedure. In the case of an inviscid fluid we derive an analogue of the Bernoulli equations. For irrotational flow we prove that the velocity field can be derived from a potential. If, in addition, the fluid is incompressible, the potential must obey the d'Alembert equation, and thus the problem is reduced to solving the d'Alembert equation with specific boundary conditions (in both space and time). The solutions indicate the existence of light velocity sound waves in an incompressible fluid (a result known from previous literature [19]). Relaxing the constraints and allowing the fluid to become linearly compressible, one can derive a wave equation from which the sound velosity can again be computed. For a stationary background flow, it has been demonstrated that the sound velocity attains its corrrect values for the incompressible and non-relatvistic limits. Finally, viscosity is introduced, bulk and shear viscosity constants are defined, and we formulate equations for the motion of a viscous fluid.