Showing papers in "arXiv: Fluid Dynamics in 2014"

Multiphase lattice Boltzmann simulations for porous media applications -- a review

Abstract: Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

Turbulent channel flow of dense suspensions of neutrally-buoyant spheres

Abstract: Dense particle suspensions are widely encountered in many applications and in environmental flows. While many previous studies investigate their rheological properties in laminar flows, little is known on the behaviour of these suspensions in the turbulent/inertial regime. The present study aims to fill this gap by investigating the turbulent flow of a Newtonian fluid laden with solid neutrally-buoyant spheres at relatively high volume fractions in a plane channel. Direct Numerical Simulation are performed in the range of volume fractions Phi=0-0.2 with an Immersed Boundary Method used to account for the dispersed phase. The results show that the mean velocity profiles are significantly altered by the presence of a solid phase with a decrease of the von Karman constant in the log-law. The overall drag is found to increase with the volume fraction, more than one would expect just considering the increase of the system viscosity due to the presence of the particles. At the highest volume fraction here investigated, Phi=0.2, the velocity fluctuation intensities and the Reynolds shear stress are found to decrease. The analysis of the mean momentum balance shows that the particle-induced stresses govern the dynamics at high Phi and are the main responsible of the overall drag increase. In the dense limit, we therefore find a decrease of the turbulence activity and a growth of the particle induced stress, where the latter dominates for the Reynolds numbers considered here.

Phoretic self-propulsion at finite P\'eclet numbers

Abstract: Phoretic self-propulsion is a unique example of force- and torque-free motion on small scales. The classical framework describing the flow field around a particle swimming by self-diffusiophoresis neglects the advection of the solute field by the flow and assumes that the chemical interaction layer is thin compared to the particle size. In this paper we quantify and characterize the effect of solute advection on the phoretic swimming of a sphere. We first rigorously derive the regime of validity of the thin-interaction layer assumption at finite values of the P\'eclet number (Pe). Within this assumption, we solve computationally the flow around Janus phoretic particles and examine the impact of solute advection on propulsion and the flow created by the particle. We demonstrate that although advection always leads to a decrease of the swimming speed and flow stresslet at high values of the P\'eclet number, an increase can be obtained at intermediate values of Pe. This possible enhancement of swimming depends critically on the nature of the chemical interactions between the solute and the surface. We then derive an asymptotic analysis of the problem at small Pe allowing to rationalize our computational results. Our computational and theoretical analysis is accompanied by a parallel study of the role of reactive effects at the surface of the particle on swimming (Damk\"ohler number).

DNS of turbulent channel flow at very low Reynolds numbers

Abstract: Direct numerical simulations (DNS) of fully-developed turbulent channel flows for very low Reynolds numbers have been performed with a larger computational box sizes than those of existing DNS. The friction Reynolds number was decreased down to 60, where the friction Reynolds number is based on the friction velocity and the channel half width. When the Reynolds number was decreased to 60 with small computational box size, the flow became laminar. Using a large box, we found that a localized turbulence was observed to sustain in the form of periodic oblique band. This type of locally disordered flow is similar to a equilibrium turbulent puff in a transitional pipe flow. Various turbulence statistics such as turbulence intensities, vorticity fluctuations and Reynolds stresses are provided. Especially, their near-wall asymptotic behavior and budget terms of turbulence kinetic energy were discussed with respect to the Reynolds-number dependence and an influence of the computational box size. Other detailed characteristics associated with the turbulence structures were also presented and discussed.

Large-scale vortices in rapidly rotating Rayleigh-B\'enard convection

Abstract: Using numerical simulations of rapidly rotating Boussinesq convection in a Cartesian box, we study the formation of long-lived, large-scale, depth-invariant coherent structures. These structures, which consist of concentrated cyclones, grow to the horizontal size of the box, with velocities significantly larger than the convective motions. We vary the rotation rate, the thermal driving and the aspect ratio in order to determine the domain of existence of these large-scale vortices (LSV). We find that two conditions are required for their formation. First, the Rayleigh number, a meaure of the thermal driving, must be several times its value at the linear onset of convection; this corresponds to Reynolds numbers, based on the convective velocity and the box depth, $\gtrsim 100$. Second, the rotational constraint on the convective structures must be strong. This requires that the local Rossby number, based on the convective velocity and the horizontal convective scale, $\lesssim 0.15$. Simulations in which certain wavenumbers are artificially suppressed in spectral space suggest that the LSV are produced by the interactions of small-scale, depth-dependent convective motions. The presence of LSV significantly reduces the efficiency of the convective heat transport.

