Showing papers in "arXiv: Fluid Dynamics in 2005"

Retraction dynamics of aquous drops upon impact on nonwetting surfaces

Abstract: We study the impact and subsequent retraction dynamics of liquid droplets upon high-speed impact on hydrophobic surfaces. Performing extensive experiments, we show that the drop retraction rate is a material constant and does not depend on the impact velocity. We show that when increasing the Ohnesorge number, $\Oh=\eta/\sqrt{\rho R_{\rm I} \gamma}$, the retraction, i.e. dewetting, dynamics crosses over from a capillaro-inertial regime to a capillaro-viscous regime. We rationalize the experimental observations by a simple but robust semi-quantitative model for the solid-liquid contact line dynamics inspired by the standard theories for thin film dewetting.

Heat transport in turbulent Rayleigh-Benard convection: Effect of finite top- and bottom-plate conductivity

Abstract: We describe three apparatus, known as the large, medium, and small apparatus, used for high-precision measurements of the Nusselt number N as a function of the Rayleigh number R for cylindrical samples of fluid and present results illustrating the influence of the finite conductivity of the top and bottom plates on the heat transport in the fluid. We used water samples at a mean temperature of 40 degrees C (Prandtl number sigma = 4.4). The samples in the large apparatus had a diameter D of 49.69 cm and heights L = 116.33, 74.42, 50.61, and 16.52 cm. For the medium apparatus we had D = 24.81 cm, and L = 90.20 and 24.76 cm. The small apparatus contained a sample with D = 9.21 cm, and L = 9.52 cm. For each aspect ratio Gamma = D/L the data covered a range of a little over a decade of R. The maximum R = 10^12 with Nusselt numbers N = 600 was reached for Gamma = 0.43. Measurements were made with both Aluminum and Copper top and bottom plates of nominally identical size and shape. For the large and medium apparatus the results with Aluminum plates fall below those obtained with Copper plates, thus confirming qualitatively the prediction by Verzicco that plates of finite conductivity diminish the heat transport in the fluid. The Nusselt number N_infinity for plates with infinite conductivity was estimated by fitting simultaneously Aluminum- and Copper-plate data sets to an effective powerlaw for N_infinity multiplied by a correction factor f(X) = 1 - exp[-(aX)^b] that depends on the ratio X of the thermal resistance of the fluid to that of the plates as suggested by Verzicco. Within their uncertainties the parameters a and b were independent of Gamma for the large apparatus and showed a small Gamma-dependence for the medium apparatus. The correction was larger for the large, smaller for the medium, and negligible for the small apparatus.

Viscous fingering of miscible slices

Abstract: Viscous fingering of a miscible high viscosity slice of fluid displaced by a lower viscosity fluid is studied in porous media by direct numerical simulations of Darcy's law coupled to the evolution equation for the concentration of a solute controlling the viscosity of miscible solutions. In contrast with fingering between two semi-infinite regions, fingering of finite slices is a transient phenomenon due to the decrease in time of the viscosity ratio across the interface induced by fingering and dispersion processes. We show that fingering contributes transiently to the broadening of the peak in time by increasing its variance. A quantitative analysis of the asymptotic contribution of fingering to this variance is conducted as a function of the four relevant parameters of the problem i.e. the log-mobility ratio R, the length of the slice l, the Peclet number Pe and the ratio between transverse and axial dispersion coefficients $\epsilon$. Relevance of the results is discussed in relation with transport of viscous samples in chromatographic columns and propagation of contaminants in porous media.

On the origin of non-Gaussian statistics in hydrodynamic turbulence

Abstract: Turbulent flows are notoriously difficult to describe and understand based on first principles One reason is that turbulence contains highly intermittent bursts of vorticity and strain-rate with highly non-Gaussian statistics Quantitatively, intermittency is manifested in highly elongated tails in the probability density functions of the velocity increments between pairs of points A long-standing open issue has been to predict the origins of intermittency and non-Gaussian statistics from the Navier-Stokes equations Here we derive, from the Navier-Stokes equations, a simple nonlinear dynamical system for the Lagrangian evolution of longitudinal and transverse velocity increments From this system we are able to show that the ubiquitous non-Gaussian tails in turbulence have their origin in the inherent self-amplification of longitudinal velocity increments, and cross amplification of the transverse velocity increments

