A Novel Operator-Splitting Technique for One-Dimensional Laminar Premixed Flames

This paper is concerned the calculation of flame structure of one-dimensional laminar premixed flames using the technique of operator-splitting. The technique utilizes an explicit method of solution with one step Euler for chemistry and a novel probabilistic scheme for diffusion. The relationship between diffusion phenomenon and Gauss-Markoff process is exploited to obtain an unconditionally stable explicit difference scheme for diffusion. The method has been applied to (a) a model problem, (b) hydrazine decomposition, (c) a hydrogen-oxygen system with 28 reactions with constant Dρ 2 approximation, and (d) a hydrogen-oxygen system (28 reactions) with trace diffusion approximation. Certain interesting aspects of behaviour of the solution with non-unity Lewis number are brought out in the case of hydrazine flame. The results of computation in the most complex case are shown to compare very favourably with those of Warnatz, both in terms of accuracy of results as well as computational time, thus showing that explicit methods can be effective in flame computations. Also computations using the Gear-Hindmarsh for chemistry and the present approach for diffusion have been carried out and comparison of the two methods is presented.

The possible influence of curvature and rotation on ocean currents

Recent laboratory investigations have shown that rotation and (streamwise) curvature can have spectacular effects on momentum transport in turbulent shear flows. A simple model that takes account of these effects (based on an analogy with buoyant flows) utilises counterparts of the Richardson number Rg and the Monin-Oboukhov length. Estimates of Rg for meanders in ocean currents like the Gulf Stream show it to be of order 1 or more, while laboratory investigations reveal strong effects even at |Rg|∼0·1. These considerations lead to the conclusion that at a cyclonic bend in the Gulf Stream, a highly unstable flow in the outer half of the jet rides over a highly stable flow in the inner half. It is conjectured that the discrepancies noticed between observation and the various theories of Gulf Stream meanders, and such phenomena as the observed detachment of eddies from the Gulf Stream, may be due to the effects of curvature and rotation on turbulent transport.

