<div>Mixing of stratified fluids due to motion of bubble swarms can happen through two major mechanisms. The first is the capture and transport of heavier liquid into the lighter layers by the bubble wake. The second is the mixing due to turbulent dispersion. Stratification also affects bubble dynamics in various ways, namely by reducing the horizontal and vertical bubble fluctuations and extent, altering the drag experienced by rising bubbles, and changing the wake dynamics. The objective of this study is to understand these explained phenomena by decoupling their effects from each other and studying them individually. CFD offers powerful capabilities to achieve the decoupling and perform in-depth analysis of the fluid flow. </div><div><br></div><div>Firstly, the study of mixing induced in stratified fluids by bubbly flow in a Hele-Shaw Cell will be performed. Simulations are run for a range of void fractions and Froude numbers. The confinement prevents turbulence production, and mixing occurs primarily due to transport of colder liquid into the hotter layers by the bubble wake. Bubbles move in a zigzag motion attributed to the periodic vortex shedding in their wake. We report the formation of horizontal clusters and establish a direct correlation between the size of clusters and the rise velocity of the bubbles. We report an increase in the buoyancy flux across the isopycnals as the void fraction increases. The fraction of energy production due to the buoyancy flux increases with the strength of stratification, giving rise to a higher mixing efficiency. At the same time, cross isopycnal diffusion is higher at weaker stratification strengths.</div><div><br></div><div>Subsequently, direct numerical simulations of up to 146 bubbles rising in unbounded stratified fluids are performed. Both the bubble dynamics and destratification effects caused by the bubble motion are analyzed. The importance of bubble deformability and bubble Reynolds numbers on the induced background mixing are studied by varying the $E\ddot{o}tv\ddot{o}s$ number in the range 1.55 to 4.95 and Reynolds number in the range 25 to 200. Highly deformable, high Reynolds number bubbles undergo path instabilities and give rise to higher levels of mixing. Liquid and bubble velocity fluctuations and pseudo-turbulence caused by the bubble motion in the unconfined setting are examined and are seen to play an important role in mixing statistics. An increase in turbulent kinetic energy (TKE) levels with void fraction is noted. TKE levels are seen to decrease slightly as the stratification strength is increased, indicating increasing stability and resistance to destratification. Regardless of the stratification strength, a kinetic energy spectrum slope value between $-3 \sim -3.25$ is reported depending on Reynolds number. The dependence of mixing parameters on the void-fraction of bubbles and stratification strength of the liquid is also presented. </div><div><br></div><div>Next, the study of buoyancy driven motion of a single air bubble in stratified liquid is undertaken. A range of parameters including Froude number, Reynolds number and Bond number are explored. The Reynolds and Bond numbers will be maintained at values where the bubble motion and wake can be assumed to be axisymmetric. Wake dynamics and drift-volumes associated with the bubble rising in the stratified fluid are analyzed. The presence of secondary and tertiary vortices, which are alternating in direction, in the wake of the bubble due to the negative buoyant force experienced by the isopycnals is reported. The isopycnals oscillate before coming back to their stable state and the frequency of oscillations increases with stratification strength. The dependence of drag coefficient, determined by an unsteady force balance, and steady state bubble velocities, on the above mentioned parameters are studied. Analysis of bubble rise in partial stratification reveals the differences between homogeneous and stratified mediums.</div><div><br></div><div>Since most stratified bubbly flows occur near the free surface, an attempt is made at modeling the bubble rise up-to the free surface and subsequent bubble bursting. A brief study of in-line bubble coalescence is also attempted and potential future work for bubbly flows with topological changes is discussed.</div>

[Truncated abstract] Internal waves redistribute energy and momentum in stratified lakes and constitute the path through which the energy that is introduced at the lake scale is cascaded down to the turbulent scales where mixing and dissipation take place. This research, based on intensive field data complemented with numerical simulations, covers several aspects of the energy flux path ranging from basin-scale waves with periods of several hours to high frequency waves with periods of few minutes. It was found that, at the basin-scale level, the horizontal shape of the lake at the level of the metalimnion controls the period and modal structure of the basin-scale natural modes, conforming to the dispersion relationship of internal waves in circular basins. The sloping bottom, in turn, produces local intensification of the wave motion due to focusing of internal wave rays over near-critical slopes, providing hot spots for the degeneration of the basin-scale waves due to shear instabilities, nonlinear processes and dissipation. Different types of high-frequency phenomena were observed in a stratified lake under different forcing conditions. The identification of the generation mechanisms revealed how these waves extract energy from the mean flow and the basin-scale waves. The changes to the stratification show that such waves contribute to mixing in different ways . . . Detailed field observations were used to develop a comprehensive description of an undocumented energy flux mechanism in which shear-instabilities with significant amplitudes away from the generation level are produced in the surface layer due to the shear generated by the wind. The vertical structure of these instabilities is such that the growing wave-related fluctuations strain the density field in the metalimnion triggering secondary instabilities. These instabilities also transport energy vertically to the thermocline where they transfer energy back to the mean flow through interaction with the background shear.