Compressibility behaviour of bentonite and organically modified bentonite slurry

Bentonite in slurry walls needs to be amended with organo-clay to control the migration of organic contaminants. Consolidation behaviour of the slurry is important because it will reduce the total effective stress owing to mobilisation of frictional force between the side wall of the trench and the slurry. Compressibility of the slurry of bentonite is expected to undergo significant changes owing to amendment with organo-clay and according to the nature of the fluid. Standard one-dimensional consolidation tests were carried out on slurries of bentonite, organo-clay and their mixtures by remoulding them to their respective liquid limit consistency with water as well as fluids of low polarity, such as carbon tetrachloride, and inundating with different fluids. Organo-clay and its mixture with bentonite when moulded with water exhibit lower compressibility than bentonite in any pore fluid, but their compressibility increases when moulded with carbon tetrachloride and inundated with the same fluid. These changes in the compressibility of bentonite amended with organo-clay are explained by particle rearrangements, changes in the development of the diffuse double layer and contribution from the water adsorbed in the inter-lamellar space of the clay.

Experimental and simulation studies on the performance of standing wave thermoacoustic prime mover for pulse tube refrigerator

The thermoacoustic prime mover (TAPM) has gained considerable attention as a pressure wave generator to drive pulse tube refrigerator (PTR) due to no moving parts, reasonable efficiency, use of environmental friendly working fluids etc. To drive PTCs, lower frequencies (f) with larger pressure amplitudes (Delta P) are essential, which are affected by geometric and operating parameters of TAPM as well as working fluids. For driving PTRs, a twin standing wave TAPM is built and studied by using different working fluids such as helium, argon, nitrogen and their binary mixtures. Simulation results of DeltaEc are compared with experimental data wherever possible. DeltaEc predicts slightly increased resonance frequencies, but gives larger Delta P and lower temperature difference Delta T across stack. High mass number working fluid leads to lower frequency with larger Delta P, but higher Delta T. Studies indicate that the binary gas mixture of right composition with lower Delta T can be arrived at to drive TAPM of given geometry. (C) 2013 Elsevier Ltd and IIR. All rights reserved.

Effect of Heat Pipe Figure of Merit on an Evaporating Thin Film

The performance of a two-phase heat transport device such as the loop heat pipe is influenced by the evaporative heat transfer coefficient in the evaporator. From previous experiments with loop heat pipes, it has been observed that fluids with a high heat pipe figure of merit have a high heat transfer coefficient. Considering an evaporating extended thin film, this paper theoretically corroborates this experimental observation by deriving a direct link between the evaporative heat flux at the interface and the fluid figures of merit (namely interline heat flow parameter and heat pipe figure of merit) in the thin film. Numerical experiments with different working fluids clearly show that a fluid with high figure of merit also has a high cumulative heat transfer in the microregion encompassing the evaporating thin film. Thus, a loop heat pipe or heat pipe that uses a working fluid with a high interline heat flow parameter and heat pipe figure of merit will lead to a high evaporative heat transfer coefficient.

Spectral response of a droplet in pulsating external flow field

A droplet introduced in an external convective flow field exhibits significant multimodal shape oscillations depending upon the intensity of the aerodynamic forcing. In this paper, a theoretical model describing the temporal evolution of normal modes of the droplet shape is developed. The fluid is assumed to be weakly viscous and Newtonian. The convective flow velocity, which is assumed to be incompressible and inviscid, is incorporated in the model through the normal stress condition at the droplet surface and the equation of motion governing the dynamics of each mode is derived. The coupling between the external flow and the droplet is approximated to be a one-way process, i.e., the external flow perturbations effect the droplet shape oscillations and the droplet oscillation itself does not influence the external flow characteristics. The shape oscillations of the droplet with different fluid properties under different unsteady flow fields were simulated. For a pulsatile external flow, the frequency spectra of the normal modes of the droplet revealed a dominant response at the resonant frequency, in addition to the driving frequency and the corresponding harmonics. At driving frequencies sufficiently different from the resonant frequency of the prolate-oblate oscillation mode of the droplet, the oscillations are stable. But at resonance the oscillation amplitude grows in time leading to breakup depending upon the fluid viscosity. A line vortex advecting past the droplet, simulated as an isotropic jump in the far field velocity, leads to the resonant excitation of the droplet shape modes if and only if the time taken by the vortex to cross the droplet is less than the resonant period of the P-2 mode of the droplet. A train of two vortices interacting with the droplet is also analysed. It shows clearly that the time instant of introduction of the second vortex with respect to the droplet shape oscillation cycle is crucial in determining the amplitude of oscillation. (C) 2014 AIP Publishing LLC.

