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.

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Hidden order in serrated flow of metallic glasses

We report results of statistical and dynamic analysis of the serrated stress-time curves obtained from compressive constant strain-rate tests on two metallic glass samples with different ductility levels in an effort to extract hidden information in the seemingly irregular serrations. Two distinct types of dynamics are detected in these two alloy samples. The stress-strain curve corresponding to the less ductile Zr65Cu15Ni10Al10 alloy is shown to exhibit a finite correlation dimension and a positive Lyapunov exponent, suggesting that the underlying dynamics is chaotic. In contrast, for the more ductile Cu47.5Zr47.5Al5 alloy, the distributions of stress drop magnitudes and their time durations obey a power-law scaling reminiscent of a self-organized critical state. The exponents also satisfy the scaling relation compatible with self-organized criticality. Possible physical mechanisms contributing to the two distinct dynamic regimes are discussed by drawing on the analogy with the serrated yielding of crystalline samples. The analysis, together with some physical reasoning, suggests that plasticity in the less ductile sample can be attributed to stick-slip of a single shear band, while that of the more ductile sample could be attributed to the simultaneous nucleation of a large number of shear bands and their mutual interactions. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Modeling and control of unified power flow controller for transient stability

FACTS controllers are emerging as viable and economic solutions to the problems of large interconnected ne networks, which can endanger the system security. These devices are characterized by their fast response, absence of inertia, and minimum maintenance requirements. Thyristor controlled equipment like Thyristor Controlled Series Capacitor (TCSC), Static Var Compensator (SVC), Thyristor Controlled Phase angle Regulator (TCPR) etc. which involve passive elements result in devices of large sizes with substantial cost and significant labour for installation. An all solid-state device using GTOs leads to reduction in equipment size and has improved performance. The Unified Power Flow Controller (UPFC) is a versatile controller which can be used to control the active and reactive power in the Line independently. The concept of UPFC makes it possible to handle practically all power flow control and transmission line compensation problems, using solid-state controllers, which provide functional flexibility, generally not attainable by conventional thyristor controlled systems. In this paper, we present the development of a control scheme for the series injected voltage of the UPFC to damp the power oscillations and improve transient stability in a power system. (C) 1998 Elsevier Science Ltd. All rights reserved.

Tensile flow behaviour of 2.25Cr-1Mo ferritic steel base metal and simulated heat affected zone structures of 2.25Cr-1Mo weld joint

Tensile tests in the temperature range 298 to 873 K have been performed on 2.25Cr-1Mo base metal and simulated heat affected zone (HAZ) structures of its weld joint, namely coarse grain bainite, fine grain bainite and intercritical structure. Tensile flow behaviour of all the microstructural conditions could be adequately described by the Hollomon equation (sigma = K-1 epsilon(n1)) at higher (> 623 K) temperatures. Deviation from the Hollomon equation was observed at low strains and lower (< 623 K) temperatures. The Ludwigson modification of Hollomon's equation, sigma = K-1 epsilon(n1) + exp (K-2 + n(2) epsilon), was found to describe the flow curve. In general, the flow parameters n(1), K-1, n(2) and K-2 were found to decrease with increase in temperature except in the intermediate temperature range (423 to 623 K). Peaks/plateaus were observed in their variation with temperature in the intermediate temperature range coinciding with the occurrence of serrated flow in the load-elongation curve. The n(1) Value increased and the K-1 value decreased with the type of microstructure in the order: coarse grain bainite, fine grain bainite, base metal and intercritical structure. The variation of nl with microstructure has been rationalized on the basis of mean free path (MFP) of dislocations which is directly related to the inter-particle spacing. Larger MFP of dislocations lead to higher strain hardening exponents n(1).

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Unsteady flow and heat transfer on a semi-infinite flat plate with an aligned magnetic field

The unsteady laminar boundary layer flow of an electrically conducting fluid past a semi-infinite flat plate with an aligned magnetic field has been studied when at time t > 0 the plate is impulsively moved with a constant velocity which is in the same or opposite direction to that of free stream velocity. The effect of the induced magnetic field has been included in the analysis. The non-linear partial differential equations have been solved numerically using an implicit finite-difference method. The effect of the impulsive motion of the surface is found to be more pronounced on the skin friction but its effect on the x-component of the induced magnetic field and heat transfer is small. Velocity defect occurs near the surface when the plate is impulsively moved in the same direction as that of the free stream velocity. The surface shear stress, x-component of the induced magnetic field on the surface and the surface heat transfer decrease with an increasing magnetic field, but they increase with the reciprocal of the magnetic Prandtl number. However, the effect of the reciprocal of the magnetic Prandtl number is more pronounced on the x-component of the induced magnetic field. (C) 1999 Elsevier Science Ltd. All rights reserved.

Aerodynamics of delta wings with flaps at hypersonic speeds

