Showing papers in "arXiv: Fluid Dynamics in 2008"

Aeroelastic instability of cantilevered flexible plates in uniform flow

Abstract: We address the flutter instability of a flexible plate immersed in an axial flow This instability is similar to flag flutter and results from the competition between destabilising pressure forces and stabilising bending stiffness In previous experimental studies, the plates have always appeared much more stable than the predictions of two-dimensional models This discrepancy is discussed and clarified in this paper by examining experimentally and theoretically the effect of the plate aspect ratio on the instability threshold We show that the two-dimensional limit cannot be achieved experimentally because hysteretical behaviour and three-dimensional effects appear for plates of large aspect ratio The nature of the instability bifurcation (sub- or supercritical) is also discussed in the light of recent numerical results

Equilibrium and traveling-wave solutions of plane Couette flow

Abstract: We present ten new equilibrium solutions to plane Couette flow in small periodic cells at low Reynolds number (Re) and two new traveling-wave solutions. The solutions are continued under changes of Re and spanwise period. We provide a partial classification of the isotropy groups of plane Couette flow and show which kinds of solutions are allowed by each isotropy group. We find two complementary visualizations particularly revealing. Suitably chosen sections of their 3D-physical space velocity fields are helpful in developing physical intuition about coherent structures observed in low Re turbulence. Projections of these solutions and their unstable manifolds from their infinite-dimensional state space onto suitably chosen 2- or 3-dimensional subspaces reveal their interrelations and the role they play in organizing turbulence in wall-bounded shear flows.

Tensorial hydrodynamic slip

Abstract: We describe a tensorial generalization of the Navier slip boundary condition and illustrate its use in solving for flows around anisotropic textured surfaces. Tensorial slip can be derived from molecular or microstructural theories or simply postulated as an constitutive relation, subject to certain general constraints on the interfacial mobility. The power of the tensor formalism is to capture complicated effects of surface anisotropy, while preserving a simple fluid domain. This is demonstrated by exact solutions for laminar shear flow and pressure-driven flow between parallel plates of arbitrary and different textures. From such solutions, the effects of rotating a texture follow from simple matrix algebra. Our results may be useful to extracting local slip tensors from global measurements, such as the permeability of a textured channel or the force required to move a patterned surface, in experiments or simulations.

A model for the symmetry breaking of the reverse Benard-von Karman vortex street produced by a flapping foil

Abstract: The vortex streets produced by a flapping foil of span-to-chord aspect ratio of 4:1 are studied in a hydrodynamic tunnel experiment. In particular, the mechanisms giving rise to the symmetry breaking of the reverse B\'enard-von K\'arm\'an vortex street that characterizes fish-like swimming and forward flapping flight are examined. Two-dimensional particle image velocimetry measurements in the mid-plane perpendicular to the span axis of the foil are used to characterize the different flow regimes. The deflection angle of the mean jet flow with respect to the horizontal observed in the average velocity field is used as a measure of the asymmetry of the vortex street. Time series of the vorticity field are used to calculate the advection velocity of the vortices with respect to the free-stream, defined as the phase velocity $U_{phase}$, as well as the circulation $\Gamma$ of each vortex and the spacing $\xi$ between consecutive vortices in the near wake. The observation that the symmetry breaking results from the formation of a dipolar structure from each couple of counter-rotating vortices shed on each flapping period serves as starting point to build a model for the symmetry breaking threshold. A symmetry breaking criterion based on the relation between the phase velocity of the vortex street and an idealized self-advection velocity of two consecutive counter-rotating vortices in the near wake is established. The predicted threshold for symmetry breaking accounts well for the deflected wake regimes observed in the present experiments and may be useful to explain other experimental and numerical observations of similar deflected propulsive vortex streets reported in the literature.

Reflexion and Diffraction of Internal Waves analyzed with the Hilbert Transform

Abstract: We apply the Hilbert transform to the physics of internal waves in two-dimensional fluids. Using this demodulation technique, we can discriminate internal waves propagating in different directions: this is very helpful in answering several fundamental questions in the context of internal waves. We focus more precisely in this paper on phenomena associated with dissipation, diffraction and reflection of internal waves.

Wall-oscillation conditions for drag reduction in turbulent channel flow