Characterization and Performance Study of Packaged Micropump for Drug Delivery

In this work, we present the characterization and performance studies of self-priming peristaltic pump for drug delivery application. Conventional materials and methods have been used to fabricate single cam mechanism based peristaltic micropump. To control the fluid flow precisely in micro liter range, a single cam mechanism has been used instead of conventional roller mechanism. The fabricated pump is suitable for liquid, gas and foam. Using water as a fluid medium, a flow rate of 12.5 mu l/rpm is achieved using a flexible silicone tube of inner diameter 1.5 mm and outer diameter 2.5 mm. Other than water, higher viscosity fluids showed a decrease in the flow rate. The designed micropump exhibits a linear dependence of flow rate in the voltage range of 2.5V to 5V. Drug delivery using micropump demands that the micropump has to pump against the blood pressure (maximum of 25kPa) with constant flow rate. Here the designed pump is able to pump the liquid with a constant flow rate of 500 mu l/min (water) up to a backpressure of 40kPa. It was observed that, by increasing the backpressure above 40kPa, flow rate of the pump gradually decreased to 125 mu l/min at 120kPa. In addition, Micropump based drug delivery demands that the micropump should be normally in closed condition in all the positions to avoid drug leakage and bleeding. Hence, micropump has been characterized for normally closed condition in all positions (0 degrees to 360 degrees). However, a minute leak of 0.14 % was found for an inlet pressure of 140kPa. Also, the normally closed region with no leak is observed up to 60kPa of pressure in all positions (0 degrees to 360 degrees).

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.

Insight into instabilities in burning droplets

The complex multiscale physics of nano-particle laden functional droplets in a reacting environment is of fundamental and applied significance for a wide variety of applications ranging from thermal sprays to pharmaceutics to modern day combustors using new brands of bio-fuels. Formation of homogenous nucleated bubbles at the superheat limit inside vaporizing droplets (with or without nanoparticles) represents an unstable system. Here we show that self-induced boiling in burning functional pendant droplets can produce severe volumetric shape oscillations. Internal pressure build-up due to ebullition activity ejects bubbles from the droplet domain causing undulations on the droplet surface and oscillations in bulk. Through experiments, we establish that the degree of droplet deformation depends on the frequency and intensity of these bubble expulsion events. In a distinct regime of single isolated bubble residing in the droplet, however, pre-ejection transient time is identified by Darrieus-Landau evaporative instability, where bubble-droplet system behaves as a synchronized driver-driven system with bulk bubble-shape oscillations being imposed on the droplet. The agglomeration of nanophase additives modulates the flow structures within the droplet and also influences the bubble inception and growth leading to different levels of instabilities. (C) 2014 AIP Publishing LLC.

Flame surface statistics of constant-pressure turbulent expanding premixed flames

In this paper we investigate the local flame surface statistics of constant-pressure turbulent expanding flames. First the statistics of local length ratio is experimentally determined from high-speed planar Mie scattering images of spherically expanding flames, with the length ratio on the measurement plane, at predefined equiangular sectors, defined as the ratio of the actual flame length to the length of a circular-arc of radius equal to the average radius of the flame. Assuming isotropic distribution of such flame segments we then convolute suitable forms of the length-ratio probability distribution functions (pdfs) to arrive at the corresponding area-ratio pdfs. It is found that both the length ratio and area ratio pdfs are near log-normally distributed and shows self-similar behavior with increasing radius. Near log-normality and rather intermittent behavior of the flame-length ratio suggests similarity with dissipation rate quantities which stimulates multifractal analysis. (C) 2014 AIP Publishing LLC.

