Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.

Fluid dynamics of the monsoon

The monsoonal regions of the world are characterized by a seasonal reversal in the direction of winds associated with the excursion of the equatorial trough (or the ITCZ) in response to the variation in the latitude of maximum insolation. This monsoonal circulation is a planetary scale phenomenon. However, the associated precipitation is critically dependent on the organization of the cumulus clouds (typically a few kilometers in horizontal extent) over the scale of synoptic vortices (typically a few hundred kilometers in horizontal extent). Thus modelling of the seasonal transitions and intraseasonal fluctuations requires an understanding of the fluid mechanics of these three scales of organizations and their interactions. The present paper is an attempt to outline the current state of understanding of these phenomena.

Long-range Coulomb interactions and nanoscale electronic inhomogeneities in correlated oxides

Electronic, magnetic, or structural inhomogeneities ranging in size from nanoscopic to mesoscopic scales seem endemic and are possibly generic to colossal magnetoresistance manganites and other transition metal oxides. They are hence of great current interest and understanding them is of fundamental importance. We show here that an extension, to include long-range Coulomb interactions, of a quantum two-fluid l-b model proposed recently for manganites [Phys. Rev. Lett. 92, 157203 (2004)] leads to an excellent description of such inhomogeneities. In the l-b model two very different kinds of electronic states, one localized and polaronic (l) and the other extended or broad band (b) coexist. For model parameters appropriate to manganites and even within a simple dynamical mean-field theory (DMFT) framework, it describes many of the unusual phenomena seen in manganites, including colossal magnetoresistance (CMR), qualitatively and quantitatively. However, in the absence of long-ranged Coulomb interaction, a system described by such a model would actually phase separate, into macroscopic regions of l and b electrons, respectively. As we show in this paper, in the presence of Coulomb interactions, the macroscopic phase separation gets suppressed and instead nanometer scale regions of polarons interspersed with band electron puddles appear, constituting a kind of quantum Coulomb glass. We characterize the size scales and distribution of the inhomogeneity using computer simulations. For realistic values of the long-range Coulomb interaction parameter V-0, our results for the thresholds for occupancy of the b states are in agreement with, and hence support, the earlier approach mentioned above based on a configuration averaged DMFT treatment which neglects V-0; but the present work has features that cannot be addressed in the DMFT framework. Our work points to an interplay of strong correlations, long-range Coulomb interaction, and dopant ion disorder, all inevitably present in transition metal oxides as the origin of nanoscale inhomogeneities rather than disorder frustrated phase competition as is generally believed. As regards manganites, it argues against explanations for CMR based on disorder frustrated phase separation and for an intrinsic origin of CMR. Based on this, we argue that the observed micrometer (meso) scale inhomogeneities owe their existence to extrinsic causes, e.g., strain due to cracks and defects. We suggest possible experiments to validate our speculation.

Radiative interactions in boundary layers

A detailed description of radiative interactions in laminar compressible boundary layers for moderate Mach numbers is presented by way of asymptotic analysis and supporting solutions. The radiation field is described by the differential approximation. While the asymptotic analysis is valid for large N (the ratio of photon mean free path to molecular mean free path) and arbitrary Boltzmann number, Bo (the ratio of convective heat flux to radiation heat flux), the solutions are obtained for Bo [double less-than sign] 1, the case of strong radiative interactions. The asymptotic analysis shows the existence of an optically thin boundary layer for large N and all Bo. For Bo [double less-than sign] 1, two outer regions are observed — one optically thin (at short distances from the leading edge) and the other optically thick (at large distances from the leading edge). An interesting feature not pointed out in the previous literature is the existence of a wall layer at large distances from the leading edge where convective heat flux can be ignored to the leading order of approximation. The radiation field in all cases can be very well approximated by a one-dimensional description. The solutions have been constructed using the ideas of matched asymptotic expansions by approximate analytical procedures and numerical methods. It is shown that, to the leading order of approximation, the radiation slip method yields exactly the same result as the more complicated matching procedure. Both the cases of linear and nonlinear radiation have been considered, the former being of interest in developing approximate methods which are subsequently generalized to handle the nonlinear problem. Detailed results are presented for both cases.

Determination of favorable inter-particle interactions for formation of substitutionally ordered solid phases from a binary mixture of oppositely charged colloidal suspensions

