Anomalous behavior of hydrodynamic modes in the two dimensional shear flow of a granular material

The growth rates of the hydrodynamic modes in the homogeneous sheared state of a granular material are determined by solving the Boltzmann equation. The steady velocity distribution is considered to be the product of the Maxwell Boltzmann distribution and a Hermite polynomial expansion in the velocity components; this form is inserted into them Boltzmann equation and solved to obtain the coeificients of the terms in the expansion. The solution is obtained using an expansion in the parameter epsilon =(1 - e)(1/2), and terms correct to epsilon(4) are retained to obtain an approximate solution; the error due to the neglect of higher terms is estimated at about 5% for e = 0.7. A small perturbation is placed on the distribution function in the form of a Hermite polynomial expansion for the velocity variations and a Fourier expansion in the spatial coordinates: this is inserted into the Boltzmann equation and the growth rate of the Fourier modes is determined. It is found that in the hydrodynamic limit, the growth rates of the hydrodynamic modes in the flow direction have unusual characteristics. The growth rate of the momentum diffusion mode is positive, indicating that density variations are unstable in the limit k--> 0, and the growth rate increases proportional to kslash} k kslash}(2/3) in the limit k --> 0 (in contrast to the k(2) increase in elastic systems), where k is the wave vector in the flow direction. The real and imaginary parts of the growth rate corresponding to the propagating also increase proportional to kslash k kslash(2/3) (in contrast to the k(2) and k increase in elastic systems). The energy mode is damped due to inelastic collisions between particles. The scaling of the growth rates of the hydrodynamic modes with the wave vector I in the gradient direction is similar to that in elastic systems. (C) 2000 Elsevier Science B.V. All rights reserved.

Novel organic porous solids with channel and layered structures from 1,3,5-triazine-2,4,6-triaminehexaacetic acid and its calcium salt

Single crystals of a symmetrically substituted molecule, 1,3,5-triazine-2,4,6-triaminehexaacetic acid, (TTHA) and its Ca2+ salt have been synthesized, the analysis of which reveals the existence of novel channel type cavities and helical packing organizations in the crystals.

Effects of electron correlations in conjugated organic molecules and solids

The discrepancies between the non-interacting models and experimental results for conjugated systems is highlighted in this brief review. The interacting model hamiltonians correctly give the forbidden singlet state below the optical gap in polyenes and also explain both the nonvanishing optical gap in polyacetylenes and the vanishing optical gap in symmetric cyanine dyes. The negative spin densities in polyene radicals is also understood in terms of a correlated picture. The role of electron-electron interactions in other strongly correlated systems, such as polydiacetylene and mixed and segregated stack charge transfer solids, are also briefly discussed.

Correlations and the Ring-Kinetic Equation in Dense Sheared Granular Flows

A formal way of deriving fluctuation-correlation relations in dense sheared granular media, starting with the Enskog approximation for the collision integral in the Chapman-Enskog theory, is discussed. The correlation correction to the viscosity is obtained using the ring-kinetic equation, in terms of the correlations in the hydrodynamic modes of the linearised Enskog equation. It is shown that the Green-Kubo formula for the shear viscosity emerges from the two-body correlation function obtained from the ring-kinetic equation.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

The stress in a slowly sheared granular column

We report measurements of the wall stress in a granular material sheared in a cylindrical Couette cell, as a function of the distance from the free surface. Our results shows that when the material is static, all components of the stress saturate to constant values within a short distance from the free surface, in conformity with earlier experiments and theoretical predictions. When the material is sheared by rotating the inner cylinder at a constant rate, the stresses are remarkably altered. The radial normal stress does not saturate, and increases even more rapidly with depth than the linear hydrostatic pressure profile. The axial shear stress changes sign on shearing, and its magnitude increases with depth. These results are discussed in the context of the predictions of the classical and Cosserat plasticity theories.

Kinetic theory for the density plateau in the granular flow down an inclined plane

The striking lack of observable variation of the volume fraction with height in the center of a granular flow down an inclined plane is analysed using constitutive relations obtained from kinetic theory. It is shown that the rate of conduction in the granular energy balance equation is O(delta(2)) smaller than the rate of production of energy due to mean shear and the rate of dissipation due to inelastic collisions, where the small parameter delta = (d/(1 - e(n))H-1/2), d is the particle diameter, en is the normal coefficient of restitution and H is the thickness of the flowing layer. This implies that the volume fraction is a constant in the leading approximation in an asymptotic analysis in small delta. Numerical estimates of both the parameter delta and its pre-factor are obtained to show that the lack of observable variation of the volume fraction with height can be explained by constitutive relations obtained from kinetic theory.

Importance of local heating in the compaction of thermally sensitive solids

