Showing papers in "arXiv: Fluid Dynamics in 2013"

The unifying theory of scaling in thermal convection: The updated prefactors

Abstract: The unifying theory of scaling in thermal convection (Grossmann & Lohse (2000)) (henceforth the GL theory) suggests that there are no pure power laws for the Nusselt and Reynolds numbers as function of the Rayleigh and Prandtl numbers in the experimentally accessible parameter regime. In Grossmann & Lohse (2001) the dimensionless parameters of the theory were fitted to 155 experimental data points by Ahlers & Xu (2001) in the regime $3\times 10^7 \le Ra \le 3 \times 10^{9}$ and $4\le Pr \le 34$ and Grossmann & Lohse (2002) used the experimental data point from Qiu & Tong (2001) and the fact that Nu(Ra,Pr) is independent of the parameter a, which relates the dimensionless kinetic boundary thickness with the square root of the wind Reynolds number, to fix the Reynolds number dependence. Meanwhile the theory is on one hand well confirmed through various new experiments and numerical simulations. On the other hand these new data points provide the basis for an updated fit in a much larger parameter space. Here we pick four well established (and sufficiently distant) Nu(Ra,Pr) data points and show that the resulting Nu(Ra,Pr) function is in agreement with almost all established experimental and numerical data up to the ultimate regime of thermal convection, whose onset also follows from the theory. One extra Re(Ra,Pr) data point is used to fix Re(Ra,Pr). As Re can depend on the definition and the aspect ratio the transformation properties of the GL equations are discussed in order to show how the GL coefficients can easily be adapted to new Reynolds number data while keeping Nu(Ra,Pr) unchanged.

On coherent structure in wall turbulence

Abstract: A new theory of coherent structure in wall turbulence is presented. The theory is the first to predict packets of hairpin vortices and other structure in turbulence, and their dynamics, based on an analysis of the Navier-Stokes equations, under an assumption of a turbulent mean profile. The assumption of the turbulent mean acts as a restriction on the class of possible structures. It is shown that the coherent structure is a manifestation of essentially low-dimensional flow dynamics, arising from a critical layer mechanism. Using the decomposition presented in McKeon & Sharma (J. Fluid Mech, 658, 2010), complex coherent structure is recreated from minimal superpositions of response modes predicted by the analysis, which take the form of radially-varying travelling waves. By way of example, simple combinations of these modes are offered that predicts hairpins and modulated hairpin packets. The phase interaction also predicts important skewness and correlation results known in the literature. It is also shown that the very large scale motions act to organise hairpin-like structures such that they co-locate with areas of low streamwise momentum, by a mechanism of locally varying the shear profile. The relationship between Taylor's hypothesis and coherence is discussed and both are shown to be the consequence of the localisation of the response modes around the critical layer. A pleasing link is made to the classical laminar inviscid theory, whereby the essential mechanism underlying the hairpin vortex is captured by two obliquely interacting Kelvin-Stuart (cat's eye) vortices. Evidence for the theory is presented based on comparison to observations of structure reported in the experimental, transitional flow and turbulent flow numerical simulation literature.

Sparsity-promoting dynamic mode decomposition

Abstract: Dynamic mode decomposition (DMD) represents an effective means for capturing the essential features of numerically or experimentally generated flow fields. In order to achieve a desirable tradeoff between the quality of approximation and the number of modes that are used to approximate the given fields, we develop a sparsity-promoting variant of the standard DMD algorithm. In our method, sparsity is induced by regularizing the least-squares deviation between the matrix of snapshots and the linear combination of DMD modes with an additional term that penalizes the $\ell_1$-norm of the vector of DMD amplitudes. The globally optimal solution of the resulting regularized convex optimization problem is computed using the alternating direction method of multipliers, an algorithm well-suited for large problems. Several examples of flow fields resulting from numerical simulations and physical experiments are used to illustrate the effectiveness of the developed method.

Visualizing multiphase flow and trapped fluid configurations in a model three-dimensional porous medium

