Este Trabalho apresenta novas soluções exatas para as equações de Navier – Stokes transientes tridimensionais para escoamentos viscosos incompressíveis. Estas soluções são obtidas por meio de Split e Transformações Auto-Bäcklund. O procedimento de Split desacopla as equações de Navier – Stokes em dois sistemas de equações diferenciais parciais, um linear e outro não-linear, ambos não-homogêneos. O sistema linear, que contém somente termos viscosos e derivadas temporais, é resolvido via Transformações Auto-Bäcklund induzidas por relações de comutação, fornecendo o campo de velocidades. Os componentes do vetor velocidade são então substituídos no sistema não-linear a fim de obter o correspondente campo de pressões. A resolução do sistema não-linear para a pressão pode ser obtida tanto numericamente (via integração direta) quanto analiticamente, empregando a equação de Helmholtz. O objetivo do presente trabalho é encontrar expressões analíticas para o campo de velocidades e obter resultados numéricos para o campo de pressão associado. O caráter híbrido das soluções proporciona uma redução significativa do tempo de processamento requerido para a simulação de escoamentos viscosos, o qual praticamente se reduz ao tempo demandado para a tarefa de pós-processamento. Com esse objetivo em mente, foi desenvolvida uma formulação tridimensional escalar para a função corrente, a fim de reduzir o tempo requerido na tarefa mais dispendiosa de pós-processamento, a saber, o traçado das linhas de corrente em torno de corpos submersos de formato arbitrário. Neste estágio de desenvolvimento, esta formulação é empregada para produzir mapas de linhas de corrente para escoamentos viscosos em torno de uma esfera para números de Reynolds elevados. / This work presents new exact solutions to the unsteady three dimensional Navier-Stokes equations for incompressible viscous flows. These solutions are obtained by means of split and auto-Bäcklund transformations. The splitting procedure decouples the Navier-Stokes equations into a linear and a nonlinear inhomogeneous system of partial differential equations. The linear system, which contains only viscous terms and time derivatives, is solved via auto-Bäcklund transformations induced by commutation relations, furnishing the velocity field. The components of the velocity vector are then replaced into the nonlinear system to obtain the corresponding pressure field. The solution of the nonlinear system for the pressure variable can be carried out either numerically (by direct integration) or analytically, using the Helmholtz equation . The aim of the proposed work is to find analytical expressions for the velocity field and to obtain numerical results to the associated pressure field. The hybrid character of the solutions provides a significant reduction on the time processing required to simulate viscous flows, which virtually reduces to the time demanded to execute post-processing tasks. Taking this fact in mind, a three dimensional scalar formulation for the streamfunction was developed in order to simplify the most time-consuming post-processing task required, e.g., plotting the streamlines around arbitrary shaped bodies. At this stage of development, this formulation is employed to produce streamline maps for viscous flows around a sphere for high Reynolds numbers.

Made available in DSpace on 2014-10-09T12:38:39Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:05:15Z (GMT). No. of bitstreams: 1 06042.pdf: 9369626 bytes, checksum: 37950a2f878d6535f671de4a025da71c (MD5) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

