Showing papers in "arXiv: Fluid Dynamics in 2016"

A Physics Informed Machine Learning Approach for Reconstructing Reynolds Stress Modeling Discrepancies Based on DNS Data

Abstract: Turbulence modeling is a critical component in numerical simulations of industrial flows based on Reynolds-averaged Navier-Stokes (RANS) equations. However, after decades of efforts in the turbulence modeling community, universally applicable RANS models with predictive capabilities are still lacking. Recently, data-driven methods have been proposed as a promising alternative to the traditional approaches of turbulence model development. In this work we propose a data-driven, physics-informed machine learning approach for predicting discrepancies in RANS modeled Reynolds stresses. The discrepancies are formulated as functions of the mean flow features. By using a modern machine learning technique based on random forests, the discrepancy functions are first trained with benchmark flow data and then used to predict Reynolds stresses discrepancies in new flows. The method is used to predict the Reynolds stresses in the flow over periodic hills by using two training flow scenarios of increasing difficulties: (1) the flow in the same periodic hills geometry yet at a lower Reynolds number, and (2) the flow in a different hill geometry with a similar recirculation zone. Excellent predictive performances were observed in both scenarios, demonstrating the merits of the proposed method. Improvement of RANS modeled Reynolds stresses enabled by the proposed method is an important step towards predictive turbulence modeling, where the ultimate goal is to predict the quantities of interest (e.g., velocity field, drag, lift) more accurately by solving RANS equations with the Reynolds stresses obtained therefrom.

A Computational Method for Sharp Interface Advection

Abstract: We devise a numerical method for passive advection of a surface, such as the interface between two incompressible fluids, across a computational mesh. The method is called isoAdvector, and is developed for general meshes consisting of arbitrary polyhedral cells. The algorithm is based on the volume of fluid (VOF) idea of calculating the volume of one of the fluids transported across the mesh faces during a time step. The novelty of the isoAdvector concept consists in two parts: First, we exploit an isosurface concept for modelling the interface inside cells in a geometric surface reconstruction step. Second, from the reconstructed surface, we model the motion of the face-interface intersection line for a general polygonal face to obtain the time evolution within a time step of the submerged face area. Integrating this submerged area over the time step leads to an accurate estimate for the total volume of fluid transported across the face. The method was tested on simple 2D and 3D interface advection problems both on structured and unstructured meshes. The results are very satisfactory both in terms of volume conservation, boundedness, surface sharpness, and efficiency. The isoAdvector method was implemented as an OpenFOAM(R) extension and is published as open source.

On the spreading of impacting drops

Abstract: The energy budget and dissipation mechanisms during droplet impact on solid surfaces are studied numerically and theoretically. We find that for high impact velocities and negligible surface friction at the solid surface (i.e. free-slip), about one half of the initial kinetic energy is transformed into surface energy, independent of the impact parameters and the detailed energy loss mechanism(s). We argue that this seemingly universal rule is related to the deformation mode of the droplet and is reminiscent of pipe flow undergoing a sudden expansion, for which the head loss can be calculated by multiplying the kinetic energy of the incoming flow by a geometrical factor. For impacts on a no-slip surface also dissipation in the shear boundary layer at the solid surface is important. In this case the geometric head loss acts as a lower bound on the total dissipation (i.e. the spreading on a no-slip surface approaches that on a free-slip surface when the droplet viscosity is send to zero). This new view on the impact problem allows for simple analytical estimates of the maximum spreading diameter of impacting drops as a function of the impact parameters and the properties of the solid surface. It bridges the gap between previous momentum balance approaches and energy balance approaches, which hitherto did not give consistent predictions in the low viscosity limit. Good agreement is found between our models and experiments, both for impacts on "slippery" or lubricated surfaces (e.g. Leidenfrost droplet impacts and head-on droplet-droplet collisions) and for impacts on no-slip surfaces.

Constrained Sparse Galerkin Regression

