Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.

Biaxial strain induced electrical inhomogenities and phase separation in the ferromagnetic metallic phase in thin films of La0.7Ca0.3MnO3: A scanning tunneling potentiometry/spectroscopy study

We have investigated the effect of biaxial strain on local electrical/electronic properties in thin films of La0.7Ca0.3MnO3 with varying degrees of biaxial strain in them. The local electrical properties were investigated as a function of temperature by scanning tunneling spectroscopy (STS) and scanning tunneling potentiometry (STP), along with the bulk probe like conductance fluctuations.The results indicate a positive correlation between the lattice mismatch biaxial strain and the local electrical/electronic inhomogenities observed in the strained sample. This is plausible since the crystal structure of the manganites interfere rather strongly with the magnetic/electronic degrees of freedom. Thus even a small imbalance (biaxial strain) can induce significant changes in the electrical properties of the system.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

Effect of screening of intermicellar interactions on the linear and nonlinear rheology of a viscoelastic gel

We report our studies of the linear and nonlinear rheology of aqueous solutions of the surfactant cetyl trimethylammonium tosylate (CTAT) with varying amounts of sodium chloride (NaCl). The CTAT concentration is fixed at 42 mM, and the salt concentration is varied between 0 and 120 mM. On increasing the salt (NaCl) concentration, we see three distinct regimes in the zero-shear viscosity and the high-frequency plateau modulus data. In regime 1, the zero-shear viscosity shows a weak increase with salt concentration due to enhanced micellar growth. The decrease in the zero-shear viscosities with salt concentration in regimes II and III can be explained in terms of intermicellar branching. The most intriguing feature of our data, however, is the anomalous behavior of the high-frequency plateau modulus in regime II (0.12 less than or equal to [NaCl]/[CTAT] less than or equal to 1.42). In this regime, the plateau modulus increases with an increase in NaCl concentration. This is highly interesting, since the correlation length of concentration fluctuations and hence the plateau modulus G(0) are not expected to change appreciably in the semidilute regime. We propose to explain the changes in regime II in terms of a possible unbinding of the organic counterions (tosylate) from the CTA(+) surfaces on the addition of NaCl. In the nonlinear flow curves of the samples with high salt content, significant deviations from the predictions of the Giesekus model for entangled micelles are observed.

Universality and scaling behavior of injected power in elastic turbulence in wormlike micellar gel

We study the statistical properties of spatially averaged global injected power fluctuations for Taylor-Couette flow of a wormlike micellar gel formed by surfactant cetyltrimethylammonium tosylate. At sufficiently high Weissenberg numbers the shear rate, and hence the injected power p(t), at a constant applied stress shows large irregular fluctuations in time. The nature of the probability distribution function (PDF) of p(t) and the power-law decay of its power spectrum are very similar to that observed in recent studies of elastic turbulence for polymer solutions. Remarkably, these non-Gaussian PDFs can be well described by a universal, large deviation functional form given by the generalized Gumbel distribution observed in the context of spatially averaged global measures in diverse classes of highly correlated systems. We show by in situ rheology and polarized light scattering experiments that in the elastic turbulent regime the flow is spatially smooth but random in time, in agreement with a recent hypothesis for elastic turbulence.

On breaks of the Indian monsoon

