Showing papers in "arXiv: Fluid Dynamics in 2012"

Hydrodynamics of self-propulsion near a boundary: predictions and accuracy of far-field approximations

Abstract: The swimming trajectories of self-propelled organisms or synthetic devices in a viscous fluid can be altered by hydrodynamic interactions with nearby boundaries. We explore a multipole description of swimming bodies and provide a general framework for studying the fluid-mediated modifications to swimming trajectories. The validity of the far-field description is probed for a selection of model swimmers of varying geometry and propulsive activity by comparison with full numerical simulations. The reduced model is then used to deliver simple but accurate predictions of hydrodynamically generated wall attraction and pitching dynamics, and may help to explain a number of experimental results.

Energy harvesting efficiency of piezoelectric flags in axial flows

Abstract: Self-sustained oscillations resulting from fluid-solid instabilities, such as the flutter of a flexible flag in axial flow, can be used to harvest energy if one is able to convert the solid energy into electricity. Here, this is achieved using piezoelectric patches attached to the surface of the flag that convert the solid deformation into an electric current powering purely resistive output circuits. Nonlinear numerical simulations in the slender-body limit, based on an explicit description of the coupling between the fluid-solid and electric systems, are used to determine the harvesting efficiency of the system, namely the fraction of the flow kinetic energy flux effectively used to power the output circuit, and its evolution with the system's parameters. The role of the tuning between the characteristic frequencies of the fluid-solid and electric systems is emphasized, as well as the critical impact of the piezoelectric coupling intensity. High fluid loading, classically associated with destabilization by damping, leads to greater energy harvesting, but with a weaker robustness to flow velocity fluctuations due to the sensitivity of the flapping mode selection. This suggests that a control of this mode selection by a careful design of the output circuit could provide some opportunities of improvement for the efficiency and robustness of the energy harvesting process.

Effect of oxidation on the Mechanical Properties of Liquid Gallium and Eutectic Gallium-Indium

Abstract: Liquid metals exhibit remarkable mechanical properties, in particular large surface tension and low viscosity. However, these properties are greatly affected by oxidation when exposed to air. We measure the viscosity, surface tension, and contact angle of gallium (Ga) and a eutectic gallium-indium alloy (eGaIn) while controlling such oxidation by surrounding the metal with an acid bath of variable concentration. Rheometry measurements reveal a yield stress directly attributable to an oxide skin that obscures the intrinsic behavior of the liquid metals. We demonstrate how the intrinsic viscosity can be obtained with precision through a scaling technique that collapses low- and high-Reynolds number data. Measuring surface tension with a pendant drop method, we show that the oxide skin generates a surface stress that mimics surface tension and develop a simple model to relate this to the yield stress obtained from rheometry. We find that yield stress, surface tension, and contact angle all transition from solid-like to liquid behavior at the same critical acid concentration, thereby quantitatively confirming that the wettability of these liquid metals is due to the oxide skin.

A Novel Model Order Reduction Approach for Navier-Stokes Equations at High Reynolds Number

Abstract: A new approach to model order reduction of the Navier-Stokes equations at high Reynolds number is proposed. Unlike traditional approaches, this method does not rely on empirical turbulence modeling or modification of the Navier-Stokes equations. It provides spatial basis functions different from the usual proper orthogonal decomposition basis function in that, in addition to optimally representing the training data set, the new basis functions also provide stable and accurate reduced-order models. The proposed approach is illustrated with two test cases: two-dimensional flow inside a square lid-driven cavity and a two-dimensional mixing layer.

Short time dynamics of viscous drop spreading

Abstract: Liquid drops start spreading directly after coming into contact with a solid sub- strate. Although this phenomenon involves a three-phase contact line, the spread- ing motion can be very fast. We experimentally study the initial spreading dy- namics, characterized by the radius of the wetted area, for viscous drops. Using high-speed imaging with synchronized bottom and side views gives access to 6 decades of time resolution. We show that short time spreading does not exhibit a pure power-law growth. Instead, we find a spreading velocity that decreases logarithmically in time, with a dynamics identical to that of coalescing viscous drops. Remarkably, the contact line dissipation and wetting effects turn out to be unimportant during the initial stages of drop spreading.

Self-propulsion in viscoelastic fluids: pushers vs. pullers

