Topic overview

physics.flu-dyn

4653 works10692 researchers

Map preview

Start with the graph, then narrow the list

4653works
10692researchers

Next steps

Use the topic as a working map

Open the full map for clusters, then return here to scan ranked papers and people.

Topic graph

See the topic as a live network

Open full explorer

Inspect nearby papers, researchers, institutions and communities without opening a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Papers in this area

24 paper(s) to start with

preprint2017arXiv

Super- and sub-rotating equatorial jets: Newtonian cooling versus Rayleigh friction

Numerical simulations of the shallow water equations on rotating spheres produce mixtures of robust vortices and alternating zonal jets, as seen in the atmospheres of the gas giant planets. However, simulations that include Rayleigh friction invariably produce a sub-rotating (retrograde) equatorial jet for Jovian parameter regimes, whilst observations of Jupiter show a super-rotating (prograde) equatorial jet that has persisted over several decades. Super-rotating equatorial jets have recently been obtained in shallow water simulations that include a Newtonian relaxation of perturbations to the layer thickness to model radiative cooling to space, and in simulations of the thermal shallow water equations that include a similar relaxation term in their temperature equation. Simulations of global quasigeostrophic forms of these different models produce equatorial jets in the same directions as the parent models, suggesting that the mechanism responsible for setting the direction lies within quasigeostrophic theory. We provide such a mechanism by calculating the effective force acting on the thickness-weighted zonal mean flow due to the decay of an equatorially trapped Rossby wave. Dec

preprint2016arXiv

A Hydrodynamic Approach to the Study of Anisotropic Instabilities in Dissipative Relativistic Plasmas

We develop a purely hydrodynamic formalism to describe collisional, anisotropic instabilities in a relativistic plasma, that are usually described with kinetic theory tools. Our main motivation is the fact that coarse-grained models of high particle number systems give more clear and comprehensive physical descriptions of those systems than purely kinetic approaches, and can be more easily tested experimentally as well as numerically. In particular, we aim at developing a theory that describes both a background non-equilibrium fluid configurations and its perturbations, to be able to account for the backreaction of the latter on the former. Our system of equations includes the usual conservation laws for the energy-momentum tensor and for the electric current, and the equations for two new tensors that encode the information about dissipation. To make contact with kinetic theory, we write the different tensors as the moments of a non-equilibrium one-particle distribution function (1pdf) which, for illustrative purposes, we take in the form of a Grad-like ansatz. Although this choice limits the applicability of the formalism to states not far from equilibrium, it retains the main fe

preprint2016arXiv

Effect of fluid inertia on the motion of a collinear swimmer

The swimming of a two-sphere system and of a three-sphere chain in an incompressible viscous fluid is studied on the basis of simplified equations of motion which take account of both Stokes friction and added mass effects. The analysis is based on an explicit expression for the asymptotic periodic swimming velocity and a corresponding evaluation of the mean rate of dissipation. The mean swimming velocity of the two-sphere system is found to be non-vanishing provided that the two spheres are not identical. The swimming of a comparable chain of three identical spheres is much more efficient.

preprint2016arXiv

Finite-time Thin Film Rupture Driven by Generalized Evaporative Loss

Rupture is a nonlinear instability resulting in a finite-time singularity as a fluid layer approaches zero thickness at a point. We study the dynamics of rupture in a generalized mathematical model of thin films of viscous fluids with evaporative effects. The governing lubrication model is a fourth-order nonlinear parabolic partial differential equation with a non-conservative loss term due to evaporation. Several different types of finite-time singularities are observed due to balances between evaporation and surface tension or intermolecular forces. Non-self-similar behavior and two classes of self-similar rupture solutions are analyzed and validated against high resolution PDE simulations.

preprint2016arXiv

Negative pressure in shear thickening band of a dilatant fluid

We perform experiments and numerical simulations to investigate spatial distribution of pressure in a sheared dilatant fluid of the Taylor-Couette flow under a constant external shear stress. In a certain range of shear stress, the flow undergoes the shear thickening oscillation around 20 Hz. We find that, during the oscillation, a localized thickened band rotates around the axis with the flow. {Based upon experiments and numerical simulations, we show that a major part of the thickened band is under negative pressure even in the case of discontinuous shear thickening}, which indicates that the thickening is caused by Reynolds dilatancy; the dilatancy causes the negative pressure in interstitial fluid, which generates contact structure in the granular medium., then frictional resistance hinders rearrangement of the structure and solidifies the medium.

