The normal injection of a liquid jet into a gaseous crossflow has many engineering applications. In this thesis, detailed numerical simulations based on the level set method are employed to understand the physical mechanism underlying the jet ``surface breakup''. The numerical observations reveal the existence of hydrodynamic instabilities on the jet periphery. The temporal growth of such azimuthal instabilities leads to the formation of interface corrugations, which are eventually sheared off of the jet surface as sheet-like structures. The sheets finally undergo disintegration into ligaments and drops during the surface breakup process. Temporal linear stability analyses are employed to understand the nature of these instabilities. To facilitate the analysis, analytical solutions for the flow fields of the jet and the crossflow are derived. We identify the ``shear instability'' as the primary destabilization mechanism in the flow. This inherently inviscid mechanism opposes the previously suggested mechanism of surface breakup (known as ``boundary layer stripping''), which is based on a viscous interpretation. The influence of the jet-to-crossflow density ratio on the flow stability are also studied. The findings show that a higher density jet leads to higher wavenumber instabilities on the jet surface and thereby subsequent smaller drops and ligaments. The stability characteristics of the most amplified modes (i.e., the wavenumber and corresponding growth rate) obtained from stability analyses and numerical simulations are in good agreement. The stability results of the jet also show that the density may have a non-monotonic stabilizing/destabilizing effect on the flow stability. To investigate such effect, the concept of wave resonance are employed to physically interpret the inviscid instability mechanism in two-phase flows with sharp interfaces and linear velocity profiles. We demonstrate that neutrally stable waves are formed due to the density jump in the flow, in addition to the well-known vorticity (Rayleigh) waves. Under certain conditions, such neutral waves are capable of resonating and generating unstable modes. The resonance of different pairs of neutral waves, therefore, results in either stabilizing or destabilizing effect of density variation. We predict similar reasoning behind the density behavior in the jet in crossflow configuration with smoothly varying velocity and density profiles.

Large-eddy simulations of turbulent flows past a circular cylinder have been performed at sub-, critical and super-critical Re using an orthogonal curvilinear grid solver, CFDship-Iowa version 6.2. An extensive verification and validation study has been carried out. Various aspects of the flow field have been investigated. The aspect ratio of the computational domain has major effects on the results. In general, large aspect ratio produced best results for the sub-critical Re. Small dependency on both aspect ratio and grid resolution was observed for the critical Re. Small aspect ratio and conservative scheme produced best results for the super-critical Re. Overall flow features and the drag crisis phenomenon have been correctly predicted. A lot of experimental and numerical studies of flow past a circular cylinder were collected and used for the validation of the present LES study. Integral and local variables were in fairly good agreement for the sub-critical Re. Sharp behavior including drag crisis was predicted for the critical Re. Although some discrepancy including early formation of turbulent separation was observed, local flow structures including separation bubble were observed for the super-critical Re. The formation of secondary vortex near the cylinder wall and its evolution into separation bubble were observed. The spectral analysis showed that the separation bubble had the instabilities close to the shear layer frequency. The proximity of shear layer to the cylinder enhanced the mixing process of boundary layer and shear layer and led to the formation of separation bubble. A snapshot POD method was used to extract flow structures in the boundary layer, shear layer and wake. In the boundary layer, the secondary vortices and separation bubble were successfully extracted. Due to the weak TKE distribution, specific flow structures were hard to find in the shear layer. Large two-dimensional flow structures representing the Karman shedding vortices were extracted for the sub- and super-critical Re.