Closed-loop separation control using machine learning

Abstract: A novel, model free, approach to experimental closed-loop flow control is implemented on a separated flow. Feedback control laws are generated using genetic programming where they are optimized using replication, mutation and cross-over of best performing laws to produce a new generation of candidate control laws. This optimization process is applied automatically to a backward-facing step flow at Re=1350, controlled by a slotted jet, yielding an effective control law. Convergence criterion are suggested. The law is able to produce effective action even with major changes in the flow state, demonstrating its robustness. The underlying physical mechanisms leveraged by the law are analyzed and discussed. Contrary to traditional periodic forcing of the shear layer, this new control law plays on the physics of the recirculation area downstream the step. While both control actions are fundamentally different they still achieve the same level of effectiveness. Furthermore the new law is also potentially easier and cheaper to implement actuator wise.

Nanoscale simulation of shale transport properties using the lattice Boltzmann method: permeability and diffusivity

Abstract: Porous structures of shales are reconstructed based on scanning electron microscopy (SEM) images of shale samples from Sichuan Basin, China. Characterization analyzes of the nanoscale reconstructed shales are performed, including porosity, pore size distribution, specific surface area and pore connectivity. The multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) fluid flow model and single-relaxation-time (SRT) LBM diffusion model are adopted to simulate the fluid flow and Knudsen diffusion process within the reconstructed shales, respectively. Tortuosity, intrinsic permeability and effective Knudsen diffusivity are numerically predicted. The tortuosity is much higher than that commonly employed in Bruggeman equation. Correction of the intrinsic permeability by taking into consideration the contribution of Knudsen diffusion, which leads to the apparent permeability, is performed. The correction factor under different Knudsen number and pressure are estimated and compared with existing corrections reported in the literature. For the wide pressure range under investigation, the correction factor is always greater than 1, indicating the Knudsen diffusion always plays a role on the transport mechanisms of shale gas in shales studied in the present study. Most of the values of correction factor are located in the transition regime, with no Darcy flow regime observed.

Sedimentation of a dilute suspension of rigid spheres at intermediate Galileo numbers: the effect of clustering upon the particle motion

Abstract: Direct numerical simulation of the gravity-induced settling of finite-size particles in triply-periodic domains has been performed under dilute conditions. For a single solid-to-fluid density ratio of 1.5 we have considered two values of the Galileo number corresponding to steady vertical motion (Ga=121) and to steady oblique motion (Ga=178) in the case of one isolated sphere. For the multi-particle system we observe strong particle clustering only in the latter case. The geometry and time scales related to clustering are determined from Voronoi tesselation and particle-conditioned averaging. As a consequence of clustering, the average particle settling velocity is increased by 12% as compared to the value of an isolated sphere; such a collective effect is not observed in the non-clustering case. By defining a local (instantaneous) fluid velocity average in the vicinity of the finite-size particles it is shown that the observed enhancement of the settling velocity is due to the fact that the downward fluid motion (with respect to the global average) which is induced in the cluster regions is preferentially sampled by the particles. It is further observed that the variance of the particle velocity is strongly enhanced in the clustering case. With the aid of a decomposition of the particle velocity it is shown that this increase is due to enhanced fluid velocity fluctuations (due to clustering) in the vicinity of the particles. Finally, we discuss a possible explanation for the observation of a critical Galileo number marking the onset of clustering under dilute conditions.

Mobilization of a trapped non-wetting fluid from a three-dimensional porous medium

Abstract: We use confocal microscopy to directly visualize the formation and complex morphologies of trapped non-wetting fluid ganglia within a model 3D porous medium. The wetting fluid continues to flow around the ganglia after they form; this flow is characterized by a capillary number, Ca. We find that the ganglia configurations do not vary for small Ca; by contrast, as Ca is increased above a threshold value, the largest ganglia start to become mobilized and are ultimately removed from the medium. By combining our 3D visualization with measurements of the bulk transport, we show that this behavior can be quantitatively understood by balancing the viscous forces exerted on the ganglia with the pore-scale capillary forces that keep them trapped within the medium. Our work thus helps elucidate the fluid dynamics underlying the mobilization of a trapped non-wetting fluid from a 3D porous medium.