Abstract: We use numerical simulations to address locomotion at zero Reynolds number in viscoelastic (Giesekus) fluids. The swimmers are assumed to be spherical, to self-propel using tangential surface deformation, and the computations are implemented using a finite element method. The emphasis of the study is on the change of the swimming kinematics, energetics, and flow disturbance from Newtonian to viscoelastic, and on the distinction between pusher and puller swimmers. In all cases, the viscoelastic swimming speed is below the Newtonian one, with a minimum obtained for intermediate values of the Weissenberg number, $We$. An analysis of the flow field places the origin of this swimming degradation in non-Newtonian elongational stresses. The power required for swimming is also systematically below the Newtonian power, and always a decreasing function of $We$. A detail energetic balance of the swimming problem points at the polymeric part of the stress as the primary $We$-decreasing energetic contribution, while the contributions of the work done by the swimmer from the solvent remain essentially $We$-independent. In addition, we observe negative values of the polymeric power density in some flow regions, indicating positive elastic work by the polymers on the fluid. The hydrodynamic efficiency, defined as the ratio of the useful to total rate of work, is always above the Newtonian case, with a maximum relative value obtained at intermediate Weissenberg numbers. Finally, the presence of polymeric stresses leads to an increase of the rate of decay of the flow velocity in the fluid, and a decrease of the magnitude of the stresslet governing the magnitude of the effective bulk stress in the fluid.

Optimal Taylor-Couette flow: direct numerical simulations

Abstract: We numerically simulate turbulent Taylor-Couette flow for independently rotating inner and outer cylinders, focusing on the analogy with turbulent Rayleigh-B\'enard flow. Reynolds numbers of $Re_i=8\cdot10^3$ and $Re_o=\pm4\cdot10^3$ of the inner and outer cylinders, respectively, are reached, corresponding to Taylor numbers Ta up to $10^8$. Effective scaling laws for the torque and other system responses are found. Recent experiments with the Twente turbulent Taylor-Couette ($T^3C$) setup and with a similar facility in Maryland at very high Reynolds numbers have revealed an optimum transport at a certain non-zero rotation rate ratio $a = -\omega_o / \omega_i$ of about $a_{opt}=0.33-0.35$. For large enough $Ta$ in the numerically accessible range we also find such an optimum transport at non-zero counter-rotation. The position of this maximum is found to shift with the driving, reaching a maximum of $a_{opt}=0.15$ for $Ta=2.5\cdot10^7$. An explanation for this shift is elucidated, consistent with the experimental result that $a_{opt}$ becomes approximately independent of the driving strength for large enough Reynolds numbers. We furthermore numerically calculate the angular velocity profiles and visualize the different flow structures for the various regimes. By writing the equations in a frame co-rotating with the outer cylinder a link is found between the local angular velocity profiles and the global transport quantities.

Coalescence of Liquid Drops: Different Models Versus Experiment

Abstract: The process of coalescence of two identical liquid drops is simulated numerically in the framework of two essentially different mathematical models, and the results are compared with experimental data on the very early stages of the coalescence process reported recently. The first model tested is the `conventional' one, where it is assumed that coalescence as the formation of a single body of fluid occurs by an instant appearance of a liquid bridge smoothly connecting the two drops, and the subsequent process is the evolution of this single body of fluid driven by capillary forces. The second model under investigation considers coalescence as a process where a section of the free surface becomes trapped between the bulk phases as the drops are pressed against each other, and it is the gradual disappearance of this `internal interface' that leads to the formation of a single body of fluid and the conventional model taking over. Using the full numerical solution of the problem in the framework of each of the two models, we show that the recently reported electrical measurements probing the very early stages of the process are better described by the interface formation/disappearance model. New theory-guided experiments are suggested that would help to further elucidate the details of the coalescence phenomenon. As a by-product of our research, the range of validity of different `scaling laws' advanced as approximate solutions to the problem formulated using the conventional model is established.

Resonant behaviour of an oscillating wave energy converter in a channel

Abstract: A mathematical model is developed to study the behaviour of an oscillating wave energy converter in a channel. During recent laboratory tests in a wave tank, peaks in the hydrodynamic actions on the converter occurred at certain frequencies of the incident waves. This resonant mechanism is known to be generated by the transverse sloshing modes of the channel. Here the influence of the channel sloshing modes on the performance of the device is further investigated. Within the framework of a linear inviscid potential-flow theory, application of the Green theorem yields a hypersingular integral equation for the velocity potential in the fluid domain. The solution is found in terms of a fast-converging series of Chebyshev polynomials of the second kind. The physical behaviour of the system is then analysed, showing sensitivity of the resonant sloshing modes to the geometry of the device, that concurs in increasing the maximum efficiency. Analytical results are validated with available numerical and experimental data.