preprint2017arXiv

Faraday waves by resonant triad interaction of surface--compression waves

The propagation of wave disturbances over a vertically oscillating liquid may form standing waves, known as Faraday waves. Here we present an alternative description of the generation and evolution of Faraday waves by nonlinear resonant triad interactions, accounting for gravity surface effects and the slight compressibility of the liquid. To this end, we analyse a triad comprising an infinitely long-crested compression mode and two oppositely propagating subharmonic surface waves.

preprint2015arXiv

Characteristics of undulatory locomotion in granular media

Undulatory locomotion is ubiquitous in nature and observed in different media, from the swimming of flagellated microorganisms in biological fluids, to the slithering of snakes on land, or the locomotion of sandfish lizards in sand. Despite the similarity in the undulating pattern, the swimming characteristics depend on the rheological properties of different media. Analysis of locomotion in granular materials is relatively less developed compared with fluids partially due to a lack of validated force models but recently a resistive force theory in granular media has been proposed and shown useful in studying the locomotion of a sand-swimming lizard. Here we employ the proposed model to investigate the swimming characteristics of a slender filament, of both finite and infinite length, undulating in a granular medium and compare the results with swimming in viscous fluids. In particular, we characterize the effects of drifting and pitching in terms of propulsion speed and efficiency for a finite sinusoidal swimmer. We also find that, similar to Lighthill's results using resistive force theory in viscous fluids, the sawtooth swimmer is the optimal waveform for propulsion speed at

preprint2016arXiv

Turbulence as a problem in non-equilibrium statistical mechanics

The transitional and well-developed regimes of turbulent shear flows exhibit a variety of remarkable scaling laws that are only now beginning to be systematically studied and understood. In the first part of this article, we summarize recent progress in understanding the friction factor of turbulent flows in rough pipes and quasi-two-dimensional soap films, showing how the data obey a two-parameter scaling law known as roughness-induced criticality, and exhibit power-law scaling of friction factor with Reynolds number that depends on the precise form of the nature of the turbulent cascade. These results hint at a non-equilibrium fluctuation-dissipation relation that applies to turbulent flows. The second part of this article concerns the lifetime statistics in smooth pipes around the transition, showing how the remarkable super-exponential scaling with Reynolds number reflects deep connections between large deviation theory, extreme value statistics, directed percolation and the onset of coexistence in predator-prey ecosystems. Both these phenomena reflect the way in which turbulence can be fruitfully approached as a problem in non-equilibrium statistical mechanics.

preprint2017arXiv

Controlling droplet spreading with topography

We present a novel experimental system that can be used to study the dynamics of picoliter droplet spreading over substrates with topographic variations. We concentrate on spreading of a droplet within a recessed stadium-shaped pixel, with applications to the manufacture of POLED displays, and find that the sloping side wall of the pixel can either locally enhance or hinder spreading depending on whether the topography gradient ahead of the contact line is positive or negative. Locally enhanced spreading occurs via the formation of thin pointed rivulets along the side walls of the pixel through a mechanism similar to capillary rise in sharp corners. We demonstrate that a thin-film model combined with an experimentally measured spreading law, which relates the speed of the contact line to the contact angle, provides excellent predictions of the evolving liquid morphologies. We also show that the spreading can be adequately described by a Cox-Voinov law for the majority of the evolution. The model does not include viscous effects and hence, the timescales for the propagation of the thin pointed rivulets are not captured. Nonetheless, this simple model can be used very effectively to

preprint2016arXiv

Random patterns in fish schooling enhance alertness: a hydrodynamic perspective

One of the most highly debated questions in the field of animal swarming and social behaviour, is the collective random patterns and chaotic behaviour formed by some animal species, in particular if there is a danger. Is such a behaviour beneficial or unfavourable for survival? Here we report on one of the most remarkable forms of animal swarming and social behaviour - fish schooling - from a hydrodynamic point of view. We found that some fish species do not have preferred orientation and they swarm in a random pattern mode, despite the excess of energy consumed. Our analyses, which includes calculations of the hydrodynamic forces between slender bodies, show that such a behaviour enhances the transfer of hydrodynamic information, and thus enhances the survivability of the school. These findings support the general hypothesis that a disordered and non-trivial collective behaviour of individuals within a nonlinear dynamical system is essential for optimising transfer of information - an optimisation that might be crucial for survival.