Waving transport and propulsion in a generalized Newtonian fluid

Abstract: Cilia and flagella are hair-like appendages that protrude from the surface of a variety of eukaryotic cells and deform in a wavelike fashion to transport fluids and propel cells. Motivated by the ubiquity of non-Newtonian fluids in biology, we address mathematically the role of shear-dependent viscosities on both the waving flagellar locomotion and ciliary transport by metachronal waves. Using a two-dimensional waving sheet as model for the kinematics of a flagellum or an array of cilia, and allowing for both normal and tangential deformation of the sheet, we calculate the flow field induced by a small-amplitude deformation of the sheet in a generalized Newtonian Carreau fluid up to order four in the dimensionless waving amplitude. Shear-thinning and shear-thickening fluids are seen to always induce opposite effects. When the fluid is shear-thinning, the rate of working of the sheet against the fluid is always smaller than in the Newtonian fluid, and the largest gain is obtained for antiplectic metachronal waves. Considering a variety of deformation kinematics for the sheet, we further show that in all cases transport by the sheet is more efficiency in a shear-thinning fluid, and in most cases the transport speed in the fluid is also increased. Comparing the order of magnitude of the shear-thinning contributions with past work on elastic effects as well as the magnitude of the Newtonian contributions, our theoretical results suggest that the impact of shear-dependent viscosities on transport could play a major biological role.

Modeling size segregation of granular materials: the roles of segregation, advection and diffusion

Abstract: Predicting segregation of granular materials composed of different-sized particles is a challenging problem. In this paper, we develop and implement a theoretical model that captures the interplay between advection, segregation, and diffusion in size bidisperse granular materials. The fluxes associated with these three driving factors depend on the underlying kinematics, whose characteristics play key roles in determining particle segregation configurations. Unlike previous models for segregation, our model uses parameters based on kinematic measures from discrete element method simulations instead of arbitrarily adjustable fitting parameters, and it achieves excellent quantitative agreement with both experimental and simulation results when applied to quasi-two-dimensional bounded heaps. The model yields two dimensionless control parameters, both of which are only functions of physically control parameters (feed rate, particle sizes, and system size) and kinematic parameters (diffusion coefficient, flowing layer depth, and percolation velocity). The P\'eclet number, $Pe$, captures the interplay of advection and diffusion, and the second dimensionless parameter, $\Lambda$, describes the interplay between segregation and advection. A parametric study of $\Lambda$ and $Pe$ demonstrates how the particle segregation configuration depends on the interplay of advection, segregation, and diffusion. The model can be readily adapted to other flow geometries.

Asymmetric steady streaming as a mechanism for acoustic propulsion of rigid bodies

Abstract: Recent experiments showed that standing acoustic waves could be exploited to induce self-propulsion of rigid metallic particles in the direction perpendicular to the acoustic wave. We propose in this paper a physical mechanism for these observations based on the interplay between inertial forces in the fluid and the geometrical asymmetry of the particle shape. We consider an axisymmetric rigid near-sphere oscillating in a quiescent fluid along a direction perpendicular to its symmetry axis. The kinematics of oscillations can be either prescribed or can result dynamically from the presence of an external oscillating velocity field. Steady streaming in the fluid, the inertial rectification of the time-periodic oscillating flow, generates steady stresses on the particle which, in general, do not average to zero, resulting in a finite propulsion speed along the axis of the symmetry of the particle and perpendicular to the oscillation direction. Our derivation of the propulsion speed is obtained at leading order in the Reynolds number and the deviation of the shape from that of a sphere. The results of our model are consistent with the experimental measurements, and more generally explains how time periodic forcing from an acoustic field can be harnessed to generate autonomous motion.

Direct numerical simulation of turbulent channel flow over porous walls