Abstract: The drag reduction properties of a turbulent channel flow modified by spanwise sinusoidal oscillations of the walls are investigated by direct numerical simulations. The work is based on the linear relation between the drag reduction and the parameter $S$, function of the maximum wall velocity and the period of the oscillation. This quantity has been found through physical arguments pertaining to the action of the oscillating Stokes layer on the near-wall turbulence dynamics. The predictive potential of the scaling parameter is exploited to gain insight into the drag-reducing effects of the oscillating wall technique. The period of oscillation which guarantees the maximum drag reduction for a given maximum wall displacement is studied for the first time. The issue of the minimum intensity of wall forcing required to produce a non-zero drag reduction effect and the dependence of the drag reduction on the Reynolds number are also addressed. The drag reduction data available in the literature are compared with the prediction given by the scaling parameter, thus attaining a comprehensive view of the state of the art.

Oscillations of the large-scale circulation in turbulent Rayleigh-B\'{e}nard convection: the off-center mode and its relationship with the torsional mode

Abstract: We report an experimental study of the large-scale circulation (LSC) in a turbulent Rayleigh-B\'{e}nard convection cell with aspect ratio unity. The temperature-extremum-extraction (TEE) method for obtaining the dynamic information of the LSC is presented. With this method, the azimuthal angular positions of the hot ascending and cold descending flows along the sidewall are identified from the measured instantaneous azimuthal temperature profile. The motion of the LSC is then decomposed into two different modes: the azimuthal mode and the translational or off-center mode. Comparing to the previous sinusoidal-fitting (SF) method, it is found that both methods give the same information about the azimuthal motion of the LSC, but the TEE method in addition can provide information about the off-center motion of the LSC, which is found to oscillate time-periodically around the cell's central vertical axis with an amplitude being nearly independent of the turbulent intensity. It is further found that the azimuthal angular positions of the hot ascending flow near the bottom plate and the cold descending flow near the top plate oscillate periodically out of phase by $\pi$, leading to the torsional mode of the LSC. These oscillations are then propagated vertically along the sidewall by the hot ascending and cold descending fluids. When they reach the mid-height plane, the azimuthal positions of the hottest and coldest fluids again oscillate out of phase by $\pi$. It is this out-of-phase horizontal positional oscillation of the hottest and coldest fluids at the same horizontal plane that produces the off-center oscillation of the LSC. A direct velocity measurement further confirms the existence of the bulk off-center mode of the flow field near cell center.

Supercritical transition to turbulence in an inertially-driven von Karman closed flow

Abstract: We study the transition from laminar flow to fully developed turbulence for an inertially-driven von Karman flow between two counter-rotating large impellers fitted with curved blades over a wide range of Reynolds number (100 - 1 000 000). The transition is driven by the destabilisation of the azimuthal shear-layer, i.e., Kelvin-Helmholtz instability which exhibits travelling/drifting waves, modulated travelling waves and chaos below the emergence of a turbulent spectrum. A local quantity -the energy of the velocity fluctuations at a given point- and a global quantity -the applied torque- are used to monitor the dynamics. The local quantity defines a critical Reynolds number Rec for the onset of time-dependence in the flow, and an upper threshold/crossover Ret for the saturation of the energy cascade. The dimensionless drag coefficient, i.e., the turbulent dissipation, reaches a plateau above this finite Ret, as expected for a "Kolmogorov"-like turbulence for Re -> infinity. Our observations suggest that the transition to turbulence in this closed flow is globally supercritical: the energy of the velocity fluctuations can be considered as an order parameter characterizing the dynamics from the first laminar time-dependence up to the fully developed turbulence. Spectral analysis in temporal domain moreover reveals that almost all of the fluctuations energy is stored in time-scales one or two orders of magnitude slower than the time-scale based on impeller frequency.

Heteroclinic connections in plane Couette flow

Abstract: Plane Couette flow transitions to turbulence for Re~325 even though the laminar solution with a linear profile is linearly stable for all Re (Reynolds number). One starting point for understanding this subcritical transition is the existence of invariant sets in the state space of the Navier Stokes equation, such as upper and lower branch equilibria and periodic and relative periodic solutions, that are quite distinct from the laminar solution. This article reports several heteroclinic connections between such objects and briefly describes a numerical method for locating heteroclinic connections. Computing such connections is essential for understanding the global dynamics of spatially localized structures that occur in transitional plane Couette flow. We show that the nature of streaks and streamwise rolls can change significantly along a heteroclinic connection.

Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