Abstract: We report an approach to fully visualize the flow of two immiscible fluids through a model three-dimensional (3D) porous medium at pore-scale resolution. Using confocal microscopy, we directly image the drainage of the medium by the non-wetting oil and subsequent imbibition by the wetting fluid. During imbibition, the wetting fluid pinches off threads of oil in the narrow crevices of the medium, forming disconnected oil ganglia. Some of these ganglia remain trapped within the medium. By resolving the full 3D structure of the trapped ganglia, we show that the typical ganglion size, and the total amount of residual oil, decreases as the capillary number Ca increases; this behavior reflects the competition between the viscous pressure in the wetting fluid and the capillary pressure required to force oil through the pores of the medium. This work thus shows how pore-scale fluid dynamics influence the trapped fluid configurations in multiphase flow through 3D porous media.

Condensation and jumping relay of droplets on lotus leaf

Abstract: Dynamic behavior of micro water droplet condensed on a lotus leaf with two-tier roughness is studied. Under laboratory environment, the contact angle of the micro droplet on single micro papilla increases smoothly from 80 deg to 160 deg during the growth of condensed water. The best-known "self-clean" phenomenon, will be lost. A striking observation is the out-of-plane jumping relay of condensed droplets triggered by falling droplets, as well as its sustained speed obtained in continuous jumping relays, enhance the automatic removal of dropwise condensation without the help from any external force. The surface tension energy dissipation is the main reason controlling the critical size of jumping droplet and its onset velocity of rebounding.

Fluid flow dynamics under location uncertainty

Abstract: We present a derivation of a stochastic model of Navier Stokes equations that relies on a decomposition of the velocity fields into a differentiable drift component and a time uncorrelated uncertainty random term. This type of decomposition is reminiscent in spirit to the classical Reynolds decomposition. However, the random velocity fluctuations considered here are not differentiable with respect to time, and they must be handled through stochastic calculus. The dynamics associated with the differentiable drift component is derived from a stochastic version of the Reynolds transport theorem. It includes in its general form an uncertainty dependent "subgrid" bulk formula that cannot be immediately related to the usual Boussinesq eddy viscosity assumption constructed from thermal molecular agitation analogy. This formulation, emerging from uncertainties on the fluid parcels location, explains with another viewpoint some subgrid eddy diffusion models currently used in computational fluid dynamics or in geophysical sciences and paves the way for new large-scales flow modelling. We finally describe an applications of our formalism to the derivation of stochastic versions of the Shallow water equations or to the definition of reduced order dynamical systems.

Model-based scaling of the streamwise energy density in high-Reynolds number turbulent channels

Abstract: We study the Reynolds number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (2010), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier-Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behavior with Reynolds number on the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it was shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. We identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_\tau \approx 10^3 - 10^{10}$).

DNS of horizontal open channel flow with finite-size, heavy particles at low solid volume fraction

Abstract: We have performed direct numerical simulation of turbulent open channel flow over a smooth horizontal wall in the presence of finite-size, heavy particles. The spherical particles have a diameter of approximately 7 wall units, a density of 1.7 times the fluid density and a solid volume fraction of 0.0005. The value of the Galileo number is set to 16.5, while the Shields parameter measures approximately 0.2. Under these conditions, the particles are predominantly located in the vicinity of the bottom wall, where they exhibit strong preferential concentration which we quantify by means of Voronoi analysis and by computing the particle-conditioned concentration field. As observed in previous studies with similar parameter values, the mean streamwise particle velocity is smaller than that of the fluid. We propose a new definition of the fluid velocity "seen" by finite-size particles based on an average over a spherical surface segment, from which we deduce in the present case that the particles are instantaneously lagging the fluid only by a small amount. The particle-conditioned fluid velocity field shows that the particles preferentially reside in the low-speed streaks, leading to the observed apparent lag. Finally, a vortex eduction study reveals that spanwise particle motion is significantly correlated with the presence of vortices with the corresponding sense of rotation which are located in the immediate vicinity of the near-wall particles.

Experimental study of parametric subharmonic instability for internal waves

