Evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source

The present work gives a comprehensive numerical study of the evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source. Using pseudospectral and predictor–corrector implicit finite difference methods, we first reproduced the known analytic results of the plane harmonic problem to a high degree of accuracy. The non-planar harmonic problems, for which the amplitude decay is faster than that for the planar case, are then treated. The results are correlated with the known asymptotic results of Scott (1981) and Enflo (1985). The constant in the old-age formula for the cylindrical canonical problem is found to be 1.85 which is rather close to 2, ‘estimated’ analytically by Enflo. The old-age solutions exhibiting strict symmetry about the maximum are recovered; these provide an excellent analytic check on the numerical solutions. The evolution of the waves for different source geometries is depicted graphically.

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.

Numerical simulation of solidification of liquid aluminum alloy flowing on cooling slope

Preparation of semisolid slurry using a cooling slope is increasingly becoming popular, primarily because of the simplicity in design and ease control of the process. In this process, liquid alloy is poured down an inclined surface which is cooled from underneath. The cooling enables partial solidification and the incline provides the necessary shear for producing semisolid slurry. However, the final microstructure of the ingot depends on several process parameters such as cooling rate, incline angle of the cooling slope, length of the slope and initial melt superheat. In this work, a CFD model using volume of fluid (VOF) method for simulating flow along the cooling slope was presented. Equations for conservation of mass, momentum, energy and species were solved to predict hydrodynamic and thermal behavior, in addition to predicting solid fraction distribution and macrosegregation. Solidification was modeled using an enthalpy approach and a volume averaged technique for the different phases. The mushy region was modeled as a multi-layered porous medium consisting of fixed columnar dendrites and mobile equiaxed/fragmented grains. The alloy chosen for the study was aluminum alloy A356, for which adequate experimental data were available in the literature. The effects of two key process parameters, namely the slope angle and the pouring temperature, on temperature distribution, velocity distribution and macrosegregation were also studied.

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.

Systematics of the magnetic-Prandtl-number dependence of homogeneous, isotropic magnetohydrodynamic turbulence

We present the results of our detailed pseudospectral direct numerical simulation (DNS) studies, with up to 1024(3) collocation points, of incompressible, magnetohydrodynamic (MHD) turbulence in three dimensions, without a mean magnetic field. Our study concentrates on the dependence of various statistical properties of both decaying and statistically steady MHD turbulence on the magnetic Prandtl number Pr-M over a large range, namely 0.01 <= Pr-M <= 10. We obtain data for a wide variety of statistical measures, such as probability distribution functions (PDFs) of the moduli of the vorticity and current density, the energy dissipation rates, and velocity and magnetic-field increments, energy and other spectra, velocity and magnetic-field structure functions, which we use to characterize intermittency, isosurfaces of quantities, such as the moduli of the vorticity and current density, and joint PDFs, such as those of fluid and magnetic dissipation rates. Our systematic study uncovers interesting results that have not been noted hitherto. In particular, we find a crossover from a larger intermittency in the magnetic field than in the velocity field, at large Pr-M, to a smaller intermittency in the magnetic field than in the velocity field, at low Pr-M. Furthermore, a comparison of our results for decaying MHD turbulence and its forced, statistically steady analogue suggests that we have strong universality in the sense that, for a fixed value of Pr-M, multiscaling exponent ratios agree, at least within our error bars, for both decaying and statistically steady homogeneous, isotropic MHD turbulence.

Acoustic performance of hoses - A parametric study

Modeling of wave propagation in hoses, unlike in rigid pipes or waveguides, introduces a coupling between the inside medium, the hose wall, and the outside medium, This alters the axial wave number and thence the corresponding effective speed of sound inside the hose resulting in sound radiation into the outside medium, also called the breakout or shell noise, The existing literature on the subject is such that a hose cannot be integrated into the,whole piping system made up of sections of hoses, pipes, and mufflers to predict the acoustical performance in terms of transmission loss (TL), The present paper seeks to fill this gap, Three one-dimensional coupled wave equations are written to account for the presence of a yielding wall with a finite lumped transverse impedance of the hose material, The resulting wave equation can readily be reduced to a transfer matrix form using an effective wave number for a moving medium in a hose section, Incorporating the effect of fluid loading due to the outside medium also allows prediction of the transverse TL and the breakout noise, Axial TL and transverse TL have been combined into net TL needed by designers, Predictions of the axial as well as transverse TL are shown to compare well with those of a rigorous 3-D analysis using only one-hundredth of the computation time, Finally, results of some parametric studies are reported for engineers involved in the acoustical design of hoses. (C) 1996 Institute of Noise Control Engineering.

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.

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.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

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.

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.