Turbulent gas-solid flows are readily encountered in many industrial and environmental processes. The development of a generic modeling technique for gas-solid turbulent flows remains a significant challenge in the field of mechanical engineering. Eulerian models are typically used to model large systems of particles. In this dissertation, a numerical analysis was carried out to assess a current state-of-the-art Eulerian two-fluid model for fully-developed turbulent gas-solid upward flow in a vertical pipe. The two-fluid formulation of Bolio et al. (1995) was adopted for the current study and the drag force was considered as the dominant interfacial force between the solids and fluid phase. In the first part of the thesis, a two-equation low Reynolds number k-ε model was used to predict the fluctuating velocities of the gas-phase which uses an eddy viscosity model. The stresses developed in the solids-phase were modeled using kinetic theory and the concept of granular temperature was used for the prediction of the solids velocity fluctuation. The fluctuating drag, i.e., turbulence modulation term in the transport equation of the turbulence kinetic energy and granular temperature was used to capture the effect of the presence of the dispersed solid particles on the gas-phase turbulence. The current study documents the performance of two popular turbulence modulation models of Crowe (2000) and Rao et al. (2011). Both models were capable of predicting the mean velocities of both the phases which were generally in good agreement with the experimental data. However, the phenomena that small particles cause turbulence suppression and large particles cause turbulence enhancement was better captured by the model of Rao et al. (2011); conversely, the model of Crowe (2000) produced turbulence enhancement in all cases. Rao et al. (2011) used a modified wake model originally proposed by Lun (2000) which is activated when the particle Reynolds number reaches 150. This enables the overall model to produce turbulence suppression and augmentation that follows the experimental trend. The granular temperature predictions of both models show good agreement with the limited experimental data of Jones (2001). The model of Rao et al. (2011) was also able to capture the effect of gas-phase turbulence on the solids velocity fluctuation for three-way coupled systems. However, the prediction of the solids volume fraction which depends on the value of the granular temperature shows noticeable deviations with the experimental data of Sheen et al. (1993) in the near-wall region. Both turbulence modulation models predict a flat profile for the solids volume fraction whereas the measurements of Sheen et al. (1993) show a significant decrease near the wall and even a particle-free region for flows with large particles. The two-fluid model typically uses a low Reynolds number k-ε model to capture the near-wall behavior of a turbulent gas-solid flow. An alternative near-wall turbulence model, i.e., the two-layer model of Durbin et al. (2001) was also implemented and its performance was assessed. The two-layer model is especially attractive because of its ability to include the effect of surface roughness. The current study compares the predictions of the two-layer model for both clear gas and gas-solid flows to the results of a conventional low Reynolds number model. The effects of surface roughness on the turbulence kinetic energy and granular temperature were also documented for gas-particle flows in both smooth and rough pipes.

The focus of the present licentiate thesis is on the effect of suction and pressure gradients on turbulent boundary-layer flows, which are investigated separately through performing numerical simulations.The first part aims at assessing history and development effects on adverse pressure-gradient (APG) turbulent boundary layers (TBL). A suitable set-up was developed to study near-equilibrium conditions for a boundary layer developingon a flat plate by setting the free-stream velocity at the top of the domain following a power law. The computational box size and the correct definition of the top-boundary condition were systematically tested. Well-resolved large-eddy simulations were performed to keep computational costs low. By varying the free-stream velocity distribution parameters, e.g. power-law exponent and virtual origin, pressure gradients of different strength and development were obtained. The magnitude of the pressure gradient is quantified in terms of the Clauser pressure-gradient parameter β. The effect of the APG is closely related to its streamwise development, hence, TBLs with non-constant and constant β were investigated. The effect was manifested in the mean flow through a much more pronounced wake region and in the Reynolds stresses through the existence of an outer peak. The terms of the turbulent kinetic energy budgets indicate the influence of the APG on the distribution of the transfer mechanism across the boundary layer. Stronger and more energetic structures were identified in boundary layers with relatively stronger pressure gradients in their development history. Due to the difficulty of determining the boundary-layer thickness in flows with strong pressure gradients or over a curvedsurface, a new method based on the diagnostic-plot concept was introduced to obtain a robust estimation of the edge of a turbulent boundary layer. In the second part, large-eddy simulations were performed on temporally developing turbulent asymptotic suction boundary layers (TASBLs). Findings from previous studies about the effect of suction could be confirmed, e.g. the reduction of the fluctuation levels and Reynolds shear stresses. Furthermore, the importance of the size of the computational domain and the time development were investigated. Both parameters were found to have a large impact on the results even on low-order statistics. While the mean velocity profile collapses in the inner layer irrespective of box size and development time, a wake region occurs for too small box sizes or early development time and vanishes once sufficiently large domains and/or integration times are chosen. The asymptotic state is charactersized by surprisingly thick boundary layers even for moderateReynolds numbers Re (based on free-stream velocity and laminar displacement thickness); for instance, Re = 333 gives rise to a friction Reynolds number Reτ = 2000. Similarly, the flow gives rise to very large structures in the outer region. These findings have important ramifications for experiments, since very large facilities are required to reach the asymptotic state even for low Reynolds numbers. / <p>QC 20160418</p>

The focus of this thesis is on the numerical study of subcritical transition to turbulence in the asymptotic suction boundary layer (ASBL). Applying constant homogeneous suction prevents the spatial growth of the boundary layer, granting access to the asymptotic dynamics. This enables research approaches which are not feasible in the spatially growing case. In a first part, the laminar–turbulent separatrix of the ASBL is investigated numerically by means of an edge-tracking algorithm. The consideration of spanwise-extended domains allows for the robust localisation of the attracting flow structures on this separatrix. The active part of the identified edge states consists of a pair of low- and high-speed streaks, which experience calm phases followed by high energy bursts. During these bursts the structure is destroyed and re-created with a shift in the spanwise direction. Depending on the streamwise extent of the domain, these shifts are either regular in direction and distance, and periodic in time, or irregular in space and erratic in time. In all cases, the same clear regeneration mechanism of streaks and vor- tices is identified, bearing strong similarities with the classical self-sustaining cycle in near-wall turbulence. Bifurcations from periodic to chaotic regimes are studied by varying the streamwise length of the (periodic) domain. The resulting bifurcation diagram contains a number of phenomena, e.g. multistability, intermittency and period doubling, usually investigated in the context of low-dimensional systems. The second part is concerned with spatio–temporal aspects of turbulent ASBL in large domains near the onset of sustained turbulence. Adiabatically decreasing the Reynolds number, starting from a fully turbulent state, we study low-Re turbulence and events leading to laminarisation. Furthermore, a robust quantitative estimate for the lowest Reynolds number at which turbulence is sustained is obtained at Re <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Capprox" /> 270. / <p>QC 20140213</p>