Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Application of partition method to vibration problems of plates

The “partition method” or “sub-domain method” consists of expressing the solution of a governing differential equation, partial or ordinary, in terms of functions which satisfy the boundary conditions and setting to zero the error in the differential equation integrated over each of the sub-domains into which the given domain is partitioned. In this paper, the use of this method in eigenvalue problems with particular reference to vibration of plates is investigated. The deflection of the plate is expressed in terms of polynomials satisfying the boundary conditions completely. Setting the integrated error in each of the subdomains to zero results in a set of simultaneous, linear, homogeneous, algebraic equations in the undetermined coefficients of the deflection series. The algebraic eigenvalue problem is then solved for eigenvalues and eigenvectors. Convergence is examined in a few typical cases and is found to be satisfactory. The results obtained are compared with existing results based on other methods and are found to be in very good agreement.

The mechanism of manganese electrodeposition

The mechanism of manganese electrodeposition from a sulphate bath on to a stainless-steel substrate has been studied by using current efficiency data to resolve the totali-E curves. A simple, two-step electron transfer mechanism:is proposed to explain the following experimentally obtained parameters: cathodic and anodic transfer coefficients, reaction order and stoichiometric number. The mechanism also explains the effect of pH oni o,Mn and on the corrosion currents.

Dendritic structures produced on solidification of multicomponent aqueous solutions

Dendrite structures of ice produced on undirectional solidification of ternary and quaternary aqueous solutions have been studied. Upon freezing, solutions containing more than one solute produce plate-shaped dendrites of ice. The spacing between dendrites increase linearly with the distance from the chill surface and the square root of local solidification time (or square root of inverse freezing rate) for any fixed composition. For fixed freezing conditions, the dendrite spacings from multicomponent aqueous solutions were a function of the concentrations and diffusion coefficients of the individual solutes. The dendrite spacing produced by freezing of a solution was changed by the addition of a solute different from those already present. If the main diffusion coefficient of the added solute is higher than that of solutes already present, the dendrite spacing is increased and vice versa. The dendrite spacing in multi-component systems increases with the total solute concentration if the constituent solutes are present in equal amounts. The dendrite spacing obtained on freezing of these dilute multicomponent solutions can be expressed by regression equations of the type Image Full-size image (2K) where L is the dendrite spacing in microns, C1, C2 and C3 are concentrations of individual solutes, Θf is the total freezing time and A1 −A8 are constants. A Yates analysis of the dendrite spacings in a factorial design of quaternary solutions indicates that there are strong interactions between individual solutes in regard to their effect on the dendrite spacings. A mass transport analysis has been used to calculate the interdendritic supersaturation ΔC of the individual solutes, the supercooling in the interdendritic liquid ΔT, and the transverse growth velocity of the dendrites, VT. In ternary solutions if two solutes are present in equal amount the supersaturation of the solute with higher main diffusion coefficient is lower, and vice versa. If a solute with higher main diffusion coefficient is added to a binary solution, the interface growth velocity, the interdendritic supersaturation of the base solute and the interdendritic supercooling increase with the quantity of solute added.

Identification of a class of non-linear systems

This paper considers the on-line identification of a non-linear system in terms of a Hammerstein model, with a zero-memory non-linear gain followed by a linear system. The linear part is represented by a Laguerre expansion of its impulse response and the non-linear part by a polynomial. The identification procedure involves determination of the coefficients of the Laguerre expansion of correlation functions and an iterative adjustment of the parameters of the non-linear gain by gradient methods. The method is applicable to situations involving a wide class of input signals. Even in the presence of additive correlated noise, satisfactory performance is achieved with the variance of the error converging to a value close to the variance of the noise. Digital computer simulation establishes the practicability of the scheme in different situations.

Growth of ionization currents in nitrogen

Townsend's primary and secondary ionization coefficients α/p and γ were determined in nitrogen over a wide range of E/p (100-1000 V cm−1 Torr−1) and p (0·4 to 12 Torr at 0 °C) using the pressure variation technique. This technique, along with the Gosseries method of evaluation of ionization coefficients, seems to be more suitable at higher values of E/p, since the errors in these coefficients could be minimized by a suitable selection of p and d, thus eliminating the non-equilibrium ionization condition.

Photon-Induced Electron Pairing