Abstract: Although major advances have been achieved over the past decades for the reduction and identification of linear systems, deriving nonlinear low-order models still is a chal- lenging task. In this work, we develop a new data-driven framework to identify nonlinear reduced-order models of a fluid by combining dimensionality reductions techniques (e.g. proper orthogonal decomposition) and sparse regression techniques from machine learn- ing. In particular, we extend the sparse identification of nonlinear dynamics (SINDy) algorithm to enforce physical constraints in the regression, namely energy-preserving quadratic nonlinearities. The resulting models, hereafter referred to as Galerkin regression models, incorporate many beneficial aspects of Galerkin projection, but without the need for a full-order or high-fidelity solver to project the Navier-Stokes equations. Instead, the most parsimonious nonlinear model is determined that is consistent with observed mea- surement data and satisfies necessary constraints. Galerkin regression models also readily generalize to include higher-order nonlinear terms that model the effect of truncated modes. The effectiveness of Galerkin regression is demonstrated on two different flow configurations: the two-dimensional flow past a circular cylinder and the shear-driven cavity flow. For both cases, the accuracy of the identified models compare favorably against reduced-order models obtained from a standard Galerkin projection procedure. Present results highlight the importance of cubic nonlinearities in the construction of accurate nonlinear low-dimensional approximations of the flow systems, something which cannot be readily obtained using a standard Galerkin projection of the Navier-Stokes equations. Finally, the entire code base for our constrained sparse Galerkin regression algorithm is freely available online.

Coalescence of sessile drops

Abstract: We present an experimental and theoretical description of the kinetics of coalescence of two water drops on a plane solid surface The case of partial wetting is considered The drops are in an atmosphere of nitrogen saturated with water where they grow by condensation and eventually touch each other and coalesce A new convex composite drop is rapidly formed that then exponentially and slowly relaxes to an equilibrium hemispherical cap The characteristic relaxation time is proportional to the drop radius R * at final equilibrium This relaxation time appears to be nearly 10 7 times larger than the bulk capillary relaxation time t b = R * $\eta$/$\sigma$, where $\sigma$ is the gas--liquid surface tension and $\eta$ is the liquid shear viscosity In order to explain this extremely large relaxation time, we consider a model that involves an Arrhenius kinetic factor resulting from a liquid--vapour phase change in the vicinity of the contact line The model results in a large relaxation time of order t b exp(L/RT) where L is the molar latent heat of vaporization, R is the gas constant and T is the temperature We model the late time relaxation for a near spherical cap and find an exponential relaxation whose typical time scale agrees reasonably well with the experiment 1 Introduction Fusion or coalescence between drops is a key process in a wide range of phenomena: phase transition in fluids and liquid mixtures or polymers, stability of foams and emulsions, and sintering in metallurgy (Eggers 1998), which is why the problem of coalescence has already received considerable attention Most of the studies of this process so far have been devoted to the coalescence of two spherical drops floating in a medium The kinetics of the process before and after the drops have touched each other is governed by the hydrodynamics inside and outside the drops and by the van der Waals forces when the drops are within mesoscopic distance from each other (Yiantsios \& Davis 1991) The composite drop that results from the coalescence of two drops relaxes to a spherical shape within a time which is dominated by the relaxation of the flow inside and outside (Nikolayev, Beysens \& Guenoun 1996; Nikolayev \& Beysens 1997) There are no studies, to our knowledge, of the coalescence of two sessile drops after they touch each other In this paper, we report a preliminary study of the dynamics and morphology of this process, in the case of hemispherical water droplets which grow slowly on a plane surface at the expense of the surrounding atmosphere, forming what is called 'dew' or 'breath figures' (Beysens et al 1991; Beysens 1995) The drops eventually touch each other and coalesce to form an elongated composite

Color of turbulence

Abstract: In this paper, we address the problem of how to account for second-order statistics of turbulent flows using low-complexity stochastic dynamical models based on the linearized Navier-Stokes equations. The complexity is quantified by the number of degrees of freedom in the linearized evolution model that are directly influenced by stochastic excitation sources. For the case where only a subset of velocity correlations are known, we develop a framework to complete unavailable second-order statistics in a way that is consistent with linearization around turbulent mean velocity. In general, white-in-time stochastic forcing is not sufficient to explain turbulent flow statistics. We develop models for colored-in-time forcing using a maximum entropy formulation together with a regularization that serves as a proxy for rank minimization. We show that colored-in-time excitation of the Navier-Stokes equations can also be interpreted as a low-rank modification to the generator of the linearized dynamics. Our method provides a data-driven refinement of models that originate from first principles and captures complex dynamics of turbulent flows in a way that is tractable for analysis, optimization, and control design.

Continuous Time Random Walks for the Evolution of Lagrangian Velocities