Hydrodynamic interactions of spherical particles in suspensions confined between two planar walls

Abstract: Hydrodynamic interactions in a suspension of spherical particles confined between two parallel planar walls are studied under creeping-flow conditions. The many-particle friction matrix in this system is evaluated using our novel numerical algorithm based on transformations between Cartesian and spherical representations of Stokes flow. The Cartesian representation is used to describe the interaction of the fluid with the walls and the spherical representation is used to describe the interaction with the particles. The transformations between these two representations are given in a closed form, which allows us to evaluate the coefficients in linear equations for the induced-force multipoles on particle surfaces. The friction matrix is obtained from these equations, supplemented with the superposition lubrication corrections. We have used our algorithm to evaluate the friction matrix for a single sphere, a pair of spheres, and for linear chains of spheres. The friction matrix exhibits a crossover from a quasi-two-dimensional behavior (for systems with small wall separation H) to the three-dimensional behavior (when the distance H is much larger than the interparticle distance L). The crossover is especially pronounced for a long chain moving in the direction normal to its orientation and parallel to the walls. In this configuration, a large pressure buildup occurs in front of the chain for small values of the gapwidth H, which results in a large hydrodynamic friction force. A standard wall superposition approximation does not capture this behavior.

Absense of slow transients, and the effect of imperfect vertical alignment, in turbulent Rayleigh-Benard convection

Abstract: We report experimental results for the influence of a tilt angle beta relative to gravity on turbulent Rayleigh-Benard convection of cylindrical samples. The measurements were made at Rayleigh numbers R up to 10^11 with two samples of height L equal to the diameter D (aspect ratio Gamma = D/L = 1). The fluid was water with a Prandtl number sigma = 4.38. In contrast to the experiences reported by Chilla et. al. (2004) for a similar sample but with Gamma = 0.5 (D = 0.5 and L = 1.0 m), we found no long relaxation times. For R = 9.4 times 10^10 we measured the Nusselt number N as a function of tilt angle beta and obtained a small beta dependence about a factor of 50 smaller than the result found by Chilla et. al. (2004) for their Gamma = 0.5 sample. We measured side-wall temperatures at eight equally spaced azimuthal locations on the horizontal mid-plane of the sample and used their cross-correlation functions to find the turn-over time of the large-scale circulation (LSC). The resulting Reynolds numbers R_e^cc were found to increase with beta. An important conclusion is that the increase of R_e^cc with beta of the LSC does not significantly influence the heat transport. Over the range 10^9 < R < 10^11 the enhancement of R_e^cc at constant beta due to the tilt could be described by a power law of R with an exponent of -1/6, consistent with a simple model that balances the additional buoyancy due to the tilt angle by the shear stress across the boundary layers. Even a small tilt angle dramatically suppressed the azimuthal meandering and the sudden reorientations characteristic of the LSC in a sample with beta = 0. The azimuthal mean of the temperature at the horizontal mid-plane within our resolution was independent of beta.

Asymptotic theory for a moving droplet driven by a wettability gradient

Abstract: An asymptotic theory is developed for a moving drop driven by a wettability gradient. We distinguish the mesoscale where an exact solution is known for the properly simplified problem. This solution is matched at both -- the advancing and the receding side -- to respective solutions of the problem on the microscale. On the microscale the velocity of movement is used as the small parameter of an asymptotic expansion. Matching gives the droplet shape, velocity of movement as a function of the imposed wettability gradient and droplet volume.

On the behavior of liquid drops on a solid surface. 1. The sliding of drops on an inclined surface

Abstract: The idea of contact angle was generalized by using the principle of minimum total energy. The problems of the shape of the two-dimensional sessile drop and the drop on an inclined surface are considered. The differential equations describing the drop shapes are found in both cases. A case of good wetting is considered to determine the profiles of the sessile and inclined drops. The generalization of the boundary conditions for the contact angle hysteresis of the drop on an inclined surface is discussed. The limit of the inclination angle at which the drop begins to run down the inclined plane is found. The models of the drop detachment from the solid and the motion of liquid inside of the drop are discussed. This is not an original text.

The hydraulic jump as a white hole

