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.

Capturability analysis of a geometric guidance law in relative velocity space

In this paper, a relative velocity approach is used to analyze the capturability of a geometric guidance law. Point mass models are assumed for both the missile and the target. The speeds of the missile and target are assumed to remain constant throughout the engagement. Lateral acceleration, obtained from the guidance law, is applied to change the path of the missile. The kinematic equations for engagements in the horizontal plane are derived in the relative velocity space. Some analytical results for the capture region are obtained for non-maneuvering and maneuvering targets. For non-maneuvering targets it is enough for the navigation gain to be a constant to intercept the target, while for maneuvering targets a time varying navigation gain is needed for interception. These results are then verified through numerical simulations.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

A sufficient criterion for Kelvin–Helmholtz instability in the magnetopause boundary‐layer region

The Kelvin–Helmholtz instability has been investigated for the magnetopause boundary‐layer region by the linearized method. The plasma in magnetosheath and magnetopause is assumed to be semi‐infinitely extended homogeneous, nondissipative, and incompressible. It is observed that, if one relation of two plasma speeds on the two sides of the magnetopause, wave number, and boundary‐layer thickness exceeds a certain threshold, the instability sets in. This new analytically sufficient criterion for excitation of instability in the three‐layer plasma flow generalizes the corresponding Chandrasekhar’s instability criterion for two‐layer plasma flow. The known results have been recovered and modified, the new results have been discovered. It is proved that the velocity threshold for the onset of instability is low when the magnitude of the magnetosheath and boundary‐layer region magnetic field and the angle between them are small. Also the threshold depends on the direction of plasma flow. The following results are observed numerically. The growth of the instability is sensitive to the magnetic field direction in the magnetosheath. A slight variation in the magnetic field direction in the second region can substantially change the relative velocity threshold for instability. When the ratio of the density of the second and third layer (magnetosphere) increases or that of the first and third layer decreases, the threshold decreases. Apart from this a necessary criterion for instability is obtained for a particular case.

Capturability of realistic generalized true proportional navigation

The capturability of a realistic generalized true proportional navigation (RGTPN) guidance law, against a nonmaneuvering target, is analyzed. The RGTPN law is obtained by relaxing the somewhat unrealistic assumption of constant closing velocity, made in all earlier studies on generalized true proportional navigation (GTPN), and incorporating the actual time-varying value in the guidance law. Closed-form solutions for the complete capture region of RGTPN is obtained in terms of both zero and acceptable non-zero miss distances. It is shown that the capture region of RGTPN in the initial relative velocity space is significantly smaller than that of GTPN, for reasonable values of navigation constant (N) and angular direction (eta) of the missile commanded latax. However, for certain values of N and eta, capturability of RGTPN is found to be better. It is also shown that if in one of the versions of GTPN, which uses constant values of both the closing velocity and the line-of-sight (LOS) angular velocity in the guidance law, the corresponding realistic time-varying quantities are used, the capture region actually expands to cover the whole of the initial relative velocity space. A number of examples are given to compare the capture performance of RGTPN with other versions of the GTPN guidance laws.

Persistent passive hopping and juggling is possible even with plastic collisions

We describe simple one-dimensional models of passive (no energy input, no control), generally dissipative, vertical hopping and one-ball juggling. The central observation is that internal passive system motions can conspire to eliminate collisions in these systems. For hopping, two point masses are connected by a spring and the lower mass has inelastic collisions with the ground. For juggling, a lower point-mass hand is connected by a spring to the ground and an upper point-mass ball is caught with an inelastic collision and then re-thrown into gravitational free flight. The two systems have identical dynamics. Despite inelastic collisions between non-zero masses, these systems have special symmetric energy-conserving periodic motions where the collision is at zero relative velocity. Additionally, these special periodic motions have a non-zero sized, one-sided region of attraction on the higher-energy side. For either very large or very small mass ratios, the one-sided region of attraction is large. These results persist for mildly non-linear springs and non-constant gravity. Although non-collisional damping destroys the periodic motions, small energy injection makes the periodic motions stable, with a two-sided region of attraction. The existence of such special energy conserving solutions for hopping and juggling points to possibly useful strategies for both animals and robots. The lossless motions are demonstrated with a table-top experiment.

Collision Cones for Quadric Surfaces

