scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Fluid Dynamics in 2009"


Journal ArticleDOI
TL;DR: The appearance of rogue waves is well known in oceanographic, optics, and cold-matter systems as mentioned in this paper, and there is even a possibility for the existence of atmospheric rogue waves.
Abstract: The appearance of rogue waves is well known in oceanographics, optics, and cold matter systems. Here we show a possibility for the existence of atmospheric rogue waves.

205 citations


Journal ArticleDOI
TL;DR: In this paper, the authors present an experimental investigation of intense turbulence generated by a class of low-blockage space-filling fractal square grids and show that the centerline streamwise variation of most of the statistical properties of the turbulent flow can be scaled by a wake interaction length-scale $x_*.
Abstract: We present an experimental investigation of intense turbulence generated by a class of low-blockage space-filling fractal square grids. We confirm the existence of a protacted production region followed by a decaying region, as first reported by Hurst & Vassilicos (Physics of Fluids, 2007). We show that the centerline streamwise variation of most of the statistical properties of the turbulent flow can be scaled by a wake interaction length-scale $x_*$. We also confirm the finding of Seoud and Vassilicos (Physics of Fluids, 2007) that the ratio of the integral length-scale $L_u$ to the Taylor micro-scale $\lambda$ remains constant in the decaying region whereas the Reynolds number $Re_\lambda$ strongly decreases. As a result the scaling $L_{u}/\lambda \sim Re_{\lambda}$ which follows from the $u'^{3}/L_u$ scaling of the dissipation rate in boundary-free shear flows and in usual grid-generated turbulence does not hold here. However, we show that the ratio $L_u/\lambda$ is an increasing function of the inlet Reynolds number $Re_0$. This extraordinary decoupling is consistent with a self-preserving single length-scale decaying homogeneous turbulence proposed by George & Wang (Physics of Fluids, 2009) with which our results are compared.

201 citations


Journal ArticleDOI
TL;DR: In this article, the influence of bending rigidity of a flexible heaving wing on its propulsive performance in a two-dimensional imposed parallel flow is investigated in the inviscid limit.
Abstract: The influence of the bending rigidity of a flexible heaving wing on its propulsive performance in a two-dimensional imposed parallel flow is investigated in the inviscid limit. Potential flow theory is used to describe the flow over the flapping wing. The vortical wake of the wing is accounted for by the shedding of point vortices with unsteady intensity from the wing's trailing edge. The trailing-edge flapping amplitude is shown to be maximal for a discrete set of values of the rigidity, at which a resonance occurs between the forcing frequency and a natural frequency of the system. A quantitative comparison of the position of these resonances with linear stability analysis results is presented. Such resonances induce maximum values of the mean developed thrust and power input. The flapping efficiency is also shown to be greatly enhanced by flexibility.

174 citations


Journal ArticleDOI
TL;DR: In this paper, the authors consider sinusoidal waves of spanwise velocity which vary in time and are modulated in space along the streamwise direction, so that the phase speed may be null, positive or negative, and they can be either stationary or traveling forward or backward in the direction of the mean flow.
Abstract: Waves of spanwise velocity imposed at the walls of a plane turbulent channel flow are studied by Direct Numerical Simulations. We consider sinusoidal waves of spanwise velocity which vary in time and are modulated in space along the streamwise direction. The phase speed may be null, positive or negative, so that the waves may be either stationary or traveling forward or backward in the direction of the mean flow. Such a forcing includes as particular cases two known techniques for reducing friction drag: the oscillating wall technique (a traveling wave with infinite phase speed) and the recently proposed steady distribution of spanwise velocity (a wave with zero phase speed). The traveling waves alter the friction drag significantly. Waves which slowly travel forward produce a large reduction of drag, that can relaminarize the flow at low values of the Reynolds number. Faster waves yield a totally different outcome, i.e. drag increase. Even faster waves produce a drag reduction effect again. Backward-traveling waves instead lead to drag reduction at any speed. The traveling waves, when they reduce drag, operate in similar fashion to the oscillating wall, with an improved energetic efficiency. Drag increase is observed when the waves travel at a speed comparable with that of the convecting near-wall turbulence structures. A diagram illustrating the different flow behaviors is presented.

142 citations