Made available in DSpace on 2015-05-08T15:00:05Z (GMT). No. of bitstreams: 1 parte1.pdf: 3576583 bytes, checksum: 2028463a4df5689828411372ce4dd847 (MD5) Previous issue date: 2006-02-24 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The present work consists of a hybrid numerical-analytical simulation of the developing and fully developed turbulent flows inside a parallel-plates channel. The Reynolds equations for the mean flow are solved through the Generalized Integral Transform Technique (GITT) in the boundary layer and streamfunction-only formulations. For turbulence closure, the respectively one-equation turbulence models, and some of their variations, developed by Secundov (1972), Baldwin e Barth (1990), Spalart e Allmaras (1992, 1994) and Menter (1997) were employed. Despite of based on the eddy viscosity concept, as the well-known one-equation k-L turbulence model previously adopted in works that made use of the present hybrid solution methodology, such models do not need any explicit length scale, and therefore are more generals, but only one transport equation for the turbulent viscosity or for a variable directly related to the eddy viscosity. In this sense, some simulations for different Reynolds numbers and different turbulent inlet conditions were developed in order to, in function of the obtained results and the convergence studies of the main potentials involved in the simulations, develop a critically and detailed discussion of the main shortcomings and capability predictions of each turbulence model adopted, such as the non-asymptotic development of the longitudinal velocity component and the friction factor fields. As a consequence of the analysis type performed over each turbulence model investigated and, in addition, due to the excellent numerical quality of the obtained results, the present work extends and consolidates the very important role that the integral transform technique may play in the computational fluid dynamics field meanwhile hybrid methodology. / O presente trabalho consiste na simulação numérico-analítica do escoamento turbulento isotérmico, em desenvolvimento e completamente desenvolvido, no interior de um canal de placas planas e paralelas, através da Técnica da Transformada Integral Generalizada (GITT). As equações médias de Reynolds, escritas segundo as hipótese de camada limite são empregadas na formulação de função corrente. Para o fechamento das equações médias da turbulência foram empregados os modelos de turbulência de uma equação de transporte (e algumas de suas variações) desenvolvidos pelos seguintes autores: Secundov (1972), Baldwin e Barth (1990), Spalart e Allmaras (1992, 1994) e Menter (1997). Apesar de baseados no conceito de viscosidade turbilhonar ( eddy viscosity ), como o tradicional modelo de turbulência de uma equação k-L, anteriormente empregado em trabalhos que fizeram uso da presente metodologia híbrida de solução, tais modelos não necessitam de uma escala explícita de comprimento, e por isso são ditos mais gerais, mas sim, de apenas uma equação de transporte para a própria viscosidade turbulenta (ou uma variável diretamente relacionada à própria viscosidade turbulenta). Nesse sentido, foram simuladas algumas situações de escoamentos (diferentes números de Reynolds e diferentes condições de entrada turbulenta) e buscou-se, em função dos resultados obtidos e de estudos de convergência dos principais potenciais envolvidos nas simulações, uma criteriosa e detalhada discussão das capacidades de previsão de certas características inerentes ao escoamento (o pico da componente longitudinal da velocidade no centro do canal e a depressão do fator de atrito, durante o desenvolvimento do escoamento, ou seja, o comportamento não assintótico dessas variáveis no escoamento) por todos os modelos de turbulência empregados. Em função do tipo de análise oferecida sobre cada modelo de turbulência utilizado, e, em adição, da excelente qualidade dos resultados obtidos com o uso da GITT, o presente trabalho extende e consolida o importante papel que a técnica da transformada integral pode desempenhar no campo da simulação computacional, enquanto metodologia híbrida.

No labirinto dos gotejadores ocorrem diferentes interações entre as partículas de argila e as características do escoamento, sendo que essas interações interferem no potencial de obstrução dos mesmos. Considerando que os principais fatores que interferem no processo de obstrução causado por partículas de argila estão relacionados à natureza das argilas, ao regime de escoamento, à força iônica e ao pH da solução, e à concentração de partículas em suspensão. Esta pesquisa apresentou como objetivo identificar como esses fatores podem influenciar o desempenho dos gotejadores. Além disso, visando a caracterização do processo de obstrução causado por partículas desta natureza, buscou-se: relacionar o potencial de agregação das partículas com o desempenho dos gotejadores; e, analisar o comportamento da deposição de partículas no interior dos labirintos dos gotejadores. Para isso, além da realização dos ensaios em bancada utilizando diferentes dispersões de argila, foram realizadas análises sobre o potencial de agregação das partículas, e, utilizando um dispositivo milifluidico, foi realizada uma análise de deposição de partículas. Durante a pesquisa foram utilizados dois tipos de argilas, caulinita e montmorilonita. A força iônica da solução e a natureza dos materiais de argila apresentaram efeito sobre o desempenho dos gotejadores, enquanto que a caulinita em solução salina de sódio proporcionou incremento máximo de vazão de 5%, a montmorilonita causou redução máxima de vazão de 15%. No geral, as partículas de argila começaram a causar redução de vazão significativas com concentrações iguais ou superiores a 1000 mg L-1. Observou-se que as partículas de argila estão sujeitas ao fenômeno de autolimpeza que ocorre no interior dos labirintos após o acionamento do sistema. Não foi possível encontrar uma relação entre o potencial de agregação das partículas de argila e o desempenho dos gotejadores. Possivelmente, o índice de agregação utilizado pela teoria DLVO clássica não seja o mais adequado para obter esta relação. Em relação a deposição de partículas, as regiões que apresentaram maiores deposições de partículas coincidem com aquelas que apresentam menores valores de intensidade de turbulência e energia cinética turbulenta. A deposição de partícula ocorre principalmente nas regiões de vórtices e estagnação, localizadas nos dois primeiros defletores dos labirintos. Não foi observado acúmulo de partículas na região do fluxo principal. Por isso, sugere-se que partículas de argila, como agente isolado de obstrução, não apresentam potencial para causar obstrução completa dos gotejadores. / In the drippers labyrinth different interactions occur between the clay particles and the flow characteristics, and these interactions interfere in the clogging potential of the particles. Considering the main factors that affect the process of clogging caused by clay particles are related to the type of the clays, the flow regime, the ionic strength and pH of the solution, and the concentration of suspended particles. This research aims to identify how these factors can influence the performance of drippers. In addition, aiming to characterize the clogging process caused by particles of clay, it was sought to relate the potential of particles aggregation with the performance of drippers; and, to analyze the behavior of particles deposition inside the labyrinths of drippers. For this, in addition to performing bench tests using different clay dispersions, analyzes were carried out on the particle aggregation potential, and a particle deposition analysis was performed using a microfluidic device. During the research two types of clays were used, kaolinite and montmorillonite. The ionic strength of the solution and the type of the clay materials had an effect on the drippers performance. Whereas the kaolinite in sodium saline provided a maximum flow rate increase of 5%, the montmorillonite caused a maximum flow rate reduction of 15%. In general, clay particles began to cause significant flow rate reduction with concentrations equal to or greater than 1000 mg L-1. It was observed that clay particles are subject to the phenomenon of self-cleaning which occurs inside the labyrinths after turn on the system. It was not possible to find a relation between the potential of clay particles aggregation and the drippers performance. Possibly, the aggregation index used by the classical DLVO theory is not the most adequate to obtain this relation. In relation to the particles deposition, the regions of greatest deposition coincide with those that present lower values of turbulence intensity and turbulent kinetic energy. Particles deposition occurs mainly in the vortex and stagnation regions, located in the first two baffles of the labyrinths. In addition, no particle accumulation occurs in the main stream region, it suggested that clay particles, as an isolated clogging agent, do not have the potential to cause total clogging of drippers.

