Showing papers in "arXiv: Fluid Dynamics in 2011"

Thermodynamically Consistent, Frame Indifferent Diffuse Interface Models for Incompressible Two-Phase Flows with Different Densities

Abstract: A new diffuse interface model for a two-phase flow of two incompressible fluids with different densities is introduced using methods from rational continuum mechanics. The model fulfills local and global dissipation inequalities and is frame indifferent. Moreover, it is generalized to situations with a soluble species. Using the method of matched asymptotic expansions we derive various sharp interface models in the limit when the interfacial thickness tends to zero. Depending on the scaling of the mobility in the diffusion equation we either derive classical sharp interface models or models where bulk or surface diffusion is possible in the limit. In the latter case a new term resulting from surface diffusion appears in the momentum balance at the interface. Finally, we show that all sharp interface models fulfill natural energy inequalities.

The Spatial Variability of Vehicle Densities as Determinant of Urban Network Capacity

Abstract: Due to the complexity of the traffic flow dynamics in urban road networks, most quantitative descriptions of city traffic so far are based on computer simulations. This contribution pursues a macroscopic (fluid-dynamic) simulation approach, which facilitates a simple simulation of congestion spreading in cities. First, we show that a quantization of the macroscopic turning flows into units of single vehicles is necessary to obtain realistic fluctuations in the traffic variables, and how this can be implemented in a fluid-dynamic model. Then, we propose a new method to simulate destination flows without the requirement of individual route assignments. Combining both methods allows us to study a variety of different simulation scenarios. These reveal fundamental relationships between the average flow, the average density, and the variability of the vehicle densities. Considering the inhomogeneity of traffic as an independent variable can eliminate the scattering of congested flow measurements. The variability also turns out to be a key variable of urban traffic performance. Our results can be explained through the number of full links of the road network, and approximated by a simple analytical formula.

Boundary layer structure in turbulent thermal convection and its consequences for the required numerical resolution

Abstract: Results on the Prandtl-Blasius type kinetic and thermal boundary layer thicknesses in turbulent Rayleigh-B\'enard convection in a broad range of Prandtl numbers are presented. By solving the laminar Prandtl-Blasius boundary layer equations, we calculate the ratio of the thermal and kinetic boundary layer thicknesses, which depends on the Prandtl number Pr only. It is approximated as $0.588Pr^{-1/2}$ for $Pr\ll Pr^*$ and as $0.982 Pr^{-1/3}$ for $Pr^*\ll\Pr$, with $Pr^*= 0.046$. Comparison of the Prandtl--Blasius velocity boundary layer thickness with that evaluated in the direct numerical simulations by Stevens, Verzicco, and Lohse (J. Fluid Mech. 643, 495 (2010)) gives very good agreement. Based on the Prandtl--Blasius type considerations, we derive a lower-bound estimate for the minimum number of the computational mesh nodes, required to conduct accurate numerical simulations of moderately high (boundary layer dominated) turbulent Rayleigh-B\'enard convection, in the thermal and kinetic boundary layers close to bottom and top plates. It is shown that the number of required nodes within each boundary layer depends on Nu and Pr and grows with the Rayleigh number Ra not slower than $\sim\Ra^{0.15}$. This estimate agrees excellently with empirical results, which were based on the convergence of the Nusselt number in numerical simulations.

Interface-resolved DNS of vertical particulate channel flow in the turbulent regime

Abstract: We have conducted a direct numerical simulation (DNS) study of dilute turbulent particulate flow in a vertical plane channel, considering thousands of finite-size rigid particles with resolved phase interfaces. The particle diameter corresponds to approximately 11 wall units and their terminal Reynolds number is set to 136. The fluid flow with bulk Reynolds number 2700 is directed upward, which maintains the particles suspended upon average. Two density ratios were simulated, differing by a factor of 4.5. The corresponding Stokes numbers of the two flow cases were O(10) in the near-wall region and O(1) in the outer flow. We have observed the formation of large-scale elongated streak-like structures with streamwise dimensions of the order of 8 channel half-widths and cross-stream dimensions of the order of one half-width. At the same time, we have found no evidence of significant formation of particle clusters, which suggests that the large structures are due to an intrinsic instability of the flow, triggered by the presence of the particles. It was found that the mean fluid velocity profile tends towards a concave shape, and the turbulence intensity as well as the normal stress anisotropy are strongly increased. The effect of varying the Stokes number while maintaining the buoyancy, particle size and volume fraction constant was relatively weak.

