A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.

Boundary conditions at a rigid wall for rough granular gases

We derive boundary conditions at a rigid wall for a granular material comprising rough, inelastic particles. Our analysis is confined to the rapid flow, or granular gas, regime in which grains interact by impulsive collisions. We use the Chapman-Enskog expansion in the kinetic theory of dense gases, extended for inelastic and rough particles, to determine the relevant fluxes to the wall. As in previous studies, we assume that the particles are spheres, and that the wall is corrugated by hemispheres rigidly attached to it. Collisions between the particles and the wall hemispheres are characterized by coefficients of restitution and roughness. We derive boundary conditions for the two limiting cases of nearly smooth and nearly perfectly rough spheres, as a hydrodynamic description of granular gases comprising rough spheres is appropriate only in these limits. The results are illustrated by applying the equations of motion and boundary conditions to the problem of plane Couette flow.

Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 2. The fluctuating-force model

A fluctuating-force model is developed for representing the effect of the turbulent fluid velocity fluctuations on the particle phase in a turbulent gas–solid suspension in the limit of high Stokes number, where the particle relaxation time is large compared with the correlation time for the fluid velocity fluctuations. In the model, a fluctuating force is incorporated in the equation of motion for the particles, and the force distribution is assumed to be an anisotropic Gaussian white noise. It is shown that this is equivalent to incorporating a diffusion term in the Boltzmann equation for the particle velocity distribution functions. The variance of the force distribution, or equivalently the diffusion coefficient in the Boltzmann equation, is related to the time correlation functions for the fluid velocity fluctuations. The fluctuating-force model is applied to the specific case of a Couette flow of a turbulent particle–gas suspension, for which both the fluid and particle velocity distributions were evaluated using direct numerical simulations by Goswami & Kumaran (2010). It is found that the fluctuating-force simulation is able to quantitatively predict the concentration, mean velocity profiles and the mean square velocities, both at relatively low volume fractions, where the viscous relaxation time is small compared with the time between collisions, and at higher volume fractions, where the time between collisions is small compared with the viscous relaxation time. The simulations are also able to predict the velocity distributions in the centre of the Couette, even in cases in which the velocity distribution is very different from a Gaussian distribution.

A sufficient criterion for Rayleigh-Taylor instability of incompressible viscous three-layer flow

Using normal mode analysis Rayleigh-Taylor instability is investigated for three-layer viscous stratified incompressible steady flow, when the top 3rd and bottom 1st layers extend up to infinity, the middle layer has a small thickness δ. The wave Reynolds number in the middle layer is assumed to be sufficiently small. A dispersion relation (a seventh degree polynomial in wave frequency ω) valid up to the order of the maximal value of all possible Kj (j less-than-or-equals, slant 0, K is the wave number) in each coefficient of the polynomial is obtained. A sufficient condition for instability is found out for the first time, pursuing a medium wavelength analysis. It depends on ratios (α and β) of the coefficients of viscosity, the thickness of the middle layer δ, surface tension ratio T and wave number K. This is a new analytical criterion for Rayleigh-Taylor instability of three-layer fluids. It recovers the results of the corresponding problem for two-layer fluids. Among the results obtained, it is observed that taking the coefficients of viscosity of 2nd and 3rd layers same can inhibit the effect of surface tension completely. For large wave number K, the thickness of the middle layer should be correspondingly small to keep the domain of dependence of the threshold wave number Kc constant for fixed α, β and T.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.

Unsteady flow and heat transfer between rotating coaxial disks

A study is made on the flow and heat transfer of a viscous fluid confined between two parallel disks. The disks are allowed to rotate with different time dependent angular velocities, and the upper disk is made to approach the lower one with a constant speed. Numerical solutions of the governing parabolic partial differential equations are obtained through a fourth-order accurate compact finite difference scheme. The normal forces and torques that the fluid exerts on the rotating surfaces are obtained at different nondimensional times for different values of the rate of squeezing and disk angular velocities. The temperature distribution and heat transfer are also investigated in the present analysis.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Pressure sensitive flow and constraint factor in amorphous materials below glass transition

The constraint factor, C (given by the hardness-yield strength ratio H/Y in the fully lastic regime of indentation), in metallic glasses, is greater than three, a reflection of the sensitivity of their plastic flow to pressure. Furthermore, C increases with increasing temperature. In this work, we examine if this is true in amorphous polymers as well, through experiments on amorphous poly(methyl methacrylate) (PMMA). Uniaxial compression as well as spherical indentation tests were conducted in the 248-348 K range to construct H/Y versus indentation strain plots at each temperature and obtain the C-values. Results show that C increases with temperature in PMMA as well. Good correlation between the loss factors, measured using a dynamic mechanical analyzer, and C, suggest that the enhanced sensitivity to pressure is possibly due to beta-relaxation. We offer possible mechanistic reasons for the observed trends in amorphous materials in terms of relaxation processes.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

Plastic response of orthotropic spherical shells under blast loading

The plastic response of a segment of a simply supported orthotropic spherical shell under a uniform blast loading applied on the convex surface of the shell is presented. The blast is assumed to impart a uniform velocity to the shell surface initially. The material of the shell is orthotropic obeying a modified Tresca yield hypersurface conditions and the associated flow rules. The deformation of the shell is determined during all phases of its motion by considering the motion of plastic hinges in different regimes of flow. Numerical results presented include the permanent deformed configuration of the shell and the total time of shell response for different degrees of orthotropy. Conclusions regarding the plastic behaviour of spherical shells with circumferential and meridional stiffening under uniform blast load are presented.

Flow-graph techniques for epicyclic gear train analysis

Flow-graph techniques are applied in this article for the analysis of an epicyclic gear train. A gear system based on this is designed and constructed for use in Numerical Control Systems.

Is there an athermal region of plastic flow in metals at intermediate temperatures?

The question of the existence or otherwise of an athermal temperature region of plastic flow in metals is examined. It is suggested that the athermal region is absent in metals with large dislocation densities. Such an explanation is provoked by a fairly recent proposition that the unzipping of attractive junctions is a plausible rate-controlling mechanism at high temperatures.

Propagation of a spherical shock in an inhomogeneous self-gravitating or nongravitating system

The propagation of a shock wave of finite strength due to an explosion into inhomogeneous nongravitating and self-gravitating systems has been considered, using similarity principles, supposing that the density varies as an inverse power of distance from the centre of explosion. A large number of systems, characterised by different density exponents and different adiabatic coefficients of the gas have been considered for different shock strengths. The numerical integration from the shock inward has been continued to the surface of singularity where density tends to infinity and which acts like a piston in the self-gravitating case and to the surface where the velocity gradient tends to infinity in the nongravitating case. The effect of variation of shock strength, density exponent and adiabatic coefficient on the location of these singularities and on the distribution of flow parameters behind the shock has been studied. The initial energy of the system and the manner of release of the explosion energy influence strongly the flow behind the shock. The results have been graphically depicted.