For over a century, the term break has been used for spells in which the rainfall over the Indian monsoon zone is interrupted. The phenomenon of 'break monsoon' is of great interest because long intense breaks are often associated with poor monsoon seasons. Such breaks have distinct circulation characteristics (heat trough type circulation) and have a large impact on rainfed agriculture. Although interruption of the monsoon rainfall is considered to be the most important feature of the break monsoon, traditionally breaks have been identified on the basis of the surface pressure and wind patterns over the Indian region. We have defined breaks (and active spells) on the basis of rainfall over the monsoon zone. The rainfall criteria are chosen so as to ensure a large overlap with the traditional breaks documented by Ramamurthy (1969) and De et al (1998). We have identified these rainbreaks for 1901-89. We have also identified active spells on the basis of rainfall over the Indian monsoon zone. We have shown that the all-India summer monsoon rainfall is significantly negatively correlated with the number of rainbreak days (correlation coefficient -0.56) and significantly positively correlated with the number of active days (correlation coefficient 0.47). Thus the interannual variation of the all-India summer monsoon rainfall is shown to be related to the number of days of rainbreaks and active spells identified here. There have been several studies of breaks (and also active spells in several cases) identified on the basis of different criteria over regions differing in spatial scales (e.g., Webster et al 1998; Krishnan et al 2000; Goswami and Mohan 2000; and Annamalai and Slingo 2001). We find that there is considerable overlap between the rainbreaks we have identified and breaks based on the traditional definition. There is some overlap with the breaks identified by Krishnan et al (2000) but little overlap with breaks identified by Webster et al (1998). Further, there are three or four active-break cycles in a season according to Webster et al (1998) which implies a time scale of about 40 days for which Goswami and Mohan (2000), and Annamalai and Slingo'(2001) have studied breaks and active minus break fluctuations. On the other hand, neither the traditional breaks (Ramamurthy 1969; and De et al 1998) nor the rainbreaks occur every year. This suggests that the 'breaks' in these studies axe weak spells of the intraseasonal variation of the monsoon, which occur every year. We have derived the OLR and circulation patterns associated with rainbreaks and active spells and compared them with the patterns associated with breaks/active minus break spells from these studies. Inspite of differences in the patterns over the Indian region, there is one feature which is seen in the OLR anomaly patterns of breaks identified on the basis of different criteria as well as the rainbreaks identified in this paper viz., a quadrapole over the Asia-west Pacific region arising from anomalies opposite (same) in sign to those over the Indian region occurring over the equatorial Indian Ocean and northern tropical (equatorial) parts of the west Pacific. Thus it appears that this quadrapole is a basic feature of weak spells of the intraseasonal,variation over the Asia-west Pacific region. Since the rainbreaks are intense weak spells, this basic feature is also seen in the composite patterns of these breaks. We find that rainbreaks (active spells) are also associated with negative (positive) anomalies over a part of the cast Pacific suggesting that the convection over the Indian region is linked to that over the east Pacific not only on the interannual scale (as evinced by the link between the Indian summer monsoon rainfall and ENSO) but on the intraseasonal scale as well.

Shear banding in a yield stress bearing Langmuir monolayer

We examine the shear-thinning behaviour of a two dimensional yield stress bearing monolayer of sorbitan tristearate at air/water interface. The flow curve consists of a linear region at low shear stresses/shear rates, followed by a stress plateau at higher values. The velocity profile obtained from particle imaging velocimetry indicates that shear banding occurs, showing coexistence of the fluidized region near the rotor and solid region with vanishing shear-rate away from the rotor. In the fluidized region, the velocity profile, which is linear at low shear rates, becomes exponential at the onset of shear-thinning, followed by a time varying velocity profile in the plateau region. At low values of constant applied shear rates, the viscosity of the film increases with time, thus showing aging behaviour like in soft glassy three-dimensional (3D) systems. Further, at the low values of the applied stress in the yield stress regime, the shear-rate fluctuations in time show both positive and negative values, similar to that observed in sheared 3D jammed systems. By carrying out a statistical analysis of these shear-rate fluctuations, we estimate the effective temperature of the soft glassy monolayer using the Galavatti-Cohen steady state fluctuation relation.

Modeling of nonreacting and reacting turbulent spray jets using a fully stochastic separated flow approach

Numerical modeling of several turbulent nonreacting and reacting spray jets is carried out using a fully stochastic separated flow (FSSF) approach. As is widely used, the carrier-phase is considered in an Eulerian framework, while the dispersed phase is tracked in a Lagrangian framework following the stochastic separated flow (SSF) model. Various interactions between the two phases are taken into account by means of two-way coupling. Spray evaporation is described using a thermal model with an infinite conductivity in the liquid phase. The gas-phase turbulence terms are closed using the k-epsilon model. A novel mixture fraction based approach is used to stochastically model the fluctuating temperature and composition in the gas phase and these are then used to refine the estimates of the heat and mass transfer rates between the droplets and the surrounding gas-phase. In classical SSF (CSSF) methods, stochastic fluctuations of only the gas-phase velocity are modeled. Successful implementation of the FSSF approach to turbulent nonreacting and reacting spray jets is demonstrated. Results are compared against experimental measurements as well as with predictions using the CSSF approach for both nonreacting and reacting spray jets. The FSSF approach shows little difference from the CSSF predictions for nonreacting spray jets but differences are significant for reacting spray jets. In general, the FSSF approach gives good predictions of the flame length and structure but further improvements in modeling may be needed to improve the accuracy of some details of the Predictions. (C) 2011 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Dark Matter Halo Properties as Deduced from the Observed H I Scale Height Data, in The Low-Frequency Radio Universe