The Significance of Simple Invariant Solutions in Turbulent Flows

Abstract: Recent remarkable progress in computing power and numerical analysis is enabling us to fill a gap in the dynamical systems approach to turbulence. One of the significant advances in this respect has been the numerical discovery of simple invariant sets, such as nonlinear equilibria and periodic solutions, in well-resolved Navier--Stokes flows. This review describes some fundamental and practical aspects of dynamical systems theory for the investigation of turbulence, focusing on recently found invariant solutions and their significance for the dynamical and statistical characterization of low-Reynolds-number turbulent flows. It is shown that the near-wall regeneration cycle of coherent structures can be reproduced by such solutions. The typical similarity laws of turbulence, i.e. the Prandtl wall law and the Kolmogorov law for the viscous range, as well as the pattern and intensity of turbulence-driven secondary flow in a square duct can also be represented by these simple invariant solutions.

Stability and turbulent transport in rotating shear flows: prescription from analysis of cylindrical and plane Couette flows data

Abstract: This paper provides a prescription for the turbulent viscosity in rotating shear flows for use e.g. in geophysical and astrophysical contexts. This prescription is the result of the detailed analysis of the experimental data obtained in several studies of the transition to turbulence and turbulent transport in Taylor-Couette flow. We first introduce a new set of control parameters, based on dynamical rather than geometrical considerations, so that the analysis applies more naturally to rotating shear flows in general and not only to Taylor-Couette flow. We then investigate the transition thresholds in the supercritical and the subcritical regime in order to extract their general dependencies on the control parameters. The inspection of the mean profiles provides us with some general hints on the mean to laminar shear ratio. Then the examination of the torque data allows us to propose a decomposition of the torque dependence on the control parameters in two terms, one completely given by measurements in the case where the outer cylinder is at rest, the other one being a universal function provided here from experimental fits. As a result, we obtain a general expression for the turbulent viscosity and compare it to existing prescription in the literature. Finally, throughout all the paper we discuss the influence of additional effects such as stratification or magnetic fields.

The decay of turbulence generated by a class of multi-scale grids

Abstract: A new experimental investigation of decaying turbulence generated by a low-blockage space-filling fractal square grid is presented. We find agreement with previous works by Seoud & Vassilicos (2007) and Mazellier & Vassilicos (2010) but also extend the length of the assessed decay region and consolidate the results by repeating the experiments with different probes of increased spatial resolution. It is confirmed that this moderately high Reynolds number Re{\lambda} turbulence (up to Re{\lambda} {\simeq}350 here) does not follow the classical high Reynolds number scaling of the dissipation rate {\epsilon} ~ u'^3/L and does not obey the equivalent proportionality between the Taylor-based Reynolds number Re{\lambda} and the ratio of integral scale L to Taylor micro-scale {\lambda}. Instead we observe an approximate proportionality between L and {\lambda} during decay. This non-classical behaviour is investigated by studying how the energy spectra evolve during decay and examining how well they can be described by self-preserving single-length scale forms. A detailed study of homogeneity and isotropy is also presented which reveals the presence of transverse energy transport and pressure transport in the part of the turbulence decay region where we take data (even though previous studies found mean flow and turbulence intensity profiles to be approximately homogeneous in much of the decay region). The exceptionally fast turbulence decay observed in the part of the decay region where we take data is consistent with the non-classical behaviour of the dissipation rate. Measurements with a regular square mesh grid as well as comparisons with active grid experiments by Mydlarski & Warhaft (1996) and Kang, Chester & Meveneau (2003) are also presented to highlight the similarities and differences between these turbulent flows and the turbulence generated by our fractal square grid.

Optimal feeding is optimal swimming for all P\'eclet numbers

Abstract: Cells swimming in viscous fluids create flow fields which influence the transport of relevant nutrients, and therefore their feeding rate. We propose a modeling approach to the problem of optimal feeding at zero Reynolds number. We consider a simplified spherical swimmer deforming its shape tangentially in a steady fashion (so-called squirmer). Assuming that the nutrient is a passive scalar obeying an advection-diffusion equation, the optimal use of flow fields by the swimmer for feeding is determined by maximizing the diffusive flux at the organism surface for a fixed rate of energy dissipation in the fluid. The results are obtained through the use of an adjoint-based numerical optimization implemented by a Legendre polynomial spectral method. We show that, to within a negligible amount, the optimal feeding mechanism consists in putting all the energy expended by surface distortion into swimming - so-called treadmill motion - which is also the solution maximizing the swimming efficiency. Surprisingly, although the rate of feeding depends strongly on the value of the Peclet number, the optimal feeding stroke is shown to be essentially independent of it, which is confirmed by asymptotic analysis. Within the context of steady actuation, optimal feeding is therefore found to be equivalent to optimal swimming for all Peclet numbers.