Abstract: On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.

Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water

Abstract: Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta ("non-Oberbeck-Boussinesqness") for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.

Fine-scale statistics of temperature and its derivatives in convective turbulence

Abstract: We study the fine-scale statistics of temperature and its derivatives in turbulent Rayleigh-Benard convection. Direct numerical simulations are carried out in a cylindrical cell with unit aspect ratio filled with a fluid with Prandtl number equal to 0.7 for Rayleigh numbers between 10^7 and 10^9. The probability density function of the temperature or its fluctuations is found to be always non-Gaussian. The asymmetry and strength of deviations from the Gaussian distribution are quantified as a function of the cell height. The deviations of the temperature fluctuations from the local isotropy, as measured by the skewness of the vertical derivative of the temperature fluctuations, decrease in the bulk, but increase in the thermal boundary layer for growing Rayleigh number, respectively. Similar to the passive scalar mixing, the probability density function of the thermal dissipation rate deviates significantly from a log-normal distribution. The distribution is fitted well by a stretched exponential form. The tails become more extended with increasing Rayleigh number which displays an increasing degree of small-scale intermittency of the thermal dissipation field for both the bulk and the thermal boundary layer. We find that the thermal dissipation rate due to the temperature fluctuations is not only dominant in the bulk of the convection cell, but also yields a significant contribution to the total thermal dissipation in the thermal boundary layer. This is in contrast to the ansatz used in scaling theories and can explain the differences in the scaling of the total thermal dissipation rate with respect to the Rayleigh number.

Relative periodic orbits in transitional pipe flow

Abstract: A dynamical system description of the transition process in shear flows with no linear instability starts with a knowledge of exact coherent solutions, among them travelling waves (TWs) and relative periodic orbits (RPOs). We describe a numerical method to find such solutions in pipe flow and apply it in the vicinity of a Hopf bifurcation from a TW which looks to be especially relevant for transition. The dominant structural feature of the RPO solution is the presence of weakly modulated streaks. This RPO, like the TW from which it bifurcates, sits on the laminar-turbulent boundary separating initial conditions which lead to turbulence from those which immediately relaminarise.

Galilean invariance of lattice Boltzmann models

Abstract: It is well-known that the original lattice Boltzmann (LB) equation deviates from the Navier-Stokes equations due to an unphysical velocity dependent viscosity. This unphysical dependency violates the Galilean invariance and limits the validation domain of the LB method to near incompressible flows. As previously shown, recovery of correct transport phenomena in kinetic equations depends on the higher hydrodynamic moments. In this Letter, we give specific criteria for recovery of various transport coefficients. The Galilean invariance of a general class of LB models is demonstrated via numerical experiments.

Direct Assessment of Vorticity Alignment with Local and Nonlocal Strain Rates in Turbulent Flows

Abstract: A direct Biot-Savart integration is used to decompose the strain rate into its local and nonlocal constituents, allowing the vorticity alignment with the local and nonlocal strain rate eigenvectors to be investigated. These strain rate tensor constituents are evaluated in a turbulent flow using data from highly-resolved direct numerical simulations. While the vorticity aligns preferentially with the intermediate eigenvector of the \textit{combined} strain rate, as has been observed previously, the present results for the first time clearly show that the vorticity aligns with the most extensional eigenvector of the \textit{nonlocal} strain rate. This in turn reveals a significant linear contribution to the vortex stretching dynamics in turbulent flows.

Turbulent Flow in Pipes and Channels as Cross-Stream "Inverse Cascades" of Vorticity

Abstract: A commonplace view of pressure-driven turbulence in pipes and channels is as "cascades" of streamwise momentum toward the viscous layer at the wall. We present in this paper an alternative picture of these flows as "inverse cascades" of spanwise vorticity, in the cross-stream direction but away from the viscous sublayer. We show that there is a constant spatial flux of spanwise vorticity, due to vorticity conservation, and that this flux is necessary to produce pressure-drop and energy dissipation. The vorticity transport is shown to be dominated by viscous diffusion at distances closer to the wall than the peak Reynolds stress, well into the classical log-layer. The Perry-Chong model based on "representative" hairpin/horsehoe vortices predicts a single sign of the turbulent vorticity flux over the whole log-layer, whereas the actual flux must change sign at the location of the Reynolds-stress maximum. Sign-reversal may be achieved by assuming a slow power-law decay of the Townsend "eddy intensity function" for wall-normal distances greater than the hairpin length-scale. The vortex-cascade picture presented here has a close analogue in the theory of quantum superfluids and superconductors, the "phase slippage" of quantized vortex lines. Most of our results should therefore apply as well to superfluid turbulence in pipes and channels. We also discuss issues about drag-reduction from this perspective.