preprint2017arXiv

Lagrangian Statistics for Navier-Stokes Turbulence under Fourier-mode reduction: Fractal and Homogeneous Decimations

We study small-scale and high-frequency turbulent fluctuations in three-dimensional flows under Fourier-mode reduction. The Navier-Stokes equations are evolved on a restricted set of modes, obtained as a projection on a fractal or homogeneous Fourier set. We find a strong sensitivity (reduction) of the high-frequency variability of the Lagrangian velocity fluctuations on the degree of mode decimation, similarly to what is already reported for Eulerian statistics. This is quantified by a tendency towards a quasi-Gaussian statistics, i.e., to a reduction of intermittency, at all scales and frequencies. This can be attributed to a strong depletion of vortex filaments and of the vortex stretching mechanism. Nevertheless, we found that Eulerian and Lagrangian ensembles are still connected by a dimensional bridge-relation which is independent of the degree of Fourier-mode decimation.

preprint2015arXiv

Continuum Nanofluidics

This paper introduces the fundamental continuum theory governing momentum transport in isotropic nanofluidic flows. The theory is an extension to the classical Navier-Stokes equation, which includes coupling between translational and rotational degrees of freedom, as well as non-local response functions that incorporates spatial correlations. The continuum theory is compared with molecular dynamics simulation data for both relaxation processes and fluid flows showing excellent agreement on the nanometer length scale. We also present practical tools to estimate when the extended theory should be used.

preprint2016arXiv

Self-similar formation of the Kolmogorov spectrum in the Leith model of turbulence

The last stage of evolution toward the stationary Kolmogorov spectrum of hydrodynamic turbulence is studied using the Leith model. This evolution is shown to manifest itself as a reflection wave in the wavenumber space propagating from the largest toward the smallest wavenumbers, and is described by a self-similar solution of a new (third) kind. This stage follows the previously studied stage of an initial explosive propagation of the spectral front from the smallest to the largest wavenumbers reaching arbitrarily large wavenumbers in a finite time, and which was described by a self-similar solution of the second kind. Nonstationary solutions corresponding to"warm cascades" characterised by a thermalised spectrum at large wavenumbers are also obtained.

preprint2016arXiv

Numerical simulations of the Princeton magneto-rotational instability experiment with conducting axial boundaries

We investigate numerically the Princeton magneto-rotational instability (MRI) experiment and the effect of conducting axial boundaries or endcaps. MRI is identified and found to reach a much higher saturation than for insulating endcaps. This is probably due to stronger driving of the base flow by the magnetically rather than viscously coupled boundaries. Although the computations are necessarily limited to lower Reynolds numbers ($\Re$) than their experimental counterparts, it appears that the saturation level becomes independent of $\Re$ when $\Re$ is sufficiently large, whereas it has been found previously to decrease roughly as $\Re^{-1/4}$ with insulating endcaps. The much higher saturation levels will allow for the first positive detection of MRI beyond its theoretical and numerical predictions.

preprint2016arXiv

Fast symmetric factorization of hierarchical matrices with applications

We present a fast direct algorithm for computing symmetric factorizations, i.e. $A = WW^T$, of symmetric positive-definite hierarchical matrices with weak-admissibility conditions. The computational cost for the symmetric factorization scales as $\mathcal{O}(n \log^2 n)$ for hierarchically off-diagonal low-rank matrices. Once this factorization is obtained, the cost for inversion, application, and determinant computation scales as $\mathcal{O}(n \log n)$. In particular, this allows for the near optimal generation of correlated random variates in the case where $A$ is a covariance matrix. This symmetric factorization algorithm depends on two key ingredients. First, we present a novel symmetric factorization formula for low-rank updates to the identity of the form $I+UKU^T$. This factorization can be computed in $\mathcal{O}(n)$ time if the rank of the perturbation is sufficiently small. Second, combining this formula with a recursive divide-and-conquer strategy, near linear complexity symmetric factorizations for hierarchically structured matrices can be obtained. We present numerical results for matrices relevant to problems in probability \& statistics (Gaussian processes), interp