The earlier work on the possibility of interband electron pairing in the presence of a strong radiation field has been further extended. Some additional terms, neglected earlier, have been taken into account and generalized to a situation where the electron-phonon coupling coefficients for the two conduction bands (valleys) are different. It is found that the pairing interaction is attractive and the strength depends on the photon density.

Application of partition method to vibration problems of plates

The “partition method” or “sub-domain method” consists of expressing the solution of a governing differential equation, partial or ordinary, in terms of functions which satisfy the boundary conditions and setting to zero the error in the differential equation integrated over each of the sub-domains into which the given domain is partitioned. In this paper, the use of this method in eigenvalue problems with particular reference to vibration of plates is investigated. The deflection of the plate is expressed in terms of polynomials satisfying the boundary conditions completely. Setting the integrated error in each of the subdomains to zero results in a set of simultaneous, linear, homogeneous, algebraic equations in the undetermined coefficients of the deflection series. The algebraic eigenvalue problem is then solved for eigenvalues and eigenvectors. Convergence is examined in a few typical cases and is found to be satisfactory. The results obtained are compared with existing results based on other methods and are found to be in very good agreement.

Cooperativity in the inhibition of oxidative phosphorylation by chlorophenoxyisobutyrate

The kinetics of inhibition of oxidative phosphorylation by the antihypercholes-terolaemic compound p-chlorophenoxyisobutyrate reveal cooperativity characteristic of allosteric interactions. Hill plots and Dixon plots give clear indication that the compound interferes with two distinct steps in the energy-transfer pathway. The values of interaction coefficients calculated from the Hill plots were two and four in the direction of ATP synthesis and one and two in the reverse direction. This could mean either that the pathways of synthesis and breakdown of ATP are different, or that if the pathways are the same, only half the inhibitor-binding sites function in the reverse direction.

On the bottleneck linear programming problem

In this paper, the results on primal methods for Bottleneck Linear Programming (BLP) problem are briefly surveyed, the primal method is presented and the degenerate case related to Bottleneck Transportation Problem (BTP) is explicitly considered. The algorithm is based on the idea of using auxiliary coefficients as is done by Garfinkel and Rao [6]. The modification presented for the BTP rectifies the defect in Hammer's method in the case of degenerate basic feasible solution. Illustrative numerical examples are also given.

Stress-optical dispersion in rubidium chloride and rubidium bromide crystals

The stress-optical coefficients C = (n3/2) (q11−q12) and C′ = (n3/2)q44 of RbCl and RbBr crystals have been measured at room temperature (26°C) over the wave length range 5750-2500 A.

Wave analysis of sound decay in rectangular rooms

The decay of sound in a rectangular room is analyzed for various boundary conditions on one of its walls. It is shown that the decay of the sound-intensity level is in general nonlinear. But for specific areas and impedances of the material it is possible to obtain a linear initial decay. It is also shown that the coefficients derived from the initial decay rates neither correspond to the predictions of Sabine's or Eyring's geometrical theories nor to the normal coefficients of Morse's wave theory. The dependence of the coefficients on the area of the material is discussed. The influence of the real and the imaginary parts of the specific acoustic impedance of the material on the coefficients is also discussed. Finally, the existence of a linear initial decay corresponding to the decay of a diffuse field in the case of a highly absorbing material partially covering a wall is explained on the basis of modal coupling.

Growth of ionization currents in nitrogen

Townsend's primary and secondary ionization coefficients α/p and γ were determined in nitrogen over a wide range of E/p (100-1000 V cm−1 Torr−1) and p (0·4 to 12 Torr at 0 °C) using the pressure variation technique. This technique, along with the Gosseries method of evaluation of ionization coefficients, seems to be more suitable at higher values of E/p, since the errors in these coefficients could be minimized by a suitable selection of p and d, thus eliminating the non-equilibrium ionization condition.

Energy losses at pipe trifurcations

In the case of pipe trifurcation, previous observations report negative energy losses in the centre branch. This causes an anomaly, because there should not be any negative energy loss due to conservation of energy principle. Earlier investigators have suggested that this may be due to the non-inclusion of kinetic energy coefficient (a) in the computations of energy losses without any experimental evidence. In the present work, through experimentally determined velocity profiles, energy loss coefficients have been evaluated. It has been found that with the inclusion of a in the computations of energy loss, there is no negative energy loss in the centre branch.