Este Trabalho apresenta novas soluções exatas para as equações de Navier – Stokes transientes tridimensionais para escoamentos viscosos incompressíveis. Estas soluções são obtidas por meio de Split e Transformações Auto-Bäcklund. O procedimento de Split desacopla as equações de Navier – Stokes em dois sistemas de equações diferenciais parciais, um linear e outro não-linear, ambos não-homogêneos. O sistema linear, que contém somente termos viscosos e derivadas temporais, é resolvido via Transformações Auto-Bäcklund induzidas por relações de comutação, fornecendo o campo de velocidades. Os componentes do vetor velocidade são então substituídos no sistema não-linear a fim de obter o correspondente campo de pressões. A resolução do sistema não-linear para a pressão pode ser obtida tanto numericamente (via integração direta) quanto analiticamente, empregando a equação de Helmholtz. O objetivo do presente trabalho é encontrar expressões analíticas para o campo de velocidades e obter resultados numéricos para o campo de pressão associado. O caráter híbrido das soluções proporciona uma redução significativa do tempo de processamento requerido para a simulação de escoamentos viscosos, o qual praticamente se reduz ao tempo demandado para a tarefa de pós-processamento. Com esse objetivo em mente, foi desenvolvida uma formulação tridimensional escalar para a função corrente, a fim de reduzir o tempo requerido na tarefa mais dispendiosa de pós-processamento, a saber, o traçado das linhas de corrente em torno de corpos submersos de formato arbitrário. Neste estágio de desenvolvimento, esta formulação é empregada para produzir mapas de linhas de corrente para escoamentos viscosos em torno de uma esfera para números de Reynolds elevados. / This work presents new exact solutions to the unsteady three dimensional Navier-Stokes equations for incompressible viscous flows. These solutions are obtained by means of split and auto-Bäcklund transformations. The splitting procedure decouples the Navier-Stokes equations into a linear and a nonlinear inhomogeneous system of partial differential equations. The linear system, which contains only viscous terms and time derivatives, is solved via auto-Bäcklund transformations induced by commutation relations, furnishing the velocity field. The components of the velocity vector are then replaced into the nonlinear system to obtain the corresponding pressure field. The solution of the nonlinear system for the pressure variable can be carried out either numerically (by direct integration) or analytically, using the Helmholtz equation . The aim of the proposed work is to find analytical expressions for the velocity field and to obtain numerical results to the associated pressure field. The hybrid character of the solutions provides a significant reduction on the time processing required to simulate viscous flows, which virtually reduces to the time demanded to execute post-processing tasks. Taking this fact in mind, a three dimensional scalar formulation for the streamfunction was developed in order to simplify the most time-consuming post-processing task required, e.g., plotting the streamlines around arbitrary shaped bodies. At this stage of development, this formulation is employed to produce streamline maps for viscous flows around a sphere for high Reynolds numbers.