Abstract: Internal waves are believed to be of primary importance as they affect ocean mixing and energy transport. Several processes can lead to the breaking of internal waves and they usually involve non linear interactions between waves. In this work, we study experimentally the parametric subharmonic instability (PSI), which provides an efficient mechanism to transfer energy from large to smaller scales. It corresponds to the destabilization of a primary plane wave and the spontaneous emission of two secondary waves, of lower frequencies and different wave vectors. Using a time-frequency analysis, we observe the time evolution of the secondary waves, thus measuring the growth rate of the instability. In addition, a Hilbert transform method allows the measurement of the different wave vectors. We compare these measurements with theoretical predictions, and study the dependence of the instability with primary wave frequency and amplitude, revealing a possible effect of the confinement due to the finite size of the beam, on the selection of the unstable mode.

Estimating Uncertainties in Statistics Computed from DNS

Abstract: Rigorous assessment of uncertainty is crucial to the utility of DNS results. Uncertainties in the computed statistics arise from two sources: finite statistical sampling and the discretization of the Navier-Stokes equations. Due to the presence of non-trivial sampling error, standard techniques for estimating discretization error (such as Richardson extrapolation) fail or are unreliable. This work provides a systematic and unified approach for estimating these errors. First, a sampling error estimator that accounts for correlation in the input data is developed. Then, this sampling error estimate is used as part of a Bayesian extension of Richardson extrapolation in order to characterize the discretization error. These methods are tested using the Lorenz equations and are shown to perform well. These techniques are then used to investigate the sampling and discretization errors in the DNS of a wall-bounded turbulent flow. For both cases, it is found that while the sampling uncertainty is large enough to make the order of accuracy difficult to determine, the estimated discretization errors are quite small. This indicates that the commonly used heuristics provide ad- equate resolution for this class of problems. However, it is also found that, for some quantities, the discretization error is not small relative to sampling error, indicating that the conventional wisdom that sampling error dominates discretization error for this class of simulations needs to be reevaluated.

Low-Reynolds number swimming in a capillary tube

Abstract: We use the boundary element method to study the low-Reynolds number locomotion of a spherical model microorganism in a circular tube. The swimmer propels itself by tangen- tial or normal surface motion in a tube whose radius is on the order of the swimmer size. Hydrodynamic interactions with the tube walls significantly affect the average swimming speed and power consumption of the model microorganism. In the case of swimming parallel to the tube axis, the locomotion speed is always reduced (resp. increased) for swimmers with tangential (resp. normal) deformation. In all cases, the rate of work nec- essary for swimming is increased by confinement. Swimmers with no force-dipoles in the far field generally follow helical trajectories, solely induced by hydrodynamic interactions with the tube walls, and in qualitative agreement with recent experimental observations for Paramecium. Swimmers of the puller type always display stable locomotion at a lo- cation which depends on the strength of their force dipoles: swimmers with weak dipoles (small {\alpha}) swim in the centre of the tube while those with strong dipoles (large {\alpha}) swim near the walls. In contrast, pusher swimmers and those employing normal deformation are unstable and end up crashing into the walls of the tube. Similar dynamics is observed for swimming into a curved tube. These results could be relevant for the future design of artificial microswimmers in confined geometries.

Kelvin wake pattern at large Froude numbers

Abstract: Gravity waves generated by an object moving at constant speed at the water surface form a specific pattern commonly known as the Kelvin wake. It was proved by Lord Kelvin that such a wake is delimited by a constant angle $\simeq 19.47^\circ$. However a recent study by Rabaud and Moisy based on the observation of airborne images showed that the wake angle seems to decrease as the Froude number $Fr$ increases, scaling as $Fr^{-1}$ for large Froude numbers. To explain such observations the authors make the strong hypothesis that an object of size $b$ cannot generate wavelengths larger than $b$. With no need of such an assumption and modelling the moving object by an axisymmetric pressure field, we analytically show that the angle corresponding to the maximum amplitude of the waves scales as $Fr^{-1}$ for large Froude numbers, whereas the angle delimiting the wake region outside which the surface is essentially flat remains constant and equal to the Kelvin angle for all $Fr$.

Generalization of the Rotne-Prager-Yamakawa mobility and shear disturbance tensors

