Fast decay of the velocity autocorrelation function in dense shear flow of inelastic hard spheres

We find in complementary experiments and event-driven simulations of sheared inelastic hard spheres that the velocity autocorrelation function psi(t) decays much faster than t(-3/2) obtained for a fluid of elastic spheres at equilibrium. Particle displacements are measured in experiments inside a gravity-driven flow sheared by a rough wall. The average packing fraction obtained in the experiments is 0.59, and the packing fraction in the simulations is varied between 0.5 and 0.59. The motion is observed to be diffusive over long times except in experiments where there is layering of particles parallel to boundaries, and diffusion is inhibited between layers. Regardless, a rapid decay of psi(t) is observed, indicating that this is a feature of the sheared dissipative fluid, and is independent of the details of the relative particle arrangements. An important implication of our study is that the non-analytic contribution to the shear stress may not be present in a sheared inelastic fluid, leading to a wider range of applicability of kinetic theory approaches to dense granular matter.

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.

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.

Immunization of male bonnet monkeys (M-radiata) with a recombinant FSH receptor preparation affects testicular function and fertility

Immunization of proven fertile adult male monkeys (n = 3) with a recombinant FSH receptor protein preparation (oFSHR-P) (representing amino acids 1-134 of the extracellular domain of the receptor Mr similar to 15KDa) resulted in production of receptor blocking antibodies. The ability of the antibody to bind a particulate FSH receptor preparation and receptors in intact granulosa cells was markedly (by 30-80%) inhibited by FSH. Serum T levels and LH receptor function following immunization remained unchanged. The immunized monkeys showed a 50% reduction (p<0.001) in transformation of spermatogonia(2C) to primary spermatocytes (4C) as determined by flow cytometry and the 4C:2C ratio showed a correlative change (R 0.81, p<0.0007) with reduction in fertility index (sperm counts X motility score). Breeding studies indicated that monkeys became infertile between 242-368 days of immunization when the fertility index was in the range of 123+/-76 to 354+/-42 (compared to a value of 1602+/-384 on day 0). As the effects observed ate near identical to that seen following immunization with FSH it is suggestive that oFSHR-P can substitute for FSH in the development of a contraceptive vaccine.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

A hydrodynamical study of the flow in renal tubules

The hydrodynamical problem of flow in proximal renal tubule is investigated by considering axisymmetric flow of a viscous, incompressible fluid through a long narrow tube of varying cross-section with reabsorption at the wall. Two cases for reabsorption have been studied (i) when the bulk flow,Q, decays exponentially with the axial distancex, and (ii) whenQ is an arbitrary function ofx such thatQ-Q 0 can be expressed as a Fourier integral (whereQ 0 is the flux atx=0). The analytic expressions for flow variables have been obtained by applying perturbation method in terms of wall parameter ε. The effects of ε on pressure drop across the tube, radial velocity and wall shear have been studied in the case of exponentially decaying bulk flow and it has been found that the results are in agreement with the existing ones for the renal tubules.

Changes in testicular function following specific deprivation of LH in the adult male rabbit

Sexually mature male rabbits actively immunized against highly purified ovine LH (oLH) were used as a model system to study the effects of endogenous LH deprivation (and therefore testosterone) on spermatogenesis as well as pituitary FSH secretion. Immunization against oLH generated antibody titres capable of cross-reacting and neutralizing rabbit LH and this resulted in a significant reduction (P<0.01) in serum testosterone levels by 2-4 weeks of immunization. A significant increase in circulating FSH concentration (from a basal level of similar to 1 ng to 60-100 ng/ml; P<0.01) was observed within 4-6 weeks of immunization, perhaps a consequence of the negative feedback effect of the lack of testosterone. The effect of LH deprivation on spermatogenesis assessed by DNA flow cytometry and histological analyses of testicular biopsy tissue revealed that lack of testosterone primarily results in a rapid reduction and complete absence of round (1C) and elongated (HC) spermatids. The immediate effect of LH/testosterone deprivation thus appears to be at the step of meiotic transformation of primary spermatocytes (4C) to 1C. A significant reduction (>80%; P<0.01) in the 4C population and a relative accumulation (>90%; P<0.01) in spermatogonia (2C) was also observed, suggesting a need for testosterone during the transformation of 2C to 1C. In all but one of the rabbits, both qualitative and quantitative recovery in spermatogenesis occurred during the recovery phase, even at a time when only a marginal increase in serum testosterone (compared with the preimmunization) levels was observed as a result of a rapid decline in the cross-reactive antibody titres. These results clearly show that LH/testosterone deprivation in addition to primarily affecting the meiotic step also regulates the conversion of 2C to 4C during spermatogenesis.