Abstract: We perform direct numerical simulations (DNS) of a turbulent channel flow over porous walls. In the fluid region the flow is governed by the incompressible Navier--Stokes (NS) equations, while in the porous layers the Volume-Averaged Navier--Stokes (VANS) equations are used, which are obtained by volume-averaging the microscopic flow field over a small volume that is larger than the typical dimensions of the pores. In this way the porous medium has a continuum description, and can be specified without the need of a detailed knowledge of the pore microstructure by indipendently assigning permeability and porosity. At the interface between the porous material and the fluid region, momentum-transfer conditions are applied, in which an available coefficient related to the unknown structure of the interface can be used as an error estimate. To set up the numerical problem, the velocity-vorticity formulation of the coupled NS and VANS equations is derived and implemented in a pseudo-spectral DNS solver. Most of the simulations are carried out at $Re_\tau=180$ and consider low-permeability materials; a parameter study is used to describe the role played by permeability, porosity, thickness of the porous material, and the coefficient of the momentum-transfer interface conditions. Among them permeability, even when very small, is shown to play a major role in determining the response of the channel flow to the permeable wall. Turbulence statistics and instantaneous flow fields, in comparative form to the flow over a smooth impermeable wall, are used to understand the main changes introduced by the porous material. A simulations at higher Reynolds number is used to illustrate the main scaling quantities.

A hyperbolic model for viscous Newtonian flows

Abstract: We discuss a pure hyperbolic alternative to the Navier-Stokes equations, which are of parabolic type. As a result of the substitution of the concept of the viscosity coefficient by a microphysics-based temporal characteristic, particle settled life (PSL) time, it becomes possible to formulate a model for viscous fluids in a form of first order hyperbolic partial differential equations. Moreover, the concept of PSL time allows the use of the same model for flows of viscous fluids (Newtonian or non-Newtonian) as well as irreversible deformation of solids. In the theory presented, a continuum is interpreted as a system of material particles connected by bonds; the internal resistance to flow is interpreted as elastic stretching of the particle bonds; and a flow is a result of bond destructions and rearrangements of particles. Finally, we examine the model for simple shear flows, arbitrary incompressible and compressible flows of Newtonian fluids and demonstrate that Newton's viscous law can be obtained in the framework of the developed hyperbolic theory as a steady-state limit. A basic relation between the viscosity coefficient, PSL time, and the shear sound velocity is also obtained.

A generalised-Lagrangian-mean model of the interactions between near-inertial waves and mean flow

Abstract: Wind forcing of the ocean generates a spectrum of inertia-gravity waves that is sharply peaked near the local inertial (or Coriolis) frequency. The corresponding near-inertial waves (NIWs) are highly energetic and play a significant role in the slow, large-scale dynamics of the ocean. To analyse this role, we develop a new model of the nondissipative interactions between NIWs and balanced motion. The model is derived using the generalised-Lagrangian-mean (GLM) framework (specifically, the glm variant of Soward & Roberts (2010)), taking advantage of the time-scale separation between the two types of motion to average over the short NIW period. We combine Salmon's (2013) variational formulation of GLM with Whitham averaging to obtain a system of equations governing the joint evolution of NIWs and mean flow. Assuming that the mean flow is geostrophically balanced reduces this system to a simple model coupling Young & Ben Jelloul's (1997) equation for NIWs with a modified quasi-geostrophic equation. In this coupled model, the mean flow affects the NIWs through advection and refraction; conversely, the NIWs affect the mean flow by modifying the potential-vorticity inversion - the relation between advected potential vorticity and advecting mean velocity - through a quadratic wave term, consistent with the GLM results of Buhler & McIntyre (1998). The coupled model is Hamiltonian and its conservation laws, for wave action and energy in particular, prove illuminating: on their basis, we identify a new interaction mechanism whereby NIWs forced at large scales extract energy from the balanced flow as their horizontal scale is reduced by differential advection and refraction so that their potential energy increases. A rough estimate suggests that this mechanism could provide a significant sink of energy for mesoscale motion and play a part in the global energetics of the ocean.

Shape-dependence of particle rotation in isotropic turbulence

Abstract: We consider the rotation of neutrally buoyant axisymmetric particles suspended in isotropic turbulence. Using laboratory experiments as well as numerical and analytical calculations, we explore how particle rotation depends upon particle shape. We find that shape strongly affects orientational trajectories, but that it has negligible effect on the variance of the particle angular velocity. Previous work has shown that shape significantly affects the variance of the tumbling rate of axisymmetric particles. It follows that shape affects the spinning rate in a way that is, on average, complementary to the shape-dependence of the tumbling rate. We confirm this relationship using direct numerical simulations, showing how tumbling rate and spinning rate variances show complementary trends for rod-shaped and disk-shaped particles. We also consider a random but non-turbulent flow. This allows us to explore which of the features observed for rotation in turbulent flow are due to the effects of particle alignment in vortex tubes.

An optimisation approach for analysing nonlinear stability with transition to turbulence in fluids as an exemplar