On the triggering of the Ultimate Regime of convection

Abstract: Rayleigh-B\'enard cells are one of the simplest systems to explore the laws of natural convection in the highly turbulent limit. However, at very high Rayleigh numbers (Ra > 1E12) and for Prandtl numbers of order one, experiments fall into two categories: some evidence a steep enhancement of the heat transfer while others do not. The origin of this apparent disagreement is presently unexplained. This puzzling situation motivated a systematic study of the triggering of the regime with an enhanced heat transfer, originally named the "Ultimate Regime" of convection. High accuracy heat transfer measurements have been conducted in convection cells with various aspect ratios and different specificities, such as altered boundary conditions or obstacles inserted in the flow. The two control parameters, the Rayleigh and Prandtl numbers have been varied independently to disentangle their relative influence. Among other results, it is found that i) most experiments reaching very high $Ra$ are not in disagreement if small differences in Prandtl numbers are taken into account, ii) the transition is not directly triggered by the large scale circulation present in the cell, iii) the sidewall of the cell have a significant influence on the transition. The characteristics of this Ultimate regime are summarized and compared with R. Kraichnan prediction for the asymptotic regime of convection.

The granular silo as a continuum plastic flow: the hour-glass vs the clepsydra

Abstract: The granular silo is one of the many interesting illustrations of the thixotropic property of granular matter: a rapid flow develops at the outlet, propagating upwards through a dense shear flow while material at the bottom corners of the container remains static. For large enough outlets, the discharge flow is continuous; however, by contrast with the clepsydra for which the flow velocity depends on the height of fluid left in the container, the discharge rate of granular silos is constant. Implementing a plastic rheology in a 2D Navier-Stokes solver (following the mu(I)-rheology or a constant friction), we simulate the continuum counterpart of the granular silo. Doing so, we obtain a constant flow rate during the discharge and recover the Beverloo scaling independently of the initial filling height of the silo. We show that lowering the value of the coefficient of friction leads to a transition toward a different behavior, similar to that of a viscous fluid, and where the filling height becomes active in the discharge process. The pressure field shows that large enough values of the coefficient of friction ($\simeq$ 0.3) allow for a low-pressure cavity to form above the outlet, and can thus explain the Beverloo scaling. In conclusion, the difference between the discharge of a hourglass and a clepsydra seems to reside in the existence or not of a plastic yield stress.

Direct Numerical Simulations of Local and Global Torque in Taylor-Couette Flow up to Re=30.000

Abstract: The torque in turbulent Taylor-Couette flows for shear Reynolds numbers Re_S up to 3x10^4 at various mean rotations is studied by means of direct numerical simulations for a radius ratio of \eta=0.71. Convergence of simulations is tested using three criteria of which the agreement of dissipation values estimated from the torque and from the volume dissipation rate turns out to be most demanding. We evaluate the influence of Taylor vortex heights on the torque for a stationary outer cylinder and select a value of the aspect ratio of \Gamma=2, close to the torque maximum. The connection between the torque and the transverse current J^\omega of azimuthal motion which can be computed from the velocity field enables us to investigate the local transport resulting in the torque. The typical spatial distribution of the individual convective and viscous contributions to the local current is analysed for a turbulent flow case. To characterise the turbulent statistics of the transport, PDF's of local current fluctuations are compared to experimental wall shear stress measurements. PDF's of instantaneous torques reveal a fluctuation enhancement in the outer region for strong counter-rotation. Moreover, we find for simulations realising the same shear Re_S>=2x10^4 the formation of a torque maximum for moderate counter-rotation with angular velocities \omega_o\approx-0.4\omega_i. In contrast, for Re_S<=4x10^3 the torque features a maximum for a stationary outer cylinder. In addition, the effective torque scaling exponent is shown to also depend on the mean rotation state. Finally, we evaluate a close connection between boundary-layer thicknesses and the torque.

Stabilizing effect of flexibility in the wake of a flapping foil