Abstract: We develop a continuous time random walk (CTRW) approach for the evolution of Lagrangian velocities in steady heterogeneous flows based on a stochastic relaxation process for the streamwise particle velocities. This approach describes persistence of velocities over a characteristic spatial scale, unlike classical random walk methods, which model persistence over a characteristic time scale. We first establish the relation between Eulerian and Lagrangian velocities for both equidistant and isochrone sampling along streamlines, under transient and stationary conditions. Based on this, we develop a space continuous CTRW approach for the spatial and temporal dynamics of Lagrangian velocities. While classical CTRW formulations have non-stationary Lagrangian velocity statistics, the proposed approach quantifies the evolution of the Lagrangian velocity statistics under both stationary and non-stationary conditions. We provide explicit expressions for the Lagrangian velocity statistics, and determine the behaviors of the mean particle velocity, velocity covariance and particle dispersion. We find strong Lagrangian correlation and anomalous dispersion for velocity distributions which are tailed toward low velocities as well as marked differences depending on the initial conditions. The developed CTRW approach predicts the Lagrangian particle dynamics from an arbitrary initial condition based on the Eulerian velocity distribution and a characteristic correlation scale.

On the correspondence between Koopman mode decomposition, resolvent mode decomposition, and invariant solutions of the Navier-Stokes equations

Abstract: The relationship between Koopman mode decomposition, resolvent mode decomposition and exact invariant solutions of the Navier-Stokes equations is clarified. The correspondence rests upon the invariance of the system operators under symmetry operations such as spatial translation. The usual interpretation of the Koopman operator is generalised to permit combinations of such operations, in addition to translation in time. This invariance is related to the spectrum of a spatio-temporal Koopman operator, which has a travelling wave interpretation. The relationship leads to a generalisation of dynamic mode decomposition, in which symmetry operations are applied to restrict the dynamic modes to span a subspace subject to those symmetries. The resolvent is interpreted as the mapping between the Koopman modes of the Reynolds stress divergence and the velocity field. It is shown that the singular vectors of the resolvent (the resolvent modes) are the optimal basis in which to express the velocity field Koopman modes where the latter are not a priori known.

Optimal undulating swimming for a single fish-like body and for a pair of interacting swimmers

Abstract: We establish through numerical simulation conditions for optimal undulatory propulsion for a single fish, and for a pair of hydrodynamically interacting fish, accounting for linear and angular recoil. We first employ systematic 2D simulations to identify conditions for minimal propulsive power of a self-propelled fish, and continue with targeted 3D simulations for a danio-like fish. We find that the Strouhal number, phase angle between heave and pitch at the trailing edge, and angle of attack are principal parameters. Angular recoil has significant impact on efficiency, while optimized body bending requires maximum bending amplitude upstream of the trailing edge. For 2D simulations, imposing a deformation based on measured displacement for carangiform swimming provides efficiency of 40%, which increases for an optimized profile to 57%; for a 3D fish, the corresponding increase is from 22% to 35%; all at Reynolds number 5000. Next, we turn to 2D simulation of two hydrodynamically interacting fish. We find that the upstream fish benefits energetically only for small distances. In contrast, the downstream fish can benefit at any position that allows interaction with the upstream wake, provided its body motion is timed appropriately with respect to the oncoming vortices. For an in-line configuration, one body length apart, the optimal efficiency of the downstream fish can increase to 66%; for an offset arrangement it can reach 81%. This proves that in groups of fish, energy savings can be achieved for downstream fish through interaction with oncoming vortices, even when the downstream fish lies directly inside the jet-like flow of an upstream fish.

Phenomenology of buoyancy-driven turbulence: recent results

Abstract: In this paper, we review the recent developments in the field of buoyancy-driven turbulence. Scaling and numerical arguments show that the stably-stratified turbulence with moderate stratification has kinetic energy spectrum $E_u(k) \sim k^{-11/5}$ and the kinetic energy flux $\Pi_u(k) \sim k^{-4/5}$, which is called Bolgiano-Obukhov scaling. The energy flux for the Rayleigh-Benard convection (RBC) however is approximately constant in the inertial range that results in Kolmorogorv's spectrum ($E_u(k) \sim k^{-5/3}$) for the kinetic energy. The phenomenology of RBC should apply to other flows where the buoyancy feeds the kinetic energy, e.g. bubbly turbulence and fully-developed Rayleigh Taylor instability. This paper also covers several models that predict the Reynolds and Nusselt numbers of RBC. Recent works show that the viscous dissipation rate of RBC scales as $\sim \mathrm{Ra}^{1.3}$, where $\mathrm{Ra}$ is the Rayleigh number.

Coherent structure coloring: identification of coherent structures from sparse data using graph theory