Control surface effectiveness is an important parameter for any aeroplane. For a hypersonic aircraft, though the power required to operate the flaps is determined by low speed flying conditions, it is imperative to know the effect of flaps at hypersonic speeds. Hence, studies have been done on this topic by aerodynamicists for over 40 years. In spite of this, only a limited data is available in the literature on this subject. This paper discusses the experimental study of the effect of sweep on the aerodynamic characteristics of thin slab delta wings with flaps at hypersonic speeds. For the purpose of this investigation, a novel special thin six-component balance, which has a thickness of 4mm and can be housed inside wings with 8mm thickness, has been designed. The wings had a sweep of 76degrees, 70degrees and 65degrees, t/c of 0.053 and flaps with 12% of wing area and 12% of wing chord. Testing were done at Mach 8.2, Re number of 2.13 x 10(6) (based on chord), from alpha = -12degrees to 12degrees and flap angle of 20degrees, 30degrees and 40degrees. Separation lengths, measured from Schlieren pictures, clearly show that there is 'no appreciable' effect of sweep on them. Also, using a simple local flow field calculation, the separation has been identified to be transitional in nature. These features of separation reflect in the force data. Because of the small separation length, the flaps (inspite of their small size) were very effective in generating additional C-N, C-M and C-l, which increased with increase in flap angle. In general, the C-N, C-M and X-CP were unaffected by sweep for symmetric flap deflection at positive incidences and asymmetric flap case, For symmetric flap case at negative incidences, only C-N was not influenced by the sweep but C-M decreased and X-CP moved upstream as the sweep is decreased, The wing with lower sweep produces higher CA and lower (L/D)(max) for both symmetric and asymmetric flaps. The rolling moment and adverse yaw increased with decrease in sweep for asymmetric flap deflection. Newtonian theory is shown to be incapable of predicting the effect of sweep on C-l, C-n and on the incremental values of C-N, C-M and C-A. In conclusion, it can be said that a small flap is generally adequate for hypersonic aeroplanes provided they operate at altitudes where transitional and turbulent separation can be expected to occur. This would make the flaps effective and thus enable ample control authority.

Experimental study and optimization of Plasma Actuators for Flow control in subsonic regime

Experimental study and optimization of Plasma Ac- tuators for Flow control in subsonic regime PRADEEP MOISE, JOSEPH MATHEW, KARTIK VENKATRAMAN, JOY THOMAS, Indian Institute of Science, FLOW CONTROL TEAM | The induced jet produced by a dielectric barrier discharge (DBD) setup is capable of preventing °ow separation on airfoils at high angles of attack. The ef-fect of various parameters on the velocity of this induced jet was studied experimentally. The glow discharge was created at atmospheric con-ditions by using a high voltage RF power supply. Flow visualization,photographic studies of the plasma, and hot-wire measurements on the induced jet were performed. The parametric investigation of the charac- teristics of the plasma show that the width of the plasma in the uniform glow discharge regime was an indication of the velocity induced. It was observed that the spanwise and streamwise overlap of the two electrodes,dielectric thickness, voltage and frequency of the applied voltage are the major parameters that govern the velocity and the extent of plasma.e®ect of the optimized con¯guration on the performance characteristics of an airfoil was studied experimentally.

Oscillatory flow and temperature fields in an open tube with temperature difference across the ends

The oscillating flow and temperature field in an open tube subjected to cryogenic temperature at the cold end and ambient temperature at the hot end is studied numerically. The flow is driven by a time-wise sinusoidally varying pressure at the cold end. The conjugate problem takes into account the interaction of oscillatory flow with the heat conduction in the tube wall. The full set of compressible flow equations with axisymmetry assumption are solved with a pressure correction algorithm. Parametric studies are conducted with frequencies of 5-15 Hz, with one end maintained at 100 K and other end at 300 K. The flow and temperature distributions and the cooldown characteristics are obtained. The frequency and pressure amplitude have negligible effect on the time averaged Nusselt number. Pressure amplitude is an important factor determining the enthalpy flow through the solid wall. The frequency of operation has considerable effect on penetration of temperature into the tube. The density variation has strong influence on property profiles during cooldown. The present study is expected to be of interest in applications such as pulse tube refrigerators and other cryocoolers, where oscillatory flows occur in open tubes. (C) 2011 Elsevier Ltd. All rights reserved.

Studies on the effects of test gas on the flow field around large angle blunt cone flying at hypersonic Mach number

The effect of the test gas on the flow field around a 120degrees apex angle blunt cone has been investigated in a shock tunnel at a nominal Mach number of 5.75. The shock standoff distance around the blunt cone was measured by an electrical discharge technique using both carbon dioxide and air as test gases. The forebody laminar convective heat transfer to the blunt cone was measured with platinum thin-film sensors in both air and carbon dioxide environments. An increase of 10 to 15% in the measured heat transfer values was observed with carbon dioxide as the test gas in comparison to air. The measured thickness of the shock layer along the stagnation streamline was 3.57 +/- 0.17 mm in air and 3.29 +/- 0.26 mm in carbon dioxide. The computed thickness of the shock layer for air and carbon dioxide were 3.98 mm and 3.02 mm, respectively. The observed increase in the measured heat transfer rates in carbon dioxide compared to air was due to the higher density ratio across the bow shock wave and the reduced shock layer thickness.

Influence of temperature and strain rate on the flow stress of an FeAl alloy

The deformation behavior of an FeAl alloy processed by hot extrusion of water atomized powder has been investigated. Compression tests are performed in the temperature range 1073–1423 K and in the strain rate range 0.001–100 s−1 up to a true plastic strain of 0.5. The flow stress has been found to be strongly dependent on temperature as well as strain rate. The stress exponent in the power law rate equation is estimated to be in the range 7.0–4.0, decreasing with temperature. The activation energy for plastic flow in the range 1073–1373 K varies from 430 kJ mol−1 at low stresses to 340 kJ mol−1 at high stresses. However, it is fairly independent of strain rate and strain. The activation area has similarly shown a stress dependence and lies in the range 160–45b2. At 1423 K and at strain rates lower than 0.1 s−1 a strain rate sensitivity of 0.3 is observed with an associated activation energy of 375 kJ mol−1. The plastic flow in the entire range of temperature and strain rate investigated appears to be controlled by a diffusion mechanism. The results have revealed that it is possible to process the alloy by superplastic forming in the range 1373–1423 K at strain rates lower than 0.1 s−1.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.