The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Model for Predicting the Mean Drop Size in Effervescent Sprays

This paper is focused on the development of a model for predicting the mean drop size in effervescent sprays. A combinatorial approach is followed in this modeling scheme, which is based on energy and entropy principles. The model is implemented in cascade in order to take primary breakup (due to exploding gas bubbles) and secondary breakup (due to shearing action of surrounding medium) into account. The approach in this methodology is to obtain the most probable drop size distribution by maximizing the entropy while satisfying the constraints of mass and energy balance. The comparison of the model predictions with the past experimental data is presented for validation. A careful experimental study is conducted over a wide range of gas-to-liquid ratios, which shows a good agreement with the model predictions: It is observed that the model gives accurate results in bubbly and annular flow regimes. However, discrepancies are observed in the transitional slug flow regime of the atomizer.

Experimental Studies on Mass Distribution of Impinging Jets

Noting that practical impinging injectors are likely to have skewness, an experimental study has been made to understand the behavior of such jets using water as the simulant. In perfectly impinging jets, a high aspect ratio ellipse-like mass distribution pattern is obtained with major axis normal to the plane of two jets whereas in skewed jets the major axis turns from its normal position. A simple analysis shows that this angle of turn is a function of skewness fraction and impingement angle only and is independent of injection velocity. Experimental data from both mass distribution and photographic technique validate this prediction.

Single Cavity Trapped Vortex Combustor Simulations

The present work describes steady and unsteady computation of reacting flow in a Trapped Vortex Combustor. The primary motivation of this study is to develop this concept into a working combustor in modern gas turbines. The present work is an effort towards development of an experimental model test rig for further understanding dynamics of a single cavity trapped vortex combustor. The steady computations with and without combustion have been done for L/D of 0.8, 1 and 1.2; also unsteady non-reacting flow simulation has been done for L/D of 1. Fuel used for the present study is methane and Eddy-Dissipation model has been used for combustion-turbulence interactions. For L/D of 0.8, combustion efficiency is maximum and pattern factor is minimum. Also, primary vortex in the cavity is more stable and symmetric for L/D of 0.8. From unsteady non-reacting flow simulations, it is found that there is no vortex shedding from the cavity but there are oscillations in the span-wise direction of the combustor.

Fluid Mechanics of Aquatic Locomotion at Large Reynolds Numbers

Abstract | There exist a huge range of fish species besides other aquatic organisms like squids and salps that locomote in water at large Reynolds numbers, a regime of flow where inertial forces dominate viscous forces. In the present review, we discuss the fluid mechanics governing the locomotion of such organisms. Most fishes propel themselves by periodic undulatory motions of the body and tail, and the typical classification of their swimming modes is based on the fraction of their body that undergoes such undulatory motions. In the angulliform mode, or the eel type, the entire body undergoes undulatory motions in the form of a travelling wave that goes from head to tail, while in the other extreme case, the thunniform mode, only the rear tail (caudal fin) undergoes lateral oscillations. The thunniform mode of swimming is essentially based on the lift force generated by the airfoil like crosssection of the fish tail as it moves laterally through the water, while the anguilliform mode may be understood using the “reactive theory” of Lighthill. In pulsed jet propulsion, adopted by squids and salps, there are two components to the thrust; the first due to the familiar ejection of momentum and the other due to an over-pressure at the exit plane caused by the unsteadiness of the jet. The flow immediately downstream of the body in all three modes consists of vortex rings; the differentiating point being the vastly different orientations of the vortex rings. However, since all the bodies are self-propelling, the thrust force must be equal to the drag force (at steady speed), implying no net force on the body, and hence the wake or flow downstream must be momentumless. For such bodies, where there is no net force, it is difficult to directly define a propulsion efficiency, although it is possible to use some other very different measures like “cost of transportation” to broadly judge performance.

Electromechanical dynamics and optimization of pectoral fin-based ionic polymer-metal composite underwater propulsor