Abstract: We present a frame-invariant method for detecting coherent structures from Lagrangian flow trajectories that can be sparse in number, as is the case in many fluid mechanics applications of practical interest. The method, based on principles used in graph coloring and spectral graph drawing algorithms, examines a measure of the kinematic dissimilarity of all pairs of fluid trajectories, either measured experimentally, e.g. using particle tracking velocimetry; or numerically, by advecting fluid particles in the Eulerian velocity field. Coherence is assigned to groups of particles whose kinematics remain similar throughout the time interval for which trajectory data is available, regardless of their physical proximity to one another. Through the use of several analytical and experimental validation cases, this algorithm is shown to robustly detect coherent structures using significantly less flow data than is required by existing spectral graph theory methods.

Wetting theory for small droplets on textured solid surfaces

Abstract: Conventional wetting theories on rough surfaces with Wenzel, Cassie-Baxter, and Penetrate modes suggest the possibility of tuning the contact angle by adjusting the surface texture. Despite decades of intensive study, there are still many experimental results that are not well understood because conventional wetting theory, which assume an infinite droplet size, has been used to explain measurements of finite-sized droplets. In this study, we suggest a wetting theory that is applicable to any droplet size based on the free energy landscape analysis of various wetting modes of finite-sized droplets on a 2D textured surface. The key finding of our study is that there are many quantized wetting angles with local free energy minima; the implication of this is remarkable. We find that the conventional theories can predict the contact angle at the global free energy minimum if the droplet size is 40 times or larger than the characteristic scale of the surface roughness. Furthermore, we confirm that the pinning origin is the local free energy minima and obtain the energy barriers of pinning as a function of geometric factors. We validate our theory against experimental results on an anisotropic rough surface. In addition, we discuss wetting on a non-uniform rough surface with a rough central region and flat edge. Our findings clarify the extent to which the conventional wetting theory is valid and expand the physical understanding of wetting phenomena of small liquid drops on rough surfaces.

Preferential concentration of inertial sub-kolmogorov particles. The roles of mass loading of particles, Stokes and Reynolds numbers

Abstract: Turbulent flows laden with inertial particles present multiple open questions and are a subject of great interest in current research. Due to their higher density compared to the carrier fluid, inertial particles tend to form high concentration regions, i.e. clusters, and low concentration regions, i.e. voids, due to the interaction with the turbulence. In this work, we present an experimental investigation of the clustering phenomenon of heavy sub-Kolmogorov particles in homogeneous isotropic turbulent flows. Three control parameters have been varied over significant ranges: $Re_{\lambda} \in [170 - 450]$, $St\in [0.1 - 5]$ and volume fraction $\phi_v\in [2\times 10^{-6} - 2\times 10^{-5}]$. The scaling of clustering characteristics, such as the distribution of Voronoi areas and the dimensions of cluster and void regions, with the three parameters are discussed. In particular, for the polydispersed size distributions considered here, clustering is found to be enhanced strongly (quasi-linearly) by $Re_{\lambda}$ and noticeably (with a square-root dependency) with $\phi_v$, while the cluster and void sizes, scaled with the Kolmogorov lengthscale $\eta$, are driven primarily by $Re_{\lambda}$. Cluster length $\sqrt{\langle A_c \rangle}$ scales up to $\approx 100 {\eta}$, measured at the highest $Re_{\lambda}$, while void length $\sqrt{\langle A_v \rangle}$ scaled also with $\eta$ is typically two times larger ($\approx 200 {\eta}$). The lack of sensitivity of the above characteristics to the Stokes number lends support to the "sweep-stick" particle accumulation scenario. The non-negligible influence of the volume fraction, however, is not considered by that model and can be connected with collective effects.

On the Quantification of Mixing in Microfluidics

Abstract: Methods for quantifying mixing in microfluidics have varied largely in the past, and various indices have been employed to represent the extent of mixing. Mixing between two or more colored liquids is usually quantified using simple mathematical functions operated over a sequence of images. The function, usually termed mixing indices, involves a measure of standard deviation. Here, we first review some mixing indices and then experimentally verify the index most representative of a mixing event. It is observed that the relative mixing index is not affected by the lighting conditions, unlike other known mixing indices. Based on this finding, the use of a relative mixing index is advocated for further use in the lab-on-a-chip community for quantifying mixing events.

Global and local statistics in turbulent convection at low Prandtl numbers