preprint2016arXiv

An efficient threshold dynamics method for wetting on rough surfaces

The threshold dynamics method developed by Merriman, Bence and Osher (MBO) is an efficient method for simulating the motion by mean curvature flow when the interface is away from the solid boundary. Direct generalization of MBO-type methods to the wetting problem with interfaces intersecting the solid boundary is not easy because solving the heat equation in a general domain with a wetting boundary condition is not as efficient as it is with the original MBO method. The dynamics of the contact point also follows a different law compared with the dynamics of the interface away from the boundary. In this paper, we develop an efficient volume preserving threshold dynamics method for simulating wetting on rough surfaces. This method is based on minimization of the weighted surface area functional over an extended domain that includes the solid phase. The method is simple, stable with $O(N \log N)$ complexity per time step and is not sensitive to the inhomogeneity or roughness of the solid boundary.

preprint2016arXiv

The effect of rotation on oscillatory double-diffusive convection (semiconvection)

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermo-compositional layers which significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and non-layered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations we then analyze the effect of rotation on properties of layered and non-layered ODDC, and study how the angle of the rotation axis with re

preprint2017arXiv

PIV in stratified gas-liquid flow in a horizontal pipe using water droplets as tracers in the gas-phase

Simultaneous Particle Image Velocimetry (PIV) measurements of stratified turbulent air/water flow in a horizontal pipe have been performed using water droplets as tracers in the gas-phase. The use of water droplets as tracers ensures that the water surface tension remains unaffected and thus allows small-scale interfacial structures, such as capillary waves to occur naturally. Experiments have been conducted in a 31 m long, 100 mm diameter PVC pipe using air and water as test fluids. For the purpose of validation of the experimental set-up and the suggested seeding technique, single-phase measurements of both air and water were compared to each other and to DNS results from literature, showing very good agreement. The two-phase measurements are presented in terms of mean- and rms-profiles. These measurements offer a qualitative demonstration of the behavior of the interfacial turbulence and its correlation with the various interfacial flow patterns. The observations made in this paper are in agreement with the conclusions drawn from the DNS study of Lakehal et al. (2003). The present results may eventually provide a better explanation to many important phenomena related to the phys

preprint2016arXiv

Effects of a small magnetic field on homoclinic bifurcations in a low-Prandtl-number fluid

Effects of a uniform magnetic field on homoclinic bifurcations in Rayleigh-Bénard convection in a fluid of Prandtl number $Pr = 0.01$ are investigated using direct numerical simulations (DNS). A uniform magnetic field is applied either in the vertical or in the horizontal direction. For a weak vertical magnetic field, the possibilities of both forward and backward homoclinic bifurcations are observed leading to a spontaneous merging of two limit cycles into one as well as a spontaneous breaking of a limit cycle into two for lower values of the Chandrasekhar's number ($Q\leq 5$). A slightly stronger magnetic field makes the convective flow time independent giving the possibility of stationary patterns at the secondary instability. For horizontal magnetic field, the $x\leftrightharpoons y$ symmetry is destroyed and neither a homoclinic gluing nor a homoclinic breaking is observed. Two low-dimensional models are also constructed: one for a weak vertical magnetic field and another for a weak horizontal magnetic field. The models qualitatively capture the features observed in DNS and help understanding the unfolding of bifurcations close to the onset of magnetoconvection.

preprint2016arXiv

Mega electrorheological phenomena in graphene nanogels

Unprecedentedly massive electrorheology has been reported for dilute graphene nanoflakes based ER fluids that have been engineered as novel, readily synthesizable polymeric gels. Polyethylene glycol based graphene gels have been synthesized and very high ER response, 25,000 percent enhancement in viscosity under influence of electric field, has been observed for low concentration systems 2 wt. percent. The gels overcome several drawbacks innate to ER fluids. The gels exhibit long term stability, high graphene packing ratio which ensures very high ER response and the microstructure of the gels ensure that fibrillation of the graphene nanoflakes under field is undisturbed by thermal fluctuations, further leading to mega ER. The gels exhibit large yield stress handling caliber with yield stress observed as high as 13 kPa at 2 wt. perc. graphene. Detailed investigations on the effects of graphene concentration, electric field strength, imposed shear resistance and transients of electric field actuation on the ER response of the gels have been performed. The gels show nearly negligible hysteresis with respect to electro viscous effects which speak of the stability and responsiveness req