Abstract: The wake of a flexible foil undergoing pitching oscillations in a low-speed hydrodynamic tunnel is used to examine the effect of chord-wise foil flexibility in the dynamical features of flapping-based propulsion. We compare the regime transitions in the wake with respect to the case of a rigid foil and show that foil flexibility inhibits the symmetry breaking of the reverse B\'enard-von K\'arm\'an wake reported in the literature. A momentum balance calculation shows the average thrust to be up to three times greater for the flexible foil than for the rigid foil. We explain both of these observations by analyzing the vortex dynamics in the very near wake.

Revealing the state space of turbulent pipe flow by symmetry reduction

Abstract: Symmetry reduction by the method of slices is applied to pipe flow in order to quotient the stream-wise translation and azimuthal rotation symmetries of turbulent flow states. Within the symmetry-reduced state space, all travelling wave solutions reduce to equilibria, and all relative periodic orbits reduce to periodic orbits. Projections of these solutions and their unstable manifolds from their $\infty$-dimensional symmetry-reduced state space onto suitably chosen 2- or 3-dimensional subspaces reveal their interrelations and the role they play in organising turbulence in wall-bounded shear flows. Visualisations of the flow within the slice and its linearisation at equilibria enable us to trace out the unstable manifolds, determine close recurrences, identify connections between different travelling wave solutions, and find, for the first time for pipe flows, relative periodic orbits that are embedded within the chaotic attractor, which capture turbulent dynamics at transitional Reynolds numbers.

Forced generation of simple and double emulsions in all-aqueous system

Abstract: We report an easy-to-implement method that allows the direct generation of water-in-water (w/w) single emulsions. The method relies on direct perturbation of the pressure that drives the flow of the dispersed phase of the emulsions. The resultant inner jet is induced to break up into droplets due to the growth of the perturbation through Rayleigh-Plateau instability [L. Rayleigh, Proc. R. Soc. London 29, 71-97 (1879)]; this leads to the formation of monodisperse droplets. By implementing this method on a modified microfluidic device, we directly generate water-in-water-in-water (w/w/w) double emulsions with good control over the size and the number of encapsulated droplets. Our approach suggests a new route to apply droplet-based microfluidics to completely water-based systems.

Viscous Marangoni propulsion

Abstract: Marangoni propulsion is a form of locomotion wherein an asymmetric release of surfactant by a body located at the surface of a liquid leads to its directed motion. We present in this paper a mathematical model for Marangoni propulsion in the viscous regime. We consider the case of a thin rigid circular disk placed at the surface of a viscous fluid and whose perimeter has a prescribed concentration of an insoluble surfactant, to which the rest of its surface is impenetrable. Assuming a linearized equation of state between surface tension and surfactant concentration, we derive analytically the surfactant, velocity and pressure fields in the asymptotic limit of low Capillary, Peclet and Reynolds numbers. We then exploit these results to calculate the Marangoni propulsion speed of the disk. Neglecting the stress contribution from Marangoni flows is seen to over-predict the propulsion speed by 50%.

Lift Evaluation of a 2D Flapping Flat Plate

Abstract: Several previous experimental and theoretical studies have shown that a leading edge vortex (LEV) on an airfoil or wing can provide lift enhancement. In this paper, unsteady 2D potential flow theory is employed to model the flow field of a flapping flat plate wing. A multi-vortices model is developed to model both the leading edge and trailing edge vortices (TEVs), which offers improved accuracy compared with using only single vortex at each separation location. The lift is obtained by integrating the unsteady Blasius equation. It is found that the motion of vortices contributes significantly to the overall aerodynamic force on the flat plate. The shedding of TEVs and the stabilization of LEVs explicitly contributes to lift enhancement. A Kutta-like condition is used to determine the vortex intensity and location at the leading edge for large angle of attack cases; however, it is proposed to relax this condition for small angle of attack cases and apply a 2D shear layer model to calculate the circulation of the new added vortex. The results of the simulation are then compared with classical numerical, modeled and experimental data for canonical unsteady flat plat problems. Good agreement with these data is observed. Moreover, these results suggested that the leading edge vortex shedding for small angles of attack should be modeled differently than that for large angles of attack. Finally, the results of vortex motion vs. lift indicate that both a motion against the streamwise direction of the LEV and a streamwise motion of the TEV contributes positive lift. This also provides the insights for future active flow control of MAVs that the formation and shedding process of LEVs and TEVs can be manipulated to provide lift enhancement.