Abstract: Statistical properties of turbulent Rayleigh-Benard convection at low Prandtl numbers (Pr), which are typical for liquid metals such as mercury, gallium or liquid sodium, are investigated in high-resolution three-dimensional spectral element simulations in a closed cylindrical cell with an aspect ratio of one and are compared to previous turbulent convection simulations in air. We compare the scaling of global momentum and heat transfer. The scaling exponents are found to be in agreement with experiments. Mean profiles of the root-mean-square velocity as well as the thermal and kinetic energy dissipation rates have growing amplitudes with decreasing Prandtl number which underlies a more vigorous bulk turbulence in the low-Pr regime. The skin-friction coefficient displays a Reynolds-number dependence that is close to that of an isothermal, intermittently turbulent velocity boundary layer. The thermal boundary layer thicknesses are larger as Pr decreases and conversely the velocity boundary layer thicknesses become smaller. We investigate the scaling exponents and find a slight decrease in exponent magnitude for the thermal boundary layer thickness as Pr decreases, but find the opposite case for the velocity boundary layer thickness scaling. A growing area fraction of turbulent patches close to the heating and cooling plates can be detected by exceeding a locally defined shear Reynolds number threshold. This area fraction is larger for lower Pr at the same Ra. Our analysis of the kurtosis of the locally defined shear Reynolds number demonstrates that the intermittency in the boundary layer is significantly increased for the lower Prandtl number and for sufficiently high Rayleigh number compared to convection in air. This complements our previous findings of enhanced bulk intermittency in low-Prandtl-number convection.

Effects of Pore-Scale Disorder on Fluid Displacement in Partially-Wettable Porous Media

Abstract: We present a systematic, quantitative assessment of the impact of pore size disorder and its interplay with flow rates and wettability on immiscible displacement of a viscous fluid. Pore-scale simulations and micromodel experiments show that reducing disorder increases the displacement efficiency and compactness, minimizing the fluid-fluid interfacial area, through (i) trapping at low rates and (ii) viscous fingering at high rates. Increasing the wetting angle suppresses both trapping and fingering, hence reducing the sensitivity of the displacement to the underlying disorder. A modified capillary number $\mathrm{Ca^*}$ that includes the impact of disorder $\lambda$ on viscous forces (through pore connectivity) is direct related to $\lambda$, in par with previous works. Our findings bear important consequences on sweep efficiency and fluid mixing and reactions, which are key in applications such as microfluidics to carbon geosequestration, energy recovery, and soil aeration and remediation.

A framework for computing effective boundary conditions at the interface between free fluid and a porous medium

Abstract: Interfacial boundary conditions determined from empirical or ad-hoc models remain the standard approach to model fluid flows over porous media, even in situations where the topology of the porous medium is known. We propose a non-empirical and accurate method to compute the effective boundary conditions at the interface between a porous surface and an overlying flow. Using multiscale expansion (homogenization) approach, we derive a tensorial generalized version of the empirical condition suggested by Beavers & Joseph (1967). The components of the tensors determining the effective slip velocity at the interface are obtained by solving a set of Stokes equations in a small computational domain near the interface containing both free flow and porous medium. Using the lid-driven cavity flow with a porous bed, we demonstrate that the derived boundary condition is accurate and robust by comparing an effective model to direct numerical simulations. Finally, we provide an open source code that solves the microscale problems and computes the velocity boundary condition without free parameters over any porous bed.

Effects of Substrate Heating and Wettability on Evaporation Dynamics and Deposition Patterns for a Sessile Water Droplet Containing Colloidal Particles

Abstract: Effects of substrate temperature, substrate wettability and particles concentration are experimentally investigated for evaporation of a sessile water droplet containing colloidal particles. Time-varying droplet shapes and temperature of the liquid-gas interface are measured using high-speed visualization and infrared thermography, respectively. The motion of the particles inside the evaporating droplet is qualitatively visualized by an optical microscope and profile of final particle deposit is measured by an optical profilometer. On a non-heated hydrophilic substrate, a ring-like deposit forms after the evaporation, as reported extensively in the literature; while on a heated hydrophilic substrate, a thinner ring with an inner deposit is reported in the present work. The latter is attributed to Marangoni convection and recorded motion of the particles as well as measured temperature gradient across the liquid-gas interface confirms this hypothesis. The thinning of the ring scales with the substrate temperature and is reasoned to stronger Marangoni convection at larger substrate temperature. In case of a non-heated hydrophobic substrate, an inner deposit forms due to very early depinning of the contact line. On the other hand, in case of a heated hydrophobic substrate, the substrate heating as well as larger particle concentration helps in the pinning of the contact line, which results in a thin ring with an inner deposit. We propose a regime map for predicting three types of deposits namely, ring, thin ring with inner deposit and inner deposit - for varying substrate temperature, substrate wettability and particles concentration. A first-order model corroborates the liquid-gas interface temperature measurements and variation in the measured ring profile with the substrate temperature.