Abstract: Rotne-Prager-Yamakawa approximation is one of the most commonly used methods of including hydrodynamic interactions in modelling of colloidal suspensions and polymer solutions. The two main merits of this approximation is that it includes all long-range terms (i.e. decaying as R^-3 or slower in interparticle distances) and that the diffusion matrix is positive definite, which is essential for Brownian dynamics modelling. Here, we extend the Rotne-Prager-Yamakawa approach to include both translational and rotational degrees of freedom, and derive the regularizing corrections to account for overlapping particles. Additionally, we show how the Rotne-Prager-Yamakawa approximation can be generalized for other geometries and boundary conditions.

Proper orthogonal decomposition closure models for fluid flows: Burgers equation

Abstract: This paper puts forth several closure models for the proper orthogonal decomposition (POD) reduced order modeling of fluid flows These new closure models, together with other standard closure models, are investigated in the numerical simulation of the Burgers equation This simplified setting represents just the first step in the investigation of the new closure models It allows a thorough assessment of the performance of the new models, including a parameter sensitivity study Two challenging test problems displaying moving shock waves are chosen in the numerical investigation The closure models and a standard Galerkin POD reduced order model are benchmarked against the fine resolution numerical simulation Both numerical accuracy and computational efficiency are used to assess the performance of the models

Drops on soft solids: Free energy and double transition of contact angles

Abstract: The equilibrium shape of liquid drops on elastic substrates is determined by minimising elastic and capillary free energies, focusing on thick incompressible substrates. The problem is governed by three length scales: the size of the drop $R$, the molecular size $a$, and the ratio of surface tension to elastic modulus $\gamma/E$. We show that the contact angles undergo two transitions upon changing the substrates from rigid to soft. The microscopic wetting angles deviate from Young's law when $\gamma/Ea \gg 1$, while the apparent macroscopic angle only changes in the very soft limit $\gamma/ER \gg 1$. The elastic deformations are worked out in the simplifying case where the solid surface energy is assumed constant. The total free energy turns out lower on softer substrates, consistent with recent experiments.

Direct Numerical Simulation of Electrokinetic Instability and Transition to Chaotic Motion