Local heating is an important parameter in compaction of thermally sensitive solids since local hot spots could conceivably raise the temperature of the system by several hundred degrees. To understand the importance of local hot spots, 20 g of ammonium perchlorate (AP) and potassium perchlorate (KP) were pressed together at 1500 kg cm−2 for 5 min. The surface structural examination of the compact revealed a secondary phase at the interfaces between the grains in the compositions ranging from 50% AP-50% KP to 10% AP-90% KP. The observation of the secondary phase only at the interfaces has been attributed to the short life times of temperature pulses present during the compaction. An interesting aspect of the investigation is the formation of a series of solid solutions of AP-KP.

Steady compressible flow of granular materials through a wedge-shaped hopper: The smooth wall, radial gravity problem

A continuum model based on the critical state theory of soil mechanics is used to generate stress and density profiles, and to compute discharge velocities for the plane flow of cohesionless materials. Two types of yield loci are employed, namely, a yield locus with a corner, and a smooth yield locus. The yield locus with a corner leads to computational difficulties. For the smooth yield locus, results are found to be relatively insensitive to the shape of the yield locus, the location of the upper traction-free surface and the density specified on this surface. This insensitivity arises from the existence of asymptotic stress and density fields, to which the solution tends to converge on moving down the hopper. Numerical and approximate analytical solutions are obtained for these fields and the latter is used to derive an expression for the discharge velocity. This relation predicts discharge velocities to within 13% of the exact (numerical) values. While the assumption of incompressibility has been frequently used in the literature, it is shown here that in some cases, this leads to discharge velocities which are significantly higher than those obtained by the incorporation of density variation.

Use of metastable equilibria for determination of Gibbs energy of solids

Attempts are made to measure activities of both components of a binary alloy (A�B) at 650 K using a solid-state galvanic cell incorporating a new composite solid electrolyte. Since the ionic conductivity of the composite solid electrolyte is three orders of magnitude higher than that of pure CaF2, the cell can be operated at lower temperatures. The alloy phase is equilibrated in separate experiments with flourides of each component and fluorine potential is measured. The mixture of the alloy (A�B) and the fluoride of the more reactive component (BF2) is stable, while (A�B) + AF2 mixture is metastable, Factors governing the possible use of metastable equilibria have been elucidated in this study. In the Co�Ni system, where the difference in Gibbs energies of formation of the fluorides is 21.4 kJ/mol, emf of the cell with metastable phases at the electrode is constant for periods ranging from 90 to 160 ks depending on alloy composition. Subsequently, the emf decreases because of the onset of the displacement reaction. In the Ni�Mn system, measurement of the activity of Ni using metastable equilibria is not fully successful at 650 K because of the large driving force for the displacement reaction (208.8 kJ/mol). Critical factors in the application of metastable equilibria are the driving force for displacement reaction and diffusion coefficients in both the alloy and fluoride solid solution.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Models for spin-state transitions in solids

The model for spin-state transitions described by Bari and Sivardiere (1972) is static and can be solved exactly even when the dynamics of the lattice are included; the dynamic model does not, however, show any phase transition. A coupling between the octahedra, on the other hand, leads to a phase transition in the dynamical two-sublattice displacement model. A coupling of the spin states to the cube of the sublattice displacement leads to a first-order phase transition. The most reasonable model appears to be a two-phonon model in which an ion-cage mode mixes the spin states, while a breathing mode couples to the spin states without mixing. This model explains the non-zero population of high-spin states at low temperatures, temperature-dependent variations in the inverse susceptibility and the spin-state population ratio, as well as the structural phase transitions accompanying spin-state transitions found in some systems.

A thermodynamic protocol for explaining subcriticality and other subtle features in self-deflagrating solids

A novel universal approach to understand the self-deflagration in solids has been attempted by using basic thermodynamic equation of partial differentiation, where burning mte depends on the initial temperature and pressure of the system. Self-deflagrating solids are rare and are reported only in few compounds like ammonium perchlorate (AP), polystyrene peroxide and tetrazole. This approach has led us to understand the unique characteristics of AP, viz. the existence of low pressure deflagration limit (LPL 20 atm), hitherto not understood sufficiently. This analysis infers that the overall surface activation energy comprises of two components governed by the condensed phase and gas phase processes. The most attractive feature of the model is the identification of a new subcritical regime I' below LPL where AP does not burn. The model is aptly supported by the thermochemical computations and temperature-profile analyses of the combustion train. The thermodynamic model is further corroborated from the kinetic analysis of the high pressure (1-30 atm) DTA thermograms which affords distinct empirical decomposition rate laws in regimes I' and 1 (20-60 atm). Using Fourier-Kirchoff one dimensional heat transfer differential equation, the phase transition thickness and the melt-layer thickness have been computed which conform to the experimental data.