Heat transport by turbulent Rayleigh-Bénard convection for $\Pra\ \simeq 0.8$ and $3\times 10^{12} \alt \Ra\ \alt 10^{15}$: Aspect ratio $Γ= 0.50$

Abstract: We report experimental results for heat-transport measurements, in the form of the Nusselt number \Nu, by turbulent Rayleigh-B\'enard convection in a cylindrical sample of aspect ratio $\Gamma \equiv D/L = 0.50$ ($D = 1.12$ m is the diameter and $L = 2.24$ m the height). The measurements were made using sulfur hexafluoride at pressures up to 19 bars as the fluid. They are for the Rayleigh-number range $3\times 10^{12} \alt \Ra \alt 10^{15}$ and for Prandtl numbers \Pra\ between 0.79 and 0.86. For $\Ra < \Ra^*_1 \simeq 1.4\times 10^{13}$ we find $\Nu = N_0 \Ra^{\gamma_{eff}}$ with $\gamma_{eff} = 0.312 \pm 0.002$, consistent with classical turbulent Rayleigh-B\'enard convection in a system with laminar boundary layers below the top and above the bottom plate. For $\Ra^*_1 < \Ra < \Ra^*_2$ (with $\Ra^*_2 \simeq 5\times 10^{14}$) $\gamma_{eff}$ gradually increases up to $0.37\pm 0.01$. We argue that above $\Ra^*_2$ the system is in the ultimate state of convection where the boundary layers, both thermal and kinetic, are also turbulent. Several previous measurements for $\Gamma = 0.50$ are re-examined and compared with the present results.

Turbulent velocity spectra in superfluid flows

Abstract: We present velocity spectra measured in three cryogenic liquid 4He steady flows: grid and wake flows in a pressurized wind tunnel capable of achieving mean velocities up to 5 m/s at temperatures above and below the superfluid transition, down to 1.7 K, and a "chunk" turbulence flow at 1.55 K, capable of sustaining mean superfluid velocities up to 1.3 m/s. Depending on the flows, the stagnation pressure probes used for anemometry are resolving from one to two decades of the inertial regime of the turbulent cascade. We do not find any evidence that the second order statistics of turbulence below the superfluid transition differ from the ones of classical turbulence, above the transition.

How to Calculate Tortuosity Easily

Abstract: Tortuosity is one of the key parameters describing the geometry and transport properties of porous media. It is defined either as an average elongation of fluid paths or as a retardation factor that measures the resistance of a porous medium to the flow. However, in contrast to a retardation factor, an average fluid path elongation is difficult to compute numerically and, in general, is not measurable directly in experiments. We review some recent achievements in bridging the gap between the two formulations of tortuosity and discuss possible method of numerical and an experimental measurements of the tortuosity directly from the fluid velocity field.

Dispersion stability and thermal conductivity of propylene glycol-based nanofluids

Abstract: The dispersion stability and thermal conductivity of propylene glycol based nanofluids containing Al2O3 and TiO2 nanoparticles were studied in the temperature range of 20 to 80 {\deg}C. Nanofluids with different concentrations of nanoparticles were formulated by the two-step method without use of dispersants. In contrast to the common belief the average particle size of nanofluids was observed to decrease with increasing temperature. The nanofluids showed excellent stability over the temperature range of interest. Thermal conductivity enhancement for both of studied nanofluids was a non-linear function of concentration while was temperature independent. Theoretical analyses were performed using existing models and comparisons were made with experimental results. The model based on the aggregation theory appears to yield the best fit. Keywords: Nanofluids, Propylene glycol, Alumina nanoparticles, Titania nanoparticles, Thermal conductivity, Dispersion stability.

Reactive-infiltration instabilities in rocks. Fracture dissolution

Abstract: A reactive fluid dissolving the surface of a uniform fracture will trigger an instability in the dissolution front, leading to spontaneous formation of pronounced well-spaced channels in the surrounding rock matrix. Although the underlying mechanism is similar to the wormhole instability in porous rocks there are significant differences in the physics, due to the absence of a steadily propagating reaction front. In previous work we have described the geophysical implications of this instability in regard to the formation of long conduits in soluble rocks. Here we describe a more general linear stability analysis, including axial diffusion, transport limited dissolution, non-linear kinetics, and a finite length system.