Transition in vortex breakdown modes in a coaxial isothermal unconfined swirling jet

This paper reports first observations of transition in recirculation pattern from an open-bubble type axisymmetric vortex breakdown to partially open bubble mode through an intermediate, critical regime of conical sheet formation in an unconfined, co-axial isothermal swirling flow. This time-mean transition is studied for two distinct flow modes which are characterized based on the modified Rossby number (Ro(m)), i.e., Ro(m) <= 1 and Ro(m) > 1. Flow modes with Ro(m) <= 1 are observed to first undergo cone-type breakdown and then to partially open bubble state as the geometric swirl number (S-G) is increased by similar to 20% and similar to 40%, respectively, from the baseline open-bubble state. However, the flow modes with Ro(m) > 1 fail to undergo such sequential transition. This distinct behavior is explained based on the physical significance associated with Ro(m) and the swirl momentum factor (xi). In essence, xi represents the ratio of angular momentum distributed across the flow structure to that distributed from central axis to the edge of the vortex core. It is observed that xi increases by similar to 100% in the critical swirl number band where conical breakdown occurs as compared to its magnitude in the S-G regime where open bubble state is seen. This results from the fact that flow modes with Ro(m) <= 1 are dominated by radial pressure gradient due to swirl/rotational effect when compared to radial pressure deficit arising from entrainment (due to the presence of co-stream). Consequently, the imparted swirl tends to penetrate easily towards the central axis causing it to spread laterally and finally undergo conical sheet breakdown. However, the flow modes with Ro(m) > 1 are dominated by pressure deficit due to entrainment effect. This blocks the radial inward penetration of imparted angular momentum thus preventing the lateral spread of these flow modes. As such these structures fail to undergo cone mode of vortex breakdown which is substantiated by a mere 30%-40% rise in xi in the critical swirl number range. (C) 2014 AIP Publishing LLC.

Multimodal shape oscillations of droplets excited by an air stream

The shape dynamics of droplets exposed to an air jet at intermediate droplet Reynolds numbers is investigated. High speed imaging and hot-wire anemometry are employed to examine the mechanism of droplet oscillation. The theory that the vortex shedding behind the droplet induces oscillation is examined. In these experiments, no particular dominant frequency is found in the wake region of the droplet. Hence the inherent free-stream disturbances prove to be driving the droplet oscillations. The modes of droplet oscillation show a band of dominant frequencies near the corresponding natural frequency, further proving that there is no particular forcing frequency involved. In the frequency spectrum of the lowest mode of oscillation for glycerol at the highest Reynolds number, no response is observed below the threshold frequency corresponding to the viscous dissipation time scale. This selective suppression of lower frequencies in the case of glycerol is corroborated by scaling arguments. The influence of surface tension on the droplet oscillation is studied using ethanol as a test fluid. Since a lower surface tension reduces the natural frequency, ethanol shows lower excited frequencies. The oscillation levels of different fluids are quantified using the droplet aspect ratio and correlated in terms of Weber number and Ohnesorge number. (C) 2014 Elsevier Ltd. All rights reserved.

Kinetics of phase separation in polymer mixtures: A molecular dynamics study

We present detailed results from a molecular dynamics (MD) simulation of phase-separation kinetics in polymer mixtures. Our MD simulations naturally incorporate hydrodynamic effects. We find that polymeric phase separation (with dynamically symmetric components) is in the same universality class as segregation of simple fluids: the degree of polymerization only slows down the segregation kinetics. For d = 2 polymeric fluids, the domain growth law is L(t) similar to t(phi) with phi showing a crossover from 1/3 -> 1/2 -> 2/3. For d = 3 polymeric fluids, we see the crossover phi = 1/3 -> 1. Our MD simulations do not yet access the inertial hydrodynamic regime (with L similar to t(2/3)) of phase separation in 3-d fluids. (C) 2014 AIP Publishing LLC.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Shear-imposed falling film

The study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469-485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.