Orientador: Jose Roberto Nunhez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-05T05:52:35Z (GMT). No. of bitstreams: 1 Galeazzo_FlavioCesarCunha_M.pdf: 7104794 bytes, checksum: 95d360083b713d4c451c26b8e6276361 (MD5) Previous issue date: 2005 / Resumo: O trabalho analisou o impacto que alterações nas condições operacionais e na geometria causam no desempenho térmico de um reator de esterificação de porte industrial da empresa M&G. As análises utilizaram um modelo numérico criado com o software FLUENT. A hipótese simplificadora de se considerar somente uma fase líquida, apesar do sistema reacional real ser trifásico, foi competente em mostrar as principais características do escoamento e da transferência de calor no interior do reator, como zonas de estagnação e área de atuação dos impelidores, além da troca de calor a partir dos tubos da serpentina axial. Os resultados levam à conclusão de que o sistema de agitação possui um papel apenas secundário no desempenho térmico do reator, frente o papel preponderante dos discos separadores de fluxo / Abstract: The work analysed the impact of operating conditions and geometric changes in the thermal performance of a esterification industrial reactor of the company M&G. The analyses used a numerical model created with the aid of the software FLUENT. The simplifying hypothesis of considering only one liquid phase instead of the real system three phases was capable of showing the main characteristics of the flow and the heat transfer inside the reactor, like stagnation zones and impellers influence zones, and the heat transfer from the axial coil. The results lead to the conclusion that the agitation system have only a secundary role in the reactor thermal performance, while the role of the flow separation discs is dominant / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Made available in DSpace on 2014-10-09T12:38:39Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:05:15Z (GMT). No. of bitstreams: 1 06042.pdf: 9369626 bytes, checksum: 37950a2f878d6535f671de4a025da71c (MD5) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

O procedimento de dimensionamento de uma linha lateral de microirrigação necessita avaliar com precisão as perdas de carga distribuídas na tubulação e as perdas de carga localizadas nas inserções dos emissores com os tubos. Estas perdas localizadas podem ser significativas quando comparadas com as perdas de carga totais, devido ao grande número de emissores instalados ao longo da linha lateral. Este trabalho reporta os resultados de um experimento sobre perda de carga distribuída, fator de atrito e perda de carga localizada em conexões de emissores &#34;on-line&#34; em tubos de polietileno de pequeno diâmetro. Foram utilizados cinco tubos com diâmetros internos de 10,0 mm, 13,0 mm, 16,3 mm, 17,4 mm e 19,7 mm. O experimento foi conduzido para números de Reynolds no intervalo de 5000 a 68000, obtidos pela variação da vazão nos tubos, a uma temperatura média da água de 20 &#177; 2 &#176;C. Os resultados foram analisados e concluiu-se que o fator de atrito f da equação de Darcy-Weisbach pode ser estimado com c = 0,300 e m = 0,25. A equação de Blasius com c = 0,316 e m = 0,25 mostrou-se conservadora na estimativa do fator de atrito, porém esse fato não constitui limitação para sua utilização em projetos de microirrigação. As análises mostraram que as duas equações proporcionam estimativas de f com pequeno desvio médio (5,1%). Para um dado conjunto tubo-conexão o coeficiente de carga cinética (KL) foi praticamente independente do número de Reynolds, para R>20000, sugerindo que cada conjunto tubo-conexão pode ser caracterizado por um valor médio de KL. Para desenvolver um procedimento de estimativa de KL, a geometria da conexão entre o emissor e o tubo foi caracterizada por um índice de obstrução IO, que depende da razão (r) entre a área da seção transversal do tubo, onde o conector está localizado, e a área da seção transversal do tubo fora do conector. Uma função potência foi ajustada aos pares experimentais (IO, KL). A seleção do modelo é consistente com o fenômeno físico uma vez que KL = 0 para r = 1 (nenhuma obstrução dentro do tubo). Para 5000<R68000 a relação foi KL = 1,23 (IO)0,51 com R2 = 0,9556 e erro padrão do ajuste igual a 0,04245. As diferenças entre os valores de KL estimados e observados são normalmente distribuídas. / Microirrigation lateral design procedure needs to accuratel evaluation of both the pipe head losses and the local losses that are due to the protrusion of emitter barbs into the flow. These local losses, in fact (in relation to the high number of emitters located along the line) can become significant compared to the overall energy loss. On this paper, the results of an experimental study on the pipe head losses, friction factor and head local losses for small-diameters polyethylene pipes are reported. The experiment was carried out using a range of Reynolds number between 5000 to 68000, obtained by varying discharge at 20 &#177; 2 &#176;C water temperature, with a internal diameter pipes of 10,0 mm, 13,0 mm, 16,3 mm, 17,4 mm and 19,7 mm. According to the results analysis and experimental conditions the friction factor (f) of the Darcy-Weisbach equation can be estimated with c = 0,300 and m = 0,25. The Blasius equation (c = 0,316 and m = 0,25) gives a conservative estimative of f, although this fact is non restrictive for microirrigationsystem design. The analysis shows that both the Blasius and the adjusted equation parameters allow accurate friction factor estimate, characterized by low mean error (5,1%). For a given pipeconnection system, the fraction KL of kinetic head was practically independent of the Reynolds number, for R>20000, which suggested that each system can be characterized by the mean value of KL. To derive an estimating procedure of KL, the geometry of the connection between the emitter and the pipe was characterized by the obstruction index IO, which is dependent on the ratio (r) between the pipe cross-section area corresponding to the section in which the emitter is located, and the pipe cross-section area. A power relationship was then fitted to the experimental IO, KL data pairs. The selection form of thr relationship is consistent with the physical phenomenon since it estimates KL = 0 for r = 1 (no obstruction into the pipe). For 5000<R<68000 the relationship was KL = 1,23 (IO)0,51 with R2 = 0,9556 and standard fit error equal to 0,04245. The differences between KL observed values and the calculated ones are normally distributed.