Journal ArticleDOI
TL;DR: Gekle et al. as mentioned in this paper combine detailed boundary-integral simulations with analytical modeling to describe the formation and break-up of Worthington jets in two common physical systems: the impact of a circular disc on a liquid surface and the release of air bubbles from an underwater nozzle.
Abstract: Helped by the careful analysis of their experimental data, Worthington (1897) described roughly the mechanism underlying the formation of high-speed jets ejected after the impact of an axisymmetric solid on a liquid-air interface. In this work we combine detailed boundary-integral simulations with analytical modeling to describe the formation and break-up of such Worthington jets in two common physical systems: the impact of a circular disc on a liquid surface and the release of air bubbles from an underwater nozzle. We first show that the jet base dynamics can be predicted for both systems using our earlier model in Gekle, Gordillo, van der Meer and Lohse. Phys. Rev. Lett. 102 (2009). Nevertheless, our main point here is to present a model which allows us to accurately predict the shape of the entire jet. Good agreement with numerics and some experimental data is found. Moreover, we find that, contrarily to the capillary breakup of liquid cylinders in vacuum studied by Rayleigh, the breakup of stretched liquid jets at high values of both Weber and Reynolds numbers is not triggered by the growth of perturbations coming from an external source of noise. Instead, the jet breaks up due to the capillary deceleration of the liquid at the tip which produces a corrugation to the jet shape. This perturbation, which is self-induced by the flow, will grow in time promoted by a capillary mechanism. We are able to predict the exact shape evolution of Worthington jets ejected after the impact of a solid object - including the size of small droplets ejected from the tip due to a surface-tension driven instability - using as the single input parameters the minimum radius of the cavity and the flow field before the jet emerges.

119 citations


Posted Content
TL;DR: It is shown that the fluid propelled increases drastically compared to synchronously beating cilia, and is unidirectional, even when the phase difference between adjacent cilia is small, which provides guidelines for the optimal design of magnetically-driven artificial cilia for microfluidic propulsion.
Abstract: Natural cilia are hair-like microtubule-based structures that are able to move fluid at low Reynolds number through asymmetric motion. In this paper we follow a biomimetic approach to design artificial cilia lining the inner surface of microfluidic channels with the goal to propel fluid. The artificial cilia consist of polymer films filled with magnetic nanoparticles. The asymmetric, non-reciprocating motion is generated by tuning an external magnetic field. To obtain the magnetic field and associated magnetization local to the cilia we solve the Maxwell equations, from which the magnetic torques can be deduced. To obtain the ciliary motion we solve the dynamic equations of motion which are then fully coupled to the fluid dynamic equations that describe fluid flow around the cilia. By doing so we show that by properly tuning the applied magnetic field, asymmetric ciliary motion can be generated that is able to propel fluid in a microchannel. The results are presented in terms of three dimensionless parameters that fully delineate the asymmetry and cycle time as a function of the relative contribution of elastic, inertial, magnetic and viscous fluid forces.

119 citations


Journal ArticleDOI
TL;DR: In this paper, the authors derived analytical formulas for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that can be used for any surface slip fraction (validated by numerical calculations).
Abstract: Super-hydrophobic array of grooves containing trapped gas (stripes), have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused on idealized cases of stick-perfect slip stripes, with limited guidance. Here, we analyze the experimentally relevant situation of a pressure-driven flow past striped slip-stick surfaces with arbitrary local slip at the gas sectors. We derive analytical formulas for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that can be used for any surface slip fraction (validated by numerical calculations). By representing eigenvalues of the slip length-tensor, they allow us to obtain the effective slip for any orientation of stripes with respect to the mean flow. Our results imply that flow past stripes is controlled by the ratio of the local slip length to texture size. In case of a large (compared to the texture period) slip at the gas areas, surface anisotropy leads to a tensorial effective slip, by attaining the values predicted earlier for a perfect local slip. Both effective slip lengths and anisotropy of the flow decrease when local slip becomes of the order of texture period. In the case of small slip, we predict simple surface-averaged, isotropic flows (independent of orientation). These results provide a framework for the rational design of super-hydrophobic surfaces and devices.

118 citations