The origin of hysteresis in the flag instability

Abstract: The flapping flag instability occurs when a flexible cantilevered plate is immersed in a uniform airflow. To this day, the nonlinear aspects of this aeroelastic instability are largely unknown. In particular, experiments in the literature all report a large hysteresis loop, while the bifurcation in numerical simulations is either supercritical or subcritical with a small hysteresis loop. In this paper, this discrepancy is addressed. First weakly nonlinear stability analyses are conducted in the slender-body and two-dimensional limits, and second new experiments are performed with flat and curved plates. The discrepancy is attributed to inevitable planeity defects of the plates in the experiments.

Using multiscale norms to quantify mixing and transport

Abstract: Mixing is relevant to many areas of science and engineering, including the pharmaceutical and food industries, oceanography, atmospheric sciences, and civil engineering In all these situations one goal is to quantify and often then to improve the degree of homogenisation of a substance being stirred, referred to as a passive scalar or tracer A classical measure of mixing is the variance of the concentration of the scalar, which can be related to the $L^2$ norm of the concentration field Recently other norms have been used to quantify mixing, in particular the mix-norm as well as negative Sobolev norms These norms have the advantage that unlike variance they decay even in the absence of diffusion, and their decay corresponds to the flow being mixing in the sense of ergodic theory General Sobolev norms weigh scalar gradients differently, and are known as multiscale norms for mixing We review the applications of such norms to mixing and transport, and show how they can be used to optimise the stirring and mixing of a decaying passive scalar We then review recent work on the less-studied case of a continuously-replenished scalar field --- the source-sink problem In that case the flows that optimally reduce the norms are associated with transport rather than mixing: they push sources onto sinks, and vice versa

Minimal seeds for shear flow turbulence: using nonlinear transient growth to touch the edge of chaos

