Showing papers in "arXiv: Fluid Dynamics in 2015"
TL;DR: In this article, the authors review the biomechanics of bacterial motility and look ahead to future challenges using hydrodynamics as an organizing framework, including the ability of cells to reorient and search their surroundings to their interactions within mechanically and chemically-complex environments.
Abstract: Bacteria predate plants and animals by billions of years Today, they are the world's smallest cells yet they represent the bulk of the world's biomass, and the main reservoir of nutrients for higher organisms Most bacteria can move on their own, and the majority of motile bacteria are able to swim in viscous fluids using slender helical appendages called flagella Low-Reynolds-number hydrodynamics is at the heart of the ability of flagella to generate propulsion at the micron scale In fact, fluid dynamic forces impact many aspects of bacteriology, ranging from the ability of cells to reorient and search their surroundings to their interactions within mechanically and chemically-complex environments Using hydrodynamics as an organizing framework, we review the biomechanics of bacterial motility and look ahead to future challenges
318 citations
TL;DR: In this article, a two-stage dynamic mode decomposition (DMD) algorithm is proposed to remove bias errors by processing snapshots asymmetrically, and the resulting unbiased and noise-aware total DMD formulation reduces to standard DMD in the absence of snapshot errors.
Abstract: The Dynamic Mode Decomposition (DMD)---a popular method for performing data-driven Koopman spectral analysis---has gained increased adoption as a technique for extracting dynamically meaningful spatio-temporal descriptions of fluid flows from snapshot measurements. Often times, DMD descriptions can be used for predictive purposes as well, which enables informed decision-making based on DMD model-forecasts. Despite its widespread use and utility, DMD regularly fails to yield accurate dynamical descriptions when the measured snapshot data are imprecise due to, e.g., sensor noise. Here, we express DMD as a two-stage algorithm in order to isolate a source of systematic error. We show that DMD's first stage, a subspace projection step, systematically introduces bias errors by processing snapshots asymmetrically. To remove this systematic error, we propose utilizing an augmented snapshot matrix in a subspace projection step, as in problems of total least-squares, in order to account for the error present in all snapshots. The resulting unbiased and noise-aware total DMD (TDMD) formulation reduces to standard DMD in the absence of snapshot errors, while the two-stage perspective generalizes the de-biasing framework to other related methods as well. TDMD's performance is demonstrated in numerical and experimental fluids examples.
180 citations
TL;DR: Spectral Proper Orthogonal Decomposition (SPOD) as discussed by the authors is based on classical POD and it can be applied to spatially and temporally resolved data, and it allows for a continuous shifting from the energetically optimal POD to the spectrally pure Fourier decomposition by changing a single parameter.
Abstract: The identification of coherent structures from experimental or numerical data is an essential task when conducting research in fluid dynamics. This typically involves the construction of an empirical mode base that appropriately captures the dominant flow structures. The most prominent candidates are the energy-ranked proper orthogonal decomposition (POD) and the frequency ranked Fourier decomposition and dynamic mode decomposition (DMD). However, these methods fail when the relevant coherent structures occur at low energies or at multiple frequencies, which is often the case. To overcome the deficit of these "rigid" approaches, we propose a new method termed Spectral Proper Orthogonal Decomposition (SPOD). It is based on classical POD and it can be applied to spatially and temporally resolved data. The new method involves an additional temporal constraint that enables a clear separation of phenomena that occur at multiple frequencies and energies. SPOD allows for a continuous shifting from the energetically optimal POD to the spectrally pure Fourier decomposition by changing a single parameter. In this article, SPOD is motivated from phenomenological considerations of the POD autocorrelation matrix and justified from dynamical system theory. The new method is further applied to three sets of PIV measurements of flows from very different engineering problems. We consider the flow of a swirl-stabilized combustor, the wake of an airfoil with a Gurney flap, and the flow field of the sweeping jet behind a fluidic oscillator. For these examples, the commonly used methods fail to assign the relevant coherent structures to single modes. The SPOD, however, achieves a proper separation of spatially and temporally coherent structures, which are either hidden in stochastic turbulent fluctuations or spread over a wide frequency range.
140 citations
TL;DR: In this article, the differences and similarities between 2D and 3D Rayleigh-Benard convection are studied using Nu(Ra) for Pr = 4.38 and Pr = 0.7 and Nu(Pr) for Ra up to $10^8.
Abstract: Two dimensional (2D) and three dimensional (3D) Rayleigh-Benard convection is compared using results from direct numerical simulations and prior experiments. The explored phase diagrams for both cases are reviewed. The differences and similarities between 2D and 3D are studied using Nu(Ra) for Pr = 4.38 and Pr = 0.7 and Nu(Pr) for Ra up to $10^8$. In the Nu(Ra) scaling at higher Pr, 2D and 3D are very similar; differing only by a constant factor up to $\text{Ra} = 10^{10}$. In contrast, the difference is large at lower Pr, due to the strong roll state dependence of Nu in 2D. The behaviour of Nu(Pr) is similar in 2D and 3D at large Pr. However, it differs significantly around Pr = 1. The Reynolds number values are consistently higher in 2D and additionally converge at large Pr. Finally, the thermal boundary layer profiles are compared in 2D and 3D.
135 citations
TL;DR: In this article, the effect of preferential sampling on the particle statistics in the presence of gravity was examined over the largest range to date (from 88-597) and it was shown that preferential sampling effects are important for low and intermediate values of the Reynolds number.
Abstract: In this study, we analyze the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number $R_\lambda$ on the particle statistics is examined over the largest range to date (from $R_\lambda = 88-597$). We first explore the effect of preferential sampling on the single-particle statistics, and use our understanding of preferential sampling to provide a physical explanation for many of the trends in the particle velocity gradients, kinetic energies, and accelerations at low $St$. As $St$ increases, inertial filtering effects become more important, causing the particle kinetic energies and accelerations to decrease.
We then consider particle-pair statistics, and focus our attention on the relative velocities and radial distribution functions (RDFs) of the particles. The relative velocity statistics indicate that preferential-sampling effects are important for $St \lesssim 0.1$ and that path-history/non-local effects become increasingly important for $St \gtrsim 0.2$. The lower-order relative velocity statistics are only weakly sensitive to changes in Reynolds number at low $St$. We find that the RDFs peak near $St$ of order unity, that they exhibit power-law scaling for low and intermediate $St$, and that they are largely independent of Reynolds number for low and intermediate $St$. We also observe that at large $St$, changes in the RDF are related to changes the scaling exponents of the relative velocity variances. The particle collision kernel is found to be largely insensitive to the flow Reynolds number, suggesting that relatively low-Reynolds-number simulations may be able to capture much of the relevant physics of droplet collisions and growth in the adiabatic cores of atmospheric clouds.
133 citations
TL;DR: The results show that the transport and spread of chemotactic cells can be independently modified and controlled by the fluid material properties, and even small amounts of polymer in solution can drastically change E. coli dynamics.
Abstract: Run-and-tumble motility is widely used by swimming microorganisms including numerous prokaryotic eukaryotic organisms. Here, we experimentally investigate the run-and-tumble dynamics of the bacterium E. coli in polymeric solutions. We find that even small amounts of polymer in solution can drastically change E. coli dynamics: cells tumble less and their velocity increases, leading to an enhancement in cell translational diffusion and a sudden decline in rotational diffusion. We show that suppression of tumbling is due to fluid viscosity while the enhancement in swimming speed is mainly due to fluid elasticity. Visualization of single fluorescently labeled DNA polymers reveals that the flow generated by individual E. coli is sufficiently strong to stretch polymer molecules and induce elastic stresses in the fluid, which in turn can act on the cell in such a way to enhance its transport. Our results show that the transport and spread of chemotactic cells can be independently modified and controlled by the fluid material properties.
122 citations
TL;DR: In this article, the authors investigate the role of hydrodynamics in fish schooling by modeling swimmers as vortex dipoles whose interactions are governed by the Biot-Savart law.
Abstract: Schooling, an archetype of collective behavior, emerges from the interactions of fish responding to visual and other informative cues mediated by their aqueous environment. In this context, a fundamental and largely unexplored question concerns the role of hydrodynamics. Here, we investigate schooling by modeling swimmers as vortex dipoles whose interactions are governed by the Biot-Savart law. When we enhance these dipoles with behavioral rules from classical agent based models we find that they do not lead robustly to schooling due to flow mediated interactions. In turn, we present dipole swimmers equipped with adaptive decision-making that learn, through a reinforcement learning algorithm, to adjust their gaits in response to non-linearly varying hydrodynamic loads. The dipoles maintain their relative position within a formation by adapting their strength and school in a variety of prescribed geometrical arrangements. Furthermore, we identify schooling patterns that minimize the individual and the collective swimming effort, through an evolutionary optimization. The present work suggests that the adaptive response of individual swimmers to flow-mediated interactions is critical in fish schooling.
95 citations
TL;DR: This work is able to show that different types of surfactants turn the droplet's surface either rigid or elastic, which alters the evaporating fluid flow, either enhancing the classical coffee-stain effect or leading to a total flow inversion.
Abstract: An evaporating droplet is a dynamic system in which flow is spontaneously generated to minimize the surface energy, dragging particles to the borders and ultimately resulting in the so-called "coffee-stain effect". The situation becomes more complex at the droplet's surface, where surface tension gradients of different nature can compete with each other yielding different scenarios. With careful experiments and with the aid of 3D particle tracking techniques, we are able to show that different types of surfactants turn the droplet's surface either rigid or elastic, which alters the evaporating fluid flow, either enhancing the classical coffee-stain effect or leading to a total flow inversion. Our measurements lead to unprecedented and detailed measurements of the surface tension difference along an evaporating droplet's surface with good temporal and spatial resolution.
85 citations
Posted Content•
TL;DR: In this article, a lattice Boltzmann scheme is presented which exhibits an increased stability and accuracy with respect to standard singleor multi-relaxation-tion-time (MRT) approaches.
Abstract: In the present paper a lattice Boltzmann scheme is presented which exhibits an increased stability and accuracy with respect to standard singleor multi-relaxa-tion-time (MRT) approaches. The scheme is based on a single-relaxation-time model where a special regularization procedure is applied. This regularization is based on the fact that, for a-thermal ows, there exists a recursive way to express the velocity distribution function at any order (in the Hermite series sense) in terms of the density, velocity, and stress tensor. A linear stability analysis is conducted which shows enhanced dispersion/dissipation relations with respect to existing models. The model is then validated on two (one 2D and one 3D) moderately high Reynolds number simulations (Re 1000) at moderate Mach numbers (Ma 0:5). In both cases the results are compared with an MRT model and an enhanced accuracy and stability is shown by the present model.
75 citations
Posted Content•
TL;DR: In this article, the authors developed an analytical solution to the linearized flow response to a Gaussian lift and drag force and used the results to establish a relationship between the local and far-field velocity required to specify lift and charge forces, focusing first on the lift force and then finding the force center location that minimize the square difference between the velocity fields induced by the Gaussian force and 2D potential flow over Joukowski airfoils.
Abstract: The actuator line model (ALM) is a commonly used method to represent lifting surfaces such as wind turbine blades within Large-Eddy Simulations (LES) In the ALM the lift and drag forces are replaced by an imposed body force which is typically smoothed over several grid points using a Gaussian kernel with some prescribed smoothing width $\epsilon$ To date, the choice of $\epsilon$ has most often been based on numerical considerations related to the grid spacing used in LES However, especially for finely resolved LES with grid spacings on the order of or smaller than the chord-length of the blade, the best choice of $\epsilon$ is not known In this work, a theoretical approach is followed to determine the most suitable value of $\epsilon$ Firstly, we develop an analytical solution to the linearized flow response to a Gaussian lift and drag force and use the results to establish a relationship between the local and far-field velocity required to specify lift and drag forces Then, focusing first on the lift force, we find $\epsilon$ and the force center location that minimize the square difference between the velocity fields induced by the Gaussian force and 2D potential flow over Joukowski airfoils We find that the optimal smoothing width $\epsilon^{\rm opt}$ is on the order of 14-25\% of the chord length of the blade, and the center of force is located at about 13-26\% downstream of the leading edge of the blade, for the cases considered These optimal values do not depend on angle of attack and depend only weakly on the type of lifting surface To represent the drag force, the optimal width of the circular Gaussian drag force field is shown to be equal to the momentum thickness of the wake
72 citations
Posted Content•
TL;DR: The discrete unified gas kinetic scheme (DUGKS) is a finite volume method for deterministic solution of the Boltzmann model equation with asymptotic preserving property as mentioned in this paper.
Abstract: The recently proposed discrete unified gas kinetic scheme (DUGKS) is a finite volume method for deterministic solution of the Boltzmann model equation with asymptotic preserving property. In DUGKS, the numerical flux of the distribution function is determined from a local numerical solution of the Boltzmann model equation using an unsplitting approach. The time step and mesh resolution are not restricted by the molecular collision time and mean free path. To demonstrate the capacity of DUGKS in practical problems, this paper extends the DUGKS to arbitrary unstructured meshes. Several tests of both internal and external flows are performed, which include the cavity flow ranging from continuum to free molecular regimes, a multiscale flow between two connected cavities with a pressure ratio of 10000, and a high speed flow over a cylinder in slip and transitional regimes. The numerical results demonstrate the effectiveness of the DUGKS in simulating multiscale flow problems.
TL;DR: In this paper, the authors simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9.
Abstract: We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results.
TL;DR: In this paper, the authors investigate the motion of a sedimenting spherical drop in the presence of an applied uniform electric field in an otherwise arbitrary direction in the limit of low surface charge convection.
Abstract: We investigate the motion of a sedimenting spherical drop in the presence of an applied uniform electric field in an otherwise arbitrary direction in the limit of low surface charge convection. We analytically solve the electric potential in and around the leaky dielectric drop, and solve for the Stokesian velocity and pressure fields. We obtain the drop velocity through perturbations in powers of the electric Reynolds number which signifies the importance of the charge relaxation time scale as compared to the convective time scale. We show that in the presence of electric field either in the sedimenting direction or orthogonal to it, there is a change in the drop velocity only in the direction of sedimentation due to an asymmetric charge distribution in the same direction. However, in the presence of an electric field applied in both the directions, and depending on the permittivities and conductivities of the two fluids, we obtain a non-intuitive lateral migration of drop in addition to the buoyancy driven sedimentation. These dynamical features can be effectively used for manipulating drops in a controlled electro-fluidic environment.
TL;DR: In this article, a set of criteria which are applied to discuss various formulations under which Lagrangian stochastic models can be found are established. But the authors do not address the question of the relative predictive capacities of different models but concentrates on their formulation since advantages and disadvantages of different formulations are not always clear.
Abstract: In this paper, we establish a set of criteria which are applied to discuss various formulations under which Lagrangian stochastic models can be found. These models are used for the simulation of fluid particles in single-phase turbulence as well as for the fluid seen by discrete particles in dispersed turbulent two-phase flows. A central issue is to put forward reliable requirements which must be met by Lagrangian stochastic models and a new element brought by the present analysis is to address the single- and two-phase flow situations from a unified point of view.This analysis does not address the question of the relative predictive capacities of different models but concentrates on their formulation since advantages and disadvantages of different formulations are not always clear. Indeed, hidden in the changes from one structure to another are some possible pitfalls which can lead to flaws in the construction of practical models and to physically-unsound numerical calculations. A first interest of the present approach is illustrated by considering some models proposed in the literature and by showing that these criteria help to assess whether these Lagrangian stochastic models can be regarded as acceptable descriptions. A second interest is to indicate how future developments can be safely built, which is also relevant for stochastic subgrid models for particle-laden flows in the context of Large Eddy Simulations.
TL;DR: In this paper, the authors examined how increasing the value of the Reynolds number (Re$) affects the ability of spanwise-forcing techniques to yield turbulent skin-friction drag reduction.
Abstract: This paper examines how increasing the value of the Reynolds number $Re$ affects the ability of spanwise-forcing techniques to yield turbulent skin-friction drag reduction. The considered forcing is based on the streamwise-travelling waves of spanwise wall velocity (Quadrio {\em et al. J. Fluid Mech.}, vol. 627, 2009, pp. 161--178). The study builds upon an extensive drag-reduction database created with Direct Numerical Simulation of a turbulent channel flow for two, 5-fold separated values of $Re$, namely $Re_\tau=200$ and $Re_\tau=1000$. The sheer size of the database, which for the first time systematically addresses the amplitude of the forcing, allows a comprehensive view of the drag-reducing characteristics of the travelling waves, and enables a detailed description of the changes occurring when $Re$ increases. The effect of using a viscous scaling based on the friction velocity of either the non-controlled flow or the drag-reduced flow is described. In analogy with other wall-based drag reduction techniques, like for example riblets, the performance of the travelling waves is well described by a vertical shift of the logarithmic portion of the mean streamwise velocity profile. Except when $Re$ is very low, this shift remains constant with $Re$, at odds with the percentage reduction of the friction coefficient, which is known to present a mild, logarithmic decline. Our new data agree with the available literature, which is however mostly based on low-$Re$ information and hence predicts a quick drop of maximum drag reduction with $Re$. The present study supports a more optimistic scenario, where for an airplane at flight Reynolds numbers a drag reduction of nearly 30\% would still be possible thanks to the travelling waves.
Posted Content•
TL;DR: In this article, a double multiple-relaxation-time lattice Boltzmann model is developed for simulating transient solid-liquid phase change problems in porous media at the representative elementary volume scale.
Abstract: In this paper, a double multiple-relaxation-time lattice Boltzmann model is developed for simulating transient solid-liquid phase change problems in porous media at the representative elementary volume scale. The model uses two different multiple-relaxation-time lattice Boltzmann equations, one for the flow field and the other for the temperature field with nonlinear latent heat source term. The model is based on the generalized non-Darcy formulation, and the solid-liquid phase change interface is traced through the liquid fraction which is determined by the enthalpy method. The model is validated by numerical simulations of conduction melting in a semi-infinite space, solidification in a semi-infinite corner, and convection melting in a square cavity filled with porous media. The numerical results demonstrate the efficiency and accuracy of the present model for simulating transient solid-liquid phase change problems in porous media.
Posted Content•
TL;DR: In this article, the transition and turbulence production in a hypersonic boundary layer was investigated in a Mach 6 quiet wind tunnel using Rayleigh-scattering visualization, fast-response pressure measurements, and particle image velocimetry.
Abstract: Transition and turbulence production in a hypersonic boundary layer is investigated in a Mach 6 quiet wind tunnel using Rayleigh-scattering visualization, fast-response pressure measurements, and particle image velocimetry. It is found that the second instability acoustic mode is the key modulator of the transition process. The second mode experiences a rapid growth and a very fast annihilation due to the effect of bulk viscosity. The second mode interacts strongly with the first vorticity mode to directly promote a fast growth of the latter and leads to immediate transition to turbulence.
TL;DR: In this paper, an effective equation of motion for the orientational dynamics of a neutrally buoyant spheroid suspended in a simple shear flow was derived, valid for arbitrary particle aspect ratios and to linear order in the shear Reynolds number.
Abstract: We derive an effective equation of motion for the orientational dynamics of a neutrally buoyant spheroid suspended in a simple shear flow, valid for arbitrary particle aspect ratios and to linear order in the shear Reynolds number. We show how inertial effects lift the degeneracy of the Jeffery orbits and determine the stabilities of the log-rolling and tumbling orbits at infinitesimal shear Reynolds numbers. For prolate spheroids we find stable tumbling in the shear plane, log-rolling is unstable. For oblate particles, by contrast, log-rolling is stable and tumbling is unstable provided that the aspect ratio is larger than a critical value. When the aspect ratio is smaller than this value tumbling turns stable, and an unstable limit cycle is born.
TL;DR: In this article, the authors conduct direct numerical simulations of spray formation in a gas-liquid mixing layer with state-of-the-art computational resources (using up to 4 billion cells and 16384 cores).
Abstract: Spray formation and atomization in a gas-liquid mixing layer is an important fundamental problem of multiphase flows. It is highly desirable to visualize the detailed atomization process and to analyze the instabilities and mechanisms involved, and massive numerical simulations are required, in addition to experiment. Rapid development of numerical methods and computer technology in the past decades now allows large-scale three-dimensional direct numerical simulations of atomization to be performed. Nevertheless, the fundamental question, whether all the physical scales involved in the primary breakup process are faithfully resolved, remains unclear. In the present study, we conduct direct numerical simulations of spray formation in a gas-liquid mixing layer with state-of-the-art computational resources (using up to 4 billion cells and 16384 cores), in order to obtain a high-fidelity numerical closeup of the detailed mechanisms of spray formation. We also aim to examine whether present computational resources are sufficient for a fully resolved direct numerical simulation of atomization.
TL;DR: This chapter presents a theoretical overview of small-scale locomotion in complex fluids with a focus on recent efforts quantifying the impact of non-Newtonian rheology on swimming microorganisms.
Abstract: Microorganisms such as bacteria often swim in fluid environments that cannot be classified as Newtonian. Many biological fluids contain polymers or other heterogeneities which may yield complex rheology. For a given set of boundary conditions on a moving organism, flows can be substantially different in complex fluids, while non-Newtonian stresses can alter the gait of the microorganisms themselves. Heterogeneities in the fluid may also be characterized by length scales on the order of the organism itself leading to additional dynamic complexity. In this chapter we present a theoretical overview of small-scale locomotion in complex fluids with a focus on recent efforts quantifying the impact of non-Newtonian rheology on swimming microorganisms.
TL;DR: In this article, the effect of varying the mass and volume fraction of a suspension of rigid spheres dispersed in a turbulent channel flow was studied and the main effect at density ratios of the order of 10 was a strong shear-induced migration towards the centerline of the channel.
Abstract: We study the effect of varying the mass and volume fraction of a suspension of rigid spheres dispersed in a turbulent channel flow. We performed several Direct Numerical Simulations using an Immersed Boundary Method for finite-size particles changing the solid to fluid density ratio, the mass fraction and the volume fraction. We find that varying the density ratio between 1 and 10 at constant volume fraction does not alter the flow statistics as much as when varying the volume fraction at constant and at constant mass fraction. Interestingly, the increase in overall drag found when varying the volume fraction is considerably higher than that obtained for increasing density ratios at same volume fraction. The main effect at density ratios of the order of 10 is a strong shear-induced migration towards the centerline of the channel. When the density ratio is further increased up to 1000 the particle dynamics decouple from that of the fluid. The solid phase behaves as a dense gas and the fluid and solid phase statistics drastically change. In this regime, the collision rate is high and dominated by the normal relative velocity among particles.
TL;DR: In this paper, Xi et al. presented new families of nonlinear three-dimensional travelling wave solutions to the Navier-Stokes equation, also known as exact coherent states, for Newtonian plane Poiseuille flow.
Abstract: Several new families of nonlinear three-dimensional travelling wave solutions to the Navier-Stokes equation, also known as exact coherent states, are computed for Newtonian plane Poiseuille flow. The symmetries and streak/vortex structures are reported and their possible connections to critical layer dynamics examined. While some of the solutions clearly display fluctuations that are localized around the critical layer (the surface on which the streamwise velocity matches the wave speed of the solution), for others this connection is not as clear. Dynamical trajectories along unstable directions of the solutions are computed. Over certain ranges of Reynolds number, two solution families are shown to lie on the basin boundary between laminar and turbulent flow. Direct comparison of nonlinear travelling wave solutions to turbulent flow in the same channel is presented. The state-space dynamics of the turbulent flow are organized around one of the newly-identified travelling wave families, and in particular the lower branch solutions of this family are closely approached during transient excursions away from the dominant behaviour. These observations provide a firm dynamical-systems foundation for prior observations that minimal channel turbulence displays time intervals of "active" turbulence punctuated by brief periods of "hibernation" (see e.g. Xi, L. and Graham, M. D., Phys. Rev. Lett., 104, 218301 (2010)). The hibernating intervals are approaches to lower branch nonlinear travelling waves. Representing these solutions on a Prandtl-von Karman plot illustrates how their bulk flow properties are related to those of Newtonian turbulence as well as the universal asymptotic state called maximum drag reduction (MDR) found in viscoelastic turbulent flow.
[...]
TL;DR: In this article, the results of a combined experimental and numerical study of the generation of internal waves using the novel internal wave generator design of Gostiaux et al. were presented.
Abstract: We present the results of a combined experimental and numerical study of the generation of internal waves using the novel internal wave generator design of Gostiaux et al. (2007). This mechanism, which involves a tunable source comprised of oscillating plates, has so far been used for a few fundamental studies of internal waves, but its full potential has yet to be realized. Our studies reveal that this approach is capable of producing a wide variety of two-dimensional wave fields, including plane waves, wave beams and discrete vertical modes in finite-depth stratifications. The effects of discretization by a finite number of plates, forcing amplitude and angle of propagation are investigated, and it is found that the method is remarkably efficient at generating a complete wave field despite forcing only one velocity component in a controllable manner. We furthermore find that the nature of the radiated wave field is well predicted using Fourier transforms of the spatial structure of the wave generator.
TL;DR: Two different approaches for calculating the effective diffusion coefficient as a function of the microstructure are presented, including a deterministic approach based on the method of multiple scales and a stochastic approach for small volume fraction of spheres based on matched asymptotic expansions.
Abstract: The problem of diffusion in a porous medium with a spatially varying porosity is considered. The particular microstructure analyzed comprises a collection of impenetrable spheres, though the methods developed are general. Two different approaches for calculating the effective diffusion coefficient as a function of the microstructure are presented. The first is a deterministic approach based on the method of multiple scales; the second is a stochastic approach for small volume fraction of spheres based on matched asymptotic expansions. We compare the two approaches, and we show good agreement between them in a number of example configurations.
TL;DR: In this article, the role of fractal clustering of particles and caustic singularities of their velocities in collisional aggregation due to turbulence has been discussed, leading to a semi-quantitative understanding on the influence of turbulence on collision rates, and point to deficiencies in current understanding of rainfall and planet formation.
Abstract: Collisions between particles suspended in a fluid play an important role in many physical processes. As an example, collisions of microscopic water droplets in clouds are a necessary step in the production of macroscopic raindrops. Collisions of dust grains are also conjectured to be important for planet formation in the gas surrounding young stars, and also to play a role in the dynamics of sand storms. In these processes, collisions are favoured by fast turbulent motions. Here we review recent advances in the understanding of collisional aggregation due to turbulence. We discuss the role of fractal clustering of particles, and caustic singularities of their velocities. We also discuss limitations of the Smoluchowski equation for modelling these processes. These advances lead to a semi-quantitative understanding on the influence of turbulence on collision rates, and point to deficiencies in the current understanding of rainfall and planet formation.
TL;DR: In this article, the early stages of drop impact onto a solid surface are considered and detailed numerical simulations and detailed asymptotic analysis of the process reveal a self-similar structure both for the velocity field and the pressure field.
Abstract: The early stages of drop impact onto a solid surface are considered. Detailed numerical simulations and detailed asymptotic analysis of the process reveal a self-similar structure both for the velocity field and the pressure field. The latter is shown to exhibit a maximum not near the impact point, but rather at the contact line. The motion of the contact line is furthermore shown to exhibit a 'tank treading' motion. These observations are apprehended at the light of a variant of Wagner theory for liquid impact. This framework offers a simple analogy where the fluid motion within the impacting drop may be viewed as the flow induced by a flat rising expanding disk. The theoretical predictions are found to be in very close agreement both qualitatively and quantitatively with the numerical observations for about three decades in time. Interestingly the inviscid self-similar impact pressure and velocities are shown to depend solely on the self-similar variables $(r/\sqrt{t},z/\sqrt{t})$. The structure of the boundary layer developing along the wet substrate is investigated as well, and is proven to be formally analogous to that of the boundary layer growing in the trail of a shockwave. Interestingly, the boundary layer structure only depends on the impact self-similar variables. This allows to construct a seamless uniform analytical solution encompassing both impact and viscous effects. The depiction of the different dynamical fields enables to quantitatively predict observables of interest, such as the evolution of the integral viscous shearing force and of the net normal force.
TL;DR: In this paper, the second order closure of the dynamics is used to study the structure and dynamics of turbulence in plane Poiseuille flow at moderately high Reynolds numbers in a closely related dynamical system, referred to as the restricted nonlinear (RNL) system.
Abstract: The perspective of statistical state dynamics (SSD) has recently been applied to the study of mechanisms underlying turbulence in various physical systems. An example implementation of SSD is the second order closure referred to as stochastic structural stability theory (S3T), which has provided insight into the dynamics of wall turbulence and specifically the emergence and maintenance of the roll/streak structure. This closure eliminates nonlinear interactions among the perturbations has been removed, restricting nonlinearity in the dynamics to that of the mean equation and the interaction between the mean and perturbation covariance. Here, this quasi-linear restriction of the dynamics is used to study the structure and dynamics of turbulence in plane Poiseuille flow at moderately high Reynolds numbers in a closely related dynamical system, referred to as the restricted nonlinear (RNL) system. RNL simulations reveal that the essential features of wall-turbulence dynamics are retained. Remarkably, the RNL system spontaneously limits the support of its turbulence to a small set of streamwise Fourier components giving rise to a naturally minimal representation of its turbulence dynamics. Although greatly simplified, this RNL turbulence exhibits natural-looking structures and statistics. Surprisingly, even when further truncation of the perturbation support to a single streamwise component is imposed the RNL system continues to produce self-sustaining turbulent structure and dynamics. RNL turbulence at the Reynolds numbers studied is dominated by the roll/streak structure in the buffer layer and similar very-large-scale structure (VLSM) in the outer layer. Diagnostics of the structure, spectrum and energetics of RNL and DNS turbulence are used to demonstrate that the roll/streak dynamics supporting the turbulence in the buffer and logarithmic layer is essentially similar in RNL and DNS.
Abstract: Compressing a porous material will decrease the volume of the pore space, driving fluid out. Similarly, injecting fluid into a porous material can expand the pore space, distorting the solid skeleton. This poromechanical coupling has applications ranging from tissue mechanics to hydrogeology. The classical theory of linear poroelasticity captures this coupling by combining Darcy's law with Terzaghi's effective stress and linear elasticity in a linearized kinematic framework. Linear poroelasticity is a good model for very small deformations, but it becomes increasingly inappropriate for moderate to large deformations, which are common in the context of phenomena such as swelling and damage, and for soft materials such as gels and tissues. The well-known theory of large-deformation poroelasticity combines Darcy's law with Terzaghi's effective stress and nonlinear elasticity in a rigorous kinematic framework. This theory has been used extensively in biomechanics to model large elastic deformations in soft tissues, and in geomechanics to model large elastoplastic deformations in soils. Here, we first provide an overview and discussion of this theory with an emphasis on the physics of poromechanical coupling. We present the large-deformation theory in an Eulerian framework to minimize the mathematical complexity, and we show how this nonlinear theory simplifies to linear poroelasticity under the assumption of small strain. We then compare the predictions of linear poroelasticity with those of large-deformation poroelasticity in the context of two uniaxial model problems: Fluid outflow driven by an applied mechanical load (the consolidation problem) and compression driven by a steady fluid throughflow. We use these problems to explore the steady and dynamical errors associated with the linear model in both situations, as well as the impact of introducing a deformation-dependent permeability.
TL;DR: In this article, stable but generally unstable solutions of the 2D Boussinesq equations are obtained for no-slip boundary conditions and Prandtl number 7, with a delicate spiral structure in the temperature field.
Abstract: Steady but generally unstable solutions of the 2D Boussinesq equations are obtained for no-slip boundary conditions and Prandtl number 7. The primary solution that bifurcates from the conduction state at Rayleigh number $Ra \approx 1708$ has been calculated up to $Ra\approx 5. 10^6$ and shows heat flux $Nu \sim 0.143\, Ra^{0.28}$ with a delicate spiral structure in the temperature field. Another solution that maximizes $Nu$ over the horizontal wavenumber has been calculated up to $Ra=10^9$ and its heat flux scales as $Nu \sim 0.115\, Ra^{0.31}$ for $10^7 < Ra \le 10^9$, quite similar to 3D turbulent data. The latter is a simple yet multi-scale coherent solution whose horizontal wavenumber scales as $0.133 \, Ra^{0.217}$ in that range. That optimum solution is unstable to larger scale perturbations and in particular to mean shear flows, yet it appears to be relevant as a backbone for turbulent solutions, possibly setting the scale, strength and spacing of elemental plumes.
TL;DR: In this article, a simulation of the Taylor-Couette (TC) problem with mild curvature was performed, and the results showed that small-gap TC is dominated by spatially fixed large-scale structures known as Taylor rolls (TRs).
Abstract: Direct numerical simulations of the Taylor-Couette (TC) problem, the flow between two coaxial and independently rotating cylinders, have been performed. The study focuses on TC flow with mild curvature (small gap) with a radius ratio of $\eta=r_i/r_o=0.909$, an aspect ratio of $\Gamma=L/d=2\pi/3$, and a stationary outer cylinder. Three inner cylinder Reynolds of $1\cdot10^5$, $2\cdot10^5$ and $3\cdot 10^5$ were simulated, corresponding to frictional Reynolds numbers between $Re_\tau\approx 1400$ and $Re_\tau \approx 4000$. An additional case with a large gap, $\eta=0.5$ and driving of $Re=2\cdot10^5$ was also performed. Small-gap TC was found to be dominated by spatially-fixed large-scale structures, known as Taylor rolls (TRs). TRs are attached to the boundary layer, and are active, i.e. they transport angular velocity through Reynolds stresses. An additional simulation with inner cylinder Reynolds number of $Re=1\cdot10^5$ and fixed outer cylinder with an externally imposed axial flow of comparable strength as the wind of the TRs was also conducted. The axial flow was found to convect the TRs without any weakening effect. For small-gap TC, evidence for the existence of logarithmic velocity fluctuations, and of an overlap layer, in which the velocity fluctuations collapse in outer units, was found. Profiles consistent with a logarithmic dependence were also found for the angular velocity in large-gap TC, albeit in a very reduced range of scales. Finally, the behaviour of both small- and large-gap TC was compared to other canonical flows. Small-gap TC has similar behaviour in the near-wall region to other canonical flows, while large-gap TC displays very different behaviour.