Journal ArticleDOI
TL;DR: In this article, 3D-Particle Tracking (3D-PTV) and Phase Sensitive Constant Temperature Anemometry in pseudo-turbulence were performed to investigate bubble clustering and to obtain the mean bubble rise velocity, distributions of bubble velocities, and energy spectra at dilute gas concentrations.
Abstract: 3D-Particle Tracking (3D-PTV) and Phase Sensitive Constant Temperature Anemometry in pseudo-turbulence--i.e., flow solely driven by rising bubbles-- were performed to investigate bubble clustering and to obtain the mean bubble rise velocity, distributions of bubble velocities, and energy spectra at dilute gas concentrations ($\alpha \leq2.2$%). To characterize the clustering the pair correlation function $G(r,\theta)$ was calculated. The deformable bubbles with equivalent bubble diameter $d_b=4-5$ mm were found to cluster within a radial distance of a few bubble radii with a preferred vertical orientation. This vertical alignment was present at both small and large scales. For small distances also some horizontal clustering was found. The large number of data-points and the non intrusiveness of PTV allowed to obtain well-converged Probability Density Functions (PDFs) of the bubble velocity. The PDFs had a non-Gaussian form for all velocity components and intermittency effects could be observed. The energy spectrum of the liquid velocity fluctuations decayed with a power law of -3.2, different from the $\approx -5/3$ found for homogeneous isotropic turbulence, but close to the prediction -3 by \cite{lance} for pseudo-turbulence.

98 citations


Journal ArticleDOI
TL;DR: In this paper, the authors present results from two 1536^3 direct numerical simulations of rotating turbulence where both energy and helicity are injected into the flow by an external forcing.
Abstract: We present results from two 1536^3 direct numerical simulations of rotating turbulence where both energy and helicity are injected into the flow by an external forcing. The dual cascade of energy and helicity towards smaller scales observed in isotropic and homogeneous turbulence is broken in the presence of rotation, with the development of an inverse cascade of energy now coexisting with direct cascades of energy and helicity. In the direct cascade range, the flux of helicity dominates over that of energy at low Rossby number. These cascades have several consequences for the statistics of the flow. The evolution of global quantities and of the energy and helicity spectra is studied, and comparisons with simulations at different Reynolds and Rossby numbers at lower resolution are done to identify scaling laws.

96 citations


Journal ArticleDOI
TL;DR: In this paper, a steady wall-based forcing is proposed to passively interact with the turbulent flow to achieve drag reduction through a suitable distribution of roughness, designed to excite a selected streamwise wavelength.
Abstract: Steady forcing at the wall of a channel flow is studied via DNS to assess its ability of yielding reductions of turbulent friction drag. The wall forcing consists of a stationary distribution of spanwise velocity that alternates in the streamwise direction. The idea behind the forcing builds upon the existing technique of the spanwise wall oscillation, and exploits the convective nature of the flow to achieve an unsteady interaction with turbulence. The analysis takes advantage of the equivalent laminar flow, that is solved analytically to show that the energetic cost of the forcing is unaffected by turbulence. In a turbulent flow, the alternate forcing is found to behave similarly to the oscillating wall; in particular an optimal wavelength is found that yields a maximal reduction of turbulent drag. The energetic performance is significantly improved, with more than 50% of maximum friction saving at large intensities of the forcing, and a net energetic saving of 23% for smaller intensities. Such a steady, wall-based forcing may pave the way to passively interacting with the turbulent flow to achieve drag reduction through a suitable distribution of roughness, designed to excite a selected streamwise wavelength.

96 citations


Journal ArticleDOI
TL;DR: In this paper, Pasquetti et al. adapted the pseudo-penalisation spectral method to integrate numerically the simultaneous dynamics of the particle and of the fluid, which gave a unique handle on the limit of validity of point-particle approximations, which are generally used in applicative situations.
Abstract: The dynamics of neutrally buoyant particles transported by a turbulent flow is investigated for spherical particles with radii of the order of the Kolmogorov dissipative scale or larger. The pseudo-penalisation spectral method that has been proposed by Pasquetti et al. (2008) is adapted to integrate numerically the simultaneous dynamics of the particle and of the fluid. Such a method gives a unique handle on the limit of validity of point-particle approximations, which are generally used in applicative situations. Analytical predictions based on such models are compared to result of very well resolved direct numerical simulations. Evidence is obtained that Faxen corrections give dominant finite-size corrections to velocity and acceleration fluctuations for particle diameters up to four times the Kolmogorov scale. The dynamics of particles with larger diameters is dominated by inertial-range physics, and is consistent with predictions obtained from dimensional analysis.