We use the HΙ scale height data along with the HΙ rotation curve as constraints to probe the shape and density profile of the dark matter halos of M31 (Andromeda) and the superthin, low surface brightness (LSB) galaxy UGC 07321. We model the galaxy as a two component system of gravitationally-coupled stars and gas subjected to the force field of a dark matter halo. For M31, we get a flattened halo which is required to match the outer galactic HΙ scale height data, with our best-fit axis ratio (0.4) lying at the most oblate end of the distributions obtained from cosmological simulations. For UGC 07321, our best-fit halo core radius is only slightly larger than the stellar disc scale length, indicating that the halo is important even at small radii in this LSB galaxy. The high value of the gas velocity dispersion required to match the scale height data can explain the low star-formation rate of this galaxy.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

Quantum destruction of spiral order in two-dimensional frustrated magnets

We study the fate of spin-1/2 spiral-ordered two-dimensional quantum antiferromagnets that are disordered by quantum fluctuations. A crucial role is played by the topological point defects of the spiral phase, which are known to have a Z(2) character. Previous works established that a nontrivial quantum spin-liquid phase results when the spiral is disordered without proliferating the Z(2) vortices. Here, we show that when the spiral is disordered by proliferating and condensing these vortices, valence-bond solid ordering occurs due to quantum Berry phase effects. We develop a general theory for this latter phase transition and apply it to a lattice model. This transition potentially provides a new example of a Landau-forbidden deconfined quantum critical point.

Explicit filtering LES for turbulent non-premixed combustion

Two models for large eddy simulation of turbulent reacting flow in homogeneous turbulence were studied. The sub-grid stress arising out of non-linearities of the Navier-Stokes equations were modeled using an explicit filtering approach. A filtered mass density function (FMDF) approach was used for closure of the sub-grid scalar fluctuations. A posteriori calculations, when compared with the results from the direct numerical simulation, indicate that the explicit filtering is adequate in representing the effect of sub-grid stress on the filtered velocity field in the absence of reaction. Discrepancies arise when reactions occur, but the FMDF approach suffices to account for sub-grid scale fluctuations of the reacting scalars, accurately.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Characterization of phase transformation behaviour and microstructural development of electroless Ni-B coating

Phase transformation behaviour of amorphous electroless Ni-B coating with a targeted composition of Ni-6wt% B is characterized in conjunction with microstructural development and hardness. Microscopic observations of the as-deposited coating display a novel microstructure which is already phase separated at multiple length scales. Spherical colonies of similar to 5 mu m consist of 2-3 mu m nodular regions which are surrounded by similar to 2-3 mu m region that contains fine bands ranging from 10 to 70 nm in width. The appearance of three crystalline phases in this binary system at different stages of heat treatment and the concomitant variation in hardness are shown to arise from nanoscale fluctuations in the as-deposited boron content from 4 to 8 wt%. High temperature annealing reveals continuous crystallization up to 430 degrees C, overlapping with the domain of B loss due to diffusion into the substrate. The implications of such a microstructure for optimal heat treatment procedures are discussed. (C) 2011 Elsevier B.V. All rights reserved.

Computation of thermodynamic properties of liquid ferrous alloys from structural measurements

Relation between X-ray scattering intensities, mean square thermal fluctuations and thermodynamic properties. High temperature X-ray diffraction study of liquid Fe-Ni and Fe-Si alloys using reflection and transmission geometries. Calculation of the structure factor as a function of wave vector. Extrapolation to zero wave vector. Calculation of the concentration-concentration correlation function defined by A. B. Bhatia and D. E. Thorton. Computation of thermodynamic quantities of mixing A G, LlH and LlS for the binary alloys. Comparison with direct thermodynamic measurements reported in the literature.