Unsteady Natural Convection Flow over a Heated Cylinder Buried in a Fluid Saturated Porous Medium

Transient natural convection flow on a heated cylinder buried in a semi-infinite liquid-saturated porous medium has been studied. The unsteadiness in the problem arises due to the cylinder which is heated (cooled) suddenly and then maintained at that temperature. The coupled partial differential equations governing the flow and heat transfer are cast into stream function-temperature formulation, and the solutions are obtained from the initial time to the time when steady state is reached. The heat transfer is found to change significantly with increasing time in a small time interval immediately after the start of the impulsive change, and steady state is reached after some time. The average Nusselt number is found to increase with Rayleigh number When the surface of the cylinder is suddenly cooled, there is a change in the direction of the heat transfer in a small time interval immediately after the start of the impulsive change in the surface temperature;however when the surface temperature is suddenly increased, no such phenomenon is observed.

Optimal multicommodity flow through the complete graph with random edge capacities

We consider a multicommodity flow problem on a complete graph whose edges have random, independent, and identically distributed capacities. We show that, as the number of nodes tends to infinity, the maximumutility, given by the average of a concave function of each commodity How, has an almost-sure limit. Furthermore, the asymptotically optimal flow uses only direct and two-hop paths, and can be obtained in a distributed manner.

Heat transfer in rarefied couette flow on the basis of discrete ordinate method

The method of discrete ordinates, in conjunction with the modified "half-range" quadrature, is applied to the study of heat transfer in rarefied gas flows. Analytic expressions for the reduced distribution function, the macroscopic temperature profile and the heat flux are obtained in the general n-th approximation. The results for temperature profile and heat flux are in sufficiently good accord both with the results of the previous investigators and with the experimental data.

Potential functions for hydrogen bond interactions - III. Empirical Potential Function for the Peptide N-H...O=C Hydrogen Bond

A careful comparison of the distribution in the (R, θ)-plane of all NH ... O hydrogen bonds with that for bonds between neutral NH and neutral C=O groups indicated that the latter has a larger mean R and a wider range of θ and that the distribution was also broader than for the average case. Therefore, the potential function developed earlier for an average NH ... O hydrogen bond was modified to suit the peptide case. A three-parameter expression of the form {Mathematical expression}, with △ = R - Rmin, was found to be satisfactory. By comparing the theoretically expected distribution in R and θ with observed data (although limited), the best values were found to be p1 = 25, p3 = - 2 and q1 = 1 × 10-3, with Rmin = 2·95 Å and Vmin = - 4·5 kcal/mole. The procedure for obtaining a smooth transition from Vhb to the non-bonded potential Vnb for large R and θ is described, along with a flow chart useful for programming the formulae. Calculated values of ΔH, the enthalpy of formation of the hydrogen bond, using this function are in reasonable agreement with observation. When the atoms involved in the hydrogen bond occur in a five-membered ring as in the sequence[Figure not available: see fulltext.] a different formula for the potential function is needed, which is of the form Vhb = Vmin +p1△2 +q1x2 where x = θ - 50° for θ ≥ 50°, with p1 = 15, q1 = 0·002, Rmin = 2· Å and Vmin = - 2·5 kcal/mole. © 1971 Indian Academy of Sciences.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

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.

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.