Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Phase transitions in triammonium hydrogen disulphate: Crystal structure of an intermediate stable phase at-90 degrees C

Triammonium hydrogen disulphate, (NH4)(3)H(SO4)(2), belongs to the family of crystal structures M3H(XO4)(2) (with M = NH4, K, Rb, Cs, and X = S, Se) which display super protonic phases at elevated temperatures, while at room temperature these are relatively poor proton conductors. The crystal structure of triammonium hydrogen disulphate has been determined by X-ray diffraction at -90 degrees C and the variation in the characteristics of the hydrogen bond is discussed in comparison with that of the structures at -110 degrees C and room temperature. It is concluded that the mechanics involving the proton migration in such systems is realised in terms of the variations in the hydrogen bond features with temperature.

Residence time distribution for a class of Gaussian Markov processes

We study the distribution of residence time or equivalently that of "mean magnetization" for a family of Gaussian Markov processes indexed by a positive parameter alpha. The persistence exponent for these processes is simply given by theta=alpha but the residence time distribution is nontrivial. The shape of this distribution undergoes a qualitative change as theta increases, indicating a sharp change in the ergodic properties of the process. We develop two alternate methods to calculate exactly but recursively the moments of the distribution for arbitrary alpha. For some special values of alpha, we obtain closed form expressions of the distribution function. [S1063-651X(99)03306-1].

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

Accelerated weathering of distribution class polymeric insulators

This paper deals with the long-term accelerated weathering of 11 kV polymeric insulators for 25000 h. Polymeric insulators were continuously subjected to accelerated weathering in a specially designed multistress-aging chamber under UV radiation, temperature and electric stress. Chemical, physical and electrical changes due to degradation have been assessed using various techniques. Some of the interesting results observed indicate that there is a significant reduction in the content of low molecular weight molecules, hydrophobicity was dynamic in nature and there is a significant increase in the surface roughness and oxidation levels with respect to the duration of the weathering.

Continuous distribution kinetics for the thermal degradation of LDPE in solution

Polymer degradation in solution has several advantages over melt pyrolysis, The degradation of low-density polyethylene (LDPE) occurs at much lower temperatures in solution (280-360degreesC) than in conventional melt pyrolysis (400-450degreesC). The thermal degradation kinetics of LDPE in solution was investigated in this work. LDPE was dissolved in liquid paraffin and degraded for 3 h at various temperatures (280-360degreesC). Samples were taken at specific times and analyzed with high-pressure liquid chromatography/gel permeation chromatography for the molecular weight distribution (MWD), The time evolution of the MWD was modeled with continuous distribution kinetics. Data indicated that LDPE followed random-chain-scission degradation. The rapid initial drop in molecular weight, observed up to 45 min, was attributed to the presence of weak links in the polymer. The rate coefficients for the breakage of weak and strong links were determined, and the corresponding average activation energies were calculated to be 88 and 24 kJ/mol, respectively. (C) 2002 John Wiley Sons, Inc.

Woody plant seedling distribution under invasive Lantana camara thickets in a dry-forest plot in Mudumalai, southern India

Lantana camara, a shrub of Central and South American origin, has become invasive across dry forests worldwide. The effect of the thicket-forming habit of L. camara as a dispersal and recruitment barrier in a community of native woody seedlings was examined in a 50-ha permanent plot located in the seasonally dry forest of Mudumalai, southern India. Sixty 100-m(2) plots were enumerated for native woody seedlings between 10-100 cm in height. Of these, 30 plots had no L. camara thickets, while the other 30 had dense thickets. The frequency of occurrence and abundance of seedlings were modelled as a function of dispersal mode (mammal, bird or mechanical) and affinities to forest habitats (dry forest, moist forest or ubiquitous) as well as presence or absence of dense L. camara thickets. Furthermore, frequency of occurrence and abundance of individual species were also compared between thickets and no L. camara. At the community level, L. camara density, dispersal mode and forest habitat affinities of species determined both frequency of occurrence and abundance of seedlings, with the abundance of dry-forest mammal-dispersed species and ubiquitous mechanically dispersed species being significantly lower under L. camara thickets. Phyllanthus emblica and Kydia calycina were found to be significantly less abundant under L. camara, whereas most other species were not affected by the presence of thickets. It was inferred that, by affecting the establishment of native tree seedlings, L. camara thickets could eventually alter the community composition of such forests.

