Index of content:
Volume 23, Issue 8, August 2011
The influence of solid walls on the Lagrangian statistics of statistically stationary two-dimensional turbulence is investigated by comparing the flow in a circular wall-bounded and in an unbounded periodic domain. A Fourier pseudo-spectral method is used, which is combined in the wall bounded case with a volume penalization technique to impose no-slip conditions. A particular emphasis is put on the acceleration of fluid particles. It is investigated to what extent the impact of the boundaries influences the shape of the probability density functions of Lagrangian velocity increments. It is shown that the influence of walls is not confined to a small near-wall region but alters the statistics in the entire flow domain. This can be explained by the vorticity generation in the turbulent boundary layer which destabilizes and leads to the formation of vortices that subsequently detach and travel into the bulk flow. The enstrophy level is thus increased with respect to the one in the unbounded periodic domain.
23(2011); http://dx.doi.org/10.1063/1.3619217View Description Hide Description
We observe the imbibition of concentrated suspensions of rigid spherical particles into capillary tubes. We find that while the meniscus position advances according to the square root of time, the pre-factor of this “diffusive” scaling is smaller than that predicted by an ad-hoc extension of classical imbibition using an effective viscosity for the imbibing fluid. The reduction in the pre-factor is traced to the presence of axial variations in particle volume fraction, which arise as a result of shear-induced migration of particles across streamlines. In addition, we develop a model that accounts for the effect of both radial and axial variations in particle volume fraction on the movement of the imbibition front which partially accounts for the main experimental observations.
23(2011); http://dx.doi.org/10.1063/1.3624790View Description Hide Description
Anisotropicviscous drag is usually believed to be a requirement for the low Reynolds numberlocomotion of slender bodies such as flagella and cilia. Here, we show that locomotion under isotropic drag is possible for extensible slender bodies. After general considerations, a two-ring swimmer and a model dinoflagellate flagellum are studied analytically to illustrate how extensibility can be exploited for self-propulsion without drag anisotropy. This new degree of freedom could be useful for some complex swimmer geometries and locomotion in complex fluid environments where drag anisotropy is weak or even absent.
23(2011); http://dx.doi.org/10.1063/1.3626028View Description Hide Description
Direct numerical simulations have been performed to study the effect of a stationary distribution of spanwise wall-velocity that oscillates in the streamwise direction on a turbulent boundary layer. For the first time, a spatially developing flow with this type of forcing is studied. The part of the boundary layer which flows over the alternating wall-velocity section is greatly affected with a drag reduction close to 50% which exhibits an oscillatory distribution with a wavenumber which is twice that of the imposed wall-velocity. The maximum in drag reduction occurs where the wall velocity is at its maximum (or minimum) and the minimum occurs where the wall velocity is zero. Comparisons of the mean spanwise velocity profiles with the analytical solution to the laminar Navier-Stokes equations show very good agreement. The streamwise velocity profile indicates a thickening of the viscous sub-layer when scaled with the local friction velocity and an upward shifting of the logarithmic region when scaled with the reference (unmanipulated) friction velocity. An estimation of the idealized power consumption shows that—with the present wall forcing magnitude—more energy is required for the spatial oscillation than what is saved by drag reduction.
- Biofluid Mechanics
23(2011); http://dx.doi.org/10.1063/1.3622319View Description Hide Description
Many physiological flows are driven by waves of muscular contractions passed along a tubular structure. This peristaltic pumping plays a role in ovum transport in the oviduct and in rapid sperm transport through the uterus. As such, flow due to peristalsis has been a central theme in classical biological fluid dynamics. Analytical approaches and numerical methods have been used to study flow in two-dimensional channels and three-dimensional tubes. In two dimensions, the effect of asymmetry due to a phase shift between the channel walls has been examined. However, in three dimensions, peristalsis in a non-axisymmetric tube has received little attention. Here, we present a computational model of peristaltic pumping of a viscous fluid in three dimensions based upon the method of regularized Stokeslets. In particular, we study the flow structure and mean flow in a three-dimensional tube whose asymmetry is governed by a single phase-shift parameter. We view this as a three-dimensional analog of the phase-shifted two-dimensional channel. We find that the maximum mean flow rate is achieved for the parameter that results in an axisymmetric tube. We also validate this approach by comparing our computational results with classical long-wavelength theory for the three-dimensional axisymmetric tube. This computational framework is easily implemented and may be adapted to more comprehensive physiological models where the kinematics of the tube walls are not specified a priori, but emerge due to the coupling of its passive elastic properties, force generating mechanisms, and the surrounding viscous fluid.
- Micro- and Nanofluid Mechanics
23(2011); http://dx.doi.org/10.1063/1.3611027View Description Hide Description
A solution is obtained for the ferrohydrodynamic problem of ferrofluidflow in the annular gap between two coaxial cylinders under the influence of the magnetic field generated by a multipole stator winding. The solution is applicable to small values of the applied magnetic field amplitude, takes into account the spatial harmonics produced by the multipole stator winding, and includes the effect of spin viscosity. The special case for which the internal diameter goes to zero, corresponding to the cylindrical “spin-up flow” geometry, is also obtained. Results show that for zero spin viscosity, no flow is predicted for any number of poles in the winding. In the case of nonzero spin viscosity, for a two-pole stator winding (m = 1) the solution predicts counter-rotating flow with respect to the magnetic field close to the inner cylinder and co-rotating flow with respect to the magnetic field near the outer cylinder, in agreement with previously reported experimental observations. For stator windings with larger number of poles , the solution predicts co-rotation of fluid and field in the annular gap. Using a four-pole stator winding , experimental measurements show that a ferrofluid co-rotates with the magnetic field in both the cylindrical and annular cases, with qualitative agreement between the theory with nonzero spin viscosity and the experiments. Non-idealities in the stator winding distribution, e.g., due to slot effects, are considered through linear superposition of higher order harmonics, again predicting zero flow in the absence of spin viscosity and counter-rotating flow with respect to the magnetic field close to the inner cylinder wall only for the case where the fundamental harmonic is dominant.
23(2011); http://dx.doi.org/10.1063/1.3623432View Description Hide Description
Processes of heat, momentum, and concentration transport in a boundary layer of a nanofluid near a flat wall were studied. The study was performed by means of numerical analysis of boundary layer equations in a self-similar form. Self-similar forms of these equations were obtained based on symmetry properties (Lie groups). In doing so, dependence of physical properties (viscosity, thermal conductivity, and diffusion coefficient) on concentration of nanofluids and temperature were taken into account. Effects of concentration of the nano-particles on velocity and temperature profiles, as well as on the relative Nusselt numbers and skin-friction coefficients, were elucidated.
An integral perturbation model of flow and momentum transport in rotating microchannels with smooth or microstructured wall surfaces23(2011); http://dx.doi.org/10.1063/1.3624599View Description Hide Description
This paper summarizes the development of an integral perturbation solution of the equations governing flow momentum transport and energy conversion in microchannels between disks of multiple-disk drag turbines such as Tesla turbines. Analysis of this type of flow problem is a key element in optimal design of Tesla drag-type turbines for geothermal or solar alternative energy technologies. In multiple-disk turbines, high speed flow enters tangentially at the outer radius of cylindrical microchannels formed by closely spaced parallel disks, spiraling through the channel to an exhaust at a small radius, or at the center of the disk. Previous investigations have generally developed models based on simplifying idealizations of the flow in these circumstances. Here, beginning with the momentum and continuity equations for incompressible and steady flow in cylindrical coordinates, an integral solution scheme is developed that leads to a dimensionless perturbation series solution that retains the full complement of momentum and viscous effects to consistent levels of approximation in the series solution. This more rigorous approach indicates all dimensionless parameters that affect flow and transport and allows a direct assessment of the relative importance of viscous, pressure, and momentum effects in different directions in the flow. The resulting lowest-order equations are solved explicitly and higher order terms in the series solutions are determined numerically. Enhancement of rotor drag in this type of turbine enhances energy conversion efficiency. We also developed a modified version of the integral perturbation analysis that incorporates the effects of enhanced drag due to surface microstructuring. Results of the model analysis for smooth disk walls are shown to agree well with experimental performance data for a prototype Tesla turbine and predictions of performance models developed in earlier investigations. Model predictions indicate that enhancement of disk drag by strategic microstructuring of the disk surfaces can significantly increase turbine efficiency. Exploratory calculations with the model indicate that turbine efficiencies exceeding 75% can be achieved by designing for optimal ranges of the governing dimensionless parameters.
Combined influence of streaming potential and substrate compliance on load capacity of a planar slider bearing23(2011); http://dx.doi.org/10.1063/1.3624615View Description Hide Description
In the present study, we investigate the combined interplay of streaming potential and substrate compliance with sliding dynamics on the load carrying capacity of a planar slider bearing. We relax previously utilized simplifying assumptions to model the electrokinetic effects and demonstrate that the streaming potential may augment the load carrying capacity of the bearing to a considerable extent. Interestingly, we also reveal that the electrokinetically augmented load carrying capacity exhibits strong dependencies on a combination of the compliance and the sliding dynamics, which have, hitherto, not been extensively explored. This rich interplay reveals certain parametric regimes of interest, which are significant from the viewpoint of practical design considerations.
23(2011); http://dx.doi.org/10.1063/1.3626022View Description Hide Description
In this paper, we further develop the notion of Eulerian indicators (EIs) for predicting Lagrangianmixing behavior. We employ a two-dimensional “blinking” Stokes flow as a model for mixing in a three-dimensional, spatially periodic channel flow. Each blinking flow alternates two distinct velocity fields that were calculated using a lid-driven cavity model. A new EI termed mobility is introduced to measure how effectively the blinking velocity fields transport fluid throughout the domain. We also calculate the transversality for these flows, which is an EI measuring how much the velocity direction at each point in the domain changes when the velocity fields blink. For the studied flows, we show that although individually the mobility and transversality do not correlate well with mixing as measured by the decay of the variance of concentration, the product of mobility and transversality does correlate well with the decay of the variance of concentration and predicts which combinations of velocity fields will produce the most effective mixing.
- Interfacial Flows
23(2011); http://dx.doi.org/10.1063/1.3615643View Description Hide Description
Using a lattice Boltzmann multiphase model, three-dimensional numerical simulations have been performed to understand dropletformation in microfluidic cross-junctions at low capillary numbers. Flow regimes, consequence of interaction between two immiscible fluids, are found to be dependent on the capillary number and flow rates of the continuous and dispersed phases. A regime map is created to describe the transition from dropletsformation at a cross-junction (DCJ), downstream of cross-junction to stable parallel flows. The influence of flow rate ratio, capillary number, and channel geometry is then systematically studied in the squeezing-pressure-dominated DCJ regime. The plug length is found to exhibit a linear dependence on the flow rate ratio and obey power-law behavior on the capillary number. The channel geometry plays an important role in droplet breakup process. A scaling model is proposed to predict the plug length in the DCJ regime with the fitting constants depending on the geometrical parameters.
23(2011); http://dx.doi.org/10.1063/1.3623425View Description Hide Description
The damping rate of the small free oscillations in a non-cylindrical, axi-symmetric liquid bridge between two circular disks is calculated and compared with some previous experimental measurements using hexadecane in a millimetric liquid bridge. Current theories, accounting for viscous damping in both the boundary layers attached to the disks and the bulk, underestimated the measured damping by a O(1) quantity; and no improvement resulted from calculations based on the full Navier-Stokes equations. These discrepancies are essentially eliminated in this paper considering the effect of the surface shear viscosity (whose value results from empirical fitness), which could be due to the presence of a contaminating monolayer. Some consequences are extracted in connection with surface wave damping in micro-fluidic devices.
23(2011); http://dx.doi.org/10.1063/1.3627150View Description Hide Description
The study of convective thermocapillary instabilities in liquid bridges [J. J. Xu and S. H. Davis, Phys. Fluids 27(5), 1102 (1984)] is revisited. A new branch of neutral mode m = 1 is found. The previously reported results are confirmed in the range of low Prandtl numbers. It is shown that for large Prandtl numbers, the flow becomes unstable at much smaller values of the Marangoni number than it was reported previously. The calculations are performed for adiabatic and heat conductivefree surface. In both cases, the critical mode is m = 1. The previously reported change of critical mode from m = 1 to m = 0 with increasing the Prandtl number is not confirmed. The corrected results provide a better agreement with the experimental data.
- Viscous and Non-Newtonian Flows
23(2011); http://dx.doi.org/10.1063/1.3615518View Description Hide Description
We study the steady, three-dimensional creeping, and viscoelasticflow around a freely rotating rigid sphere subject to simple shear flow imposed at infinity. The viscoelasticity of the ambient fluid is modeled using the second-order-fluid model, the Upper Convected Maxwell, the exponential affine Phan-Thien-Tanner, and the Giesekus constitutive equations. A spherical coordinate system with origin at the center of the sphere is used to describe the flow field. The solution of the governing equations is expanded as a series for small values of the Deborah number. The resulting sequence of differential equations is solved analytically up to second order and numerically up to fourth order in Deborah number by employing fully spectral representations for all the primary variables. In particular, Chebyshev polynomials are used in the radial coordinate and the double Fourier series in the longitudinal and latitudinal coordinates. The numerical results up to second-order agree within machine accuracy with the available analytical solutions clearly indicating the correctness and accuracy of the numerical method developed here. Analytical expressions for the angular velocity of the rigid sphere up to fourth order, which show the slowdown of the rotation of the sphere with respect to the Newtonian creeping case, are also derived. For small Deborah numbers, these expressions, along with those presented in a recent letter [Housiadas and Tanner, Phys. Fluids 23, 051702 (2011)] are in agreement with the few available experimental data and numerical results.
23(2011); http://dx.doi.org/10.1063/1.3624697View Description Hide Description
In this paper, the steady axisymmetric stagnation point flow of a viscous and incompressible fluid over a shrinking circular cylinder with mass transfer (suction) is studied. The flow is induced by a cylinder shrinking with a linear velocity distribution from the stagnation line. The fluid flowsolution is an exact solution of the Navier-Stokes equation, which is reduced to a nonlinear ordinary differential equation. This equation is solved numerically for some values of the governing parameters that involves a Reynolds numberR, a shrinking parameter λ, and a suction parameter γ. The effects of these governing parameters on the velocity and temperature profiles, skin friction coefficient, Nusselt number, as well as the distributions of the streamlines are investigated. The obtained results for the case of fixed cylinder are compared with those from the open literature and it is shown that they are in excellent agreement. The solutions are non-unique for some values of λ . The streamlines show that the flow structures become complicated due to the shrinking and suction effects.
- Particulate, Multiphase, and Granular Flows
23(2011); http://dx.doi.org/10.1063/1.3623275View Description Hide Description
Although the segregation of dry granular surfaceflow has been widely studied, cases where the air is completely replaced by a liquid still remain unexplored. In this study, we report on experiments performed to investigate the phenomena of particle segregation and flowing behavior in a rotating drum with liquids of different viscosities and different filling degrees. The experimental results indicate that the viscosity of the interstitial fluid has a significant effect on the granular flow in the slurrygranular flow. The segregation index and angle of repose are shown to decrease with increased liquid viscosity. When the liquid viscosity is the same, the increase in filling degree causes the segregation index to increase, while the net rate of mixing seems to decrease. A new dimensionless flow variable is used to distinguish the flow regimes. We find that the flow regime changes from rolling regime to cascading regime when the dimensionless flow variable is below a critical value. Furthermore, the change of segregation index occurs during the transition of the granular flow regime.
Roles of particle-wall and particle-particle interactions in highly confined suspensions of spherical particles being sheared at low Reynolds numbers23(2011); http://dx.doi.org/10.1063/1.3613972View Description Hide Description
The roles of particle-wall and particle-particle interactions are examined for suspensions of spherical particles in a viscous fluid being confined and sheared at low Reynolds numbers by two parallel walls moving with equal but opposite velocities. Both particle-wall and particle-particle interactions are shown to decrease the rotational velocity of the spheres, so that in the limit of vanishingly small gaps between the spheres and the walls, the spheres acquire a rotational slip relative to the walls. The presence of the walls also increases the particle stresslet and, therefore, the total viscous dissipation. In the limit of vanishingly small gaps, the increased viscous dissipation in the gaps between pairs of spheres aligned in the flow direction is largely compensated by the reduction in the dissipation in the gaps between the spheres and the walls due to a reduction in the rotational velocity of the spheres. As a result, the effect of short-range particle interactions on the stresslet is generally insignificant. On the other hand, the channel-scale particle interactions in the shear flow induced by the moving walls decrease the particle stresslet, primarily because the fraction of pairs of spheres that are aligned parallel to the flow (the presence of which in a shear flow reduces the stresslet) is relatively higher than in unbounded suspensions. Expressions are also derived for the total stress in dilute random suspensions that account for both the particle-wall and the channel-scale particle-particle interactions in determining the rotational velocities and stresses. The latter are shown to be consistent with recent numerical [Y. Davit and P. Peyla, Europhys. Lett. 83, 64001 (2008)] and experimental [P. Peyla and C. Verdier, Europhys. Lett. 94, 44001 (2011)] findings according to which, for a range of sphere radius to gap width ratios, the effect of particle-particle interactions is to decrease the total dissipation.
First-order virial expansion of short-time diffusion and sedimentation coefficients of permeable particles suspensions23(2011); http://dx.doi.org/10.1063/1.3626196View Description Hide Description
For suspensions of permeable particles, the short-time translational and rotational self-diffusion coefficients, and collective diffusion and sedimentation coefficients are evaluated theoretically. An individual particle is modeled as a uniformly permeable sphere of a given permeability, with the internal solvent flow described by the Debye-Bueche-Brinkman equation. The particles are assumed to interact non-hydrodynamically by their excluded volumes. The virial expansion of the transport properties in powers of the volume fraction is performed up to the two-particle level. The first-order virial coefficients corresponding to two-body hydrodynamic interactions are evaluated with very high accuracy by the series expansion in inverse powers of the inter-particle distance. Results are obtained and discussed for a wide range of the ratio, x, of the particle radius to the hydrodynamic screening length inside a permeable sphere. It is shown that for , the virial coefficients of the transport properties are well-approximated by the hydrodynamic radius (annulus) model developed by us earlier for the effective viscosity of porous-particle suspensions.
- Laminar Flows
23(2011); http://dx.doi.org/10.1063/1.3614482View Description Hide Description
Flow separation control by a non-thermal plasma actuator is considered for a NACA 0015 airfoil at a chord Reynolds number of 1.9 × 105. Static hysteresis in the lift coefficient is demonstrated for increasing and then decreasing sinusoidal voltage amplitude supplying a typical single dielectric barrier discharge actuator at the leading edge of the model. In addition to these open-loop experiments, unsteady surface pressure signals are examined for transient processes involving forced reattachment and natural separation. The results show that strong pressure oscillations in the relatively slow separation process, compared to reattachment, precede the ultimate massive flow separation. To enhance the contrast between the parts of the signal related to the attached flow and those related to the incipient separation, RMS estimate of filtered values of is used to define a flow separation predictor that is implemented in feedback control. Two simple controllers are proposed, one based on a predefined threshold of the unsteady Cp and another that utilizes the flow separation predictor to identify incipient separation. The latter effectively leverages the hysteresis in the post-stall regime to reduce the electrical power consumed by the actuator while maintaining continuously attached flow.
23(2011); http://dx.doi.org/10.1063/1.3623422View Description Hide Description
Direct numerical simulations are performed to examine roles of streamwise dynamics in the spreading and mixing of flows in a two-step rectangular sudden expansion channel. The configuration is observed to facilitate higher entrainment by virtue of developed inflow type streamwise dynamics, and the system-generated passive forcing provided the necessary impetus for sustainable growth/evolution of the vortices. In addition, with the introduction and proper placement of two tiny rectangular “tabs” over the inlet walls, the downstream growth/dynamics of the vortices could be effectively modified. Through their inflow/outflow type dynamics, the streamwise vortices are found to decisively dictate the transverse jet spreading. The physical process thereby either led to a quick axis switching of the jet, or stopped axis switching altogether. However, the stretching of the azimuthal vortices is found to remain directly linked to the streamwise continuation of jet’s azimuthal curvature variation. A pressure analysis as presented herein, and the simulated nature of dynamics of two azimuthal components of vortices further reveal that the transverse pressure gradient skewing is a dominant source of streamwise vorticity in such flows. With the knowledge of simulated transverse pressure distribution, we formulate here a mechanism which efficiently predicts the inception/dynamics of all the streamwise vortices. Moreover, a novel mathematical foundation in support of the pressure analysis has been outlined here, which ensures further broader universal scope of applicability of our proposed pressure analysis.
- Instability and Transition
Sidewall and thermal boundary condition effects on the evolution of longitudinal rolls in Rayleigh-Bénard-Poiseuille convection23(2011); http://dx.doi.org/10.1063/1.3605698View Description Hide Description
Experimental and numerical studies of steady longitudinal convection rolls that develop in a Poiseuille air flow in a rectangular channel heated from below and cooled from the top are conducted in the range 3500 ≤ Ra ≤ 6000 and 20 ≤ Re ≤ 200. The effect of the lateral vertical walls on the onset and development of the convection cells is investigated by changing the transverse aspect ratio of the channel from 4.7 to 18.4. The influence of the entrance temperature and adiabatic or conductive thermal boundary conditions at the side and top walls of the channel is also investigated. The scenario of the roll formation is described in details. It results in a symmetric pattern in the form of steady longitudinal rolls with an even number of rolls that depends not only on the aspect ratio but possibly on the inlet temperature of the flow. It is shown that the fully developed pattern is determined by the two rolls nearby each vertical side wall that are triggered just at the entrance of the channel due to the presence of velocity boundary layers adjacent to the walls. It is also shown that the heat conduction in the top horizontal wall of the experimental channel must be taken into account in the numerical simulations so that the experimental wavenumber can be properly depicted.