Ionic polymer-metal composites are soft artificial muscle-like bending actuators, which can work efficiently in wet environments such as water. Therefore, there is significant motivation for research on the development and design analysis of ionic polymer-metal composite based biomimetic underwater propulsion systems. Among aquatic animals, fishes are efficient swimmers with advantages such as high maneuverability, high cruising speed, noiseless propulsion, and efficient stabilization. Fish swimming mechanisms provide biomimetic inspiration for underwater propulsor design. Fish locomotion can be broadly classified into body and/or caudal fin propulsion and median and/or paired pectoral fin propulsion. In this article, the paired pectoral fin-based oscillatory propulsion using ionic polymer-metal composite for aquatic propulsor applications is studied. Beam theory and the concept of hydrodynamic function are used to describe the interaction between the beam and water. Furthermore, a quasi-steady blade element model that accounts for unsteady phenomena such as added mass effects, dynamic stall, and the cumulative Wagner effect is used to obtain hydrodynamic performance of the ionic polymer-metal composite propulsor. Dynamic characteristics of ionic polymer-metal composite fin are analyzed using numerical simulations. It is shown that the use of optimization methods can lead to significant improvement in performance of the ionic polymer-metal composite fin.

Minimum power consumption for drag reduction on a circular cylinder by tangential surface motion

We investigate the effect of a prescribed tangential velocity on the drag force on a circular cylinder in a spanwise uniform cross flow. Using a combination of theoretical and numerical techniques we make an attempt at determining the optimal tangential velocity profiles which will reduce the drag force acting on the cylindrical body while minimizing the net power consumption characterized through a non-dimensional power loss coefficient (C-PL). A striking conclusion of our analysis is that the tangential velocity associated with the potential flow, which completely suppresses the drag force, is not optimal for both small and large, but finite Reynolds number. When inertial effects are negligible (R e << 1), theoretical analysis based on two-dimensional Oseen equations gives us the optimal tangential velocity profile which leads to energetically efficient drag reduction. Furthermore, in the limit of zero Reynolds number (Re -> 0), minimum power loss is achieved for a tangential velocity profile corresponding to a shear-free perfect slip boundary. At finite Re, results from numerical simulations indicate that perfect slip is not optimum and a further reduction in drag can be achieved for reduced power consumption. A gradual increase in the strength of a tangential velocity which involves only the first reflectionally symmetric mode leads to a monotonic reduction in drag and eventual thrust production. Simulations reveal the existence of an optimal strength for which the power consumption attains a minima. At a Reynolds number of 100, minimum value of the power loss coefficient (C-PL = 0.37) is obtained when the maximum in tangential surface velocity is about one and a half times the free stream uniform velocity corresponding to a percentage drag reduction of approximately 77 %; C-PL = 0.42 and 0.50 for perfect slip and potential flow cases, respectively. Our results suggest that potential flow tangential velocity enables energetically efficient propulsion at all Reynolds numbers but optimal drag reduction only for Re -> infinity. The two-dimensional strategy of reducing drag while minimizing net power consumption is shown to be effective in three dimensions via numerical simulation of flow past an infinite circular cylinder at a Reynolds number of 300. Finally a strategy of reducing drag, suitable for practical implementation and amenable to experimental testing, through piecewise constant tangential velocities distributed along the cylinder periphery is proposed and analysed.

Analytical theory and stability analysis of an elongated nanoscale object under external torque

We consider the rotational motion of an elongated nanoscale object in a fluid under an external torque. The experimentally observed dynamics could be understood from analytical solutions of the Stokes equation, with explicit formulae derived for the dynamical states as a function of the object dimensions and the parameters defining the external torque. Under certain conditions, multiple analytical solutions to the Stokes equations exist, which have been investigated through numerical analysis of their stability against small perturbations and their sensitivity towards initial conditions. These experimental results and analytical formulae are general enough to be applicable to the rotational motion of any isolated elongated object at low Reynolds numbers, and could be useful in the design of non-spherical nanostructures for diverse applications pertaining to microfluidics and nanoscale propulsion technologies.