A public turbulence database cluster and applications to study Lagrangian evolution of velocity increments in turbulence

Abstract: A public database system archiving a direct numerical simulation (DNS) data set of isotropic, forced turbulence is described in this paper. The data set consists of the DNS output on $1024^3$ spatial points and 1024 time-samples spanning about one large-scale turn-over timescale. This complete $1024^4$ space-time history of turbulence is accessible to users remotely through an interface that is based on the Web-services model. Users may write and execute analysis programs on their host computers, while the programs make subroutine-like calls that request desired parts of the data over the network. The users are thus able to perform numerical experiments by accessing the 27 Terabytes of DNS data using regular platforms such as laptops. The architecture of the database is explained, as are some of the locally defined functions, such as differentiation and interpolation. Test calculations are performed to illustrate the usage of the system and to verify the accuracy of the methods. The database is then used to analyze a dynamical model for small-scale intermittency in turbulence. Specifically, the dynamical effects of pressure and viscous terms on the Lagrangian evolution of velocity increments are evaluated using conditional averages calculated from the DNS data in the database. It is shown that these effects differ considerably among themselves and thus require different modeling strategies in Lagrangian models of velocity increments and intermittency.

The critical layer in pipe flow at high Reynolds number

Abstract: We report the computation of a family of traveling wave solutions of pipe flow up to $Re=75000$. As in all lower-branch solutions, streaks and rolls feature prominently in these solutions. For large $Re$, these solutions develop a critical layer away from the wall. Although the solutions are linearly unstable, the two unstable eigenvalues approach 0 as $Re\to\infty$ at rates given by $Re^{-0.41}$ and $Re^{-0.87}$ -- surprisingly, the solutions become more stable as the flow becomes less viscous. The formation of the critical layer and other aspects of the $Re\to\infty$ limit could be universal to lower-branch solutions of shear flows. We give implementation details of the GMRES-hookstep and Arnoldi iterations used for computing these solutions and their spectra, while pointing out the new aspects of our method.

MHD turbulence in a channel with spanwise magnetic field

Abstract: The effect of a uniform spanwise magnetic field on a turbulent channel flow is investigated for the case of low magnetic Reynolds number. DNS and LES computations are performed for two values of the hydrodynamic Reynolds number (10^4 and 2\times 10^4) and the Hartmann number varying in a wide range. It is shown that the main effect of the magnetic field is the suppression of turbulent velocity fluctuations and momentum transfer in the wall-normal direction. This leads to drag reduction and transformation of the mean flow profile. The centerline velocity grows, the mean velocity gradients near the wall decrease, and the typical horizontal dimensions of the coherent structures enlarge upon increasing the Hartmann number. Comparison between LES and DNS results shows that the dynamic Smagorinsky model accurately reproduces the flow transformation.

Phase-field model for the Rayleigh--Taylor instability of immiscible fluids

Abstract: The Rayleigh--Taylor instability of two immiscible fluids in the limit of small Atwood numbers is studied by means of a phase-field description. In this method the sharp fluid interface is replaced by a thin, yet finite, transition layer where the interfacial forces vary smoothly. This is achieved by introducing an order parameter (the phase field) whose variation is continuous across the interfacial layers and is uniform in the bulk region. The phase field model obeys a Cahn--Hilliard equation and is two-way coupled to the standard Navier--Stokes equations. Starting from this system of equations we have first performed a linear analysis from which we have analytically rederived the known gravity-capillary dispersion relation in the limit of vanishing mixing energy density and capillary width. We have performed numerical simulations and identified a region of parameters in which the known properties of the linear phase (both stable and unstable) are reproduced in a very accurate way. This has been done both in the case of negligible viscosity and in the case of nonzero viscosity. In the latter situation only upper and lower bounds for the perturbation growth-rate are known. Finally, we have also investigated the weakly-nonlinear stage of the perturbation evolution and identified a regime characterized by a constant terminal velocity of bubbles/spikes. The measured value of the terminal velocity is in perfect agreement with available theoretical prediction. The phase-field approach thus appears to be a valuable tecnhique for the dynamical description of the stages where hydrodynamic turbulence and wave-turbulence enter into play.