Estudos do escoamento em vertedouros em degraus foram realizados experimentalmente por diversos autores, mas a utilização de ferramentas de computação numérica para a simulação destes casos ainda é muito escassa. Neste trabalho realizou-se um estudo da simulação numérica do escoamento sobre vertedores em degraus, utilizando-se um programa de fluidodinâmica computacional (CFD). As configurações geométricas e hidráulicas de alguns autores da literatura foram reproduzidas através das simulações em CFD. Com as simulações testaram-se as equações propostas experimentalmente ou compararam-se diretamente os resultados experimentais com os numéricos. As características dos escoamentos nos vertedouros em degraus que foram analisadas e comparadas são as seguintes: perda ou dissipação de energia, distribuição da velocidade, distribuição da concentração de ar, perfis de pressão do degrau, um estudo da resistência de atrito e um estudo da cavitação. Foram obtidas boas comparações entre os resultados experimentais e numéricos, desta forma, todas estas análises serviram, tanto para validar a utilização das ferramentas de fluidodinâmica computacional, como para caracterizar os escoamentos para uma possível proposta de melhora no desempenho desta estrutura hidráulica. Conhecendo-se o problema da cavitação, buscou-se uma solução para minimizar a ocorrência deste fenômeno no escoamento em vertedouro em degraus. Assim, foi realizado um estudo do escoamento em um aerador de fundo. Neste estudo do aerador de fundo foram caracterizados numericamente os perfis de concentração de ar, os perfis de pressão sob o salto do aerador, os perfis de velocidade e ainda foi realizado um estudo da vazão de ar que alimenta a aeração sob o jato do aerador. Após este estudo do aerador de fundo, uma nova de configuração geométrica foi proposta para o vertedouro em degrau, com a construção de um aerador de fundo no início do escoamento de forma a minimizar o risco de cavitação ou ainda aumentar as vazões máximas de serviço. / Several experimental studies about stepped spillways flow were carried out by several authors, no ever the use of numeric computation tools is still very scarce. This work study numeric simulation in stepped spillways flow was carried out using a computation fluid dynamic (CFD) software. Several hydraulics conditions from the literature were reproduced through simulations by CFD tools. Those simulation results were compared with experimental results. The characteristics of the flow in stepped spillways were analyzed and compared focusing the following aspects: energy dissipation, velocity distribution, air concentration distribution, pressure profiles in the steps, friction resistance and cavitation. The numerical simulation results had a good approximation to the experimental results. All these analyses could be used to validate the use of the computational fluid dynamic tools, to characterize the flow for a possible improvement proposal of this hydraulic structure. The cavitation problem is the great challenge in thigh velocity flows. One solution to reduce the cavitation in the stepped spillways is the construction of a bottom aerator. A numerical study of a bottom aerator was carried out, focusing out of air concentrations, pressures under jet and velocity profiles. It was studied the discharge of air under the jet. As final result, this work proposes a new geometric configuration for stepped spillways with a bottom aerator. This bottom aerator was suggested to be set in the beginning of the spillway chute in order to reduce cavitation risk or to increase the maximum discharge in spillway.