Conformal Cytocompatible Ferrite Coatings Facilitate the Realization of a Nanovoyager in Human Blood

Controlled motion of artificial nanomotors in biological environments, such as blood, can lead to fascinating biomedical applications, ranging from targeted drug delivery to microsurgery and many more. In spite of the various strategies used in fabricating and actuating nanomotors, practical issues related to fuel requirement, corrosion, and liquid viscosity have limited the motion of nanomotors to model systems such as water, serum, or biofluids diluted with toxic chemical fuels, such as hydrogen peroxide. As we demonstrate here, integrating conformal ferrite coatings with magnetic nanohelices offer a promising combination of functionalities for having controlled motion in practical biological fluids, such as chemical stability, cytocompatibility, and the generated thrust. These coatings were found to be stable in various biofluids, including human blood, even after overnight incubation, and did not have significant influence on the propulsion efficiency of the magnetically driven nanohelices, thereby facilitating the first successful ``voyage'' of artificial nanomotors in human blood. The motion of the ``nanovoyager'' was found to show interesting stick-slip dynamics, an effect originating in the colloidal jamming of blood cells in the plasma. The system of magnetic ``nanovoyagers'' was found to be cytocompatible with C2C12 mouse myoblast cells, as confirmed using MTT assay and fluorescence microscopy observations of cell morphology. Taken together, the results presented in this work establish the suitability of the ``nanovoyager'' with conformal ferrite coatings toward biomedical applications.

Independent Positioning of Magnetic Nanomotors

There is considerable interest in powering and maneuvering nanostructures remotely in fluidic media using noninvasive fuel-free methods, for which small homogeneous magnetic fields are ideally suited. Current strategies include helical propulsion of chiral nanostructures, cilia-like motion of flexible filaments, and surface assisted translation of asymmetric colloidal doublets and magnetic nanorods, in all of which the individual structures are moved in a particular direction that is completely tied to the characteristics of the driving fields. As we show in this paper, when we use appropriate magnetic field configurations and actuation time scales, it is possible to maneuver geometrically identical nanostructures in different directions, and subsequently position them at arbitrary locations with respect to each other. The method reported here requires proximity of the nanomotors to a solid surface, and could be useful in applications that require remote and independent control over individual components in microfluidic environments.

扩压段对高超声速推进风洞起动的影响

高超声速推进风洞是进行超燃冲压发动机模型地面模拟的重要试验设备.其中,扩压段的设计非常重要,它不是孤立的,而与主流系统和引射系统密切相关.通过对不同几何形状的扩压段在不同安装位置进行调试,研究在风洞起动过程中整个系统流动状态的演变规律,探索扩压段的优化设计.实验表明,扩压段应安装在距主喷管较近的位置,并具有与主喷管相匹配的入口面积和形状.扩压段的形状应设计成能使主气流通过一系列斜激波串减速的通道.扩压段应具有一定的长度以保证主气流减速到一定程度.

Computational study of viscous effects on lobed mixer flow features and performance

This article describes a computational study of viscous effects on lobed mixer flowfields. The computations, which were carried out using a compressible, three-dimensional, unstructured-mesh Navier-Stokes solver, were aimed at assessing the impacts on mixer performance of inlet boundary-layer thickness and boundary-layer separation within the lobe. The geometries analyzed represent a class of lobed mixer configurations used in turbofan engines. Parameters investigated included lobe penetration angles from 22 to 45 deg, stream-to-stream velocity ratios from 0.5 to 1.0, and two inlet boundary-layer displacement thicknesses. The results show quantitatively the increasing influence of viscous effects as lobe penetration angle is increased. It is shown that the simple estimate of shed circulation given by Skebe et al. (Experimental Investigation of Three-Dimensional Forced Mixer Lobe Flow Field, AIAA Paper 88-3785, July, 1988) can be extended even to situations in which the flow is separated, provided an effective mixer exit angle and height are defined. An examination of different loss sources is also carried out to illustrate the relative contributions of mixing loss and of boundary-layer viscous effects in cases of practical interest.