The solid phase formed by a binary mixture of oppositely charged colloidal particles can be either substitutionally ordered or substitutionally disordered depending on the nature and strength of interactions among the particles. In this work, we use Monte Carlo molecular simulations along with the Gibbs-Duhem integration technique to map out the favorable inter-particle interactions for the formation of substitutionally ordered crystalline phases from a fluid phase. The inter-particle interactions are modeled using the hard core Yukawa potential but the method can be easily extended to other systems of interest. The study obtains a map of interactions depicting regions indicating the type of the crystalline aggregate that forms upon phase transition.

Note: Effect of fluid phase compositions on the formation of substitutionally ordered solid phases from a binary mixture of oppositely charged colloidal suspensions

A binary mixture of oppositely charged colloidal particles can self-assemble into either a substitutionally ordered or substitutionally disordered crystalline phase depending on the nature and strength of interactions among the particles. An earlier study had mapped out favorable inter-particle interactions for the formation of substitutionally ordered crystalline phases from a fluid phase using Monte Carlo molecular simulations along with the Gibbs-Duhem integration technique. In this paper, those studies are extended to determine the effect of fluid phase composition on formation of substitutionally ordered solid phases.

Nonlinear gravity-wave interactions in stratified turbulence

To investigate the dynamics of gravity waves in stratified Boussinesq flows, a model is derived that consists of all three-gravity-wave-mode interactions (the GGG model), excluding interactions involving the vortical mode. The GGG model is a natural extension of weak turbulence theory that accounts for exact three-gravity-wave resonances. The model is examined numerically by means of random, large-scale, high-frequency forcing. An immediate observation is a robust growth of the so-called vertically sheared horizontal flow (VSHF). In addition, there is a forward transfer of energy and equilibration of the nonzero-frequency (sometimes called ``fast'') gravity-wave modes. These results show that gravity-wave-mode interactions by themselves are capable of systematic interscale energy transfer in a stratified fluid. Comparing numerical simulations of the GGG model and the full Boussinesq system, for the range of Froude numbers (Fr) considered (0.05 a parts per thousand currency sign Fr a parts per thousand currency sign 1), in both systems the VSHF is hardest to resolve. When adequately resolved, VSHF growth is more vigorous in the GGG model. Furthermore, a VSHF is observed to form in milder stratification scenarios in the GGG model than the full Boussinesq system. Finally, fully three-dimensional nonzero-frequency gravity-wave modes equilibrate in both systems and their scaling with vertical wavenumber follows similar power-laws. The slopes of the power-laws obtained depend on Fr and approach -2 (from above) at Fr = 0.05, which is the strongest stratification that can be properly resolved with our computational resources.

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.

Anomalous Viscosity Reduction and Hydrodynamic Interactions of Polymeric Nanocolloids in Polymers

One of the central dogmas of fluid physics is the no-slip boundary condition, whose validity has come under scrutiny, especially in the fields of micro and nanofluidics. Although various studies show the violation of the no-slip condition its effect on flow of colloidal particles in viscous media has been rarely explored. Here we report unusually large reduction of effective viscosity experienced by polymeric nano colloids moving through a highly viscous and confined polymer, well above its glass transition temperature. The extent of reduction in effective interface viscosity increases with decreasing temperature and polymer film thickness. Concomitant with the reduction in effective viscosity we also observe apparent divergence of the wave vector dependent hydrodynamic interaction function of these colloids with an anomalous power law exponent of similar to 2 at the lowest temperatures and film thickness studied. Such strong hydrodynamic interactions are not expected for polymeric colloidal motion in polymer melts. We suggest hydrodynamics, especially slip present at the colloid-polymer interface which determines the observed reduction in interface viscosity and presence of strong hydrodynamic interactions.

Intracellular concentrations of Ca2+ modulate the strength of signal and alter the outcomes of cytotoxic T-lymphocyte antigen-4 (CD152)-CD80/CD86 interactions in CD4(+) T lymphocytes