Bubbly turbulent flows occur in a variety of industrial, naval and geophysical problems. In these flows, the bubbles in the flow interact with turbulence and/or vortical structures present in the continuous phase, resulting in bubble motion and deformation, and at the same time modifying the turbulence and/or vortical structures. Despite the fact that this has been a subject of interest for some time, mechanisms of bubble break-up due to turbulence and turbulence modulation due to bubbles are not well understood. To help understand this two-way coupled problem, we study in this thesis, the interaction of single and multiple bubbles with vortical structures; the thesis being broadly divided in to three parts. In the first part, we study the interaction of a single bubble with a single vortical structure, namely a vortex ring, formed in the continuous phase (water). This may be thought of as a simplified case of the interaction of bubbles with vortical structures in any turbulent flow. We then increase the complexity and study the interaction of a single bubble with naturally occurring vortical structures present in a fully developed turbulent channel flow, and then finally to the case of a large number of bubbles injected in to a fully developed turbulent channel. The bubble motions and deformations in all three cases are directly imaged using high speed visualizations, while the flow field information is obtained using time-resolved Particle-Image Velocimetry (PIV) in the first two cases, and from pressure drop measurements within the channel in the latter case. The interaction of a single vortex ring with a bubble has been studied for a large range of vortex ring strengths, represented in terms of a Weber number (We). We find that in all cases, the bubble is first captured by the low pressure within the core of the ring, then stretched azimuthally within the core, and gradually broken up in to a number of smaller bubbles. Along with these bubble deformations, the vorticity within the core of the ring is also modified significantly due to bubble capture. In particular, at low We, we find that the core of the ring fragments as a result of the interaction resulting in a large reduction in the enstrophy of the ring and its convection speed. In the second part of the thesis, interaction of a single bubble with naturally occurring vortical structures present in a fully developed turbulent channel is studied. In this case, single bubbles of different sizes are injected either from bottom or top wall into a channel at Reynolds number of about 60,000. We study the trajectories of the single bubble, and also investigate the effect that such bubbles have on the naturally occurring vortical structures present in these flows. The injected bubble is found to have three broadly different types of bubble paths when injected from the bottom wall, which are sliding along the wall, bouncing motions and vertical escape from the vicinity of the wall. Even at the same bubble diameter Db and channel flow Re, we find that different realizations show considerable variations, with all three bubble paths being possible. PIV measurements of a bubble captured by a naturally occurring vortical structure in the flow, shows a more rapid decrease in enstrophy compared to naturally occurring structures in the absence of bubbles, as seen in the interaction of a bubble with a vortex ring. We also find that the bubble can interact with multiple vortical structures, depending on their strength and spatial distribution in the flow, resulting in a complex bouncing bubble motion. In the third part of the study, a large number of bubbles are injected in to the channel through porous plates fixed at the top and bottom channel walls. The main parameters here are the channel Re, bubble void fraction (α) and the orientation of injection. In this case, in addition to bubble visualizations, the pressure drop through the channel is measured at different vertical locations. These measurements show large vertical variations in the measured pressure drop due to the presence of bubbles. The overall drag reduction in these cases is obtained from an integral of the pressure drop variation along the vertical direction. The visualizations show a number of bubble dynamics regimes depending on the parameters, with possibilities of both increased and decreased drag compared to the reference no bubble case. From simultaneous measurements, we relate the variations in drag reduction to the different bubble dynamics regimes. We find that at the same void fraction (α), the drag reduction obtained can be very different due to changes in bubble dynamics regimes caused by changes in other parameters. Top wall injection is observed to give good drag reductions over a wide range of flow Re and α, but is seen to saturate beyond a threshold α. In contrast, the bottom wall injection case shows that drag reduction continuously increases with αat high Re. The present study shows a maximum of about 60% increase and a similar 60% reduction in wall drag over the entire range of conditions investigated.