Journal ArticleDOI
TL;DR: In this paper, the authors examined the relationship between the buoyancy and mixing in a wind-and buoyancy-driven thermally stratified ocean to determine whether these constrain the mechanical sources of stirring.
Abstract: A key issue in stratified turbulence theory concerns the nature of the link between D(APE), the dissipation rate of available potential energy APE, and W_{r,turbulent}, the turbulent rate of change of background gravitational potential energy GPE_r, which are both controlled by molecular diffusion. For Boussinesq fluids with a linear equation of state, this link is simply W_{r,turbulent}=D(APE), widely interpreted as implying that GPE_r increases at the expense of APE, in contrast with the laminar case where GPE_r increases at the expense of internal energy (IE). This idea is revisited here by regarding IE as the sum of three distinct subcomponents: available internal energy (AIE), exergy (IE_{exergy}), and dead internal energy (IE_0). In this new view, D(APE) is the dissipation rate of APE into IE_0, while both W_{r,laminar} and W_{r,turbulent} convert IE_{exergy} into GPE_r. The equality W_{r,turbulent}=D(APE) thus states that IE_{exergy} is converted into GPE_r at the same rate as APE is dissipated into IE_0. For non-Boussinesq fluids, the equality D(APE)=W_{r,turbulent} is at best a good approximation, for W_{r,turbulent} is generally smaller than D(APE), and sometimes even negative for a strongly nonlinear equation of state. In a second step, the link between stirring and mixing is examined for a wind-and buoyancy-driven thermally stratified ocean to determine whether these constrain the mechanical sources of stirring, as recently advocated. It is established that the coupling between stirring and mixing cannot refute the traditional buoyancy-driven view of the so-called meridional overturning circulation, in contrast to recent claims. In fact, the buoyancy forcing appears to be as important as the mechanical forcing in stirring and driving the large-scale ocean circulation.

Journal ArticleDOI
TL;DR: In this paper, it was shown that high values of the magnetic Prandtl number have a strong influence on convection-driven dynamos in rotating spherical shells filled with electrically conducting fluids, while low values promote dynamo action through the shear provided by differential rotation, while the generation of magnetic fields is more difficult to sustain in high-Prandtl-number fluids.
Abstract: The value of the Prandtl number $P$ exerts a strong influence on convection-driven dynamos in rotating spherical shells filled with electrically conducting fluids. Low Prandtl numbers promote dynamo action through the shear provided by differential rotation, while the generation of magnetic fields is more difficult to sustain in high-Prandtl-number fluids where higher values of the magnetic Prandtl number $P_m$ are required. The magnetostrophic approximation often used in dynamo theory appears to be valid only for relatively high values of $P$ and $P_m$. Dynamos with a minimum value of $P_m$ seem to be most readily realizable in the presence of convection columns at moderately low values of $P$. The structure of the magnetic field varies strongly with $P$ in that dynamos with a strong axial dipole field are found for high values of $P$ while the energy of this component is exceeded by that of the axisymmetric toroidal field and by that of the non-axisymmetric components at low values of $P$. Some conclusions are discussed in relation to the problem of the generation of planetary magnetic fields by motions in their electrically conducting liquid cores.

Journal ArticleDOI
TL;DR: In this article, the relative dispersion of pairs of inertial particles in incompressible, homogeneous, and isotropic turbulence is studied by means of direct numerical simulations at two values of the Taylor-scale Reynolds number $Re_{\lambda} \sim 200$ and 400.
Abstract: The relative dispersion of pairs of inertial particles in incompressible, homogeneous, and isotropic turbulence is studied by means of direct numerical simulations at two values of the Taylor-scale Reynolds number $Re_{\lambda} \sim 200$ and 400. The evolution of both heavy and light particle pairs is analysed at varying the particle Stokes number and the fluid-to-particle density ratio. For heavy particles, it is found that turbulent dispersion is schematically governed by two temporal regimes. The first is dominated by the presence, at large Stokes numbers, of small-scale caustics in the particle velocity statistics, and it lasts until heavy particle velocities have relaxed towards the underlying flow velocities. At such large scales, a second regime starts where heavy particles separate as tracers particles would do. As a consequence, at increasing inertia, a larger transient stage is observed, and the Richardson diffusion of simple tracers is recovered only at large times and large scales. These features also arise from a statistical closure of the equation of motion for heavy particle separation that is proposed, and which is supported by the numerical results. In the case of light particles with high density ratios, strong small-scale clustering leads to a considerable fraction of pairs that do not separate at all, although the mean separation increases with time. This effect strongly alters the shape of the probability density function of light particle separations.