Stability of stratified two-phase flows in horizontal channels

Abstract: Linear stability of stratified two-phase flows in horizontal channels to arbitrary wavenumber disturbances is studied. The problem is reduced to Orr-Sommerfeld equations for the stream function disturbances, defined in each sublayer and coupled via boundary conditions that account also for possible interface deformation and capillary forces. Applying the Chebyshev collocation method, the equations and interface boundary conditions are reduced to the generalized eigenvalue problems solved by standard means of numerical linear algebra for the entire spectrum of eigenvalues and the associated eigenvectors. Some additional conclusions concerning the instability nature are derived from the most unstable perturbation patterns. The results are summarized in the form of stability maps showing the operational conditions at which a stratified-smooth flow pattern is stable. It is found that for gas-liquid and liquid-liquid systems the stratified flow with a smooth interface is stable only in confined zone of relatively low flow rates, which is in agreement with experiments, but is not predicted by long-wave analysis. Depending on the flow conditions, the critical perturbations can originate mainly at the interface (so-called "interfacial modes of instability") or in the bulk of one of the phases (i.e., "shear modes"). The present analysis revealed that there is no definite correlation between the type of instability and the perturbation wavelength.

Drops bouncing off macro-textured superhydrophobic surfaces

Abstract: Recent experiments with droplets impacting a macro-textured superhydrophobic surfaces revealed new regimes of bouncing with a remarkable reduction of the contact time. We present here a comprehensive numerical study that reveals the physics behind these new bouncing regimes and quantify the role played by various external and internal forces that effect the dynamics of a drop impacting a complex surface. For the first time, three-dimensional simulations involving macro-textured surfaces are performed. Aside from demonstrating that simulations reproduce experiments in a quantitative manner, the study is focused on analyzing the flow situations beyond current experiments. We show that the experimentally observed reduction of contact time extends to higher Weber numbers, and analyze the role played by the texture density. Moreover, we report a non-linear behavior of the contact time with the increase of the Weber number for application relevant imperfectly coated textures, and also study the impact on tilted surfaces in a wide range of Weber numbers. Finally, we present novel energy analysis techniques that elaborate and quantify the interplay between the kinetic and surface energy, and the role played by the dissipation for various Weber numbers.

Searching for turbulence models by artificial neural network

Abstract: Artificial neural network (ANN) is tested as a tool for finding a new subgrid model of the subgrid-scale (SGS) stress in large-eddy simulation. ANN is used to establish a functional relation between the grid-scale (GS) flow field and the SGS stress without any assumption of the form of function. Data required for training and test of ANN are provided by direct numerical simulation (DNS) of a turbulent channel flow. It is shown that ANN can establish a model similar to the gradient model. The correlation coefficients between the real SGS stress and the output of ANN are comparable to or larger than similarity models, but smaller than a two-parameter dynamic mixed model.

Properties of the turbulent/non-turbulent interface in boundary layers

Abstract: The turbulent/non-turbulent interface is analysed in a direct numerical simulation of a boundary layer in the range $Re_\theta=2800-6600$, with emphasis on the behaviour of the relatively large-scale fractal intermittent region. This requires the introduction of a new definition of the distance between a point and a general surface, which is compared with the more usual vertical distance to the top of the layer. Interfaces are obtained by thresholding the enstrophy field and the magnitude of the rate-of-strain tensor, and it is concluded that, while the former are physically relevant features, the latter are not. By varying the threshold, a topological transition is identified as the interface moves from the free stream into the turbulent core. A vorticity scale is defined that collapses that transition for different Reynolds numbers, roughly equivalent to the root-mean-squared vorticity at the edge of the boundary layer. Conditionally averaged flow variables are analysed as functions of the new distance, both within and outside the interface. It is found that the interface contains a nonequilibrium layer whose thickness scales well with the Taylor microscale, enveloping a self-similar layer spanning a fixed fraction of the boundary-layer thickness. Interestingly, the straining structure of the flow is similar in both regions. Irrotational pockets within the turbulent core are also studied. They form a self-similar set whose size decreases with increasing depth, presumably due to break-up by the turbulence, but the rate of viscous diffusion is independent of the pocket size. The raw data used in the analysis are freely available from our web page (this http URL).

Coherent structures and extreme events in rotating multiphase turbulent flows