Abstract: In the geometry of the circular hydraulic jump, the velocity of the liquid in the interior region exceeds the speed of capillary-gravity waves (ripplons), whose spectrum is `relativistic' in the shallow water limit. The velocity flow is radial and outward, and thus the relativistic ripplons cannot propagating into the interior region. In terms of the effective 2+1 dimensional Painleve-Gullstrand metric appropriate for the propagating ripplons, the interior region imitates the white hole. The hydraulic jump represents the physical singularity at the white-hole horizon. The instability of the vacuum in the ergoregion inside the circular hydraulic jump and its observation in recent experiments on superfluid 4He by E. Rolley, C. Guthmann, M.S. Pettersen and C. Chevallier in physics/0508200 are discussed.

Fluctuation of magnetic induction in von K\'arm\'an swirling flows

Abstract: Studies of magnetic induction in von K\'arm\'an swirling flows have so far linked the time-averaged induced magnetic field to the structure of the mean flow. They have evidenced the Omega and Parker mechanism generated respectively by the flow differential rotation and helicity, which underly the Duddley and James \cite{DuddleyJames} dynamos. Using an array of Hall probes we study here the dynamical regime. In the experimental flow, turbulence is fully developed and large fluctuations are observed in the magnetic induction processes. We find that the large scale turbulent fluctuations have different characteristics when induction results from the differential rotation or from the dynamics of the stagnation point in the mid plane of the von K\'arm\'an flow. Symmetry considerations indicate that the dynamical flow spends half of its time away from the time-averaged structure. The consequences of these observations for dynamo experiments are discussed.

Scaling collapse and structure functions: Identifying self-affinity in finite length time series

Abstract: Empirical determination of the scaling properties and exponents of time series presents a formidable challenge in testing, and developing, a theoretical understanding of turbulence and other out-of-equilibrium phenomena. We discuss the special case of self affine time series in the context of a stochastic process. We highlight two complementary approaches to the differenced variable of the data: i) attempting a scaling collapse of the Probability Density Functions which should then be well described by the solution of the corresponding Fokker-Planck equation and ii) using structure functions to determine the scaling properties of the higher order moments. Using the example of a Levy flight, we demonstrate a method for conditioning the time series that recovers the underlying self affine scaling in a finite length time series.

Formation of singularities on the surface of a liquid metal in a strong electric field

Abstract: The nonlinear dynamics of the free surface of an ideal conducting liquid in a strong external electric field is studied. It is establish that the equations of motion for such a liquid can be solved in the approximation in which the surface deviates from a plane by small angles. This makes it possible to show that on an initially smooth surface for almost any initial conditions points with an infinite curvature corresponding to branch points of the root type can form in a finite time.

A model of laminated wave turbulence

Abstract: A model of laminated wave turbulence is presented. This model consists of two co-existing layers - one with continuous waves' spectra, covered by KAM theory and Kolmogorov-like power spectra, and one with discrete waves' spectra, covered by discrete classes of waves and Clipping method. Some known laboratory experiments and numerical simulations are explained in the frame of this model.

Lattice Boltzmann Method with regularized non-equilibrium distribution functions

Abstract: A new lattice Boltzmann (LB) model is introduced, based on a regularization of the pre-collision distribution functions in terms of the local density, velocity, and momentum flux tensor. The model dramatically improves the precision and numerical stability for the simulation of fluid flows by LB methods. This claim is supported by simulation results of some 2D and 3D flows.

Liquid interfaces in viscous straining flows: Numerical studies of the selective withdrawal transition

Abstract: This paper presents a numerical analysis of the transition from selective withdrawal to viscous entrainment. In our model problem, an interface between two immiscible layers of equal viscosity is deformed by an axisymmetric withdrawal flow, which is driven by a point sink located some distance above the interface in the upper layer. We find that steady-state hump solutions, corresponding to selective withdrawal of liquid from the upper layer, cease to exist above a threshold withdrawal flux, and that this transition corresponds to a saddle-node bifurcation for the hump solutions. Numerical results on the shape evolution of the steady-state interface are compared against previous experimental measurements. We find good agreement where the data overlap. However, the numerical results' larger dynamic range allows us to show that the large increase in the curvature of the hump tip near transition is not consistent with an approach towards a power-law cusp shape, an interpretation previously suggested from inspection of the experimental measurements alone. Instead the large increase in the curvature at the hump tip reflects a logarithmic coupling between the overall height of the hump and the curvature at the tip of the hump.

Turbulent-laminar patterns in plane Couette flow

Abstract: Regular patterns of turbulent and laminar fluid motion arise in plane Couette flow near the lowest Reynolds number for which turbulence can be sustained. We study these patterns using an extension of the minimal flow unit approach to simulations of channel flows pioneered by Jimenez and Moin. In our case, computational domains are of minimal size in only two directions. The third direction is taken to be large. Furthermore, the long direction can be tilted at any prescribed angle to the streamwise direction. We report on different patterned states observed as a function of Reynolds number, imposed tilt, and length of the long direction. We compare our findings to observations in large aspect-ratio experiments.

Instability of Taylor-Couette Flow between Concentric Rotating Cylinders

Abstract: The energy gradient theory is used to study the instability of Taylor-Couette flow between concentric rotating cylinders. This theory has been proposed in our previous works. In our previous studies, the energy gradient theory was demonstrated to be applicable for wall-bounded parallel flows. It was found that the critical value of the energy gradient parameter Kmax at turbulent transition is about 370-389 for wall-bounded parallel flows (which include plane Poiseuille flow, pipe Poiseuille flow and plane Couette flow) below which no turbulence occurs. In this paper, the detailed derivation for the calculation of the energy gradient parameter in the flow between concentric rotating cylinders is provided. The calculated results for the critical condition of primary instability (with semi-empirical treatment) are found to be in very good agreement with the experiments in the literature. A possible mechanism of spiral turbulence generation observed for counter-rotation of two cylinders can also be explained using the energy gradient theory. The energy gradient theory can serve to relate the condition of transition in Taylor-Couette flow to that in plane Couette flow. The latter reasonably becomes the limiting case of the former when the radii of cylinders tend to infinity. It is our contention that the energy gradient theory is possibly fairly universal for analysis of flow instability and turbulent transition, and is found valid for both pressure and shear driven flows in parallel and rotating flow configurations.

Inertial Range Scaling, Karman-Howarth Theorem and Intermittency for Forced and Decaying Lagrangian Averaged MHD in 2D

Abstract: We present an extension of the Karman-Howarth theorem to the Lagrangian averaged magnetohydrodynamic (LAMHD-alpha) equations. The scaling laws resulting as a corollary of this theorem are studied in numerical simulations, as well as the scaling of the longitudinal structure function exponents indicative of intermittency. Numerical simulations for a magnetic Prandtl number equal to unity are presented both for freely decaying and for forced two dimensional MHD turbulence, solving directly the MHD equations, and employing the LAMHD-alpha equations at 1/2 and 1/4 resolution. Linear scaling of the third-order structure function with length is observed. The LAMHD-alpha equations also capture the anomalous scaling of the longitudinal structure function exponents up to order 8.

Rough-conduit flows and the existence of fully developed turbulence

Abstract: It is widely believed that at high Reynolds number (Re) all turbulent flows approach a state of "fully developed turbulence" defined by a unique, Re-independent statistics of the velocity fluctuations. Yet direct measurements of the velocity fluctuations have failed to yield clear-cut empirical evidence of the existence of fully developed turbulence. Here we relate the friction coefficient ($f$) of rough-conduit flows to the statistics of the velocity fluctuations. In light of experimental measurements of $f$, our results yield unequivocal evidence of the existence of fully developed turbulence.

A geostrophic-like model for large Hartmann number flows

Abstract: A flow of electrically conducting fluid in the presence of a steady magnetic field has a tendency to become quasi two-dimensional, i.e. uniform in the direction of the magnetic field, except in thin so-called Hartmann boundary layers. The condition for this tendency is that of a strong magnetic field, corresponding to large values of the dimensionless Hartmann number (Ha >> 1). This is analogous to the case of low Ekman number rotating flows, with Ekman layers replacing Hartmann layers. This has been at the origin of the homogeneous model for flows in a rotating frame of reference, with its rich structure: geostrophic contours and shear layers of Stewartson, Munk and Stommel. In magnetohydrodynamics, the characteristic surfaces introduced by Kulikovskii play a role similar to the role of the geostrophic contours. However, a general theory for quasi two-dimensional magnetohydrodynamics is lacking. In this paper, a model is proposed which provides a general framework for quasi two-dimensional magnetohydrodynamic flows. Not only can this model account for otherwise disconnected past results, but it is also used to predict a new type of shear layer, of typical thickness Ha^{-1/4}.

Transport coefficients for electrolytes in arbitrarily shaped nano and micro-fluidic channels

Abstract: We consider laminar flow of incompressible electrolytes in long, straight channels driven by pressure and electro-osmosis. We use a Hilbert space eigenfunction expansion to address the general problem of an arbitrary cross section and obtain general results in linear-response theory for the hydraulic and electrical transport coefficients which satisfy Onsager relations. In the limit of non-overlapping Debye layers the transport coefficients are simply expressed in terms of parameters of the electrolyte as well as the geometrical correction factor for the Hagen-Poiseuille part of the problem. In particular, we consider the limits of thin non-overlapping as well as strongly overlapping Debye layers, respectively, and calculate the corrections to the hydraulic resistance due to electro-hydrodynamic interactions.

Magnetofluid Unification in the Yang-Mills Lagrangian

Abstract: The Yang-Mills magnetofluid unification is constructed using lagrangian approach by imposing certain gauge symmetry to the matter inside the fluid. The model provides a general description for relativistic fluid interacting with Abelian or non-Abelian gauge field. The differences with the hybrid magnetofluid model are discussed, and few physical consequences of this formalism are worked out.

Autocatalytic reaction front in a pulsative periodic flow

Abstract: Autocatalytic reaction fronts between reacted and unreacted species may propagate as solitary waves, namely at a constant front velocity and with a stationary concentration profile, resulting from a balance between molecular diffusion and chemical reaction A velocity field in the supporting medium may affect the propagation of such fronts through different phenomena: convection, diffusion enhancement, front shape changes We report here on an experimental study and lattice BGK numerical simulations of the effect of an oscillating flow on the autocatalytic reaction between iodate and arsenous acid in a Hele-Shaw cell In the low frequency range covered by the experiments, the front behavior is controlled by the flow across the gap and can be reproduced with 2D numerical simulations An analytical solution is also derived under the condition of a weak flow velocity, and is found to be in reasonable agreement with our data

Vortex Stretching and Reconnection in a Compressible Fluid

Abstract: Vortex stretching in a compressible fluid is considered. Two-dimensional and axisymmetric cases are considered separately. The flows associated with the vortices are perpendicular to the plane of the uniform straining flows. Externally-imposed density build-up near the axis leads to enhanced compactness of the vortices - "dressed" vortices (in analogy to "dressed" charged particles in a dielectric system). The compressible vortex flow solutions in the 2D as well as axisymmetric cases identify a length scale relevant for the compressible case which leads to the Kadomtsev-Petviashvili spectrum for compressible turbulence. Vortex reconnection process in a compressible fluid is shown to be possible even in the inviscid case - compressibility leads to defreezing of vortex lines in the fluid.

A Flux-Conservative Formalism for Convective and Dissipative Multi-Fluid Systems, with Application to Newtonian Superfluid Neutron Stars

Abstract: We develop a flux-conservative formalism for a Newtonian multi-fluid system, including dissipation and entrainment (i.e. allowing the momentum of one fluid to be a linear combination of the velocities of all fluids). Maximum use is made of mass, energy, and linear and angular momentum conservation to specify the equations of motion. Also used extensively are insights gleaned from a convective variational action principle, key being the distinction between each velocity and its canonically conjugate momentum. Dissipation is incorporated to second order in the ``thermodynamic forces'' via the approach pioneered by Onsager. An immediate goal of the investigation is to understand better the number, and form, of independent dissipation terms required for a consistent set of equations of motion in the multi-fluid context. A significant, but seemingly innocuous detail, is that one must be careful to isolate ``forces'' that can be written as total gradients, otherwise errors can be made in relating the net internal force to the net externally applied force. Our long-range aim is to provide a formalism that can be used to model dynamical multi-fluid systems both perturbatively and via fully nonlinear 3D numerical evolutions. To elucidate the formalism we consider the standard model for a heat-conducting, superfluid neutron star, which is believed to be dominated by superfluid neutrons, superconducting protons, and a highly degenerate, ultra-relativistic gas of normal fluid electrons. We determine that in this case there are, in principle, 19 dissipation coefficients in the final set of equations.

Gauge invariant fluid lagrangian and its application to cosmology

Abstract: A lagrangian for relativistic fluid systems with matters inside is developed using gauge principle. In the model, the gauge boson represents the fluid field in a form $A_\mu \equiv \epsilon_\mu \phi$, where $\epsilon_\mu$ contains the fluid kinematics and $\phi$ is an auxiliary field representing the fluid distribution. This leads to a new relativistic equation of motion for fluid, but which further coincides to the classical Euler equation at non-relativistic limit. The lagrangian is applied to model homogeneous universe as a bulk pure fluid system. Taking the simplest case of fluid with radial velocity and uniform distribution, the free energy density is calculated and its behaviour around Hubble distance is discussed.

Equilibrium Conditions for the Floating of Multiple Interfacial Objects

Abstract: We study the effect of interactions between objects floating at fluid interfaces, for the case in which the objects are primarily supported by surface tension. We give conditions on the density and size of these objects for equilibrium to be possible and show that two objects that float when well-separated may sink as the separation between the objects is decreased. Finally, we examine the equilbrium of a raft of strips floating at an interface, and find that rafts of sufficiently low density may have infinite spatial extent, but that above a critical raft density, all rafts sink if they are sufficiently large. We compare our numerical and asymptotic results with some simple table-top experiments, and find good quantitative agreement.

Gravity currents in a porous medium at an inclined plane

Abstract: We consider the release from a point source of relatively heavy fluid into a porous saturated medium above an impermeable slope. We consider the case where the volume of the resulting gravity current increases with time like $t^\alpha$ and show that for $\alpha 3$, this situation is reversed with spreading occurring predominantly downslope for short times. The governing equations admit similarity solutions whose scaling behaviour we determine, with the full similarity form being evaluated by numerical computations of the governing partial differential equation. We find that the results of these analyses are in good quantitative agreement with a series of laboratory experiments. Finally, we briefly discuss the implications of our work for the sequestration of carbon dioxide in aquifers with a sloping, impermeable cap.

Lagrangian dynamics of the Navier-Stokes equation

Abstract: Most researches on fluid dynamics are mostly dedicated to obtain the solutions of Navier-Stokes equation which governs fluid flow with particular boundary conditions and approximations. We propose an alternative approach to deal with fluid dynamics using the lagrangian. We attempt to develop a gauge invariant lagrangian which reconstructs the Navier-Stokes equation through the Euler-Lagrange equation. The lagrangian consists of gauge boson field $\A_\mu$ with appropriate content describing the fluid dynamics, i.e. $\A_\mu = (\Phi, -\uv)$. An example of applying the lagrangian to the interaction of fluid in a solitonic medium is also given.

Determining the Critical Condition for Flow Transition in a Full-Developed Annulus Flow

Abstract: Axial flow in an annulus between two concentric cylinders is commonly seen in various flow devices used in chemical processing industries and petroleum science and engineering. The flow state in the annulus strongly influences the performance of fluid transportation in the devices. Therefore, the determination of flow state which is laminar flow or turbulent flow is an important task to predict the performance of the flow devices. In previous works, we have proposed an energy gradient method for studying the flow instability and turbulent transition. In this method, it is shown that the flow instability and turbulent transition in wall bounded shear flows depend on the relative magnitude of the gradient of the total mechanical energy in transverse direction and the rate of loss of the total mechanical energy along the streamwise direction for a given imposed disturbance. For pipe and plane Poiseuille flows, it has been demonstrated that the transition to turbulence for these wall bounded parallel flows occurs at a consistent value of the energy gradient parameter (Kmax). In present study, the critical condition for turbulent transition in annulus flow is calculated with the energy gradient method for various radius ratios. The critical flow rate and critical Reynolds number are given for various radius ratios. Then, the analytical results are compared with the experiments in the literature. Finally, the implication of the result is discussed in terms of the drag reduction and mixing as well as heat transfer in practical industrial applications of various fluid delivery devices.