Kinetics of Extraction of Iron(II1) by Bis(2-ethylhexyl Phosphate in a Laminar Liquid-Liquid Jet Reactor

The kinetics of iron(II1) extraction by bis(Zethylhexy1) phosphate (HDEHP, HA) in kerosene from sulfuric acid solutions has been studied in a liquid-liquid laminar jet reactor. The contact time of the interface in this reacting device is of the same order of magnitude as the surface renewal time in dispersion mixing and much less than that obtained in the relatively quiescent condition of the Lewis cell. Yet the analysis of the data in this study suggested a rate-controlling step involving surface saturation quite in conformity with that obtained in the Lewis cell and not with that in dispersion mixing as reported in the literature. Further, the mechanism suggested a weaker dependence of the rate on hydrogen ion concentration which was reported by other workers.

Development of pseudo-two-dimensional turbulent wakes

Pseudotwo-dimensional wakes are generated by introducing spanwise cellular structures into an otherwise plane turbulent wake by means of the castellated blunt trailing edges of different configurations. The transverse growths of these coflowing cellular wakes are found to be independent of each other without any noticeable spanwise interaction. This wake growth is examined in the light of the plane equilibrium wake analysis. Though these wakes are not found to be exactly self-similar, their growth shows a nonmonotonous approach toward the asymptotic state appropriate to that of a plane wake. The dye emission in the wakes illustrated a coherent vortical structure in the transverse plane, similar to that of the usual two-dimensional wake, in spite of the initial spanwise irregularities.

Breakage of viscoelastic drops in turbulent stirred dispersions

The existing models of drop breakage in stirred turbulent dispersions are applicable only to purely viscous dispersed phases. In their present form, they are found to underpredict the diameters of the largest stable drops formed when a viscoelastic fluid is dispersed into a Newtonian liquid. In purely viscous fluids, the turbulent stresses are opposed both by the stresses due to interfacial tension and the viscous stresses generated as the drop deforms. In viscoelastic fluids, drop deformation produces additional retractive elastic stresses which also oppose turbulent stresses. As the deformation rates are large, the retractive stresses can be large in magnitude. Assuming that these additional stresses decay with time, a model of viscoelastic drop breakage in turbulent stirred dispersions has been developed. The new model quantitatively predicts the dmax of viscoelastic fluids. The model, however, does not predict the observation that when the time constant of the fluid becomes large (λ > 0.5 s), the fluid can not be dispersed into droplets up to agitator speeds of about 10 rps in our equipment.

Wall Shear Fluctuations in a Turbulent Boundary Layer

IN the last two decades, the instantaneous structure of a turbulent boundary layer has been examined by many in an effort to understand the dynamics of the flow. Distinct and well-defined flow patterns that seem to have great relevance to the turbulence production mechanism have been observed in the wall region.1'2 The flow near the wall is intermittent with periodic eruptions of the fluid, a phenomenon generally termed "bursting process." Earlier investigations in this field were limited to liquid flows at low speeds and the entire flowpattern was observed using flow visualization techniques.Study was later extended to boundary-layer flows in windtunnels at higher speeds and Reynolds numbers using hot-wiresignals for the analysis of the bursting phenomenon.

Turbulent spot growth in favorable pressure gradients

Experimental results are presented on the lateral growth of turbulent spots in a series of flows with favorable pressure gradients. It is shown that the wedge angle increases slowly with the Reynolds number and that a favorable pressure gradient inhibits the growth of turbulent spots and, in general, results in a nonlinear turbulent wedge. As soon as the pressure gradient decreases to the point where the flow becomes supercritical, however, spot growth increases rapidly and the associated turbulent wedge becomes linear.

Relaminarization in highly accelerated turbulent boundary layers

The mean flow development in an initially turbulent boundary layer subjected to a large favourable pressure gradient beginning at a point x0 is examined through analyses expected a priori to be valid on either side of relaminarization. The ‘quasi-laminar’ flow in the later stages of reversion, where the Reynolds stresses have by definition no significant effect on the mean flow, is described by an asymptotic theory constructed for large values of a pressure-gradient parameter Λ, scaled on a characteristic Reynolds stress gradient. The limiting flow consists of an inner laminar boundary layer and a matching inviscid (but rotational) outer layer. There is consequently no entrainment to lowest order in Λ−1, and the boundary layer thins down to conserve outer vorticity. In fact, the predictions of the theory for the common measures of boundary-layer thickness are in excellent agreement with experimental results, almost all the way from x0. On the other hand the development of wall parameters like the skin friction suggests the presence of a short bubble-shaped reverse-transitional region on the wall, where neither turbulent nor quasi-laminar calculations are valid. The random velocity fluctuations inherited from the original turbulence decay with distance, in the inner layer, according to inverse-power laws characteristic of quasi-steady perturbations on a laminar flow. In the outer layer, there is evidence that the dominant physical mechanism is a rapid distortion of the turbulence, with viscous and inertia forces playing a secondary role. All the observations available suggest that final retransition to turbulence quickly follows the onset of instability in the inner layer.It is concluded that reversion in highly accelerated flows is essentially due to domination of pressure forces over the slowly responding Reynolds stresses in an originally turbulent flow, accompanied by the generation of a new laminar boundary layer stabilized by the favourable pressure gradient.

Turbulent Wall-Jets in Conical Diffusers

Abstract is not available.

Drag reduction by counterflow supersonic jet for a blunt cone in high enthalpy flows

Counterflow supersonic jet is used as a drag reduction device during the experiments in free piston driven shock tunnel, HST3. Accelerometer based force balance is employed to measure the drag force experienced by the 60-degree apex angle blunt cone model without and with the supersonic jet opposing the hypersonic flow. It is observed that the drag force decreases with increase in injection pressure ratio until the critical injection pressure is reached. Maximum reduction in drag force of 44 percent is recorded at the critical injection pressure ratio 22.36. Further increase in injection pressure ratio has reduced the percentage drag reduction. Change in nature of the flowfield around the model has also been observed across the critical injection pressure ratio.

Experimental and numerical investigation of jet injection in a wall bounded supersonic flow

Numerical and experimental studies of a supersonic jet (Helium) inclined at 45 degrees to a oncoming Mach 2 flow have been carried out. The numerical study has been used to arrive at a geometry that could reduce an oncoming Mach 5.75 flow to Mach 2 flow and in determining the jet parameters. Experiments are carried out in the IISc. hypersonic shock tunnel HST2 at similar conditions obtained from numerical studies. Flow visualization studies carried out using Schlieren technique clearly show the presence of the bow shock in front of the jet exposed to supersonic cross flow. The jet Mach number is experimentally found to be approximate to 3. Visual observations show that the jet has penetrated up to 60% of the total height of the chamber.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.

Laminar and turbulent natural convection in rectangular cavities

"This chapter discusses laminar and turbulent natural convection in rectangular cavities. Natural convection in rectangular two-dimensional cavities has become a standard problem in numerical heat transfer because of its relevance in understanding a number of problems in engineering. Current research identified a number of difficulties with regard to the numerical methods and the turbulence modeling for this class of flows. Obtaining numerical predictions at high Rayleigh numbers proved computationally expensive such that results beyond Ra ∼ 1014 are rarely reported. The chapter discusses a study in which it was found that turbulent computations in square cavities can't be extended beyond Ra ∼ O (1012) despite having developed a code that proved very efficient for the high Ra laminar regime. As the Rayleigh number increased, thin boundary layers began to form next to the vertical walls, and the central region became progressively more stagnant and highly stratified. Results obtained for the high Ra laminar regime were in good agreement with existing studies. Turbulence computations, although of a preliminary nature, indicated that a second moment closure model was capable of predicting the experimentally observed flow features."--Publisher Summary