Large scale flow around turbulent spots

Abstract: Numerical simulations of a model of plane Couette flow focusing on its in-plane spatio-temporal properties are used to study the dynamics of turbulent spots. While the core of a spot is filled with small scale velocity fluctuations, a large scale flow extending far away and occu- pying the full gap between the driving plates is revealed upon filtering out small scales. It is characterized by streamwise inflow towards the spot and spanwise outflow from the spot, giving it a quadrupolar shape. A correction to the base flow is present within the spot in the form of a spanwise vortex with vorticity opposite in sign to that of the base flow. The Reynolds stresses are shown to be at the origin of this recirculation, whereas the quadrupolar shape of the in-plane flow results from the transport of this recirculation by the base flow that pumps it towards the spot in the streamwise direction and flushes it in the spanwise direction to insure mass conservation. These results shed light on earlier observations in plane Couette flow or other wall flows experiencing a direct transition to turbulence by spot nucleation.

Self-Similarity and Universality in Rayleigh-Taylor, Boussinesq Turbulence

Abstract: We report and discuss case study simulations of the Rayleigh-Taylor instability in the Boussinesq, incompressible regime developed to turbulence. Our main focus is on a statistical analysis of density and velocity fluctuations inside of the already developed and growing in size mixing zone. Novel observations reported in the manuscript concern self-similarity of the velocity and density fluctuations spectra inside of the mixing zone snapshot, independence of the spectra of the horizontal slice level, and universality showing itself in a virtual independence of the internal structure of the mixing zone, measured in the re-scaled spatial units, of the initial interface perturbations.

Analytical solution of Stokes flow inside an evaporating sessile drop: Spherical and cylindrical cap shapes

Abstract: Exact analytical solutions are derived for the Stokes flows within evaporating sessile drops of spherical and cylindrical cap shapes. The results are valid for arbitrary contact angle. Solutions are obtained for arbitrary evaporative flux distributions along the free surface as long as the flux is bounded at the contact line. The field equations, E^4(Psi)=0 and Del^4(Phi)=0, are solved for the spherical and cylindrical cap cases, respectively. Specific results and computations are presented for evaporation corresponding to uniform flux and to purely diffusive gas phase transport into an infinite ambient. Wetting and non-wetting contact angles are considered with the flow patterns in each case being illustrated. For the spherical cap with evaporation controlled by vapor phase diffusion, when the contact angle lies in the range 0

Direct Numerical Simulation of Turbulent Heat Transfer Modulation in Micro-Dispersed Channel Flow

Abstract: The object of this paper is to study the influence of dispersed micrometer size particles on turbulent heat transfer mechanisms in wall-bounded flows. The strategic target of the current research is to set up a methodology to size and design new-concept heat transfer fluids with properties given by those of the base fluid modulated by the presence of dynamically-interacting, suitably-chosen, discrete micro- and nano- particles. We run Direct Numerical Simulation (DNS) for hydrodynamically fully-developed, thermally-developing turbulent channel flow at shear Reynolds number Re=150 and Prandtl number Pr=3, and we tracked two large swarms of particles, characterized by different inertia and thermal inertia. Preliminary results on velocity and temperature statistics for both phases show that, with respect to single-phase flow, heat transfer fluxes at the walls increase by roughly 2% when the flow is laden with the smaller particles, which exhibit a rather persistent stability against non-homogeneous distribution and near-wall concentration. An opposite trend (slight heat transfer flux decrease) is observed when the larger particles are dispersed into the flow. These results are consistent with previous experimental findings and are discussed in the frame of the current research activities in the field. Future developments are also outlined.

Osmotically driven pipe flows and their relation to sugar transport in plants