Changes in turbulent dissipation in a channel flow with oscillating walls

Abstract: Harmonic oscillations of the walls of a turbulent plane channel flow are studied by direct numerical simulations to improve our understanding of the physical mechanism for skin-friction drag reduction. The simulations are carried out at constant pressure gradient in order to define an unambiguous inner scaling: in this case, drag reduction manifests itself as an increase of mass flow rate. Energy and enstrophy balances, carried out to emphasize the role of the oscillating spanwise shear layer, show that the viscous dissipations of the mean flow and of the turbulent fluctuations increase with the mass flow rate, and the relative importance of the latter decreases. We then focus on the turbulent enstrophy: through an analysis of the temporal evolution from the beginning of the wall motion, the dominant, oscillation-related term in the turbulent enstrophy is shown to cause the turbulent dissipation to be enhanced in absolute terms, before the slow drift towards the new quasi-equilibrium condition. This mechanism is found to be responsible for the increase in mass flow rate. We finally show that the time-average volume integral of the dominant term relates linearly to the drag reduction.

A nonlinear Schr\"odinger equation for water waves on finite depth with constant vorticity

Abstract: A nonlinear Schr\"odinger equation for the envelope of two dimensional surface water waves on finite depth with non zero constant vorticity is derived, and the influence of this constant vorticity on the well known stability properties of weakly nonlinear wave packets is studied. It is demonstrated that vorticity modifies significantly the modulational instability properties of weakly nonlinear plane waves, namely the growth rate and bandwidth. At third order we have shown the importance of the coupling between the mean flow induced by the modulation and the vorticity. Furthermore, it is shown that these plane wave solutions may be linearly stable to modulational instability for an opposite shear current independently of the dimensionless parameter kh, where k and h are the carrier wavenumber and depth respectively.

On the Three-dimensional Central Moment Lattice Boltzmann Method

Abstract: A three-dimensional (3D) lattice Boltzmann method based on central moments is derived. Two main elements are the local attractors in the collision term and the source terms representing the effect of external and/or self-consistent internal forces. For suitable choices of the orthogonal moment basis for the three-dimensional, twenty seven velocity (D3Q27), and, its subset, fifteen velocity (D3Q15) lattice models, attractors are expressed in terms of factorization of lower order moments as suggested in an earlier work; the corresponding source terms are specified to correctly influence lower order hydrodynamic fields, while avoiding aliasing effects for higher order moments. These are achieved by successively matching the corresponding continuous and discrete central moments at various orders, with the final expressions written in terms of raw moments via a transformation based on the binomial theorem. Furthermore, to alleviate the discrete effects with the source terms, they are treated to be temporally semi-implicit and second-order, with the implicitness subsequently removed by means of a transformation. As a result, the approach is frame-invariant by construction and its emergent dynamics describing fully 3D fluid motion in the presence of force fields is Galilean invariant. Numerical experiments for a set of benchmark problems demonstrate its accuracy.

Shell Models of Magnetohydrodynamic Turbulence

Abstract: Shell models of hydrodynamic turbulence originated in the seventies. Their main aim was to describe the statistics of homogeneous and isotropic turbulence in spectral space, using a simple set of ordinary differential equations. In the eighties, shell models of magnetohydrodynamic (MHD) turbulence emerged based on the same principles as their hydrodynamic counter-part but also incorporating interactions between magnetic and velocity fields. In recent years, significant improvements have been made such as the inclusion of non-local interactions and appropriate definitions for helicities. Though shell models cannot account for the spatial complexity of MHD turbulence, their dynamics are not over simplified and do reflect those of real MHD turbulence including intermittency or chaotic reversals of large-scale modes. Furthermore, these models use realistic values for dimensionless parameters (high kinetic and magnetic Reynolds numbers, low or high magnetic Prandtl number) allowing extended inertial range and accurate dissipation rate. Using modern computers it is difficult to attain an inertial range of three decades with direct numerical simulations, whereas eight are possible using shell models. In this review we set up a general mathematical framework allowing the description of any MHD shell model. The variety of the latter, with their advantages and weaknesses, is introduced. Finally we consider a number of applications, dealing with free-decaying MHD turbulence, dynamo action, Alfven waves and the Hall effect.