Abstract: We propose a general strategy for determining the minimal finite amplitude isturbance to trigger transition to turbulence in shear flows. This involves constructing a variational problem that searches over all disturbances of fixed initial amplitude, which respect the boundary conditions, incompressibility and the Navier--Stokes equations, to maximise a chosen functional over an asymptotically long time period. The functional must be selected such that it identifies turbulent velocity fields by taking significantly enhanced values compared to those for laminar fields. We illustrate this approach using the ratio of the final to initial perturbation kinetic energies (energy growth) as the functional and the energy norm to measure amplitudes in the context of pipe flow. Our results indicate that the variational problem yields a smooth converged solution providing the amplitude is below the threshold amplitude for transition. This optimal is the nonlinear analogue of the well-studied (linear) transient growth optimal. At and above this threshold, the optimising search naturally seeks out disturbances that trigger turbulence by the end of the period, and convergence is then practically impossible. The first disturbance found to trigger turbulence as the amplitude is increased identifies the `minimal seed' for the given geometry and forcing (Reynolds number). We conjecture that it may be possible to select a functional such that the converged optimal below threshold smoothly converges to the minimal seed at threshold. This seems at least approximately true for our choice of energy growth functional and the pipe flow geometry chosen here.

Isotropisation at small scales of rotating helically-driven turbulence

Abstract: We present numerical evidence of how three-dimensionalization occurs at small scale in rotating turbulence with Beltrami (ABC) forcing, creating helical flow. The Zeman scale $\ell_{\Omega}$ at which the inertial and eddy turn-over times are equal is more than one order of magnitude larger than the dissipation scale, with the relevant domains (large-scale inverse cascade of energy, dual regime in the direct cascade of energy $E$ and helicity $H$, and dissipation) each moderately resolved. These results stem from the analysis of a large direct numerical simulation on a grid of $3072^3$ points, with Rossby and Reynolds numbers respectively equal to 0.07 and $2.7\times 10^4$. At scales smaller than the forcing, a helical wave-modulated inertial law for the energy and helicity spectra is followed beyond $\ell_{\Omega}$ by Kolmogorov spectra for $E$ and $H$. Looking at the two-dimensional slow manifold, we also show that the helicity spectrum breaks down at $\ell_{\Omega}$, a clear sign of recovery of three-dimensionality in the small scales.

Highly focused supersonic microjets

Abstract: The paper describes the production of thin, focused microjets with velocities up to 850 m/s by the rapid vaporization of a small mass of liquid in an open liquid-filled capillary. The vaporization is caused by the absorption of a low-energy laser pulse. A likely explanation of the observed phenomenon is based on the impingement of the shock wave caused by the nearly-instantaneous vaporization on the free surface of the liquid. An experimental study of the dependence of the jet velocity on several parameters is conducted, and a semi-empirical relation for its prediction is developed. The coherence of the jets, their high velocity and good reproducibility and controllability are unique features of the system described. A possible application is to the development of needle-free drug injection systems which are of great importance for global health care.

Force and torque acting on particles in a transitionally rough open channel flow

Abstract: Direct numerical simulation of open channel flow over a geometrically rough wall has been performed at a bulk Reynolds number of approximately 2900. The wall consisted of a layer of spheres in a square arrangement. Two cases have been considered. In the first case the spheres are small (with diameter equivalent to 10.7 wall units) and the limit of the hydraulically smooth flow regime is approached. In the second case the spheres are more than three times larger (49.3 wall units) and the flow is in the transitionally rough flow regime. Special emphasis is given on the characterisation of the force and torque acting on a particle due to the turbulent flow. It is found that in both cases the mean drag, lift and spanwise torque are to a large extent produced at the top region of the particle surface. The intensity of the particle force fluctuations is significantly larger in the large-sphere case, while the trend differs for the fluctuations of the individual components of the torque. A simplified model is used to show that the torque fluctuations might be explained by the spheres acting as a filter with respect to the size of the flow scales which can effectively generate torque fluctuations. Fluctuations of both force and torque are found to exhibit strongly non-Gaussian probability density functions with particularly long tails, an effect which is more pronounced in the small-sphere case. Some implications of the present results for sediment erosion are briefly discussed.

Optimal Taylor-Couette turbulence

Abstract: Strongly turbulent Taylor-Couette flow with independently rotating inner and outer cylinders with a radius ratio of \eta = 0.716 is experimentally studied. From global torque measurements, we analyse the dimensionless angular velocity flux Nu_\omega(Ta, a) as a function of the Taylor number Ta and the angular velocity ratio a = -\omega_o/\omega_i in the large-Taylor-number regime 10^{11} \lesssim Ta \lesssim 10^{13}. We analyse the data with the common power-law ansatz for the dimensionless angular velocity transport flux Nu_\omega(Ta, a) = f(a)Ta^\gamma, with an amplitude f(a) and an exponent \gamma. The data are consistent with one effective exponent \gamma = 0.39\pm0.03 for all a. The amplitude of the angular velocity flux f(a) = Nu_\omega(Ta, a)/Ta^0.39 is measured to be maximal at slight counter-rotation, namely at an angular velocity ratio of a_opt = 0.33\pm0.04. This value is theoretically interpreted as the result of a competition between the destabilizing inner cylinder rotation and the stabilizing but shear-enhancing outer cylinder counter-rotation. With the help of laser Doppler anemometry, we provide angular velocity profiles and identify the radial position r_n of the neutral line. While for moderate counter-rotation -0.40 \omega_i \lesssim \omega_o < 0, the neutral line still remains close to the outer cylinder and the probability distribution function (p.d.f.) of the bulk angular velocity is observed to be monomodal. For stronger counter-rotation the neutral line is pushed inwards towards the inner cylinder; in this regime the p.d.f. of the bulk angular velocity becomes bimodal, reflecting intermittent bursts of turbulent structures beyond the neutral line into the outer flow domain, which otherwise is stabilized by the counter-rotating outer cylinder. Finally, a hypothesis is offered allowing a unifying view for all these various results.

The end of nanochannels

Abstract: Current theories of nanochannel flow impose no upper bound on flow rates, and predict friction through nanochannels can be vanishingly small. We reassess neglecting channel entry effects in extremely long channels and find violations at the nanoscale. Even in frictionless nanochannels, end effects provide a finite amount of friction. Hence, the speed at which nanochannels transport liquids is limited. Flow-rate and slip length measurements are reevaluated using calculations which include end-effect friction. End effects are critical for the design of new technological devices and to understand biological transport.

Experimental and numerical study of mean zonal flows generated by librations of a rotating spherical cavity

Abstract: We study both experimentally and numerically the steady zonal flow generated by longitudinal librations of a spherical rotating container. This study follows the recent weakly nonlinear analysis of Busse (2010), developed in the limit of small libration frequency - rotation rate ratio, and large libration frequency - spin-up time product. Using PIV measurements as well as results from axisymmetric numerical simulations, we confirm quantitatively the main features of Busse's analytical solution: the zonal flow takes the form of a retrograde solid body rotation in the fluid interior, which does not depend on the libration frequency nor on the Ekman number, and which varies as the square of the amplitude of excitation. We also report the presence of an unpredicted prograde flow at the equator near the outer wall.

Angular momentum transport and turbulence in laboratory models of Keplerian flows

Abstract: We present angular momentum transport (torque) measurements in two recent experimental studies of the turbulent flow between independently rotating cylinders. In addition to these studies, we reanalyze prior torque measurements to expand the range of control parameters for the experimental Taylor-Couette flows. We find that the torque may be described as a product of functions that depend only on the Reynolds number, which describes the turbulent driving intensity, and the rotation number, which characterizes the effects of global rotation. For a given Reynolds number, the global angular momentum transport for Keplerian-like flow profiles is approximately 14% of the maximum achievable transport rate. We estimate that this level of transport would produce an accretion rate of $\dot{M}/\dot{M_0} \sim 10^{-3}$ in astrophysical disks. We argue that this level of transport from hydrodynamics alone could be significant.

Hydrodynamic models of self-organized dynamics: derivation and existence theory

Abstract: This paper is concerned with the derivation and analysis of hydrodynamic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. The starting point is the kinetic model considered in earlier work of Degond & Motsch with the addition of an attraction-repulsion interaction potential. Introducing different scalings than in Degond & Motsch, the non-local effects of the alignment and attraction-repulsion interactions can be kept in the hydrodynamic limit and result in extra pressure, viscosity terms and capillary force. The systems are shown to be symmetrizable hyperbolic systems with viscosity terms. A local-in-time existence result is proved in the 2D case for the viscous model and in the 3D case for the inviscid model. The proof relies on the energy method.

Particle Jet Formation During Explosive Dispersal

Abstract: Previous experimental studies have shown that when a layer of solid particles is explosively dispersed, the particles often develop a non-uniform spatial distribution. The instabilities within the particle bed and at the particle layer interface likely form on the timescale of the shock propagation through the particles. The mesoscale perturbations are manifested at later times in experiments by the formation of coherent clusters of particles or jet-like particle structures, which are aerodynamically stable. A number of different mechanisms likely contribute to the jet formation including shock fracturing of the particle bed and particle-particle interactions in the early stages of the dense gas-particle flow. Aerodynamic wake effects at later times contribute to maintaining the stability of the jets. The experiments shown in this fluid dynamics video were carried out in either spherical or cylindrical geometry and illustrate the formation of particle jets during the explosive dispersal process. The number of jet-like structures that are generated during the dispersal of a dry powder bed is compared with the number formed during the dispersal of the same volume of water. The liquid dispersal generates a larger number of jets, but they fragment and dissipate sooner. When the particle bed is saturated with water and explosively dispersed, the number of particle jets formed is larger than both the dry powder and pure water charges. More extensive experiments that explore the effect of particle size, density and the mass ratio of explosive to particles on the susceptibility for jet formation are reported in Frost et al. (Proc. of 23rd ICDERS, Irvine, CA, 2011).

On the velocity distribution in homogeneous isotropic turbulence: correlations and deviations from Gaussianity

Abstract: We investigate the single-point probability density function of the velocity in three-dimensional stationary and decaying homogeneous isotropic turbulence. To this end we apply the statistical framework of the Lundgren-Monin-Novikov hierarchy combined with conditional averaging, identifying the quantities that determine the shape of the probability density function. In this framework the conditional averages of the rate of energy dissipation, the velocity diffusion and the pressure gradient with respect to velocity play a key role. Direct numerical simulations of the Navier-Stokes equation are used to complement the theoretical results and assess deviations from Gaussianity.

Microfluidic propulsion by the metachronal beating of magnetic artificial cilia: a numerical analysis

Abstract: In this work we study the effect of metachronal waves on the flow created by magnetically-driven plate-like artificial cilia in microchannels using numerical simulations. The simulations are performed using a coupled magneto-mechanical solid-fluid computational model that captures the physical interactions between the fluid flow, ciliary deformation and applied magnetic field. When a rotating magnetic field is applied to super-paramagnetic artificial cilia, they mimic the asymmetric motion of natural cilia, consisting of an effective and recovery stroke. When a phase-difference is prescribed between neighbouring cilia, metachronal waves develop. Due to the discrete nature of the cilia, the metachronal waves change direction when the phase difference becomes sufficiently large, resulting in antiplectic as well as symplectic metachrony. We show that the fluid flow created by the artificial cilia is significantly enhanced in the presence of metachronal waves and that the fluid flow becomes unidirectional. Antiplectic metachrony is observed to lead to a considerable enhancement in flow compared to symplectic metachrony, when the cilia spacing is small. Obstruction of flow in the direction of the effective stroke for the case of symplectic metachrony was found to be the key mechanism that governs this effect.

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

Abstract: We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow v.s. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque $G=T/ u^2L$, where $L$ is the cylinder length, scales with Reynolds number $R$ and gap width $\eta$, $G=1.46 \eta^{3/2} (1-\eta)^{-7/4}R^{3/2}$. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: $$ G=0.50\frac{\eta^{2}}{(1-\eta)^{3/2}}\frac{R^{2}}{\ln[\eta^2(1-\eta)R^ 2/10^4]^{3/2}}.$$ These predictions are found to be in excellent agreement with available experimental data. Predictions for scaling of velocity fluctuations are also provided.

Surface nanobubbles as a function of gas type

Abstract: We experimentally investigate the nucleation of surface nanobubbles on PFDTS-coated silicon as a function of the specific gas dissolved in the water. In each case we restrict ourselves to equilibrium conditions ($c=100%$, $T_{liquid} = T_{substrate}$). Not only is nanobubble nucleation a strong function of gas type, but there also exists an optimal system temperature of $\sim 35-40\mathrm{^oC}$ where nucleation is maximized, which is weakly dependent on gas type. We also find that contact angle is a function of nanobubble radius of curvature for all gas types investigated. Fitting this data allows us to describe a line tension which is dependent on the type of gas, indicating that the nanobubbles are sat on top of adsorbed gas molecules. The average line tension was $\tau \sim -0.8 \mathrm{nN}$.

Finite volume and pseudo-spectral schemes for the fully nonlinear 1D Serre equations

Abstract: After we derive the Serre system of equations of water wave theory from a generalized variational principle, we present some of its structural properties. We also propose a robust and accurate finite volume scheme to solve these equations in one horizontal dimension. The numerical discretization is validated by comparisons with analytical, experimental data or other numerical solutions obtained by a highly accurate pseudo-spectral method.

Stokes flow near the contact line of an evaporating drop

Abstract: The evaporation of sessile drops in quiescent air is usually governed by vapour diffusion. For contact angles below $90^\circ$, the evaporative flux from the droplet tends to diverge in the vicinity of the contact line. Therefore, the description of the flow inside an evaporating drop has remained a challenge. Here, we focus on the asymptotic behaviour near the pinned contact line, by analytically solving the Stokes equations in a wedge geometry of arbitrary contact angle. The flow field is described by similarity solutions, with exponents that match the singular boundary condition due to evaporation. We demonstrate that there are three contributions to the flow in a wedge: the evaporative flux, the downward motion of the liquid-air interface and the eigenmode solution which fulfils the homogeneous boundary conditions. Below a critical contact angle of $133.4^\circ$, the evaporative flux solution will dominate, while above this angle the eigenmode solution dominates. We demonstrate that for small contact angles, the velocity field is very accurately described by the lubrication approximation. For larger contact angles, the flow separates into regions where the flow is reversing towards the drop centre.

Local transient rheological behavior of concentrated suspensions

Abstract: This paper reports experiments on the shear transient response of concentrated non-Brownian suspensions. The shear viscosity of the suspensions is measured using a wide-gap Couette rheometer equipped with a Particle Image Velocimetry (PIV) device that allows measuring the velocity field. The suspensions made of PMMA particles (31$\mu$m in diameter) suspended in a Newtonian index- and density-matched liquid are transparent enough to allow an accurate measurement of the local velocity for particle concentrations as high as 50%. In the wide-gap Couette cell, the shear induced particle migration is evidenced by the measurement of the time evolution of the flow profile. A peculiar radial zone in the gap is identified where the viscosity remains constant. At this special location, the local particle volume fraction is taken to be the mean particle concentration. The local shear transient response of the suspensions when the shear flow is reversed is measured at this point where the particle volume fraction is well defined. The local rheological measurements presented here confirm the macroscopic measurements of Gadala-Maria and Acrivos (1980). After shear reversal, the viscosity undergoes a step-like reduction, decreases slower and passes through a minimum before increasing again to reach a plateau. Upon varying the particle concentration, we have been able to show that the minimum and the plateau viscosities do not obey the same scaling law with respect to the particle volume fraction. These experimental results are consistent with the scaling predicted by Mills and Snabre (2009) and with the results of numerical simulation performed on random suspensions [Sierou and Brady (2001)]. The minimum seems to be associated with the viscosity of an isotropic suspension, or at least of a suspension whose particles do not interact through non-hydrodynamic forces, while the plateau value would correspond to the viscosity of a suspension structured by the shear where the non-hydrodynamic forces play a crucial role.

High-resolution truncated plurigaussian simulations for the characterization of heterogeneous formations

Abstract: Integrating geological concepts, such as relative positions and proportions of the different lithofacies, is of highest importance in order to render realistic geological patterns. The truncated plurigaussian simulation method provides a way of using both local and conceptual geological information to infer the distributions of the facies and then those of hydraulic parameters. The method (Le Loc'h and Galli 1994) is based on the idea of truncating at least two underlying multi-Gaussian simulations in order to create maps of categorical variable. In this manuscript we show how this technique can be used to assess contaminant migration in highly heterogeneous media. We illustrate its application on the biggest contaminated site of Switzerland. It consists of a contaminant plume located in the lower fresh water Molasse on the western Swiss Plateau. The highly heterogeneous character of this formation calls for efficient stochastic methods in order to characterize transport processes.

The role of Stewartson and Ekman layers in turbulent rotating Rayleigh-B\'enard convection

Abstract: When the classical Rayleigh-Benard (RB) system is rotated about its vertical axis roughly three regimes can be identified. In regime I (weak rotation) the large scale circulation (LSC) is the dominant feature of the flow. In regime II (moderate rotation) the LSC is replaced by vertically aligned vortices. Regime III (strong rotation) is characterized by suppression of the vertical velocity fluctuations. Using results from experiments and direct numerical simulations of RB convection for a cell with a diameter-to-height aspect ratio equal to one at $Ra \sim 10^8-10^9$ ($Pr=4-6$) and $0 \lesssim 1/Ro \lesssim 25$ we identified the characteristics of the azimuthal temperature profiles at the sidewall in the different regimes. In regime I the azimuthal wall temperature profile shows a cosine shape and a vertical temperature gradient due to plumes that travel with the LSC close to the sidewall. In regime II and III this cosine profile disappears, but the vertical wall temperature gradient is still observed. It turns out that the vertical wall temperature gradient in regimes II and III has a different origin than that observed in regime I. It is caused by boundary layer dynamics characteristic for rotating flows, which drives a secondary flow that transports hot fluid up the sidewall in the lower part of the container and cold fluid downwards along the sidewall in the top part.

Damping by branching: a bioinspiration from trees

Abstract: Man-made slender structures are known to be sensitive to high levels of vibration, due to their flexibility, which often cause irreversible damage. In nature, trees repeatedly endure large amplitudes of motion, mostly caused by strong climatic events, yet with minor or no damage in most cases. A new damping mechanism inspired by the architecture of trees is here identified and characterized in the simplest tree-like structure, a Y-shape branched structure. Through analytical and numerical analyses of a simple two-degree-of-freedom model, branching is shown to be the key ingredient in this protective mechanism that we call damping-by-branching. It originates in the geometrical nonlinearities so that it is specifically efficient to damp out large amplitudes of motion. A more realistic model, using flexible beam approximation, shows that the mechanism is robust. Finally, two bioinspired architectures are analyzed, showing significant levels of damping achieved via branching with typically 30% of the energy being dissipated in one oscillation. This concept of damping-by-branching is of simple practical use in the design of slender flexible structures.