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>