Emergence and equilibration of jets in beta-plane turbulence: applications of Stochastic Structural Stability Theory

Abstract: Stochastic Structural Stability Theory (S3T) provides analytical methods for understanding the emergence and equilibration of jets from the turbulence in planetary atmospheres based on the dynamics of the statistical mean state of the turbulence closed at second order. Predictions for formation and equilibration of turbulent jets made using S3T are critically compared with results of simulations made using the associated quasi-linear and nonlinear models. S3T predicts the observed bifurcation behavior associated with the emergence of jets, their equilibration and their breakdown as a function of parameters. Quantitative differences in bifurcation parameter values between predictions of S3T and results of nonlinear simulations are traced to modification of the eddy spectrum which results from two processes: nonlinear eddy-eddy interactions and formation of discrete non-zonal structures. Remarkably, these non-zonal structures, which substantially modify the turbulence spectrum, are found to arise from S3T instability. Formation as linear instabilities and equilibration at finite amplitude of multiple equilibria for identical parameter values in the form of jets with distinct meridional wavenumbers is verified as is the existence at equilibrium of finite amplitude non-zonal structures in the form of nonlinearly modified Rossby waves. When zonal jets and nonlinearly modified Rossby waves coexist at finite amplitude the jet structure is generally found to dominate even if it is linearly less unstable. The physical reality of the manifold of S3T jets and non-zonal structures is underscored by the existence in nonlinear simulations of jet structure at subcritical S3T parameter values which are identified with stable S3T jet modes excited by turbulent fluctuations.

Fluid transport by individual microswimmers

Abstract: We discuss the path of a tracer particle as a microswimmer moves past on an infinite straight trajectory. If the tracer is sufficiently far from the path of the swimmer it moves in a closed loop. As the initial distance between the tracer and the path of the swimmer $\rho$ decreases, the tracer is displaced a small distance backwards (relative to the direction of the swimmer velocity). For much smaller tracer-swimmer separations, however, the tracer displacement becomes positive and diverges as $\rho \to 0$. To quantify this behaviour we calculate the Darwin drift, the total volume swept out by a material sheet of tracers, initially perpendicular to the swimmer path, during the swimmer motion. We find that the drift can be written as the sum of a {\em universal} term which depends on the quadrupolar flow field of the swimmer, together with a non-universal contribution given by the sum of the volumes of the swimmer and its wake. The formula is compared to exact results for the squirmer model and to numerical calculations for a more realistic model swimmer.

Highly focused supersonic microjets: numerical simulations

Abstract: By focusing a laser pulse inside a capillary partially filled with liquid, a vapour bubble is created which emits a pressure wave. This pressure wave travels through the liquid and creates a fast, focused axisymmetric microjet when it is reflected at the meniscus. We numerically investigate the formation of this microjet using axisymmetric boundary-integral simulations, where we model the pressure wave as a pressure pulse applied on the bubble. We find a good agreement between the simulations and experimental results in terms of the time evolution of the jet and on all parameters that can be compared directly. We present a simple analytical model that accurately predicts the velocity of the jet after the pressure pulse and its maximum velocity.

Aggregation and settling in aqueous polydisperse alumina nanoparticle suspensions

Abstract: Nanoparticle suspensions (also called nanofluids) are often polydisperse and tend to settle with time. Settling kinetics in these systems are known to be complex and hence challenging to understand. In this work, polydisperse spherical alumina (Al2O3) nanoparticles in the size range of ~10-100nm were dispersed in water and examined for aggregation and settling behaviour near its isoelectric point (IEP). A series of settling experiments were conducted and the results were analysed by photography and by Small Angle X-ray Scattering (SAXS). The settling curve obtained from standard bed height measurement experiments indicated two different types of behaviour, both of which were also seen in the SAXS data. But the SAXS data were remarkably able to pick out the rapid settling regime as a result of the high temporal resolution (10s) used. By monitoring the SAXS intensity, it was further possible to record the particle aggregation process for the first time. Optical microscopy images were produced on drying and dried droplets extracted from the suspension at various times. Dried deposits showed the rapid decrease in the number of very large particles with time which qualitatively validates the SAXS prediction, and therefore its suitability as a tool to study unstable polydisperse colloids. Keywords: Nanoparticles, nanofluids, polydisperse, aggregation, settling, alumina, microscopy, SAXS