Journal ArticleDOI
TL;DR: Several sets of measurements of the hydrodynamic resistance of droplets, expressed in terms of a "resistance length" are presented to reveal its dependence on a number of parameters, such as the viscosity, the volume of droplet, their velocity as well as the spacing between them.
Abstract: The overall traffic of droplets in a network of microfluidic channels is strongly influenced by the liquid properties of the moving droplets. In particular, the effective hydrodynamic resistance of individual droplets plays a key role in their global behavior. We here propose two simple and low-cost experimental methods for measuring this parameter by analyzing the dynamics of a regular sequence of droplets injected into an "asymmetric loop" network. The choice of a droplet taking either route through the loop is influenced by the presence of previous droplets which modulate the hydrodynamic resistance of the branches they are sitting in. We propose to extract the effective resistance of a droplet from easily observable time series, namely from the choices the droplets make at junctions and from the inter-droplet distances. This becomes possible when utilizing a recently proposed theoretical model, based on a number of simplifying assumptions. We here present several sets of measurements of the hydrodynamic resistance of droplets, expressed in terms of a "resistance length". The aim is twofold, (1) to reveal its dependence on a number of parameters, such as the viscosity, the volume of droplets, their velocity as well as the spacing between them. At the same time (2), by using a standard measurement technique, we compare the limitations of the proposed methods. As an important result of this comparison we obtain the range of validity of the simplifying assumptions made in the theoretical model.

Journal ArticleDOI
TL;DR: In this paper, critical balance -a scale-by-scale balance between the linear propagation and nonlinear interaction time scales -can be used as a universal scaling conjecture for determining the spectra of strong turbulence in anisotropic wave systems.
Abstract: It is proposed that critical balance - a scale-by-scale balance between the linear propagation and nonlinear interaction time scales - can be used as a universal scaling conjecture for determining the spectra of strong turbulence in anisotropic wave systems. Magnetohydrodynamic (MHD), rotating and stratified turbulence are considered under this assumption and, in particular, a novel and experimentally testable energy cascade scenario and a set of scalings of the spectra are proposed for low-Rossby-number rotating turbulence. It is argued that in neutral fluids, the critically balanced anisotropic cascade provides a natural path from strong anisotropy at large scales to isotropic Kolmogorov turbulence at very small scales. It is also argued that the kperp^{-2} spectra seen in recent numerical simulations of low-Rossby-number rotating turbulence may be analogous to the kperp^{-3/2} spectra of the numerical MHD turbulence in the sense that they could be explained by assuming that fluctuations are polarised (aligned) approximately as inertial waves (Alfven waves for MHD).

Journal ArticleDOI
TL;DR: The study of flow of non-Newtonian fluids in porous media is very important and serves a wide variety of practical applications in processes such as enhanced oil recovery from underground reservoirs, filtration of polymer solutions and soil remediation through the removal of liquid pollutants.
Abstract: The study of flow of non-Newtonian fluids in porous media is very important and serves a wide variety of practical applications in processes such as enhanced oil recovery from underground reservoirs, filtration of polymer solutions and soil remediation through the removal of liquid pollutants. These fluids occur in diverse natural and synthetic forms and can be regarded as the rule rather than the exception. They show very complex strain and time dependent behavior and may have initial yield-stress. Their common feature is that they do not obey the simple Newtonian relation of proportionality between stress and rate of deformation. Non-Newtonian fluids are generally classified into three main categories: time-independent whose strain rate solely depends on the instantaneous stress, time-dependent whose strain rate is a function of both magnitude and duration of the applied stress and viscoelastic which shows partial elastic recovery on removal of the deforming stress and usually demonstrates both time and strain dependency. In this article the key aspects of these fluids are reviewed with particular emphasis on single-phase flow through porous media. The four main approaches for describing the flow in porous media are examined and assessed. These are: continuum models, bundle of tubes models, numerical methods and pore-scale network modeling.