Abstract: This article introduces, and reviews recent work using, a simple optimisation technique for analysing the nonlinear stability of a state in a dynamical system. The technique can be used to identify the most efficient way to disturb a system such that it transits from one stable state to another. The key idea is introduced within the framework of a finite-dimensional set of ordinary differential equations (ODEs) and then illustrated for a very simple system of 2 ODEs which possesses bistability. Then the transition to turbulence problem in fluid mechanics is used to show how the technique can be formulated for a spatially-extended system described by a partial differential equation (the well-known Navier-Stokes equation). Within that context, the optimisation technique bridges the gap between (linear) optimal perturbation theory and the (nonlinear) dynamical systems approach to fluid flows. The fact that the technique has now been recently shown to work in this very high dimensional setting augurs well for its utility in other physical systems.

Continuum thermodynamics of chemically reacting fluid mixtures

Abstract: We consider viscous, heat conducting mixtures of molecularly miscible chemical species forming a fluid in which the constituents can undergo chemical reactions. Assuming a common temperature for all components, we derive a closed system of partial mass and partial momentum balances plus a mixture balance of internal energy. This is achieved by careful exploitation of the entropy principle and requires appropriate definitions of absolute temperature and chemical potentials, based on an adequate definition of thermal energy excluding diffusive contributions. The resulting interaction forces split into a thermo-mechanical and a chemical part, where the former turns out to be symmetric in case of binary interactions. For chemically reacting systems and as a new result, the chemical interaction force is a contribution being non-symmetric outside of chemical equilibrium. The theory also provides a rigorous derivation of the so-called generalized thermodynamic driving forces, avoiding the use of approximate solutions to the Boltzmann equations. Moreover, using an appropriately extended version of the entropy principle and introducing cross-effects already before closure as entropy invariant couplings between principal dissipative mechanisms, the Onsager symmetry relations become a strict consequence. With a classification of the factors in the binary products of the entropy production according to their parity--instead of the classical partition into so-called fluxes and driving forces--the apparent anti-symmetry of certain couplings is thereby also revealed. If the diffusion velocities are small compared to the speed of sound, the Maxwell-Stefan equations follow in the case without chemistry, thereby neglecting wave phenomena in the diffusive motion. This results in a reduced model with only mass being balanced individually. In the reactive case ...

On the physics of fizzing: How bubble bursting controls droplets ejection

Abstract: Bubbles at a free surface surface usually burst in ejecting myriads of droplets. Focusing on the bubble bursting jet, prelude for these aerosols, we propose a simple scaling for the jet velocity and we unravel experimentally the intricate roles of bubble shape, capillary waves, gravity and liquid properties. We demonstrate that droplets ejection unexpectedly changes with liquid properties. In particular, using damping action of viscosity, self-similar collapse can be sheltered from capillary ripples and continue closer to the singular limit, therefore producing faster and smaller droplets.These results pave the road to the control of the bursting bubble aerosols.

Direct numerical simulation of turbulent channel flow up to $Re_\tau \approx 5200$

Abstract: A direct numerical simulation of incompressible channel flow at $Re_\tau$ = 5186 has been performed, and the flow exhibits a number of the characteristics of high Reynolds number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Karman constant $\kappa = 0.384 \pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in $k$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bi-modal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.

Toward compressed DMD: spectral analysis of fluid flows using sub-Nyquist-rate PIV data

Abstract: Dynamic mode decomposition (DMD) is a powerful and increasingly popular tool for performing spectral analysis of fluid flows. However, it requires data that satisfy the Nyquist-Shannon sampling criterion. In many fluid flow experiments, such data are impossible to capture. We propose a new approach that combines ideas from DMD and compressed sensing. Given a vector-valued signal, we take measurements randomly in time (at a sub-Nyquist rate) and project the data onto a low-dimensional subspace. We then use compressed sensing to identify the dominant frequencies in the signal and their corresponding modes. We demonstrate this method using two examples, analyzing both an artificially constructed test dataset and particle image velocimetry data collected from the flow past a cylinder. In each case, our method correctly identifies the characteristic frequencies and oscillatory modes dominating the signal, proving the proposed method to be a capable tool for spectral analysis using sub-Nyquist-rate sampling.

The Numerical Simulation of Ship Waves using Cartesian Grid Methods

