A Comprehensive Theory of Erosive Burning in Solid Rocket Propeilants

The theory of erosive burning has been constructed front first principles using turbulent boundary layer concepts. It is shown that the problem constitutes one of solution of flame propagation equation for turbulent flow. The final approximate solution for the case of single step overall kinetics reveals the combined effects of fluid mechanics and chemical kinetics. The results obtained from this theory are compared with earlier experimental results. The dependence of erosive burning characteristics on various parameters has been elucidated.

Uniformly valid analytical solution to the problem of a decaying shock wave

An explicit representation of an analytical solution to the problem of decay of a plane shock wave of arbitrary strength is proposed. The solution satisfies the basic equations exactly. The approximation lies in the (approximate) satisfaction of two of the Rankine-Hugoniot conditions. The error incurred is shown to be very small even for strong shocks. This solution analyses the interaction of a shock of arbitrary strength with a centred simple wave overtaking it, and describes a complete history of decay with a remarkable accuracy even for strong shocks. For a weak shock, the limiting law of motion obtained from the solution is shown to be in complete agreement with the Friedrichs theory. The propagation law of the non-uniform shock wave is determined, and the equations for shock and particle paths in the (x, t)-plane are obtained. The analytic solution presented here is uniformly valid for the entire flow field behind the decaying shock wave.

Cooling of a composite slab in a two-fluid medium

A mixed boundary value problem associated with the diffusion equation that involves the physical problem of cooling of an infinite parallel-sided composite slab in a two-fluid medium, is solved completely by using the Wiener-Hopf technique. An analytical solution is derived for the temperature distribution at the quench fronts being created by two different layers of cold fluids having different cooling abilities moving on the upper surface of the slab at constant speedv. Simple expressions are derived for the values of the sputtering temperatures of the slab at the points of contact with the respective layers, assuming the front layer of the fluid to be of finite width and the back layer of infinite extent. The main problem is solved through a three-part Wiener-Hopf problem of a special type and the numerical results under certain special circumstances are obtained and presented in the form of a table.

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Some observations of tip-vortex cavitation

Cavitation has been observed in the trailing vortex system of an elliptic planform hydrofoil. A complex dependence on Reynolds number and gas content is noted at inception. Some of the observations can be related to tension effects associated with the lack of sufficiently large-sized nuclei. Inception measurements are compared with estimates of pressure in the vortex obtained from LDV measurements of velocity within the vortex. It is concluded that a complete correlation is not possible without knowledge of the fluctuating levels of pressure in tip-vortex flows. When cavitation is fully developed, the observed tip-vortex trajectory shows a surprising lack of dependence on any of the physical parameters varied, such as angle of attack, Reynolds number, cavitation number, and dissolved gas content.

Study Of Fluid-Flow And Residence-Time Distribution In A Continuous Slab Casting Mold

The fluid-flow pattern and residence-time distribution (r.t.d.) of the fluid in a continuous casting mould have been studied using a water model. The two recirculating zones below the discharge ports have been found to be asymmetric. The effect of casting speed, discharge port diameter, shroud well depth and the immersion depth on r.t.d. have been investigated. The r.t.d. curve has been well represented by a model of two backmix cells of equal volume in series. The exist of the fluid has been found to be non-uniform across the cross-section of the mould. The fluid-flow pattern has been observed to change with time in a random fashion. Dead volume of upto 31.8% has been found with smaller discharge ports.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

Cooling of an infinite slab in a two-fluid medium

A mixed boundary-valued problem associated with the diffusion equation, that involves the physical problem of cooling of an infinite slab in a two-fluid medium, is solved completely by using the Wiener-Hopf technique. An analytical solution is derived for the temperature distribution at the quench fronts being created by two different layers of cold fluids having different cooling abilities moving on the upper surface of the slab at constant speed. Simple expressions are derived for the values of the sputtering temperatures of the slab at the points of contact with the respective layers, assuming one layer of the fluid to be of finite extent and the other of infinite extent. The main problem is solved through a three-part Wiener - Hopf problem of a special type, and the numerical results under certain special circumstances are obtained and presented in the form of a table.

Direct numerical simulations of statistically steady, homogeneous, isotropic fluid turbulence with polymer additives

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, and isotropic fluid turbulence. We find clear manifestations of dissipation-reduction phenomena: on the addition of polymers to the turbulent fluid, we obtain a reduction in the energy dissipation rate; a significant modification of the fluid-energy spectrum, especially in the deep-dissipation range; and signatures of the suppression of small-scale structures, including a decrease in small-scale vorticity filaments. We also compare our results with recent experiments and earlier DNS studies of decaying fluid turbulence with polymer additives.

Experimental studies of a distorted turbulent spot in a three-dimensional flow

We report here on the results of a series of experiments carried out on a turbulent spot in a distorted duct to study the effects of a divergence with straight streamlines preceded by a short stretch of transverse streamline curvature, both in the absence of any pressure gradient. It is found that the distortion produces substantial asymmetry in the spot: the angles at which the spot cuts across the local streamlines are altered dramatically (in contradiction of a hypothesis commonly made in transition zone modelling), and the Tollmien-Schlichting waves that accompany the wing tips of the spot are much stronger on the outside of the bend than on the inside. However there is no strong effect on the internal structure of the spot and the eddies therein, or on such propagation characteristics as overall spread rate and the celerities of the leading and trailing edges. Both lateral streamline curvature and non-homogeneity of the laminar boundary layer into which the spot propagates are shown to be strong factors responsible for the observed asymmetry. It is concluded that these factors produce chiefly a geometric distortion of the coherent structure in the spot, but do not otherwise affect its dynamics in any significant way.

The sputtering temperature of a cooling cylindrical rod without and with an insulated core in a two-fluid medium

Utilising Jones' method associated with the Wiener-Hopf technique, explicit solutions are obtained for the temperature distributions on the surface of a cylindrical rod without an insulated core as well as that inside a cylindrical rod with an insulated inner core when the rod, in either of the two cases, is allowed to enter, with a uniform speed, into two different layers of fluid with different cooling abilities. Simple expressions are derived for the values of the sputtering temperatures of the rod at the points of entry into the respective layers, assuming the upper layer of the fluid to be of finite depth and the lower of infinite extent. Both the problems are solved through a three-part Wiener-Hopf problem of special type and the numerical results under certain special circumstances are obtained and presented in tabular forms.

Peristaltic transport of two immiscible viscous fluids in a circular tube

Peristaltic motion of two immiscible viscous incompressible fluids in a circular tube is studied in pumping and copumping ranges under long-wavelength and low-Reynolds-number assumptions. The effect of the peripheral-layer viscosity on the time-averaged flux and the mechanical efficiency is studied. The formation and growth of the trapping zone in the core and the peripheral layer are explained. It is observed that the bolus volume in the peripheral layer increases with an increase in the viscosity ratio. The limits of the time-averaged flux (Q) over bar for trapping in the core are obtained. The trapping observed in the peripheral layer decreases in size with an increase in (Q) over bar but never disappears. The development of the complete trapping of the core fluid by the peripheral-layer fluid with an increase in the time-averaged flux is demonstrated. The effect of peripheral-layer viscosity on the reflux layer is investigated. It is also observed that the reflux occurs in the entire pumping range for all viscosity ratios and it is absent in the entire range of copumping.

Stability of spatially developing boundary layers in pressure gradients

A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R(-1) in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R(-1)) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R(-1)) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R(-1) except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism.

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.