The problem of collision prediction in dynamic environments appears in several diverse fields, which include robotics, air vehicles, underwater vehicles, and computer animation. In this paper, collision prediction of objects that move in 3-D environments is considered. Most work on collision prediction assumes objects to be modeled as spheres. However, there are many instances of object shapes where an ellipsoidal or a hyperboloid-like bounding box would be more appropriate. In this paper, a collision cone approach is used to determine collision between objects whose shapes can be modeled by general quadric surfaces. Exact collision conditions for such quadric surfaces are obtained in the form of analytical expressions in the relative velocity space. For objects of arbitrary shapes, exact representations of planar sections of the 3-D collision cone are obtained.

Capturability of Augmented Pure Proportional Navigation Guidance Against Time-Varying Target Maneuvers

This paper proposes a variation of the pure proportional navigation guidance law, called augmented pure proportional navigation, to account for target maneuvers, in a realistic nonlinear engagement geometry, and presents its capturability analysis. These results are in contrast to most work in the literature on augmented proportional navigation laws that consider a linearized geometry imposed upon the true proportional navigation guidance law. Because pure proportional navigation guidance law is closer to a realistic implementation of proportional navigation than true proportional navigation law, and any engagement process is predominantly nonlinear, the results obtained in this paper are more realistic than any available in the literature. Sufficient conditions on speed ratio, navigation gain, and augmentation parameter for capturability, and boundedness of lateral acceleration, against targets executing piecewise continuous maneuvers with time, are obtained. Further, based on a priori knowledge of the maximum maneuver capability of the target, a significant simplification of the guidance law is proposed in this paper. The proposed guidance law is also shown to require a shorter time of interception than standard pure proportional navigation and augmented proportional navigation. To remove chattering in the interceptor maneuver at the end phase of the engagement, a hybrid guidance law using augmented pure proportional navigation and pure proportional navigation is also proposed. Finally, the guaranteed capture zones of standard and augmented pure proportional navigation guidance laws against maneuvering targets are analyzed and compared in the normalized relative velocity space. It is shown that the guaranteed capture zone expands significantly when augmented pure proportional navigation is used instead of pure proportional navigation. Simulation results are given to support the theoretical findings.

A strongly coupled anisotropic fluid from dilaton driven holography

We consider a system consisting of 5 dimensional gravity with a negative cosmological constant coupled to a massless scalar, the dilaton. We construct a black brane solution which arises when the dilaton satisfies linearly varying boundary conditions in the asymptotically AdS(5) region. The geometry of this black brane breaks rotational symmetry while preserving translational invariance and corresponds to an anisotropic phase of the system. Close to extremality, where the anisotropy is big compared to the temperature, some components of the viscosity tensor become parametrically small compared to the entropy density. We study the quasi normal modes in considerable detail and find no instability close to extremality. We also obtain the equations for fluid mechanics for an anisotropic driven system in general, working upto first order in the derivative expansion for the stress tensor, and identify additional transport coefficients which appear in the constitutive relation. For the fluid of interest we find that the parametrically small viscosity can result in a very small force of friction, when the fluid is enclosed between appropriately oriented parallel plates moving with a relative velocity.

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.

Unsteady compressible flow due to relative rotation and translation of a viscous fluid on a disk

The modification of the axisymmetric viscous flow due to relative rotation of the disk or fluid by a translation of the boundary is studied. The fluid is taken to be compressible, and the relative rotation and translation velocity of the disk or fluid are time-dependent. The nonlinear partial differential equations governing the motion are solved numerically using an implicit finite difference scheme and Newton's linearisation technique. Numerical solutions are obtained at various non-dimensional times and disk temperatures. The non-symmetric part of the flow (secondary flow) describing the translation effect generates a velocity field at each plane parallel to the disk. The cartesian components of velocity due to secondary flow exhibit oscillations when the motion is due to rotation of the fluid on a translating disk. Increase in translation velocity produces an increment in the radial skin friction but reduces the tangential skin friction.

An Explicit Guidance Scheme For A Low-Thrust Terminal Rendezvous

An explicit near-optimal guidance scheme is developed for a terminal rendezvous of a spacecraft with a passive target in circular orbit around the earth. The thrust angle versus time profile for the continuous-thrust, constant-acceleration maneuver is derived, based on the assumption that the components of inertial acceleration due to relative position and velocity are negligible on account of the close proximity between the two spacecraft. The control law is obtained as a ''bilinear tangent law'' and an analytic solution to the state differential equations is obtained by expanding a portion of the integrand as an infinite series in time. A differential corrector method is proposed, to obtain real-time updates to the guidance parameters at regular time intervals. Simulation of the guidance scheme is carried out using the Clohessy-Wiltshire equations of relative motion as well as the inverse-square two-body equations of motion. Results for typical examples are presented.

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.