Journal ArticleDOI
TL;DR: In this paper, the turbulence transition of plane Couette flow in large domains where localised perturbations are observed to generate growing turbulent spots was studied and invariant structures intermediate between laminar and turbulent flow were determined.
Abstract: We study the turbulence transition of plane Couette flow in large domains where localised perturbations are observed to generate growing turbulent spots. Extending previous studies on the boundary between laminar and turbulent dynamics we determine invariant structures intermediate between laminar and turbulent flow. In wide but short domains we find states that are localised in spanwise direction, and in wide and long domains the states are also localised in downstream direction. These localised states act as critical nuclei for the transition to turbulence in spatially extended domains.

Posted Content
TL;DR: In this article, a more systematic approach based on studying discrete equilibrium structure factors for a broad class of explicit linear finite-volume schemes is introduced, which provides a better characterization of the accuracy of a spatio-temporal discretization as a function of wavenumber and frequency, allowing us to distinguish between behavior at long wavelengths and short wavelengths, where stability concerns are of greater importance.
Abstract: This paper describes the development and analysis of finite-volume methods for the Landau-Lifshitz Navier-Stokes (LLNS) equations and related stochastic partial differential equations in fluid dynamics. The LLNS equations incorporate thermal fluctuations into macroscopic hydrodynamics by the addition of white-noise fluxes whose magnitudes are set by a fluctuation-dissipation relation. Originally derived for equilibrium fluctuations, the LLNS equations have also been shown to be accurate for non-equilibrium systems. Previous studies of numerical methods for the LLNS equations focused primarily on measuring variances and correlations computed at equilibrium and for selected non-equilibrium flows. In this paper, we introduce a more systematic approach based on studying discrete equilibrium structure factors for a broad class of explicit linear finite-volume schemes. This new approach provides a better characterization of the accuracy of a spatio-temporal discretization as a function of wavenumber and frequency, allowing us to distinguish between behavior at long wavelengths, where accuracy is a prime concern, and short wavelengths, where stability concerns are of greater importance. We use this analysis to develop a specialized third-order Runge Kutta scheme that minimizes the temporal integration error in the discrete structure factor at long wavelengths for the one-dimensional linearized LLNS equations. Together with a novel method for discretizing the stochastic stress tensor in dimension larger than one, our improved temporal integrator yields a scheme for the three-dimensional equations that satisfies a discrete fluctuation-dissipation balance for small time steps and is also sufficiently accurate even for time steps close to the stability limit.

Journal ArticleDOI
TL;DR: In this article, the transition from dynamos dominated by non-axisymmetric components of the magnetic field to those dominated by the axishemetric components depends on the magnetic Prandtl number as well as on the ordinary PrandTL number for higher values of the rotation parameter.
Abstract: For the understanding of planetary and stellar dynamos an overview of the major parameter dependences of convection driven dynamos in rotating spherical fluid shells is desirable. Although the computationally accessible parameter space is limited, earlier work is extended with emphasis on higher Prandtl numbers and uniform heat flux condition at the outer boundary. The transition from dynamos dominated by non-axisymmetric components of the magnetic field to those dominated by the axisymmetric components depends on the magnetic Prandtl number as well as on the ordinary Prandtl number for higher values of the rotation parameter $\tau$. The dependence of the transition on the latter parameter is also discussed. A variety of oscillating dynamos is presented and interpreted in terms of dynamo waves, standing oscillation or modified relaxation oscillations.