Simulation and experimental study on compositional evolution of Li-Co in LiCoO(2) thin films during sputter deposition

The compositional evolution in sputter deposited LiCoO(2) thin films is influenced by process parameters involved during deposition. The electrochemical performance of these films strongly depends on their microstructure, preferential orientation and stoichiometry. The transport process of sputtered Li and Co atoms from the LiCoO(2) target to the substrate, through Ar plasma in a planar magnetron configuration, was investigated based on the Monte Carlo technique. The effect of sputtering gas pressure and the substrate-target distance (d(st)) on Li/Co ratio, as well as, energy and angular distribution of sputtered atoms on the substrate were examined. Stable Li/Co ratios have been obtained at 5 Pa pressure and d(st) in the range 5 11 cm. The kinetic energy and incident angular distribution of Li and Co atoms reaching the substrate have been found to be dependent on sputtering pressure. Simulations were extended to predict compositional variations in films prepared at various process conditions. These results were compared with the composition of films determined experimentally using x-ray photoelectron spectroscopy (XPS). Li/Co ratio calculated using XPS was in moderate agreement with that of the simulated value. The measured film thickness followed the same trend as predicted by simulation. These studies are shown to be useful in understanding the complexities in multicomponent sputtering. (C) 2011 American Institute of Physics. doi:10.1063/1.3597829]

Continuous distribution kinetics for ultrasonic degradation of poly(methyl methacrylate)

Ultrasonic degradation of poly(methyl methacrylate) (PMMA) was carried out in several solvents and some mixtures of solvents. The time evolution of molecular weight distribution (MWD), determined by gel permeation chromatography, is analysed by continuous distribution kinetics. The rate coefficients for polymer degradation are determined for each solvent. The variation of rate coefficients is correlated with the vapour pressure of the solvent, kinematic viscosity of the solution and solvent-polymer interaction parameters. The vapour pressure and the kinematic viscosity of the solution are found to be more critical than other parameters (such as the Huggins and Flory-Huggins constants) in determining the degradation rates. (C) 2001 Society of Chemical Industry.

Wear of metals: a consequence of stable/unstable material response

Wear of metals in dry sliding is dictated by the material response to traction. This is demonstrated by considering the wear of aluminium and titanium alloys. In a regime of stable homogeneous deformation the material approaching the surface from the bulk passes through microprocessing zones of flow, fracture, comminution and compaction to generate a protective tribofilm that retains the interaction in the mild wear regime. If the response leads to microstructural instabilities such as adiabatic shear bands, the near-surface zone consists of stacks of 500 nm layers situated parallel to the sliding direction. Microcracks are generated below the surface to propagate normally away from the surface though microvoids situated in the layers, until it reaches a depth of 10-20 mum. A rectangular laminate debris consisting of a 20-40 layer stack is produced, The wear in this mode is severe.

An adaptive cascade compensator for a wide class of stable minimum or nonminimum phase plants through time moment estimation

Model Reference Adaptive Control (MRAC) of a wide repertoire of stable Linear Time Invariant (LTI) systems is addressed here. Even an upper bound on the order of the finite-dimensional system is unavailable. Further, the unknown plant is permitted to have both minimum phase and nonminimum phase zeros. Model following with reference to a completely specified reference model excited by a class of piecewise continuous bounded signals is the goal. The problem is approached by taking recourse to the time moments representation of an LTI system. The treatment here is confined to Single-Input Single-Output (SISO) systems. The adaptive controller is built upon an on-line scheme for time moment estimation of a system given no more than its input and output. As a first step, a cascade compensator is devised. The primary contribution lies in developing a unified framework to eventually address with more finesse the problem of adaptive control of a large family of plants allowed to be minimum or nonminimum phase. Thus, the scheme presented in this paper is confined to lay the basis for more refined compensators-cascade, feedback and both-initially for SISO systems and progressively for Multi-Input Multi-Output (MIMO) systems. Simulations are presented.