Abstract: By using direct numerical simulations (DNS) at unprecedented resolution we study turbulence under rotation in the presence of simultaneous direct and inverse cascades. The accumulation of energy at large scale leads to the formation of vertical coherent regions with high vorticity oriented along the rotation axis. By seeding the flow with millions of inertial particles, we quantify -for the first time- the effects of those coherent vertical structures on the preferential concentration of light and heavy particles. Furthermore, we quantitatively show that extreme fluctuations, leading to deviations from a normal-distributed statistics, result from the entangled interaction of the vertical structures with the turbulent background. Finally, we present the first-ever measurement of the relative importance between Stokes drag, Coriolis force and centripetal forces along the trajectories of inertial particles. We discover that vortical coherent structures lead to unexpected diffusion properties for heavy and light particles in the directions parallel and perpendicular to the rotation axis.

Synchronised Swimming of Two Fish

Abstract: We study the fluid dynamics of two fish-like bodies with synchronised swimming patterns. Our studies are based on two-dimensional simulations of viscous incompressible flows. We distinguish between motion patterns that are externally imposed on the swimmers and self-propelled swimmers that learn manoeuvres to achieve certain goals. Simulations of two rigid bodies executing pre-specified motion indicate that flow-mediated interactions can lead to substantial drag reduction and may even generate thrust intermittently. In turn we examine two self-propelled swimmers arranged in a leader-follower configuration, with a-priori specified body-deformations. We find that the swimming of the leader remains largely unaffected, while the follower experiences either an increase or decrease in swimming speed, depending on the initial conditions. Finally, we consider a follower that synchronises its motion so as to minimise its lateral deviations from the leader's path. The leader employs a steady gait while the follower uses a reinforcement learning algorithm to adapt its swimming-kinematics. We find that swimming in a synchronised tandem can yield up to about 30% reduction in energy expenditure for the follower, in addition to a 20% increase in its swimming-efficiency. The present results indicate that synchronised swimming of two fish can be energetically beneficial.

A nonlinear model for rotationally constrained convection with Ekman pumping

Abstract: It is a well established result of linear theory that the influence of differing mechanical boundary conditions, i.e., stress-free or no-slip, on the primary instability in rotating convection becomes asymptotically small in the limit of rapid rotation. This is accounted for by the diminishing impact of the viscous stresses exerted within Ekman boundary layers and the associated vertical momentum transport by Ekman pumping. By contrast, in the nonlinear regime recent experiments and supporting simulations are now providing evidence that the efficiency of heat transport remains strongly influenced by Ekman pumping in the rapidly rotating limit. In this paper, a reduced model is developed for the case of low Rossby number convection in a plane layer geometry with no-slip upper and lower boundaries held at fixed temperatures. A complete description of the dynamics requires the existence of three distinct regions within the fluid layer: a geostrophically balanced interior where fluid motions are predominately aligned with the axis of rotation, Ekman boundary layers immediately adjacent to the bounding plates, and thermal wind layers driven by Ekman pumping in between. The reduced model uses a classical Ekman pumping parameterization to alleviate the need for spatially resolving the Ekman boundary layers. Results are presented for both linear stability theory and a special class of nonlinear solutions described by a single horizontal spatial wavenumber. It is shown that Ekman pumping allows for significant enhancement in the heat transport relative to that observed in simulations with stress-free boundaries. Without the intermediate thermal wind layer the nonlinear feedback from Ekman pumping would be able to generate a heat transport that diverges to infinity. This layer arrests this blowup resulting in finite heat transport at a significantly enhanced value.

Geometric control of active collective motion

Abstract: Recent experimental studies have shown that confinement can profoundly affect self-organization in semi-dilute active suspensions, leading to striking features such as the formation of steady and spontaneous vortices in circular domains and the emergence of unidirectional pumping motions in periodic racetrack geometries. Motivated by these findings, we analyze the two-dimensional dynamics in confined suspensions of active self-propelled swimmers using a mean-field kinetic theory where conservation equations for the particle configurations are coupled to the forced Navier-Stokes equations for the self-generated fluid flow. In circular domains, a systematic exploration of the parameter space casts light on three distinct states: equilibrium with no flow, stable vortex, and chaotic motion, and the transitions between these are explained and predicted quantitatively using a linearized theory. In periodic racetracks, similar transitions from equilibrium to net pumping to traveling waves to chaos are observed in agreement with experimental observations and are also explained theoretically. Our results underscore the subtle effects of geometry on the morphology and dynamics of emerging patterns in active suspensions and pave the way for the control of active collective motion in microfluidic devices.

Role of natural convection in the dissolution of sessile droplets