Journal ArticleDOI
TL;DR: In this paper, the authors studied the influence of a generic slip boundary condition on wake dynamics and showed that slip significantly delays the onset of recirculation and shedding in the wake behind a circular cylinder.
Abstract: We study numerically the two-dimensional flow past a circular cylinder as a prototypical transitional flow, and investigate the influence of a generic slip boundary condition on the wake dynamics. We show that slip significantly delays the onset of recirculation and shedding in the wake behind the cylinder. As expected, the drag on the cylinder decreases with slip, with an increased drag sensitivity for large Reynolds numbers. We also show that past the critical shedding Reynolds number, slip decreases the vorticity intensity in the wake, as well as the lift forces on the cylinder, but increases the shedding frequency. We further provide evidence that the shedding transition can be interpreted as a critical accumulation of surface vorticity, similarly to related studies on wake instability of axisymmetric bodies. Finally, we propose that our results could be used as a passive method to infer the effective friction properties of slipping surfaces.

Journal ArticleDOI
TL;DR: In this paper, a modified Bautista-Manero model was used for modeling the flow of viscoelastic materials which also show thixotropic attributes, and an algorithm, originally proposed by Philippe Tardy, that employs this model to simulate steadystate time-dependent flow was implemented in a non-Newtonian flow simulation code using pore-scale modeling.
Abstract: In this article, the extensional flow and viscosity and the converging-diverging geometry were examined as the basis of the peculiar viscoelastic behavior in porous media. The modified Bautista-Manero model, which successfully describes shearthinning, elasticity and thixotropic time-dependency, was used for modeling the flow of viscoelastic materials which also show thixotropic attributes. An algorithm, originally proposed by Philippe Tardy, that employs this model to simulate steadystate time-dependent flow was implemented in a non-Newtonian flow simulation code using pore-scale modeling and the initial results were analyzed. The findings are encouraging for further future development.

Posted Content
TL;DR: In this paper, the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows was theoretically analyzed and numerically evaluated.
Abstract: In this work, we have theoretically analyzed and numerically evaluated the accuracy of high-order lattice Boltzmann (LB) models for capturing non-equilibrium effects in rarefied gas flows. In the incompressible limit, the LB equation is proved to be equivalent to the linearized Bhatnagar-Gross-Krook (BGK) equation. Therefore, when the same Gauss-Hermite quadrature is used, LB method closely assembles the discrete velocity method (DVM). In addition, the order of Hermite expansion for the equilibrium distribution function is found not to be correlated with the approximation order in terms of the Knudsen number to the BGK equation, which was previously suggested by \cite{2006JFM...550..413S}. Furthermore, we have numerically evaluated the LB models for a standing-shear-wave problem, which is designed specifically for assessing model accuracy by excluding the influence of gas molecule/surface interactions at wall boundaries. The numerical simulation results confirm that the high-order terms in the discrete equilibrium distribution function play a negligible role. Meanwhile, appropriate Gauss-Hermite quadrature has the most significant effect on whether LB models can describe the essential flow physics of rarefied gas accurately. For the same order of the Gauss-Hermite quadrature, the exact abscissae will also modestly influence numerical accuracy. Using the same Gauss-Hermite quadrature, the numerical results of both LB and DVM methods are in excellent agreement for flows across a broad range of the Knudsen numbers, which confirms that the LB simulation is similar to the DVM process. Therefore, LB method can offer flexible models suitable for simulating continuum flows at Navier Stokes level and rarefied gas flows at the linearized Boltzmann equation level.

Journal ArticleDOI
TL;DR: In this article, the effect of boundary slip by using lattice Boltzmann simulations is reviewed and applied to investigate fluid flow in hydrophobic and rough microchannels as well as over surfaces covered by nano- or microscale gas bubbles.
Abstract: In this contribution we review recent efforts on investigations of the effect of (apparent) boundary slip by utilizing lattice Boltzmann simulations. We demonstrate the applicability of the method to treat fundamental questions in microfluidics by investigating fluid flow in hydrophobic and rough microchannels as well as over surfaces covered by nano- or microscale gas bubbles.