The costimulatory receptors CD28 and cytotoxic T-lymphocyte antigen (CTLA)-4 and their ligands, CD80 and CD86, are expressed on T lymphocytes; however, their functional roles during T cell-T cell interactions are not well known. The consequences of blocking CTLA-4-CD80/CD86 interactions on purified mouse CD4(+) T cells were studied in the context of the strength of signal (SOS). CD4(+) T cells were activated with phorbol 12-myristate 13-acetate (PMA) and different concentrations of a Ca2+ ionophore, Ionomycin (I), or a sarcoplasmic Ca2+ ATPase inhibitor, Thapsigargin (TG). Increasing concentrations of I or TG increased the amount of interleukin (IL)-2, reflecting the conversion of a low to a high SOS. During activation with PMA and low amounts of I, intracellular concentrations of calcium ([Ca2+](i)) were greatly reduced upon CTLA-4-CD80/CD86 blockade. Further experiments demonstrated that CTLA-4-CD80/CD86 interactions reduced cell cycling upon activation with PMA and high amounts of I or TG (high SOS) but the opposite occurred with PMA and low amounts of I or TG (low SOS). These results were confirmed by surface T-cell receptor (TCR)-CD3 signalling using a low SOS, for example soluble anti-CD3, or a high SOS, for example plate-bound anti-CD3. Also, CTLA-4-CD80/CD86 interactions enhanced the generation of reactive oxygen species (ROS). Studies with catalase revealed that H2O2 was required for IL-2 production and cell cycle progression during activation with a low SOS. However, the high amounts of ROS produced during activation with a high SOS reduced cell cycle progression. Taken together, these results indicate that [Ca2+](i) and ROS play important roles in the modulation of T-cell responses by CTLA-4-CD80/CD86 interactions.

Spectroscopic studies of the interactions of the pyrethroid insecticide fenvalerate with gramicidin

Fenvalerate is a pyrethroid insecticide which interacts with ionic channels. Using circular dichroism technique we have studied the interaction of fenvalerate with gramicidin, a model channel peptide which transports ions. In most organic solvents, gramicidin exists as a double helix except in trifluoroethanol where it exists as a channel forming single stranded beta(6.3) helical monomer. In model lipid membranes, under certain experimental conditions, gramicidin exists as a channel forming single stranded beta(6.3) helical dimer. Our results show that fenvalerate interacts more with the single stranded beta(6.3) helical monomer or dimer than with the double helical form of gramicidin. This was further confirmed by an increase in the rate of gramicidin mediated proton transport in liposomes by fenvalerate, using the pH sensitive fluorophore, pyranine.

Unsteady-Flow Of A Viscous-Fluid Between 2 Parallel Disks With A Time-Varying Gap Width And A Magnetic-Field

The unsteady incompressible viscous fluid flow between two parallel infinite disks which are located at a distance h(t*) at time t* has been studied. The upper disk moves towards the lower disk with velocity h'(t*). The lower disk is porous and rotates with angular velocity Omega(t*). A magnetic field B(t*) is applied perpendicular to the two disks. It has been found that the governing Navier-Stokes equations reduce to a set of ordinary differential equations if h(t*), a(t*) and B(t*) vary with time t* in a particular manner, i.e. h(t*) = H(1 - alpha t*)(1/2), Omega(t*) = Omega(0)(1 - alpha t*)(-1), B(t*) = B-0(1 - alpha t*)(-1/2). These ordinary differential equations have been solved numerically using a shooting method. For small Reynolds numbers, analytical solutions have been obtained using a regular perturbation technique. The effects of squeeze Reynolds numbers, Hartmann number and rotation of the disk on the flow pattern, normal force or load and torque have been studied in detail

Experiments on a plume with off-source heating: Implications for cloud fluid dynamics

We report here on a series of laboratory experiments on plumes, undertaken with the object of simulating the effect of the heat release that occurs in clouds on condensation of water vapor. The experimental technique used for this purpose relies on ohmic heating generated in an electrically conducting plume fluid subjected to a suitable alternating voltage across specified axial stations in the plume flow [Bhat et al., 1989]. The present series of experiments achieves a value of the Richardson number that is toward the lower end of the range that characteristics cumulus clouds. It is found that the buoyancy enhancement due to heating disrupts the eddy structures in the flow and reduces the dilution owing to entrainment of ambient fluid that would otherwise have occurred in the central region of the plume. Heating also reduces the spread rate of the plume, but as it accelerates the flow as well, the overall specific mass flux in the plume does not show a very significant change at the heat input employed in the experiment. However, there is some indication that the entrainment rate (proportional to the streamwise derivative of the mass flux) is slightly higher immediately after heat injection and slightly lower farther downstream. The measurements support a previous proposal for a cloud scenario [Bhat and Narasimha, 1996] and demonstrate how fresh insights into certain aspects of the fluid dynamics of clouds may be derived from the experimental techniques employed here.

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.