preprint2016arXiv

Shock-driven transition to turbulence: emergence of power-law scaling

We consider two cases of interaction between a planar shock and a cylindrical density interface. In the first case (planar normal shock), the axis of the gas cylinder is parallel to the shock front, and baroclinic vorticity deposited by the shock is predominantly two-dimensional (directed along the axis of the cylinder). In the second case, the cylinder is tilted, resulting in an oblique shock interaction, and a fully three-dimensional shock-induced vorticity field. The statistical properties of the flow for both cases are analyzed based on images from two orthogonal visualization planes, using structure functions of the intensity maps of fluorescent tracer pre-mixed with the heavy gas. At later times, these structure functions exhibit power-law-like behavior over a considerable range of scales. Manifestation of this behavior is remarkably consistent in terms of dimensionless time defined based on Richtmyer's linear theory within the range of Mach numbers from 1.1 to 2.0 and the range of gas cylinder tilt angles with respect to the plane of the shock front (0 to 30 degrees).

preprint2016arXiv

A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions

A new parallel, computationally efficient immersed boundary method for solving three-dimensional, viscous, incompressible flows on unbounded domains is presented. Immersed surfaces with prescribed motions are generated using the interpolation and regularization operators obtained from the discrete delta function approach of the original (Peskin's) immersed boundary method. Unlike Peskin's method, boundary forces are regarded as Lagrange multipliers that are used to satisfy the no-slip condition. The incompressible Navier-Stokes equations are discretized on an unbounded staggered Cartesian grid and are solved in a finite number of operations using lattice Green's function techniques. These techniques are used to automatically enforce the natural free-space boundary conditions and to implement a novel block-wise adaptive grid that significantly reduces the run-time cost of solutions by limiting operations to grid cells in the immediate vicinity and near-wake region of the immersed surface. These techniques also enable the construction of practical discrete viscous integrating factors that are used in combination with specialized half-explicit Runge-Kutta schemes to accura

preprint2017arXiv

Sharp Interface Limit for a Stokes/Allen-Cahn System

We consider the sharp interface limit of a coupled Stokes/Allen-Cahn system, when a parameter $\varepsilon>0$ that is proportional to the thickness of the diffuse interface tends to zero, in a two dimensional bounded domain. For sufficiently small times we prove convergence of the solutions of the Stokes/Allen-Cahn system to solutions of a sharp interface model, where the interface evolution is given by the mean curvature equation with an additional convection term coupled to a two-phase Stokes system with an additional contribution to the stress tensor, which describes the capillary stress. To this end we construct a suitable approximation of the solution of the Stokes/Allen-Cahn system, using three levels of the terms in the formally matched asymptotic calculations, and estimate the difference with the aid of a suitable refinement of a spectral estimate of the linearized Allen-Cahn operator. Moreover, a careful treatment of the coupling terms is needed.

preprint2016arXiv

Stabilizing effect of tip splitting on the interface motion

Pattern-forming processes, such as electrodeposition, dielectric breakdown or viscous fingering are often driven by instabilities. Accordingly, the resulting growth patterns are usually highly branched, fractal structures. However, in some of the unstable growth processes the envelope of the structure grows in a highly regular manner, with the perturbations smoothed out over the course of time. In this paper, we show that the regularity of the envelope growth can be connected to small-scale instabilities leading to the tip splitting of the fingers at the advancing front of the structure. Whenever the growth velocity becomes too large, the finger splits into two branches. In this way it can absorb an increased flux and thus damp the instability. Hence, somewhat counterintuitively, the instability at a small scale results in a stability at a larger scale. The quantitative analysis of these effects is provided by means of the Loewner equation, which one can use to reduce the problem of the interface motion to that of the evolution of the conformal mapping onto the complex plane. This allows an effective analysis of the multifingered growth in a variety of different geometries. We show

People in this topic

12 visible researcher(s)