Journal ArticleDOI
TL;DR: In this paper, a sloshing mode has been found in experiments on the large-scale circulation of turbulent Rayleigh-Benard convection in a cylindrical sample of aspect ratio one.
Abstract: In agreement with a recent experimental discovery by Xia et. al. (2009), we also find a sloshing mode in experiments on the large-scale circulation (LSC) of turbulent Rayleigh-Benard convection in a cylindrical sample of aspect ratio one. The sloshing mode has the same frequency as the torsional oscillation discovered by Funfschilling and Ahlers (2004). We show that both modes can be described by an extension of a model developed previously [Brown and Ahlers (2008)] which consists of permitting a lateral displacement of the LSC circulation plane away from the vertical center line of the sample as well as a variation in displacements with height (such displacements had been excluded in the original model). Pressure gradients produced by the side wall of the container on average center the plane of the LSC so that it prefers to reach its longest diameter. If the LSC is displaced away from this diameter, the walls provide a restoring force. Turbulent fluctuations drive the LSC away from the central alignment, and combined with the restoring force they lead to oscillations. These oscillations are advected along with the LSC. This model predicts the correct wavenumber and phase of the oscillations, as well as estimates of the frequency, amplitude, and probability distributions of the displacements.

Journal ArticleDOI
TL;DR: In this paper, it was shown that the non-relativistic limit of the conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry.
Abstract: The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schroedinger group, which also involves, in addition, Schroedinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally-invariant relativistic theory, the recently discussed Conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.

Journal ArticleDOI
TL;DR: In this article, the scale invariance and non-intermittent properties of the cascade of energy to small scales of rotating turbulence have been investigated in the context of helical rotating turbulence.
Abstract: directnumerical simulations of helical rotating turbulence. Symmetric and anti-symmetric velocity incre-ments are examined, as well as probability density functions of the velocity field and of the helicitydensity. It is found that the direct cascade of energy to small scales is scale invariant and non-intermittent, whereas the direct cascade of helicity is highly intermittent. Furthermore, the studyof structure functions of different orders allows us to identify a recovery of isotropy of strong eventsat very small scales in the flow. Finally, we observe the juxtaposition in space of strong laminar andpersistent helical columns next to time-varying vortex tangles, the former being associated with theself-similarity of energy and the latter with the intermittency of helicity.I. INTRODUCTION

Journal ArticleDOI
TL;DR: In this article, an analytical method to calculate the effective slip length of a mesh-like super-hydrophobic surface is presented. But the authors focus on the case of isolated holes in an otherwise continuous no-slip surface.
Abstract: When a liquid droplet is located above a super-hydrophobic surface, it only barely touches the solid portion of the surface, and therefore slides very easily on it. More generally, super-hydrophobic surfaces have been shown to lead to significant reduction of viscous friction in the laminar regime, so it is of interest to quantify their effective slipping properties as a function of their geometric characteristics. Most previous studies have considered flows bounded by arrays of either long grooves, or isolated solid pillars on an otherwise flat solid substrate, and for which therefore the surrounding air constitutes the continuous phase. Here we consider instead the case where the super-hydrophobic surface is made of isolated holes in an otherwise continuous no-slip surface, and specifically focus on the mesh-like geometry recently achieved experimentally. We present an analytical method to calculate the friction of such a surface in the case where the mesh is thin. The results for the effective slip length of the surface are computed, compared to simple estimates, and a practical fit is proposed displaying a logarithmic dependance on the area fraction of the solid surface.

Journal ArticleDOI
TL;DR: Whereas the pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle height.
Abstract: We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods - Computational Fluid Dynamics (CFD), Lattice-Boltzmann Equations (LBE) and Molecular Dynamics (MD) The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles theta > 45 While for theta = 30, theta = 40 (no flow) and theta = 60 (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 the numerical experiments provide different results Whilst pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle heigth

Journal ArticleDOI
TL;DR: In this article, the dynamics of an oscillated sessile drop of incompressible liquid with the focus on the contact line hysteresis is under theoretical consideration, where the solid substrate is subject to transverse oscillations, which are assumed small amplitude and high frequency.
Abstract: The dynamics of an oscillated sessile drop of incompressible liquid with the focus on the contact line hysteresis is under theoretical consideration The solid substrate is subject to transverse oscillations, which are assumed small amplitude and high frequency The dynamic boundary condition that involves an ambiguous dependence of the contact angle on the contact line velocity is applied: the contact line starts to move only when the deviation of the contact angle exceeds a certain critical value As a result, the stick-slip dynamics can be observed The frequency response of surface oscillations on the substrate and at the pole of the drop are analyzed It is shown that novel features such as the emergence of antiresonant frequency bands and nontrivial competition of different resonances are caused by contact line hysteresis