Separating and reattaching flows are important in a large number of engineering configurations. The flow through a sudden expansion (backward-facing step) represents a conceptually simple case of this class of flows and hence has been the subject of numerous studies. The present study focuses on the effect of the expansion ratio (defined as the ratio of downstream channel height to upstream channel height) on the unsteady flow features in the reattachment region and further downstream. It is known that this flow demonstrates two different instabilities; the Kelvin-Helmholtz shear layer instability, which scales with the shear layer thickness, and the instability associated with the separation bubble, which scales with the step height and has similarities to K´arman vortex shedding behind a cylinder.In addition to these, there is also a possibility of the presence of the ‘preferred’ mode of the jet issuing from the inlet channel of the sudden expansion, especially at high expansion ratios, where the flow resembles a wall jet. The aim of this study is to investigate experimentally the changes in the instability of the separation bubble, as the expansion ratio is changed, and its possible interactions with the other instabilities in the flow.One might expect some changes in the flow with expansion ratio, as at low expansion ratios, the configuration represents a simple backward-facing step geometry, while at high expansion ratios, the geometry approaches that of a wall jet. A variable expansion ratio backward-facing step facility has been developed in an open circuit wind tunnel.This facility permits continuous variation of the expansion ratio from 1 to around 6. Attention is focused on the turbulent regime of the flow, where the flow structure has been found in previous studies to be relatively insensitive to the Reynolds number. The inlet conditions have been kept constant with a thin turbulent boundary layer at the step, the boundary layer thickness at separation being approximately 14 % of the inlet channel height. The Reynolds number based on the inlet channel height, H, is kept constant at Re=48,000 and the expansion ratio is varied by changing the channel height downstream of the step. Detailed hot wire measurements have been made to characterize the spatial variation of the dominant frequencies in the flow at different expansion ratios. The expansion ratio has been varied from a low value of 1.14 to a high value of 3.25, and detailed measurements are obtained for five expansion ratios of 1.14, 1.3, 1.5, 2.0 and 3.0. Further, to elucidate the dominant vortical structures in the flow, Particle Image Velocimetry measurements have been undertaken simultaneously with hot wire measurements for the case of expansion ratio 1.5, which have permitted the conditional averaging of vorticity fields.These investigations have brought forth some interesting features of the flow over a backward-facing step. Results for the time-mean properties of the flow indicate that the shear layer separating from the step deviates from a free mixing layer behaviour away from the step, possibly due to its interaction with the wall and the recirculation region underneath it. At any given streamwise location, the shear layer momentum thickness, θ, is seen to increase with the expansion ratio. Further, upto reattachment, the momentum thickness of the shear layer is seen to scale with the step height, h, independent of its initial thickness at separation, θo, as long as the boundary layer at separation is sufficiently thin as compared to the step height. Investigations for the unsteady flow features show that the frequency of the dominant peak in the velocity spectrum, supposed to represent the passage frequency (Strouhal number, S, based on the step height, h, and the inlet velocity, U) of the vortical structures, varies in the cross stream (y) direction, in addition to its expected variation in the streamwise (x) direction. The variation of the Strouhal number in the cross stream direction is seen to scale with the local momentum thickness of the shear layer, except for locations very close to the step. To characterize the development of the dominant frequency in the streamwise direction, the maximum value of the Strouhal number at a streamwise location is taken to be the representative value for that streamwise location. The Strouhal number is seen to decrease in the streamwise direction, from a very high value near the step, to a value of approximately 0.08 in the reattachment region, and remains constant further downstream. This value, supposed to represent the large scale structures shed from the reattachment region, is seen to remain very close to 0.08 for all Expansion ratios investigated. Conditional averaging of the vorticity fields in the reattachment region is done for an expansion ratio of 1.5, to get a detailed picture of the unsteady flow field. The hot wire signal at the outer edge of the shear layer in the reattachment region, which represents the non-dimensional structure passage frequency of S=0.08, is used as the conditioning signal. Results seem to indicate that the recirculation region, or the ‘bubble’ divides into two cells, and sheds the downstream cell quasi-periodically. The passage of these structures through the reattachment region seems to be concomitant With a local vertical motion of the shear layer. Further, the streamwise development of the local Strouhal number, Sθ, based on the local momentum thickness of the shear layer, and the local free stream velocity, Umax, indicates a possibility of a coupling between the shear layer and the structures shed from the reattachment region.