Abstract: In plants, osmotically driven flows are believed to be responsible for translocation of sugar in the pipe-like phloem cell network, spanning the entire length of the plant. In this paper, we present an experimental and theoretical study of transient osmotically driven flows through pipes with semipermeable walls. We extend the experimental work of Eschrich, Evert and Young \cite[]{Eschrich:1972} by providing a more accurate version of their experiment allowing for better comparison with theory. In the experiments we measure the dynamics and structure of a "sugar front", i.e. the transport and decay of a sudden loading of sugar in a pipe which is closed in both ends. We include measurements of pressure inside the membrane tube allowing us to compare the experiments directly with theory and, in particular, to confirm quantitatively the exponential decay of the front in a closed this http URL a novel setup we are able to measure the entire concentration profile as the sugar front moves. In contrast to previous studies we find very good agreement between experiment and theory. In the limit of low axial resistance (valid in our experiments as well as in many cases in plants) we show that the equations can be solved exactly by the method of characteristics yielding, in general, an implicit solution. Further we show that under more general conditions the equations of motion can be rewritten as a single integro-differential equation, which can be readily solved numerically. The applicability of our results to plants is discussed and it is shown that it is probable that the pressure-flow hypothesis can account for short distance transport of sugar in plants.

Evidence of thin-film precursors formation in hydrokinetic and atomistic simulations of nano-channel capillary filling

Abstract: We present hydrokinetic Lattice Boltzmann and Molecular Dynamics simulations of capillary filling of high-wetting fluids in nano-channels, which provide clear evidence of the formation of thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor films is found to obey the Lucas-Washburn law as the main capillary front, z2(t) proportional to t, although with a larger prefactor, which we find to take the same value for both geometries under inspection. Both hydrokinetic and Molecular Dynamics approaches indicate a precursor film thickness of the order of one tenth of the capillary diameter. The quantitative agreement between the hydrokinetic and atomistic methods indicates that the formation and propagation of thin precursors can be handled at a mesoscopic/hydrokinetic level, thereby opening the possibility of using hydrokinetic methods to space-time scales and complex geometries of direct experimental relevance.

Rotational dynamics of a soft filament: wrapping transition and propulsive forces

Abstract: We analyze experimentally the shape of a long elastic filament rotating in a viscous liquid. We identify a continuous but sharp transition from a straight to an helical shape, resulting from the competition between viscous stresses and elastic forces. This induced helicity generates a propulsive force along the axis of rotation. In addition, we show that the shape transition is associated with an unstable branch in the force-torque relation. A linearized model of the fluid-structure interaction is proposed to account for all the features of the non-linear filament dynamics.

Electro-viscous effects in capillary filling of nanochannels

Abstract: We theoretically examine the widespread hypothesis of an electro-viscous origin of the increase in apparent viscosity observed in recent experiments on capillary filling of nanochannels. Including Debye-layer corrections to the hydraulic resistance we find that the apparent viscosity reaches a maximum in the mesoscopic regime where the channel height (or more generally the hydraulic radius) is comparable to the screening length. However, for realistic estimates of central parameters, we find that the electro-viscous contribution to the apparent viscosity is at most a 1% effect.

Wetting failure and contact line dynamics in a Couette flow

Abstract: Liquid-liquid wetting failure is investigated in a two-dimensional Couette system with two immiscible fluids of arbitrary viscosity. The problem is solved exactly using a sharp interface treatment of hydrodynamics (lubrication theory) as a function of the capillary number, viscous ratio and separation of scale, i.e. slip length versus macroscopic scale of the system. The existence of critical velocities, above which no stationary solutions are found, is analyzed in detail in terms of the relevant parameters of the system. Comparisons with existing analysis for other geometries are also carried out. A numerical method of analysis is also presented, based on diffuse interface models obtained from multiphase extensions of the lattice Boltzmann equation (LBE). Sharp interface and diffuse interface models are quantitatively compared face to face indicating the correct limit of applicability of the diffuse interface models.

Stable Manifolds and the Transition to Turbulence in Pipe Flow

Abstract: Lower-branch traveling waves and equilibria computed in pipe flow and other shear flows appear intermediate between turbulent and laminar motions. We take a step towards connecting these lower-branch solutions to transition by deriving a numerical method for finding certain special disturbances of the laminar flow in a short pipe. These special disturbances cause the disturbed velocity field to approach the lower-branch solution by evolving along its stable manifold. If the disturbance were slightly smaller, the flow would relaminarize, and if slightly larger, it would transition to a turbulent state.