Abstract: The dissolution process of small (initial (equivalent) radius $R_0 Ra_t$, where $Ra_t = 12$ is the transition Ra-number as extracted from the data. For $Ra < Ra_t$ and smaller, convective transport is progressively overtaken by diffusion and the above scaling relations break down.

Actuator line modeling of vertical-axis turbines

Abstract: To bridge the gap between high and low fidelity numerical modeling tools for vertical-axis (or cross-flow) turbines (VATs or CFTs), an actuator line model (ALM) was developed and validated for both a high and a medium solidity vertical-axis turbine at rotor diameter Reynolds numbers $Re_D \sim 10^6$. The ALM is a combination of classical blade element theory and Navier--Stokes based flow models, and in this study both $k$--$\epsilon$ Reynolds-averaged Navier--Stokes (RANS) and Smagorinsky large eddy simulation (LES) turbulence models were tested using the open-source OpenFOAM computational fluid dynamics framework. The RANS models were able to be run on coarse grids while still providing good convergence behavior in terms of the mean power coefficient, and also approximately four orders of magnitude reduction in computational expense compared with 3-D blade-resolved RANS simulations. Submodels for dynamic stall, end effects, added mass, and flow curvature were implemented, resulting in reasonable performance predictions for the high solidity rotor, more discrepancies for the medium solidity rotor, and overprediction for both cases at high tip speed ratio. The wake results showed that the ALM was able to capture some of the important flow features that contribute to VAT's relatively fast wake recovery---a large improvement over the conventional actuator disk model. The mean flow field was better realized with the LES, which still represented a computational savings of two orders of magnitude compared with 3-D blade-resolved RANS, though vortex breakdown and subsequent turbulence generation appeared to be underpredicted, which necessitates further investigation of optimal subgrid scale modeling.

Puffing and micro-explosion behavior in combustion of butanol/Jet A-1 and acetone-butanol-ethanol (A-B-E)/Jet A-1 fuel droplets

Abstract: The present investigation deals with the puffing and micro-explosion characteristics in the combustion of a single droplet comprising butanol/Jet A-1, acetone-butanol-ethanol (A-B-E)/Jet A-1 blends, and A-B-E. The onset of nucleation, growth of vapor bubble and subsequent breakup of droplet for various fuel blends have been analyzed from the high-speed images. Puffing was observed to be the dominant phenomenon in 30% butanol blend, while micro-explosion was found to be the dominant one in other fuel blends (blend with 50% butanol or 30% A-B-E or 50% A-B-E). It was observed that puffing always preceded the micro-explosion. The probability of micro-explosion in droplets with A-B-E blends was found to be higher than that of butanol blends. Although the rate of bubble growth was almost similar for all butanol and A-B-E blends, the final bubble diameter before the droplet breakup was found to be higher for 50/50 blends than that of 30/70 blends. The occurrence of micro-explosion shortened the droplet lifetime, and this effect appeared to be stronger for droplets with 50/50 composition. Micro-explosion led to the ejection of both larger and smaller secondary droplets; however, puffing resulted in relatively smaller secondary droplets compared to micro-explosion. Puffing/micro-explosion were also observed in the secondary droplets.

A nonlinear small-deformation theory for transient droplet electrohydrodynamics

Abstract: The deformation of a viscous liquid droplet suspended in another liquid and subject to an applied electric field is a classic multiphase flow problem best described by the Melcher-Taylor leaky dielectric model The main assumption of the model is that any net charge in the system is concentrated on the interface between the two liquids as a result of the jump in Ohmic currents from the bulk Upon application of the field, the drop can either attain a steady prolate or oblate shape with toroidal circulating flows both inside and outside arising from tangential stresses on the interface due to action of the field on the surface charge distribution Since the pioneering work of \cite{taylor1966}, there have been numerous computational and theoretical studies to predict the deformations measured in experiments Most existing theoretical models, however, have either neglected transient charge relaxation or nonlinear charge convection by the interfacial flow In this work, we develop a novel small-deformation theory accurate to second order in electric capillary number ${O}(Ca_E^2)$ for the complete Melcher-Taylor model that includes transient charge relaxation, charge convection by the flow, as well as transient shape deformation The main result of the paper is the derivation of coupled evolution equations for the induced electric multipoles and for the shape functions describing the deformations on the basis of spherical harmonics Our results, which are consistent with previous models in the appropriate limits, show excellent agreement with fully nonlinear numerical simulations based on an axisymmetric boundary-element formulation and with existing experimental data in the small-deformation regime