Abstract: A new type of instability - electrokinetic instability - and an unusual transition to chaotic motion near a charge-selective surface was studied by numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. Two kinds of initial conditions were considered: (a) white noise initial conditions to mimic "room disturbances" and subsequent natural evolution of the solution; (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution; two-dimensional steady electroconvective vortices (stationary point in a proper phase space); unsteady vortices aperiodically changing their parameters (homoclinic contour); periodic motion (limit cycle); and chaotic motion. The transition to chaotic motion did not include Hopf bifurcation. Numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt-current characteristics.

Alignment of vorticity and rods with Lagrangian fluid stretching in turbulence

Abstract: Stretching in continuum mechanics is naturally described using the Cauchy-Green strain tensors. These tensors quantify the Lagrangian stretching experienced by a material element, and provide a powerful way to study processes in turbulent fluid flows that involve stretching such as vortex stretching and alignment of anisotropic particles. Analyzing data from a simulation of isotropic turbulence, we observe preferential alignment between anisotropic particles and vorticity. We show that this alignment arises because both of these quantities independently tend to align with the strongest Lagrangian stretching direction, as defined by the maximum eigenvector of the left Cauchy-Green strain tensor. In particular, anisotropic particles approach almost perfect alignment with the strongest stretching direction. The alignment of vorticity with stretching is weaker, but still much stronger than previously observed alignment of vorticity with the eigenvectors of the Eulerian strain rate tensor. The alignment of strong vorticity is almost the same as that of rods that have experienced the same stretching.

Boundary layer dynamics at the transition between the classical and the ultimate regime of Taylor-Couette flow

Abstract: Direct numerical simulations of turbulent Taylor-Couette flow are performed up to inner cylinder Reynolds numbers of {Re_i=10^5} for a radius ratio of {\eta=r_i/r_o=0.714} between the inner and outer cylinder. With increasing {Re_i}, the flow undergoes transitions between three different regimes: (i) a flow dominated by large coherent structures, (ii) an intermediate transitional regime, and (iii) a flow with developed turbulence. In the first regime the large--scale rolls completely drive the meridional flow while in the second one the coherent structures recover only on average. The presence of a mean flow allows for the coexistence of laminar and turbulent boundary layer dynamics. In the third regime the mean flow effects fade away and the flow becomes dominated by plumes. The effect of the local driving on the azimuthal and angular velocity profiles is quantified, in particular we show when and where those profiles develop.

Physics of Rheologically-Enhanced Propulsion: Different Strokes in Generalized Stokes

Abstract: Shear-thinning is an important rheological property of many biological fluids, such as mucus, whereby the apparent viscosity of the fluid decreases with shear. Certain microscopic swimmers have been shown to progress more rapidly through shear-thinning fluids, but is this behavior generic to all microscopic swimmers, and what are the physics through which shear-thinning rheology affects a swimmer's propulsion? We examine swimmers employing prescribed stroke kinematics in two-dimensional, inertialess Carreau fluid: shear-thinning "Generalized Stokes" flow. Swimmers are modeled, using the method of femlets, by a set of immersed, regularized forces. The equations governing the fluid dynamics are then discretized over a body-fitted mesh and solved with the finite element method. We analyze the locomotion of three distinct classes of microswimmer: (1) conceptual swimmers comprising sliding spheres employing both one- and two-dimensional strokes, (2) slip-velocity envelope models of ciliates commonly referred to as "squirmers" and (3) monoflagellate pushers, such as sperm. We find that morphologically identical swimmers with different strokes may swim either faster or slower in shear-thinning fluids than in Newtonian fluids. We explain this kinematic sensitivity by considering differences in the viscosity of the fluid surrounding propulsive and payload elements of the swimmer, and using this insight suggest two reciprocal sliding sphere swimmers which violate Purcell's Scallop theorem in shear-thinning fluids. We also show that an increased flow decay rate arising from shear-thinning rheology is associated with a reduction in the swimming speed of slip-velocity squirmers. For sperm-like swimmers, a gradient of thick to thin fluid along the flagellum alters the force it exerts upon the fluid, flattening trajectories and increasing instantaneous swimming speed.

Spatially localized solutions of shear flows

Abstract: We present several new spatially localized equilibrium and traveling-wave solutions of plane Couette and channel flows. The solutions exhibit strikingly concentrated regions of vorticity that are flanked on either side by high-speed streaks. For several traveling-wave solutions of channel flow, the concentrated vortex structures are confined to the near-wall region and form particularly isolated and elemental coherent structures in the near-wall region of shear flows. The solutions are constructed by a variety of methods: application of windowing functions to previously known spatially periodic solutions, continuation from plane Couette to channel flow conditions, and from initial guesses obtained from turbulent simulation data. We show how the symmetries of localized solutions derive from the symmetries of their periodic counterparts, analyze the exponential decay of their tails, examine the scale separation and scaling of their streamwise Fourier modes, and show that they develop critical layers for large Reynolds numbers.

Mixing and transport by ciliary carpets: a numerical study

Abstract: We use a 3D computational model to study the fluid transport and mixing due to the beating of an infinite array of cilia. In accord with recent experiments, we observe two distinct regions: a fluid transport region above the cilia and a fluid mixing region below the cilia tip. The metachronal wave due to phase differences between neighboring cilia is known to enhance the fluid transport above the ciliary tip. In this work, we show that the metachronal wave also enhances the mixing rates in the sub-ciliary region, often simultaneously with the flow rate enhancement. Our results suggest that this simultaneous enhancement in transport and mixing is due to an enhancement in shear flow. As the flow above the cilia increases, shear rate in the fluid increases and such shear enhances stretching, which is an essential ingredient for mixing. Estimates of the mixing time scale indicate that, compared to diffusion, the mixing due to the cilia beat may be significant and sometimes dominates chemical diffusion.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

Abstract: We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates.

Slipping motion of large neutrally-buoyant particles in turbulence

Abstract: Direct numerical simulations are used to investigate the individual dynamics of large spherical particles suspended in a developed homogeneous turbulent flow. A definition of the direction of the particle motion relative to the surrounding flow is introduced and used to construct the mean fluid velocity profile around the particle. This leads to an estimate of the particle slipping velocity and its associated Reynolds number. The flow modifications due to the particle are then studied. The particle is responsible for a shadowing effect that occurs in the wake up to distances of the order of its diameter: the particle pacifies turbulent fluctuations and reduces the energy dissipation rate compared to its average value in the bulk. Dimensional arguments are presented to draw an analogy between particle effects on turbulence and wall flows. Evidence is obtained on the presence of a logarithmic sublayer at distances between the thickness of the viscous boundary layer and the particle diameter $\Dp$. Finally, asymptotic arguments are used to relate the viscous sublayer quantities to the particle size and the properties of the outer turbulence. It is shown in particular that the skin-friction Reynolds number behaves as $Re_\tau \propto (\Dp/\eta)^{4/3}$.

Leidenfrost temperature increase for impacting droplets on carbon-nanofiber surfaces

Abstract: Droplets impacting on a superheated surface can either exhibit a contact boiling regime, in which they make direct contact with the surface and boil violently, or a film boiling regime, in which they remain separated from the surface by their own vapor. The transition from the contact to the film boiling regime depends not only on the temperature of the surface and kinetic energy of the droplet, but also on the size of the structures fabricated on the surface. Here we experimentally show that surfaces covered with carbon-nanofibers delay the transition to film boiling to much higher temperature compared to smooth surfaces. We present physical arguments showing that, because of the small scale of the carbon fibers, they are cooled by the vapor flow just before the liquid impact, thus permitting contact boiling up to much higher temperatures than on smooth surfaces. We also show that, as long as the impact is in the film boiling regime, the spreading factor of impacting droplets follows the same $\We^{3/10}$ scaling (with $\We$ the Weber number) found for smooth surfaces, which is caused by the vapor flow underneath the droplet.

Orientation dynamics of small, triaxial-ellipsoidal particles in isotropic turbulence

Abstract: The orientation dynamics of small anisotropic tracer particles in turbulent flows is studied using direct numerical simulation (DNS) and results are compared with Lagrangian stochastic models. Generalizing earlier analysis for axisymmetric ellipsoidal particles (Parsa et al. 2012), we measure the orientation statistics and rotation rates of general, triaxial ellipsoidal tracer particles using Lagrangian tracking in DNS of isotropic turbulence. Triaxial ellipsoids that are very long in one direction, very thin in another, and of intermediate size in the third direction exhibit reduced rotation rates that are similar to those of rods in the ellipsoid's longest direction, while exhibiting increased rotation rates that are similar to those of axisymmetric discs in the thinnest direction. DNS results differ significantly from the case when the particle orientations are assumed to be statistically independent from the velocity gradient tensor. They are also different from predictions of a Gaussian process for the velocity gradient tensor, which does not provide realistic preferred vorticity-strain-rate tensor alignments. DNS results are also compared with a stochastic model for the velocity gradient tensor based on the recent fluid deformation approximation (RFDA). Unlike the Gaussian model, the stochastic model accurately predicts the reduction in rotation rate in the longest direction of triaxial ellipsoids since this direction aligns with the flow's vorticity, with its rotation perpendicular to the vorticity being reduced. For disc-like particles, or in directions perpendicular to the longest direction in triaxial particles, the model predicts {noticeably} smaller rotation rates than those observed in DNS, a behavior that can be understood based on the probability of vorticity orientation with the most contracting strain-rate eigen-direction in the model.

Energy spectra in turbulent bubbly flows

Abstract: We conduct experiments in a turbulent bubbly flow to study the nature of the transition between the classical $-$5/3 energy spectrum scaling for a single-phase turbulent flow and the $-$3 scaling for a swarm of bubbles rising in a quiescent liquid and of bubble-dominated turbulence. The bubblance parameter, which measures the ratio of the bubble-induced kinetic energy to the kinetic energy induced by the turbulent liquid fluctuations before bubble injection, is often used to characterise the bubbly flow. We vary the bubblance parameter from $b = \infty$ (pseudo-turbulence) to $b = 0$ (single-phase flow) over 2-3 orders of magnitude ($0.01 - 5$) to study its effect on the turbulent energy spectrum and liquid velocity fluctuations. The probability density functions (PDFs) of the liquid velocity fluctuations show deviations from the Gaussian profile for $b > 0$, i.e. when bubbles are present in the system. The PDFs are asymmetric with higher probability in the positive tails. The energy spectra are found to follow the $-$3 scaling at length scales smaller than the size of the bubbles for bubbly flows. This $-$3 spectrum scaling holds not only in the well-established case of pseudo-turbulence, but surprisingly in all cases where bubbles are present in the system ($b > 0$). Therefore, it is a generic feature of turbulent bubbly flows, and the bubblance parameter is probably not a suitable parameter to characterise the energy spectrum in bubbly turbulent flows. The physical reason is that the energy input by the bubbles passes over only to higher wave numbers, and the energy production due to the bubbles can be directly balanced by the viscous dissipation in the bubble wakes as suggested by Lance $\&$ Bataille (1991). In addition, we provide an alternative explanation by balancing the energy production of the bubbles with viscous dissipation in the Fourier space.

Current-driven flow instabilities in large-scale liquid metal batteries, and how to tame them

Abstract: The use of liquid metal batteries is considered as one promising option for electric grid stabilisation. While large versions of such batteries are preferred in view of the economies of scale, they are susceptible to various magnetohydrodynamic instabilities which imply a risk of short-circuiting the battery due to the triggered fluid flow. Here we focus on the current driven Tayler instability and give critical electrical currents for its onset as well as numerical estimates for the appearing flow structures and speeds. Scaling laws for different materials, battery sizes and geometries are found. We further discuss and compare various means for preventing the instability.

Mechanical energy dissipation induced by sloshing and wave breaking in a fully coupled angular motion system. Part I: Theoretical formulation and Numerical Investigation

Abstract: A dynamical system involving a driven pendulum filled with liquid, is analyzed in the present paper series. The study of such a system is conducted in order to understand energy dissipation resulting from the shallow water sloshing and induced wave breaking. This analysis is relevant for the design of Tuned Liquid Damper devices. The complexity and violence of the flow generated by the roll motion results in the impossibility of using an analytical approach, requiring in turn the use of a suitable numerical solver. In Part I, the coupled dynamical system is thoroughly described, revealing its nonlinear features associated with the large amplitude of the forcing, both in terms of mechanical and fluid dynamical aspects. A smoothed particle hydrodynamics (SPH) model, largely validated in literature, is used to calculate the frequency behavior of the whole system. For small rotation angles, a semi-analytical model of the energy dissipated by the fluid, based on a hydraulic jump solution, is developed; the energy transfer is numerically calculated in order to extend the analysis to large oscillation angles. The experimental part of the investigation is carried out in Part II of this work.

Counter-gradient heat transport in two-dimensional turbulent Rayleigh-B\'{e}nard convection

Abstract: We present high-resolution numerical investigations of heat transport by two-dimensional (2D) turbulent Rayleigh-B\'{e}nard (RB) convection over the Rayleigh number range $10^8 \leqslant Ra\leqslant 10^{10}$ and the Prandtl number range $0.7\leqslant Pr \leqslant10$. We find that there exist strong counter-gradient local heat flux with magnitude much larger than the global Nusselt number $Nu$ of the system. Two mechanisms for generating counter-gradient heat transport are identified: one is due to the bulk dynamics and the other is due to the competitions between the corner-flow rolls and the large-scale circulation (LSC). While the magnitude of the former is found to increase with increasing Prandtl number, that of the latter maximizes at medium $Pr$. We further reveal that the corner-LSC competitions lead to the anomalous $Nu$-$Pr$ relation in 2D RB convection, i.e. $Nu(Pr)$ minimizes, rather than maximizes as in three-dimensional cylindrical case, at $Pr\approx2\sim3$ for moderate $Ra$.