Abstract: Two different cartesian-grid methods are used to simulate the flow around the DDG 5415. The first technique uses a "coupled level-set and volume-of-fluid" (CLS) technique to model the free-surface interface. The no-flux boundary condition on the hull is imposed using a finite-volume technique. The second technique uses a level-set technique (LS) to model the free-surface interface. A body-force technique is used to impose the hull boundary condition. The predictions of both numerical techniques are compared to whisker-probe measurements of the DDG 5415. The level-set technique is also used to investigate the breakup of a two-dimensional spray sheet.

Ship waves in the presence of uniform vorticity

Abstract: Lord Kelvin's result that waves behind a ship lie within a half-angle 19 deg 28' is perhaps the most famous and striking result in the field of surface waves. We solve the linear ship wave problem in the presence of a shear current of constant vorticity S, and show that the Kelvin angles (one each side of wake) as well as other aspects of the wake depend closely on the "shear Froude number" Frs=VS/g (based on length g/S^2 and the ship's speed V), and on the angle between current and the ship's line of motion. In all directions except exactly along the shear flow there exists a critical value of Frs beyond which no transverse waves are produced, and where the full wake angle reaches 180 deg. Such critical behaviour is previously known from waves at finite depth. For side-on shear, one Kelvin angle can exceed 90 deg. On the other hand, the angle of maximum wave amplitude scales as 1/Fr (Fr based on size of ship) when Fr >> 1, a scaling virtually unaffected by the shear flow.

A framework for studying the effect of compliant surfaces on wall turbulence

Abstract: This paper extends the resolvent formulation proposed by McKeon & Sharma (2010) to consider turbulence-compliant wall interactions. Under this formulation, the turbulent velocity field is expressed as a linear superposition of propagating modes, identified via a gain-based decomposition of the Navier-Stokes equations. Compliant surfaces, modeled as a complex wall-admittance linking pressure and velocity, affect the gain and structure of these modes. With minimal computation, this framework accurately predicts the emergence of the quasi-2D propagating waves observed in recent direct numerical simulations. Further, the analysis also enables the rational design of compliant surfaces, with properties optimized to suppress flow structures energetic in wall turbulence. It is shown that walls with unphysical negative damping are required to interact favorably with modes resembling the energetic near-wall cycle, which could explain why previous studies have met with limited success. Positive-damping walls are effective for modes resembling the so-called very large-scale motions (VLSMs), indicating that compliant surfaces may be better suited for application at higher Reynolds number. Unfortunately, walls that suppress structures energetic in natural turbulence are also predicted to have detrimental effects elsewhere in spectral space. Consistent with previous experiments and simulations, slow-moving spanwise-constant structures are particularly susceptible to further amplification. Mitigating these adverse effects will be central to the development of compliant coatings that have a net positive influence on the flow.

Undulatory swimming in shear-thinning fluids: Experiments with C. elegans

Abstract: The swimming behaviour of microorganisms can be strongly influenced by the rheology of their fluid environment. In this manuscript, we experimentally investigate the effects of shear-thinning viscosity on the swimming behaviour of an undulatory swimmer, the nematode Caenorhabditis elegans. Tracking methods are used to measure the swimmer's kinematic data (including propulsion speed) and velocity fields. We find that shear-thinning viscosity modifies the velocity fields produced by the swimming nematode but does not modify the nematode's speed and beating kinematics. Velocimetry data show significant enhancement in local vorticity and circulation and an increase in fluid velocity near the nematode's tail compared to Newtonian fluids of similar effective viscosity. These findings are compared to recent theoretical and numerical results.

Surface structure determines dynamic wetting

Abstract: Liquid wetting of a surface is omnipresent in nature and the advance of micro-fabrication and assembly techniques in recent years offers increasing ability to control this phenomenon. Here, we identify how surface roughness influences the initial dynamic spreading of a partially wetting droplet by studying the spreading on a solid substrate patterned with microstructures just a few micrometers in size. We reveal that the roughness influence can be quantified in terms of a line friction coefficient for the energy dissipation rate at the contact line, and that this can be described in a simple formula in terms of the geometrical parameters of the roughness and the line-friction coefficient of the planar surface. We further identify a criterion to predict if the spreading will be controlled by this surface roughness or by liquid inertia. Our results point to the possibility of selectively controlling the wetting behavior by engineering the surface structure.

The effect of the Basset history force on particle clustering in Homogeneous and Isotropic Turbulence

Abstract: We study the effect of the Basset history force on the dynamics of small particles transported in homogeneous and isotropic turbulence and show that this term, often neglected in previous numerical studies, reduces the small-scale clustering typical of inertial particles. The contribution of this force to the total particle acceleration is, on average, responsible for about 10% of the total acceleration and particularly relevant during rare strong events. At moderate density ratios, i.e. sand or metal powder in water, its presence alters the balance of forces determining the particle acceleration.

Heterogeneous Computing on Mixed Unstructured Grids with PyFR

Abstract: PyFR is an open-source high-order accurate computational fluid dynamics solver for mixed unstructured grids that can target a range of hardware platforms from a single codebase. In this paper we demonstrate the ability of PyFR to perform high-order accurate unsteady simulations of flow on mixed unstructured grids using heterogeneous multi-node hardware. Specifically, after benchmarking single-node performance for various platforms, PyFR v0.2.2 is used to undertake simulations of unsteady flow over a circular cylinder at Reynolds number 3 900 using a mixed unstructured grid of prismatic and tetrahedral elements on a desktop workstation containing an Intel Xeon E5-2697 v2 CPU, an NVIDIA Tesla K40c GPU, and an AMD FirePro W9100 GPU. Both the performance and accuracy of PyFR are assessed. PyFR v0.2.2 is freely available under a 3-Clause New Style BSD license (see this http URL).

On the kurtosis of deep-water gravity waves

Abstract: In this paper, we revisit Janssen's (2003) formulation for the dynamic excess kurtosis of weakly nonlinear gravity waves at deep water. For narrowband directional spectra, the formulation is given by a sixfold integral that depends upon the Benjamin-Feir index and the parameter $R=\sigma_{\theta}^{2}/2 u^{2}$, a measure of short-crestedness for the dominant waves with $ u$ and $\sigma_{\theta}$} denoting spectral bandwidth and angular spreading. Our refinement leads to a new analytical solution for the dynamic kurtosis of narrowband directional waves described with a Gaussian type spectrum. For multidirectional or short-crested seas initially homogenous and Gaussian, in a focusing (defocusing) regime dynamic kurtosis grows initially, attaining a positive maximum (negative minimum) at the intrinsic time scale \[ \tau_{c}= u^{2}\omega_{0}t_{c}=1/\sqrt{3R},\qquad\mathrm{or}\qquad t_{c}/T_{0}\approx0.13/ u\sigma_{\theta}, \] where $\omega_{0}=2\pi/T_{0}$ denotes the dominant angular frequency. Eventually the dynamic excess kurtosis tends monotonically to zero as the wave field reaches a quasi-equilibrium state characterized with nonlinearities mainly due to bound harmonics. Quasi-resonant interactions are dominant only in unidirectional or long-crested seas where the longer-time dynamic kurtosis can be larger than that induced by bound harmonics, especially as the Benjamin-Feir index increases. Finally, we discuss the implication of these results on the prediction of rogue waves.

Local boundary layer scales in turbulent Rayleigh-Benard convection

Abstract: We compute fully local boundary layer scales in three-dimensional turbulent Rayleigh-Benard convection. These scales are directly connected to the highly intermittent fluctuations of the fluxes of momentum and heat at the isothermal top and bottom walls and are statistically distributed around the corresponding mean thickness scales. The local boundary layer scales also reflect the strong spatial inhomogeneities of both boundary layers due to the large-scale, but complex and intermittent, circulation that builds up in closed convection cells. Similar to turbulent boundary layers, we define inner scales based on local shear stress which can be consistently extended to the classical viscous scales in bulk turbulence, e.g. the Kolmogorov scale, and outer scales based on slopes at the wall. We discuss the consequences of our generalization, in particular the scaling of our inner and outer boundary layer thicknesses and the resulting shear Reynolds number with respect to Rayleigh number. The mean outer thickness scale for the temperature field is close to the standard definition of a thermal boundary layer thickness. In the case of the velocity field, under certain conditions the outer scale follows a similar scaling as the Prandtl-Blasius type definition with respect to Rayleigh number, but differs quantitatively. The friction coefficient c_epsilon scaling is found to fall right between the laminar and turbulent limits which indicates that the boundary layer exhibits transitional behavior. Additionally, we conduct an analysis of the recently suggested dissipation layer thickness scales versus Rayleigh number and find a transition in the scaling. We also performed one study of aspect ratio equal to three in the case of Ra=1e+8.